Abstract
Silicon-based composite electrodes in lithium ion batteries attract increasing attention because of their high theoretical capacity. Here, numerical simulations are used to better understand the interplay between electrochemical and morphological behavior of the silicon-graphite (1:2.7) composite electrode during galvanostatic cycling. Finite element methodology is used to solve a one-dimensional model based on the porous electrode and concentrated solution theory. Porosity changes in the silicon electrode and solid electrolyte interphase layer growth are also included in the model. The simulation results show that at lower rates, the electrode with high initial porosity is being fully utilized before the lower cut-off potential is reached. When comparing the computational results with experimental observations, it can be seen that the main reason for capacity fade is the increase in tortuosity in the diffusion pathway of lithium ions due to cracking of the silicon composite electrode upon electrochemical cycling.