基于电化学-热耦合模型的锂离子电池组件产热分析

doi:10.19799/j.cnki.2095-4239.2023.0082

基于电化学-热耦合模型的锂离子电池组件产热分析

杨佳兴^,, 张恒运^,, 徐屹东

上海工程技术大学机械与汽车工程学院，上海 201620

Heat generation analysis for lithium-ion battery components using electrochemical and thermal coupled model

YANG Jiaxing^,, ZHANG Hengyun^,, XU Yidong

School of Mechanical & Automotive Engineering, Shanghai University of Engineering Science, Shanghai 201620, China

通讯作者: 张恒运，教授，研究方向为动力电池热管理，E-mail：zhanghengyun@sues.edu.cn。

收稿日期: 2023-02-20 修回日期: 2023-03-14

基金资助:

国家自然科学基金. 51876113
上海市自然科学基金项目. 21ZR1426300

Received: 2023-02-20 Revised: 2023-03-14

作者简介 About authors

杨佳兴（1996—），男，硕士研究生，研究方向为动力电池热管理，E-mail：yoogating@163.com； E-mail：yoogating@163.com。

摘要

研究电池电化学过程产热对锂离子电池的热管理至关重要。本工作建立了三元NMC锂离子电池的电化学-热耦合模型，首先通过对该电池进行不同倍率的放电与温度实验测试，验证了该模型在电压和温度变化预测准确性。然后针对不同温度下的表现进行模拟仿真研究。在室温下，无论倍率大小，负极产热总是小于正极产热，虽然负极的极化热高于正极，但其可逆吸热较大，导致产热水平低于正极。而随着放电倍率的增加，正极产热所占比例减小，负极所占比例先增加后减小，而集流体产热所占比例持续增加。然而，低温条件下的电池放电表现出与室温情况不同的产热特性，首先，低温导致低倍率负极产热率比例大大增加，负极可逆热为总可逆热的主要贡献热。而高倍率负极产热率减少，正极则呈相反趋势。其次在低温下放电时间随倍率增加呈现不同趋势，高倍率下放电电压快速降低导致放电不完全，在低倍率0.5～1 C放电运行时出现了电压反弹现象但基本放电完全，这是由于低温限制了负极颗粒内部锂离子及时向外扩散，造成电阻增加与电压快速降低，同时大量产热导致自身温升，从而在低倍率下获得电压反弹并保持持续放电的能力。

关键词： 锂离子电池 ; 电化学热耦合模型 ; 产热 ; 低温特性 ; 电压反弹

Abstract

Understanding the thermal aspects of the battery electrochemical process is essential for managing the heat generated in lithium-ion batteries. Herein, an electrochemical-thermal coupled model was developed for an NMC lithium-ion battery. Experimental tests were conducted to validate the accuracy of the model in predicting voltage and temperature variations at different discharge rates and temperatures. Simulation studies were then conducted to estimate battery performance at different temperatures. At normal temperatures, regardless of the discharge rate, the NE heat generation was always smaller than the PE heat generation. While the NE polarization heat was higher than the PE polarization heat, the reversible heat absorption of NE was larger than that of PE. Consequently, the NE heat generation remained lower than that of PE. As the discharge rate increased, the proportion of PE heat generation decreased, whereas the proportion of NE heat generation increased first and then decreased. Simultaneously, the proportion of collector heat generation continued to rise. However, the heat generation characteristics differed when the battery discharge was discharged at subzero temperature compared to the normal temperature case. Firstly, subzero temperatures considerably increased the proportion of NE generation at low discharge rates. The reversible heat from NE became the primary contributor to the total reversible heat. Conversely, at high discharge rates, the NE heat generation rate decreased, while the PE heat generation rate exhibited the opposite trend. Secondly, the discharge time at subzero temperatures showed different trends with increasing discharge rates. At high discharge rates, the voltage rapidly decreased, leading to incomplete discharge and a phenomenon known as voltage rebound. At low discharge rates of 0.5—1 C discharge operation, the basic discharge process was completed, but the voltage rebound still occurred. The subzero temperature restricted the timely diffusion of lithium ions from the NE particles, leading to increased resistance and a rapid voltage drop. Meanwhile, a large amount of heat generation contributed to the temperature rise in the battery. Therefore, the battery exhibited the ability to rebound and maintain a continuous discharge at low rates.

Keywords： lithium-ion battery ; electrochemical-thermal coupled model ; heat generation ; subzero temperature characteristics ; voltage rebound

PDF (4212KB) 元数据多维度评价相关文章导出 EndNote| Ris| Bibtex 收藏本文

本文引用格式

杨佳兴, 张恒运, 徐屹东. 基于电化学-热耦合模型的锂离子电池组件产热分析[J]. 储能科学与技术, 2023, 12(8): 2615-2625

YANG Jiaxing. Heat generation analysis for lithium-ion battery components using electrochemical and thermal coupled model[J]. Energy Storage Science and Technology, 2023, 12(8): 2615-2625

锂离子电池由于其高能量密度、快速充电和放电能力、无记忆效应和低自放电率而被认为是最受欢迎的动力源。然而，锂电池的性能与工作条件和环境温度密切相关。温度过低或者过高都对电池的使用寿命和内部组成产生不利影响^[1]。

