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

Residual Control and InitialResidual: How does it exactly work?

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

Like Tree11Likes
  • 2 Post By hrushi.397
  • 7 Post By alexeym
  • 2 Post By chriss85

Reply
 
LinkBack Thread Tools Search this Thread Display Modes
Old   March 17, 2015, 12:15
Default Residual Control and InitialResidual: How does it exactly work?
  #1
Member
 
Hrushi
Join Date: Jan 2013
Posts: 58
Rep Power: 13
hrushi.397 is on a distinguished road
Dear all,

Could someone please clarify to me how residualControl in fvSolution works? As far as I know, residualControl helps you control the flow of the simulation based on the InitialResidual. If the initialResidual during successive corrector steps falls below the level as specified in the residualControl, then the corrector should stop, irrespective of number of corrector steps specified in nCorrectors specified. Is my interpretation correct?

Regards,

Hrushi
kcavatar and JD_PM like this.
hrushi.397 is offline   Reply With Quote

Old   March 18, 2015, 04:09
Default
  #2
Senior Member
 
Agustín Villa
Join Date: Apr 2013
Location: Alcorcón
Posts: 313
Rep Power: 14
agustinvo is on a distinguished road
Hi,

in each variable

Code:
    U
    {
        solver          GAMG;
        ...
    
       tolerance       1e-12;
       relTol          0.;
    }
tolerance will stop you iterations for that equation in function of the Final Residual.
Then the residualControl is working with the InitialResidual, as you said.

I think this is what you were asking.
agustinvo is offline   Reply With Quote

Old   March 18, 2015, 06:12
Default
  #3
Member
 
Hrushi
Join Date: Jan 2013
Posts: 58
Rep Power: 13
hrushi.397 is on a distinguished road
Hi Augstin,

Thank you for the reply. I understand the tolerance and how it works. And if I understand what you are saying, then the residualControl in fvSolution specifies the control on the initialResidual.

In my solver, I tried to set residualControl for pressure at 0.1. But I could see that even after such a relaxed residualControl, my simulations ran for 15 correctors that I had specified earlier. I do not understand why this is happening. Any suggestions?

Regards,

Hrushikesh
hrushi.397 is offline   Reply With Quote

Old   March 18, 2015, 06:39
Default
  #4
Senior Member
 
Alexey Matveichev
Join Date: Aug 2011
Location: Nancy, France
Posts: 1,930
Rep Power: 38
alexeym has a spectacular aura aboutalexeym has a spectacular aura about
Send a message via Skype™ to alexeym
Hi,

resildualControl limits number of outer correctors (nOuterCorrectors). Just take a look at pimpleFoam source:

Code:
        // --- Pressure-velocity PIMPLE corrector loop
        while (pimple.loop())
        {
            ...
            // --- Pressure corrector loop
            while (pimple.correct())
            {
                #include "pEqn.H"
            }
            ...
        }
You have got two while loops, inner stops when pimple.correct() returns false, outer stops when pimple.loop() returns false. So take a look at pimpleControl sources:

pimpleControl.C
Code:
void Foam::pimpleControl::read()
{
    ...
    nCorrPIMPLE_ = pimpleDict.lookupOrDefault<label>("nOuterCorrectors", 1);
    nCorrPISO_ = pimpleDict.lookupOrDefault<label>("nCorrectors", 1);
    ...
}
pimpleControlI.H
Code:
inline bool Foam::pimpleControl::correct()
{
    corrPISO_++;

    if (debug)
    {
        Info<< algorithmName_ << " correct: corrPISO = " << corrPISO_ << endl;
    }

    if (corrPISO_ <= nCorrPISO_)
    {
        return true;
    }
    else
    {
        corrPISO_ = 0;
        return false;
    }
}
So pimple.correct will return false only when number of corrector iterations exceed nCorrectors. While pimple.loop

Code:
bool Foam::pimpleControl::loop()
{
    ...

    if (corr_ == nCorrPIMPLE_ + 1)
    {
        if ((!residualControl_.empty()) && (nCorrPIMPLE_ != 1))
        {
            Info<< algorithmName_ << ": not converged within "
                << nCorrPIMPLE_ << " iterations" << endl;
        }

        corr_ = 0;
        mesh_.data::remove("finalIteration");
        return false;
    }

    bool completed = false;
    if (converged_ || criteriaSatisfied())
    {
        if (converged_)
        {
            Info<< algorithmName_ << ": converged in " << corr_ - 1
                << " iterations" << endl;

            mesh_.data::remove("finalIteration");
            corr_ = 0;
            converged_ = false;

            completed = true;
        }
        else
        {
            Info<< algorithmName_ << ": iteration " << corr_ << endl;
            storePrevIterFields();

            mesh_.data::add("finalIteration", true);
            converged_ = true;
        }
    }
    ...
    return !completed;
}
will continue till number of iterations exceed nOutetCorrectors or this condition becomes true

Code:
converged_ || criteriaSatisfied()
and it is criteriaSatisfied() that checks initial residuals:

Code:
bool Foam::pimpleControl::criteriaSatisfied()
{
    ...
            const bool absCheck = residual < residualControl_[fieldI].absTol;
            bool relCheck = false;

            scalar relative = 0.0;
            if (!storeIni)
            {
                const scalar iniRes =
                    residualControl_[fieldI].initialResidual
                  + ROOTVSMALL;

                relative = residual/iniRes;
                relCheck = relative < residualControl_[fieldI].relTol;
            }

            achieved = achieved && (absCheck || relCheck);
    ...
    return checked && achieved;
}
Though IIRC it was already explained several times on the forum.
ilhado, granzer, cryabroad and 4 others like this.

Last edited by alexeym; March 18, 2015 at 10:03.
alexeym is offline   Reply With Quote

Old   March 18, 2015, 08:53
Default
  #5
Member
 
Hrushi
Join Date: Jan 2013
Posts: 58
Rep Power: 13
hrushi.397 is on a distinguished road
Hi Alexey,

Thank you very much for your response. It is pretty clear to me now.

Regards,

Hrushi
hrushi.397 is offline   Reply With Quote

Old   March 18, 2015, 13:10
Default
  #6
Senior Member
 
Join Date: Oct 2013
Posts: 397
Rep Power: 18
chriss85 will become famous soon enough
Also see http://www.cfd-online.com/Forums/blo...hm-part-i.html and http://www.cfd-online.com/Forums/blo...m-part-ii.html for a good explanation.
Tobi and dalschouten like this.
chriss85 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



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