applications/solvers/multiphase/interPhaseChangeFoam/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::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()))
  16     );
  17
  18     UEqn.relax();
  19
  20     if (momentumPredictor)
  21     {
  22         solve
  23         (
  24             UEqn
  25          ==
  26             fvc::reconstruct
  27             (
  28                 fvc::interpolate(rho)*(g & mesh.Sf())
  29               + (
  30                     fvc::interpolate(interface.sigmaK())*fvc::snGrad(alpha1)
  31                   - fvc::snGrad(p)
  32                 ) * mesh.magSf()
  33             )
  34         );
  35     }