|
[Sponsors] |
Continuity. Stop critertion of PISO and SIMPLE in OpenFOAM solvers. |
![]() |
|
LinkBack | Thread Tools | Search this Thread | Display Modes |
![]() |
![]() |
#1 |
Senior Member
Claus Meister
Join Date: Aug 2009
Location: Wiesbaden, Germany
Posts: 241
Rep Power: 18 ![]() |
Hallo Foamers!
I have a question related to OpenFOAM solvers with SIMPLE or PISO loop. Why was as stop criterion the number of iteration used or implemented,respectively, and not a number of accuracy, e.g. 1e-10 (sum,global,cumulative)? Has this a specific reason? Cheers! |
|
![]() |
![]() |
![]() |
![]() |
#2 |
Super Moderator
Niklas Nordin
Join Date: Mar 2009
Location: Stockholm, Sweden
Posts: 693
Rep Power: 30 ![]() ![]() |
||
![]() |
![]() |
![]() |
![]() |
#3 |
Senior Member
Claus Meister
Join Date: Aug 2009
Location: Wiesbaden, Germany
Posts: 241
Rep Power: 18 ![]() |
I did not know that! Honestly, I am a bit surprised. How can I set up such criterion? In guides and documents, I have nothing read about this, anywhere! Where I can find it?
|
|
![]() |
![]() |
![]() |
![]() |
#4 |
Super Moderator
Niklas Nordin
Join Date: Mar 2009
Location: Stockholm, Sweden
Posts: 693
Rep Power: 30 ![]() ![]() |
if you are using simpleFoam it is already there, so you just have to add it to the SIMPLE dictionary.
Code:
SIMPLE { nNonOrthogonalCorrectors 1; convergence 1.0e-8; } where you will find initConvergenceCheck.H Code:
// initialize values for convergence checks scalar eqnResidual = 1, maxResidual = 0; scalar convergenceCriterion = 0; simple.readIfPresent("convergence", convergenceCriterion); and then the convergenceCheck.H which is the last include of the iteration-loop Code:
// check convergence if (maxResidual < convergenceCriterion) { Info<< "reached convergence criterion: " << convergenceCriterion << endl; runTime.writeAndEnd(); Info<< "latestTime = " << runTime.timeName() << endl; } |
|
![]() |
![]() |
![]() |
![]() |
#5 |
Senior Member
Claus Meister
Join Date: Aug 2009
Location: Wiesbaden, Germany
Posts: 241
Rep Power: 18 ![]() |
Thank you for you posted. I will looking for it!
|
|
![]() |
![]() |
![]() |
Thread Tools | Search this Thread |
Display Modes | |
|
|
![]() |
||||
Thread | Thread Starter | Forum | Replies | Last Post |
continuityErrs.h, coninuity equation, automatic stop are reaching continuity | idrama | OpenFOAM | 0 | September 20, 2010 15:09 |
SIMPLE, SIMPLER or PISO? | Gem Leku | Main CFD Forum | 7 | October 24, 2007 11:20 |
SIMPLE and PISO | Aditya | Main CFD Forum | 5 | April 1, 2006 18:52 |
diffrence beetwen SIMPLE and PISO | grzes | FLUENT | 0 | January 29, 2006 05:57 |
PISO vs. SIMPLE | benedikt flurl | Main CFD Forum | 2 | April 14, 2005 06:54 |