90°弯管内流动的理论模型及流动特性的数值研究110精.docVIP

第21卷第3期 2004年6月　计算力学学报　 ChineseJournalofComputationalMechanics Vol.21,No.3 June2004 文章编号:100724708(2004)0320314208 90°弯管内流动的理论模型及流动特性的数值研究 丁　珏,　翁培奋3 (上海大学应用数学和力学研究所,上海200072) 摘　要:从三维不可压缩雷诺时均Navier2Stokes方程出发,对90°弯曲管道内湍流流动进行数值模拟。网格划分采用六面体网格,湍流模型为RNGk2Ε模型,在近壁区采用两层壁面模型进行修正,流场的计算结果与实验数据吻合较好。在此基础上,本文数值研究了来流方向对流场结构和流动特性的影响。得出在弯管流场中发生了分离现象,且随着来流侧滑角的增大,分离区范围增大。此外,随着来流从同一侧滑角变换至同一攻角时,横截面的二次流图像中也从具有两个对称主涡变成只具有一个主涡的现象。关键词:弯曲管道;Navier2Stokes方程;湍流模型;流动特性中图分类号:O357.5　　　文献标识码:A 1　引　言 弯曲管道广泛地应用于工业、农业等机械设备上,如压缩机、泵及各种类型的热交换器,甚至核动力的管道系统中。这些形形色色的弯管实现了流体输运和热量、质量交换等功能。一般情况下,影响流体流动性质的因素有很多,诸如弯管的弯曲程度,流体的来流马赫数Ma,流体运动方向等。在纵多因素的影响下,弯曲管道内的流场呈现出十分复杂的流动特性。如一定来流条件下,在管壁附近形成分离区,管道横截面上产生二次流动,这些现象不仅造成流体总压和能量的损失,而且形成的局部障碍区域也使流动系统的阻力增大,降低了热量、质量的交换效率。因此,弯管内的流动一直受到内流研究者的关注 [123] 较大成绩。史峰[4]等人利用标准k2Ε湍流模型研究了Rc 弯道内的湍流流动。计算发D=1.69的90°现利用k2Ε模型得到的流场时均速度分布与实验数据相比,在弯段Η=45°截面前,二者符合较好,即计算误差较小;当弯段Η45°,存在较大的误差。文献[5]采用湍流大涡模型,对90°弯管进行了数值研究,拓宽了湍流大涡模型在具有强曲率弯曲管道内流体流动问题的应用。 基于以往研究的基础,本文引进Yahhot和Orszag [7] 新近应用重整化群方法提出的RNGk2Ε 湍流模型,建立90°大曲率弯管内部流场的Navier2Stokes方程,应用交错网格系统下的SIMPLE算 法求解,并将计算结果与文献[1]中的实验数据进行对比。此外,针对以往研究较少的流体运动方向对流场结构及流动特性影响,本文着重进行了分析。给出了几种来流侧滑角和攻角下流场内重要参量的分布情况,揭示出流场中出现的物理现象和规律。 。 早期的实验研究以Taylor为代表,他采用激光多谱勒测速仪对90°方形截面弯管内雷诺数分别为790和40000的来流条件进行了系统的研究,给出了层流和湍流两种流态下的时均速度分布,以及弯曲段壁面的压力分布等重要实验结果[1]。文献[2]应用LDV测量了弯曲度为Rc D=1.69(Rc,D分别代表曲率的平均半径与弯管的水动力直径)的矩形截面弯道内的湍流。而文献[3]则通过热线测速系统对三维90°弯管内流场的气流特性进行了实验研究。在理论研究方面[426],国内研究者已取得 收稿日期:2002209229;修改稿收到日期:20032032311 基金项目:教育部青年教师教学科研奖励基金;上海市曙光 计划项目及上海市重点学科资助项目1 作者简介:丁　珏(19732),女,副研究员; 翁培奋3(19652),男,教授,博士生导师1 2　数学模型 2.1　控制方程 通过对连续方程和瞬间Navier2Stokes方程时均化 ,得到直角坐标系下定常条件、不可压缩流体流动遵循的方程。　　连续方程 =0　(i=1,2,3;j=1,2,3)xi - (1) 　　N2S方程 ? 1995-2006 Tsinghua Tongfang Optical Disc Co., Ltd. All rights reserved. 　第3期 -- 丁　珏,等:90°弯管内流动的理论模型及流动特性的数值研究 - 31 5 --′(Θuj=-+-Θu′gΦΛΘiuj)-i2 xjxixjxj - 　　　5=u,v,w,k,Ε 其中S5,#5分别代表通用变量5的源项和有效扩散系数。J是Jacobi变换矩阵。u,v和w分别代表曲线坐标下的速度分量。 Φi2= -- 1,0, i=2i≠2 (2) 其中　ui,p分别代表流体的平均流速和平均压力。雷诺应力张量-Θuiuj使用由Boussinesgue等提出的涡粘性模型来计算。2.2　RNGk-Ε湍流模式[7] RNGk2Ε湍流模型既适应于高雷诺数情