From 75e6ec24cf5c54d3d06619c82232cf1fec43fb94 Mon Sep 17 00:00:00 2001
From: Thomas Fetzer <thomas.fetzer@iws.uni-stuttgart.de>
Date: Tue, 23 Feb 2016 13:52:28 +0100
Subject: [PATCH] [components] Implement thermal conductivities for components

---
 dumux/material/components/air.hh        |  6 +++++-
 dumux/material/components/co2.hh        | 14 ++++++++++++++
 dumux/material/components/mesitylene.hh | 13 +++++++++++++
 dumux/material/components/n2.hh         | 16 ++++++++++++++++
 dumux/material/components/o2.hh         | 16 ++++++++++++++++
 dumux/material/components/xylene.hh     | 14 ++++++++++++++
 6 files changed, 78 insertions(+), 1 deletion(-)

diff --git a/dumux/material/components/air.hh b/dumux/material/components/air.hh
index a5bb1f9fa5..1a5c751808 100644
--- a/dumux/material/components/air.hh
+++ b/dumux/material/components/air.hh
@@ -234,15 +234,19 @@ public:
 
     /*!
      * \brief Thermal conductivity \f$\mathrm{[[W/(m*K)]}\f$ of air.
+     *
      * Isobaric Properties for Nitrogen in: NIST Standard \cite NIST <BR>
      * evaluated at p=.1 MPa, T=20°C <BR>
      * Nitrogen: 0.025398 <BR>
      * Oxygen: 0.026105 <BR>
      * lambda_air is approximately 0.78*lambda_N2+0.22*lambda_O2
+     *
+     * \param temperature absolute temperature in \f$\mathrm{[K]}\f$
+     * \param pressure of the phase in \f$\mathrm{[Pa]}\f$
      */
     static Scalar gasThermalConductivity(Scalar temperature, Scalar pressure)
     {
-        return 0.0255535; // conductivity of pure air [W/(m K)]
+        return 0.0255535;
     }
 };
 
diff --git a/dumux/material/components/co2.hh b/dumux/material/components/co2.hh
index c0ad2c2e3e..c217e12c9c 100644
--- a/dumux/material/components/co2.hh
+++ b/dumux/material/components/co2.hh
@@ -347,6 +347,20 @@ public:
         // no difference for supercritical CO2
         return gasViscosity(temperature, pressure);
     }
+
+    /*!
+     * \brief Thermal conductivity \f$\mathrm{[[W/(m*K)]}\f$ of CO2.
+     *
+     * Thermal conductivity of CO2 at T=20°C, see:
+     * http://www.engineeringtoolbox.com/carbon-dioxide-d_1000.html
+     *
+     * \param temperature absolute temperature in \f$\mathrm{[K]}\f$
+     * \param pressure of the phase in \f$\mathrm{[Pa]}\f$
+     */
+    static Scalar gasThermalConductivity(Scalar temperature, Scalar pressure)
+    {
+        return 0.087;
+    }
 };
 
 template <class Scalar, class CO2Tables>
diff --git a/dumux/material/components/mesitylene.hh b/dumux/material/components/mesitylene.hh
index 32e355d32f..76ee24cee0 100644
--- a/dumux/material/components/mesitylene.hh
+++ b/dumux/material/components/mesitylene.hh
@@ -317,6 +317,19 @@ public:
         return (C6H5 + 3*CH3 - 2*H)/molarMass(); // J/(mol K) -> J/(kg K)
     }
 
+    /*!
+     * \brief Thermal conductivity \f$\mathrm{[[W/(m*K)]}\f$ of mesitylene
+     *
+     * see: http://pubs.acs.org/doi/pdf/10.1021/ci000139t
+     *
+     * \param temperature absolute temperature in \f$\mathrm{[K]}\f$
+     * \param pressure of the phase in \f$\mathrm{[Pa]}\f$
+     */
+    static Scalar liquidThermalConductivity( Scalar temperature,  Scalar pressure)
+    {
+        return 0.1351;
+    }
+
 protected:
     /*!
      * \brief The molar density of pure mesitylene at a given pressure and temperature
diff --git a/dumux/material/components/n2.hh b/dumux/material/components/n2.hh
index 4d5cb57be0..9df5eddd46 100644
--- a/dumux/material/components/n2.hh
+++ b/dumux/material/components/n2.hh
@@ -272,6 +272,22 @@ public:
         // convertion from micro poise to Pa s
         return mu/1e6 / 10;
     }
+
+    /*!
+     * \brief Thermal conductivity \f$\mathrm{[[W/(m*K)]}\f$ of nitrogen.
+     *
+     * Isobaric Properties for Nitrogen and Oxygen in: NIST Standard
+     * Reference Database Number 69, Eds. P.J. Linstrom and
+     * W.G. Mallard evaluated at p=.1 MPa, does not
+     * change dramatically with p and can be interpolated linearly with temperature
+     *
+     * \param temperature absolute temperature in \f$\mathrm{[K]}\f$
+     * \param pressure of the phase in \f$\mathrm{[Pa]}\f$
+     */
+    static Scalar gasThermalConductivity(Scalar temperature, Scalar pressure)
+    {
+        return 6.525e-5 * (temperature - 273.15) + 0.024031;
+    }
 };
 
 } // end namespace
diff --git a/dumux/material/components/o2.hh b/dumux/material/components/o2.hh
index df3ce0b4a8..28df1fcbd5 100644
--- a/dumux/material/components/o2.hh
+++ b/dumux/material/components/o2.hh
@@ -221,6 +221,22 @@ public:
         // convertion from micro poise to Pa s
         return mu/1e6 / 10;
     }
+
+    /*!
+     * \brief Thermal conductivity \f$\mathrm{[[W/(m*K)]}\f$ of nitrogen.
+     *
+     * Isobaric Properties for Nitrogen and Oxygen in: NIST Standard
+     * Reference Database Number 69, Eds. P.J. Linstrom and
+     * W.G. Mallard evaluated at p=.1 MPa, does not
+     * change dramatically with p and can be interpolated linearly with temperature
+     *
+     * \param temperature absolute temperature in \f$\mathrm{[K]}\f$
+     * \param pressure of the phase in \f$\mathrm{[Pa]}\f$
+     */
+    static Scalar gasThermalConductivity(Scalar temperature, Scalar pressure)
+    {
+        return 8.044e-5 * (temperature - 273.15) + 0.024486;
+    }
 };
 
 } // end namespace
diff --git a/dumux/material/components/xylene.hh b/dumux/material/components/xylene.hh
index 81ed4d1a83..785685281e 100644
--- a/dumux/material/components/xylene.hh
+++ b/dumux/material/components/xylene.hh
@@ -354,6 +354,20 @@ public:
 
         return r; // [Pa s]
     }
+
+    /*!
+     * \brief Thermal conductivity \f$\mathrm{[[W/(m*K)]}\f$ of xylene
+     *
+     * Thermal Conductivity of p-Xylene taken from the Dortmund Data Bank, see:
+     * http://www.ddbst.de/en/EED/PCP/TCN_C176.php
+     *
+     * \param temperature absolute temperature in \f$\mathrm{[K]}\f$
+     * \param pressure of the phase in \f$\mathrm{[Pa]}\f$
+     */
+    static Scalar liquidThermalConductivity( Scalar temperature,  Scalar pressure)
+    {
+        return 0.13;
+    }
 };
 
 } // end namespace
-- 
GitLab