populationBalanceModel: Moved access function for continuous phase turbulence

Patch contributed by Institute of Fluid Dynamics, Helmholtz-Zentrum
Dresden - Rossendorf (HZDR)

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<phaseCompressibleTurbulenceModel>
-        (
-            IOobject::groupName
-            (
-                turbulenceModel::propertiesName,
-                popBal_.continuousPhase().name()
-            )
-        );
-}
-
-
 void Foam::diameterModels::binaryBreakupModel::correct()
 {}
 

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();
 

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)
                 )
             )
         );

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<phaseCompressibleTurbulenceModel>
-        (
-            IOobject::groupName
-            (
-                turbulenceModel::propertiesName,
-                popBal_.continuousPhase().name()
-            )
-        );
-}
-
-
 void Foam::diameterModels::breakupModel::correct()
 {}
 

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<daughterSizeDistributionModel>& dsdPtr() const
         {

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())))

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;

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<phaseCompressibleTurbulenceModel>
-        (
-            IOobject::groupName
-            (
-                turbulenceModel::propertiesName,
-                popBal_.continuousPhase().name()
-            )
-        );
-}
-
-
 void Foam::diameterModels::coalescenceModel::correct()
 {}
 

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();
 

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<phaseCompressibleTurbulenceModel>
+        (
+            IOobject::groupName
+            (
+                turbulenceModel::propertiesName,
+                continuousPhase_.name()
+            )
+        );
+}
+
+
 void Foam::diameterModels::populationBalanceModel::solve()
 {
     const dictionary& solutionControls = mesh_.solverDict(name_);

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();
 };