首页 > 行业资讯 > 《储能科学与技术》推荐|赵浩东 等：增设导流孔及翅片化通道墙强化液冷板散热性能的新策略

《储能科学与技术》推荐|赵浩东 等：增设导流孔及翅片化通道墙强化液冷板散热性能的新策略

时间：2023-10-25 来源： 浏览：

《储能科学与技术》推荐|赵浩东 等：增设导流孔及翅片化通道墙强化液冷板散热性能的新策略

原创 赵浩东 张甫仁 等 储能科学与技术
储能科学与技术

esst2012

中文核心、科技核心和cscd核心期刊，化学工业出版社和中国化工学会主办，主编黄学杰研究员。投稿及下载官网：http://esst.cip.com.cn/CN/2095-4239/home.shtml；欢迎给公众号投稿

收录于合集
#2023年第10期 11 个
#研究 302 个

作者： 赵浩东( ), 张甫仁( ), 杜柏林, 李雪, 黄郅凯, 孙世政    

单位： 重庆交通大学机电与车辆工程学院

引用： 赵浩东, 张甫仁, 杜柏林, 等. 增设导流孔及翅片化通道墙强化液冷板散热性能的新策略[J]. 储能科学与技术, 2023, 12(10): 3108-3119.

DOI： 10.19799/j.cnki.2095-4239.2023.0364

摘 要   为了增强液冷板的散热性能，通过在液冷通道内添加翅片形成局部扰流是一种主要方法，但是该方法会导致压降增大。基于此，本文设计了一种通道内增设隔板并带有导流孔，同时对通道墙翅片化的新型液冷板结构。在以降低压降和平均温度为目标的情况下，首先通过单因素法分析和讨论了导流孔和导流翅片个数的影响，得出当导流孔和导流翅片的个数分别为4和11时，液冷板的综合散热性能最佳。为了进一步优化液冷板的散热性能，采用多目标优化方法对不同导流孔之间的距离( )以及导流孔与隔板起始处的距离( )进行优化，结果显示多目标优化后模型的综合性能有了进一步提升。其次通过正交实验，讨论了翅片的倾斜角度以及开口宽度对平均温度与压降的影响，优化结果表明，翅片间的开口宽度影响最大。通过正交实验优化后平均温度降低了0.869 ℃(2.33%)，压降降低了18.257 Pa(71.62%)。最后，讨论了雷诺数在100～400的范围内变动时，不同液冷板结构的努塞尔数、压降和综合评价指标的变化情况。本研究有助于推动电池热管理散热的应用。
关键词   翅片；导流孔；多目标优化；正交实验
目前随着社会的进步，汽车数量不断上升，导致有害气体和颗粒物的排放量增加，对环境和人体健康造成影响。因此一些国家逐渐停止对燃油汽车的生产，推动新能源汽车的制造，解决汽车污染物排放问题。纯电动汽车、混合动力汽车、燃料电池汽车是新能源汽车的典型代表。
对于电动汽车而言，电池扮演着重要的角色。目前电池热管理系统中常用的冷却方式为液冷与风冷。液体冷却方式所采用的介质相比于空气而言传热系数和热容量更高，因此液体冷却方式是更适合电池模组的冷却方式。除此之外，冷却方式还有相变材料冷却、热管冷却、混合冷却等。
液冷系统可以通过改变液冷通道的结构形式、翅片的形状或结构参数来改变流场的分布，从而改善散热情况。Xie等为了提高传热系数在液冷通道中增加了挡板，并分析了挡板高度和数量对散热性能的影响。得到的结果是增加挡板能够提高壁面温度，增加了冷却液和电池间的热交换。为了提升传热系数和散热性能，可以在液冷通道中增加挡板改变通道结构，但是通道结构型液冷板往往会导致压力损失(压降)增大。相较之下，翅片型液冷板在相近的散热表现下，压力损失通常更低。Khoshvaght-Aliabadi等发现与挡板相比，采用翅片的液冷系统具有更好的散热性能。Law等提出了一种具有翅片的MPFHS(micro-pin-fin heat sink，微型针翅式散热器)，翅片可以形成第二通道，以提高散热性能。
可见，在液冷板内增设翅片可以有效地增强液冷板的散热性能和温度均匀性。该方法的实质就是通过在局部位置增设翅片形成局部扰流，从而提升系统的散热性能，但压降增大是该方法难以避免的问题。基于此，为了能提升液冷板的散热性能并实现压降的降低，本文提出在通道内增设隔板，形成局部扰流来强化散热性能；通过在隔板上增设导流孔，并将通道墙进行翅片化，以实现压降的降低。

