FreeFOAM The Cross-Platform CFD Toolkit
compressibleContinuityErrs.H
Go to the documentation of this file.
1 {
2  #include <finiteVolume/rhoEqn.H>
3 }
4 {
5  scalar sumLocalContErr =
6  sum
7  (
8  mag(rho.internalField() - (psi*p)().internalField())
9  )/sum(rho.internalField());
10 
11  scalar globalContErr =
12  sum(rho.internalField() - (psi*p)().internalField())
13  /sum(rho.internalField());
14 
16 
17  Info<< "time step continuity errors : sum local = " << sumLocalContErr
18  << ", global = " << globalContErr
19  << ", cumulative = " << cumulativeContErr << endl;
20 }