13 fvc::interpolate(
rho)*fvc::flux(
HbyA)
21 while (
simple.correctNonOrthogonal())
34 if (
simple.finalNonOrthogonalIter())
45 #include "compressibleContinuityErrs.H"
48 U.correctBoundaryConditions();
56 Info<<
"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 dimensionedScalar rhoMin
const dimensionedScalar rhoMax
tmp< volScalarField > rAU
const dictionary & simple