pEqn.H File Reference

Include dependency graph for pEqn.H:

Go to the source code of this file.

Functions
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))
tUEqn	clear ()
surfaceScalarField	phiHbyA ("phiHbyA", fvc::interpolate(rho) *fvc::flux(HbyA))
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 (simple.transonic())
p_rgh	relax ()
U	correctBoundaryConditions ()
fvOptions	correct (U)
pressureControl	limit (p)
	if (closedVolume &&!thermo.incompressible())
	if (!simple.transonic())

Variables
	rho = thermo.rho()
bool	closedVolume = simple.transonic() ? false : adjustPhi(phiHbyA, U, p_rgh)
	phiHbyA = phig
	else
	p = p_rgh + rho*gh + pRef
	U = HbyA + rAU*fvc::reconstruct((phig + p_rghEqn.flux())/rhorAUf)

Function Documentation

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)

)

clear()

tUEqn clear

(

)

phiHbyA()

surfaceScalarField phiHbyA	(	"phiHbyA"	,
		fvc::interpolate(rho) *fvc::flux(HbyA)
	)

makeRelative()

MRF makeRelative	(	fvc::interpolate(rho)	,
		phiHbyA
	)

phig()

surfaceScalarField phig

(

-rhorAUf *ghf *

fvc::snGradrho) *mesh.magSf(

)

constrainPressure()

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

p_rghEqn()

fvScalarMatrix p_rghEqn	(	p_rgh	,
		dimMass/	dimTime
	)

if() [1/3]

if

(

simple.

transonic()

)

Definition at line 28 of file pEqn.H.

References Foam::fvc::div(), fvOptions, Foam::fvc::interpolate(), Foam::fvc::laplacian(), p_rgh, p_rghEqn(), phiHbyA, pressureControl, psi, rho, rhorAUf(), and simple().

Here is the call graph for this function:

relax()

p_rgh relax

(

)

correctBoundaryConditions()

U correctBoundaryConditions

(

)

correct()

fvOptions correct

(

U

)

limit()

pressureControl limit

(

p

)

if() [2/3]

if

(

closedVolume &&!thermo.

incompressible()

)

Definition at line 97 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() [3/3]

if

(

!simple.

transonic()

)

Definition at line 107 of file pEqn.H.

References rho.

Variable Documentation

rho

rho = thermo.rho()

Definition at line 1 of file pEqn.H.

closedVolume

bool closedVolume = simple.transonic() ? false : adjustPhi(phiHbyA, U, p_rgh)

Definition at line 17 of file pEqn.H.

phiHbyA

phiHbyA = phig

Definition at line 21 of file pEqn.H.

else

else

Initial value:

{
    while (simple.correctNonOrthogonal())
    {
        p_rghEqn =
            fvc::div(phiHbyA)
          - fvm::laplacian(rhorAUf, p_rgh)
          ==
            fvOptions(psi, p_rgh, rho.name());

        p_rghEqn.setReference
        (
            pressureControl.refCell(),
            pressureControl.refValue()
        );

        p_rghEqn.solve();
    }
}

phi = phiHbyA + p_rghEqn.flux()

Definition at line 59 of file pEqn.H.

p

p = p_rgh + rho*gh + pRef

Definition at line 80 of file pEqn.H.

U

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

Definition at line 89 of file pEqn.H.