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}
scalar sumLocalContErr
scalar globalContErr
volScalarField & p
const volScalarField & psi
const volScalarField & p0
Definition: EEqn.H:36
dimensioned< Type > sum(const DimensionedField< Type, GeoMesh > &df)
dimensioned< typename typeOfMag< Type >::type > mag(const dimensioned< Type > &dt)
scalar rho0