Simulation stops at Starting time loop

JMaverick93 · March 11, 2018, 19:13

Hey guys.

So, I'm a begginer at OpenFoam and I'm trying simulate the flow in a de Laval nozzle but when I try to use rhoCentralFoam it just stops at the message of starting time loop, like so:

Code:

Create time

Create mesh for time = 0

Reading thermophysical properties

Selecting thermodynamics package
{
    type            hePsiThermo;
    mixture         pureMixture;
    transport       const;
    thermo          hConst;
    equationOfState perfectGas;
    specie          specie;
    energy          sensibleInternalEnergy;
}

Reading field U

Creating turbulence model

Selecting turbulence model type laminar
Selecting laminar stress model Stokes
fluxScheme: Kurganov

Starting time loop

What could be causing it and how can I fix it?

Code:

BlockMeshDict
convertToMeters    0.001;

vertices
    (
        (0        0       -0.5)    // 0
        (20.96    0       -0.5)    // 1
        (86.82    0       -0.5)    // 2
        (0        18.1    -0.5)    // 3
        (20.96    6       -0.5)    // 4
        (86.82    20      -0.5)    // 5
        (0        0        0.5)    // 6
        (20.96    0        0.5)    // 7
        (86.82    0        0.5)    // 8
        (0        18.1     0.5)    // 9
        (20.96    6        0.5)    // 10
        (86.82    20       0.5)    // 11
    );

edges
    (
    );

blocks
    (
        hex (7 6 9 10 1 0 3 4)     // 0
        (10 10 1)                  
        simpleGrading (1 1 1)
        hex (8 7 10 11 2 1 4 5)     // 1
        (10 10 1)                  
        simpleGrading (1 1 1)
    );

boundary
    (
        inlet
            {
                type patch;
                faces
                (
                    (6 0 3 9)
                );
            }
        bottom
             {
                 type symmetryPlane;
                 faces
                 (
                     (6 0 1 7)
                     (7 1 2 8)
                 );
             }
         outlet
            {
                type patch;
                faces
                (
                    (8 2 5 11)
                );
            }
         top
             {
                 type wall;
                 faces
                 (
                     (3 9 10 4)
                     (4 10 11 5)
                 );
             }
    );

mergePatchPairs
(
);

Code:

ControlDict
application     rhoCentralFoam;

startFrom       latestTime;

startTime       0;

stopAt          endTime;

endTime         1;

deltaT          0.00000001;

writeControl    adjustableRunTime;

writeInterval   0.001;

purgeWrite      0;

writeFormat     ascii;

writePrecision  6;

writeCompression off;

timeFormat      general;

timePrecision   6;

runTimeModifiable true;

adjustTimeStep  yes;

maxCo           0.1;

maxDeltaT       1;

Code:

p
dimensions      [1 -1 -2 0 0 0 0];

internalField   uniform 6895000;

boundaryField
{
    inlet
    {
        type            totalPressure;
        p0              uniform 6895000;
        gamma           1.133;
        value           $internalField;
    }

    outlet
    {
        type            fixedValue;
        value           uniform 101325;
    }

    bottom
    {
        type            symmetryPlane;
    }

    top
    {
        type            zeroGradient;
    }

    defaultFaces
    {
        type            empty;
    }
}

Code:

T
dimensions      [0 0 0 1 0 0 0];

internalField   uniform 1722;

boundaryField
{
    inlet
    {
        type            fixedValue;
        value           uniform 1722;
    }

    outlet
    {
        type            zeroGradient;
    }

    bottom
    {
        type            symmetryPlane;
    }

    top
    {
        type            zeroGradient;
    }

    defaultFaces
    {
        type            empty;
    }
}

Code:

U
Dimensions      [0 1 -1 0 0 0 0];

internalField   uniform (0 0 0);

boundaryField
{
    inlet
    {
        type            pressureInletVelocity;
        value           uniform (0 0 0);
    }

    outlet
    {
        type            zeroGradient;
    }

    bottom
    {
        type            symmetryPlane;
    }

    top
    {
        type            zeroGradient;
    }

    defaultFaces
    {
        type            empty;
    }
}

Thanks in advance!

shereez234 · March 12, 2018, 10:16

Quote:

Originally Posted by JMaverick93

Hey guys.

So, I'm a begginer at OpenFoam and I'm trying simulate the flow in a de Laval nozzle but when I try to use rhoCentralFoam it just stops at the message of starting time loop, like so:

Code:

Create time

Create mesh for time = 0

Reading thermophysical properties

