applications/solvers/surfaceTracking/include/scalePhi.H

   1 if(interface.twoFluids() && pNeedRef)
   2 {
   3     phi.boundaryField()[interface.bPatchID()] =
   4         (
   5             U.boundaryField()[interface.bPatchID()]
   6           & mesh.Sf().boundaryField()[interface.bPatchID()]
   7         );
   8
   9     scalarField weights =
  10         mag(phi.boundaryField()[interface.bPatchID()]);
  11
  12     if(mag(gSum(weights)) > VSMALL)
  13     {
  14         weights /= gSum(weights);
  15     }
  16
  17     phi.boundaryField()[interface.bPatchID()] -=
  18         weights*gSum(phi.boundaryField()[interface.bPatchID()]);
  19
  20 //     Info << "Free surface shadow flow rate: "
  21 //         << gSum(phi.boundaryField()[interface.bPatchID()]) << endl;
  22
  23     phi.boundaryField()[interface.bPatchID()] +=
  24         p.boundaryField()[interface.bPatchID()].snGrad()
  25        *mesh.magSf().boundaryField()[interface.bPatchID()]
  26        /AU.boundaryField()[interface.bPatchID()];
  27 }
  28