Issue 3339 dead lithium

CHANGELOG.md

@@ -2,6 +2,7 @@
 
 ## Bug fixes
 
+- The irreversible plating model now increments `f"{Domain} dead lithium concentration [mol.m-3]"`, not `f"{Domain} lithium plating concentration [mol.m-3]"` as it did previously. ([#3485](https://github.com/pybamm-team/PyBaMM/pull/3485))
 - Fixed a bug where the JaxSolver would fails when using GPU support with no input parameters ([#3423](https://github.com/pybamm-team/PyBaMM/pull/3423))
 
 # [v23.9rc0](https://github.com/pybamm-team/PyBaMM/tree/v23.9rc0) - 2023-10-31

pybamm/models/submodels/interface/lithium_plating/plating.py

@@ -115,18 +115,26 @@ def set_rhs(self, variables):
             ]
             L_sei = variables[f"{Domain} total SEI thickness [m]"]
 
-        # In the partially reversible plating model, coupling term turns reversible
-        # lithium into dead lithium. In other plating models, it is zero.
         lithium_plating_option = getattr(self.options, domain)["lithium plating"]
-        if lithium_plating_option == "partially reversible":
+        if lithium_plating_option == "reversible":
+            # In the reversible plating model, there is no dead lithium
+            dc_plated_Li = -a_j_stripping / self.param.F
+            dc_dead_Li = pybamm.Scalar(0)
+        elif lithium_plating_option == "irreversible":
+            # In the irreversible plating model, all plated lithium is dead lithium
+            dc_plated_Li = pybamm.Scalar(0)
+            dc_dead_Li = -a_j_stripping / self.param.F
+        elif lithium_plating_option == "partially reversible":
+            # In the partially reversible plating model, the coupling term turns
+            # reversible lithium into dead lithium over time.
             dead_lithium_decay_rate = self.param.dead_lithium_decay_rate(L_sei)
             coupling_term = dead_lithium_decay_rate * c_plated_Li
-        else:
-            coupling_term = pybamm.Scalar(0)
+            dc_plated_Li = -a_j_stripping / self.param.F - coupling_term
+            dc_dead_Li = coupling_term
 
         self.rhs = {
-            c_plated_Li: -a_j_stripping / self.param.F - coupling_term,
-            c_dead_Li: coupling_term,
+            c_plated_Li: dc_plated_Li,
+            c_dead_Li: dc_dead_Li,
         }
 
     def set_initial_conditions(self, variables):