zhyang-dev/CFDEMcoupling-PFM
1
0
Fork 0
Code Issues Pull Requests Actions Packages Projects Releases Wiki Activity

pressure eqn and bc updates

Browse Source
This commit is contained in:
s126103
2019-07-23 09:13:56 +02:00
parent c28480a802
commit 39f6e7d056
2 changed files with 6 additions and 6 deletions
Download Patch File Download Diff File Expand all files Collapse all files

4
applications/solvers/cfdemSolverPiso/pEqn.H
View File

@ -21,7 +21,7 @@ phi = phiHbyA + rAUf*(fvc::interpolate(Ksl/rho) * phiS);
MRF.makeRelative(phi);
// adjustment of phi (only in cases w.o. p boundary conditions) not tested yet
adjustPhi(phi, U, p);
//adjustPhi(phi, U, p);
// Update the pressure BCs to ensure flux consistency
constrainPressure(p, Uvoidfraction, phiHbyA, rAUvoidfraction, MRF);
@ -55,4 +55,4 @@ else
U = HbyA - voidfraction*rAU*fvc::grad(p) + Ksl/rho*Us*rAU;
U.correctBoundaryConditions();
fvOptions.correct(U);
fvOptions.correct(U);

8
src/finiteVolume/fields/fvPatchFields/derived/uniformFixedValueTube/uniformFixedValueTubeFvPatchField.C
View File

@ -89,7 +89,7 @@ uniformFixedValueTubeFvPatchField<Type>::uniformFixedValueTubeFvPatchField
:
fixedValueFvPatchField<Type>(p, iF),
uniformValue_(Function1<Type>::New("uniformValue", dict)),
pName_("p"), //JOKER
pName_("p_rgh"), //JOKER
phiName_("phi"), //JOKER
velocityFieldName_("U"),
densityFieldName_("rho"),
@ -110,7 +110,7 @@ uniformFixedValueTubeFvPatchField<Type>::uniformFixedValueTubeFvPatchField
:
fixedValueFvPatchField<Type>(ptf),
uniformValue_(ptf.uniformValue_().clone().ptr()),
pName_("p"), //JOKER
pName_("p_rgh"), //JOKER
phiName_("phi"), //JOKER
velocityFieldName_("U"),
densityFieldName_("rho"),
@ -132,7 +132,7 @@ uniformFixedValueTubeFvPatchField<Type>::uniformFixedValueTubeFvPatchField
:
fixedValueFvPatchField<Type>(ptf, iF),
uniformValue_(ptf.uniformValue_().clone().ptr()),
pName_("p"), //JOKER
pName_("p_rgh"), //JOKER
phiName_("phi"), //JOKER
velocityFieldName_("U"),
densityFieldName_("rho"),
@ -181,7 +181,7 @@ void uniformFixedValueTubeFvPatchField<Type>::updateCoeffs()
//calc cell velocity from flux phip
//scalar vel = phip/this->patch().magSf();
// some relaxation might be useful?
// calc pressure drop
// fvPatchField<Type>::operator==(0.5*zeta*density*mag(velocity)*mag(velocity)*uniformValue_->value(t));
//dP = zeta * l/d * rho u^2 / 2