From ffc354252f005609e0082b1944105052093737a8 Mon Sep 17 00:00:00 2001
From: Henry Weller <http://openfoam.org>
Date: Thu, 18 Oct 2018 12:10:28 +0100
Subject: [PATCH] equationOfState::perfectFluid: Added enthalpy departure
 function

Updated steamInjection to solve for internal energy for the water.
---
 .../specie/equationOfState/perfectFluid/perfectFluidI.H       | 4 +++-
 .../steamInjection/constant/thermophysicalProperties.water    | 4 +---
 2 files changed, 4 insertions(+), 4 deletions(-)

diff --git a/src/thermophysicalModels/specie/equationOfState/perfectFluid/perfectFluidI.H b/src/thermophysicalModels/specie/equationOfState/perfectFluid/perfectFluidI.H
index 1ec512663e..b05343df32 100644
--- a/src/thermophysicalModels/specie/equationOfState/perfectFluid/perfectFluidI.H
+++ b/src/thermophysicalModels/specie/equationOfState/perfectFluid/perfectFluidI.H
@@ -95,7 +95,9 @@ inline Foam::scalar Foam::perfectFluid<Specie>::rho(scalar p, scalar T) const
 template<class Specie>
 inline Foam::scalar Foam::perfectFluid<Specie>::H(scalar p, scalar T) const
 {
-    return 0;
+    const scalar RT = this->R()*T;
+
+    return -rho0_*RT/(rho0_ + p/RT);
 }
 
 
diff --git a/tutorials/multiphase/reactingTwoPhaseEulerFoam/laminar/steamInjection/constant/thermophysicalProperties.water b/tutorials/multiphase/reactingTwoPhaseEulerFoam/laminar/steamInjection/constant/thermophysicalProperties.water
index 5471a07c5e..b9341f6376 100644
--- a/tutorials/multiphase/reactingTwoPhaseEulerFoam/laminar/steamInjection/constant/thermophysicalProperties.water
+++ b/tutorials/multiphase/reactingTwoPhaseEulerFoam/laminar/steamInjection/constant/thermophysicalProperties.water
@@ -23,11 +23,9 @@ thermoType
     thermo          hRefConst;
     equationOfState perfectFluid;
     specie          specie;
-    energy          sensibleEnthalpy;
+    energy          sensibleInternalEnergy;
 }
 
-dpdt yes;
-
 pressureWorkAlphaLimit 0;
 
 mixture