1 模型

1.1　几何模型

本文根据传统的液冷板结构温度较高，提出一种新型液冷板结构。传统结构和新型结构的对比如图1所示。由于方形电池的长、宽分别为140 mm、65 mm，因此液冷板整体尺寸为140 mm×65 mm×3 mm；通道内的宽度为2 mm；隔板尺寸为51 mm×1 mm；冷却液入口和出口尺寸均为5 mm×2 mm。其他相关参数见表1。

图1     液冷板几何结构及尺寸

表1   液冷板相关尺寸参数

1.2　数值模型

本文中的模型采用ANSYS Workbench进行网格划分，在ANSYS Fluent中进行瞬态共轭问题的求解。控制方程通过压力分解器进行求解。采用基于压力的分离算法和二阶迎风格式对控制方程进行离散。在冷却液质量流量为0.1～1 g/s的基础上计算出雷诺数小于2300，按规定为层流状态。有关方程及计算如下所示。
动量方程：
(1)
能量方程：
(2)
连续性方程：
(3)
其中， P 代表微元体受到的压力； T 为温度； 为比热容； 为流体的耗散项； k 为传热系数； ρ 表示流体的密度； u ， v ， w 分别表示流体流动的速度在 x ， y ， z 方向上的分量。

1.3　边界条件

液冷板入口采用质量流量入口，流量为0.5 g/s，温度为25 ℃，出口采用压力出口，其值为0 Pa，冷板与环境间的传热系数为5 W/(m·℃)，液冷板加热面的热通量为3412 W/m 2 ，为了便于计算作出以下假设：
（1）设置在电池与液冷板之间不存在接触热阻；
（2）冷却剂是稳定的、不可压缩的流体；
（3）冷却剂在流动过程中处于稳定状态；
（4）锂电池在工作过程中受热均匀；
（5）液冷板和冷却剂的热材料性能不随温度变化。

1.4　电池产热的综合评价指标

本次需要计算出热阻、传热系数、努塞尔数等参数值来对散热能力和泵功率等进行评价。计算方法如下所述。
雷诺数(Reynolds number)是一种无量纲常数，表示流动情况，流体的状态是层流还是湍流可以通过雷诺数的大小区分。其表达式如下：
(4)
其中， 表示水力直径； 表示动力黏度。
水力直径 可表示如下：
(5)
其中， 表示液冷通道的横截面积； p 表示通道的周长； w 和 h 分别表示液冷通道横截面的宽和高。
努塞尔数的计算如下：
(6)
式中， h 为对流换热系数； k 为液冷板的传热系数。对流换热系数 h 的表达式为
(7)
式中， Q 为热流密度。
平均热阻的计算公式如下：
(6)
T W 表示液冷板加热面的平均温度； q 为冷板表面热通量。
外界水泵给液冷板输送冷却液需要的能耗用泵功率 表示，对应的计算公式如下：
(8)
其中， 表示压降； 表示体积流量。
(9)
其中， 是进口处压力； 为出口的压力。
鉴于液冷板换热性能的提高可能会导致其压降升高，因此需要综合考虑这两种性能对液冷板的影响。本文选择综合评价指标 FOM 作为液冷板性能的综合评价指标，以综合考虑压降和换热性能的影响。其表达式如下：
(10)
式中， 表示初始液冷板的传热系数； 表示优化液冷板的传热系数。

1.5　网格独立性

为了平衡计算精度和时间，合适的网格数就显得尤为重要。这里，使用ANSYS Workbench对5种网格工况进行独立性分析。局部网格的展示如图2(a)所示。由图2(b)可知，当网格的数量超过887725时，平均温度 和最高温度 的差值均处于5%的误差范围内，因此本文采用的网格数为887725。

图2     网格展示与独立性分析

2 数值模拟验证

本文使用图3(a)所示的液冷板实验平台来验证CFD方法的有效性，该实验平台分为3个部分：发热系统，冷却系统，数据采集系统。将质量流量设为0.5 g/s、1 g/s、1.5 g/s和2 g/s，然后进行实验研究。实验中所采用的加工仪器以及液冷板如图3(b)所示，在仿真模型基础上，在基础液冷板尺寸上增加7 mm宽度来添加螺栓孔，因此实际液冷板尺寸是154 mm×79 mm×6 mm。

