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.     

Turbulent separated flow over and downstream of a two-element airfoil

Flow characteristics in the vicinity of the flap of a single-slotted airfoil are presented and analysed. The flow remained attached over the model surfaces except in the vicinity of the flap trailing edge where a small region of boundary-layer separation extended over the aft 7% of flap chord. The airfoil configuration was tested at a Mach number of 0.09 and a chord Reynolds number of 1.8 × 10⁶ in the NASA Ames Research Center 7- by 10-Foot Wind Tunnel. The flow was complicated by the presence of a strong, initially inviscid, jet, emanating from the slot between airfoil and flap, and a gradual merging of the main airfoil wake and flap suction-side boundary layer.Research Engineer, NRC Research AssociateAerospace Engineer     

Time resolved PIV analysis of flow over a NACA 0015 airfoil with Gurney flap

A NACA 0015 airfoil with and without a Gurney flap was studied in a wind tunnel with Re _c = 2.0 × 10⁵ in order to examine the evolving flow structure of the wake through time-resolved PIV and to correlate this structure with time-averaged measurements of the lift coefficient. The Gurney flap, a tab of small length (1–4% of the airfoil chord) that protrudes perpendicular to the chord at the trailing edge, yields a significant and relatively constant lift increment through the linear range of the C _L versus α curve. Two distinct vortex shedding modes were found to exist and interact in the wake downstream of flapped airfoils. The dominant mode resembles a Kàrmàn vortex street shedding behind an asymmetric bluff body. The second mode, which was caused by the intermittent shedding of fluid recirculating in the cavity upstream of the flap, becomes more coherent with increasing angle of attack. For a 4% Gurney flap at α = 8°, the first and second modes corresponded with Strouhal numbers based on flap height of 0.18 and 0.13. Comparison of flow around ‘filled’ and ‘open’ flap configurations suggested that the second shedding mode was responsible for a significant portion of the overall lift increment.     

EXPERIMENTAL INVESTIGATION ON LIFT INCREMENT OF A PLASMA CIRCULATION CONTROL AIRFOIL

利用等离子体激励器发展了新型的环量增升技术,并对二维NACA0012翼型绕流实施控制。由于NACA0012翼型为尖后缘构型,环量增升装置由2个非对称型介质阻挡放电等离子体激励器构成。一个等离子体激励器贴附于翼型吸力面靠近后缘处,其诱导的壁面射流沿来流方向指向下游;另一个等离子体激励器贴附于翼型压力面靠近后缘处,其诱导的壁面射流与来流方向相反指向上游。在风洞中通过时间解析二维PIV系统对翼型绕流流场进行了测量,基于翼型弦长的雷诺数Re=20 000。结果表明在等离子体激励器的控制下,翼型压力面靠近后缘处可以形成一个定常回流区,从而起到虚拟气动外形的作用,因此翼型吸力面的流场得到加速,压力面的流场得到减速,使得翼型压力面的吸力以及压力面的压力都得到增加,进而增加了翼型的环量。风洞天平测力实验进一步验证了该环量增升技术的有效性。在整个攻角范围内,施加控制的翼型的升力系数相比没有控制的工况有明显的提高。     

Unsteady flow structure and vorticity convection over the airfoil oscillating at high reduced frequency

Unsteady vortex structures and vorticity convection over the airfoil (NACA 0012), oscillating in the uniform inflow, are studied by flow visualization and velocity measurements. The airfoil, pivoting at one-third of the chord, oscillates periodically near the static stalling angle of attack (AOA) at high reduced-frequency. The phase-triggering and modified phase-averaged techniques are employed to reconstruct the pseudo instantaneous velocity field over the airfoil. During the down stroke cycle, the leading-edge separation vortex is growing and the vortex near the trailing edge begins to shed into the wake. During the upstroke cycle, the leading-edge separation vortex is matured and moves downstream, and the counter clockwise vortex is forming near the trailing edge. Convection speeds and wavelength of the unsteady vortex structure over the airfoil equal to that of the counter clockwise vortex shed into the wake. This kind of vortex structure is termed as “synchronized shedding” type. The wavelength of unsteady vortex structure over the airfoil is significantly different from that at low reduced-frequency. Consistent convection speeds of the leading-edge separation vortex are acquired from the spatial-temporal variations of local circulation and local surface vorticity generation, and equals that predicted from flow visualization. Spatial-temporal variations of the local surface vorticity generation clearly reveal the formation and passage of the leading-edge separation vortex only in the region where the flow does not separate completely from the surface. Significant amounts of the surface vorticity are generated within the leading-edge region of the airfoil during the upstroke cycle. Only negligible amount of surface vorticity is produced within the region of complete flow separation. During the down stroke cycle, the surface vorticity generation is mild along the airfoil surface, except the leading-edge region where a small scale leading-edge separation vortex is forming and growing.     

