Flow field characteristics study of a flapping airfoil using computational fluid dynamics

The flow field of a flapping airfoil in Low Reynolds Number (LRN) flow regime is associated with complex nonlinear vortex shedding and viscous phenomena. The respective fluid dynamics of such a flow is investigated here through Computational Fluid Dynamics (CFD) based on the Finite Volume Method (FVM). The governing equations are the unsteady, incompressible two-dimensional Navier-Stokes (N-S) equations. The airfoil is a thin ellipsoidal geometry performing a modified figure-of-eight-like flapping pattern. The flow field and vortical patterns around the airfoil are examined in detail, and the effects of several unsteady flow and system parameters on the flow characteristics are explored. The investigated parameters are the amplitude of pitching oscillations, phase angle between pitching and plunging motions, mean angle of attack, Reynolds number (Re), Strouhal number (St) based on the translational amplitudes of oscillations, and the pitching axis location (x/c). It is shown that these parameters change the instantaneous force coefficients quantitatively and qualitatively. It is also observed that the strength, interaction, and convection of the vortical structures surrounding the airfoil are significantly affected by the variations of these parameters.     

低雷诺数俯仰振荡翼型等离子体流动控制

针对低雷诺数翼型气动性能差的特点, 通过介质阻挡放电(dielectric barrier discharge, DBD)等离子体激励控制的方法, 提高翼型低雷诺数下的气动特性,改善其流场结构. 采用二维准直接数值模拟方法求解非定常不可压Navier-Stokes方程,对具有俯仰运动的NACA0012翼型的低雷诺数流动展开数值模拟.同时将介质阻挡放电激励对流动的作用以彻体力源项的形式加入Navier-Stokes方程,通过数值模拟探究稳态DBD等离子体激励对俯仰振荡NACA0012翼型气动特性和流场特性的影响.为了进行流动控制, 分别在上下表面的前缘和后缘处安装DBD等离子体激励器,并提出四种激励器的开环控制策略,通过对比研究了这些控制策略在不同雷诺数、不同减缩频率以及激励位置下的控制效果.通过流场结构和动态压强分析了等离子体进行流场控制的机理. 结果表明,前缘DBD控制中控制策略B(负攻角时开启上表面激励器,正攻角时开启下表面激励器)效果最好,后缘DBD控制中控制策略C(逆时针旋转时开启上表面激励器,顺时针旋转时开启下表面激励器)效果最好,前缘DBD控制效果会随着减缩频率的增大而下降, 同时会导致阻力增大.而后缘DBD控制可以减小压差阻力, 优于前缘DBD控制,对于计算的所有减缩频率(5.01~11.82)都有较好的增升减阻效果.在不同雷诺数下, DBD控制的增升效果较为稳定, 而减阻效果随着雷诺数的降低而变差,这是由流体黏性效应增强导致的.      

High fidelity numerical simulation of airfoil thickness and kinematics effects on flapping airfoil propulsion

High-fidelity numerical simulations with the spectral difference (SD) method are carried out to investigate the unsteady flow over a series of oscillating NACA 4-digit airfoils. Airfoil thickness and kinematics effects on the flapping airfoil propulsion are highlighted. It is confirmed that the aerodynamic performance of airfoils with different thickness can be very different under the same kinematics. Distinct evolutionary patterns of vortical structures are analyzed to unveil the underlying flow physics behind the diverse flow phenomena associated with different airfoil thickness and kinematics and reveal the synthetic effects of airfoil thickness and kinematics on the propulsive performance. Thickness effects at various reduced frequencies and Strouhal numbers for the same chord length based Reynolds number (=1200) are then discussed in detail. It is found that at relatively small Strouhal number (=0.3), for all types of airfoils with the combined pitching and plunging motion (pitch angle 20°, the pitch axis located at one third of chord length from the leading edge, pitch leading plunge by 75°), low reduced frequency (=1) is conducive for both the thrust production and propulsive efficiency. Moreover, relatively thin airfoils (e.g. NACA0006) can generate larger thrust and maintain higher propulsive efficiency than thick airfoils (e.g. NACA0030). However, with the same kinematics but at relatively large Strouhal number (=0.45), it is found that airfoils with different thickness exhibit diverse trend on thrust production and propulsive efficiency, especially at large reduced frequency (=3.5). Results on effects of airfoil thickness based Reynolds numbers indicate that relative thin airfoils show superior propulsion performance in the tested Reynolds number range. The evolution of leading edge vortices and the interaction between the leading and trailing edge vortices play key roles in flapping airfoil propulsive performance.     

