37 if (
mesh.nInternalFaces())
45 CoNum = 0.5*
gMax(sumPhi/
mesh.V().field())*runTime.deltaTValue();
48 0.5*(
gSum(sumPhi)/
gSum(
mesh.V().field()))*runTime.deltaTValue();
56 )().primitiveField()/
mesh.V().field()
57 )*
sqrt(magg).
value()*runTime.deltaTValue();
60 Info<<
"Courant number mean: " << meanCoNum
61 <<
" max: " << CoNum <<
endl;
63 Info<<
"Gravity wave Courant number max: " << waveCoNum
dimensionedScalar sqrt(const dimensionedScalar &ds)
Ostream & endl(Ostream &os)
Add newline and flush stream.
Type gSum(const FieldField< Field, Type > &f)
const Type & value() const
Return const reference to value.
volScalarField scalarField(fieldObject, mesh)
Type gMax(const FieldField< Field, Type > &f)
static tmp< GeometricField< Type, fvsPatchField, surfaceMesh > > interpolate(const GeometricField< Type, fvPatchField, volMesh > &tvf, const surfaceScalarField &faceFlux, Istream &schemeData)
Interpolate field onto faces using scheme given by Istream.
volScalarField & h
Planck constant.
dimensioned< scalar > mag(const dimensioned< Type > &)
tmp< GeometricField< Type, fvPatchField, volMesh > > surfaceSum(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)