增设导流孔及翅片化通道墙强化液冷板散热性能的新策略

摘要

为了增强液冷板的散热性能，通过在液冷通道内添加翅片形成局部扰流是一种主要方法，但是该方法会导致压降增大。基于此，本文设计了一种通道内增设隔板并带有导流孔，同时对通道墙翅片化的新型液冷板结构。在以降低压降和平均温度为目标的情况下，首先通过单因素法分析和讨论了导流孔和导流翅片个数的影响，得出当导流孔和导流翅片的个数分别为4和11时，液冷板的综合散热性能最佳。为了进一步优化液冷板的散热性能，采用多目标优化方法对不同导流孔之间的距离()以及导流孔与隔板起始处的距离()进行优化，结果显示多目标优化后模型的综合性能有了进一步提升。其次通过正交实验，讨论了翅片的倾斜角度以及开口宽度对平均温度与压降的影响，优化结果表明，翅片间的开口宽度影响最大。通过正交实验优化后平均温度降低了0.869℃(2.33%)，压降降低了18.257Pa(71.62%)。最后，讨论了雷诺数在100～400的范围内变动时，不同液冷板结构的努塞尔数、压降和综合评价指标的变化情况。本研究有助于推动电池热管理散热的应用。关键词

翅片；导流孔；多目标优化；正交实验目前随着社会的进步，汽车数量不断上升，导致有害气体和颗粒物的排放量增加，对环境和人体健康造成影响。因此一些国家逐渐停止对燃油汽车的生产，推动新能源汽车的制造，解决汽车污染物排放问题。纯电动汽车、混合动力汽车、燃料电池汽车是新能源汽车的典型代表。对于电动汽车而言，电池扮演着重要的角色。目前电池热管理系统中常用的冷却方式为液冷与风冷。液体冷却方式所采用的介质相比于空气而言传热系数和热容量更高，因此液体冷却方式是更适合电池模组的冷却方式。除此之外，冷却方式还有相变材料冷却、热管冷却、混合冷却等。液冷系统可以通过改变液冷通道的结构形式、翅片的形状或结构参数来改变流场的分布，从而改善散热情况。为了提高传热系数在液冷通道中增加了挡板，并分析了挡板高度和数量对散热性能的影响。得到的结果是增加挡板能够提高壁面温度，增加了冷却液和电池间的热交换。为了提升传热系数和散热性能，可以在液冷通道中增加挡板改变通道结构，但是通道结构型液冷板往往会导致压力损失(压降)增大。相较之下，翅片型液冷板在相近的散热表现下，压力损失通常更低。发现与挡板相比，采用翅片的液冷系统具有更好的散热性能。提出了一种具有翅片的MPFHS(micro-pin-finheatsink，微型针翅式散热器)，翅片可以形成第二通道，以提高散热性能。可见，在液冷板内增设翅片可以有效地增强液冷板的散热性能和温度均匀性。该方法的实质就是通过在局部位置增设翅片形成局部扰流，从而提升系统的散热性能，但压降增大是该方法难以避免的问题。基于此，为了能提升液冷板的散热性能并实现压降的降低，本文提出在通道内增设隔板，形成局部扰流来强化散热性能；通过在隔板上增设导流孔，并将通道墙进行翅片化，以实现压降的降低。1模型1.1几何模型本文根据传统的液冷板结构温度较高，提出一种新型液冷板结构。传统结构和新型结构的对比如图1所示。由于方形电池的长、宽分别为140mm、65mm，因此液冷板整体尺寸为140mm×65mm×3mm；通道内的宽度为2mm；隔板尺寸为51mm×1mm；冷却液入口和出口尺寸均为5mm×2mm。其他相关参数见表1。图1

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

