In this section, we organize the papers based on their focus on the structural characteristics of the airfoil. This classification takes into account whether the airfoil is designed as a single or multi-element body as well as whether it is intended to be rigid or flexible. In conventional aircraft wings, airfoils are designed to be optimal for a specific flight condition. However, when analyzing flight conditions for different regimes such as takeoff, landing, or cruise flight, researchers often investigate multi-element geometries with components such as flaps or slats [
3,
4,
19] (
Figure 3a) In recent years, the analysis of flexible geometries has gained significant interest among researchers. This innovative design concept aims to enhance aerodynamic and energy efficiency during flight by actively adjusting the shape of the wing in real-time. In an article published by Terrence A. Weisshaar [
49] the author highlights that the term "morphing aircraft" has come to be associated with exotic smart materials, which is too restrictive. Morphing objectives include efficient, multipoint adaptability achieved by state change, including geometric reconfiguration (
Figure 3b).
3.1.1. Paper Review
A summary of the papers analyzed is presented below:
Effective directional control of an airplane demands adjustments to the angle of attack of the wings. However, at high angles of attack, the torque required to operate the mobile control surface increases significantly, limiting the functionality of these systems. One solution is to add a separate rigid element (tab) to the mobile control surface to improve its capabilities (
Figure 5a). However, this approach creates a gap between the control surface and the tab, which can disrupt the surrounding flow. To address this issue, Yu-Tai Lee et al. [
8] developed a concept that uses a single element geometry made of flexible material (such as shape memory alloy) that utilizes electrical power to manipulate the control surface, eliminating the need for hydraulics and the gap between the flap and the tab. (
Figure 5b).
Orthotropic materials display unique mechanical properties in different directions, making them ideal for developing innovative structures with intriguing possibilities. These materials have properties that vary along the orthogonal axes that define their section or volume. The concept of orthotropy serves as the foundation for developing composite materials. Paolo Bettini et al.[
7] have developed complex composite cellular structures that are capable of undergoing large displacements and withstanding limited stress on their components. Numerical simulations are used to investigate and compare the advantages of using orthotropic materials, with the results being compared to those obtained from metallic prototypes (
Figure 6).
Yokozeki et al. [
15] developed a morphing airfoil with variable geometry by utilizing corrugated structures and wire tension actuation, which was tested in a wind tunnel. The flexible seamless flap was actuated by a wire system in the leading part of the body, using C-shaped circular sections and vertically straight sections for the corrugated structure with uniform curves. The wire was attached to the trailing edge region and passed through the corrugated structure near the lower wing surface, with the other end connected to a servomotor mounted inside the wing. (
Figure 7).
Roelof Vos et al. [
9] have developed a warp-controlled wing twisting prototype based on the concept of controlled wing twisting. By introducing variable torsion through rotating ribs around a fixed aluminum rod, they were able to create a lightweight and robust structure without relying on extremely powerful actuators. The wing itself is made of carbon fiber skin. This study demonstrated that active wing twist can effectively control the rolling motion of an aircraft, maximizing the lift-drag ratio under various flight conditions. While the weight of the wings is lighter compared to those with hinged morphing mechanisms made of aluminum or steel, the weight of the skin and driving mechanisms in the whole morphing structure is still too heavy for a compliant wing aircraft.
In a study by Hongda Li et al. [
17], an innovative structural solution was developed for an adaptive variable camber compatible wing. The proposed solution features a skin that can change thickness and a trailing edge transformation mechanism designed based on a new type of artificial muscles (
Figure 8). Although artificial muscles, particularly Pneumatic Artificial Muscles (PAMs), are not new in aerospace applications, they have been studied extensively for use in driving flaps or morphing mechanisms.
As mentioned in the previous paragraph, the weight of the skin is a significant contributor to the overall weight of morphing wing structures. Consequently, the development of a lightweight, flexible skin is a key challenge in creating successful morphing wings. In a study by Maurizio Arena et al. [
25], an innovative structural concept was developed and tested on a two-bay trailing edge segment consisting of three articulated finger-like ribs and a multi-box arrangement covered by an elastomer-based morphing skin. The primary objective of this study was to understand and evaluate how the morphing skin could affect the dynamic response of the adaptive structural system. (
Figure 9).
Topology optimization is a mathematical method that optimizes material layout within a specific design space to maximize system performance for a given set of loads, boundary conditions, and constraints. In a study by Pedro Gomes et al. [
26], a strategy for density-based topology optimization of fluid-structure interaction problems is proposed to address some of the shortcomings associated with non-stiffness-based designs. The objective of this study is to improve the passive aerodynamic shape adaptation of highly compliant airfoils at multiple operating points. The authors found that better results were obtained using topology optimization in addition to shape optimization because the optimized material layout allowed the airfoil to better adapt to the different fluid loads at different speeds.
In the study conducted by S. A. Meguid et al. [
45] the authors provide a comprehensive review of the underlying technologies related to shape morphing wing concepts and describes the design, analysis, prototyping, and flight testing of a novel morphing wing system for small unmanned aerial vehicles (UAVs). The wing system was designed to be lightweight and flexible, using a camber-adjustable airfoil shape and a flexible-rib system (
Figure 10) driven by servomotor-actuated mechanisms. The designed prototype was manufactured and tested in compliance with aircraft design standards. The success of the project, including flight testing of the flexible rib assembly, is summarized in the paper.
Q. Chanzy et al. in [
46] describes the process of designing, testing, and optimizing a novel set of morphing wings for a 15kg maximum take-off weight Unmanned Aerial Vehicle (
Figure 11). A comparator set of conventional wings was used throughout the study to assess the advantages of adopting the morphing technology. Recent advancements in morphing wing technologies have been driven by the use of composite materials that offer precise control over the stiffness properties and shape changes of the wings. In this study, a careful combination of fiber and foam structural elements was used to control the wing shape changes. The wings have an open slot at the trailing edge, allowing the upper and lower trailing edge surfaces to slide over each other and change their camber. The results demonstrate that the morphing wings provide a controllable aircraft while reducing drag by 40% compared to wings with conventional ailerons in a deflected position.
In an article presented by U. Fasel et al. [
48] the authors discuss the potential of combining continuous carbon fiber composite additive manufacturing and morphing mechanisms to improve the flight performance of aerospace structures. The high design freedom of the manufacturing process allows for tailored fiber placement, fully exploiting the anisotropy and strength of the composite material. Morphing mechanisms, which use the elastic properties of the material to achieve shape changes, are prohibitively expensive to manufacture with conventional processes due to their complex internal topologies. The article presents a proof of concept of the combination of these two technologies, where a complete composite primary- and morphing structure of a fixed-wing drone (
Figure 12)is additively manufactured and successfully flight-tested.
In the papers we have reviewed thus far, the designs have primarily been either structure-based or material-based to achieve desired shape changes. However, some works have taken a combined approach, incorporating both the structural characteristics and material properties into their designs. Some examples involve using material properties to morph specific parts of the section, while employing distinctive structural elements to morph other areas. Additionally, some works utilize smart materials as actuators within an internally structured design, based on a topology optimization approach. This is the case of the work developed by Yaqing Zhang et. al. [
50] The paper describes the design of a variable camber morphing wing using compliant mechanisms and hyperelastic structure topology optimization (
Figure 13). A laminated leading-edge skin is designed to fit the curvature changing phenomenon of the leading edge during deformation. A morphing wing demonstrator was manufactured to testify its deformation capability. The proposed design can realize larger deflection of leading and trailing edges, leading to improved aerodynamic performance and sufficient strength to resist aerodynamic and structural loadings.
Another example of a hybrid design can be found in the research conducted by B. K. S. Woods et al. [
51] The paper presents a new airfoil morphing structure called the Fishbone Active Camber (FishBAC) that uses compliant mechanisms to create significant and continuous changes in airfoil camber and section aerodynamic properties (
Figure 14). The structure consists of a chordwise bending beam spine with stringers branching off to connect it to a pre-tensioned Elastomeric Matrix Composite (EMC) skin surface. Actuators mounted in the D-spar induce bending moments on the spine through tendons, similar to natural musculature systems. The paper focuses on a trailing edge morph configuration, where the compliant spine connects a rigid leading edge D-spar to a solid trailing edge strip. A prototype device is built to explore various aspects of manufacturing this concept and to prove the large deflection capability of the FishBAC.
There are other papers that are also worth commenting on:
The article published by H. L. Şahin et al. [
52] discusses the integration of multiple scientific disciplines in the design and analysis of morphing wings. It presents a new scissor-structural mechanism (SSM) for the trailing edge of morphing wings. Using a novel kinematic synthesis concept, the SSM achieves minimal design error in tracking the airfoil camber.
In [
53] A. V. Ivchenko et al. discusses two proposed design solutions for creating an adaptive compliant aircraft wing frame using a reconfigurable honeycomb core with hidden topological surfaces. The first solution involves lever actuators integrated into the honeycomb structure, while the second solution utilizes mechatronic nodes as actuators. The research also explores the potential of smart materials and artificial muscles.
The potential of combining structural material and fluid pressure to create a morphing structure with integrated shape-change and stiffness functions, is shown in the work presented by S. Vasista et al. [
54] by utilizing a cellular arrangement, distributed actuation and precise control. The moving isosurface threshold topology optimization method was used to design the cell shapes.
We will now refer to some recently published review articles related to this topic:
The first article, authored by S. Ameduri et al. [
23], provides a review of the progress made in morphing research within the European Project "AIRGREEN2" in the Clean-Sky 2 GRA platform. This paper provides an overview of the scientific achievements in morphing research by three Italian research organizations: the Italian Aerospace Research Centre, the University of Naples "Federico II," and the Politecnico of Milano. The development plans followed by each organization are discussed, highlighting their roles and accomplishments. Examples of their projects include the development of wind tunnel models, variable chamber trailing edges, adaptive wing systems, and high-lift generation architectures. The paper also outlines the specific contributions of each organization in various European-level projects. The ongoing activities and upcoming milestones, such as flight tests and critical design reviews, are described. The overall goal of the project is to realize adaptive winglet, wingtip, and control systems for enhanced performance in climb, cruise, and landing.
Another article presented by J. Sun et al. [
55] discusses the advantages of morphing technologies for aircraft and reviews recent developments in the application of smart materials and structures in morphing aircraft. Traditional aircraft are optimized for only one or two flight conditions, not for the entire flight envelope, while morphing aircraft can adapt their shape to optimize performance across a wide range of conditions. The use of "smart" materials and structures, which offer benefits such as high energy density, ease of control, variable stiffness, and strain tolerance, provides new possibilities for designing morphing aircraft. The article focuses on four categories of applications: actuators, sensors, controllers, and structures.
Finally we quote the article, presented by Tuba Majid et al. [
28] provides an overview of the current state-of-the-art in camber morphing mechanisms in the context of design and implementation. It highlights the importance of thoroughly investigating morphing mechanisms and their design principles and methodologies from the implementation and test flight aspects. The article categorizes the camber morphing mechanisms into three categories based on their design features and synthesis methodology: structure-based, material-based, and hybrid. The review provides researchers with choices of materials, actuators, and design methodologies for the implementation of morphing mechanisms with practical aspects of flight performance and applications. However, a majority of these works primarily focuses on structural integrity and aerodynamic performance, lacking comprehensive data on aerodynamic benefits in different flight conditions, energy requirements, morphing control, weight, and fatigue failure.
3.1.2. Numerical Approach and Discretization Criterion
In this section, we will primarily concentrate on the numerical approach utilized to analyze and characterize the fluid’s behavior, along with its parameterization.
In aerodynamics, the flow regime is often described in terms of Mach number. The flow is subsonic when Mach < 1, transonic when Mach ≅ 1, supersonic when Mach > 1 and hipersonic when Mach
1 [
31]. We will focus only in subsonic and transonic regimes.
In the subsonic regime, which is the primary focus of airfoil performance analysis, compressibility effects are generally negligible. Under these conditions, the fluid density can be assumed to be constant, and the flow is considered incompressible. However, as the Mach number exceeds a critical value at certain points on the airfoil surface, the flow velocity surpasses the speed of sound, leading to the formation of shock waves. These shock waves typically form perpendicular to the flow direction and cause the velocity to drop below the speed of sound, transitioning the flow back to the subsonic regime. During typical cruising conditions in the transonic regime, the presence of shock waves becomes a significant contributor to drag. Therefore, from a numerical perspective, a non-viscous flow solver capable of calculating the drag induced by weak shock waves is an appropriate tool for simulating the fluid’s physical behavior, as long as there is no extensive flow separation.
As it was described in section 2.3, the Navier-Stokes equations represent a set of coupled, time-dependent, and nonlinear equations. Analytically solving these equations is often impractical or impossible in many cases. Therefore, numerical analysis becomes necessary to obtain an approximate solution. The Reynolds-averaged Navier-Stokes (RANS) equations provide a reduced representation of the general Navier-Stokes equations. These equations separate the steady-state solution from the time-varying fluctuations within the system. The derivation of the RANS equations employs a mathematical technique called Reynolds decomposition, which separates the expected value of a quantity from its fluctuations. In the context of RANS equations, the solution is divided into a time-independent mean flow velocity and time-varying fluctuations about the mean:
, were
x is the position vector. Following the development of RANS equations, Boussinesq introduced the concept of eddy viscosity, which allows for modeling the momentum transfer caused by turbulent eddies. This concept simplifies the Reynolds stress tensor into a more manageable form [
32]. This is in analogy with how the momentum transfer caused by the molecular motion in a gas can be described by a molecular viscosity.
In 1987, Mark Drela and Michael B. Giles [
1] introduced and validated a viscous-inviscous analysis method called ISES (Implicit Sequential Euler Solver). This method proved to be well-suited for analyzing transonic and low Reynolds number airfoils. They opted for the ISES method over Navier-Stokes analysis, citing the latter’s computational inefficiency for routine design work and lack of proven superiority in terms of accuracy compared to faster zone approaches. The ISES code provides a rapid and dependable solution for addressing strong viscous-inviscid interactions commonly encountered in transonic and low Reynolds number airfoil flows. The presented results demonstrate the capability of the ISES code to accurately predict transitioning separation bubbles and their associated performance losses, as well as the rapid degradation of airfoil performance at decreasing Reynolds numbers.
The aerospace industry heavily relies on computational fluid dynamics (CFD) tools for designing and analyzing flying vehicles, leveraging the exponential growth of computing capabilities. Over time, significant expertise has been gained in using the Reynolds-averaged Navier-Stokes (RANS) equations and coupled integral boundary layer (IBL) methods to predict the flow fields and design trends of aerospace vehicles with attached flow. However, the accuracy of these simulation tools diminishes when analyzing regions with flow separation, such as multi-element configurations, which poses challenges in modeling the phenomenon precisely. In a comprehensive survey conducted by Christopher L. Rumsey et al. [
33], various CFD methods applied to high-lift multi-element computations are examined. The study highlights the crucial need for aircraft designers to have reliable CFD tools capable of accurately predicting the effects resulting from small configuration changes under high-lift conditions, emphasizing the importance of employing both IBL and RANS methods. This conclusion underscores the significance of advancing CFD techniques to enhance the understanding and optimization of complex aerodynamic phenomena in aerospace design.
The design of high-lift systems holds significant importance in aerodynamics, directly impacting the performance of takeoff and landing. However, this design process is complex due to various factors and constraints that must be considered by designers. Evaluating multiple candidate configurations is essential during the initial design phase before finalizing the optimal design. Low-fidelity tools, such as potential flow theory, are commonly employed in the conceptual phase due to their cost-effectiveness and ability to analyze numerous configurations. Nonetheless, relying solely on low-fidelity tools can lead to unexpected outcomes when transitioning to high-fidelity computational fluid dynamics (CFD) analysis in subsequent phases. Hence, it is crucial to enhance the accuracy of low-fidelity tools to avoid surprises when designs undergo CFD assessment. Addressing this concern, Partha Mondal et al. [
19] developed a discrete-vortex-method (DVM)-based model for evaluating functions within a genetic algorithm (GA)-based flapposition-optimization strategy. DVM is renowned for effectively capturing camber effects and lifting surface interference, making it a suitable approach for optimizing high-lift systems. By integrating DVM into the optimization process, the study aims to enhance the fidelity of low-fidelity tools and expand their applicability in high-lift system optimization.
Mesh generation is a fundamental step in the field of computational fluid dynamics and shape optimization techniques, particularly when dealing with search spaces that involve vastly different geometries. Two common types of meshing schemes are unstructured and structured. Unstructured meshing allows for fast variation of element size but lacks flexibility in handling domains with complex geometries. On the other hand, structured meshing can handle any geometry but may produce low-quality meshes and distorted finite elements for irregular contours. To address these challenges, D. López et al. [
6] propose a structured meshing procedure that creates subdomains to generate high-quality structured meshes in critical flow regions around aerodynamic sections. They demonstrate the effectiveness of this procedure by applying it to a straightforward shape optimization problem in fluid dynamics. The method involves dividing the mesh domain into blocks that vary according to the airfoil’s geometry. This approach ensures the generation of quality structured meshes for any geometry within the search space.
An alternative approach to address the meshing challenge in shape optimization with highly varied geometries is to employ mesh deformation techniques, particularly when working with gradient-based search methods where geometry changes are gradual. S. Jajobsson et al. [
34] introduced a method in their paper that utilizes an interpolant of boundary node displacements to extrapolate these displacements throughout the entire mesh domain, enabling deformation of both structured and unstructured meshes in 2D or 3D. The proposed method demonstrates its effectiveness by achieving deformed meshes with qualities comparable to the initial mesh in numerical examples. Furthermore, it offers the advantage of low computational cost and minimal memory requirements, making it suitable for situations that involve repeated deformations of the flow domain. By utilizing this approach, the mesh can be adapted and deformed to accommodate varying geometries efficiently and effectively in shape optimization processes.
When discussing aircraft design, the parameterization of geometry plays a key role. According to Daniel P. Raymer [
37], this design process can be divided into three phases: conceptual design, preliminary design, and detailed design. Wenbin Song et al. [
36] conducted a study on parameterization methods for airfoil shape optimization within a CAD-based design optimization framework. In their research, different parameterization methods were explored. The NACA airfoils, for instance, are defined by a series of y-coordinates at specified chord locations. Non-uniform Rational B-splines (NURBS) utilize control points to define the curvature of airfoils. B-splines and Bezier curves have been widely employed to fit airfoil shapes using interpolation techniques. Another parameterization method mentioned is PARSEC, developed by Sobieczky [
35], which uses geometric parameters such as leading edge radius, thickness-to-chord ratio, and maximum thickness to define the airfoil shape.
Scientists and engineers employ computer simulations to examine the relationships between input parameters and model outputs. However, conducting comprehensive parameter studies can be challenging or even impractical when simulations are computationally expensive and models have numerous inputs. In such cases, engineers seek ways to reduce the dimensionality of the input parameter space to enable more feasible studies. One promising approach to address this challenge is the use of active subspaces, which are a set of emerging dimension reduction tools. These tools help identify the most significant directions within the parameter space. By uncovering the active subspace, engineers can focus their analysis on the most influential factors, effectively reducing the dimensionality of the problem. In the book "Active Subspaces: Emerging Ideas for Dimension Reduction in Parameter Studies" by P. Constantine [
38], various techniques for discovering a model’s active subspace are discussed. The author also proposes methods for leveraging this reduced dimensionality to facilitate parameter studies that would otherwise be impractical or unfeasible. By utilizing active subspaces, engineers can gain valuable insights into the underlying structure of the parameter space, enabling them to streamline and optimize their simulations. This approach opens up new possibilities for efficient parameter studies and enhances the understanding of complex systems.
In their research paper, Zachary J. Grey et al. [
39] delve into the investigation of two commonly employed shape parameterizations: a) parametric sections (PARSEC) [
35] and b) class-shape transformation (CST) [
40]. Their primary focus is on estimating the active subspaces associated with airfoil performance metrics, specifically lift and drag. These metrics are computed using computational fluid dynamics considering various airfoil parameters. The study demonstrates the consistency of design-oriented insights derived from active subspaces across the two parameterizations. This consistency highlights the robustness and reliability of active subspaces as a dimension reduction technique. By reducing the dimensionality of the problem, active subspaces enable the use of simple visualization tools that provide valuable insights into the lift/drag Pareto front in multiobjective optimization for airfoil design.
Haipeng Li et al. [
5] present a very interesting paper that explores the application of an exergy-based method in aerodynamic designs. They investigate the impact of exergy on both two-dimensional airfoils and various three-dimensional wing planforms. To achieve the desired shapes of airfoils and wings, the researchers couple a numerical code with an exergy analysis, a shape-shifting algorithm, and an optimization code. The concept of exergy is rooted in the second law of thermodynamics, which states that exergy is not conserved during real processes. Instead, it is consumed to meet system objectives and some is irreversibly destroyed. The destruction of exergy corresponds to the generation of entropy. By minimizing entropy generation, one can better utilize exergy. Consequently, entropy generation minimization (EGM) methods have proven valuable in many thermodynamic optimization problems. In this paper, the authors combine entropy generation minimization with a RANS-based numerical flow solver, a shape-shifting algorithm, and an optimizer, resulting in a novel design methodology. This integrated approach enables the creation of optimized designs that maximize exergy utilization and minimize entropy generation..
3.1.3. Shape Optimization Approach
Optimization techniques have emerged as indispensable tools in the overall process of numerical analysis aimed at aerodynamic shape design. Due to the vast amount of literature available on this subject and limited space constraints, we can only highlight a few notable articles in this section that specifically focus on shape optimization processes. It is important to acknowledge that the mentioned papers represent only a fraction of the extensive body of literature on shape optimization. Nonetheless, they serve as valuable references for those seeking to explore this topic further.
In the past, optimization algorithms were primarily employed to minimize the difference between calculated and specified surface flow data, but this approach was not directly linked to optimizing the aerodynamic shape. In contrast, modern aerodynamic optimization methods integrate conventional analysis techniques with optimization algorithms, enabling iterative modifications to the geometry with the goal of minimizing a specific "cost" function, such as drag or the deviation from desired surface pressure distributions. While these aerodynamic optimization methods offer significant advantages, it is important to acknowledge that they can incur a high computational cost. The iterative nature of the optimization process, combined with the complexity of the underlying flow simulations, often demands substantial computational resources and time. By addressing the computational cost associated with aerodynamic optimization, researchers strive to make these methods more viable for real-world applications, enabling engineers to design more efficient and high-performance aerospace systems.
In 1992, G. S. Dulikravich conducted a study on Aerodynamic Shape Design [
2], which highlighted the potential of sensitivity analysis formulations as a promising approach for evaluation and optimization under off-design conditions. However, it was noted that sensitivity analysis can be computationally expensive, requiring a large number of analysis runs. To address this challenge, the concept of single cycle optimization was introduced and welcomed. This approach allows for a stable iterative algorithm, where an optimizer is applied to each updated configuration even before the flowfield has fully converged to the new geometry. This early integration of optimization within the iterative process shows promise in reducing computational expenses. Through a comparative analysis of optimization-based approaches, it was concluded that more economical methods can be achieved with gradient-based algorithms. These algorithms leverage the gradient information of the cost function to guide the optimization process efficiently. The availability of a wide variety of gradient-based optimization algorithms provides flexibility in choosing the most suitable one for a given problem. Additionally, it is worth noting that remarkable results were obtained using an evolution-type algorithm, which demonstrated reduced sensitivity to local minima. This alternative approach offers a different perspective by exploring the search space in a more comprehensive manner, potentially avoiding convergence to suboptimal solutions.
The Genetic Algorithm (GA) [
58], as an evolutionary computing technique, has gained significant popularity among researchers due to its effectiveness in finding global optima for complex functions. However, its computational time requirements pose a challenge in aerodynamic optimization problems that involve the use of computational fluid dynamics for fitness function computation. To address this limitation, researchers have explored various approaches to improve the optimum design and enhance the convergence rate of the optimization process, with a particular focus on proper population dispersion methods. Among the different dispersion methods, those incorporating mutation adaptation have proven to be particularly powerful. In a paper by M. Ebrahimi et al. [
44], an efficient method for scattering GA individuals in the design space is proposed to accelerate airfoil shape optimization. This method leverages the variation of the mutation rate for each gene of the chromosomes, guided by feedback from the current population. The proposed method introduces diversity into the population, allowing for a more thorough exploration of the design space.
In the optimization process involving genetic algorithms, fitness evaluations play a crucial role, particularly in shape optimization analyses of multi-element airfoils where viscous turbulent effects dominate the gap region of the control surfaces. To effectively utilize evolutionary algorithms, the integration of a mesh regeneration process into the design optimization cycle becomes necessary, as significant shape changes need to be incorporated. In their paper, Yu-Tai Lee et al. [
8] addressed this challenge by presenting a novel approach that combines a computational fluid dynamics process with a GA-based optimization procedure to evaluate the fitness function during the design/optimization iterations. This integration was considered essential due to the presence of strong viscous and turbulent effects in the gap region between the stabilizer and flap, particularly at large flap angles. By employing a multi-element unstructured Reynolds-Averaged Navier-Stokes (RANS) framework within the CFD process, the researchers were able to accurately capture and analyze the intricate flow phenomena associated with the multi-element airfoil configuration. The GA-based optimization procedure then utilized the insights gained from the CFD simulations to iteratively refine the airfoil shape and enhance its performance.
In a study related to the aforementioned topic, Ernesto Benini et al. [
10] presented a comprehensive approach for the design of multicomponent airfoils consisting of slat, main, and flap components. Their research focused on employing a multi-objective genetic algorithm in conjunction with an Euler/boundary-layer flow solver to drive the design procedure.
As discussed earlier in section 3.1.1, the accurate prediction of the behavior of high-lift multicomponent airfoil devices using numerical simulation poses a significant challenge for designers. The flow characteristics in these cases are complex, and transition phenomena occur within each aerodynamic element. To tackle these complexities, a multiobjective design approach for high-lift multicomponent airfoils becomes an attractive option for the optimization process. This approach allows designers to consider multiple conditions and requirements simultaneously, leading to more comprehensive and robust designs. In particular, direct multi-objective optimization proves to be valuable when dealing with conflicting requirements that need to be satisfied concurrently. By employing this method, designers can explore trade-offs between various objectives and identify design solutions that achieve a balance between competing criteria. D. Quagliarella et al. [
4] have contributed to this field by developing a multi-point design optimization method specifically designed for transonic and high-lift airfoil configurations. Their approach incorporates a multiobjective genetic algorithm that is coupled with two distinct flow field solvers. This integration enables the system to consider both transonic and high-lift requirements simultaneously through multiobjective optimization. By leveraging the strengths of this methodology, designers can explore a broader design space and identify optimal solutions that satisfy diverse performance criteria. The development of a multi-point design optimization method tailored to transonic and high-lift airfoil configurations is a significant advancement in the field.
In many design processes, especially in the context of shape optimization, understanding the underlying relationship between design variables and output variables can be a challenging task. In such cases, surrogate models provide a valuable tool to statistically relate input data to output data, enabling the approximation of complex, higher-order models in a simplified manner. Alfredo Arias-Montaño et al. [
12] presented a noteworthy study that applied a surrogate-based multiobjective evolutionary optimization approach to optimize airfoil aerodynamic design. The researchers employed multiple surrogate models, operating in parallel, with the objective of leveraging their respective strengths when tackling a computationally expensive multi-objective optimization problem. The team was able to create efficient approximations of the underlying aerodynamic models, reducing the computational burden associated with repeated evaluations.
Particle Swarm Optimization (PSO) [
59], although part of the emerging methods within evolutionary computation, differs from the previous ones in that each particle or individual in the swarm has memory and is capable of remembering the best position in the search space it has visited. Swarm algorithms have emerged as a powerful class of optimization methods capable of addressing both single and multi-objective optimization problems, with or without constraints. Unlike conventional evolutionary methods, swarm algorithms simulate social behavior, where each individual in a swarm adjusts its flight based on personal flight experience and the collective experience of its companions. The success of this strategy hinges on the effective exchange of information between individuals. In a paper presented by T. Ray and H. M. Tsai [
43], a swarm algorithm is applied to airfoil design optimization studies. The algorithm is grounded in a sociobehavioral model, providing designers with the desired flexibility to tackle various forms of the airfoil shape optimization problem, whether it involves constraints or not, and encompasses both single and multi-objective optimization. The study highlights that while the socialist behavior exhibited by the swarm algorithm may lead to an increased number of costly evaluations, it also offers the potential for further exploration and makes it a valuable tool for challenging airfoil design optimization problems.
While optimization methods based on evolutionary algorithms offer great potential, it is worth noting that for shape optimization processes involving small changes in the final geometry, simpler methods can often sufficient and even prove to be more robust than other optimization approaches. A noteworthy conclusion supporting this viewpoint is presented in the work by Matthias Schramm et al. [
42]. In their study, the researchers highlight that in scenarios where the modifications to the geometry are relatively minor, gradient-based optimization using the method of steepest descent can yield satisfactory results. This approach involves iteratively adjusting the design variables in the direction of the steepest descent, aiming to locate the minimum of the objective function. The advantage of the steepest descent method lies in its simplicity and computational efficiency compared to more complex optimization techniques.