Separation behaviour in front of a two-dimensional fence

The wall-bounded turbulent shear flow in front of a two-dimensional fence was investigated experimentally. In this prototype of a rapidly separating flow the assumptions for a first-order boundary-layer theory cease to apply. This is caused by the streamline curvature, the ensuing pressure gradient normal to the wall and the large vertical velocity component v in front of the fence. For the present experiment, where the ratio of the boundary layer thickness δ₀ measured without the fence and the fence height h is 0.82, the time mean separation length upstream of the fence l_f is 0.65h. However, instantaneous reverse flow events can be detected up to 4 mean separation lengths l_f upstream of the fence. The maximum value of the reverse flow factor χ_w is 95% indicating a strong reverse flow region. The experiments were performed by LDA and a wall pulsed-wire skin-friction meter. They show the limits of first-order boundary-layer theory and provide the first comprehensive data set of mean and fluctuating velocities and of wall shear-stress for this type of separating flow.     

An experimental study of a turbulent wing-body junction and wake flow

Extensive measurements were conducted in an incompressible turbulent flow around the wing-body junction formed by a 3∶2 semi-elliptic nose/NACA 0020 tail section and a flat plate. Mean and fluctuating velocity measurements were performed adjacent to the wing and up to 11.56 chord lengths downstream. The appendage far wake region was subjected to an adverse pressure gradient. The authors' results show that the characteristic horseshoe vortex flow structure is elliptically shaped, with ? (W)/?Y forming the primary component of the streamwise vorticity. The streamwise development of the flow distortions and vorticity distributions is highly dependent on the geometry-induced pressure gradients and resulting flow skewing directions. The primary goal of this research was to determine the effects of the approach boundary layer characteristics on the junction flow. To accomplish this goal, the authors' results were compared to several other junction flow data sets obtained using the same body shape. The trailing vortex leg flow structure was found to scale on T. A parameter known as the momentum deficit factor (MDF = (Re _T)² (θ/T)) was found to correlate the observed trends in mean flow distortion magnitudes and vorticity distribution. Changes in δ/T were seen to affect the distribution of u′, with lower ratios producing well defined local turbulence maxima. Increased thinning of the boundary layer near the appendage was also observed for small values of δ/T.     

On the generation of yawing moment using active flow control

The generation of control moments without moving control surfaces is of great practical importance. Following a successful flight demonstration of creating roll motion without ailerons using differential, lift oriented, flow control the current study is a first step towards generating yawing motion via differential flow controlled drag.A wind tunnel study was conducted on a 21% thick Glauert type airfoil. The upper surface flow is partially separated from the two-thirds chord location and downstream on this airfoil at all incidence angles. An array of mass-less Piezo-fluidic actuators, located at x/c = 0.65, are capable of fully reattaching the flow in a gradual, controlled manner. The actuators are individually operated such that the boundary layer could be controlled in a 3D fashion.Several concepts for creating yaw motion without moving control surface are examined. The ultimate goal is to generate the same lift on both wings, while decreasing the drag on one wing and increasing the drag on the other, therefore creating a yawing moment. Decreased drag is created by effective part-span separation delay while increased drag can be created by enhanced generation of vortex shedding or by highly localized 3D actuation.Detailed measurements of 3D surface pressure distributions and wake data with three velocity and streamwise vorticity components are presented and discussed along with surface flow visualization images. The data provide evidence that yawing moments can be generated with AFC.     

Stochastic estimation of flow near the trailing edge of a NACA0012 airfoil

A stochastic estimation technique has been applied to simultaneously acquired data of velocity and surface pressure as a tool to identify the sources of wall-pressure fluctuations. The measurements have been done on a NACA0012 airfoil at a Reynolds number of Re _c  = 2 × 10⁵, based on the chord of the airfoil, where a separated laminar boundary layer was present. By performing simultaneous measurements of the surface pressure fluctuations and of the velocity field in the boundary layer and wake of the airfoil, the wall-pressure sources near the trailing edge (TE) have been studied. The mechanisms and flow structures associated with the generation of the surface pressure have been investigated. The “quasi-instantaneous” velocity field resulting from the application of the technique has led to a picture of the evolution in time of the convecting surface pressure generating flow structures and revealed information about the sources of the wall-pressure fluctuations, their nature and variability. These sources are closely related to those of the radiated noise from the TE of an airfoil and to the vibration issues encountered in ship hulls for example. The NACA0012 airfoil had a 30 cm chord and aspect ratio of 1.     

