Simulation of Aerodynamic Performances of Flexible Flapping Wing Airfoils北大核心CSCD

与固定翼相比,在低速、小Reynolds数条件下,扑翼飞行具有显著的气动性能优势,受到越来越多的重视。然而,目前对扑翼翼型的研究以刚性翼型为主,对柔性翼型气动性能认识还不清楚。该文建立了柔性椭圆翼型的流固耦合仿真模型,分析了不同风速、迎角下柔性椭圆翼型的周围流场、变形以及气动性能。仿真结果表明,较刚性翼型,柔性翼型延缓了尾涡脱落时间,有效降低升力扰动振荡频率;柔性翼型显著抑制了尾流流场的扰动,降低升力扰动振荡幅值,合适的弹性模量翼型使得扰动振荡完全消除。研究结果可为软飞行器气动设计提供参考。     

Multi-fidelity design optimization of transonic airfoils using physics-based surrogate modeling and shape-preserving response prediction

A computationally efficient design methodology for transonic airfoil optimization has been developed. In the optimization process, a numerically cheap physics-based low-fidelity surrogate (the transonic small-disturbance equation) is used in lieu of an accurate, but computationally expensive, high-fidelity (the compressible Euler equations) simulation model. Correction of the low-fidelity model is achieved by aligning its corresponding airfoil surface pressure distribution with that of the high-fidelity model using a shape-preserving response prediction technique. The resulting method requires only a single high-fidelity simulation per iteration of the design process. The method is applied to airfoil lift maximization in two-dimensional inviscid transonic flow, subject to constraints on shock-induced pressure drag and airfoil cross-sectional area. The results showed that more than a 90% reduction in high-fidelity function calls was achieved when compared to direct high-fidelity model optimization using a pattern-search algorithm.     

A comparison of the flow around a flexible and a non-flexible airfoil for flapping wing propulsion

The flow around flapping airfoils is evaluated using the Stereoscopic Particle Image Velocimetry (PIV) technique. Resulting from the measurements is the distribution of the turbulent shear stress. This serves as an indicator for the position of transition from laminar to turbulent flow. The distributions were compared for both a flapping non–flexible airfoil and a flexible airfoil. It was determined that there are differences between the flow around the flexible and the non–flexible airfoils. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Numerical simulation of airfoil vibrations induced by turbulent flow

The subject of this paper is the numerical simulation of the interaction between two-dimensional incompressible viscous flow and a vibrating airfoil. A solid elastically supported airfoil with two degrees of freedom, which can rotate around the elastic axis and oscillate in the vertical direction, is considered. The numerical simulation consists of the stabilized finite element solution of the Reynolds averaged Navier–Stokes equations with algebraic models of turbulence, coupled with the system of ordinary differential equations describing the airfoil motion. Since the computational domain is time dependent and the grid is moving, the Arbitrary Lagrangian–Eulerian (ALE) method is used. The developed method was applied to the simulation of flow-induced airfoil vibrations.     

The new airfoil model NACA0015, modal analysis and flutter properties

The experimental airfoil model NACA0015 was used to study aeroelastic phenomena during self-excited profile vibration. It provides data for control of aeroelastic calculation programs at subsonic speeds of the stream. The model movability is two-dimensional with two-degree of freedom dynamic system, one in pitch and the second in plunge, and is proposed to be a dynamic system having two near corresponding eigenfrequencies. To quantitatively evaluate flow field using interferometry, a special test section design and profile was constructed. It utilized a large visual field for the optical system together with the option of changing support stiffness for both degrees of freedom. In this paper experimental results from the range of Reynolds numbers Re = (2.63–2.83) 10⁵ are published. The identified eigenvalues and eigenmodes for zero flow velocity are compared with measured flutter properties (frequency, modes and time evolutions) of the airfoil.     

Numerical bifurcation analysis of static stall of airfoil and dynamic stall under unsteady perturbation

