applications/solvers/multiphase/lesInterFoam/UEqn.H

   1     surfaceScalarField muEff
   2     (
   3         "muEff",
   4         twoPhaseProperties.muf()
   5       + fvc::interpolate(rho*turbulence->nuSgs())
   6     );
   7
   8     fvVectorMatrix UEqn
   9     (
  10         fvm::ddt(rho, U)
  11       + fvm::div(rhoPhi, U)
  12       - fvm::laplacian(muEff, U)
  13       - (fvc::grad(U) & fvc::grad(muEff))
  14     //- fvc::div(muEff*(fvc::interpolate(dev(fvc::grad(U))) & mesh.Sf()))
  15     );
  16
  17     if (momentumPredictor)
  18     {
  19         solve
  20         (
  21             UEqn
  22          ==
  23             fvc::reconstruct
  24             (
  25                 (
  26                     fvc::interpolate(interface.sigmaK())*fvc::snGrad(gamma)
  27                   - ghf*fvc::snGrad(rho)
  28                   - fvc::snGrad(pd)
  29                 ) * mesh.magSf()
  30             )
  31         );
  32     }