TEqn.H

Go to the documentation of this file.

 
    alphat = turbulence->nut()/Prt;
    alphat.correctBoundaryConditions();
    //Info << "input alphat is : " << alphat << "m2/s"<<endl;
 
    volScalarField alphaTimeAverage("alpha", turbulence->nu()/Pr);
    //Info << "input alpha is : " << alphaTimeAverage << "m2/s"<<endl;
 
    volScalarField alphaEff("alphaEff", turbulence->nu()/Pr + alphat);
    //Info << "input alphaEff is : " << alphaEff << "m2/s"<<endl;
 
//  Info << "turbulence->nu() is " << turbulence->nu() << endl;
//  Info << "prant number Pr is " << Pr << endl;
//  Info << "turbulence->nut() is " << turbulence->nut() << endl;
//  Info << "prant number Prt is " << Prt << endl;
 
    fvScalarMatrix TEqn
    (
        fvm::div(phi, T)
      - fvm::laplacian(alphaEff, T)
     ==
        radiation->ST(rhoCpRef, T)
      + fvOptions(T)
    );
 
    TEqn.relax();
 
    fvOptions.constrain(TEqn);
 
    TEqn.solve();
 
    radiation->correct();
 
    //Info << "radiation model is " <<endl;
 
    //Info << radiation << endl;
 
    fvOptions.correct(T);
 
    rhok = 1.0 - beta*(T - TRef);
 
    Tcoeff =T - TRef ;