Aerodynamic Design Optimization of a Morphing Leading Edge and Trailing Edge Airfoil–Application on the UAS ‐ S45

: This work presents an aerodynamic optimization method for a Droop Nose Leading Edge (DNLE) and Morphing Trailing Edge (MTE) of a UAS ‐ S45 root airfoil by using Bezier ‐ PARSEC parameterization. The method is performed using a hybrid optimization technique based on a Particle Swarm Optimization (PSO) algorithm combined with a Pattern Search algorithm. This is needed to provide an efficient exploitation of the potential configurations obtained by the PSO algorithm. The drag minimization and the endurance maximization were investigated for these configurations individually as two single ‐ objective optimization functions. The aerodynamic calculations in the optimization framework were performed using the XFOIL solver with flow transition estimation criteria, and these results were next validated with a Computational Fluid Dynamics solver using the Transition (cid:4666)𝛾(cid:3398)𝑅𝑒 (cid:3087) (cid:4667) Shear Stress Transport (SST) turbulence model. The optimization was conducted at different flight conditions. Both the DNLE and MTE optimized airfoils showed a significant improvement in the overall aerodynamic performance, and MTE airfoils increased the efficiency of C L3/2 /C D by 10.25%, indicating better endurance performance. Therefore, both DNLE and MTE configurations show promising results in enhancing the aerodynamic efficiency of the UAS ‐ S45 airfoil.


Outline of the Research
The goal of reducing global fuel consumption and fuel-related emissions has placed tremendous pressure on the aviation industry. In 2018, the aviation industry was responsible for 895 million tons of CO2 emissions globally emitted into the atmosphere [1]. Reducing fuel consumption will benefit both the world environment and the air transport industry. According to the Aviation Transport Action Group (ATAG), a reduction in fuel burning will have a significant impact on the aviation industry because the largest operating cost of this industry is fuel [1]. Drastic measures are still required even though some steps have already been taken to reduce these emissions. Various solutions have been implemented by the aviation industry, including smart material technology, laminar flow technology, air traffic management technologies, advanced propulsion techniques and sustainable fuels [2][3][4]. Morphing wing technology is one of the technologies showing high potential in decreasing aircraft fuel consumption [5][6][7][8].
Even though there is no settled definition for "morphing", this term is borrowed in Aviation Technology from avian flight to describe the ability to modify maneuvers at certain flight characteristics in order to obtain the best possible performance. This type of application requires morphing structures capable of adapting to changing flight conditions. Morphing systems include several wing shape modifications, span or sweep changes, changes in twist or dihedral angle [9] and variations of the airfoil camber [7] or of the thickness distribution [4]. In the design phase, both military and civilian aircraft traditionally fly at a single or few optimum flight conditions; morphing, however, is envisioned to increase the number of optimum operational points for a given aircraft [5].
The morphing strategy scope is broad in the sense that an optimal solution can be generated in various ways with respect to a wide spectrum of mission profiles and types, and flight regimes. Novel technologies need to be part of recent efforts to produce environmentally sustainable and efficient aircraft. To meet environmental guidelines, including those of the International Civil Aviation Organization's (ICAO's) Committee on Aviation Environmental Protection (CAEP), innovative technologies such as "morphing" are fascinating as they provide advantages over conventional wing configurations [10]. Additionally, "morphing" is more practically applied in Unmanned Aerial Vehicles (UAVs) because of their reduced scale and lower complexity in terms of wing design structure and energy consumption expressed in terms of actuation power [6,11,12], in addition to the advantages that morphing designs are lighter and less noisy than their conventional designs. Some morphing opportunities with potential benefits for increasing aerodynamic performance of the UAS-S45 are shown in Figure 1. Aircraft are designed to achieve optimal aerodynamic performance (maximum lift-to-drag ratio) for a mission specification and for an extended range of flight conditions. Nevertheless, these mission specifications change continuously throughout the different flight phases, and an aircraft often flies at non-optimal flight conditions. Although the conventional hinged high lifting devices and trailing edge surfaces (discrete control surfaces) in aircraft are effective in controlling the airflow for different flight conditions, they create surface discontinuities, which increase the drag [13,14]. The disadvantages of these hinged surfaces are found in both their deployed and retracted configurations [15]. When deployed, the gaps between the high lifting surface and the wing can cause noise and turbulence, and therefore can produce a turbulent boundary layer and thus increase drag. Even when retracted, the trailing edge hinges still produce a turbulent boundary layer. Numerous researchers believe that the laminar flow technology has the maximum potential to reduce drag and to avoid flow separation [16,17]. For this goal, the wings require the design of thin airfoils, seamless high-quality surfaces and variable droop leading edges. The use of these high continuous smooth lifting devices in UAVs will further broaden their flight envelope and extend their endurance. These challenges can be addressed by using morphing technology in the current aviation industry.
The aim of this study is thus to investigate the aerodynamic optimization by employing constrained shape parameterization for UAS-S45 root airfoil using Droop Nose Leading Edge (DNLE) and Morphing Trailing Edge (MTE) technologies. The aerodynamic design of continuous morphing wing surfaces employed to replace traditional discrete control surfaces, such as flaps, ailerons or sometimes slats, to adjust the camber of the wing, are discussed in this paper along with the parameterization strategy and the optimization algorithm.

