3 alphat.correctBoundaryConditions();
autoPtr< compressible::turbulenceModel > turbulence
dimensionedScalar Pr("Pr", dimless, laminarTransport)
fvMatrix< scalar > fvScalarMatrix
fvScalarMatrix TEqn(fvm::ddt(T)+fvm::div(phi, T)-fvm::laplacian(alphaEff, T)==radiation->ST(rhoCpRef, T)+fvOptions(T))
tmp< GeometricField< Type, fvPatchField, volMesh > > div(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
tmp< GeometricField< Type, fvPatchField, volMesh > > laplacian(const GeometricField< Type, fvPatchField, volMesh > &vf, const word &name)
tmp< GeometricField< Type, fvPatchField, volMesh > > ddt(const dimensioned< Type > dt, const fvMesh &mesh)
GeometricField< scalar, fvPatchField, volMesh > volScalarField
virtual void correct()
Main update/correction routine.
dimensionedScalar Prt("Prt", dimless, laminarTransport)
virtual tmp< fvScalarMatrix > ST(const dimensionedScalar &rhoCp, volScalarField &T) const
Temperature source term.
dimensionedScalar rhoCpRef("rhoCpRef", dimDensity *dimEnergy/dimMass/dimTemperature, 1.0)
dimensionedScalar beta("beta", dimless/dimTemperature, laminarTransport)
dimensionedScalar TRef("TRef", dimTemperature, laminarTransport)
volScalarField alphaEff("alphaEff", turbulence->nu()/Pr+alphat)
autoPtr< radiation::radiationModel > radiation(radiation::radiationModel::New(T))