EEqnRhoSimple.H

Go to the documentation of this file.

{
 
    fvSourceEnergy = fvSource & U;
 
    volScalarField& he = thermo.he();
 
    fvScalarMatrix EEqn(
        fvm::div(phi, he)
        + (he.name() == "e"
               ? fvc::div(phi, volScalarField("Ekp", 0.5 * magSqr(U) + p / rho))
               : fvc::div(phi, volScalarField("K", 0.5 * magSqr(U))))
        - fvm::laplacian(turbulencePtr_->alphaEff(), he)
        - fvSourceEnergy
        - fvOptions(rho, he));
 
    EEqn.relax();
 
    fvOptions.constrain(EEqn);
 
    // get the solver performance info such as initial
    // and final residuals
    SolverPerformance<scalar> solverE = EEqn.solve();
 
    DAUtility::primalResidualControl(solverE, printToScreen_, "he", daGlobalVarPtr_->primalMaxRes);
 
    // bound he
    DAUtility::boundVar(allOptions, he, printToScreen_);
 
    thermo.correct();
}