Literature Survey
The potential in the field of aircraft morphing is evident as hundreds of research groups worldwide have been dedicated to studying the different aspects of morphing. Two of the earliest groups were launched by NASA, and they worked within the "Active Flexible Wing program" and the "Mission Adaptive Wing program" [18,19]. In these programs, flutter suppression, load alleviation and load control during rapid roll maneuvers were investigated.
One of the potential consequences of flexible morphing structures is that both dynamic properties and aerodynamic loads of the wing are affected. Hence, the aeroelastic effects arising from interactions between the configuration-varying aerodynamics and the morphing structure are significant. Li et al. have presented important morphing studies in terms of aeroelastic, control and optimization aspects [20].
The nonlinear aeroelastic behaviour of a composite wing with morphing trailing-edge has been studied [21]. Moreover, the nonlinear aeroelastic behaviour of a composite wing with morphing trailing-edge has been studied.
Aero-elastically stable configurations of a morphing wing trailing edge driven by electromechanical actuators have been investigated [22]. The effects induced by trailing-edge actuators stiffnesses on the aeroelastic behaviour of the wing were simulated using different approaches. The results revealed that flutter could be avoided if sufficient stiffness was provided by these actuators.
An analytical sensitivity calculation platform for flexible wings has been studied [23,24]. This platform has been used to perform wing aeroelastic optimization and stability analysis.
In another study, a low-fidelity model of an active camber aero-elastic morphing wing was developed to investigate the critical speed values by varying its chord-wise dimension [25]. A wide range of configurations was explored to predict the dynamic behaviour of these active camber morphing wings.
The combined flutter behaviour and gust response of a series of flexible airfoils has been investigated [26]. Rayleigh's beam equation was used to model this series of flexible airfoils. The effect of chordwise flexibility of a compliant airfoil was investigated numerically to demonstrate its dynamical stability [27,28]. An actuated two-dimensional membrane airfoil was investigated experimentally and numerically, and it was found that membrane flexibility decreased the drag and delayed the stall.
These aeroelastic studies are significant in understanding the flexible morphing wing and its potentially relevant effects. However, our study does not consider the structural effects of morphing, and therefore the optimization is aimed to improve the aerodynamic drag performance and endurance in the given range of angles of attack.
Other research projects in the US and Canada include the "Smart Materials and Structures Demonstration" program initiated to develop new affordable smart materials to establish the performance gains by investigating smart rotors, the smart aircraft and the marine propulsion system; the "Aircraft Morphing" program [29]; the "Advanced Fighter Technology Integration" (AFTI) program, aimed to develop and demonstrate in flight a smooth variable camber wing and flight control system capable of adjusting the wing's shape in response to flight conditions to maximize aerodynamic efficiency; the "Active Aeroelastic Wing" program, demonstrated to improve aircraft roll control through aerodynamically induced wing twist on a full-scale high-performance aircraft at transonic and supersonic speeds; the "Morphing Aircraft Structures" program [30]; the "Mission Optimized Smart Structures" (MOSS), in which aerodynamic optimization was performed with the stretching of the leading edge using a novel skin material developed with a nanocomposite at the National Research Council of Canada [31]; and the "Controller Design and Validation for Laminar Flow Improvement on a Morphing Research Wing-Validation of Numerical Studies with Wind Tunnel Tests-CRIAQ MDO 7.1" [32,33]. The CRIAQ 7.1 project took place at our Laboratory of Applied Research in Active Controls, Avionics and AeroServoElasticity LARCASE.
The European Union has conducted several projects, including the "Active Aeroelastic Aircraft Structures" (3AS) project, to develop novel active aeroelastic control strategies to improve aircraft performance (structural weight, better control effectiveness) by controlling structural deformations to modulate the desired aerodynamic deformations [34]; the "Aircraft Wing Advanced Technology Operations" (AWIATOR) project, in which novel fixed wing configurations were introduced aimed at reducing the vortex hazard by implementing larger winglets and further improving aircraft efficiency and reducing farfield impact [35]; the "New Aircraft Concepts Research" (NACRE) project, in which Powered Tails and Advanced Wings were studied to obtain high environmental performance (noise and CO2 emissions)-both high aspect ratio low-sweep wings and forward-swept wings (with natural laminar flow) contributed to achieving good fuel efficiency [36]; the "Smart Leading Edge Device" (SmartLED) project [13]; the "Smart High Lift Devices for Next Generation Wings" (SADE) project [37]. Both SmartLED and SADE were aimed to develop and investigate the morphing high lift devices, "smart leading edge" to enable seamless high lift devices and therefore enable laminar wings and "smart single-slotted flap", for the next generation aircraft of high surface quality for drag reduction. Projects by the European Union also included the "Smart Fixed Wing Aircraft" (SFWA) project [38] and the "Smart Intelligent Aircraft Structures" (SARISTU) project to address the physical integration of smart intelligent structural concepts. The SARISTU included a series of research collaborations to addresses aircraft weight and operational cost reductions as well as an improvement in the flight profile specific aerodynamic performance [39]. Other projects by the European Union were the "Novel Air Vehicle Configurations" (NOVEMOR) project, in which morphing wing solutions (span and camber strategies and wing-tip devices) were proposed to enhance lift capabilities and maneuvering [40]; the "Clean-Sky 1 and 2" projects [41,42]; the "Combined Morphing Assessment Software Using Flight Envelope Data" (CHANGE) project, which developed a modular software architecture capable of determining and achieving optimum wing shape [43]; and the "Sustainable and Energy Efficient Aviation" (SE2A) project aimed to investigate the Morphing structures for the 1g-wing to exploit the nonlinear structural behavior of wing design components to achieve passive load alleviation [44].
The "Morphing Architectures and Related Technologies for Wing Efficiency Improvement-CRIAQ MDO 505" was also realized at the LARCASE in the continuation of the CRIAQ 7.1 project mentioned above. The achievements of the international Canadian-Italian CRIAQ MODO 505 project are mentioned in various publications [45,46].
The aircraft optimization process has evolved dramatically over the past decades. The design of new intelligent algorithms and computational solvers has substantially impacted the overall design process, including the "morphing aircraft optimization". The verification and validation of aircraft design using optimization techniques has reduced its huge experimental costs and has been achieved by the use of efficient algorithms and computational solvers [47,48]. An optimization process is initiated by minimizing or maximizing an objective function for the target design concept with respect to design variables subjected to given constraints, which were applied to limit the search space and therefore to yield physically feasible optimization results. Both gradient-based and gradient-free intelligent algorithms have been used in these optimization processes, and their results are dependent on the type of optimization problem to be solved [47]. These algorithms are discussed in the following sections. Many optimization techniques based on bioinspired processes and Surrogate-Assisted natural methods have been formulated, and mathematical and statistical analysis of these algorithms has been performed [49][50][51][52]. The comparative performance study based on convergence trends was performed with different geometry parameterization techniques. Both advantages and disadvantages can be found in these optimization algorithms for complex shape functions. However, the goal of this study is to perform aerodynamic design optimization of morphing airfoil using Particle Swarm Optimization (PSO) algorithm combined with Pattern Search technique.
Similarly, various computational solvers are employed based on their desired solution accuracy and computational cost [53]. Researchers have conducted various aerodynamic and structural optimization studies, including design studies of different surfaces (components) of the wing, such as its upper surface, trailing edge or leading edge. Adjoint-based airfoil optimization has been studied using Euler equations [54], while other researchers have implemented extensive gradient-based aerodynamic shape optimization methodologies for transonic wing design using high-fidelity RANS solvers [55]. Many other optimization studies have been performed using different optimization techniques. Some of the most notable works are summarized in Table 1. The study presented in this paper analyses the aerodynamic optimization of the Hydra Technologies S45 Balaam UAS airfoil. This unmanned aerial surveillance (UAS) system is needed to provide security and surveillance capabilities for the Mexican Air Forces, as well as civilian protection in dangerous situations. This study of Morphing Wing Technology has the aim to enhance the aerodynamic efficiency and the effective range of the UAS S45, as well as to extend the flight time. The next section of this paper presents the optimization framework of the overall methodology, which includes two objective functions' formulations ("drag minimization" and "endurance maximization"). The parameterization technique is also presented in Section 2 and has the aim to obtain the "optimal aerodynamic shapes"; the computational solvers for calculating aerodynamic coefficients and the optimization algorithms employed are also presented. The results obtained by these solvers and algorithms for the UAS-S45 optimized designs are compared with the baseline UAS-S45 model results.

