结构参数对不锈钢板式热沉流动传热性能影响的仿真研究.pdf

结构参数对不锈钢板式热沉流动传热结构参数对不锈钢板式热沉流动传热性能影响的仿真研究性能影响的仿真研究何鸿辉1，龚洁2*，周盈2，杨雪2，金玮安2（1.中国空间技术研究院；2.北京卫星环境工程研究所：北京100094）摘要：为研究不锈钢板式热沉的结构参数对其流动传热性能的影响，进一步提升其换热效率，利用ANSYS 软件建立了与真实结构较为一致而流道高度、焊点直径及焊点排布方式各不相同的热沉几何模型，并利用 Fluent 软件进行数值模拟。结果表明：在流速相等的条件下，较大的流道高度 H 和较小的焊点直径 D 会使热沉的流阻下降，温度均匀性提升，从而改善热沉的传热性能；当供液压力足够时，对于焊点排布方式分别为矩形和菱形的热沉，尽管后者的流阻相对更大，但同时也拥有更高的温度均匀性。以上研究结果可为不锈钢板式热沉结构的设计与优化提供参考。关键词：板式热沉；结构参数；流动传热；焊点排布；数值模拟中图分类号：V416.5文献标志码：A文章编号：1673-1379(2023)04-0374-06DOI:10.12126/see.2023023Simulation of the effect of structural parameters on flow andheat transfer of stainless steel plate heat sinkHEHonghui1,GONGJie2*,ZHOUYing2,YANGXue2,JINWeian2(1.ChinaAcademyofSpaceTechnology;2.BeijingInstituteofSpacecraftEnvironmentEngineering:Beijing100094,China)Abstract:Inordertostudytheeffectofstructuralparametersofthestainlesssteelplateheatsinkonitsflowand heat transfer performance to improve the heat exchange efficiency,the software ANSYS was used toestablishageometricmodeloftheheatsinkthatwasconsistentwiththerealstructurebutwithdifferentflowpathheight,solderjointdiameterandsolderjointarrangement.ThesoftwareFluentwasusedtoconductthenumericalsimulation.Theresultsshowthat,undertheequalflowvelocity,largerflowpathheightHandsmallersolderjointdiameterDwouldreducetheflowresistanceoftheheatsinkandimprovethetemperatureuniformity,therebyimprovingtheheattransferoftheheatsink.Whentheliquidsupplypressureissufficient,fortheheatsinks whose solder joint arrangements are rectangular and rhombic,although the latter s flow resistance isrelativelylarger,ithadhighertemperatureuniformity.Theresearchmayprovideareferenceforthestructuraldesignandoptimizationofthestainlesssteelplateheatsink.Keywords:plate heat sink;structural parameters;flow and heat transfer;solder joint arrangement;numericalsimulation收稿日期：2023-03-02；修回日期：2023-07-17引用格式：何鸿辉,龚洁,周盈,等.结构参数对不锈钢板式热沉流动传热性能影响的仿真研究J.航天器环境工程,2023,40(4):374-379HE H H,GONG J,ZHOU Y,et al.Simulation of the effect of structural parameters on flow and heat transfer of stainless steel plateheat sinkJ.Spacecraft Environment Engineering,2023,40(4):374-379Vol.40,No.4航天器环境工程第40卷第4期374SPACECRAFTENVIRONMENTENGINEERING2023年8月https:/E-mail:Tel:[phone] 引言引言热沉是空间热环境模拟设备的核心部件，用于模拟冷黑环境，其换热性能的优劣直接影响地面环境模拟试验的效果1。板式热沉是将 2 层不锈钢板边界进行激光焊接，然后按一定规律在板间点焊，最后经充压胀形等工序使不锈钢板片胀起并在焊点处形成蜂窝状凹面构造出工质流道，从而实现热沉表面直接与工质的换热2-5。蜂窝结构不仅改变了热沉流道截面积，也兼具扰流作用，有效提升了热沉的换热效率和承受热负荷的能力6-7。张磊等4,6研究了工质流速对板式热沉换热的影响，发现增加流速使热沉换热效率提高的同时也将增大流体压力损失；同时发现等流量下热沉的流道深度与焊点间距是影响其换热性能的重要因素，需选择合理的流道结构设计参数。孙华强8和单巍巍等9对比分析了液氮在传统管板式热沉与不锈钢板式热沉流道中的传热效果和阻力特性，发现后者的换热效率及温度均匀性皆优于前者。冯宁等10对板式调温热沉壁面温度分布情况进行模拟，发现循环氮气入口速度越大，热沉温度分布越均匀。然而影响板式热沉流动传热性能的因素较多，上述文献均未考虑等流速条件下热沉流道高度、焊点直径与焊点排布方式的影响。为更全面地了解不锈钢板式热沉传热性能的影响因素，本文建立了与真实结构接近而结构参数不同的多个不锈钢板式热沉几何模型，并利用 Fluent软件对热沉内工质的流动传热进行数值模拟，比较不同流道高度、焊点直径以及焊点排布方式下热沉传热性能的优劣，以期为热沉结构的优化设计提供参考。1 数值模拟模型数值模拟模型1.1几何模型本文研究的不锈钢板式热沉几何结构参数主要包括流道高度 H、焊点直径 D 和焊点排布方式，如图 1 所示。板式热沉的结构比较复杂，目前已有的数值模拟研究对其几何模型均有不同程度的简化。为了得到与板式热沉实际结构较为一致的几何模型，本文利用 ANSYS 软件进行建模，根据热沉的成形特点将胀形过程视为准静态过程，同时考虑材料回弹，并忽略胀形速度、加速度及进出口位置等对不锈钢板应力分布的影响，主要设置包括：1）本构关系模型采用考虑弹性问题冷加工的线性强化模型；2）屈服准则采用考虑中间主应力影响的 Von-Mises 准则；3）有限元计算方法采用静力隐式算法，模型求解采用 Newton-Raphson 迭代法；4）对板件内表面施加垂直向外方向的均匀压力，开启大变形开关。DH(b)焊点矩形排布(c)焊点菱形排布(a)流道高度 H与焊点直径 D图1板式热沉的几何结构参数Fig.1Geometricstructuralparametersoftheplateheatsink通过上述建模步骤得到的板式热沉几何结构单元如图 2 所示；同时，以流道高度 7mm、焊点间距 75mm、焊点直径 14mm 的热沉实物为例，对其两相邻焊点间的胀起高度进行多次测量，得到实物与仿真结果的胀起高度对比（见图 3），可以看到三维仿真模型与实物结构有着较高的相合性。图2板式热沉的几何结构单元Fig.2Geometricstructuralunitoftheplateheatsink9876543210-1-2胀板高度/mm水平方向坐标/mm0510 15 20 25 30 35 40 45 50 55 60 65 70 75仿真计算实物测量图3热沉实物与模型的胀起高度对比Fig.3Comparisonofbulgeheightbetweentheheatsinkanditsmodel板式热沉模型的平面示意见图 4，长 1200mm、宽 400mm，流体自左侧流入、右侧流出。为研究热沉充分发展段的流动换热规律，避免进出口结构的影响，选择在长度方向上多建立几排焊点，以区分第4期何鸿辉等：结构参数对不锈钢板式热沉流动传热性能影响的仿真研究375进出口对焊点的影响段与充分发展段，同时在最后一排焊点后设置 400mm 的无焊点延长段。1200400流入流出图4模型平面示意图Fig.4Planediagramofthemodel1.2数学模型流体流动所遵循的控制方程包括质量守恒方程、动量守恒方程及能量守恒方程11。质量守恒方程为t+div(u)=0，(1)式中：为流体密度，kg/m3；t 为时间，s；u 为流体速度矢量，m/s。动量守恒方程为t(u)+div(uu)=-px+div(gradu)+Sut(v)+div(vu)=-py+div(gradv)+Svt(w)+div(wu)=-pz+div(gradw)+Sw，(2)式中：u、v、w 分别为速度矢量在 x、y、z 方向上的分量，m/s；p 为流体静压，Pa；为流体分子黏度，m2/s；Su、Sv、Sw为动量方程的 3 个广义源项。能量守恒方程为t(T)+div(Tu)=div(cpgradT)+ST，(3)式中：T 为流体温度，K；为流体导热系数，W/(mK)；cp为流体比热容，J/(kgK)；ST为黏性耗散项。1.3边界条件与数值求解方法根据板式热沉的结构特点，对几何模型进行非结构化网格划分和网格无关性验证后，将网格文件导入 Fluent 软件进行仿真计算。计算边界条件设置如下：1）入口边界条件采用速度入口，入口温度设为300K，入口流体速度为 1m/s；2）出口边界条件采用压力出口，环境压力设为标准大气压；3）上、下壁面和焊点周围及侧边焊缝均设为无滑移边界条件，下壁面温度分布服从恒热流条件，热流密度设为 6000W/m2，其他壁面绝热。热沉内工质选择水，湍流模型选用 RNGk-模型12，近壁面处理选用增强壁面函数，压力速度耦合选用 Coupled 算法，各参数的离散均为二阶精度迎风格式，解收敛判断标准是相对残差不大于110-5。2 仿真结果与讨论仿真结果与讨论2.1流道高度的影响对流道高度为 4mm、7mm、10mm，焊点间距为 75mm，焊点直径为 14mm，焊点排布方式为矩形的不锈钢板式热沉进行模拟计算，并从每排焊点圆心所在截面，绘制不同流道高度下截面平均压力随截面位置变化的曲线（图 5），同时得到热沉的温度分布云图（图 6）及不同截面平均压力分布云图（图 7）。00.20.40.60.81.01.2010002000300040005000600070008000截面平均圧力/Pa截面与入口面的距离/mH=10 mmH=7 mmH=4 mm图5不同流道高度下截面平均压力随截面位置的变化Fig.5Variationofaveragepressureofthecross-sectionwithitspositionofthecross-sectionfordifferentflowpathheights309.9308.9308.0307.0306.0305.0304.0303.0302.1301.1300.1温度/K(a)H=4 mm温度/K308.9308.0307.2306.3305.4304.5303.6302.7301.0301.9300.1(b)H=7 mm温度/K(c)H=mm376航天器环境工程第40卷温度/K(a)H=mm温度/K(b)H=mm温度/K309.8308.8307.8306.9305.9304.9304.0303.0302.0301.1300.1(c)H=10 mm图6不同流道高度下的热沉温度分布云图Fig.6Heatsinktemperaturedistributioncontourfordifferentflowpathheights8626.27763.26900.25174.24311.33448.32585.31722.3859.3-3.66037.2压力/Pa(a)H=4 mm压力/Pa5587.75028.24468.83909.43349.92790.52231.11671.61112.2552.8-6.7(b)H=7 mm压力/Pa4490.94040.93590.93140.92690.92240.91790.91340.9890.9440.9-9.1(c)H=10 mm图7不同流道高度下的热沉压力分布云图Fig.7Heat sink pressure distribution contour for differentflowpathheights由图 5 可知，工质经过前 2 排焊点后进入充分发展阶段，且由焊点段过渡到无焊点段时的压力变化并不均匀，对比焊点段和无焊点段可以明显看出焊点会造成更大的压降（流阻）。由图 6 和图 7 可知：焊点附近的温度与压力梯度变化均十分明显；在相同流速下，当流道高度 H 处于 410mm之间，H 越大则压降越小，且换热面温差也将减小，表明此时热沉表面整体结构的温度均匀性较好。这是由于 H 的增大使得流道尺寸相应增加，焊点对工质流动的阻碍作用减弱，同时工质与热沉内表面的接触面积变大，从而降低了流阻并使换热更为充分。因此在实际泵驱冷却循环过程中，应尽量选用 H 较大的胀板结构，通过较低压降实现更大的循环流量和温度均匀性，以改善热沉的流动传热性能。2.2焊点直径的影响对焊点直径为 6mm、10mm、14mm，流道高度为 10mm，焊点间距为 75mm，焊点排布方式为矩形的不锈钢板式热沉结构进行模拟计算，绘制不同焊点直径下截面平均压力随截面位置变化的曲线（图 8），同时得到热沉温度分布云图（图 9）及不同截面平均压力分布云图（图 10）。00.20.40.60.81.01.20500100015002000250030003500截面平均圧力/Pa截面与入口面的距离/mD=14 mmD=10 mmD=6 mm图8不同焊点直径下截面平均压力随截面位置的变化Fig.8Variationofaveragepressureofthecross-sectionwithitspositionfordifferentsolderjointdiameters306.3305.7305.1304.4303.8303.2302.6302.0301.3300.7300.1温度/K(a)D=6 mm307.2306.5305.8305.0304.3303.6302.9302.2301.5300.8300.1温度/K(b)D=10 mm309.8308.8307.8306.9305.9304.9304.0303.0302.0301.1300.1温度/K(c)D=14 mm图9不同焊点直径下的热沉温度分布云图Fig.9Heatsinktemperaturedistributioncontourfordifferentsolderjointdiameters3371.53033.42695.42357.42019.41681.41343.41005.4667.3329.3-8.7压力/Pa(a)D=6 mm-6.2压力/Pa(b)D=mm-9.1压力/Pa(c)D=mm第4期何鸿辉等：结构参数对不锈钢板式热沉流动传热性能影响的仿真研究377-8.7压力/Pa(a)D=mm3851.53465.73079.92694.22308.41922.61536.91151.1765.3379.6-6.2压力/Pa(b)D=10 mm4490.94040.93590.93140.92690.92240.91790.91340.9890.9440.9-9.1压力/Pa(c)D=14 mm图10不同焊点直径下的热沉压力分布云图Fig.10Heatsinkpressuredistributioncontourfordifferentsolderjointdiameters由图 8图 10 可知，工质经过前 3 排焊点后进入充分发展阶段，对比 3 种焊点直径下的压力变化可以发现，当焊点直径 D 处于 614mm 之间，D 较大时，压降有所上升，且换热面温差增大，即此时热沉的温度均匀性较差。这是由于焊点尺寸增大将使热沉流道截面积缩小、扰流作用增强，工质流动时的压力损失增大；另一方面，焊点内部没有流体冷却，只能通过不锈钢导热，而不锈钢的导热性能相对较差，故使焊点位置局部温度较高。因此，实际加工时，应在结构强度允许的前提下，尽量选用焊点尺寸较小的胀板结构。2.3焊点排布方式的影响对焊点排布方式分别为矩形和菱形，流道高度为 7mm，焊点间距为 75mm，焊点直径为 14mm的不锈钢板式热沉结构进行模拟计算，绘制不同焊点排布方式下截面平均压力随截面位置变化的曲线（图 11），同时得到热沉的温度分布云图（图 12）及不同截面平均压力分布云图（图 13）。00.20.40.60.81.01.201000200030004000500060007000矩形排布菱形排布截面与入口面的距离/m截面平均圧力/Pa图11不同焊点排布方式下截面平均压力随截面位置的变化Fig.11Variation of average pressure of the cross-sectionwithitspositionfordifferentsolderjointarrangements308.9308.0307.2306.3305.4304.5303.6302.7301.9301.0300.1温度/K(a)矩形排布308.1307.3306.5305.7304.9304.1303.3302.5301.7300.9300.1温度/K(b)菱形排布图12不同焊点排布方式下的热沉温度分布云图Fig.12Heat sink temperature distribution contour fordifferentsolderjointarrangements5587.75028.24468.83909.43349.42790.52231.11671.61112.2552.8-6.7压力/Pa(a)矩形排布7741.87011.46181.05350.64520.13689.72859.32028.91198.5368.1-462.4压力/Pa(b)菱形排布图13不同焊点排布方式下的热沉压力分布云图Fig.13Heatsinkpressuredistributioncontourfordifferentsolderjointarrangements由图 11图 13 可知，在相同流速下，相比焊点矩形排布，菱形排布会造成更大的压力损失，但换热面温差也更小，即焊点菱形排布可以改善热沉的传热性能，提升热沉表面的温度均匀性。这是由于菱形交错排列的焊点使工质流动的速度大小和方向改变得更为明显，增强了其扰动程度，减少了流动死区，增大压力损失的同时亦强化了传热。因此，在系统供液压力足够的前提下，可尽量选用菱形焊点排布结构以提高板式热沉的温度均匀性。3 结束语结束语本文利用 ANSYS 软件建立了与不锈钢板式热沉实际结构相近的三维模型，并使用 Fluent 软件进行数值模拟，分析了等流速时流道高度 H、焊点直378航天器环境工程第40卷径 D 以及焊点排布方式这 3 种结构参数对热沉流动传热性能的影响。结果发现：在结构强度允许的前提下，应尽量选用 H 较大、D 较小的胀板结构来降低热沉工质流阻，提升热沉表面的温度均匀性；在系统供液压力足够的前提下，选用焊点为菱形排布的胀板结构可改善热沉的传热性能。以上研究结果可为不锈钢板式热沉结构的设计与优化提供理论参考。参考文献（References）黄本诚,童靖宇.空间环境工程学M.北京:中国科学技术出版社,2010:124-1271黄本诚,刘波涛,李志胜.空间模拟器设计技术M.北京:中国宇航出版社,2017:154-1552葛佩红,罗清海,涂敏,等.蜂窝换热板研究进展J.制冷与空调,2019,19(12):5-10GEPH,LUOQH,TUM,etal.ResearchadvanceofhoneycombheatexchangeplateJ.RefrigerationandAir-conditioning,2019,19(12):5-103张磊,刘敏,刘波涛.流速及进出液口形式对板式热沉换热性能影响J.航天器环境工程,2012,29(5):566-570ZHANGL,LIUM,LIUBT.Effectsofvelocitiesandinlet-outlet forms on heat transfer of plate-style heatsinkJ.SpacecraftEnvironmentEngineering,2012,29(5):566-5704蒋志广,许贞龙,张立建.不锈钢板式热沉液压胀形工艺参数仿真分析J.航天器环境工程,2022,39(1):33-39JIANGZG,XUZL,ZHANGLJ.Simulationofhydro-bulging process of stainless-steel plate heat sinkJ.5SpacecraftEnvironmentEngineering,2022,39(1):33-39张磊,刘然,刘敏.蜂窝流道热沉强化传热数值模拟J.装备环境工程,2013,10(4):121-125ZHANGL,LIUR,LIUM.NumericalsimulationofheattransferenhancementinhoneycombshroudJ.EquipmentEnvironmentEngineering,2013,10(4):121-1256金玮安,张磊.大热流背景下蜂窝板结构参数的多目标优化J.航天器环境工程,2023,40(2):128-133JIN W A,ZHANG L.Multi-objective optimization ofstructuralparametersofhoneycombplateunderlargeheatflowJ.SpacecraftEnvironmentEngineering,2023,40(2):128-1337孙华强.小型热真空环境模拟装置热沉的热力学分析D.沈阳:东北大学,2013:59-708单巍巍,丁文静,王紫娟,等.不锈钢胀板式热沉性能分析及其数值模拟C/第十届全国低温工程大会暨中国航天低温专业信息网 2011 年度学术交流会论文集,2011:557-5619冯宁,刘智勇,颜昌林,等.流速对鱼骨式与板式调温热沉性能的影响J.兰州交通大学学报,2022,41(5):115-122FENG N,LIU Z Y,YAN C L,et al.Effect of flowvelocity on the performance of fishbone type and platetype temperature regulating heat sinkJ.Journal ofLanzhouJiaotongUniversity,2022,41(5):115-12210陶文铨.传热学M.5 版.北京:高等教育出版社,2019:34-3811YAKHOT V,ORZAG S A.Renormalization groupanalysis of turbulence:basic theoryJ.Journal ofScientificComputation,1996(1):3-1112（编辑：冯露漪）一作简介：何鸿辉，高级工程师，主要从事空间环境模拟设备研制等工作。*通信作者：龚洁，研究员，主要从事空间环境模拟技术的研究。第4期何鸿辉等：结构参数对不锈钢板式热沉流动传热性能影响的仿真研究379