[rans] Move Prandtl and Schmidt numbers to problem on RANS level

dumux/freeflow/rans/problem.hh

@@ -371,6 +371,28 @@ public:
     const Scalar karmanConstant() const
     { return 0.41; }
 
+    /*!
+     * \brief Return the turbulent Prandtl number \f$ [-] \f$ which is used to convert
+     *        the eddy viscosity to an eddy thermal conductivity
+     */
+    Scalar turbulentPrandtlNumber() const
+    {
+        static const Scalar turbulentPrandtlNumber
+            = getParamFromGroup<Scalar>(this->paramGroup(), "RANS.TurbulentPrandtlNumber", 1.0);
+        return turbulentPrandtlNumber;
+    }
+
+    /*!
+     * \brief Return the turbulent Schmidt number \f$ [-] \f$ which is used to convert
+     *        the eddy viscosity to an eddy diffusivity
+     */
+    Scalar turbulentSchmidtNumber() const
+    {
+        static const Scalar turbulentSchmidtNumber
+            = getParamFromGroup<Scalar>(this->paramGroup(), "RANS.TurbulentSchmidtNumber", 1.0);
+        return turbulentSchmidtNumber;
+    }
+
 public:
     std::vector<unsigned int> wallElementID_;
     std::vector<Scalar> wallDistance_;

dumux/freeflow/rans/twoeq/kepsilon/volumevariables.hh

@@ -155,10 +155,8 @@ public:
     template<class Problem>
     void calculateEddyDiffusivity(const Problem& problem)
     {
-        static const auto turbulentSchmidtNumber
-            = getParamFromGroup<Scalar>(problem.paramGroup(),
-                                        "RANS.TurbulentSchmidtNumber", 1.0);
-        eddyDiffusivity_ = RANSParentType::kinematicEddyViscosity() / turbulentSchmidtNumber;
+        eddyDiffusivity_ = RANSParentType::kinematicEddyViscosity()
+                           / problem.turbulentSchmidtNumber();
     }
 
     /*!

dumux/freeflow/rans/twoeq/komega/volumevariables.hh

@@ -161,10 +161,8 @@ public:
     template<class Problem>
     void calculateEddyDiffusivity(const Problem& problem)
     {
-        static const auto turbulentSchmidtNumber
-            = getParamFromGroup<Scalar>(problem.paramGroup(),
-                                        "RANS.TurbulentSchmidtNumber", 1.0);
-        eddyDiffusivity_ = RANSParentType::kinematicEddyViscosity() / turbulentSchmidtNumber;
+        eddyDiffusivity_ = RANSParentType::kinematicEddyViscosity()
+                           / problem.turbulentSchmidtNumber();
     }
 
     /*!

dumux/freeflow/rans/twoeq/lowrekepsilon/volumevariables.hh

@@ -144,10 +144,8 @@ public:
     template<class Problem>
     void calculateEddyDiffusivity(const Problem& problem)
     {
-        static const auto turbulentSchmidtNumber
-            = getParamFromGroup<Scalar>(problem.paramGroup(),
-                                        "RANS.TurbulentSchmidtNumber", 1.0);
-        eddyDiffusivity_ = RANSParentType::kinematicEddyViscosity() / turbulentSchmidtNumber;
+        eddyDiffusivity_ = RANSParentType::kinematicEddyViscosity()
+                           / problem.turbulentSchmidtNumber();
     }
 
     /*!

dumux/freeflow/rans/volumevariables.hh

@@ -215,14 +215,10 @@ public:
     template<class Problem>
     void calculateEddyThermalConductivity(const Problem& problem)
     {
-        static const auto turbulentPrandtlNumber
-            = getParamFromGroup<Scalar>(problem.paramGroup(),
-                                        "RANS.TurbulentPrandtlNumber", 1.0);
-
         eddyThermalConductivity_ = ParentType::kinematicEddyViscosity()
                                    * ParentType::asImp_().density()
                                    * ParentType::asImp_().heatCapacity()
-                                   / turbulentPrandtlNumber;
+                                   / problem.turbulentPrandtlNumber();
     }
 
     /*!

dumux/freeflow/rans/zeroeq/volumevariables.hh

@@ -183,10 +183,8 @@ public:
     template<class Problem>
     void calculateEddyDiffusivity(const Problem& problem)
     {
-        static const auto turbulentSchmidtNumber
-            = getParamFromGroup<Scalar>(problem.paramGroup(),
-                                        "RANS.TurbulentSchmidtNumber", 1.0);
-        eddyDiffusivity_ = RANSParentType::kinematicEddyViscosity() / turbulentSchmidtNumber;
+        eddyDiffusivity_ = RANSParentType::kinematicEddyViscosity()
+                           / problem.turbulentSchmidtNumber();
     }
 
     /*!