Method
This optimization framework study performs the aerodynamic optimization of the leading edge and trailing edge of a morphing airfoil by using an intelligent and iterative process based on user-defined aerodynamics and constraints. The methodology involves the objective function formulation integrated with a geometrical shape parameterization model, an aerodynamic flow solver and the optimization algorithm.
For this study, an optimization framework is designed to allow the integration of shape generation using the direct manipulation of the airfoil shape variables while respecting the geometrical constraints. The optimization procedure framework is described in Figure 2, and it consists of a parameterization block based on the Bezier-Parsec (BP) parameterization airfoil shape technique and an aerodynamic solver, such as the panel solver XFOIL (Version 6.99, MIT, USA, 2013 ). The validation of its results is done using high-fidelity solver in Fluent based on Reynolds-averaged Navier-Stokes (RANS) equations in the Transition SST model. A hybrid optimizer based on Particle Swarm Optimization (PSO) algorithm combined with the Pattern Search technique is also used. The optimization procedure is shown in Figure 2 and is described in detail in Sections 2.1-2.4.

Objective Function Definition
Numerous researchers have stated that there is always a trade-off between aerodynamic coefficients and performance parameters in aerodynamic design problems [66]. A difficult balance exists between increasing aircraft performance and meeting design constraints efficiently in practical aerodynamic design problems. For example, the objective of aerodynamic optimization is lift maximization, drag minimization or maximization of the lift-to-drag ratio while keeping CL and/or CD constant as constraints. Likewise, various airfoil shapes are optimized for their morphing by considering an increased lift-to-drag ratio and, at the same time, for meeting the given constraints. Other airfoil shape optimizations may involve those of multiple aerodynamic coefficients based on performance parameters such as range, loiter and control. Such optimization techniques can be applied to solve complex aerodynamic design problems for a broad range of flight conditions.
In addition, the design problem shown in this paper must consider the extended range of flight operating conditions that a UAV is expected to encounter within its flight envelope. As mentioned above, this study considers the optimization of the airfoil shape for drag minimization and for endurance maximization as individual target objective functions.
The mathematical formulation of the drag minimization objective function is the following. ; (1) Su The lift coefficients of the UAS-S45 are calculated for each angle of attack using design parameters mentioned in a previous paper of our team [67]. The constraints implemented are both maximum lift coefficient and lift coefficient of baseline airfoil at each angle of attack. The optimization runs are executed with a minimum lift coefficient at each angle of attack. Constraints are applied by penalty functions to the objective function to ensure optimum performance is obtained at each given angle of attack.
For the design of a long-endurance UAV, Equations (2) and (3), which are needed to determine the endurance (E) of a propeller-powered UAV, are well established using the Breguet formula [68].
where is the specific fuel consumption, is the propeller efficiency, is the free stream density, is the wing planform area, is the initial weight of the aircraft, is the final weight of the aircraft and and are the lift and drag coefficients, respectively.
The parameters that do not affect the airfoil shape design in the endurance performance can be neglected in Equation (2), and the parameters that can affect the airfoil shape design are only the two aerodynamic coefficients and ; therefore, Equation (2) can be reformulated as (3) Therefore, an optimized airfoil maximizes the lift and minimizes the drag coefficients in the above Equation (3); therefore, the overall endurance aerodynamic efficiency is maximized. The weighted average of the flight envelope is used to find out the time that the UAV would have to spend in different flight regimes. Since the present study concerns the surveillance UAS-S45 optimization, it is assumed that 80% of total flight time is spent in loiter flight conditions, which require maximum endurance. Therefore, the objective function to maximize the endurance aerodynamic efficiency with the same constrained optimization is given by the following equation: where ∈ , ℎ ℎ

