Calendar ageing variable plotting

[Akila1993]

Hi All,

I am using following code for calendar ageing experiment with degradation models. Suppose that I want to plot the rate of change of SEI layer growth or Voltage drop (η_SEI) due to the SEI growth against time how I can modify the code.

Could anyone support me in this regard.

Thanks !

model = pybamm.lithium_ion.DFN(
    {
        "thermal": "lumped", #lumped thermal submodel
        "surface form": "differential", #Adding capacitance to the model
        "SEI":"solvent-diffusion limited", #"reaction limited",
        "SEI porosity change": "true",
        #"lithium plating": "partially reversible",
        #"lithium plating porosity change": "true",  # alias for "SEI porosity change"
        #"particle mechanics":"swelling only",
        #"SEI on cracks": "true",
        #"loss of active material": "stress-driven",

    }
)
pybamm.set_logging_level("INFO")

updated_parameters["Current function [A]"] = 0

sim = pybamm.Simulation(model, parameter_values=updated_parameters,solver=pybamm.CasadiSolver(mode='safe without grid'),var_pts=var_pts)

years = 1
days = years * 365
hours = days * 24
minutes = hours * 60
seconds = minutes * 60

t_eval = np.linspace(0, seconds, 100)

sol = sim.solve(t_eval=t_eval)