TEqn.H
Go to the documentation of this file.
1 {
2  radiation->correct();
3  rhoCp = rho*fluid.Cp();
4 
5  const surfaceScalarField rhoCpPhi(fvc::interpolate(fluid.Cp())*rhoPhi);
6 
7  const volScalarField kappaEff
8  (
9  "kappaEff",
10  fluid.kappa() + fluid.Cp()*turbulence->mut()/fluid.Prt()
11  );
12 
13  fvScalarMatrix TEqn
14  (
15  fvm::ddt(rhoCp, T)
16  + fvm::div(rhoCpPhi, T, "div(phi,T)")
17  - fvm::Sp(fvc::ddt(rhoCp) + fvc::div(rhoCpPhi), T)
18  - fvm::laplacian(kappaEff, T, "laplacian(kappa,T)")
19  ==
20  fluid.heatTransfer(T)
21  + radiation->ST(T)
22  + fvOptions(rhoCp, T)
23  );
24 
25  TEqn.relax();
26 
27  fvOptions.constrain(TEqn);
28 
29  TEqn.solve();
30 
31  fvOptions.correct(T);
32 
33  fluid.correct();
34 
35  Info<< "min/max(T) = "
36  << min(T).value() << ", " << max(T).value() << endl;
37 }
fluid
twoPhaseSystem & fluid
Definition: setRegionFluidFields.H:1
Foam::fvScalarMatrix
fvMatrix< scalar > fvScalarMatrix
Definition: fvMatricesFwd.H:37
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::interpolate
static tmp< GeometricField< Type, faePatchField, edgeMesh > > interpolate(const GeometricField< Type, faPatchField, areaMesh > &tvf, const edgeScalarField &faceFlux, Istream &schemeData)
Interpolate field onto faces using scheme given by Istream.
Foam::fac::div
tmp< GeometricField< Type, faPatchField, areaMesh > > div(const GeometricField< Type, faePatchField, edgeMesh > &ssf)
Definition: facDiv.C:43
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)
fvOptions
fv::options & fvOptions
Definition: setRegionFluidFields.H:21
Foam::volScalarField
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Definition: volFieldsFwd.H:81
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
T
const volScalarField & T
Definition: createFieldRefs.H:2
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))
Foam::surfaceScalarField
GeometricField< scalar, fvsPatchField, surfaceMesh > surfaceScalarField
Definition: surfaceFieldsFwd.H:64
radiation
autoPtr< radiation::radiationModel > radiation(radiation::radiationModel::New(T))
Sp
zeroField Sp
Definition: alphaSuSp.H:2