Even though highly relevant, the coupled effect of current flow, heat generation and capacity fade have not been adequately examined in the literature. While several works have simulated the temperature distribution in large format cells [1,2] or capacity fade for cells [3] under isothermal conditions, very few have examined their coupled effect. Studies considering the coupled effect of capacity fade with heat generation have incorporated a pseudo-2D electrochemical degradation model [4,5]. The thermal model ranges from one to three-dimensional, with higher dimensional thermal models considering uniform heat generation within the cell. However, in a practical scenario, a large-format cell suffers from a non-uniform temperature distribution, leading to non-uniform electrochemical reactions and degradation.
Hence, a detailed coupled thermo-electrochemical, capacity-fade model is required to understand cell degradation and temperature rise during its operational life. In a step towards this goal, a two-dimensional physics-based, coupled thermo-electrochemical model with capacity degradation will be demonstrated for cylindrical lithium-ion cells. The electrochemical model with capacity fade is based on the porous-electrode and concentrated solution theories [6]. The thermal model considers the effect of ohmic heat in various cell components and the reversible and irreversible heat of reactions in the electrodes. The contribution of side reactions in heat generation is incorporated. The effect of changing porosity and thermal and electronic impedance with ageing on the cell is considered in the model. The results will give a better understanding of (a) the safe operation of the cell as the internal temperature of the cell changes with ageing, and of (b) the dependence of cell ageing on temperature. Fig.1 shows the spatial distribution of temperature at the end of discharge for an LMO cell under different convective heat transfer coefficients. The cells were subjected to a 5C current at an ambient temperature of 25°C. While the temperature reached for a convective heat transfer coefficient h=5 W/m2.K is much higher, the temperature gradient developed is more significant for h=50 W/m2.K. This can be attributed to the increased Biot number for the cell under forced cooling. A difference of 8°C for h=50 W/m2.K between the internal and external cell temperature show that the capacity-fade model needs to consider the local temperature distribution within the cell for better prediction of cell degradation with cycling.
References:
[1] Pan, Y. W., Hua, Y., Zhou, S., He, R., Zhang, Y., Yang, S., Liu, X., Lian, Y., Yan, X. & Wu, B. (2020). Journal of Power Sources, 459, 228070.
[2] Zhao, Y., Diaz, L. B., Patel, Y., Zhang, T., & Offer, G. J. (2019). Journal of The Electrochemical Society, 166(13), A2849.
[3] Atalay, S., Sheikh, M., Mariani, A., Merla, Y., Bower, E., & Widanage, W. D. (2020). Journal of Power Sources, 478, 229026.
[4] Xu, M., Wang, R., Zhao, P., & Wang, X. (2019). Journal of Power Sources, 438, 227015.
[5] Jiang, G., Zhuang, L., Hu, Q., Liu, Z., & Huang, J. (2020). Applied Thermal Engineering, 171, 115080.
[6] Rashid, M., & Gupta, A. (2014). ECS Electrochemistry Letters, 3(10), A95.