对转开式转子非定常气动干扰特性分析

采用动态面搭接技术,求解非定常雷诺平均Navier-Stokes方程实现对转开式转子的非定常气动数值模拟.对10片桨叶单个转子和10×10对转开式转子构型分别进行模拟,对比分析前后两个转子间的气动干扰和滑流流动干扰.结果表明:与单个转子相比,对转开式转子前转子的拉力系数和功率系数减小,后转子的拉力系数和功率系数则都增大,并在一个旋转周期内都呈现20次周期性波动.前后转子拉力系数的频谱分析显示振荡发生在桨叶通过频率的偶数倍,且在2倍时后转子拉力的振幅最大.前转子桨尖涡与后转子的周期性干扰,引起后转子桨叶力分布的变化和非定常性.与单个转子相比,前转子后滑流轴向速度偏大,周向速度偏小.后转子对滑流有二次加速及旋涡恢复作用.      

跨声叶栅角区激波/附面层干扰的开槽控制措施

提出了使用叶根槽作为一种被动控制手段来控制跨声叶栅的角区分离问题。在压力面与吸力面的压差作用下,叶根槽可产生自发射流,为叶栅吸力面侧角区注入高能流体,从而控制跨声叶栅的角区分离问题。通过数值模拟的方法分析了在不同攻角下叶根槽对压气机叶栅性能的影响及作用机理。结果表明:在小攻角下,叶根槽射流可破坏角区环形涡,从而有效减小跨声叶栅角区分离,提高叶栅的流通能力,改善叶栅性能;在大攻角下,叶根槽射流已不能破坏角区环形涡,但仍能为角区低能流体充能,减弱角区分离,从而拓宽叶栅工作范围。在0°攻角下总压损失系数可降低11.6%,同时叶栅攻角裕度由2°拓宽为3°。      

马赫数1.95超声速欠膨胀喷流声辐射数值分析

针对马赫数为1.95的欠膨胀超声速喷流声辐射特性,采用高精度计算格式的大涡模拟（LES）方法进行数值研究。通过对喷流平均流及湍流脉动统计结果的对比,确定数值方法的精确性。细致分析了超声速喷流流动特征,特别关注了其中激波胞格和湍流相互作用现象。基于流动与声场辐射间的关联性,分析了超声速喷流马赫波、宽频激波噪声辐射特性及形成机理。数值结果表明超声速喷流的剪切层演化及激波胞格与湍流间的相互作用构成欠膨胀超声速喷流的主要噪声源。     

转静干涉对转子叶片颤振特性的影响

在不同气动工况、不同几何模型下,采用自行开发程序对全环多排的高压压气机进行了流固耦合数值模拟,分析了进口导流叶片(IGV,以下简称导叶)对转子叶片颤振稳定性的影响。通过对气动弹性标准算例4的数值模拟验证了程序在颤振领域的有效性。针对导叶-转子模型和单转子模型,考察了近堵塞点、近设计点和近失速点3个工况下,节径变化对叶片颤振稳定性的影响,给出了气动弹性最不稳定状态对应的叶片振动形式。通过对比发现,导叶作用随工况而异,近堵塞点导叶使得转子1阶弯曲模态气动阻尼提高130.63%。研究表明:导叶引起的非定常压力波反射增强了转子叶片的非定常压力扰动幅值,使得弯曲振型的颤振稳定性增强。基于单转子模型的颤振分析给出了不准确的气动阻尼值。      

超声速气流中液体横向射流的气液相互作用过程数值研究

采用Euler-Lagrange方法对来流马赫数为1.94的超声速气流中液体横向射流的气液相互作用过程进行数值研究。计算给出的射流穿透深度、液滴Sauter平均直径(SMD)及液滴速度分布均与实验吻合较好。仿真结果较详细地揭示液体射流喷雾与气流之间的强烈相互作用过程。受液雾影响,射流前形成较强激波,气流依次经过激波及液雾区域,气流速度存在两次下降过程。计算结果揭示,超声速来流可以与射流的液滴轨迹相交,气流经液雾前沿进入液雾区域后,流向往壁面偏折。本文首次发现并提出,由于气液相互作用诱导形成两组反向反转漩涡对,这对于理解两相混合过程具有重要意义。气液相对滑移速度的分析表明,液滴在穿透自由来流并开始转向时受气流作用最为显著,完成转向后气液相互作用逐渐减弱。     

