OpenFOAM-2.4.x/applications/solvers/multiphase/twoPhaseEulerFoam/CourantNos.H
2014-02-17 16:00:08 +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);
}