环形叶栅分离旋涡频谱特性研究

施加非定常气流脉动激励将大幅度减少叶栅分离损失。而激励的效果与本身的气流脉动脉动频率以及强度等因素有关。为了考察非定常气流脉动激励对叶栅气流分离的作用，以及激励作用效果最佳时激励脉动频率与叶栅分离旋涡脱落频率之间的关系，进行了叶栅分离旋涡脱落频率的测量。实验在环形叶栅实验台上进行，利用IFA300热线风速仪等测试仪器进行实验，分析了不同实验条件下环形扩压叶栅后流动脉动的分布情况。结果表明，在不同迎角以及不同测量位置处(叶中和端壁压)，所测环形扩压叶栅后的流动分离均存在某一特征频率。此外对导流叶栅后和环形扩压叶栅后的流动情况的对比说明了在叶轮机械非定常流中流动分离频率的锁定(lock on)现象^[1]。     

亚声速扩压平面叶栅尾迹动态压力场测量与分析

为研究叶片在不同攻角下引起的气流分离对叶栅出口气流紊流度的影响,借助动态压力测量设备和测试技术,完成了某扩压平面叶栅在进口马赫数为0.677,攻角分别为0°、-10°和+8°三种典型工况下,尾迹非定常流动的测量。通过测量尾迹区域沿栅距方向和轴向的尾迹动态压力,并对动态压力数据进行时域分析,得到叶栅尾迹非定常流动的时均总压和压力脉动云图,揭示出尾迹区流动过程,同时还与叶栅的气动性能和气流稳定性进行了关联。     

进口气流角对三维振荡涡轮叶片非定常流动影响的数值模拟研究

进一步发展了非定常雷诺平均N—S方程求解程序,并基于影响系数法,在三种不同振型下就进口气流角对三维涡轮振荡叶片绕流的影响进行了数值模拟研究。数值模拟结果表明,所发展的程序对振荡叶栅流动模拟具有较好的精度;进口气流角对振荡叶栅内部非定常流动以及叶片表面的非定常气动力有着较为重要的影响,且在不同模态下,进口气流角对振荡叶片非定常流动的影响规律表现不同。     

压气机叶栅角区分离中速度PDF双峰现象的DES数值模拟

为了研究在实验中发现的角区分离区域内速度概率密度分布(PDF)双峰现象并探索其机理,利用脱体涡模拟(DES)对压气机叶栅中角区分离流动进行数值模拟,成功模拟出双峰现象,并讨论了双峰现象的可能机理。DES模拟结果中,双峰点位置以及内外峰随双峰点位置变化而此消彼长的发展趋势与实验结果一致。对来流攻角脉动特性的分析表明,角区分离中双峰现象与叶栅前缘附近非定常来流和角区内分离流动的相互作用有关。     

基于DDES方法的叶栅分离旋涡的非定常流动数值研究

为了研究叶栅分离旋涡的非定常流动特性,选取亚声速叶栅为研究对象,进行了基于延迟分离涡模拟方法 (DDES)的非定常数值模拟。针对进口Ma0.3的情况,对比分析了进口攻角、叶栅稠度、叶片弯角、叶片中弧线弯度分布和叶片厚度等因素对于叶栅分离旋涡运动的影响。结果表明,当攻角增大时,尾迹附近流场会逐渐由定常模式转化为非定常对涡脱落模式,同时对涡脱落频率不断减小;随着攻角的进一步增大,叶背分离强度不断增加,叶背分离涡开始单独脱落,并主导流场旋涡运动。强烈的叶背分离涡脱落引起了叶片气动力的剧烈脉动,其幅值是小攻角状态下由尾迹对涡脱落引起叶片气动力脉动的30倍以上。叶栅几何参数的改变对于流场旋涡运动同样有很大影响,流场同样呈现出"定常尾迹"、"对涡脱落"、"叶背分离涡脱落"这三种主要的旋涡运动形式之一。叶背分离点的相对位置则是影响旋涡非定常运动形式和旋涡脱落频率的主要内在因素。     

三维振荡叶栅内部非定常流动数值模拟研究

进一步发展了非定常雷诺平均N-S方程求解程序,并基于影响系数法对三种不同振型下的三维涡轮振荡叶片绕流问题进行了数值模拟研究.数值模拟结果表明,所发展的程序对振荡叶栅流动模拟具有较好的精度;振型对振荡叶栅内部非定常流动以及叶片表面的非定常气动力有着较为重要的影响,其中叶片运动导致叶栅各通道中流量的再分配以及出口气流角的变化等对叶片的负荷分布会产生较为显著的影响.     

振动叶栅中尾迹/激波相干的数值研究

