上游尾迹与涡轮叶栅通道涡相互作用研究

采用三维粘性非定常数值模拟方法研究了上游尾迹与涡轮叶栅通道涡的相互作用,对定常、非定常时均以及瞬时时刻流动机理进行了分析.结果表明:上游尾迹的非定常作用一方面增强了叶栅通道涡的径向涡,使得流动损失增大;另一方面能够一定程度上抑制通道涡中流向涡的发展,对控制损失起到正面作用,端区的综合非定常效应取决于两者之间的平衡.在本文计算条件下,上述两方面综合影响使得通道涡的非定常损失增大.      

基于分离涡模拟的起落架气动噪声研究

用基于SA湍流模式的非定常RANS、分离涡模拟(DES)和延迟分离涡模拟(DDES)分别对四轮基本起落架模型进行了数值模拟,并根据所得的非定常流场计算了表面声压级分布和声压谱。三种方法所用的时间成本大致相同,URANS略低于其它两种,而DDES由于在附着流动区更好地保持RANS特性,故时间成本略低于DES。计算显示,非定常RANS在附着流动区能够得到合理的结果,但不能准确刻画分离流动的流动形态。DES和DDES都能较好地刻画起落架绕流的定常和非定常特性,DDES的结果略好于DES。因为起落架流动属于大分离流动,所以DES也能够得到相对正确的结果。研究结果验证了分离涡模拟在起落架大分离非定常流动预测与噪声预测中的可行性,对减小起落架噪声的方法研究具有一定的意义。     

细长旋成体背风面非对称流场的DES模拟

利用脱体涡模拟DES方法对细长旋成体非对称绕流进行数值模拟,比较了基于S-A湍流模型的DES与RANS方法对背风面非对称分离涡的数值模拟能力,发现在大迎角非对称涡未破裂情况下,DES与RANS均能模拟出与实验相符的非对称绕流;当迎角增大,背风面非对称涡发生破裂时,RANS无法准确捕捉到背风面流场的非定常性,而DES能准确预测非对称涡的飘起与破裂,并与实验值接近。计算结果表明:与RANS方法相比,DES方法具有更好的模拟大攻角,大分离流动的能力,尤其在非对称流场的分离涡破裂模拟方面具有明显优势,能够更真实模拟出细长旋成体背风面分离涡破裂之后的非定常流动特征。     

悬停共轴刚性旋翼直升机涡尾迹高精度数值模拟

针对共轴刚性旋翼直升机流场涡尾迹的高精度数值模拟需求，结合结构化运动嵌套网格技术、高精度数值格式和非定常RANS方法，建立了共轴刚性旋翼直升机的非定常涡流场高精度数值模拟方法及高效并行化自主软件平台。运用所发展的方法和软件平台，进行了类“X-2”构型共轴刚性旋翼直升机悬停状态涡尾迹结构的高精度数值模拟。结果表明，所发展的数值模拟方法具有较高的计算效率和计算精度，能够捕捉到精细的涡尾迹结构及其非定常生成、发展和演化过程。     

液氢绕回转体非定常空化的流动特性

为研究航空发动机内部低温燃料的非定常空化流动特性,采用大涡模拟(LES)对液氢绕回转体的非定常流动特性进行了分析。通过与实验结果的对比可知,所采用的数值计算方法能够有效地模拟液氢绕回转体流动的非定常过程。讨论了空化演化过程和旋涡动力特性,研究结果表明:空化的非定常演变过程大致可以分为3个阶段:附着空穴生长阶段、大尺度空泡发展阶段和小尺度空泡发展阶段,且反向射流是造成空穴不稳定及脱落的主要原因;探讨了空化与旋涡之间的交互作用,明确了空穴前端和尾端处的涡量源于旋涡伸长项,空穴内部涡量源于旋涡扩张项,空穴交界面、反向射流头部区域涡量源于斜压扭矩项。      

瞬时畸变极值预测的随机涡方法

对于湍流型进气畸变压力场,文献[1]提出了一种随机涡模拟方法。将非定常N-S方程与统计学理论相结合,将随机涡尺寸与畸变指数峰值相关联。对此,作者在文献[2]中作了详细的推导、分析和预测,发现在原模型的应用中,存在着两个明显困难:合理确定瞬时畸变分布;迅速和准确确定涡核半径a。实验和模拟结果表明畸变分布介于正态分布与Weibull分布之间变化。      

微型飞行器低雷诺数矩形扑翼非定常气动特性的数值模拟

采用双时间步方法求解三维可压缩非定常N-S方程,数值模拟了微型飞行器低雷诺数矩形扑翼的非定常绕流,首先将得到的结果与文献进行了对比,数据间具有较好的一致性.然后针对不同的展弦比、减缩频率及初始攻角,计算了矩形扑翼的非定常气动特性及表面流态和动态压力分布,并分析了翼尖涡对扑翼非定常气动特性的影响.     

基于非定常动量源方法的全机耦合气动分析

基于非定常动量源方法开展了直升机全机耦合气动分析研究。首先,以定常动量源模型为基础给出了高效的求解方法并验证其准确性。然后,将动量源项随时间变化添加至虚拟桨叶的所在位置,建立了能够体现流场非定常特性的非定常动量源模型。在此基础上,应用所建立的模拟方法深入开展了直升机全机耦合气动计算与分析。最后,结合飞行平衡模型分析了飞行参数对直升机气动特性的影响。研究结果表明,所建立的非定常动量源方法可以准确地模拟出桨叶的涡尾迹,可用于直升机全机耦合的气动分析研究。     

短舱对螺旋桨滑流影响的IDDES数值模拟

基于非结构重叠网格技术,对短舱与螺旋桨滑流间的相互作用进行了非定常数值模拟研究。为了更好地捕捉螺旋桨尾涡的细节信息,计算采用基于Spalart-Allmaras模型的改进延迟脱体涡模拟(IDDES)方法,并在非定常计算过程中运用网格自适应技术以提高流场特征的空间分辨率。研究结果表明:IDDES方法获得的拉力系数计算值与实验值吻合良好,短舱的存在会增大螺旋桨的拉力系数;短舱对螺旋桨桨毂涡的结构影响较大,但对桨尖涡的螺旋结构影响较小;对单独螺旋桨算例来说,桨尖涡与桨毂涡的失稳发展过程都具有周期性,且在有/无短舱情况下桨尖涡的失稳位置相同,失稳后桨尖涡之间配对融合过程一致,从而说明桨毂涡对桨尖涡的失稳没有影响。     

低雷诺数翼型层流分离现象大涡模拟方法

为准确预测翼型低雷诺数条件下出现的层流分离流动现象,发展了基于结构化拼接网格的大涡模拟方法。控制方程为Favre滤波Naver-Stokes方程,并选用多种亚格子模型。空间离散采用AUSM格式以及高阶WENO格式,时间推进采用显式方法和隐式方法。以SD7003翼型在雷诺数60 000及4°迎角下的层流分离流动为研究对象,对比分析了数值格式、亚格子模型、网格尺度对流场预测结果的影响。数值结果表明:划分的计算网格能够有效解析小尺度流动结构,基于隐式亚格子模型、采用AUSM格式和双时间步推进的大涡模拟方法能够准确预测流动分离、转捩、再附等复杂流动现象,计算得到的平均压强系数与雷诺应力分布与文献数据吻合较好。与转捩模式计算结果对比进一步表明:发展的大涡模拟方法能够准确预测翼型低雷诺数层流分离流动非定常演化过程,为下一步的研究工作建立了有效的数值模拟手段。     

模拟厚度翼型非定常运动的涡方法

本文研究二维厚度翼型的大攻角非定常运动。用沿翼型表面分布的源和面涡来模拟翼型的厚度效应和升力效应,应用离散涡方法模拟前缘分离涡层和尾涡层。计算了大攻角翼型作俯仰运动的情形,与实验结果符合较好。     

高负荷涡轮端区非定常流动相互作用研究

采用三维黏性非定常数值模拟方法研究了某型高压涡轮端区非定常流动相互作用,着重研究了上游静子尾迹与转子二次流的非定常作用机制,同时还分析了负荷分布、激波等对端区非定常流动的影响。结果表明,静子尾迹的非定常作用一定程度减小了转子轮毂二次流的径向涡量;尾迹对流向涡量的影响取决于尾迹沿叶高的分布,当吸力面一侧的尾迹具有与二次流方向相反的流向涡量时,二次流的流向涡量减小;非定常效应还使得转子叶片根部负荷略为减小,也一定程度上抑制了转子轮毂二次流的发展。此外,受静子尾缘激波的影响,转子叶片表面负荷分布发生明显的周期性变化,导致叶片表面承受较强的非定常力,在涡轮设计中必须考虑。另外,通过计算涡轮级中的熵增和熵产,定量地分析了端区非定常相互作用产生的损失,并得到了一些有意义的结论。     

三角翼前缘脱体涡和尾涡的卷起及其相互作用

本文应用二维非定常比拟和恪子涡(Vortex-in-Cell)方法,数值模拟了三角翼前缘涡层的卷起以及与尾涡的相互作用。由于使用了数千个点涡和较小的空间网格,获得了前缘涡层小涡配对出现的较大的离散涡,模拟了紧卷涡层的不稳定现象。这是以前类似的工作未曾实现的。     

基于DDES模拟的叶顶泄漏流与尾迹非定常干涉机理

