一种基于动模态分解的翼型流动转捩预测新方法

考虑自由转捩的定常/非定常流动Navier-Stokes方程数值求解,对于翼型流动细节的精确模拟和气动力的精确预测均具有十分重要的意义。采用动模态分解(DMD)方法进行流动稳定性分析,再结合e N方法,提出了一套适用于翼型绕流的转捩预测新方法,称为DMD/e N方法。相比于传统的线性稳定性分析方法,DMD方法不需要求解附面层方程和线性稳定性方程,也没有引入平行流假设,具有更好的理论适用性和算法鲁棒性。开展了NLF0416、S809和SD7003等翼型的转捩预测数值验证研究,通过与实验结果以及与传统的基于线性稳定性分析的e N方法的比较,验证了本文所发展的转捩预测新方法在预测翼型的定常流动和非定常流动转捩方面的正确性,也表明了该方法具有解决含层流分离泡的翼型绕流转捩预测的能力。     

Self-sustained shock oscillations on airfoils at transonic speeds

Self-sustained shock wave oscillations on airfoils at transonic flow conditions are associated with the phenomenon of buffeting. The physical mechanisms of the periodic shock motion are not yet fully understood even though experiments performed over fifty years ago have demonstrated the presence of oscillatory shock waves on the airfoil surfaces at high subsonic speeds. The unsteady pressure fluctuations generated by the low-frequency large-amplitude shock motions are highly undesirable from the structural integrity and aircraft maneuverability point of view. For modern supercritical wing design with thick profiles, the shock-induced fluctuations are particularly severe and methods to reduce the shock wave amplitudes to lower values or even to delay the oscillations to higher Mach numbers or incidence angles will result in expanding the buffet boundary of the airfoil. This review begins with a recapitulation of the classical work on shock-induced bubble separation and trailing edge separation of a turbulent boundary layer. The characteristics of the unsteady pressure fluctuations are used to classify the types of shock-boundary layer interaction. The various modes of shock wave motion for different flow conditions and airfoil configurations are described. The buffet boundaries obtained using the standard trailing edge pressure divergence technique and an alternative approach of measuring the divergence of normal fluctuating forces are compared to show the equivalence. The mechanisms of self-sustained shock oscillations are discussed for symmetrical circular-arc airfoils at zero incidence and for supercritical airfoils at high incidence angles with fully separated flows. The properties of disturbances in the wake are examined from linear stability analysis of two-dimensional compressible flows. The advances in high-speed computing make predictions of buffeting flows possible. Navier–Stokes solvers and approximate boundary layer-inviscid flow interaction methods are shown to give good correlation of frequencies and other unsteady flow characteristics with experiments. Finally, passive and active methods of shock oscillation control show promising results in delaying buffet onset to higher Mach numbers or incidence angles, thus enhancing the transonic performance of airfoils.     

Interaction of wake vortices with the ground

The aim of this study is to describe wake vortices near the ground, including rebound. A calculation code solving the incompressible Navier–Stokes equations has been developed; it allows us to compute unsteady two-dimensional flows and steady three-dimensional flows having a preferential direction. The two-dimensional version of the code allows the capture of the rebound of the vortices near the wall as well as the creation of secondary vortices. In the three-dimensional case, the results are compared with experiments. The importance of the axial velocity is also analyzed.     

振荡扰流片非定常流场的N-S方程数值模拟

 通过有限差分法求解不可压缩流体的二维非定常Navier-Stokes方程,研究扰流片作振荡运动引起的强迫非定常分离流特性。着重研究了雷诺数、缩减频率等因素的影响。得到涡在扰流片上生成、发展、脱落以及脱落后的演化过程;并给出扰流片和平板表面的压力分布。计算结果和实验结果符合良好。     

大攻角非对称涡特性及消除侧向力和力矩方法

本文提出了一个消除非对称力的头部外形方案。通过风洞实验,论证了方案是有效和稳定的。实验发现由对称涡向非对称涡的转变有一个不稳定区,建议用动态测力,得到与激光蒸汽屏流场相一致的结果,并研究其相互关系。     

上游尾迹与涡轮转子泄漏流相互作用数值模拟

叶轮机内部流动本质上是周期性非定常的,研究涡轮转子叶尖区域的非定常相互作用机理,对提高小展弦比高负荷涡轮性能具有重要意义.利用数值模拟方法研究了上游静子尾迹与涡轮转子叶尖泄漏流的非定常相互作用,分析了定常结果、时间平均结果以及瞬时时刻结果的流动图画.结果表明:上游静子尾迹与涡轮转子尖区二次流的相互作用能明显影响泄漏涡和机匣通道涡的时空演化规律,从而改变转子尖区的损失分布.上游尾迹在转子通道中传播时,诱导泄漏涡和通道涡区域出现周期性的扰动涡对,扰动涡对沿着泄漏涡和通道涡的轨迹向下游运动,使得转子尖区二次流结构呈现周期性变化.      

不同轴长比椭球体的绕流计算与分离形态分析

