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 (limitedp)
p	relax ()
U	correctBoundaryConditions ()
fvOptions	correct (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(), and while().

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

)

rhorAtUf()

tmp<surfaceScalarField> rhorAtUf

(

pimple.

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

)

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

)

Definition at line 52 of file pEqn.H.

References Foam::correction(), Foam::fvc::ddt(), Foam::fvc::div(), fvOptions, Foam::fvc::grad(), HbyA, Foam::fvc::interpolate(), Foam::fvc::laplacian(), mesh, phi, phiHbyA, pimple, pressureControl, psi, rAU, rho, rhorAUf(), and Foam::fvc::snGrad().

Here is the call graph for this function:

while()

while

(

pimple.

correctNonOrthogonal()

)

Definition at line 118 of file pEqn.H.

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

Here is the call graph for this function:

if() [4/9]

if

(

mesh.

steady()

)

Definition at line 137 of file pEqn.H.

References limitedp, pressureControl, psi, and thermo.

correctRho()

thermo correctRho

(

psi *p -

psip0

)

if() [5/9]

if

(

limitedp

)

Definition at line 148 of file pEqn.H.

References rho, and thermo.

relax()

p relax

(

)

correctBoundaryConditions()

U correctBoundaryConditions

(

)

correct()

fvOptions correct

(

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

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())
    {
        adjustMass = adjustPhi(phiHbyA, U, p);
    }

    if (pimple.consistent())
    {
        phiHbyA += (rhorAAtUf - rhorAUf)*fvc::snGrad(p)*mesh.magSf();
        HbyA += (rAAtU - rAU)*fvc::grad(p);
    }

    fvScalarMatrix pDDtEqn
    (
        fvc::ddt(rho) + psi*correction(fvm::ddt(p))
      + fvc::div(phiHbyA)
     ==
        fvOptions(psi, p, rho.name())
    )

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.