小流量下转静系盘腔传热特性

采用热敏液晶(TLC)瞬态传热实验技术对小流量下的中心进气转静系盘腔传热特性进行研究,介绍了同步触发拍摄系统,考察了流量与转速变化对动盘面传热的影响.结果表明:在所研究工况范围内,流量变化对动盘面中低半径处传热影响较大,对高半径处影响较小;转速升高使动盘面传热效果整体增强;当湍流流动参数较小时,动盘面中高半径区域努塞尔数随半径增大而增大;采用湍流流动参数对动盘面平均努塞尔数进行归纳,实验结果与雷诺相似分析结果的偏差小于15%,一定程度上验证了雷诺相似对转静盘腔传热问题的适用性.      

对转涡轮盘腔内各参数变化对转盘壁面摩擦力矩的影响

对转涡轮盘腔内壁面摩擦力矩的大小直接影响盘腔内气体流动结构和换热效率,进而影响对转涡轮发动机的性能。应用RNG k-ε湍流模型对对转涡轮盘腔内气体流动进行了数值模拟,研究了对转涡轮盘腔间距、中心进气流量、转盘转速的变化对盘腔内壁面摩擦力矩的影响。数值模拟结果表明,盘腔间距在计算范围内的变化对转盘壁面摩擦力矩的影响不明显;转盘壁面摩擦力矩随中心进气流量和转盘转速的增大而增大。     

旋转盘腔紊流流动的数值研究

为了比较近壁紊流模型对旋转盘腔内流动特征模拟精度的影响，采用在核心区选用标准k-ε模型、在近壁区选用单方程和Launder-Sharma两方程两种低雷诺数模型的区域紊流模化方法，对具有径向出流的转-静盘腔和具有径向出流的共旋盘腔结构的流动进行了数值计算研究，计算结果表明：在计算的两个算例中，前者采用近壁单方程模型要稍优于两方程模型，后者采用近壁两方程模型要明显优于单方程模型，因此近壁区紊流特性的处理对旋转盘腔内的流场模拟精度具有重要的影响。     

低位预旋进气转静系盘腔流动特性研究

采用实验和数值模拟的方法对低位预旋进气转静系盘腔的流动特性进行研究,考察了流量与动盘转速变化对腔室内流动的影响,结果表明腔内径向静压差随流量增大而增大,随转速升高而增大;引入腔内有效角动量的概念,用角动量系数表征其与盘腔入口角动量之比,通过分析来流角动量与盘面转矩的相对大小,推导得到影响角动量系数的两个主要无量纲参数为修正的湍流流动参数与预旋来流雷诺数;对于给定结构形式的预旋系统,当改变转速时,角动量系数随修正湍流流动参数增大而减小,当修正的湍流流动参数大于0.0646时,角动量系数变化缓慢,腔内有效角动量约为入口角动量的70%～90%;当改变来流速度时,角动量系数随修正湍流流动参数增大呈现先减小后增大的变化趋势。     

中心进气旋转盘流动与换热的数值研究

本文采用混合长度模型及共轭计算的方法对中心进气转静系旋转盘、腔内的流动与换热进行了数值计算, 并与实验结果进行了比较, 证明混合长度模型对转静系旋转盘、腔的流动与换热的计算是可行的。计算结果显示: 冲击区努赛尔特数最大值不在驻点, 而是偏离驻点一段距离; 随着转速和冷气流量的增加, 盘面换热总体加强。      

有预旋进气转静盘腔中的流动和换热特性数值研究

运用RNG湍流模型对有静盘外缘预旋进气和轴向中心进气、转盘外缘轴向出气和外围屏径向出气的转-静盘腔中流动和传热过程进行了数值研究。在盘腔结构中,静止内隔片和旋转内隔片将盘腔分为内转-静盘腔室和外预旋腔室。由于盘腔流动结构具有周向周期性,取盘腔的1/30作为计算域。研究发现:旋转雷诺数Re_θ、预旋喷嘴进气无量纲流动速度C_w,_p、静盘中心轴向进气的无量纲流动速度C_w,_d和预旋比β_p都不同程度地影响盘腔内流动和换热。预旋和隔片的协同作用可以很好地起到阻止热燃气入侵转-静盘腔和提供低静温冷却气体流入涡轮叶片冷却通道。      

轮缘封严气体对涡轮转子性能影响的非定常数值研究(英文)

