2.3. Battery Thermal Runaway Modeling Method
- (1)
Modeling of side-effect mechanisms
The battery temperature continues rising and triggers the self-exothermic side reactions, accelerating the evolution of the thermal runaway. This paper mainly considers the four main side reactions: the SEI film decomposition reaction, the negative electrode decomposition reaction, the positive electrode decomposition reaction, and the electrolyte oxidation reaction. Moreover, the exothermic rate of each reaction is described according to Arrhenius’ formula. The cumulative heat generation rates induced by the side reactions are computed using the following equation:
where
,
,
and
correspond to the heat generation rate per unit volume of the four side reactions mentioned above, in W/m
−3.
SEI film is formed during the first cycle of the battery to avoid the negative electrode directly reacting with the electrolyte, and it is easy to decompose and exothermic at 80~120 °C [
12]. The following equation can describe the decomposition reaction process and its heat generation rate:
where
can represent sei (SEI film), ne (negative reaction), pe (positive reaction), and e (electrolyte reaction);
is the temperature of the cell in K;
is a dimensionless quantity, which represents the concentration of the reactants;
,
,
, and
represent the frequency factor, the activation energy of the reaction, the relationship between the reaction rate and the concentration of the reactants, and the correction term of the reaction, respectively;
denotes the exothermic heat of the reaction per unit of weight in J/kg;
is the mass per unit volume of reactants in kg/m; and R is the ideal gas constant.
Then, the electrolyte reacts with the negative electrode when the temperature reaches around 120 °C, generating heat of
[
13]. Then, the positive electrode decomposes partially at around 160 °C with the heat of
[
14]. The electrolyte can decompose exothermically at elevated temperatures (>200 °C), as
[
15].
- (2)
Thermal runaway trigger modeling
The thermal runaway side-reaction heat generation described above occurs significantly only when the battery temperature surpasses 80~100 °C. Consequently, it becomes imperative to model the catalytic heat, , associated with thermal runaway before the temperature rises from ambient conditions to initiate the side reaction.
The thermal runaway of batteries can be classified into three types of misuse based on their triggering mechanisms: mechanical abuse (e.g., needle penetration and compression), electrical abuse (including overcharge, internal short circuit, and external short circuit), and thermal abuse (such as exposure to high temperatures). These triggering mechanisms are interconnected and can influence one another. For instance, in the case of mechanical abuse, deformation of the battery may result in the rupture of internal separators. This deformation can lead to contact between positive and negative electrodes, causing an internal short circuit and electrical abuse. The progression of an internal short-circuit generates abnormal heat, leading to a rapid temperature rise and eventual thermal abuse. The heightened temperature induced by various forms of misuse can ultimately trigger thermal runaway.
In the context of TRP within a battery module, cells undergoing thermal runaway generate a significant amount of abnormal heat. This abnormal heat is predominantly transferred to adjacent cells through thermal conduction, heating normal cells and subsequently triggering thermal runaway. In this study, the analysis of battery thermal propagation behavior categorizes triggering events leading to thermal runaway into two modes based on the abnormal heat generation rate of cells in thermal runaway: the slow heat-generation mode associated with fault evolution and the rapid heat-generation mode associated with sudden abuse.
- a.
Fault evolution thermal runaway-trigger modeling method
The abnormal heat production rate of the battery module is comparatively low, leading to a gradual initiation of battery thermal runaway, characterized by an extended duration. The thermal runaway is instigated by heat sources, encompassing internal chemical reactions within the battery and the heat generated during the charging and discharging cycles.
where
is the chemical reaction heat generation in the fault evolution, which is related to the battery temperature and fitted via the experiment; and
is the heat generation in the charging and discharging cycle.
Take overcharge abuse as an example. The chemical reaction occurring within the battery is the reaction between plated lithium and electrolyte induced by overcharging. The rate of reaction can be mathematically expressed as follows (6) [
16]. The reaction rate is represented by the following:
where
is the amount of the deposited lithium on the anode,
is the frequency factor of the reaction,
is the activation energy of the reaction,
is the reaction heat per unit mass of lithium plating,
is the reaction mass of lithium plating, and
determines whether or not the reaction occurs between lithium and electrolyte. When
>> 1,
is equal to 1, the reaction occurs; when
= 0,
is equal to 0, then the reaction does not occur [
13].
Bernardi’s heat generation rate equation can calculate the general case of heat production in charge/discharge cycles. Considering that the overcharge process is mainly battery ohmic heat generation, it can be simplified by (7):
where
is the charging current, and
is the DC internal resistance obtained experimentally.
- b.
Sudden abuse thermal runaway-trigger modeling method
In the case of thermal runaway triggered by sudden abuse, the battery module exhibits a relatively high abnormal heat production rate, resulting in a rapid initiation of battery thermal runaway. The time required for thermal runaway triggered by sudden abuse is short. The thermal runaway trigger in this context can be simplified as follows:
where
is the heat generation of sudden abuse, time-dependent, fitted by experiment.
Therefore, the total heat generated by the thermal runaway of the battery is as follows:
where
represents the total additional heat.