UEqn.H
Go to the documentation of this file.
1 MRF.correctBoundaryVelocity(U);
2
3 fvVectorMatrix UEqn
4 (
5 fvm::ddt(rho, U) + fvm::div(phi, U)
6 + MRF.DDt(rho, U)
7 + turbulence->divDevRhoReff(U)
8 ==
10 );
11
12 UEqn.relax();
13
14 fvOptions.constrain(UEqn);
15
16 if (pimple.momentumPredictor())
17 {
18 solve
19 (
20 UEqn
21 ==
22 fvc::reconstruct
23 (
24 (
25 - ghf*fvc::snGrad(rho)
26 - fvc::snGrad(p_rgh)
27 )*mesh.magSf()
28 )
29 );
30
31 fvOptions.correct(U);
32 K = 0.5*magSqr(U);
33 }
CGAL::Exact_predicates_exact_constructions_kernel K
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
CEqn solve()