applications/solvers/multiphase/twoPhaseEulerFoam/CourantNos.H

   1 #   include "CourantNo.H"
   2
   3 {
   4     scalar UrCoNum = 0.5*gMax
   5     (
   6         fvc::surfaceSum(mag(phia - phib))().internalField()/mesh.V().field()
   7     )*runTime.deltaTValue();
   8
   9     Info<< "Max Ur Courant Number = " << UrCoNum << endl;
  10
  11     CoNum = max(CoNum, UrCoNum);
  12 }