By the finite element method combined with Arbitrary-Lagrangian-Eulerian (ALE) frame and explicit Characteristic Based Split Scheme (CBS), the complex flows around stationary and sinusoidal pitching airfoil are studied numerically. In particular, the static and dynamic stalls are analyzed in detail, and the natures of the static stall of NACA0012 airfoil are given from viewpoint of bifurcations. Following the bifurcation in Map, the static stall is proved to be the result from saddle-node bifurcation which involves both the hysteresis and jumping phenomena, by introducing a Map and its Floquet multiplier, which is constructed in the numerical simulation of flow field and related to the lift of the airfoil. Further, because the saddle-node bifurcation is sensitive to imperfection or perturbation, the airfoil is then subjected to a perturbation which is a kind of sinusoidal pitching oscillation, and the flow structure and aerodynamic performance are studied numerically. The results show that the large-scale flow separation at the static stall on the airfoil surface can be removed or delayed feasibly, and the ensuing lift could be enhanced significantly and also the stalling incidence could be delayed effectively. As a conclusion, it can be drawn that the proper external excitation can be considered as a powerful control strategy for the stall. As an unsteady aerodynamic behavior of high angle of attack, the dynamic stall can be investigated from viewpoint of nonlinear dynamics, and there exists a rich variety of nonlinear phenomena, which are related to the lift enhancement and drag reduction.     

A modified discrete-vortex method algorithm with shedding criterion for aerodynamic coefficients prediction at high angle of attack

Low-order methods require less computing power than classical computational fluid dynamics and can be implemented on a laptop computer, which is needed for engineering tasks. Discrete vortex methods are such low order methods that can describe the unsteady separated flow around an airfoil. After a presentation of the leading edge suction parameter discrete vortex method, a modified algorithm is proposed, in order to reduce the computing cost, and compared with the previous one. Several reference unsteady airfoil motions are discussed in terms of gain in the computation time with comparisons between the previous scheme and the present one. The accuracy of the new method is demonstrated through aerodynamic coefficients. The application of the present discrete vortex method to a transient pitching motion of an airfoil is also presented, in order to understand the leading edge vortex formation, and its implication in terms of lift and drag coefficients. The method is not limited to unsteady or transient motions but can also simulate the flow around a constant angle of attack airfoil. In that case, an original method of fast summation of the vortices located far away from the airfoil, allows a linear dependence of the computation time versus the number of vortices shed, which is a great improvement over the quadratic dependence observed in the classical discrete vortex methods. The development of the aerodynamic coefficients with angle of attack, from values ranging between −10° and 90°, is obtained for a purely two-dimensional flow. In particular, the shape of the lift coefficient of the airfoil in the fully detached flow region is established. Comparisons with relevant experimental or computational fluid dynamics data are discussed in order to grasp the influence of upstream turbulence level and three-dimensional effects in the measured data in the fully detached flow region.     

Application of finite element method in aeroelasticity

The subject of this paper is the numerical simulation of the interaction of two-dimensional incompressible viscous flow and a vibrating airfoil, which can rotate around the elastic axis and oscillate in the vertical direction. The numerical simulation consists of the finite element approximation of the Navier–Stokes equations coupled with the system of ordinary differential equations describing the airfoil motion. The arbitrary Lagrangian–Eulerian (ALE) formulation of the Navier–Stokes equations, stabilization the finite element discretization and coupling of both models is discussed. Moreover, the Reynolds averaged Navier–Stokes (RANS) system of equations together with the Spallart–Almaras turbulence model is also discussed. The computational results of aeroelastic calculations are presented and compared with the NASTRAN code solutions.     

Computational fluid dynamics simulation of aerodynamic performance and flow separation by single element and slatted airfoils under rainfall conditions

