applications/solvers/multiphase/compressibleInterDyMFoam/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(rhoPhi, U)
  12       - fvm::laplacian(muEff, U)
  13       - (fvc::grad(U) & fvc::grad(muEff))
  14     //- fvc::div(muf*(mesh.Sf() & fvc::interpolate(fvc::grad(U)().T())))
  15     );
  16
  17     UEqn.relax();
  18
  19     if (momentumPredictor)
  20     {
  21         solve
  22         (
  23             UEqn
  24          ==
  25             fvc::reconstruct
  26             (
  27                 fvc::interpolate(rho)*(g & mesh.Sf())
  28               + (
  29                     fvc::interpolate(interface.sigmaK())*fvc::snGrad(alpha1)
  30                   - fvc::snGrad(p)
  31                 ) * mesh.magSf()
  32             )
  33         );
  34     }