研究了尾迹/激波相干对振动叶栅非定常流动的影响.数值方法利用了先进的谐波方法.该方法基于傅立叶变换,把叶栅非定常流变量分解为一时间平均量和若干谐波分量,继而,时域非定常N-S方程被转化为频域定常傅立叶系数方程.定常傅立叶系数方程的求解采用了伪时间推进方法.为加速数值解收敛,采用了本地时间步长和多重网格技术.通过NASA67转子算例结果表明:上游尾迹对振动叶栅非定常流动具有显著影响.当存在尾迹相干时,叶片通道内的激波在轴向方向上表现出较大的漂移运动,叶片所受的非定常力幅值变大,同时失去了时间周期性.      

非定常环境中叶栅总压恢复系数的数值研究

利用以Ni格式和隐式双时间步法求解二维非定常N-S方程为基础的数值模拟方法研究了流动的非定常性对涡轮叶栅总压恢复系数的影响,并与文中提出的定常模型和准定常模型进行了对比。结果发现,非定常流动能够提高叶栅的总压恢复系数,但这是建立在与准定常流相比较的层面上,而和定常流比较时,非定常流只会使总压恢复系数下降。      

一种快速模拟振荡叶栅非定常流的数值方法

为了提高叶轮机内周期性非定常流数值计算效率,发展了一种计算振荡叶栅非线性非定常流的谐波平衡方法,将周期性非定常流动参数采用傅里叶级数表示,从而将传统的定常流动数值计算方法应用于非定常流问题求解中,大大提高了非定常流计算效率。计算结果表明,发展的谐波平衡方法与线性方法、非线性时域法计算结果吻合较好,且其计算效率比时域法高1个数量级。另外,系统研究了谐波阶次对计算效率与计算精度的影响,认为在气流未出现严重分离时,应采用低阶谐波平衡方法,以减少计算耗时。此外,从压气机振荡叶栅非定常流场计算结果可以看出,叶片吸力面诱导涡的生成频率与振荡频率并不一致,涡生成频率要高于振荡频率。     

振荡叶栅气动性能实验研究(英文)

为了提供双凸翼型叶栅在所有叶间相位角及实际折合频率下高亚音和跨音定常和振荡气动性能的基本数据,在 NASA路易斯研究中心的跨音振荡叶栅试验设备上进行了一系列试验。为了进行该试验,研制出并使用了非定常气动影响系数法,即叶栅中一次只有一个翼型振荡,测出该振荡翼型与邻近静止翼型上非定常压力的矢量和,藉此就能确定在特定叶间相位角下相当于全部叶片振荡的叶栅的非定常气动性能。      

不同振动形式下的翼型失速特性

通过求解雷诺平均Navier-Stokes方程,研究了翼型在强迫振动和自然结构振动下的大迎角流场特性,尤其是失速迎角附近的流场和气动特性.研究结果表明:在接近颤振临界速度情况下结构自然振动可以引起翼型大尺度的分离,导致失速分离涡提前出现;强迫性的沉浮运动和俯仰运动在一定幅度下也可以引起失速性质的大分离,而且沉浮和俯仰振动的频率和振幅都是影响翼型大尺度分离的重要因素.     

考虑压气机动静干涉效应的非设计工况性能数值模拟

1引言非定常流动是叶轮机械内部固有的本质属性,而现有的设计体系中不论是在设计或分析时基本上采用定常流动的假设,从而忽略了非定常流动对叶轮机械在气动、结构强度等方面的影响。目前,国内外对非定常流动的研究比较多,但其中采用非定常流动数值模拟方法预估压气机非设计工况     

翼型大迎角振动对气动性能影响的实验研究与初步分析

通过在NF-3低速风洞专门研制的翼型模型及相应的俯仰和沉浮振动机构,选用NACA0012翼型进行大迎角下不同频率的振动实验,研究了模型振动平均状态下对其气动力特性的影响情况,并在N-S方程基础上对振动流场进行了初步分析。实验与计算研究的结果表明：在临近定常失速迎角的大迎角条件下,翼型的振动可以引起旋涡分离,导致翼型升力减小和失速迎角的提前。就所涉及的两种振动模式而言,俯仰振动的影响大于沉浮振动,所以,模型设计和加工时要特别注意加强机翼弦向的扭转刚度。     

用N－S方程计算翼型非定常粘性大攻角绕流

本文给出了用ＮＳ方程对翼型非定常粘性大攻角绕流的数值模拟。控制方程为二级时均可压缩完全ＮＳ方程；湍流模型采用双层代数涡粘性模型￣［１］。使用近似因式分解ＡＤＩ差分格式离散求解，网格是用保角变换方法生成的相对翼型固定的Ｃ型网格。本文给出两类典型非定常绕流数值模拟结果：翼型过失速常攻角周期流动和大攻角强迫俯仰谐振非定常绕流。并与国外的实验和计算结果进行了比较，表明了本方法准确、高效的特点。     

大展弦比复合材料机翼失速颤振分析

