TEqn.H

Go to the documentation of this file.

//Solve the Momentum equation
//Info << "sovleScalarT = " << solveScalarT << endl;
if (solveScalarT)
{
 //   Info << " go into TEqn" << endl;
//dimensionedScalar Pr("Pr", dimless, laminarTransport);
//
//dimensionedScalar Prt("Prt", dimless, laminarTransport);
 
        // Thermal expansion coefficient [1/K]
    dimensionedScalar Pr(laminarTransport.lookup("Pr"));
 
    // Turbulent Prandtl number
    dimensionedScalar Prt(laminarTransport.lookup("Prt"));
//Info <<"pt and ptr is :" << Pr <<"  " << Prt << endl;
 
    alphat = turbulence->nut()/Prt;
 
    alphat.correctBoundaryConditions();
 
    volScalarField alphaEff("alphaEff", turbulence->nu()/Pr + alphat); 
 
    fvScalarMatrix TEqn
    (
//      fvm::ddt(T)
         fvm::div(phi,T)
        - fvm::laplacian(DT,T)
    
    ==  
        fvOptions(T)
    );
 
    TEqn.relax();
    fvOptions.constrain(TEqn);
    solve(TEqn);
    fvOptions.correct(T);
}