变接触面液冷系统的电池模组温度一致性研究

干年妃 孙长乐 刘东旭 雷富强

摘   要：为了使电池系统单体间的温度具有较好的一致性，设计了一种在单体间隙中填充铝柱的液冷热管理系统，建立了单体的电化学-热耦合模型，对比研究了不同入口流速、不同截面边长和高度组合的铝柱液冷系统对电池热性能的影响. 研究结果表明，截面边长和高度梯度变化的铝柱液冷系统的冷却性能和单体间的温度一致性都要优于截面边长和高度为定值的系统. 在放电倍率为3 C，入口流速为0.10 m/s时，此组合方式下电池模块中单体间的最大温差保持在3.72 ℃以内，满足电池系统热管理的需求. 最后，针对此模型对不同放电倍率进行了仿真验证，结果显示单体间温度具有较好的一致性.

关键词：锂离子电池;电池热管理;液冷;变接触面;温差

Abstract：It is far more difficult to reduce the temperature differences among cells in a battery module than to reduce the maximum temperature of module. In order that the cells can have better temperature distributing conformity， this paper designed a liquid cooling thermal management system by packing aluminum column in the gaps among  the cells and set up an electrochemical-thermal coupling model for the cells. Comparative study of the impacts on the thermal properties of the cells exerted by liquid cooling system of aluminum columns under the combination modes of different entrance velocity，different section length and height was also conducted. The results indicated that the cooling performance of the aluminum column liquid cooling system with gradient variation of cross-section length and height and the temperature distributing conformity among the cells were better than those with the fixed values of section length and height. Under the condition that the discharge rate was 3C and the entrance velocity was 0.1 m/s，the maximum temperature difference among the cells in the module group based on the proposed coupling model was kept within 3.72 ℃，meeting the requirements of thermal management of the battery system. Finally， different discharge rates were verified through the simulation based on this model， indicating that the cells had better temperature distributing conformity.

Key words：lithium ion battery;battery thermal management;liquid cooling;variable contact surface;temperature difference

伴随着全球石油短缺和环境污染的巨大压力，新能源汽车由于在节能减排和环保方面相对于传统车辆有明显的优势，近些年已呈现爆发式的增长. 发展新能源汽车在全球也已成不可逆转之势，世界主要汽车生产国家纷纷制定了将全面禁售燃油车的时间表. 以锂离子电池为动力源的纯电动汽车作为新能源汽车的主要发展方向之一受到了各国政府和主要汽车制造商的大力扶持和推动[1]，锂离子电池具有较高的能量密度和较长的循环寿命，是新能源汽车，尤其是纯电动汽车的理想动力源. 温度对电池寿命和工作性能产生重要影响，其中最高温度和温度一致性是影响电池寿命和工作性能的最重要的两个因素[2]. 在充放电过程中，过高的温度会导致电池破裂、失火甚至爆炸，而在低温环境下，动力电池内部的内阻增大，电池的充放电容量急剧降低，在低温环境下使用电池，會增加电池的耗损，降低电池的使用寿命[3]. 因此，为了将电池温度控制在可接受范围内，并在充放电过程中保持单体间均匀的温度分布，就需要对电池进行热管理（BTM）[4].先前的学者对电池热管理方面已经做了大量的研究，包括电池单体生热模型的建立[5-8]、冷却结构的布置、冷却材料的选择等等[9-10]. 就目前而言，根据冷却介质，电池热管理可分为空冷、液冷、相变材料冷却和热管冷却，以及多种方式的耦合冷却. 由于空气冷却的冷却性能有限，不能满足充放电倍率较大的电池系统的散热要求，而相变冷却虽然冷却性能良好，但其复杂的结构和成本限制了其在实际中的应用[11]. 因此，对于高充放电率、高产热率的电池组而言，液体冷却具有明显的优势[12]. 方形电池的液冷设计往往通过冷却板或直接浸入液体介质中来实现[13]. 微通道液体冷却被广泛应用到方形电池的热管理中[14-17]. Qian等人[18]研究了不同微通道数量及入口方向对电池冷却性能的影响;Xie等人[19]提出了一种Y形微通道冷却系统;Salimpour等人[20]设计了一种环形的冷却结构，并对其进行了优化;Deng等人[21]研究了蛇形微通道的管道数和温度一致性的关系. 相对于方形电池，圆柱形电池模组液冷系统的研究相对较少. Zhao等人[22]在充放电过程中，采用带有液体流动的波状通道对圆柱形电池模组进行冷却;胡兴军等人[23]针对圆柱形电池模组设计了8种间接接触的液冷结构并对比了它们的冷却效果;Basu等人[24]设计的铝片导电元件与圆柱形锂离子电池组紧密接触;Zhao等人[25]对于42110圆柱形电池设计了一种基于套筒的微通道液冷系统. 之后Rao等人[26]设计了一种通过改变铝块与电池接触面积来调整单体间温差的液冷系统.

