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

Go to the source code of this file.

Functions

 if ((!mesh.steady() &&!pimple.simpleRho())||pimple.consistent())
 
const volScalarField psip0 (psi *p)
 
const volScalarField rAU ("rAU", 1.0/UEqn.A())
 
const surfaceScalarField rhorAUf ("rhorAUf", fvc::interpolate(rho *rAU))
 
tmp< volScalarField > rAtU (pimple.consistent() ? volScalarField::New("rAtU", 1.0/(1.0/rAU - UEqn.H1())) :tmp< volScalarField >(nullptr))
 
tmp< surfaceScalarField > rhorAtUf (pimple.consistent() ? surfaceScalarField::New("rhoRAtUf", fvc::interpolate(rho *rAtU())) :tmp< surfaceScalarField >(nullptr))
 
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, rhorAAtUf, MRF)
 
 if (pimple.transonic())
 
 while (pimple.correctNonOrthogonal())
 
 if (mesh.steady())
 
thermo correctRho (psi *p - psip0)
 
 if (constrained)
 
p 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

const volScalarField & rAAtU = pimple.consistent() ? rAtU() : rAU
 
const surfaceScalarField & rhorAAtUf
 
bool adjustMass = false
 
 else
 
 U = HbyA - rAAtU*fvc::grad(p)
 
 K = 0.5*magSqr(U)
 

Function Documentation

◆ if() [1/9]

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

Definition at line 1 of file pEqn.H.

References Foam::fvc::interpolate(), Foam::compressible::New(), p, pimple, psi, psip0(), rAtU(), rAU, rho, rhorAtUf(), rhorAUf(), thermo, and UEqn.

Here is the call graph for this function:

◆ psip0()

const volScalarField psip0 ( psi p)

Referenced by if().

Here is the caller graph for this function:

◆ rAU()

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

◆ rhorAUf()

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

◆ rAtU()

tmp<volScalarField> rAtU ( pimple.  consistent) ? volScalarField::New("rAtU", 1.0/(1.0/rAU - UEqn.H1())) :tmp< volScalarField >(nullptr)

Referenced by if(), and while().

Here is the caller graph for this function:

◆ rhorAtUf()

tmp<surfaceScalarField> rhorAtUf ( pimple.  consistent) ? surfaceScalarField::New("rhoRAtUf", fvc::interpolate(rho *rAtU())) :tmp< surfaceScalarField >(nullptr)

Referenced by if().

Here is the caller graph for this function:

◆ HbyA()

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

◆ if() [2/9]

if ( pimple.nCorrPiso()<=  1)

Definition at line 32 of file pEqn.H.

References Foam::fvc::ddtCorr(), Foam::fvc::flux(), HbyA, 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 ,
phiHbyA   
)

◆ constrainPressure()

constrainPressure ( p  ,
rho  ,
U  ,
phiHbyA  ,
rhorAAtUf  ,
MRF   
)

◆ if() [3/9]

if ( pimple.  transonic())

◆ while()

while ( pimple.  correctNonOrthogonal())

Definition at line 118 of file pEqn.H.

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

Here is the call graph for this function:

◆ if() [4/9]

if ( mesh.  steady())

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

References rho, and thermo.

◆ relax()

p relax ( )

◆ correctBoundaryConditions()

U correctBoundaryConditions ( )

◆ constrain()

fvConstraints constrain ( U  )

◆ if() [6/9]

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

Definition at line 172 of file pEqn.H.

References Foam::fvc::domainIntegrate(), initialMass, psi, and thermo.

Here is the call graph for this function:

◆ if() [7/9]

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

Definition at line 179 of file pEqn.H.

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

Here is the call graph for this function:

◆ if() [8/9]

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

Definition at line 187 of file pEqn.H.

References rho.

◆ if() [9/9]

if ( thermo.  dpdt())

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

◆ rAAtU

const volScalarField& rAAtU = pimple.consistent() ? rAtU() : rAU

Definition at line 26 of file pEqn.H.

◆ rhorAAtUf

const surfaceScalarField& rhorAAtUf
Initial value:
=
pimple.consistent() ? rhorAtUf() : rhorAUf
pimpleNoLoopControl & pimple
surfaceScalarField rhorAUf("rhorAUf", fvc::interpolate(rho *rAU))
tmp< surfaceScalarField > rhorAtUf(pimple.consistent() ? surfaceScalarField::New("rhoRAtUf", fvc::interpolate(rho *rAtU())) :tmp< surfaceScalarField >(nullptr))

Definition at line 27 of file pEqn.H.

◆ adjustMass

bool adjustMass = false

Definition at line 47 of file pEqn.H.

◆ else

else
Initial value:
{
if (mesh.steady())
{
}
if (pimple.consistent())
{
}
(
==
fvModels.source(psi, p, rho.name())
)
tmp< fvMatrix< Type > > correction(const fvMatrix< Type > &)
Return the correction form of the given matrix.
rAU
Definition: pEqn.H:1
tmp< GeometricField< typename outerProduct< vector, Type >::type, fvPatchField, volMesh >> grad(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
Definition: fvcGrad.C:52
fvMatrix< scalar > fvScalarMatrix
Definition: fvMatricesFwd.H:42
pimpleNoLoopControl & pimple
surfaceScalarField rhorAUf("rhorAUf", fvc::interpolate(rho *rAU))
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
dynamicFvMesh & mesh
rho
Definition: pEqn.H:1
const volScalarField & rAAtU
Definition: pEqn.H:26
bool adjustMass
Definition: pEqn.H:47
phiHbyA
Definition: pEqn.H:32
Foam::fvModels & fvModels
adjustPhi(phiHbyA, Urel, p)
U
Definition: pEqn.H:72
const surfaceScalarField & rhorAAtUf
Definition: pEqn.H:27
volVectorField & HbyA
Definition: pEqn.H:13
const volScalarField & psi
p
Definition: pEqn.H:101
tmp< GeometricField< Type, fvsPatchField, surfaceMesh > > snGrad(const GeometricField< Type, fvPatchField, volMesh > &vf, const word &name)
Definition: fvcSnGrad.C:45

Definition at line 98 of file pEqn.H.

◆ U

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

Definition at line 160 of file pEqn.H.

◆ K

K = 0.5*magSqr(U)

Definition at line 163 of file pEqn.H.