Parameterization Strategy
Parameterization of the airfoil shapes for the UAS design and optimization requires a mathematical formulation, and it becomes an important part of the overall optimization process. Studies have shown that the choice of shape parameterization technique has a strong influence on the solution accuracy, robustness and computational time of the overall optimization process [69]. In order to obtain an optimal aerodynamic solution, the shape function must be directly related to the airfoil geometry, have a flexible design space, and be robust for the control of all parameters. Some of the well-known methods are the Discrete Based Approach, the Bezier curves, the B-Spline or NURBS curves, the Cubic Spline and the Free-Form Representation [69,70]. The drawbacks of these methods are that they do not use airfoil shape parameters, they require a large number of design variables and they often provide inaccurate shapes for the leading edge and trailing edge of an airfoil.
These limitations have been partly removed by using the polynomial-based functions, such as PARSEC and CST. "PARSEC" is characterized by eleven design coefficients that control airfoil shape parameters, and they are the: (a) edge radius, (b) upper crest abscissa, (c) upper crest ordinate, (d) upper crest curvature, (e) lower crest abscissa, (f) lower crest ordinate, (g) lower crest curvature, (h) trailing edge ordinate, (i) trailing edge thickness, (j) trailing edge direction and (k) trailing edge wedge angle. Several studies have shown PARSEC method's superior performance with respect to the Discrete Based Approach, based on their convergence rate, flexibility and epistasis, which is due to the nonlinear dependency of the objective function on the design parameters [71]. However, it has also been found that the leading edge and trailing edge design coefficients in this parameterization method do not provide accurate airfoil shape; therefore, implementing morphing design using only these PARSEC design coefficients might further increase the airfoil shape complexity.
Therefore, a combined approach, the Bezier-PARSEC (BP) parameterization method, was developed [72] that included the advantages of Bezier with PARSEC combination, and it was found to be suitable for this study. The BP parameterization was implemented in this paper with the aim to increase the solution accuracy, flexibility and efficiency of morphing leading edge and trailing edge optimizations. More significantly, this technique provides the large search space needed for morphing leading edge and trailing edge design, and it reduces the computational time.
The BP method has been further divided into two parameterization sub-methods: BP3333 and BP3434. Four third-order Bezier curves are used to represent the airfoil shape in the BP3333 parameterization sub-method. The BP3434 parameterization has a 3rd-degree edge thickness curve, a 4th-degree trailing edge thickness curve, a 3rd-degree leading edge camber curve and a 4th-degree trailing edge camber curve. Due to the smaller number of degrees of freedom in the BP3333 method with respect to those of the BP3434 method, specifically at the trailing edge and at the leading edge, the BP3434 method was chosen in this research.
A fourth-degree Bezier curve is given by: and The parameterization is controlled by 15 parameters: 10 aerodynamic and 5 Bezier parameters. Figure 3 shows the graphical representations of these parameters. The control points definitions for the leading and trailing edge Bezier curves are given in Table 2 for the thickness profile curve, and in Table 3 for the camber profile curve [73]. The upper and lower bound values of the parameters used in the BP parameterization optimization method calculated in this study for the UAS-S45 are presented in Table 4.

Description of Flow Solvers
Various solvers have been employed in aerodynamic optimization studies, such as the Fully Potential Flow, Coupled Boundary-Layer, Euler and Viscous Navier-Stokes solvers [74]. The type of solver chosen by a user depends upon the type of optimization problem. Results obtained with the flow solver must be consistently accurate in order for the optimization to be considered highly adequate.
For this aerodynamic research, incompressible flow solvers are needed as they can perform fast computation, and also can provide high accuracy results. Therefore, in this paper, two types of aerodynamic solvers are used for the aerodynamic optimization: a Panel-based method and a Reynolds-Averaged Navier-Stokes equations (RANS) solver. The choice to keep the computational cost low for an airfoil analysis leads directly to the use of the XFOIL code, which is known for its very good combination of execution speed and accuracy of results. Furthermore, in order to obtain higher accuracy results for solving the viscous boundary layer and the flow separation than the XFoil solver, a well-established high-fidelity Computational Fluid Dynamics (CFD) solver, Ansys Fluent was also used in this morphing airfoil analysis.

XFoil
XFoil is a code for airfoil design and analysis consisting of inviscid, inverse and viscous formulations [75]. The viscous formulation can be used to calculate free or imposed forced flow transition to handle transitional separation bubbles, to calculate aerodynamic coefficients and to cope with moderately trailing edge separation. It uses an approximate e N envelope method to calculate the flow transition. The turbulence level settings are kept as default with free transition features. Both the boundary layer and the wake parameters are calculated with a two-equation integral boundary layer formulation as shown in Equations (8) and (9), respectively [76].
where shape factors * and * * are defined in the following equations * At high angles of attack where stall occurs, XFoil has difficulties in giving a converging solution for the airfoil analysis. To obtain the convergence of the solution, non-converged solutions are given a high penalty by setting the fitness function arbitrarily large via a penalty function, and thus they are eventually eliminated during the design process. The solver tolerance value used for this penalty function calculation is 0.002.
Many researchers have employed XFoil software in morphing airfoil studies. In one of these studies, the results obtained using XFoil were compared with those obtained using the Reynolds-Averaged Navier-Stokes (RANS) solvers for camber morphing, and only a slight difference between the results was observed [77].

