|
[Sponsors] |
tutorial buoyantCavity - derived 2D case - convergence |
![]() |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
![]() |
![]() |
#1 |
New Member
Herbert Dobernig
Join Date: Apr 2020
Posts: 4
Rep Power: 6 ![]() |
I derived a 2D case from the tutorial buoyantCavity by changing boundary type of patch "frontAndBack" to "empty" and reducing grid to 1 in z-direction.
I set "endTime" to 10000 iterations buoyantSimpleFoam stops after 3350 iterations: log.buoyantSimpleFoam: "SIMPLE solution converged in 3350 iterations" After the 3350 itartions heat transfer has not balanced to a sufficient degree and thus temerature ist still rising significantly. What are recommended ways to prevent "buoyantSimpleFoam" from stopping after 3350 iterations? |
|
![]() |
![]() |
![]() |
![]() |
#2 |
Member
Lukas
Join Date: Sep 2021
Posts: 36
Rep Power: 4 ![]() |
Good Morning Herbert,
you have two ways to prevent stopping , 1: Increase Residual stop criterium (as example from 1e-4 to 1e-8) 2: remove residual stop criterium completely and monitor the temperature until the profile converges(and stop the simulation manually after 10000 its) |
|
![]() |
![]() |
![]() |
![]() |
#3 |
Senior Member
Gerhard Holzinger
Join Date: Feb 2012
Location: Austria
Posts: 341
Rep Power: 28 ![]() ![]() |
The residuals "only" tell you how good or accurate the flow is being approximated by the simulation. This is independent whether your case has reached equilibrium. In a transient simulation, the flow will always change no matter how small the residuals are.
Thus, if you want to stop the simulation when your case has reached thermal equilibrium, you need to construct some criterion which describes exactly this: is the thermal equilibrium reached? My first try would be to create a codedFunctionObject which describes this, and which writes a file if your condition for thermal equilibrium is met. Then, you can use the stopAtFile functionObject to stop the simulation when it detects the existence of this special file. |
|
![]() |
![]() |
![]() |
![]() |
#4 |
New Member
Herbert Dobernig
Join Date: Apr 2020
Posts: 4
Rep Power: 6 ![]() |
From steady-state solvers like buoyantSimpleFoam I would expect stabilization of results as a mandatory stop criterium.
Thank you for the "work-around" using a codingFunctionObject. Is there good reason, why stabilization of results as stop criterium is not implemented by default in steady-state solvers? |
|
![]() |
![]() |
![]() |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
SU2 Unsteady NACA0012 Tutorial case | Gj419 | SU2 | 1 | November 23, 2021 11:55 |
Convergence order of a simple laplacianFoam case | daven | OpenFOAM Verification & Validation | 0 | September 9, 2021 09:50 |
HronTurekFsi tutorial case using icoFsiElasticNonLinULSolidFoam | MarcelK | OpenFOAM Running, Solving & CFD | 9 | February 25, 2019 02:25 |
Turning the Tutorial propeller case into a flow driven case with 6Dof | efirvida | OpenFOAM | 0 | March 7, 2017 09:45 |
Free surface boudary conditions with SOLA-VOF | Fan | Main CFD Forum | 10 | September 9, 2006 12:24 |