|
[Sponsors] |
February 20, 2020, 17:20 |
interFoam not proceeding after certain time
|
#1 |
Senior Member
Josh McCraney
Join Date: Jun 2018
Posts: 220
Rep Power: 8 |
Hi FOAMers!
I'm simulating a 1m bubble in space. At 163 seconds (about 43 seconds of clocktime) the simulation does not proceed. It is still working (I see output logs on the terminal screen) but it will not proceed beyond the time of 163 seconds. I let the clocktime go to 75 seconds, but still nothing. This is odd, because if I change the initial conditions a little, I can get simulations for 1000 seconds with clocktimes at only a few minutes. Any idea what is happening? The log output early on (when it runs well) is: Code:
smoothSolver: Solving for alpha.water, Initial residual = 0.000270691, Final residual = 5.43977e-11, No Iterations 2 Phase-1 volume fraction = 0.3852 Min(alpha.water) = 0 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 0.3852 Min(alpha.water) = -1.10156e-79 Max(alpha.water) = 1 DICPCG: Solving for p_rgh, Initial residual = 0.00525566, Final residual = 0.000171449, No Iterations 3 time step continuity errors : sum local = 2.74885e-06, global = 6.78163e-07, cumulative = 5.46983e-05 DICPCG: Solving for p_rgh, Initial residual = 0.00101003, Final residual = 4.74852e-05, No Iterations 49 time step continuity errors : sum local = 7.61139e-07, global = 1.35427e-07, cumulative = 5.48337e-05 DICPCG: Solving for p_rgh, Initial residual = 0.000129676, Final residual = 8.93693e-08, No Iterations 96 time step continuity errors : sum local = 1.41385e-09, global = 7.0341e-11, cumulative = 5.48338e-05 ExecutionTime = 5.9 s ClockTime = 6 s Courant Number mean: 0.0220595 max: 0.184886 Interface Courant Number mean: 0.00260337 max: 0.107418 deltaT = 0.102041 Time = 24.8776 Code:
Courant Number mean: 0.0119459 max: 0.210334 Interface Courant Number mean: 0.000120731 max: 0.00966653 deltaT = 7.99842e-77 Time = 292.591 smoothSolver: Solving for alpha.water, Initial residual = 2.16354e-06, Final residual = 4.96374e-10, No Iterations 1 Phase-1 volume fraction = 0.384993 Min(alpha.water) = -1.23499e-106 Max(alpha.water) = 1 MULES: Correcting alpha.water MULES: Correcting alpha.water Phase-1 volume fraction = 0.384993 Min(alpha.water) = -1.13528e-11 Max(alpha.water) = 1 DICPCG: Solving for p_rgh, Initial residual = 0.00155701, Final residual = 7.37893e-05, No Iterations 71 time step continuity errors : sum local = 2.35383e-07, global = 2.0591e-10, cumulative = 0.000198153 DICPCG: Solving for p_rgh, Initial residual = 0.00247107, Final residual = 0.000107513, No Iterations 70 time step continuity errors : sum local = 1.39834e-07, global = -3.65717e-09, cumulative = 0.000198149 DICPCG: Solving for p_rgh, Initial residual = 0.00122908, Final residual = 8.98552e-08, No Iterations 105 time step continuity errors : sum local = 2.12212e-10, global = -6.87585e-12, cumulative = 0.000198149 ExecutionTime = 133.33 s ClockTime = 150 s Last edited by joshmccraney; February 20, 2020 at 18:46. |
|
February 20, 2020, 19:09 |
|
#2 |
Senior Member
Josh McCraney
Join Date: Jun 2018
Posts: 220
Rep Power: 8 |
The test case runs much better (still freezes) when changing the file fvSolution, the p_rgh had relTol set to 0.05. When I changed this to 0, the case ran well.
Have no idea why though. |
|
February 21, 2020, 02:09 |
|
#3 |
Senior Member
Herpes Free Engineer
Join Date: Sep 2019
Location: The Home Under The Ground with the Lost Boys
Posts: 932
Rep Power: 12 |
Is it possible to replicate the issue with a tutorial?
Or is it possible to attach the case causing the issue?
__________________
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 |
|
February 21, 2020, 11:25 |
|
#4 |
Senior Member
Josh McCraney
Join Date: Jun 2018
Posts: 220
Rep Power: 8 |
It seems to be working now. I changed the tolerances in fvSolution. Nothing serious, though. I think the big change was the relTol in my second post. At least for future people who have this "freezing" issue, this is how I fixed it.
|
|
Thread Tools | Search this Thread |
Display Modes | |
|
|
Similar Threads | ||||
Thread | Thread Starter | Forum | Replies | Last Post |
negative min.(alpha1) value interFoam | Arjun Jayakumar | OpenFOAM | 11 | December 21, 2019 11:59 |
laplacianFoam with source term | Herwig | OpenFOAM Running, Solving & CFD | 17 | November 19, 2019 14:47 |
Inconsistencies in reading .dat file during run time in new injection model | Scram_1 | OpenFOAM | 0 | March 23, 2018 23:29 |
pressure in incompressible solvers e.g. simpleFoam | chrizzl | OpenFOAM Running, Solving & CFD | 13 | March 28, 2017 06:49 |
InterFoam negative alpha | karasa03 | OpenFOAM | 7 | December 12, 2013 04:41 |