锂离子电池是一个典型的电化学储能装置，包含了各种电化学反应、离子输运以及热传递过程。Bernardi等^[2]和Thomas等^[3]提出了电化学系统及多孔插入电极的通用能量平衡方程，包括不可逆热源项、可逆热源项以及相变和热容量变化项。国内外学者也基于不同产热模型对锂离子电池的电热特性进行了分析^[4-5]，但这些模型在使用过程有所简化，导致其内部各组件产热特性无法评估。Newman等^[6]提出了基于多孔电极理论的准二维模型(P2D模型)，考虑了液态电解液中的浓溶液和固态电极的传输动力学，并利用Fick扩散定律来描述固相材料中锂离子的插入和脱离，自此该模型得到了广泛应用^[7-11]。Lyu等^[8]利用P2D模型研究了高镍NMC储能锂离子电池的热特性。Du等^[9]建立了一个电化学-热耦合模型来考虑锂离子电池内部不可逆热量的产生和演变。张志超等^[10]利用P2D模型研究了层叠式锂离子电池在多尺度的电化学与热行为。

需要指出的是，目前的研究大都按照电池各部分热量产生机制进行求解。而在实际应用中，基于电池内部分层组件产热进行区分更加合理，故需要对电池分层组件产热机制进行研究。因此本工作旨在研究NMC锂离子电池分层组件在不同温度、不同倍率下的产热特性。为此，基于P2D模型分析微观层面的电化学浓度和电荷传输并获得电池各层组件的产热变化规律。研究了不同放电倍率下常温与低温产热特性。其中正、负极的极化热，正、负极的可逆热，以及电解液、颗粒的欧姆热，都是通过数值求解确定的。

1 电化学-热耦合模型

选用额定容量为2.7 Ah、以三元材料为正极材料的N18650CL-29圆柱形锂离子电池进行模型分析和实验测试，其基本参数见表1。电化学-热耦合模型的流程如图1所示，它可以分解为以下几个步骤：①拆解真实电池，确定集流体以及电极尺寸，通过建立集流体的等效电阻计算集流体产热量。②输入电极模型的参数，并根据电极模型计算电极产热。③用电池的总产热耦合进三维热模型中计算温升，并将电池温度再度耦合至电极模型中重新计算产热率，实现双向耦合直至迭代收敛。

表1 实验中使用的电池的基本参数

Table 1 Basic parameters of the battery used in the experiment

参数	数值
标准容量/mAh	2700
标称电压/V	3.60
充电截止电压/V	4.20
放电截止电压/V	2.50

新窗口打开| 下载CSV

图1

新窗口打开| 下载原图ZIP| 生成PPT

图1 电化学-热耦合流程

Fig. 1 Electrochemical thermal coupled flow diagram

图2是包含负极、隔膜和正极的锂离子电池一维模型，在宏观方向上沿电极厚度方向，在微观方向沿活性材料的固体颗粒的半径方向。

图2

新窗口打开| 下载原图ZIP| 生成PPT

图2 电极和颗粒尺度的电化学模型

Fig. 2 Schematic diagram of the electrochemical model in both electrode scale and particle scale

1.1　电化学模型

负极由人造石墨制成，电解液为LiPF₆溶液(1∶1∶1 EC/EMC/DMC)，正极为Li[Ni_0.5Mn_0.3Co_0.2]O₂。P2D模型用表2中的5个控制方程表示，首先对于正负极电极颗粒而言，在电极颗粒内部的物质传输分别用Fick第二定律与欧姆定律表示，Fick第二定律将固相球形颗粒径向上的任意一点锂离子浓度随时间t的变化率转化为方程右侧中该点锂离子浓度随空间的变化率。而用欧姆定律描述了电子在固相颗粒沿电极厚度方向的移动。在电极颗粒内部的电荷守恒用式(3)表示。在电解液中的输运通常用浓溶液理论修正的Nernst-Planck方程来描述，而基于浓溶液理论，给出液相中锂离子的物质守恒方程。等效参数由下式进行描述：

表2 P2D模型控制方程

Table 2 P2D model control equations

位置		控制方程		边界条件
正负电极颗粒	物质传输	$\frac{\partial c_{s}}{\partial t} = \frac{1}{r^{2}} \frac{\partial}{\partial r} (D_{s} r^{2} \frac{\partial c_{s}}{\partial r})$	(1)	$\frac{\partial c_{s}}{\partial t} \|_{r = 0} = 0$ , $- F D_{s} \frac{\partial c_{s}}{\partial r} \|_{r = r_{p}} = j_{P}$ , $c_{s} \|_{t = 0} = c_{s, 0}$
	物质传输	$i_{s} = - σ_{e f f} \nabla φ_{s}$	(2)	$i_{s} \|_{x = 0} = i_{s} \|_{x = l_{n e} + l_{S} + l_{p e}} = I$ , $i_{s} \|_{x = l_{n e}} = i_{s} \|_{x = l_{n e} + l_{S}} = 0$
	电荷守恒	$i_{l} + i_{s} = I$ ， $\nabla \cdot i_{l} + \nabla \cdot i_{s} = 0$	(3)	$\nabla \cdot i_{l} = α_{P} j_{P}$ ， $\nabla \cdot i_{s} = - α_{P} j_{P}$
