基于IDDES方法的翼型结冰失速分离流动数值模拟

应用基于k-ωSST湍流模型的IDDES(Improved Delayed Detached Eddy Simulation)方法,就失速点附近翼型前缘典型双角状积冰导致的复杂分离流动进行了数值模拟研究.通过与风洞试验结果进行对比,表明对于此类分离流动问题,IDDES方法能够在壁面附近取得良好的速度预测结果,有效解析分离区域内的中小尺度湍流结构,较为准确地描述大尺度时均分离泡的再附位置和形态特征,适用于翼型结冰后复杂流动的精细分析.同时计算结果显示当此带冰翼型位于失速点附近时,角状冰后方脱落剪切层内部的旋涡不稳定析出和输运过程促进了外部流动与回流区域流动间的掺混,将导致流动发生非定常再附现象.     

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

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

三维旋转效应对桨叶气动特性影响的计算

在旋转的作用下,流动分离向下游延迟,结果大大改变了风力机桨叶的气动特性,这种由旋转所带来的影响被称为三维旋转效应。本文以桨叶表面边界层理论分析的结果为基础,给出了三维旋转效应对气流分离影响的解析关系,将其与非定常气动和动态失速模型相结合,并适当考虑旋转条件下桨叶翼型载荷的分布特性,使其能计算三维旋转效应下风力机桨叶的非定常气动特性.通过算例表明,当这一三维旋转效应模型与非定常动态失速模型耦合后再计算桨叶非定常气动特性时,能大大改善计算结果。     

翼型动失速问题研究

本文对俯仰振荡和快速上仰翼型的非定常动失速问题进行了计算研究。对绕流流场采用了高效率的分区精确数值模拟的研究方法,保证了计算结果的精度,大大节省了对非定常粘性流动至关重要的计算机时。文中单独划出前缘驻点区,对形成前缘动失速旋涡的前缘附近区域的流场进行精确计算。算例显示了动失速流场的形成、变化过程及其对气动特性的影响,计算结果和实验的比较是满意的。     

基于后缘小翼的旋翼翼型动态失速控制分析

针对后缘小翼(TEF)的典型运动参数对旋翼翼型动态失速特性的控制进行了研究。发展了一套适用于带有后缘小翼控制的旋翼翼型非定常流动特性模拟的高效、高精度CFD方法。通过求解Poisson方程生成围绕旋翼翼型的黏性贴体和正交网格,为保证后缘小翼附近的网格生成质量,建立了基于翼型点重构的方法来描述后缘小翼的偏转运动;为克服变形网格方法可能导致网格畸变的不足,发展了一套适用于带有后缘小翼控制的旋翼翼型运动嵌套网格方法。基于非定常雷诺平均Navier-Stokes(URANS)方程、双时间法、Spalart-Allmaras(S-A)湍流模型和Roe-Monotone Upwind-centered Scheme for Conservation Laws(Roe-MUSCL)插值格式,发展了旋翼翼型非定常气动特性分析的高精度数值方法,并采用Lower-Upper Symmetric Gauss-Seidel(LU-SGS)隐式时间推进方法及并行技术提高计算效率。以有试验结果验证的HH-02翼型和SC1095翼型为算例,精确捕捉了动态失速状态下的气动力迟滞效应,验证了本文方法的有效性。着重针对SC1095旋翼翼型的动态失速状态开展后缘小翼的控制分析,提出了可以体现翼型升力、阻力及力矩综合特性的关系式Po和Pc,揭示了后缘小翼振荡频率、相位差和偏转幅值对动态失速特性影响的规律。研究结果表明:当后缘小翼偏转的相对运动频率为1.0,且小翼运动规律与翼型振荡规律之间的相位差为0°时,后缘小翼能够更好地抑制翼型动态失速现象;在此状态下,当偏转幅值为10°时,SC1095翼型最大阻力系数和最大力矩系数可以分别降低19%和27%。     

一种用于分离流动的网格自适应算法

