3 volScalarField rUA = 1.0/UEqn.A();
8 surfaceScalarField phid
13 (fvc::interpolate(U) & mesh.Sf())
14 + fvc::ddtPhiCorr(rUA, rho, U, phi)
18 for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)
22 betav*fvm::ddt(psi, p)
24 - fvm::laplacian(rho*invA, p)
29 if (nonOrth == nNonOrthCorr)
38 fvc::interpolate(rho)*
40 (fvc::interpolate(U) & mesh.Sf())
41 + fvc::ddtPhiCorr(rUA, rho, U, phi)
44 for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)
48 betav*fvm::ddt(psi, p)
50 - fvm::laplacian(rho*invA, p)
55 if (nonOrth == nNonOrthCorr)
63 #include "continuityErrs.H"
65 U -= invA & (betav*fvc::grad(p));
66 U.correctBoundaryConditions();
68 DpDt = fvc::DDt(surfaceScalarField("phiU", phi/fvc::interpolate(rho)), p);