EaEqn.H
Go to the documentation of this file.
1 {
2  volScalarField& hea = thermo.he();
3 
5  (
6  betav*fvm::ddt(rho, hea) + mvConvection->fvmDiv(phi, hea)
8  + (
9  hea.name() == "ea"
10  ? fvc::div
11  (
13  p,
14  "div(phiv,p)"
15  )
16  : -betav*dpdt
17  )
18  - fvm::laplacian(Db, hea)
19  + betav*fvOptions(rho, hea)
20  );
21 
22  EaEqn.relax();
23 
24  fvOptions.constrain(EaEqn);
25 
26  EaEqn.solve();
27 
28  fvOptions.correct(hea);
29 
30  thermo.correct();
31 }
fvMatrix< scalar > fvScalarMatrix
Definition: fvMatricesFwd.H:42
fvScalarMatrix EaEqn(betav *fvm::ddt(rho, hea)+mvConvection->fvmDiv(phi, hea)+betav *fvc::ddt(rho, K)+fvc::div(phi, K)+(hea.name()=="ea" ? fvc::div(phi/fvc::interpolate(rho), p, "div(phiv,p)") :-betav *dpdt) - fvm::laplacian(Db, hea)+betav *fvOptions(rho, hea))
fv::options & fvOptions
surfaceScalarField & phi
tmp< GeometricField< Type, fvPatchField, volMesh > > div(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
Definition: fvcDiv.C:47
rhoReactionThermo & thermo
Definition: createFields.H:28
tmp< GeometricField< Type, fvPatchField, volMesh > > laplacian(const GeometricField< Type, fvPatchField, volMesh > &vf, const word &name)
Definition: fvcLaplacian.C:45
CGAL::Exact_predicates_exact_constructions_kernel K
tmp< GeometricField< Type, fvPatchField, volMesh > > ddt(const dimensioned< Type > dt, const fvMesh &mesh)
Definition: fvcDdt.C:45
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Definition: volFieldsFwd.H:52
volScalarField & dpdt
const volScalarField & betav
static tmp< GeometricField< Type, fvsPatchField, surfaceMesh > > interpolate(const GeometricField< Type, fvPatchField, volMesh > &tvf, const surfaceScalarField &faceFlux, Istream &schemeData)
Interpolate field onto faces using scheme given by Istream.
volScalarField Db("Db", turbulence->muEff())
tmp< fv::convectionScheme< scalar > > mvConvection(fv::convectionScheme< scalar >::New(mesh, fields, phi, mesh.divScheme("div(phi,Yi_h)")))
volScalarField & p