Selecting thermodynamics package
{
    type            hePsiThermo;
    mixture         pureMixture;
    transport       const;
    thermo          hConst;
    equationOfState perfectGas;
    specie          specie;
    energy          sensibleInternalEnergy;
}

Reading field U

Creating turbulence model

Selecting turbulence model type laminar
Selecting laminar stress model Stokes
fluxScheme: Kurganov

Starting time loop

What could be causing it and how can I fix it?

Code:

BlockMeshDict
convertToMeters    0.001;

vertices
    (
        (0        0       -0.5)    // 0
        (20.96    0       -0.5)    // 1
        (86.82    0       -0.5)    // 2
        (0        18.1    -0.5)    // 3
        (20.96    6       -0.5)    // 4
        (86.82    20      -0.5)    // 5
        (0        0        0.5)    // 6
        (20.96    0        0.5)    // 7
        (86.82    0        0.5)    // 8
        (0        18.1     0.5)    // 9
        (20.96    6        0.5)    // 10
        (86.82    20       0.5)    // 11
    );

edges
    (
    );

blocks
    (
        hex (7 6 9 10 1 0 3 4)     // 0
        (10 10 1)                  
        simpleGrading (1 1 1)
        hex (8 7 10 11 2 1 4 5)     // 1
        (10 10 1)                  
        simpleGrading (1 1 1)
    );

boundary
    (
        inlet
            {
                type patch;
                faces
                (
                    (6 0 3 9)
                );
            }
        bottom
             {
                 type symmetryPlane;
                 faces
                 (
                     (6 0 1 7)
                     (7 1 2 8)
                 );
             }
         outlet
            {
                type patch;
                faces
                (
                    (8 2 5 11)
                );
            }
         top
             {
                 type wall;
                 faces
                 (
                     (3 9 10 4)
                     (4 10 11 5)
                 );
             }
    );

mergePatchPairs
(
);

Code:

ControlDict
application     rhoCentralFoam;

startFrom       latestTime;

startTime       0;

stopAt          endTime;

endTime         1;

deltaT          0.00000001;

writeControl    adjustableRunTime;

writeInterval   0.001;

purgeWrite      0;

writeFormat     ascii;

writePrecision  6;

writeCompression off;

timeFormat      general;

timePrecision   6;

runTimeModifiable true;

adjustTimeStep  yes;

maxCo           0.1;

maxDeltaT       1;

Code:

p
dimensions      [1 -1 -2 0 0 0 0];

internalField   uniform 6895000;

boundaryField
{
    inlet
    {
        type            totalPressure;
        p0              uniform 6895000;
        gamma           1.133;
        value           $internalField;
    }

    outlet
    {
        type            fixedValue;
        value           uniform 101325;
    }

    bottom
    {
        type            symmetryPlane;
    }

    top
    {
        type            zeroGradient;
    }

    defaultFaces
    {
        type            empty;
    }
}

Code:

T
dimensions      [0 0 0 1 0 0 0];

internalField   uniform 1722;

boundaryField
{
    inlet
    {
        type            fixedValue;
        value           uniform 1722;
    }

    outlet
    {
        type            zeroGradient;
    }

    bottom
    {
        type            symmetryPlane;
    }

    top
    {
        type            zeroGradient;
    }

    defaultFaces
    {
        type            empty;
    }
}

Code:

U
Dimensions      [0 1 -1 0 0 0 0];

internalField   uniform (0 0 0);

boundaryField
{
    inlet
    {
        type            pressureInletVelocity;
        value           uniform (0 0 0);
    }

    outlet
    {
        type            zeroGradient;
    }

    bottom
    {
        type            symmetryPlane;
    }

    top
    {
        type            zeroGradient;
    }

    defaultFaces
    {
        type            empty;
    }
}

Thanks in advance!

maxDeltaT = 1 (change)

JMaverick93 · March 12, 2018, 10:56

Hello shereez234, thanks for replying!

Change to what exactly? I've already tried to run with a differente value (0.001) but the results were the same.

harsha_kulkarni · March 18, 2018, 07:08

Quote:

Originally Posted by JMaverick93

Hey guys.

So, I'm a begginer at OpenFoam and I'm trying simulate the flow in a de Laval nozzle but when I try to use rhoCentralFoam it just stops at the message of starting time loop, like so:

Code:

Create time

Create mesh for time = 0

Reading thermophysical properties

Selecting thermodynamics package
{
    type            hePsiThermo;
    mixture         pureMixture;
    transport       const;
    thermo          hConst;
    equationOfState perfectGas;
    specie          specie;
    energy          sensibleInternalEnergy;
}

Reading field U

Creating turbulence model

