applications/solvers/multiphase/twoPhaseEulerFoam/CourantNos.H

   1 #   include "CourantNo.H"
   2
   3 {
   4     scalar UrCoNum = max
   5     (
   6         mesh.surfaceInterpolation::deltaCoeffs()*mag(phia - phib)
   7        /mesh.magSf()
   8     ).value()*runTime.deltaT().value();
   9
  10     Info<< "Max Ur Courant Number = " << UrCoNum << endl;
  11
  12     CoNum = max(CoNum, UrCoNum);
  13 }