TEqn.H
Go to the documentation of this file.
1 {
2  tmp<volScalarField> tcp(thermo->Cp());
3  const volScalarField& cp = tcp();
4 
5  const dimensionedScalar Cp1 = thermo->Cp1();
6  const dimensionedScalar Cp2 = thermo->Cp2();
7 
8  rhoCp = rho*cp;
9 
10  kappaEff = thermo->kappa() + rho*cp*turbulence->nut()/Prt;
11 
12  const surfaceScalarField rhoCpPhi
13  (
14  "rhoCpPhi",
15  rhoPhi*(Cp1 - Cp2) + phi*rho2*Cp2
16  );
17 
18  fvScalarMatrix TEqn
19  (
20  fvm::ddt(rhoCp, T)
21  + fvm::div(rhoCpPhi, T)
22  - fvm::Sp(fvc::ddt(rhoCp) + fvc::div(rhoCpPhi), T)
23  - fvm::laplacian(kappaEff, T)
24  + mixture->TSource()
25  );
26 
27 
28  TEqn.relax();
29  TEqn.solve();
30 
31  Info<< "min/max(T) = " << min(T).value() << ", "
32  << max(T).value() <<endl;
33 }
Foam::fvScalarMatrix
fvMatrix< scalar > fvScalarMatrix
Definition: fvMatricesFwd.H:37
cp
const volScalarField & cp
Definition: TEqn.H:3
Foam::max
label max(const labelHashSet &set, label maxValue=labelMin)
Find the max value in labelHashSet, optionally limited by second argument.
Definition: hashSets.C:40
Foam::fac::div
tmp< GeometricField< Type, faPatchField, areaMesh > > div(const GeometricField< Type, faePatchField, edgeMesh > &ssf)
Definition: facDiv.C:43
mixture
Info<< "Creating temperaturePhaseChangeTwoPhaseMixture\"<< endl;autoPtr< temperaturePhaseChangeTwoPhaseMixture > mixture
Definition: createFields.H:39
Foam::endl
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:531
Foam::fac::ddt
tmp< GeometricField< Type, faPatchField, areaMesh > > ddt(const dimensioned< Type > dt, const faMesh &mesh)
Definition: facDdt.C:40
kappaEff
kappaEff
Definition: TEqn.H:10
rhoCpPhi
const surfaceScalarField rhoCpPhi(fvc::interpolate(fluid.Cp()) *rhoPhi)
Foam::volScalarField
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Definition: volFieldsFwd.H:81
thermo
psiReactionThermo & thermo
Definition: createFields.H:28
Foam::fac::laplacian
tmp< GeometricField< Type, faPatchField, areaMesh > > laplacian(const GeometricField< Type, faPatchField, areaMesh > &vf, const word &name)
Definition: facLaplacian.C:40
Foam::min
label min(const labelHashSet &set, label minValue=labelMax)
Find the min value in labelHashSet, optionally limited by second argument.
Definition: hashSets.C:26
rhoCp
rhoCp
Definition: TEqn.H:3
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
rho2
volScalarField & rho2
Definition: setRegionFluidFields.H:28
T
const volScalarField & T
Definition: createFieldRefs.H:2
Prt
dimensionedScalar Prt("Prt", dimless, laminarTransport)
Cp1
const dimensionedScalar Cp1
Definition: TEqn.H:5
Foam::dimensionedScalar
dimensioned< scalar > dimensionedScalar
Dimensioned scalar obtained from generic dimensioned type.
Definition: dimensionedScalarFwd.H:33
tcp
tmp< volScalarField > tcp
Definition: setRegionSolidFields.H:6
rhoPhi
rhoPhi
Definition: rhoEqn.H:10
Foam::Info
messageStream Info
Information stream (stdout output on master, null elsewhere)
TEqn
fvScalarMatrix TEqn(fvm::ddt(T)+fvm::div(phi, T) - fvm::laplacian(alphaEff, T)==radiation->ST(rhoCpRef, T)+fvOptions(T))
Cp2
const dimensionedScalar Cp2
Definition: TEqn.H:6
Foam::surfaceScalarField
GeometricField< scalar, fvsPatchField, surfaceMesh > surfaceScalarField
Definition: surfaceFieldsFwd.H:64
Sp
zeroField Sp
Definition: alphaSuSp.H:2