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