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())
const volScalarField	psip0 (psi *p)
const volScalarField	rAU ("rAU", 1.0/UEqn.A())
const surfaceScalarField	rhorAUf ("rhorAUf", fvc::interpolate(rho *rAU))
volVectorField	HbyA (constrainHbyA(rAU *UEqn.H(), U, p_rgh))
	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)
surfaceScalarField	phig (-rhorAUf *ghf *fvc::snGrad(rho) *mesh.magSf())
	constrainPressure (p_rgh, rho, U, phiHbyA, rhorAUf, MRF)
fvScalarMatrix	p_rghEqn (p_rgh, dimMass/dimTime)
	if (pimple.transonic())
	while (pimple.correctNonOrthogonal())
	if (mesh.steady())
thermo	correctRho (psi *p - psip0)
	if (constrained)
p_rgh	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
bool	adjustMass = mesh.steady() && adjustPhi(phiHbyA, U, p_rgh)
	phiHbyA = phig
	else
	phi = phiHbyA + p_rghEqn.flux()
	p = p_rgh + rho*gh + pRef
	U = HbyA + rAU*fvc::reconstruct((phig + p_rghEqn.flux())/rhorAUf)
	K = 0.5*magSqr(U)

Function Documentation

if() [1/9]

if

(

!mesh.

steady) &&!pimple.simpleRho(

)

Definition at line 1 of file pEqn.H.

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

Here is the call graph for this function:

psip0()

const volScalarField psip0

(

psi *

p

)

rAU()

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

rhorAUf()

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

HbyA()

volVectorField HbyA

(

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

)

if() [2/9]

if

(

pimple.nCorrPiso()<=

1

)

Definition at line 14 of file pEqn.H.

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

phig()

surfaceScalarField phig

(

-rhorAUf *ghf *

fvc::snGradrho) *mesh.magSf(

)

Referenced by mixtureKEpsilon< BasicMomentumTransportModel >::correct(), Foam::hydrostaticInitialisation(), if(), and while().

Here is the caller graph for this function:

constrainPressure()

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

p_rghEqn()

fvScalarMatrix p_rghEqn	(	p_rgh	,
		dimMass/	dimTime
	)

Referenced by if(), and while().

Here is the caller graph for this function:

if() [3/9]

if

(

pimple.

transonic()

)

Definition at line 41 of file pEqn.H.

References Foam::correction(), Foam::fvc::ddt(), Foam::fvc::div(), fvModels, Foam::fvc::interpolate(), Foam::fvc::laplacian(), p_rgh, p_rghEqn(), phiHbyA, pimple, pressureReference, psi, and rho.

Here is the call graph for this function:

while()

while

(

pimple.

correctNonOrthogonal()

)

Definition at line 85 of file pEqn.H.

References Foam::fvc::laplacian(), p_rgh, p_rghEqn(), and pressureReference.

Here is the call graph for this function:

if() [4/9]

if

(

mesh.

steady()

)

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

References p_rgh, rho, and thermo.

relax()

p_rgh relax

(

)

correctBoundaryConditions()

U correctBoundaryConditions

(

)

constrain()

fvConstraints constrain

(

U

)

if() [6/9]

if

(

adjustMass &&!thermo.

incompressible()

)

Definition at line 140 of file pEqn.H.

References Foam::fvc::domainIntegrate(), gh, initialMass, p_rgh, pRef, psi, rho, and thermo.

Here is the call graph for this function:

if() [7/9]

if

(

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

adjustMass

)

Definition at line 148 of file pEqn.H.

References rho, and thermo.

if() [8/9]

if

(

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

)

Definition at line 153 of file pEqn.H.

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

Here is the call graph for this function:

if() [9/9]

if

(

thermo.

dpdt()

)

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

adjustMass

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

Definition at line 28 of file pEqn.H.

phiHbyA

phiHbyA = phig

Definition at line 32 of file pEqn.H.

Referenced by Foam::constrainPhiHbyA(), if(), NullMRF::relative(), and while().

else

else

Initial value:

{
    fvScalarMatrix p_rghDDtEqn
    (
        fvc::ddt(rho) + psi*correction(fvm::ddt(p_rgh))
      + fvc::div(phiHbyA)
     ==
        fvModels.source(psi, p_rgh, rho.name())
    )

Definition at line 76 of file pEqn.H.

phi

phi = phiHbyA + p_rghEqn.flux()

Definition at line 99 of file pEqn.H.

Referenced by forAllConstIter(), if(), and while().

p

p = p_rgh + rho*gh + pRef

Definition at line 101 of file pEqn.H.

Referenced by if(), and while().

U

U = HbyA + rAU*fvc::reconstruct((phig + p_rghEqn.flux())/rhorAUf)

Definition at line 128 of file pEqn.H.

K

K = 0.5*magSqr(U)

Definition at line 131 of file pEqn.H.