Experimental investigation of a blunt trailing edge flow field with application to sound generation

The unsteady lift generated by turbulence at the trailing edge of an airfoil is a source of radiated sound. The objective of the present research was to measure the velocity field in the near wake region of an asymmetric beveled trailing edge in order to determine the flow mechanisms responsible for the generation of trailing edge noise. Two component velocity measurements were acquired using particle image velocimetry. The chord Reynolds number was 1.9 × 10⁶. The data show velocity field realizations that were typical of a wake flow containing an asymmetric periodic vortex shedding. A phase average decomposition of the velocity field with respect to this shedding process was utilized to separate the large scale turbulent motions that occurred at the vortex shedding frequency (i.e., those responsible for the production of tonal noise) from the smaller scale turbulent motions, which were interpreted to be responsible for the production of broadband sound. The small scale turbulence was found to be dependent on the phase of the vortex shedding process implying a dependence of the broadband sound generated by the trailing edge on the phase of the vortex shedding process.     

Influences of wind tunnel parameters on airfoil characteristics at high subsonic speeds

The influences of several wind tunnel parameters on airfoil characteristics are experimentally investigated in a transonic wind tunnel. Quantified as Mach number errors they show decisive effects and confirm that they have to be taken into consideration in a valuation of test results. Differences in data gained in measurements of several wind tunnels can be partly explained thereby.List of symbols C chord length - C _d drag coefficient - C ₀ d drag coefficient at zero lift - C ₁ lift coefficient - C _{1 max} lift coefficient at maximum lift - C _p pressure coefficient - D divergence - M Mach number - P ₀ stagnation pressure - Re Reynolds number - S _h suction along the horizontal walls - S _v suction along the vertical walls - X axial coordinate - X _ref position of the reference point - X _s shock position - angle of attack     

Pressure field due to drag reducing outer layer devices in turbulent boundary layers

The wall static pressure in the vicinity of drag reducing outer layer devices in flat wall turbulent boundary layers has been measured and compared with an inviscid theory. Symmetric and cambered airfoil devices have been examined at small angles of attack and very low chord Reynolds numbers. Airfoil devices impose a sequence of strong favorable and adverse pressure gradients on the boundary layer whose drag is to be reduced. At very small angles of attack (± 2°), this pressure field extends up to about three chord lengths downstream of the trailing edge of an airfoil device. Also examined are the pressures on the upper and lower surfaces of a symmetric airfoil device in the freestream and near the wall. The freestream pressure distribution around an airfoil section is altered by the wall proximity. The relevance of lift enhancement caused by wall proximity to drag reduction has been discussed. The pressure distributions on the flat wall beneath the symmetric airfoil devices are predicted well by the inviscid theory. However, the remaining pressure distributions are predicted only qualitatively, presumably because of strong viscous effects.     

Symmetry groups and the pseudo-Riemann spacetimes for mixed-hardening elastoplasticity

The constitutive postulations for mixed-hardening elastoplasticity are selected. Several homeomorphisms of irreversibility parameters are derived, among which X_a⁰ and X_c⁰ play respectively the roles of temporal components of the Minkowski and conformal spacetimes. An augmented vector X_a:=(YQ_a^t,YQ_a⁰)^t is constructed, whose governing equations in the plastic phase are found to be a linear system with a suitable rescaling proper time. The underlying structure of mixed-hardening elastoplasticity is a Minkowski spacetime $\text{M}{}^{\mathrm{n+1}}$ on which the proper orthochronous Lorentz group SO_o(n,1) left acts. Then, constructed is a Poincaré group ISO_o(n,1) on space X:=X_a+X_b, of which X_b reflects the kinematic hardening rule in the model. We also find that the space (Q_a^t,q₀^a) is a Robertson–Walker spacetime, which is conformal to X_a through a factor Y, and conformal to X_c:=(ρQ_a^t,ρQ_a⁰)^t through a factor ρ as given by ρ(q₀^a)=Y(q₀^a)/[1−2ρ₀Q_a⁰(0)+2ρ₀Y(q₀^a)Q_a⁰(q₀^a)]. In the conformal spacetime the internal symmetry is a conformal group.     

A study on the mechanism of high-lift generation by an airfoil in unsteady motion at low reynolds number

