Concurrent wing and high-lift system aerostructural optimization

A method is presented for concurrent aerostructural optimization of wing planform, airfoil and high lift devices. The optimization is defined to minimize the aircraft fuel consumption for cruise, while satisfying the field performance requirements. A coupled adjoint aerostructural tool, that couples a quasi-three-dimensional aerodynamic analysis method with a finite beam element structural analysis is used for this optimization. The Pressure Difference Rule is implemented in the quasi-three-dimensional analysis and is coupled to the aerostructural analysis tool in order to compute the maximum lift coefficient of an elastic wing. The proposed method is able to compute the maximum wing lift coefficient with reasonable accuracy compared to high-fidelity CFD tools that require much higher computational cost. The coupled aerostructural system is solved using the Newton method. The sensitivities of the outputs of the developed tool with respect to the input variables are computed through combined use of the chain rule of differentiation, automatic differentiation and coupled-adjoint method. The results of a sequential optimization, where the wing shape and high lift device shape are optimized sequentially, is compared to the results of simultaneous wing and high lift device optimization.


Introduction
Although knowledge of the physics of high-lift devices (HLD) has come a long way since the fundamental paper of A.M.O Smith on high-lift aerodynamics in 1975 (Smith 1975), analysis and optimization of high-lift devices still proves to be a difficult subject. Through the use of Computational Fluid Dynamics (CFD) and increased computing capabilities, extensive research on the subject has become possible. In early days, this research mainly focused on achieving high-lift requirements to satisfy take-off and landing performance requirements. However, over the past years the focus has switched to reducing weight and complexity (van Dam 2002) as aircraft manufacturers tend to use less complex high-lift devices (Reckzeh 2003). The importance of weight and aerodynamic performance of high-lift devices in aircraft design is illustrated by Meredith (1993). According to Meredith, an increase of 0.1 in lift coefficient at constant angle of attack results in a reduction of approach attitude by about one degree, reducing landing gear length and thereby saving up to 1400 lb. Moreover, an increase of 1.5% in maximum lift coefficient (C L max ) may result in an extra 6600 lb payload at fixed approach speed while an 1% increase in take-off lift over drag ratio (L/D) is equal to a 2800 lb increase in payload or a 150 nm range increase.
Even though numerous semi-empirical methods exist to predict the wing weight, drag and lift of multi-element wings (Raymer 2012;Torenbeek 1982;Roskam 2000;Pepper et al. 1996), the accuracy of these methods does not yield the level of accuracy required by the industry, requiring e.g. a drag prediction accuracy of one drag count (van Dam 2003). To achieve the required accuracy, more physics based methods are required such as Computational Fluid Dynamics (CFD) and Finite Element Methods (FEM) tools. Example of application of such high-fidelity analysis for wing optimization can be found in the work of Martins et al. (2004), Kennedy and Martins (2014) and Barcelos and Maute (2008). The downside of these tools is that they require the use of high performance computational resources, making optimization problems in some cases too costly to solve. An alternative to the high-fidelity 3D aerodynamic solvers is the quasi-three-dimensional (Q3D) analysis methodology, which combines two-dimensional viscous airfoil data with inviscid three-dimensional wing aerodynamic data. This methodology requires only a portion of the computational power required for high-fidelity tools while generating sufficiently accurate results. Examples of using the Q3D method for aerodynamic analysis was presented by van Dam (2002), Elham (2015) and Mariens et al. (2014). Elham and Van Tooren developed a coupled-adjoint aerostructural analysis and optimization tool by coupling a Q3D method to a FEM (Elham and van Tooren 2016a). This tool has been validated for drag prediction and twist deformation. Using the coupled adjoint method, the tool is able to compute the derivatives of the outputs with respect to the inputs analytically enabling gradient based optimization.
In the traditional design methodology, the design of wing shape and HLD is done sequentially. The wing planform and airfoil shapes are designed (or optimized) first and then the HLD shape is determined (Flaig and Hilbig 1993;Nield 1995). It is known that sequential design and optimization may result in a sub-optimal design. In this paper a method for concurrent aerostructural optimization of wing and HLD is presented. In such a method the shape of the wing planform, airfoil, HLD as well as the wingbox structure is optimized simultaneously to minimize the aircraft mission fuel weight and satisfy the aircraft field performance requirements, that are the main drivers for HLD design.
The structure of this paper is as follows: First, the basic framework of the aerostructural analysis and optimization is described. Then the modifications applied to the aerostructural tool are explained, followed by a description of the method for predicting maximum lift. Then the method of coupling the modified methods is explained, followed by a validation of the extended model. Finally, a test case optimization is presented for a Fokker 100 class wing.

