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(mag(rho - rho0 - psi*(p - p0)))/sum(rho)).value();
7 
8  scalar globalContErr = (sum(rho - rho0 - psi*(p - p0))/sum(rho)).value();
9 
11 
12  Info<< "time step continuity errors : sum local = " << sumLocalContErr
13  << ", global = " << globalContErr
14  << ", cumulative = " << cumulativeContErr << endl;
15 }
16 
17 // ************************ vim: set sw=4 sts=4 et: ************************ //