扇形气膜孔冷却效率的数值模拟
- 格式:pdf
- 大小:4.85 MB
- 文档页数:5
前缘气膜孔对涡轮静叶冷却效果影响的数值模拟
杨凡;郑洪涛;李智明
【期刊名称】《热能动力工程》
【年(卷),期】2006(21)4
【摘要】采用全三维数值模拟技术,利用k-ε双方程湍流模型和SIMPLE算法,通过求解三维粘性可压缩Favre平均Navier-Stokes方程,对某新型燃气轮机第一级气膜冷却静叶叶栅的三维湍流流场进行了数值模拟。
分别通过改变燃气轮机前缘气膜孔的参数,计算出叶片外表面的温度分布和冷却空气流量大小。
结果表明,前缘气膜孔的直径、数目以及射流方向对叶片表面冷却效果的影响是非常显著的。
从而提出了一种叶片前缘气膜冷却设计的新方案,为工程设计提供了有价值的参考。
【总页数】5页(P345-349)
【关键词】燃气轮机;前缘气膜冷却;数值模拟;涡轮;第一级静叶
【作者】杨凡;郑洪涛;李智明
【作者单位】哈尔滨工程大学动力与核能工程学院
【正文语种】中文
【中图分类】TK474.7
【相关文献】
1.涡轮静叶复合角度气膜冷却孔排布置优化的数值研究 [J], 张玲;汪山入;董海瑞;张宏洋
2.叶栅前缘单排冷却孔气膜冷却效果的数值研究 [J], 袁野;万剑峰
3.涡轮静叶前缘气膜冷却数值模拟 [J], 杨凡;曹辉;郑洪涛;李智明
4.涡轮叶栅前缘气膜冷却数值模拟 [J], 颜培刚;王松涛;韩万金;王仲奇
5.涡轮叶栅前缘槽缝气膜冷却的数值模拟 [J], 王晓东;康顺
因版权原因,仅展示原文概要,查看原文内容请购买。
旋转对复合角度气膜冷却叶片的数值模拟
李少华;李知骏;王梅丽;郭婷婷
【期刊名称】《汽轮机技术》
【年(卷),期】2011(053)003
【摘要】采用Realizable k-ε紊流模型并结合SIMPLEC算法,对前缘复合角度α=30°、β=45°,α=90°、β=45°的动叶栅在不同旋转状速度下的气膜冷却效率进行计算.分析了不同转速、吹风比、叶片前缘射流角度对气膜冷却效率的影响.计算结果表明:旋转导致冷却射流向叶顶偏移,转速越高气膜冷却效率越低;高转速时叶盆区域有回流涡旋形成;高吹风比使得冷却射流在吸力面的贴壁性变差;比较两种前缘冷气喷射角度的计算结果可以看出,前缘冷却气流喷射角度较小时的气膜冷却效果较好.
【总页数】3页(P164-166)
【作者】李少华;李知骏;王梅丽;郭婷婷
【作者单位】东北电力大学,吉林132012;东北电力大学,吉林132012;东北电力大学,吉林132012;北京国电龙源环保工程有限公司,北京102206
【正文语种】中文
【中图分类】TK262
【相关文献】
1.旋转对涡轮叶片气膜冷却影响的数值模拟 [J], 袁锋;吴亚东;竺晓程;杜朝辉
2.旋流燃烧器内二次风叶片旋转角度对流场影响的数值模拟 [J], 车娟
3.旋转状态下曲率对叶片气膜冷却特性的影响 [J], 张玲;张璐琦;白博升
4.燃气轮机旋转状态下的动叶气膜冷却效果数值模拟研究 [J], 李录平;唐学智;张浩;黄章俊
5.旋转涡轮叶片端部气膜冷却的数值计算 [J], 杨汇涛;曹玉璋
因版权原因,仅展示原文概要,查看原文内容请购买。
连续收缩扩张气膜孔排冷却特性的数值模拟周骏飞;王新军;张峰【摘要】为研究连续收缩扩张孔的冷却特性,在C3X静叶片上分别建立了连续收缩扩张气膜孔冷却模型、圆柱气膜孔冷却模型和展向扩张气膜孔冷却模型,连续收缩扩张气膜孔每排23个、孔间距为20 mm,展向扩张孔每排19个、孔间距为24 mm,圆柱孔每排19个、孔间距为24 mm.同时,在叶片前部开设了一个U形冷却通道,尾部开设了一个直冷却通道,冷气通过这2个内部冷却通道进入气膜孔.利用ANSYS-ICEM商用软件对上述3种模型进行了结构化网格划分,采用ANSYS-CFX 商用软件和SST湍流模型进行了数值计算和分析比较,结果表明:连续收缩扩张孔的气膜冷却效率高于圆柱孔和展向扩张孔,在孔口附近和高吹风比下的优势最明显;连续收缩扩张孔使冷气射流在相邻两孔的交汇处形成了类似反肾形涡结构,该涡的强度不大,但具有良好的延续性和较大的冷气覆盖面积;复合冷却时冷气射流脱离壁面的现象更明显,孔口附近总冷却效率低于绝热冷却效率.在连续收缩扩张孔的实际应用中选择偏大的吹风比和更小的入射角可以提高气膜冷却效率.【期刊名称】《西安交通大学学报》【年(卷),期】2016(050)003【总页数】7页(P38-44)【关键词】收缩扩张孔;气膜冷却;内部冷却;静叶片;数值计算【作者】周骏飞;王新军;张峰【作者单位】西安交通大学能源与动力工程学院,710049,西安;西安交通大学能源与动力工程学院,710049,西安;西安交通大学能源与动力工程学院,710049,西安【正文语种】中文【中图分类】TK263燃气轮机高温部件必须受到冷却保护,气膜冷却是在透平叶片表面形成低温气体保护层,使叶片能够承受更高的燃气温度。
传统圆柱形气膜孔流动损失大,气膜覆盖率低,在高吹风比下容易产生射流而脱离壁面,且无法满足设计者的要求。
Goldstein 率先提出了一种带有展向扩张的成型孔结构[1]。
Bunker综述了4种典型扩张孔结构并进行了系统分析[2],表明与传统的圆柱孔相比,扩张孔的气膜冷却效率更高,气膜覆盖面积更大,且在较高吹风比下有着明显的优势。
文章编号:1000-033X (2017)05-0089-05收稿日期:2016-12-15基金项目:国家科技支撑计划项目(2013BAF07B04);唐山市重点汽车实验室建设项目(12130201A-2)作者简介:刘佳鑫(1983-),男,吉林桦旬人,博士后,讲师,研究方向为工程车辆节能降噪技术.基于CFD 数值仿真的工程机械冷却风扇性能分析刘佳鑫1,2,王宝中1,邢梦龙1,秦四成3,蒋炎坤2,龙海洋1(1.华北理工大学机械工程学院,河北唐山 063009;2.华中科技大学能源与动力工程学院,湖北武汉 430074;3.吉林大学机械科学与工程学院,吉林长春 130022)摘 要:为了提升冷却风扇的性能,使流经散热器空气的流动状态得到改善,基于国内某工程机械用冷却风扇建立几何模型,利用CFD 数值仿真对冷却风扇进行模拟.结果表明:原风扇仿真模型正确;孤立翼型法与合理的特征控制可用于风扇的几何开发与建模;当体积流量达到10.24m 3㊃s -1时,新风扇全压效率提升了约0.81%,轴功率降低了0.12kW .关键词:数值仿真;冷却风扇;风扇设计;工程机械中图分类号:U415.5 文献标志码:BAnal y sis on Performance of Coolin g Fan for ConstructionMachiner y Based on CFD Numerical SimulationLIU Jia-xin 1,2,WANG Bao-zhon g 1,XING Men g -lon g 1,QIN Si-chen g 3,JIANG Yan-kun 2,LONG Hai-y an g 1(1.School of Mechanical En g ineerin g ,North China Universit y of Science and Technolo gy ,Tan g shan 063009,Hebei ,China ;2.School of Ener gy and Power En g ineerin g ,Huazhon g Universit y of Science and Technolo gy ,Wuhan 430074,Hubei ,China ;3.School of Mechanical Science and En g ineerin g ,Jilin Universit y ,Chan g chun 130022,Jilin ,China )Abstract :In order to increase the p erformance of the coolin g fan so as to im p rove the air flow throu g h the radiator ,the g eometric model was established based on the domestic coolin g fan for construction machiner y .The p erformance of the coolin g fan was simulated b y CFD numericalsimulation.The results show that the ori g inal simulation model is correct ;the isolated airfoilmethod and reasonable characteristic control can be used for g eometric develo p ment and modelin g of the fan ;when the volume flow reaches 10.24m 3㊃s -1,full p ressure efficienc y of the new fan increases b y about 0.81%,and shaft p ower reduces b y 0.12kW.Ke y words :numerical simulation ;coolin g fan ;fan desi g n ;construction machiner y0 引 言冷却风扇作为发动机冷却系统的一个重要部件,在为动力舱创造空气流动环境的同时,也为散热器组提供充足的冷却风.通过冷却风扇控制空气流量,可以合理地实现对冷却液温度的调控,保证各部件在合理的温度区间,从而保持发动机高效率的工作状态.由此可见,冷却风扇性能直接影响着设备总体的性能和工作可靠性.针对冷却风扇性能,国内外专家学者已从冷却风扇设计方法入手积极开展研究:上海理工大学王企鲲运用 等密流型 与 变密流型 2种设计方法进行扭叶片改型设计,并利用CFD (Com p utationalFluid D y namics )仿真对比两者差异,且对2种设计方法作出了评价[1];刘全忠通过数值模拟,以汽轮发电机为研究对象,获得了安装角等参数与效率之间的关系,并对其展开性能优化[2];高红霞等使用控制涡设计理论对直升机轴流风扇气动性能进行改进提升;M Henner 利用CFD 数值模拟的方法,通过调整结构设计参数,实现对风扇性能的改进[3];Zhou Jianhui 则采用理论设计与CFD 仿真验证相结合的方式对CPU 冷却风扇进行了重新设计[4].基于以上已取得的研究成果,本文利用CFD 数值仿真对用于某国内工程机械上的冷却风扇进行性能分析,通过试验数据验证仿真模型;利用孤立翼型法对其重新设计,并对新风扇进行性能仿真,对比仿真结果,以确认性能改进和提升的有效性.1 原风扇仿真与试验验证1.1 风扇与风道三维模型根据图纸建立原风扇模型,如图1所示.风扇直径为780mm ;轮毂比为0.33;弯掠角为0ʎ;轮毂直径为260mm ;安装角为26ʎ;翼型为等厚度弯板.图1 原风扇模型建立风道三维几何模型,如图2所示.风道采用圆形截面,直径等于其风扇当量直径,风道的入口长度为4倍当量直径,出口长度为6倍当量直径.在出口风道2~4倍风扇当量直径处建立整流栅.图2 风道模型1.2 网格划分与边界设定利用Gambit 对模型进行网格划分与边界设定.采用结构性网格对风扇表面进行网格划分;对旋转域内部以及进出口处使用结构性与非结构性网格进行网格混合加密,加密后无负网格和扭曲网格,如图3所示.图3 网格模型将风扇㊁风道表面㊁整流栅设定为wall ;将风扇旋转域与风道进出口公共面设定为interior ;将风道进口设定为velocit y -inlet ;风道出口设定为p res-sure-outlet ;将旋转域以及进口风道㊁出口风道设定为fluid ;风扇转速为2000r ㊃min -1;选用标准k-ε湍流模型[5-8].图4 原风扇进口处空气状态1.3 结果分析与试验验证如图4所示,在原风扇空气进口处,低压区主要集中在风扇中上部分,呈阶梯状分布.这是由于扇叶扫掠过的区域空气体积骤降,周围的空气无法迅速填充而造成的;随后在压差作用下,低压区空气体积不断被填充,该区域逐步回升到环境压强;随着转速的增加,空气持续流通,阶梯状分布的低压区逐渐变成环状低压区.图5(a )中,原风扇空气出口端轮毂处出现较大的低压区;高压区则呈环状分布在0.3倍叶高至叶尖之间,并在叶尖处压强达到最大,这是风扇对气流做功的结果.图5(b )中,气流在叶尖处流速最大,并由此逐渐向两边递减;最低气流速度出现在轮毂处,这种现象可能会导致回流的产生.图5 原风扇出口处空气状态在风扇转速一定的条件下,对7组不同流量值进行试验仿真,在5倍当量直径处提取全压值,计算全压效率,与试验数据进行对比,绘制两者全压及全压效率性能曲线,如图6所示.从图6(a )可以看出:当风量降低时,二者偏差逐渐增大,最大误差为4.5%;当风量大于6m 3㊃s -1时,原风扇仿真值与试验值逐渐趋于吻合,最大误差约2.38%.从图6(b )可以看出,随着风量的增加,全压效率都呈现出上升趋势,且均在风量为10.24m 3㊃s -1时出现最大值,原风扇仿真值与试验值最大相差约3%.2 风扇设计与数值仿真2.1 风扇设计选用孤立翼型设计法对冷却风扇进行设计.设计分为3个阶段:选定设计参数㊁预备设计㊁正图6 原风扇试验与仿真性能对比式设计.(1)设计参数.参照原风扇性能和结构参数,设定全压为550Pa ,体积流量为10.24m 3㊃s -1,新风扇采用Clark-y 翼型.(2)预备设计.预备设计计算如下Λ=V a ΩR λb =ΛX bεSb =K thλb2C L δ=2εSb cos βm ìîíïïïïïïïï(1)式中:Λ为叶尖流速系数;V a 为通过叶轮的平均轴向速度;Ω为叶轮叶片的旋转角速度;R 为叶轮半径;λb 为叶片根部流速系数;X b 为轮毂比;εSb 为根部旋流系数;K th 为理论全压系数;C L 为升力系数;δ为实度;βm 为叶片根部合速度与轴向速度的夹角.预备设计筛选依据为:流速系数不大于1.4,旋流系数不大于1,叶片根部载荷因子不大于1.(3)正式设计.正式设计的计算如下n ㊃C R =2πσr-φ=90ʎ-βm +αRe =CV a v cos βmȡ2ˑ105ìîíïïïïïï(2)式中:n 为叶片数;C 为弦长;σ为实度比;r-为相对半径;φ为安装角;α为气流攻角;Re 为雷诺数;v 为气体动力黏度.2.2 三维成型依据理论设计,建立新风扇三维模型控制方程,即三维空间坐标转换[9]x =r sin αy =r cos αz =x 0ìîíïïï(3)式中:r 为柱坐标系下圆柱面半径;α为翼型上对应在柱坐标上一点与圆柱截面圆心的连线与zo y 平面之间的夹角;x ㊁y ㊁z 分别为柱坐标系下各点坐标.为扩大冷却风扇的稳定工作范围,提升其工作效率,采用周向前弯布局,如图7所示.前弯角推导公式为[10]θr -=θmax +arcsinr -2R +X 2b R2r-M 2+X 2bR2-arcsinR +X 2bR2M 2+X 2b R2(4)式中:r -为圆弧上任意点相对半径;θmax 为叶尖重心处前弯角最大值;M 为圆弧圆心横坐标.图7 前弯角示意图依据以上控制方程建立新风扇模型,网格划分㊁边界设置与1.2相同,如图8所示.2.3 仿真结果分析由图9(a )压强分布可知,相比原风扇,新风扇每片叶片的后端都出现了大范围均匀分布的低压区,贯穿整片叶片.由图9(b )速度分布可以看出,新风扇最大速度出现在叶尖处,这是由导风罩对气流的挤压作用与叶尖处较高的线速度共同引起的;同时新风扇表现出较为明显的梯度特征,随着风扇的转动,其梯度分布也趋于向风扇转动方向移动.图8 新风扇三维模型图9 新风扇进口处空气状态如图10(a )所示,新风扇出口端低压区分布较小;压强最大值出现在叶片中上部,高压区则均匀分布在大部分出口截面上.如图10(b )所示,速度呈现出均匀环状分布,最大速度出现在0.7倍叶高以上范围内;由于轮毂处空气无法流通,因此速度最小.3 仿真结果对比对原风扇和新风扇进行仿真性能对比,结果如图11所示.由图11(a )可以看出,随着风量的增大,两风扇仿真值吻合度逐渐提升.从图11(b )可以看出,两风扇全压效率都呈现出上升趋势,新风扇整体略高于原风扇,在体积流量达到10.24m 3㊃s -1时,新风扇较原风扇全压效率提升了0.81%.在图11(c )中,新风扇轴功率整体低于原风扇,当流量为10.24m 3㊃s -1时,新风扇比原风扇轴功率降低了0.12kW .图10 新风扇出口处空气状态图11风扇仿真性能对比4结语本文基于国内某工程机械用冷却风扇,使用三维软件建立原风扇几何模型,利用CFD在虚拟风道中进行性能模拟,并将试验数据与之进行对比;随后,使用孤立翼型法重新设计,通过推导成型控制方程建立三维模型,在相同边界下进行仿真,将两风扇的仿真结果相对比,最终得到以下结论. (1)与虚拟风道相结合的CFD数值仿真可用于模拟和预测冷却风扇性能数据.(2)孤立翼型法与翼型控制方程相结合,能更好地进行冷却风扇设计与三维建模.(3)在大部分的流量区间内,新风扇轴功率等性能参数均优于原风扇.当流量为10.24m3㊃s-1时,新风扇全压效率提升了0.81%,轴功率降低了0.12kW,证明改进有效.参考文献:[1]王企鲲,陈康民.轴流风扇两种扭叶片设计方法及其气动性能的比较[J].流体机械,2010,38(9):24-30.[2]刘全忠,宫汝志,王洪杰,等.汽轮发电机冷却风扇的数值模拟及优化[J].哈尔滨工业大学学报,2010,42(3):442-445.[3]高红霞,余建祖,谢永奇.直升机用高转速㊁小流量轴流风扇设计[J].航空动力学报,2006,21(1):119-124.[4]JIANHUI Z,CHUNXIN Y.Desi g n and Simulation of the CPU Fan and Heat Sinks[J].IEEE Transactions on Com p o-nents and Packin g Technolo g ies,2008,31(4):890-903.[5]刘佳鑫,秦四成,徐振元,等.虚拟风洞下的车辆散热器模块性能改进[J].吉林大学学报:工学版,2014,44(2):330-334.[6]刘佳鑫,秦四成,徐振元,等.基于CFD仿真的车辆散热器性能对比分析[J].华南理工大学学报:自然科学版,2012,40(5):24-29.[7]刘佳鑫,秦四成,徐振元,等.工程车辆散热器模块散热性能数值仿真[J].西南交通大学学报,2012,47(4):623-628.[8]刘佳鑫,秦四成,孔维康,等.虚拟风洞下车辆散热器模块传热性能数值仿真[J].吉林大学学报:工学版,2012,42(4):834-839.[9]廖庚华,刘庆平,陈坤,等.基于CATIA的轴流风机叶片仿生参数化建模[J].吉林大学学报:工学版,2012,42(2):403-406.[10]刘佳鑫.工程机械散热模块传热性能研究[D].长春:吉林大学,2013.[责任编辑:杜卫华]。
第36卷第2期2019年6月苏州科技大学学报(自然科学版)Journal of Suzhou University of Science and Technology (Natural Science)
Vol.36 No.2
Jun. 2019
doi : 10.12084/jassn.2096-3289.2019.02.012
扇形气膜孔冷却效率的数值模拟
王晓春,李娟(苏州科技大学环境科学与工程学院,江苏苏州215009)
摘要:为了改进气膜孔结构,提高航空发动机涡轮叶片表面气膜冷却效率,采用商业软件ANSYSCFX,对扇形圆 角扩张气膜孔的气膜冷却效率进行了数值模拟研究。根据绝热壁面假设下的温度分布云图可以发现,冷气在气膜孔 的出口区域分叉成两股,冷气分叉使最期望被冷却的孔中心的下游区域没有得到充分冷却。随后又根据流固耦合条 件下的模拟结果,提出了定量描述气膜冷却的气膜冷却效率随位置和毕渥数的拟合公式。现有文献关于气膜冷却的 研究往往采用冷却效率云图,而采用数学公式定量描述冷却效率的方法不仅更加简洁,而且更便于定量比较不同结 构气膜孔的冷却效率。关键词:气膜孔;冷却效率;吹风比;数值模拟;扇形孔中图分类号:TK124 文献标志码:A 文章编号:2096-3289(2019)02-0061-05
为了提高战斗机的性能,战斗机不断地向大推力、高速度、高机动性能的方向发展。因此,燃气透平的进 口温度需要不断上升。目前,先进透平燃气进口温度已经远远超过现有透平材料的熔点。为了保证透平高温 部件的安全运行及合理寿命,现有航空发动机端壁中广泛采用了气膜冷却技术。气膜冷却特性的控制受到 很多因素的影响,例如气膜孔的几何参数、叶片的几何参数、孔的气动参数等。在气膜孔几何形状设计方面,Goldstein等人m可能是最先关注不同气膜孔形状导致的冷却特性差异,随 后人们研究了多种扩张气膜孔结构。Wittig等人'Thole等人Pl、GritSch等人W先后比较了圆柱形、出口横向扩 张、横向向前扩张孔的冷却特性。张玲等人[3研究了气膜孔轴向倾角和展向倾角对端壁冷却效率的影响。在 众多气膜孔结构中,扇形扩张气膜孔结构的气膜覆盖面积较大,因而受到广泛的关注和研究[&"\笔者此前的 研究中也比较了三种扇形气膜孔结构的冷却特性™。为清晰的显示扇形气膜孔的气膜覆盖范围,对绝热壁面 假设下的气膜冷却效率进行了数值模拟,结果表明,冷却流在气膜孔下游逐渐分叉成两股,导致覆盖面积变 小,气膜孔下游不能被很好地冷却。在气膜孔的冷却效率研究中,广泛采用实验获得表面冷却效率云图。例如:翟颖妮等人[9]采用热色液晶 技术测量了倾斜角圆柱孔、水滴孔以及曲面簸箕孔的气膜冷却特性;常艳等人[1°]利用红外热像测温技术研究 了带有突片结构的气膜孔的气膜冷却特性;李广超等人[11]也利用红外技术测量了双出口气膜孔的冷却效率。 随着计算机性能的提升,数值模拟方法因其经济且高效的特性在众多领域得到广泛应用[12_13]。笔者通过数值 模拟的方法获得了气膜孔表面冷却效率云图,并以此为基础,尝试通过定量的方式来描述冷却效率,并提出 了精度较高的拟合公式,用公式定量描述冷却效率的方法更便于定量比较不同结构气膜孔的冷却效率。
1研究模型及方法
1.1气膜孔结构图1展示了数值模拟所采用扇形气膜孔结构。该气膜孔由纯圆柱进口和侧向扩张出口两部分组成。纯 圆柱部分的特征直径为0.6 mm,圆柱轴线与气膜孔出口平面的夹角为30°,侧向扩张角为14°。
[收稿日期]2018-03-16[基金项目]江苏省自然科学基金青年基金资助项目(BK20170382);苏州市科技发展计划项目(SNG2018045)[作者简介]王晓春(1987-),男,安徽滁州人,讲师,博士,研究方向:航空发动机叶片流动及换热的实验和数值模拟。62苏州科技大学学报(自然科学版)2019 年
(a)外形轮廓 (b)几何尺寸
图1扇形气膜孔几何结构1.2气膜冷却的数值模型
数值模拟的研究区域如图2所示。图2(a)显示了固壁绝热条件下的研究区域:主流域、冷却流域和气膜 孔。图2(b)展示了流固耦合条件下的模拟区域,比图2(a)多了固体壁面域。图2(c)给出模拟区域的几何尺 寸。研究区域的稳态流动和换热模型可以用可压缩雷诺时均的NaVier-St〇kes(RNS)方程求解。该方程可用张 量表示成如下形式dp | d{PUj ) =Q
dt d%i
d{puj ) | d{pUflii ) _dp l d r I dUj | duf 2 g dui U dj-pu/u/) ⑴dt dxj d%i d%i ^ \ dxj dxj 3 ij d%i IJ dxj
^f1+i^[Ui{pE+p )]=i^l(k+^r采用商业软件ANSYSCFX14.0对该方程进行数值计算。选择RNGk-s揣流模型。该软件计算时利用有限容积法,模拟区域划分和尺寸介绍如图2所示。计算网格利用商业软件ANSYS ICEM 14.0生成。气膜孔 采用六面体0型网格,其他区域采用普通六面体网格。对壁面附近网格加密,使边界层部分第一层网格的y+ 值全部小于1。通过对网格加密来验证无关性,加密前后,气膜孔附近的平均温度偏差小于3%,由此认为加 密后的网格已经符合要求,作为后面模拟研究中采用的网格。
(a)主流和冷却流区域;(b)流体和固体区域;(c)计算模型尺寸
图2数值计算模型
2结果及讨论
2.1扇形气膜孔出口的冷却流分叉现象为了研究单纯由气膜孔引起的冷却效率差异,采用绝热壁面假设,忽略固壁导热的影响,采用商业软件 ANSYS CFX对图2(a)所示的研究区域进行数值模拟。主流和冷却流均采用理想气体。模拟工况的边界条件第2期王晓春,等:扇形气膜孔冷却效率的数值模拟63
表面温度/K ■ 900 ’
850799749698648丨5981 547 BR=0.5孔出口区域■ 900 844789733678622I 567■511 朋=1.0
900 ^9 ■■ 900 ^9843 842 841 |785 728 670 613 ■ 556■4985/2=1.5784726668610j 552* 4941 BR=2.0782723664605546A488BR=3.0
\J
图3不同吹风比下扇形气膜孔冷却效果比较2.2定量描述气膜冷却效率为了研究固体壁面加人气膜冷却后的综合冷却效果,选择图2(b)所示模型为研究对象。主流和冷却流 工质仍采用理想气体。由于固体导热对综合冷却效果影响很大[8],所以选取7种不同导热率it(60.5、36.0、 15.1、3.98、0.58、0.15、0.058 WJn^K)-1)的工况进行数值模拟。模拟边界条件的设置如下:主流入口速度 设为12〇111_8-1,温度7^为900 1^主流入口雷诺数办为11 988;冷气入口温度7;为486 1^,速度7(!为 1.024 ml-1(相应的吹风比为2);主流和冷却流通道侧壁面设为绝热。为定量描述气膜冷却效率,在固壁表面选取了 6个特征位置,如图4所示。固壁表面气膜孔出口下游方 向设为%方向,气膜孔的圆柱轴线与固壁表面的交点设为x
轴的原点,6个特征位置分别在%轴上的坐标为1Z)、4D、8Z)、12D、16Z)、2(W)处。通过数值模拟,得到7种导热率下,固壁表面的温度分 布7;,进而可以算出6个特征位置的冷却效率。这里的冷却 效率为综合冷却效率,公式如下
设置如下:主流入口速度设为120 m_s'温度L为900 K,主流入口雷诺数为11 988(特征长度选主 流通道入口的水力直径);冷气入口温度7;为486 K,吹风比(57?)分别为0.5、1.0、1.5、2.0、3.0,相应的冷气入 口速度F。分别为0.256、0.512、0.768、1.024、1.536 m*S-1;除进出口外,研究区域所有表面均设置为绝热表面。 吹风比的定义如下br=M^_ (2)P〇〇V 00
冷却后的壁面温度分布如图3所示。图3中标尺最大最小值为各自温度分布图的最大最小值。可以发 现:(1)在气膜孔的出口区域,冷气有两个最小值区域,说明冷气在出口处已经分叉成两股。同时,随着吹风 比的增加这种分叉更加明显。(2)在气膜孔出口下游区域,冷气明显分成了两股,吹风比越大,这种分叉程度 越显著。最期望被冷却的孔中心的下游区域,反而变成局部高温。(3)吹风比从0.5增加到1.5时,气膜孔下游 区域温度明显降低。此时增加冷气量,会起到比较好的冷却效果。吹风比从1.5增加到2.0时,温度变化不明 显。吹风比从2.0增加到3.0时可以发现,气膜孔下游区域温度反而上升。因而,扇形气膜孔在吹风比为1.5 时具有最佳的冷却效果。
冈媒
图4孔出口下游6个特征点64苏州科技大学学报(自然科学版)2019 年固壁温度7;与研究位置%和固壁导热率A;有关,研究中,常采用毕渥数(所)描述导热率的变化,因而, 冷却效率可以表示成7]=rj(x,Bi)(4)
如果选择固定位置的%,那么气膜冷却效率将会是毕渥数的函数即(5)这里执定义如下
Bi=h;1k(6)
其中[表示固体壁面的厚度,1.8 mm#为固体表面的对流换热系数,由ANSYSCFX计算获得。笔者对6个特征位置的冷却效率进行拟合,公式(7)在6个特征点均可以得到比较理想的拟合结果。BUD, -Bi /D2 /r,x
V=A~Cie -C2e (7)
式中d、^、仏、/^、^是拟合方程的拟合系数。图5给出了 6个特征位置的拟合曲线。
T] = 0.878 - 0.458e ^ - 0.298e ^ ^ = 〇-632-0.316e 3239 -0.196e
1 nr7) 0.7 r7?
Bix=ADBi;;=0.470-0.234e 2 747 -0.118e 0238
0.5「乃
0.4
0.30.2
0.1c=&D
Bi
7] = 0.388 -0.196e - 0.076e
0.4 0.3 0.2
°'1h t=12D
10.4
0.30.2
= 0.340-0.175^ 2427 -0.052e 0228 77 = 0.307-0.158^"^ -0.039^"^r7? 0.35 rV
Bi^=16DBi
0.300.250.200.150.100.0512120.0
x=20D 1.5 5.0Bi
"lo.o
图5特征位置冷却效率与毕渥数的关系此外利用标准差分进化算法,对7种导热率工况下6个特征位置的数值模拟结果进行拟合,可以得到如 下的拟合公式7y=0.264+0.655e-_/6.628-0.150e-_.493)-0.612e-_/6628 • e_(肌493) (8)根据拟合公式(7)和(8),当端壁导热率趋于无穷大,即毕渥数趋于0时,气膜冷却效率具有最小值;当端 壁导热率趋于0,即毕渥数趋于无穷大时,气膜冷却效率具有最大值。变化趋势与实际相符合。同时,利用公 式(7)(如图5所示)和公式(8)所得7种导热率工况下6个特征位置的计算值,与采用ANSYS CFX模拟值 的最大偏差小于5.2%。