23 fvc::interpolate(
rho)*fvc::flux(
HbyA)
41 surfaceScalarField
phid
55 fvc::ddt(
rho) +
psi*correction(fvm::ddt(
p))
61 while (
pimple.correctNonOrthogonal())
68 pEqn.solve(
mesh.solver(
p.select(
pimple.finalInnerIter())));
70 if (
pimple.finalNonOrthogonalIter())
83 fvc::ddt(
rho) +
psi*correction(fvm::ddt(
p))
89 while (
pimple.correctNonOrthogonal())
93 pEqn.solve(
mesh.solver(
p.select(
pimple.finalInnerIter())));
95 if (
pimple.finalNonOrthogonalIter())
103#include "compressibleContinuityErrs.H"
109U.correctBoundaryConditions();
115 p.correctBoundaryConditions();
130 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
tmp< fvVectorMatrix > tUEqn(fvm::ddt(rho, U)+fvm::div(phi, U)+MRF.DDt(rho, U)+turbulence->divDevRhoReff(U)==fvOptions(rho, U))
volScalarField rhorAtU("rhorAtU", rho *rAtU)
volScalarField rAtU(1.0/(1.0/rAU - UEqn.H1()))
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))))
autoPtr< surfaceVectorField > rhoUf
const dimensionedScalar rhoMin
const dimensionedScalar rhoMax
const pressureControl & pressureControl
tmp< volScalarField > rAU