pEqn.H File Reference

Include dependency graph for pEqn.H:

Go to the source code of this file.

Functions
surfaceScalarField	rAUf ("rAUf", fvc::interpolate(rAU))
volVectorField	HbyA (constrainHbyA(rAU *UEqn.H(), U, p_rgh))
surfaceScalarField	phiHbyA ("phiHbyA", fvc::flux(HbyA)+fvc::interpolate(rho *rAU) *fvc::ddtCorr(U, Uf))
MRF	makeRelative (phiHbyA)
surfaceScalarField	phig ((mixture.surfaceTensionForce() - ghf *fvc::snGrad(rho)) *rAUf *mesh.magSf())
	constrainPressure (p_rgh, U, phiHbyA, rAUf, MRF)
	if (pimple.transonic())
volScalarField	p_rgh_0 (p_rgh)
	while (pimple.correctNonOrthogonal())
surfaceVectorField	n (mesh.Sf()/mesh.magSf())
mixture	thermo1 ().correctRho(psi1 *(p_rgh - p_rgh_0))
mixture	thermo2 ().correctRho(psi2 *(p_rgh - p_rgh_0))
p_rgh	correctBoundaryConditions ()

Variables
	phiHbyA = phig
tmp< fvScalarMatrix >	p_rghEqnComp1
tmp< fvScalarMatrix >	p_rghEqnComp2
	else
	Uf = n*(fvc::absolute(phi, U)/mesh.magSf() - (n & Uf))
	rho = alpha1*rho1 + alpha2*rho2
	p_rgh = p - rho*gh
	K = 0.5*magSqr(U)

Function Documentation

rAUf()

surfaceScalarField rAUf	(	"rAUf"	,
		fvc::interpolate(rAU)
	)

HbyA()

volVectorField HbyA

(

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

)

phiHbyA()

surfaceScalarField phiHbyA	(	"phiHbyA"	,
		fvc::flux(HbyA)+fvc::interpolate(rho *rAU) *fvc::ddtCorr(U, Uf)
	)

makeRelative()

MRF makeRelative

(

phiHbyA

)

phig()

surfaceScalarField phig

(

(mixture.surfaceTensionForce() - ghf *fvc::snGrad(rho)) *rAUf *mesh.magSf()

)

constrainPressure()

constrainPressure	(	p_rgh	,
		U	,
		phiHbyA	,
		rAUf	,
		MRF
	)

if()

if

(

pimple.

transonic()

)

Definition at line 32 of file pEqn.H.

References Foam::correction(), Foam::fvc::ddt(), Foam::deleteDemandDrivenData(), Foam::fvc::div(), Foam::fvc::interpolate(), p_rgh, phi, psi1, psi2, rho1, rho2, and Sp.

Here is the call graph for this function:

p_rgh_0()

volScalarField p_rgh_0

(

p_rgh

)

Referenced by forAllConstIter(), and while().

Here is the caller graph for this function:

while()

while

(

pimple.

correctNonOrthogonal()

)

Definition at line 87 of file pEqn.H.

References alpha1, alpha2, Foam::fvc::div(), fvOptions, gh, HbyA, Foam::fvc::interpolate(), Foam::fvc::laplacian(), Foam::max(), mesh, n(), p, p_rgh, p_rghEqnComp1, p_rghEqnComp2, phi, phig(), phiHbyA, pimple, pMin(), rAUf(), Foam::fvc::reconstruct(), rho1, rho2, solve(), U, and Uf.

Here is the call graph for this function:

n()

surfaceVectorField n

(

mesh.

Sf)/mesh.magSf(

)

thermo1()

mixture thermo1

(

)

thermo2()

mixture thermo2

(

)

correctBoundaryConditions()

p_rgh correctBoundaryConditions

(

)

Variable Documentation

phiHbyA

phiHbyA = phig

Definition at line 21 of file pEqn.H.

p_rghEqnComp1

tmp<fvScalarMatrix> p_rghEqnComp1

Definition at line 29 of file pEqn.H.

Referenced by while().

p_rghEqnComp2

p_rghEqnComp2

Initial value:

=
            pos(alpha2)
           *(
               (
                   fvc::ddt(alpha2, rho2) + fvc::div(alphaPhi2*rho2f)
                 - (fvOptions(alpha2, mixture.thermo2().rho())&rho2)
               )/rho2
             - fvc::ddt(alpha2) - fvc::div(alphaPhi2)
             + (alpha2*psi2/rho2)*correction(fvm::ddt(p_rgh))
            )

Definition at line 30 of file pEqn.H.

Referenced by while().

else

else

Initial value:

{


        p_rghEqnComp1 =
            pos(alpha1)
           *(
                (
                    fvc::ddt(alpha1, rho1) + fvc::div(alphaPhi1*rho1f)
                  - (fvOptions(alpha1, mixture.thermo1().rho())&rho1)
                )/rho1
              - fvc::ddt(alpha1) - fvc::div(alphaPhi1)
              + (alpha1*psi1/rho1)*correction(fvm::ddt(p_rgh))
            )

Definition at line 58 of file pEqn.H.

Uf

Uf = n*(fvc::absolute(phi, U)/mesh.magSf() - (n & Uf))

Definition at line 124 of file pEqn.H.

rho

rho = alpha1*rho1 + alpha2*rho2

Definition at line 131 of file pEqn.H.

p_rgh

p_rgh = p - rho*gh

Definition at line 134 of file pEqn.H.

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

K

K = 0.5*magSqr(U)

Definition at line 137 of file pEqn.H.