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

Setting up the PIMPLE algorithm correctly

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

Like Tree3Likes
  • 2 Post By clapointe
  • 1 Post By clapointe

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   August 28, 2018, 09:32
Default Setting up the PIMPLE algorithm correctly
  #1
Member
 
benoit favier
Join Date: Jun 2017
Posts: 64
Rep Power: 9
Benben is on a distinguished road
Hello,

I am trying to setup a case of "reactingFoam" on openfoam 6.
I read this tutorial on the pimple algorythm. If i understand correctly, the solver is supposed to stop doing "Outercorrector" loops once the residual reaches a certain point.

However, in my case, it always does the maximum number of "Outercorrector" loops regardless of my residual tolerance.

My "fvSolution" file is the following.

Code:
solvers
{
    "rho.*"
    {
        solver          diagonal;
    }
    p
    {
        solver          PCG;
        preconditioner  DIC;
        tolerance       1e-6;
        relTol          0.;
    }
    
    pFinal
    {
        $p;
        tolerance        1e-6;
        relTol           0.;
    }
    "(U|h|k|omega)"
    {
        solver          PBiCGStab;
        preconditioner  DILU;
        tolerance       1e-6;
        relTol          0.;
    }
    
    "(U|h|k|epsilon)Final"
    {
        $U;
        relTol          0;
    }
    Yi
    {
        solver          PBiCGStab;
        preconditioner  DILU;
        tolerance       1e-6;
        relTol          0.;
    }
}
PIMPLE
{
    nNonOrthogonalCorrectors 0;
    nCorrectors          2;
    nOuterCorrectors    50;
 
    OuterCorrectorResidualControl
    {
        "(U|p)" 
         {
              tolerance       1.E-4; 
              relTol          0.; 
         } 
    }
}

relaxationFactors
{
    fields
    {
        p               0.4;
    }
    equations
    {
        U               0.4;
        k               0.9;
        e               0.9;
        ".*"            0.9;
    }
}
What did'nt I understand here ?
Benben is offline   Reply With Quote

Old   August 28, 2018, 13:03
Default
  #2
Senior Member
 
Join Date: Aug 2015
Posts: 494
Rep Power: 15
clapointe is on a distinguished road
It should be outerCorrectorResidualControl -- with the lowercase "o". See this tutorial for reference : https://github.com/OpenFOAM/OpenFOAM...tem/fvSolution.

Caelan
mechy and Benben like this.
clapointe is offline   Reply With Quote

Old   August 29, 2018, 02:30
Default
  #3
Member
 
benoit favier
Join Date: Jun 2017
Posts: 64
Rep Power: 9
Benben is on a distinguished road
Thank you.
OpenFoam should really warn when there are useless elements in a dictionnary.

But then it tells me some strange things at the begining of a simulation :
Code:
PIMPLE: No convergence criteria found

PIMPLE: Corrector convergence criteria found
        "(U|p|h)": tolerance 0.001, relTol 0
        Calclations will do 5 corrections if the convergence criteria are not met
Benben is offline   Reply With Quote

Old   August 29, 2018, 10:25
Default
  #4
Senior Member
 
Join Date: Aug 2015
Posts: 494
Rep Power: 15
clapointe is on a distinguished road
For what it's worth I see the same double message at the top of my log file. The convergence monitoring works just fine though.

Caelan
Benben likes this.
clapointe is offline   Reply With Quote

Old   August 29, 2018, 14:58
Default
  #5
Member
 
benoit favier
Join Date: Jun 2017
Posts: 64
Rep Power: 9
Benben is on a distinguished road
Ok, thank you for your help
Benben 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
pimple algorithm details for transient, incompressible flow cfdopenfoam Main CFD Forum 0 November 15, 2015 23:31
Reducing div(U) in PIMPLE algorithm hrushi.397 OpenFOAM Running, Solving & CFD 0 April 19, 2015 05:12
A question on the PIMPLE algorithm GerhardHolzinger OpenFOAM Running, Solving & CFD 4 February 13, 2015 06:49
Explanation of PIMPLE algorithm Tobi OpenFOAM 7 November 8, 2014 07:01
why PIMPLE doesn't converge at each time step but no diverging? immortality OpenFOAM Running, Solving & CFD 7 May 19, 2013 14:16


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