The aerodynamic force and flow structure of NACA 0012 airfoil performing an unsteady motion at low Reynolds number (Re=100) are calculated by solving Navier-Stokes equations. The motion consists of three parts: the first translation, rotation and the second translation in the direction opposite to the first. The rotation and the second translation in this motion are expected to represent the rotation and translation of the wing-section of a hovering insect. The flow structure is used in combination with the theory of vorticity dynamics to explain the generation of unsteady aerodynamic force in the motion. During the rotation, due to the creation of strong vortices in short time, large aerodynamic force is produced and the force is almost normal to the airfoil chord. During the second translation, large lift coefficient can be maintained for certain time period and , the lift coefficient averaged over four chord lengths of travel, is larger than 2 (the corresponding steady-state lift coefficient is only 0.9). The large lift coefficient is due to two effects. The first is the delayed shedding of the stall vortex. The second is that the vortices created during the airfoil rotation and in the near wake left by previous translation form a short “vortex street” in front of the airfoil and the “vortex street” induces a “wind”; against this “wind” the airfoil translates, increasing its relative speed. The above results provide insights to the understanding of the mechanism of high-lift generation by a hovering insect. The project supported by the National Natural Science Foundation of China (19725210)     

Flows around two airfoils performing fling and subsequent translation and translation and subsequent clap

The aerodynamic forces and flow structures of two airfoils performing “fling and subsequent translation“ and “translation and subsequent clap“ are studied by numerically solving the Navier-Stokes equations in moving overset grids. These motions are relevant to the flight of very small insects. The Reynolds number, based on the airfoil chord length c and the translation velocity U, is 17. It is shown that: (1) For two airfoils performing fling and subsequent translation, a large lift is generated both in the fling phase and in the early part of the translation phase. During the fling phase,a pair of leading edge vortices of large strength is generated; the generation of the vortex pair in a short period results in a large time rate of change of fluid impulse, which explains the large lift in this period. During the early part of the translation, the two leading edge vortices move with the airfoils;the relative movement of the vortices also results in a large time rate of change of fluid impulse, which explains the large lift in this part of motion. (In the later part of the translation, the vorticity in the vortices is diffused and convected into the wake.) The time averaged lift coefficient is approximately 2.4 times as large as that of a single airfoil performing a similar motion. (2) For two airfoils performing translation and subsequent clap, a large lift is generated in the clap phase. During the clap, a pair of trailing edge vortices of large strength are generated; again, the generation of the vortex pair in a short period (which results in a large time rate of change of fluid impulse) is responsible for the large lift in this period. The time averaged lift coefficient is approximately 1.6 times as large as that of a single airfoil performing a similar motion. (3) When the initial distance between the airfoils (in the case of clap, the final distance between the airfoils) varies from 0.1 to 0.2c, the lift on an airfoil decreases only slightly but the torque decreases greatly. When the distance is about lc， the interference effects between the two airfoils become very small.     

Effect of velocity ratio on the streamwise vortex structures in the wake of a stack

The time-averaged velocity and streamwise vorticity fields within the wake of a stack were investigated in a low-speed wind tunnel using a seven-hole pressure probe. The experiments were conducted at a Reynolds number, based on the stack external diameter, of Re_D=2.3×10⁴. The stack, of aspect ratio AR=9, was mounted normal to a ground plane and was partially immersed in a flat-plate turbulent boundary layer, where the ratio of the boundary layer thickness to the stack height was δ/H≈0.5. The jet-to-cross-flow velocity ratio was varied from R=0 to 3, which covered the downwash, crosswind-dominated and jet-dominated flow regimes. In the downwash and crosswind-dominated flow regimes, two pairs of counter-rotating streamwise vortex structures were identified within the stack wake. The tip vortex pair located close to the free end of the stack, and the base vortex pair located close to the ground plane within the flat-plate boundary layer, were similar to those found in the wake of a finite circular cylinder, and were associated with the upwash and downwash flow fields within the stack wake, respectively. In the jet-dominated flow regime, a third pair of streamwise vortex structures was observed, referred to as the jet-wake vortex pair, which occurred within the jet-wake region above the free end of the stack. The jet-wake vortex pair had the same orientation as the base vortex pair and was associated with the jet rise. The peak vorticity and strength of the streamwise vortex structures were functions of the jet-to-cross-flow velocity ratio. For the tip vortex structures, their peak vorticity and strength reduced as the jet-to-cross-flow velocity ratio increased.     

Experimental study of pitching and plunging airfoils at low Reynolds numbers

