36 if (
mesh.nInternalFaces())
43 CoNum = 0.5*
gMax(sumPhi/
mesh.V().field())*runTime.deltaTValue();
46 0.5*(
gSum(sumPhi)/
gSum(
mesh.V().field()))*runTime.deltaTValue();
48 acousticCoNum = 0.5*
gMax 53 )().primitiveField()/
mesh.V().field()
54 )*runTime.deltaTValue();
57 Info<<
"phi Courant Number mean: " << meanCoNum
59 <<
" acoustic max: " << acousticCoNum
dimensionedScalar sqrt(const dimensionedScalar &ds)
Ostream & endl(Ostream &os)
Add newline and flush stream.
tmp< surfaceScalarField > interpolate(const RhoType &rho)
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)
const volScalarField & psi