27 while (
simple.correctNonOrthogonal())
29 fvScalarMatrix p_rghEqn
38 if (
simple.finalNonOrthogonalIter())
49 U.correctBoundaryConditions();
56 #include "continuityErrs.H"
64 p += dimensionedScalar
74 /fvc::domainIntegrate(
psi);
82 Info<<
"rho min/max : " << min(
rho).value() <<
" " <<
max(
rho).value()
const surfaceScalarField & ghf
const volScalarField & gh
Basic thermodynamics type based on the use of fitting functions for cp, h, s obtained from the templa...
const volScalarField & psi
surfaceScalarField phig("phig", -rhorAUf *ghf *fvc::snGrad(rho) *mesh.magSf())
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))
adjustPhi(phiHbyA, U, p_rgh)
dimensionedScalar compressibility
tmp< volScalarField > rAU
GeometricField< scalar, fvPatchField, volMesh > volScalarField
const dictionary & simple
dimensionedScalar initialMass