From f6d24cf29f449a7446806c6c077c92920d009899 Mon Sep 17 00:00:00 2001 From: Will Bainbridge Date: Wed, 16 May 2018 21:36:34 +0100 Subject: [PATCH] populationBalanceModel: Moved access function for continuous phase turbulence Patch contributed by Institute of Fluid Dynamics, Helmholtz-Zentrum Dresden - Rossendorf (HZDR) --- .../binaryBreakupModel/binaryBreakupModel.C | 15 --------------- .../binaryBreakupModel/binaryBreakupModel.H | 4 ---- .../LaakkonenAlopaeusAittamaa.C | 7 ++++--- .../breakupModels/breakupModel/breakupModel.C | 15 --------------- .../breakupModels/breakupModel/breakupModel.H | 4 ---- .../CoulaloglouTavlaridesCoalescence.C | 5 +++-- .../PrinceBlanch/PrinceBlanch.C | 5 +++-- .../coalescenceModel/coalescenceModel.C | 15 --------------- .../coalescenceModel/coalescenceModel.H | 4 ---- .../populationBalanceModel.C | 16 ++++++++++++++++ .../populationBalanceModel.H | 4 ++++ 11 files changed, 30 insertions(+), 64 deletions(-) diff --git a/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/populationBalanceModel/binaryBreakupModels/binaryBreakupModel/binaryBreakupModel.C b/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/populationBalanceModel/binaryBreakupModels/binaryBreakupModel/binaryBreakupModel.C index da326223b..9d1c0060c 100644 --- a/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/populationBalanceModel/binaryBreakupModels/binaryBreakupModel/binaryBreakupModel.C +++ b/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/populationBalanceModel/binaryBreakupModels/binaryBreakupModel/binaryBreakupModel.C @@ -78,21 +78,6 @@ Foam::diameterModels::binaryBreakupModel::binaryBreakupModel // * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * // -const Foam::phaseCompressibleTurbulenceModel& -Foam::diameterModels::binaryBreakupModel::continuousTurbulence() const -{ - return - popBal_.mesh().lookupObject - ( - IOobject::groupName - ( - turbulenceModel::propertiesName, - popBal_.continuousPhase().name() - ) - ); -} - - void Foam::diameterModels::binaryBreakupModel::correct() {} diff --git a/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/populationBalanceModel/binaryBreakupModels/binaryBreakupModel/binaryBreakupModel.H b/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/populationBalanceModel/binaryBreakupModels/binaryBreakupModel/binaryBreakupModel.H index 5ed5385e5..01e4b5b72 100644 --- a/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/populationBalanceModel/binaryBreakupModels/binaryBreakupModel/binaryBreakupModel.H +++ b/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/populationBalanceModel/binaryBreakupModels/binaryBreakupModel/binaryBreakupModel.H @@ -38,7 +38,6 @@ SourceFiles #define binaryBreakupModel_H #include "populationBalanceModel.H" -#include "phaseCompressibleTurbulenceModel.H" // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // @@ -136,9 +135,6 @@ public: // Member Functions - //- Return a reference to the turbulence model for the continuous phase - const phaseCompressibleTurbulenceModel& continuousTurbulence() const; - //- Correct diameter independent expressions virtual void correct(); diff --git a/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/populationBalanceModel/breakupModels/LaakkonenAlopaeusAittamaa/LaakkonenAlopaeusAittamaa.C b/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/populationBalanceModel/breakupModels/LaakkonenAlopaeusAittamaa/LaakkonenAlopaeusAittamaa.C index 088817961..3e77f16fa 100644 --- a/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/populationBalanceModel/breakupModels/LaakkonenAlopaeusAittamaa/LaakkonenAlopaeusAittamaa.C +++ b/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/populationBalanceModel/breakupModels/LaakkonenAlopaeusAittamaa/LaakkonenAlopaeusAittamaa.C @@ -81,7 +81,7 @@ Foam::diameterModels::breakupModels::LaakkonenAlopaeusAittamaa::setBreakupRate const sizeGroup& fi = *popBal_.sizeGroups()[i]; breakupRate = - C1_*cbrt(continuousTurbulence().epsilon()) + C1_*cbrt(popBal_.continuousTurbulence().epsilon()) *erfc ( sqrt @@ -89,12 +89,13 @@ Foam::diameterModels::breakupModels::LaakkonenAlopaeusAittamaa::setBreakupRate C2_*popBal_.sigmaWithContinuousPhase(fi.phase()) /( continuousPhase.rho()*pow(fi.d(), 5.0/3.0) - *pow(continuousTurbulence().epsilon(), 2.0/3.0) + *pow(popBal_.continuousTurbulence().epsilon(), 2.0/3.0) ) + C3_*continuousPhase.mu() /( sqrt(continuousPhase.rho()*fi.phase().rho()) - *cbrt(continuousTurbulence().epsilon())*pow(fi.d(), 4.0/3.0) + *cbrt(popBal_.continuousTurbulence().epsilon()) + *pow(fi.d(), 4.0/3.0) ) ) ); diff --git a/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/populationBalanceModel/breakupModels/breakupModel/breakupModel.C b/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/populationBalanceModel/breakupModels/breakupModel/breakupModel.C index 51dc97b50..af11fec39 100644 --- a/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/populationBalanceModel/breakupModels/breakupModel/breakupModel.C +++ b/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/populationBalanceModel/breakupModels/breakupModel/breakupModel.C @@ -81,21 +81,6 @@ Foam::diameterModels::breakupModel::breakupModel // * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * // -const Foam::phaseCompressibleTurbulenceModel& -Foam::diameterModels::breakupModel::continuousTurbulence() const -{ - return - popBal_.mesh().lookupObject - ( - IOobject::groupName - ( - turbulenceModel::propertiesName, - popBal_.continuousPhase().name() - ) - ); -} - - void Foam::diameterModels::breakupModel::correct() {} diff --git a/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/populationBalanceModel/breakupModels/breakupModel/breakupModel.H b/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/populationBalanceModel/breakupModels/breakupModel/breakupModel.H index 5f41e1cf1..e50a732da 100644 --- a/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/populationBalanceModel/breakupModels/breakupModel/breakupModel.H +++ b/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/populationBalanceModel/breakupModels/breakupModel/breakupModel.H @@ -38,7 +38,6 @@ SourceFiles #include "populationBalanceModel.H" #include "daughterSizeDistributionModel.H" -#include "phaseCompressibleTurbulenceModel.H" // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // @@ -148,9 +147,6 @@ public: return popBal_; } - //- Return a reference to the turbulence model for the continuous phase - const phaseCompressibleTurbulenceModel& continuousTurbulence() const; - //- Return const-reference to daughter size distribution pointer const autoPtr& dsdPtr() const { diff --git a/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/populationBalanceModel/coalescenceModels/CoulaloglouTavlaridesCoalescence/CoulaloglouTavlaridesCoalescence.C b/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/populationBalanceModel/coalescenceModels/CoulaloglouTavlaridesCoalescence/CoulaloglouTavlaridesCoalescence.C index 147827b56..8e56d6c79 100644 --- a/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/populationBalanceModel/coalescenceModels/CoulaloglouTavlaridesCoalescence/CoulaloglouTavlaridesCoalescence.C +++ b/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/populationBalanceModel/coalescenceModels/CoulaloglouTavlaridesCoalescence/CoulaloglouTavlaridesCoalescence.C @@ -25,6 +25,7 @@ License #include "CoulaloglouTavlaridesCoalescence.H" #include "addToRunTimeSelectionTable.H" +#include "phaseCompressibleTurbulenceModel.H" // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * // @@ -79,11 +80,11 @@ addToCoalescenceRate coalescenceRate += C1_*(pow(fi.x(), 2.0/3.0) + pow(fj.x(), 2.0/3.0)) *sqrt(pow(fi.x(), 2.0/9.0) + pow(fj.x(), 2.0/9.0)) - *cbrt(continuousTurbulence().epsilon())/(1 + popBal_.alphas()) + *cbrt(popBal_.continuousTurbulence().epsilon())/(1 + popBal_.alphas()) *exp ( - C2_*continuousPhase.mu()*continuousPhase.rho() - *continuousTurbulence().epsilon() + *popBal_.continuousTurbulence().epsilon() /sqr(popBal_.sigmaWithContinuousPhase(fi.phase())) /pow3(1 + popBal_.alphas()) *pow4(cbrt(fi.x())*cbrt(fj.x())/(cbrt(fi.x()) + cbrt(fj.x()))) diff --git a/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/populationBalanceModel/coalescenceModels/PrinceBlanch/PrinceBlanch.C b/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/populationBalanceModel/coalescenceModels/PrinceBlanch/PrinceBlanch.C index fca600503..030c689de 100644 --- a/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/populationBalanceModel/coalescenceModels/PrinceBlanch/PrinceBlanch.C +++ b/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/populationBalanceModel/coalescenceModels/PrinceBlanch/PrinceBlanch.C @@ -26,6 +26,7 @@ License #include "PrinceBlanch.H" #include "addToRunTimeSelectionTable.H" #include "mathematicalConstants.H" +#include "phaseCompressibleTurbulenceModel.H" // * * * * * * * * * * * * * * Static Data Members * * * * * * * * * * * * * // @@ -95,7 +96,7 @@ addToCoalescenceRate /(16.0*popBal_.sigmaWithContinuousPhase(fi.phase())) ) *log(h0_/hf_) - *cbrt(continuousTurbulence().epsilon())/pow(rij, 2.0/3.0) + *cbrt(popBal_.continuousTurbulence().epsilon())/pow(rij, 2.0/3.0) ); if (turbulentCollisions_) @@ -103,7 +104,7 @@ addToCoalescenceRate coalescenceRate += ( C1_*pi*sqr(fi.d() + fj.d()) - *cbrt(continuousTurbulence().epsilon()) + *cbrt(popBal_.continuousTurbulence().epsilon()) *sqrt(pow(fi.d(), 2.0/3.0) + pow(fj.d(), 2.0/3.0)) ) *collisionEfficiency; diff --git a/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/populationBalanceModel/coalescenceModels/coalescenceModel/coalescenceModel.C b/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/populationBalanceModel/coalescenceModels/coalescenceModel/coalescenceModel.C index 799467ecf..0c9f9d3ab 100644 --- a/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/populationBalanceModel/coalescenceModels/coalescenceModel/coalescenceModel.C +++ b/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/populationBalanceModel/coalescenceModels/coalescenceModel/coalescenceModel.C @@ -78,21 +78,6 @@ Foam::diameterModels::coalescenceModel::coalescenceModel // * * * * * * * * * * * * * * * Member Functions * * * * * * * * * * * * * // -const Foam::phaseCompressibleTurbulenceModel& -Foam::diameterModels::coalescenceModel::continuousTurbulence() const -{ - return - popBal_.mesh().lookupObject - ( - IOobject::groupName - ( - turbulenceModel::propertiesName, - popBal_.continuousPhase().name() - ) - ); -} - - void Foam::diameterModels::coalescenceModel::correct() {} diff --git a/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/populationBalanceModel/coalescenceModels/coalescenceModel/coalescenceModel.H b/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/populationBalanceModel/coalescenceModels/coalescenceModel/coalescenceModel.H index 9c4d30817..d836d819c 100644 --- a/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/populationBalanceModel/coalescenceModels/coalescenceModel/coalescenceModel.H +++ b/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/populationBalanceModel/coalescenceModels/coalescenceModel/coalescenceModel.H @@ -36,7 +36,6 @@ SourceFiles #define coalescenceModel_H #include "populationBalanceModel.H" -#include "phaseCompressibleTurbulenceModel.H" // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // @@ -134,9 +133,6 @@ public: // Member Functions - //- Return a reference to the turbulence model for the continuous phase - const phaseCompressibleTurbulenceModel& continuousTurbulence() const; - //- Correct diameter independent expressions virtual void correct(); diff --git a/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/populationBalanceModel/populationBalanceModel/populationBalanceModel.C b/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/populationBalanceModel/populationBalanceModel/populationBalanceModel.C index f5a889a86..6d269cd03 100644 --- a/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/populationBalanceModel/populationBalanceModel/populationBalanceModel.C +++ b/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/populationBalanceModel/populationBalanceModel/populationBalanceModel.C @@ -34,6 +34,7 @@ License #include "fvcDdt.H" #include "fvmSup.H" #include "fvcDiv.H" +#include "phaseCompressibleTurbulenceModel.H" // * * * * * * * * * * * * Private Member Functions * * * * * * * * * * * * // @@ -1161,6 +1162,21 @@ Foam::diameterModels::populationBalanceModel::sigmaWithContinuousPhase } +const Foam::phaseCompressibleTurbulenceModel& +Foam::diameterModels::populationBalanceModel::continuousTurbulence() const +{ + return + mesh_.lookupObject + ( + IOobject::groupName + ( + turbulenceModel::propertiesName, + continuousPhase_.name() + ) + ); +} + + void Foam::diameterModels::populationBalanceModel::solve() { const dictionary& solutionControls = mesh_.solverDict(name_); diff --git a/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/populationBalanceModel/populationBalanceModel/populationBalanceModel.H b/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/populationBalanceModel/populationBalanceModel/populationBalanceModel.H index 679361071..98886c0be 100644 --- a/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/populationBalanceModel/populationBalanceModel/populationBalanceModel.H +++ b/applications/solvers/multiphase/reactingEulerFoam/phaseSystems/populationBalanceModel/populationBalanceModel/populationBalanceModel.H @@ -154,6 +154,7 @@ SourceFiles #include "sizeGroup.H" #include "phasePair.H" #include "pimpleControl.H" +#include "phaseCompressibleTurbulenceModelFwd.H" // * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * * // @@ -419,6 +420,9 @@ public: const phaseModel& dispersedPhase ) const; + //- Return reference to turbulence model of the continuous phase + const phaseCompressibleTurbulenceModel& continuousTurbulence() const; + //- Solve the population balance equation void solve(); };