1. 引言
锂离子电池作为一种典型的电化学储能系统,涵盖多种电化学反应、质量传递、电荷转移以及热量传递等过程。很多学者对复杂的电化学系统进行了研究。加州大学伯克利分校的Bernardi等人 [1] 对电化学系统提出了一个通用的能量平衡方程,该方程包括不可逆热源项,可逆热源项和相变与热容变化项。2003年Newman等人 [2] 基于热力学基本参数,推导出一种多孔插入式电极的能量平衡方程,该方程包含不可逆热、可逆热、热容、电池中的化学反应热和混合热。随着计算机数值仿真技术的发展,已经能够建立更完善的数学物理模型,全面系统地分析计算各物理场的相互作用机理,探究其内部规律。
张志超等人 [3] 对锂离子电池充放电过程进行了产热研究,研究从电池的材料、电池循环寿命、充放电的电流大小以及环境温度和荷电状态五个维度来展开,全面概述各种因素对动力电池的影响。云凤玲 [4] 研究了高比能的三元锂离子电池的过电位与熵系数,也发现电池在放电过程中的温度变化源于不可逆热与可逆热,其中不可逆热与电池极化过电位与欧姆产热有关,而可逆热与电池熵系数有关,同时,对电池进行循环并测量其熵系数,发现循环前后熵系数变化趋势一致。张凯兰 [5] 研究了锂离子电池在极端情况下的热失控机制,并利用产热模型进行了电化学–热–力的耦合。吴彬 [6] 提出了基于时域和频域的两种热参数估计方法,与交流阻抗类似,他把电学模型应用到热学,推导出了热阻抗谱模型,并用该模型对动力电池的导热系数进行参数辨识。时天禄 [7] 对锂离子电池进行了电化学–热耦合并进一步进行电池热管理研究,电化学–热模型从微观角度分析了电极中锂离子浓度过电位等的变化同时结合冷却结构,完成了电化学–热耦合下的动力电池热管理的研究。He等人 [8] 为不同尺寸的圆柱形锂离子电池建立了电化学–热耦合模型。Hamza等人 [9] 利用多尺度电化学–热学模型对一个电动双层电容器进行建模,研究其产热和温度分布。
上述研究者运用电化学模型研究了电池的产热情况,但并没有从锂离子电池机理出发,研究不同电极内部参数大小对产热的影响,也缺少对多孔电极中电极参数对容量的影响,而多孔电极粒径大小会直接影响电极材料的设计与制备,本文将依靠电化学模型完成多孔电极的不同粒径分布的仿真研究。
2. 模型与参数选择
2.1. 数学模型
选用额定容量为2.7 Ah、以人造石墨为正极材料的锂离子电池进行模型分析。其基本参数见表1。
Table 1. Basic parameters of the battery
表1. 电池的基本参数
我们建立了一个微观尺度的电化学–热耦合模型,以研究多孔电极不同尺寸粒径下的极化热。
P2D模型在这里被用来模拟离子传输并在建模中使用了以下假设。
1) 在负极和正极中的传输过程用同一套传输方程来描述。
2) 隔膜的多孔结构通过引入与孔隙度相关的平均曲率系数来描述。
3) 等效传质系数,如负极和正极的等效固相电子传导率和等效液相离子传导率与各自的孔隙率有关。
4) 在模拟中不考虑附属的物理成分,如导电剂和粘合剂。
5) 固体和液体两相共存,在任何一点上都被认为是体积平均理论框架下的混合均匀相。
基于P2D模型上述假设,该模型用五个控制方程表示,首先对于电极颗粒而言,在电极颗粒内部的物质传输分别用Fick第二定律与欧姆定律来表示,其中Fick第二定律将固相球形颗粒径向方向上的任意一点锂离子浓度随时间(t)的变化率转化为方程右侧中该点锂离子浓度随空间的变化率。而用欧姆定律描述了电子在固相颗粒沿电极厚度方向的移动。在电极颗粒内部的电荷守恒用Buter-Volmer方程表示。其次对于电解液而言,在电解液中的输运通常用浓溶液理论修正的Nernst-Planck方程来描述,而基于浓溶液理论,给出液相中锂离子的物质守恒方程。
极化热的来源主要与固液项界面的电化学反应有关,而电化学反应是锂离子电池产生持续充放电能力的基础,通常利用Buter-Volmer方程来描述电池内部的电极反应动力学。
(1)
其中,j0代表交换电流密度,αa和αc为传递系数,表示电极电位对反应活化能的影响程度,h是过电位。B-V方程将反应电流密度与过电位联系起来,表明在两个电极上均存在着可逆的电化学反应,即同时存在着氧化和还原反应电流,两个电化学反应相反。
(2)
其中,k为反应速率常数,cl,ref为液相锂离子参考浓度,cs,surf为固相颗粒表面锂离子浓度,cs,max为固相颗粒最大锂离子浓度。过电位用方程定义:
(3)
其中电极的极化热为:
(4)
2.2. 参数选择
本部分利用P2D电极模型进一步分析电极设计问题,即负极颗粒粒径分布对电极产热的影响,从而为负极材料粒度选型提供指导意见。由于负极颗粒的粒径分布主要影响负极区域的曲折系数以及活性颗粒内部的离子扩散过程。颗粒粒径分布函数可以通过扫描电子显微镜或者微米X光计算机成像获得的重构图像进行统计得到或者利用粒度仪等仪器获得。本文选取四个维度粒径进行仿真分析。如表2所示,每个维度下包含不同的粒径大小与固液相的体积分数。
Table 2. Different particle size distribution parameters
表2. 不同颗粒粒径分布参数
本模型采用锂离子电池模块进行求解,模型基于P2D模型进行计算,在进行产热及容量计算时,使用解析的方式求出该模型的瞬态解。
3. 不同粒径仿真分析
通过计算不同粒径完成了对该电池的极化热分析。从图1中可以看出,随着粒径分布的增加,电池极化热也增加。在同一粒径分布下,可以看到随着放电时间的增加,电池的极化热呈现先降低后增加的状态。关于粒径分布,可以发现,当粒径分布越大,对该电池的放电时间也存在显著的影响。当粒径分布为D1时,电池可以放3600 s的电量,而当粒径分布为D4时,电池放电时间出现衰减,没法在原有规定的放电时间下完成放电。由此可以知道,提高与完善负极颗粒粒径分布集中度、降低颗粒粒径大小可有效改善电池内部离子在电解液以及活性颗粒内部的有效传输速率,从而造成了电池极化热的降低,并且提高电池的可用容量。
Figure 1. Polarization heat of batteries with different particle size distributions
图1. 不同粒径分布下的电池极化热
图2采用D3粒径分布分离出了锂离子电池的正负极极化热与过电势,从图2(a)中可以看出,在放电前期和中期,正负极极化热持续稳定输出,负极极化热大于正极极化热。而在放电末期,正极极化热基本保持不变,负极极化热迅速增加。而对该产热进行深入分析,可以看到图2(b)中的过电势随时间的变化,负极过电势在放电末期的过电位迅速增加,锂离子电池放电结束时负极锂离子消耗过大,导致“低锂状态”使其负极内部极化加剧,从而极化热产热较大。
Figure 2. Positive and negative: (a) Polarization heat, (b) Overpotential
图2. 正负极:(a) 极化热,(b) 过电势
4. 结论
本文利用电化学模型研究了锂离子电池在不同粒径分布下的产热与容量情况。通过对锂离子电池进行不同粒径分布的1C放电倍率的仿真研究,发现电池的颗粒粒径整体大小会影响电池的极化产热,结果表示颗粒粒径越大电池的极化热越大。而在同一粒径下,电池极化热随着放电过程呈现先降低后增加的趋势。由于粒径变动影响孔隙率等体积分数,进而影响了电极颗粒传输参数,从而造成了电池放电时间的衰减。