用拟压缩性方法求解不可压Ｎ－Ｓ方程,对４∶１旋成体椭球和其他轴长比椭球的层流绕流进行数值模拟。分别用二阶中心差分格式和三阶迎风紧致差分格式来逼近无粘项。从计算所得的表面摩擦应力线、空间流线和横截面流场特性,分析了１０°≤α≤７０°迎角范围内的分离和脱体涡系形态,算出了相应的气动特性。计算结果表明,三阶迎风紧致差分比二阶中心差分可更好地模拟分离涡结构。迎角直到４０°时,算出的流态仍然对称,而迎角为７０°时,算出的流态为非对称非定常,与实验流态基本相符。椭球体横截面形状越扁平,开式分离就越提前产生。椭球体横截面形状横向越窄,分离区就越小。由于改变椭球体横截面形状具有改变分离和脱体涡系的效果,本文的研究结果有助于改善三维物体的气动力特性,为改进机体设计提供理论依据     

圆锥脉冲等离子体诱导流场特性分析

对等离子体诱导流场特性进行研究,有利于解决双稳态非对称分离涡带来的连续比例控制困难的问题。在封闭光学玻璃箱体内,应用介质阻挡放电等离子体对20°顶角圆锥附近静止大气进行了定常和脉冲循环控制,对等离子体诱导的圆锥截面绕流速度场进行了二维PIV测量,对定常控制和脉冲循环控制下最大绕流速度及最大轴向涡量进行了比较分析。实验结果表明：相对于定常控制模式,脉冲循环控制下沿垂直于圆锥截面对称面径线分布的时间平均切向速度和轴向涡量范围较广;在脉冲循环控制下,动量传递的主要表现在离散涡的形成而不是气流的加速。     

基于RANS/LES混合方法的分离流动模拟

飞行器在大迎角、快速俯仰机动时,流场中含有大尺度、非定常的涡结构,传统雷诺平均Navier-Stokes （RANS）模型不能准确模拟流场结构,根据国际上相关研究的发展趋势,需要采用混合RANS/大涡模拟（LES）模型来对复杂分离流动进行准确模拟。本文对基于分区混合与湍流尺度混合的双重RANS/LES混合计算模型进行发展与应用。通过典型简化模型的静、动态湍流大分离流动,测试和验证所采用的脱体涡模拟（DES）类方法,重点研究改进的延迟DES （IDDES）模型在动态问题应用中的正确性和有效性,并对所采用的数值模拟方法和相应的计算软件的可靠性、鲁棒性以及精度进行了考核验证。典型算例包括超声速圆柱底部流动、跨声速方腔流动、NACA0015机翼深失速分离涡模拟等。计算表明：发展的IDDES类混合计算模型可有效解决对数层不匹配的问题;对于定态非定常分离流动,DES、DDES、IDDES等模型计算结果差别不大,随着流动的非定常特性增强,IDDES模型的优势逐渐显现;对于动态非定常分离流动,则需要采用IDDES类模型。     

超声速对流推力矢量喷管的数值模拟

利用二次流对流剪切激励主射流方式来产生射流对流推力矢量方法是一种效率很高的矢量化喷流新方法,又称为对流推力矢量法。本文运用三阶精度MUSCL TVD有限体积格式求解二维非定常、可压缩、雷诺平均N-S方程和Spalart-Allmaras湍流模型,计算了对流推力矢量喷管的超声速射流问题,获得几种不同二次流吸气压强下对狭缝和喷管肩部壁面压强分布、剪切层诱导旋涡的结构特征和主射流矢量角等结果,计算结果与实验数据对比取得良好的一致性,并对计算结果进行了详细分析和讨论。     

A novel high-order scheme for numerical simulation of wake flow over helicopter rotors in hover

Accurate prediction of tip vortices is crucial for predicting the hovering performance of a helicopter rotor. A new high-order scheme(we call it WENO-K) proposed by our research group is employed to minimize numerical dissipation and extended to numerical simulation of unsteady compressible viscous flows dominated by tip vortices over hovering rotors. WENO-K is referred to as an adaptively optimized WENO scheme with Gauss-Kriging reconstruction, and its advantage is to reduce dissipation in smoo...     

突然起动圆柱的非定常绕流计算

本文用分区方法定量模拟高Re数情况下突然起动圆柱的初期对称流动。绕圆柱的二维不可压流场可分为无粘区、附着粘性区和分离区。无粘区有解析解，用非定常层流边界层方程解附着粘性区，用离散涡方法处理随时间变化的分离区。试用了混合的Rankine和Lamb涡模型及二次离散涡。给出了Re=300、9500两个算例，基本符合有关实验结果。     

脱体涡流型翼面非定常气动力计算方法

本文研究了细长三角翼在脱体涡流型时变攻角所产生的非定常气动力的机理,并采用了一种时间历程的,考虑翼面前、后缘涡面随局部气流拖出并由时间发展逐步形成的数值方法,称之为涡面生成法,计算结果表明该方法是稳定和有效的。     

