v3/a02994_source.html

     // Solve the momentum equation

     MRF.correctBoundaryVelocity(U);

     tmp<fvVectorMatrix> UEqn
     (
         fvm::div(phi, U)
       + MRF.DDt(U)
       + turbulence->divDevReff(U)
      ==
         fvOptions(U)
     );

     UEqn().relax();

     fvOptions.constrain(UEqn());

     if (simple.momentumPredictor())
     {
         solve
         (
             UEqn()
           ==
             fvc::reconstruct
             (
                 (
                   - ghf*fvc::snGrad(rhok)
                   - fvc::snGrad(p_rgh)
                 )*mesh.magSf()
             )
         );

         fvOptions.correct(U);
     }
Foam::fvc::reconstruct
tmp< GeometricField< typename outerProduct< vector, Type >::type, fvPatchField, volMesh >> reconstruct(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
Definition: fvcReconstruct.C:54

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

p_rgh
volScalarField & p_rgh
Definition: setRegionFluidFields.H:31

turbulence
autoPtr< compressible::turbulenceModel > turbulence
Definition: createFields.H:23

mesh
dynamicFvMesh & mesh
Definition: createDynamicFvMesh.H:18

rhok
rhok
Definition: TEqn.H:27

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

MRF
IOMRFZoneList & MRF
Definition: setRegionFluidFields.H:15

UEqn
fvVectorMatrix UEqn(fvm::ddt(rho, U)+fvm::div(phi, U)+turbulence->divDevRhoReff(U)==fvOptions(rho, U))

Foam::fvc::snGrad
tmp< GeometricField< Type, fvsPatchField, surfaceMesh > > snGrad(const GeometricField< Type, fvPatchField, volMesh > &vf, const word &name)
Definition: fvcSnGrad.C:45

phi
surfaceScalarField & phi
Definition: setRegionFluidFields.H:8

simple
const dictionary & simple
Definition: readFluidMultiRegionSIMPLEControls.H:1

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

U
U
Definition: pEqn.H:82

ghf
const surfaceScalarField & ghf
Definition: setRegionFluidFields.H:35