图3     实验平台与液冷板模型
实验过程如图4所示，实验仪器主要有温度采集仪(LD5200)、电源(SS-1003)、恒温箱(SPX-150B)、流量计(LZB-4WB)、蠕动泵(550)、恒温水浴锅(HH-4)、计算机等。恒温箱控制温度，恒温水浴锅提供恒定温度的冷却液。智能蠕动泵将冷却液带入转子流量计。质量流量通过转子流量计与液冷板模型进口相连接来改变其大小。此次采用的冷却液无杂质。另外，本实验采用铝块和加热棒来替代电池。使用5个PT100热电阻来测量温度，图5展示了温度采集点的布置形式。通过温度采集仪对得到的相关数据进行记录和保存。为了保证25 ℃的恒温环境，将实验装置放入恒温箱。

图4     实验平台的搭建及原理图

1—电源；2—培养箱；3—流量计；4—蠕动泵；5—恒温水浴锅；6—计算机；7—温度控制器；8—铝块和液冷板

图5     液冷板上 5 个测温点的分布图
通过表2得出不同流量下实验结果与仿真结果间的误差小于5%，因此仿真结果可靠。

表2   实验数据与仿真数据对比

3 结果与讨论

3.1　导流孔与导流翅片个数的讨论

液冷板是对称的结构，因此只在一侧隔板上增设导流孔，另一侧则做相同变化。可以通过在隔板上增设导流孔的方式来减小冷却液的流动距离，从而达到降低压降的效果。将隔板命名为隔板A、隔板B、隔板C，如图6所示。在隔板上增加1～2个长度为2 cm，宽度与隔板相同的导流孔，共8种组合情况，如表3所示。对得到的8种组合进行仿真得出每个组合的 、 和 FOM 值如图7所示。通过图7得到当组合方式为 时综合性能最佳。

图6     液冷板内隔板名称

表3   8种组合形式

图7     8 种组合方式的液冷板仿真结果
在上述基础上，将铝块进行分割，如图8所示，得到多个导流翅片，使冷却液二次分流，增加液冷板与冷却液的热交换面积。观察图9得知，当翅片数量增加时， 略有起伏， 则不断增大，而有11个翅片时，综合评价指标 FOM 最大，此时综合性能最好，故组合方式为 ，其 为36.641 ℃， 为11.644 Pa，相较于初始结构， 和 分别降低了0.682 ℃和13.846 Pa。

图8     导流翅片的分割形式

图9     导流翅片个数的影响
3种液冷板的温度云图和流速云图，如图10所示。从温度云图得知，相对于传统液冷板结构和新型液冷板结构，优化液冷板结构具有更均匀的温度分布。从流速云图可以得知，这种结果是因为传统液冷板结构的流量分布不均，存在一些死流区域，而新型液冷板结构添加了隔板以减缓冷却液的流速所致。优化液冷板结构由于添加了导流孔和导流翅片，死流区域明显减少，温度分布更均匀。

图10     液冷板优化前后的温度云图和流速云图

3.2　基于多目标优化方法的导流孔的结构优化

基于目前的优化液冷板结构为 ，本节将对4个导流孔进行多目标优化。以隔板上4个导流孔(分别命名为m、n、p、q)作为优化变量，从而通过对4个导流孔的位置分布的讨论来优化液冷板的平均温度和压降。
本研究所采用的方法为最优拉丁超立抽样设计，以应用于探究图11中包括4个导流孔在内的设计变量。其中，m距离隔板A起始处长度为 ，n距离隔板B起始处长度为 ，p与n之间的距离为 ，q距离隔板C起始处长度为 。上述变量的值范围可参见表4。

图11     设计变量参数示意图

表4   变量的名称及其取值范围

在变量空间中随机选取52个样本点，分布图如图12所示。对样本点进行仿真模拟可以得到每个样本点对应的响应值平均温度 和压降 ，对应关系见表5。 、 的计算公式为

图12     样本点分布图

表5   样本点与其对应的响应值

