A New Modified Harmony Search Optimization Algorithm for Evaluating Airfoil Shape Parameterization Methods and Aerodynamic Optimization

Document Type : Research Paper


1 M.Sc. Student, University of Birjand, Birjand, South Khorasan, 9718776391, Iran

2 Corresponding Author, Associate Professor, Department of Mechanical Engineering, University of Birjand, Birjand, South Khorasan, 97175/615, Iran, mmjafarian@birjand.ac.ir

3 Associate Professor, Department of Mechanical Engineering, University of Birjand, Birjand, South Khorasan, 97175/615, Iran


In this work, a modification has been made to increase the efficiency and convergence of the harmony search algorithm. Then, the capability of this amendment was investigated by applying it to the following aerodynamic problems for the first time. First, the methods of airfoil shape parametrization (Bezier curves, Parsec method, and NACA 4-digit airfoil) were investigated using an inverse optimization design by the present modified harmony search optimization algorithm. Then, inverse and direct optimization of an airfoil were carried out by the modified algorithm. Aerodynamic analysis of the problem was obtained using compressible Reynolds-Averaged Navier-Stokes (RANS) equations along with the Spalart-Allmaras turbulence model. Results showed that the Bezier curves and the Parsec method have higher flexibility than the NACA 4-digit airfoil. The Parsec method was introduced as the best approach, because of fewer control parameters. The inverse optimization results showed that the present airfoil shape optimization set can obtain the target shape with high accuracy. The Direct optimization with a maximum lift to drag ratio target function revealed that the shock waves significantly weaken at the optimum airfoil. Generally, the results obtained verify that using the modified harmony search algorithm together with the Parsec method provides a powerful tool for direct and inverse aerodynamic optimization.


Main Subjects