Aerostructural analysis and optimization framework
The present research expands the methodtool developed by Elham and van Tooren (2016a) for coupled-adjoint aerostructural analysis and optimization of lifting surfaces to include the analysis and optimization of HLD. This tool employs a Q3D method to compute the aerodynamic characteristics of a wing. Using the Q3D method reduces the computational effort to perform aerodynamic analysis but attains a high level of accuracy. In the Q3D methodology proposed by Elham and van Tooren, the total wing drag is then decomposed into three components: profile drag, induced drag and wave drag. The induced drag is calculated from a Vortex Lattice Method (VLM) analysis using Trefftz plane analysis based on the method of Katz and Plotkin (1991). The VLM analysis is also used to compute the lift distribution over the wing. Profile drag and wave drag are then computed using the viscous 2D solver MSES (Drela and Giles 1987;Drela 2013) following the strip theory. In this approach, the wing is divided into a number of spanwise sections (or strips), for which the aerodynamic forces and moments are computed using the effective flow properties, which are determined from the free stream flow properties taking into account the effects of sweep ( ) and downwash (α i ). The Q3D method is coupled to the structural analysis tool FEMWET (Elham and van Tooren 2016b), which simulates the wingbox structure using equivalent panels and computes the wing deformation using a FEM. The coupled aerostructural system is formulated using 4 governing equations R 1 to R 4 as follows: ⎡ (1) R 1 and R 2 are the governing equations of respectively the VLM and FEM. R 3 relates to the condition where the total lift is equal to the weight times the design load factor. R 4 states that the sectional viscous lift coefficient should be the same as the lift distribution from the VLM analysis, corrected for sweep. The system is solved using the Newton method for iteration where the updates on the state variables are determined using (2). ⎡ In order to facilitate gradient based optimization, the coupled-adjoint method is used Kenway et al. (2014). The sensitivity of any function of interest I with respect to any design variable x is computed by: where λ is the adjoint vector computed using the following equation: For a complete description of this coupled-adjoint aerostructural analysis and optimization method and derivation of the coupled-adjoint method, the reader is referred to Elham and van Tooren (2016a).

Maximum lift prediction
Since MSES more often than not fails to converge at high angles of attack, the Pressure Difference Rule (PDR), developed by Valarezo and Chin (1994) is used for the estimation of C L max . The PDR states that for a given chord Reynolds number and free stream Mach number, there exist a relation between the wing stall and the pressure difference between the suction peak and trailing edge pressure. While Valarezo and Chin made use of a higher-order panel method to obtain the pressure difference at several spanwise section, any other reliable method may be used such as the Q3D method described in this paper, due to the fact that empirical data is used in the analysis which takes viscous effects into account. Furthermore it was identified that this rule can be used for 3D wing analysis even though it relies on 2D sectional data. This is due to the fact that at the critical stall section, the suction peak of the 3D wing will be equal to that of the 2D flow for the respective airfoil section.
The PDR is implemented as follows: First, the effective Reynolds number is computed at several wing stations based on the local clean chord, taking into account sweep and downwash. Then, using the free stream Mach number, the critical pressure difference ( C p crit ) for a number of spanwise sections is computed from Fig. 1.
The effective pressure distribution over each specified spanwise section is then computed from a 2D linear strength vortex panel method, based on the method of Katz and Plotkin (1991), using the effective flow properties as described in Section 2. To analyze the airfoil using the panel method code the effective angle of attack and the effective Mach number are required. These values are obtained from the global angle of attack and free stream Mach number by adjusting for sweep effects and downwash (see Elham and van Tooren (2016a) for more details).
Besides producing the given outputs, the panel method is able to produce the derivatives of the outputs with respect to the inputs using a combination of the chain rule of differentiation and Automatic Differentiation (AD) in  Valarezo and Chin (1994) reverse mode using the Matlab AD toolbox Intlab (Rump 1999).
In order to account for the decambering effect of the boundary layer and wakes on aft segments of a multielement wing, Valarezo and Chin incorporated a reduction in effective flap deflection in their research (Valarezo and Chin 1994). The same flap reduction method has been applied to the Fokker 100 wing aerodynamic analysis (Obert 1986). The flap reduction angles used in both researches are shown in Fig. 2. Since the flap reduction angles of both researches are matching up to moderate flap angles of 20 • , it is assumed that the same reduction angles can be used for any flap configuration in this research.  Valarezo and Chin (1994) and Obert (1986)  Based on the PDR method the wing stall occurs when at one of the spanwise stations C p 2d = C p crit . Here C p 2d is the 2D pressure difference computed using the Q3D method, see Fig. 3.

