17 bool closedVolume =
false;
29 while (
simple.correctNonOrthogonal())
31 tmp<fvScalarMatrix> tpEqn;
58 if (
simple.finalNonOrthogonalIter())
78 U.correctBoundaryConditions();
93 Info<<
"rho max/min : "
tmp< volTensorField > trTU
tmp< GeometricField< Type, fvPatchField, volMesh > > laplacian(const GeometricField< Type, fvPatchField, volMesh > &vf, const word &name)
tmp< surfaceScalarField > interpolate(const RhoType &rho)
fv::IOoptionList & fvOptions
Ostream & endl(Ostream &os)
Add newline and flush stream.
fvVectorMatrix UEqn(fvm::ddt(rho, U)+fvm::div(phi, U)+turbulence->divDevRhoReff(U)==fvOptions(rho, U))
GeometricField< scalar, fvPatchField, volMesh > volScalarField
dimensioned< Type > max(const dimensioned< Type > &, const dimensioned< Type > &)
GeometricField< scalar, fvsPatchField, surfaceMesh > surfaceScalarField
Switch pressureImplicitPorosity(false)
GeometricField< vector, fvPatchField, volMesh > volVectorField
adjustPhi(phiHbyA, U, p_rgh)
dimensionedScalar initialMass
PtrList< dimensionedScalar > rhoMin(fluidRegions.size())
dimensioned< Type > domainIntegrate(const GeometricField< Type, fvPatchField, volMesh > &vf)
const dictionary & simple
tmp< GeometricField< Type, fvPatchField, volMesh > > div(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
tmp< GeometricField< typename outerProduct< vector, Type >::type, fvPatchField, volMesh >> grad(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
psiReactionThermo & thermo
dimensioned< Type > min(const dimensioned< Type > &, const dimensioned< Type > &)
PtrList< dimensionedScalar > rhoMax(fluidRegions.size())
tmp< volScalarField > trAU
Calculates and prints the continuity errors.