Description

Can we add the inverse kinetics method for "intercalation kinetics": "linear"?

Motivation

Running linear kinetics along with "surface form": "false".

Possible Implementation

I think it should be as simple as just sub-classing the InverseButlerVolmer class and deleting the arcsinh function from

  def _get_overpotential(self, j, j0, ne, T, u):
      return (2 * (self.param.R * T) / self.param.F / ne) * pybamm.arcsinh(
          j / (2 * j0 * u)
      )

plus some logic for handling model options.

Additional context

No response