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(rhoPhi, U)
6 - fvm::Sp(fvc::ddt(rho) + fvc::div(rhoPhi), U)
7 + MRF.DDt(rho, U)
8 + turbulence->divDevRhoReff(rho, U)
9 ==
11 );
12
13 UEqn.relax();
14
15 fvOptions.constrain(UEqn);
16
17 if (pimple.momentumPredictor())
18 {
19 solve
20 (
21 UEqn
22 ==
23 cellMask*fvc::reconstruct
24 (
25 (
26 mixture.surfaceTensionForce()
27 - ghf*fvc::snGrad(rho)
28 - fvc::snGrad(p_rgh)
29 ) * mesh.magSf()
30 )
31 );
32
33 fvOptions.correct(U);
34 }
rhoPhi
Definition: rhoEqn.H:10
volScalarField & p_rgh
fv::options & fvOptions
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()
Info<< "Creating temperaturePhaseChangeTwoPhaseMixture\n"<< endl;autoPtr< temperaturePhaseChangeTwoPhaseMixture > mixture
Definition: createFields.H:39