根据先前学者的研究可知，圆柱形电池液冷系统主要通过改变结构形式来达到降低最高温度和温差的目的. 相对于降低模组的最高温度，调节单体间温度的一致性更加困难. 根据Rao等人[26]对于改变铝块长度来实现降低电池单体间温差的研究，本文在此基础上提出了一种新的结构，将铝柱放置在18650电池单体间的空隙中，分别研究了铝柱截面边长和铝柱的高度对电池温度的影响，并通过改变铝柱的高度和截面边长来增加与电池单体的接触面梯度以降低温差，并研究了不同放电倍率下单体的温差，为进一步提高电池热管理系统的工作性能提供参考和指导.

1   模型与方法

1.1   单体模型

电池单体采用松下NCR18650PF，最大容量为2 900 mA·h，标称电压为3.7 V. 锂离子电池因其不同的长度尺度和几何复杂性，各层在层法向的尺寸通常为数十微米，但在电池片方向达数十厘米，且电池内部为螺旋绕制结构，所以电池单体三维模型的建立较为复杂. 但是就锂离子电池的传热过程而言，只存在一个不同方向的传热速率的差异，因此可以将电池三维模型中的活性电池材料简化成具有各向异性导热特征的均质材料. 为了方便仿真，采用快速仿真的电化学-热耦合模型，其耦合过程如图1所示. 一维单电池模型用于对电池化学性质进行建模以及计算电化学模型的产热，并将热量耦合至三维模型，三维模型用于对温度进行建模，并将温度耦合至一维模型. 此模型大大减少了计算时间，适用于研究电池包热管理.

1.2   冷却系统结构

图2展示了电池模块液冷散热系统的示意图，电池模块共由21个18650电池单体组成，沿x方向有3个单体，沿y方向有7个单体，单体间的间距都为1 mm. 铝柱放置在4个单体形成的空隙中，铝柱与单体的接触面均为单体的弧形状态，如图2（c）所示，铝柱截面形状为一个正方形减去单体弧形所占的面积，将正方形边长a记为铝柱的截面边长，铝柱的高度记为h，管道分布如图2（b）所示，管道竖直通过铝柱，管道直径为3 mm. 为节省空间，提高铝柱质量，最外侧的4个铝柱侧边与管道中心的距离定为2 mm. 此模型具有对称性，且由于中间一列单体两侧都布置有铝柱，所以外侧两列单体更具有参考性，如图2（c）所示，将外侧一列单体标号为1～7，以便于后续的研究. 考虑导热系数和黏度，采用液态水作为冷却介质. 电池、铝柱和冷却液的热物理性质见表1.

1.3   平衡方程

锂离子电池生热量包括副反应热、反应热、极化热和焦耳热. 由于过充和过放、电解质分解及自放电产生的副反应热几乎为零[27]，故只考虑剩下的三部分热量. 在进行充放电时，锂离子和电子发生嵌入和脱嵌，运动过程中产生热量如式（1）所示，两者遵循电荷守恒定律. 反应热Q1为：

1.4   边界条件与网格验证

