2 volScalarField rUA("rUA", 1.0/UEqn.A());
3 surfaceScalarField rUAf("(1|A(U))", fvc::interpolate(rUA));
7 surfaceScalarField phiU
9 (fvc::interpolate(U) & mesh.Sf())
10 + fvc::ddtPhiCorr(rUA, U, phi)
13 phi = phiU + rUAf*fvc::interpolate(rhok)*(g & mesh.Sf());
15 for (int nonOrth=0; nonOrth<=nNonOrthCorr; nonOrth++)
19 fvm::laplacian(rUAf, p) == fvc::div(phi)
22 if (corr == nCorr-1 && nonOrth == nNonOrthCorr)
24 pEqn.solve(mesh.solver(p.name() + "Final"));
28 pEqn.solve(mesh.solver(p.name()));
31 if (nonOrth == nNonOrthCorr)
37 U += rUA*fvc::reconstruct((phi - phiU)/rUAf);
38 U.correctBoundaryConditions();
40 #include "continuityErrs.H"