液冷板相关尺寸参数1.2数值模型本文中的模型采用ANSYSWorkbench进行网格划分，在ANSYSFluent中进行瞬态共轭问题的求解。控制方程通过压力分解器进行求解。采用基于压力的分离算法和二阶迎风格式对控制方程进行离散。在冷却液质量流量为0.1～1g/s的基础上计算出雷诺数小于2300，按规定为层流状态。有关方程及计算如下所示。动量方程：(1)能量方程：(2)连续性方程：(3)其中，P代表微元体受到的压力；T为温度；为比热容；为流体的耗散项；k为传热系数；ρ表示流体的密度；u，v，w分别表示流体流动的速度在x，y，z方向上的分量。1.3边界条件液冷板入口采用质量流量入口，流量为0.5g/s，温度为25℃，出口采用压力出口，其值为0Pa，冷板与环境间的传热系数为5W/(m·℃)，液冷板加热面的热通量为3412W/m2，为了便于计算作出以下假设：（1）设置在电池与液冷板之间不存在接触热阻；（2）冷却剂是稳定的、不可压缩的流体；（3）冷却剂在流动过程中处于稳定状态；（4）锂电池在工作过程中受热均匀；（5）液冷板和冷却剂的热材料性能不随温度变化。1.4电池产热的综合评价指标本次需要计算出热阻、传热系数、努塞尔数等参数值来对散热能力和泵功率等进行评价。计算方法如下所述。雷诺数(Reynoldsnumber)是一种无量纲常数，表示流动情况，流体的状态是层流还是湍流可以通过雷诺数的大小区分。其表达式如下：(4)其中，表示水力直径；表示动力黏度。水力直径可表示如下：(5)其中，表示液冷通道的横截面积；p表示通道的周长；w和h分别表示液冷通道横截面的宽和高。努塞尔数的计算如下：(6)式中，h为对流换热系数；k为液冷板的传热系数。对流换热系数h的表达式为(7)式中，Q为热流密度。平均热阻的计算公式如下：(6)TW表示液冷板加热面的平均温度；q为冷板表面热通量。外界水泵给液冷板输送冷却液需要的能耗用泵功率表示，对应的计算公式如下：(8)其中，表示压降；表示体积流量。(9)其中，是进口处压力；为出口的压力。鉴于液冷板换热性能的提高可能会导致其压降升高，因此需要综合考虑这两种性能对液冷板的影响。本文选择综合评价指标FOM作为液冷板性能的综合评价指标，以综合考虑压降和换热性能的影响。其表达式如下：(10)式中，表示初始液冷板的传热系数；表示优化液冷板的传热系数。1.5网格独立性为了平衡计算精度和时间，合适的网格数就显得尤为重要。这里，使用ANSYSWorkbench对5种网格工况进行独立性分析。局部网格的展示如图2(a)所示。由图2(b)可知，当网格的数量超过887725时，平均温度和最高温度的差值均处于5%的误差范围内，因此本文采用的网格数为887725。图2

网格展示与独立性分析2数值模拟验证本文使用图3(a)所示的液冷板实验平台来验证CFD方法的有效性，该实验平台分为3个部分：发热系统，冷却系统，数据采集系统。将质量流量设为0.5g/s、1g/s、1.5g/s和2g/s，然后进行实验研究。实验中所采用的加工仪器以及液冷板如图3(b)所示，在仿真模型基础上，在基础液冷板尺寸上增加7mm宽度来添加螺栓孔，因此实际液冷板尺寸是154mm×79mm×6mm。图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个长度为2cm，宽度与隔板相同的导流孔，共8种组合情况，如表3所示。对得到的8种组合进行仿真得出每个组合的、和FOM值如图7所示。通过图7得到当组合方式为时综合性能最佳。图6

液冷板内隔板名称表3

8种组合形式图7

8种组合方式的液冷板仿真结果在上述基础上，将铝块进行分割，如图8所示，得到多个导流翅片，使冷却液二次分流，增加液冷板与冷却液的热交换面积。观察图9得知，当翅片数量增加时，略有起伏，则不断增大，而有11个翅片时，综合评价指标FOM最大，此时综合性能最好，故组合方式为，其为36.641℃，为11.644Pa，相较于初始结构，和分别降低了0.682℃和13.846Pa。图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.907mm，=3.276mm，=0mm，=19.709mm。此时=36.55℃，=10.31Pa。对上述得到的优化液冷板结构进行模拟得到仿真结果，然后与ASA算法得到的预测值进行对比，结果如表7。由计算得知预测结果与仿真结果之间的相对误差值均小于2%，说明了预测结果的可靠性。本次优化得到的液冷板结构与初始液冷板相比降低0.728℃(1.95%)，降低14.992Pa(58.82%)。图14展示出了初始结构与优化结构的温度云图与压力云图。观察温度云图得知，优化后液冷板的高温区域减少，温度分布更加合理；观察压力云图得知，优化后液冷板的压降降低，通过改变导流孔位置可以适当地降低压降。表7

预测值与仿真值的对比图14

多目标优化前后温度云图变化3.3基于正交实验的液冷板优化结构在上述最优的液冷板结构基础上，本节将对导流翅片的角度以及开口宽度进行讨论。由于液冷板是对称结构，本次只对一侧的翅片进行变化，另一侧做相同的变化即可。由图15所示，将4排翅片分别命名为组D、组E、组F、组G，4组翅片与垂直方向的夹角为、、、，夹角的大小设置为70°、80°、90°、100°、110°共5个等级。其次在每个翅片的中间开口，开口宽度为W，W的取值设置为0mm、2mm、4mm、6mm、8mm。然后列出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.257Pa(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综合性能的分析本节将讨论