Next Article in Journal
Adaptive Robust Time-Varying Formation Control of Quadrotors under Switching Topologies: Theory and Experiment
Next Article in Special Issue
Extended Hierarchical Kriging Method for Aerodynamic Model Generation Incorporating Multiple Low-Fidelity Datasets
Previous Article in Journal
A Comparison of the Damage Tolerance of AA7075-T6, AA2024-T3, and Boeing Space, Intelligence, and Weapons Systems AM-Built LPBF Scalmalloy
Previous Article in Special Issue
Aero-Engine Preliminary Design Optimization and Operability Studies Supported by a Compressor Mean-Line Design Module
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Efficient Global Aerodynamic Shape Optimization of a Full Aircraft Configuration Considering Trimming

1
National Key Laboratory of Science and Technology on Aerodynamic Design and Research, Xi’an 710072, China
2
Institute of Aerodynamic and Multidisciplinary Design Optimization, School of Aeronautics, Northwestern Polytechnical University, Xi’an 710072, China
*
Authors to whom correspondence should be addressed.
Aerospace 2023, 10(8), 734; https://doi.org/10.3390/aerospace10080734
Submission received: 19 July 2023 / Revised: 16 August 2023 / Accepted: 19 August 2023 / Published: 21 August 2023
(This article belongs to the Special Issue Aerodynamic and Multidisciplinary Design Optimization)

Abstract

:
Most existing aerodynamic shape optimization (ASO) studies do not take the balanced pitching moment into account and thus the optimized configuration has to be trimmed to ensure zero pitching moment, which causes additional drag and reduces the benefit of ASO remarkably. This article proposes an efficient global ASO method that directly enforces a zero pitching moment constraint. A free-form deformation (FFD) parameterization combing Laplacian smoothing method is implemented to parameterize a full aircraft configuration and ensure sufficiently smooth aerodynamic shapes. Reynolds-averaged Navier–Stokes (RANS) equations are solved to simulate transonic viscous flows. A surrogate-based multi-round optimization strategy is used to drive ASO towards the global optimum. To verify the effectiveness of the proposed method, we adopt two design optimization strategies for the NASA Common Research Model (CRM) wing–body–tail configuration. The first strategy is to optimize the configuration without considering balance of pitching moment, and then manually trim the optimized configuration by deflecting the horizontal tail. The second one is to directly enforce the zero pitching moment constraint in the optimization model and take the deflection angle of the horizontal tail as an additional design variable. Results show that: (1) for the first strategy, about 4-count drag-reducing benefits would be lost when manually trimming the optimal configuration; (2) the second strategy can achieve 3.2-count more drag-reducing benefits than the first strategy; (3) compared with gradient-based optimization (GBO), surrogate-based optimization (SBO) is more efficient than GBO for ASO problems with around 80 design variables, and the benefit of ASO achieved by SBO is comparable to that obtained by GBO.

1. Introduction

Following the continuous advancement of computer performance and numerical computation methods, computational fluid dynamics (CFD) has been developed to an extent that it has become an important tool to predict the aerodynamic performance of a full aircraft configuration and replace expensive wind tunnel experiments partially. As a result, CFD-driven aerodynamic shape optimization (ASO) has partly replaced the traditional “cut and try” method and inverse design method, and is becoming the important way to improve the aerodynamic performance of an aircraft. However, it is still difficult to perform ASO of complex aircraft configurations such as the wing–body–tail configuration, since high-fidelity CFD simulation methods such as the Reynolds-averaged Navier–Stokes (RANS) equations are computationally expensive, and a large number of design variables and constraints can cause the optimization algorithm to suffer from numerical difficulties.
Surrogate-based optimization (SBO) has become one of the most attractive methods for ASO with high-fidelity flow simulations. The solution of expensive cost functions can be approximated by cheap-to-evaluate surrogate models to drive the addition and evaluation of new sample point(s) towards the global optimum [1,2,3]. SBO requires much less expensive function evaluations to reach the global optimum compared with a genetic algorithm (GA), especially when expensive high-fidelity CFD is employed [4,5]. It is free of gradient information and can be used to solve problems with highly nonlinear and multi-modal design space [6,7,8]. Therefore, SBO has gained much attention in the area of engineering design optimization. Despite the popularity of SBO, it is still suffering from the difficulty associated with the “curse of dimensionality”, which means that the computational cost grows exponentially as the number of design variables increases. To improve the SBO efficiency, some advanced technologies have been proposed, such as evaluating the responses of the samples in parallel [9], variable-fidelity surrogate modeling [10,11,12,13,14,15,16], surrogate modeling incorporating cheap gradients [17,18,19,20], deep-learning-based optimal sampling [21], and applying dimension reduction [22,23] and constraint aggregation during optimization [24,25].
However, the design optimization of complex aircraft configurations such as the wing–body–tail configuration brings a large challenge to SBO. First, both design of experiment (DoE) and sub-optimization within the framework of SBO are pure numerical methods designed to generate samples without considering feasibility. Thus, many samples associated with unsmooth aerodynamic shapes would be generated via DoE and sub-optimization during optimization, which results in an unpractical optimal aerodynamic shape [22]. Second, the size of the design space directly affects the optimization results, and an inappropriate design space makes it difficult for SBO to obtain the global optimal solution [26,27]. Third, although fast design can be achieved by using a low-fidelity CFD solver in the optimization, the benefit of drag reduction cannot be realized completely when verifying the optimal solution by the high-fidelity CFD simulation such as solving the RANS equations [28]. Finally, it is unsuitable to design the wing of a complex aircraft configuration only since the benefit of drag reduction achieved would decrease significantly due to the trim drag induced by trimming the full aircraft configuration [29]. Recently, some design works considering trimming have been performed via the gradient-based optimization (GBO) [30,31,32,33,34]; however, the adjoint equations for solving the derivatives of the added constraints, such as the zero pitching moment constraint, must be derived, and solving the added adjoint equations can be expected to increase the computational cost of GBO significantly.
This article aims to propose an efficient global ASO method using surrogate models to design complex aircraft configurations, in which trimming is considered by adding a zero pitching moment constraint during optimization. The multi-round optimization strategy is used to drive SBO towards the global optimum by adjusting the design space in each round according to the optimum from the previous round. The Laplacian smoothing method is conducted to avoid obtaining unsmooth optimal aerodynamic shapes. To validate the proposed method, two optimization strategies are applied to design the NASA CRM wing–body–tail configuration. The first strategy is to optimize the wing only and then manually trim the optimal configuration by deflecting the horizontal tail as a solid body. The second strategy is to consider trimming during optimization. Then, all optimization results obtained via SBO are compared with those obtained via GBO. Two goals are involved in this work; one is to validate SBO for solving ASO of a complex aircraft configuration, since there is hardly any research focus on high-fidelity surrogate-based ASO of such a complex aircraft configuration. The other one is to compare the optimization results achieved through different optimization strategies, thus providing some useful guidance for future work.
This article proceeds with Section 2, where the numerical tools employed in this study are described. Section 3 illustrates the problem formulation, while Section 4 presents ASO of the NASA CRM wing–body–tail configuration and corresponding results. In Section 5, conclusions are presented.

