Repository: pybamm-team/PyBaMM

Issue:
Add parameters and functions required for thermal modelling
As title


Discussion:

Patch:
diff --git a/docs/source/parameters/index.rst b/docs/source/parameters/index.rst
index 3212b4ef0..bce29d859 100644
--- a/docs/source/parameters/index.rst
+++ b/docs/source/parameters/index.rst
@@ -6,6 +6,7 @@ Parameters
   parameter_values
   geometric_parameters
   electrical_parameters
+  thermal_parameters
   standard_parameters_lithium_ion
   standard_parameters_lead_acid
   print_parameters
diff --git a/docs/source/parameters/thermal_parameters.rst b/docs/source/parameters/thermal_parameters.rst
new file mode 100644
index 000000000..d06ef55ed
--- /dev/null
+++ b/docs/source/parameters/thermal_parameters.rst
@@ -0,0 +1,4 @@
+Thermal Parameters
+=====================
+
+.. automodule:: pybamm.parameters.thermal_parameters
diff --git a/input/parameters/lithium-ion/graphite_entropic_change_Moura.py b/input/parameters/lithium-ion/graphite_entropic_change_Moura.py
new file mode 100644
index 000000000..1ca0626f5
--- /dev/null
+++ b/input/parameters/lithium-ion/graphite_entropic_change_Moura.py
@@ -0,0 +1,33 @@
+import autograd.numpy as np
+
+
+def graphite_entropic_change_Moura(sto, c_n_max):
+    """
+        Graphite entropic change in open circuit potential (OCP) at a temperature of
+        298.15K as a function of the stochiometry taken from Scott Moura's FastDFN code
+        [1].
+
+        References
+        ----------
+        .. [1] https://github.com/scott-moura/fastDFN
+
+          Parameters
+          ----------
+          sto: double
+               Stochiometry of material (li-fraction)
+
+    """
+
+    du_dT = (
+        -1.5 * (120.0 / c_n_max) * np.exp(-120 * sto)
+        + (0.0351 / (0.083 * c_n_max)) * ((np.cosh((sto - 0.286) / 0.083)) ** (-2))
+        - (0.0045 / (0.119 * c_n_max)) * ((np.cosh((sto - 0.849) / 0.119)) ** (-2))
+        - (0.035 / (0.05 * c_n_max)) * ((np.cosh((sto - 0.9233) / 0.05)) ** (-2))
+        - (0.0147 / (0.034 * c_n_max)) * ((np.cosh((sto - 0.5) / 0.034)) ** (-2))
+        - (0.102 / (0.142 * c_n_max)) * ((np.cosh((sto - 0.194) / 0.142)) ** (-2))
+        - (0.022 / (0.0164 * c_n_max)) * ((np.cosh((sto - 0.9) / 0.0164)) ** (-2))
+        - (0.011 / (0.0226 * c_n_max)) * ((np.cosh((sto - 0.124) / 0.0226)) ** (-2))
+        + (0.0155 / (0.029 * c_n_max)) * ((np.cosh((sto - 0.105) / 0.029)) ** (-2))
+    )
+
+    return du_dT
diff --git a/input/parameters/lithium-ion/lico2_entropic_change_Moura.py b/input/parameters/lithium-ion/lico2_entropic_change_Moura.py
new file mode 100644
index 000000000..246221dc0
--- /dev/null
+++ b/input/parameters/lithium-ion/lico2_entropic_change_Moura.py
@@ -0,0 +1,30 @@
+import autograd.numpy as np
+
+
+def lico2_entropic_change_Moura(sto, c_p_max):
+    """
+        Lithium Cobalt Oxide (LiCO2) entropic change in open circuit potential (OCP) at
+        a temperature of 298.15K as a function of the stochiometry. The fit is taken
+        from Scott Moura's FastDFN code [1].
+
+        References
+        ----------
+        .. [1] https://github.com/scott-moura/fastDFN
+
+          Parameters
+          ----------
+          sto: double
+               Stochiometry of material (li-fraction)
+
+    """
+
+    du_dT = (
+        0.07645 * (-54.4806 / c_p_max) * ((1.0 / np.cosh(30.834 - 54.4806 * sto)) ** 2)
+        + 2.1581 * (-50.294 / c_p_max) * ((np.cosh(52.294 - 50.294 * sto)) ** (-2))
+        + 0.14169 * (19.854 / c_p_max) * ((np.cosh(11.0923 - 19.8543 * sto)) ** (-2))
+        - 0.2051 * (5.4888 / c_p_max) * ((np.cosh(1.4684 - 5.4888 * sto)) ** (-2))
+        - (0.2531 / 0.1316 / c_p_max) * ((np.cosh((-sto + 0.56478) / 0.1316)) ** (-2))
+        - (0.02167 / 0.006 / c_p_max) * ((np.cosh((sto - 0.525) / 0.006)) ** (-2))
+    )
+
+    return du_dT
diff --git a/input/parameters/lithium-ion/mcmb2528_lif6-in-ecdmc_lico2_parameters_Dualfoil.csv b/input/parameters/lithium-ion/mcmb2528_lif6-in-ecdmc_lico2_parameters_Dualfoil.csv
index 38ed6a4aa..0be21766c 100644
--- a/input/parameters/lithium-ion/mcmb2528_lif6-in-ecdmc_lico2_parameters_Dualfoil.csv
+++ b/input/parameters/lithium-ion/mcmb2528_lif6-in-ecdmc_lico2_parameters_Dualfoil.csv
@@ -53,7 +53,35 @@ Reference OCP vs SHE in the positive electrode [V],,,
 Charge transfer coefficient,0.5,Scott Moura FastDFN,
 Double-layer capacity [F.m-2],0.2,,
 ,,,