根据冷却液的性质、管道直径及入口流速，可计算雷诺数来确定是采用层流模型还是湍流模型. 本文最大入口流速为0.2 m/s，雷诺数为671.36，由此可知，本文采用层流模型. 在所有仿真中，入口冷却液温度、电池模组初始温度及环境温度均设为25 ℃. 冷却液的入口边界条件采用速度边界条件，具体速度值在仿真过程中确定;出口边界条件采用压力边界条件，法向压力设置为0. 电池模块外层有保持架、集流片、固定架等覆盖件，且整个电池包也由外壳所包覆，电池单体、冷却通道和铝块的对流传热系数设为2 W/（m2·K）.

利用有限元软件COMSOL求解耦合传热问题，COMSOL具有强大的多物理场耦合功能. 由于控制方程的高度非线性和模型中几何尺度的不同，计算精度和计算时间依赖于网格和求解器. 模型都采用自由四面体网格，以管道出口边界的温度为参考，测试了几种网格密度对其影响，以保证解的网格无关性.

2   结果与讨论

由于电池的热量都是通过铝柱传递到冷却管道，因此铝柱和电池单体的接触面积与冷却性能有很大的关系. 铝柱与电池的接触面积可以通过改变铝柱的截面大小和铝柱的高度来进行调节. 在接下来的研究中，我们将分别讨论铝柱截面大小和铝柱高度这两个变量对电池冷却性能的影响.

2.1   铝柱截面大小对温度的影响

选定铝柱高度为定值40 mm，只改變截面边长以达到改变与电池接触面积大小的目的，截面边长a以1 mm的间距从7 mm变化至14 mm;冷却液入口流速分别设置为0.05、0.10、0.15和0.20 m/s，电池以3 C倍率放电至终止电压，放电时长为1 200 s. 放电结束后整个电池模块的最高温度和最大温差如图4所示.

由图4可知，冷却液入口流速对温度影响显著，流速越快，最高温度及最大温差越低. 随着入口流速从0.05 m/s增加到0.20 m/s，铝柱截面边长为7 ～ 14 mm的电池模组最高温度降低了4.72～5.69 ℃;最大温差从7.17～9.16 ℃全都降低到5 ℃以内. 铝柱截面大小对温度同样也有所影响，截面越大，最高温度及最大温差越低. 随着截面边长从7 mm增加到14 mm，在流速分别为0.05、0.10、0.15和0.20 m/s下模块的最高温度分别降低1.97、1.41、1.14和0.99 ℃，最大温差分别降低1.99、1.36、1.09和0.95 ℃.

过高的入口流速虽然能有效降低电池的温度和温差，但会造成水泵更高的功率损耗. 图5是入口流速为0.10 m/s时，在不同铝柱截面边长下1～7号电池单体的最高温度分布.

由图5可以看出，离管道入口越近的单体温度越低，且存在一个明显的梯度变化，通过调整铝柱与单体的接触面积来缓和这种梯度变化，调整接触面积的方法为调整铝柱截面边长. 每个铝柱的截面边长a由边长梯度因子α确定：

图7分别显示了α = 0、α = 0.5及α = 1的情况下，电池以3 C放电倍率放电至终止电压后的温度云图. 在不同边长梯度因子α下的电池模组最高温度都出现模组外侧的离出口最近的两个电池单体表面，最低温度出现在模组中间一列离管道入口最近的电池单体表面，当α从0增加到1时，电池模组离管道入口较近的单体温度有所上升，其温度一致性表现更好.

图8为α = 0、α = 0.5及α = 1时1～7号单体的最高温度曲线. 由图8可见，当截面边长梯度增大时，离管道入口越近的单体最高温度呈现上升趋势，当α = 1时单体间的最高温度差相比于α = 0时下降了0.45 ℃. 可得出结论，改变铝柱截面大小对改善单体间温差有一定的影响，但由于截面边长的调整范围有限，导致这种影响并不显著. 因此还需要用其他方式去缓解温度梯度.

2.2   铝柱高度对温度的影响研究

