Extend pybamm.ParameterValues.create_from_bpx() to automate initial concentrations for blended electrodes

[ejfdickinson]

Description

Calling pybamm.ParameterValues.create_from_bpx() with a BPX file containing a blended electrode returns a warning as follows:

Initial concentrations cannot be set using stoichiometry limits for blend electrodes. Please set the initial concentrations manually.

The method could automate this calculation by solving the equations specified in the BPX standard.

Motivation

No response

Possible Implementation

No response

Additional context

No response

[rtimms]

Tagging #2682 as related. Currently the electrode state of health models in pybamm don't support blended electrodes.

[ejfdickinson]

@rtimms Off-topic I know, but for this reason I find it odd that calc_esoh=True is the default in pybamm.Simulation.solve().

[rtimms]

Yes, this is simply because calc_esoh predates blended electrodes. I’m open to changing this.

[valentinsulzer]

The proper fix, which we are planning to do, is to do "just in time" calculations of esoh values only when they are requested