+# Density,,,
+Negative current collector density [kg.m-3], 8954,, 
+Negative electrode density [kg.m-3], 1657,, 
+Separator density [kg.m-3], 397,, 
+Positive electrode density [kg.m-3], 3262,, 
+Positive current collector density [kg.m-3], 2707,, 
+,,,
+# Specific heat capacity
+Negative current collector specific heat capacity [J.kg-1.K-1], 385,,
+Negative electrode specific heat capacity [J.kg-1.K-1], 700,,
+Separator specific heat capacity [J.kg-1.K-1], 700,,
+Negative electrode specific heat capacity [J.kg-1.K-1], 700,,
+Positive current collector specific heat capacity [J.kg-1.K-1], 897,,
+,,,
+# Thermal conductivity
+Negative current collector thermal conductivity [W.m-1.K-1], 401,, 
+Negative electrode thermal conductivity [W.m-1.K-1], 1.7,, 
+Separator thermal conductivity [W.m-1.K-1], 0.16,, 
+Positive electrode thermal conductivity [W.m-1.K-1], 2.1,, 
+Positive current collector thermal conductivity [W.m-1.K-1], 237,, 
+,,,
+# Thermal parameters
+Heat transfer coefficient [W.m-2.K-1], 10,,
+Typical temperature variation [K], 2.4,, 
+Lumped effective thermal density [J.K-1.m-3], 1.8116E6,,
+Effective thermal conductivity [W.m-1.K-1], 59.3964,,
+
 # Initial Conditions,,,
 Initial concentration in electrolyte [mol.m-3],1000,Scott Moura FastDFN,
 Initial concentration in negative electrode [mol.m-3],19986.609595075,Scott Moura FastDFN,
 Initial concentration in positive electrode [mol.m-3],30730.7554385565,Scott Moura FastDFN,
+Initial temperature [K], 298.15,
diff --git a/pybamm/__init__.py b/pybamm/__init__.py
index ec9c2f948..91ab9e252 100644
--- a/pybamm/__init__.py
+++ b/pybamm/__init__.py
@@ -163,6 +163,7 @@ from .parameters.parameter_values import ParameterValues
 from .parameters import standard_current_functions
 from .parameters import geometric_parameters
 from .parameters import electrical_parameters
+from .parameters import thermal_parameters
 from .parameters import standard_parameters_lithium_ion, standard_parameters_lead_acid
 from .parameters.print_parameters import print_parameters, print_evaluated_parameters
 
diff --git a/pybamm/models/base_models.py b/pybamm/models/base_models.py
index f885f3143..811bc4f14 100644
--- a/pybamm/models/base_models.py
+++ b/pybamm/models/base_models.py
@@ -411,6 +411,12 @@ class StandardBatteryBaseModel(BaseModel):
                 "Positive electrode diffusivity": os.path.join(
                     input_path, "lico2_diffusivity_Dualfoil.py"
                 ),
+                "Negative electrode OCV entropic change": os.path.join(
+                    input_path, "graphite_entropic_change_Moura.py"
+                ),
+                "Positive electrode OCV entropic change": os.path.join(
+                    input_path, "lico2_entropic_change_Moura.py"
+                )
             },
         )
 
