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

Go to the source code of this file.


 if (!pimple.simpleRho())
const volScalarField psip0 (psi *p)
volScalarField rAU (1.0/UEqn.A())
surfaceScalarField rhorAUf ("rhorAUf", fvc::interpolate(rho *rAU))
volVectorField HbyA (constrainHbyA(rAU *UEqn.H(), U, p))
 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)
 constrainPressure (p, rho, U, phiHbyA, rhorAUf, MRF)
 if (pimple.transonic())
 while (pimple.correctNonOrthogonal())
thermo correctRho (psi *p - psip0)
 if (limitedp)
p relax ()
U correctBoundaryConditions ()
fvOptions correct (U)
 if (thermo.dpdt())


bool limitedp = pressureControl.limit(p)
 U = HbyA - rAU*fvc::grad(p)
 K = 0.5*magSqr(U)

Function Documentation

◆ if() [1/5]

if ( !pimple.  simpleRho())

Definition at line 1 of file pEqn.H.

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

Here is the call graph for this function:

◆ psip0()

const volScalarField psip0 ( psi p)

◆ rAU()

volScalarField rAU ( 1.0/UEqn.  A())

◆ rhorAUf()

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

◆ HbyA()

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

◆ if() [2/5]

if ( pimple.nCorrPiso()<=  1)

Definition at line 14 of file pEqn.H.

References constrainPressure(), Foam::fvc::ddtCorr(), Foam::fvc::flux(), Foam::fvc::interpolate(), makeRelative(), 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 ,

◆ constrainPressure()

constrainPressure ( p  ,
rho  ,
U  ,
phiHbyA  ,
rhorAUf  ,

◆ if() [3/5]

if ( pimple.  transonic())

Definition at line 32 of file pEqn.H.

References Foam::correction(), Foam::fvc::ddt(), Foam::fvc::div(), fvOptions, Foam::fvc::interpolate(), Foam::fvc::laplacian(), pDDtEqn(), phi, phiHbyA, pimple, psi, and rho.

Here is the call graph for this function:

◆ while()

while ( pimple.  correctNonOrthogonal())

Definition at line 75 of file pEqn.H.

References Foam::fvc::laplacian(), phi, phiHbyA, and pimple.

Here is the call graph for this function:

◆ correctRho()

thermo correctRho ( psi *p psip0)

◆ if() [4/5]

if ( limitedp  )

Definition at line 93 of file pEqn.H.

References rho, and thermo.

◆ relax()

p relax ( )

◆ correctBoundaryConditions()

U correctBoundaryConditions ( )

◆ correct()

fvOptions correct ( U  )

◆ if() [5/5]

if ( thermo.  dpdt())

Definition at line 117 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

◆ else

Initial value:
fvOptions(psi, p,
tmp< fvMatrix< Type > > correction(const fvMatrix< Type > &)
Return the correction form of the given matrix.
Definition: pEqn.H:22
fvMatrix< scalar > fvScalarMatrix
Definition: fvMatricesFwd.H:42
fvScalarMatrix pDDtEqn(fvc::ddt(rho)+psi *correction(fvm::ddt(p))+fvc::div(phiHbyA)==fvOptions(psi, p,
fv::options & fvOptions
Definition: pEqn.H:1
Definition: pEqn.H:50
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
const volScalarField & psi

Definition at line 66 of file pEqn.H.

◆ limitedp

bool limitedp = pressureControl.limit(p)

Definition at line 88 of file pEqn.H.

◆ U

U = HbyA - rAU*fvc::grad(p)

Definition at line 104 of file pEqn.H.

◆ K

K = 0.5*magSqr(U)

Definition at line 107 of file pEqn.H.