Efficiency of an auto-propelled flapping airfoil

The present study deals with an investigation of the flow aerodynamic characteristics and the propulsive velocity of a system equipped with a nature inspired propulsion system. In particular, the study is aimed at studying the effect of the flapping frequency on the flow behavior. We consider a NACA0014 airfoil undergoing a vertical sinusoidal flapping motion. In contrast to nearly all previous studies in the literature, the present work does not impose any velocity on the inlet flow. During each iteration the outer flow velocity is computed after having determined the forces exerted on the airfoil. Forward motion may only be produced by flapping motion of the airfoil. This is more consistent with the physical phenomenon. The non-stationary viscous flow around the flapping airfoil is simulated using Ansys-Fluent 12.0.7. The airfoil movement is achieved using the deformable mesh technique and an in-house developed User Define Function (UDF). Our results show the influence of flapping frequency and amplitude on both the airfoil velocity and the propulsive efficiency. The resulting motion is contrasts to the applied forces. In the present study, the frequency ranges from 0.1 to 20 Hz while the airfoil amplitude values considered are: 10%, 17.5%, 25% and 40%.     

Power extraction efficiency improvement of a fully-activated flapping foil: With the help of an auxiliary rotating foil

The power extraction efficiency improvement of a fully-activated flapping foil with the help of an auxiliary rotating foil is numerically examined in this work. A NACA0015 airfoil is placed in a two-dimensional laminar flow and synchronously performs the imposed pitching and plunging motions. An auxiliary smaller foil, which rotates about its center, is arranged below the flapping foil. As a consequence, the vortex interaction between the flapping foil and the rotating foil occurs. At a Reynolds number of 1100 and the position of the pitching axis at one-third chord, the effects of the distance between the flapping foil and the auxiliary foil, the phase difference between the rotating motion and the flapping motion as well as the frequency of flapping motion on the power extraction performance are systematically investigated. It is found that compared to the single flapping foil, the efficiency of power extraction for the flapping foil with an auxiliary device can be improved. Based on the numerical analysis, it is indicated that the enhanced plunging contribution, which is caused by the increased lift force owing to the vortex interaction, directly helps the efficiency improvement.     

Numerical investigation of transonic flow over deformable airfoil with plunging motion

In this article, the transonic inviscid flow over a deformable airfoil with plunging motion is studied numerically. A finite volume method based on the Roe scheme developed in a generalized coordinate is used along with an arbitrary Lagrangian-Eulerian method and a dynamic mesh algorithm to track the instantaneous position of the airfoil.The effects of different governing parameters such as the phase angle, the deformation amplitude, the initial angle of attack, the flapping frequency, and the Mach number on the unsteady flow field and aerodynamic coefficients are investigated in detail. The results show that maneuverability of the airfoil under various flow conditions is improved by the deformation. In addition, as the oscillation frequency of the airfoil increases, its aerodynamic performance is significantly improved.     

Computational study of a pitching bio-inspired corrugated airfoil

The phenomenon of insect flight has been of scientific interest for many years and is more recently inspiring modern engineering devices such as Micro Aerial Vehicles (MAVs). Insect flight is characterized by unsteady fluid dynamics at low Reynolds numbers. The importance of viscous effects to the successful flapping flight of insects has been identified and with the current state of computing and Computational Fluid Dynamics (CFD) these effects can now be studied in detail. The present work attempts to simplify this complex phenomenon by considering symmetric oscillating rotational motion of a wing (pitching). What is of interest in this study is how the shape of a corrugated idealized insect wing affects the performance and flow characteristics around the pitching wing. Two dimensional CFD on an oscillating wing has been performed and reported. Measurements were taken to ensure the accuracy of the computational solution and the results validated against experimental PIV results. A range of frequencies and rotational amplitudes have been investigated. Lift and drag coefficients have been analyzed for all cases to quantify the effects of unsteady flow features on the performance of the oscillating wing. It was found that the wing shape used in this study resulted in the viscous features formed on the top of the wing exhibiting high sensitivity to the oscillating conditions and these influenced the performance of the wing. The flow features formed on the bottom of the wing remained similar throughout the cases tested. In the pitching regime this wing profile did not perform as well as published results for smooth airfoils in terms of thrust and propulsion efficiency. However this may be due to reduced frequency effects becoming important at our high pitching amplitude which need to be investigated further. There may be other oscillatory regimes that more accurately represent flapping flight in which the corrugated foil outperforms a smooth counterpart but these are yet to be investigated. Further research in this area may help answer the question as to how evolutionarily significant other benefits of a corrugated wing, such as being light and strong, are compared to its aerodynamic properties, the present results seem to favor the former.     