2. Methodology

2.1. Geometric Parametrization

The free-form deformation (FFD) method [35] based on the non-uniform rational B-spline (NURBS) is applied to parametrize the NASA CRM wing–body–tail configuration, which can change the geometry by perturbing the FFD control point. As shown in Figure 1, the FFD volume has a total of 25 FFD volumetric blocks, one FFD volumetric block with 90 control points is used to parametrize the wing, and another one with 66 control points is used to parametrize the horizontal tail.

2.2. Laplacian Smoothing of Aerodynamic Shape

In a typical progression of surrogate-based ASO, a DoE method such as Latin hypercube sampling (LHS) [36] is used to generate initial sample points, and the sub-optimization defined by a certain infill-sampling criteria is used to select new samples in each SBO iteration. However, both DoE and sub-optimization are pure numerical methods designed to generate samples without considering feasibility. As a result, an increasing number of samples associated with unsmooth aerodynamic shapes are generated during optimization, which results in an unpractical optimal aerodynamic shape.
To address this issue, the Laplacian smoothing method is introduced in SBO to smooth the aerodynamic shapes corresponding to the samples generated by DoE and sub-optimization. The Laplacian smoothing method is defined as
y i = θ y i + 1 θ n n e i g h b o r j = 1 n n e i g h b o r y j ,
where yi denotes the y-direction deformation of the i-th FFD control point used to parameterize the wing and where θ is used as a weight coefficient, which can influence the smoothing operation speed. Li et al. found that a practical aerodynamic shape can be obtained via an efficient smoothing operation when θ = 0.5 [22]. Therefore, θ is set as 0.5 in this article. In addition, nneighbor denotes the number of FFD control points adjacent to the i-th control point, and it is set as 4 for our studies.
The pseudo code of the Laplacian smoothing method is given in Algorithm 1, which describes the smoothing progress in detail. As shown in Figure 2, the aerodynamic shape of a typical sample obtained randomly via DoE is smoothed by using the Laplacian smoothing method, resulting in a smooth and practical aerodynamic shape.
Algorithm 1 Laplacian smoothing algorithm
Input yinitial /*initial y-axis direction deformations of FFD control point* /
Output ysmoothing /*y-axis direction deformations of FFD control points obtained via Laplacian smoothing* /
Procedure Smoothing (yinitial)
ysmoothing = yinitial
ε g = 0.035 /*convergence criterion is set as 0.035* /
ε g = 0
vsmoothing = FED(yinitial)
/ * perform FFD parameterization and get initial y - axis deirection coordinates of the wing surface mesh * /
while ε < ε g do
   for each FFD control point i do
      y s m o o t h i n g i = θ y s m o o t h i n g i + 1 θ n n e i g h b o r j = 1 n n e i g h b o r y s m o o t h i n g j / *Laplacian smoothing function * / *θ = 0.5 in this article, nneighbor = 4 for the wing* /
   end for
   Vsmoothing = FFD(ysmoothing) / *update the y - axis direction coordinates of the wing surface mesh * /
    ϵ = V s m o o t h i n g V i n i t i a l V i n i t i a l / * calculate the relative error * /
end while
end Procedure

2.3. Flow Solver

We use the open-source code “Adflow” as the flow solver [37]. It is a finite-volume, cell-centered multi-block solver for RANS equations, and can be employed to calculate the gradients required for the optimizations by using the discrete adjoint method.

2.3.1. Governing Equation of Flow and Discretization Scheme

For the control volume V and surface element dS, Navier–Stokes equations in integral form can be defined as
Ω Q t d V + Ω F n d S = Ω F v n d S ,
where Q = ρ , ρ u , ρ v , ρ w , ρ E T denotes the conservation variables. Ω denotes the control volume boundary, n is the surface normal of Ω . q = ( u , v , w ) T is the Cartesian velocity vector. ρ is the density. E is the specific total energy. F , F v denote the inviscid and viscous flux vector, respectively.
The Jameson–Schmidt–Turkel (JST) scheme with second-order accuracy [38] is used for the spatial discretization. The spatial discretization of Equation (2) for each grid cell i can be expressed as
Q i V i t + R ¯ i = 0 ,
where V i is the volume of current grid cell i. R ¯ i = R i c R i v D i is the total residual of Equation (3), R i , R i v denote the residual of the inviscid and the viscous flux respectively, and D i is the artificial dissipation term added to prevent non-physical oscillations.
The first-order backward difference is performed on the time derivative term, Equation (3) can be written into the following implicit form
V i Q i k + 1 Q i k Δ t + R ¯ i k + 1 = 0 ,
where R ¯ i k + 1 can be written as
R ¯ i k + 1 = R i c , k + 1 R i v , k D i k ,
where the superscript k represents the number of the time step. The diagonally dominant alternating direction implicit (DD-ADI) method [39] is employed for the implicit solution. Furthermore, the Spalart–Allmaras (SA) one-equation turbulence model is used to compute turbulence viscosity coefficient.

2.3.2. Discrete Adjoint Method

For ASO problems, derivatives of the cost functions F with respect to the design variables x can be obtained by
d F d x = F x + ψ T R x ,
where ψ denotes the adjoint variables, and R / x represents partial derivatives of the residual of the governing equations with respect to the design variables. ψ can be achieved by solving adjoint equations.
R w T ψ = F w T ,
where w denotes flowfield variables, such as ρ and u , v , w . The preconditioned generalized minimal residual (GMRES) [40] solver from PETSc [41] is used for solving the adjoint equations.

2.4. Mesh Perturbation Method

Since the deformations of the FFD control volumes can change the embedded geometry during optimization, an in-house mesh perturbation code “GridDeform” is used to perturb the CFD mesh according to the new geometry. The core idea of the mesh perturbation method used in this article is to interpolate the deformation of the surface mesh to the far field by using the radial basis function (RBF). In addition, the Intel Math Kernel Library (MKL) is used to solve large linear systems under the parallel computing environment, which makes the in-house mesh perturbation code capable of perturbing the CFD mesh with large-scale cells efficiently in parallel.

2.5. SBO Algorithm

To find the global optimum, a surrogate-based multi-round optimization strategy is introduced into our study. The main loop of this method is shown in Figure 3. Its core idea is to adjust the design space in each round according to the optimal aerodynamic shape obtained from the previous round.
The in-house SBO code “SurroOpt” [8,42,43,44] is used in our studies. Figure 4 sketches the optimization framework of SurroOpt. The optimization progress can be described as
(1)
LHS is selected to generate 100 initial samples. Then, the aerodynamic shapes generated by these samples are smoothed by using the Laplacian smoothing method, and corresponding responses are solved by Adflow;
(2)
Initial kriging models are built based on the samples and corresponding responses;
(3)
The parallel infill-sampling criteria, combining the expected improvement (EI) and minimizing surrogate prediction (MSP), are used to generate two new samples. The new aerodynamic shapes are smoothed by using the Laplacian smoothing method, and corresponding responses are solved by Adflow;
(4)
Kriging models are updated based on these smoothed samples and their responses.
(5)
Steps 3 and 4 are repeated until one of the termination conditions is satisfied.

3. Problem Formulation

3.1. Baseline Geometry

The NASA CRM wing–body–tail configuration is treated as the baseline geometry, as shown in Figure 5. The reference point is located at 25% of mean aerodynamic chord. The nose is located at the origin of coordinates and the span-wise direction is set parallel to the z-axis direction. The reference point is located at (x, y, z) = (31.32972, −0.50934, 0.007669) m. The reference length is 7.00532 m, and the reference area is 191.8448 m2.

3.2. Mesh Convergence Study

Three coarsening levels of multi-block CFD structured grids are generated by ANSYS ICEM CFD and shown in Figure 6. L2, L1, and L0 represent the coarse grid, medium grid, and fine grid, respectively, and each grid has 273 blocks. The first cell height ∆1, mesh size, and y+ values are listed in Table 1. A grid convergence study is performed based on these meshes. The flow condition is the lift coefficient C L = 0.5, the freestream Mach number Ma = 0.85, and the Reynolds number Re = 5 × 106. In addition, the drag coefficient CD of the infinite-cell-grid is evaluated via Richardson extrapolation.
Table 2 lists the CFD simulation results of the baseline meshes. It is obvious that the drag coefficient of the baseline geometry decreases with increasing mesh size, while the moment coefficient increases. The same trend can be seen in Figure 7, which shows the drag and moment convergence plots with respect to 1/MeshSize2/3. The drag coefficient of the L0 grid is only 1.6 counts more than the zero grid spacing drag, which proves that the accuracy of the L0 mesh is close to the infinite-cell-grid.
Figure 8 compares sectional pressure distributions from the CFD simulation results with the experimental data from the National Transonic Facility (NTF) test. We can see that the CFD results get closer to the experimental data with increasing mesh size, and the L1 and L0 grids have sufficient accuracy. After considering the trade-off between computational time and accuracy, L1 mesh is applied to optimizations.

3.3. Optimization Problem Formulation

To ensure the lift constraint CL = 0.5 is satisfied, the angle of attack is treated as a design variable. As shown in Figure 1, the wing FFD volumetric block has a total of 90 control points. To fix the leading and trailing edge of the wing, the two control points at these locations of each control section share one design variable to move the same distance along the opposite directions. We also choose the y-axis direction displacements of the other 70 control points as design variables. Therefore, there are 80 wing shape variables in total.
To ensure that a practical design can be obtained, a volume constraint and a large number of thickness constraints are added to the optimization. The volume is enforced to be no less than that of the baseline wing. As shown in Figure 9, we impose 525 thickness constraints on the wing with 25 spanwise and 21 chordwise locations. These thicknesses are enforced to be no less than 85% of those of the baseline wing. To handle so many geometric constraints efficiently, the constant-parameter Kreisselmeier–Steinhauser (KS) method [45] is adopted to lump all the thickness constraints into one during optimization. Thus, the cost of optimization is greatly reduced. The KS method can be defined as
K S g x = g max x + 1 ρ k s ln j = 1 n e ρ g j x g max x
where gmax(x) denotes the maximum constraint value, and ρ k s is the aggregation parameter.
Two optimization strategies are adopted in this study. The first strategy optimizes the wing–body–tail configuration without considering trimming and then manually trims the optimized configuration by deflecting the horizontal tail as a solid body. The second strategy enforces the balanced pitching moment as a constraint of the optimization, and the trimming is realized automatically by treating the tail deflection angle as a design variable. The tail rotation axis is parallel to the z-axis and is located at 40% of the tail root chord. Figure 10 shows the −6° deflection of the horizontal tail relative to the original configuration. The design variables and constraints used for both strategies are listed in Table 3.

4. Optimization Results

The ASO of the NASA CRM wing–body–tail configuration is performed via SBO and GBO in this section. For SBO, a multi-round optimization strategy is used to drive ASO towards the global optimum. The LHS method is used to select initial samples. EI and MSP are used to select 2 new samples, corresponding responses are solved simultaneously to update the surrogate models. For GBO, the gradient-based optimizer based on sequential quadratic programming (SQP) is used in our studies, and the gradients are provided by Adflow. All optimization works are conducted by the advanced TH-1A computing system at the National Supercomputer Center (NSCC) in Tianjin. Each computing node in this cluster has two 2.6 GHz Intel Xeon CPUs (e5-2690 V4) and 128 GB RAM in total.

4.1. Strategy I: Untrimmed Optimization + Manual Trimming

This strategy optimizes the wing only and then manually trims the optimized configuration by deflecting the horizontal tail as a solid body.
Convergence histories of the objectives are shown in Figure 11. The time costs of the optimizations are listed in Table 4; we can see that SBO is more efficient than GBO for the current example with 81 design variables, since the responses of the two new samples generated by EI and MSP are solved simultaneously in each iteration of SBO. The optimal objective function value obtained via SBO is comparable to that obtained via GBO. Figure 12 shows the pressure distributions of the baseline and optimal configurations, the shock wave on the suction side is notably weakened after the optimization. The drag coefficient drops from 283.77 counts to 277.48 counts after two rounds of optimizations within the framework of SBO, which is only 0.65 counts more than that (276.83 counts) obtained via GBO.
The optimal configurations obtained via Strategy I are then trimmed manually. Figure 13 depicts the pressure distributions of the baseline and trimmed configurations. The tail horizontal deflection angles are −0.835° and −0.899° for SBO and GBO, respectively, which results in zero pitching moment coefficient for the optimal configurations, leaving the total drag coefficient at 281.77 counts and 281.4 counts due to the trim drag. The trim drag coefficient is about 4 counts.

4.2. Strategy II: Trimmed Optimization

For this strategy, the tail deflection angle is treated as a design variable to trim the aircraft, and the balanced pitching moment constraint is enforced to ensure that the pitch moment coefficient is zero.
Convergence histories of the objectives obtained via SBO and GBO are plotted in Figure 14. The time costs of the optimizations are listed in Table 5, we can see that SBO is more efficient than GBO for the design of complex aircraft configurations with 82 design variables, since the responses of the two new samples generated by EI and MSP are solved simultaneously in each iteration of SBO. Figure 15 depicts the pressure distributions of the baseline and trimmed configurations. The shock wave on the suction side is weakened obviously after optimization. The drag coefficient drops from 283.77 counts to 278.61 counts after two rounds of optimizations within the framework of SBO, which is only 0.77 counts different from that (279.38 counts) obtained via GBO.

4.3. Comparison of the Results Obtained by Using Different Optimization Strategies

The optimization results obtained by all optimization strategies are compared, as listed in Table 6. It is obvious that better optimal configurations can be obtained by considering trimming during optimization via Strategy II compared with Strategy I. A total of 3.16-drag-count drag-reducing benefits is achieved by considering trimming during optimization rather than trimming the optimized configuration manually. In addition, the optimization results obtained by GBO and SBO are compared, as listed in Table 7. One can see that the benefit of ASO achieved by SBO is comparable to that obtained by GBO, since the difference between the results obtained by these two optimization algorithms is less than 1 drag count.

5. Conclusions

In this article, the ASO of a full aircraft configuration considering trimming was investigated within the framework of a surrogate-based efficient global optimization. Two optimization strategies were adopted to optimize the NASA CRM wing–body–tail configuration. The first strategy optimizes the wing directly and then manually trims the optimized configuration by deflecting the horizontal tail. The second strategy optimizes the configuration considering trimming by adding a balanced pitching moment constraint during optimization. The optimization results are compared with those obtained via gradient-based optimization with gradients computed by the adjoint method, and the following conclusions can be drawn:
(1) Surrogate-based optimization combining the multi-round optimization strategy and Laplacian smoothing offers great potential in ASO of a full aircraft configuration. For the current test examples, it is more efficient than gradient-based optimization for design problems with around 80 design variables, and the drag reduction is comparable possibly due to the unimodal feature of the NASA CRM wing ASO [46,47]. As the number of design variables increases, gradient-based optimization becomes more efficient than surrogate-based optimization.
(2) The benefit of drag minimization can be expected to lose approximately 4 drag counts when trimming the optimal configuration obtained via untrimmed optimization. Therefore, when optimizing the wing–body–tail configuration, it is necessary to consider the influence of the trim drag on the optimization results.
(3) Compared with the strategy that manually trims the optimal configuration obtained from an untrimmed optimization, a trimmed optimization can achieve 3.16-drag-count drag-reducing benefits by enforcing the balanced pitching moment constraint during optimization.
Beyond the scope of this article, gradient-enhanced surrogate model will be investigated for the ASO of a full aircraft configuration.

Author Contributions

Conceptualization, K.W. and Z.H.; methodology, K.W. and Z.H.; software, K.W. and Z.H.; validation, Z.H., K.Z. and W.S.; formal analysis, K.W.; investigation, K.W.; resources, K.W., Z.H. and W.S.; data curation, K.W.; writing—original draft preparation, K.W.; writing—review and editing, Z.H., K.Z. and W.S.; visualization, K.W.; supervision, Z.H. and K.Z.; project administration, K.W.; funding acquisition, Z.H. and K.Z. All authors have read and agreed to the published version of the manuscript.

Funding

This research is funded by the National Natural Science Foundation of China (grant No. U20B2007 and 11972305), as well as the Shaanxi Science Fund for Distinguished Young Scholars (grant No. 2020JC-13) and Natural Science Fund of Shaanxi Province (grant No. 2020JM-127).

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Queipo, N.V.; Haftka, R.T.; Shyy, W.; Goela, T.; Vaidyanathan, R.; Tucker, P.K. Surrogate-based analysis and optimization. Prog. Aerosp. Sci. 2005, 45, 1–28. [Google Scholar] [CrossRef]
  2. Forrester, A.I.J.; Keane, A.J. Recent advances in surrogate-based optimization. Prog. Aerosp. Sci. 2009, 45, 50–79. [Google Scholar] [CrossRef]
  3. Haftka, R.T.; Villanueva, D.; Chaudhuri, A. Parallel surrogate-assisted global optimization with expensive functions-a survey. Struct. Multidiscip. Optim. 2016, 54, 3–13. [Google Scholar] [CrossRef]
  4. Parr, J.M.; Keane, A.J.; Forrester, A.I.J.; Holden, C.M.E. Infill sampling criteria for surrogate-based optimization with constraint handling. Eng. Optimiz. 2012, 44, 1147–1166. [Google Scholar] [CrossRef]
  5. Han, Z.H. Kriging surrogate model and its application to design optimization: A review of recent progress. Acta Aeronaut. Astronaut. Sin. 2016, 37, 3197–3225. (In Chinese) [Google Scholar]
  6. Viana, F.A.C.; Simpson, T.W.; Balabanov, V.; Toropov, V. Metamodeling in multidisciplinary design optimization: How far have we really come? AIAA J. 2014, 52, 670–690. [Google Scholar] [CrossRef]
  7. Han, Z.H.; Goertz, S.; Zimmermann, R. Improving variable-fidelity surrogate modeling via gradient enhanced Kriging and a generalized hybrid bridge function. Aerosp. Sci. Tech. 2013, 25, 177–189. [Google Scholar] [CrossRef]
  8. Han, Z.H. SurroOpt: A generic surrogate-based optimization code for aerodynamic and multidisciplinary design. In Proceedings of the ICAS 2016, Daejeon, Republic of Korea, 25–30 September 2016. [Google Scholar]
  9. Liu, J.; Song, W.P.; Han, Z.H.; Zhang, Y. Efficient aerodynamic shape optimization of transonic wings using a parallel infilling strategy and surrogate models. Struct. Multidiscip. Optim. 2017, 55, 925–943. [Google Scholar] [CrossRef]
  10. Han, Z.H.; Xu, C.Z.; Zhang, L.; Zhang, Y.; Zhang, K.S.; Song, W.P. Efficient aerodynamic shape optimization using variable-fidelity surrogate models and multilevel computational grids. Chin. J Aeronaut. 2020, 33, 31–47. [Google Scholar] [CrossRef]
  11. Courrier, N.; Boucard, P.A.; Soulier, B. Variable-fidelity modeling of structural analysis of assemblies. J. Global. Optim. 2016, 64, 577–613. [Google Scholar] [CrossRef]
  12. Ha, H.; Oh, S.; Yee, K. Feasibility study of hierarchical kriging model in the design optimization process. J. Korean Soc. Aero. Spac. Sci. 2014, 42, 108–118. [Google Scholar] [CrossRef]
  13. Jo, Y.; Choi, S.; Lee, D.J.; Choi, D.Z. Adaptive variable-fidelity analysis and design using dynamic fidelity indicators. AIAA J. 2016, 54, 3564–3579. [Google Scholar] [CrossRef]
  14. Leifsson, L.; Koziel, S. Aerodynamic shape optimization by variable-fidelity computational fluid dynamics models: A review of recent progress. J. Comput. Sci. 2015, 10, 45–54. [Google Scholar] [CrossRef]
  15. Leifsson, L.; Koziel, S.; Tesfahunegn, Y.A. Multiobjective aerodynamic optimization by variable-fidelity models and response surface surrogate models. AIAA J. 2016, 54, 531–541. [Google Scholar] [CrossRef]
  16. Zhang, Y.; Han, Z.H.; Zhang, K.S. Variable-fidelity expected improvement for efficient global optimization of expensive functions. Struct. Multidiscip. Optim. 2018, 58, 1431–1451. [Google Scholar] [CrossRef]
  17. Han, Z.H.; Zhang, Y.; Song, C.X.; Zhang, K.S. Weighted gradient-enhanced kriging for high-dimensional surrogate modeling and design optimization. AIAA J. 2017, 55, 4330–4346. [Google Scholar] [CrossRef]
  18. Laurenceau, J.; Meaux, M.; Montagnac, M.; Sagaut, P. Comparison of gradient-based and gradient-enhanced response-surface-based optimizers. AIAA J. 2012, 48, 981–994. [Google Scholar] [CrossRef]
  19. Mortished, C.; Ollar, J.; Jones, R.; Benzie, P.; Toropov, V.; Sienz, J. Aircraft wing optimization based on computationally efficient gradient-enhanced Kriging. In Proceedings of the 57th AIAA/ASCE/AHS/ASC Structures, Structural Dynamics, and Materials Conference, San Diego, CA, USA, 4–8 January 2016. [Google Scholar]
  20. Bouhlel, M.A.; Martins, J.R.R.A. Gradient-enhanced kriging for high-dimensional problems. Eng. Comput. Ger. 2019, 35, 157–173. [Google Scholar] [CrossRef]
  21. Li, J.C.; Zhang, M. Adjoint-free Aerodynamic Shape Optimization of the Common Research Model Wing. AIAA J. 2021, 59, 1990–2000. [Google Scholar] [CrossRef]
  22. Li, J.C.; Cai, J.S.; Qu, K. Surrogate-based aerodynamic shape optimization with the active subspace method. Struct. Multidiscip. Optim. 2019, 59, 403–419. [Google Scholar] [CrossRef]
  23. Bouhlel, M.A.; Bartoli, N.; Regis, R.G.; Otsmane, A.; Morlier, J. Efficient global optimization for high-dimensional constrained problems by using the Kriging models combined with the partial least squares method. Eng. Optimiz. 2018, 50, 2038–2053. [Google Scholar] [CrossRef]
  24. Zhang, K.S.; Han, Z.H.; Gao, Z.J.; Wang, Y. Constraint aggregation for large number of constraints in wing surrogate-based optimization. Struct. Multidiscip. Optim. 2019, 59, 421–438. [Google Scholar] [CrossRef]
  25. Wang, K.; Han, Z.H.; Zhang, K.S.; Song, W.P. An Efficient Geometric Constraint Handling Method for Surrogate-Based Aerodynamic Shape Optimization. Eng. Appl. Comp. Fluid 2023, 17, e2153173. [Google Scholar] [CrossRef]
  26. Zhang, Y.; Han, Z.H.; Shi, L.X.; Song, W.P. Multi-round Surrogate-based Optimization for Benchmark Aerodynamic Design Problems. In Proceedings of the 54th AIAA Aerospace Sciences Meeting, San Diego, CA, USA, 4–8 January 2016. [Google Scholar]
  27. Li, J.C.; Zhang, M. Drag Reduction of Transonic Wings with Surrogate-Based Optimization. In Proceedings of the 2018 Asia-Pacific International Symposium on Aerospace Technology, Chengdu, China, 15–19 October 2018. [Google Scholar]
  28. Han, Z.H.; Zhang, Y.; Xu, C.Z.; Wang, K.; Wu, M.M.; Zhu, Z.; Song, W.P. Aerodynamic optimization design of large civil aircraft wings using surrogate-based model. Acta Aeronaut. Astronaut. Sin. 2019, 40, 522398. (In Chinese) [Google Scholar]
  29. Chen, S.; Lyu, Z.J.; Kenway, G.K.W.; Martins, J.R.R.A. Aerodynamic shape optimization of common research model wing-body-tail configuration. J. Aircr. 2016, 53, 276–293. [Google Scholar] [CrossRef]
  30. Lyu, Z.J.; Martins, J.R.R.A. Aerodynamic design optimization studies of a blended-wing-body aircraft. J. Aircr. 2014, 51, 1604–1617. [Google Scholar] [CrossRef]
  31. Koo, D.; Zingg, D.W. Progress in Aerodynamic Shape Optimization Based on the Reynolds-Averaged Navier-Stokes Equations. In Proceedings of the 54th AIAA Aerospace Sciences Meeting, San Diego, CA, USA, 4–8 January 2016. [Google Scholar]
  32. Merle, A.; Stück, D.W.; Rempke, A. An Adjoint-based Aerodynamic Shape Optimization Strategy for Trimmed Aircraft with Active Engines. In Proceedings of the 35th AIAA Applied Aerodynamics Conference, Denver, CO, USA, 5–9 June 2017. [Google Scholar]
  33. Kenway, G.K.W.; Secco, N.; Martins, J.R.R.A.; Mishra, A.; Duraisamy, K. An Efficient Parallel Overset Method for Aerodynamic Shape Optimization. In Proceedings of the 58th AIAA/ASCE/AHS/ASC Structures, Structural Dynamics, and Materials Conference, Grapevine, TX, USA, 9–13 January 2017. [Google Scholar]
  34. Lei, R.W.; Bai, J.Q.; Xu, D.Y. Aerodynamic optimization of civil aircraft with wing-mounted engine jet based on adjoint method. Aerosp. Sci. Tech. 2019, 93, 105285. [Google Scholar] [CrossRef]
  35. Sederberg, T.W.; Parry, S.R. Free-form deformation of solid geometric models. ACM Siggraph Comput. Graph. 1986, 20, 151–160. [Google Scholar] [CrossRef]
  36. Giunta, A.A.; Wojtkiewicz, J.S.F.; Eldred, M.S. Overview of Modern Design of Experiments Methods for Computa-tional Simulations. In Proceedings of the 41st AIAA Aerospace Sciences Meeting and Exhibit, Reno, NV, USA, 6–9 January 2003. [Google Scholar]
  37. Mader, C.A.; Kenway, G.K.W.; Yildirim, A.; Martins, J.R.R.A. ADflow: An Open-Source Computational Fluid Dynamics Solver for Aerodynamic and Multidisciplinary Optimization. J. Aerosp. Inf. Syst. 2020, 17, 508–527. [Google Scholar] [CrossRef]
  38. Jameson, A.; Schmidt, W.; Turkel, E. Numerical Solution of the Euler Equations by Finite Volume Methods Using Runge–Kutta Time Stepping Schemes. In Proceedings of the 14th Fluid and Plasma Dynamics Conference, Palo Alto, CA, USA, 23–25 June 1981. [Google Scholar]
  39. Klopfer, G.; Hung, C.; Wijngaart, R.V.D.; Onufer, J. A diagonalized diagonal dominant alternating direction implicit (D3ADI) scheme and subiteration correction. In Proceedings of the 29th AIAA, Fluid Dynamics Conference, Albuquerque, NM, USA, 15–18 June 1998. [Google Scholar]
  40. Saad, Y.; Schultz, M.H. GMRES: A Generalized Minimal Residual Algorithm for Solving Nonsymmetric Linear Systems. Siam J. Sci. Comput. 1986, 7, 856–869. [Google Scholar] [CrossRef]
  41. Balay, S.; Brown, J.; Buschelman, K.; Eijkhout, V.; Gropp, W.D.; Kaushik, D.; Knepley, M.G.; McInnes, L.C.; Smith, B.F.; Zhang, H. PETSc Users Manual, Technical Report ANL-95/11—Revision 3.4; Argonne National Laboratory: Lemont, IL, USA, 2013. [Google Scholar]
  42. Han, Z.H.; Chen, J.; Zhang, K.S.; Xu, Z.M.; Zhu, Z.; Song, W.P. Aerodynamic Shape Optimization of Natural-Laminar-Flow Wing Using Surrogate-Based Approach. AIAA J. 2018, 56, 2579–2593. [Google Scholar] [CrossRef]
  43. Bu, Y.P.; Song, W.P.; Han, Z.H.; Zhang, Y.; Zhang, L. Aerodynamic/aeroacoustic variable-fidelity optimization of helicopter rotor based on hierarchical Kriging model. Chin. J. Aeronaut. 2020, 33, 476–492. [Google Scholar] [CrossRef]
  44. Han, Z.H.; Xu, C.Z.; Qiao, J.L.; Liu, F.; Chi, J.B.; Meng, G.Y.; Zhang, K.S.; Song, W.P. Recent progress of efficient global aerodynamic shape optimization using surrogate-based approach. Acta Aeronaut. Astronaut. Sin. 2020, 41, 623344. (In Chinese) [Google Scholar]
  45. Kreisselmeier, G.; Steinhauser, R. Systematic control design by optimizing a vector performance index. In Proceedings of the IFAC Symposium on computer Aided Design of Control Systems, Zurich, Switzerland, 29–31 August 1979. [Google Scholar]
  46. He, X.L. Multi-Component Adaptive Aerodynamic Shape Optimization Methods for Civil Transport Aircrafts. Ph.D. Thesis, Northwestern Polytechnical University, Wuhan, China, 2019. (In Chinese). [Google Scholar]
  47. Yu, L.; Lyu, Z.J.; Xu, Z.L.; Martins, J.R.R.A. On the influence of optimization algorithm and initial design on wing aerodynamic shape optimization. Aerosp. Sci. Tech. 2018, 75, 183–199. [Google Scholar] [CrossRef]
Figure 1. FFD parameterization for the NASA CRM wing–body–tail configuration.
Figure 1. FFD parameterization for the NASA CRM wing–body–tail configuration.
Aerospace 10 00734 g001
Figure 2. Comparison of wing–body–tail configurations before and after Laplacian smoothing.
Figure 2. Comparison of wing–body–tail configurations before and after Laplacian smoothing.
Aerospace 10 00734 g002
Figure 3. Flowchart of multi-round optimization strategy.
Figure 3. Flowchart of multi-round optimization strategy.
Aerospace 10 00734 g003
Figure 4. Framework of SurroOpt [8,42,43,44].
Figure 4. Framework of SurroOpt [8,42,43,44].
Aerospace 10 00734 g004
Figure 5. NASA CRM wing–body–tail configuration.
Figure 5. NASA CRM wing–body–tail configuration.
Aerospace 10 00734 g005
Figure 6. O-H multiblock structured meshes for the mesh convergence study.
Figure 6. O-H multiblock structured meshes for the mesh convergence study.
Aerospace 10 00734 g006
Figure 7. Mesh convergence for drag and moment coefficients of the baseline geometry.
Figure 7. Mesh convergence for drag and moment coefficients of the baseline geometry.
Aerospace 10 00734 g007
Figure 8. Comparison of sectional pressure distributions of the baseline geometry.
Figure 8. Comparison of sectional pressure distributions of the baseline geometry.
Aerospace 10 00734 g008
Figure 9. The 525 (21 chordwise × 25 spanwise) geometric constraints enforced on the wing.
Figure 9. The 525 (21 chordwise × 25 spanwise) geometric constraints enforced on the wing.
Aerospace 10 00734 g009
Figure 10. Tail deflection resulting from the deflection of a sub-FFD block.
Figure 10. Tail deflection resulting from the deflection of a sub-FFD block.
Aerospace 10 00734 g010
Figure 11. Convergence history of the total drag coefficient (Strategy I).
Figure 11. Convergence history of the total drag coefficient (Strategy I).
Aerospace 10 00734 g011
Figure 12. Comparison of pressure distributions of the baseline with untrimmed optimal configurations from Strategy I.
Figure 12. Comparison of pressure distributions of the baseline with untrimmed optimal configurations from Strategy I.
Aerospace 10 00734 g012
Figure 13. Comparison of pressure distributions of the baseline with optimal configurations trimmed manually from Strategy I.
Figure 13. Comparison of pressure distributions of the baseline with optimal configurations trimmed manually from Strategy I.
Aerospace 10 00734 g013
Figure 14. Convergence history of the total drag coefficient (Strategy II).
Figure 14. Convergence history of the total drag coefficient (Strategy II).
Aerospace 10 00734 g014
Figure 15. Comparison of pressure distributions of the baseline and optimal configurations from Strategy II.
Figure 15. Comparison of pressure distributions of the baseline and optimal configurations from Strategy II.
Aerospace 10 00734 g015
Table 1. Details of baseline CFD grids.
Table 1. Details of baseline CFD grids.
Number1Y+Number of Cells (Million)
L23.69 × 10−51.00.887
L12.46 × 10−52/36.776
L01.64 × 10−54/959.355
Table 2. CFD simulation results of the mesh convergence study.
Table 2. CFD simulation results of the mesh convergence study.
Mesh LevelsCLCDCMZ
L000.50.02720/
L00.50.02736−0.0422
L10.50.02838−0.0463
L20.50.03351−0.0640
Table 3. Summary of all optimization strategies.
Table 3. Summary of all optimization strategies.
Optimization StrategyDesign VariablesObjectiveConstraints
DescriptionQuantityFunctionDescriptionQuantity
Strategy IAngle of attack1Min. CDCL = 0.5Lift constraint1
t−0.85t0 ≥ 0Geometric constraints525
vv0Volume constraint1
Wing shape80 Δ y L E upper = Δ y LE lower Fix the leading edge5
Δ y TE upper = Δ y TE lower Fix the trailing edge5
Strategy IIAngle of attack1Min. CDCL =0.5 1
CMZ= 0zero pitching moment constraint1
Wing shape80t−0.85t0 ≥ 0 525
vv0 1
Tail deflection angle1 Δ y LE upper = Δ y LE lower 5
Δ y TE upper = Δ y TE lower 5
Table 4. Time cost of optimization via Strategy I.
Table 4. Time cost of optimization via Strategy I.
Optimization AlgorithmComputational Cost (Hours)
GBO182.27
SBO165.58
Table 5. Time cost of optimization via Strategy II.
Table 5. Time cost of optimization via Strategy II.
Optimization AlgorithmComputational Cost (Hours)
GBO218.14
SBO200.25
Table 6. Comparison of the results of all optimization strategies.
Table 6. Comparison of the results of all optimization strategies.
Optimization StrategyCLCD
(Counts)
CMZTail Deflection AngleAngle of Attack
Strategy Iuntrimmed0.5277.48−0.056492.2432°
trimmed manually0.5281.770−0.8350°2.3487°
Strategy II0.5278.610−0.3996°2.3604°
Table 7. Comparison of the results obtained via different optimization methods.
Table 7. Comparison of the results obtained via different optimization methods.
Optimization MethodCLCD
(Counts)
CMZTail Deflection AngleAngle of Attack
GBOStrategy I: trimmed manually0.5281.400−0.8990°2.4087°
Strategy II0.5279.380−0.4065°2.3543°
SBOStrategy I: trimmed manually0.5281.770−0.8350°2.3487°
Strategy II0.5278.610−0.3996°2.3604°
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Wang, K.; Han, Z.; Zhang, K.; Song, W. Efficient Global Aerodynamic Shape Optimization of a Full Aircraft Configuration Considering Trimming. Aerospace 2023, 10, 734. https://doi.org/10.3390/aerospace10080734

AMA Style

Wang K, Han Z, Zhang K, Song W. Efficient Global Aerodynamic Shape Optimization of a Full Aircraft Configuration Considering Trimming. Aerospace. 2023; 10(8):734. https://doi.org/10.3390/aerospace10080734

Chicago/Turabian Style

Wang, Kai, Zhonghua Han, Keshi Zhang, and Wenping Song. 2023. "Efficient Global Aerodynamic Shape Optimization of a Full Aircraft Configuration Considering Trimming" Aerospace 10, no. 8: 734. https://doi.org/10.3390/aerospace10080734

APA Style

Wang, K., Han, Z., Zhang, K., & Song, W. (2023). Efficient Global Aerodynamic Shape Optimization of a Full Aircraft Configuration Considering Trimming. Aerospace, 10(8), 734. https://doi.org/10.3390/aerospace10080734

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop