add reaction heating (pybamm-team#4557)

* add reaction heating * style: pre-commit fixes --------- Co-authored-by: pre-commit-ci[bot] <66853113+pre-commit-ci[bot]@users.noreply.github.com> Co-authored-by: Eric G. Kratz <[email protected]>
kratman · Nov 4, 2024 · f4e4955 · f4e4955
1 parent 106c249
commit f4e4955
Showing 1 changed file with 6 additions and 2 deletions.
diff --git a/src/pybamm/models/submodels/thermal/base_thermal.py b/src/pybamm/models/submodels/thermal/base_thermal.py
@@ -143,8 +143,12 @@ def _get_standard_coupled_variables(self, variables):
             phase_names = ["primary ", "secondary "]
 
         if self.options.electrode_types["negative"] == "planar":
-            Q_rxn_n = pybamm.FullBroadcast(
-                0, ["negative electrode"], "current collector"
+            i_n = variables["Lithium metal total interfacial current density [A.m-2]"]
+            eta_r_n = variables["Lithium metal interface reaction overpotential [V]"]
+            Q_rxn_n = pybamm.PrimaryBroadcast(
+                i_n * eta_r_n / self.param.n.L,
+                ["negative electrode"],
+                "current collector",
             )
             Q_rev_n = pybamm.FullBroadcast(
                 0, ["negative electrode"], "current collector"