Running a set of simulations

[MatHunt]

So I want to see how the terminal voltage and the overpotential change when you change the separator porosity. So I thought I would do this whole thing in a loop. So I did the following(for the terminal voltage:
`import numpy as np
import pybamm
model = pybamm.lithium_ion.DFN()
parameter_values = pybamm.ParameterValues("Chen2020")
por = [0.1,0.2,0.3,0.4,0.5,0.6,0.7,0.8,0.9,1.0]

for i in range(len(por)):
parameter_values["Separator porosity"] = por[i]
sim = pybamm.Simulation(model, parameter_values=parameter_values)
sim.solve([0, 3600])
sim.plot("Terminal voltage [V]")`

The program falls over when asked to plot, and I have no idea why, as it seems to work for the actual simulation. The second part, is how do I obtain the overpotentials?

[rtimms]

You need square brackets around the list of variables you want to plot, e.g. sim.plot(["Terminal voltage [V]"]).

It would be faster to do this using InputParameters, then you only need to set the simulation up once

import pybamm

model = pybamm.lithium_ion.DFN()

parameter_values = pybamm.ParameterValues("Chen2020")
parameter_values["Separator porosity"] = "[input]"

sim = pybamm.Simulation(model, parameter_values=parameter_values)

pors = [0.1, 0.3, 0.5, 0.7, 0.9]
sols = []
for por in pors:
    sol = sim.solve([0, 3600], inputs={"Separator porosity": por})
    sols.append(sol)

pybamm.dynamic_plot(sols, ["Terminal voltage [V]"], labels=pors)

[rtimms]

To find other variables you can do model.variables.search("blah")

[MatHunt]

It was the first thing I did, but I couldn't find the variable overpotential, so I assumed that it was called something else.