An investigation into the effects of unsteady parameters on the aerodynamics of a low Reynolds number pitching airfoil

The growing applications of low Reynolds number (LRN) operating vehicles impose the need for accurate LRN flow solutions. These applications usually involve complex unsteady phenomena, which depend on the kinematics of the vehicle such as pitching, plunging, and flapping of a wing. The objective of the present study is to address the issues related to LRN aerodynamics of a harmonically pitching NACA0012 airfoil. To this end, the influence of unsteady parameters, namely, amplitude of oscillation, d, reduced frequency, k, and Reynolds number, Re, on the aerodynamic performance of the model is investigated. Computational fluid dynamics (CFD) is utilized to solve Navier–Stokes (N–S) equations discretized based on the Finite Volume Method (FVM). The resulting instantaneous lift coefficients are compared with analytical data from Theodorsen’s method. The simulation results reveal that d, k, and Re are of great importance in the aerodynamic performance of the system, as they affect the maximum lift coefficients, hysteresis loops, strength, and number of the generated vortices within the harmonic motion, and the extent of the so-called figure-of-eight phenomenon region. Thus, achieving the optimum lift coefficients demands a careful selection of these parameters.     

Numerical analysis on transitions and symmetry-breaking in the wake of a flapping foil

Flying and marine animals often use flapping wings or tails to generate thrust. In this paper, we will use the simplest flapping model with a sinusoidal pitching motion over a range of frequency and amplitude to investigate the mechanism of thrust generation. Previous work focuses on the Karman vortex street and the reversed Karman vortex street but the transition between two states remains unknown. The present numerical simulation provides a complete scenario of flow patterns from the Karman vortex street to reversed Karman vortex street via aligned vortices and the ultimate state is the deflected Karman vortex street, as the parameters of flapping motions change. The results are in agreement with the previous experiment. We make further discussion on the relationship of the observed states with drag and thrust coefficients and explore the mechanism of enhanced thrust generation using flapping motions.     

飞行器纵向俯仰气动特性计算研究

应用有限体积方法求解三维可压缩雷诺平均N-S方程,计算了巡航导弹外形飞行器作小振幅俯仰运动时的动态绕流流场和空气动力特性,开展了导弹绕不同转轴、以不同频率和在不同迎角范围内进行俯仰运动的非定常气动力迟滞特性研究。计算结果表明,当导弹作快速俯仰运动时,在上仰和下俯过程中的同一迎角瞬间,绕导弹流场流动明显不同,表现出明显的非定常迟滞特性。导弹的非定常气动力迟滞特性随俯仰运动频率的增大明显增强,且气动力迟滞曲线随着俯仰轴位置的变化而变化。在同一减缩频率下,导弹在不同迎角范围内作周期俯仰运动时,相同的运动相位角所对应的升力系数对迎角的导数是一致的,而不同减缩频率下升力系数对迎角的导数随运动相位角变化曲线明显不同。     

Numerical study of large amplitude,nonsinusoidal motion and camber effects on pitching airfoil propulsion

