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

interFoam not proceeding after certain time

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

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   February 20, 2020, 17:20
Default interFoam not proceeding after certain time
  #1
Senior Member
 
Josh McCraney
Join Date: Jun 2018
Posts: 220
Rep Power: 8
joshmccraney is on a distinguished road
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

Last edited by joshmccraney; February 20, 2020 at 18:46.
joshmccraney is offline   Reply With Quote

Old   February 20, 2020, 19:09
Default
  #2
Senior Member
 
Josh McCraney
Join Date: Jun 2018
Posts: 220
Rep Power: 8
joshmccraney is on a distinguished road
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.
joshmccraney is offline   Reply With Quote

Old   February 21, 2020, 02:09
Default
  #3
HPE
Senior Member
 
HPE's Avatar
 
Herpes Free Engineer
Join Date: Sep 2019
Location: The Home Under The Ground with the Lost Boys
Posts: 932
Rep Power: 12
HPE is on a distinguished road
Is it possible to replicate the issue with a tutorial?

Or is it possible to attach the case causing the issue?
HPE is offline   Reply With Quote

Old   February 21, 2020, 11:25
Default
  #4
Senior Member
 
Josh McCraney
Join Date: Jun 2018
Posts: 220
Rep Power: 8
joshmccraney is on a distinguished road
Quote:
Originally Posted by HPE View Post
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.
joshmccraney is offline   Reply With Quote

Reply

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
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


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