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
rho	relax ()
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)
const 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.schemes().steady())
thermo	correctRho (psi *p - psip0)
	if (constrained)
U	correctBoundaryConditions ()
fvConstraints	constrain (U)
	if (adjustMass &&!thermo.incompressible())
	if (mesh.schemes().steady()||pimple.simpleRho()||adjustMass)
	if (mesh.schemes().steady()||pimple.simpleRho())
	if (thermo.dpdt())

Variables
	rho = thermo.rho()
bool	adjustMass = mesh.schemes().steady() && adjustPhi(phiHbyA, U, p_rgh)
	phiHbyA = phig
	else
	phi = phiHbyA + p_rghEqn.flux()
const bool	constrained = fvConstraints.constrain(p)
	p = p_rgh + rho*gh + pRef
	U = HbyA + rAU*fvc::reconstruct((phig + p_rghEqn.flux())/rhorAUf)
	K = 0.5*magSqr(U)

Function Documentation

relax()

rho relax

(

)

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() [1/8]

if

(

pimple.nCorrPiso()<=

1

)

Definition at line 12 of file pEqn.H.

References Foam::fvc::ddtCorr(), Foam::fvc::flux(), HbyA, Foam::fvc::interpolate(), MRF(), phi, phiHbyA, rho, rhorAUf(), 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()

const 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() [2/8]

if

(

pimple.

transonic()

)

Definition at line 39 of file pEqn.H.

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

Here is the call graph for this function:

while()

while

(

pimple.

correctNonOrthogonal()

)

Definition at line 84 of file pEqn.H.

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

Here is the call graph for this function:

if() [3/8]

if

(

mesh.

schemes).steady(

)

Definition at line 100 of file pEqn.H.

References gh, p, p_rgh, pRef, and rho.

correctRho()

thermo correctRho

(

psi *p -

psip0

)

if() [4/8]

if

(

constrained

)

Definition at line 113 of file pEqn.H.

References p_rgh, rho, and thermo.

correctBoundaryConditions()

U correctBoundaryConditions

(

)

constrain()

fvConstraints constrain

(

U

)

if() [5/8]

if

(

adjustMass &&!thermo.

incompressible()

)

Definition at line 145 of file pEqn.H.

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

Here is the call graph for this function:

if() [6/8]

if

(

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

adjustMass

)

Definition at line 153 of file pEqn.H.

References rho, and thermo.

if() [7/8]

if

(

mesh.

schemes).steady()||pimple.simpleRho(

)

Definition at line 158 of file pEqn.H.

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

Here is the call graph for this function:

if() [8/8]

if

(

thermo.

dpdt()

)

Definition at line 166 of file pEqn.H.

References Foam::fvc::ddt(), Foam::fvc::div(), dpdt, mesh, Foam::fvc::meshPhi(), p, rho, and U.

Here is the call graph for this function:

Variable Documentation

rho

rho = thermo.rho()

Definition at line 1 of file pEqn.H.

adjustMass

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

Definition at line 26 of file pEqn.H.

phiHbyA

phiHbyA = phig

Definition at line 30 of file pEqn.H.

Referenced by Foam::constrainPhiHbyA(), if(), 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 75 of file pEqn.H.

phi

phi = phiHbyA + p_rghEqn.flux()

Definition at line 98 of file pEqn.H.

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

constrained

const bool constrained = fvConstraints.constrain(p)

Definition at line 108 of file pEqn.H.

Referenced by fvConstraints::constrain(), limitPressure::constrain(), indexedCell< Gt, Cb >::featurePointExternalCell(), indexedCell< Gt, Cb >::featurePointInternalCell(), indexedCell< Gt, Cb >::hasConstrainedPoint(), and if().

p

p = p_rgh + rho*gh + pRef

Definition at line 125 of file pEqn.H.

Referenced by if(), and while().

U

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

Definition at line 129 of file pEqn.H.

K

K = 0.5*magSqr(U)

Definition at line 132 of file pEqn.H.