fixup after rebase

OPM · Jan 15, 2025 · de11bb3 · de11bb3
1 parent 6284483
commit de11bb3
Show file tree

Hide file tree

Showing 4 changed files with 11 additions and 17 deletions.
diff --git a/opm/models/blackoil/blackoilenergymodules.hh b/opm/models/blackoil/blackoilenergymodules.hh
@@ -477,19 +477,11 @@ public:
         }
     }
 
-    void updateTemperature_([[maybe_unused]] const ElementContext& elemCtx,
-                            [[maybe_unused]] unsigned dofIdx,
-                            [[maybe_unused]] unsigned timeIdx)
+    void updateTemperature_(const ElementContext& elemCtx,
+                            unsigned dofIdx,
+                            unsigned timeIdx)
     {
-        auto& fs = asImp_().fluidState_;
-        const Scalar T = elemCtx.problem().temperature(elemCtx, dofIdx, timeIdx);
-        // if TEMP is used we solve for temperature seperatly
-        if constexpr (enableTemperature) {
-            const Evaluation TE = Evaluation::createVariable(T, Indices::temperatureIdx);
-            fs.setTemperature(TE);
-        } else {
-            fs.setTemperature(T);
-        }
+        updateTemperature_(elemCtx.problem(), elemCtx.globalSpaceIndex(dofIdx, timeIdx), timeIdx);
     }
 
     template<class Problem>
@@ -500,7 +492,7 @@ public:
                             [[maybe_unused]] const LinearizationType& lintype
         )
     {
-        updateTemperature_(problem,globalDofIdx, timeIdx);
+        updateTemperature_(problem, globalDofIdx, timeIdx);
     }
 
     void updateEnergyQuantities_(const Problem& problem, unsigned globalDofIdx, unsigned timeIdx)

diff --git a/opm/models/blackoil/blackoilintensivequantities.hh b/opm/models/blackoil/blackoilintensivequantities.hh
@@ -175,9 +175,9 @@ public:
 
     void updateTempSalt(const ElementContext& elemCtx, unsigned dofIdx, unsigned timeIdx)
     {
-        if constexpr (enableTemperature || enableEnergy) {
-            asImp_().updateTemperature_(elemCtx, dofIdx, timeIdx);
-        }
+        // the temperature is updated even if the energy equations are not solved
+        // to allow for temperature dependent properites
+        asImp_().updateTemperature_(elemCtx, dofIdx, timeIdx);
 
         if constexpr (enableBrine) {
             asImp_().updateSaltConcentration_(elemCtx, dofIdx, timeIdx);

diff --git a/opm/simulators/flow/BlackoilModelProperties.hpp b/opm/simulators/flow/BlackoilModelProperties.hpp
@@ -67,7 +67,7 @@ struct EnableSolvent<TypeTag, TTag::FlowProblem>
 
 template<class TypeTag>
 struct EnableTemperature<TypeTag, TTag::FlowProblem>
-{ static constexpr bool value = true; };
+{ static constexpr bool value = false; };
 
 template<class TypeTag>
 struct EnableEnergy<TypeTag, TTag::FlowProblem>

diff --git a/opm/simulators/flow/TemperatureModel.hpp b/opm/simulators/flow/TemperatureModel.hpp
@@ -157,6 +157,8 @@ class TemperatureModel : public GenericTemperatureModel<GetPropType<TypeTag, Pro
         const unsigned int numCells = simulator_.model().numTotalDof();
         for (unsigned globI = 0; globI < numCells; ++globI) {
             intQuants_[globI] = simulator_.model().intensiveQuantities(globI, /*timeIdx*/ 0);
+            intQuants_[globI].updateTemperature_(simulator_.problem(), globI, /*timeIdx*/ 0);
+            intQuants_[globI].updateEnergyQuantities_(simulator_.problem(), globI, /*timeIdx*/ 0);
         }
         advanceTemperatureFields();
     }