[obiscolly50]

I'm experiencing a problem of compressibleInterFoam solver exiting when i try to restart from a non-zero time (latestTime). Starting the simulation from 0 works just fine. But say i stop at 0.1 and try to start from latestTime, the solver exits. This is the same for both series and parallel runs. Anyone has an idea of what's going on? Here's the code i'm getting:

obi@Obi:/mnt/d/linux/rootfs/opt/openfoam6/tutorials/multiphase/compressibleInterFoam/laminar/depth$ compressibleInterFoam
/*---------------------------------------------------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  6
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
Build  : 6-1a0c91b3baa8
Exec   : compressibleInterFoam
Date   : Sep 22 2018
Time   : 15:48:37
Host   : "Obi"
PID    : 447
I/O    : uncollated
Case   : /mnt/d/linux/rootfs/opt/openfoam6/tutorials/multiphase/compressibleInterFoam/laminar/depth
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0.01


PIMPLE: No convergence criteria found


PIMPLE: No corrector convergence criteria found
        Calclations will do 3 corrections

Reading field p_rgh

Reading field U

Reading/calculating face flux field phi

Constructing twoPhaseMixtureThermo

Selecting surfaceTensionModel liquidProperties
Selecting thermodynamics package
{
    type            heRhoThermo;
    mixture         pureMixture;
    properties      liquid;
    energy          sensibleInternalEnergy;
}

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

Reading thermophysical properties


Reading g

Reading hRef
Calculating field g.h

Restarting alpha
Selecting turbulence model type laminar
Selecting laminar stress model Stokes
Creating field kinetic energy K

No MRF models present

No finite volume options present
Courant Number mean: 0.11987808 max: 0.45361393

Starting time loop

Courant Number mean: 0.11987808 max: 0.45361393
Interface Courant Number mean: 0.0054768835 max: 0.41423052
deltaT = 0.00033333333
Time = 0.01033333333

PIMPLE: Iteration 1
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.48394679  Min(alpha.water) = -3.8593614e-13  Max(alpha.water) = 1
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.48400956  Min(alpha.water) = -3.6033325e-13  Max(alpha.water) = 1
smoothSolver:  Solving for T, Initial residual = 0.0024525484, Final residual = 6.7194677e-10, No Iterations 5
obi@Obi:/mnt/d/linux/rootfs/opt/openfoam6/tutorials/multiphase/compressibleInterFoam/laminar/depth$

[pbrady2013]

Hey,


Can you post a small sample to reproduce on another machine, OS, version?


Failing that, I don't use compressible much but I do inter, I tend to run my restarts in a new case folder and copy the 0.1 fields into the 0 folder of my new case and have see no weirdness like you describe.


Cheers,
-pete

[obiscolly50]

Originally Posted by pbrady2013

Hey,


Can you post a small sample to reproduce on another machine, OS, version?


Failing that, I don't use compressible much but I do inter, I tend to run my restarts in a new case folder and copy the 0.1 fields into the 0 folder of my new case and have see no weirdness like you describe.


Cheers,
-pete

Hi Pete,

Thanks for the reply. I'm just test running the depthCharge2D tutorial available under compressibleInterFoam. I haven't made any changes to the parameters, except controlDict where i specify the start and end times.

When doing the restarts as you mentioned, is it necessary to delete any files from the last time step?

[Tobi]

Hi,

Originally Posted by obiscolly50

..., is it necessary to delete any files from the last time step?

No you dont have to delete any file. That would be funny otherwise
Based on your output, I cannot see any abort nor error/warning. Probably you are doing anything wrong.

- Which FOAM version are you using
- What did you do? I guess starting from 0 to 0.1 (endTime)
- Restart with 0.1 (startTime)

[obiscolly50]

Originally Posted by Tobi

Hi,



No you dont have to delete any file. That would be funny otherwise
Based on your output, I cannot see any abort nor error/warning. Probably you are doing anything wrong.

- Which FOAM version are you using
- What did you do? I guess starting from 0 to 0.1 (endTime)
- Restart with 0.1 (startTime)

Hi Tobi,

I'm currently using OpenFoam 6. In the past, I've used v1706 and restarted simulations from the last time step with no issues.
As you said, I ran from 0 to 0.1s, but restarting with 0.1 as new startTime and 0.2 as end time takes me back to the prompt without completing a full iteration or producing any errors or warning. But running from 0 to 0.2 s produces no problem which is quite weird.

[Tobi]

Give me one moment. I will check it.

[Tobi]

Here we go. No problems:

Create time

Create mesh for time = 0.1


PIMPLE: No convergence criteria found


PIMPLE: No corrector convergence criteria found
        Calclations will do 3 corrections

Reading field p_rgh

Reading field U

Reading/calculating face flux field phi

Constructing twoPhaseMixtureThermo

Selecting surfaceTensionModel liquidProperties
Selecting thermodynamics package 
{
    type            heRhoThermo;
    mixture         pureMixture;
    properties      liquid;
    energy          sensibleInternalEnergy;
}

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

Reading thermophysical properties


Reading g

Reading hRef
Calculating field g.h

Restarting alpha
Selecting turbulence model type laminar
Selecting laminar stress model Stokes
Creating field kinetic energy K

No MRF models present

No finite volume options present
Courant Number mean: 0.10664387 max: 0.48430098

Starting time loop

Courant Number mean: 0.10664387 max: 0.48430098
Interface Courant Number mean: 0.0085487643 max: 0.47101548
deltaT = 0.00034482759
Time = 0.1003448276

PIMPLE: Iteration 1
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.4845822  Min(alpha.water) = -1.4207953e-05  Max(alpha.water) = 1
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.48454596  Min(alpha.water) = -1.3415277e-05  Max(alpha.water) = 1
smoothSolver:  Solving for T, Initial residual = 0.00011731011, Final residual = 3.5544292e-09, No Iterations 2
GAMGPCG:  Solving for p_rgh, Initial residual = 0.030081878, Final residual = 1.7224384e-08, No Iterations 5
PIMPLE: Iteration 2
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.48461754  Min(alpha.water) = -1.4333223e-05  Max(alpha.water) = 1
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.48461936  Min(alpha.water) = -1.3581991e-05  Max(alpha.water) = 1
smoothSolver:  Solving for T, Initial residual = 4.0596051e-05, Final residual = 9.4069755e-10, No Iterations 2
GAMGPCG:  Solving for p_rgh, Initial residual = 0.004244985, Final residual = 5.497686e-08, No Iterations 4
PIMPLE: Iteration 3
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.48461802  Min(alpha.water) = -1.4279122e-05  Max(alpha.water) = 1
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.48461882  Min(alpha.water) = -1.3483091e-05  Max(alpha.water) = 1
smoothSolver:  Solving for T, Initial residual = 8.1976307e-06, Final residual = 4.3512811e-10, No Iterations 2
GAMGPCG:  Solving for p_rgh, Initial residual = 0.00082954307, Final residual = 8.5665535e-08, No Iterations 3
ExecutionTime = 0.27 s


Courant Number mean: 0.10793185 max: 0.49006626
Interface Courant Number mean: 0.0086214425 max: 0.478485
deltaT = 0.00034482759
Time = 0.1006896552

- Starting from 0 to 0.1s
- StartTime = 0.1 s to endTime 0.2 s
- Output as above.

[obiscolly50]

Originally Posted by Tobi

Here we go. No problems:

- Starting from 0 to 0.1s
- StartTime = 0.1 s to endTime 0.2 s
- Output as above.

Thanks for looking into it. I'm now suspecting there might be a problem with my installation. I'm not sure reinstalling will help. What version do you use?

[Tobi]

I use the dev branch of the Foundation version. I did not have time to check weather there was a commit on this or not. However, I guess it was not a code problem. It should be related to your installation. But the funny thing is, that you get back to the prompt without any error/warning or even an output of foam. Would be interesting to see in which peace of the code it leaves the program.

You can also check other solvers such as rhoPimpleFoam. If the behavior is the same, it should be related to you installation.

[imech]

Originally Posted by obiscolly50

I'm experiencing a problem of compressibleInterFoam solver exiting when i try to restart from a non-zero time (latestTime). Starting the simulation from 0 works just fine. But say i stop at 0.1 and try to start from latestTime, the solver exits. This is the same for both series and parallel runs. Anyone has an idea of what's going on? Here's the code i'm getting:

obi@Obi:/mnt/d/linux/rootfs/opt/openfoam6/tutorials/multiphase/compressibleInterFoam/laminar/depth$ compressibleInterFoam
/*---------------------------------------------------------------------------*\
  =========                 |
  \\      /  F ield         | OpenFOAM: The Open Source CFD Toolbox
   \\    /   O peration     | Website:  https://openfoam.org
    \\  /    A nd           | Version:  6
     \\/     M anipulation  |
\*---------------------------------------------------------------------------*/
Build  : 6-1a0c91b3baa8
Exec   : compressibleInterFoam
Date   : Sep 22 2018
Time   : 15:48:37
Host   : "Obi"
PID    : 447
I/O    : uncollated
Case   : /mnt/d/linux/rootfs/opt/openfoam6/tutorials/multiphase/compressibleInterFoam/laminar/depth
nProcs : 1
sigFpe : Enabling floating point exception trapping (FOAM_SIGFPE).
fileModificationChecking : Monitoring run-time modified files using timeStampMaster (fileModificationSkew 10)
allowSystemOperations : Allowing user-supplied system call operations

// * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * //
Create time

Create mesh for time = 0.01


PIMPLE: No convergence criteria found


PIMPLE: No corrector convergence criteria found
        Calclations will do 3 corrections

Reading field p_rgh

Reading field U

Reading/calculating face flux field phi

Constructing twoPhaseMixtureThermo

Selecting surfaceTensionModel liquidProperties
Selecting thermodynamics package
{
    type            heRhoThermo;
    mixture         pureMixture;
    properties      liquid;
    energy          sensibleInternalEnergy;
}

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

Reading thermophysical properties


Reading g

Reading hRef
Calculating field g.h

Restarting alpha
Selecting turbulence model type laminar
Selecting laminar stress model Stokes
Creating field kinetic energy K

No MRF models present

No finite volume options present
Courant Number mean: 0.11987808 max: 0.45361393

Starting time loop

Courant Number mean: 0.11987808 max: 0.45361393
Interface Courant Number mean: 0.0054768835 max: 0.41423052
deltaT = 0.00033333333
Time = 0.01033333333

PIMPLE: Iteration 1
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.48394679  Min(alpha.water) = -3.8593614e-13  Max(alpha.water) = 1
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.48400956  Min(alpha.water) = -3.6033325e-13  Max(alpha.water) = 1
smoothSolver:  Solving for T, Initial residual = 0.0024525484, Final residual = 6.7194677e-10, No Iterations 5
obi@Obi:/mnt/d/linux/rootfs/opt/openfoam6/tutorials/multiphase/compressibleInterFoam/laminar/depth$

