CFD Online Discussion Forums

CFD Online Discussion Forums (https://www.cfd-online.com/Forums/)
-   OpenFOAM Running, Solving & CFD (https://www.cfd-online.com/Forums/openfoam-solving/)
-   -   interFoam not proceeding after certain time (https://www.cfd-online.com/Forums/openfoam-solving/224469-interfoam-not-proceeding-after-certain-time.html)

joshmccraney February 20, 2020 16:20

interFoam not proceeding after certain time
 
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

The log output at the end (when it freezes) is:

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


joshmccraney February 20, 2020 18:09

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.

HPE February 21, 2020 01:09

Is it possible to replicate the issue with a tutorial?

Or is it possible to attach the case causing the issue?

joshmccraney February 21, 2020 10:25

Quote:

Originally Posted by HPE (Post 758989)
Is it possible to replicate the issue with a tutorial?

Or is it possible to attach the case causing the issue?

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.


All times are GMT -4. The time now is 10:32.