35 scalar sumLocalContErr =
runTime.deltaTValue()*
36 mag(contErr)().weightedAverage(
mesh.V()).value();
38 scalar globalContErr =
runTime.deltaTValue()*
39 contErr.weightedAverage(
mesh.V()).value();
42 Info<<
"time step continuity errors : sum local = " << sumLocalContErr
43 <<
", global = " << globalContErr
volScalarField::Internal contErr((fvc::ddt(rho)+fvc::div(rhoPhi) -(fvModels.source(alpha1, mixture.thermo1().rho())&rho1) -(fvModels.source(alpha2, mixture.thermo2().rho())&rho2))())
tmp< GeometricField< Type, fvPatchField, volMesh > > div(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
Ostream & endl(Ostream &os)
Add newline and flush stream.
tmp< GeometricField< Type, fvPatchField, volMesh > > ddt(const dimensioned< Type > dt, const fvMesh &mesh)
GeometricField< scalar, fvPatchField, volMesh > volScalarField
dimensioned< scalar > mag(const dimensioned< Type > &)