CourantNos.H
Go to the documentation of this file.
1 #include "CourantNo.H"
2 
3 {
4  scalar UrCoNum = 0.5*gMax
5  (
6  fvc::surfaceSum(mag(phi1 - phi2))().internalField()/mesh.V().field()
7  )*runTime.deltaTValue();
8 
9  Info<< "Max Ur Courant Number = " << UrCoNum << endl;
10 
11  CoNum = max(CoNum, UrCoNum);
12 }
dimensioned< scalar > mag(const dimensioned< Type > &)
tmp< GeometricField< Type, fvPatchField, volMesh > > surfaceSum(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
messageStream Info
surfaceScalarField & phi2
Definition: createFields.H:24
dynamicFvMesh & mesh
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:251
dimensioned< Type > max(const dimensioned< Type > &, const dimensioned< Type > &)
surfaceScalarField & phi1
Definition: createFields.H:19
scalar CoNum
Type gMax(const FieldField< Field, Type > &f)
conserve internalField()+