|
[Sponsors] |
compressibleInterFoam exits without warning or error |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
September 22, 2018, 18:14 |
compressibleInterFoam exits without warning or error
|
#1 |
Member
Obi
Join Date: Jul 2016
Location: Canada
Posts: 45
Rep Power: 10 |
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$ |
|
September 23, 2018, 20:44 |
Workaround and Sample
|
#2 |
Member
Peter Brady
Join Date: Apr 2014
Location: Sydney, NSW, Australia
Posts: 54
Rep Power: 12 |
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 |
|
September 23, 2018, 20:58 |
|
#3 | |
Member
Obi
Join Date: Jul 2016
Location: Canada
Posts: 45
Rep Power: 10 |
Quote:
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? |
||
September 24, 2018, 05:27 |
|
#4 | |
Super Moderator
Tobias Holzmann
Join Date: Oct 2010
Location: Bad Wörishofen
Posts: 2,711
Blog Entries: 6
Rep Power: 52 |
Hi,
Quote:
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 |
||
September 24, 2018, 10:57 |
|
#5 | |
Member
Obi
Join Date: Jul 2016
Location: Canada
Posts: 45
Rep Power: 10 |
Quote:
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. |
||
September 24, 2018, 12:42 |
|
#7 |
Super Moderator
Tobias Holzmann
Join Date: Oct 2010
Location: Bad Wörishofen
Posts: 2,711
Blog Entries: 6
Rep Power: 52 |
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 |
|
September 24, 2018, 12:59 |
|
#8 |
Member
Obi
Join Date: Jul 2016
Location: Canada
Posts: 45
Rep Power: 10 |
||
September 24, 2018, 13:13 |
|
#9 |
Super Moderator
Tobias Holzmann
Join Date: Oct 2010
Location: Bad Wörishofen
Posts: 2,711
Blog Entries: 6
Rep Power: 52 |
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 |
|
September 11, 2020, 11:10 |
|
#10 | |
New Member
lpp
Join Date: Sep 2020
Posts: 2
Rep Power: 0 |
Quote:
|
||
June 16, 2023, 05:37 |
CompressibleInterFoam stopping without giving any error
|
#11 |
New Member
Shanu
Join Date: Jun 2023
Posts: 3
Rep Power: 3 |
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 |
|
Tags |
compressibleinterfoam |
Thread Tools | Search this Thread |
Display Modes | |
|
|
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 |