为了研究涡轮转静盘腔的轮缘封严气体对涡轮转子性能的影响,对其进行了详细的三维非定常数值模拟。通过与无盘腔模型比较,非定常结果显示交替进行的主流燃气入侵和封严气体出流改变了转子叶片进口的气流角度,拉伸了端壁附近的马蹄涡,因此引起涡轮出口流动条件的较大变化。存在涡轮转静盘腔但没有冷气时,主流导叶尾迹的部分流体入侵盘腔,对涡轮转子性能有正面的影响。封严气体与上游导叶尾迹的相互作用引起转子通道内熵增,当冷气量占主流流量的1.37%时,涡轮效率降低2.1%。证实了轮缘封严气体对涡轮转子性能有很大的影响,在涡轮设计中必须给予考虑。     

不同进口条件下向心流转静系盘腔流动特性研究

在典型中小型航空发动机空气系统中,用于高压涡轮冷却、封严的空气一般需要经过离心压气机叶轮背腔,因此掌握以离心压气机叶轮背腔为代表的向心入流转静系盘腔内流动特点及压力分布是保证空气系统各项功能实现的关键。采用数值模拟方法对带有向心入流的转静系盘腔流动开展研究,研究不同来流条件下不同间距比的转静系盘腔流动特点及盘腔内压力分布。结果表明:在间距比G=0.01~0.2内,不同进口条件下盘腔内的流动均为Batchelor流型,即转盘与静盘具有独立边界层,边界层之间为核心区;当径向罗斯比数远小于1时,在核心区内流动满足径向平衡方程,此时盘腔内旋转比分布决定了盘腔内压力分布;对于满足径向平衡方程的此类盘腔,盘腔内流动由进口旋转比β_0、紊流参数λ_T、间距比G决定;进一步的,得到了不同β_0,λ_T,G下盘腔出口旋转比及核心区内旋转比变化规律,分析发现小间距比工况下核心区内旋转比满足5/7幂指数关系;大间距比工况下旋转比满足修正5/7幂指数关系,通过得到的旋转比关联式可以计算出盘腔内的压力分布。     

中心进气转静盘腔的流量分配特性试验

采用实验方法, 对具有中心进气、两个出口的盘腔系统的流动特性进行了研究.获得了静盘与转盘间隙以及转盘外缘小孔的冷气流量.结果表明:盘面出流孔的流量随总流量或转速的增加而增加, 该流量与总流量的流量比随总流量的增加而降低, 随着转速增加而增加.盘面出流孔数量为原出流孔数量的一半时, 流出流量也近似是原流出流量的一半.      

直接供气预旋转静系流动特性的实验

用实验方法研究直接供气转静系预旋腔内的流动特性,得到了静盘表面的静压分布、预旋孔的排气系数以及盘缘出口流量的分配规律.结果表明:流量改变对静盘压力分布形状的影响较大,而转速改变的影响则较小;在几何条件不变以及预旋进气和中心进气流量相等的前提下,预旋孔排气系数主要由预旋进口流量系数决定;对于该实验模型和研究工况来说,盘缘出口流量系数比在0.7以上.      

转静盘腔流动瞬态响应特性研究

以航空发动机转静盘腔内非稳态流动为研究对象,基于用户自定义函数编程设置进口压力按斜坡函数变化时的非稳态边界条件,利用数值模拟方法,在滞后效应、流场演化和物性参数方面,研究了不同无量纲跃升幅值下转静盘腔流动瞬态响应机理。结果表明:无量纲跃升幅值由1.05增至1.2,盘腔内部的当地响应时间延长了81.1%,相对总温的超调量增加了157.0%,其峰值时间增加了31.3%,且进口压力波动的传播方向总是由进口指向出口;瞬态过程中,转静盘腔的容积效应和旋转效应共同作用使出口流量明显滞后进口流量,且其滞留质量流量与相对总温的超调现象密切相关;冲击射流、旋转贴壁射流和静子壁面侧的涡系构成了盘腔流场演化的主体,无量纲跃升幅值越大,流场演化过程越剧烈,响应时间越长;声速和热扩散率与盘腔响应速度密切相关。     

湍流边界层厚度对三维空腔流动的影响

采用脱体涡模拟(DES)方法开展了不同湍流边界层厚度(TTBL)下的三维空腔非定常流动数值计算。空腔长、宽、深比例为5:1:1,来流马赫数为0.85,雷诺数为13.47×10⁶ m^-1,各工况湍流边界层厚度比值为1:2:4:8。研究结果表明,湍流边界层厚度对自由剪切层的发展、空腔底部静态压力分布、脉动压力及空腔流动类型均有重要影响,且随着边界层厚度的增大,下游剪切层覆盖的范围会增大,但是剪切层增长率降低;空腔前后静态压力压差减小、压力梯度下降;腔内局部测点的脉动压力声压级下降,各阶声压峰值频率向低频方向偏移;空腔流动类型往开式流动方向转换。     

