Simulation of potential and species distribution in a Li||Bi liquid metal battery using coupled meshes

Abstract

In this work a 1D finite volume based model using coupled meshes is introduced to capture potential and species distribution throughout the discharge process in a lithium–bismuth liquid metal battery while neglecting hydrodynamic effects, focusing on the electrochemical properties of the cell and the mass transport in electrolyte and cathode. Interface reactions in the electrical double layer are considered through the introduction of a discrete jump of the potential modelled as periodic boundary condition to resolve interfacial discontinuities in the cell potential. A balanced-force like approach is implemented to ensure consistent calculation at the interface level. It is found that mass transport and concentration gradients have a significant effect on the cell overpotentials and thus on cell performance and cell voltage. By quantifying overvoltages in the Li Bi cell with a mixed cation electrolyte, it is possible to show that diffusion and migration current density could have counteractive effects on the cell voltage. Furthermore, the simulated limiting current density is observed to be much lower than experimentally measured, which can be attributed to convective effects in the electrolyte that need to be addressed in future simulations. The solver is based on the open source library OpenFOAM and thoroughly verified against the equivalent system COMSOL multiphysics and further validated with experimental results.