Aerostructural coupling
In order to enable wing aerostructural analysis and optimization at high-lift conditions, the PDR described in the previous section is coupled with FEMWET. This further enhances the accuracy of predicting C L max as it takes into account aeroelastic effects in the analysis. The coupled system is defined following set of equations: Here the third governing equation is the Kreisselmeier-Steinhauser (KS) aggregation function (Wrenn 1989). The KS function is essential in establishing a single governing equation for the PDR condition, which allows for coupling this condition with the FEMWET analysis, effectively achieving an aerostructural prediction of C L max : where f k = C p crit k − C p 2d k . In this equation a value of 80 is used for ρ KS . The fourth equation in (5) states that the inviscid 2D lift computed by the panel method is equal to the VLM lift distribution, corrected for sweep. The system is solved using the same Newton method for iteration described in Section 2.
As mentioned before, the coupling of the PDR with FEMWET allows for taking aeroelastic effects into account in determining C L max . A representation of the coupling is presented in Fig. 4.

Sensitivity analysis
In order to use the Newton method for iteration, the partial derivatives of the governing equations with respect to the state variables (matrix J in (2)) are required. Additionally, to perform gradient based optimization, the sensitivities of any function of interest with respect to the design variables such as the wing planform or airfoil shape are required. The present tool computes all of the required derivatives through a combination of AD, chain rule of differentiation and the aforementioned coupled-adjoint method.
In Table 1, the partial derivatives of the governing equations of the coupled system (5) are presented. The first row is computed with relative ease where the partial derivative of R 1 with respect to is the Aerodynamic Influence Matrix (AI C) matrix and the partial derivatives of R 1 with respect to U , α and α i are computed through AD. Moving to the second row, the derivative of R 2 with respect to is computed using AD and the derivative of R 2 with respect to U is the stiffness matrix K.
Computing the third and fourth row requires more attention. Starting with row 3, the partial derivative of f k in (6) with respect to α and α i are computed analytically through: Moving on to the fourth row in Table 1, the partial derivative of R 4 with respect to may be computed using AD. The remaining terms in the fourth rows are a bit more challenging however. These terms are computed using the combined use of AD, chain rule of differentiation and the adjoint method within MSES.

Verification and validation
While the aerostructural tool developed by Elham and Van Tooren has been validated for wing drag and wing deformation (Elham and van Tooren 2016a), the enhanced method needs to be validated for maximum wing lift coefficient prediction and computation of wing lift over drag ratios in high-lift conditions. Finally, the sensitivities are verified through Finite Differencing (FD).

Maximum lift coefficient
For validating the accuracy of the C L max computation using the PDR, the 3D Royal Aircraft Establishment (RAE) experimental database (Lovell 1977) is used. For the current research, the basic body-off model was used without wing extension and with trailing edge Fowler flaps extending from 0.142 of the half span to the wing tip, see Figs. 5 and 6. The analysis was conducted in the RAE 11.5-by 8.5-ft low-speed wind tunnel at a nominal Mach number of 0.223, corresponding to a Reynolds number of 1.35e6 based on the mean wing chord. Furthermore, in the PDR analysis it was assumed that the wing was rigid and thus no wing deflection occurs. Below, in Fig. 7a the computed results are shown together with experimental values for C L max (Lovell 1977).
The maximum error between the computed and experimental results is 4.38% at a flap deflection of 10 • . A second test case has been performed to validate C L max of the Fokker 100 class wing, for which the geometry and flow parameters are described in Section 7. Using the PDR, the clean C L max was predicted to be 1.71 (See Fig. 7b). Compared to the actual value of 1.72 (Obert 2009), this is an error of 0.58%. The largest error is 9.17% at a flap deflection of 24 • . Considering the fact that in Valarezo and Chin (1994) the PDR has been validated against experimental data for numerous multi-element wing combinations up to flap angles of 40  degrees, the inaccuracy of the Fokker 100 C L max at high flap deflections may well be attributed to the fact that the exact flap geometry of the Fokker 100 wing was unavailable for this research.