diff --git a/pybamm/parameters/electrical_parameters.py b/pybamm/parameters/electrical_parameters.py
index a450f8f0c..6fb44f63f 100644
--- a/pybamm/parameters/electrical_parameters.py
+++ b/pybamm/parameters/electrical_parameters.py
@@ -1,9 +1,6 @@
 #
-# Geometric Parameters
+# Standard electrical parameters
 #
-"""
-Standard electrical parameters
-"""
 import pybamm
 import numpy as np
 
@@ -16,8 +13,7 @@ def abs_non_zero(x):
 
 
 # --------------------------------------------------------------------------------------
-"Dimensional Parameters"
-# Electrical
+# Dimensional Parameters
 I_typ = pybamm.Parameter("Typical current [A]")
 Q = pybamm.Parameter("Cell capacity [A.h]")
 C_rate = abs(I_typ / Q)
diff --git a/pybamm/parameters/standard_parameters_lithium_ion.py b/pybamm/parameters/standard_parameters_lithium_ion.py
index baf30e699..b183566f5 100644
--- a/pybamm/parameters/standard_parameters_lithium_ion.py
+++ b/pybamm/parameters/standard_parameters_lithium_ion.py
@@ -84,6 +84,9 @@ c_p_init_dimensional = pybamm.Parameter(
     "Initial concentration in positive electrode [mol.m-3]"
 )
 
+# thermal
+Delta_T = pybamm.thermal_parameters.Delta_T
+
 # --------------------------------------------------------------------------------------
 "2. Dimensional Functions"
 
@@ -108,14 +111,28 @@ def D_p_dimensional(c_p):
     return pybamm.FunctionParameter("Positive electrode diffusivity", c_p)
 
 
-def U_n_dimensional(c_s_n):
+def U_n_dimensional(sto):
     "Dimensional open-circuit voltage in the negative electrode [V]"
-    return pybamm.FunctionParameter("Negative electrode OCV", c_s_n)
+    return pybamm.FunctionParameter("Negative electrode OCV", sto)
+
+
+def U_p_dimensional(sto):
+    "Dimensional open-circuit voltage in the positive electrode [V]"
+    return pybamm.FunctionParameter("Positive electrode OCV", sto)
 
 
-def U_p_dimensional(c_s_p):
-    "Dimensional open-circuit voltage of of the positive electrode [V]"
-    return pybamm.FunctionParameter("Positive electrode OCV", c_s_p)
+def dUdT_n_dimensional(sto):
+    "Dimensional entropic change of the negative electrode open circuit voltage [V.K-1]"
+    return pybamm.FunctionParameter(
+        "Negative electrode OCV entropic change", sto, c_n_max
+    )
+
+
+def dUdT_p_dimensional(sto):
+    "Dimensional entropic change of the positive electrode open circuit voltage [V.K-1]"
+    return pybamm.FunctionParameter(
+        "Positive electrode OCV entropic change", sto, c_p_max
+    )
 
 
 # can maybe improve ref value at some stage
@@ -152,6 +169,8 @@ tau_diffusion_e = L_x ** 2 / D_e_dimensional(c_e_typ)
 tau_diffusion_n = R_n ** 2 / D_n_dimensional(c_n_max)
 tau_diffusion_p = R_p ** 2 / D_p_dimensional(c_p_max)
 
+# Thermal diffusion timescale
+tau_th_yz = pybamm.thermal_parameters.tau_th_yz
 
 # --------------------------------------------------------------------------------------
 "4. Dimensionless Parameters"
@@ -161,8 +180,9 @@ C_p = tau_diffusion_p / tau_discharge
 C_e = tau_diffusion_e / tau_discharge
 C_r_n = tau_r_n / tau_discharge
 C_r_p = tau_r_p / tau_discharge
+C_th = tau_th_yz / tau_discharge
 
-# Concentrtion ratios
+# Concentration ratios
 gamma_e = c_e_typ / c_n_max
 gamma_p = c_p_max / c_n_max
 
@@ -217,10 +237,27 @@ voltage_high_cut = (
     voltage_high_cut_dimensional - (U_p_ref - U_n_ref)
 ) / potential_scale
 
