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

Merge branch 'master' of github.com-OpenFOAM:OpenFOAM/OpenFOAM-dev

Browse Source
This commit is contained in:
Henry Weller
2020-07-03 15:45:18 +01:00
parent 2f30c1364b 3f5df8fb78
commit 3cc2ded6d0
1 changed files with 12 additions and 4 deletions
Download Patch File Download Diff File Expand all files Collapse all files

16
applications/solvers/multiphase/reactingEulerFoam/phaseCompressibleMomentumTransportModels/kineticTheoryModels/kineticTheoryModel/kineticTheoryModel.C
View File

@ -391,6 +391,7 @@ void Foam::RASModels::kineticTheoryModel::correct()
{
// Local references
volScalarField alpha(max(alpha_, scalar(0)));
const phaseSystem& fluid = phase_.fluid();
const phaseModel& continuousPhase = this->continuousPhase();
const volScalarField& rho = phase_.rho();
const surfaceScalarField& alphaRhoPhi = alphaRhoPhi_;
@ -437,10 +438,17 @@ void Foam::RASModels::kineticTheoryModel::correct()
);
// Drag
const dragModel& drag =
phase_.fluid().lookupSubModel<dragModel>(phase_, continuousPhase);
volScalarField beta(drag.K());
volScalarField beta
(
fluid.foundSubModel<dragModel>(phase_, continuousPhase)
? fluid.lookupSubModel<dragModel>(phase_, continuousPhase).K()
: volScalarField::New
(
"beta",
phase_.mesh(),
dimensionedScalar(dragModel::dimK, 0)
)
);
// Eq. 3.25, p. 50 Js = J1 - J2
volScalarField J1("J1", 3*beta);