Wing weight
The total wing weight is then computed using the following equation: where W FEM wingbox is the optimum wingbox weight computed using FEM. The factor 2 counts for non-optimum weight components. Elham and van Tooren (2016a) used a factor of 1.5 for the non-optimum weights. As reported in Elham and van Tooren (2016a) using the factor of 1.5 for non-optimum weights for A320 class wings produced remarkably accurate results. Initial investigation in the wing weight prediction of the Fokker 100 wing showed that this method underestimates the actual wing weight by 12%, which is likely due to the fact that the Fokker 100 design is a relatively old design and therefore far from optimal. After increasing the nonoptimal weight factor from 1.5 to 2, results agreed much better with actual wing weight data of the Fokker 100 wing (Paul 1993). The method developed by Torenbeek (1992) is used to compute the wing secondary weights. This method includes an empirical method for flap weight estimation. The primary and secondary wing weight of Fokker 100 is estimated to be 2853 kg and 1543 kg respectively using the mentioned method. It makes the total wing structural weight equal to 4369 kg. Comparing to the actual wing weight of Fokker 100 equal to 4343 kg, the weight estimation methods seems to be accurate enough for this research.

Airfield performance
Airfield performance is governed by the take-off and landing distances. Regulations for take-off and landing performance are presented in FAR 25, which defines the take-off distance as the ground covered from standstill to the point where the aircraft is at 35 ft above the ground. The landing distance is defined as the ground covered from a point where the aircraft is at 50 ft above the ground to standstill. Takeoff distance is the sum of the ground roll, rotation, transition and climbing distances. Airfield performance is approximated following the method from Raymer (2012) which has been adjusted for current regulations. This method takes the wing design parameters, stall speed and design weight, together with the computed drag and lift coefficient from the Q3D analysis at the characteristic take-off and landing velocities and calculates the respective take-off and landing distances: The ground run distance is computed from: Here V i is the initial velocity taken as 0 for the ground run, and V f is the velocity at the start of rotation which may be no less than 1.1 V s-1g according to Raymer. Coefficients K A and K T are defined as: HereT is the average thrust taken to be equal to the thrust at V = 1/ √ 2V f . The rolling friction coefficient μ r is taken to be 0.03. Weight is equal to MTOW and the density is taken to be 1.225kg/m 3 . For large aircraft, the rotation time to liftoff attitude may be assumed to be three seconds. Therefore s R can be approximated to be 3.3V s-1g . During transition, the aircraft accelerates from 1.1 V s-1g to 1.13 V s-1g . The average speed during transition is therefore 1.12 V s-1g . The average lift coefficient during transition may be assumed to be 90% of C L max | TO . The average load factor during transition is equal to 1.2, which gives a rotation arc radius of R = 0.205V s-1g . The climb angle at the end of transition is determined from: The distance traveled and altitude gained during transition is computed by: The value of h TR needs to be checked against the 35 ft obstacle height. If the obstacle height is cleared before the end of the transition segment, the following equation is to be used to compute the transition distance: Finally, unless the obstacle height is cleared during transition, the horizontal distance traveled during climb to clear the 35 ft obstacle height is found from: Landing distance is computed analogously to take-off distance, taking into account that approach speed V A must be at least 1.23 times higher than V s-1g , the approach angle should not be steeper than 3 deg and the touch down velocity V T D is assumed to be 1.15 times V s-1g according to Raymer. This results in an average flaring velocity V F L of 1.19 times V S 0 . The load factor during landing can be taken as 1.2 and the rolling friction coefficient due to deployed brakes during the ground run can be taken to be 10 times higher than during take-off. It should be noted that typically, the aircraft rolls free for 1 to 3 seconds before the pilot applies the brakes.
Validation of this method is performed for the Fokker 100 aircraft of which details are described in Section 7. The computed take-off and landing distance are compared to Fokker 100 performance data in Table 2.

