Files
openfoam/applications/solvers/multiphase/twoPhaseEulerFoam/CourantNos.H
2012-02-29 16:39:53 +00:00

13 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);
}