A novel method for alleviating numerical stiffness in Li-ion thermal abuse models
Numerical modeling of thermal runaway in Lithium-ion batteries has become a critical tool for designing safer battery systems. Significant progress has been made in developing kinetic mechanisms for decomposition reactions and including additional physics such as venting and combustion. However, the...
Main Authors: | , , |
---|---|
Format: | Article |
Language: | English |
Published: |
Elsevier
2023-08-01
|
Series: | Journal of Power Sources Advances |
Subjects: | |
Online Access: | http://www.sciencedirect.com/science/article/pii/S266624852300015X |
_version_ | 1827866600684912640 |
---|---|
author | Jason Ostanek Mohammad Parhizi Judith Jeevarajan |
author_facet | Jason Ostanek Mohammad Parhizi Judith Jeevarajan |
author_sort | Jason Ostanek |
collection | DOAJ |
description | Numerical modeling of thermal runaway in Lithium-ion batteries has become a critical tool for designing safer battery systems. Significant progress has been made in developing kinetic mechanisms for decomposition reactions and including additional physics such as venting and combustion. However, the governing heat conduction equation and decomposition reaction equations become numerically stiff during thermal runaway, which limits the utility of thermal abuse models to low-dimensional formulations. The present work introduces a new solution strategy, which switches from the full, 3D transient heat conduction formulation to an adiabatic, 0D lumped body formulation only during the stiff portion of the simulation, i.e., only during thermal runaway. To test the new solver, a 3D thermal abuse model was configured to simulate an oven test of an 18650-format cell. The new solver was exercised for scenarios of varying degrees of stiffness, and the results were compared with a baseline solver using typical integration methods. For an extremely stiff scenario, computation speed was increased by a factor of 183x relative to the baseline solver, with little impact on solution accuracy, thus effectively alleviating the numerical stiffness issue. The new solution strategy addresses the poor scalability of high-dimensional models, such as 3D-CFD-based thermal abuse models, and improves their practicality for industrial use. |
first_indexed | 2024-03-12T15:07:52Z |
format | Article |
id | doaj.art-0a0f8ac54fdb4d8c94183ce68cfcc2e9 |
institution | Directory Open Access Journal |
issn | 2666-2485 |
language | English |
last_indexed | 2024-03-12T15:07:52Z |
publishDate | 2023-08-01 |
publisher | Elsevier |
record_format | Article |
series | Journal of Power Sources Advances |
spelling | doaj.art-0a0f8ac54fdb4d8c94183ce68cfcc2e92023-08-12T04:34:51ZengElsevierJournal of Power Sources Advances2666-24852023-08-0123100123A novel method for alleviating numerical stiffness in Li-ion thermal abuse modelsJason Ostanek0Mohammad Parhizi1Judith Jeevarajan2School of Engineering Technology, Purdue University, 401 N. Grant St., West Lafayette, IN, 47907, USA; Corresponding author.Electrochemical Safety Research Institute, UL Research Institutes, 5000 Gulf Fwy, UHTB, Houston, TX, 77204, USAElectrochemical Safety Research Institute, UL Research Institutes, 5000 Gulf Fwy, UHTB, Houston, TX, 77204, USANumerical modeling of thermal runaway in Lithium-ion batteries has become a critical tool for designing safer battery systems. Significant progress has been made in developing kinetic mechanisms for decomposition reactions and including additional physics such as venting and combustion. However, the governing heat conduction equation and decomposition reaction equations become numerically stiff during thermal runaway, which limits the utility of thermal abuse models to low-dimensional formulations. The present work introduces a new solution strategy, which switches from the full, 3D transient heat conduction formulation to an adiabatic, 0D lumped body formulation only during the stiff portion of the simulation, i.e., only during thermal runaway. To test the new solver, a 3D thermal abuse model was configured to simulate an oven test of an 18650-format cell. The new solver was exercised for scenarios of varying degrees of stiffness, and the results were compared with a baseline solver using typical integration methods. For an extremely stiff scenario, computation speed was increased by a factor of 183x relative to the baseline solver, with little impact on solution accuracy, thus effectively alleviating the numerical stiffness issue. The new solution strategy addresses the poor scalability of high-dimensional models, such as 3D-CFD-based thermal abuse models, and improves their practicality for industrial use.http://www.sciencedirect.com/science/article/pii/S266624852300015XLi-ion cellsNumerical simulationThermal runaway |
spellingShingle | Jason Ostanek Mohammad Parhizi Judith Jeevarajan A novel method for alleviating numerical stiffness in Li-ion thermal abuse models Journal of Power Sources Advances Li-ion cells Numerical simulation Thermal runaway |
title | A novel method for alleviating numerical stiffness in Li-ion thermal abuse models |
title_full | A novel method for alleviating numerical stiffness in Li-ion thermal abuse models |
title_fullStr | A novel method for alleviating numerical stiffness in Li-ion thermal abuse models |
title_full_unstemmed | A novel method for alleviating numerical stiffness in Li-ion thermal abuse models |
title_short | A novel method for alleviating numerical stiffness in Li-ion thermal abuse models |
title_sort | novel method for alleviating numerical stiffness in li ion thermal abuse models |
topic | Li-ion cells Numerical simulation Thermal runaway |
url | http://www.sciencedirect.com/science/article/pii/S266624852300015X |
work_keys_str_mv | AT jasonostanek anovelmethodforalleviatingnumericalstiffnessinliionthermalabusemodels AT mohammadparhizi anovelmethodforalleviatingnumericalstiffnessinliionthermalabusemodels AT judithjeevarajan anovelmethodforalleviatingnumericalstiffnessinliionthermalabusemodels AT jasonostanek novelmethodforalleviatingnumericalstiffnessinliionthermalabusemodels AT mohammadparhizi novelmethodforalleviatingnumericalstiffnessinliionthermalabusemodels AT judithjeevarajan novelmethodforalleviatingnumericalstiffnessinliionthermalabusemodels |