CFD Online Logo CFD Online URL
www.cfd-online.com
[Sponsors]
Home > Forums > Software User Forums > OpenFOAM > OpenFOAM Running, Solving & CFD

compressibleInterFoam exits without warning or error

Register Blogs Members List Search Today's Posts Mark Forums Read

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   September 22, 2018, 18:14
Default compressibleInterFoam exits without warning or error
  #1
Member
 
Obi
Join Date: Jul 2016
Location: Canada
Posts: 45
Rep Power: 10
obiscolly50 is on a distinguished road
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:


Code:
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$
obiscolly50 is offline   Reply With Quote

Old   September 23, 2018, 20:44
Default Workaround and Sample
  #2
Member
 
Peter Brady
Join Date: Apr 2014
Location: Sydney, NSW, Australia
Posts: 54
Rep Power: 12
pbrady2013 is on a distinguished road
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
pbrady2013 is offline   Reply With Quote

Old   September 23, 2018, 20:58
Default
  #3
Member
 
Obi
Join Date: Jul 2016
Location: Canada
Posts: 45
Rep Power: 10
obiscolly50 is on a distinguished road
Quote:
Originally Posted by pbrady2013 View Post
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?
obiscolly50 is offline   Reply With Quote

Old   September 24, 2018, 05:27
Default
  #4
Super Moderator
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Bad Wörishofen
Posts: 2,711
Blog Entries: 6
Rep Power: 52
Tobi has a spectacular aura aboutTobi has a spectacular aura aboutTobi has a spectacular aura about
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Hi,

Quote:
Originally Posted by obiscolly50 View Post
..., 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)
__________________
Keep foaming,
Tobias Holzmann
Tobi is offline   Reply With Quote

Old   September 24, 2018, 10:57
Default
  #5
Member
 
Obi
Join Date: Jul 2016
Location: Canada
Posts: 45
Rep Power: 10
obiscolly50 is on a distinguished road
Quote:
Originally Posted by Tobi View Post
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.
obiscolly50 is offline   Reply With Quote

Old   September 24, 2018, 12:39
Default
  #6
Super Moderator
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Bad Wörishofen
Posts: 2,711
Blog Entries: 6
Rep Power: 52
Tobi has a spectacular aura aboutTobi has a spectacular aura aboutTobi has a spectacular aura about
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Give me one moment. I will check it.
__________________
Keep foaming,
Tobias Holzmann
Tobi is offline   Reply With Quote

Old   September 24, 2018, 12:42
Default
  #7
Super Moderator
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Bad Wörishofen
Posts: 2,711
Blog Entries: 6
Rep Power: 52
Tobi has a spectacular aura aboutTobi has a spectacular aura aboutTobi has a spectacular aura about
Send a message via ICQ to Tobi Send a message via Skype™ to Tobi
Here we go. No problems:
Code:
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.
__________________
Keep foaming,
Tobias Holzmann
Tobi is offline   Reply With Quote

Old   September 24, 2018, 12:59
Default
  #8
Member
 
Obi
Join Date: Jul 2016
Location: Canada
Posts: 45
Rep Power: 10
obiscolly50 is on a distinguished road
Quote:
Originally Posted by Tobi View Post
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?
obiscolly50 is offline   Reply With Quote

Old   September 24, 2018, 13:13
Default
  #9
Super Moderator
 
Tobi's Avatar
 
Tobias Holzmann
Join Date: Oct 2010
Location: Bad Wörishofen
Posts: 2,711
Blog Entries: 6
Rep Power: 52
Tobi has a spectacular aura aboutTobi has a spectacular aura aboutTobi has a spectacular aura about
Send a message via ICQ to Tobi Send a message via Skype™ to 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.
__________________
Keep foaming,
Tobias Holzmann
Tobi is offline   Reply With Quote

Old   September 11, 2020, 11:10
Default
  #10
New Member
 
lpp
Join Date: Sep 2020
Posts: 2
Rep Power: 0
imech is on a distinguished road
Quote:
Originally Posted by obiscolly50 View Post
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:


Code:
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?
imech is offline   Reply With Quote

Old   June 16, 2023, 05:37
Default CompressibleInterFoam stopping without giving any error
  #11
New Member
 
Shanu
Join Date: Jun 2023
Posts: 3
Rep Power: 3
chaudhad is on a distinguished road
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
chaudhad is offline   Reply With Quote

Reply

Tags
compressibleinterfoam

Thread Tools Search this Thread
Search this Thread:

Advanced Search
Display Modes

Posting Rules
You may not post new threads
You may not post replies
You may not post attachments
You may not edit your posts

BB code is On
Smilies are On
[IMG] code is On
HTML code is Off
Trackbacks are Off
Pingbacks are On
Refbacks are On


Similar Threads
Thread Thread Starter Forum Replies Last Post
Compile calcMassFlowC aurore OpenFOAM Programming & Development 13 March 23, 2018 08:43
error compiling modified applications yvyan OpenFOAM Programming & Development 21 March 1, 2016 05:53
Ansys Fluent 13.0 UDF compilation problem in Window XP (32 bit) Yogini Fluent UDF and Scheme Programming 7 October 3, 2012 08:24
Version 15 on Mac OS X gschaider OpenFOAM Installation 113 December 2, 2009 11:23
user subroutine error CFDUSER CFX 2 December 9, 2006 07:31


All times are GMT -4. The time now is 17:14.