Loading...
Searching...
No Matches
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))().primitiveField()/mesh.V().field()
7 )*runTime.deltaTValue();
8
9 Info<< "Max Ur Courant Number = " << UrCoNum << endl;
10
11 CoNum = Foam::max(CoNum, UrCoNum);
12}
surfaceScalarField & phi2
surfaceScalarField & phi1
dynamicFvMesh & mesh
engineTime & runTime
label max(const labelHashSet &set, label maxValue=labelMin)
Find the max value in labelHashSet, optionally limited by second argument.
Definition hashSets.C:40
messageStream Info
Information stream (stdout output on master, null elsewhere).
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition Ostream.H:519
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
Type gMax(const FieldField< Field, Type > &f)
scalar CoNum