applications/solvers/multiphase/twoLiquidMixingFoam/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(muEff*(fvc::interpolate(dev(fvc::grad(U))) & mesh.Sf()))
  15     );
  16
  17     if (oCorr == nOuterCorr-1)
  18     {
  19         if (mesh.solutionDict().relax("UFinal"))
  20         {
  21             UEqn.relax(mesh.solutionDict().relaxationFactor("UFinal"));
  22         }
  23         else
  24         {
  25             UEqn.relax(1);
  26         }
  27     }
  28     else
  29     {
  30         UEqn.relax();
  31     }
  32
  33     if (momentumPredictor)
  34     {
  35         solve
  36         (
  37             UEqn
  38          ==
  39             fvc::reconstruct
  40             (
  41                 (- ghf*fvc::snGrad(rho) - fvc::snGrad(p_rgh))*mesh.magSf()
  42             ),
  43             mesh.solutionDict().solver(oCorr == nOuterCorr-1 ? "UFinal" : "U")
  44         );
  45     }