35 if (
mesh.nInternalFaces())
51 CoNum = 0.5*
gMax(sumPhi/
mesh.V().field())*runTime.deltaTValue();
54 0.5*(
gSum(sumPhi)/
gSum(
mesh.V().field()))*runTime.deltaTValue();
57 Info<<
"Courant Number mean: " << meanCoNum
58 <<
" max: " << CoNum <<
endl;
#define forAll(list, i)
Loop across all elements in list.
dimensioned< Type > max(const dimensioned< Type > &, const dimensioned< Type > &)
multiphaseSystem::phaseModelList & phases
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)