Sensitivities Verification
Finally, the sensitivities computed by the presented analysis tool are verified through Finite Differencing (FD). In their paper, Elham and van Tooren (2016a) have already verified the sensitivities of several outputs of FEMWET. For this reason, only the sensitivities of landing distance with respect to a number of input variables are verified for a Fokker 100 class wing. From the research by Elham and van Tooren, it has been observed that while MSES computes sensitivities with a large number of digits, it only reports them up to 4 digits. The reported sensitivity of drag is therefore limited to an order of 10 − 3, which is not sufficient to perform verification through finite differencing. In order to still be able to verify the computed sensitivities, an empirical aerodynamic analysis tool has been developed, which does not rely on sensitivities computed by MSES but instead computes them using AD. The results of the sensitivity verification is shown below in Table 3.

Test case application
As a test case, the aerostructural optimization of a Fokker 100 class aircraft wing is considered. The aerostructural optimization is performed using SNOPT (Gill et al. 2005) and is formulated as follows: In the design vector X, the variables T, P and P f represent respectively the equivalent panel thickness of the wingbox, the planform geometry variables and the high-lift device planform geometry variables. The last two group of design variables are shown in Fig. 8.
The fourth group of design variables, G j , is used to define the wing airfoils' shapes. At 8 spanwise positions, the airfoil geometry is perturbed normal to its surface by a value ( n) which is determined based on basis functions g j , the Chebyshev polynomials, and the mode amplitudes G j as follows: Here s is the fractional arc length of the segment that the perturbation is applied to. In order to send a single set of shape design variables to the optimizer which applies to all analysis cases, the airfoil surface is split in sections as shown in Fig. 9. Through this method, the effect of perturbations on the cruise wing with respect to e.g. cruise drag can be combined with the effect of the same perturbation on C L max in landing configuration. A total of 9 shape variables are used to perturb respectively sections A-C and A-B and 1 shape variable is used for respectively sections C-G and B-G in order to prevent extreme discontinuities at intersections B and C. A total of 20 shape variables are thus defined per airfoil section and are active in both clean and highlift wing configuration analysis. It should be noted that as a result of splitting the airfoil sections, points B and C remain

