Simulations were conducted in the Python Battery Mathematical Modelling (PyBaMM) software package (using v0.4.0). The equations for the thermal-electrochemical model are summarised in Table and Table . In order to solve the model a finite volume scheme was used, with 30 grid points for each electrode and the separator, and 150 grid points for each particle; resulting in a system of 9092 ODEs and 150 algebraic equations. An exponential mesh was used to help with the convergence of the solver. In order to solve the system, a CasADI solver was used. Each simulation of discharge plus relaxation takes 10 to 20 minutes using an Intel Core i7-7660U (2.50GHz) processor and 16 GB RAM. This is because the nonlinear diffusion takes the solver many iterations to converge. The computational time could be reduced by using more sophisticated numerical methods, but this is out of the scope of this work. It should be noted that reduced models such as the Thermal Single Particle Model with electrolyte (TSPMe) yield very similar results with a significant reduction of the computational time.