# Models

**URL:** https://pybamm.discourse.group/c/pybamm-qa/models/7.md

[Latest](https://pybamm.discourse.group/latest.md) · [Categories](https://pybamm.discourse.group/categories.md) · [Tags](https://pybamm.discourse.group/tags.md)

---

## [About the Models category](https://pybamm.discourse.group/t/about-the-models-category/14)

<div class="topic-metadata">

**Author:** [@pybamm](https://pybamm.discourse.group/u/pybamm)\
**Replies:** 0\
**Last updated:** [6 November 2024 14:08 UTC](https://pybamm.discourse.group/t/about-the-models-category/14 "2024-11-06T14:08:48Z")

</div>

Subcategory to ask questions about PyBaMM built-in models. For models not included into PyBaMM (i.e. community models) please refer to the relevant post in the “Model sharing” category.

---

## [Dimensionless form of the DFN model II](https://pybamm.discourse.group/t/dimensionless-form-of-the-dfn-model-ii/487)

<div class="topic-metadata">

**Author:** [@isaac\_basil](https://pybamm.discourse.group/u/isaac_basil)\
**Replies:** 0\
**Last updated:** [8 October 2026 08:55 UTC](https://pybamm.discourse.group/t/dimensionless-form-of-the-dfn-model-ii/487 "2026-10-08T08:55:34Z")

</div>

Hi, There’s a previous thread on the dimensionless form of the DFN. It says that the dimensionless form is no longer used. However, in the fickian\_diffusion.py file (which I believe is the default submodel for the DFN) …

---

## [Simulations with O’Regan 2022 parameter set not working at 2C 0°C, 2C 10°C and WLTP 0°C](https://pybamm.discourse.group/t/simulations-with-o-regan-2022-parameter-set-not-working-at-2c-0-c-2c-10-c-and-wltp-0-c/456)

<div class="topic-metadata">

**Author:** [@ashkal](https://pybamm.discourse.group/u/ashkal)\
**Replies:** 2\
**Last updated:** [17 April 2026 15:32 UTC](https://pybamm.discourse.group/t/simulations-with-o-regan-2022-parameter-set-not-working-at-2c-0-c-2c-10-c-and-wltp-0-c/456 "2026-04-17T15:32:44Z")

</div>

Hello, I am working with the O’Regan parameter set to get temperature prediction for LGM50 to compare the predictions with the model that I have developed. I am using this code as the reference: https://docs.pybamm.org…

---

## [Discrepancies in OCV-SOC Curves: Comparing Full-Cell Discharge Method vs. Electrode OCP Reconstruction](https://pybamm.discourse.group/t/discrepancies-in-ocv-soc-curves-comparing-full-cell-discharge-method-vs-electrode-ocp-reconstruction/325)

<div class="topic-metadata">

**Author:** [@pkalbhor](https://pybamm.discourse.group/u/pkalbhor)\
**Replies:** 4\
**Last updated:** [5 March 2026 13:31 UTC](https://pybamm.discourse.group/t/discrepancies-in-ocv-soc-curves-comparing-full-cell-discharge-method-vs-electrode-ocp-reconstruction/325 "2026-03-05T13:31:40Z")

</div>

I am calculating OCV-SOV curve for full cell from two different methods Using slow discharge expriment Using OCP curves of each electrodes Both are very different as shown in below image. What am I missing here?

---

## [Modifying DFN model to add convection term](https://pybamm.discourse.group/t/modifying-dfn-model-to-add-convection-term/445)

<div class="topic-metadata">

**Author:** [@kate2525](https://pybamm.discourse.group/u/kate2525)\
**Replies:** 0\
**Last updated:** [28 January 2026 21:59 UTC](https://pybamm.discourse.group/t/modifying-dfn-model-to-add-convection-term/445 "2026-01-28T21:59:42Z")

</div>

For a cell with electrolyte flow through (axial flow) and connected to a well mixed external tank is it even possible to modify the existing DFN model to introduce the species balance for the external tank and modify th…

---

## [State vector entries for composite electrode](https://pybamm.discourse.group/t/state-vector-entries-for-composite-electrode/405)

<div class="topic-metadata">

**Author:** [@DanielNicolae](https://pybamm.discourse.group/u/DanielNicolae)\
**Replies:** 1\
**Last updated:** [4 November 2025 14:52 UTC](https://pybamm.discourse.group/t/state-vector-entries-for-composite-electrode/405 "2025-11-04T14:52:13Z")

</div>

What do the last two entries in the state vector represent for an SPM/SPMe with a composite anode? var\_pts = { "x\_n": 20, "x\_s": 20, "x\_p": 20, "r\_n": 20, "r\_n\_prim": 20, "r\_n\_sec": 20, "r\_p…

---

## [How to disable already existing events in a model?](https://pybamm.discourse.group/t/how-to-disable-already-existing-events-in-a-model/201)

<div class="topic-metadata">

**Author:** [@alighato](https://pybamm.discourse.group/u/alighato)\
**Replies:** 1\
**Last updated:** [26 September 2025 14:22 UTC](https://pybamm.discourse.group/t/how-to-disable-already-existing-events-in-a-model/201 "2025-09-26T14:22:39Z")

</div>

I am trying to run a Single Particle Model (SPM) simulation and want to allow the cell voltage to go beyond the usual operating limits (i.e., I do not want the simulation to stop once the voltage crosses the typical lowe…

---

## [Correct Model and Solver pipeline assembly for better performance](https://pybamm.discourse.group/t/correct-model-and-solver-pipeline-assembly-for-better-performance/394)

<div class="topic-metadata">

**Author:** [@StephenV](https://pybamm.discourse.group/u/StephenV)\
**Replies:** 0\
**Last updated:** [9 September 2025 06:55 UTC](https://pybamm.discourse.group/t/correct-model-and-solver-pipeline-assembly-for-better-performance/394 "2025-09-09T06:55:19Z")

</div>

Hello Dear Colleagues, First of all, thank you for PyBaMM! I really enjoy working with it. I have a question regarding performance optimization and the proper assembly of a pipeline (model + solver). I’ve been followin…

---

## [Hysteresis in ECM model?](https://pybamm.discourse.group/t/hysteresis-in-ecm-model/387)

<div class="topic-metadata">

**Author:** [@lefile](https://pybamm.discourse.group/u/lefile)\
**Replies:** 0\
**Last updated:** [26 August 2025 11:54 UTC](https://pybamm.discourse.group/t/hysteresis-in-ecm-model/387 "2025-08-26T11:54:20Z")

</div>

Hi Everyone! Just a quick question - can I use hysteresis sub-model within equivalent-circuit model? I’d like to define charge and discharge OCVs. Best regards!

---

## [Modifying experiment without rebuilding sim](https://pybamm.discourse.group/t/modifying-experiment-without-rebuilding-sim/327)

<div class="topic-metadata">

**Author:** [@senzhan](https://pybamm.discourse.group/u/senzhan)\
**Replies:** 0\
**Last updated:** [15 May 2025 12:48 UTC](https://pybamm.discourse.group/t/modifying-experiment-without-rebuilding-sim/327 "2025-05-15T12:48:42Z")

</div>

I read that you can modify the current function without rebuilding the simulation again. This saves a lot of time when you need to run a series of simulations. Is there a way to modify experiment without rebuilding simul…

---

## [Accelerating SEI Formation](https://pybamm.discourse.group/t/accelerating-sei-formation/319)

<div class="topic-metadata">

**Author:** [@AlasdairCrawford](https://pybamm.discourse.group/u/AlasdairCrawford)\
**Replies:** 0\
**Last updated:** [1 May 2025 23:31 UTC](https://pybamm.discourse.group/t/accelerating-sei-formation/319 "2025-05-01T23:31:54Z")

</div>

Hello, I want to be able to simulate one cycle to represent multiple ageing cycles (For example to simulate a battery aging over 800 cycles, I might only simulate 8 cycles with ageing factor of 100). For reference, this…

---

## [Compatibility Issue with SEI on Cracks in Composite Model](https://pybamm.discourse.group/t/compatibility-issue-with-sei-on-cracks-in-composite-model/307)

<div class="topic-metadata">

**Author:** [@Hesam\_Tavakoli](https://pybamm.discourse.group/u/Hesam_Tavakoli)\
**Replies:** 0\
**Last updated:** [19 April 2025 19:06 UTC](https://pybamm.discourse.group/t/compatibility-issue-with-sei-on-cracks-in-composite-model/307 "2025-04-19T19:06:14Z")

</div>

Hi, I’ve encountered an issue while working with PyBaMM, specifically with the SEI on cracks model when used in conjunction with a composite model. The following error was raised: KeyError: "'X-averaged negative partic…

---

## [Lumped Thermal Model - Can I directly enter a cell level effective heat capacity?](https://pybamm.discourse.group/t/lumped-thermal-model-can-i-directly-enter-a-cell-level-effective-heat-capacity/234)

<div class="topic-metadata">

**Author:** [@ChaddKiggins](https://pybamm.discourse.group/u/ChaddKiggins)\
**Replies:** 1\
**Last updated:** [4 April 2025 00:32 UTC](https://pybamm.discourse.group/t/lumped-thermal-model-can-i-directly-enter-a-cell-level-effective-heat-capacity/234 "2025-04-04T00:32:59Z")

</div>

Hi all, I have a question regarding the lumped thermal model. Currently, it is formulated by defining the thickness and densities of the anode (film & current collector), separator, & cathode (film & current collector).…

---

## [Importing previous simulation output into a new simulation using an interpolant](https://pybamm.discourse.group/t/importing-previous-simulation-output-into-a-new-simulation-using-an-interpolant/236)

<div class="topic-metadata">

**Author:** [@isaac\_basil](https://pybamm.discourse.group/u/isaac_basil)\
**Replies:** 5\
**Last updated:** [31 March 2025 06:57 UTC](https://pybamm.discourse.group/t/importing-previous-simulation-output-into-a-new-simulation-using-an-interpolant/236 "2025-03-31T06:57:19Z")

</div>

Hi, I am trying to solve 2 simulations sequentially, and use an output from the first simulation as an input to the second. This output is a function of both space and time, j(x,t). I have 2 questions: I’m trying to …

---

## [Particle phases default option](https://pybamm.discourse.group/t/particle-phases-default-option/218)

<div class="topic-metadata">

**Author:** [@fmaletic](https://pybamm.discourse.group/u/fmaletic)\
**Replies:** 1\
**Last updated:** [11 March 2025 00:29 UTC](https://pybamm.discourse.group/t/particle-phases-default-option/218 "2025-03-11T00:29:23Z")

</div>

Hi, I’ve noticed that there is a difference in simulation results when using "particle phases": ("1", "1") and with the default option, which, from what I understand from the source code, is just "particle phases": "…

---

## [How to disable already existing events in a model?](https://pybamm.discourse.group/t/how-to-disable-already-existing-events-in-a-model/203)

<div class="topic-metadata">

**Author:** [@alighato](https://pybamm.discourse.group/u/alighato)\
**Replies:** 2\
**Last updated:** [6 March 2025 18:38 UTC](https://pybamm.discourse.group/t/how-to-disable-already-existing-events-in-a-model/203 "2025-03-06T18:38:42Z")

</div>

Hello everyone, I am trying to run a Single Particle Model (SPM) simulation and want to allow the cell voltage to go beyond the usual operating limits (i.e., I do not want the simulation to stop once the voltage crosses…

---

## [Modify the plating.py issue](https://pybamm.discourse.group/t/modify-the-plating-py-issue/197)

<div class="topic-metadata">

**Author:** [@travis\_wu](https://pybamm.discourse.group/u/travis_wu)\
**Replies:** 0\
**Last updated:** [21 February 2025 02:24 UTC](https://pybamm.discourse.group/t/modify-the-plating-py-issue/197 "2025-02-21T02:24:10Z")

</div>

I want to set j\_stripping to 0 when eta\_plating is greater than 0, while keeping the rest of the calculations unchanged. So, I modified the code as follows: I have confirmed that the logic is correct, but I am unable…

---

## [Question on the DFN surface formulation model equations](https://pybamm.discourse.group/t/question-on-the-dfn-surface-formulation-model-equations/170)

<div class="topic-metadata">

**Author:** [@tosh07](https://pybamm.discourse.group/u/tosh07)\
**Replies:** 0\
**Last updated:** [13 February 2025 12:42 UTC](https://pybamm.discourse.group/t/question-on-the-dfn-surface-formulation-model-equations/170 "2025-02-13T12:42:56Z")

</div>

Hi Everyone, This is just a question about the model equations for the ‘surface formulation’ of the DFN model. I am running some simple simulations to compare this formulation with other model types. I was just curious…

---

## [Discussing/understanding the Marquis2019 SPMe Asymptotic Derivation](https://pybamm.discourse.group/t/discussing-understanding-the-marquis2019-spme-asymptotic-derivation/88)

<div class="topic-metadata">

**Author:** [@amundmr](https://pybamm.discourse.group/u/amundmr)\
**Replies:** 1\
**Last updated:** [7 February 2025 15:48 UTC](https://pybamm.discourse.group/t/discussing-understanding-the-marquis2019-spme-asymptotic-derivation/88 "2025-02-07T15:48:20Z")

</div>

Hi all! Hope you had a great vacation. I hope I am not off-topic by asking math & model derivation questions here, but I am trying to get to the crux of Marquis et al 2019’s paper which derives the SPMe model that is us…

---

## [Clarification on Reproducing Composite Battery Degradation Model in PyBaMM (Bonkile et al. work)](https://pybamm.discourse.group/t/clarification-on-reproducing-composite-battery-degradation-model-in-pybamm-bonkile-et-al-work/150)

<div class="topic-metadata">

**Author:** [@DineshAdepu](https://pybamm.discourse.group/u/DineshAdepu)\
**Replies:** 0\
**Last updated:** [6 February 2025 10:18 UTC](https://pybamm.discourse.group/t/clarification-on-reproducing-composite-battery-degradation-model-in-pybamm-bonkile-et-al-work/150 "2025-02-06T10:18:29Z")

</div>

Dear all, I am attempting to reproduce the composite battery degradation model presented in Bonkile et al. \[1\]. The paper states that the modeling work has been integrated into the main PyBaMM repository. However, whil…

---

## [Positive electrode porosity](https://pybamm.discourse.group/t/positive-electrode-porosity/136)

<div class="topic-metadata">

**Author:** [@algo-eng](https://pybamm.discourse.group/u/algo-eng)\
**Replies:** 0\
**Last updated:** [31 January 2025 07:48 UTC](https://pybamm.discourse.group/t/positive-electrode-porosity/136 "2025-01-31T07:48:10Z")

</div>

Selected SEI solvent diffusion limited, film resistance, porosity change true For LFP cell After first cycle it is showing positive electrode porosity 1 Why it is?

---

## [Error plotting the different contributions to LLI against throughput capacity in coupled degradation model](https://pybamm.discourse.group/t/error-plotting-the-different-contributions-to-lli-against-throughput-capacity-in-coupled-degradation-model/76)

<div class="topic-metadata">

**Author:** [@AkilaJ](https://pybamm.discourse.group/u/AkilaJ)\
**Replies:** 1\
**Last updated:** [9 January 2025 10:36 UTC](https://pybamm.discourse.group/t/error-plotting-the-different-contributions-to-lli-against-throughput-capacity-in-coupled-degradation-model/76 "2025-01-09T10:36:55Z")

</div>

Hello Everyone, I used the following code to run the degradation submodels suggested by @DrSOKane, However, when I tried to obtain the plotting as per the example code I got an empty plot. Also, I cannot see the followi…

---

## [Extracting the heat source without running the thermal model](https://pybamm.discourse.group/t/extracting-the-heat-source-without-running-the-thermal-model/28)

<div class="topic-metadata">

**Author:** [@brosaplanella](https://pybamm.discourse.group/u/brosaplanella)\
**Replies:** 1\
**Last updated:** [13 November 2024 14:39 UTC](https://pybamm.discourse.group/t/extracting-the-heat-source-without-running-the-thermal-model/28 "2024-11-13T14:39:31Z")

</div>

This topic has been transferred from Github discussions. I’m running an electrochemical model (SPM) for now and am looking to calculate the heat source Q described in the thermal lumped model. I don’t need to calculate …

---

## [Can I create a custom ECM model?](https://pybamm.discourse.group/t/can-i-create-a-custom-ecm-model/21)

<div class="topic-metadata">

**Author:** [@brosaplanella](https://pybamm.discourse.group/u/brosaplanella)\
**Replies:** 1\
**Last updated:** [11 November 2024 15:51 UTC](https://pybamm.discourse.group/t/can-i-create-a-custom-ecm-model/21 "2024-11-11T15:51:59Z")

</div>

This topic has been transferred from Github discussions. I know pybamm offers a 1RC and 2RC Thevenin model. However is it possible to use pybamm’s Base Model class to create something simple like a Rint model or higher …