CFD Fluent Solver
The validation of solvers and mesh settings is needed before carrying out the optimization studies. The validation studies were performed at a Reynolds number of 2.4 × 10 6 and a Mach number of 0.10. The range of angles of attack was considered from 0° to 16° to include the stall angle. The aerodynamic coefficients and the pressure distribution were obtained and further used for solution validation. A C-shaped computational domain with structured grids, as seen in Figure 4, was used in domain discretization after the performance of the mesh convergence test. The turbulence model used in this study was the Transition SST. This model uses a combination of SST K-ω coupled with intermittency γ and transition onset Reynolds number.
is the critical Reynolds number where the intermittency starts [78]. Four transport equations of the Transition SST model are given below: The spatial discretization used in the solution is of the second order for all the turbulence statistical parameters, including the momentum, the turbulent kinetic energy, the specific dissipation rate, the intermittency and the Reynolds momentum thickness.
ANSYS Fluent was used in this study to carry out the numerical analysis. The meshing was done using ANSYS Mesher, as shown in Figure 4. The rectangular computational domain was designed around the airfoil geometries. The outlet length was selected to be 30c, which was adequate to allow for the full development of the wake flow, based on the size of the computational domain adopted in previous studies. The distance from the domain inlet to the airfoil was 10c to prevent the inlet boundary from unphysically impacting the induction field upstream of the airfoil. The high-quality grids were generated with dense grids in the boundary layer of the airfoil with a gradual decrease of grid cells away from the airfoil surface. The first layer's thickness was calculated based on inlet velocity and adopted accordingly to capture the laminar and transitional boundary layer regions. Furthermore, the grid sensitivity analysis was performed according to the number of cells in the computational domain.
A freestream velocity of 34 m/s with a turbulence rate of 0.01 percent was imposed at the inlet. The surface of the airfoil was modelled as a zero-roughness no-slip wall, and its values were executed on the grid domain, which prescribed its motion relative to the rest of the computational domain. A symmetry boundary condition was applied to achieve a parallel flow at the top and bottom of the domain by assuming zero normal velocity and zero normal gradients of the flow quantities. A zero static gauge pressure was applied at the exit of the domain. Iterations were finalized when all scaled residuals were below 1 × 10 6 .

Optimization Algorithm
An intelligent search algorithm is essential for the direct numerical optimization of airfoil design. This algorithm operates iteratively and utilizes the inputs of the shape parameterization method to define the airfoil shape, and then it uses the flow solver to calculate the aerodynamic coefficients. Its efficiency is assessed according to its ability to provide a global optimal solution with reasonable computational resources. Therefore, the choice of such an algorithm influences the solution convergence and its feasibility.
Researchers have implemented various algorithms, such as the Gradient-Based, Adjoint-Based and Evolutionary algorithms with the aim to investigate and solve various optimization problems [74]. These algorithms have been considered to solve different aerodynamic optimization problems based on their advantages and disadvantages. The hybrid optimizer based on the Particle Swarm Optimization (PSO) algorithm was used in this paper for the UAS-S45 optimization, and it was further combined with the Pattern Search algorithm in order to enhance the solution convergence and its refinement.
Evolutionary algorithms have attracted the attention of shape designers due to their characteristics, such as their suitability for modelling discontinuous shape functions, obtaining global optimal solutions, ease of parallel computing, etc. During morphing airfoil optimization, gradient-based algorithms converge fast, but they only cover a small area out of a large airfoil search space, and thus a local minimum solution instead of global minimum solution is found. It is also difficult to find the gradient of the non-linear flow fields for gradient-based optimizers. Therefore, a PSO algorithm is used, based on a simplified social behavior closely related to the swarming theory, where the solutions are represented by a set of particles that heuristically navigate through a design space [79]. The efficiency of PSO algorithms over genetic algorithms is due to their independence from parameters, such as crossovers and mutations; instead, the solution is updated by sharing its information amongst its populations of particles.

PSO Algorithm
In the PSO algorithm, each particle is a solution to a given optimization problem and is composed of two vectors: "position" and "velocity". A position vector is used to store the positioning of the particles in the given dimensional space. The velocity vector is updated using the following equation for each iteration k: where = position vector; = randomly generated velocity vector of the particles at search initialization; = vector representing the best solution achieved by the particle; = vector representing the "global best solution" collectively achieved by the swarm; and = stochastic acceleration terms that pull each particle toward and positions respectively; and = random numbers in the uniform range [0, 1]; The update of the position vector at the iteration is given by the following equation: 1 The PSO procedure can be outlined as follows: 1. Generate an initial swarm with arbitrary values of the particle position and arbitrary initial velocities in the n-dimensional search space; 2. Evaluate the fitness of each particle of the swarm by searching a new velocity vector applied to each individual particle; the new velocity vector is influenced by its best position, the swarm's best position and its previous velocity; 3. Update the velocity of the particles as shown in Equation (14) and the position of the particles as shown in Equation (15) and evaluate the new fitness function; 4. Repeat steps 2 and 3 until the maximum number of iterations is reached or until the results no longer improve for a given number of iterations.
Pattern Search (PS) is a direct search optimization method that is used to minimize a function by comparing its value at each iteration with its value at its previous iteration in a finite set of trial points. The method is outlined as follows: 1. Select the initial solution; 2. Explore the direction search by evaluating each objective function value, one at a time; 3. Evaluate the objective function. If the objective function is minimized, then update it; otherwise, do not update it; 4. Evaluate the next values of the objective function following the previous steps 2 to 3; 5. After exploring the search of all possible values of the objective function and the minimization of the overall objective function, the PS method follows a pattern move; and 6. If the overall objective function has not been minimized, repeat the process.

