1 surfaceScalarField muEff
4 twoPhaseProperties.muf()
5 + fvc::interpolate(rho*turbulence->nut())
12 - fvm::laplacian(muEff, U)
13 - (fvc::grad(U) & fvc::grad(muEff))
14 //- fvc::div(muEff*(fvc::interpolate(dev(fvc::grad(U))) & mesh.Sf()))
17 if (oCorr == nOuterCorr-1)
26 if (momentumPredictor)
34 fvc::interpolate(rho)*(g & mesh.Sf())
35 - mesh.magSf()*fvc::snGrad(p)
37 mesh.solver(oCorr == nOuterCorr-1 ? "UFinal" : "U")