选定铝柱截面边长为定值10 mm，铝柱高度以5 mm为间隔从20 mm变化至60 mm;冷却液入口流速分别设为0.05、0.10、0.15和0.20 m/s，电池以3 C倍率放电至终止电压，放电时长为1 200 s. 放电结束后整个模块的最高温度和最大温差分别如图9所示. 由图9可见，在不同铝柱高度下，冷却液的入口流速依然是对温度影响最大的因素，随着入口流速从0.05 m/s增加到0.20 m/s，铝柱高度为20 ～ 60 mm的电池模组最高温度降低了4.67～5.71 ℃;最大温差从7.33～9.29 ℃全都降低到5 ℃以内. 铝柱高度对温度同样也有所影响，铝柱越高，最高温度及最大温差越低. 随着铝柱高度从20 mm增加到60 mm，在流速分别为0.05、0.10、0.15和0.20 m/s下模块的最高温度分别降低2.51、2.03、1.72和1.56 ℃，最大温差分别降低1.96、1.78、1.56和1.44 ℃.

在不同铝柱高度的情况下，离管道入口越近的单体温度同样越低，且存在一个明显的梯度变化，通过调整铝柱与单体的接触面积来缓和这种梯度变化. 调整接触面积的方式为调整铝柱高度，每个铝柱高度h由高度梯度因子β确定：

由于管道结构的限制，铝柱高度最高为60 mm，如果铝柱过短，会造成单体两端温度过高，因此高度梯度因子β的取值范围为0 ≤ β ≤ 6. 图11为β = 0、β = 4、β = 5以及β = 6时的冷却系统结构及温度分布云图. 在不同高度梯度因子β下的电池模组最高温度都出现在模组外侧两列离出口最近的两个电池单体表面，最低温度出现在模组中间一列离管道入口最近的电池单体表面，但相比于β = 0，β逐渐从4增加到6时，铝柱高度梯度越大，离管道入口越近的单体温度上升越明显，而离管道出口相近的单体温度有略微下降的趋势. 当β = 6时，高度梯度最大，模组的温度一致性表现最好.

图12为β = 0、β = 4、β = 5以及β = 6时1～7号单体最高温度分布点线图. 通过图12可以看出，梯度因子越大，离管道入口越近的单体温度上升越明显，整个曲线表现得愈平缓. 当β = 0、β = 4、β = 5以及β = 6时1～7号单体最高温度差分别为3.17、2.59、2.24和1.78 ℃;相比于β = 0，当β = 4、β = 5和β =6时，单体最高温度差分别降低0.58、0.93和1.39 ℃. 可得出结论，在相对较低的入口速度下，通过调节铝柱高度的梯度可以有效缓和单体间的温度差异.

2.3   铝柱高度与截面的组合方式对温度的影响

由于对铝柱截面边长和高度的单一参数调整并不能得到单体间足够小的温度梯度，为了使铝柱与电池单体的接触面积梯度足够大，采用铝柱高度与截面的组合方法对模块温度进行研究，取最大截面边长梯度因子为α = 1，取最大高度梯度因子β = 6，这样可以使铝柱的电池单體的接触面积梯度达到最大. 图13为铝柱高度与截面的组合方式下的冷却系统结构图及温度分布云图. 由图13可见，铝柱与电池的接触面积梯度已经达到最大，单体间温度分布也较为均匀.

图14显示的是当取最大截面边长梯度因子α=1、最大高度梯度因子β=6时与截面边长取定值14 mm、高度取定值60 mm时1～7号电池单体最高温度对比图. 可以看出，当截面边长梯度因子α=1与高度梯度因子β=6组合时，曲线较为平缓，相比于铝柱截面边长和高度为定值时1～7单体最高温度差降低1.83 ℃，整个模组的最大温差仅为3.72 ℃.

2.4    不同放电倍率下的温度分布研究

由于前文中的模型只是针对放电倍率为3 C的极限工况而言的，此模型在中低放电倍率下的温度分布一致性仍需要观察. 利用前文中的模型，观察放电倍率分别为0.5 C、1 C和2 C下的温度分布. 图15为不同放电倍率下模组的温度云图及单体最高温度分布图.

