1 surfaceScalarField muEff
4 twoPhaseProperties->muf()
5 + fvc::interpolate(rho*turbulence->nut())
12 - fvm::Sp(fvc::ddt(rho) + fvc::div(rhoPhi), U)
13 - fvm::laplacian(muEff, U)
14 - (fvc::grad(U) & fvc::grad(muEff))
15 //- fvc::div(muEff*(fvc::interpolate(dev2(fvc::grad(U))) & mesh.Sf()))
20 if (momentumPredictor)
29 fvc::interpolate(interface.sigmaK())*fvc::snGrad(alpha1)
30 - ghf*fvc::snGrad(rho)