pEqn.H File Reference

Include dependency graph for pEqn.H:

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

relax()

rho relax

(

)

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

)

HbyA()

volVectorField HbyA

(

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

)

if() [1/8]

if

(

pimple.nCorrPiso()<=

1

)

Definition at line 30 of file pEqn.H.

References constrainPressure(), Foam::fvc::ddtCorr(), Foam::fvc::flux(), HbyA, Foam::fvc::interpolate(), makeRelative(), MRF(), p, 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
	)

constrainPressure()

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

if() [2/8]

if

(

pimple.

transonic()

)

Definition at line 51 of file pEqn.H.

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

Here is the call graph for this function:

while()

while

(

pimple.

correctNonOrthogonal()

)

Definition at line 119 of file pEqn.H.

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

Here is the call graph for this function:

if() [3/8]

if

(

mesh.

schemes).steady(

)

Definition at line 138 of file pEqn.H.

References fvConstraints::constrain(), constrained, fvConstraints, p, psi, psip0(), and thermo.

Here is the call graph for this function:

correctRho()

thermo correctRho

(

psi *p -

psip0

)

if() [4/8]

if

(

constrained

)

Definition at line 149 of file pEqn.H.

References p, rho, and thermo.

correctBoundaryConditions()

U correctBoundaryConditions

(

)

constrain()

fvConstraints constrain

(

U

)

if() [5/8]

if

(

adjustMass &&!thermo.

incompressible()

)

Definition at line 173 of file pEqn.H.

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

Here is the call graph for this function:

if() [6/8]

if

(

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

adjustMass

)

Definition at line 180 of file pEqn.H.

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

Here is the call graph for this function:

if() [7/8]

if

(

mesh.

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

)

Definition at line 188 of file pEqn.H.

References rho.

if() [8/8]

if

(

thermo.

dpdt()

)

Definition at line 193 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.

rAAtU

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

Definition at line 24 of file pEqn.H.

rhorAAtUf

const surfaceScalarField& rhorAAtUf

Initial value:

=
    pimple.consistent() ? rhorAtUf() : rhorAUf

Definition at line 25 of file pEqn.H.

adjustMass

bool adjustMass = false

Definition at line 49 of file pEqn.H.

else

else

Initial value:

{
    if (mesh.schemes().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)
     ==
        fvModels.source(psi, p, rho.name())
    )

Definition at line 99 of file pEqn.H.

U

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

Definition at line 161 of file pEqn.H.

K

K = 0.5*magSqr(U)

Definition at line 164 of file pEqn.H.