EEqn.H
Go to the documentation of this file.
1{
3
4 fvScalarMatrix EEqn
5 (
6 fvm::div(phi, he)
7 + (
8 he.name() == "e"
9 ? fvc::div(phi, volScalarField("Ekp", 0.5*magSqr(U) + p/rho))
10 : fvc::div(phi, volScalarField("K", 0.5*magSqr(U)))
11 )
12 - fvm::laplacian(turbulence->alphaEff(), he)
13 ==
14 rho*(U&g)
15 + radiation->Sh(thermo, he)
16 + fvOptions(rho, he)
17 );
18
19 EEqn.relax();
20
21 fvOptions.constrain(EEqn);
22
23 EEqn.solve();
24
25 fvOptions.correct(he);
26
27 thermo.correct();
28 radiation->correct();
29}
volScalarField & he
Definition: YEEqn.H:52
const uniformDimensionedVectorField & g
fv::options & fvOptions
surfaceScalarField & phi
Basic thermodynamics type based on the use of fitting functions for cp, h, s obtained from the templa...
U
Definition: pEqn.H:72
volScalarField & p
fvScalarMatrix EEqn(fvm::ddt(rho, he)+mvConvection->fvmDiv(phi, he)+fvc::ddt(rho, K)+fvc::div(phi, K)+(he.name()=="e" ? fvc::div(fvc::absolute(phi/fvc::interpolate(rho), U), p, "div(phiv,p)") :-dpdt) - fvm::laplacian(turbulence->alphaEff(), he)==Qdot+fvOptions(rho, he))
autoPtr< radiation::radiationModel > radiation(radiation::radiationModel::New(T))
compressible::turbulenceModel & turbulence
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Definition: volFieldsFwd.H:82