benchmark_isothermal_potential_pair.py

#
# Benchmark a "1+1+1D DFN" isothermal potential pair pouch cell model in PyBaMM
#
import pybamm
from benchmark import benchmark

# file to save times to
filename = "isothermal_potential_pair.pkl"

# number of finite volumes in current colector domain (must be 1 or >2)
npts = [1, 3, 4, 8, 16, 32, 64, 128, 256, 512]

# number of repeated solves
repeats = 3

# load 1+1D model (standard DFN)
dfn = pybamm.lithium_ion.DFN(
    options={},
    name="1+1D DFN",
)

# load 1+1+1D model
dfn_pouch = pybamm.lithium_ion.DFN(
    options={
        "current collector": "potential pair",
        "dimensionality": 1,
    },
    name="1+1+1D DFN",
)

# benchmark
benchmark(dfn, dfn_pouch, filename, npts, repeats)