Time Spectral Adjoint Based Design for Flutter and Limit Cycle Oscillation Suppression

When designing aircraft wings shapes, it is important to ensure that the flight envelope does not overlap with regions of flutter or Limit Cycle Oscillation (LCO). A quick assessment of these dynamic aeroelastic for various design candidates is key to successful design. Flutter based design requires the sensitivity of flutter parameters to be known with the respect of design parameters. Traditionally, frequency domain based methods have been used to predict flutter characteristics and its sensitivity. However, this approach is only applicable for linear or linearized models and cannot be applied to systems undergoing LCO or other nonlinear effects. Though the time accurate approach can be implemented to overcome this problem, it is computationally expensive. Also, the unsteady adjoint formulation for sensitivity analysis, requires the state and adjoint variables to be stored at every time step, which prohibitively increases the memory requirement. In this work, these problems have been overcome by implementing a time spectral method based approach to compute flutter onset, LCOs and their design sensitivities in a computationally efficient manner. The time spectral based formulation approximates the solution as a discrete Fourier series and directly solves for the periodic steady state, leading to a steady formulation. This can lead to the time spectral approach to be faster than the time accurate approach. More importantly, the steady formulation of the time spectral method also eliminates the memory issues faced by the unsteady adjoint formulation. The time spectral based flutter/LCO prediction method was used to predict flutter and LCO characteristics of the AGARD 445.6 wing and pitch/plunge airfoil section with NACA 64A010 airfoil. Furthermore, the adjoint based sensitivity analysis was used to carry out aerodynamic shape optimization, with an objective of maximizing the flutter velocity with and without constraints on the drag coefficient. The resulting designs show significant increase in the flutter velocity and the corresponding LCO velocity profile. The resulting airfoils display a greater sensitivity to the transonic shock which in turn leads to greater aerodynamic damping and hence leading to an increase in flutter velocity.