13 if (
p.needReference())
24 while (
pimple.correctNonOrthogonal())
28 fvm::laplacian(alphacf*
rAUcf,
p)
30 fvc::ddt(alphac) + fvc::div(alphacf*
phiHbyA)
35 pEqn.solve(
mesh.solver(
p.select(
pimple.finalInnerIter())));
37 if (
pimple.finalNonOrthogonalIter())
45 Uc.correctBoundaryConditions();
surfaceScalarField phic(mixture.cAlpha() *mag(alphaPhic/mesh.magSf()))
surfaceScalarField rAUcf("Dp", fvc::interpolate(rAUc))
volScalarField rAUc(1.0/UcEqn.A())
surfaceScalarField phicForces(fvc::flux(rAUc *cloudVolSUSu/rhoc)+rAUcf *(g &mesh.Sf()))
fvVectorMatrix UcEqn(fvm::ddt(alphac, Uc)+fvm::div(alphaPhic, Uc) - fvm::Sp(fvc::ddt(alphac)+fvc::div(alphaPhic), Uc)+continuousPhaseTurbulence->divDevRhoReff(Uc)==(1.0/rhoc) *cloudSU)
adjustPhi(phiHbyA, U, p_rgh)
GeometricField< vector, fvPatchField, volMesh > volVectorField