Multi-fidelity wing aerostructural optimization using a trust region filter-SQP algorithm

A trust region filter-SQP method is used for wing multi-fidelity aerostructural optimization. Filter method eliminates the need for a penalty function, and subsequently a penalty parameter. Besides, it can easily be modified to be used for multi-fidelity optimization. A low fidelity aerostructural analysis tool is presented, that computes the drag, weight and structural deformation of lifting surfaces as well as their sensitivities with respect to the design variables using analytical methods. That tool is used for a mono-fidelity wing aerostructral optimization using a trust region filter-SQP method. In addition to that, a multi-fidelity aerostructural optimization has been performed, using a higher fidelity CFD code to calibrate the results of the lower fidelity model. In that case, the lower fidelity tool is used to compute the objective function, constraints and their derivatives to construct the quadratic programming subproblem. The high fidelity model is used to compute the objective function and the constraints used to generate the filter. The results of the high fidelity analysis are also used to calibrate the results of the lower fidelity tool during the optimization. This method is applied to optimize the wing of an A320 like aircraft for minimum fuel burn. The results showed about 9 % reduction in the aircraft mission fuel burn.


Introduction
According to industry criteria for aircraft design, the drag prediction accuracy using numerical methods should be within one drag count (one ten thousandth of the drag coefficient) (van Dam 2003). To confirm a need for such a level of accuracy, Meredith (1993) showed that changing the drag coefficient by 0.0001 has an effect equivalent to the weight of one passenger for a long-haul aircraft. Similarly, one percent error in wing structural weight estimation of the same class aircraft is equal to the weight of 4 to 5 passengers. The need for such a high level of accuracy forces the designers to use high fidelity, physics based analysis for aerostructural analysis, design and optimization of aircraft. The traditional design methods based on empirical, statistic based methods do not satisfy the required level of accuracy.
On the other hand, using high fidelity methods for aerostructural optimization requires massive computational power Kennedy and Martins 2014), that is a serious barrier against using high fidelity aerostructural optimization in early design stages. As a solution multi-fidelity optimization techniques are used to keep the level of accuracy similar to the results of the high fidelity analysis methods, while reduce the computational cost of the optimization. Alexandrov et al. (2001) presented a model management framework for multi-fidelity aerodynamic shape optimization of lifting surfaces based on a trust region algorithm. In that model a lower fidelity tool is used for the optimization, while a higher fidelity tool is occasionally, but systematically, called to monitor the optimization process. March and Willcox (2010) suggested a multi-fidelity optimization framework based on a trust region algorithm, in which the gradient of the objective function is computed using the low fidelity model, but the algorithm is provably convergent to the solution of the high fidelity model. The same algorithm is used by Elham (2015) for aerodynamic shape optimization of lifting surfaces, where an adjoint quasi-three-dimensional (Q3D) model is used for prediction of the wing drag and its derivatives, and a three-dimensional CFD tool is used to calibrate the results of the Q3D model. Applications of multi-fidelity model-based strategy for aerodynamic and aerostructural optimization are presented by Rajnarayan et al. (2008) and Berci et al. (2014).
Besides the framework for model management in a multifidelity optimization, the choice of a proper algorithm for numerical optimization is important. Aerostructural optimization of lifting surfaces, or the whole aircraft in general, involves hundreds to thousands of design variables, and tens to hundreds of constraints. Besides, computing the objective function and the constraints required execution of CFD and FEM analysis, which takes considerable amount of time. Therefore the gradient based optimization algorithms are the most efficient algorithms for solving such problems (Martins et al. 2005). In an optimization using a gradient based algorithm, in order to achieve a quadratic rate of convergence, an underlying Newton iteration is required, which is the basics of Sequential Quadratic Programming (SQP). The SQP algorithms are based on iteratively solving a quadratic model of the objective function and linear models of the constraints. The SQP approach has been used in both line search and trust region frameworks (Nocedal and Wright 2000). Many of the SQP methods use a penalty function, which combines the objective function and the constraints. An 1 penalty function and the augmented Lagrangian function are popular choices for the penalty function. Selection of a proper penalty function, and the method used for updating the penalty parameter is a challenge (Nocedal and Wright 2000).  proposed a method to eliminate the need for a penalty function in an SQP algorithm. They suggested a filter that rejects the unacceptable solutions. In that so-called trust region filter-SQP algorithm, no penalty function is constructed and the filter is applied to the objective function and the (norm of the) constraints separately.
In this research a modified version of the filter method proposed by Fletcher and Leyffer is applied to multi-fidelity aerostructural optimization of lifting surfaces. In the next section, the aerostructural analysis tool used in this research is briefly explained. Then in Section 3, the trust region filter-SQP method is discussed in detail. Eventually in Section 4 the filter method is applied to both a mono-fidelity and a multi-fidelity wing aerostructural optimization.