Constraint Violation
Combined Sequential 1 Sequential 2 Fig. 11 History of wing aerostructural optimizations stationary during the optimization. Although it is possible to accommodate for horizontal and vertical perturbation of these points through the use of finite differencing, this is beyond the scope of the present research. The fifth group of design variables is used to perturb the flap's position using the mode amplitudes G t k (see Fig. 9). The mode amplitudes consist of two translational modes. G t 1 controls the horizontal translation of the flap and G t 2 the vertical translation. The third mode amplitude G t 3 controls the flap deflection.
The final group of variables are used to avoid unnecessary iterations for aeroelastic analysis. The optimization problem is subject to a number of constraints including constraints on structural failure and aileron effectiveness as described in Elham and van Tooren (2016a). In the same research, Elham and Van Tooren included a constraint on wing loading to take airfield performance into account. In the present research, the wing loading constraint is replaced by constraints on take-off and landing distance.
In total, a number of 835 inequality and 2 equality constraints are defined, and a total of 214 design variables are used for optimization (Table 4).
For aerodynamic analysis, it is assumed that the high-lift system consists of one continuous single slotted flap with a maximum deflection of 20 • . It is assumed that the takeoff is performed with flaps retracted as the Fokker 100 has been certified to take-off without flaps (Obert 2009). So the flap is optimized only for landing requirement. However the take-off performance was still taken into account since modifying the wing planform geometry affects the take-off performance. The wingbox internal structure is initiated by performing an aeroelastic optimization, aimed to minimize structural weight while satisfying failure constraints such as  Table 5 and a constraint on aileron effectiveness. The structural analysis is performed for the load cases listed in Table 5. Fatigue is simulated by limiting the stress in the wing box lower panel to 42% of the maximum allowable stress of the material in a 1.3g gust load case (Hürlimann et al. 2011). The aircraft rolling moment due to aileron deflection (L δ ) in the critical roll case is limited to be higher or equal to L δ of the initial aircraft.
The mission fuel weight (W fuel ) is computed based on the method of Roskam (2003). The required fuel use for cruise is computed using the Bréguet range equation, while statistical factors are used to determine the fuel weight of the remaining segments of the mission. The total aircraft drag is assumed to be the sum of the wing drag and the drag of the rest of the aircraft based on Fokker 100 aircraft data (Obert 2009). The drag of the aircraft minus wing is kept constant during the optimization. The aircraft range, cruise Mach number, altitude and engine parameters are determined based on aircraft data. The aircraft Maximum  take-off Weight (MTOW) is assumed to be equal to the payload weight, the aircraft fuel weight, the wing structural weight and a weight components that is called the rest weight. which is the operational empty weight minus the wing structural weight. The rest weight of the aircraft is computed from the Fokker 100 weight data and is kept constant during optimization.
The extended design structure matrix of this optimization is shown in Fig. 10.
In order to compare the proposed optimization procedure to conventional optimization schemes, two optimizations were performed for the Fokker 100 wing: -Wing A (Concurrent optimization): The wing planform and airfoil shape and the high lift device geometry are optimized simultaneously for minimizing the aircraft mission fuel weight and satisfying the field performance constraints, (7). -Wing B (Sequential optimization): Wing planform and airfoils are first optimized for minimum fuel weight (based on the cruise condition) without any airfield performance constraint. Then the high-lift devices are optimized for satisfying the field performance requirements.
The optimization histories of both wings A and B are shown in Fig. 11. Both optimizations have been performed using 8 parallel 3.50 GHz processors and 63 GB of RAM. The optimization of wing A finished in 12 hours after 18 iterations with a total number of 109 objective function evaluations. The initial optimization of Wing B finished in 7 hours after 15 iterations with a total number of 73 function evaluations. The subsequent high-lift device optimization finished in 12.5 hours after 4 iterations and 86 function evaluations. The computational cost breakdown of the tool is as follows: Wing aerostructural analysis is performed in 159.6 seconds, of which 133.5 seconds may be attributed to the coupled adjoint sensitivity analysis. The high-lift wing analysis, which includes the computation of C L max , takes 349.1 seconds to complete of which 268.7 seconds can be attributed to the coupled-adjoint analysis. In Fig. 12 the contours of the wingbox and fuel tank (both shown in blue) of wing A is shown. The front and rear spar position as function of local chord length were kept constant during the optimization. In Fig. 13 the optimized planforms of cases A and B are overlaid with the initial planform design. It can be seen that both optimized wings have a higher aspect ratio of respectively 10.65 and 9.74 compared to 8.26 of the initial wing. Quarter chord sweep angles of cases A and B were reduced to respectively 14.9 • and 16.9 • compared to the initial sweep angle of 18.5 • ( Table 6). The difference in the planforms originate from the different optimization schemes as a pure cruise wing design tends to have an increased sweep angle and larger wing loading to reduce pressure drag than a wing design which takes into account airfield performance. The fuel weight reduction of wing A is 9.65% and 9.20% for wing B while the wing weight is only marginally reduced by respectively 3.6% and 3.3% for wings A and B respectively. From Table 7 it can be deducted that the reduced fuel weight is primarily achieved through reducing the total drag by respectively 28.7% and 25.0% for wings A and B. The induced drag of respectively wings A and B is reduced by 24.4% and 14.1% by increasing the aspect ratio. The largest reduction has been in pressure drag which has seen a reduction of respectively 57.1% and 60.3%. This reduction in pressure drag is a result of the optimized airfoil shapes which can be seen in Fig. 14. Despite the increased normal Mach number due to the reduced sweep angle, the optimizer was able to reduce the strength of the shock wave and even at several spanwise positions completely remove the shock wave on the airfoils at cruise speed, see Fig. 15. The lift distribution on the initial as well as the optimized wings are plotted in Fig. 16.
Notably, the structural weight of both wings A and B has increased to a small degree. While a higher aspect ratio typically results in a larger structural weight near the root in order to withstand the increased bending moment, the weight penalty has been partially countered by reducing wing sweep and increasing wing flexibility. This increase in flexibility has been achieved by reducing the equivalent panel thickness of the outboard wing sections, which becomes visible in Fig. 17. The wing tip vertical and twist deformation under 1g load of the initial wing are respectively 0.39m and −0.84 • . For a 2.5g pull up these values are 1.12m and −2.24 • . The vertical and twist deformation of the tip of the wing A under 1g load are 0.94m and −1.29 • and for wing B 0.75m and −1.25 • . For a 2.5g pull up these deformations are 2.37m and −3.02 • for wing A and 1.93m and −2.99 • for wing B. The wing deformation under 2.5g pull up is shown graphically in (Fig. 18).
In Table 8 the flap designs for the initial and optimized designs are listed. Because wing B has a smaller wing span and wing area than wing A, it is expected to require a larger flap span to maintain airfield performance. Indeed, from Table 8 it follows that the flap span for the wing A is reduced by 11.8% and for wing B only by 7.23%. The flap weight of the wings A and B is also reduced by respectively 22.9% and 17.9%. Flap overlap (h f ), gap (g f ) and deflection angle (δ f ) have all increased in landing configuration in order to increase flap efficiency to counter the reduced flap span and increased wing loading as can be seen in Fig. 19. It is also   Table 9. It can be seen that both wings A and B have increased C L max in take-off and landing configuration as a result of both a reduced sweep angle and increased flap deflection. While both wings have a reduced wing area, the stall speed V s 1 g at MTOW and MLW remain similar to the initial design as a result of the increased C L max values.
While wing A meets both airfield performance requirements, wing B is not able to achieve take-off performance with flaps retracted. Although this seems to be an infeasible solution, it is a result of the optimization formulation, in which take-off is always performed with flaps retracted. So the wing planform was not design to satisfy the field performance and the flap optimization was not performed for take-off condition. When considering the design of wing B, it can be seen that its wing loading has increased due to the small wing area. While this is beneficial for cruise flight, it is undesirable for maintaining airfield performance.