The effects of large amplitude and nonsinusoidal motion on pitching airfoil aerodynamics for thrust generation were numerically studied with a 2-D NACA0012 airfoil used, and various 2-D NACA asymmetric airfoils were applied for camber effect study. The large amplitude effect study has been undertaken over a wide range of reduced frequency k (from 6 to 18) and pitching amplitude θ (from 5° to 30°) at Re=1.35×10⁴ with sinusoidal pitching profile used. It is shown that the large pitching amplitude results in much more thrust generated than that at low pitching amplitude and the increase of thrust with amplitude becomes slow when the amplitude reaches some degree. However, the propulsive efficiency noticeably decreases with the increase of θ at a fixed k.An adjustable parameter K was employed to realize various nonsinusoidal motions and the effect of nonsinusoidal motion was investigated with various unsteady parameters (θ, k) applied. The results reveal that nonsinusoidal motion has a noticeable effect on the aerodynamic performance, as it affects the instantaneous force coefficients, maximum thrust coefficients and flow structures. An increase in K results in a better thrust generation performance at fixed θ and k, especially for K>0. It is also shown that the larger K noticeably influences the wake pattern and induces a stronger reverse von Karman vortex street in the wake, which in turn leads to the increased thrust. The camber study was performed on various 2-D NACA airfoils with different cambers and camber locations undergoing sinusoidal pitching motion at θ=5° and Re=1.35×10⁴. It is found that varying camber offers little improvement in thrust generation performance.     

Thrust augmentation of flapping airfoils in low Reynolds number flow using a flexible membrane

The unsteady aerodynamic thrust and aeroelastic response of a two-dimensional membrane airfoil under prescribed harmonic motion are investigated computationally with a high-order Navier–Stokes solver coupled to a nonlinear membrane structural model. The effects of membrane prestress and elasticity are examined parametrically for selected plunge and pitch–plunge motions at a chord-based Reynolds number of 2500. The importance of inertial membrane loads resulting from the prescribed flapping is also assessed for pure plunging motions. This study compares the period-averaged aerodynamic loads of flexible versus rigid membrane airfoils and highlights the vortex structures and salient fluid–membrane interactions that enable more efficient flapping thrust production in low Reynolds number flows.     

Effect of flexure on aerodynamic propulsive efficiency of flapping flexible airfoil

The aim of present study is to investigate the effect of chord-wise flexure amplitude on unsteady aerodynamic characteristics for a flapping airfoil with various combinations of Reynolds number and reduced frequency. Unsteady, viscous flows over a single flexible airfoil in plunge motion are computed using conformal hybrid meshes. The dynamic mesh technique is applied to illustrate the deformation modes of the flexible flapping airfoil. In order to investigate the influence of the flexure amplitude on the aerodynamic performance of the flapping airfoil, the present study considers eight different flexure amplitudes (a₀) ranging from 0 to 0.7 in intervals of 0.1 under conditions of Re=10⁴, reduced frequency k=2, and dimensionless plunge amplitude h₀=0.4. The computed unsteady flow fields clearly reveal the formation and evolution of a pair of leading edge vortices along the body of the flexible airfoil as it undergoes plunge motion. Thrust-indicative wake structures are generated when the flexure amplitude of the airfoil is less than 0.5 of the chord length. An enhancement in the propulsive efficiency is observed for a flapping airfoil with flexure amplitude of 0.3 of the chord length. This study also calculates the propulsive efficiency and thrust under various Reynolds numbers and reduced frequency conditions. The results indicate that the propulsive efficiency has a strong correlation with the reduced frequency. It is found that the flow conditions which yield the highest propulsive efficiency correspond to Strouhal number St of 0.255.     

Numerical investigations into the asymmetric effects on the aerodynamic response of a pitching airfoil

The effects of asymmetric sinusoidal motion on pitching airfoil aerodynamics were studied by numerical simulations for 2-D flow around a NACA0012 airfoil at Re=1.35×10⁵. Various unsteady parameters (amplitude of oscillation, d; reduced frequency, k) were applied to investigate the effect of asymmetry parameter S on the instantaneous force coefficients and flow patterns. The results reveal that S has a noticeable effect on the aerodynamic performance, as it affects the instantaneous force coefficient, maximum lift and drag coefficient, hysteresis loops and the flow structures.     

Comparative Analysis of the Characteristics of the Vortex Wake behind a Flapping Wing Performing Oscillations of Different Types

