二次空气系统盘腔一维非定常计算方法

通过对盘腔结构的一维算法展开研究,提出适用于具有非标准结构、包含多个进出口的复杂盘腔的一维建模规则,并采用特征线法,实现盘腔结构的一维非定常流动与换热的计算,通过将广义的变截面、有摩擦、考虑对外热交换的盘腔一维可压缩非定常流动与换热的偏微分方程组转换为常微分方程组,进而获得盘腔内部沿流动方向压力、流量、温度等系统参数的特征线数值解,采用商业计算流体软件CFX对其进行验证.结果表明:所提出的盘腔建模与计算规则能够直观且准确地描述盘腔一维非定常流动及换热特征,仿真结果与CFX计算结果最大偏差约为9.4%,可应用于具有非标准结构的多接口复杂盘腔的模拟,可以作为二次空气系统盘腔结构一维非定常计算研究的有效手段.      

二次空气系统一维非定常计算方法

研究了燃气轮机二次空气系统的非定常计算方法.通过将整个二次空气系统划分为两部分,特征"长"元件以及局部损失元件.特征线法用于单个特征"长"元件的求解,通过将连续方程、动量方程、能量方程构成的偏微分方程组转换为常微分方程组并联立求解,可以充分考虑内部流动与换热的相互影响;基于压力修正的流体网络法则用于由多个局部损失元件组成的局部流体网络,首先对动量方程组进行求解,进而对各自的能量方程进行求解.如何在采用这两种方法进行求解的上述两类元件之间进行数据交互将被讨论.此外,所提出的时间步长定义方法以及改进的元件计算模型使得二次空气系统的非定常分析得以实现.计算程序经与某系统定常实验数据及单个元件的非定常计算数据进行对比,在10s的计算时间内,基于特征线法得到结果与CFD数据偏差达到10.2%,结果表明该计算程序可以满足二次空气系统非定常计算的要求.     

气壁镀镍和冷却剂入口对再生冷却的影响

考虑冷却剂入口模型及气壁镀镍,对液体火箭发动机推力室再生冷却通道和冷却剂进行三维流动与传热耦合计算.采用经验公式计算燃气侧对流及辐射换热,冷却剂为甲烷,考虑其物性随温度和压力的变化.所得结果表明:冷却剂入口二次流及突扩流场的叠加,使局部压力损失变大,影响进口下游较大区域流动传热状况;气壁镀镍能提高气壁的耐高温性,降低气壁锆铜的温度.      

涡轮叶片凹槽状叶顶非定常流动传热特性的数值研究

为了研究涡轮叶片凹槽状叶顶的非定常流动传热特性,以GE-E3第一级静叶和动叶为研究对象,采用ANSYS-CFX数值求解三维Reynolds-Averaged Navier-Stokes(RANS)和标准k-ω紊流模型。数值预测的叶顶换热系数分布与实验数据吻合良好,从而验证了数值方法的可靠性。数值计算结果表明:静叶尾迹对动叶顶部的流动和换热特性影响显著。压力面侧前缘区域和吸力面中间位置的流场受动静干涉影响显著。叶顶表面的换热系数脉动主要出现在靠近前缘的凹槽底部表面和再附着线附近及吸力面侧肩壁。靠近动叶尾缘区域的换热系数脉动同时受动静干涉作用和下游流场的影响。定常计算得到的换热系数在前缘冲击区和分离线附近高于非定常时均值,在压力面侧肩壁附近小于非定常时均值。定常计算得到的平均换热系数要高出非定常计算结果3.5%。     

入口带倒圆结构的旋转长孔计算模型

通过对旋转长孔的各种复杂工况开展CFD仿真,对现有经验关联式进行修正,并基于流体网络法对修正后的计算模型进行开发,开发完成的模型采用计算流体力学商业软件Simerics进行验证,结果表明:该模型考虑了入口倒圆结构对流动的影响,对于旋转长孔的模拟准确度较高,在低雷诺数及超声速等各种流动状态下均具有较广的适用性,并且只需基本的几何特征参数即可进行一维数值模拟,因此可以作为二次空气系统一维数值仿真的标准化模型使用.      

航空发动机整机空气系统流动与传热数值模拟

