SEqn.H

Go to the documentation of this file.

//Solve the Momentum equation
if (solveScalar)
{
    fvScalarMatrix SEqn
    (
        fvm::ddt(S)
        + fvm::div(phi,S)
        - fvm::laplacian(DT,S)
        - fvm::laplacian(
                turbulence->nut()/Sct,S
                )
    ==  
        fvOptions(S)
    );
 
    SEqn.relax();
    fvOptions.constrain(SEqn);
    solve(SEqn);
    fvOptions.correct(S);
}