The wake characteristics of a custom-designed airfoil performing pitching oscillations, heaving oscillations, and a combination of pitch and heave oscillations are compared in this study. The influence of flapping parameters are investigated at a constant Reynolds number Re\(_{c} = 2640\) and is presented for the Strouhal numbers based on the oscillation amplitude, St_A, varying in the \(0.1 \leqslant {\text{S}}{{{\text{t}}}_{A}} \leqslant 0.4\) range. The generation of vorticity above and below the airfoil depends on the airfoil’s initial direction of motion and remains the same for all types of flapping oscillations investigated. The evolution of the leading-edge and trailing-edge vortices is presented. The heaving oscillations of the airfoil are found to have a greater influence on the characteristics of the leading edge vortex. The wake behind the combined pitch-heave oscillations appears to be governed by pitching oscillations below \({\text{S}}{{{\text{t}}}_{A}} = 0.24\), whereas it is driven by heaving oscillations above \({\text{S}}{{{\text{t}}}_{A}} = 0.24\). The force computations indicate that the mere existence of the reverse von Kármán street is not sufficient to develop the thrust on the airfoil. The periodic component of velocity fluctuations significantly influences the wake characteristics. The anisotropic stress field developed around the airfoil due to the periodic fluctuations of the velocity is presented. The coherent structures developed in the wake are identified using the proper orthogonal decomposition and a qualitative comparison of the structures for different flapping oscillations is presented. The energy transfer from the flapping airfoil to the fluid for different flapping oscillations is highest for heaving oscillations followed by combined pitch-heave oscillations and pitching oscillations.

     

How a flexible tail improves the power extraction efficiency of a semi-activated flapping foil system: A numerical study

The improvement of power extraction of a semi-activated flapping foil system via the use of a flexible tail is numerically investigated in this work. A NACA0015 airfoil arranged in a two-dimensional laminar flow synchronously executes a forced pitching motion and an induced plunging motion. A flat plate attached to the trailing edge of the foil is utilized to model a tail, and thereby they are considered as a unit for the purpose of power extraction. The tail is either rigid or deformable. At a Reynolds number of 1100 and the position of the pitching axis at third chord, the effects of the mass and flexibility of the tail as well as the frequency of pitching motion on the net power extraction are systematically examined. It is found that compared to the foil with a rigid tail, the efficiency of net power extraction for the foil with a deformable tail can be improved. Based on the numerical analysis, it is indicated that the enhanced power extraction, which is caused by the increased lift force, directly contributes to the net efficiency improvement. Moreover, owing to high enhancement of power extraction, a flexible tail with high flexibility is recommended in the semi-activated flapping foil based power extraction system.     

An analytic approach to theoretical modeling of highly unsteady viscous flow excited by wing flapping in small insects

Numerous studies on the aerodynamics of insect wing flapping were carried out on different approaches of flight investigations, model experiments, and numerical simulations, but the theoretical modeling remains to be explored. In the present paper, an analytic approach is presented to model the flow interactions of wing flapping in air for small insects with the surrounding flow fields being highly unsteady and highly viscous. The model of wing flapping is a 2-D flat plate, which makes plunging and pitching oscillations as well as quick rotations reversing its positions of leading and trailing edges, respectively, during stroke reversals. It contains three simplified aerodynamic assumptions: (i) unsteady potential flow; (ii) discrete vortices shed from both leading and trailing edges of the wing; (iii) Kutta conditions applied at both edges. Then the problem is reduced to the solution of the unsteady Laplace equation, by using distributed singularities, i.e., sources/sinks, and vortices in the field. To validate the present physical model and analytic method proposed via benchmark examples, two elemental motions in wing flapping and a case of whole flapping cycles are analyzed, and the predicted results agree well with available experimental and numerical data. This verifies that the present analytical approach may give qualitatively correct and quantitatively reasonable results. Furthermore, the total fluid-dynamic force in the present method can be decomposed into three parts: one due to the added inertial (or mass) effect, the other and the third due to the induction of vortices shed from the leading-and the trailing-edge and their images respectively, and this helps to reveal the flow control mechanisms in insect wing flapping. The project supported by the National Natural Science Foundation of China (10072066) and the Chinese Academy of Sciences (KJCX-SW-LO4, KJCX2-SW-L2)     

Experimental analysis of passive bio-inspired covert feathers for stall and post-stall performance enhancement