Hybrid Optimization Scheme
A hybrid optimizer is used in this study, as shown in Figure 5. In this figure, the hybrid optimization implies that the PSO-based optimization process is applied firstly to its variables to obtain their optimized values. Since the first optimization process may not be accurate as it would result in low variations of values to be optimized, a second optimization process would then be applied to the results of the first optimization process. In this case, PSO, a non-elitist algorithm, is then followed by Pattern Search, using hybrid optimization. The PSO and Pattern Search can then be used to find a local optimum in a limited search area; the use of both PSO and PS algorithms through hybrid optimization ensures better results than a single optimization process. The independent variables given as inputs to the optimization process are also limited due to their upper and lower bounds. These bounds ensure that the search space for both the PSO and Pattern Search algorithms remains within them. During the PSO algorithm execution, if a solution is found outside the bounds of the search space, any such solution outside these bounds will have its value replaced by that of its given upper and/or lower bounds, without any influence on other solutions. The hybrid optimizer's procedure is shown in Figure 5.

Discussion and Results
The aerodynamic optimization considered in this study applies to both Droop Nose Leading Edge Morphing (DNLE) and Morphing Trailing Edge (MTE) configurations. To utilize the full potential of the "morphing" concept, a baseline shape optimization was carried out followed by its morphing shape design at various flight conditions given in Table 5. These cases are based on various altitudes and Reynolds numbers that are computed for the airspeed of 34 m/s. These given flight conditions are chosen because of the fact that the UAS-S45 can reach altitudes of 20,000 ft and a stall speed of 34 m/s. The design framework uses endurance maximization as an optimization function to enhance the performance of the UAS-S45 over a large part of its flight regimes. Figure 6 shows the flight envelope of UAS-S45. From a physical standpoint, every flight condition is defined by a specific Mach number, altitude (Reynolds number and temperature are calculated according to ISO standard atmosphere) and a fixed lift coefficient. The present study considered the optimization of the leading edge and trailing edge of the UAS S45 conventional wing that could be both morphed. The difference of performances obtained between the baseline airfoil and the morphed one demonstrate the overall benefits of morphing airfoil design. Figure 7 Figure 8 shows the aerodynamic coefficients CL and CD variations of the baseline UAS-S45 root airfoil for two Reynolds numbers Re = 1.2 × 10 6 and Re = 2.4 × 10 6 , which were calculated using the XFoil and the Transition SST in Ansys/Fluent software. As seen in Figure 8, the lift coefficients' variations with the angle of attack were predicted accurately, while small differences were found in drag coefficient variations with the angle of attack. These differences in the drag coefficients' variations are found largely at angles of attack higher than 4° because the pressure-based drag coefficient becomes the dominant component of the total drag coefficient and flow separation conditions result in a decreased friction-based drag coefficient. This result implies that drag forces are very sensitive to turbulence effects and the Transition SST model can solve the turbulence problems more accurately than the XFoil solver. Figure 9a shows the friction coefficient variations with the chord for the angles of attack of 0°, 2° and 6° obtained by use of the Transition SST model; this figure indicates very good accuracy in predicting the flow separation phenomena as the skin friction coefficient increases slightly with the increase in angle of attack. Figure 9a,b also shows the initial flow separation, the separated region and its reattachment. It is seen in Figure 9a that the laminar-to-turbulent transition region is located at 18% of the chord (0.18c) at an angle of attack of 6° and at 43% of the chord (0.43c) at an angle of attack of 2°.The variations of the skin friction coefficient along the chord were assessed at different angles of attack using the Transition SST model in Figure 9. The velocity magnitude contours plots at three angles of attack (4°, 10°, 16°) at Re = 2.4 × 10 6 are shown in Figure 10. The turbulent intensities play an important role in understanding the flow behavior over an airfoil, and they are shown in Figure 10 for the angles of attack of 4° and 10°.
The contour lines remain attached to the airfoil at low angles of attack such as 4°, as shown in Figure 10a; the flow starts to separate at the trailing edge of the airfoil at an angle of attack of 10° as shown in Figure 10b and changes into a fully separated flow at the angle of attack of 16° as shown in Figure 10c. It is important to understand that boundary layer separation takes place at an angle of attack of approximately 10°. The separation region shown in Figure 10c causes increased drag as it induces a large wake that completely changes the flow downstream of the point of separation. It is clear that the turbulent intensity is minimal at the airfoil leading edge at a low angle of attack of 4°, as seen in Figure 11a, while at the angle of attack of 10°, the turbulent intensity increases over the airfoil surface, including at its leading edge.