电解液	物质传输	$i_{l} = - κ_{e f f} \nabla φ_{l} + \frac{2 κ_{e f f} R T}{F} (1 - t_{+}) (1 + \frac{d l n f_{\pm}}{d l n c_{l}}) \nabla l n c_{l}$	(4)	$i_{l} \|_{x = 0} = i_{l} \|_{x = l_{n e} + l_{S} + l_{p e}} = 0$
电解液	物质守恒	$ε \frac{\partial c_{l}}{\partial t} = \nabla \cdot (D_{e f f} \nabla c_{l}) + (1 - t_{+}) \frac{\nabla \cdot i_{l}}{F}$	(5)	$i_{l} \|_{x = l_{n e}} = i_{l} \|_{x = l_{n e} + l_{S}} = I$

新窗口打开| 下载CSV

\{\begin{array}{l} σ_{e f f} = σ {ε_{i}}^{B r u g}, i = n e, p e \\ κ_{e f f} = κ {ε_{l}}^{B r u g} \\ D_{e f f} = D_{l} {ε_{l}}^{B r u g} \end{array}

(6)

式中， σ 为固相电导率；κ为液相电导率；D_l为锂离子液相扩散系数；Brug为Bruggeman系数。

电化学反应是锂离子电池产生持续充放电能力的基础，利用Buter-Volmer方程来描述电池内部的电极反应动力学。

j_{P} = j_{0} [e x p (\frac{α_{a} F}{R T} η) - e x p (- \frac{α_{c} F}{R T} η)]

(7)

式中，j₀为交换电流密度； α_a和 α_c为传递系数； η 为过电位。B-V方程将反应电流密度与过电位联系起来。

j_{0} = F k {(\frac{c_{l}}{c_{l, r e f}})}^{α_{a}} \cdot c_{s, s u r f}^{α_{c}} {(c_{s, m a x} - c_{s, s u r f})}^{α_{a}}

(8)

式中，k为电极界面反应速率常数；c_{l, ref}为液相锂离子参考浓度；c_{s, surf}为固相颗粒表面锂离子浓度；c_{s, max}为固相颗粒最大锂离子浓度。过电位定义为 η = φ_s- φ_l-E。 φ_s为固相电势， φ_l为液相电势，E为电极固相材料的平衡电势，与电极SOC和温度(T)有关，相应的方程为

E = E_{r e f} + \frac{\partial E}{\partial T} (T - T_{r e f})

(9)

全电池的熵系数通过测量不同温度下的开路电压得到。采用文献[12]中负极熵系数，并利用实测的全电池熵系数计算得到正极熵系数^[13-14]，如图3所示。

图3

新窗口打开| 下载原图ZIP| 生成PPT

图3 正负电极熵系数随锂离子化学计量的变化

Fig. 3 NE and PE entropy coefficients with lithium-ion stoichiometry

这里电极SOC为锂离子的化学计量浓度，其公式为：

x, y = \frac{c_{s}}{c_{s, m a x}}

(10)

1.2　产热模型

锂离子电池的产热率包括来自负极和正极、集流体与电解液的产热

q_{t o t} = q_{n e} + q_{p e} + q_{C} + q_{e l}

(11)

其中单个正、负电极包含3个主要的热源项，即极化热(q_pol)、欧姆热(q_Ohm)和可逆热(q_rev)，由式(12)～(14)表示，集流体中的热量由式(15)的欧姆热表示，电解液欧姆热由式(16)计算。产热公式汇总在表3中。

表3 产热方程

Table 3 Heat generation equation

位置	方程
正负电极颗粒	$q_{p o l} = α_{P} j_{P} η$	(12)
	$q_{O h m} = σ_{e f f} \nabla φ_{s} \cdot \nabla φ_{s} + κ_{e f f} \nabla φ_{l} \cdot \nabla φ_{l} + \frac{2 κ_{e f f} R T}{F} (t_{+} - 1) (1 + \frac{d l n f_{\pm}}{d l n c_{l}}) \nabla l n c_{l} \cdot \nabla φ_{l}$	(13)
	$q_{r e v} = α_{P} j_{P} T \frac{\partial E}{\partial T}$	(14)
集流体	$q_{C} = (R_{C, C u} + R_{C, A l}) I^{2}$	(15)
电解液	$q_{e l} = κ_{e f f} \nabla φ_{l} \cdot \nabla φ_{l} + \frac{2 κ_{e f f} R T}{F} (t_{+} - 1) (1 + \frac{d l n f_{\pm}}{d l n c_{l}}) \nabla l n c_{l} \cdot \nabla φ_{l}$	(16)

新窗口打开| 下载CSV

1.3　模型参数

P2D模型参数包括正极和负极颗粒的平均半径、正极和负极的平衡电位、正极和负极的熵系数、电极活性材料的最大锂离子浓度、电极的固相扩散系数、电解液扩散系数、电极的固相电导率、电解液电导率等。本工作中NMC锂电池的参数来自于文献[8, 15-22]，室温下电化学参数见表4，变温动态参数见表5。

表4 室温下电化学模型参数

Table 4 Electrochemical model parameters at normal temperature

参数	负极	隔膜	正极
厚度/μm	72.5	32	69.5
颗粒半径r/μm	7		5
固相体积分数 ε_s	0.5386		0.6348
电解液体积分数 ε_l	0.2807	0.43	0.3061
固相最大锂离子浓度c_{s, max}/(mol/m³)	28607		49388
固相初始锂离子浓度c_{s, 0}/(mol/m³)	0.84c_{s, max, ne}		0.279c_{s, max, pe}
电解液参考锂离子浓度c_{l, 0}=c_ref/(mol/m³)	1000	1000	1000
传递系数 α_ne, α_pe	0.5		0.5
Bruggeman因子Brug	1.5	1.5	1.5
锂离子固相扩散活化能E_aR/(J/mol)	68025.7		9976.8
锂离子固相参考扩散系数D_{s, ref}/(m²/s)	1.4523 $\times$ 10^-13		5 $\times$ 10^-13
等效固相电导率 σ /(S/m)	100		3.8
转移数t₊	0.363
电池径向热导率 λ_xy /[W/(m·K)]	1.63
电池轴向热导率 λ_z /[W/(m·K)]	36.96
电池密度 ρ /(kg/m³)	2510
比定压热容c_p /[J/(kg·K)]	1028

新窗口打开| 下载CSV

表5 变温电化学动态参数

Table 5 Variable temperature electrochemical dynamic parameters

参数	动态表达式
等效电解液电导率 к /(S/m)	$κ (c_{r e f}, T) = 10^{- 4} c_{r e f} {(\begin{array}{l} - 8.2488 + 0.053248 T - 0.00002987 T^{2} + \\ 0.26235 \times 10^{- 3} c_{r e f} - 0.0093063 \times 10^{- 3} c_{r e f} T + \\ 0.000008069 \times 10^{- 3} c_{r e f} T^{2} + \\ 0.22002 \times 10^{- 6} c_{r e f}^{2} - 0.0001765 \times 10^{- 6} c_{r e f}^{2} T \end{array})}^{2}$
锂离子固相扩散系数D_s/(m²/s)	$D_{s} = {D_{s}}_{, r e f} e x p [\frac{E_{a R}}{R} (\frac{1}{T_{r e f}} - \frac{1}{T})]$
锂离子液相扩散系数D_l/(m²/s)	$l g [D_{l} (c_{r e f}, T)] = {(- 4.43 - \frac{54}{T - 229 - 0.005 c_{r e f}} - 0.00022 c_{r e f})}^{- 4}$
活性系数ν(f_±, t₊)	$ν (c_{r e f}, T) = \frac{1 - 0.363}{1 - 0.399} \{\begin{array}{l} 0.601 - 0.24 \times 10^{- 2} c_{r e f}^{0.5} + \\ 0.982 \cdot [1 - 0.0052 (T - 293.15)] \cdot (10^{- 3} c_{r e f})^{1.5} \end{array}\}$

新窗口打开| 下载CSV

1.4　实验测试验证

图4(a)展示了电池热测试实验图。使用K型热电偶和温度采集装置(Agilent 34972A)监测并记录电池温度，时间间隔为3 s。为了更准确地测量电池的放电温度，电池表面去除绝缘膜并焊接热电偶，直接测量电池的金属表面温度。在验证研究中，在室温25 ℃，电池外表面的对流传热系数确定为13.4 W/(m²·K)^[23]。同时记录电压和电流。根据电池内部集流体结构构建集流体等效电路并得到负极集流体电阻为9.19 mΩ，正极集流体电阻为3.4 mΩ，该部分产热不能忽略，需要考虑在内。

图4

新窗口打开| 下载原图ZIP| 生成PPT

图4 (a) 18650锂离子热测试电池；(b) 基于极耳与集流体关系的负极集流体ECM；(c) 基于极耳与集流体关系的正极集流体ECM

Fig. 4 Schematic of (a) a 18650 Li-ion thermal test battery; (b) NE collector ECM based on lug-collector relationship; (c) PE collector ECM based on lug-collector relationship

图5为室温(25 ℃)下不同放电倍率0.5 C、1 C和3 C下仿真与实验的对比。其中图5(a)为电压对比而图5(b)是温升对比，可以看出该模型在电压和温度上都有很好的一致性。

图5

新窗口打开| 下载原图ZIP| 生成PPT

图5 (a) 显示0.5 C、1 C和3 C的放电电压与时间的关系；(b) 显示0.5 C、1 C和3 C的放电温度与时间的关系。

Fig. 5 (a) shows 0.5 C, 1 C and 3 C discharge voltage versus time; (b) shows 0.5 C, 1 C and 3 C discharge temperature versus time

2 结果与讨论

2.1　室温放电的产热分析

本节分析在室温和不同放电倍率下(0.5 C、1 C、2 C、3 C)电池内各组件的产热，不同倍率下的电池平均总产热符合文献[24]产热与倍率的拟合关系，并着重对0.5 C放电倍率下各组件的瞬时产热进行了分析对比。图6显示了室温下不同放电倍率的负极、电解液、集流体和正极的产热率。由图6(a)可以看出，尽管负极的极化热高于正极，但由于负极的可逆产热存在吸热阶段，因此负极产热量低，而正极产热率大于负极，容易局部升温导致安全问题。而随着放电倍率的增加，正极、负极、电解液与集流体的产热均有所增加。由表6可以看出，随着放电倍率的增加，负极产热比例先增加后减小，正极产热的比例在减小而集流体产热的比例在增加，这是由于集流体欧姆产热与电流平方成正比，导致产热比例快速上升。而从低温下的电池各组件产热比例可以看出，随着放电倍率的增加，正极和负极占比在减小，而集流体产热的比例在增加。

图6

新窗口打开| 下载原图ZIP| 生成PPT

图6 室温下不同放电倍率电池内部组件产热率

Fig. 6 Heat comparison of different heat generation components at 25 ℃

表6 室温与低温下不同放电倍率各组件产热比例

Table 6 Heat production ratio of each component with different discharge rates at normal and subzero temperature

电池组件	室温				低温
电池组件	0.5 C	1 C	2 C	3 C	0.5 C	1 C	2 C	3 C
负极	19.56%	29.02%	26.3%	23.73%	69.34%	65.13%	36.58%	29.96%
正极	59.26%	45.96%	34.08%	28.36%	20.98%	20.81%	38.73%	39.85%
电解液	0.95%	0.93%	0.52%	0.58%	0.27%	0.92%	0.78%	0.96%
集流体	20.23%	24.08%	39.11%	47.33%	9.41%	13.13%	23.91%	29.23%

新窗口打开| 下载CSV

2.1.1　不可逆热分析

图7显示了室温0.5 C下负极、正极的极化热与过电位的比较。可以看出，负极的极化热总是比正极的极化热大。放电结束时，负极的极化热快速增长，而正极的极化热几乎保持不变。这主要是因为在接近放电结束时，大量的锂离子在负极被消耗掉了，导致负极内部严重缺锂的极化状态。负极过电位快速增加作为对比，正极过电位基本不变。

图7

新窗口打开| 下载原图ZIP| 生成PPT

图7 0.5 C下25 ℃时，正负极极化热与过电位随时间变化的曲线

Fig. 7 Plot of polarization heat and overpotential of PE and NE with time at 0.5 C at 25 ℃

欧姆热分为集流体欧姆热、负极和正极颗粒欧姆热以及电解液欧姆热，其中集流体欧姆热是主要贡献热。在电极内部，电解液欧姆热是主要的贡献热。如图8所示，电解液欧姆热大于颗粒欧姆热。随着放电时间的增加，电解液中的热量缓慢减少，固相材料的欧姆热保持不变。正极的欧姆热大于负极的欧姆热，两者受温度的影响较小。

图8

新窗口打开| 下载原图ZIP| 生成PPT

图8 在0.5 C下25 ℃时，集流体欧姆热与不同组件欧姆热随时间变化的曲线

Fig. 8 Collector ohmic heat versus different component ohmic heat at 0.5 C at 25 ℃

2.1.2　可逆热分析

图9显示了正极与负极在室温与低温时的可逆热比较。在室温下正极的可逆热全程放热，而负极在放电后期才转为放热，而在低温下，0.5 C倍率的负极可逆放热在放电中期超过了正极可逆热，总可逆热的主要贡献来自负极可逆放热。

图9

新窗口打开| 下载原图ZIP| 生成PPT

图9 正负极在0.5 C放电时的可逆热比较

Fig. 9 Comparison of reversible heat for PE and NE

2.2　低温放电的产热分析

低温限制了锂离子在液体中的移动和在固体中的扩散，使得锂离子在电池内的传输动力不足，导致内阻和产热增加。另外，通常在低温下发现的电压反弹现象^[25]与内部电化学状态的关系也将在本节阐述。

2.2.1　电池内部热源分析

图10显示了低温下(-15 ℃)不同放电倍率的负极、电解液、集流体和正极产热率。与室温情况相比，低温低倍率(0.5～1 C)下负极的产热率明显大于正极，而在2～3 C高倍率，正极产热较负极大。图10(a)为0.5 C放电倍率下的产热，其发生在负极的极化热和可逆热达到了69.34%，而正极的产热量下降到20.98%，其余9.68%的热量发生在电解液和集流体的欧姆热。随着放电倍率的增加，尽管负极极化热一直贡献了最大产热，但在3 C下，集流体产热超过负极极化热。

图10

新窗口打开| 下载原图ZIP| 生成PPT

图10 -15 ℃下不同放电倍率电池内部组件产热

Fig. 10 Comparison of different heat generation components at -15 ℃

图11展示了不同放电倍率下的正负极表面SOC(锂离子的化学计量浓度)的变化。其中小倍率放电时间较长，并且正极表面SOC变化较大，说明此种倍率下放电是完全的。而在2 C与3 C高倍率时，正极表面SOC升高很小，在此倍率下电池放电不完全。

图11

新窗口打开| 下载原图ZIP| 生成PPT

图11 在低温-15 ℃下不同放电倍率正负极表面SOC

Fig. 11 Surface SOC of NE and PE at different discharge rates at subzero temperature (-15 ℃)

2.2.2　0.5 C放电机理分析

图12显示了不同组件在-15 ℃下的产热与电压的关系，以-15 ℃的0.5 C为例。放电3575 s为电压反弹点，同时也对应于最大产热率。该温度下负极的总产热率高于正极，此时负极的极化热远高于正极的极化热，这与室温情况不同。这是因为低温增加了锂离子从负极颗粒向外扩散的阻力，从而增加了产热率。

图12

新窗口打开| 下载原图ZIP| 生成PPT

图12 在0.5 C放电下-15 ℃时，每个组件的端电压与产热率的关系

Fig. 12 Terminal voltage versus heat generation for each component at 0.5 C and -15 ℃

图13显示了0.5 C在-15 ℃下，4个不同时间节点的负极颗粒锂离子浓度。在10 s时负极颗粒内外都有足够的锂离子，但当放电过程进行到3575 s时，负极颗粒表面锂离子浓度降至低点(199 mol/m³)，此刻对应产热曲线的最高点。由于电池产热增加温度升高，加速了锂离子从电极内部向外部的扩散速度，于是负极颗粒表面锂离子浓度先是增加，而在放电后期浓度耗尽(6405 s)，导致放电趋于终止。

图13

新窗口打开| 下载原图ZIP| 生成PPT

图13 在0.5 C放电下-15 ℃时锂离子在负极颗粒内的浓度随时间变化的关系

Fig. 13 Lithium-ion concentration within the NE particle as a function of time at 0.5 C at -15 ℃

3 结论

构建电化学-热耦合模型来研究18650圆柱形锂离子电池在不同工作条件下的产热情况，分析了电池各产热部分在室温和低温下的演变，并以0.5 C为例，分析了低温下放电电压反弹内在机制，归纳出以下结论。

（1）分析了室温下电池不同组件的各种产热情况。尽管负极的极化热高于正极，但可逆热吸热导致其产热率低于正极。而随着放电倍率的增加，负极产热在整个电池中所占比例先增大后减小，正极产热占比减小，而集流体产热与电流平方正相关，导致高倍率占比增加。负极极化热在放电末期快速增长的原因是由于缺锂导致负极过电位增加。

（2）在-15 ℃的低温产热与室温工况下有明显不同。首先其在低倍率下放电负极的产热量超过了正极，而在高倍率下负极产热低于正极的产热特性。在低倍率下放电时，负极的极化热与可逆热占比较高，而高倍率下放电，正极极化热与可逆热占比较高，而不论何种倍率，集流体产热总是随放电倍率增加。

（3）同时在-15 ℃电池0.5 C放电，会出现电压反弹，反弹时间与负极颗粒表面SOC引发的产热极大值和锂离子浓度极小值相对应。电池低锂态运行使得温度升高并加速了锂离子在负极颗粒中的向外运动，导致放电曲线中的电压反弹。

符号说明

c	锂离子浓度，mol/m³
*c_p*	定压比热容，J/(kg·K)
D	扩散系数，m²/s
E	电压，V
E_aR	扩散系数活化能，J/mol
F	法拉第常数，C/mol
h_c	对流换热系数，W/(m²·K)
i	局部电流密度，A/m²
j₀	交换电流密度，A/m²
k	反应速率常数，m^2.5/(mol^0.5·s)
L	厚度，m
q_C	集流体欧姆热，W/m³
q_el	电解液欧姆热，W/m³
q_tot	总产热率，W/m³
q_ne, q_pe	负极、正极产热率，W/m³
R	理想气体常数，J/(mol·K)
r_p	颗粒半径，μm
T	温度， ℃
t	时间，s
t₊	转移数
U	开路电压，V
ν	活性因子
希腊符号
α_n, α_p	负极或正极传递系数
ε	体积分数
к	液相电导率，S/m
η	过电位，V
λ	热导率，W/(m·K)
ρ	电池密度，kg/m³
σ	固相电导率，S/m
φ	电势，V
角标
Brug	Bruggeman因子
C	集流体
eff	有效值
el	电解液
max	最大值
min	最小值
ne	负极
pe	正极
ref	参考值
S	隔膜
s	固相
surf	表面
0	初始值

缩略词

ECM	等效电路模型
NE	负极
PE	正极

参考文献

原文顺序

文献年度倒序

文中引用次数倒序

被引期刊影响因子

[1]

FENG X N, OUYANG M G, LIU X, et al. Thermal runaway mechanism of lithium ion battery for electric vehicles: A review[J]. Energy Storage Materials, 2018, 10: 246-267.

[本文引用: 1]

[2]

BERNARDI D, PAWLIKOWSKI E, NEWMAN J. A general energy balance for battery systems[J]. Journal of the Electrochemical Society, 1985, 132(1): 5-12.

[本文引用: 1]

[3]

THOMAS K E, NEWMAN J. Thermal modeling of porous insertion electrodes[J]. Journal of the Electrochemical Society, 2003, 150(2): doi: 10.1149/1.1531194.

[本文引用: 1]

[4]

HU X S, LI S B, PENG HUEI. A comparative study of equivalent circuit models for Li-ion batteries[J]. Journal of Power Sources, 2012, 198: 359-367.

[本文引用: 1]

[5]

BABA N, YOSHIDA H, NAGAOKA M, et al. Numerical simulation of thermal behavior of lithium-ion secondary batteries using the enhanced single particle model[J]. Journal of Power Sources, 2014, 252: 214-228.

[本文引用: 1]

[6]

DOYLE M, FULLER T F, NEWMAN J. Modeling of galvanostatic charge and discharge of the lithium/polymer/insertion cell[J]. Journal of the Electrochemical Society, 1993, 140(6): 1526-1533.

[本文引用: 1]

[7]

李夔宁, 谢运成, 谢翌, 等. 基于电化学热耦合模型的富镍锂离子电池产热分析[J]. 储能科学与技术, 2021, 10(3): 1153-1162.

[本文引用: 1]

LI K N, XIE Y C, XIE Y, et al. Analysis of heat production of nickel-rich lithium-ion battery based on electrochemical thermal coupling model[J]. Energy Storage Science and Technology, 2021, 10(3): 1153-1162.

[本文引用: 1]

[8]

LYU P Z, HUO Y T, QU Z G, et al. Investigation on the thermal behavior of Ni-rich NMC lithium-ion battery for energy storage[J]. Applied Thermal Engineering, 2020, 166: doi: 10.1016/j.applthermaleng.2019.114749.

[本文引用: 2]

[9]

DU S L, LAI Y Q, AI L A, et al. An investigation of irreversible heat generation in lithium ion batteries based on a thermo-electrochemical coupling method[J]. Applied Thermal Engineering, 2017, 121: 501-510.

[本文引用: 1]

[10]

张志超, 郑莉莉, 杜光超, 等. 基于多尺度锂离子电池电化学及热行为仿真实验研究[J]. 储能科学与技术, 2020, 9(1): 124-130.

[本文引用: 1]

ZHANG Z C, ZHENG L L, DU G C, et al. Electrochemical and thermal behavior simulation experiments based on multiscale lithium ion batteries[J]. Energy Storage Science and Technology, 2020, 9(1): 124-130.

[本文引用: 1]

[11]

韦雪晴, 邓海鹏, 周宇, 等. 锂离子电池组的三维电化学-热耦合仿真分析[J]. 储能科学与技术, 2022, 11(12): 3965-3977.

[本文引用: 1]

WEI X Q, DENG H P, ZHOU Y, et al. Three-dimensional electrochemical-thermal coupling simulation analysis of lithium ion battery pack[J]. Energy Storage Science and Technology, 2022, 11(12): 3965-3977.

[本文引用: 1]

[12]

FANG R Q, GE H, WANG Z H, et al. A two-dimensional heterogeneous model of lithium-ion battery and application on designing electrode with non-uniform porosity[J]. Journal of the Electrochemical Society, 2020, 167(13): doi: 10.1149/1945-7111/abb83a.

[本文引用: 1]

[13]

LIU S B, ZHANG H Y, XU X B. A study on the transient heat generation rate of lithium-ion battery based on full matrix orthogonal experimental design with mixed levels[J]. Journal of Energy Storage, 2021, 36: doi: 10.1016/j.est.2021.102446.

[本文引用: 1]

[14]

FARAG M, SWEITY H, FLECKENSTEIN M, et al. Combined electrochemical, heat generation, and thermal model for large prismatic lithium-ion batteries in real-time applications[J]. Journal of Power Sources, 2017, 360: 618-633.

[本文引用: 1]

[15]

LI J E, CHENG Y, AI L H, et al. 3D simulation on the internal distributed properties of lithium-ion battery with planar tabbed configuration[J]. Journal of Power Sources, 2015, 293: 993-1005.

[本文引用: 1]

[16]

LI H H, SAINI A, LIU C Y, et al. Electrochemical and thermal characteristics of prismatic lithium-ion battery based on a three-dimensional electrochemical-thermal coupled model[J]. Journal of Energy Storage, 2021, 42: doi: 10.1016/j.est.2021.102976.

[17]

云凤玲. 高比能量锂离子动力电池热性能及电化学-热耦合行为的研究[D]. 北京: 北京有色金属研究总院, 2016.

[18]

BAHIRAEI F, GHALKHANI M, FARTAJ A, et al. A pseudo 3D electrochemical-thermal modeling and analysis of a lithium-ion battery for electric vehicle thermal management applications[J]. Applied Thermal Engineering, 2017, 125: 904-918.

[19]

DONG T, PENG P, JIANG F M. Numerical modeling and analysis of the thermal behavior of NCM lithium-ion batteries subjected to very high C-rate discharge/charge operations[J]. International Journal of Heat and Mass Transfer, 2018, 117: 261-272.

[20]

VALO̸EN L O, REIMERS J N. Transport properties of LiPF₆-based Li-ion battery electrolytes[J]. Journal of the Electrochemical Society, 2005, 152(5): doi: 10.1149/1.1872737.

[21]

DOYLE M, NEWMAN J, GOZDZ A S, et al. Comparison of modeling predictions with experimental data from plastic lithium ion cells[J]. Journal of the Electrochemical Society, 1996, 143(6): 1890-1903.

[22]

WU W, XIAO X R, HUANG X S. The effect of battery design parameters on heat generation and utilization in a Li-ion cell[J]. Electrochimica Acta, 2012, 83: 227-240.

[本文引用: 1]

[23]

XU X B, ZHANG H Y, LIU S B, et al. Surrogate models for lithium-ion battery heat generation based on orthogonal experiments by eliminating external wire connection effect[J]. Applied Thermal Engineering, 2022, 213: doi: 10.1016/j.applthermaleng.2022.118655.

[本文引用: 1]

[24]

吴青余, 张恒运, 李俊伟. 校准量热法测量锂电池比热容和生热率[J]. 汽车工程, 2020, 42(1): 59-65.

[本文引用: 1]

WU Q Y, ZHANG H Y, LI J W. Calibrated calorimetry for measuring the specific heat capacity and heat generation rate of lithium-ion battery[J]. Automotive Engineering, 2020, 42(1): 59-65.

[本文引用: 1]

[25]

JI Y, ZHANG Y C, WANG C Y. Li-ion cell operation at low temperatures[J]. Journal of the Electrochemical Society, 2013, 160(4): doi: 10.1149/2.047304jes.

[本文引用: 1]

... 锂离子电池由于其高能量密度、快速充电和放电能力、无记忆效应和低自放电率而被认为是最受欢迎的动力源.然而，锂电池的性能与工作条件和环境温度密切相关.温度过低或者过高都对电池的使用寿命和内部组成产生不利影响^[1]. ...

... 锂离子电池是一个典型的电化学储能装置，包含了各种电化学反应、离子输运以及热传递过程.Bernardi等^[2]和Thomas等^[3]提出了电化学系统及多孔插入电极的通用能量平衡方程，包括不可逆热源项、可逆热源项以及相变和热容量变化项.国内外学者也基于不同产热模型对锂离子电池的电热特性进行了分析^[4-5]，但这些模型在使用过程有所简化，导致其内部各组件产热特性无法评估.Newman等^[6]提出了基于多孔电极理论的准二维模型(P2D模型)，考虑了液态电解液中的浓溶液和固态电极的传输动力学，并利用Fick扩散定律来描述固相材料中锂离子的插入和脱离，自此该模型得到了广泛应用^[7-11].Lyu等^[8]利用P2D模型研究了高镍NMC储能锂离子电池的热特性.Du等^[9]建立了一个电化学-热耦合模型来考虑锂离子电池内部不可逆热量的产生和演变.张志超等^[10]利用P2D模型研究了层叠式锂离子电池在多尺度的电化学与热行为. ...

... P2D模型参数包括正极和负极颗粒的平均半径、正极和负极的平衡电位、正极和负极的熵系数、电极活性材料的最大锂离子浓度、电极的固相扩散系数、电解液扩散系数、电极的固相电导率、电解液电导率等.本工作中NMC锂电池的参数来自于文献[8, 15-22]，室温下电化学参数见表4，变温动态参数见表5. ...

... 全电池的熵系数通过测量不同温度下的开路电压得到.采用文献[12]中负极熵系数，并利用实测的全电池熵系数计算得到正极熵系数^[13-14]，如图3所示. ...

... 图4(a)展示了电池热测试实验图.使用K型热电偶和温度采集装置(Agilent 34972A)监测并记录电池温度，时间间隔为3 s.为了更准确地测量电池的放电温度，电池表面去除绝缘膜并焊接热电偶，直接测量电池的金属表面温度.在验证研究中，在室温25 ℃，电池外表面的对流传热系数确定为13.4 W/(m²·K)^[23].同时记录电压和电流.根据电池内部集流体结构构建集流体等效电路并得到负极集流体电阻为9.19 mΩ，正极集流体电阻为3.4 mΩ，该部分产热不能忽略，需要考虑在内. ...

... 本节分析在室温和不同放电倍率下(0.5 C、1 C、2 C、3 C)电池内各组件的产热，不同倍率下的电池平均总产热符合文献[24]产热与倍率的拟合关系，并着重对0.5 C放电倍率下各组件的瞬时产热进行了分析对比.图6显示了室温下不同放电倍率的负极、电解液、集流体和正极的产热率.由图6(a)可以看出，尽管负极的极化热高于正极，但由于负极的可逆产热存在吸热阶段，因此负极产热量低，而正极产热率大于负极，容易局部升温导致安全问题.而随着放电倍率的增加，正极、负极、电解液与集流体的产热均有所增加.由表6可以看出，随着放电倍率的增加，负极产热比例先增加后减小，正极产热的比例在减小而集流体产热的比例在增加，这是由于集流体欧姆产热与电流平方成正比，导致产热比例快速上升.而从低温下的电池各组件产热比例可以看出，随着放电倍率的增加，正极和负极占比在减小，而集流体产热的比例在增加. ...

... 低温限制了锂离子在液体中的移动和在固体中的扩散，使得锂离子在电池内的传输动力不足，导致内阻和产热增加.另外，通常在低温下发现的电压反弹现象^[25]与内部电化学状态的关系也将在本节阐述. ...

〈

〉