Conclusion
An enhanced coupled-adjoint aerostructural analysis and optimization tool has been presented which enables the optimization of high-lift devices from the start of the design process. The semi-empirical Pressure Difference Rule has been implemented in an existing quasi-three-dimensional aerodynamic analysis and coupled to the structural solver FEMWET to compute the wing C L max taking into account aeroelastic effects. The coupled system is solved using the Newton method.
The modified aerostructural tool is able to compute the derivatives of the outputs with respect to the inputs using a combination of the chain rule of differentiation, automatic differentiation and coupled-adjoint method. Wing weight is computed using the empirical method of Torenbeek which enables optimization taking into account the specific weight of high-lift devices. Airfield performance is determined using the method of Raymer which includes lift and drag terms at take-off and landing. Validation of the modifications showed good levels of accuracy for L/D, C L max , wing weight and airfield performance. The tool was then used for gradient based aerostructural wing and high-lift system optimization for a Fokker 100 class wing. Design variables include flap span and settings, wing planform, airfoil geometry and wingbox structure. The optimization was performed for minimizing the fuel weight, while satisfying constraints on structural failure, C L max in take-off and landing configuration and the roll requirement. The flaps were retracted in the take-off configuration as the Fokker 100 is certified to perform take-off with flaps retracted. Two types of optimizations were performed. The proposed concurrent optimization scheme where the highlift system was optimized from the start of the process and a more conventional sequential optimization, in which the planform was first optimized for cruise performance after which the high-lift system was sized to minimize fuel weight for the fixed optimized planform taking into account airfield performance.
The concurrent optimization resulted in a fuel weight reduction of 9.65%, while sequential optimization reduced fuel weight by 9.20%. The reduced fuel weight was attributed to a reduction in pressure drag resulting from modified airfoil shapes, reducing shock waves over the wing. The optimized wings both have an increased aspect ratio and reduced sweep angle, which resulted in a reduction of induced drag. To counter the weight penalty due to these modifications, the structural stiffness was reduced near the wing tip. The optimizer reduced flap weight of both optimized wings by respectively 22.9% and 17.9% by reducing flap span. It can be concluded that the proposed method of combining the optimization of high-lift and cruise wing design is promising and provides ample opportunities for more research.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http:// creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.