EauEqn.H
Go to the documentation of this file.
1 if (ign.ignited())
2 {
3  volScalarField& heau = thermo.heu();
4 
5  fvScalarMatrix heauEqn
6  (
7  betav*fvm::ddt(rho, heau) + mvConvection->fvmDiv(phi, heau)
8  + (betav*fvc::ddt(rho, K) + fvc::div(phi, K))*rho/thermo.rhou()
9  + (
10  heau.name() == "eau"
11  ? mvConvection->fvcDiv(fvc::absolute(phi, rho, U), p/rho)
12  *rho/thermo.rhou()
13  : -betav*dpdt*rho/thermo.rhou()
14  )
15  - fvm::laplacian(Db, heau)
16 
17  // These terms cannot be used in partially-premixed combustion due to
18  // the resultant inconsistency between ft and heau transport.
19  // A possible solution would be to solve for ftu as well as ft.
20  //- fvm::div(muEff*fvc::grad(b)/(b + 0.001), heau)
21  //+ fvm::Sp(fvc::div(muEff*fvc::grad(b)/(b + 0.001)), heau)
22  ==
23  betav*fvModels.source(rho, heau)
24  );
25 
26  fvConstraints.constrain(heauEqn);
27 
28  heauEqn.solve();
29 
31 }
fvMatrix< scalar > fvScalarMatrix
Definition: fvMatricesFwd.H:42
fluidReactionThermo & thermo
Definition: createFields.H:28
tmp< GeometricField< Type, fvPatchField, volMesh > > div(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
Definition: fvcDiv.C:47
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
tmp< fv::convectionScheme< scalar > > mvConvection(fv::convectionScheme< scalar >::New(mesh, fields, phi, mesh.divScheme("div(phi,ft_b_ha_hau)")))
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Definition: volFieldsFwd.H:57
volScalarField & dpdt
Foam::fvConstraints & fvConstraints
phi
Definition: correctPhi.H:3
bool constrain(fvMatrix< Type > &eqn) const
Apply constraints to an equation.
Foam::fvModels & fvModels
volScalarField Db("Db", turbulence->muEff())
tmp< surfaceScalarField > absolute(const tmp< surfaceScalarField > &tphi, const volVectorField &U)
Return the given relative flux in absolute form.
Definition: fvcMeshPhi.C:188
U
Definition: pEqn.H:72
volScalarField & p