UEqn.H
Go to the documentation of this file.
1 // Solve the momentum equation
2
3 MRF.correctBoundaryVelocity(U);
4
5 fvVectorMatrix UEqn
6 (
7 fvm::ddt(U) + fvm::div(phi, U)
8 + MRF.DDt(U)
9 + turbulence->divDevReff(U)
10 ==
12 );
13
14 UEqn.relax();
15
16 fvOptions.constrain(UEqn);
17
18 if (pimple.momentumPredictor())
19 {
20 solve
21 (
22 UEqn
23 ==
24 fvc::reconstruct
25 (
26 (
27 - ghf*fvc::snGrad(rhok)
28 - fvc::snGrad(p_rgh)
29 )*mesh.magSf()
30 )
31 );
32
33 fvOptions.correct(U);
34 }
volScalarField & p_rgh
fv::options & fvOptions
surfaceScalarField & phi
const surfaceScalarField & ghf
IOMRFZoneList & MRF
pimpleControl & pimple
U
Definition: pEqn.H:72
fvVectorMatrix & UEqn
Definition: UEqn.H:13
dynamicFvMesh & mesh
compressible::turbulenceModel & turbulence
rhok
Definition: TEqn.H:27
CEqn solve()