TEqn.H

Go to the documentation of this file.

{
    alphat = turbulence->nut()/Prt;
    alphat.correctBoundaryConditions();
 
    volScalarField alphaEff("alphaEff", turbulence->nu()/Pr + alphat);
 
    fvScalarMatrix TEqn
    (
        fvm::ddt(T)
      + fvm::div(phi, T)
      - fvm::laplacian(alphaEff, T)
     ==
        radiation->ST(rhoCpRef, T)
      + fvOptions(T)
    );
 
    TEqn.relax();
 
    fvOptions.constrain(TEqn);
 
    TEqn.solve();
 
    radiation->correct();
 
    fvOptions.correct(T);
 
    rhok = 1.0 - beta*(T - TRef);
}