22 fvc::interpolate(
rho)*fvc::flux(
HbyA)
34 surfaceScalarField
phid
44 fvc::ddt(
rho) +
psi*correction(fvm::ddt(
p))
50 while (
pimple.correctNonOrthogonal())
57 pEqn.solve(
mesh.solver(
p.select(
pimple.finalInnerIter())));
59 if (
pimple.finalNonOrthogonalIter())
69 fvc::ddt(
rho) +
psi*correction(fvm::ddt(
p))
75 while (
pimple.correctNonOrthogonal())
79 pEqn.solve(
mesh.solver(
p.select(
pimple.finalInnerIter())));
81 if (
pimple.finalNonOrthogonalIter())
93U.correctBoundaryConditions();
99 p.correctBoundaryConditions();
105#include "compressibleContinuityErrs.H"
118 dpdt -= fvc::div(fvc::meshPhi(
rho,
U),
p);
CGAL::Exact_predicates_exact_constructions_kernel K
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)
fvScalarMatrix pDDtEqn(fvc::ddt(rho)+psi *correction(fvm::ddt(p))+fvc::div(phiHbyA)==fvOptions(psi, p, rho.name()))
surfaceScalarField phid("phid", fvc::interpolate(psi) *(fvc::flux(HbyA)+MRF.zeroFilter(rhorAUf *fvc::ddtCorr(rho, U, phi)/fvc::interpolate(rho))))
const dimensionedScalar rhoMin
const dimensionedScalar rhoMax
const pressureControl & pressureControl
tmp< volScalarField > rAU