+# Thermal
+rho_cn = pybamm.thermal_parameters.rho_cn
+rho_n = pybamm.thermal_parameters.rho_n
+rho_s = pybamm.thermal_parameters.rho_s
+rho_p = pybamm.thermal_parameters.rho_p
+rho_cp = pybamm.thermal_parameters.rho_cp
+
+lambda_cn = pybamm.thermal_parameters.lambda_cn
+lambda_n = pybamm.thermal_parameters.lambda_n
+lambda_s = pybamm.thermal_parameters.lambda_s
+lambda_p = pybamm.thermal_parameters.lambda_p
+lambda_cp = pybamm.thermal_parameters.lambda_cp
+
+Theta = pybamm.thermal_parameters.Theta
+h = pybamm.thermal_parameters.h
+
 # Initial conditions
 c_e_init = c_e_init_dimensional / c_e_typ
 c_n_init = c_n_init_dimensional / c_n_max
 c_p_init = c_p_init_dimensional / c_p_max
+T_init = pybamm.thermal_parameters.T_init
 
 
 # --------------------------------------------------------------------------------------
@@ -240,25 +277,37 @@ def kappa_e(c_e):
     return kappa_e_dimensional(c_e_dimensional) / kappa_scale
 
 
-def D_n(c_n):
+def D_n(c_s_n):
     "Dimensionless negative particle diffusivity"
-    c_n_dimensional = c_n * c_n_max
-    return D_n_dimensional(c_n_dimensional) / D_n_dimensional(c_n_max)
+    c_s_n_dimensional = c_s_n * c_n_max
+    return D_n_dimensional(c_s_n_dimensional) / D_n_dimensional(c_n_max)
 
 
-def D_p(c_p):
+def D_p(c_s_p):
     "Dimensionless positive particle diffusivity"
-    c_p_dimensional = c_p * c_p_max
-    return D_p_dimensional(c_p_dimensional) / D_p_dimensional(c_p_max)
+    c_s_p_dimensional = c_s_p * c_p_max
+    return D_p_dimensional(c_s_p_dimensional) / D_p_dimensional(c_p_max)
 
 
-def U_n(c_n):
-    "Dimensionless open-circuit sp.potential in the negative electrode"
-    sto = c_n
+def U_n(c_s_n):
+    "Dimensionless open-circuit potential in the negative electrode"
+    sto = c_s_n
     return (U_n_dimensional(sto) - U_n_ref) / potential_scale
 
 
-def U_p(c_p):
-    "Dimensionless open-circuit sp.potential in the positive electrode"
-    sto = c_p
+def U_p(c_s_p):
+    "Dimensionless open-circuit potential in the positive electrode"
+    sto = c_s_p
     return (U_p_dimensional(sto) - U_p_ref) / potential_scale
