[blockMesh] Velocity,Pressure locked in Blocks

Tobi94 · December 30, 2019, 16:01

Hello,

i want to do a simulation like the tutorial/incompressible/icoFoam/elbow but with my own Geometry.

The trouble by my simulation is the inlets in the System are locked in the Blocks who have the Inletpatches.
And so the hole system reaches its convergenz after s few seconds and no flow through the system.

If anyone know a possible reason, i would be thankful.

Pls ignore the commentaries in the Codes, btw, they are German.

Thankfully
Newbie Tobi

BlockMeshDict

Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v1906                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    object      blockMeshDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

scale   1.0;    // convertToMeters: 1m * scale = [Einheit]

vertices
(
    (0 0 0) 	//0  - unten, rechtes Ufer, Startpunkt
    (20 0 0)	//1  - unten, rechtes Ufer, gegenüber ende Auslasswand1 (7)
    (35 0 0)	//2  - unten, rechtes Ufer, gegenüber ende Auslasswand2 (8)
    (50 0 0)	//3  - unten, rechtes Ufer, gegenüber Buhnenkopf1 (14)
    (60 0 0)	//4  - unten, rechtes Ufer, gegenüber Buhnenkopf2 (15)
    (200 0 0)	//5  - unten, rechtes Ufer, Endpunkt
    (0 180 0)   //6  - unten, linkes Ufer, Startpunkt
    (20 180 0)	//7  - unten, linkes Ufer, ende Auslasswand1
    (35 180 0)	//8  - unten, linkes Ufer, ende Auslasswand2
    (50 180 0)	//9  - unten, linkes Ufer, Buhnenwurzel1
    (60 180 0)	//10 - unten, linkes Ufer, Buhnenwurzel2 
    (200 180 0)	//11 - unten, linkes Ufer, Endpunkt
    (20 200 0)	//12 - unten, anfang Auslasswand1
    (35 200 0)	//13 - unten, anfang Auslasswand2
    (50 120 0)	//14 - unten, Buhnenkopf1
    (60 120 0) 	//15 - unten, Buhnenkopf2

    (0 0 1) 	//16 - oben, rechtes Ufer, Startpunkt
    (20 0 1)	//17 - oben, rechtes Ufer, gegenüber ende Auslasswand1 (23)
    (35 0 1)	//18 - oben, rechtes Ufer, gegenüber ende Auslasswand2 (24)
    (50 0 1)	//19 - oben, rechtes Ufer, gegenüber Buhnenkopf1 (30)
    (60 0 1)	//20 - oben, rechtes Ufer, gegenüber Buhnenkopf2 (31)
    (200 0 1)	//21 - oben, rechtes Ufer, Endpunkt
    (0 180 1)   //22 - oben, linkes Ufer, Startpunkt
    (20 180 1)	//23 - oben, linkes Ufer, ende Auslasswand1
    (35 180 1)	//24 - oben, linkes Ufer, ende Auslasswand2
    (50 180 1)	//25 - oben, linkes Ufer, Buhnenwurzel1
    (60 180 1)	//26 - oben, linkes Ufer, Buhnenwurzel2 
    (200 180 1)	//27 - oben, linkes Ufer, Endpunkt
    (20 200 1)	//28 - oben, anfang Auslasswand1
    (35 200 1)	//29 - oben, anfang Auslasswand2
    (50 120 1)	//30 - oben, Buhnenkopf1
    (60 120 1) 	//31 - oben, Buhnenkopf2

    (35 120 0)  //32 - unten, Buhne links
    (35 120 1)  //33 - oben, Buhne links
    (200 120 0) //34 - unten, Buhne rechts
    (200 120 1) //35 - oben, Buhne rechts
    (20 120 0) 	//36 - unten, Flussmitte
    (20 120 1)	//37 - oben, Flussmitte
    (0 120 0)   //38 - unten, mitte Einlass
    (0 120 1)   //39 - oben, mitte Einlass
);

blocks
(
    hex (0 1 36 38 16 17 37 39)	(20 120 1) simpleGrading (1 1 1) //  1 Haupt
    hex (38 36 7 6 39 37 23 22) (20 60 1) simpleGrading (1 1 1)   //  2 Haupt
    hex (1 2 32 36 17 18 33 37) (15 120 1) simpleGrading (1 1 1) //  3 Haupt
    hex (36 32 8 7 37 33 24 23) (15 60 1) simpleGrading (1 1 1)   //  4 Haupt
    hex (7 8 13 12 23 24 29 28) (15 20 1) simpleGrading (1 1 1)  //  5 Auslass
    hex (2 3 14 32 18 19 30 33) (15 120 1) simpleGrading (1 1 1)  //  6 Haupt
    hex (32 14 9 8 33 30 25 24) (15 60 1) simpleGrading (1 1 1)   //  7 Haupt
    hex (3 4 15 14 19 20 31 30) (10 120 1) simpleGrading (1 1 1)  //  8 Haupt
 // hex (14 15 10 9 30 31 26 25) (10 60 1) simpleGrading (1 1 1)  //  9 Buhne
    hex (4 5 34 15 20 21 35 31) (140 120 1) simpleGrading (1 1 1) // 10 Haupt
    hex (15 34 11 10 31 35 27 26) (140 60 1) simpleGrading (1 1 1)// 11 Haupt
);