对于可压缩粘性流动,提出利用流场速度的紊乱度作为指示变量进行网格自适应.Jameson中心格式有限体积法、五步Runge-Kutta时间推进法/双时间推进法求解定常/非定常N-S方程.基于雷诺平均N-S方程模拟紊流,选用SA一方程模型.在数值求解二维静态失速和动态失速问题过程中,加入网格自适应算法,提高数值模拟对流动分离特性的捕捉和分辨能力.算例结果表明在流场发生失速后,运用本文的自适应算法能够在增加少量网格单元的情况下明显提高计算精度.     

Review of research into the concept of the microblowing technique for turbulent skin friction reduction

A new technology for reducing turbulent skin friction, called the Microblowing Technique (MBT), is presented. Results from proof-of-concept experiments show that this technology could potentially reduce turbulent skin friction by more than 50% of the skin friction of a solid flat plate for subsonic and supersonic flow conditions. The primary purpose of this review paper is to provide readers with information on the turbulent skin friction reduction obtained from many experiments using the MBT. Although the MBT has a penalty for obtaining the microblowing air associated with it, some combinations of the MBT with suction boundary layer control methods are an attractive alternative for a real application. Several computational simulations to understand the flow physics of the MBT are also included. More experiments and computational fluid dynamic (CFD) computations are needed for the understanding of the unsteady flow nature of the MBT and the optimization of this new technology.     

Aerodynamic shape optimization for alleviating dynamic stall characteristics of helicopter rotor airfoil

In order to alleviate the dynamic stall effects in helicopter rotor, the sequential quadratic programming(SQP) method is employed to optimize the characteristics of airfoil under dynamic stall conditions based on the SC1095 airfoil. The geometry of airfoil is parameterized by the class-shape-transformation(CST) method, and the C-topology body-fitted mesh is then automatically generated around the airfoil by solving the Poisson equations. Based on the grid generation technology, the unsteady Reynolds-averaged Navier-Stokes(RANS) equations are chosen as the governing equations for predicting airfoil flow field and the highly-efficient implicit scheme of lower–upper symmetric Gauss–Seidel(LU-SGS) is adopted for temporal discretization. To capture the dynamic stall phenomenon of the rotor more accurately, the Spalart–Allmaras turbulence model is employed to close the RANS equations. The optimized airfoil with a larger leading edge radius and camber is obtained. The leading edge vortex and trailing edge separation of the optimized airfoil under unsteady conditions are obviously weakened, and the dynamic stall characteristics of optimized airfoil at different Mach numbers, reduced frequencies and angles of attack are also obviously improved compared with the baseline SC1095 airfoil. It is demonstrated that the optimized method is effective and the optimized airfoil is suitable as the helicopter rotor airfoil.     

预测轴流压气机工作稳定性的非线性模型

建立了预测轴流压气机工作稳定性的非线性模型,采用有限差分数值模拟单转子和一台三级轴流压气机旋转失速的形成和发展过程,预测失速特征参数。由于合理地考虑了失速流动的非定常和非线性因素的影响,大多数失速特征参数的预测结果与试验测量值相符。     

Experimental investigation on SPS casing treatment with bias flow

Generally, casing treatment（CT） is a passivity method to enhance the stall margin of fan/compressor. A novel casing treatment based on the small disturbance theory and vortex and wave interaction suggestion is a method combining passive control and active control, which has been proved effective at enhancing the stall margin of fan/compressor in experiment. In order to investigate the mechanism of this kind of casing treatment, an experimental investigation of a stall precursor-suppressed（SPS） casing treatment with air suction or blowing air is conducted in the present paper. The SPS casing treatment is designed to suppressing stall precursors to realize stall margin enhancement in turbomachinery. The experimental results show that the casing treatment with blowing air of small quantity can improve the stall margin by about 8% with about 1% efficiency loss. By contrast, the SPS casing treatment with micro-bias flow does not improve the stall margin much more than that without bias flow, even worse. Meanwhile, the present investigation has also attempted to reveal the mechanism of stall margin improvement with the casing treatment.It is found that the stall margin improvements vary with the modification of the unsteady shedding flow and the unsteady wall boundary impedance. The experimental results agree fairly well with the theoretical prediction using a flow stability model of rotating stall.     

