compressibleContinuityErrs.H
Go to the documentation of this file.
1 {
2  scalar sumLocalContErr =
3  (sum(mag(rho - rho0 - psi*(p - p0)))/sum(rho)).value();
4 
5  scalar globalContErr = (sum(rho - rho0 - psi*(p - p0))/sum(rho)).value();
6 
8 
9  Info<< "time step continuity errors : sum local = " << sumLocalContErr
10  << ", global = " << globalContErr
11  << ", cumulative = " << cumulativeContErr << endl;
12 }
dimensioned< scalar > mag(const dimensioned< Type > &)
scalar rho0
messageStream Info
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:251
volScalarField & p
Definition: createFields.H:51
const volScalarField & psi
Definition: createFields.H:24
dimensioned< Type > sum(const DimensionedField< Type, GeoMesh > &df)
scalar globalContErr