UEqn.H
Go to the documentation of this file.
1  // Solve the Momentum equation
2 
3  MRF.correctBoundaryVelocity(U);
4 
5  UEqn =
6  (
7  fvm::div(phi, U)
8  + MRF.DDt(rho, U)
9  + turb.divDevRhoReff(U)
10  ==
11  fvOptions(rho, U)
12  );
13 
14  UEqn.relax();
15 
16  fvOptions.constrain(UEqn);
17 
18  if (momentumPredictor)
19  {
20  solve
21  (
22  UEqn
23  ==
24  fvc::reconstruct
25  (
26  (
27  - ghf*fvc::snGrad(rho)
28  - fvc::snGrad(p_rgh)
29  )*mesh.magSf()
30  )
31  );
32 
33  fvOptions.correct(U);
34  }
momentumPredictor
const bool momentumPredictor
Definition: readFluidMultiRegionSIMPLEControls.H:6
solve
solve(UEqn==-fvc::grad(p))
Foam::fac::div
tmp< GeometricField< Type, faPatchField, areaMesh > > div(const GeometricField< Type, faePatchField, edgeMesh > &ssf)
Definition: facDiv.C:43
turb
compressible::turbulenceModel & turb
Definition: setRegionFluidFields.H:8
fvOptions
fv::options & fvOptions
Definition: setRegionFluidFields.H:21
ghf
const surfaceScalarField & ghf
Definition: setRegionFluidFields.H:16
mesh
dynamicFvMesh & mesh
Definition: createDynamicFvMesh.H:6
MRF
IOMRFZoneList & MRF
Definition: setRegionFluidFields.H:20
Foam::fvc::reconstruct
tmp< GeometricField< typename outerProduct< vector, Type >::type, fvPatchField, volMesh >> reconstruct(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
Definition: fvcReconstruct.C:49
U
U
Definition: pEqn.H:72
UEqn
fvVectorMatrix & UEqn
Definition: UEqn.H:13
p_rgh
volScalarField & p_rgh
Definition: setRegionFluidFields.H:13
Foam::fvc::snGrad
tmp< GeometricField< Type, fvsPatchField, surfaceMesh > > snGrad(const GeometricField< Type, fvPatchField, volMesh > &vf, const word &name)
Definition: fvcSnGrad.C:40