进口探针支杆尾迹诱发压气机转子叶片共振的数值研究

为了探究包含支杆直径10mm进口探针的压气机罕见地出现第一级转子叶片动应变超限报警的产生机理,以更换直径8mm+3mm的探针后报警解除的情况作为对比方案,开展了进口探针支杆诱发压气机转子叶片振动的流动机理研究。以试验压气机为研究对象,利用ANSYS APDL命令流语言使用载荷步文件法对该压气机转子叶片进行了单向瞬态流固耦合数值计算。结果表明,两类探针支杆尾迹影响下,压气机转子叶片都发生了共振;引起叶片共振的激振力频率来源于支杆尾迹诱导频率与支杆通过频率的叠加,其中支杆尾迹诱导频率是由支杆尾迹涡脱落与转子叶片通过非同相位所诱导产生;支杆直径减小,支杆尾迹引起的叶片表面激振力减小,使叶片振动应变相应降低。     

Bump进气道中鼓包诱导的激波/边界层干扰特性

为了探索Bump进气道中鼓包诱导的锥形激波和机身发展而来的湍流边界层干扰问题,分析其气动优势,首先选取了半锥和半棱锥这两种与鼓包的流场结构具有一定相似性的构型作为参照,采用数值仿真方法,分别对这三类典型的三维激波/湍流边界层干扰问题进行了流场分析。在此基础之上,设计了三个不同马赫数的鼓包,并研究了设计马赫数对鼓包流场特性的影响。结果表明:当三类构型的无黏激波强度相等时,半锥诱导产生的旋涡强度最强,鼓包次之,半棱锥最弱。尽管鼓包诱导的流场非常复杂,其干扰流场却呈现出准锥形相似的特性。虽然半锥对边界层的排移能力最强,但是综合考虑边界层排移能力及进气道出口流场畸变下,鼓包最具优势,这也是其被选为超声速进气道前缘压缩面的重要原因之一。此外,在设计状态下,适当增加设计马赫数能改善鼓包排移边界层的能力,但设计马赫数太高,边界层排移能力基本不变,反而使得进气道总压损失急剧增加。      

压气机动静叶干扰的非定常流动特性

采用大涡模拟(large eddy simulation,LES)方法计算动叶尾迹对静叶干扰的流场信息。利用涡量分布揭示动叶尾迹在静叶通道内的演化过程,利用压力梯度识别激波结构及波振源,运用动力学模态分解(dynamic mode decomposition,DMD)方法对静叶通道内流场的时空结构进行模态分解。结果表明:流场中存在3处波振源,分别位于动叶尾缘、静叶前缘和静叶尾缘处;发现静叶通道内流场的频谱具有多峰值现象,模态分解的第1阶流动代表动叶尾迹在通道内随时间迁移,对应频率为动叶通过频率(blade passing frequency,BPF)是通道内旋涡非定常波动的主导频率;第2阶流动是动叶通过频率的2倍频流动,旋涡的空间尺度为1阶模态的1/2,为更小尺度的扰动。      

Unsteady Flow Variability Driven by Rotor-stator Interaction at Rotor Exit

Numerical investigation of the unsteady flow variability driven by rotorstator interaction in a transonic axial compressor is performed. Two models with close and far axial gap between rotor and stator rows are studied in the simulation. Particular attention is attached to the analysis of mechanisms involved in driving rotor wake oscillation, rotor wake skewing and flow angle fluctuation at rotor exit. The results show that smaller axial gap is favorable to enhance the interaction in the region between two adjacent rows, and the fluctuation of the static pressure difference between two sides of rotor wake is improved by potential field from down stator, which is the driving force for rotor wake oscillation. The interaction between rotor and stator is weakened by increasing axial distance, rotor wake shifts to suction side of rotor blade with 5%-10% of rotor pitch, the absolute value of flow angle at rotor exit is less than that in the case of close interspace for every time step, and the fluctuation amplitude is also decreased.     

多块协调变形的网格变形技术及其应用

基于多块网格协调变形的理论,提出了适用于叶轮机内流多物边界的更具一般性的网格变形技术,先获得块边界处的变形量,然后插值得到块内网格各点的变形量。通过二维振荡叶栅和带叶尖间隙的三维跨音风扇叶片的非定常绕流算例,证明网格变形技术在叶轮机气弹稳定性数值模拟中的应用是有效的,能够在保证网格高精度的条件下,提高网格变形速度。