Optimization of Morphing Leading Edge
The Droop Nose Leading Edge (DNLE) optimization was performed to satisfy the objective functions defined in Equation (1) for drag minimization and in Equation (5) for endurance maximization. The optimization processes were carried out for two cases of DNLE control point deflections, defined as "Opt. Case I" and "Opt. Case II".
In "Opt. Case I", a morphing DNLE is a smooth single structural layout with a morphing leading edge starting location fixed at 22% of the chord. The trailing-edge is not allowed to move, while the leading edge points are only allowed to move in the vertical direction to effectively adapt the chord-length variation of the baseline airfoil.
In "Opt. Case II", the morphing DNLE is free to change its shape anywhere from 0 to 22% of the chord. The starting point and the end point of the leading edge are used to represent the morphing DNLE shape. The parameterization of the DNLE shape was obtained by using the Bezier-PARSEC (BP) method, in which the leading edge radius of the airfoil and 5 Bezier control points were used to maintain the slope continuity within the airfoil shape.
The XFoil software was used to efficiently investigate the aerodynamics of the optimized morphing geometries. The aerodynamic performance of the DNLE was evaluated to determine its best airfoil design shape for each flight condition. The computations were performed for the three flight conditions listed in Table 5. The results were obtained in terms of aerodynamic lift and drag coefficients variations with angle of attack and pressure distributions versus the airfoil chord locations for the optimized airfoil geometries. Figure 12a illustrates the initial (reference) versus the final (optimized or morphed) airfoil shape, while Figure 12b shows the convergence trend for the minimum drag optimization function for an optimized airfoil shape with DNLE at an angle of attack of 50. The algorithms PSO with Pattern Search converge to the optimal solution after 90 iterations. The x-axis indicates the number of function evaluations, and the y-axis shows the value of the fitness. The fitness value is the objective function value of each particle of the algorithm, expressing the performance of any airfoil shapes. These results for the droop nose leading edge (DNLE) optimization are shown in Figures 13 and 14, and the UAS-S45 baseline airfoil aerodynamic coefficients were compared to those of the optimized airfoils. Figure 13a Figure 13a-c show that higher lift coefficients were obtained for Opt. Cases I and II than those of the baseline airfoil. Figure 13d shows that the DNLE designs for Opt. Case I at three different flight conditions. It is found that Opt. Case I DNLE design is the best at Flight Condition III. These optimal shapes revealed that maximum increment of the lift coefficient with respect to the baseline airfoil of up to almost 21% was found for Opt. Case I design seen in Figure 13c. In addition, an increase of 9.6% was obtained for the maximum lift coefficient, associated with a stall angle delay of 3° as seen in Figure 13c. Therefore, it was found that both Opt. Case I and Opt. Case II results are better than the baseline airfoil results. These better results can be explained by the pressure distribution change around the airfoil. The DNLE design does not cause an adverse pressure gradient as the flap configuration, but instead, it decreases the maximum speed at the airfoil leading edge, which is similar to the findings [80]. Figure 14 presents the evaluation of the lift vs. drag polar of the baseline airfoil and of the optimized airfoils for three different flight conditions. Therefore, the increased aerodynamic performance of optimized airfoils in terms of their lift versus drag variation for all three flight conditions can be seen in Figure 14a The pressure coefficients' variations with the chord of optimized airfoils for three flight conditions are shown in Figure 15. It is clear that the major changes in pressures take place near the leading edge at its upper surface, and that the pressure peaks then remain smooth on the rest of the airfoil for all configurations in Figure 15a for an angle of attack of 4° and in Figure 15b for an angle of attack of 10°. The chord-wise pressure distribution reveals that increased performance was obtained for DNLE configurations. By comparing the chord-wise pressure distributions for different configurations at angles of attack of 10° for the baseline versus optimized airfoils, it is found that a kink in the suction peak is followed by a constant pressure distribution. The DNLE-optimized configuration at flight condition III undergoes flow separation, and reattaches shortly on the airfoil surface, as seen in Figure 15b.   Figure 17 shows that the optimization process led to an increase in endurance maximization characterized by CL 3/2 /CD in the optimized DNLE airfoils for three different flight conditions. The maximum endurance capability given by CL 3/2 /CD for the same lift coefficient of the baseline airfoil to the lift coefficient of the morphing DNLE configuration Opt. Case I was found to be 1.17 and 1.21, respectively. Furthermore, the values of CL 3/2 /CD increased from 117 to 132, thus indicating better endurance performance for the UAS-S45 DNLE airfoil configurations than for its baseline airfoil at three flight conditions.
The flow transition behavior can be seen in the velocity contour plots of one of the DNLE airfoils, as shown in Figure 18a for an angle of attack of 4° and in Figure 18b for an angle of attack of 10°. The comparison of baseline airfoil and DNLE configurations at a 10° angle of attack shows that higher gradients are formed for the DNLE airfoil in comparison with those of the baseline airfoil, illustrated in Figure 10. However, the trailing edge flow separation region remains the same for both airfoils. In addition, the Turbulent Kinetic Energy (TKE) for the baseline airfoil and for the DNLE airfoils by use of the Transition SST turbulence model are shown in Figure 18c,d. It is evident that the leading edge is propagating the energy towards the downstream; thus a considerable decrease is observed in the wake turbulence. The TKE occurring in the baseline airfoil, which starts from the upper surface and moves towards the trailing edge was not observed in the DNLE-optimized airfoil at an angle of attack of 4°, as seen in Figure 18c.

