applications/solvers/surfaceTracking/bubbleInterTrackFoam/scaleSpacePhi.H

   1 if(runTime.timeIndex() > 1)
   2 {
   3     fvsPatchScalarField& spacePhi = phi.boundaryField()[spacePatchID];
   4
   5     scalar inletFlux = gSum(neg(spacePhi)*spacePhi);
   6
   7     scalar outletFlux = gSum(pos(spacePhi)*spacePhi);
   8
   9     if(outletFlux < VSMALL)
  10     {
  11         outletFlux = VSMALL;
  12     }
  13
  14     scalar outflowScaling = -inletFlux/outletFlux;
  15
  16     spacePhi += pos(spacePhi)*spacePhi*(outflowScaling - 1.0);
  17
  18 //     Info<< "Space patch continuity: "
  19 //         << gSum(phi.boundaryField()[spacePatchID]) << endl;
  20
  21 //     U.boundaryField()[spacePatchID] ==
  22 //         U.boundaryField()[spacePatchID]
  23 //       + pos(spacePhi)*U.boundaryField()[spacePatchID]*(outflowScaling - 1.0);
  24 }