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, phi, 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())
 
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))
 
p_rgh correctBoundaryConditions ()
 

Variables

 phiHbyA = phig
 
tmp< fvScalarMatrix > p_rghEqnComp1
 
tmp< fvScalarMatrix > p_rghEqnComp2
 
 else
 
 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)+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   
)

◆ if() [1/2]

if ( pimple.  transonic())

Definition at line 32 of file pEqn.H.

References alpha1(), alpha2, alphaPhi1, alphaPhi2(), Foam::correction(), Foam::fvc::ddt(), Foam::fvc::div(), fvOptions, Foam::fvc::interpolate(), mixture, p_rgh, phi, 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 96 of file pEqn.H.

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

Here is the call graph for this function:

◆ ref()

p_rghEqnComp1 ref ( )

◆ p_rgh_0()

volScalarField p_rgh_0 ( p_rgh  )

◆ while()

while ( pimple.  correctNonOrthogonal())

◆ 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

p_rghEqnComp1
Initial value:
=
*(
(
- (fvOptions(alpha1, mixture.thermo1().rho())&rho1)
)/rho1
)
tmp< fvMatrix< Type > > correction(const fvMatrix< Type > &)
Return the correction form of the given matrix.
fv::options & fvOptions
surfaceScalarField rho1f(fvc::interpolate(rho1))
tmp< GeometricField< Type, fvPatchField, volMesh > > div(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
Definition: fvcDiv.C:47
p_rgh
Definition: pEqn.H:140
const surfaceScalarField & alphaPhi1
tmp< GeometricField< Type, fvPatchField, volMesh > > ddt(const dimensioned< Type > dt, const fvMesh &mesh)
Definition: fvcDdt.C:45
dimensionedScalar pos(const dimensionedScalar &ds)
volScalarField & alpha1(mixture.alpha1())
phaseChangeTwoPhaseMixture & mixture
Definition: createFields.H:38
const dimensionedScalar & rho1
Definition: createFields.H:43

Definition at line 29 of file pEqn.H.

◆ p_rghEqnComp2

p_rghEqnComp2
Initial value:
=
*(
(
- (fvOptions(alpha2, mixture.thermo2().rho())&rho2)
)/rho2
)
tmp< fvMatrix< Type > > correction(const fvMatrix< Type > &)
Return the correction form of the given matrix.
fv::options & fvOptions
tmp< GeometricField< Type, fvPatchField, volMesh > > div(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
Definition: fvcDiv.C:47
p_rgh
Definition: pEqn.H:140
alpha2
Definition: alphaEqn.H:115
tmp< GeometricField< Type, fvPatchField, volMesh > > ddt(const dimensioned< Type > dt, const fvMesh &mesh)
Definition: fvcDdt.C:45
dimensionedScalar pos(const dimensionedScalar &ds)
surfaceScalarField rho2f(fvc::interpolate(rho2))
surfaceScalarField alphaPhi2("alphaPhi2", phi - alphaPhi1)
const dimensionedScalar & rho2
Definition: createFields.H:44
phaseChangeTwoPhaseMixture & mixture
Definition: createFields.H:38

Definition at line 30 of file pEqn.H.

◆ else

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

Definition at line 70 of file pEqn.H.

◆ rho

Definition at line 154 of file pEqn.H.

◆ p_rgh

p_rgh = p - rho*gh

Definition at line 157 of file pEqn.H.

◆ K

K = 0.5*magSqr(U)

Definition at line 160 of file pEqn.H.