茹卡乌斯卡斯横掠顺排管束的数值模拟与公式验证分析

茹卡乌斯卡斯横掠顺排管束的数值模拟与公式验证摘要：茹卡乌斯卡斯实验关联式在换热器的设计中有着广泛的应用。针对茹卡乌斯卡斯研究流体横掠顺排管束流动与换热特性的实验段为原型，经过适当的简化，建立二维模型，运用大型CFD软件Fluent对该模型内流体的流动与换热特性进行了数值模拟研究。将数值模拟结果与公式计算进行对比，并考察了因素进口速度、进口温度、换热管壁温及换热管外径对换热性能（努塞尔数）和使用茹卡乌斯卡斯公式进行计算的误差的影响。关键词：数值模拟、横掠顺排管束、公式验证、Nu数Abstract:Zhukauskasexperimentalcorrelationhasbeenwidelyusedinheatexchangerdesign.FluidforZhukauskasACROSSLINEBUNDLEtestsectionflowandheattransfercharacteristicsoftheprototype,afterappropriatesimplification,toestablishatwo-dimensionalmodel,usingtheFluentCFDsoftwareinlargemodelfluidflowandheattransfercharacteristicsofnumericalsimulation.Thenumericalresultsarecomparedwiththeformula,andexaminethefactorsofinletvelocity,inlettemperature,theheattransfertubewalltemperatureandtheoutsidediameteroftheheattransferperformanceofheattransfertubes(Nusseltnumber)anduseZhukauskasformulaconductingimpacterrorofcalculation.Keywords:numericalsimulation,ACROSSLINEBUNDLEformula,Nunumber流体外掠管束的流动比较复杂,其换热和流动特性主要采用实验数据和实验关联式进行预测。在众多应用于换热器设计的实验关联式中,茹卡乌斯卡斯公式有着非常广泛的应用[1]。本文研究Re数在65102102，流体外掠顺排管束为20排的实验模型，根据计算流体力学理论，运用大型通用流体力学计算软件Fluent对实验模型内流体的换热和流动特性进行数值模拟，并将数值模拟结果与公式25.0)Pr/(PrPrRe033.036.08.0wffffNu计算结果进行误差分析。1模型描述1.1模型与参数模型为干空气横掠顺排管束，示意图见图1；几何尺寸及操作条件见表1。图1模型正面示意图表1几何尺寸及操作条件1.2控制方程本问题所涉及的流动为二维稳态定常流动,采用SIMPLE方法求解N-S方程，湍流模型采用标准k-ε模型。[1-4]方程式如下[2]-[5]：连续性方程0xujj动量方程)]xuxu(xu[xxpx)uu(t)u(ijjijijijjiik方程)]xuxu(xuxk)[(xx)ku(t)k(ijjijitjktjjjε方程kc)xuxu(xukc]x)[(xx)u(t)(22ijjijit1jtjjj能量方程]xT)c[(xx)Tuc(t)Tc(jttpjjjpp其中2tkc式中：x为坐标；u为速度分量；i，j为张量指标；t为时间；k为湍流脉动动能；ε为脉动动能耗散率；T为温度；ρ为流体密度；μ为动力粘度；cｐ为定压比热容；μｔ为湍流粘度；σ几何尺寸操作条件干空气50℃物性参数模型长度L/mm520进口速度V/(m/s)200密度ρ/(kg/m3)1.093模型宽度H/mm104进口温度T/℃50热容cp/(J/(kg*K))1005换热管外径D/mm20管壁温度Tw/℃20热导率λ/(W/(m*K))0.0283中心距S1、S2/mm26流动介质空气动力粘度v106/(m2/s)17.95ｔ为能量普朗特数，一般取0.85；σk为脉动动能普朗特数，一般取1.0；σε为脉动动能耗散率普朗特数，一般取1.3；模型系数cμ，c1和c2分别取值0.09，1.44和1.92。1.3网格划分对计算区域采用非结构化壳网格划分，以四边形网格为主，允许一部分三角形网格单元的存在；在为了准确的计算流体在管束壁面边界层内的变化，在管子外划分4层的边界层，底层网格高度为0.5mm，增长率为1.1。全局网格如图2所示，部分网格如图3所示。图2全局网格显示图3部分网格显示1.4边界条件模型计算边界条件设定如下：管束壁面温度为恒温20℃；流体进口温度为50℃；流体入口采用速度进口；出口采用自由流动；为避免两侧壁面的对体流动产生影响，对轴向中间截面采用对称边界。1.5计算方法在进行模拟时,湍流模型选用标准k-ε模型，壁面处流体无滑移。控制方程采用有限体积法离散，压力速度耦合采用SIMPLE算法,压力插补格式选择standard格式,其他项采用一阶迎风格式[6]。2结果与讨论2.1网格独立性分析为了检查数值计算的准确性，分别对网格数目为5380、8804、16061、36541、49499进行独立性检验。模拟计算结果如图4所示。从图中可看出当网格数达到36541时，再增大网格数对Nu数影响已经很小了，考虑到计算的准确性和花费时间，最后选取36541作为计算网格数目[7]-[8]。0100002000030000400005000015001600170018001900200021002200Nu数网格数图4Nu与网格数之间的关系2.2与经验关联式的对比流体横掠顺排管束（≥16）的平均表面换热系数由茹卡乌斯卡斯的如下公式计算：25.0)Pr/(PrPrRe033.036.08.0wffffNu（5102Re6102）由fluent模拟结果：出口温度取右侧出口处的温度为35.16℃，截面最小处速度取g1、g2、g3、g4截面处的最大速度为862.81m/s。当时，Prf=0.702。计算雷诺数为：tw=50°C，Prw=0.703；由茹卡乌斯卡斯关联式得努赛尔数：计算质量流率对数平均温差由能量平衡方程，计算表面传热系数有误差561061.91095.1702.0200ReudCtttfff58.42216.35502'97.1770)699.0703.0(703.0)1061.9(033.025.036.08.06fNu)/(73.22200093.1104.0104.0smkguqm74.2116.35205020ln5016.35mt)/(77.3103)74.21(02.014.3204)5016.35(100573.22)'(2KmWtAttcqhmffpm66.21931024.2802.077.31033fhdNu193.066.219397.177066.2193e2.3模拟结果的云图分析利用fluent14.0进行数值模拟计算，计算域温度变化如图5所示；计算域压力变化如图6所示；计算域速度变化如图7所示。由温度云图可知，流体温度从入口到出口逐渐降低，是由于流体与管壁进行了热量交换；由压力云图可知，流体压力从入口到出口逐渐降低；由速度云图可知，流体速度从入口到出口逐渐降低；由局部速度矢量图可知，流体在绕管束流动时会在近壁面产生漩涡。图5温度云图图6压力云图图7速度云图图8局部速度矢量图2.4不同影响因素的分析为验证公式的适用性，考察因素进口速度、进口温度、管壁温度，换热管外径对指标误差的影响。因素一进口速度变化时，进口速度与模拟结果Nu数、公式计算Nu数关系如图9所示；因素二进口温度变化时，进口温度与模拟结果Nu数、公式计算Nu数关系如图10所示；因素三管壁温度变化时，管壁温度与模拟结果Nu数、公式计算Nu数关系如图11所示；因素四换热管外径变化时，换热管外径与模拟结果Nu数、公式计算Nu数关系如图12所示。由图9可知Nu数随着进口速度的增大而增大，进口速度在60m/s时公式计算与数值模拟的Nu相比误差为9.48％，误差随着速度的增大而降低；由图10可知Nu数随着进口温度的增大而增大，误差基本在19％左右；由图11可知管壁温度对Nu数几乎没有影响，误差处于19％的水平；由图12可知Nu数随着换热管外径的增大而减小，换热管外径在12mm是误差较小，为11.64％，误差随着换热管的外径增大而增大。0501001502002503003504005001000150020002500300035004000Nu数进口速度m/s数值模拟公式验证5060708090150016501800195021002250Nu数进口温度/℃数值模拟公式验证图9进口速度对Nu数的影响图10进口温度对Nu数的影响-20-100102010001200140016001800200022002400260028003000Nu数管壁温度/℃数值模拟公式验证121620242832350030002500200015001000Nu数换热管外径/mm数值模拟公式验证图11管壁温度对Nu数的影响图12换热管外径对Nu数的影响2.5误差分析从软件数值模拟与茹卡乌斯卡斯实验关联式比较，知数值模拟会产生一定的误差，分析产生误差的原因如下：1）湍流模型的选择以及湍流强度的计算；2）FLUENT模拟不出管壁的粗糙度；3）求解方法的选择，一阶迎风格式或者二阶迎风格式。3结论应用Fluent软件,以茹卡乌斯卡斯研究横掠管束的模型进行数值模拟研究。将数值模拟结果与茹卡乌斯卡斯实验关联式计算的结果进行对比，表明数值模拟产生的误差是可控的。通过单因素分析得到如下结论：1）1）将三维模型简化为二维模型大大减少计算时间，将两侧设为对称边界条件，减少壁面对流体流动状态产生影响，采用对称边界获得的数值结果是可靠的；2）在四个因素中，进口温度对误差影响较大；3）管壁温度对Nu数无影响；4）进口速度越大，换热效果越好；换热管外径越大，换热效果越差。参考文献：[1]董其伍,欧阳克,刘敏珊,张丽娜.茹卡乌斯卡斯横掠错排管束实验模型的数值模拟[J].压力容器,2010,01:21-26.[2]余徐飞,王治云,李起耘,黄文和,陆廷安,杨茉.横掠周期性密集管束流动换热的数值模拟[J].上海理工大学学报,2015,06:563-567+576.[3]刘敏珊,杨帆,董其伍,谢建.流体横掠管束模拟中壁面函数影响研究[J].热能动力工程,2010,05:497-500+573.[4]杨茉,赵明,叶剑军,章立新,余敏,卢玫.非稳态横掠管束周期性充分发展对流换热的数值模拟[J].工程热物理学报,2003,06:988-991.[5]邓斌,陶文铨.管壳式换热器壳侧湍流流动的数值模拟及实验研究[J].西安交通大学学报,2003,09:889-893+924.[6]邓斌,陶文铨.管壳式换热器壳侧湍流流动与换热的三维数值模拟[J].化工学报,2004,07:1053-1059.[7]潘维,池作和,斯东波,岑可法.匀速流体横掠管束的流场数值模拟[J].浙江大学学报(工学版),2004,08:110-113.[8]康细洋.格子Boltzmann方法研究外掠圆管及管束的流动和换热特性[D].重庆大学,2009.