在0.5 C放电倍率下，电池模组的最高温度为25.33 ℃，最大温差为0.31 ℃;在1 C放电倍率下，电池模组的最高温度为25.87 ℃，最大温差为0.82 ℃;在2 C放电倍率下，电池模组的最高温度为27.48 ℃，最大温差为2.05 ℃;在0.5 C和1 C放电倍率下，1～7号电池单体的最高温度曲线表现较为平缓，在2 C放电倍率下，曲线表现出上升的趋势，但也在较小的区间内，说明此模型能够有效调节单体间的温差分布.

3   结   论

本文设计了一种在圆柱形锂离子电池间隙中填充铝柱的液冷热管理系统，并对系统的散热性能进行了讨论. 研究了铝柱截面边长和铝柱高度对电池模块散热性能的影响. 针对冷却系统，首先选定一个铝柱的高度为定值，通过改变铝柱的截面边长大小来观察系统的散热性能. 结果表明，随着截面边长的增大，电池模块的最高温度和最大温差均减小. 为改善温度分布的一致性，采用变截面边长来缓和单体间的温度梯度. 当α = 1时，单体间的最高温度差降低了0.45 ℃. 之后选定一个铝柱的截面边长为定值，通过改变铝柱的高度来观察系统的散热性能. 结果显示，随着铝柱高度的增加，电池模块的最高温度及最大温差均减小，但仍存在一个温度梯度，采用变铝柱高度来缓和这种温度梯度. 研究了不同β值对温度一致性的影响，当β = 6时单体间的温度一致性表现最好，相比于β = 0，当β = 6时单体间的最高温度差降低了1.39 ℃. 为使单体间的温度梯度达到最小，采用α = 1和β = 6的组合方式，经过观察，此组合方式下，模组的最大温差可降至3.72 ℃，相比于铝柱截面边长和高度为定值时单体间最高温度差降低了1.83 ℃.

文章的最后，研究了模型對于中低放电倍率下的温度分布. 在0.5 C和1 C放电倍率下，单体间最大温差均在1 ℃以内;在2 C放电倍率下，单体间最大温差为2.05 ℃，说明此电池冷却系统对调节电池温度一致性表现优异.

参考文献

[1]    SHAFIEI E，LEAVER J，DAVIDSDOTTIR B. Cost-effectiveness analysis of inducing green vehicles to achieve deep reductions in greenhouse gas emissions in New Zealand[J]. Journal of Cleaner Production，2017，150（1）：339—351.

[2]    XIA G D，CAO L，BI G L. A review on battery thermal management in electric vehicle application[J]. Journal of Power Sources，2017，367（1）：90—105.

[3]    PING P，WANG Q S，HUANG P F，et al. Thermal behaviour analysis of lithium-ion battery at elevated temperature using deconvolution method[J]. Applied Energy，2014，129（15）：261—273.

[4]    RAVDEL B，ABRAHAM K M，GITZENDANNER R，et al. Thermal stability of lithium-ion battery electrolytes[J]. Journal of Power Sources，2003，119/121（1）：805—810.

[5]    FENG X，LU L，OUYANG M，et al. A 3D thermal runaway propagation model for a large format lithium ion battery module[J]. Energy，2016，115（1）：194—208.

[6]    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.

[7]    CHIEW J，CHIN C S，TOH W D，et al. A pseudo three-dimensional electrochemical-thermal model of a cylindrical LiFePO4 /graphite battery[J]. Applied Thermal Engineering，2019，147（25）：450—463.

[8]    YE Y H，SHI Y X，CAI N S，et al. Electro-thermal modeling and experimental validation for lithium ion battery[J]. Journal of Power Sources，2012，199（1）：227—238.

[9]    DENG Y W，FENG C L，E J Q，et al. Effects of different coolants and cooling strategies on the cooling performance of the power lithium ion battery system：a review[J]. Applied Thermal Engineering，2018，142：10—29.

[10]  WANG Y，GAO Q，WANG G H，et al. A review on research status and key technologies of battery thermal management and its enhanced safety[J]. International Journal of Energy Research，2018，42（13）：4008—4033.

