Refactor "Inner SEI lithium interstitial diffusivity [m2.s-1]" in Arrhenius form

[Mingzefei]

Description

I want to refactor "Inner SEI lithium interstitial diffusivity [m2.s-1]" in Arrhenius form. I have found the reference value for the activation energy of this reaction and intend to refactor this param as a function of temperature.

The reason I have identified is that "Inner SEI lithium interstitial diffusivity [m2.s-1]" is hardcoded and cannot be rewritten as a function.

Motivation

All other diffusivity params are in function form, but the SEI diffusivity is in parameter form, which is somewhat strange.
And I am very keen to achieve this functionality.

Possible Implementation

I am considering delving into the local pybamm code for modifications.
Is there a better way to achieve this? Or can someone provide some guidance?

Additional context

No response

[kratman]

@Mingzefei There are a few PRs working on SEI and other degradation mechanisms at the moment. You might want to start by taking a look at some of them to see how your idea fits into the current progress:

• make SEI single layer only #4470
• Domain degradation parameters #4463
• Composite porosity #4417

These three PRs are all linked to some of the current degradation work. Once you see what is changing, then you can start to refactor and put up a PR.

PyBaMM also has a slack channel and some other forums where you could discuss it further, if needed