Technologies inspired by the functioning and behavior of biological beings are commonly developed for aircraft flight. Among the bio-inspired approaches that have grown in interest, particularly for unmanned aerial vehicle flight, is based on the behavior of bird’s cover feathers under higher angles of attack. The covert feathers, when activated by separated flows, promote lift increment that helps in certain maneuvers. This work investigates the benefit in the stall and post-stall performance of employing bio-inspired covert feathers devices attached to an airfoil’s upper surface. To fill the gaps in the recent technical literature, experimental analysis of an SD7003 airfoil was executed in a wind tunnel with the application of bio-inspired covert feathers of different shapes and tapes in three chordwise positions. The bio-inspired devices were conceived to resemble the feathers’ lightness and discrete-distribution along with the wing model. Experiments were carried out measuring the aerodynamic forces and moment at Reynolds number around 170,000 for static and dynamic ramp-up and hold pitching motion. It has been confirmed that the use of bio-inspired covert feathers brought benefits to the stall and post-stall behavior of the airfoil. The maximum lift has increased, and the transition from attached to stalled flow around the airfoil tends to be smoother when the devices were used. Four shapes for the bio-inspired devices and three positions in chordwise direction were considered. The best performance among the case was encountered for a jagged bio-inspired device taped at a quarter-chord position. Indeed, the most forward position for all the devices resulted in higher maximum lift and increment to the respective angle of attack. Ramp-up and hold wind tunnel tests also confirmed the best performance of jagged bio-inspired devices nearer the leading edge. The aerodynamic response to the pitching motion showed that the stall and post-stall regime occur much smoother, indicating that the approach presents good potential for dynamic stall or gust response passive control.

     

Influence of heat transfer on the aerodynamic performance of a plunging and pitching NACA0012 airfoil at low Reynolds numbers

The unsteady low Reynolds number aerodynamics phenomena around flapping wings are addressed in several investigations. Elsewhere, airfoils at higher Mach numbers and Reynolds numbers have been treated quite comprehensively in the literature. It is duly noted that the influence of heat transfer phenomena on the aerodynamic performance of flapping wings configurations is not well studied. The objective of the present study is to investigate the effect of heat transfer upon the aerodynamic performance of a pitching and plunging NACA0012 airfoil in the low Reynolds number flow regime with particular emphasis upon the airfoil's lift and drag coefficients. The compressible Navier–Stokes equations are solved using a finite volume method. To consider the variation of fluid properties with temperature, the values of dynamic viscosity and thermal diffusivity are evaluated with Sutherland's formula and the Eucken model, respectively. Instantaneous and mean lift and drag coefficients are calculated for several temperature differences between the airfoil surface and freestream within the range 0–100 K. Simulations are performed for a prescribed airfoil motion schedule and flow parameters. It is learnt that the aerodynamic performance in terms of the lift C_L and drag C_D behavior is strongly dependent upon the heat transfer rate from the airfoil to the flow field. In the plunging case, the mean value of C_D tends to increase, whereas the amplitude of C_L tends to decrease with increasing temperature difference. In the pitching case, on the other hand, the mean value and the amplitude of both C_D and C_L decrease. A spectral analysis of C_D and C_L in the pitching case shows that the amplitudes of both C_D and C_L decrease with increasing surface temperature, whereas the harmonic frequencies are not affected.     

横摆振荡翼型动态气动掠效应试验研究

大型风力机设计对获取翼型更加全面、准确的动态载荷提出更高要求, 研究翼型横摆振荡动态气动特性具有重要意义. 借助"电子凸轮"技术和动态数据同步采集手段, 针对翼型动态“掠效应”首次开展了横摆振荡风洞试验研究, 研究表明: 横摆振荡翼型的气动曲线存在明显迟滞效应, 吸力面压力周期性波动是主要诱因, 且随着振荡频率、初始迎角和振幅的增大, 气动迟滞特性均增强; 升力和压差阻力随横摆角变化的迟滞回线呈"W"形, 俯仰力矩迟滞回线呈"M"形, 升力差量迟滞回线呈"$\infty$"形; 负行程下翼型气动力相对于正行程下的更高, 且负行程下翼型气动力随振荡频率的增大而略有增大, 正行程下则明显减小; 升力系数功率谱密度分布在振荡频率倍频处的能量集中的幅值随着振荡频率增大有增大趋势; 吸力面1.2%和40%弦长处压力的滞回特性较强, 是由于翼面剪切层涡和动态分离涡周期性发展、运动、破裂和重建; 振幅为$10^{\circ}$时, 升力迟滞曲线呈"$^{\wedge}$"形, 振幅为$30^{\circ}$ 时, 升力迟滞曲线呈"$^{\wedge\wedge\wedge}$"形.