From e5bf7d89d23454909b2e35f807304af380584285 Mon Sep 17 00:00:00 2001
From: sergio <s.ferraris@opencfd.co.uk>
Date: Tue, 12 Apr 2022 13:45:55 -0700
Subject: [PATCH] ENH: Adding thermo types to Lee model

---
 .../InterfaceCompositionModels.C              | 53 +++++++++++++++++++
 .../solidSpecie/Make/options                  |  3 +-
 .../include/solidThermoPhysicsTypes.H         | 15 ++++++
 3 files changed, 70 insertions(+), 1 deletion(-)

diff --git a/src/phaseSystemModels/multiphaseInter/phasesSystem/InterfaceCompositionModel/InterfaceCompositionModels.C b/src/phaseSystemModels/multiphaseInter/phasesSystem/InterfaceCompositionModel/InterfaceCompositionModels.C
index 4c934613f53..e60470d49e7 100644
--- a/src/phaseSystemModels/multiphaseInter/phasesSystem/InterfaceCompositionModel/InterfaceCompositionModels.C
+++ b/src/phaseSystemModels/multiphaseInter/phasesSystem/InterfaceCompositionModel/InterfaceCompositionModels.C
@@ -299,6 +299,20 @@ namespace Foam
             tabulatedThermoPhysics
         );
 
+        // From pure solid tabulated to pure tabulated flow
+        makeInterfacePureType
+        (
+            Lee,
+            heSolidThermo,
+            solidThermo,
+            pureMixture,
+            hTabulatedPolyIcoSolidThermoPhysics,
+            heRhoThermo,
+            rhoThermo,
+            pureMixture,
+            tabulatedThermoPhysics
+        );
+
         // From pure phase (exp-Transp, hPower solidThermo) to phase (ico-rho)
         makeInterfacePureType
         (
@@ -427,6 +441,45 @@ namespace Foam
             hConstSolidThermoPhysics
         );
 
+        makeInterfacePureType
+        (
+            Lee,
+            heRhoThermo,
+            rhoThermo,
+            pureMixture,
+            constRhoHThermoPhysics,
+            heSolidThermo,
+            solidThermo,
+            pureMixture,
+            hTabulatedPolyIcoSolidThermoPhysics
+        );
+
+        makeInterfacePureType
+        (
+            Lee,
+            heRhoThermo,
+            rhoThermo,
+            pureMixture,
+            tabulatedThermoPhysics,
+            heSolidThermo,
+            solidThermo,
+            pureMixture,
+            hTabulatedPolyIcoSolidThermoPhysics
+        );
+
+        makeInterfacePureType
+        (
+            Lee,
+            heRhoThermo,
+            rhoThermo,
+            pureMixture,
+            icoPoly8HThermoPhysics,
+            heSolidThermo,
+            solidThermo,
+            pureMixture,
+            hTabulatedPolyIcoSolidThermoPhysics
+        );
+
 
 
         // interfaceHeatResistance model definitions
diff --git a/src/thermophysicalModels/solidSpecie/Make/options b/src/thermophysicalModels/solidSpecie/Make/options
index bdc8b074c03..26cf19b5aca 100644
--- a/src/thermophysicalModels/solidSpecie/Make/options
+++ b/src/thermophysicalModels/solidSpecie/Make/options
@@ -1,5 +1,6 @@
 EXE_INC = \
-    -I$(LIB_SRC)/thermophysicalModels/specie/lnInclude
+    -I$(LIB_SRC)/thermophysicalModels/specie/lnInclude \
+    -I$(LIB_SRC)/thermophysicalModels/thermophysicalProperties/lnInclude
 
 LIB_LIBS = \
     -lspecie
diff --git a/src/thermophysicalModels/solidSpecie/include/solidThermoPhysicsTypes.H b/src/thermophysicalModels/solidSpecie/include/solidThermoPhysicsTypes.H
index d607e82af9a..525c34d7452 100644
--- a/src/thermophysicalModels/solidSpecie/include/solidThermoPhysicsTypes.H
+++ b/src/thermophysicalModels/solidSpecie/include/solidThermoPhysicsTypes.H
@@ -40,10 +40,12 @@ Description
 #include "hConstThermo.H"
 #include "hPolynomialThermo.H"
 #include "hPowerThermo.H"
+#include "hTabulatedThermo.H"
 #include "constIsoSolidTransport.H"
 #include "constAnIsoSolidTransport.H"
 #include "exponentialSolidTransport.H"
 #include "polynomialSolidTransport.H"
+#include "tabulatedSolidTransport.H"
 
 #include "sensibleEnthalpy.H"
 #include "thermo.H"
@@ -121,6 +123,19 @@ namespace Foam
             >,
             8
         > hPolyTranspPolyIcoSolidThermoPhysics;
+
+     typedef
+        tabulatedSolidTransport
+        <
+            species::thermo
+            <
+                hTabulatedThermo
+                <
+                    icoPolynomial<specie, 8>
+                >,
+                sensibleEnthalpy
+            >
+        > hTabulatedPolyIcoSolidThermoPhysics;
 }
 
 
-- 
GitLab