skip to main content


Title: Extension of the Unsteady Vortex Lattice Method: viscous Oseen vortices and thickness effects
In this study, the standard unsteady vortex lattice method (UVLM) is modified by different approaches to accommodate the viscous effects on the flow field and lift force history of a moving airfoil. Firstly, the conventional inviscid vortices are replaced with ones that satisfy Oseen’s approximation of the Navier–Stokes equations. These are known as viscous vortices and are a function of the Reynolds number. Compared to the standard UVLM, the results show a different behavior for the wake deformation at lower Reynolds numbers where the viscous effects are more dominant. The second modification is implemented by an elegant alteration in the conservation of circulation condition based on the fact that the whole domain composed of the fluid and solid can be treated as a single kinematic entity. The no-slip boundary condition enables considering the moving solid vorticity in the Kelvin condition, where the results for the frequency response indicate more phase lag compared to the standard UVLM or potential flow theories (i.e., Theodorsen), which increases at higher reduced frequencies. Furthermore, the effect of the position of the shed vortices is analyzed and seen to affect the amplitude of the resulting lift history. Finally, the no-slip boundary condition is also explicitly considered by implementing source singularities along the surface of a thick airfoil, resulting in a better prediction of the lift magnitude compared to CFD results. The combination of all the aforementioned modifications has the potential to enhance the performance of the UVLM in case of viscous flows at relatively high reduced frequencies.  more » « less
Award ID(s):
1635673
NSF-PAR ID:
10130335
Author(s) / Creator(s):
; ;
Date Published:
Journal Name:
AIAA SciTech
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. null (Ed.)
    In this paper, we started by summarizing our recently developed viscous unsteady theory based on coupling potential flow with the triple deck boundary layer theory. This approach provides a viscous extension of potential flow unsteady aerodynamics. As such, a Reynolds- number-dependence could be determined. We then developed a finite-state approximation of such a theory, presenting it in a state space model. This novel nonlinear state space model of the viscous unsteady aerodynamic loads is expected to serve aerodynamicists better than the classical Theodorsen’s model, as it captures viscous effects (i.e., Reynolds number dependence) as well as nonlinearity and additional lag in the lift dynamics; and allows simulation of arbitrary time-varying airfoil motions (not necessarily harmonic). Moreover, being in a state space form makes it quite convenient for simulation and coupling with structural dynamics to perform aeroelasticity, flight dynamics analysis, and control design. We then proceeded to develop a linearization of such a model, which enables analytical results. So, we derived an analytical representation of the viscous lift frequency response function, which is an explicit function of, not only frequency, but also Reynolds number. We also developed a state space model of the linearized response. We finally simulated the nonlinear and linear models to a non- harmonic, small-amplitude pitching maneuver at 100 , 000 Reynolds number and compared the resulting lift and pitching moment with potential flow, in reference to relatively higher fidelity computations of the Unsteady Reynolds-Averaged Navier-Stokes equations. 
    more » « less
  2. Large-eddy simulations (LES) of the fluid flow over a NACA0018 airfoil at AOA =20 degrees angle of attack are performed to investigate the effect of surface morphing oscillations on the aerodynamic performance of the airfoil over a wide range of Reynolds numbers (Re = 5,000 to 500,000). These oscillations are in the form of low amplitude backward (opposite to the airfoil's forward motion) traveling wave actuations on the upper surface of the airfoil. The sharp interface curvilinear immersed boundary (CURVIB) method is used to handle the moving surface of the airfoil. The nondimensional amplitude is a*=0.001 (a*=a/L; a: amplitude, L: chord length of the airfoil) and reduced frequency (f*= fL/U; f is the frequency and U is the freestream velocity) is chosen to match the leading edge vortex shedding frequency. The results of the simulations at the post-stall angle of attack (AOA =20 degrees) show that the lift coefficient increases more than 20% and the drag coefficient decreases more than 40% within the Reynolds number range of Re = 50,000-500,000 for traveling wave actuation of amplitude, a*=0.001, and frequency, f*=8. However, the lift and drag coefficients of the actuated airfoil were similar to the baseline airfoil for Re = 5,000. 
    more » « less
  3. Abstract

    Bio-inspired flow control strategies can provide a new paradigm of efficiency and adaptability to overcome the operational limitations of traditional flow control. This is particularly useful to small-scale uncrewed aerial vehicles since their mission requirements are rapidly expanding, but they are still limited in terms of agility and adaptability when compared to their biological counterparts, birds. One of the flow control strategies that birds implement is the deployment of covert feathers. In this study, we investigate the performance characteristics and flow physics of torsionally hinged covert-inspired flaps mounted on the suction side of a NACA2414 airfoil across different Reynolds numbers, specifically 200,000 and 1,000. These two Reynolds numbers are representative of different avian flight regimes where covert feathers have been observed to deploy during flight, namely cruising and landing/perching. We performed experiments and simulations where we varied the flap location, the hinge stiffness, and the moment of inertia of the flap to investigate the aerodynamic performance and describe the effects of the structural parameters of the flap on the aerodynamic lift improvements. Results of the study show up to 12% lift improvement post-stall for the flapped cases when compared to the flap-less baseline. The post-stall lift improvement is sensitive to the flap’s structural properties and location. For instance, the hinge stiffness controls the mean deflection angle of the flap, which governs the resulting time-averaged lift improvements. The flap moment of inertia, on the other hand, controls the flap dynamics, which in turn controls the flap’s lift-enhancing mechanism and how the flap affects the instantaneous lift. By examining the time-averaged and instantaneous lift measurement, we uncover the mechanisms by which the covert-inspired flap improves lift and highlights similarities and differences across Reynolds numbers. This article highlights the feasibility of using covert-inspired flaps as flow control across different flight missions and speeds.

     
    more » « less
  4. This study presents the first 3D two-way coupled fluid structure interaction (FSI) simulation of a hybrid anechoic wind tunnel (HAWT) test section with modeling all important effects, such as turbulence, Kevlar wall porosity and deflection, and reveals for the first time the complete 3D flow structure associated with a lifting model placed into a HAWT. The Kevlar deflections are captured using finite element analysis (FEA) with shell elements operated under a membrane condition. Three-dimensional RANS CFD simulations are used to resolve the flow field. Aerodynamic experimental results are available and are compared against the FSI results. Quantitatively, the pressure coefficients on the airfoil are in good agreement with experimental results. The lift coefficient was slightly underpredicted while the drag was overpredicted by the CFD simulations. The flow structure downstream of the airfoil showed good agreement with the experiments, particularly over the wind tunnel walls where the Kevlar windows interact with the flow field. A discrepancy between previous experimental observations and juncture flow-induced vortices at the ends of the airfoil is found to stem from the limited ability of turbulence models. The qualitative behavior of the flow, including airfoil pressures and cross-sectional flow structure is well captured in the CFD. From the structural side, the behavior of the Kevlar windows and the flow developing over them is closely related to the aerodynamic pressure field induced by the airfoil. The Kevlar displacement and the transpiration velocity across the material is dominated by flow blockage effects, generated aerodynamic lift, and the wake of the airfoil. The airfoil wake increases the Kevlar window displacement, which was previously not resolved by two-dimensional panel-method simulations. The static pressure distribution over the Kevlar windows is symmetrical about the tunnel mid-height, confirming a dominantly two-dimensional flow field. 
    more » « less
  5. null (Ed.)
    A closed-loop control algorithm for the reduction of turbulent flow separation over NACA 0015 airfoil equipped with leading-edge synthetic jet actuators (SJAs) is presented. A system identification approach based on Nonlinear Auto-Regressive Moving Average with eXogenous inputs (NARMAX) technique was used to predict nonlinear dynamics of the fluid flow and for the design of the controller system. Numerical simulations based on URANS equations are performed at Reynolds number of 106 for various airfoil incidences with and without closed-loop control. The NARMAX model for flow over an airfoil is based on the static pressure data, and the synthetic jet actuator is developed using an incompressible flow model. The corresponding NARMAX identification model developed for the pressure data is nonlinear; therefore, the describing function technique is used to linearize the system within its frequency range. Low-pass filtering is used to obtain quasi-linear state values, which assist in the application of linear control techniques. The reference signal signifies the condition of a fully re-attached flow, and it is determined based on the linearization of the original signal during open-loop control. The controller design follows the standard proportional-integral (PI) technique for the single-input single-output system. The resulting closed-loop response tracks the reference value and leads to significant improvements in the transient response over the open-loop system. The NARMAX controller enhances the lift coefficient from 0.787 for the uncontrolled case to 1.315 for the controlled case with an increase of 67.1%. 
    more » « less