跨声速压气机转子叶尖非定常流动的数值模拟

为了进一步加深对转子叶尖区域非定常流动现象的认识,针对某一跨声速单转子轴流压气机,采用三维非定常数值方法开展了详细研究。对单转子在不同进气条件以及不同工作流量下分别进行非定常模拟,来研究该转子叶尖区域复杂流动结构以及叶尖非定常流动的发展过程。结果表明:不同进气条件下,转子叶尖区域流场结构形式表现不同.当进口非轴向进气时,叶尖非定常流动呈单通道周期形式;而轴向进气时,叶尖区域出现了类似于“旋转不稳定”的沿周向传播的非定常流动现象,且随着工作流量的减小,叶尖区域流场结构也相应发生改变.      

旋转帽罩电加热防冰瞬态过程研究

基于欧拉两相流理论,对旋转帽罩水滴撞击特性进行了数值模拟,提出了旋转部件表面对流换热计算方法,并基于改进Messinger结冰模型开发了旋转帽罩电热防冰计算程序,对旋转帽罩瞬态防冰过程进行了数值模拟研究.结果表明:旋转运动对旋转帽罩表面对流换热起到增强作用,且在供给相同加热热流密度时转速越大,防冰表面温度越低;防冰系统启动阶段,旋转帽罩表面会发生结冰和冰脱落现象,考虑结冰过程后系统的响应时间缩短;当电加热功率相同时,周期电加热防冰方式更为节能;当电加热能耗相同时,周期电加热方式系统响应更快.     

小型垂直轴风力机叶片结冰风洞试验与数值计算

利用风洞试验和数值模拟相结合的手段,研究了小型垂直轴风力机叶片在旋转状态下的结冰特性以及结冰后翼型与风力机的气动特性变化,以期为建立较复杂的大型水平轴风力机叶片旋转试验系统和研究其结冰机理、防除冰技术提供参考。试验在东北农业大学自行设计的利用自然低温的结冰风洞中进行,获得了采用NACA0018翼型的小型2叶片垂直轴风力机风轮在5种尖速比下的结冰分布：风力机叶片结冰遍布叶片整个表面,随着尖速比的增大,结冰形状出现不对称性。同时,数值模拟结果表明：叶片结冰后,随着尖速比的增加和结冰量的增多,升力系数降低阻力系数增大的趋势明显,风力机的功率系数也随之下降。分析发现,叶片结冰导致不同旋转角下叶片翼型周围的压力场和速度场发生了不同程度的变化,从而气动特性发生变化,影响了风力机性能。     

旋转部件复杂表面水滴撞击计算

为模拟旋转体结冰问题的过冷水滴运动及撞击过程,基于欧拉方法及单旋转坐标系模型,建立了三维旋转水滴运动模型,并提出了相应的数值求解方法。采用单旋转坐标系对空气及水滴两相流动过程进行处理,通过引入惯性力,将惯性系下的周期转动边界转换为定常流动边界;利用欧拉方法,使用单向耦合形式描述空气—水滴流场;在单旋转坐标系下,向控制方程内引入科里奥利加速度及牵连加速度,进行非惯性系下欧拉方程的修正,从而描述水滴运动过程;采用有限容积求解器对空气及水滴运动的控制方程组进行求解,通过引入源项定义单旋转坐标系下的惯性力,得到空气流场及水滴场的速度、体积分数分布,进而得到表面水滴撞击特性。采用上述方法对旋转帽罩与叶片模型进行算例分析,结果表明所建立的旋转水滴计算方法有效,对比静止状态表面的结果,旋转对帽罩的水滴撞击特性影响甚小,而对桨叶存在显著影响;由于帽罩具有中心对称的特性,因而旋转带来的切向速度变化对其水滴撞击特性影响不明显;桨叶表面水滴收集系数随旋转角速度增大而增大,同时收集系数在表面的分布会向迎风方向偏移,较大的角速度对应了更为显著的收集系数增幅与偏移现象。     

含碰摩故障的转子-盘片-机匣系统动力学特性分析

以转子-盘片-机匣耦合系统为研究对象,基于ANSYS有限元软件建立有限元模型,通过实验测试得到的固有频率验证了模型的有效性;采用点点接触来模拟叶片和机匣之间的碰摩故障,并通过测试的响应数据,进一步验证了所采用模型的正确性和可行性.最后,仿真分析了定转速情况下不同侵入量、定侵入量和定间隙下不同转速对系统动力学响应的影响.研究结果表明:叶尖碰摩主要激发高倍频成分,且在系统俯仰刚体模态频率附近会出现倍频幅值放大现象;定侵入量条件下,高转速碰摩情况时系统的频率成分减少,系统的法向碰摩力也随着转速的增大而变化;定间隙条件下,在高转速碰摩情况时倍频幅值会超过转频,在一次连续碰摩过程中会出现多次碰撞反弹现象.      

