skip to main content

Title: Leading-edge-vortex tailoring on unsteady airfoils using an inverse aerodynamic approach
In this paper, we present an approach to obtain a desired leading-edge vortex (LEV) shedding pattern from unsteady airfoils through the execution of suitable motion kinematics. Previous research revealed that LEV shedding is associated with the leading-edge suction parameter (LESP) exceeding a maximum threshold. A low-order method called LESP-modulated discrete vortex method (LDVM) was also developed to predict the onset and termination of LEV shedding from an airfoil undergoing prescribed motion kinematics. In the current work, we present an inverse-aerodynamic formulation based on the LDVM to generate the appropriate motion kinematics to achieve a prescribed LESP variation, and thus, the desired LEV shedding characteristics from the airfoil. The algorithm identifies the kinematic state of the airfoil required to attain the target LESP value through an iterative procedure performed inside the LDVM simulation at each time step. Several case studies are presented to demonstrate design scenarios such as tailoring the duration and intensity of LEV shedding, inducing LEV shedding from the chosen surface of the airfoil, promoting or suppressing LEV shedding during an unsteady motion on demand, and achieving similar LEV shedding patterns using different maneuvers. The kinematic profiles generated by the low-order formulation are also simulated using a high-fidelity unsteady more » Reynolds-averaged Navier–Stokes method to confirm the accuracy of the low-order model. « less
; ; ;
Award ID(s):
Publication Date:
Journal Name:
Physics of Fluids
Page Range or eLocation-ID:
Sponsoring Org:
National Science Foundation
More Like this
  1. The rising global trend to reduce dependence on fossil fuels has provided significant motivation toward the development of alternative energy conversion methods and new technologies to improve their efficiency. Recently, oscillating energy harvesters have shown promise as highly efficient and scalable turbines, which can be implemented in areas where traditional energy extraction and conversion are either unfeasible or cost prohibitive. Although such devices are quickly gaining popularity, there remain a number of hurdles in the understanding of their underlying fluid dynamics phenomena. The ability to achieve high efficiency power output from oscillating airfoil energy harvesters requires exploitation of the complexities of the event of dynamic stall. During dynamic stall, the oncoming flow separates at the leading edge of the airfoil to form leading ledge vortex (LEV) structures. While it is well known that LEVs play a significant role in aerodynamic force generation in unsteady animal flight (e.g. insects and birds), there is still a need to further understand their spatiotemporal evolution in order to design more effective energy harvesting enhancement mechanisms. In this work, we conduct extensive experimental investigations to shed-light on the flow physics of a heaving and pitching airfoil energy harvester operating at reduced frequencies of k =more »fc=U1 = 0.06-0.18, pitching amplitude of 0 = 75 and heaving amplitude of h0 = 0:6c. The experimental work involves the use of two-component particle image velocimetry (PIV) measurements conducted in a wind tunnel facility at Oregon State University. Velocity fields obtained from the PIV measurements are analyzed qualitatively and quantitatively to provide a description of the dynamics of LEVs and other flow structures that may be present during dynamic stall. Due to the difficulties of accurately measuring aerodynamic forces in highly unsteady flows in wind tunnels, a reduced-order model based on the vortex-impulse theory is proposed for estimating the aerodynamic loadings and power output using flow field data. The reduced-order model is shown to be dominated by two terms that have a clear physical interpretation: (i) the time rate of change of the impulse of vortical structures and (ii) the Kutta-Joukowski force which indirectly represents the history effect of vortex shedding in the far wake. Furthermore, the effects of a bio-inspired flow control mechanism based on deforming airfoil surfaces on the flow dynamics and energy harvesting performance are investigated. The results show that the aerodynamic loadings, and hence power output, are highly dependent on the formation, growth rate, trajectory and detachment of the LEV. It is shown that the energy harvesting efficiency increases with increasing reduced frequency, peaking at 25% when k = 0.14, agreeing very well with published numerical results. At this optimal reduced frequency, the time scales of the LEV evolution and airfoil kinematics are matched, resulting in highly correlated aerodynamic load generation and airfoil motion. When operating at k > 0:14, it is shown that the aerodynamic moment and airfoil pitching motion become negatively correlated and as a result, the energy harvesting performance is deteriorated. Furthermore, by using a deforming airfoil surface at the leading and trailing edges, the peak energy harvesting efficiency is shown to increase by approximately 17% and 25% relative to the rigid airfoil, respectively. The performance enhancement is associated with enhanced aerodynamic forces for both the deforming leading and trailing edges. In addition, The deforming trailing edge airfoil is shown to enhance the correlation between the aerodynamic moment and pitching motion at higher reduced frequencies, resulting in a peak efficiency at k = 0:18 as opposed to k = 0:14 for the rigid airfoil.« less
  2. This paper presents a state-variable formulation to model and simulate the 2D unsteady aerodynamics of an airfoil undergoing arbitrary motion kinematics. The model builds upon a large-angle unsteady aerodynamic formulation in which the airfoil is represented using a lumped vortex element (LVE) model. The airfoil is divided into several panels, with a bound vortex placed on each panel. At any time instant, the bound-vortex strengths are determined by employing zero-normal-flow conditions at the control points located on each panel. The vorticity shed from the trailing edge of the airfoil is modeled using discrete vortices that move freely in the flow field. The required state variables are first identified, and all the time derivative terms of the state variables are then derived to form the final state-variable representation. Trailing-edge vortex shedding is incorporated using the Kelvin condition. The final state variable equation can be solved as an ordinary differential equation using any standard ODE-solving algorithm. Three case studies are presented here to evaluate the predictions of the model. In the cases considered here, the airfoil undergoes various unsteady plunge motions. The aerodynamic load history and the wake patterns are compared against the results from the low-order model developed by Narsipur etmore »al. [1] in previous research. The comparison shows that the current state-variable formulation captures the unsteady flow characteristics and the aerodynamic load in good agreement with the reference results.« less
  3. This study explores the feasibility of using the vortex impulse approach, based on experimen- tally generated velocity elds to estimate the energy harvesting performance of a sinusoidally apping foil. Phase-resolved, two-component particle image velocimetry measurements are conducted in a low-speed wind tunnel to capture the ow eld surrounding the apping foil at reduced frequencies of k = fc=U1 = 0.06 - 0.16, pitching amplitude of 0 = 70 and heaving amplitude of h0=c = 0:6. The model results show that for the conditions tested, a maximum energy harvesting eciency of 25% is attained near k = 0:14, agreeing very well with published numerical and experimental results in both accuracy and general trend. The vortex impulse method identi es key contributions to the transient power production from both linear and angular momentum e ects. The eciency reduction at larger values of reduced frequencies is shown to be a result of the reduced power output from the angular momentum. Further, the impulse formulation is decomposed into contributions from posi- tive and negative vorticity in the ow and is used to better understand the uid dynamic mechanisms responsible for producing a peak in energy harvesting performance at k = 0:14. At the largermore »k values, there is a reduction of the advective time scales of the leading edge vortex (LEV) formation. Consequently, the LEV that is shed during the previous half cycle interacts with the foil at the current half cycle resulting in a large negative pitching power due to the reversed direction of the kinematic motion. This vortex capture process signif- icantly decreases the total cycle averaged power output and energy harvesting eciency. These results show the link between the kinematic motion and LEV time scales that a ect the overall power production.« less
  4. MDPI (Ed.)
    Through computational fluid dynamics (CFD) simulations of a model manta ray body, the hydrodynamic role of manta-like bioinspired flapping is investigated. The manta ray model motion is reconstructed from synchronized high-resolution videos of manta ray swimming. Rotation angles of the model skeletal joints are altered to scale the pitching and bending, resulting in eight models with different pectoral fin pitching and bending ratios. Simulations are performed using an in-house developed immersed boundary method-based numerical solver. Pectoral fin pitching ratio (PR) is found to have significant implications in the thrust and efficiency of the manta model. This occurs due to more optimal vortex formation and shedding caused by the lower pitching ratio. Leading edge vortexes (LEVs) formed on the bottom of the fin, a characteristic of the higher PR cases, produced parasitic low pressure that hinders thrust force. Lowering the PR reduces the influence of this vortex while another LEV that forms on the top surface of the fin strengthens it. A moderately high bending ratio (BR) can slightly reduce power consumption. Finally, by combining a moderately high BR = 0.83 with PR = 0.67, further performance improvements can be made. This enhanced understanding of manta-inspired propulsive mechanics fills a gapmore »in our understanding of the manta-like mobuliform locomotion. This motivates a new generation of manta-inspired robots that can mimic the high speed and efficiency of their biological counterpart« less
  5. 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 formore »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.

    « less