#729 started experimental tortuosity class

Author: Scottmar93

pybamm/models/submodels/tortuosity/experimental_tortuosity.py

@@ -0,0 +1,37 @@
+#
+# Class for Experimentally supplied tortuosity
+#
+import pybamm
+
+from .base_tortuosity import BaseModel
+
+
+class Experimental(BaseModel):
+    """Submodel for experimental tortuosity
+
+    **Extends:** :class:`pybamm.tortuosity.BaseModel`
+    """
+
+    def __init__(self, param, phase, set_leading_order=False):
+        super().__init__(param, phase)
+        self.set_leading_order = set_leading_order
+
+    def get_coupled_variables(self, variables):
+        param = self.param
+
+        if self.phase == "Electrolyte":
+            eps_n, eps_s, eps_p = variables["Porosity"].orphans
+            tor = pybamm.Concatenation(
+                eps_n ** param.b_e_n, eps_s ** param.b_e_s, eps_p ** param.b_e_p
+            )
+        elif self.phase == "Electrode":
+            eps_n, eps_s, eps_p = variables["Active material volume fraction"].orphans
+            tor = pybamm.Concatenation(
+                eps_n ** param.b_s_n, eps_s ** param.b_s_s, eps_p ** param.b_s_p
+            )
+
+        variables.update(
+            self._get_standard_tortuosity_variables(tor, self.set_leading_order)
+        )
+
+        return variables