applications/solvers/multiphase/cavitatingFoam/UEqn.H

   1     surfaceScalarField muEff
   2     (
   3         "muEff",
   4         twoPhaseProperties.muf()
   5       + fvc::interpolate(rho*turbulence->nut())
   6     );
   7
   8     fvVectorMatrix UEqn
   9     (
  10         fvm::ddt(rho, U)
  11       + fvm::div(phi, U)
  12       - fvm::laplacian(muEff, U)
  13     //- (fvc::grad(U) & fvc::grad(muf))
  14       - fvc::div(muEff*(fvc::interpolate(dev(fvc::grad(U))) & mesh.Sf()))
  15     );
  16
  17     UEqn.relax();
  18
  19     if (momentumPredictor)
  20     {
  21         solve(UEqn == -fvc::grad(p));
  22     }