+
+
+def dUdT_n(c_s_n):
+    "Dimensionless entropic change in negative open-circuit potential"
+    sto = c_s_n
+    return dUdT_n_dimensional(sto) * Delta_T / potential_scale
+
+
+def dUdT_p(c_s_p):
+    "Dimensionless entropic change in positive open-circuit potential"
+    sto = c_s_p
+    return dUdT_p_dimensional(sto) * Delta_T / potential_scale
diff --git a/pybamm/parameters/thermal_parameters.py b/pybamm/parameters/thermal_parameters.py
new file mode 100644
index 000000000..02f92b9c9
--- /dev/null
+++ b/pybamm/parameters/thermal_parameters.py
@@ -0,0 +1,72 @@
+#
+# Standard thermal parameters
+#
+import pybamm
+
+# --------------------------------------------------------------------------------------
+# Dimensional parameters
+
+# Reference temperature
+T_ref = pybamm.Parameter("Reference temperature [K]")
+
+# Density
+rho_cn_dim = pybamm.Parameter("Negative current collector density [kg.m-3]")
+rho_n_dim = pybamm.Parameter("Negative electrode density [kg.m-3]")
+rho_s_dim = pybamm.Parameter("Separator density [kg.m-3]")
+rho_p_dim = pybamm.Parameter("Positive electrode density [kg.m-3]")
+rho_cp_dim = pybamm.Parameter("Positive current collector density [kg.m-3]")
+
+# Specific heat capacity
+c_p_cn_dim = pybamm.Parameter(
+    "Negative current collector specific heat capacity [J.kg-1.K-1]"
+)
+c_p_n_dim = pybamm.Parameter("Negative electrode specific heat capacity [J.kg-1.K-1]")
+c_p_s_dim = pybamm.Parameter("Separator specific heat capacity [J.kg-1.K-1]")
+c_p_p_dim = pybamm.Parameter("Negative electrode specific heat capacity [J.kg-1.K-1]")
+c_p_cp_dim = pybamm.Parameter(
+    "Positive current collector specific heat capacity [J.kg-1.K-1]"
+)
+
+# Thermal conductivity
+lambda_cn_dim = pybamm.Parameter(
+    "Negative current collector thermal conductivity [W.m-1.K-1]"
+)
+lambda_n_dim = pybamm.Parameter("Negative electrode thermal conductivity [W.m-1.K-1]")
+lambda_s_dim = pybamm.Parameter("Separator thermal conductivity [W.m-1.K-1]")
+lambda_p_dim = pybamm.Parameter("Positive electrode thermal conductivity [W.m-1.K-1]")
+lambda_cp_dim = pybamm.Parameter(
+    "Positive current collector thermal conductivity [W.m-1.K-1]"
+)
+
+# Thermal parameters
+h_dim = pybamm.Parameter("Heat transfer coefficient [W.m-2.K-1]")
+Delta_T = pybamm.Parameter("Typical temperature variation [K]")
+rho_eff_dim = pybamm.Parameter("Lumped effective thermal density [J.K-1.m-3]")
+lambda_eff_dim = pybamm.Parameter("Effective thermal conductivity [W.m-1.K-1]")
+
+# Initial temperature
+T_init_dim = pybamm.Parameter("Initial temperature [K]")
+
+# --------------------------------------------------------------------------------------
+# Timescales
+tau_th_yz = rho_eff_dim * (pybamm.geometric_parameters.L_z ** 2) / lambda_eff_dim
+
+# --------------------------------------------------------------------------------------
+# Dimensionless parameters
+
+rho_cn = rho_cn_dim * c_p_cn_dim / rho_eff_dim
+rho_n = rho_n_dim * c_p_n_dim / rho_eff_dim
+rho_s = rho_s_dim * c_p_s_dim / rho_eff_dim
+rho_p = rho_p_dim * c_p_p_dim / rho_eff_dim
+rho_cp = rho_cp_dim * c_p_cp_dim / rho_eff_dim
+
+lambda_cn = lambda_cn_dim / lambda_eff_dim
+lambda_n = lambda_n_dim / lambda_eff_dim
+lambda_s = lambda_s_dim / lambda_eff_dim
+lambda_p = lambda_p_dim / lambda_eff_dim
+lambda_cp = lambda_cp_dim / lambda_eff_dim
+
+Theta = Delta_T / T_ref
+h = h_dim * pybamm.geometric_parameters.L_x / lambda_eff_dim
+
+T_init = (T_init_dim - T_ref) / Delta_T


Test patch:
diff --git a/tests/unit/test_parameters/test_dimensionless_parameter_values_lithium_ion.py b/tests/unit/test_parameters/test_dimensionless_parameter_values_lithium_ion.py
index 8da1dd582..458718ffc 100644
--- a/tests/unit/test_parameters/test_dimensionless_parameter_values_lithium_ion.py
+++ b/tests/unit/test_parameters/test_dimensionless_parameter_values_lithium_ion.py
@@ -1,10 +1,8 @@
 #
-# Tests for the Base Parameter Values class
+# Tests lithium ion parameters load and give expected values
 #
 import pybamm
 
-import os
-
 import unittest
 import numpy as np
 