为了研究航空发动机整机空气系统内流动与传热特性,本文以涡扇发动机整机空气系统为对象,主要包括压气机盘腔、涡轮盘腔、低压轴前后腔、前后轴承腔等旋转盘腔和大量的流阻换热单元以及相邻的结构部件。所研究的是一个多腔相连,多进口,多出口,流固耦合传热的复杂问题,分析了发动机真实状况下空气系统内的流动与传热特性。采用Mixture多相流模型进行数值计算得到了该系统的速度场、压力场和温度场。结果表明：该系统内的流场是复杂的多涡流场,多个出口出现燃气倒灌;整个系统的腔压大致从前到后逐渐升高;由于主流通道燃气的入侵,导致高、低压涡轮盘的温度较高,后轴承腔内滑油的冷却作用较为明显。本文的工作使得对发动机整机空气系统的研究不再局限于一维计算,为空气系统的设计提供了理论依据。     

固体火箭发动机反向喷管推力终止过程的压力冲击问题

本文介绍一个反向喷管推力终止的冷模拟实验.实验表明,突然打开固体火箭发动机反向喷管后,燃烧室中可能出现强烈的压力冲击,这种压力冲击来自喷管喉部节流和形成准稳态流动之前的非定常过程.本文对这一效应做了估算.用特征线法及激波拟合法进行了一维非定常流动的数值分析.预估的压力-时间曲线与实验结果相当一致.     

旋转状态下矩形微小通道流动与换热试验研究

为了研究微小通道结构在航空发动机涡轮叶片中应用的前景和可行性,以空气为冷却介质,在Re＝1000～3000、转速为0～500 r/min、Ro＝0～3.5×10-3条件下,对水力直径为1 mm的旋转微小通道组的流动和换热特性进行试验研究。结果表明：微小通道流阻系数呈现粗糙壁通道特征,通道临界Re≈2350,流阻系数以及临界Re随转速增加未见明显改变。在静止状态下,通道组综合换热系数随Re增大而增大,换热系数分布沿流动方向逐渐减小;在旋转状态下,通道组平均综合换热系数略有增大,旋转对换热特性的影响随着流动的发展而增大     

1.5级大子午扩张涡轮非定常流动和传热特性分析

大子午扩张涡轮端壁二次流与热负荷之间的关系,对于端壁冷却非常重要。本文采用几何约化法,对某1.5级大子午扩张涡轮进行数值模拟,研究了大子午扩张涡轮上端壁非定常流动和传热特性。计算结果表明：大子午扩张涡轮上通道涡尺度较大且位置发生改变,沿径向向下移动约20%叶高;R1出口泄漏涡、通道涡和尾迹是造成S2流动和传热非定常性的主要因素;传热与二次流密切相关,对传热研究必须与流动相结合。研究结果将有助于提高对大子午扩张涡轮端壁非定常流动和传热特性的认识。     

2级涡轮内部流动定常与非定常计算差异研究

为获取进而认识涡轮内部流动状态,以某2级约化形式的动力涡轮为研究对象,分别对其进行定常和非定常数值计算和分析.研究表明:定常计算与非定常计算对涡轮内部流动的模拟结果,如叶片表面的压力分布、叶排进出口的气流角、叶片通道中的二次流流向涡、展向涡、叶片通道中的损失等,均存在差异;流动的非定常性越强,定常与非定常计算结果的差异越大,且该差异大小对于静叶与动叶呈相反的展向分布规律.     

旋转带肋U型通道内流动与换热的数值模拟

为了获得能够模拟高参数涡轮叶片内冷通道换热效果的模型，数值模拟了旋转状态下U型通道内的流场和温度场，比较了数值模拟与实验的结果。结果表明：所采用的计算程序和模型与实验结果吻合。旋转状态下，通道内各面换热的变化是和通道内流场的变化密切相关的；哥氏力在垂直于旋转半径截面上的不均匀分布引起流动较大变化，对通道内各面换热的影响比较大。哥氏力的作用较大幅度强化指向面换热，小幅强化两侧面换热，而弱化背向面换热。对于带肋通道，总体上阻力系数随着旋转速度的增加而升高。     

轴向通流旋转盘腔内换热的数值模拟

