v3/a03035_source.html

 {
     volScalarField& hea = thermo.he();

     solve
     (
         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)
     );

     thermo.correct();
 }
dpdt
volScalarField & dpdt
Definition: setRegionFluidFields.H:12

K
CGAL::Exact_predicates_exact_constructions_kernel K
Definition: CGALTriangulation3DKernel.H:56

betav
const volScalarField & betav
Definition: setRegionSolidFields.H:30

Foam::fvc::laplacian
tmp< GeometricField< Type, fvPatchField, volMesh > > laplacian(const GeometricField< Type, fvPatchField, volMesh > &vf, const word &name)
Definition: fvcLaplacian.C:45

Foam::MULES::interpolate
tmp< surfaceScalarField > interpolate(const RhoType &rho)
Definition: IMULESTemplates.C:40

Foam::fvc::ddt
tmp< GeometricField< Type, fvPatchField, volMesh > > ddt(const dimensioned< Type > dt, const fvMesh &mesh)
Definition: fvcDdt.C:45

rho
rho
Definition: readInitialConditions.H:96

solve
solve(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))

Foam::volScalarField
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Definition: volFieldsFwd.H:52

p
volScalarField & p
Definition: createFields.H:51

Db
volScalarField Db("Db", turbulence->muEff())

phi
surfaceScalarField & phi
Definition: setRegionFluidFields.H:8

mvConvection
tmp< fv::convectionScheme< scalar > > mvConvection(fv::convectionScheme< scalar >::New(                       mesh,                       fields,                       phi,                       mesh.divScheme("div(phi,Yi_h)")       ))

Foam::fvc::div
tmp< GeometricField< Type, fvPatchField, volMesh > > div(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
Definition: fvcDiv.C:47

thermo
psiReactionThermo & thermo
Definition: createFields.H:32