have you solved this problem?

[chaudhad]

Hi guys, I know it is an old post but I am stuck for a very long time. Actually, I am simulating a 2phase flow in a cylinder where initially air is filled and water is moved inside the cylinder with a constant velocity to compress the air. Simulations have already been done with blaockMesh utility to generate a mesh and they have come up with a result that matches the experimental results. Now we want to do some modifications in the geometry like adding some cooling pipes and then simulating the flow. Since the geometry is a bit complicated, I decided to use snappyHexMesh to mesh my geometry. but it seems like whenever I try to make my mesh fine by adding more boundary layers like 15 in this case, the simulation stops without giving any warning. I am still unable to validate the previous simulation done with blockMesh when I am using snappyHExMesh on the same geometry. Coarse Mesh gives results differently and whenever I try to make my mesh refined the simulation stops. I am using the LES turbulence model for my case. I have attached everything in the post. Please if someone can help, it will be highly appreciated. I have attached log files of the simulation as well.

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

error



Courant Number mean: 0.0014648472 max: 0.18931605
Interface Courant Number mean: 1.6251368e-05 max: 0.1644763
deltaT = 0.00028510116
Time = 0.00335071

MULES: Solving for alpha.water
Phase-1 volume fraction = 0.00012129885 Min(alpha.water) = -0.077692331 Max(alpha.water) = 1
MULES: Solving for alpha.water
Phase-1 volume fraction = 0.00012668907 Min(alpha.water) = -0.081284766 Max(alpha.water) = 1





// ************************************************** *********************** //
checkMesh

Checking geometry...
Overall domain bounding box (-0.025900001 -0.025892001 -5.74917e-18) (0.025900001 0.025892001 0.90600002)
Mesh has 3 geometric (non-empty/wedge) directions (1 1 1)
Mesh has 3 solution (non-empty) directions (1 1 1)
Boundary openness (-5.0241235e-18 -4.1763932e-17 -3.4894248e-17) OK.
Max cell openness = 8.8521418e-16 OK.
Max aspect ratio = 20.398317 OK.
Minimum face area = 2.7732229e-07. Maximum face area = 3.0123228e-05. Face area magnitudes OK.
Min volume = 6.6540004e-10. Max volume = 1.1361717e-07. Total volume = 0.0019039706. Cell volumes OK.
Mesh non-orthogonality Max: 38.107286 average: 4.9532135
Non-orthogonality check OK.
Face pyramids OK.
Max skewness = 0.43586216 OK.
Coupled point location match (average 0) OK.

Mesh OK.

End

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

0 folder

alphat

boundaryField
{
"(cylinderWalls|topX|bottomX|topY|bottomY|topZ|bot tomZ)"
{
type compressible::alphatWallFunction;
value uniform 0;
}

pipeWalls
{
type compressible::alphatWallFunction;
value uniform 0;
}

top
{
type compressible::alphatWallFunction;
value uniform 0;
}

inlet
{
type calculated;
value uniform 0;
}
}

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

nut

boundaryField
{


"(cylinderWalls|topX|bottomX|topY|bottomY|topZ|bot tomZ)"
{
type nutUWallFunction;
value uniform 0;
}

pipeWalls
{
type nutUWallFunction;
value uniform 0;
}

top
{
type nutUWallFunction;
value uniform 0;
}

inlet
{
type calculated;
value uniform 0;
}
}

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



simulation.txt

fvScheme.txt

fvsolution.txt

snappyHex.txt