[1] Avinash, G.S., and Anil Lal, S., "Inverse Design of Airfoil using Vortex Element Method", International Journal of Fluid Machinery and Systems, Vol. 11(2), pp. 163-170, (2018).
[2] Sobester, A., and Thomas, B.T., "The Quest for a Truly Parsimonious Airfoil Parameterization Scheme", ICAS Congress Including the 8th AIAA ATIO Conference, Sep. 14-19, Anchorage, Alaska, pp. 1-9, (2008).
[3] Shahrokhi, A., and Jahangirian, A., "The Effects of Shape Parameterization on the Efficiency of Evolutionary Design Optimization for Viscous Transonic Airfoils", Journal of Aerospace Science and Technology (JAST), Vol. 5, pp. 35-43, (2008).
[4] Boehm, W., "Bezier Presentation of Airfoils", Computer Aided Geometric Design, Vol. 4, pp. 17-22, (1987).
[5] Jahangirian, A., and Shahrokhi, A., "Inverse Design of Transonic Airfoils using Genetic Algorithm and a New Parametric Shape Method", Inverse Problems in Science and Engineering, Vol. 17, pp. 681-699, (2009).
[6] H´ajek, J., "Aerodynamic Optimization of Airfoils and Wings using Fast Solvers", Ph.D. Thesis, Charles University, Prague, Czech, (2009).
[7] Wauquiez, C., "Shape Optimization of Low-speed Airfoils using MATLAB and Automatic Differentiation", Licentiate’s Thesis, Department of Numerical Analysis and Computing Science University of Royal Institute of Technology, Stockholm, Sweden (2000).
[8] Geem, Z.W., "Music-inspired Harmony Search Algorithm", Studies in Computational Intelligence, Vol. 191, pp. 1-14,  (2009).
[9] Soemarwoto, B.I., "The Variational Method for Aerodynamic Optimization using the Navier-stokes Equations", Icase Report, pp. 71-97, (1997).
[10] Gardner, B.A., and Selig, M.S., "Airfoil Design using a Genetic Algorithm and an Inverse Method", 41st Aerospace Sciences Meeting and Exhibit, 6-9 January, Reno, Nevada, (2003).
[11] De’ Michieli Vitturi, M., and Beux, F., "A Discrete Gradient-based Approach for Aerodynamic Shape Optimization in Turbulent Viscous Flow", Finite Elements in Analysis and Design, Vol. 43, pp. 68–80, (2006).
[12] Shahrokhi, A., and Jahangirian, A., "Airfoil Shape Parameterization for Optimum Navier-stokes Design with Genetic Algorithm", Aerospace Science and Technology, Vol. 11, pp. 443-450, (2007).
[13] Xu, Z., and Xia, J., "Aerodynamic Optimization Based on Continuous Adjoint Method for a Flexible Wing", International Journal of Aerospace Engineering, Article ID 4706925, (2016).
[14] Alves, R., Goncalves, L., and Aguiar, J., "Parsec Parameterization Methodology for Enhancing Airfoils Geometry using PSO Algorithm", CILAMCSE 2016, Proceedings of the XXXVII Iberian Latin-American Congress on Computational Methods in Engineering, Braslia, DF, Brazil, November 6-9, (2016).
[15] Ebrahimi, M., and Jahangirian, A.R., "Airfoil Shape Optimization with Adaptive Mutation Genetic Algorithm", Journal of Aerospace Science and Technology (JAST), Vol. 11(1), pp. 47-59, (2017).
[16] Yu, Y., Lyu, Z., Xu, Z., and Martins, J., "On the Influence of Optimization Algorithm and Initial Design on Wing Aerodynamic Shape Optimization", Aerospace Science and Technology, Vol. 75, pp. 183-199, (2018).
[17] Fesanghary, M., Damangir, E., and Soleimani, I., "Design Optimization of Shell and Tube Heat Exchangers using Global Sensitivity Analysis and Harmony Search Algorithm", Applied Thermal Engineering, Vol. 29(5–6), pp. 1026-1031, (2009).
[18] Kaveh, A., Shakouri, A., and Abadi, M., "Cost Optimization of a Composite Floor System using an Improved Harmony Search Algorithm", Journal of Constructional Steel Research, Vol. 66(5), pp. 664-669, (2010).
[19] Yousefi, M., Enayatifar, R., Darus, A.N., and Abdullah, A.H., "Optimization of Plate-fin Heat Exchangers by an Improved Harmony Search Algorithm", Applied Thermal Engineering, Vol. 50(1), pp. 877-885, (2013).
[20] Boryczka, U., and Szwarc, K., "The Adaptation of the Harmony Search Algorithm to the ATSP with the Evaluation of the Influence of the Pitch Adjustment Place on the Quality of Results", Journal of Information and Telecommunication, Vol. 3(1), pp. 2-18, (2019).
[21] Yarmohamadi, H., Zhang, Q., Jiao, P., and Alavi, A.H., "An Enhanced Adaptive Global-best Harmony Search Algorithm for Continuous Optimization Problems", Engineering Reports Journal, Vol. 2(11), pp. 1-19, (2020).
[22] Abu Doush, I., and Santos, E., "Best Polynomial Harmony Search with Best β-hill Climbing Algorithm", Journal of Intelligent Systems, Vol. 30(1), pp. 1-17, (2020).
[23]      Zhang, Y., Li, J., and Li, L., "A Reward Population-based Differential Genetic Harmony Search Algorithm", Algorithms, Vol. 15(23), pp. 1-19, (2022).
[24] Iglesias, A., “Bezier Curves and Surfaces", Department of Applied Mathematics and Computational Sciences University of Cantabria UC—CAGD Group”, https://ocw.unican.es/pluginfile.php/1949/course/section/2243/01-Bezier.pdf, (2001).
[25] Endo, M., "Wind Turbine Airfoil Optimization by Particle Swarm Method", M.Sc. Thesis, Department of Mechanical and Aerospace Engineering University of Case Western Reserve, Cleveland, Ohio, USA, (2011).
[26] Spalart, P.R., and Allmaras, S.R., "A One-equation Turbulence Model for Aerodynamic Flows", La Recherche Aerospatiale, Vol. 1, pp. 5-21, (1994).
[27] Jameson, A., and Schmidt, W., "Numerical Solution of the Euler Equations by Finite Volume Methods using Runge-Kutta Time-stepping Schemes", American Institute of Aeronautics and Astronautics (AIAA Journal), Vol. 81, pp. 1259, (1981).
[28] Jameson, A., "The Present Status, Challenges, and Future Developments in Computational Fluid Dynamics", Twelfth Australasian Fluid Mechanics Conference, Dec. 7-10,  The University of Sydney, Sydney, Australia, pp. 423-468,  (1995).
[29] Lee, K.S., and Geem, Z.W., "A New Meta-heuristic Algorithm for Continuous Engineering Optimization: Harmony Search Theory and Practice", Computer Methods in Applied Mechanics and Engineering, Vol. 194, pp. 3902–3933, (2005).
[30] Pan, Q.K., Suganthan, P.N., Tasgetiren, M.F., and Liang, J.J., "A Self-adaptive Global Best Harmony Search Algorithm for Continuous Optimization Problems", Applied Mathematics and Computation, Vol. 216, pp. 830–848, (2010).