This study investigated the effects of rainfall on flow separation and the aerodynamic performance of single element and slatted NACA 0012 airfoils by using a mathematical model developed with the commercial computational fluid dynamics solver ANSYS FLUENT 18.2. A two-way momentum coupled Eulerian–Lagrangian multiphase approach was used to simulate the formation of the water film layer on the airfoil's surface. According to the results, very low values of the lift-to-drag ratio at low angles of attack reflected severe degradation of the aerodynamic performance of the airfoil in the presence of water accumulated on its surface. The impact of rain droplets on the leading-edge slat surface led to less water accumulating on the main section of the airfoil. In particular, the maximum water film mass concentrated on the airfoil surface decreased from 15 g to 1 g compared with the single element airfoil. Hence, the thickness of the water film layer was not sufficiently large to significantly affect the aerodynamic coefficients of the slatted airfoil, especially the maximum lift coefficient, compared with the thicker water film layer on the single element airfoil. In addition, the use of slats clearly enhanced the aerodynamic coefficients and increased the stall angle from 13° to 22° in dry conditions, and from 16° to 24° in rainy conditions. Slats also significantly decreased the boundary layer thickness and delayed the separation at higher angles of attack.     

Aerodynamic Optimization of Airfoils Using Adaptive Parameterization and Genetic Algorithm

A new method for airfoil shape parameterization is presented, and its influences on the optimum design and convergence of the evolutionary optimization process are investigated. An online adaptive method is used that alters the airfoil parametric function during the process of optimization. A geometric inverse design is carried out, and the capability of the method for producing general airfoil shapes is assessed. The performance of the method is then evaluated by aerodynamic shape optimization. The result indicates that the proposed method improves the optimum design airfoil significantly. In addition, it reduces the total number of flow solver calls, which consequently reduces the required computational time.     

An insight into the rotational stall delay

Blade element momentum (BEM) theory which is based on the two-dimensional (2D) aerodynamic properties of airfoil blade element is the most common computational engineering method for the prediction of loads and power curves of wind turbines. Although most BEM models yield acceptable results for high tip-speed ratios where the local angles of attack are small, no generally accepted model exists up to date that consistently predicts the loads and power in stall regime for stall-controlled turbines. Understanding of the stall delay phenomenon on wind turbines remains, to this day, incomplete. The lack of a conceptual model for the complex three-dimensional (3D) flow field on the rotor blade, where stall is begins, how it progresses and where stall is practically terminated, has hindered the finding of a unanimously accepted solution. The paper aims at giving a better understanding of the delayed stall events and a reasonably simple correction model that complements the 2D airfoil characteristics used to a BEM method. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

High-Order Multioperator Compact Schemes for Numerical Simulation of Unsteady Subsonic Airfoil Flow

On the basis of high-order schemes, the viscous gas flow over the NACA2212 airfoil is numerically simulated at a free-stream Mach number of 0.3 and Reynolds numbers ranging from 10³ to 10⁷. Flow regimes sequentially varying due to variations in the free-stream viscosity are considered. Vortex structures developing on the airfoil surface are investigated, and a physical interpretation of this phenomenon is given.     

Solving the Inverse Mixed Boundary-Value Problem of Aerohydrodynamics in a New Statement

We consider the following inversemixed boundary-value problem of aerohydrodynamics. It is required to find a form of an airfoil circulated by a potential flow of an incompressible non-viscous liquid. A part of the profile is known; it is a broken line, convex upwards. The other part is found by the values of velocity potential, given as a function of parameter which can be either the abscissa, or the ordinate of an airfoil point. On a small part of the airfoil, containing the leading edge, the parameter is the ordinate, and for other points of the unknown part it is the abscissa.     

Turbulent Flow with Embedded Vortical Structures Induced by Vortex Generators in a Cascade

