2312/twoPhaseEulerFoam_2CourantNos_8H_source.html

     #include "CourantNo.H"

 {
     scalar UrCoNum = 0.5*gMax
     (
         fvc::surfaceSum(mag(phi1 - phi2))().primitiveField()/mesh.V().field()
     )*runTime.deltaTValue();

     Info<< "Max Ur Courant Number = " << UrCoNum << endl;

     CoNum = max(CoNum, UrCoNum);
 }
Foam::mag
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)

Foam::max
label max(const labelHashSet &set, label maxValue=labelMin)
Find the max value in labelHashSet, optionally limited by second argument.
Definition: hashSets.C:40

runTime
engineTime & runTime
Definition: createEngineTime.H:13

Foam::endl
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:531

phi1
surfaceScalarField & phi1
Definition: setRegionFluidFields.H:10

mesh
dynamicFvMesh & mesh
Definition: createDynamicFvMesh.H:6

phi2
surfaceScalarField & phi2
Definition: setRegionFluidFields.H:14

CoNum
scalar CoNum
Definition: compressibleCourantNo.H:25

Foam::gMax
Type gMax(const FieldField< Field, Type > &f)
Definition: FieldFieldFunctions.C:573

Foam::Info
messageStream Info
Information stream (stdout output on master, null elsewhere)

Foam::fvc::surfaceSum
tmp< GeometricField< Type, fvPatchField, volMesh > > surfaceSum(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
Definition: fvcSurfaceIntegrate.C:128