FreeFOAM The Cross-Platform CFD Toolkit
continuityErrs.H
Go to the documentation of this file.
1 {
2  volScalarField thermoRho = psi*p + (1.0 - gamma)*rhol0;
3 
5 
6  scalar sumLocalContErr =
7  (
8  fvc::domainIntegrate(mag(rho - thermoRho))/totalMass
9  ).value();
10 
11  scalar globalContErr =
12  (
13  fvc::domainIntegrate(rho - thermoRho)/totalMass
14  ).value();
15 
17 
18  Info<< "time step continuity errors : sum local = " << sumLocalContErr
19  << ", global = " << globalContErr
20  << ", cumulative = " << cumulativeContErr
21  << endl;
22 }