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.

Functions

 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 (constrained)
 
p_rgh relax ()
 
U correctBoundaryConditions ()
 
fvConstraints constrain (U)
 
 if (adjustMass &&!thermo.incompressible())
 
 if (mesh.steady()||pimple.simpleRho()||adjustMass)
 
 if ((mesh.steady()||pimple.simpleRho()) &&!pimple.transonic())
 
 if (thermo.dpdt())
 

Variables

bool adjustMass = mesh.steady() && adjustPhi(phiHbyA, U, p_rgh)
 
 phiHbyA = phig
 
 else
 
 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 ,
phiHbyA   
)

◆ phig()

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

Referenced by mixtureKEpsilon< BasicMomentumTransportModel >::correct(), Foam::hydrostaticInitialisation(), if(), and while().

Here is the caller graph for this function:

◆ constrainPressure()

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

◆ 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())

Definition at line 41 of file pEqn.H.

References Foam::correction(), Foam::fvc::ddt(), Foam::fvc::div(), fvModels, Foam::fvc::interpolate(), Foam::fvc::laplacian(), p_rgh, p_rghEqn(), phiHbyA, pimple, pressureReference, psi, and rho.

Here is the call graph for this function:

◆ while()

while ( pimple.  correctNonOrthogonal())

Definition at line 85 of file pEqn.H.

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

Here is the call graph for this function:

◆ if() [4/9]

if ( mesh.  steady())

Definition at line 103 of file pEqn.H.

References fvConstraints::constrain(), fvConstraints, psi, and thermo.

Here is the call graph for this function:

◆ correctRho()

thermo correctRho ( psi *p psip0)

◆ if() [5/9]

if ( constrained  )

Definition at line 114 of file pEqn.H.

References p_rgh, rho, and thermo.

◆ relax()

p_rgh relax ( )

◆ correctBoundaryConditions()

U correctBoundaryConditions ( )

◆ constrain()

fvConstraints constrain ( U  )

◆ if() [6/9]

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

Definition at line 140 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 148 of file pEqn.H.

References rho, and thermo.

◆ if() [8/9]

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

Definition at line 153 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 161 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 = mesh.steady() && 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

else
Initial value:
{
fvScalarMatrix p_rghDDtEqn
(
==
fvModels.source(psi, p_rgh, rho.name())
)
tmp< fvMatrix< Type > > correction(const fvMatrix< Type > &)
Return the correction form of the given matrix.
fvMatrix< scalar > fvScalarMatrix
Definition: fvMatricesFwd.H:42
p_rgh
Definition: pEqn.H:159
tmp< GeometricField< Type, fvPatchField, volMesh > > div(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
Definition: fvcDiv.C:47
tmp< GeometricField< Type, fvPatchField, volMesh > > ddt(const dimensioned< Type > dt, const fvMesh &mesh)
Definition: fvcDdt.C:45
rho
Definition: pEqn.H:1
phiHbyA
Definition: pEqn.H:32
Foam::fvModels & fvModels
const volScalarField & psi

Definition at line 76 of file pEqn.H.

◆ phi

phi = phiHbyA + p_rghEqn.flux()

Definition at line 99 of file pEqn.H.

Referenced by forAllConstIter(), if(), and while().

◆ p

p = p_rgh + rho*gh + pRef

Definition at line 101 of file pEqn.H.

Referenced by if(), and while().

◆ U

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

Definition at line 128 of file pEqn.H.

◆ K

K = 0.5*magSqr(U)

Definition at line 131 of file pEqn.H.