edges
(
);

boundary
(
    Wall_right
    {
        type wall;
        faces
        (
            (0 1 17 16) // rechtes Ufer, 1 Haupt
            (1 2 18 17) // rechtes Ufer, 2 Haupt
            (2 3 19 18) // rechtes Ufer, 3 Haupt
            (3 4 20 19) // rechtes Ufer, 4 Haupt
            (4 5 21 20) // rechtes Ufer, 5 Haupt
        );
    }


    Wall_left
    {
        type wall;
        faces
        (
            (7 6 22 23) // linkes Ufer, 1 Haupt
            (12 7 23 28) // Auslasswand1
            (8 13 29 24) // Auslasswand2
            (9 8 24 25) // linkes Ufer, 3 Haupt
            (14 9 25 30) // Buhnenwand1
            (15 14 30 31) // Buhnenkopf
            (10 15 31 26) // Buhnenwand2
            (11 10 26 27) // linkes Ufer, 5 Haupt
        );
    }


    inletOne
    {
        type patch;
        faces
        (
            (13 12 28 29) // Kraftwerk Auslass
        );
    }
    inletTwo
    {
        type patch;
        faces
        (
            (6 38 39 22) // Fluss, Block 1
            (38 0 16 39) // Fluss, Block 2
        );
    }
    outlet
    {
        type patch;
        faces
        (
            (5 34 35 21)  // Fluss, Block 10
            (34 11 27 35) // Fluss, Block 11
        );
    }


    frontAndBackPlanes
    {
        type empty;
        faces
        (
            (16 17 37 39) //  1 Block, oben
            (39 37 23 22) //  2 Block, oben
            (17 18 33 37) //  3 Block, oben
            (37 33 24 23) //  4 Block, oben
            (23 24 29 28) //  5 Block, oben
            (18 19 30 33) //  6 Block, oben
            (33 30 25 24) //  7 Block, oben
            (19 20 31 30) //  8 Block, oben
          //(30 31 26 25) //  9 Block, oben
            (20 21 35 31) // 10 Block, oben
            (31 35 27 26) // 11 Block, oben

            (1 0 38 36)   //  1 Block, unten
            (36 38 6 7)   //  2 Block, unten
            (2 1 36 32)   //  3 Block, unten
            (32 36 7 8)   //  4 Block, unten
            (8 7 12 13)   //  5 Block, unten
            (3 2 32 14)   //  6 Block, unten
            (14 32 8 9)   //  7 Block, unten
            (4 3 14 15)   //  8 Block, unten
          //(15 14 9 10)  //  9 Block, unten
            (5 4 15 34)   // 10 Block, unten
            (34 15 10 11) // 11 Block, unten
        );
    }

);

mergePatchPairs
(
);

// ************************************************************************* //

controlDict

Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v1906                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      controlDict;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

application     icoFoam;

startFrom       latestTime;

startTime       0;

stopAt          endTime;

endTime         250;

deltaT          1.0;

writeControl    timeStep;

writeInterval   1;

purgeWrite      0;

writeFormat     ascii;

writePrecision  6;

writeCompression off;

timeFormat      general;

timePrecision   6;

runTimeModifiable true;


// ************************************************************************* //

p

Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v1906                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volScalarField;
    object      p;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 2 -2 0 0 0 0];

internalField   uniform 0;

boundaryField
{
    Wall_right
    {
        type            zeroGradient;
    }
    Wall_left
    {
        type            zeroGradient;
    }
    inletOne
    {
        type            zeroGradient;
    }

    inletTwo
    {
        type            zeroGradient;
    }

    outlet
    {
        type            fixedValue;
        value           uniform 0;
    }

    frontAndBackPlanes
    {
        type            empty;
    }
}

// ************************************************************************* //

U

Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v1906                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       volVectorField;
    object      U;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

dimensions      [0 1 -1 0 0 0 0];

internalField   uniform (0 0 0);

