YEqn.H
Go to the documentation of this file.
1 tmp<fv::convectionScheme<scalar>> mvConvection
2 (
3  fv::convectionScheme<scalar>::New
4  (
5  mesh,
6  fields,
7  phi,
8  mesh.schemes().div("div(phi,Yi_h)")
9  )
10 );
11 
12 reaction->correct();
13 
14 forAll(Y, i)
15 {
16  if (composition.solve(i))
17  {
18  volScalarField& Yi = Y[i];
19 
20  fvScalarMatrix YiEqn
21  (
22  fvm::ddt(rho, Yi)
23  + mvConvection->fvmDiv(phi, Yi)
24  + thermophysicalTransport->divj(Yi)
25  ==
26  reaction->R(Yi)
27  + fvModels.source(rho, Yi)
28  );
29 
30  YiEqn.relax();
31 
32  fvConstraints.constrain(YiEqn);
33 
34  YiEqn.solve("Yi");
35 
36  fvConstraints.constrain(Yi);
37  }
38 }
39 
40 composition.normalise();
Foam::fvScalarMatrix
fvMatrix< scalar > fvScalarMatrix
Definition: fvMatricesFwd.H:42
Foam::compressible::New
autoPtr< CompressibleMomentumTransportModel > New(const volScalarField &rho, const volVectorField &U, const surfaceScalarField &phi, const viscosity &viscosity)
Definition: compressibleMomentumTransportModelTemplates.C:34
composition
basicSpecieMixture & composition
Definition: createFieldRefs.H:3
reaction
combustionModel & reaction
Definition: setRegionFluidFields.H:3
mesh
fvMesh & mesh
Definition: createMeshesPostProcess.H:9
Foam::fvc::ddt
tmp< GeometricField< Type, fvPatchField, volMesh > > ddt(const dimensioned< Type > dt, const fvMesh &mesh)
Definition: fvcDdt.C:45
Foam::volScalarField
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Definition: volFieldsFwd.H:58
fvConstraints
Foam::fvConstraints & fvConstraints
Definition: setRegionFluidFields.H:37
phi
phi
Definition: correctPhi.H:3
thermophysicalTransport
fluidReactionThermophysicalTransportModel & thermophysicalTransport
Definition: setRegionFluidFields.H:16
Foam::fvConstraints::constrain
bool constrain(fvMatrix< Type > &eqn) const
Apply constraints to an equation.
Definition: fvConstraintsTemplates.C:29
fvModels
Foam::fvModels & fvModels
Definition: setRegionFluidFields.H:36
fields
Info<< "Reading field p_rgh\"<< endl;volScalarField p_rgh(IOobject("p_rgh", runTime.timeName(), mesh, IOobject::MUST_READ, IOobject::AUTO_WRITE), mesh);pressureReference pressureReference(p, p_rgh, pimple.dict(), thermo.incompressible());mesh.schemes().setFluxRequired(p_rgh.name());hydrostaticInitialisation(p_rgh, p, rho, U, gh, ghf, pRef, thermo, pimple.dict());Info<< "Creating field dpdt\"<< endl;volScalarField dpdt(IOobject("dpdt", runTime.timeName(), mesh), mesh, dimensionedScalar(p.dimensions()/dimTime, 0));Info<< "Creating field kinetic energy K\"<< endl;volScalarField K("K", 0.5 *magSqr(U));dimensionedScalar initialMass=fvc::domainIntegrate(rho);multivariateSurfaceInterpolationScheme< scalar >::fieldTable fields
Definition: createFields.H:131
mvConvection
tmp< fv::convectionScheme< scalar > > mvConvection(fv::convectionScheme< scalar >::New(mesh, fields, phi, mesh.schemes().div("div(phi,Yi_h)")))
Y
PtrList< volScalarField > & Y
Definition: createFieldRefs.H:4
forAll
forAll(Y, i)
Definition: YEqn.H:14