EEqn.H
Go to the documentation of this file.
1{
3
4 fvScalarMatrix EEqn
5 (
6 fvm::ddt(rho, he) + mvConvection->fvmDiv(phi, he)
7 + fvc::ddt(rho, K) + fvc::div(phi, K)
8 + (
9 he.name() == "e"
10 ? fvc::div
11 (
12 fvc::absolute(phi/fvc::interpolate(rho), U),
13 p,
14 "div(phiv,p)"
15 )
16 : -dpdt
17 )
18 - fvm::laplacian(turbulence->alphaEff(), he)
19 ==
20 rho*(U&g)
21 + parcels.Sh(he)
22 + surfaceFilm.Sh()
23 + radiation->Sh(thermo, he)
24 + Qdot
25 + fvOptions(rho, he)
26 );
27
28 EEqn.relax();
29
30 fvOptions.constrain(EEqn);
31
32 EEqn.solve();
33
34 fvOptions.correct(he);
35
36 thermo.correct();
37 radiation->correct();
38
39 Info<< "T gas min/max = " << min(T).value() << ", "
40 << max(T).value() << endl;
41}
CGAL::Exact_predicates_exact_constructions_kernel K
Y[inertIndex] max(0.0)
volScalarField & he
Definition: YEEqn.H:52
tmp< fv::convectionScheme< scalar > > mvConvection(fv::convectionScheme< scalar >::New(mesh, fields, phi, mesh.divScheme("div(phi,Yi_h)")))
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
const volScalarField & T
regionModels::surfaceFilmModel & surfaceFilm
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))
volScalarField & dpdt
compressible::turbulenceModel & turbulence
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Definition: volFieldsFwd.H:82
scalar Qdot
Definition: solveChemistry.H:2