Presented work is the next step after several experimental examinations of vortex generator influence on a flow separation occurring on a model of the NACA 63A421 airfoil with deflected simple flap. In this stage of research the vortices produced by vortex generators (VGs) were studied using Particle Image Velocimetry technique (PIV) and numerical simulations. Vane type VGs with two spacings among VGs pairs in straight channel with turbulent flow were tested. The average velocity flow field, peak of vorticity and circulation decay downstream of VGs were evaluated. (© 2012 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

A two-algebraic turbulence model for compressible flow in turbomachinery cascade

The objective of this study is to examine the Baldwin–Lomax turbulence model in a finite volume solver to introduce a computer code for complex two–dimensional flows in turbomachinery. The turbulent model was tested with investigating the turbulent flow over a flat plate and other test cases. Three test cases are presented and the computed results are compared with experimental data. The calculated velocity profile agreed well with the experimental data in plate test case and the solver is validated in test case of flow over a semi NACA–0012 airfoil. The solver is used for flow through a multi–blade cascade of an axial compressor in design condition to show its capability of multi–block solution. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Application in Aerohydrodynamics of the Solution of an Inverse Boundary Value Problem for Analytic Functions

We consider a modified inverse boundary value problem of aerohydrodynamics in which it is required to find the shape of an airfoil streamlined by a potential flow of an incompressible nonviscous fluid, when the distribution of the velocity potential on one section of the airfoil is given as a function of the abscissa, while, on other sections of the airfoil, as a function of the ordinate of the point. The velocity of the undisturbed flow streamlining the sought-for airfoil is determined in the process of solving the problem. It is shown that, under rather general conditions on the initially set functions, the sought-for contour is closed unlike the inverse problem in the case when, on the unknown contour, the velocity distribution is given as a function of the arc abscissa of the contour point. We also consider the case when, on the entire desired contour, the distribution of the velocity potential is given as a function of one and the same Cartesian coordinate of the contour point.     

风力机叶型增升及叶尖流动控制研究

研究了两种改善风力机叶型气动性能的流动控制技术,分别对风力机专用S809翼型和较大升阻比的FX 60-100翼型进行应用研究．首先,通过在叶型前缘加装流动偏转器,研究流动偏转器对叶型流动分离的控制效果．并采用多岛基因算法,对流动偏转器进行多参数优化．结果表明:流动偏转器可以有效控制叶型的失速特性,推迟失速攻角和增加升力;基因优化算法能更大地提升流动偏转器的控制效果．其次,基于对风力机叶尖旋涡和尾涡特征以及叶片表面压力分布的分析,在叶片尖部加装不同倾斜角的旋涡扩散器控制叶尖涡．结果表明:涡扩散器能够提高叶尖涡涡核的总压,削弱其旋涡强度,使风力机尾流旋涡耗散更快,从而可以减小噪声,提高叶片效率．     

翼涡干扰前缘开孔被动控制数值研究

开孔方法是一种简单的流动被动控制方法.为找到一种有效降低桨涡干扰效应的被动控制方法,以NACA 0012翼型作为研究对象,建立了4种前缘开孔的模型.在不同来流速度、涡的强度和干扰距离条件下,对4种前缘开孔模型和无孔的基准翼型进行了二维平行桨涡干扰（翼涡干扰）数值模拟,对比了升力系数的变化.结果表明:前缘开孔可以降低翼涡干扰效应,但对翼型升力系数有一定的影响;宽度为2.5%弦长的直孔能在翼型升力系数损失较小的情况下有效地降低翼涡干扰效应,且适用范围较广.     

Design of a subsonic airfoil with upstream blowing

The problem is solved of designing a symmetric airfoil with upstream blowing opposite to subsonic irrotational steady flow of an inviscid incompressible fluid. The solution relies on Sedov’s idea of a stagnation region developing in the neighborhood of the stagnation point. An iterative solution process is developed, and examples of airfoils are constructed. The numerical results are analyzed, and conclusions are drawn about the effect of blowing parameters on the airfoil geometry and the resultant force acting on the airfoil.     

On Transonic Marginal Separation

Laminar boundary‐layer separation near the leading edge of a thin airfoil is one of the principal factors that limits the lift force acting on the airfoil. Marginal separation, in particular, denotes the onset of separation which is accompanied by the formation of a short separation bubble. Using asymptotic analysis this effect is studied in the limit of high Reynolds number and for transonic external flow conditions. It is assumed that the fluid under consideration is a perfect gas and the airfoil surface is taken to be thermally insulated. Results to be presented include the analytical investigation of the emerging three layer structure, the associated transonic far field and the calculation of representative wall shear stress distributions. (© 2004 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)