Aerostructural analysis
The FEMWET tool presented by Elham and van Tooren (2016b) is used for wing aerostructural analysis. FEMWET is based on a Q3D aerodynamic analysis method, which is combined with a finite beam element model of the structure. In the Q3D approach an inviscid vortex lattice method (VLM) is combined with a viscous 2D airfoil analysis code for prediction of the wing total viscous drag. The idea behind the Q3D approach is to avoid using a high fidelity 3D CFD solver, but still predict drag with the same level of accuracy. In the proposed method the wing drag is decomposed into the induced drag and the parasite drag. To compute the wing total drag, first a VLM is executed to compute the lift distribution over the wing as well as the wing induced drag using a Trefftz plane analysis. The VLM code is connected to a FEM to automatically deform the VLM mesh based on the wing structural deformation. Then several sections along the span are analyzed using a 2D airfoil analysis CFD code. The parasite drag is computed based on the pressure and friction drag of the 2D sections. In order to perform the 2D analysis, several corrections, including the corrections due to sweep, induced angle of attack and wing structural deformation are applied to the section geometry as well as the flow conditions. Details of this method are presented in reference (Elham and van Tooren 2016b).
The wing box structure is modeled using four equivalent panels: the upper panel, including the wing box upper skin, stringers and spars caps; the lower panel including the lower skin,stringers and spar caps; and two vertical panels including the front and the rear spars webs. The stiffened panel efficiency method presented by Niu (1997) is used to model the effect of stringers on the thickness of the equivalent panels. For finite element analysis of the wing box structure an equivalent Timoshenko beam is placed at the shear center of the wing box, see Fig. 1. Using this FEM model different failure criteria, referred to tension, compression, Euler buckling and shear buckling in several wing box elements, located in (y e , z e ) distance from the shear center (see Fig. 1), are calculated. For more details one can refer to Elham and van Tooren (2016b). When the Q3D aerodynamic solver is combined with the finite beam element model, four governing equations appear as follows: The first equation is the governing equation of the VLM and the second equation is the governing equation of the FEM. The third equation in fact indicates that the lift should be equal to the design weight multiplied by the load factor. The last equation is required to guarantee that the lift predicted by the VLM is the same as the lift computed by 2D section analysis at effective angle of attack. The effective angle of attack is the angle of attack that a 2D local section feels. Therefore the local downwash angle is required to compute the effective angle of attack. In such a system four sets of state variables are presented: the strengths of vortices in the VLM ( ), the displacements in FEM (U ), the global angle of attack (α) and the local downwash angles at each 2D section (α i ). For a given vector of design variables, X, the system of governing equations are solved using the Newton method: To perform the Newton iteration, the matrix of the partial derivatives J is required. All the elements of that matrix are computed by a combined use of analytical methods and Automatic Differentiation (AD). The Matlab toolbox Intlab (Rump 1999) is used for AD. More details of the sensitivity analysis are presented in Elham and van Tooren (2016b). In order to compute the sensitivity of the output (e.g. wing drag or structural failure loads) the coupled adjoint method ) is used, where the total derivative of any function of interest I , with respect to a design variable x is presented as follows: in which λ = [λ 1 λ 2 λ 3 λ 4 ] T is the adjoint vector and computed using the following equation: The FEMWET tool is also able to compute the aileron effectiveness and its sensitivity with respect to the design variables. The aileron effectiveness is defined as L δ elastic /L δ rigid . The parameter L δ is the derivative of the wing rolling moment w.r.t the aileron deflection angle. This parameter is an important requirement for aircraft performance and is strongly affected by the wing stiffness. Low wing (mainly torsional) stiffness may result in a poor roll performance or even aileron reversal. Therefore for a wing aerostructural optimization a constraint on the aileron effectiveness seems necessary.
Elham and van Tooren (2016b) performed some analysis and aerostructural optimization to verify the results of the FEMWET tool. In order to verify the accuracy of the wing drag prediction, the results of Q3D aerodynamic solver were compared to the results of a higher fidelity CFD code called MATRICS-V (van der Wees et al. 1993). The MATRICS-V flow solver is based on fully conservative full potential outer flow in quasi-simultaneous interaction with an integral boundary layer method on the wing. The MATRICS-V tool was developed by NLR and has been validated using wind tunnel test as well as the flight test results for Fokker 100 aircraft, see Figs. 2 and 3. Therefore the drag of the Fokker 100 wing predicted by the Q3D solver was compared to the drag predicted by MATRICS-V, see Fig. 4. From this figure one can observe a high level of accuracy for drag prediction using the Q3D solver.
In order to validate the wing weight and structural deformation predicted by FEMWET, Elham and van Tooren (2016b) performed an aeroelastic wingbox optimization of Airbus A320-200 wing. The optimization included five different load cases to evaluate the structural failure with respect to tensile and compressive loads, buckling, fatigue and aileron effectiveness. The total wing weight predicted by FEMWET is equal to 8791kg, which is very close to the actual wing weight of A320-200 equal to 8801kg (Obert 2009). The wing twist deformation predicted by FEMWET was compared to the actual wing twist deformation under 1g load (obtained from Obert (2009)). Figure 5 shows the the wing actual twist distribution (including the wing incidence angle at the wing root) compared to the computed values. The results show an acceptable level of accuracy for the tool.

