pEqn.H File Reference
Include dependency graph for pEqn.H:
This graph shows which files directly or indirectly include this file:

Go to the source code of this file.


 if (!mesh.steady() &&!pimple.simpleRho())
const volScalarField psip0 (psi *p)
const volScalarField rAU ("rAU", 1.0/UEqn.A())
const surfaceScalarField rhorAUf ("rhorAUf", fvc::interpolate(rho *rAU))
volVectorField HbyA (constrainHbyA(rAU *UEqn.H(), U, p_rgh))
 if (pimple.nCorrPiso()<=1)
surfaceScalarField phiHbyA ("phiHbyA", fvc::interpolate(rho) *fvc::flux(HbyA)+MRF.zeroFilter(rhorAUf *fvc::ddtCorr(rho, U, phi, rhoUf)))
MRF makeRelative (fvc::interpolate(rho), phiHbyA)
surfaceScalarField phig (-rhorAUf *ghf *fvc::snGrad(rho) *mesh.magSf())
 constrainPressure (p_rgh, rho, U, phiHbyA, rhorAUf, MRF)
fvScalarMatrix p_rghEqn (p_rgh, dimMass/dimTime)
 if (pimple.transonic())
 while (pimple.correctNonOrthogonal())
 if (mesh.steady())
thermo correctRho (psi *p - psip0)
 if (limitedp)
p_rgh relax ()
U correctBoundaryConditions ()
fvOptions correct (U)
 if (adjustMass &&!thermo.incompressible())
 if (mesh.steady()||pimple.simpleRho()||adjustMass)
 if ((mesh.steady()||pimple.simpleRho()) &&!pimple.transonic())
 if (thermo.dpdt())


bool adjustMass = pimple.transonic() ? false : adjustPhi(phiHbyA, U, p_rgh)
 phiHbyA = phig
 phi = phiHbyA + p_rghEqn.flux()
 p = p_rgh + rho*gh + pRef
 U = HbyA + rAU*fvc::reconstruct((phig + p_rghEqn.flux())/rhorAUf)
 K = 0.5*magSqr(U)

Function Documentation

◆ if() [1/9]

if ( !mesh.  steady) &&!pimple.simpleRho()

Definition at line 1 of file pEqn.H.

References Foam::constrainHbyA(), HbyA, Foam::fvc::interpolate(), p, p_rgh, psi, psip0(), rAU, rho, rhorAUf(), thermo, U, and UEqn.

Here is the call graph for this function:

◆ psip0()

const volScalarField psip0 ( psi p)

◆ rAU()

