UEqn.H
Go to the documentation of this file.
1MRF.correctBoundaryVelocity(U);
2
3fvVectorMatrix UEqn
4(
5 fvm::ddt(alphacRho, U)
6 + MRF.DDt(alphacRho, U)
7 - fvm::Sp(fvc::ddt(rho) + fvc::div(rhoPhi), U)
8 + fvm::div(rhoPhi, U)
9 + turbulence->divDevRhoReff(U)
10 ==
12 + cloudSU
13);
14
15UEqn.relax();
16fvOptions.constrain(UEqn);
17
18volScalarField rAUc(1.0/UEqn.A());
19surfaceScalarField rAUcf(fvc::interpolate(rAUc));
20
21
22surfaceScalarField phicForces
23(
24 (fvc::interpolate(rAUc*cloudVolSUSu) & mesh.Sf())
25);
26
27
28if (pimple.momentumPredictor())
29{
30 solve
31 (
32 UEqn
33 ==
34 fvc::reconstruct
35 (
37 +
38 (
39 fvc::interpolate
40 (
41 mixture.sigmaK()
42 )*fvc::snGrad(alpha1)
43 - ghf*fvc::snGrad(rho)
44 - fvc::snGrad(p_rgh)
45 ) * mesh.magSf()
46 )
47 );
48
49 fvOptions.correct(U);
50}
rhoPhi
Definition: rhoEqn.H:10
volScalarField & p_rgh
fv::options & fvOptions
const surfaceScalarField & ghf
IOMRFZoneList & MRF
pimpleControl & pimple
const volScalarField & alpha1
U
Definition: pEqn.H:72
fvVectorMatrix & UEqn
Definition: UEqn.H:13
dynamicFvMesh & mesh
compressible::turbulenceModel & turbulence
surfaceScalarField phicForces((fvc::interpolate(rAUc *cloudVolSUSu) &mesh.Sf()))
volScalarField rAUc(1.0/UEqn.A())
surfaceScalarField rAUcf(fvc::interpolate(rAUc))
CEqn solve()
Info<< "Creating temperaturePhaseChangeTwoPhaseMixture\n"<< endl;autoPtr< temperaturePhaseChangeTwoPhaseMixture > mixture
Definition: createFields.H:39