pEqn.H File Reference

Include dependency graph for pEqn.H:

Go to the source code of this file.

Functions
volVectorField	HbyA (constrainHbyA(rAU() *UEqn.H(), U, p_rgh))
surfaceScalarField	phiHbyA ("phiHbyA", fvc::flux(HbyA)+MRF.zeroFilter(fvc::interpolate(rho *rAU()) *fvc::ddtCorr(U, phi, Uf)))
MRF	makeRelative (phiHbyA)
surfaceScalarField	phig ((mixture.surfaceTensionForce() - ghf *fvc::snGrad(rho)) *rAUf *mesh.magSf())
	constrainPressure (p_rgh, U, phiHbyA, rAUf, MRF)
fvScalarMatrix	Sp_rgh (phaseChange.Sp_rgh(rho, gh, p_rgh))
	if (pimple.transonic())
surfaceScalarField	rho2f (fvc::interpolate(rho2))
	if (mesh.moving())
p_rghEqnComp1	ref () *
volScalarField	p_rgh_0 (p_rgh)
	while (pimple.correctNonOrthogonal())
mixture	thermo1 ().correctRho(psi1 *(p_rgh - p_rgh_0))
mixture	thermo2 ().correctRho(psi2 *(p_rgh - p_rgh_0))
mixture	correct ()
p_rgh	correctBoundaryConditions ()

Variables
	phiHbyA = phig
tmp< fvScalarMatrix >	p_rghEqnComp1
tmp< fvScalarMatrix >	p_rghEqnComp2
	else
	p_rgh = p - rho*gh
	K = 0.5*magSqr(U)

Function Documentation

HbyA()

volVectorField HbyA

(

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

)

phiHbyA()

surfaceScalarField phiHbyA	(	"phiHbyA"	,
		fvc::flux(HbyA)+MRF.zeroFilter(fvc::interpolate(rho *rAU()) *fvc::ddtCorr(U, phi, 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
	)

Sp_rgh()

fvScalarMatrix Sp_rgh

(

phaseChange.

Sp_rghrho, gh, p_rgh

)

Referenced by while().

Here is the caller graph for this function:

if() [1/2]

if

(

pimple.

transonic()

)

Definition at line 44 of file pEqn.H.

References alpha1(), alpha2, alphaPhi2(), Foam::correction(), Foam::fvc::ddt(), Foam::fvc::div(), Foam::fvc::interpolate(), p_rgh, phi, psi1, psi2, rho1, rho1f(), rho2, rho2f(), and Sp.

Here is the call graph for this function:

rho2f()

surfaceScalarField rho2f

(

fvc::interpolate(rho2)

)

Referenced by if().

Here is the caller graph for this function:

if() [2/2]

if

(

mesh.

moving()

)

Definition at line 94 of file pEqn.H.

References alpha1(), alpha2, Foam::fvc::div(), fvModels, mesh, p_rgh, p_rgh_0(), pimple, Foam::pos(), rho1, and rho2.

Here is the call graph for this function:

ref()

p_rghEqnComp1 ref

(

)

p_rgh_0()

volScalarField p_rgh_0

(

p_rgh

)

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

Here is the caller graph for this function:

while()

while

(

pimple.

correctNonOrthogonal()

)

Definition at line 117 of file pEqn.H.

References Foam::fvc::absolute(), alpha1(), alpha2, fvConstraints::constrain(), Foam::fvc::correctUf(), dgdt(), Foam::fvc::div(), fvConstraints, gh, HbyA, Foam::fvc::laplacian(), Foam::max(), p, p_rgh, p_rghEqnComp1, p_rghEqnComp2, phi, phig(), phiHbyA, pimple, pMin(), psi1, psi2, rAU, rAUf(), Foam::fvc::reconstruct(), rho1, rho2, solve(), Sp_rgh(), U, and Uf.

Here is the call graph for this function:

thermo1()

mixture thermo1

(

)

thermo2()

mixture thermo2

(

)

correct()

mixture correct

(

)

correctBoundaryConditions()

p_rgh correctBoundaryConditions

(

)

Variable Documentation

phiHbyA

phiHbyA = phig

Definition at line 30 of file pEqn.H.

p_rghEqnComp1

p_rghEqnComp1

Initial value:

=
            (
                (fvc::ddt(alpha1, rho1) + fvc::div(alphaPhi1*rho1f))/rho1
              - fvc::ddt(alpha1) - fvc::div(alphaPhi1)
              + (alpha1*psi1/rho1)*correction(fvm::ddt(p_rgh))
            )

Definition at line 41 of file pEqn.H.

Referenced by while().

p_rghEqnComp2

p_rghEqnComp2

Initial value:

=
            (
               (fvc::ddt(alpha2, rho2) + fvc::div(alphaPhi2*rho2f))/rho2
             - fvc::ddt(alpha2) - fvc::div(alphaPhi2)
             + (alpha2*psi2/rho2)*correction(fvm::ddt(p_rgh))
            )

Definition at line 42 of file pEqn.H.

Referenced by while().

else

else

Initial value:

{
# 1 "/home/ubuntu/OpenFOAM-9/applications/solvers/multiphase/compressibleInterFoam/rhofs.H" 1
surfaceScalarField rho1f(fvc::interpolate(rho1))

Definition at line 76 of file pEqn.H.

p_rgh

p_rgh = p - rho*gh

Definition at line 159 of file pEqn.H.

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

K

K = 0.5*magSqr(U)

Definition at line 162 of file pEqn.H.