Trust region filter-SQP method
A general optimization problem is defined as follows: Only equality constraints are considered here, however inequality constraints can be defined as equality constraints using slack variables and adding simple bounds, see Nocedal and Wright (2000). Using an SQP approach, (8) is solved by solving the following quadratic problem (QP) iteratively: where H is the Hessian matrix of the Lagrangian function, g and A are the gradient of the objective function and the constraints respectively. Most of the available SQP algorithms use the solution of the QP as a search direction, then find a step length that minimizes a one-dimensional problem, which results in a sufficient decrease of a penalty function. The penalty function combines the objective function and the constraints in a single function. One-dimensional optimization algorithms, (such as polynomial interpolation) (Vanderplaats 2007) or line search algorithms (Nocedal and Wright 2000) are used to find the step length. An alternative to this approach is the use of trust region algorithms (Conn et al. 2000). Trust region methods define a region around the current point, where the approximations of the objective function and constraints are trusted. The radius of the trust region plays the role of the step length, so the trust region algorithms find the search direction and the step length simultaneously. However the need for a penalty function is still there. As mentioned earlier, definition of a penalty function and consequently a method for updating an associated penalty parameter is a challenge. Some aspects of the difficulties associated with the choice of the penalty function and the penalty parameter are discussed in . The filter method presented by  eliminates the need for a penalty function in an SQP algorithm. The concept of the filter is based on the two goals of a constrained optimization problem: minimizing the objective function and minimizing the constraint violation. So in the method proposed by Fletcher and Leyffer a filter is used to only accept the solutions that are not dominated by the Pareto front between the objective function and the constraint violation, see Fig. 6. If ϑ is the 2-norm of the (equality) constraints, the pair (f 1 , ϑ 1 ) is said to dominate (f 2 , ϑ 2 ) if and only if both f 1 ≤ f 2 and ϑ 1 ≤ ϑ 2 . Defining F as the filter, that includes a set of pairs (f j , ϑ j ) such that no pair dominate any other, a pair (f, ϑ) is acceptable to F if it is not dominated by any pair in the filter. In a trust region filter-SQP method, the QP shown in (9) is solved within a trust region, then the solution is checked by the filter. If the filter rejects the solution, the radius of the trust region is reduced, see  for more details. A common problem in trust region methods, is that the QP may have no feasible solution if the radius of the trust region is small. In such cases a constraint restoration is required. The idea of constraint restoration is to minimize ϑ(x) starting from the current iteration. A basic trust region filter-SQP algorithm in shown in Algorithm 1.
In order to prove the convergence of the trust region filter-SQP algorithm, a small envelope is required around the current filter, in which no point is accepted. This envelope in fact enforces a sufficient decrease in the objective function and the constraint. According to  a point is acceptable to the filter if: either f ≤ f j − γ ϑ j or ϑ ≤ βϑ j for all j ∈ F (10) The proof of convergence of such an algorithm is given in . A more refined trust region filter-SQP method is presented by Conn et al. (2000), although in that algorithm a composite step optimization is used, where first in a normal step the norm of the constraints is minimized within the trust region and then in a tangential step the objective function is reduced. In this research the algorithm presented in Conn et al. (2000) is modified to use the original SQP method presented by  instead of a composite step optimization. This algorithm is presented in Algorithm 2. Conn et al. (2000) suggested the following values for the constants in Algorithm 2: γ 0 = 0.5 γ 1 = 2 η 1 = 0.01 η 2 = 0.9 γ ϑ = 10 −4 k ϑ = 10 −4 In some of the trust region methods, as in the one suggested by Conn et al. (2000), a 2 trust region subproblem is used. However in the filter method suggested by  a ∞ trust region subproblem is used to ensure that the subproblem remains tractable as a QP. In this research instead of 2 or ∞ , the trust region is defined to keep the absolute value of s lower than the trust region radius, see Algorithm 2. It has an advantage when the algorithm is applied to a wing aerostructural optimization. This advantage is explained in Section 4. In such an approach the trust region for a 2D case is a rectangle instead of a circle for 2 or a square for ∞ .
To check the algorithm an analytical optimization test case is used as follows:

