AUTHOR=Li Weisi , León Quiroga Vanessa , Crompton K. R. , Ostanek Jason K. 

TITLE=High Resolution 3-D Simulations of Venting in 18650 Lithium-Ion Cells

JOURNAL=Frontiers in Energy Research

VOLUME=9

YEAR=2021

URL=https://www.frontiersin.org/journals/energy-research/articles/10.3389/fenrg.2021.788239

DOI=10.3389/fenrg.2021.788239

ISSN=2296-598X

ABSTRACT=<p>High temperature gases released through the safety vent of a lithium-ion cell during a thermal runaway event contain flammable components that, if ignited, can increase the risk of thermal runaway propagation to other cells in a multi-cell pack configuration. Computational fluid dynamics (CFD) simulations of flow through detailed geometric models of four vent-activated commercial 18650 lithium-ion cell caps were conducted using two turbulence modeling approaches: Reynolds-averaged Navier-Stokes (RANS) and scale-resolving simulations (SRS). The RANS method was compared with independent experiments of discharge coefficient through the cap across a range of pressure ratios and then used to investigate the ensemble-averaged flow field for the four caps. At high pressure ratios, choked flow occurs either at the current collector plate when flow through the current collector plate is more restrictive or the positive terminal vent holes when flow through the current collector plate is less restrictive. Turbulent mixing occurred within the vent cap assembly, in the jets emerging from the vent holes, and in recirculating zones directly above the vent cap assembly. The global maximum turbulent viscosity ratio (<inline-formula id="inf1"><mml:math id="m1" xmlns:mml="http://www.w3.org/1998/Math/MathML"><mml:mrow><mml:msub><mml:mi>μ</mml:mi><mml:mi>T</mml:mi></mml:msub><mml:mo>/</mml:mo><mml:mi>μ</mml:mi></mml:mrow></mml:math></inline-formula>) of the MTI, LG MJ1, K2, and LG M36 caps at pressure ratio of <italic>P</italic><sub><italic>1</italic></sub><italic>/P</italic><sub><italic>2</italic></sub> = 7 were 4,575, 3,360, 3,855, and 2,993, respectively. SRS and RANS simulations showed that both velocity magnitude and fluctuating velocity magnitude were lower for vent holes which are obstructed by the burst disk. SRS showed high levels of fluctuating velocity in the jets, up to 48.5% of the global maximum velocity. The present CFD models and the resulting insights provide the groundwork for future studies to investigate how jet structure and turbulence levels influence combustion and heat transfer in propagating thermal runaway scenarios.</p>