This repository has been archived on 2023-11-20. You can view files and clone it, but cannot push or open issues or pull requests.
foam-extend4.1-coherent-io/applications/solvers/multiphase/settlingFoam/compressibleContinuityErrs.H

22 lines
599 B
C
Raw Normal View History

scalar sumLocalContErr =
runTime.deltaT().value()*
mag
(
fvc::ddt(rho)
+ fvc::div(phi)
)().weightedAverage(rho*mesh.V()).value();
scalar globalContErr =
runTime.deltaT().value()*
(
fvc::ddt(rho)
+ fvc::div(phi)
)().weightedAverage(rho*mesh.V()).value();
cumulativeContErr += globalContErr;
Info<< "time step continuity errors : sum local = " << sumLocalContErr
<< ", global = " << globalContErr
<< ", cumulative = " << cumulativeContErr
<< endl;