basicThermo: Added volScalarField overload of THE

applications/solvers/multiphase/compressibleInterFoam/twoPhaseMixtureThermo/twoPhaseMixtureThermo.C

@@ -262,6 +262,18 @@ Foam::tmp<Foam::volScalarField> Foam::twoPhaseMixtureThermo::hc() const
 }
 
 
+Foam::tmp<Foam::volScalarField> Foam::twoPhaseMixtureThermo::THE
+(
+    const volScalarField& h,
+    const volScalarField& p,
+    const volScalarField& T0
+) const
+{
+    NotImplemented;
+    return T0;
+}
+
+
 Foam::tmp<Foam::scalarField> Foam::twoPhaseMixtureThermo::THE
 (
     const scalarField& h,

applications/solvers/multiphase/compressibleInterFoam/twoPhaseMixtureThermo/twoPhaseMixtureThermo.H

@@ -217,6 +217,14 @@ public:
             //- Enthalpy of formation [J/kg]
             virtual tmp<volScalarField> hc() const;
 
+            //- Temperature from enthalpy/internal energy
+            virtual tmp<volScalarField> THE
+            (
+                const volScalarField& h,
+                const volScalarField& p,
+                const volScalarField& T0    // starting temperature
+            ) const;
+
             //- Temperature from enthalpy/internal energy for cell-set
             virtual tmp<scalarField> THE
             (

applications/solvers/multiphase/compressibleMultiphaseInterFoam/multiphaseMixtureThermo/multiphaseMixtureThermo.C

@@ -430,6 +430,18 @@ Foam::tmp<Foam::volScalarField> Foam::multiphaseMixtureThermo::hc() const
 }
 
 
+Foam::tmp<Foam::volScalarField> Foam::multiphaseMixtureThermo::THE
+(
+    const volScalarField& h,
+    const volScalarField& p,
+    const volScalarField& T0
+) const
+{
+    NotImplemented;
+    return T0;
+}
+
+
 Foam::tmp<Foam::scalarField> Foam::multiphaseMixtureThermo::THE
 (
     const scalarField& h,

applications/solvers/multiphase/compressibleMultiphaseInterFoam/multiphaseMixtureThermo/multiphaseMixtureThermo.H

@@ -339,6 +339,14 @@ public:
             //- Enthalpy of formation [J/kg]
             virtual tmp<volScalarField> hc() const;
 
+            //- Temperature from enthalpy/internal energy
+            virtual tmp<volScalarField> THE
+            (
+                const volScalarField& h,
+                const volScalarField& p,
+                const volScalarField& T0    // starting temperature
+            ) const;
+
             //- Temperature from enthalpy/internal energy for cell-set
             virtual tmp<scalarField> THE
             (

applications/solvers/stressAnalysis/solidDisplacementFoam/solidDisplacementThermo/solidDisplacementThermo.C

@@ -360,6 +360,18 @@ Foam::tmp<Foam::volScalarField> Foam::solidDisplacementThermo::hc() const
 }
 
 
+Foam::tmp<Foam::volScalarField> Foam::solidDisplacementThermo::THE
+(
+    const volScalarField& h,
+    const volScalarField& p,
+    const volScalarField& T0
+) const
+{
+    NotImplemented;
+    return tmp<volScalarField>(nullptr);
+}
+
+
 Foam::tmp<Foam::scalarField> Foam::solidDisplacementThermo::THE
 (
     const scalarField& he,

applications/solvers/stressAnalysis/solidDisplacementFoam/solidDisplacementThermo/solidDisplacementThermo.H

@@ -245,6 +245,14 @@ public:
             //- Enthalpy of formation [J/kg]
             virtual tmp<volScalarField> hc() const;
 
+            //- Temperature from enthalpy/internal energy
+            virtual tmp<volScalarField> THE
+            (
+                const volScalarField& h,
+                const volScalarField& p,
+                const volScalarField& T0    // starting temperature
+            ) const;
+
             //- Temperature from enthalpy/internal energy for cell-set
             virtual tmp<scalarField> THE
             (

src/thermophysicalModels/basic/basicThermo/basicThermo.H

@@ -343,6 +343,14 @@ public:
             //- Enthalpy of formation [J/kg]
             virtual tmp<volScalarField> hc() const = 0;
 
+            //- Temperature from enthalpy/internal energy
+            virtual tmp<volScalarField> THE
+            (
+                const volScalarField& h,
+                const volScalarField& p,
+                const volScalarField& T0    // starting temperature
+            ) const = 0;
+
             //- Temperature from enthalpy/internal energy for cell-set
             virtual tmp<scalarField> THE
             (

src/thermophysicalModels/basic/heThermo/heThermo.C

@@ -677,6 +677,28 @@ Foam::heThermo<BasicThermo, MixtureType>::CpByCpv() const
 }
 
 
+template<class BasicThermo, class MixtureType>
+Foam::tmp<Foam::volScalarField> Foam::heThermo<BasicThermo, MixtureType>::THE
+(
+    const volScalarField& h,
+    const volScalarField& p,
+    const volScalarField& T0
+) const
+{
+    return volScalarFieldProperty
+    (
+        "T",
+        dimTemperature,
+        &MixtureType::cellThermoMixture,
+        &MixtureType::patchFaceThermoMixture,
+        &MixtureType::thermoMixtureType::THE,
+        h,
+        p,
+        T0
+    );
+}
+
+
 template<class BasicThermo, class MixtureType>
 Foam::tmp<Foam::scalarField> Foam::heThermo<BasicThermo, MixtureType>::THE
 (

src/thermophysicalModels/basic/heThermo/heThermo.H

@@ -277,6 +277,14 @@ public:
             //- Enthalpy of formation [J/kg]
             virtual tmp<volScalarField> hc() const;
 
+            //- Temperature from enthalpy/internal energy
+            virtual tmp<volScalarField> THE
+            (
+                const volScalarField& h,
+                const volScalarField& p,
+                const volScalarField& T0    // starting temperature
+            ) const;
+
             //- Temperature from enthalpy/internal energy for cell-set
             virtual tmp<scalarField> THE
             (