const volScalarField rAU ( "rAU"  ,
1.0/UEqn.  A() 

◆ rhorAUf()

const surfaceScalarField rhorAUf ( "rhorAUf"  ,
fvc::interpolate(rho *rAU  

◆ HbyA()

volVectorField HbyA ( constrainHbyA(rAU *UEqn.H(), U, p_rgh )

◆ if() [2/9]

if ( pimple.nCorrPiso()<=  1)

Definition at line 14 of file pEqn.H.

References Foam::fvc::ddtCorr(), Foam::fvc::flux(), Foam::fvc::interpolate(), MRF, phi, phiHbyA, rho, rhoUf, tUEqn(), and U.

Here is the call graph for this function:

◆ phiHbyA()

surfaceScalarField phiHbyA ( "phiHbyA"  ,
fvc::interpolate(rho) *fvc::flux(HbyA)+MRF.zeroFilter(rhorAUf *fvc::ddtCorr(rho, U, phi, rhoUf))   

◆ makeRelative()

MRF makeRelative ( fvc::interpolate(rho ,

◆ phig()

surfaceScalarField phig ( -rhorAUf *ghf fvc::snGradrho) *mesh.magSf()

◆ constrainPressure()

constrainPressure ( p_rgh  ,
rho  ,
U  ,
phiHbyA  ,
rhorAUf  ,

◆ p_rghEqn()

fvScalarMatrix p_rghEqn ( p_rgh  ,
dimMass/  dimTime 

Referenced by if(), and while().

Here is the caller graph for this function:

◆ if() [3/9]

if ( pimple.  transonic())

◆ while()

while ( pimple.  correctNonOrthogonal())

Definition at line 90 of file pEqn.H.

References Foam::fvc::laplacian(), p_rgh, p_rghEqn(), and pressureControl.

Here is the call graph for this function:

◆ if() [4/9]

if ( mesh.  steady())

Definition at line 108 of file pEqn.H.

References limitedp, pressureControl, psi, and thermo.

◆ correctRho()

thermo correctRho ( psi *p psip0)

◆ if() [5/9]

if ( limitedp  )

Definition at line 119 of file pEqn.H.

References p_rgh, rho, and thermo.

◆ relax()

p_rgh relax ( )

◆ correctBoundaryConditions()

U correctBoundaryConditions ( )

◆ correct()

fvOptions correct ( U  )

◆ if() [6/9]

if ( adjustMass &&!thermo.  incompressible())

Definition at line 145 of file pEqn.H.

References Foam::fvc::domainIntegrate(), gh, initialMass, p_rgh, pRef, psi, rho, and thermo.

Here is the call graph for this function:

◆ if() [7/9]

if ( mesh.steady()||pimple.simpleRho()||  adjustMass)

Definition at line 153 of file pEqn.H.

References rho, and thermo.

◆ if() [8/9]

if ( (mesh.steady()||pimple.simpleRho()) &&!pimple.transonic()  )

Definition at line 158 of file pEqn.H.

References Foam::fvc::correctRhoUf(), phi, rho, rhoUf, and U.

Here is the call graph for this function:

◆ if() [9/9]

if ( thermo.  dpdt())

Definition at line 166 of file pEqn.H.

References Foam::fvc::ddt(), Foam::fvc::div(), dpdt, mesh, Foam::fvc::meshPhi(), rho, and U.

Here is the call graph for this function:

Variable Documentation

◆ adjustMass

bool adjustMass = pimple.transonic() ? false : adjustPhi(phiHbyA, U, p_rgh)

Definition at line 28 of file pEqn.H.

◆ phiHbyA

phiHbyA = phig

Definition at line 32 of file pEqn.H.

Referenced by Foam::constrainPhiHbyA(), if(), NullMRF::relative(), and while().

◆ else

Initial value:
if (mesh.steady())
fvOptions(psi, p_rgh,
tmp< fvMatrix< Type > > correction(const fvMatrix< Type > &)
Return the correction form of the given matrix.
fvScalarMatrix p_rghDDtEqn(fvc::ddt(rho)+psi *correction(fvm::ddt(p_rgh))+fvc::div(phiHbyA)==parcels.Srho()+surfaceFilm.Srho()+fvOptions(psi, p_rgh,
fvMatrix< scalar > fvScalarMatrix
Definition: fvMatricesFwd.H:42
fv::options & fvOptions
Definition: pEqn.H:1
tmp< GeometricField< Type, fvPatchField, volMesh > > div(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
Definition: fvcDiv.C:47
Definition: pEqn.H:140
tmp< GeometricField< Type, fvPatchField, volMesh > > ddt(const dimensioned< Type > dt, const fvMesh &mesh)
Definition: fvcDdt.C:45
dynamicFvMesh & mesh
bool adjustMass
Definition: pEqn.H:47
Definition: pEqn.H:32
adjustPhi(phiHbyA, Urel, p)
Definition: pEqn.H:72
const volScalarField & psi

Definition at line 76 of file pEqn.H.

◆ phi

phi = phiHbyA + p_rghEqn.flux()

Definition at line 104 of file pEqn.H.

Referenced by ISAT< CompType, ThermoType >::add(), effectivenessHeatExchangerSource::addSup(), waveTransmissiveFvPatchField< Type >::advectionSpeed(), advectiveFvPatchField< Type >::advectionSpeed(), Implicit< CloudType >::cacheFields(), COxidationIntrinsicRate< CloudType >::calculate(), chemPointISAT< CompType, ThermoType >::checkSolution(), VariableHardSphere< CloudType >::collide(), LarsenBorgnakkeVariableHardSphere< CloudType >::collide(), relaxation::correct(), kEqn< BasicMomentumTransportModel >::correct(), kOmega< BasicMomentumTransportModel >::correct(), dynamicKEqn< BasicMomentumTransportModel >::correct(), kEpsilon< compressible::momentumTransportModel >::correct(), LaunderSharmaKE< BasicMomentumTransportModel >::correct(), RNGkEpsilon< BasicMomentumTransportModel >::correct(), realizableKE< BasicMomentumTransportModel >::correct(), kEpsilonLopesdaCosta< BasicMomentumTransportModel >::correct(), mixtureKEpsilon< BasicMomentumTransportModel >::correct(), v2f< BasicMomentumTransportModel >::correct(), kOmegaSST< eddyViscosity< RASModel< BasicMomentumTransportModel > >, BasicMomentumTransportModel >::correct(), waxSolventEvaporation::correctModel(), Foam::fvc::ddt(), Foam::fvc::ddtCorr(), Stokes5::elevation(), fluxCorrectedVelocityFvPatchVectorField::evaluate(), age::execute(), PatchFlowRateInjection< CloudType >::flowRate(), Foam::fvc::flux(), for(), forAll(), forAllConstIter(), chemPointISAT< CompType, ThermoType >::grow(), if(), Foam::incGammaRatio_Q(), chemPointISAT< CompType, ThermoType >::inEOA(), Square< Type >::integrate(), interpolationCellPointFace< Type >::interpolate(), Foam::MULES::limit(), magSqr< Type >::magSqr(), backwardDdtScheme< Type >::meshPhi(), CrankNicolsonDdtScheme< Type >::meshPhi(), MomentumTransportModel< volScalarField, geometricOneField, incompressibleMomentumTransportModel, TransportModel >::MomentumTransportModel(), extendedUpwindCellToFaceStencil::neiStencil(), viscosityModel::New(), laminarModel< BasicMomentumTransportModel >::New(), RASModel< phaseCompressibleMomentumTransportModel >::New(), LESModel< BasicMomentumTransportModel >::New(), null< Type >::operator()(), magSqr< Type >::operator()(), rhoMagSqr< Type >::operator()(), dynamicRefineFvMesh::refine(), compressibleMomentumTransportModel::rho(), rhoMagSqr< Type >::rhoMagSqr(), thermalBaffle::solveEnergy(), EulerCoordinateRotation::transformVector(), dynamicRefineFvMesh::unrefine(), ETAB< CloudType >::update(), TAB< CloudType >::update(), totalFlowRateAdvectiveDiffusiveFvPatchScalarField::updateCoeffs(), outletMappedUniformInletFvPatchField< Type >::updateCoeffs(), mappedFlowRateFvPatchVectorField::updateCoeffs(), pressureInletVelocityFvPatchVectorField::updateCoeffs(), swirlFlowRateInletVelocityFvPatchVectorField::updateCoeffs(), pressureInletOutletParSlipVelocityFvPatchVectorField::updateCoeffs(), pressureNormalInletOutletVelocityFvPatchVectorField::updateCoeffs(), waveSurfacePressureFvPatchScalarField::updateCoeffs(), pressureDirectedInletVelocityFvPatchVectorField::updateCoeffs(), pressureDirectedInletOutletVelocityFvPatchVectorField::updateCoeffs(), porousBafflePressureFvPatchField< Type >::updateCoeffs(), syringePressureFvPatchScalarField::updateCoeffs(), plenumPressureFvPatchScalarField::updateCoeffs(), Sine< Type >::value(), Square< Type >::value(), Implicit< CloudType >::velocityCorrection(), Airy::vi(), while(), Gulders::~Gulders(), and GuldersEGR::~GuldersEGR().

◆ p

p = p_rgh + rho*gh + pRef

Definition at line 106 of file pEqn.H.

◆ U

U = HbyA + rAU*fvc::reconstruct((phig + p_rghEqn.flux())/rhorAUf)

Definition at line 133 of file pEqn.H.

◆ K

K = 0.5*magSqr(U)

Definition at line 136 of file pEqn.H.