Efficient aeroelastic wing optimization through a compact aerofoil decomposition approach

Efficient optimization of an aeroelastic wing is presented through multi-disciplinary analysis using low-dimensional modal design variables. Much work in wing optimization has concentrated on high-fidelity surface control, therefore utilising often hundreds of design variables. However, whilst fine surface control can be useful, problems can arise such as large disparities in design variable values when planform variables are introduced, slow convergence speeds, and lack of compatibility with global algorithms. Therefore, the focus of this paper is to filter the design space of this problem to reduce the dimensionality and complexity of the problem. Orthogonal geometric design variables are derived in the geometric space via singular value decomposition. Orthogonality of design variables leads to a well-conditioned design space and ensures effective optimizer convergence. These variables are applied in a sectional fashion for fixed planform drag minimization of a flexible transonic wing, using a gradient-based optimizer. Shock-free solutions are demonstrated when optimizing a rigid wing, indicating suitability of the aerofoil modes for sectional-based wing optimization. However, it is shown that these wings have poor performance when subsequently deformed under flight loads, hence optimisation including full aeroelastic performance is performed. Encouragingly, shock-free solutions are again computed. Loading is shifted outboard, leading to increased tip deflection. Monotonic improvement in the objective function (drag) with increase in dimensionality is also proven. Furthermore, applying these sectional deformation modes globally across the wing with only 10 design variable leads to a 28% drag reduction, which is within 7 drag counts of when the modes are applied locally through 82 design variables. This therefore opens the possibility of introducing global optimization algorithms to high-fidelity aeroelastic wing optimization.


