1. Introduction
The formulation of design requirements is assumed critical for effective optimization in all design disciplines [
1], since they characterize the design process through design space exploration for the best solution via evaluation and analysis [
2,
3,
4,
5]. This paper is focused on providing new perspective from objective space via formulating problem-oriented design requirement in the field of aerodynamic design such as Natural Laminar Flow optimization of a transonic airfoil and aero engine compressor blade.
The destructive impact of carbon emissions and the enduring bad economic situation of airlines makes it necessary for aircraft manufacturers to minimize fuel burn in aircraft design, one solution to which has been drag reduction [
6,
7]. The skin friction is important part of drag composition (it occupies nearly 50% of total aircraft drag in some cases [
6]). Therefore, Natural Laminar Flow (NLF) optimization has been worked on for skin friction drag reduction by extending laminar area over aircraft surface, on account of lower skin friction drag magnitude in laminar than in turbulence [
6,
8,
9].
In many aerodynamic designs, numerical optimization algorithms based on geometrical parameterization are used for the samples in design space to be iterated towards a pre-defined design requirement. The classification of skin friction drag and pressure drag comes from their different generation mechanisms and the fact that they can be easily calculated without involving the other’s information. The formulation of design requirement in NLF optimization usually reduces to the summation of pressure drag
and laminar-related performance (i.e., skin friction drag
) as in Equation (
1)
where
and
are the assigned weights that represents the importance attached to the aerodynamic performance behind variables given by designers. Wang et al. and Iuliano et al. put
under major consideration (solely minimize
) in design requirement formulation with
being kept under a passive constant constraint [
8,
10]. Chen and Tang and Tang et al. added wave drag (a subset of pressure drag) in design requirement formulation where they put much consideration on shock wave at transonic regime [
11,
12]. Han et al. and Rashad and Zingg simply put total drag as the design requirement [
13,
14,
15]. Wang et al. set 80% to
in Equation (
1) for increasing the effect of laminarization
[
16].
The above settings represent different perspectives in design process. From industrial view, making total drag of a geometry alone as the design requirement is ultimate, but this would hinder the effort of methodological manipulation in emphasizing NLF’s influence [
17,
18,
19,
20]. Usually, it is not how much importance should be adequately given to laminarization that matters in such cases, but whether laminarization is considered in drag reduction design. In fact, these studies are more of verification of transition models in the prediction of laminar transition with coupling of ordinary drag force computation. Adjusting the corresponding weights in Equation (
1) is a solution to the challenge of reasonable design requirement formulation. The problem, however, is that this is usually done according to experience. For example, simply singling out laminarization in the design requirement formulation ignores the fact that the change of geometry is also linked to the change of pressure drag force as well. Simple collocation of skin friction drag and pressure drag in multi-objective problem of min
is also not a perfect solution as well in that it ignores the factor of impact of each objective. However, it can be used as a good measurement of optimization progress, and it will be used in the first case in this paper.
This paper proposes to formulate design requirement of NLF optimization from the perspective of objective space via data manipulation approach. For an aerodynamic design with given algorithm and geometric representation method (i.e., a given design space), the objective space is the projection from design space via aerodynamics or CFD (Computational Fluid Dynamics) models. This paper believes that objective space contains the intrinsic information of the ‘optimization potentials’ of samples in terms of aerodynamic performance: thus, it can be injected in the design requirement formulation as an aide to capture the nature of the problem at hand. In the two proof-of-concept studies in this paper, the improved formulation of concerned design requirement injected with objective space information is compared with traditional design requirement formulation with evenly distributed weight distribution in terms of design efficiency and comprehensive aerodynamic performance, respectively.
The solution proposed in this paper is done via Principal Component Analysis (PCA). PCA is used to strengthen preprocessing of data to improve manipulation efficiency [
21,
22], dimension reduction [
23], and in the sample selection of initialization [
24], etc. The methods for data manipulation have been utilized in aerodynamic designs, but mostly in design space than in objective space. In fact, apart from extending the subjects of aerodynamic design to complex aircraft components [
25,
26,
27] and strengthening the link between geometry and aerodynamic performance [
21,
28,
29,
30], aerodynamic designers deployed data preparation in design space for preprocessing, e.g., dimensional reduction by parameter filtering according to physical intuitiveness(e.g., cutting out geometry parameters with less impact on ultimate aerodynamic design goal [
31]) and feature extracting from flowfield for economical use of data [
32,
33], etc. It is noted that the intention of PCA in this paper is not simply for variable reduction but for extracting the main modes from objective space considered to be the essential information that may help establish a more reasonable formulation of design requirement in aerodynamic designs.
The rest of the paper is organized as follows:
Section 2 studies the CFD solver used in the paper, which is validated by wind-tunnel experiment. First, grid sampling is applied to prepare for the creation of objective space, which offers data to be processed by PCA so that a more problem-oriented formulation of design requirement can be obtained. The eigenvectors of the result are projected onto the principal axis generated via PCA and thus forms the formulation of design requirement. Secondly, the approach of objective space analysis is applied to an NLF optimization of transonic airfoil under fixed lift via Differential Evolution algorithm, the result of which shows that the information of objective space help improve the optimization route in aerodynamic design (cf.
Section 3). Thirdly, a reformulated design requirement is also applied to an NLF aero engine compressor blade design, showing that a third aerodynamic performance objective’s concern can be accessorily reflected in a two-objective problem if an objective space of three aerodynamic performances (i.e., the pressure drag and laminar length related drag as well as pressure loss in the case of blade design) is created and influence design requirement formulation (cf.
Section 4).
3. Proof-of-Concept Study 1: A Differential Evolution Optimization with Conventional and Improved Design Requirement Formulation
3.1. NLF Optimization of RAE2822 Airfoil in Transonic Flow under Fixed Lift
The AIAA Aerodynamic Design Optimization Discussion Group (ADODG) has defined benchmark optimization problem of airfoil RAE2822 in viscous (Reynolds number is 6.5 million), transonic flow (freestream Mach number is 0.734), subject to lift (0.824), pitching moment (not smaller than −0.092), and area (decrement is not allowed) constraints [
36]. The problem has been studied by many scholars, which has significance with respect to engineering applications [
37,
38,
39,
40].
This paper intends to introduce NLF concerns (i.e., skin friction drag) to the original problem. For better illustration of the effect of laminar design, the fixed value of lift force is therefore modified form 0.824 to 0.5 to prevent the angle of attack being too high which goes against the favorable direction of laminar design. In recognition of the above, the authors formulate an optimization problem as follows:
where
,
,
are the lift, pressure drag and skin friction drag applied on the airfoil samples;
is the angle of attack;
y is the vertical coordinates of airfoil;
A is the area occupied by the airfoil samples while
is the area occupied by RAE2822. The farfield flow comes at Mach number 0.734. The lift coefficient is fixed (hence the angle of attack is determined by dichotomy method in this paper).
and are the weights assigned to problem. This paper will set up two optimizations with different design requirement. In the basic case, the weights will be assigned as and for and ; in the case with modified design requirement formulation, the weights will be calculated via analysis in the following section with consideration of objective space information.
3.2. Design Variables in Geometry Description
The parametric geometric representation approach Class-Shape Transformation (CST) conquered many aerodynamic designs with its smoothness and convenience in shape description [
31,
41,
42,
43]. Conventional CST (Equation (
5)) imposes shape functions (Equation (
3)) on specific class function (Equation (
4)) with accumulated Bernstein polynomials assuring shape smoothness [
44,
45,
46], where
n is the scale of the method;
u is the non-dimensional position,
is the trailing edge thickness correction (equal to zero in a sharp edge case);
are the parameters. For ordinary airfoil surfaces, two parameters:
and
take the value of 1 and 0.5, respectively.
The difference in form of
-norm of vertical coordinates between
Class_Shape(
u) function and original data can be minimized until the parameters of
are obtained via least square method:
A CST method of 11 scale is applied to get airfoil RAE2822’s shape variation, the Bernstein polynomial decomposition of which is shown in
Figure 4. The x-y plot represents coordinates of different Bernstein polynomials. For given x position, the summation of local values of Bernstein polynomials is the actual y coordinate of airfoil. The difference between the original RAE2822 airfoil surface and fitted one via CST is shown in
Figure 5, which meets the demand of accuracy of model tolerance of
m recommended by [
35].
The design space is created based on the variation of the first 4 CST parameters representing top, front region of the airfoil RAE2822. Each valid airfoil sample needs to meet detailed geometric constraints: not only the parameter is variated within bounding limits, the generated shapes have to be verified later in terms of curvature distribution, area check, and monotonicity check (cf.
Figure 6):
The curvature () distribution has to be within ±0.1 based on that of the original airfoil surface.
The monotonicity must be consistent with that of the original airfoil surface to avoid twisted shape which is definitely harmful to the approaching towards the target via Differential Evolution algorithm.
And constraint of maximum area of 125% of the initial area is enforced to avoid shape with too high drag before CFD calculation.
These constraints affect the data domain and its distribution in design space as well as objective space, which has impact on the final formulation of design requirement.
Principal Component Analysis as Data Manipulator
PCA is introduced to reduce the number of variables or design space dimensions on the premise of maintaining/capturing intrinsic information (e.g., the dominant physics of a system) in the original empirical data, when the problem requires expensive high-fidelity simulations [
47,
48,
49]. PCA is also referred to as Karhunen–Loève transformation, Proper Orthogonal Decomposition, singular systematic analysis, singular value decomposition under different circumstances, which has other variations in methodologies, e.g., kernel PCA, probabilistic nonlinear PCA/Gaussian process latent variable models, etc. [
50]. The method is proposed against traditional variable screening that reduce dimensionality by simply removing variables; the latter often fails because there may be correlations between each variable in high-dimensional data [
51]. PCA can solve that problem by providing new linear combinations of the original variables representing orthogonal eigenvectors of principal component axes in space, each of which with different variances [
23,
24]. The first principal component is associated with the largest variance, and the following ones are perpendicular to the previous ones [
52].
In some aerodynamic designs, the application of PCA is mainly on the data processing, e.g., classifier. Some examples can be seen in reduced-order models for simulations balancing computation cost and flowfield accuracy in numerous fields [
22,
53,
54,
55,
56].
The data-oriented technique has been studied to accelerate the process of capturing the dominant modes of a system [
57]. This type of PCA applications is focused on the idea that if samples of similar features are clustered, then the methods can be applied to them more effectively.
Beyond that, PCA is not often applied to extract useful information linked with objective space. One similar example is found in [
58] where PCA is used to guide the application of evolution operators in Genetic Algorithm where the front of non-dominated solutions is divided into sub-fronts. Its information was obtained from the calculated aerodynamic performance distribution; then, the result was fed back to the design space to further manipulate the samples. This guidance from objective space is no longer simply an aerodynamic performance progress of individual samples as in many traditional genetic algorithm; on the contrary, the extracted global information of objective space alleviates the dependence of randomness in the design space, which brings benefits to their study.
In such spirit, this paper is also focused on studying objective space information by PCA via the formulation of design requirement, which constitutes the novelty of this paper.
3.3. Information Extracted from Objective Space
If there is some information that represents some features to be deployed in optimization process, then it should produce result that is irrelevant of sampling strategy. With the approach of enumeration, a grid strategy is used for sampling of different scales. An index
n indicates the evenly distribution of sampling seeds in every CST parameter, or sampling scale, as is shown in
Table 3. The total number of samples easily exceeds reasonable calculation limitation if the value of
n is large enough because of iteration through all the parameters. This traditional sampling strategy is for revealing the intrinsic nature hidden in objective space. The huge amount of calculation corresponding to different values of
n is not necessary in the real practice: if the relevant data is enough, then the eigenvectors can be estimated within approximation.
For airfoil samples corresponding to given
n value, their aerodynamic performance of
n airfoils obtained through CFD constitute a database represented by array (
7), where the two columns are skin friction drag and pressure drag force, respectively.
The data from objective space systematically forms the covariance matrix
in Equation (
8). The spectral decomposition is shown in Equation (9) where
represents diagonal matrix of eigenvalues and
represents agglomeration of eigenvectors.
The objective space and eigenvectors are obtained for each
n in
Table 3, the result of which is shown in
Figure 7 where green, solid lines represent first principal component and the red, dashed lines represent second principal component. Since the design requirement in this case involves two aerodynamic performances, only two components are required. The boundary of the objective space can be clearly seen, which is caused by the various geometrical constraint posed by the problem. Similar to the ‘convergence’ of mesh independence check, the directions of eigenvectors are also ‘converged’ to a stable direction with the increasing of
n (i.e., the directions of eigenvectors for objective space for
n equals 4 and 5 are close), which shows the intrinsic information in the objective space determined by the problem description can be captured. The paper thinks that the understanding of the objective space and eigenvectors is directly associated with a more reasonable and problem-oriented formulation of design requirement. With the sampling, a comprehensive pattern that the data is arranged in objective space can be perceived, where the eigenvectors can be a quantitative description of it.
Therefore, if the formulation of design requirement adopts such information, the optimization can be better guided by leading every sample to move in the direction that is parallel to the best direction decided by the nature of the problem. Also, the eigenvectors has engineering meaning in that a better optimization criterion can be proposed in industrial aerodynamic design process as a balanced solution to the ultimate effect of total drag reduction and the perspective of emphasizing part of drag reduction via methodological guidance.
In the next section, a newly formulated design requirement is proposed to compete with conventional design requirement in the Differential Evolution NLF optimization to the transonic airfoil RAE2822. A more reasonable, problem-oriented weight is determined by the PCA analysis in this section, where is 0.878 and is 0.122, while the conventional design requirement is 0.5 and 0.5.
3.4. Evolutionary Optimization and Its Significance
The problem description and data manipulation in the previous sections are the preparation for the Differential Evolution optimization in this section. As a variant of an evolutionary algorithm, Differential Evolution (DE) iteratively improves a candidate solution based on floating-point value representation rather than the binary strings [
59]. It has been deployed in many aerodynamic designs [
8,
60,
61].
The flow chart in
Figure 8 is a possible application of extracting information from objective space in current evolutionary types of optimization. A better formulation of design requirement encourages initialized samples to move in the directions that are intrinsic to the nature of the problem.
3.5. Realization of Algorithm
Since the purpose of this paper is to study the impact of objective formulation on the result of aerodynamic design, the realization of optimization algorithm in different cases must be kept in similar configuration. The algorithm starts with a randomly populated initial generation of D-dimensional parameter vectors representing candidates with considerations of aforementioned constraints in the design and objective space, where is the number of parents in a population and is the number of parameters in this paper. After initializing the population, each target vector of candidate in generation undergoes mutation and crossover operation before producing a sample to be tested. The perturbations on the parameters are the factors to be multiplied by the original values (instead of the values to be treated as profile description parameters) thus assuring the uniformity in random seeding. Only when the newly generated candidate is within the various constraint and has better aerodynamic performance measured by conventional/improved optimization criterion (i.e., design requirement) than its parent, then the old candidate is updated with the new one. In this paper, surrogate models are not used to avoid inaccuracy in description of design and objective space.
The two algorithms with different design requirement are realized in two cases. As is shown in Algorithm 1,
x represents candidate,
v represents guided possible marching step,
represent random candidate in current generation; candidates, constant
F represents scaling factor in amplifying differential steps.
Algorithm 1: Differential Evolution in this paper. |
|
The randomness in the evolutionary marching of both cases is minimized to illustrate the comparison of the effect of design requirement, which is the core of this paper. Therefore, the mutation ratio and crossover operation is kept the same for the same serial number of candidates in the same generation in both cases; the selection of samples for obtaining differential vector is the same in terms of the serial number of samples (e.g., Nos. 3, 7, 11 are picked for generating differential vector in both cases, although they represent totally different locations in the objective space.).
3.6. Result Analysis
The optimization process took about 120 h with each CFD simulation occupying approximately 45 min, considering the fact that the evolution has to be processed in non-parallel way for each candidate of each generation. When the iteration with improved design requirement comes to generation no. 79, there has been no updates of candidates in objective space for the last 10 generations. Therefore, a ‘convergence’ of Differential Evolution is considered to be obtained. The result in
Figure 9 shows that the improvement in design requirement formulation takes effect in enhancing the optimization by applying the intrinsic nature of the problem via objective space analysis. With the improved design requirement formulation, the Pareto front can be obtained faster than conventional configuration. The cost for extracting information from objective space brings rewards in the final optimization process.
The Pareto front represented as thick, purple dashed line shows that due to the various constraint in the problem description, most of the candidates are in the region relatively lower than the corner of the line. Therefore, most of the optimization is done mainly in the aspect of pressure drag reduction(e.g., candidate no. 2), while some in the aspect of laminar length increasing (e.g., candidate no. 12). The candidate no. 5 is selected as the optimized airfoil to be compared with original airfoil RAE2822 because it is located in the corner of the obtained Pareto front (therefore possessing good performance in two aspects). The comparison of aerodynamic performance is displayed in
Figure 10 in terms of pressure coefficient contour and
Figure 11 in terms of two-dimensional comparison of pressure distribution and laminar length in form of skin friction drag coefficient, respectively. From
Figure 10, a shock with less strength and more rearward position is shown in optimized case that gives explanation of the better pressure drag performance. The quantitative improvement in laminar-related performance is also shown in
Figure 11.
5. Conclusions
The paper illustrates an approach that uses the importance in extracting information from objective space and applying it to improve optimization process via the formulation of design requirement. NLF optimization often meets with challenge of the treatment to more than one objective in terms of formulation of design requirement, which has been rarely done in a problem-oriented way by many studies. The study of new weight’ combination has its direct significance in NLF-related engineering problems. It is equivalently the problem of distributing weight among all the concerned aerodynamic performances. The analysis is mainly on the comparison between improved weight distribution and traditional weight distribution (i.e., 50-50 format of design requirement).
The method proposed by this paper uses PCA as the data manipulator that verifies the eigenvectors of objective space that is irrelevant to sampling strategy and has significance in understanding the nature of the problem at hand. Application of this information is finished in two cases of the paper:
proof-of-concept study 1 analyzes the objective space of skin friction and pressure drag force and compared the final effect brought by conventional and improved design requirement in terms of Pareto front formation. proof-of-concept study 2 analyzes the objective space of skin friction, pressure drag, and a third aerodynamic performance of pressure loss. The case illustrates that the introduction of a third aerodynamic performance into the creation of objective space can guide the formulation of two-objective optimization design requirement and get accessory improvement in the comprehensive aerodynamic performance.
In conclusion, the reward for the extraction and use of objective space information not only influences the acceleration of Pareto front marching (as in proof-of-concept study (1)), but also provides a formulation of design requirement that reflects the nature of the problem at hand (as in proof-of-concept study (2)). The reconsideration of design requirement formulation in objective space analysis in this paper has potentials to be extended to other multi-objective optimizations. The significance of the paper is that the information of objective space can be injected in the optimization and its result will inspire the discovery of certain pattern hidden behind the nature of the problem.