YEqn.H
Go to the documentation of this file.
1tmp<fv::convectionScheme<scalar>> mvConvection
2(
3 fv::convectionScheme<scalar>::New
4 (
5 mesh,
6 fields,
7 phi,
8 mesh.divScheme("div(phi,Yi_h)")
9 )
10);
11
12{
13 combustion->correct();
14 Qdot = combustion->Qdot();
15 volScalarField Yt(0.0*Y[0]);
16
18 {
19 if (i != inertIndex && composition.active(i))
20 {
21 volScalarField& Yi = Y[i];
22
23 fvScalarMatrix YEqn
24 (
25 mvConvection->fvmDiv(phi, Yi)
26 - fvm::laplacian(turbulence->muEff(), Yi)
27 ==
28 parcels.SYi(i, Yi)
29 + combustion->R(Yi)
30 + fvOptions(rho, Yi)
31 );
32
33 YEqn.relax();
34
35 fvOptions.constrain(YEqn);
36
37 YEqn.solve(mesh.solver("Yi"));
38
39 fvOptions.correct(Yi);
40
41 Yi.max(0.0);
42 Yt += Yi;
43 }
44 }
45
46 Y[inertIndex] = scalar(1) - Yt;
47 Y[inertIndex].max(0.0);
48}
fv::options & fvOptions
surfaceScalarField & phi
basicSpecieMixture & composition
PtrList< volScalarField > & Y
Qdot
Definition: YEqn.H:14
volScalarField Yt(0.0 *Y[0])
tmp< fv::convectionScheme< scalar > > mvConvection(fv::convectionScheme< scalar >::New(mesh, fields, phi, mesh.divScheme("div(phi,Yi_h)")))
dynamicFvMesh & mesh
label inertIndex
compressible::turbulenceModel & turbulence
Info<< "Creating combustion model\n"<< endl;autoPtr< CombustionModel< psiReactionThermo > > combustion(CombustionModel< psiReactionThermo >::New(thermo, turbulence()))
multivariateSurfaceInterpolationScheme< scalar >::fieldTable fields
Definition: createFields.H:97
#define forAll(list, i)
Loop across all elements in list.
Definition: stdFoam.H:333