solveFluid.H

Go to the documentation of this file.

if (finalIter)
{
    mesh.data::add("finalIteration", true);
}
 
if (frozenFlow)
{
    #include "EEqns.H"
}
else
{
    fluid.solve();
    fluid.correct();
 
    #include "YEqns.H"
 
    if (faceMomentum)
    {
        #include "pUf/UEqns.H"
        #include "EEqns.H"
        #include "pUf/pEqn.H"
    }
    else
    {
        #include "pU/UEqns.H"
        #include "EEqns.H"
        #include "pU/pEqn.H"
    }
 
    fluid.correctKinematics();
 
    // Update alpha's for new U
    fluid.correctTurbulence();
}
 
if (finalIter)
{
    mesh.data::remove("finalIteration");
}