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

TL;DR

This paper presents a 1D finite volume model using coupled meshes to simulate potential and species distribution in a Li||Bi liquid metal battery, highlighting the effects of mass transport on cell performance and voltage.

Contribution

A novel coupled mesh finite volume model for Li||Bi liquid metal batteries that captures electrochemical and mass transport effects without hydrodynamics.

Findings

Mass transport significantly affects overpotentials and cell voltage.

Diffusion and migration currents can counteract each other, influencing voltage.

Simulated limiting current density is lower than experimental, indicating convective effects are important.

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. It is openly available at https://doi.org/10.14278/rodare.2313.