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())
 
volScalarField rAU ("rAU", 1.0/UEqn.A())
 
surfaceScalarField rhorAUf ("rhorAUf", fvc::interpolate(rho *rAU))
 
volVectorField HbyA (constrainHbyA(rAU *UEqn.H(), U, p_rgh))
 
 if (pimple.nCorrPiso()<=1)
 
surfaceScalarField phig (-rhorAUf *ghf *fvc::snGrad(rho) *mesh.magSf())
 
surfaceScalarField phiHbyA ("phiHbyA", fvc::flux(rho *HbyA)+MRF.zeroFilter(rhorAUf *fvc::ddtCorr(rho, U, phi)))
 
MRF makeRelative (fvc::interpolate(rho), phiHbyA)
 
 constrainPressure (p_rgh, rho, U, phiHbyA, rhorAUf, MRF)
 
 if (pimple.transonic())
 
tmp< volScalarField > psip0 (mesh.steady() ? tmp< volScalarField >() :psi *p)
 
 while (pimple.correctNonOrthogonal())
 
pressureControl limit (p)
 
 if (!mesh.steady())
 
 if (thermo.dpdt())
 

Variables

const bool closedVolume = mesh.steady() && adjustPhi(phiHbyA, U, p_rgh)
 
const bool adjustMass = closedVolume && !thermo.incompressible()
 
 phiHbyA = phig
 
 p = p_rgh + rho*gh
 
 else
 

Function Documentation

◆ if() [1/5]

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

Definition at line 1 of file pEqn.H.

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

Here is the call graph for this function:

◆ rAU()

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

◆ rhorAUf()

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

◆ HbyA()

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

◆ if() [2/5]

if ( pimple.nCorrPiso()<=  1)

Definition at line 9 of file pEqn.H.

References Foam::fvc::ddtCorr(), Foam::fvc::flux(), ghf, Foam::fvc::interpolate(), mesh, MRF, phi, phig(), phiHbyA, rho, Foam::fvc::snGrad(), tUEqn(), and U.

Here is the call graph for this function:

◆ phig()

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

◆ phiHbyA()

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

◆ makeRelative()

MRF makeRelative ( fvc::interpolate(rho ,
phiHbyA   
)

◆ constrainPressure()

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

◆ if() [3/5]

if ( pimple.  transonic())

Definition at line 41 of file pEqn.H.

References Foam::fvc::div(), Foam::fvc::interpolate(), mesh, p, p_rgh, phiHbyA, psi, psip0(), and rho.

Here is the call graph for this function:

◆ psip0()

tmp<volScalarField> psip0 ( mesh.steady() ? tmp< volScalarField >() :psi p)

◆ while()

while ( pimple.  correctNonOrthogonal())

Definition at line 58 of file pEqn.H.

References Foam::fvc::div(), fvOptions, Foam::fvc::laplacian(), Foam::magSqr(), p_rgh, phi, phig(), phiHbyA, pimple, pressureControl, Foam::fvc::reconstruct(), and U.

Here is the call graph for this function:

◆ limit()

pressureControl limit ( p  )

◆ if() [4/5]

if ( !mesh.  steady())

Definition at line 99 of file pEqn.H.

References psi, and thermo.

◆ if() [5/5]

if ( thermo.  dpdt())

Definition at line 106 of file pEqn.H.

References Foam::fvc::ddt(), dpdt, and mesh.

Here is the call graph for this function:

Variable Documentation

◆ closedVolume

const bool closedVolume = mesh.steady() && adjustPhi(phiHbyA, U, p_rgh)

Definition at line 25 of file pEqn.H.

◆ adjustMass

const bool adjustMass = closedVolume && !thermo.incompressible()

Definition at line 26 of file pEqn.H.

◆ phiHbyA

phiHbyA = phig

Definition at line 28 of file pEqn.H.

◆ p

p = p_rgh + rho*gh

Definition at line 94 of file pEqn.H.

◆ else

else
Initial value:
{
}
{
p_rgh = p - rho*gh;
p.correctBoundaryConditions();
}
if (adjustMass || mesh.steady() || pimple.simpleRho())
{
rho = thermo.rho();
}
if (mesh.steady() && !pimple.transonic())
{
rho.relax();
}
Info<< "Min/max rho:" << min(rho).value() << ' '
<< max(rho).value() << endl
rho
Definition: pEqn.H:1
pimpleNoLoopControl & pimple
dimensionedScalar initialMass
Definition: createFields.H:57
p
Definition: pEqn.H:50
dimensioned< Type > max(const dimensioned< Type > &, const dimensioned< Type > &)
p_rgh
Definition: pEqn.H:140
Ostream & endl(Ostream &os)
Add newline and flush stream.
Definition: Ostream.H:256
dimensioned< Type > domainIntegrate(const GeometricField< Type, fvPatchField, volMesh > &vf)
const bool adjustMass
Definition: pEqn.H:26
rhoReactionThermo & thermo
Definition: createFields.H:28
dynamicFvMesh & mesh
dimensioned< Type > min(const dimensioned< Type > &, const dimensioned< Type > &)
const volScalarField & gh
messageStream Info
const volScalarField & psi

Definition at line 118 of file pEqn.H.