Loading...
Searching...
No Matches
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)))
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))
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)  )

References rAU.

◆ HbyA()

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

References p_rgh, rAU, U, and UEqn.

◆ phiHbyA()

surfaceScalarField phiHbyA ( "phiHbyA" ,
fvc::flux(HbyA)+MRF.zeroFilter(fvc::interpolate(rho *rAU) *fvc::ddtCorr(U, phi))  )

References HbyA, MRF, phi, rAU, rho, and U.

◆ makeRelative()

MRF makeRelative ( phiHbyA )

References MRF, and phiHbyA.

◆ phig()

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

References ghf, mesh, mixture, rAUf(), and rho.

Here is the call graph for this function:

◆ constrainPressure()

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

References MRF, p_rgh, phiHbyA, rAUf(), and U.

Here is the call graph for this function:

◆ if()

if ( pimple. transonic())

Definition at line 29 of file pEqn.H.

References alpha1, alpha2, alphaPhi1, alphaPhi2, fvOptions, mixture, p_rgh, p_rghEqnComp1, p_rghEqnComp2, phi, pimple, psi1, psi2, rho1, rho1f(), rho2, and rho2f().

Here is the call graph for this function:

◆ rho2f()

surfaceScalarField rho2f ( fvc::interpolate(rho2) )

References alpha1, alpha2, alphaPhi1, alphaPhi2, mixture, p_rgh, p_rghEqnComp1, p_rghEqnComp2, psi1, psi2, rho1, rho1f(), rho2, and rho2f().

Referenced by if(), if(), if(), and rho2f().

Here is the call graph for this function:
Here is the caller graph for this function:

◆ p_rgh_0()

volScalarField p_rgh_0 ( p_rgh )

References p_rgh.

◆ while()

while ( pimple. correctNonOrthogonal())

Definition at line 100 of file pEqn.H.

References alpha1, alpha2, fvOptions, gh, HbyA, 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:

◆ thermo1()

mixture thermo1 ( )

Definition at line 21 of file setRegionFluidFields.H.

Referenced by correct(), Cpv1(), E1Eqn(), for(), if(), and validate().

Here is the caller graph for this function:

◆ thermo2()

mixture thermo2 ( )

Definition at line 22 of file setRegionFluidFields.H.

Referenced by correct(), Cpv2(), E2Eqn(), for(), if(), and validate().

Here is the caller graph for this function:

◆ correctBoundaryConditions()

p_rgh correctBoundaryConditions ( )

References p_rgh.

Variable Documentation

◆ phiHbyA

phiHbyA = phig

Definition at line 21 of file pEqn.H.

◆ p_rghEqnComp1

p_rghEqnComp1
Initial value:
=
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))
)
volScalarField & p_rgh
fv::options & fvOptions
const volScalarField & alpha1
const auto & alphaPhi1
const volScalarField & psi1
volScalarField & rho1
surfaceScalarField rho1f(fvc::interpolate(rho1))
Info<< "Creating temperaturePhaseChangeTwoPhaseMixture\n"<< endl;autoPtr< temperaturePhaseChangeTwoPhaseMixture > mixture

Definition at line 26 of file pEqn.H.

◆ 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))
)
const auto & alphaPhi2
const volScalarField & psi2
volScalarField & rho2
const volScalarField & alpha2
surfaceScalarField rho2f(fvc::interpolate(rho2))

Definition at line 27 of file pEqn.H.

◆ else

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

Definition at line 70 of file pEqn.H.

◆ rho

Definition at line 138 of file pEqn.H.

◆ p_rgh

p_rgh = p - rho*gh

Definition at line 141 of file pEqn.H.

◆ K

K = 0.5*magSqr(U)

Definition at line 144 of file pEqn.H.