12 fvc::interpolate(
rho)*fvc::flux(
HbyA)
33 if (
simple.finalNonOrthogonalIter())
44#include "compressibleContinuityErrs.H"
47U.correctBoundaryConditions();
55Info<<
"p min/max = " << min(
p).value() <<
", " <<
max(
p).value() << endl;
Basic thermodynamics type based on the use of fitting functions for cp, h, s obtained from the templa...
const volScalarField & psi
surfaceScalarField rhorAUf("rhorAUf", fvc::interpolate(rho *rAU))
tmp< fvVectorMatrix > tUEqn(fvm::ddt(rho, U)+fvm::div(phi, U)+MRF.DDt(rho, U)+turbulence->divDevRhoReff(U)==fvOptions(rho, U))
const volScalarField psip0(psi *p)
const dimensionedScalar rhoMin
const dimensionedScalar rhoMax
tmp< volScalarField > rAU
const dictionary & simple