Adjoint-assisted Pareto front tracing in aerodynamic and conjugate heat transfer shape optimization. (English) Zbl 1521.76286
Summary: In this paper, a prediction-correction algorithm, built on the method proposed in [S. Schmidt and V. Schulz, Pac. J. Optim. 4, No. 2, 243–258 (2008; Zbl 1163.90742)], uses the adjoint method to trace the Pareto front. The method is initialized by a point on the Pareto front obtained by considering one of the objectives only. During the prediction and correction steps, different systems of equations are derived and solved by treating the Karush-Kuhn-Tucker (KKT) optimality conditions in two different ways. The computation of second derivatives of the objective functions (Hessian matrix) which appear in the equations solved to update the design variables is avoided. Instead, two approaches are used: (a) the computation of Hessian-vector products driving a Krylov subspace solver and (b) the approximations of the Hessian via Quasi-Newton methods. Three different variants of the prediction-correction method are developed, applied to 2D aerodynamic shape optimization problems with geometrical constraints and compared in terms of computational cost. It is shown that the inclusion of the prediction step in the algorithm and the use of Quasi-Newton methods with Hessian approximations in both steps has the lowest computational cost. This method is, then, used to compute the Pareto front in a 3D conjugate heat transfer shape optimization problem, with the total pressure losses and max. solid temperature as the two contradicting objectives.
MSC:
| 76G25 | General aerodynamics and subsonic flows |
| 76N25 | Flow control and optimization for compressible fluids and gas dynamics |
Keywords:
Pareto front; continuous adjoint method; aerodynamic shape optimization; conjugate heat transfer optimizationCitations:
Zbl 1163.90742References:
| [1] | Schmidt, S.; Schulz, V., Pareto-curve continuation in multi-objective optimization., Pac J Optim, 4, 2, 243-257 (2008) · Zbl 1163.90742 |
| [2] | Fonseca, C.; Fleming, P., An overview of evolutionary algorithms in multiobjective optimization, Evol Comput, 3, 1, 1-16 (1995) |
| [3] | Kapsoulis, D.; Tsiakas, K.; Trompoukis, X.; Asouti, V.; Giannakoglou, K., A PCA-assisted hybrid algorithm combining EAs and adjoint methods for CFD-based optimization, Appl Soft Comput, 73, 520-529 (2018) |
| [4] | Pironneau, O., On optimum profiles in stokes flow, J Fluid Mesh, 59, 117-128 (1973) · Zbl 0274.76022 |
| [5] | Jameson, A., Aerodynamic design via control theory, J Sci Comput, 3, 233-260 (1988) · Zbl 0676.76055 |
| [6] | Anderson, W.; Venkatakrishnan, V., Aerodynamic design optimization on unstructured grids with a continuous adjoint formulation, Comput Fluids, 28, 4-5, 443-480 (1999) · Zbl 0968.76074 |
| [7] | Giles, M.; Pierce, N., Adjoint equations in CFD: duality, boundary conditions and solution behaviour, AIAA Paper 1997-1850, 13th fluid dynamics conference. New Orleans, LA (1997) |
| [8] | Peter, J.; Dwight., R., Numerical sensitivity analysis for aerodynamic optimization: a survey of approaches, Comput Fluids, 39, 3, 373-391 (2010) · Zbl 1242.76301 |
| [9] | Papoutsis-Kiachagias, E.; Giannakoglou, K., Continuous adjoint methods for turbulent flows, applied to shape and topology optimization: industrial applications, Arch Comput Methods Eng, 23, 2, 255-299 (2016) · Zbl 1348.76054 |
| [10] | Zingg, D.; Nemec, M.; Pulliam, T., A comparative evaluation of genetic and gradient- based algorithms applied to aerodynamic optimization, Eur J Comput Mech, 17, 103-126 (2012) · Zbl 1292.76062 |
| [11] | Das, I.; Dennis, J., A closer look at the drawbacks of minimizing weighted sums of objectives for Pareto set generation in multicriteria optimization problems, Struct Optim, 14, 63-69 (1997) |
| [12] | Kim, I.; de Weck, O., Adaptive weighted-sum method for bi-objective optimization: Pareto front generation, Struct Multidiscip Optim, 29, 149-158 (2005) |
| [13] | Shankaran, S.; Barr, B., Efficient gradient-based algorithms for the construction of Pareto fronts, ASME Paper GT2011-45069. Vancouver, Canada (2011) |
| [14] | Fliege, J.; Svaiter, B., Steepest descent methods for multicriteria optimization, Math Methods Oper Res, 51, 479-494 (2000) · Zbl 1054.90067 |
| [15] | Das, I.; Dennis, J. E., Normal-boundary intersection: a new method for generating the Pareto surface in nonlinear multicriteria optimization problems, SIAM J Optim, 8, 3, 631-657 (1998) · Zbl 0911.90287 |
| [16] | Gembicki, F.; Haimes, Y., Approach to performance and sensitivity multiobjective optimization: the goal attainment method, IEEE Trans Automat Contr, 20, 6, 769-771 (1975) |
| [17] | Schy, A.; Giesy, D., Multiobjective insensitive design of airplane control systems with uncertain parameters, AIAA Paper 81-1818, guidance and control conference. Albuquerque, NM, USA (1981) |
| [18] | Schy, A.; Giesy, D., Tradeoff studies in multiobjective insensitive design of airplane control systems, AIAA Paper 83-2273, guidance and control conference. Gatlinburg, TN, USA (1983) |
| [19] | Nocedal, J.; Wright, S., Numerical optimization. (1999), Springer: Springer New York · Zbl 0930.65067 |
| [20] | Vasilopoulos, I.; Asouti, V.; Giannakoglou, K.; Meyer, M., Gradient-based Pareto front approximation applied to turbomachinery shape optimization, Eng Comput (2019) |
| [21] | Gkaragkounis, K.; Papoutsis-Kiachagias, E.; Asouti, V.; Giannakoglou, K., Adjoint-based Pareto front tracing in aerodynamic shape optimization, ICCFD10-2018-322. Barcelona, Spain (2018) |
| [22] | Saad, Y.; Schultz, M., GMRES: A generalized minimal residual algorithm for solving nonsymmetric linear systems, SIAM J Sci StatComput, 7, 3, 856-869 (1986) · Zbl 0599.65018 |
| [23] | Ghavami Nejad, M.; Papoutsis-Kiachagias, E.; Giannakoglou, K., Aerodynamic shape optimization using the truncated Newton method and continuous adjoint, ECCOMAS Congress 2016, VII European congress on computational methods in applied sciences and engineering. Crete Island, Greece (2016) · Zbl 1348.76054 |
| [24] | Papadimitriou, D.; Giannakoglou, K., Aerodynamic design using the truncated Newton algorithm and the continuous adjoint approach, Int J Numer Methods Fluids, 68, 6, 724-739 (2012) · Zbl 1452.76209 |
| [25] | Papadimitriou, D.; Giannakoglou, K., Direct, adjoint and mixed approaches for the computation of Hessian in airfoil design problems, Int J Numer Methods Fluids, 56, 10, 1929-1943 (2008) · Zbl 1141.76058 |
| [26] | Gkaragkounis, K.; Papoutsis-Kiachagias, E.; Giannakoglou, K., The continuous adjoint method for shape optimization in conjugate heat transfer problems with turbulent incompressible flows, Appl Therm Eng, 140, 351-362 (2018) |
| [27] | Papoutsis-Kiachagias, E.; Magoulas, N.; Mueller, J.; Othmer, C.; Giannakoglou, K., Noise reduction in car aerodynamics using a surrogate objective function and the continuous adjoint method with wall functions, Comput Fluids, 122, 223-232 (2015) · Zbl 1390.76826 |
| [28] | Papadimitriou, D.; Giannakoglou, K., A continuous adjoint method for the minimization of losses in cascade viscous flows, 44th AIAA Aerospace sciences meeting and exhibit (2012) |
| [29] | Patankar, S.; Spalding, D., Calculation procedure for heat, mass and momentum transfer in three-dimensional parabolic flows, Int J Heat Mass Transf, 15, 1787-1806 (1972) · Zbl 0246.76080 |
| [30] | Deb, K.; Pratap, A.; Agarwal, S.; Meyarivan, T., A fast and elitist multiobjective genetic algorithm: NSGA-II, IEEE Trans Evol Comput, 6, 2, 182-197 (2002) |
| [31] | Spalart, P.; Allmaras, S., A one-equation turbulence model for aerodynamic flows., AIAA Paper 1992-0439 (1992) |
| [32] | Tucker, P., Differential equation-based wall distance computation for DES and RANS, J Comput Phys, 190, 229-248 (2003) · Zbl 1236.76028 |
| [33] | Zymaris, A.; Papadimitriou, D.; Giannakoglou, K.; Othmer, C., Continuous adjoint approach to the Spalart-Allmaras turbulence model for incompressible flows, Comput Fluids, 38, 8, 1528-1538 (2009) · Zbl 1242.76064 |
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.