高负荷吸附式压气机叶栅数值与实验研究

针对无通道激波单纯由强逆压梯度诱导的附面层分离进行了吸附式数值与实验研究.研究对象为某大转折角高负荷吸附式压气机叶栅,利用准三维叶栅通道计算程序（MISES）进行抽吸流场数值模拟,在确定抽吸位置后进行了风洞实验验证.结果表明:抽吸后总压损失系数大幅度降低,对于单纯由强逆压梯度诱导的附面层分离,最佳抽吸位置应该位于附面层分离之后尚未充分发展之处;在确定抽吸位置时可以根据设计状态的分离状况进行;实际中需要的抽吸流量小于计算值;在数值计算中,具体的抽吸模型还有待进一步改进和修正,以使数值模拟更加准确.      

航空发动机加力燃油总管的数值计算分析

采用商业软件对航空发动机加力燃油总管进行了数值计算,研究了喷油孔的附面层与质量流量的关系,同时分别对喷油孔允许的最大和最小值进行数值模拟,就其流量值与设计值进行了对比。结果表明:喷油孔的加工质量对加力燃油总管流量的影响很大,喷油孔必须按照公差等级Ⅰ的规定进行加工,以保证其质量流量满足设计要求;对于小孔径结构的数值计算,只有充分考虑附面层的因素,才能提高数值计算的准确性。     

预旋进气转-静盘腔中的流动和传热数值研究

采用计算流体动力学(CFD)软件CFX5.7对预旋进气旋转盘腔中的流动和传热过程进行了全三维数值模拟.对4种工况进行了计算,计算结果与试验结果吻合较好,验证了所采用方法的计算精度,并分析了旋转雷诺数、无量纲切向速度、预旋比对流场结构、传热的影响.      

Numerical investigation of windage heating within shrouded rotor-stator cavity system with central inflow

The rotating disk surface temperature rise due to windage heating effect by numerically modeling the turbulent flow within a rotor-stator cavity which is available with a peripheral shroud and imposed through airflow was dealt with. The windage heating may be defined as viscous friction heating caused by relative velocity differences across the boundary layers between the fluid and the rotating disk surface. The kinetic energy dissipation process could transform the rotating shaft power into thermal heating. Commercial finite volume based solver, ANSYS/CFX was employed to numerically simulate this physical process by using the shear stress transport (SST) turbulence model. CFD results include the rotating disk surface temperature axial distribution and tangential velocity distribution of the fluid domain. The velocity difference between the result obtained by particle image velocimetry (PIV) experiments and CFD simulation are within 5%. The adiabatic disk temperature rise can be calculated by the tangential velocity of disk and fluid in large gap ratio and turbulent parameter. CFD temperature distribution results and those estimated via velocity differences are within 10%.     

Safety Analysis of Flow Parameters in a Rotor-stator Cavity

In order to ensure the safety of engine life limited parts (ELLP) according to airworthiness regulations, a numerical approach integrating one-way fluid structure interaction (FSI) and probabilistic risk assessment (PRA) is developed, by which the variation of flow parameters in a rotor-stator cavity on the safety of gas turbine disks is investigated. The results indicate that the flow parameters affect the probability of fracture of a gas turbine disk since they can change the distribution of stress and temperature of the disk. The failure probability of the disk rises with increasing rotation Reynolds number and Chebyshev number, but descends with increasing inlet Reynolds number. In addition, a sampling based sensitivity analysis with finite difference method is conducted to determine the sensitivities of the safety with respect to the flow parameters. The sensitivity estimates show that the rotation Reynolds number is the dominant variable in safety analysis of a rotor-stator cavity among the flow parameters.     

离心压缩机密封及空腔流动的一维/三维耦合建模研究

为了在保证数值预测精度的同时降低离心压缩机的建模复杂度,以一维密封泄漏计算替代完整的密封及空腔建模,开展了离心压缩机盖盘密封及空腔流动的源项建模方法研究。结果表明:基于源项方法的密封及空腔简化建模仅需在主流道建模基础上施加泄漏质量流量、总温和摩擦损失源项边界条件,即可较高精度地还原完整的盖盘密封及空腔建模的性能和关键流动细节。与完整建模相比,在全工况下源项建模预测的等熵效率的绝对偏差小于0.65%,总压比和等熵效率的相对性能偏差小于1%。该方法可在满足下游部件匹配设计精度需求的同时,减少50%以上的网格数量,并显著提高计算的经济性和收敛性。泄漏流量的计算精确度是影响源项建模法性能预测准确性的最关键因素,耦合建模可以采用一维计算中推荐的密封泄漏系数值。