Surrogate-based aerodynamic shape optimization for delaying airfoil dynamic stall using Kriging regression and infill criteria

Abstract

The dynamic stall phenomenon is characterized by the formation of a leading-edge vortex, which is responsible for adverse aerodynamic forces and moments adversely impacting the structural strength and life of a system. Aerodynamic shape optimization (ASO) provides a cost-effective approach to delay or mitigate the dynamic stall characteristics. Unfortunately, ASO requires multiple evaluations of accurate but time-consuming computational fluid dynamics (CFD) simulations to produce optimum designs rendering the optimization process computationally costly. The current work proposes a surrogate-based optimization (SBO) technique to alleviate the computational burden of ASO to delay and mitigate the deep dynamic stall characteristics of airfoils. In particular, the Kriging regression surrogate model is used for approximating the objective and constraint functions. The airfoil geometry is parametrized using six PARSEC parameters. The objective and constraint functions are evaluated with the unsteady Reynolds-averaged Navier-Stokes equations with a C-grid mesh topology and Menter's shear stress transport turbulence model. The approach is demonstrated on a vertical axis wind turbine airfoil at a Reynolds number of 135,000 and a Mach number of 0.1 undergoing a sinusoidal oscillation with a reduced frequency of 0.05. The surrogate model is constructed with 60 initial samples and further refined with 20 infill samples using expected improvement. The generated surrogate model is validated with the normalized root mean square error based on 20 test data samples. The refined surrogate model is utilized for finding the optimal design using multi-start gradient-based search. The optimal airfoil has a higher thickness, larger leading-edge radius, and an aft camber compared to the baseline. These geometric shape changes delay the dynamic stall angle by over 3∘ and reduces the severity of the pitching moment coefficient fluctuation. Finally, global sensitivity analysis is conducted on the optimal design using Sobol' indices revealing the most influential shape variables and their interaction effects impacting the airfoil dynamic stall characteristics.