Gradient-based Pareto front approximation applied to turbomachinery shape optimization

Multi-objective optimization has been rising in popularity, especially within an industrial environment, where several cost functions often need to be considered during the design phase. Traditional gradient-free approaches, such as evolutionary algorithms, can be employed to compute a front of equally suitable solutions a designer can choose from, with a high computational cost though, particularly in high-dimensional design spaces. In this paper, a gradient-based algorithm is developed for efficiently tracing the Pareto front in bi-objective aerodynamic shape optimization problems, where an adjoint method is used for the computation of the objective functions’ gradients with respect to the design variables. After obtaining a starting point on the front, a prediction–correction approach is employed to compute new Pareto points. Satisfying the Karush–Kuhn–Tucker conditions provides a prediction for the next point, which is, then, corrected by solving a minimum distance problem. The prediction step, though, requires the costly computation of the Hessian matrix. This is avoided here using the BFGS (Nocedal and Wright in Numerical optimization, Springer, New York, 2006 ) technique. The proposed method is first demonstrated in a less expensive lift/drag optimization of an isolated airfoil and then applied to the bi-objective optimization of a 3D compressor stator. The extension of the proposed method to cases with more than two objectives is straightforward, on condition that an algorithm is found to coordinate the way the Pareto front is swept.

[1]  Marcus Meyer,et al.  Linking Parametric CAD with Adjoint Surface Sensitivities , 2016 .

[2]  Stephen J. Wright,et al.  Numerical Optimization , 2018, Fundamental Statistical Inference.

[3]  Kyriakos C. Giannakoglou,et al.  Efficient PCA-driven EAs and metamodel-assisted EAs, with applications in turbomachinery , 2014 .

[4]  P. Spalart A One-Equation Turbulence Model for Aerodynamic Flows , 1992 .

[5]  Marcus Meyer,et al.  CAD-Based Aerodynamic Optimization of a Compressor Stator Using Conventional and Adjoint-Driven Approaches , 2017 .

[6]  Kyriakos C. Giannakoglou,et al.  A Two–Step Mesh Adaptation Tool Based on RBF with Application to Turbomachinery Optimization Loops , 2018, Computational Methods in Applied Sciences.

[7]  Pierre Moinier,et al.  Edge-Based Multigrid and Preconditioning for Hybrid Grids , 1999 .

[8]  K. Giannakoglou,et al.  Aerodynamic design using the truncated Newton algorithm and the continuous adjoint approach , 2012 .

[9]  Marcus Meyer,et al.  CAD-Based Adjoint Shape Optimisation of a One-Stage Turbine With Geometric Constraints , 2015 .

[10]  Leigh Lapworth,et al.  PADRAM: Parametric Design and Rapid Meshing System for Turbomachinery Optimisation , 2003 .

[11]  T. Pulliam,et al.  A comparative evaluation of genetic and gradient-based algorithms applied to aerodynamic optimization , 2008 .

[12]  B. Launder,et al.  The numerical computation of turbulent flows , 1990 .

[13]  Marius Swoboda,et al.  A Full Parametric Model for Turbomachinery Blade Design and Optimisation , 2004, DAC 2004.

[14]  John E. Dennis,et al.  Normal-Boundary Intersection: A New Method for Generating the Pareto Surface in Nonlinear Multicriteria Optimization Problems , 1998, SIAM J. Optim..

[15]  I. Y. Kim,et al.  Adaptive weighted-sum method for bi-objective optimization: Pareto front generation , 2005 .

[16]  Les A. Piegl,et al.  The NURBS Book , 1995, Monographs in Visual Communication.

[17]  Matthias Ehrgott,et al.  Multicriteria Optimization (2. ed.) , 2005 .

[18]  S. Shankaran,et al.  Efficient Gradient-Based Algorithms for the Construction of Pareto Fronts , 2011 .

[19]  Jasbir S. Arora,et al.  Survey of multi-objective optimization methods for engineering , 2004 .

[20]  Kyriakos C. Giannakoglou,et al.  Unsteady CFD computations using vertex‐centered finite volumes for unstructured grids on Graphics Processing Units , 2011 .

[21]  Michael Dellnitz,et al.  On the hierarchical structure of Pareto critical sets , 2019, J. Glob. Optim..

[22]  Michael B. Giles,et al.  On the use of Runge-Kutta time-marching and multigrid for the solution of steady adjoint equations , 2000 .

[23]  Ulf Schlichtmann,et al.  A Successive Approach to Compute the Bounded Pareto Front of Practical Multiobjective Optimization Problems , 2009, SIAM J. Optim..

[24]  R. K. Ursem Multi-objective Optimization using Evolutionary Algorithms , 2009 .

[25]  Shahrokh Shahpar,et al.  Aerodynamic Optimization of High-Pressure Turbines for Lean-Burn Combustion System , 2013 .

[26]  Kyriakos C. Giannakoglou,et al.  Shape Optimization of Wind Turbine Blades Using the Continuous Adjoint Method and Volumetric NURBS on a GPU Cluster , 2019 .

[27]  A. Strauß Theory Of Wing Sections Including A Summary Of Airfoil Data , 2016 .

[28]  A. Peirce Computer Methods in Applied Mechanics and Engineering , 2010 .

[29]  Matthias Ehrgott,et al.  Multicriteria Optimization , 2005 .

[30]  Stefan Volkwein,et al.  POD-based error control for reduced-order bicriterial PDE-constrained optimization , 2017, Annu. Rev. Control..

[32]  Marcus Meyer,et al.  Numerical and Experimental Investigations on Optimized 3D Compressor Airfoils , 2018, Volume 2A: Turbomachinery.

[33]  Michael Dellnitz,et al.  Multiobjective Optimal Control Methods for the Navier-Stokes Equations Using Reduced Order Modeling , 2018, Acta Applicandae Mathematicae.