TEqn.H
Go to the documentation of this file.
1{
2 alphat = turbulence->nut()/Prt;
3 alphat.correctBoundaryConditions();
4
5 volScalarField alphaEff("alphaEff", turbulence->nu()/Pr + alphat);
6
7 fvScalarMatrix TEqn
8 (
9 fvm::ddt(T)
10 + fvm::div(phi, T)
11 - fvm::laplacian(alphaEff, T)
12 ==
13 radiation->ST(rhoCpRef, T)
14 + fvOptions(T)
15 );
16
17 TEqn.relax();
18
19 fvOptions.constrain(TEqn);
20
21 TEqn.solve();
22
23 radiation->correct();
24
25 fvOptions.correct(T);
26
27 rhok = 1.0 - beta*(T - TRef);
28}
fv::options & fvOptions
surfaceScalarField & phi
const volScalarField & T
dimensionedScalar rhoCpRef("rhoCpRef", dimDensity *dimEnergy/dimMass/dimTemperature, 1.0)
autoPtr< radiation::radiationModel > radiation(radiation::radiationModel::New(T))
compressible::turbulenceModel & turbulence
rhok
Definition: TEqn.H:27
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))
dimensionedScalar TRef("TRef", dimTemperature, laminarTransport)
dimensionedScalar beta("beta", dimless/dimTemperature, laminarTransport)
dimensionedScalar Prt("Prt", dimless, laminarTransport)
dimensionedScalar Pr("Pr", dimless, laminarTransport)