Morphing Trailing Edge Optimization
The optimization process defined in Section 3.1 was implemented on the MTE airfoil by using the same two optimization functions defined by Equation (1) for drag minimization and by Equation (5) for endurance maximization. The optimization of the Morphing Trailing Edge (MTE) airfoil shape for the three flight conditions given in Table  5 was performed.
The optimization processes were carried out for two types of MTE control points. In the first type, referred to as Opt. Case I, MTE is a smooth single shape with a morphing starting point and trailing edge end point. The leading edge is not allowed to move, while the trailing edge points are only allowed to move in the vertical direction with respect to the baseline airfoil.
In the Opt. Case II, the MTE is a three-segmented finger-like configuration. The points are outlined to represent the morphing trailing edge starting location, the length of the first segment and the length of the second segment. The morphing location was kept at x/c = 0.60, so that enough space was provided for the actuation mechanism, and which is where its upper and lower limits were given as constraints. The upper bounds and the lower bounds were chosen to prevent the generation of morphing airfoil shapes with no likelihood of having very good aerodynamic efficiency. XFoil was used as an aerodynamic solver and therefore as a powerful tool that allowed obtaining fast convergence solutions. The results were later validated by the use of a RANS models in Ansys/Fluent solver. Figure 19a illustrates the initial (reference) versus the final (optimized or morphed) airfoil shape, while Figure 19b illustrates the convergence trend for the minimum drag optimization function for an optimized airfoil shape with MTE at an angle of attack of 50. It is seen that the graph in Figure 19b trends to a constant after some iteration steps, which shows convergence to its optimal solution.   The pressure coefficients of the MTE airfoil cases at three different flight conditions are shown in Figure 22. It can be observed that the suction peaks of the optimized airfoils are higher than those of the baseline airfoil. The MTE airfoil results expressed in terms of "effective camber change" lead to a significant increase in the negative pressure value at the trailing edge surface. The low performance of the baseline airfoil compared to that of the optimized airfoils can be related to the unfavorable pressure gradient on its upper surface near the leading edge, thereby causing an earlier laminar-turbulent transition of the boundary layer by slightly increasing the drag.   Figure 23a illustrates the initial (reference) versus the final (optimized or morphed) airfoil shape, while Figure 23b presents the convergence graph of the aerodynamic endurance optimization for the given MTE optimized airfoil while utilizing the PSO algorithm with pattern search. Figure 23b shows that the algorithm performs well and improves the airfoil shape to maximize the endurance. The optimization process using the endurance maximization objective function led to improved endurance performance in the optimized airfoils with respect to the endurance performance of the baseline airfoil. The performance measured in terms of CL 3/2 /CD variation with CL is shown for three different flight conditions in Figure 24. The endurance performance of each optimized airfoil shows a higher value in the MTE Opt. Case I with respect to its baseline airfoil configurations. For the same endurance performance defined by the CL 3/2 /CD ratio, the lift coefficients were found to be 1.17 for the baseline airfoil and 1.45 for the MTE airfoil configuration for the optimized flight condition II, respectively. Furthermore, as seen in Figure 24, the efficiency of CL 3/2 /CD increased by approximately 10.25% for flight condition III; this finding indicates the improved endurance performance of the UAS-S45 airfoil. The difference between the results obtained by the two objective functions of drag minimization and endurance maximization to obtain optimal airfoil configurations was determined. The optimized airfoil geometrical shapes were obtained for drag minimization and endurance maximization as shown in Figure 25. The main difference observed in the airfoil geometries, as shown in Figure 25, is that the trailing edge obtained with the drag minimization objective function gave a smaller deflection than the baseline airfoil. Likewise, the endurance maximization objective function requires higher airfoil deflection with a continuous trailing edge. The contour plots of the velocity magnitude and of the turbulent kinetic energy are visualized for the baseline airfoil and for the airfoils optimized by the Transition SST turbulence model in Figure 26. These contour plots reveal that for a given MTE deflection, larger TE separation regions are found at 4° angles of attack than at 10° angles of attack. Similarly, the vertex of turbulent kinetic energy (TKE) originating from the MTE at an angle of attack of 40° has more strength than the vertex an the angle of attack of 10°. The TKE contour plot of the baseline airfoil at an angle of attack of 4° is approximately comparable to the TKE contour plot of the MTE airfoil at an angle of attack of 10°. Therefore, the MTE-optimized airfoils can provide increased aerodynamics performance with respect to the performance of the original airfoils.

Conclusions
This study was conducted to perform aerodynamic optimization of Morphing Trailing Edge (MTE) and Droop Nose Leading Edge (DNLE) airfoils for the UAS-S45 at different flight conditions. A new hybrid optimization technique was chosen by coupling the Particle Swarm Optimization (PSO) with the Pattern Search (PS) algorithms. The optimization function was designed to minimize drag with respect to given constraints such as airfoil lower and upper bounds, as well as to increase endurance in given flight conditions. The Bezier-PARSEC (BP) technique was used to parameterize the baseline airfoil shape, as well as to obtain its various optimized morphing configurations by using different constraints on the morphing of both the Morphing Trailing Edge (MTE) and the Droop Nose Leading Edge (DNLE) of the wing.
Within the aerodynamics optimization framework, the low-fidelity solver XFoil and the high-fidelity CFD solver Ansys Fluent were used. The results obtained using both solvers were compared for their validation. Specifically, the flow transition was predicted using Menter's Transition SST turbulence model. In addition, the optimization framework was done from an aerodynamics perspective, and therefore, no structural studies were taken into consideration.
Both DNLE and MTE airfoil optimizations were aimed to increase the aerodynamic performance of the UAS-S45 for a wide range of angles of attack. For each of the three considered flight conditions, enhanced aerodynamic efficiency was obtained by optimizing the morphing airfoil design with respect to its UAS-S45 baseline airfoil. The optimization results have shown an increase of lift coefficients in DNLE airfoils until the stall angle of attack and thus a delay in the stall angle. An improvement in the lift coefficient was produced for the DNLE airfoils with respect to the baseline airfoil of up to 21%. In addition, an increase of 9.6% was obtained for the maximum lift coefficient, and the stall angle was also delayed by 3 degrees. Similarly, the aerodynamic performance showed a significant improvement for the MTE configurations. An increase in the maximum lift coefficient of up to 8.13% and of the efficiency of CL 3/2 /CD by 10.25% was obtained, thus indicating an increased endurance performance for the MTE airfoils of the UAS-S45. The delay of the leading edge separation by use of DNLE airfoils was another interesting result.
The flight conditions were chosen from the UAS-S45 manufacturer Hydra Technology user manual. These conditions do not cover aircraft maneuvers or gust responses. UAS-S45 maneuvers and gust response studies could be done in the future. Other objectives of this morphing optimization study will include the three-dimensional analysis of a UAS-S45 wing with a combined droop nose leading edge (DNLE) and continuous morphing trailing edge (MTE). The improvements at the performance level, such as the typical fuel savings with the optimized morphing concepts for a given flight mission, will also be studied in the future. Based on the aero-structural studies, various configurations will be analyzed, and the internal actuation mechanism will be implemented. Wind tunnel and flight tests could be further performed to validate the optimized UAS-S45 model.

Data Availability Statement:
The data presented in this study are available on request from the corresponding author.
Acknowledgments: Special thanks are due to the Natural Sciences and Engineering Research Council of Canada (NSERC) for the Canada Research Chair Tier 1 in Aircraft Modelling and Simulation Technologies funding. We would also like to thank Odette Lacasse for her support at the ETS, as well as to Hydra Technologies' team members Mr. Carlos Ruiz, Eduardo Yakin and Alvaro Gutierrez Prado in Mexico.

Conflicts of Interest:
The authors declare no conflict of interest.