#729 updated diffusion

Author: Scottmar93

examples/scripts/ecker_set.py

@@ -3,7 +3,7 @@
 pb.set_logging_level("INFO")
 
 model = pb.lithium_ion.DFN()
-model.convert_to_format = "python"
+# model.convert_to_format = "python"
 
 chemistry = pb.parameter_sets.Ecker2015
 # chemistry = pb.parameter_sets.Marquis2019
@@ -13,4 +13,4 @@
 
 solver = pb.CasadiSolver(mode="safe")
 sim.solve(solver=solver)
-sim.plot(["Discharge capacity [A.h]", "Terminal voltage [V]"])
+sim.plot()

input/parameters/lithium-ion/anodes/graphite_Ecker2015/graphite_electrolyte_reaction_rate_Ecker2015.py

@@ -30,6 +30,11 @@ def graphite_electrolyte_reaction_rate_Ecker2015(T, T_inf, E_r, R_g):
     :`numpy.Array`
         Reaction rate
     """
+
+    # from part I (just after Figure 11)
+    k_0_dm = 4.58 * 10 ** (-9)  # [dm4.5.s-1.mol-0.5]
+    k_0_m = k_0_dm * (0.1) ** 4.5  # [m4.5.s-1.mol-0.5]
+
     m_ref = 2.79 * 10 ** (-6)
     arrhenius = exp(E_r / R_g * (1 / T_inf - 1 / T))
 

input/parameters/lithium-ion/cathodes/LiNiCoO2_Ecker2015/nco_diffusivity_Ecker2015.py

@@ -1,5 +1,4 @@
 from pybamm import exp
-from pybamm import FunctionParameter
 
 
 def nco_diffusivity_Ecker2015(sto, T, T_inf, E_D_s, R_g):