反叶片角向缝机匣处理影响轴流压气机性能的机理

采用试验和非定常数值模拟方法研究了反叶片角向缝机匣处理对亚声速轴流压气机性能的影响.试验与数值计算结果均指出反叶片角向缝机匣处理后,转子在降低12%左右最大等熵效率的基础上获得30.1%左右的失速裕度改进量.详细的转子叶顶流场分析表明反叶片角向缝机匣处理改善了叶顶间隙泄漏流动,消除了实体壁机匣时叶顶通道低速泄漏流形成的堵塞,因此转子稳定性得以提高.同时转子叶片与处理缝的相对位置变化会影响处理缝的扩稳能力.处理缝中回流、处理缝形成的喷射流与叶顶通道主流的相互作用都造成较大的流动损失,进而使转子等熵效率降低.反叶片角向缝机匣处理前移约55%叶顶轴向弦长的数值研究结果表明,与反叶片角向缝机匣处理比较,反叶片角向缝机匣处理前移获得的失速裕度改进量降低约10.4%,但等熵效率的恶化程度降低近64%.      

动力效应对民机起飞构型气动特性影响的数值研究

采用在点对点多块结构化网格系统上求解三维可压缩雷诺平均N-S（Navier-Stokes）方程的方法,研究了发动机动力效应对民机起飞构型气动特性的影响.首先,通过涡轮动力模拟器风洞试验模型及民机高升力构型标准模型,对研究方法进行了验证,计算和试验结果的良好吻合说明研究方法用于模拟发动机动力效应及预测民机高升力构型气动特性是可行的.其次,针对分别安装通气短舱和动力短舱的某翼吊涡扇发动机民机起飞构型,研究了发动机动力效应对其气动特性的影响规律.结果表明:在通气短舱基础上预测的起飞最大升力系数、失速迎角以及阻力和力矩,在考虑动力效应后会产生明显的不同.动力效应可使外形最大升力系数增大并延迟机翼失速,但同时也带来了低头力矩增大和升阻比降低的不利影响.建议在对机翼及其高升力装置设计结果的数值校核中使用考虑进排气的动力短舱代替目前常用的通气短舱,这能为设计改进提供更可靠的依据.      

低速轴流压气机旋转失速演化机制研究

为深入探讨轴流压气机旋转失速机理,利用全环非定常数值模拟在出口边界添加节流阀模型模拟压气机节流过程,通过试验测量与数值模拟结果对比分析,澄清了低速轴流压气机失速先兆及失速团的产生和演化过程。研究结果表明：压气机在转子叶片安装角存在几何偏差的情况下,spike型失速先兆是由2阶模态的整阶扰动演化而来,并最终直接形成2个大尺寸全叶高的失速团;与轴对称转子模型相比,考虑存在几何偏差的转子模型的数值模拟结果更为接近试验测量结果。     

旋转盘腔盘罩间隙比的敏感性分析

为保证寿命限制件之一的涡轮盘满足适航性要求,采用单向流固耦合(fluid structure interaction,简称FSI)数值方法,研究转静系旋转盘腔盘罩间隙比的变化对转盘安全性的影响机理.并且,由流阻、换热效果和应力分布三方面构成的工程评价体系对盘罩间隙比的影响进行评价以及敏感性分析.结果表明:盘罩间隙比的变化能够影响旋转盘腔内流动结构的强度,从而改变盘面换热效果和转盘温度分布,导致与温度梯度相关的热应力也发生变化.随着盘罩间隙比的增加,旋转盘腔的流阻损失基本不变,转盘迎风面平均换热效果呈现先增强后减弱的趋势,转盘整体应力水平上升,以及出现在盘心处的最大等效应力值增大.盘罩间隙比的变化能够从材料许用应力和实际使用载荷两方面影响涡轮盘的失效概率,因此,在涡轮盘腔设计阶段,需考虑盘罩间隙比对涡轮盘安全性的影响.      

Analysis of the buoyancy-induced heat transfer in a rotating cavity

In order to get a better knowledge of the heat transfer in compressor cavities of aero-engines, the simplified rotating cavity with two-plane discs, a shaft and a cylindrical rim has been investigated numerically and compared with the available measurements. The numerical results in agreement with the available experiments show large-scale instabilities. The disk local Nusselt numbers show mainly radial rising distributions for the heated disks with radial rising temperature profiles. In the present work, at the Reynolds number of 20000, the disk local Nusselt numbers are the correlations of the local Grashof number to the power of 1.89~2.6, and the value of the power is increased as the rotational Reynolds number goes up. At the rotational Reynolds number of 800000, the local Nusselt numbers are the correlations of the local Grashof number to the power of 0.68~2.6, and the value of the power is decreased as Reynolds number goes up. The area-averaged disk Nusselt number is the correlation of the Reynolds number to the power of 0.479 and the rotational Grashof number to the power of 0.12.