pEqn.H File Reference
Include dependency graph for pEqn.H:

Go to the source code of this file.


const volScalarField psip0 (psi *p)
volScalarField rAU (1.0/UEqn.A())
mesh interpolate (rAU)
surfaceScalarField rhorAUf ("rhorAUf", fvc::interpolate(rho *rAU))
volVectorField HbyA (constrainHbyA(rAU *UEqn.H(), U, p_rgh))
surfaceScalarField phig (-rhorAUf *ghf *fvc::snGrad(rho) *mesh.magSf())
surfaceScalarField phiHbyA ("phiHbyA", fvc::flux(rho *HbyA)+phig)
 if (ddtCorr)
MRF makeRelative (fvc::interpolate(rho), phiHbyA)
 constrainPressure (p_rgh, rho, U, phiHbyA, rhorAUf, MRF)
fvScalarMatrix p_rghDDtEqn (fvc::ddt(rho)+psi *correction(fvm::ddt(p_rgh))+fvc::div(phiHbyA)==fvOptions(psi, p_rgh,
 while (pimple.correctNonOrthogonal())
 if (p_rgh.needReference())
 if (thermo.dpdt())
surfaceScalarField faceMask (localMin< scalar >(mesh).interpolate(cellMask))


dimensionedScalar compressibility = fvc::domainIntegrate(psi)
bool compressible = (compressibility.value() > SMALL)
 rho = thermo.rho()
 p = p_rgh + rho*gh

Function Documentation

◆ psip0()

const volScalarField psip0 ( psi p)

◆ rAU()

volScalarField rAU ( 1.0/UEqn.  A())

◆ interpolate()

mesh interpolate ( rAU  )

◆ rhorAUf()

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

◆ HbyA()

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

◆ phig()

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

◆ phiHbyA()

surfaceScalarField phiHbyA ( "phiHbyA"  ,
fvc::flux(rho *HbyA)+  phig 

◆ if() [1/3]

if ( ddtCorr  )

Definition at line 24 of file pEqn.H.

References ddtCorr, interpolate(), mesh, MRF, phi, phiHbyA, rho, rhorAUf(), and U.

Here is the call graph for this function:

◆ makeRelative()

MRF makeRelative ( fvc::interpolate(rho ,

◆ constrainPressure()

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

◆ p_rghDDtEqn()

fvScalarMatrix p_rghDDtEqn ( fvc::ddt(rho)+psi *correction(fvm::ddt(p_rgh))+fvc::div(phiHbyA = =fvOptions(psip_rgh,

◆ while()

while ( pimple.  correctNonOrthogonal())

Definition at line 48 of file pEqn.H.

References fvOptions, HbyA, Foam::fac::laplacian(), Foam::magSqr(), mesh, p_rgh, p_rghDDtEqn(), phi, phig(), phiHbyA, pimple, rAU(), Foam::fvc::reconstruct(), rhorAUf(), and U.

Here is the call graph for this function:

◆ if() [2/3]

if ( p_rgh.  needReference())

Definition at line 84 of file pEqn.H.

References compressibility, compressible, Foam::fvc::domainIntegrate(), Foam::getRefCellValue(), gh, initialMass, p, p_rgh, pRefCell, pRefValue, psi, psip0(), and rho.

Here is the call graph for this function:

◆ if() [3/3]

if ( thermo.  dpdt())

Definition at line 115 of file pEqn.H.

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

Here is the call graph for this function:

◆ faceMask()

surfaceScalarField faceMask ( localMin< scalar >  mesh).interpolate(cellMask)

Variable Documentation

◆ compressibility

dimensionedScalar compressibility = fvc::domainIntegrate(psi)

Definition at line 1 of file pEqn.H.

Referenced by if().

◆ compressible

bool compressible = (compressibility.value() > SMALL)

Definition at line 2 of file pEqn.H.

Referenced by turbulenceFields::execute(), for(), if(), and while().

◆ rho

rho = thermo.rho()

Definition at line 4 of file pEqn.H.

◆ p

p = p_rgh + rho*gh

Definition at line 79 of file pEqn.H.

◆ else

Initial value:
thermo.correctRho(psi*p - psip0)

Definition at line 105 of file pEqn.H.

const volScalarField psip0(psi *p)
Basic thermodynamics type based on the use of fitting functions for cp, h, s obtained from the templa...
Definition: pEqn.H:50
const volScalarField & psi
Definition: createFieldRefs.H:1