研究了大展弦比复合材料机翼在较大迎角状态下的失速颤振特性,探讨了结构几何非线性和由复合材料剪裁产生的刚度耦合效果对机翼失速颤振特性的影响.首先,将复合材料机翼建模为转角和位移均可为有限值的非线性薄壁单闭室截面Euler梁,并在综合考虑结构几何非线性、气动非线性和材料各向异性对机翼运动状态的影响的基础上,建立机翼的运动微分方程.然后,使用小扰动分析的方法得到机翼在平衡位置附近的振动方程,采用ONERA半经验的非定常失速气动力模型,获得机翼在平衡位置附近的非线性失速颤振分析方程.最后,利用谐波平衡法求解并判定机翼颤振稳定性.通过算例,首先验证了算法的正确性,然后研究了几何非线性对失速颤振的影响,并讨论不同的复合材料铺层方式导致机翼失速特性的改变.     

Numerical Investigation of the Effects of Wind Tunnel Model Vibrations on the Measured Aerodynamic Properties of Airfoils and Wings

This paper presents a numerical prediction of the unsteady flow field around oscillating airfoils at high angles of attack by solving unsteady Reynolds-averaged Navier-Stokes equations with SST turbulence model in order to simulate the effects of wind tunnel model vibrations on the aerodynamic properties of airfoils,especially high-aspect-ratio wings in a wind tunnel.The effects of the phase lagging between different modes of oscillations,i.e.,the airfoil plunging oscillation mode,the pitching oscillation mode,and the forward-backward oscillation mode,are also studied.It is shown that the vibrations (oscillations) of airfoils can cause the unsteady shedding of large-size separated vortex to precede the stationary stall incidence,hence lead to a stall onset at some earlier (lower) incidence than that in the steady sense.The different phase lagging has different effect on the flow field.When the pitching oscillation mode has small phase lagging behind the plunging oscillation mode,the effect of vibrations is large.Besides,if the amplitude of the oscillations is large enough,and the different modes of vibrations match or combine appropriately,the unsteady stall may occur 2° earlier in angle of attack than the case where airfoils keep stationary.     

等离子体激励控制激波与边界层干扰流动分离数值研究

针对高超声速进气道激波与边界层干扰流动分离控制问题,提出了一种低功率重频非定常激励方式,并基于雷诺平均Navier-Stokes(N-S)方程,从唯象学的角度出发,将等离子激励简化为功率密度源项,对比研究了定常与低功率重频非定常等离子体气动激励的作用机理与控制效果。结果表明:定常激励的能量沉积作用对于激波控制非常有效,并可诱导出斜激波,但是对于流动分离控制而言,其能量沉积显然过于强大,反而会使流动分离更加严重,无法满足控制要求;当采用低功率重频非定常激励方式时,对于不同功率密度的情况均存在最佳激励时长与频率,当功率密度为5.0×109W/m3时,最大射流速度可以达到895m/s,并且可以在一定程度上减弱激波与边界层干扰流动分离。      

Stall flutter prediction based on multi-layer GRU neural network

The modeling of dynamic stall aerodynamics is essential to stall flutter, due to the flow separation in a large-amplitude pitching oscillation process. A newly neural network based Reduced Order Model (ROM) framework for predicting the aerodynamic forces of an airfoil undergoing large-amplitude pitching oscillation at various velocities is presented in this work. First, the dynamic stall aerodynamics is calculated by solving RANS equations and the transitional SST-γ model. Afterwards, the stall flutter bifurcation behavior is calculated by the above CFD solver coupled with structural dynamic equation. The critical flutter speed and limit-cycle oscillation amplitudes are consistent with those obtained by experiments. A newly multi-layer Gated Recurrent Unit (GRU) neural network based ROM is constructed to accelerate the calculation of aerodynamic forces. The training and validation process are carried out upon the unsteady aerodynamic data obtained by the proposed CFD method. The well-trained ROM is then coupled with the structure equation at a specific velocity, the Limit-Cycle Oscillation (LCO) of stall flutter under this flow condition is predicted precisely and more quickly. In order to predict both the critical flutter velocity and LCO amplitudes after bifurcation at different velocities, a new ROM with GRU neural network considering the variation of flow velocities is developed. The stall flutter results predicted by ROM agree well with the CFD ones at different velocities. Finally, a brief sensitivity analysis of two structural parameters of ROM is carried out. It infers the potential of the presented modeling method to depict the nonlinearity of dynamic stall and stall flutter phenomenon.     

双级风扇三维粘性非定常流动的数值解

针对一台双级跨音速轴流风扇动静叶片排相互干扰的三维非定常流动,进行了N-S方程的数值求解。首先计算风扇定常流场,以定常流收敛结果为起始场计算非定常流动。计算结果表明,非定常流动对各叶片排进口气流角有较大的影响,其中静叶的变化大于动叶,与定常流结果相比,非定常的瞬态攻角增大的量大于减少的量,最多时增大近10°左右。静子叶片受到的非定常气动力变化幅度也大于转子叶片,叶片排所处的轴向位置不同,非定常流动的影响会有较大的差别。