35 if (
mesh.nInternalFaces())
40 /
rho.primitiveField()
43 CoNum = 0.5*
gMax(sumPhi/
mesh.V().field())*runTime.deltaTValue();
46 0.5*(
gSum(sumPhi)/
gSum(
mesh.V().field()))*runTime.deltaTValue();
49 Info<<
"Courant Number mean: " << meanCoNum
50 <<
" max: " << CoNum <<
endl;
Ostream & endl(Ostream &os)
Add newline and flush stream.
Type gSum(const FieldField< Field, Type > &f)
volScalarField scalarField(fieldObject, mesh)
Type gMax(const FieldField< Field, Type > &f)
dimensioned< scalar > mag(const dimensioned< Type > &)
tmp< GeometricField< Type, fvPatchField, volMesh > > surfaceSum(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)