[11]  CHEN D F，JIANG J C，KIM G H，et al. Comparison of different cooling methods for lithium ion battery cells[J]. Applied Thermal Engineering，2016，94（5）：846—854.

[12]  E J Q，HAN D D，QIU A，et al. Orthogonal experimental design of liquid-cooling structure on the cooling effect of a liquid-cooled battery thermal management system[J]. Applied Thermal Engineering，2018，132（5）：508—520.

[13]  VANGILS R W，DANILOV D，NOTTEN P H L，et al. Battery thermal management by boiling heat-transfer[J]. Energy Conversion and Management，2014，79：9—17.

[14]  AN Z，SHAH K，JIA L，et al. A parametric study for optimization of minichannel based battery thermal management system[J]. Applied Thermal Engineering，2019，154（25）：593—601.

[15]  PANCHAL S，KHASOW R，DINCER I，et al. Thermal design and simulation of mini-channel cold plate for water cooled large sized prismatic lithium-ion battery[J]. Applied Thermal Engineering，2017，122（25）：80—90.

[16]  HUANG Y Q，MEI P，LU Y J，et al. A novel approach for lithium-ion battery thermal management with streamline shape mini channel cooling plates[J]. Applied Thermal Engineering，2019，157（5）：113623.

[17]  XU J，LAN C J，QIAO Y，et al. Prevent thermal runaway of lithium-ion batteries with minichannel cooling[J]. Applied Thermal Engineering，2017，110（5）：883—890.

[18]  QIAN Z，LI Y M，RAO Z H. Thermal performance of lithium-ion battery thermal management system by using mini-channel cooling[J]. Energy Conversion and Management，2016，126（15）：622—631.

[19]  XIE G N，SHEN H，WANG C C. Parametric study on thermal performance of microchannel heat sinks with internal vertical Y-shaped bifurcations[J]. International Journal of Heat and Mass Transfer，2015，90：948—958.

[20]  SALIMPOUR M R，MENBARI A. Analytical optimization of constructal channels used for cooling a ring shaped body based on minimum flow and thermal resistances[J]. Energy，2015，81（1）：645—651.

[21]  DENG T，ZHANG G D，RAN Y. Study on thermal management of rectangular Li-ion battery with serpentine-channel cold plate[J]. International Journal of Heat and Mass Transfer，2018，125：143—152.

[22]  ZHAO C R，CAO W J，DONG T，et al. Thermal behavior study of discharging/charging cylindrical lithium-ion battery module cooled by channeled liquid flow[J]. International Journal of Heat and Mass Transfer，2018，120：751—762.

[23]  胡興军，惠政，郭鹏，等. 锂离子电池组间接触液冷散热结构研究[J]. 湖南大学学报（自然科学版），2019，46（2）：49—55.

HU X J，HUI Z，GUO P，et al. Research on the indirect contact liquid cooling heat dissipation structure of lithium-ion battery pack[J]. Journal of Hunan University（Natural Sciences），2019，46（2）：49—55.（In Chinese）

[24]  BASU S，HARIHARAN K S，KOLAKE S M，et al. Coupled electrochemical thermal modelling of a novel Li-ion battery pack thermal management system[J]. Applied Energy，2016，181（1）：1—13.

[25]  ZHAO J T，RAO Z H，LI Y M. Thermal performance of mini-channel liquid cooled cylinder based battery thermal management for cylindrical lithium-ion power battery[J]. Energy Conversion and Management，2015，103：157—165.

[26]  RAO Z H，QIAN Z，KUANG Y，et al. Thermal performance of liquid cooling based thermal management system for cylindrical lithium-ion battery module with variable contact surface[J]. Applied Thermal Engineering，2017，123：1514—1522.

[27]  SAW L H，YE Y H，TAY A A O. Electrochemical-thermal analysis of 18650 lithium iron phosphate cell[J]. Energy Conversion and Management，2013，75：162—174.

[28]  CHEN J，KANG J Q，TAN Z X. Analysis of thermal performance of 18650 Li-Ion battery based on an electrochemical-thermal coupling model[J]. Hans Journal of Chemical Engineering and Technology，2018，8（2）：97—107.