1 tmp<fv::convectionScheme<scalar> > mvConvection
3 fv::convectionScheme<scalar>::New
8 mesh.divScheme("div(phi,Yi_hs)")
13 combustion->correct();
14 dQ = combustion->dQ();
15 label inertIndex = -1;
16 volScalarField Yt(0.0*Y[0]);
20 if (Y[i].name() != inertSpecie)
22 volScalarField& Yi = Y[i];
23 fvScalarMatrix R(combustion->R(Yi));
28 + mvConvection->fvmDiv(phi, Yi)
29 - fvm::laplacian(turbulence->alphaEff(), Yi)
37 YiEqn.solve(mesh.solver("Yi"));
48 Y[inertIndex] = scalar(1) - Yt;
49 Y[inertIndex].max(0.0);
54 + mvConvection->fvmDiv(phi, hs)
55 - fvm::laplacian(turbulence->alphaEff(), hs)
59 + radiation->Shs(thermo)
69 Info<< "min/max(T) = " << min(T).value() << ", " << max(T).value() << endl;