(11)
(12)
其中， 和 分别表示 和 关于自变量( ， ， ， )的函数； 和 分别表示 和 的误差。
为了判断拟合精度，计算计算模型的样本决定系数 ：如果 <0.9，表示该模型不合格；如果 ≥0.9，则表示该模型合格。由图13所示，通过函数关系式得到平均温度 与压降 的样本决定系数 值分别为0.95644和0.99757，满足要求，因此该模型可以用于下一步的算法优化。表6为实验值与预测值的相关参数。

图13     响应值 和 的拟合精度

表6   预测值与实际值的相关参数表

由ASA算法(基于模拟退火算法改进的自适应模拟退火算法)通过对原始模型产生扰动得到新模型，然后根据Metropolis接受准则来确定新模型的接受概率，接受概率表达式如下：
(13)
式中， 为新解目标函数值； 表示初始解目标函数值； 表示温度值。
经过ASA算法寻优后，得到最优设计变量的参数为： =26.907 mm， =3.276 mm， =0 mm， =19.709 mm。此时 =36.55 ℃， =10.31 Pa。
对上述得到的优化液冷板结构进行模拟得到仿真结果，然后与ASA算法得到的预测值进行对比，结果如表7。由计算得知预测结果与仿真结果之间的相对误差值均小于2%，说明了预测结果的可靠性。本次优化得到的液冷板结构与初始液冷板相比 降低0.728 ℃ (1.95%)， 降低14.992 Pa(58.82%)。图14展示出了初始结构与优化结构的温度云图与压力云图。观察温度云图得知，优化后液冷板的高温区域减少，温度分布更加合理；观察压力云图得知，优化后液冷板的压降降低，通过改变导流孔位置可以适当地降低压降。

表7   预测值与仿真值的对比

图14     多目标优化前后温度云图变化

3.3　基于正交实验的液冷板优化结构

在上述最优的液冷板结构基础上，本节将对导流翅片的角度以及开口宽度进行讨论。由于液冷板是对称结构，本次只对一侧的翅片进行变化，另一侧做相同的变化即可。由图15所示，将4排翅片分别命名为组D、组E、组F、组G，4组翅片与垂直方向的夹角为 、 、 、 ，夹角的大小设置为70°、80°、90°、100°、110°共5个等级。其次在每个翅片的中间开口，开口宽度为 W ， W 的取值设置为0 mm、2 mm、4 mm、6 mm、8 mm。然后列出5因素5水平的正交表如表8所示，共25组实验数据，以 与 为评价指标对每组实验进行仿真模拟，最后通过极差分析得到最优结构。

图15     各因素示意图

表8   25组实验数据及仿真结果

通过表8内的数据，根据下列公式可得到评价指标的极差分析结果：
(14)
(15)
(16)
(17)
(18)
其中， i 表示水平， j 代表因素； 是某一实验结果； 是相同水平的实验结果之和； n 表示水平的个数； 则为方差。
对液冷板平均温度与压降的极差与方差分析，由表9可知对液冷板平均温度 影响从大到小依次为 ，对应的最佳组合为 W (3) α 1 (1) α 2 (1) α 3 (2) α 4 (1)，命名为U；当考虑到压力损失时，由表10可知各因素对压降影响程度的顺序为 ，对应的最佳组合为 W (5) α 1 (2) α 2 (3) α 3 (3) α 4 (2)，命名为V。通过计算得到组合U和组合V的 FOM 值分别为1.794和1.859，因此选择组合V。
组合V与组合U的温度云图与压力云图如图16所示，压力云图直观地反映出了组合V的压降明显低于组合U。对比组合V与初始液冷板，组合V的 与 分别降低了0.869 ℃(2.33%)和18.257 Pa(71.62%)。

图16     2 种组合的温度云图与压力云图

3.4　冷板特性分析

根据图17，将4种不同结构的液冷板分别命名为case0(传统液冷板)、case1(初始液冷板)、case2(3.1节中得到的优化液冷板)、case3(3.3节中得到的最优液冷板)，然后讨论4种液冷板结构的 、 和 FOM  3个指标随雷诺数的变化情况。

图17     4 种液冷板结构
3.4.1　不同液冷板结构的热特性分析
是对流热量和传导热量的比值，它是无量纲常数， 越大，传热性能越好。图18展示了4种液冷板的 值随雷诺数 Re (100～400)变化的趋势。

