Add parameter Conserve inner energy and use it

Author: Tor Harald Sandve

opm/simulators/flow/FlowProblemBlackoil.hpp

@@ -334,6 +334,13 @@ class FlowProblemBlackoil : public FlowProblem<TypeTag>
 
         this->initFluidSystem_();
 
+        // conserve inner energy instead of enthalpy if TEMP is used
+        // or THERMAL and parameter ConserveInnerEnergyThermal is true (default false)
+        bool isThermal = eclState.getSimulationConfig().isThermal();
+        bool isTemp = eclState.getSimulationConfig().isTemp();
+        bool conserveInnerEnergy = isTemp || (isThermal && Parameters::Get<Parameters::ConserveInnerEnergyThermal>());
+        FluidSystem::setEnergyEqualEnthalpy(conserveInnerEnergy);
+
         if (FluidSystem::phaseIsActive(FluidSystem::oilPhaseIdx) &&
             FluidSystem::phaseIsActive(FluidSystem::gasPhaseIdx)) {
             this->maxOilSaturation_.resize(this->model().numGridDof(), 0.0);

opm/simulators/flow/FlowProblemParameters.cpp

@@ -72,6 +72,10 @@ void registerFlowProblemParameters()
     Parameters::Register<Parameters::NumSatfuncConsistencySamplePoints>
         ("Maximum number of reported failures for each individual saturation function consistency check");
 
+    Parameters::Register<Parameters::ConserveInnerEnergyThermal>
+        ("Conserve inner energy and not enthalpy "
+         "even if THERMAL is used.");
+
     // By default, stop it after the universe will probably have stopped
     // to exist. (the ECL problem will finish the simulation explicitly
     // after it simulated the last episode specified in the deck.)
@@ -90,6 +94,8 @@ void registerFlowProblemParameters()
     // the default for the allowed volumetric error for oil per second
     Parameters::SetDefault<Parameters::NewtonTolerance<Scalar>>(1e-2);
     Parameters::SetDefault<Parameters::EnableGravity>(true);
+
+    Parameters::SetDefault<Parameters::ConserveInnerEnergyThermal>(false);
 }
 
 template void registerFlowProblemParameters<double>();

opm/simulators/flow/FlowProblemParameters.hpp

@@ -56,6 +56,9 @@ struct OutputMode { static constexpr auto value = "all"; };
 // between writing restart files
 struct RestartWritingInterval { static constexpr int value = 0xffffff; }; // disable
 
+// Conserve inner energy instead of enthalpy even if THERMAL is used
+struct ConserveInnerEnergyThermal { static constexpr bool value = false; };
+
 } // namespace Opm::Parameters
 
 namespace Opm {