Aerodynamic Shape Optimization with Time Spectral Flutter Adjoint

Flutter onset characteristic are an important consideration in commercial airliner design. Previous work with high-fidelity aerostructural optimization has shown a tendency for optimization algorithms to produce unrealistically high span designs, particularly when maximizing range or minimizing fuel burn, in an effort to maximize aerodynamic efficiency. In order to constraint this tendency, we propose to include a flutter constraint to these optimization problems. To be meaningful for this class of optimization problem, the flutter model should be able to predict the strong nonlinearity in the transonic regime, as commercial airliners largely operate in this regime. In addition, these problems typically include a large number of design variables, therefore we use gradient-based optimization algorithms, requiring the constraint formulation to be continuous, differentiable and have an efficient method for computing the gradients of the constraint. In this paper, we apply Euler time-spectral computational fluid dynamics methods to model the flutter constraint and we propose a coupled adjoint method to calculate the constraint sensitivity with respect to the design variables. The coupled adjoint method has the advantage that the gradient evaluation time is independent of the number of design variables. In the literature, the harmonic balance method has been used to model flutter constraints. A two-level adjoint formulation has been proposed to evaluate the sensitivity of flutter onset velocity with respect to the design variables. Such method requires solving the harmonic balance aerodynamic adjoint O(NCSD) times where NCSD is the structural degree of freedom multiplied by number of time instances. On the contrary, we propose the coupled adjoint method which directly deals with the whole aeroelastic system and solves only 1 adjoint equation. We verify the adjoint sensitivity computation with the finite difference method. Finally, we present a MDO problem for the classic Isogai case in which we maximize the flutter velocity index with respect to aerodynamic shape design variables. We gain a 10.9% increase of the flutter velocity index through this optimization.

[1]  Joaquim R. R. A. Martins,et al.  Computational Modeling of Flutter Constraint for High-Fidelity Aerostructural Optimization , 2019, AIAA Scitech 2019 Forum.

[2]  J. Martins,et al.  Buffet-Onset Constraint Formulation for Aerodynamic Shape Optimization , 2017 .

[3]  Antony Jameson,et al.  Time Spectral Method for Periodic Unsteady Computations over Two- and Three- Dimensional Bodies , 2005 .

[4]  Joaquim R. R. A. Martins,et al.  High-fidelity aerostructural optimization considering buffet onset , 2015 .

[5]  Qiqi Wang,et al.  Unsteady Aerostructure Coupled Adjoint Method for Flutter Suppression , 2015 .

[6]  Carol D. Wieseman,et al.  Aeroelastic Tailoring of Transport Wings Including Transonic Flutter Constraints , 2015 .

[7]  Qiqi Wang,et al.  A Monte Carlo method for solving unsteady adjoint equations , 2008, J. Comput. Phys..

[8]  Seulgi Yi,et al.  High Fidelity Prediction of Flutter/LCO using Time Spectral Method , 2018 .

[9]  John T. Hwang,et al.  Review and Unification of Methods for Computing Derivatives of Multidisciplinary Computational Models , 2013 .

[10]  Jeffrey P. Thomas,et al.  Discrete Adjoint Method for Nonlinear Aeroelastic Sensitivities for Compressible and Viscous Flows , 2013 .

[11]  Earl H. Dowell,et al.  A Fixed Point Iteration Approach for Harmonic Balance Based Aeroelastic Computations , 2018 .

[12]  Michael A. Saunders,et al.  SNOPT: An SQP Algorithm for Large-Scale Constrained Optimization , 2002, SIAM J. Optim..

[13]  Jeffrey P. Thomas,et al.  Discrete Adjoint Method for Aeroelastic Design Optimization , 2014 .

[14]  Antony Jameson,et al.  The computational efficiency of non-linear frequency domain methods , 2006, J. Comput. Phys..

[15]  Joaquim R. R. A. Martins,et al.  A coupled Newton-Krylov time spectral solver for flutter prediction , 2018 .

[16]  C. Mader,et al.  Derivatives for Time-Spectral Computational Fluid Dynamics Using an Automatic Differentiation Adjoint , 2012 .

[17]  Eli Livne,et al.  Design-Oriented Computational Fluid Dynamics-Based Unsteady Aerodynamics for Flight-Vehicle Aeroelastic Shape Optimization , 2015 .

[18]  Hang Li,et al.  A novel approach for flutter prediction of pitch–plunge airfoils using an efficient one-shot method , 2018, Journal of Fluids and Structures.

[19]  Graeme J. Kennedy,et al.  Scalable Parallel Approach for High-Fidelity Steady-State Aeroelastic Analysis and Adjoint Derivative Computations , 2014 .

[20]  Jeffrey P. Thomas,et al.  Nonlinear Inviscid Aerodynamic Effects on Transonic Divergence, Flutter, and Limit-Cycle Oscillations , 2001 .

[21]  Simao Marques,et al.  Prediction of Transonic LCOs using an Aeroelastic Harmonic Balance Method , 2014 .

[22]  Joaquim R. R. A. Martins,et al.  pyOpt: a Python-based object-oriented framework for nonlinear constrained optimization , 2011, Structural and Multidisciplinary Optimization.

[23]  Bret Stanford,et al.  Economical Unsteady High-Fidelity Aerodynamics for Structural Optimization with a Flutter Constraint , 2017 .

[24]  Joaquim R. R. A. Martins,et al.  Multipoint High-Fidelity Aerostructural Optimization of a Transport Aircraft Configuration , 2014 .

[25]  Eric Blades,et al.  A fast mesh deformation method using explicit interpolation , 2012, J. Comput. Phys..

[26]  Joaquim R. R. A. Martins,et al.  Flutter and post-flutter constraints in aircraft design optimization , 2019, Progress in Aerospace Sciences.

[27]  Juan J. Alonso,et al.  Fully-implicit time-marching aeroelastic solutions , 1994 .

[28]  Xian-Da Zhang,et al.  Matrix Analysis and Applications , 2017 .

[29]  Jeffrey P. Thomas,et al.  Computation of Unsteady Nonlinear Flows in Cascades Using a Harmonic Balance Technique , 2002 .

[30]  Graeme J. Kennedy,et al.  Application of a Time-Accurate Aeroelastic Coupling Framework to Flutter-Constrained Design Optimization , 2018, 2018 Multidisciplinary Analysis and Optimization Conference.

[31]  Joaquim R. R. A. Martins,et al.  The complex-step derivative approximation , 2003, TOMS.

[32]  J. Alonso,et al.  ADjoint: An Approach for the Rapid Development of Discrete Adjoint Solvers , 2006 .

[33]  Joaquim R. R. A. Martins,et al.  High-Fidelity Aerostructural Design Optimization of a Supersonic Business Jet , 2002 .

[34]  Joaquim R. R. A. Martins,et al.  Towards gradient-based design optimization of flexible transport aircraft with flutter constraints , 2014 .