1 surfaceScalarField muEff
4 twoPhaseProperties.muf()
5 + fvc::interpolate(rho*turbulence->nut())
12 - fvm::laplacian(muEff, U)
13 //- (fvc::grad(U) & fvc::grad(muf))
14 - fvc::div(muEff*(fvc::interpolate(dev(fvc::grad(U))) & mesh.Sf()))
19 if (momentumPredictor)
21 solve(UEqn == -fvc::grad(p));
24 Info<< "max(U) " << max(mag(U)).value() << endl;