OpenFOAM-5.x/applications/solvers/multiphase/twoPhaseEulerFoam/CourantNos.H
2014-12-10 22:40:10 +00:00

12 lines
266 B
C

# include "CourantNo.H"
{
scalar UrCoNum = 0.5*gMax
(
fvc::surfaceSum(mag(phi1 - phi2))().internalField()/mesh.V().field()
)*runTime.deltaTValue();
Info<< "Max Ur Courant Number = " << UrCoNum << endl;
CoNum = max(CoNum, UrCoNum);
}