How to quickly plot parameter as a function during a given input and also during simulation?

[RuiheLi]

Dear PyBaMM team,

I wonder are there any quick ways available to plot out what the parameter as a function defined in the PyBaMM parameter sets looks like for a given input (like c_e = np.linspace(0, 1000, 100))? Also is it possible to plot the in situ parameter, let's say, diffusivity, during the simulation? Currently I can only do this by copy paste the function from the input parameter set and manually call that function, which is not very convenient though.

Big thanks in advance!

[YannickNoelStephanKuhn]

Dear Ruihe,

from pybamm.parameters.lithium_ion_parameters import LithiumIonParameters
import numpy as np

D_e_symbol = LithiumIonParameters().D_e_dimensional
c_e_eval = np.linspace(0, 1000, 100)
parameter_set = pybamm.ParameterValues("Marquis2019")

D_e_eval = [
    parameter_set.process_symbol(D_e_symbol(c_e_point, 298.15)).evaluate()
    for c_e_point in c_e_eval
]

Same for c_e_eval being the array of electrolyte concentration solutions. Put that into a convenience function and it should be just what you need.

[RuiheLi]

Dear Yannick,

Thanks so much for the kind reply. That is really really helpful!

[rtimms]

You can also do

import pybamm

parameter_set = pybamm.ParameterValues("Marquis2019")

c_e = pybamm.linspace(0, 1000, 100)
D_e = parameter_set["Electrolyte diffusivity [m2.s-1]"]

pybamm.plot(c_e, D_e(c_e, 298.15))

to plot.

And

import pybamm 

model = pybamm.lithium_ion.SPMe()
c_e = model.variables["Electrolyte concentration [mol.m-3]"]
T = model.variables["Cell temperature [K]"]
D_e = model.param.D_e_dimensional
model.variables["Electrolyte diffusivity [m2.s-1]"] = D_e(c_e, T)

sim = pybamm.Simulation(model)
sim.solve([0, 3600])

sim.plot(["Electrolyte diffusivity [m2.s-1]", "Terminal voltage [V]"])

to add the diffusivity as a variable.

[RuiheLi]

Thanks so much for the nice answer Robert!