Measurements of the unsteady flow structure and force time history of pitching and plunging SD7003 and flat plate airfoils at low Reynolds numbers are presented. The airfoils were pitched and plunged in the effective angle of attack range of 2.4°–13.6° (shallow-stall kinematics) and ?6° to 22° (deep-stall kinematics). The shallow-stall kinematics results for the SD7003 airfoil show attached flow and laminar-to-turbulent transition at low effective angle of attack during the down stroke motion, while the flat plate model exhibits leading edge separation. Strong Re-number effects were found for the SD7003 airfoil which produced approximately 25 % increase in the peak lift coefficient at Re = 10,000 compared to higher Re flows. The flat plate airfoil showed reduced Re effects due to leading edge separation at the sharper leading edge, and the measured peak lift coefficient was higher than that predicted by unsteady potential flow theory. The deep-stall kinematics resulted in leading edge separation that led to formation of a large leading edge vortex (LEV) and a small trailing edge vortex (TEV) for both airfoils. The measured peak lift coefficient was significantly higher (~50 %) than that for the shallow-stall kinematics. The effect of airfoil shape on lift force was greater than the Re effect. Turbulence statistics were measured as a function of phase using ensemble averages. The results show anisotropic turbulence for the LEV and isotropic turbulence for the TEV. Comparison of unsteady potential flow theory with the experimental data showed better agreement by using the quasi-steady approximation, or setting C(k) = 1 in Theodorsen theory, for leading edge–separated flows.     

A Model Problem of the Design and Aerodynamical Calculation of the Gliding Wing Airfoil of a Wing-in-Ground Vehicle

The problem of the design and aerodynamical calculation of a wing airfoil whose trailing edge glides over a plane horizontal surface is formulated and solved. The known airfoil undersurface is a rectilinear segment forming a given angle with the ground surface, while the upper surface is sought on the basis of a preassigned velocity distribution. This distribution is taken from a class of hydrodynamically advantageous distributions ensuring separationless flow around the airfoil within the framework of the mathematical flow model adopted. The problem is reduced to a mixed boundary value problem in a half-plane and solved analytically. For calculating the lift coefficient the assumption of a thin jetlet flowing between the horizontal region of the airfoil contour and the ground surface is introduced. The effect of the law of pressure decrease along the jetlet, from the stagnation to the exit value, on the lift coefficient is studied. On the basis of the calculations performed a conclusion is drawn concerning the influence of the rectilinear region slope on the airfoil shape; it is also shown how the slope and the value of the maximum velocity on the airfoil affect its shape and the lift coefficient.     

An experimental study of a bio-inspired corrugated airfoil for micro air vehicle applications

An experimental study was conducted to investigate the aerodynamic characteristics of a bio-inspired corrugated airfoil compared with a smooth-surfaced airfoil and a flat plate at the chord Reynolds number of Re _C  = 58,000–125,000 to explore the potential applications of such bio-inspired corrugated airfoils for micro air vehicle designs. In addition to measuring the aerodynamic lift and drag forces acting on the tested airfoils, a digital particle image velocimetry system was used to conduct detailed flowfield measurements to quantify the transient behavior of vortex and turbulent flow structures around the airfoils. The measurement result revealed clearly that the corrugated airfoil has better performance over the smooth-surfaced airfoil and the flat plate in providing higher lift and preventing large-scale flow separation and airfoil stall at low Reynolds numbers (Re _C  < 100,000). While aerodynamic performance of the smooth-surfaced airfoil and the flat plate would vary considerably with the changing of the chord Reynolds numbers, the aerodynamic performance of the corrugated airfoil was found to be almost insensitive to the Reynolds numbers. The detailed flow field measurements were correlated with the aerodynamic force measurement data to elucidate underlying physics to improve our understanding about how and why the corrugation feature found in dragonfly wings holds aerodynamic advantages for low Reynolds number flight applications.     

Low Reynolds number airfoil aerodynamic loads determination via line integral of velocity obtained with particle image velocimetry

The small magnitude lift forces generated by both a NACA 0012 airfoil and a thin flat plate at Re?=?29,000 and 54,000 were determined through the line integral of velocity, obtained with particle image velocimetry, via the application of the Kutta–Joukowsky theorem. Surface pressure measurements of the NACA0012 airfoil were also obtained to validate the lift coefficient C _l. The bound circulation was found to be insensitive to the size and aspect ratio of the rectangular integration loop for pre-stall angles. The present C _l data were also found to agree very well with the surface pressure-determined lift coefficient for pre-stall conditions. A large variation in C _l with the loop size and aspect ratio for post-stall conditions was, however, observed. Nevertheless, the present flat-plate C _l data were also found to collectively agree with the published force-balance measurements at small angles of attack, despite the large disparity exhibited among the various published data at high angles. Finally, the ensemble-averaged wake velocity profiles were also used to compute the drag coefficient and, subsequently, the lift-to-drag ratio.