第44卷第19期包装工程
2023年10月PACKAGING ENGINEERING·273·
刘欢欢1,2，任晓龙2，张泽斌1*
（1.郑州大学机械与动力工程学院，郑州450001；2.广东佛山联创工程研究生院，广东佛山528300）
摘要：目的将锂离子电池在5C放电倍率、环境温度298.15 K的工作条件下的最大温度和温差控制在适宜范围内。方法在锂电池两侧插入带通道的仿生树状液冷板，建立电池模块液冷散热模型。 首先验证液冷数值模型的准确性，然后通过正交试验设计研究分支流道角度、入口流速、分支流道宽度和液冷板厚度对冷板散热性能的影响，并通过极差分析对4种影响因素进行重要程度排序，确定最佳的参数组合。结果经优化后，锂电池的最高温度为302.4 K、温差为3.4 K，满足设计要求，但冷却剂在流动循环过程中存在泄漏风险。结论提出的树状液冷板满足锂电池热管理性能的要求，优化后电池模块的表面温度在安全范围内，同时电池的温度均匀性显著提高。应选择导热系数高的冷板材料，并注重冷板结构设计，防止冷却剂泄漏。 关键词：锂离子电池；液冷板；树状结构；正交试验
中图分类号：TK172；TM912.9文献标识码：A 文章编号：1001-3563(2023)19-0273-10
DOI：10.19554/jki.1001-3563.2023.19.035
Numerical Optimization of Bionic Tree-like Channel Liquid Cold Plate on
Lithium-ion Battery
LIU Huan-huan1,2, REN Xiao-long2, ZHANG Ze-bin1*
(1. School of Mechanical and Power Engineering, Zhengzhou University, Zhengzhou 450001, China;
2. Guangdong Shunde Innovative Design Institute, Guangdong Foshan 528300, China)
ABSTRACT: The work aims to achieve the goal that the maximum temperature and temperature difference of lithium-ion battery are within the suitable temperature range when the discharge rate is 5C and the ambient temperature is 298.15 K.
Bionic tree-like liquid cold plates with channels were inserted on both sides of the lithium-ion battery to establish the liquid cold plate heat dissipation model of the battery module. Firstly, the accuracy of the
numerical model of liquid cold plate was verified, and then the effect of the angle of branch runner, the inlet velocity, the width of branch runner and the thickness of liquid cold plate on the cooling performance of cold plate was studied by orthogonal experimental design.
The importance of the four factors was ranked by range analysis, and the optimal parameter combination was determined.
The maximum temperature of the optimized scheme battery was 302.4 K and the temperature difference was 3.4 K, which met the design requirements. However, there was a risk of coolant leakage during the flow cycle. The tree-like liquid cold plate can meet the requirements of battery thermal management performance. The optimized battery module has a safe range of surface temperature, and the battery temperature uniformity is significantly improved. Cold plate materials with high thermal conductivity should be selected, and attention should be paid to the structural design of cold plate to prevent coolant leakage.
KEY WORDS: lithium-ion battery; liquid cold plate; bionic tree-like structure; orthogonal experiment
收稿日期：2023-04-14
·274·包装工程2023年10月
锂离子电池具有能量密度高、循环寿命长等特点，在储能系统中得到了广泛应用。储能系统中的锂电池在放电、充电过程中会产生热，并且难以及时散发热量，易引起锂电池组之间的热量聚集，出现运行温差过大等现象，过高的温度会导致锂电池起火甚至爆炸[1]。分析近年来储能系统的安全事故发现，火灾隐患多由锂电池热失控引起[2]。由此可见，开发高效的锂电池热管理系统对于锂电池模块的性能和使用寿命至关重要，锂电池热管理系统应保证锂电池模块之间的温差在5 K 以内，最高温度保持在293.15~313.15 K的最佳范围内[3]。目前，锂电池热管理系统有4种常用的冷却方法：空气冷却、液体冷却、相变材料冷却和热管冷却[4]。研究人员还提出了一种混合热管理系统，通过结合2种常用的冷却方法来提高锂电池模块的冷却效率[5]。
与其他冷却系统相比，液冷系统需要水泵、管道、冷板等配件，导致其重量和成本增加，同时冷却剂存在泄漏风险。由于液体冷却系统具有较高的传热系数和较好的导热性能，目前已成为锂电池热管理研究的重点。刘霏霏等[6]分析了并行流道冷板的入口流速、通道数量、通道宽度等参数对锂电池最高温度和温差的影响，结果表明，通道数量对冷板散热性能的影响最显著。付平等[7]通过比较传统蛇形结构冷板、外置式双进出水口蛇形冷板和内嵌式双进出水口蛇形冷板的散热性能后发现，双进出水口的内置式结构是一种较理想的冷板结构。Li等[8]采用仿生波纹并行通道冷板，研究了放电速率、流道形状、冷却液质量流量等因素对锂电池最高温度和温差的影响，发现具有仿生波纹流道的冷板在冷却效率上优于并行直流道冷板。Wang等[9]采用正交试验设计方法研究了蛛网式通道冷板对锂电池散热性能的影响。在
蛛网形通道的结构参数中，通道宽度对冷板冷却性能的影响最大，通道角度对冷板冷却性能的影响最小。杜柏林等[10]建立了相变材料/液冷复合的锂电池热管理系统，并结合多目标优化算法得出了锂电池组最高温度和温差最小的参数组合。
传统的蛇形流道冷板存在功耗大的缺点，并行通道冷板存在温度梯度大的缺点，因此在设计冷板结构时需要综合考虑散热效果和功耗情况。自然系统存在的流动结构（如树冠、叶脉、血管等）可为设计新的冷板提供思路，这些流道系统在很大程度上遵循默里定律[11]。默里定律给出了树状分支结构在2个连续分支层次上分支直径的最优关系。
文中提出一种树状流道结构冷板，首先通过单因素分析分支流道角度、宽度、冷板厚度和冷却液流速对冷板散热性能的影响，并采用正交试验设计、极差分析法对各因素进行优化分析，以得到各因素影响冷板散热性能的重要程度及最佳参数组合。
1 模型和方法
1.1 几何模型
选取8个方形锂电池和9个液冷板组成锂电池组模块。锂电池的长度、宽度、高度分别为118、63、13 mm，锂电池组模块结构如图1a所示。在2个相邻的锂电池之间存在铝制树状液冷板，液冷板的长度和宽
度分别等于锂电池的长度和宽度。由于锂电池组采用几何对称排列方式，因此采用简化的锂电池模块。选择2个1/2锂电池和1个液冷板作为计算域，如图1b所示，以降低计算成本，将顶部和底部表面设置为对称面。树状冷板的二维结构如图1c所示。模型中使用的材料热物性参数如表1所示。
图1 锂电池热管理方案示意图
Fig.1 Schematic diagram of the lithium-ion battery thermal management system
表1 材料热物性参数
Tab.1 Thermo-physical properties of materials
热物性参数密度/(kg·m-3) 导热系数/(W·m−1·K−1)比热容/(J·kg−1·K−1) 动力黏度/(Pa·s) 锂电池 2 450 3.9 1 108
铝 2 719 202.4 871
水998.2 0.6 4 182 0.001 003
第44卷 第19期 刘欢欢，等：锂离子电池仿生树状通道液冷板数值优化 ·275·
1.2 控制方程
冷却剂的质量、动量、能量守恒方程[12]见式（1）—（3）。
l ()0u ρ∇⋅= (1)
2l ()u u p u ρμ⋅∇=−∇+∇
(2)
2l pl l l l ()c u t k t ρ⋅∇=∇
(3)
式中：u 为冷却剂速度，m/s ；l ρ为冷却剂的密度，kg/m 3；p 为冷却剂静压，Pa ；μ为冷却剂动力黏
度，Pa·s ；c pl 为冷却剂的比热容，J/(kg·K)；t 1为冷却剂的温度，℃；k l 为冷却剂的导热系数，W/(m·K)。
冷却剂雷诺数见式（4）。
l uD Re ρμ
= (4)
式中：D 为冷却通道水力直径，m 。 液冷剂的压降计算见式（5）。
2
2
22i i
L u u
P f
D ρρξΔ=+
(5)
式中：f 为摩擦因数；L 为流道的长度，m ；ξi
为局部损失系数。
液冷剂与锂电池的换热量的计算见式（6）。
tra l w ()q hA t t =−
(6)
式中：h 为对流换热系数，W/(m 2·K)；A 为换热面积，m 2；t 1、t w 分别为冷却剂、壁面的温度，℃。
1.3 初始和边界条件
采用Ansys Fluent 求解器，通过计算可知冷却剂流动入口的雷诺数均小于2 300，因此将冷却剂流动设置为层流模型条件。将环境温度设定为298.15 K ，锂电池表面为自然对流，对流传热系数为5 W/(m 2·K)。入口边界为速度条件，出口边界为零压力条件。
1.4 产热模型验证
采用Rao 等[13]建立的锂电池产热模型，以恒定电流放电720 s 来计算体积发热率，放电速率为5C （C 表示锂电池的放电倍率），锂电池的体积发热率为240 000 W/m 3。从图2可知，得到的锂电池最高温度与文献[13]中锂电池的最高温度吻合较好。由此可见，这里建立的锂电池发热模型合理，后续研究将以该模型为基础。
2 数值仿真
2.1 条件假设
1）冷却剂流动是不可压缩的。
2）锂电池是一个均匀的产热部件。
3）将液冷板的上下表面视为无厚度表面[14]。
图2 产热模型的验证 Fig.2 Verification of the heat
generation model
2.2 网格数量无关性验证
文中在对模型进行网格划分时采用四面体网格，为了检验网格数量的无关性，采用6种网格数量（分别为25万、33万、42万、56万、90万、127万）进行仿真分析。作为测定网格数量无关性的初始方案，将试验条件设置为冷却剂入口流速0.1 m/s 、分支流道角度45°、分支流道宽度2 mm 、冷板厚度2 mm ，后面在进行仿真计算时若无特别说明，则试验条件不变。计算结果如图3所示，当网格数量从56万增至127万时，温度最大变化量约为0.1 K ，相对
误差在0.03%内。
为了在确保计算精度足够的同时减少时间成本，选择网格数量56万进行后面的仿真计算。
图3 锂电池最大温度的网格独立性试验结果 Fig.3 Grid independence test results of the maximum lithium-ion battery temperature
2.3 液冷数值模型验证
为了研究数值模型的准确性，对铝制微通道散热器进行了仿真模拟，Fazeli 等[15]进行了相应实验，模型几何结构如图4a 所示。在散热器底面施加恒定的
·276· 包 装 工 程 2023年10月
图4 液冷数值方法验证
Fig.4 Validation of the liquid cold plate numerical method
热流，进水流速为0.071 m/s ，进口冷却液温度为27 ℃。将模拟结果沿底面中心线的温度与实验数据进行了比较，结果如图4b 所示。可以看出，最大相对误差小于4%，因此文中的模拟结果可靠。
3 仿真结果分析
3.1 分支流道角度的影响
锂电池的温度和温差随时间的变化曲线如图5所示。从图5a 可以看出，在前40 s 内，各方案中锂电池的最高温度基本相同，升温速率在初始时刻最大，之后逐渐减小，这是因为冷却剂流过冷板通道的每个角落并与锂电池进行对流换热需要一定的时间。在40 s 后，各方案中锂电池的最高温度开始出现差异，即分支流道角度对冷板散热性能产生的影响开始显现。
在分支流道角度为30°时，锂电池的最高温度、温差分别为306.2、6.7 K ，此时冷板的散热性能最优。在不同分支流道角度下，锂电池的最高温度均在理想工作范围内，但温差大于5 K 。当分支流道角度从30°增至60°时，锂电池的最高温度增加了1.5 K ，温差增加了1.7 K 。在分支流道角度大于60°后，锂电池的最高温度、温差分别仅增加了0.1、0.2 K ，可见在分支流道角度大于60°后对树状冷板散热性能的影响较小。不同分支流道角度下锂电池温度分布云图如图6所示。可以看出，锂电池温度沿着冷却剂的流动方向不断升高，这是因为冷却剂不断被锂电池加热，导致换热速率不断减小。从图5、图6、表2中也可看出，随着分支流道角度的增加，最高温度、温差和压降都呈增大的趋势。这是因为流动换热面积随着分支角的增大而减小，阻碍了冷却剂与锂电池的对流换热。
图5 不同分支流道角度下锂电池的最高温度和温差随时间的变化 Fig.5 Evolution of maximum temperature and temperature difference in the
lithium-ion battery under different branch runner angles
第44卷 第19期 刘欢欢，等：锂离子电池仿生树状通道液冷板数值优化 ·277·
图6 不同分支流道角度下锂电池温度分布云图
Fig.6 Temperature distribution cloud map of the lithium-ion
battery under different branch runner angles
表2 不同分支流道角度下冷却剂的压降
Tab.2 Pressure drop of coolant under different branch
runner angles 角度/(°) 30 40 50 60 70 80 压降/Pa
40.4
42.9
43.944.2
44.7
44.9
3.2 冷却剂入口流速的影响
在不同流速下，锂电池的最高温度、温差随时间
的变化过程如图7所示，
锂电池的温度分布云图如图8所示。可以看出，当进口速度为0.04 m/s 时，锂电池的最高温度达到313.9
K ，超出了锂电池的理想工作温度范围。在其他5个速度下锂电池的最高温度均低于308.1 K ，只有当流速大于0.20 m/s 时，锂电池的温差才能保持在5 K 以下。随着流速的增加，锂电池的最高温度不断降低，这是因为流速的增加导致对流传热系数增大，在相同时间内冷却剂可以带走更多
的热量。
当流速从0.04 m/s 增至0.16 m/s 时，最高温度和温差显著降低，最高温度降低了9.2 K ，温差降低了8.4 K 。当入口流速从0.16 m/s 变化到0.24 m/s 时，最高温度和温差仅分别降低了1.3、1.1 K ，锂电池的最高温度和温差缓慢下降。这是因为冷却剂在流道入口段的局部传热系数较高，沿着流动方向逐渐降低，并趋于稳定。入口段的长度与雷诺数呈正相关，当入口流速从0.04 m/s 增至0.16 m/s 时，与入口流速从0.16 m/s 变化到0.24 m/s 时相比速度的变化范围大，这意味着雷诺数增加的值较大，因此平均传热系数增加得较多。由式（6）可知，冷却剂带出的热量较多，因此最高温度和温差变化明显。由此可见，增大冷却剂流速可以显著改善锂电池模块的冷却效果，但是进一步增大入口流速会减缓锂电池模块的冷却效果，同时会增大冷却剂的压降。如表3所示，当冷却剂的入口
流速从0.04 m/s 增至0.16 m/s 时，
速度增加了0.12 m/s ，
图7 不同入口流速下锂电池的最高温度和温差随时间的变化
Fig.7 Evolution of maximum temperature and temperature difference in the lithium-ion
battery under different inlet flow rates
豫公网安备41160202000603
站长QQ:729038198
关于我们
投诉建议