skip to main content


Title: Quasi-Steady versus Navier–Stokes Solutions of Flapping Wing Aerodynamics
Various tools have been developed to model the aerodynamics of flapping wings. In particular, quasi-steady models, which are considerably faster and easier to solve than the Navier–Stokes equations, are often utilized in the study of flight dynamics of flapping wing flyers. However, the accuracy of the quasi-steady models has not been properly documented. The objective of this study is to assess the accuracy of a quasi-steady model by comparing the resulting aerodynamic forces against three-dimensional (3D) Navier–Stokes solutions. The same wing motion is prescribed at a fruit fly scale. The pitching amplitude, axis, and duration are varied. Comparison of the aerodynamic force coefficients suggests that the quasi-steady model shows significant discrepancies under extreme pitching motions, i.e., the pitching motion is large, quick, and occurs about the leading or trailing edge. The differences are as large as 1.7 in the cycle-averaged lift coefficient. The quasi-steady model performs well when the kinematics are mild, i.e., the pitching motion is small, long, and occurs near the mid-chord with a small difference in the lift coefficient of 0.01. Our analysis suggests that the main source for the error is the inaccuracy of the rotational lift term and the inability to model the wing-wake interaction in the quasi-steady model.  more » « less
Award ID(s):
1761618
NSF-PAR ID:
10095723
Author(s) / Creator(s):
; ; ; ;
Date Published:
Journal Name:
Fluids
Volume:
3
Issue:
4
ISSN:
2311-5521
Page Range / eLocation ID:
81
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. Abstract

    A limiting factor in the design of smaller size uncrewed aerial vehicles is their inability to navigate through gust-laden environments. As a result, engineers have turned towards bio-inspired engineering approaches for gust mitigation techniques. In this study, the aerodynamics of a red-tailed hawk’s response to variable-magnitude discrete transverse gusts was investigated. The hawk was flown in an indoor flight arena instrumented by multiple high-speed cameras to quantify the 3D motion of the bird as it navigated through the gust. The hawk maintained its flapping motion across the gust in all runs; however, it encountered the gust at different points in the flapping pattern depending on the run and gust magnitude. The hawk responded with a downwards pitching motion of the wing, decreasing the wing pitch angle to between −20and −5, and remained in this configuration until gust exit. The wing pitch data was then applied to a lower-order aerodynamic model that estimated lift coefficients across the wing. In gusts slower than the forward flight velocity (low gust ratio), the lift coefficient increases at a low-rate, to a maximum of around 2–2.5. In gusts faster than the forward flight velocity (high gust ratio), the lift coefficient initially increased rapidly, before increasing at a low-rate to a value around 4–5. In both regimes, the hawk’s observed height change due to gust interaction was similar (and small), despite larger estimated lift coefficients over the high gust regime. This suggests another mitigation factor apart from the wing response is present. One potential factor is the tail pitching response observed here, which prior work has shown serves to mitigate pitch disturbances from gusts.

     
    more » « less
  2. Fluid–structure interaction (FSI) plays a significant role in the deformation of flapping insect wings. However, many current FSI models are high-order and rely on direct computational methods, thereby limiting parametric studies as well as insights into the physics governing wing dynamics. We develop a novel flapping wing FSI framework that accommodates general wing geometry and fluid loading. We use this framework to study the unilaterally coupled FSI of an idealized hawkmoth forewing considering two fluid models: Reynolds-averaged Navier–Stokes computational fluid dynamics (RANS CFD) and blade element theory (BET). We first compare aerodynamic modal forces estimated by the low-order BET model to those calculated via high fidelity RANS CFD. We find that for realistic flapping kinematics, BET estimates modal forces five orders of magnitude faster than CFD within reasonable accuracy. Over the range flapping kinematics considered, BET and CFD estimated modal forces vary maximally by 350% in magnitude and approximately π/2 radians in phase. The large reduction in computational time offered by BET facilitates high-dimensional parametric design of flapping-wing-based technologies. Next, we compare the contributions of aerodynamic and inertial forces to wing deformation. Under the unilateral coupling assumption, aerodynamic and inertial-elastic forces are on the same order of magnitude—however, inertial-elastic forces primarily excite the wing’s bending mode whereas aerodynamic forces primarily excite the wing’s torsional mode. This suggests that, via conscientious sensor placement and orientation, biological wings may be able to sense independently inertial and aerodynamic forces. 
    more » « less
  3. Abstract

    An experimental study was undertaken to evaluate the power extraction of an airfoil undergoing large amplitude pitching and heaving using a trailing edge flapping motion for the application of energy harvesting for steady flow over the airfoil. The airfoil was a NACA0015 design, pitching at the 1/3 chord position, with an actively controlled trailing edge flap hinged at the 2/3 chord location (chord length of c = 150mm and aspect ratio AR = 2, however end plates were used to simulate a two-dimensional airfoil). Data were obtained over a range of wind speeds corresponding to Reynolds numbers in the 30,000–60,000 range in a low-speed wind tunnel with turbulence intensities below 2%. The results are characterized using the reduced frequency, k = fc/U∞ over the range of 0.04–0.08, where f is the oscillating frequency in Hz, and U∞ is the freestream velocity. The pitching and heaving amplitudes are θ0 = 70° and h0 = 0.6c respectively, with a phase delay of 90°. Two trailing edge motion profiles are presented, examining the relative phase of trailing edge flap to the pitching phase. For each motion, a positive and negative case are considered. This is a total of 4 trailing edge motion profiles. Trailing edge motion amplitudes of 20° and 40° are compared and results contrasted. Direct transient force measurements were used to obtain the cycle variation of induced aerodynamic loads (lift coefficient) as well as the power output and efficiency. Results are used to identify the influence of trailing edge flap oscillations on the overall performance for energy harvesting, with a maximum efficiency increase of 21.3% and corresponding cycle averaged heaving power coefficient increase of 29.9% observed as a result of trailing edge motion.

     
    more » « less
  4. An analytical model is developed for the lift force produced by unsteady rotating wings; this configuration is a simple representation of a flapping wing. Modelling this is important for the aerodynamic and control-system design for bio-inspired drones. Such efforts have often been limited to being two-dimensional, semi-empirical, sometimes computationally expensive, or quasi-steady. The current model is unsteady and three-dimensional, yet simple to implement, requiring knowledge of only the wing kinematics and geometry. Rotating wings produce a vortex loop consisting of the root vortex, leading-edge vortex, tip vortex and trailing-edge vortex, which grows with time. This is modelled as a tilted planar loop, geometrically specified by the wing size, orientation and motion. By equating the angular impulse of the vortex loop to that of the fluid volume driven by the wing, the circulatory lift force is derived. Potential flow theory gives the fluid-inertial lift. Adding these two contributions yields the total lift formula. The model shows good agreement with a range of experimental and computational cases. Also, a steady-state lift model is developed that compares well with previous work for various angles of attack. 
    more » « less
  5. null (Ed.)
    Abstract

    Flapping insect wings experience appreciable deformation due to aerodynamic and inertial forces. This deformation is believed to benefit the insect’s aerodynamic force production as well as energetic efficiency. However, the fluid-structure interaction (FSI) models used to estimate wing deformations are often computationally demanding and are therefore challenged by parametric studies. Here, we develop a simple FSI model of a flapping wing idealized as a two-dimensional pitching-plunging airfoil. Using the Lagrangian formulation, we derive the reduced-order structural framework governing wing’s elastic deformation. We consider two fluid models: quasi-steady Deformable Blade Element Theory (DBET) and Unsteady Vortex Lattice Method (UVLM). DBET is computationally economical but does not provide insight into the flow structure surrounding the wing, whereas UVLM approximates flows but requires more time to solve. For simple flapping kinematics, DBET and UVLM produce similar estimates of the aerodynamic force normal to the surface of a rigid wing. More importantly, when the wing is permitted to deform, DBET and UVLM agree well in predicting wingtip deflection and aerodynamic normal force. The most notable difference between the model predictions is a roughly 20° phase difference in normal force. DBET estimates wing deformation and force production approximately 15 times faster than UVLM for the parameters considered, and both models solve in under a minute when considering 15 flapping periods. Moving forward, we will benchmark both low-order models with respect to high fidelity computational fluid dynamics coupled to finite element analysis, and assess the agreement between DBET and UVLM over a broader range of flapping kinematics.

     
    more » « less