Problem formulation
An A320 like aircraft wing is considered as a test case, see Fig. 8. The aircraft mission fuel weight is considered as the objective function. It is computed using the method presented by Roskam (1986). In that method the fuel weight of the cruise phase of the flight is computed using the Breguet method, and the fuel weights of the other segments of the mission, e.g. take-off, climb, etc. are computed using some statistical factors. In order to use the Breguet equation, the aircraft lift over drag ratio during the cruise is required. The lift and drag of the elastic wing during the cruise is computed using the FEMWET for an average aircraft all up weight equal to W ave = MT OW × (MT OW − W f uel ) as suggested by Torenbeek (2013). The drag of the other aircraft components such as fuselage, tail etc. is assumed to be constant. The design variables are categorized into four groups. The first group includes the design variables describing the wing planform geometry. Six design variables are used for that purpose: the wing root chord C r , span b, taper ratio λ, leading edge sweep angle , twist angle at kink k and twist angle at tip t . The wing aileron geometry was defined using three parameters. The start and end position of the aileron was fixed at 75 % and 95 % of the wing span. The aileron chord was fixed as 25 % of the local wing chord. The second group of design variables defines the wing airfoil shape. In order to reduce the number of design variables and guarantee a smooth shape for the airfoil, the airfoil shape is parametrized using the Chebychev polynomials. The Chebychev polynomials, g i , are defined as: Using the Chebychev polynomials the original airfoil shape is perturbed as: where n is the airfoil perturbation normal to its current surface, s is the fractional arc length of each side of the airfoil and G j are the mode amplitudes, that are defined as design variables. 160 design variables are used to control the airfoils shapes at 8 wing spanwise positions. The third group includes the design variables defining the wing box structure. The thickness of the wing box four equivalent panels in 10 spanwise positions are defined as design variables, 40 in total. As mentioned earlier the aircraft fuel weight is defined as the objective function. The fuel weight is a function of the aircraft total weight, which is a function of the fuel weight. Also for wing structural analysis the aircraft MTOW is required which is a function of the fuel weight and wing structural weight. So in order to avoid any iteration during the optimization two surrogate variables for aircraft fuel weight and maximum take-off weight are added to the design vector as the fourth group of the design variables. The aerostructural optimization is subject to several constraints. A series of constraints are used to avoid any structural failure. These constraints are define based on different load cases, which are determined using the aircraft flight envelope. Reference (Dillinger 2014) shows the flight envelope of the A320 aircraft. Using those data six load cases are selected for aircraft structural analysis and optimization as shown in Table 1. Load cases number 1 to 3 are used for analyzing the wing box failure. Depending on the load case (positive or negative load factor) the upper and lower panels can be under tension or compression. and the spars webs are  History of the analytical test case optimization experiencing shear loads. All these different scenarios are taken into account. Load case number 4 is used to simulate fatigue in the wing lower panel and the load case number 5 is used to simulate the aileron effectiveness. Load case 6 is used for wing drag prediction during the cruise.
In order to reduce the total number of the constraints on structural failure, the Kresselmeier-Steinhauser (KS) function (Kreisselmeier and Steinhauser 1980) is used to aggregate the constraints, as follows: where g max is the maximum of all constraints. The parameter ρ KS has been set to 50 as suggested by Poon and Martins (2007). All the failure constraints due to the 5 load cases are aggregated into 22 constraints using the KS function. A constraint is defined to keep the aircraft roll Fig. 8 Planform and wing box dimensions moment due to aileron deflection (L δ = dL/dδ) higher or equal to the L δ of the original wing. Another constraint is used to keep the wing loading lower or equal to the wing loading of the initial wing. Finally two consistency constraints are defined for the two surrogate design variables. The wing aerostructural optimization is formulated as follows: min W * f uel (X) X = [C r , b, λ, , kink , tip , G i , t u j , t l j , t f s j , t rs j ,

Aerostructural optimization using the trust region filter-SQP method
In the first step a mono-fidelity optimization was performed.
In that approach the aerostructural analysis method presented in Section 2 is used to predict the elastic wing drag and deformation. All the inequality constraints in (15) are changed to equality constraints by using slack variables. In total 24 slack variablesare added to the design vector. Besides, a series of bounds are defined to keep the values of these variables positive or zero.
In the second step a multi-fidelity optimization has been performed. March and Willcox (2012) modified the composite step filter method presented by Conn et al. (2000) to be used in a multi-fidelity optimization. The same approach as suggested by March and Willcox is used here to modify Algorithm 2 to be used for multi-fidelity optimization. In the modified algorithm, the Q3D aerodynamic analysis (connected with the FEM) is used as the low fidelity model and the MATRICS-V CFD code is used as the high fidelity tool. MATRICS-V is a 3D CFD code, which provides more accurate results than the Q3D method, but the computational time of running MATERICS-V is higher than Q3D. Besides, no analytical sensitivity analysis method is implemented in MATRICS-V. The MATRICS-V has been used for aerodynamic optimization using finite differencing for sensitivity analysis (Elham and van Tooren 2014) with limited number of design variables. However increasing the number of design variables and coupling the aerodynamic solver with a FEM for fully coupled aerostructural optimization, makes the use of finite differencing almost impossible. Therefore in the current research the low-fidelity aerostructural analysis tool is used to generate the TRQP, since that tool can compute the required sensitivities using analytical methods. Then the filter is applied based on the results of the MATRICS-V code. The drag predicted using the low fidelity model is calibrated using the results of the high fidelity model. Three calibration factors are defined for three different drag components (the induced drag, C D i , the pressure drag, C D p ,and the friction drag, C D f ) as follows: After each iteration, if the new point is acceptable to the filter, the calibration factors are updated and the next iteration is performed. This guarantees that at the end of the optimization the drag predicted by the low fidelity model is the same as the drag predicted by the high fidelity model. This procedure is summarized in Algorithm 4. Since a surrogate variable is used for the aircraft fuel weight, which is the objective function, the variable level of fidelity does not affect the objective function. The value of the objective function is always the value of the surrogate fuel weight in the design vector. However the actual fuel weight is computed inside the constraint function in order to generate the consistency constraint on the fuel weight. Therefore two different values of the consistency constraint on the fuel weight and also the MTOW (which is a function of the fuel weight) are computed, one using the drag computed by the low fidelity model and one by using the drag predicted by the high fidelity model. The trust region filter-SQP method algorithm used for such a multi-fidelity optimization is shown in Algorithm 3.
As mentioned before in the TRQP, instead of 2 or ∞ , a bound around the absolute value of s is defined. So in TRQP is a vector with the same length as s. It helps to define different values of for different design variables. As one can see in (15) different groups of design variables are used for a wing aerostructural optimization. Although all the design variables are normalized to have the same order of magnitude, defining the same trust region radius  Fig. 11 The shape of the airfoils is several spanwise positions for all of them does not seem logical. For example assuming the trust region radius allows 10 percent change in the design variables ( = 0.1 for a design vector normalized with the initial values of the design variables). The effect of 10 % change in the thickness of the wing box skin at wing tip on the aircraft fuel weight is not the same as the effect of 10 % change in the aircraft MTOW on the aircraft fuel weight. By defining as a vector, some design variables are allowed to have a larger change and some are more tightly limited. Selecting the proper values of for different design variables is a challenge by itself. As mentioned, the value of should be determined based on the sensitivity of the objective function and constraints with respect to the design variables. For those variables that have higher influence on the objective function a lower value for should be used. In this case the design variables are categorized into four groups. The first group includes the design variables defining the wingbox structure. The second group includes the variables defining the airfoils geometry. The variables defining the wing planform geometry are categorized into the third group. And finally the variables defining the aircraft MTOW and the aircraft fuel weight are in the fourth  group. Defining the value of , the trust region radios of the first and second group of the variables are set equal to , the trust region radios of the third group is set to /2, and the trust region radios of the fourth group is set to /3.
The history of the both mono-fidelity and multi-fidelity optimization is shown in Fig. 9. Figure 10 shows the planform of the initial wing compared to the planform of the optimized wing using both the mono-fidelity and the multifidelity optimizations. The shape of the airfoils in several spanwise position and the pressure distribution on those airfoils are shown in Figs. 11 and 12 respectively. The characteristics of the optimized aircraft are summarized in Table 2.
From Table 2 one can observe that the multi-fidelity optimization resulted in slightly less fuel weight reduction (about 0.7 % less than mono-fidelity optimization). The drag of the wing optimized using the multi-fuidelity approach is slightly higher than the drag of the wing optimized using only the low fidelity model. It can be explained  Fig. 4, showing that the low fidelity model slightly underestimates the drag comparing to the high fidelity model. The wing optimized using the multi-fidelity approach has slightly higher sweep than the one optimized using the mono-fidelity method, therefore it has slightly higher wing structural weight. In general the results of the multi-fidelity optimization are more realistic, since a more accurate drag analysis is used. However the results of the mono-fidelity optimization are quite similar to the results of the multi-fidelity optimization, that indicates a good level of accuracy of the low fidelity model. It should also be noted that the mono-fidelity optimization was about 4 times faster than the multi-fidelity optimization.
In both cases the optimizer moved toward a larger wing span and a lower leading edge sweep. The larger span results in a lower induced drag, but a higher wing structural weight. Lower sweep, on the other hand, results in a lower structural weight, but may increase the wave drag. However the optimizer managed to modify the airfoil shapes to eliminate the shock wave from the surface of the optimized wings, see Fig. 12.
The constraint on the aileron effectiveness is usually an active constraint in wing aeroelastic optimization. In fact to achieve higher aileron effectiveness, and consequently a higher value of L δ , a higher wing stiffness (mainly torsional stiffness) is required for a given wing and aileron geometry. Increasing the torsional stiffness results in a higher structural weight. The study of Elham and van Tooren (2016a) showed that the wing structural weight increases quadratically by increasing the required value for aileron effectiveness. The optimizer in this research, both in the mono-fidelity and multi-fidelity cases, moved toward a more flexible wing to reduce the wing structural weight. The initial wing has a tip vertical deflection of 1.48m and tip twist of -3.8 degrees under 2.5g pull up load, while the optimized wing (using multi-fidelity method) has a tip vertical deflection of 1.77m and tip twist of -3.9 degrees. The aircraft roll requirement was satisfied by increasing the aileron arm and the aileron surface (both were resulted from a larger span). The larger aileron area and arm allowed to keep the value of L δ higher than the required (4.04 × 10 4 for the optimized wing vs 3.80 × 10 4 for the initial wing) with a lower value of the aileron effectiveness (0.43 for the optimized wing vs 0.53 for the initial wing), that resulted from a lower wing stiffness.

Conclusions
In this research a trust region filter-SQP method is used for wing multi-fidelity aerostructural optimization. This algorithm allows to combine a lower fidelity model, that predicts the sensitivity of the objective function and the constraints, with a higher fidelity model, that is more accurate but more expensive to be executed. The low fidelity model is used to generate the TRQP subproblem. The high fidelity model is used to generate the filter and also to calibrate the results of the low fidelity model. Using that approach a high fidelity CFD tool that does not provide the sensitivities, can be used for a gradient based optimization. In addition to that, the aerodynamic solver is coupled with a structural solver for a fully coupled aerostructural optimization.
This research showed that a high fidelity aerodynamic solver can be used in an optimization without a need for high performance computer clusters. This reduces the cost and time of the optimization and makes high fidelity optimization feasible using ordinary computers.
A mono-fidelity as well as a multi-fidelity wing aerostructural optimization have been performed using the proposed algorithm. The results showed about 9 % reduction in the aircraft fuel weight. The optimizer found the optimum planform shape, airfoil shape as well as the wing box structure to achieve that amount of reduction in the fuel weight.
In this study only a high fidelity aerodynamic solver was combined with the low fidelity aerostructural analysis tool. As the next step a high fidelity FEM can be combined with the low fidelity tool as well. In that approach both the aerodynamic and structural analysis can be performed using two different levels of fidelity. There will be no need for modifying the proposed algorithm. The same algorithm can be used, but this time the constraints related to the wing structure will be computed using two different tools.
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.