液冷板性能多参数优化和性能分析.doc

            液冷板性能多参数优化和性能分析 　　摘 要：为了提高液冷板的综合性能，对影响其散热效果及流阻的重要因素――冷板厚度、流道的宽度及并行流道的数目进行了多参数数值优化，并对优化前后冷板的流动特性及热特性进行了分析。结果表明：冷板厚度为9 mm， 宽度为5.4 mm，并行流道数目为12时，冷板的性能最好，与优化前相比，冷板表面的最高温度T和流阻P分别降低了1.9 %和13.8；流道的宽度对冷板性能影响最为显著，并行流道数目对冷板性能影响最小。本文旨在为液冷板后续优化研究提供指导。 　　关键词：液冷板；流道参数；数值优化；性能分析 　　DOI：10.16640/j.cnki.37-1222/t.2018.06.122 　　随着电子技术的迅速发展，电子元件日趋向微小型化发展，集成度不断提高，使得其热流密度迅速增加，传统的风冷已无法满足散热需求。与风冷散热相比，水冷散热换热效率较好、均热性强、成本低、可靠性高，成为目前应用较为广泛的电子设备散热方式之一。但是由于电子元件发热功率的不断提高，如何实现电子设备的高效散热，保证电子元件可靠性是目前的研究热点之一，因此提高液冷板的热效率，改善其流动特性显得尤为重要[1-5]。 　　液冷板的优化研究目前集中在流道结构独立参数对冷板散热性能或流阻特性的影响上，而对流道结构多参数对液冷板综合性能的影响研究较少。由于实验研究电子设备散热系统的散热效果周期长且耗资巨大，因此，本文以提高液冷板的综合性能为目标，采用数值方法，基于RBF神经网络模型和多目标遗传算法对液冷板流道结构多参数进行优化求解及分析，为今后液冷板性能的研究提供指导。 　　1 研究对象 　　本文以1台内部流道为并联形式的液冷板为计算模型，其主要参数为：流量Q=40 mL/s，液冷板的长度a=200 mm，b=200 mm，厚度h=6.5 mm，流倒宽度d=6 mm，并行流道数目n=15，冷却工质为水。液冷板三维模型如图1所示。 　　2 研究方法 　　2.1 优化变量 　　综合考虑到冷板的实际使用要求及成本，保持原冷板长宽不变，流量不变，因此选取冷板流道的3个关键参数作为优化变量，分别为：冷板厚度h，流道宽度d，并行流道数目n，其中冷板两侧到流道上下底面的距离均为1.5 mm。 　　2.2 优化方法 　　冷板流道的多参数优化以相同工况下冷板表面最高温度最低和冷板进出口压差即流阻最小为目标函数，基于RBF神经网络构造冷板的性能预测模型，通过BBD中心组合试验和CFD数值计算获得一系列样本点，并对构建的神经网络预测模型进行训练，然后采用NSGA-II多目标遗传算法求得一组最优解，具体优化步骤如下： 　　a.确定所研究的优化变量； 　　b.对原模型进行CFD数值计算，获得原模型性能参数； 　　c.基于原模型流道参数初始值，确定优化变量取值，建立BBD中心组合试验表，按试验表进行数值计算，得到温度、流阻值； 　　d.将每组试验作为一个样本点，训练RBF神经网络预测模型，并确定模型准确性； 　　e.基于NSGA-II多目标进化遗传算法，求得性能最优的一组参数值，并计算验证。 　　3 数值计算方法 　　计算域包括冷板、水体、进口延长段和出口延长段共4部分。系统计算域模型如图2所示。所有流场均设为静止坐标系。 　　利用CFX15.0对冷板进行稳态数值计算，湍流模型采用k-ε湍流模型。计算域进口边界条件设为速度进口，为2 m/s，出口边界条件为压力出口，静压为0 Pa；工质为水，进口水温303.15 K；对冷板上底面进行持续加热，热流密度为3000 W/m2；冷板外壁面及进出口延长段壁面均设为绝热；使用PRESTO格式离散压力项，二阶迎风格式离散对流项，其余采用二阶中心差分格式；均采用无滑移壁面，粗糙度设为0.05 mm；网格关联采用GGI方式。 　　采用ICEM CFD对上述各部件进行四面体非结构网格划分，并对固液交界面进行了网格加密，网格质量均在0.2以上。为进行网格无关性分析，对原模型共划分5套不同网格数的网格，均采用CFX进行前处理，并保证边界条件、湍流模型等设置一致，以冷板表面的最高温度为评判标准，结果列于表1。从表1可看出，随网格数的增加，冷板表面的最高温度的误差在2%以内。综合计算经济性等各方面因素，后续研究均采用方案4的网格尺寸进行数值计算。 　　流道内流体与冷板的换热过程属于强对流换热，根据牛顿冷却定律，对流换热的换热量Q公式如下： 　　（1） 式中：hi为对流换热系数；Ai为固体壁面对流换热表面积；Twi和Tfi分别为固体壁面温度和流体温度。对原模型的数值计算结果进行分析，发现冷板的最高温度为313.21 K，流阻为5.83 kPa。 　　由式（1）可看出，对流换热面积直接影响冷板的散热效果