v4/lagrangian_2reactingParcelFoam_2simpleReactingParcelFoam_2UEqn_8H_source.html

     MRF.correctBoundaryVelocity(U);

     tmp<fvVectorMatrix> tUEqn
     (
         fvm::div(phi, U)
       + MRF.DDt(rho, U)
       + turbulence->divDevRhoReff(U)
      ==
         rho()*g
       + parcels.SU(U)
       + fvOptions(rho, U)
     );
     fvVectorMatrix& UEqn = tUEqn.ref();

     UEqn.relax();

     fvOptions.constrain(UEqn);

     solve(UEqn == -fvc::grad(p));

     fvOptions.correct(U);
turbulence
autoPtr< compressible::turbulenceModel > turbulence
Definition: createFields.H:23

phi
surfaceScalarField & phi
Definition: setRegionFluidFields.H:8

Foam::fvc::grad
tmp< GeometricField< typename outerProduct< vector, Type >::type, fvPatchField, volMesh >> grad(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
Definition: fvcGrad.C:52

U
U
Definition: pEqn.H:83

solve
solve(UEqn==-fvc::grad(p))

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

tUEqn
tmp< fvVectorMatrix > tUEqn(fvm::ddt(rho, U)+fvm::div(phi, U)+MRF.DDt(rho, U)+turbulence->divDevRhoReff(U)==fvOptions(rho, U))

fvOptions
fv::options & fvOptions
Definition: setRegionFluidFields.H:16

MRF
IOMRFZoneList & MRF
Definition: setRegionFluidFields.H:15

rho
rho
Definition: readInitialConditions.H:96

g
const dimensionedVector & g
Definition: setRegionFluidFields.H:33

Foam::fvVectorMatrix
fvMatrix< vector > fvVectorMatrix
Definition: fvMatricesFwd.H:45

p
volScalarField & p
Definition: createFieldRefs.H:4