v6/heatTransfer_2buoyantBoussinesqSimpleFoam_2TEqn_8H_source.html

 {
     alphat = turbulence->nut()/Prt;
     alphat.correctBoundaryConditions();

     volScalarField alphaEff("alphaEff", turbulence->nu()/Pr + alphat);

     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();

     fvOptions.correct(T);

     rhok = 1.0 - beta*(T - TRef);
 }
Pr
dimensionedScalar Pr("Pr", dimless, laminarTransport)

Foam::fvScalarMatrix
fvMatrix< scalar > fvScalarMatrix
Definition: fvMatricesFwd.H:42

fvOptions
fv::options & fvOptions
Definition: setRegionFluidFields.H:51

phi
surfaceScalarField & phi
Definition: setRegionFluidFields.H:28

Foam::fvc::div
tmp< GeometricField< Type, fvPatchField, volMesh > > div(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
Definition: fvcDiv.C:47

Foam::fvc::laplacian
tmp< GeometricField< Type, fvPatchField, volMesh > > laplacian(const GeometricField< Type, fvPatchField, volMesh > &vf, const word &name)
Definition: fvcLaplacian.C:45

Foam::volScalarField
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Definition: volFieldsFwd.H:52

turbulence
Info<< "Reading field U\"<< endl;volVectorField U(IOobject("U", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);volScalarField rho(IOobject("rho", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::AUTO_WRITE), thermo.rho());volVectorField rhoU(IOobject("rhoU", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), rho *U);volScalarField rhoE(IOobject("rhoE", runTime.timeName(), mesh, IOobject::NO_READ, IOobject::NO_WRITE), rho *(e+0.5 *magSqr(U)));surfaceScalarField pos(IOobject("pos", runTime.timeName(), mesh), mesh, dimensionedScalar("pos", dimless, 1.0));surfaceScalarField neg(IOobject("neg", runTime.timeName(), mesh), mesh, dimensionedScalar("neg", dimless, -1.0));surfaceScalarField phi("phi", fvc::flux(rhoU));Info<< "Creating turbulence model\"<< endl;autoPtr< compressible::turbulenceModel > turbulence(compressible::turbulenceModel::New(rho, U, phi, thermo))
Definition: createFields.H:94

T
const volScalarField & T
Definition: createFieldRefs.H:2

Prt
dimensionedScalar Prt("Prt", dimless, laminarTransport)

rhoCpRef
dimensionedScalar rhoCpRef("rhoCpRef", dimDensity *dimEnergy/dimMass/dimTemperature, 1.0)

beta
dimensionedScalar beta("beta", dimless/dimTemperature, laminarTransport)

TRef
dimensionedScalar TRef("TRef", dimTemperature, laminarTransport)

TEqn
fvScalarMatrix TEqn(fvm::ddt(T)+fvm::div(phi, T) - fvm::laplacian(alphaEff, T)==radiation->ST(rhoCpRef, T)+fvOptions(T))

alphaEff
volScalarField alphaEff("alphaEff", turbulence->nu()/Pr+alphat)

rhok
rhok
Definition: TEqn.H:27

radiation
autoPtr< radiation::radiationModel > radiation(radiation::radiationModel::New(T))