Selecting turbulence model type laminar
Selecting laminar stress model Stokes
fluxScheme: Kurganov

Starting time loop

What could be causing it and how can I fix it?

Code:

BlockMeshDict
convertToMeters    0.001;

vertices
    (
        (0        0       -0.5)    // 0
        (20.96    0       -0.5)    // 1
        (86.82    0       -0.5)    // 2
        (0        18.1    -0.5)    // 3
        (20.96    6       -0.5)    // 4
        (86.82    20      -0.5)    // 5
        (0        0        0.5)    // 6
        (20.96    0        0.5)    // 7
        (86.82    0        0.5)    // 8
        (0        18.1     0.5)    // 9
        (20.96    6        0.5)    // 10
        (86.82    20       0.5)    // 11
    );

edges
    (
    );

blocks
    (
        hex (7 6 9 10 1 0 3 4)     // 0
        (10 10 1)                  
        simpleGrading (1 1 1)
        hex (8 7 10 11 2 1 4 5)     // 1
        (10 10 1)                  
        simpleGrading (1 1 1)
    );

boundary
    (
        inlet
            {
                type patch;
                faces
                (
                    (6 0 3 9)
                );
            }
        bottom
             {
                 type symmetryPlane;
                 faces
                 (
                     (6 0 1 7)
                     (7 1 2 8)
                 );
             }
         outlet
            {
                type patch;
                faces
                (
                    (8 2 5 11)
                );
            }
         top
             {
                 type wall;
                 faces
                 (
                     (3 9 10 4)
                     (4 10 11 5)
                 );
             }
    );

mergePatchPairs
(
);

Code:

ControlDict
application     rhoCentralFoam;

startFrom       latestTime;

startTime       0;

stopAt          endTime;

endTime         1;

deltaT          0.00000001;

writeControl    adjustableRunTime;

writeInterval   0.001;

purgeWrite      0;

writeFormat     ascii;

writePrecision  6;

writeCompression off;

timeFormat      general;

timePrecision   6;

runTimeModifiable true;

adjustTimeStep  yes;

maxCo           0.1;

maxDeltaT       1;

Code:

p
dimensions      [1 -1 -2 0 0 0 0];

internalField   uniform 6895000;

boundaryField
{
    inlet
    {
        type            totalPressure;
        p0              uniform 6895000;
        gamma           1.133;
        value           $internalField;
    }

    outlet
    {
        type            fixedValue;
        value           uniform 101325;
    }

    bottom
    {
        type            symmetryPlane;
    }

    top
    {
        type            zeroGradient;
    }

    defaultFaces
    {
        type            empty;
    }
}

Code:

T
dimensions      [0 0 0 1 0 0 0];

internalField   uniform 1722;

boundaryField
{
    inlet
    {
        type            fixedValue;
        value           uniform 1722;
    }

    outlet
    {
        type            zeroGradient;
    }

    bottom
    {
        type            symmetryPlane;
    }

    top
    {
        type            zeroGradient;
    }

    defaultFaces
    {
        type            empty;
    }
}

Code:

U
Dimensions      [0 1 -1 0 0 0 0];

internalField   uniform (0 0 0);

boundaryField
{
    inlet
    {
        type            pressureInletVelocity;
        value           uniform (0 0 0);
    }

    outlet
    {
        type            zeroGradient;
    }

    bottom
    {
        type            symmetryPlane;
    }

    top
    {
        type            zeroGradient;
    }

    defaultFaces
    {
        type            empty;
    }
}

Thanks in advance!

This case is working on my machine

Similar Threads
Thread	Thread Starter	Forum	Replies	Last Post
[Other] refineWallLayer Error	Yuby	OpenFOAM Meshing & Mesh Conversion	2	November 11, 2021 12:04
pressure in incompressible solvers e.g. simpleFoam	chrizzl	OpenFOAM Running, Solving & CFD	13	March 28, 2017 06:49
Extrusion with OpenFoam problem No. Iterations 0	Lord Kelvin	OpenFOAM Running, Solving & CFD	8	March 28, 2016 12:08
Stuck in a Rut- interDyMFoam!	xoitx	OpenFOAM Running, Solving & CFD	14	March 25, 2016 08:09
simpleFoam error - "Floating point exception"	mbcx4jc2	OpenFOAM Running, Solving & CFD	12	August 4, 2015 03:20

March 12, 2018, 10:56		#3
JMaverick93 New Member Jonathan Maverick Join Date: Mar 2018 Posts: 2 Rep Power: 0	Hello shereez234, thanks for replying! Change to what exactly? I've already tried to run with a differente value (0.001) but the results were the same.