Introduction
The design of an aircraft is inherently an exercise in the coupling and subsequent compromise of multiple, and often independent, disciplines. Nowhere is this more apparent than in the design of aircraft wings, where the aerodynamic performance typically dominates, but where the design is dictated by the behaviour of many other disciplines. Multidisciplinary analysis (MDA) is the framework that proceeds when analysis of each discipline is integrated into an overall coupled system. Furthermore, beyond analysis, the overall multi-disciplinary system may be linked to an optimization process leading to the multi-disciplinary optimization (MDO) framework (Martins and Lambe 2013).
Commonly, engineering optimization problems embed high-fidelity computational analysis into the process, which is used to determine various metrics against which to optimize. The problem then reduces to that of finding the optimum design (often a geometric problem) that optimizes a cost function subject to restrictions imposed by constraint functions. The cost and constraint functions typically rely on outputs from the analysis. An optimization algorithm is then required to link the analysis to the design process, where the algorithm uses values of the cost and constraints to determine a vector of design variables that link to changes in the design to hopefully improve the cost. A parameterization is the vehicle that governs how the design variables affect the design. An example single-discipline problem is finding the optimum aerodynamic aerofoil or wing shape to minimise drag subject to a constraint on lift, where a computational fluid dynamics (CFD) solver is used to determine the aerodynamic performance; see Hicks and Henne (1978) Qin et al. (2004), Morris et al. (2008), Allen and Rendall (2013), Poole et al. (2017Poole et al. ( , 2018b, Destarac et al. 2018) for example.
In aircraft wing design, the aero-structural-coupled system dominates the performance. When determining aerodynamic performance of the wing in flight, the rigid wing shape is deflected as aerodynamic loads are transferred to the structure and an equilibrium state results. Simulation of this system commonly requires coupling of separate CFD and computational structural mechanics (CSM) solvers in a partitioned manner (Geuzaine et al. 2003;Taylor et al. 2004;Woodgate et al. 2005;. This solution approach has the advantage of requiring little development of new codes, with the main development requirements being the solution coupling, and the force and displacement interpolation and mesh deformation needed due to non-coincident fluid and structural meshes. The coupling requires that an outer iteration loop be added to march the solution to equilibrium. The interpolation and mesh deformation are often handled separately, but a unified approach using radial basis function (RBF) interpolation was shown to be highly effective . A similar approach, but using B-spline-based methods, has also recently been presented (Zhang et al. 2017).
An alternative solution approach is the monolithic method, where a single solver that solves the governing equations synchronously is developed (Michler et al. 2004;Hübner et al. 2004). Whilst accuracy and convergence tends to be better with the monolithic approach, the development costs are a barrier. For example, all four codes that submitted results to the aero-structural benchmarking case of the Sixth AIAA CFD Drag Prediction Workshop used a partitioned solution approach (Keye and Mavriplis 2018).
Once a suitable aero-structural solver is developed, full MDO of the wing that optimizes for the performance of the coupled aerodynamic and structural responses can be realised. There exists a deep history of performing ever increasing fidelity of aero-structural optimization (Sobieszczanski-Sobieski and Haftka 1997). Due to the cost associated with aero-structural optimization, local gradient-based optimization algorithms have usually been the tool of choice. The problem then exists of determining the gradients of the cost and constraints with respect to the design variables. The adoption of the adjoint approach (Jameson 1988) permits all gradients of flow quantities with respect to surface changes to be evaluated for a computing cost that is similar to a single flow solution. Using this, coupled aero-structural adjoint methods have been presented (Maute et al. 2001;Martins et al. 2005;, which has permitted highfidelity, large-scale aero-structural optimization. Furthermore, recent developments include the ongoing work to couple more disciplines into the process, including coupling aerodynamics, structures and acoustics in a fully coupled unsteady optimization of rotorcraft (Mavriplis et al. 2017).
Large-scale, high-fidelity MDO lends itself to detailed shape optimization. Such problems often have very large numbers of design variables. However, the quantity of design variables often has a significant effect on the convergence of the optimizer, as well as the definition of multimodality in the problem. Furthermore, whilst gradient-based methods are the most popular optimization algorithms in MDO, if alternatives (such as population-based) methods are to be exploited, then lower-dimensional design spaces become key, particularly as it been shown that the number of function evaluations increases quadratically for increasing numbers of design variables (Lyu et al. 2014). An alternative approach can also be to employ multi-start gradient (Chernukhin and Zingg 2013). Dimensionality reduction techniques have become a useful approach for reducing the dimensionality associated with aerodynamic shape optimization. Of these, singular value decomposition-based approaches take a training matrix of data (for example a number of aerofoil shapes (Toal et al. 2010;Ghoman et al. 2012)), and project a reduced-order basis approximation of the original data. Work by the authors has shown that this is a very efficient approach for producing a reduced set of aerofoil deformation modes (Poole et al. 2015a) that are suitable for aerodynamic optimization of aerofoils  and wings . Furthermore, comprehensive experiments (Masters et al. 2017) have shown singular value decomposition (SVD) modes to be the most efficient approach at representing a generic aerofoil compared to most other commonly used parameterization methods, and have shown to produce single figure drag values for the AIAA ADODG Test Case 1 (Kedward et al. 2020b), when combined with a surface gradient limiting constraint (Kedward et al. 2020a).
The work presented in this paper extends the application of these modes further by considering high-fidelity aerostructural wing optimization. The overall goal of this work is to make high-fidelity global multi-disciplinary optimization of wings a reality, and the use of a compact aerofoil decomposition for the design variables is a pillar to this goal. Whilst fixed planform optimization of rigid aerodynamic wings is generally considered to be unimodal (Yu et al. 2018), introducing planform changes along with sectional changes introduces clear multimodality (Poole et al. 2018b;Streuber and Zingg 2021). Furthermore, similar behaviour is also likely to be seen in aeroelastic wing optimization so there is a clear need to consider how global methods could be introduced into high-fidelity wing optimization. Specifically, the objectives of this work are to first determine the suitability of the aerofoil decomposition design variables for aeroelastic optimization and ensure no adverse effects are found when structural responses are included in the problem. Second, this work aims to determine whether a compact wing deformation approach with only a handful of sectional variables, can be successful in high-fidelity elastic wing optimization.
The remainder of the paper is organised as follows: the compact modal decomposition and how this is applied to wing optimization is outlined in Sect. 2; the solver and datum solution is presented in Sect. 3; the optimization problem and chosen optimizer is presented in Sect. 4; optimization of the rigid wing and aeroelastic wing are presented in Sects. 5 and 6 , respectively; finally, conclusions are given in Sect. 7.

Decomposition for wing design variables
The design variables used in the optimization process are the weightings of various sectional deformations. These sectional deformations come about by performing a matrix decomposition that uses SVD on a training library of aerofoils. The resulting modes represent a reduced basis projection of the full-basis aerofoil design space. The modes are also orthogonal, which leads to a well conditioned design space, and no coupling of variables in the geometric space. This aids optimizer convergence considerably. The surface deformation modes are defined in the two-dimensional surface space, so need mapping onto both the three-dimensional surface space and the volume space. This is achieved in a unified manner through radial basis function (RBF) interpolation via a set of control points. As RBF interpolation only needs a point cloud and no connectivity, it is general. This also leads to aerofoil surface modes being general.

SVD for sectional deformations
To obtain aerofoil deformation modes, a training library of aerofoils needs to first be collated. The selection of the training library is one of the most important steps in this process since the characteristics of the library map to the characteristics of the deformation modes, therefore, careful selection of the training library is key to allowing a wide design space to be represented and therefore allowing a variety of aerofoil shapes to be explored. The authors have previously published work (Poole et al. 2015a) that used a metric-based filtering approach to select the training library. In this work, the library of aerofoils is as previously used by the authors for transonic aerofoil optimization ). The aerofoil data was obtained from the UIUC database 1 , and was subsequently smoothed and re-parameterised to ensure a consistent fractional arc-length distribution of points.
A training library contains M aerofoils each parameterised with N surface points, where the i-th surface point of the m-th aerofoil has coordinates (x i m , z i m ) . To obtain aerofoil deformation modes, the vector difference between each surface point of all aerofoils is computed. The vector difference of the i-th surface point between the m-th and n-th aerofoils is given as (Δx i m,n , Δz i m,n ) . It is worth noting that one option is to consider vertical deformations only (in this case, the aerofoils would need to have equal horizontal point distributions) but this can be restrictive in areas of high curvature, such as the leading edge so both horizontal and vertical deformations are considered. More detail is given in (Poole et al. 2015a).
All deformations are collated into a single deformation matrix that has 2N rows and M def = M(M − 1)∕2 columns: The deformation matrix has an SVD given by: where is an orthonormal 2N × 2N matrix, is a diagonal matrix with min{2N, M def } diagonal entries arranged in descending order, and T is an orthonormal M def × M def matrix. The columns of contain the aerofoil deformation modes, and these are extracted column-wise and used for the optimization. The diagonal entries of contain the singular values, which may be used to determine now much 'energy' each mode has in the overall system.
The use of the SVD here is for geometric filtration of the high-degree training matrix into a low-degree representation using deformation modes. To do this, the first D modes (i.e. first D columns of ) may be extracted and a D-rank approximation of the original matrix may be obtained by: where the tilde is used to denote reduced forms of the SVD matrices. Once a low-rank approximation is found through SVD, the following is true (Eckart and Young 1936): where is any matrix of rank D and ‖ ⋅ ‖ F is the Frobenius norm. Hence, the error between the low-rank approximation and the full rank approximation will always be at least as (2) good as the error between any other k-rank matrix and the full-rank matrix. In this sense, the SVD produces an optimal low order projection of the higher-dimensional space into the lower-dimensional one.
The first four deformation modes are shown in Fig. 1 (these are superimposed on NACA0012 for visualisation purposes). It is clear to see that the first mode represents a thickness change, then the second represents a camber change, indicating that for this library of aerofoils, these are the two most important aerofoil design parameters respectively.
Given the required number of deformation modes, D, in the optimization, the design variables are the weightings of each deformation mode. The overall deformation is then a linear superposition of each mode: where i is the design variable relating to the i-th mode and i is the i-th mode, which is the i-th column of . Since the modes are independent of one another, the design space of lower-dimensional problems is entirely contained within the design space of higher dimensional problems. As such, adding further modes can only lead to at least a better optimum objective function. For aerofoil drag optimization, experience of the authors shows that there are normally few benefits of going above 10 modes, and that if there is an improvement in the objective function, this is outweighed by the additional cost due to increases in optimizer iterations and gradient evaluations. This is therefore an extremely efficient method of performing aerofoil optimization.

Application to wing deformation
The aerofoil deformation modes are surface deformations, however, to ensure body-fitted meshes are retained in the optimization the mesh also needs to deform. Furthermore, the surface deformation modes are applied sectionally, which is defined here.
The sectional deformations are applied at a fixed number of spanwise stations, i.e. Eq. 4 is applied at these stations locally. However, the sectional deformations can also be applied in a global nature. The sectional deformations are applied using the RBF control point approach, where a set of control points are defined in the fluid domain and global volume interpolation translates deformation of the control points to deformation of the aerodynamic mesh. Hence, the modal deformations are used to drive deformation of the control points that subsequently deform the wing surface and mesh. These deformations are decoupled, so the control point modal deformations are determined off-line and then applied in the optimization process. This has the added advantage of producing modal deformations that are not restricted to the original discretisation (i.e. number of points -rows of Eq. 1 -and distribution of points) used to derive them.
At the heart of this technique in an RBF interpolation developed originally for aero-structure coupling and mesh motion by Rendall and Allen (2008). An RBF interpolation, s, is a linear combination of basis functions, whose argument is the Euclidean distance, ‖ ⋅ ‖ between the point to be interpolated in the domain, , and the N points in the known data set. Therefore, the influence that a known point has is controlled by a function, , that depends on the distance from the interpolated point: Control points decouple the shape deformations from the surface mesh and provide a unified framework for surface and mesh deformation. Given n c control points, a global RBF interpolation of this nature provides exact recovery of data at known sites, and interpolation of that data away from the sites. In the case of optimization, the data to be interpolated is deformation of the control points, hence a deformation field is created. The position of the aerodynamic mesh points in the field therefore defines the deformation of those points. Since exact recovery of data at the know sites (in this case the position of the control points, which for the j-th control point is defined as (x c j , y c j , z c j ) ) is specified, the interpolation takes the form: Analogous definitions hold for the y and z coordinates. The radial basis function i,j = (‖ c i − c j ‖) can take a number of forms, but the radially decaying functions of Wendland (2005) are a good choice for the mesh deformation problem to give the interpolation a local character and ensure deformation is contained in a region near the moving body. The C 2 function is used here. Once the linear system is solved, the resulting deformation field can be evaluated at the location of each mesh point. The deformation of an aerodynamic mesh point is given by: with analogous definitions for y and z.
Using RBF interpolation has the advantage of being able to specify the level of control since control points can be placed arbitrarily in, on the boundary, or outside the fluid domain. However, irrespective of the location of the control points, modal deformations defined on the aerofoil surface may not be coincident with the control points so control point deformations must be defined. The authors (Poole et al. 2015b) showed a number of techniques for achieving this, and an inverse RBF interpolation is an effective approach. The set-up of control points around the wing used in this paper is shown in Fig. 2.
Modal deformations are applied sectionally at ten spanwise stations (on the red control points), whilst intermediate control points (shown in blue) are used to permit smooth spanwise deformations between the deformation slices   where i = (y i − y m )∕(y m+1 − y m ) . Hence, it can be seen that if the main deformations are equivalent of two neighbouring sections ( Δ m = Δ m+1 ) then any intermediate slice will have the same deformation. Further, it is therefore simple to also define global deformations, where each deformation slice is deformed by the same mode and scaled by the local chord. This means that wing optimization can be performed using only a handful of design variables. The global approach is given in more detail by Allen et al. (2018).
In addition to the sectional deformations, a global twist deformation is introduced. A linear twist is applied from zero at the root to the twist angle at the tip. Hence, the local angle of rotation, , of a main control point slice located at y m is given by (y m ) = y m ∕s.
The control point cage is constructed around the local wing coordinates. For example, Fig. 2 shows the cage around both a rigid wing and a deformed wing. Hence, the sequence of operations is given a design variable vector, , twist deformations of the main sections occurs, followed by sectional deformations at main spanwise stations (if a global deformation is used, then each main spanwise section is deformed by the same mode scaled by the local chord

Aero-structural solver and case definition
The aero-structural solver is based on the structured multiblock solver of Allen (2002Allen ( , 2006, with the coupled solver added through the work of Rendall and Allen (2008). In the case considered in this work, the flow is governed by the Euler equations, which are solved using finite volume integration with the Jameson-Schmidt-Turkel scheme (Jameson et al. 1981). Convergence acceleration is achieved through multigrid (Allen 2001). A modal structural solver is applied where coupling between the aerodynamic and structural grids, and subsequent mesh motion, is achieved through RBF coupling ) that uses a reduced point cloud between aerodynamic and structural neighbour nodes (Rendall and Allen 2010). Strong coupling of the two systems is employed for dynamic calculations. Newmark temporal integration (Newmark 1959) is used to march the solution.
To demonstrate using the compact aerofoil decomposition for aero-structural optimization, the MDO wing (Allwright 1997;Haase et al. 2002) is used, where the structural model and associated structural mode shapes are defined. It should be noted that the geometry referred to through this paper as rigid, is defined in Sec 8, Chapter 4 of (Haase et al. 2002) as 'jig-shape'. This had been back calculated from the designed flight shape for the purposes of the aeroelastic computations of the UNSI project. The aeroeastic solutions in this paper represent the calculated flight shape at the specified loading condition. Figure 3 shows the MDO wing planform. Also shown is the aerodynamic surface mesh and the structural grid. A 2515 node wing-box structural model is used with the modes defined by Haase et al. (2002).
Throughout this work, a trimmed cruise condition of C L = 0.4 at M ∞ = 0.85 is used. Unless otherwise stated, all wings are trimmed to this condition. The MDO wing has 18 defined structural modes and all are used for this work.
To ensure sufficient aerodynamic resolution, a mesh dependence study is presented. A family of eight-block structured C-meshes (block structure is given in Fig. 4 was generated using the methods of Allen (2008) to give high quality meshes. These range in size from 2.1 million to 0.13 million cells, and are designated L1 (2.1M), L2 (1.1M), L3 (580k), L4 (260k) and L5 (130k); sizes were chosen with approximately a two-times scaling between mesh levels, and to maximise the number of multigrid levels for each mesh.
Each mesh was used to produce both a rigid and aeroelastic solution and the final force coefficients of each run are given in Table 1. The structural deflection is calculated based on the deformation of the structural node closest to the aerodynamic tip. Also given are the run-times relative to the run-time of the coarsest mesh, with the runs being performed in serial to obtain comparable figures. All wings were run at a trimmed cruise C L of 0.4. Figure 5 gives the surface pressure coefficients of the coarsest and finest meshes. Clearly, there is little difference in the flowfields, and the structural deformations are almost identical. The L4 mesh, whilst being relatively coarse, appears to be a good compromise between run-time and accuracy whilst providing accurate capture of the structural deformation field.
The primary mesh used for the study is the L4 mesh, shown in Fig. 6, which contains 273k nodes. This mesh has a 97 × 57 surface mesh, 21 nodes on either side of the wake, and 25 nodes between the inner and outer boundary. Figure 7 shows the solution of the rigid MDO wing and the aeroelastic wing (i.e. in flight shape). and Fig. 8 gives the span-wise loading distribution. There is a strong stock present on the rigid wing, that when in flight shape is no  longer present outboard. The loading is shifted inboard and the distribution is almost perfectly triangular in flight shape. This results in a lower root bending moment at the expense of a drag penalty.

Optimization problem and algorithm
A generic single-objective optimization problem optimizes a cost function, J, which is a function of a vector of D design variables, , subject to a vector of inequality constraints, , and equality, , constraints. The problem considered in this paper is aerodynamic shape optimization. The objective is drag minimization subject to constraints on lift coefficient, C L and internal volume V. The problem is given by: The design variables of the problem are aerofoil deformation modes applied at ten equally spaced sections across the span of the wing. Furthermore, to allow induced drag to reduce,  Table 2. Using structural modal analysis is an effective means to determine the coupled aero-structural performance of the wing. However, when used in an optimization process, the primary assumption is that modifying the aerodynamic shape has minimal effect on the underlying structural model so the same structural modes can be used. This assumption holds assuming the structure is not changed, and since a wing-box is modelled, no planform changes occur but only minor surface changes, this is reasonable.
In the case considered in this paper, the overall goal of this work is to permit large-scale global aeroelastic shape optimization. Global optimization is particularly useful when the design space is known to be multimodal, and the overall goal is to locate the globally optimal solution. In the first instance, only sectional changes are considered and since fixed planform wing optimization is generally considered to be unimodal (to within numerical tolerances) (Yu et al. 2018), the gradient-based optimization algorithm, feasible sequential quadratic programming (FSQP) algorithm as implemented in version 3.7 (Zhou et al. 1997), is used. FSQP is based on the sequential quadratic programming (SQP) approach, but modified to avoid the so-called 'Maratos' effect (restriction of a step size due to the requirement of feasibility) (Maratos 1978), by combining a search along an arc (Mayne and Polack 1982) with a non-monotone procedure for that search (Grippo et al. 1986). The non-monotone line search relaxes the rule that a step must result in an improvement in the objective function at each major iteration, and instead requires an improvement over the best of the last four iterations. This means that the objective function can occasionally worsen from one iteration to the next. The Hessian is updated via the Broyden-Fletcher-Goldfarb-Shanno (BFGS) scheme, with identity initialisation. As is common with gradient-based optimization, the convergence is not based on a fixed number of iterations, rather on ensuring either convergence of the Karush-Kuhn-Tucker (KKT) conditions, or no change in the objective function within the non-monotone procedure, whichever occurs first. The FSQP algorithm is fully described and analysed in (Panier and Tits 1991;Bonnans et al. 1992).
The gradient-based optimizer requires the sensitivities of the cost and constraint functions with respect to each design variable at each major iteration. For this work, a secondorder finite difference stencil is used so for each design variable, two extra flow solutions are required (one each for the positive and negative perturbations) to evaluate the sensitivities. For computational efficiency, a parallel decomposition of the gradient evaluation is employed such that each design variable sensitivity is assigned to its own CPU, which

Rigid wing optimization results
Initially, to determine both the optimization set-up and a suitable number of SVD modal design variables, optimization of the rigid MDO wing is performed at the trimmed cruise flight condition ( C L = 0.4 at M ∞ = 0.85 ). This is performed using the local modal deformation cases only (i.e. 42, 62 and 82 variables), with global modes (i.e. 6, 8 and 10 variables) applied to the aeroelastic case to demonstrate the flexibility of this approach. Table 3 gives the optimization results of the optimizations, whilst Fig. 9 gives the convergence of the optimizer. For all three of the design variable combinations tested, both constraints are active and the drag has reduced substantially. Of particular interest is the monotonic decrease in the acquired objective function with increasing number of design variables. In the geometric space, the modal design variables are perfectly orthogonal, which has the implication that the design space of lower numbers of design variables is contained within the design space of higher numbers. In the aerodynamic space, this is not necessarily the case; though the authors have shown this to be the case for aerofoil optimization . Albeit, the monotonic decrease indicates that geometrically orthogonal aerofoil design variables exhibit similar performance when applied to wing optimization. Figure 10 shows the surface pressures of the baseline wing and 82 variables optimized wing, and demonstrates that the result is shock-free at the design condition. Since this problem is inviscid, the two primary sources of drag are wave drag and induced drag. As the problem is shock-free, wave drag has been substantially reduced.

Aeroelastic solution of optimised rigid wing
The rigid optimizations have produced the expected shockfree result at the design condition, however, it is interesting to see how this wing performs under aeroelastic loading. Hence, the optimized wing is run with the structural model to determine the aeroelastic solution. Table 4 gives the force coefficients of the aeroelastic MDO wing and the aeroelastic solution of the optimized rigid wing, whilst Fig. 11 shows the surface pressure coefficients. It should be noted that all wings are trimmed to C L = 0.4 . Clearly the aeroelastic solution of the optimized rigid wing has a shock structure, with a single shock at the inboard and outboard sections, and a double shock midspan and this has resulted in an overall drag increase. Also, the extra outboard load of the rigid optimized wing has resulted in a lower tip deflection.
It is well known that single-point drag minimization produces highly optimized point-design solutions with poor off-design performance (Poole et al. 2018a), so an increase in drag for the aeroelastic solution of a rigid wing is not surprising. The changes in loading that occur due to shape changes in optimization, lead to a change in the deflected  shape which leads to a further change in the loading. It is therefore not enough to consider optimizing a wing statically, but that the full aero-structural solution must be taken into account to obtain any meaningful results.

Aeroelastic wing optimization results
As a result of the experiment presented above, it is clear that wing optimizations should contain wing shape responses also. In that light, aerodynamic shape optimizations are presented of the aeroelastic MDO wing. Table 5 shows the optimization results for varying numbers of locally applied modes, and Fig. 12 gives the optimizer convergence. As in the rigid case, monotonic decrease in the objective function is observed, which is to be expected with the modal design variables, but is encouraging to observe. Both the lift and volume constraints are active for all the cases. The tip deflection has increased substantially for the optimized wing indicating that transferring load from root to tip (the load is almost perfectly triangular in the MDO wing, but close to elliptical in the optimized wing) and allowing a more flexible solution (something observed in modern transonic aircraft) permits shock strength reduction. This is demonstrated in the surface pressure shown in Fig. 13.

Local modal deformations
Comparing the final optimized shapes for both the rigid and aeroelastic optimizations, Fig. 14 shows the surface changes made by the optimizer to the wing shape for both the upper and lower surfaces. It is clear that there are substantial differences in the final optimized shapes. The rigid wing case has pushed the trailing edge up, whereas the  aeroelastic case has pushed this down. Furthermore, the magnitude of these deformations are much larger in the aeroelastic case than the rigid one. It is therefore not surprising that a rigid optimization subsequently run elastic, produces a completely different solution to an optimization that is run using the elastic case. It is therefore clear that to perform wing optimizations, the coupling between aerodynamic forces and structural response must be included.

Global modal deformations
Finally, a study is presented on an approach for performed extremely efficient aeroelastic wing optimization. The use of the framework developed in this paper means that aerofoil modal deformations can be applied across the whole wing (scaled by the local chord). This results in extremely low-dimensional wing optimization problems, with an order   of magnitude fewer design variables than including local deformations. Such lower number of variables does mean that very localised deformations of the surface are not permissible, but it is interesting to consider whether performing fixed planform, sectional optimization of elastic wings is permissible with only a few design variables. Table 6 gives the final optimization results for the global sectional deformation runs, whilst the optimizer convergence is given in Fig. 15. Like the local runs, a monotonic decrease in the objective function (drag) is observed, which is expected since use of these modal parameters leads to the design spaces of lower numbers of modes being entirely enclosed within that of higher numbers of modes. Furthermore, it is extremely encouraging to see that large drag reductions are still possible with a very small number of design variables (especially for full 3-D wing optimization). The optimizer converges rapidly, though like in the local cases, a small 'tail' of the convergence is seen. This can be common in aerodynamic optimization, where the optimizer is tweaking the final surface to reduce the small amount of  Table 7 provides a summary of the total number of CFD calls required during the optimization. The global application naturally has a much reduced total cost across all three cases tested. This comes about primarily due to requiring fewer gradient evaluations but also generally a reduction in the number of iterations the optimiser required to converge. This is particularly prevalent in the 8 modes cases, where global application of the modes shows a 70% reduction in the cost of the optimization, with only a small drag penalty.    Figure 16 compares the surface pressures for the 8 aerofoil modes applied locally (82 variable case) and globally (10 variable case). Broadly, the pressure distributions are very similar, which is encouraging. Figure 17 compares the loading of these wings showing similar shapes, though with the local modes having the fidelity to shift loading outboard (and therefore also have a slightly higher tip deflection). Due to the fidelity of the local deformations, these have successfully managed to eliminate the shock, whilst in the global case, this is not the case. However, comparing the final drag values, the difference between the two cases is only seven drag counts, which is 4% of the initial drag. This is a performance penalty over the local deformations, though it should be noted that this is inviscid flow so if viscous drag was included, the percentage penalty would be different. Still, it shows that a highly compact design space is still explorable with these modes. The low number of modes has not had an appreciably negative effect on the aeroelastic solution demonstrating that extremely low-dimensional flexible wing optimization is a reality. Having such a low number of design variables for elastic wing optimization opens up the possibility of introducing global optimization methods (where a very low dimensional design space is necessary) to high-fidelity wing optimization, particularly where coupling may occur between sectional and planform changes.
A final comparison of all the wings considered here is shown in Fig. 18, emphasising the different in flight shapes of the optimised wings.

Conclusions
Efficient multi-disciplinary analysis and subsequent optimization have been presented for aeroelastic wing optimization using an efficient low-dimensional surface deformation scheme. Whilst much aero-structural optimization work to date has considered high-fidelity, high-dimensional problems, little work has considered the possibility of reducing the dimensionality of the problem. As such, the authors have previously presented an efficient method for determining aerofoil deformations via singular value decomposition. The resulting design variables are geometrically orthogonal resulting in a well-conditioned design space. The work presented in this paper has developed an approach to apply these to elastic wings in such a way that the aerofoil deformations can be applied locally to each section or globally (scaled by the local chord). This flexibility allows the user to perform localised optimization using many design variables, but also low-dimensional optimization through global sectional deformations.
The compact aerofoil design variables were applied sectionally for inviscid drag minimization subject to lift and internal volume constraints. The MDO wing was considered, with optimizations presented on the rigid wing and elastic wing. A gradient-based optimization approach with parallel gradient evaluation was utilised. Shock-free solutions were demonstrated for the rigid wing, indicating suitability of the aerofoil modes for sectional-based wing optimization. However, it was demonstrated that a wing designed in rigid shape had poor performance when loaded in flight shape. Hence, the optimization framework was applied to the aeroelastic wing to produce shock-free solutions. The resulting shapes had a much increased tip deflection. The orthogonality  behaviour of the modes is preserved through to the optimization and lead to monotonic improvement in objective with increase in dimensionality. No adverse optimization performance of these modes were found when structural response was included, so this approach is suitable for such optimization problems. Further, utilising the flexibility of the approach, global sectional optimization was performed and it was shown that for an equivalent set of modes to localised deformation (10 design variables for global deformations and 82 for local), that there was only a 4% difference in the final optimization result. The global deformations therefore show great promise when it comes to being utilised with a global optimization algorithm, where low dimensionality is a key since cost scales quadratically with design variables for these algorithms. This opens up the possibility of performing high-fidelity aeroelastic wing shape optimization with global optimization algorithms, which could have particularly use when considering problems that change planform as well as section.