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)+MRF.zeroFilter(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 = max(p_rgh + rho*gh, pMin)
 
 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)+MRF.zeroFilter(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 alpha1, alpha2, alphaPhi1, alphaPhi2, fvOptions, mixture, p_rgh, p_rghEqnComp1, p_rghEqnComp2, phi, psi1, psi2, rho1, rho1f(), rho2, and rho2f().

Here is the call graph for this function:

◆ p_rgh_0()

volScalarField p_rgh_0 ( p_rgh  )

Referenced by while().

Here is the caller graph for this function:

◆ while()

while ( pimple.  correctNonOrthogonal())

Definition at line 87 of file pEqn.H.

References alpha1, alpha2, fvOptions, gh, HbyA, max(), mesh, p, p_rgh, p_rghEqnComp1, p_rghEqnComp2, phi, phig(), phiHbyA, pimple, pMin, rAU, rAUf(), rho1, rho2, solve(), and U.

Here is the call graph for this function:

◆ n()

surfaceVectorField n ( mesh.  Sf)/mesh.magSf()

◆ thermo1()

mixture thermo1 ( )

Definition at line 21 of file setRegionFluidFields.H.

◆ thermo2()

mixture thermo2 ( )

Definition at line 22 of file setRegionFluidFields.H.

◆ 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 if(), and while().

◆ p_rghEqnComp2

p_rghEqnComp2
Initial value:
=
fvc::ddt(rho2) + psi2*correction(fvm::ddt(p_rgh))
+ fvc::div(phi, rho2) - fvc::Sp(fvc::div(phi), rho2)
volScalarField & p_rgh
surfaceScalarField & phi
const volScalarField & psi2
volScalarField & rho2

Definition at line 30 of file pEqn.H.

Referenced by if(), and while().

◆ else

else
Initial value:
{
fvc::ddt(rho1) + psi1*correction(fvm::ddt(p_rgh))
+ fvc::div(phi, rho1) - fvc::Sp(fvc::div(phi), rho1)
const volScalarField & psi1
volScalarField & rho1
tmp< fvScalarMatrix > p_rghEqnComp1
Definition: pEqn.H:29

Definition at line 73 of file pEqn.H.

◆ Uf

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

Definition at line 128 of file pEqn.H.

◆ rho

Definition at line 135 of file pEqn.H.

◆ p

p = max(p_rgh + rho*gh, pMin)

Definition at line 138 of file pEqn.H.

◆ p_rgh

p_rgh = p - rho*gh

Definition at line 139 of file pEqn.H.

◆ K

K = 0.5*magSqr(U)

Definition at line 144 of file pEqn.H.