Challenges in high-alpha vehicle dynamics

The evolution of aerospace vehicles towards ever-increasing maneuverability, including flight at high angles of attack and vehicle motions of large amplitudes and high angular rates, has led to the need for prediction of vehicle aerodynamics that are dominated by unsteady separated flow effects. The existing data base is reviewed to determine to what degree the following critical issues are understood: 1. Cause and effect of asymmetric forebody flow separation with associated vortices. 2. Cause of slender wing rock. 3. Effect of vehicle motion on dynamic airfoil stall. 4. Extrapolation from subscale tests to full-scale free flight. To extend the present knowledge to include the coupling existing between novel aerodynamic controls and the vehicle dynamics is the challenge facing designers of future agile aircraft operating at high angles of attack.     

基于CFD/CSD紧耦合的旋翼翼型气弹分析

采用CFD/CSD（计算流体力学/计算结构力学）紧耦合的方法,以Fluent软件作为主控平台,通过UDF（用户自定义函数）及I/O（输入/输出）文件读写的方式实现结构响应和气动载荷的数据交换,耦合求解了旋翼桨叶剖面的气动力和振动响应.在此基础上研究旋翼桨叶剖面在变距、沉降（挥舞）和周期交变来流条件下的气动特性和振动响应特性.结果表明:桨叶剖面在轻失速情况下,气动载荷周期性比较好,表现出光滑的迟滞环曲线,结构沉降响应也表现出光滑的周期性现象,扭转响应出现局部轻微振荡.深失速情况下,气动载荷及结构响应都表现出强烈的非线性振荡,高频成分较为明显.      

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.     

Numerical investigation of aerodynamic flow actuation produced by surface plasma actuator on 2D oscillating airfoil

Numerical simulation of unsteady flow control over an oscillating NACA0012 airfoil is investigated. Flow actuation of a turbulent flow over the airfoil is provided by low current DC surface glow discharge plasma actuator which is analytically modeled as an ion pressure force produced in the cathode sheath region. The modeled plasma actuator has an induced pressure force of about 2 k Pa under a typical experiment condition and is placed on the airfoil surface at 0% chord length and/or at 10% chord length. The plasma actuator at deep-stall angles(from 5° to 25°) is able to slightly delay a dynamic stall and to weaken a pressure fluctuation in down-stroke motion. As a result, the wake region is reduced. The actuation effect varies with different plasma pulse frequencies, actuator locations and reduced frequencies. A lift coefficient can increase up to 70% by a selective operation of the plasma actuator with various plasma frequencies and locations as the angle of attack changes. Active flow control which is a key advantageous feature of the plasma actuator reveals that a dynamic stall phenomenon can be controlled by the surface plasma actuator with less power consumption if a careful control scheme of the plasma actuator is employed with the optimized plasma pulse frequency and actuator location corresponding to a dynamic change in reduced frequency.     

高前进比旋翼气动特性分析方法研究

提出一种前进比高达0.8的旋翼气动特性分析方法,该方法针对高前进比旋翼前行桨叶压缩性、后行桨叶失速效应严重以及桨叶偏流作用和反流区大的特点,建立了旋翼气动力模型以及与之相适应的旋翼诱导速度时变非均布模型与桨叶非定常挥舞运动模型,然后根据高前进比旋翼气动力、旋翼诱导速度和桨叶挥舞运动三者之间的内在耦合关系提出了高前进比旋翼气动特性的动态响应计算方法,最后以H-34旋翼为例计算了该旋翼高前进比状态的气动特性,并将计算结果与风洞试验数据进行对比验证,结果合理。     

Modern helicopter rotor aerodynamics