图18     4 种液冷板结构的 值随雷诺数变化趋势
由图18可知，4种液冷板结构均呈现出随着雷诺数的增加，努塞尔数也逐渐增加的趋势。相比之下，传统液冷板结构case0的 值最小，这表示在对流换热方面传统结构还存在不足。其原因包括传统直通道液冷板结构会导致冷却液在液冷板内滞留时间短，使得换热时间较少；还包括传统液冷板结构的液冷通道面积小，造成换热面积减少。结构case1的 值随着雷诺数增加不断增大，因为case1液冷板在传统液冷板基础上添加了隔板，这就使得冷却液在液冷板内滞留的时间变长，延长了换热时间。结构case2、case3是在case1的基础上通过增设导流孔和翅片的方式来进一步优化，优化后使冷却液在液冷通道中形成二次分流，增加换热面积，所以case2和case3的 值比case1的大。从图18可知，当雷诺数的取值大于250时，优化结构case3的 值最大。
3.4.2　压降特性分析
液冷板所需要的泵功率能够通过压降的大小来反映，泵功率越小，所需的压降越小，代表外部给液冷板供给的能量越小。如图19所示，展示了4种液冷板结构中压降随雷诺数(100～400)的变化。

图19     4 种结构的压降随雷诺数的变化趋势
根据图19得知，当雷诺数增加时，4种结构(case0、case1、case2、case3)的压降均呈现出上升的势头。这是因为随着冷却液质量流量增加，会导致通道内的翅片对冷却液造成阻流，使流动阻力变大，所以压降增加。相比之下，结构case1压降最大，因为隔板的添加导致冷却液流动距离加长。结构case2、case3的压降比case1低，是因为在降板上添加了导流孔并优化了翅片结构，使得部分冷却液流动距离大大减小。当雷诺数取值在100～300时，结构case3的压降最小。
3.4.3　综合性能的分析
本节将讨论在雷诺数变化的情况下，不同微通道结构的综合性能。本次讨论规定在不同的雷诺数下case1均为基础的工况，综合评价指标 FOM 的取值为1，在此基础上计算出其他3种工况的综合评价指标 FOM 的值。结果如图20所示。

图20     不同雷诺数下对应的 FOM 值
根据图20可知，雷诺数在100～400内时，相较于case1，其他3种工况的 FOM 取值均大于1，这是因为其他3种工况提供的热对流性能比产生的流动阻力大。随着雷诺数不断变大，case0的综合评价指标 FOM 值会减小，这表示与热对流性能相比，压降对结构case1影响更大；case2和case3的综合评价指标 FOM 值则出现了先增后减的情况，并且 FOM 的值远大于1，因此case2和case3的热性能更好。由图20可知， FOM 取得最大值时，雷诺数为150，液冷板结构为case3，相对于case1提升91.05%。

4 结论

本文提出了一种带隔板的新型液冷板结构，分析并讨论了一些参数(导流孔和导流翅片的数量、导流孔位置、翅片倾斜角度，以及翅片开口尺寸等)对液冷板性能的影响。得到的结论如下：
（1）经过对导流孔和翅片讨论后，得知导流孔个数为4，导流翅片个数为11时，液冷板的散热性能最好。
（2）经过对4个导流孔位置进行多目标优化后，液冷板的 降低了0.728 ℃ (1.95%)， 降低了14.992 Pa(58.82%)。
（3）通过正交实验，对翅片与垂直线的倾斜角度与翅片中间开口宽度进行优化后，液冷板的 降低了0.869 ℃(2.34%)， 降低了18.257 Pa(71.62%)。
（4）在热特性分析后，可知当雷诺数变大后，经过优化的液冷板具有更大的 值。尤其是雷诺数为150时，结构case3的 FOM 值最大，相比于case1提升了91.05%。

第一作者： 赵浩东（1997—），男，硕士研究生，研究方向为电池热管理，E-mail：；

通讯作者： 张甫仁，博士后，教授，研究方向为电池热管理，E-mail：zfr@cqjtu.edu.cn。

本刊推荐

（点击图片可跳转到相应文章合集）

邮发代号：80-732

联系热线：010-64519601/9602/9643

投稿网址：http://esst.cip.com.cn/CN/2095-4239/home.shtml

版权：如无特殊注明，文章转载自网络，侵权请联系cnmhg168#163.com删除！文件均为网友上传，仅供研究和学习使用，务必24小时内删除。
相关推荐
热门信息