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
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}
Y[inertIndex] max(0.0)
rhoPhi
Definition: rhoEqn.H:10
surfaceScalarField & phi
volScalarField & rho2
Basic thermodynamics type based on the use of fitting functions for cp, h, s obtained from the templa...
const volScalarField & T
compressible::turbulenceModel & turbulence
fvScalarMatrix TEqn(fvm::ddt(T)+fvm::div(phi, T) - fvm::laplacian(alphaEff, T)==radiation->ST(rhoCpRef, T)+fvOptions(T))
const surfaceScalarField rhoCpPhi(fvc::interpolate(fluid.Cp()) *rhoPhi)
rhoCp
Definition: TEqn.H:3
kappaEff
Definition: TEqn.H:10
const dimensionedScalar Cp1
Definition: TEqn.H:5
const dimensionedScalar Cp2
Definition: TEqn.H:6
const volScalarField & cp
Definition: TEqn.H:3
tmp< volScalarField > tcp
dimensionedScalar Prt("Prt", dimless, laminarTransport)
Info<< "Creating temperaturePhaseChangeTwoPhaseMixture\n"<< endl;autoPtr< temperaturePhaseChangeTwoPhaseMixture > mixture
Definition: createFields.H:39