boundaryField
{
    Wall_right
    {
        type            noSlip;
    }
    Wall_left
    {
        type            noSlip;
    }

    inletOne
    {
        type            fixedValue;
        value           uniform (6 0 0);
    }

    inletTwo
    {
        type            fixedValue;
        value           uniform (0 -10 0);
    }

    outlet
    {
        type            zeroGradient;
    }



    frontAndBackPlanes
    {
        type            empty;
    }
}

// ************************************************************************* //

fvSchemes

Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v1906                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      fvSchemes;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

ddtSchemes
{
    default         Euler;
}

gradSchemes
{
    default         Gauss linear;
}

divSchemes
{
    default         none;
    div(phi,U)      Gauss limitedLinearV 1;
}

laplacianSchemes
{
    default         Gauss linear corrected;
}

interpolationSchemes
{
    default         linear;
}

snGradSchemes
{
    default         corrected;
}


// ************************************************************************* //

fvSolution

Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v1906                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "system";
    object      fvSolution;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

solvers
{
    p
    {
        solver          PCG;
        preconditioner  DIC;
        tolerance       1e-06;
        relTol          0.05;
    }

    pFinal
    {
        $p;
        relTol          0;
    }

    U
    {
        solver          smoothSolver;
        smoother        symGaussSeidel;
        tolerance       1e-05;
        relTol          0;
    }
}

PISO
{
    nCorrectors     3;
    nNonOrthogonalCorrectors 3;
}


// ************************************************************************* //

transportProperties

Code:

/*--------------------------------*- C++ -*----------------------------------*\
| =========                 |                                                 |
| \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox           |
|  \\    /   O peration     | Version:  v1906                                 |
|   \\  /    A nd           | Web:      www.OpenFOAM.com                      |
|    \\/     M anipulation  |                                                 |
\*---------------------------------------------------------------------------*/
FoamFile
{
    version     2.0;
    format      ascii;
    class       dictionary;
    location    "constant";
    object      transportProperties;
}
// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //

nu              0.01;


// ************************************************************************* //

HPE · December 30, 2019, 18:17

Hi,

Is there any reason why 'inlet2' has -10 m/s inflow speed?

Tobi94 · December 31, 2019, 09:27

Hi HPE,

with a low inlet like -0.5 the flow don't leave the first Block so i was trying what happends, if i make the flow faster.

But also, sadly i switch the inlets between blockMeshDic and U. I corrected U but the Problem ist still subsits.

With a low inlet flows the convergenz is reached after 1sec:

Code:

    inletOne
    {
        type            fixedValue;
        value           uniform (0 -1 0);
    }

    inletTwo
    {
        type            fixedValue;
        value           uniform (2 0 0);
    }

December 30, 2019, 18:17		#2
HPE Senior Member Herpes Free Engineer Join Date: Sep 2019 Location: The Home Under The Ground with the Lost Boys Posts: 931 Rep Power: 13	Hi, Is there any reason why 'inlet2' has -10 m/s inflow speed? __________________ The OpenFOAM community is the biggest contributor to OpenFOAM: User guide/Wiki-1/Wiki-2/Code guide/Code Wiki/Journal Nilsson/Guerrero/Holzinger/Holzmann/Nagy/Santos/Nozaki/Jasak/Primer Governance Bugs/Features: OpenFOAM (ESI-OpenCFD-Trademark) Bugs/Features: FOAM-Extend (Wikki-FSB) Bugs: OpenFOAM.org How to create a MWE New: Forkable OpenFOAM mirror

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
[Commercial meshers] fluent3DMeshToFoam	bego	OpenFOAM Meshing & Mesh Conversion	31	August 16, 2023 10:04
[Commercial meshers] COnvert FLuent MEsh to openfoam with interface	manuc	OpenFOAM Meshing & Mesh Conversion	1	July 25, 2017 04:13
[Commercial meshers] converting Fluent mesh to openfoam standard mesh	deepesh	OpenFOAM Meshing & Mesh Conversion	31	March 29, 2017 06:59
dsmcInitialise - dsmcFoam	archymedes	OpenFOAM Pre-Processing	94	July 15, 2016 17:14
[Other] How to create an MRF zone ?	aminem	OpenFOAM Meshing & Mesh Conversion	2	December 8, 2014 11:45

December 31, 2019, 09:27		#3
Tobi94 New Member Tobias R Join Date: Dec 2019 Posts: 3 Rep Power: 7	Hi HPE, with a low inlet like -0.5 the flow don't leave the first Block so i was trying what happends, if i make the flow faster. But also, sadly i switch the inlets between blockMeshDic and U. I corrected U but the Problem ist still subsits. With a low inlet flows the convergenz is reached after 1sec: Code: inletOne { type fixedValue; value uniform (0 -1 0); } inletTwo { type fixedValue; value uniform (2 0 0); }