8 phi = (fvc::interpolate(U) & mesh.Sf())
9 + fvc::ddtPhiCorr(rUA, U, phi);
13 // Non-orthogonal pressure corrector loop
14 for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)
19 fvm::laplacian(rUA, p) == fvc::div(phi)
22 pEqn.setReference(pRefCell, pRefValue);
28 && nonOrth == nNonOrthCorr
31 pEqn.solve(mesh.solver("pFinal"));
38 if (nonOrth == nNonOrthCorr)
44 #include <finiteVolume/continuityErrs.H>
46 // Explicitly relax pressure for momentum corrector except for last corrector
47 if (oCorr != nOuterCorr-1)
52 U -= rUA*fvc::grad(p);
53 U.correctBoundaryConditions();