The helicopter rotor wake is among the most complex fluid dynamic structures being three dimensional and in many cases unsteady. The wake begins at the blade(s) where the flow can be transonic near the blade tip and undergo compressible dynamic stall. Farther down in the wake, the flow is essentially incompressible. Moreover, the rotor blades undergo complex unsteady motions because of the necessity to balance moments; they are elastic as well. In this paper, the fundamental aeromechanics of the wake and the flow on the blade is discussed and the primary methods of analysis, computation, and experiment employed to uncover the physics of the rotor wake are described.     

一种槽道式处理机匣的准定常流动模型Ⅰ: 建模方法及验证

为降低对槽道式处理机匣流动进行模拟的计算量,特别地,为了提供一种可用于处理机匣工程设计的快速评估手段,基于对带有槽道式处理机匣的压气机非定常流动物理的理解,提出了一种处理机匣的准定常流动模型,并利用带有处理机匣的跨声压气机转子实验测量结果,以及非定常数值模拟结果对所提出的模型进行了验证.结果表明:处理机匣流动模型很好地预测了由于采用处理机匣所取得的转子失速裕度的提高,并且与非定常模拟结果比较,流动模型对转子主流以及处理机匣内部流动模拟的准确性也能得到保证.采用处理机匣流动模型所需的计算量仅约为非定常模拟的1%,这保证了提出的模型可用于快速评估槽道式处理机匣的气动特性,从而为其工程优化设计提供了有效的手段.      

多片后缘小翼对直升机旋翼桨叶动态失速及桨毂振动载荷的控制

 减缓直升机后行桨叶动态失速发生、降低直升机桨毂振动载荷是提高直升机飞行速度、改进直升机飞行性能的重要途径。本文研究了直升机在高速高载情况下利用多片受控的桨叶后缘小翼对直升机的后行桨叶动态失速和桨毂振动载荷同时进行控制的有效方法。建立了弹性桨叶和后缘刚性小翼的结构动力学模型。桨叶剖面气动载荷采用Leishman-Beddoes 二维非定常动态失速模型计算,后缘小翼剖面气动载荷采用Hariharan-Leishman二维亚声速非定常气动模型计算。采用伽辽金和数值积分相结合的方法求解旋翼系统的气弹响应。建立了有效的多片后缘小翼控制策略和控制方法,分析了3片后缘小翼的运动规律及对后行桨叶动态失速和桨毂振动载荷的控制效果,结果表明利用多片小翼的运动是控制桨叶动态失速和桨毂振动载荷的有效方法。     

新一代战斗机非定常流动数值研究综述

新一代战斗机强调超机动能力和强隐身性，其中大攻角下的静态失速、动态失速及内埋弹仓绕流是与高机动和强隐身密切相关的、极具挑战的几类典型的非定常流动，它们对数值方法提出了极高的要求。为了高精度地仿真流场、清楚地揭示流动机理，有效地控制非定常流动，非常有必要发展高精度且高效率的RANS-LES混合方法体系，包含RANS-LES混合方法本身、与RANS-LES混合方法匹配的高精度自适应耗散格式、基准湍流模式、高质量计算网格、高精度时间推进方法、非定常量的统计方法等，具有极强的紧迫性。提出、发展、验证并应用该类方法数值仿真新一代战斗机（包括单独部件、组合体、甚至全机）的非定常流动，数值预测结果与风洞实验数据吻合良好；此类方法可为新型战斗机设计提供理论依据和分析手段。     

带气腔的离心压缩机旋转失速的三维数值模拟

使用商业CFD计算软件CFX求解三维雷诺平均的Navier-Stokes方程组, 结合出口气腔模型对某带无叶扩压器的离心压缩机的旋转失速现象进行数值模拟.首先使用定常计算得到了该离心压缩机的稳态性能曲线, 并和实验测量值进行了比较.然后引入出口气腔模型, 模拟离心压缩机内的旋转失速流动.在小流量下, 从矢量图和径向速度等值线图中观察得到了离心压缩机内部流场的非定常流动现象.还研究了气腔模型不同参数对失速流动的影响, 结果表明气腔体积越大, 计算得到的失速频率越低.