modeling and simulation of the thermal runaway of ... · u(t), i(t), soc t(t), q(t) l. cai and r....
TRANSCRIPT
![Page 1: Modeling and Simulation of the Thermal Runaway of ... · U(t), I(t), SOC T(t), Q(t) L. Cai and R. E. White, Journal of Power Sources 196 (2011) 5985–5989 Thermal conductivity Heat](https://reader035.vdocument.in/reader035/viewer/2022081404/5f064cfe7e708231d4174dd8/html5/thumbnails/1.jpg)
Institute for Applied Materials-Applied Materials Physics (IAM-AWP)
www.kit.edu
Thermal
runaway
Heat
KIT – The Research University in the Helmholtz Association
A. Melcher, C. Ziebert, B. Lei, W. Zhao, M. Rohde, H. J. Seifert
Modeling and Simulation of the Thermal Runaway
of Cylindrical Li-Ion Cells
![Page 2: Modeling and Simulation of the Thermal Runaway of ... · U(t), I(t), SOC T(t), Q(t) L. Cai and R. E. White, Journal of Power Sources 196 (2011) 5985–5989 Thermal conductivity Heat](https://reader035.vdocument.in/reader035/viewer/2022081404/5f064cfe7e708231d4174dd8/html5/thumbnails/2.jpg)
KIT, IAM-AWP2
Increase of safety and reliability of lithium-ion batteries for EV/HEV and stationary applications
Aim: Improvement of TMS and BMS by determination of quantitative data using battery calorimetry in combination with modelling and simulation
www.techatplay.com
Motivation
http://insideevs.comwww.eawg.ch
→ For improving battery management system (BMS) and thermal management system (TMS) electrochemical and thermal behavior of the cells have to be thoroughly studied
→ Overheating
→ Overcharge
→ Overdischarge
→ Short Circuit
→ Accident
Possible Safety Impacts
Dr. C. Ziebert – COMSOL Conference 2016 Munich, 12.10.-14.10.2016
![Page 3: Modeling and Simulation of the Thermal Runaway of ... · U(t), I(t), SOC T(t), Q(t) L. Cai and R. E. White, Journal of Power Sources 196 (2011) 5985–5989 Thermal conductivity Heat](https://reader035.vdocument.in/reader035/viewer/2022081404/5f064cfe7e708231d4174dd8/html5/thumbnails/3.jpg)
KIT, IAM-AWP3
B Dunn et al. Science 2011;334:928-935
Negative Electrode
Oxidation
𝐿𝑖 → 𝐿𝑖+ + 𝑒−
Positive Electrode
Reduction
𝐿𝑖+ + 𝑒−→ 𝐿𝑖
Working principle of a Lithium-ion cell
Discharge process
Dr. C. Ziebert – COMSOL Conference 2016 Munich, 12.10.-14.10.2016
![Page 4: Modeling and Simulation of the Thermal Runaway of ... · U(t), I(t), SOC T(t), Q(t) L. Cai and R. E. White, Journal of Power Sources 196 (2011) 5985–5989 Thermal conductivity Heat](https://reader035.vdocument.in/reader035/viewer/2022081404/5f064cfe7e708231d4174dd8/html5/thumbnails/4.jpg)
KIT, IAM-AWP4
Tomography images of typical cylindrical 18650 Cells
Thin layers of active materials and separators are rolled to a cylindrical cell
Positive Electrode: e.g. LiCoO2 or LiMn2O4 on 10-25 mm thick Al foil
Negative Electrode: graphite on 10-12 mm thick Cu foil
Separator: 16-25 mm thick polyolefine membrane (PE, PP, PE/PP)
Dr. C. Ziebert – COMSOL Conference 2016 Munich, 12.10.-14.10.2016
![Page 5: Modeling and Simulation of the Thermal Runaway of ... · U(t), I(t), SOC T(t), Q(t) L. Cai and R. E. White, Journal of Power Sources 196 (2011) 5985–5989 Thermal conductivity Heat](https://reader035.vdocument.in/reader035/viewer/2022081404/5f064cfe7e708231d4174dd8/html5/thumbnails/5.jpg)
5 KIT, IAM-AWP
General scheme of the modeling of
the thermal runaway of a Li-ion cell
Thermal Model
Electrochemical
Model
Exothermic Model
(new contribution)
Dr. C. Ziebert – COMSOL Conference 2016 Munich, 12.10.-14.10.2016
![Page 6: Modeling and Simulation of the Thermal Runaway of ... · U(t), I(t), SOC T(t), Q(t) L. Cai and R. E. White, Journal of Power Sources 196 (2011) 5985–5989 Thermal conductivity Heat](https://reader035.vdocument.in/reader035/viewer/2022081404/5f064cfe7e708231d4174dd8/html5/thumbnails/6.jpg)
6 KIT, IAM-AWP
Heat transport equation Initial and boundary conditions
Thermal modeling based on partial differential
equations
𝜌: density of the cell cp: heat capacity
𝜅: thermal conductivity Tenv: environmental temperature
T0: initial temperature profile inside the cell at t = 0s
n: outward pointing normal vector h: heat transfer coefficient
𝜀: emissivity 𝜎: Stefan-Boltzmann constant
Dr. C. Ziebert – COMSOL Conference 2016 Munich, 12.10.-14.10.2016
![Page 7: Modeling and Simulation of the Thermal Runaway of ... · U(t), I(t), SOC T(t), Q(t) L. Cai and R. E. White, Journal of Power Sources 196 (2011) 5985–5989 Thermal conductivity Heat](https://reader035.vdocument.in/reader035/viewer/2022081404/5f064cfe7e708231d4174dd8/html5/thumbnails/7.jpg)
7 KIT, IAM-AWP
Constant fuel model
Newman model
Modeling the heat sources of a Li-ion cell
• Electrochemical contributions
• Exothermic contributions
Electrochemical heat source:
Electrochemical heat sources from Newman model
I: applied current U: cell voltage Ueq: equilibrium voltage of the cell
Dr. C. Ziebert – COMSOL Conference 2016 Munich, 12.10.-14.10.2016
• Reversible heat:
• Irreversible heat:
![Page 8: Modeling and Simulation of the Thermal Runaway of ... · U(t), I(t), SOC T(t), Q(t) L. Cai and R. E. White, Journal of Power Sources 196 (2011) 5985–5989 Thermal conductivity Heat](https://reader035.vdocument.in/reader035/viewer/2022081404/5f064cfe7e708231d4174dd8/html5/thumbnails/8.jpg)
8 KIT, IAM-AWP
Exothermic heat sources for thermal abuse
Simplification:
Constant Fuel Model
𝑄𝑖 𝑥, 𝑡 = 𝑞𝑖𝑅𝑖 𝑥, 𝑡
𝑅𝑖 𝑥, 𝑡 = 𝐴𝑖𝑐𝑖,0𝑒𝑥𝑝 −𝐸𝑎,𝑖
𝑅𝑇 𝑥, 𝑡
i ∈ {sei, pe, ne,ele}
ci: concentration of Li-ions
qi: reaction enthalpy in Jg-1
Ri: reaction rate in 1/s
Ai: frequency factor in 1/s
ci,0: initial concentration
Ea,i: activation energy in Jmol-1
R: universal gas constant
Dr. C. Ziebert – COMSOL Conference 2016 Munich, 12.10.-14.10.2016
![Page 9: Modeling and Simulation of the Thermal Runaway of ... · U(t), I(t), SOC T(t), Q(t) L. Cai and R. E. White, Journal of Power Sources 196 (2011) 5985–5989 Thermal conductivity Heat](https://reader035.vdocument.in/reader035/viewer/2022081404/5f064cfe7e708231d4174dd8/html5/thumbnails/9.jpg)
9 KIT, IAM-AWP
Simulation of a single 18650 cell with LiCoO2 cathode in
Battery and Fuel Cell Module of COMSOL Multiphysics 5.2
MeshingGeometry
Activ
eB
atte
ryM
ate
rial
Man
dre
l
Batte
ryC
asin
gAdaptive Triangulation
Quadratic Basis
Functions
Simulated Cell:
• Geometry: cylindrical 18650 cell
• Cathode: LiCoO2
• Anode: LixC6
• Electrolyte: 1:1 EC : DEC with LiPF6 salt
Source: Melcher, C. Ziebert, M. Rohde, B. Lei, H.J. Seifert, Modeling and
Simulation of the Thermal Runaway Behavior of Cylindrical Li-Ion Cells —
Computing of Critical Parameters, Energies 9 (2016) 292, doi:10.3390/en9040292.
What can be simulated / evaluated:
• Temperature field and gradients in time and space
• Mean temperature of cell, surface, single points
• Temperature profile along the coordinate axis
• Cell voltage and concentration of Li-ions
Dr. C. Ziebert – COMSOL Conference 2016 Munich, 12.10.-14.10.2016
![Page 10: Modeling and Simulation of the Thermal Runaway of ... · U(t), I(t), SOC T(t), Q(t) L. Cai and R. E. White, Journal of Power Sources 196 (2011) 5985–5989 Thermal conductivity Heat](https://reader035.vdocument.in/reader035/viewer/2022081404/5f064cfe7e708231d4174dd8/html5/thumbnails/10.jpg)
KIT, IAM-AWP10
1d – electrochemical 2d – thermal (axial-symmetric)
Diffusion of Li-ions in the electrodes
Diffusion of Li-ions in the electrolyte
Ohmic Losses
U(t), I(t), SOC
T(t), Q(t)
L. Cai and R. E. White, Journal of Power Sources 196 (2011) 5985–5989
Thermal conductivity
Heat capacity
Heat transport – Temperature distribution
Kim et.al, Journal of The Electrochemical Society, 158 (2011), A955-A969
Multi-scale multi-domain top-down approach of the Multi-Scale Multi-Dimensional (MSMD) model
Micro- Meso- Macroscale
COMSOL Multiphysics®
The Multi-Scale Multi-Dimensional (MSMD) model
Dr. C. Ziebert – COMSOL Conference 2016 Munich, 12.10.-14.10.2016
![Page 11: Modeling and Simulation of the Thermal Runaway of ... · U(t), I(t), SOC T(t), Q(t) L. Cai and R. E. White, Journal of Power Sources 196 (2011) 5985–5989 Thermal conductivity Heat](https://reader035.vdocument.in/reader035/viewer/2022081404/5f064cfe7e708231d4174dd8/html5/thumbnails/11.jpg)
11 KIT, IAM-AWP
Simulation of an electric load under adiabatic conditions
Dr. C. Ziebert – COMSOL Conference 2016 Munich, 12.10.-14.10.2016
Cell temperature profileLoad profile
![Page 12: Modeling and Simulation of the Thermal Runaway of ... · U(t), I(t), SOC T(t), Q(t) L. Cai and R. E. White, Journal of Power Sources 196 (2011) 5985–5989 Thermal conductivity Heat](https://reader035.vdocument.in/reader035/viewer/2022081404/5f064cfe7e708231d4174dd8/html5/thumbnails/12.jpg)
12 KIT, IAM-AWP
Temperature profile radial
Dr. C. Ziebert – COMSOL Conference 2016 Munich, 12.10.-14.10.2016
![Page 13: Modeling and Simulation of the Thermal Runaway of ... · U(t), I(t), SOC T(t), Q(t) L. Cai and R. E. White, Journal of Power Sources 196 (2011) 5985–5989 Thermal conductivity Heat](https://reader035.vdocument.in/reader035/viewer/2022081404/5f064cfe7e708231d4174dd8/html5/thumbnails/13.jpg)
13 KIT, IAM-AWP
Contour lines and radial profilesCell operates in normal mode
Source: A. Melcher, C. Ziebert, B. Lei, W.J. Zhao, M. Rohde, H.J. Seifert, Modellierung und Simulation des thermischen
Runaways in zylindrischen Li-Ionen Batterien, in: D. Tikhomiriv, H.-Th. Mammen, Th. Pawletta, Hrsg. ARGESIM Report 51,
ASIM Mitteilung AM 158, S. 8-28, ARGESIM Verlag Wien, Hochschule Hamm-Lippstadt 2016, ISBN 978-3-901608-48-3.
Dr. C. Ziebert – COMSOL Conference 2016 Munich, 12.10.-14.10.2016
![Page 14: Modeling and Simulation of the Thermal Runaway of ... · U(t), I(t), SOC T(t), Q(t) L. Cai and R. E. White, Journal of Power Sources 196 (2011) 5985–5989 Thermal conductivity Heat](https://reader035.vdocument.in/reader035/viewer/2022081404/5f064cfe7e708231d4174dd8/html5/thumbnails/14.jpg)
14 KIT, IAM-AWP
Cell goes into thermal runaway
Contour lines and radial profiles
Dr. C. Ziebert – COMSOL Conference 2016 Munich, 12.10.-14.10.2016
![Page 15: Modeling and Simulation of the Thermal Runaway of ... · U(t), I(t), SOC T(t), Q(t) L. Cai and R. E. White, Journal of Power Sources 196 (2011) 5985–5989 Thermal conductivity Heat](https://reader035.vdocument.in/reader035/viewer/2022081404/5f064cfe7e708231d4174dd8/html5/thumbnails/15.jpg)
15 KIT, IAM-AWP
Classification of thermal runaway
Consideration of the phase space
Two-dimensional projections:
Dr. C. Ziebert – COMSOL Conference 2016 Munich, 12.10.-14.10.2016
![Page 16: Modeling and Simulation of the Thermal Runaway of ... · U(t), I(t), SOC T(t), Q(t) L. Cai and R. E. White, Journal of Power Sources 196 (2011) 5985–5989 Thermal conductivity Heat](https://reader035.vdocument.in/reader035/viewer/2022081404/5f064cfe7e708231d4174dd8/html5/thumbnails/16.jpg)
16 KIT, IAM-AWP
All over !!!!!
Cell ist dead
Warning:
T too big
dT/dt to big
Thermal Runaway starts
Thermal
Runaway
Classification of thermal runaway: T-dT/dt-Plane
T/[K]
dT/dt
[K/min]
10
100
400 470
All OK !!!!
Dr. C. Ziebert – COMSOL Conference 2016 Munich, 12.10.-14.10.2016
![Page 17: Modeling and Simulation of the Thermal Runaway of ... · U(t), I(t), SOC T(t), Q(t) L. Cai and R. E. White, Journal of Power Sources 196 (2011) 5985–5989 Thermal conductivity Heat](https://reader035.vdocument.in/reader035/viewer/2022081404/5f064cfe7e708231d4174dd8/html5/thumbnails/17.jpg)
KIT, IAM-AWP17
Modeling and Simulation of the Thermal Runaway of Cylindrical Li-Ion Cells
Summary and Outlook
Results:
Extension of the classical model allows to simulate thermal runaway for cycling under adiabatic conditions
Thermal runaway can be classified using the phase space
Improvement of Thermal Management and Safety
Outlook:
Comparison with experiments in battery calorimeters
Refinement of model parameters: e.g. solid fuel model, microstructure-based modelling of the electrodes instead of porous electrode theory
Additional effects: e.g. temperature dependent separator performance, venting
Transfer to different cell geometries and chemistries
Dr. C. Ziebert – COMSOL Conference 2016 Munich, 12.10.-14.10.2016
![Page 18: Modeling and Simulation of the Thermal Runaway of ... · U(t), I(t), SOC T(t), Q(t) L. Cai and R. E. White, Journal of Power Sources 196 (2011) 5985–5989 Thermal conductivity Heat](https://reader035.vdocument.in/reader035/viewer/2022081404/5f064cfe7e708231d4174dd8/html5/thumbnails/18.jpg)
18 KIT, IAM-AWP
Thank You
For Your AttentionThis R&D project is partially funded by the Federal Ministry for Education and
Research (BMBF) within the framework “IKT 2020 Research for Innovations”
under the grant 16N12515 and is supervised by the Project Management Agency
VDI | VDE | IT
Supervised by
Dr. C. Ziebert – COMSOL Conference 2016 Munich, 12.10.-14.10.2016
![Page 19: Modeling and Simulation of the Thermal Runaway of ... · U(t), I(t), SOC T(t), Q(t) L. Cai and R. E. White, Journal of Power Sources 196 (2011) 5985–5989 Thermal conductivity Heat](https://reader035.vdocument.in/reader035/viewer/2022081404/5f064cfe7e708231d4174dd8/html5/thumbnails/19.jpg)
19 KIT, IAM-AWP
Backup Slides
Dr. C. Ziebert – COMSOL Conference 2016 Munich, 12.10.-14.10.2016
![Page 20: Modeling and Simulation of the Thermal Runaway of ... · U(t), I(t), SOC T(t), Q(t) L. Cai and R. E. White, Journal of Power Sources 196 (2011) 5985–5989 Thermal conductivity Heat](https://reader035.vdocument.in/reader035/viewer/2022081404/5f064cfe7e708231d4174dd8/html5/thumbnails/20.jpg)
KIT, IAM-AWP20
Simulation parameters for electrochemical-thermal model
Source: Melcher, C. Ziebert, M. Rohde, B. Lei, H.J. Seifert, Modeling and Simulation of the Thermal Runaway Behavior of
Cylindrical Li-Ion Cells —Computing of Critical Parameters, Energies 9 (2016) 292, doi:10.3390/en9040292.
Dr. C. Ziebert – COMSOL Conference 2016 Munich, 12.10.-14.10.2016
![Page 21: Modeling and Simulation of the Thermal Runaway of ... · U(t), I(t), SOC T(t), Q(t) L. Cai and R. E. White, Journal of Power Sources 196 (2011) 5985–5989 Thermal conductivity Heat](https://reader035.vdocument.in/reader035/viewer/2022081404/5f064cfe7e708231d4174dd8/html5/thumbnails/21.jpg)
KIT, IAM-AWP21
Simulation parameters for exothermic model
Source: Melcher, C. Ziebert, M. Rohde, B. Lei, H.J. Seifert, Modeling and Simulation of the Thermal Runaway Behavior of
Cylindrical Li-Ion Cells —Computing of Critical Parameters, Energies 9 (2016) 292, doi:10.3390/en9040292.
Dr. C. Ziebert – COMSOL Conference 2016 Munich, 12.10.-14.10.2016
![Page 22: Modeling and Simulation of the Thermal Runaway of ... · U(t), I(t), SOC T(t), Q(t) L. Cai and R. E. White, Journal of Power Sources 196 (2011) 5985–5989 Thermal conductivity Heat](https://reader035.vdocument.in/reader035/viewer/2022081404/5f064cfe7e708231d4174dd8/html5/thumbnails/22.jpg)
KIT, IAM-AWP22
Battery
Monitoring
Battery
Management
Software
Battery
Monitoring
Automotive uC
IKEBA Project
Integrated Components and Integrated Design of Energy Efficient Battery Systems
5 cooperating partners Duration: 05/2013 – 04/2016 Budget: 7 Million Euro
Dr. C. Ziebert – COMSOL Conference 2016 Munich, 12.10.-14.10.2016