solveFluid.H

Go to the documentation of this file.

if (finalIter)
{
    mesh.data::add("finalIteration", true);
}
 
if (frozenFlow)
{
    #include "EEqn.H"
}
else
{
    if (oCorr == 0)
    {
        #include "rhoEqn.H"
    }
 
    #include "UEqn.H"
    #include "YEqn.H"
    #include "EEqn.H"
 
    if (!coupled)
    {
        Info<< "\nSolving for fluid region " << fluidRegions[i].name() << endl;
 
        // --- PISO loop
        for (int corr=0; corr<nCorr; corr++)
        {
            #include "pEqn.H"
        }
 
        turbulence.correct();
 
        rho = thermo.rho();
    }
}
 
if (finalIter)
{
    mesh.data::remove("finalIteration");
}