以数值模拟的方法,采用旋转坐标系稳态方程,研究了轴向通流旋转盘腔内的换热。主要讨论了流动对换热的作用以及旋转系下各力对换热的影响,给出了盘腔内的换热随各无量纲参数的变化规律。研究结果表明:哥氏力的增大削弱了轴向通流旋转盘腔内的换热,惯性力和浮升力的增大增强了换热;反映在无量纲参数上,随着进口雷诺数的提高,盘腔内的换热增强;随着瑞利数的提高(提高转速),盘腔内的换热经历一个缓慢变化-突增-缓慢变化的过程,换热的突然增强是冷气流穿透盘腔所致。      

瞬态法测量高转速旋转盘表面传热系数

将实际发动机涡轮盘冷却系统简化为中心进气旋转盘,用瞬态实验的方法对该结构的换热特性进行了研究.以高转速旋转换热实验台为基础,建立了转静系盘腔中瞬态换热的实验流程和数据处理方法,得到了转盘表面的努塞尔数,并研究了流量系数、旋转雷诺数和出气间隙比对努塞尔数的影响.研究结果表明:对于中心进气冷却结构,在转盘低半径处,转盘表面的表面传热主要由冷气冲击控制,随着半径的的增大冲击对换热的影响减弱.转盘表面的努塞尔数随冷气流量系数的增大而增大,随旋转雷诺数的增大而增大,随出气间隙比的增大而减小.      

旋转U型通道流动与换热特性的数值研究

应用k-ω SST（shear stress transport）湍流模型,计算分析旋转U型通道在不同进口雷诺数（10000~60000）和高旋转数（0~2.013）范围内的流动与换热特性.结果表明:在静止和旋转状态下,进口雷诺数越大,努塞尔数越大.相比于同一工况下的静止状态,旋转显著增强了径向外流直通道的换热强度,径向内流直通道换热强度增大不明显.旋转数对U型通道换热的影响主要通过改变哥氏力和浮升力的大小.受哥氏力的影响,径向外流直通道后缘面换热增强,前缘面换热减弱.浮升力诱发了近壁面的流动分离,使得径向外流直通道前缘面不同位置处的换热强度随旋转数的增加而先减小后增大,计算得到的临界旋转数变化规律与实验测量结果保持一致,即无量纲距离参数与临界旋转数的乘积为定值.      

高压涡轮后腔流阻特性与瞬态换热试验研究

在航空发动机二次流系统中,涡轮盘腔的流动和换热问题伴随着复杂的几何、流动及热边界条件,为探究其流场和换热特性对发动机设计的重要作用,对一具有预旋进气孔和高、低半径出口的高压涡轮后腔内的流阻特性和转盘盘面的换热特性进行了试验研究,主要应用瞬态液晶测试技术对转盘表面的对流换热特性进行了测量。在试验中,旋转雷诺数Re_ф变化范围为8×10~5~1.0×10~6,无量纲二次流量(流量系数)C_W变化范围为5.29×10~3~1.19×10~4。试验结果表明:腔内压力及流阻特性受进气流量C_W和转盘转速Re_ф的影响;转盘表面的换热随着半径的递增以及预旋比β_p的增大而增强;出口湍流参数λ_T对换热特性影响很小。     

外缘预旋进气的旋转空腔主盘局部换热及流阻特性研究

某型发动机的实际涡轮盘腔冷却结构被简化成外缘预旋进气的旋转空腔模型,以实验方法研究了旋转雷诺数（Ｒｅω）、冷气雷诺数（Ｒｅｚ）、哥拉晓夫数（Ｇｒ＊）对主盘局部换热及流阻特性的影响。实验发现,冷气雷诺数、旋转雷诺数的增大,使主盘局部换热系数在盘缘附近迅速增大,而在接近盘心时出现负值;阻力系数随旋转雷诺数的增大而减小,随冷气雷诺数的增大而增大。     

预旋进气小尺寸涡轮叶片冷却的流动换热

为了获得一种具有直通式冷气预旋进气系统的小型燃气轮机涡轮叶片的流动与换热特性,采用气热耦合计算方法进行数值研究,分析了总压损失、冷却效果和涡轮效率随预旋角、冷气雷诺数和无量纲质量流量的变化规律。结果表明,涡轮叶片预旋进气冷却的总压损失随冷气雷诺数和无量纲质量流量的增大而增大,但基本不受预旋角大小的影响;涡轮叶片的冷却效果随预旋角的减小、冷气雷诺数或无量纲质量流量的增大而增强,但不会改变其表面的温度分布特征;预旋进气冷却时的涡轮效率随冷气雷诺数的增大、预旋角或无量纲质量流量的减小而提高。