pEqn.H
Go to the documentation of this file.
1 {
2  volScalarField rAU("rAU", 1.0/UEqn.A());
6  (
7  "phiHbyA",
9  + MRF.zeroFilter(fvc::interpolate(rho*rAU)*fvc::ddtCorr(U, phi))
10  );
11  MRF.makeRelative(phiHbyA);
12 
14  (
15  (
16  mixture.surfaceTensionForce()
18  )*rAUf*mesh.magSf()
19  );
20 
22 
23  // Update the pressure BCs to ensure flux consistency
25 
26  tmp<fvScalarMatrix> p_rghEqnComp1;
27  tmp<fvScalarMatrix> p_rghEqnComp2;
28 
29  if (pimple.transonic())
30  {
31  #include "rhofs.H"
32 
34  surfaceScalarField phid2("phid2", fvc::interpolate(psi2)*phi);
35 
36  p_rghEqnComp1 =
37  pos(alpha1)
38  *(
39  (
41  - (fvOptions(alpha1, mixture.thermo1().rho())&rho1)
42  )/rho1
44  + (alpha1/rho1)
45  *correction
46  (
48  + fvm::div(phid1, p_rgh) - fvm::Sp(fvc::div(phid1), p_rgh)
49  )
50  );
51  p_rghEqnComp1.ref().relax();
52 
53  p_rghEqnComp2 =
54  pos(alpha2)
55  *(
56  (
58  - (fvOptions(alpha2, mixture.thermo2().rho())&rho2)
59  )/rho2
61  + (alpha2/rho2)
62  *correction
63  (
65  + fvm::div(phid2, p_rgh) - fvm::Sp(fvc::div(phid2), p_rgh)
66  )
67  );
68  p_rghEqnComp2.ref().relax();
69  }
70  else
71  {
72  #include "rhofs.H"
73 
74  p_rghEqnComp1 =
75  pos(alpha1)
76  *(
77  (
79  - (fvOptions(alpha1, mixture.thermo1().rho())&rho1)
80  )/rho1
83  ) - surfaceFilm.Srho()/rho1;
84 
85  p_rghEqnComp2 =
86  pos(alpha2)
87  *(
88  (
90  - (fvOptions(alpha2, mixture.thermo2().rho())&rho2)
91  )/rho2
94  );
95  }
96 
97  // Cache p_rgh prior to solve for density update
99 
100  while (pimple.correctNonOrthogonal())
101  {
102  fvScalarMatrix p_rghEqnIncomp
103  (
106  );
107 
108  solve
109  (
110  p_rghEqnComp1() + p_rghEqnComp2() + p_rghEqnIncomp
111  );
112 
113  if (pimple.finalNonOrthogonalIter())
114  {
115  p = max(p_rgh + (alpha1*rho1 + alpha2*rho2)*gh, pMin);
116  p_rgh = p - (alpha1*rho1 + alpha2*rho2)*gh;
117 
118  dgdt =
119  (
120  alpha1*(p_rghEqnComp2 & p_rgh)
121  - alpha2*(p_rghEqnComp1 & p_rgh)
122  );
123 
124  phi = phiHbyA + p_rghEqnIncomp.flux();
125 
126  U = HbyA
127  + rAU*fvc::reconstruct((phig + p_rghEqnIncomp.flux())/rAUf);
128  U.correctBoundaryConditions();
129  fvOptions.correct(U);
130  }
131  }
132 
133  // Update densities from change in p_rgh
134  mixture.thermo1().correctRho(psi1*(p_rgh - p_rgh_0));
135  mixture.thermo2().correctRho(psi2*(p_rgh - p_rgh_0));
136 
138 
139  // Correct p_rgh for consistency with p and the updated densities
140  p_rgh = p - rho*gh;
141  p_rgh.correctBoundaryConditions();
142 
143  K = 0.5*magSqr(U);
144 }
tmp< fvMatrix< Type > > correction(const fvMatrix< Type > &)
Return the correction form of the given matrix.
rAU
Definition: pEqn.H:1
phiHbyA
Definition: pEqn.H:22
psi2
Definition: TEqns.H:35
fvMatrix< scalar > fvScalarMatrix
Definition: fvMatricesFwd.H:42
const surfaceScalarField & alphaPhi2
fv::options & fvOptions
rho
Definition: pEqn.H:1
pimpleNoLoopControl & pimple
p
Definition: pEqn.H:50
IOMRFZoneList & MRF
surfaceScalarField rho1f(fvc::interpolate(rho1))
dimensioned< Type > max(const dimensioned< Type > &, const dimensioned< Type > &)
surfaceScalarField rAUf("rAUf", fvc::interpolate(rAU))
volScalarField p_rgh_0(p_rgh)
tmp< GeometricField< Type, fvPatchField, volMesh > > div(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
Definition: fvcDiv.C:47
surfaceScalarField phig("phig", -rhorAUf *ghf *fvc::snGrad(rho) *mesh.magSf())
tmp< GeometricField< typename flux< Type >::type, fvsPatchField, surfaceMesh > > ddtCorr(const GeometricField< Type, fvPatchField, volMesh > &U, const GeometricField< Type, fvsPatchField, surfaceMesh > &Uf)
Definition: fvcDdt.C:170
p_rgh
Definition: pEqn.H:140
const surfaceScalarField & alphaPhi1
GeometricField< vector, fvPatchField, volMesh > volVectorField
Definition: volFieldsFwd.H:55
tmp< GeometricField< Type, fvPatchField, volMesh > > laplacian(const GeometricField< Type, fvPatchField, volMesh > &vf, const word &name)
Definition: fvcLaplacian.C:45
CGAL::Exact_predicates_exact_constructions_kernel K
tmp< GeometricField< Type, fvPatchField, volMesh > > ddt(const dimensioned< Type > dt, const fvMesh &mesh)
Definition: fvcDdt.C:45
GeometricField< scalar, fvPatchField, volMesh > volScalarField
Definition: volFieldsFwd.H:52
dimensionedScalar pos(const dimensionedScalar &ds)
surfaceScalarField rho2f(fvc::interpolate(rho2))
dynamicFvMesh & mesh
const volScalarField & alpha1
rhoEqn solve()
tmp< volVectorField > constrainHbyA(const tmp< volVectorField > &tHbyA, const volVectorField &U, const volScalarField &p)
Definition: constrainHbyA.C:33
alpha2
Definition: alphaEqn.H:115
tmp< fvScalarMatrix > p_rghEqnComp2
Definition: pEqn.H:27
dimensioned< scalar > magSqr(const dimensioned< Type > &)
regionModels::surfaceFilmModel & surfaceFilm
const surfaceScalarField & ghf
dimensionedScalar pMin("pMin", dimPressure, fluid)
static tmp< GeometricField< Type, fvsPatchField, surfaceMesh > > interpolate(const GeometricField< Type, fvPatchField, volMesh > &tvf, const surfaceScalarField &faceFlux, Istream &schemeData)
Interpolate field onto faces using scheme given by Istream.
tmp< GeometricField< typename outerProduct< vector, Type >::type, fvPatchField, volMesh >> reconstruct(const GeometricField< Type, fvsPatchField, surfaceMesh > &ssf)
psi1
Definition: TEqns.H:34
U
Definition: pEqn.H:72
volVectorField & HbyA
Definition: pEqn.H:13
fvVectorMatrix & UEqn
Definition: UEqn.H:13
const volScalarField & gh
phi
Definition: pEqn.H:18
constrainPressure(p_rgh, rho, U, phiHbyA, rhorAUf, MRF)
const dimensionedScalar & rho2
Definition: createFields.H:44
tmp< surfaceScalarField > flux(const volVectorField &vvf)
Return the face-flux field obtained from the given volVectorField.
Definition: fvcFlux.C:32
GeometricField< scalar, fvsPatchField, surfaceMesh > surfaceScalarField
phaseChangeTwoPhaseMixture & mixture
Definition: createFields.H:38
tmp< fvScalarMatrix > p_rghEqnComp1
Definition: pEqn.H:26
const dimensionedScalar & rho1
Definition: createFields.H:43
tmp< GeometricField< Type, fvsPatchField, surfaceMesh > > snGrad(const GeometricField< Type, fvPatchField, volMesh > &vf, const word &name)
Definition: fvcSnGrad.C:45
zeroField Sp
Definition: alphaSuSp.H:2