@@ -15,41 +13,7 @@ class TestDimensionlessParameterValues(unittest.TestCase):
         correctly for the specific set of parameters for LCO from dualfoil. The values
         are those converted from those in Scott's transfer which previous versions of
         the DFN work with. A 1C rate corresponds to a 24A/m^2 current density"""
-        input_path = os.path.join(os.getcwd(), "input", "parameters", "lithium-ion")
-        values = pybamm.ParameterValues(
-            os.path.join(
-                input_path, "mcmb2528_lif6-in-ecdmc_lico2_parameters_Dualfoil.csv"
-            ),
-            {
-                "Typical current [A]": 1,
-                "Current function": os.path.join(
-                    os.getcwd(),
-                    "pybamm",
-                    "parameters",
-                    "standard_current_functions",
-                    "constant_current.py",
-                ),
-                "Electrolyte diffusivity": os.path.join(
-                    input_path, "electrolyte_diffusivity_Capiglia1999.py"
-                ),
-                "Electrolyte conductivity": os.path.join(
-                    input_path, "electrolyte_conductivity_Capiglia1999.py"
-                ),
-                "Negative electrode OCV": os.path.join(
-                    input_path, "graphite_mcmb2528_ocp_Dualfoil.py"
-                ),
-                "Positive electrode OCV": os.path.join(
-                    input_path, "lico2_ocp_Dualfoil.py"
-                ),
-                "Negative electrode diffusivity": os.path.join(
-                    input_path, "graphite_mcmb2528_diffusivity_Dualfoil.py"
-                ),
-                "Positive electrode diffusivity": os.path.join(
-                    input_path, "lico2_diffusivity_Dualfoil.py"
-                ),
-            },
-        )
-
+        values = pybamm.LithiumIonBaseModel().default_parameter_values
         param = pybamm.standard_parameters_lithium_ion
 
         c_rate = param.i_typ / 24  # roughly for the numbers I used before
@@ -229,6 +193,81 @@ class TestDimensionlessParameterValues(unittest.TestCase):
             values.process_symbol(param.sigma_p * c_rate).evaluate(None, None), 47.57, 1
         )
 
+    def test_thermal_parameters(self):
+        values = pybamm.LithiumIonBaseModel().default_parameter_values
+        param = pybamm.standard_parameters_lithium_ion
+        c_rate = param.i_typ / 24
+
+        # Density
+        np.testing.assert_almost_equal(
+            values.process_symbol(param.rho_cn).evaluate(), 1.9019, 2
+        )
+        np.testing.assert_almost_equal(
+            values.process_symbol(param.rho_n).evaluate(), 0.6403, 2
+        )
+        np.testing.assert_almost_equal(
+            values.process_symbol(param.rho_s).evaluate(), 0.1535, 2
+        )
+        np.testing.assert_almost_equal(
+            values.process_symbol(param.rho_p).evaluate(), 1.2605, 2
+        )
+        np.testing.assert_almost_equal(
+            values.process_symbol(param.rho_cp).evaluate(), 1.3403, 2
+        )
+
+        # Thermal conductivity
+        np.testing.assert_almost_equal(
+            values.process_symbol(param.lambda_cn).evaluate(), 6.7513, 2
+        )
+        np.testing.assert_almost_equal(
+            values.process_symbol(param.lambda_n).evaluate(), 0.0296, 2
+        )
+        np.testing.assert_almost_equal(
+            values.process_symbol(param.lambda_s).evaluate(), 0.0027, 2
+        )
+        np.testing.assert_almost_equal(
+            values.process_symbol(param.lambda_p).evaluate(), 0.0354, 2
+        )
+        np.testing.assert_almost_equal(
+            values.process_symbol(param.lambda_cp).evaluate(), 3.9901, 2
+        )
+
+        # other thermal parameters
+
+        # note: in paper this is 0.0534 * c_rate which conflicts with this
+        # if we do C_th * c_rate we get 0.0534 so probably error in paper
+        np.testing.assert_almost_equal(
+            values.process_symbol(param.C_th / c_rate).evaluate(), 0.0253, 2
+        )
+
+        np.testing.assert_almost_equal(
+            values.process_symbol(param.Theta / c_rate).evaluate(), 0.008, 2
+        )
+
+        np.testing.assert_almost_equal(
+            values.process_symbol(param.h).evaluate(), 3.7881 * 10 ** (-5), 7
+        )
+
+        np.testing.assert_equal(values.process_symbol(param.T_init).evaluate(), 0)
+
+    def test_parameter_functions(self):
+        values = pybamm.LithiumIonBaseModel().default_parameter_values
+        param = pybamm.standard_parameters_lithium_ion
+
+        c_test = pybamm.Scalar(0.5)
+
+        values.process_symbol(param.U_n(c_test)).evaluate()
+        values.process_symbol(param.U_p(c_test)).evaluate()
+        values.process_symbol(param.dUdT_n(c_test)).evaluate()
+        values.process_symbol(param.dUdT_p(c_test)).evaluate()
+
+        values.process_symbol(param.D_p(c_test)).evaluate()
+        values.process_symbol(param.D_n(c_test)).evaluate()
+
+        c_e_test = pybamm.Scalar(1)
+        values.process_symbol(param.D_e(c_e_test)).evaluate()
+        values.process_symbol(param.kappa_e(c_e_test)).evaluate()
+
 
 if __name__ == "__main__":
     print("Add -v for more debug output")