压气机线性振荡叶栅气弹试验研究（一）：非定常气动响应

 为了加强对叶轮机三维气动弹性机理的理解,为数值模拟方法提供试验验证数据,开展了压气机线性叶栅气动弹性试验研究。本研究分为两部分,本文为第1部分,研究压气机线性振荡叶栅的气弹稳定性,第2部分重点辨别叶尖间隙对此振荡叶栅气弹稳定性的影响。组建了线性叶栅试验台,其试验段由7个真实可控扩散叶型组成,中间叶片由凸轮机构驱动模拟三维一阶弯曲振动。使用外置压力传感器在6个叶高处测量叶片表面的定常压力分布和3个折合频率下叶片表面非定常压力分布。由管传递函数修正连接管造成的非定常信号失真。运用影响系数法构造协调叶栅任意叶片间相位角下的非定常气动响应结果。试验结果表明：非定常气动响应具有完全的三维特征;折合频率和叶片间相位角对叶栅的气动弹性稳定性具有决定性影响。     

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

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

低背鳍对细长平板三角翼流场影响的

 对细长平板三角翼及其对称面上加低背鳍组合体在低速风洞进行了二维粒子图像测速（PIV）实验,三角翼后掠角为82.5°,背鳍当地高度与三角翼当地半展长的比值为0.6,实验迎角为30°,无侧滑角,基于三角翼根弦长的雷诺数为2.33×10⁶。实验结果表明：单独细长平板三角翼分离涡流场对称、定常;加上背鳍后,组合体分离涡流场变得定常、非对称和非锥型。实验结果证实了低高度背鳍对细长平板三角翼分离涡的稳定性起着削弱和破坏的作用,初步验证了前人关于细长锥体分离涡的稳定性理论,并给出了30°迎角下分离涡失稳后的具体表现特性。     

双三角翼前缘剖面形状对涡运动的影响

采用数值模拟和理论分析相结合的方法,研究了前缘剖面形状对双三角翼涡运动的影响,分析了前缘剖面形状对三角翼、双三角翼涡运动影响的不同机理 :对三角翼,尖前缘可以形成组织最好的涡结构,但对于双三角翼,圆前缘生成的旋涡结构较靠近翼面,涡结构紧密,诱导能力较强,可以形成有利的涡涡干扰,使内翼涡通过剪切层向外翼涡输入涡量更加容易,合并涡变得更加稳定,推迟了涡破裂,而且由于涡较靠近翼面,因而可以产生较高的非线性涡升力,这同传统的认识是不一致的。     

具有叶顶间隙涡轮转子叶栅流动的拓扑及旋涡结构观测

为了了解具有顶部间隙的涡轮转子叶栅流道内及间隙内的流动状况,采用五孔球头测针和五孔微型束状测针分别对叶栅流道和间隙进行了测量,同时对端壁及叶片壁面进行了流动显示,采用拓扑分析方法对显示结果进行了详细分析,探讨了间隙存在时叶栅各种旋涡的形成机理。测量及显示结果表明：由于顶部间隙的存在,在叶栅顶部形成如泄漏涡等复杂的涡系结构,这些涡系之间及它们与上通道涡之间发生强烈的相互作用,明显增大了叶栅的顶部损失;在叶栅尾缘附近存在着部分回流区域。     

跨声速压气机转子叶尖流场旋转不稳定现象的数值研究

对某跨声速压气机转子在不同工作流量下的叶尖非定常流场进行了数值研究.结果显示:大流量状态下,该转子叶尖流场几乎不发生振荡.此时,叶尖流场可以按定常流场进行分析;小流量状态下,叶尖泄漏涡大幅振荡,相邻叶片通道内的叶尖泄漏流之间也存在周期性相互干涉.其结果是在稳定状态时出现由于叶尖泄漏涡的振荡及其周向传播造成的"旋转不稳定"现象."旋转不稳定"流场结构主模态旋涡个数大约为40%的叶片通道个数;其周向尺度占据2~3个栅距.     

轴流压气机转子尖部三维复杂流动Ⅱ——数值模拟研究

在实验研究和理论分析的基础上,对亚声速压气机转子尖部复杂流动做了数值模拟研究,旨在进一步深入研究叶尖间隙和进口总压分布对转子叶尖复杂流动的影响。首先通过对与实验条件相同的转子尖部复杂流动的数值模拟研究,校验了数值模拟结果,并分析了转子尖部复杂流动速度场、压力场和涡量场分布特性。然后通过改变叶尖间隙尺度和进口总压分布,研究了二者对近叶尖复杂流动的影响。研究结果表明:当叶尖间隙小于1%叶片弦长时角区旋涡的发展是导致转子失速的关键;而当叶尖间隙大于2%叶片弦长时叶尖泄漏涡的发展是导致转子失速的关键;改变进口总压分布可以合理地组织转子尖部流动并扩大转子工作裕度。此外,通过观测近叶片吸力面二维涡线的发展趋势可以判断叶尖复杂流动的发展状态。