8 + alphacf*fvc::interpolate(
rho*
rAUc)*fvc::ddtCorr(
U,
phi)
28 while (
pimple.correctNonOrthogonal())
30 surfaceScalarField Dp(
"Dp", alphacf*
rAUcf);
32 fvScalarMatrix p_rghEqn
34 fvm::laplacian(Dp,
p_rgh)
36 fvc::ddt(alphac) + fvc::div(alphacf*
phiHbyA)
43 if (
pimple.finalNonOrthogonalIter())
51 +
rAUc*fvc::reconstruct((
phig - p_rghEqn.flux()/alphacf)/
rAUcf);
53 U.correctBoundaryConditions();
65 p += dimensionedScalar
surfaceScalarField rAUcf("Dp", fvc::interpolate(rAUc))
volScalarField rAUc(1.0/UcEqn.A())
surfaceScalarField phicForces(fvc::flux(rAUc *cloudVolSUSu/rhoc)+rAUcf *(g &mesh.Sf()))
const surfaceScalarField & ghf
const volScalarField & gh
const volScalarField & alpha1
surfaceScalarField phig("phig", -rhorAUf *ghf *fvc::snGrad(rho) *mesh.magSf())
adjustPhi(phiHbyA, U, p_rgh)
GeometricField< vector, fvPatchField, volMesh > volVectorField
Info<< "Creating temperaturePhaseChangeTwoPhaseMixture\n"<< endl;autoPtr< temperaturePhaseChangeTwoPhaseMixture > mixture