叶顶泄漏流是涡轮叶片内部气动损失和非定常性的重要来源之一，其包含多个空间和时间尺度流场结构。采用延迟脱体涡高精度湍流模拟方法（DDES）研究了某高压涡轮动叶叶顶泄漏流，分析了泄漏涡结构及其非定常特性，进而研究了叶顶泄漏流与尾迹间的非定常干涉现象及其机制，最后基于非定常熵输运方程分解并研究了泄漏流区域的损失机制。研究发现：DDES方法可以精细捕捉泄漏流和尾迹区的多尺度流场结构及其相互干涉作用；泄漏流与尾迹间的相互干涉会导致尾迹涡轨迹发生明显偏移，并伴随着涡破碎现象的产生，具有很强的非定常效应；基于非定常熵输运方程的损失分析表明，涡破碎现象造成的损失在脉动损失中占比可达23.3%，是非常重要的损失来源之一。     

涡对与物体的相互作用研究

本文从非定常形式的Ｎ－Ｓ方程出发，采用有限差分法数值模拟了入射涡对与圆柱及翼型相互作用的流场结构，旨在分析讨论其复杂的非定常流动现象，如二次涡及三次涡的形成，多涡配对，非定常分离及剪切层的演化等，以及不同的物形及涡对与物体间不同的相对位置对影响。此外，文中还通过分析物面压力分布随时间的变化来探讨非常分离的流动特性。     

某涡轮级三维流场非定常数值研究

基于耦合Lantry-Menter转捩模型的SST(Shear stress transport)湍流模型,对某涡轮级在3个雷诺数工况下的全三维流场分别进行了定常、非定常数值模拟,目的在于对比分析两种计算方法所得结果之间的差异。计算结果表明,进口雷诺数为4×104时和5×105时,定常与非定常计算的涡轮级效率的差异较大,最大相差0.6个百分点,而当雷诺数为8×104时,两者差异很小;非定常计算能够较好的模拟叶尖泄漏涡的发展过程、低雷诺数下上游尾迹与下游转子叶片吸力面分离边界层干涉诱导卷起涡的形成、输运及其所导致的叶片吸力面表面压力大幅波动等非定常流动现象。最后,给出了转子出口截面上周向与径向的非定常压力脉动的分布,为基于Lighthill声类比方法的计算气动声学研究奠定了基础。     

超声速小攻角下的圆锥转捩攻角效应

采用高阶差分格式求解三维可压缩滤波Navier-Stokes方程,对超声速零攻角和小攻角圆锥边界层转捩进行了空间大涡模拟研究.在三维Tollmien-Schlichting(T-S)波扰动作用下,计算得到了流场非定常转捩过程.计算结果清晰地展现了非定常大尺度结构演化特征和攻角效应引起的非对称转捩现象,流场发卡涡、横流涡等转捩主控大尺度结构与实验一致.     

扑翼升力特性的非定常涡格法计算研究

针对西北工业大学航空学院微型飞行器研究室研制的PY-1型扑翼机,建立了扑翼的运动模型,并采用非定常涡格法模拟了该扑翼的升力特性,分析了迎角、飞行速度对升力特性的影响。计算结果与实际飞行试验吻合,验证了非定常涡格法模拟扑翼非定常气动特性的正确性和有效性。研究结果可为扑翼机的气动设计提供参考。     

轴向倾斜缝对压气机转子叶尖泄漏涡的非定常控制数值研究

为研究轴向倾斜缝对转子叶尖泄漏涡的非定常控制作用,采用瞬时叶片排数值模拟方法对带不同数量轴向倾斜缝的Rotor 67转子进行非定常数值模拟。结果表明,轴向倾斜缝能有效扩宽转子的稳定工作范围;轴向倾斜缝缝数应使其通过频率与预测的转子叶尖泄漏涡非定常脉动频率相近为宜;当加入轴向倾斜缝这一结构激励源后,叶尖泄漏涡不再表现出自身的脉动频率,而只以轴向倾斜缝的激励频率为主频进行脉动;在小流量工作区间,越接近失速点,轴向倾斜缝对于叶尖泄漏涡的控制作用越强,叶尖泄漏涡的压力脉动也越剧烈。     

基于RANS-LES混合方法的翼型大迎角非定常分离流动研究

使用雷诺平均Navier-Stokes方程-大涡模拟（RANS-LES）混合方法中的延迟分离涡模拟（DDES）方法,模拟了NACA 0015翼型在大迎角下的静态绕流和强迫振荡运动并与实验值进行了比较。在大迎角静态翼型大分离流动模拟中,DDES方法捕获了非定常RANS计算未能获得的机翼背风面的涡脱落现象。在所采用的RANS和DDES模型中,基于剪切应力输运（SST）湍流模型的SST-DDES混合方法给出的时均压力系数分布与实验吻合得最好。在大迎角强迫振荡翼型绕流模拟中,DDES方法得到的非定常气动载荷与实验值吻合得很好,正确地反映了最大迎角处阻力和俯仰力矩的阶跃性突变,而非定常RANS计算则给出了完全错误的趋势。