From 9a71022d93b148aebe863a7f1a6dafbf5cb88a61 Mon Sep 17 00:00:00 2001 From: Will Bainbridge Date: Wed, 4 Mar 2020 14:40:02 +0000 Subject: [PATCH] reacting*EulerFoam/.../populationBalanceModel: Corrected mass-transfer source terms Patch contributed by Institute of Fluid Dynamics, Helmholtz-Zentrum Dresden - Rossendorf (HZDR) --- .../sizeGroup/shapeModels/fractal/fractal.C | 11 ++++++++++- .../populationBalanceModel/populationBalanceModel.C | 7 ++++++- 2 files changed, 16 insertions(+), 2 deletions(-) diff --git a/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/diameterModels/velocityGroup/sizeGroup/shapeModels/fractal/fractal.C b/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/diameterModels/velocityGroup/sizeGroup/shapeModels/fractal/fractal.C index eaa5ef2301..62fadc35fd 100644 --- a/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/diameterModels/velocityGroup/sizeGroup/shapeModels/fractal/fractal.C +++ b/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/diameterModels/velocityGroup/sizeGroup/shapeModels/fractal/fractal.C @@ -27,6 +27,7 @@ License #include "sinteringModel.H" #include "fvmDdt.H" #include "fvmDiv.H" +#include "fvcDiv.H" #include "fvmSup.H" #include "fvcSup.H" #include "fvcDdt.H" @@ -214,7 +215,15 @@ void Foam::diameterModels::shapeModels::fractal::correct() fvc::ddt(alpha, rho, fi)*kappa_.oldTime() + alpha*rho*fi*fvm::ddt(kappa_) + fvm::div(fAlphaRhoPhi, kappa_) - + fvm::SuSp(-phase.continuityError()*fi, kappa_) + + fvm::SuSp + ( + fi + *( + fi.VelocityGroup().dmdt() + - (fvc::ddt(alpha, rho) + fvc::div(phase.alphaRhoPhi())) + ), + kappa_ + ) == - sinteringModel_->R() + fvc::Su(Su_*rho, kappa_) diff --git a/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/populationBalanceModel/populationBalanceModel/populationBalanceModel.C b/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/populationBalanceModel/populationBalanceModel/populationBalanceModel.C index dd176dafaa..c64d72eb54 100644 --- a/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/populationBalanceModel/populationBalanceModel/populationBalanceModel.C +++ b/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/populationBalanceModel/populationBalanceModel/populationBalanceModel.C @@ -1304,7 +1304,12 @@ void Foam::diameterModels::populationBalanceModel::solve() phase.alphaRhoPhi(), fi ) - + fvm::SuSp(-phase.continuityError(), fi) + + fvm::SuSp + ( + fi.VelocityGroup().dmdt() + - (fvc::ddt(alpha, rho) + fvc::div(phase.alphaRhoPhi())), + fi + ) == fvc::Su(Su_[i]*rho, fi) - fvm::SuSp(SuSp_[i]*rho, fi)