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 
8  (
9  fvm::ddt(T)
10  + fvm::div(phi, T)
12  ==
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 }
dimensionedScalar Pr("Pr", dimless, laminarTransport)
fvScalarMatrix TEqn(fvm::ddt(T)+fvm::div(phi, T)-fvm::laplacian(alphaEff, T)==radiation->ST(rhoCpRef, T)+fvOptions(T))
volScalarField alphaEff("alphaEff", turbulence->nu()/Pr+alphat)
dimensionedScalar beta("beta", dimless/dimTemperature, laminarTransport)
tmp< GeometricField< Type, fvPatchField, volMesh > > laplacian(const GeometricField< Type, fvPatchField, volMesh > &vf, const word &name)
Definition: fvcLaplacian.C:45
autoPtr< compressible::turbulenceModel > turbulence
Definition: createFields.H:23
tmp< GeometricField< Type, fvPatchField, volMesh > > ddt(const dimensioned< Type > dt, const fvMesh &mesh)
Definition: fvcDdt.C:45
dimensionedScalar Prt("Prt", dimless, laminarTransport)
virtual void correct()
Main update/correction routine.
const volScalarField & T
Definition: createFields.H:25
rhok
Definition: TEqn.H:27
fv::IOoptionList & fvOptions
fvMatrix< scalar > fvScalarMatrix
Definition: fvMatricesFwd.H:42
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Definition: volFieldsFwd.H:52
dimensionedScalar rhoCpRef("rhoCpRef", dimDensity *dimEnergy/dimMass/dimTemperature, 1.0)
surfaceScalarField & phi
autoPtr< radiation::radiationModel > radiation(radiation::radiationModel::New(T))
dimensionedScalar TRef("TRef", dimTemperature, laminarTransport)
tmp< GeometricField< Type, fvPatchField, volMesh > > div(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
Definition: fvcDiv.C:47
virtual tmp< fvScalarMatrix > ST(const dimensionedScalar &rhoCp, volScalarField &T) const
Temperature source term.