Research Article  Open Access
Jing Li, Cong Wang, Huiting Bian, "Aerodynamic Design Optimization of Transonic NaturalLaminarFlow Airfoil at Low Reynolds Number", Mathematical Problems in Engineering, vol. 2020, Article ID 5812129, 22 pages, 2020. https://doi.org/10.1155/2020/5812129
Aerodynamic Design Optimization of Transonic NaturalLaminarFlow Airfoil at Low Reynolds Number
Abstract
The position and size of laminar separation bubble on airfoil surfaces exert a profound impact on the efficiency of transonic naturallaminarflow airfoil at low Reynolds number. Based on the particle swarm algorithm, an optimization methodology in the current work would be established with the aim of designing a high and robust performance transonic naturallaminarflow airfoil at low Reynolds number. This methodology primarily includes two design processes: a traditional deterministic optimization at ondesign point and a multiobjective of uncertaintybased optimization. First, a multigroup cooperative particle swarm optimization was used to obtain the optimal deterministic solution. The crowing distance multiobjective particle swarm optimization and the nonintrusive polynomial chaos expansion method were then adopted to determinate the Paretooptimal front of uncertaintybased optimization. Additionally, the transition model was employed to predict the laminarturbulent transition. Regarding to the established optimization methodology, a propeller tip airfoil of solar energy unmanned aerial vehicle was finally designed. During optimization processes, the minimized pressure drag was particularly chosen as the optimization objective, while the friction drag increment served as a constraint condition. The deterministic results indicate that the optimized airfoil has a good ability to control the separation and reattachment positions, and the pressure drag can be greatly reduced when the laminar separation bubble is weakened. The multiobjective results show that the uncertaintybased optimized airfoil possesses a significant robust performance by considering the uncertainty of Mach number. The findings evidently demonstrate that the proposed optimization methodology and mathematical model are valuable tools to design a highefficiency airfoil for the propeller tip.
1. Introduction
The solar powered unmanned aerial vehicle (UAV), generally flying at an altitude of 10–40 km, where it is with low atmospheric density, has attracted wide attention due to its significant practical engineering application in recent years [1, 2]. Nevertheless, note that the low atmospheric density existing in high altitude environment could result in the Reynolds number range of the propeller airfoils falling within 10^{5}–10^{6} [3]. Subsequently, the aerodynamic performance of airfoil could be deteriorated at such low Reynolds number (LRN), and then, the propeller efficiency would be plummeted [4, 5]. As such, inspired by the importance and urgent need of lowspeed flight at LRN, comprehensive studies have been dedicated to improve the aerodynamic efficiency of its propeller [6, 7]. Zhang et al. carried out the optimization design of E387 and PSU 94–097 airfoils at and , respectively [8]. The designed airfoils showed better aerodynamic characteristics predicted through the Xfoil code than the original ones. Ma et al. established a hierarchical multiobjective optimization platform based on Xfoil code [9]. Their results indicated that three new airfoils optimized from LRN highlift airfoils have excellent aerodynamic characteristics in both ondesign and offdesign states. Wang et al. accomplished the multiobjective optimization design for SD7073 airfoil and proposed a significant design concept that is to change boundary layer properties via controlling the transition positions [10]. Chen et al. discovered the low lifttodrag ratio for the obtained tandem airfoil at Re = 220,000 that could be overcome under the lower lift coefficient, and it particularly demonstrates a manifestation at a high angle of attack [11].
These above new optimized naturallaminarflow (NLF) airfoils presented the outstanding performance in different aspects for lowspeed flight at LRN, but they can scarcely be implemented in the propeller tip. The reason for this is that the velocity nearby the blade tip would be high subsonic or transonic, while the velocity of its root region is rather low [12, 13]. Up to now, there are very few transonic LRN airfoils that are available for the blade tip of solar powered UAV. Undoubtedly, it gives rise to more difficulties in designing transonic LRN airfoil compared with lowspeed NLF one. Moreover, LSB will inevitably appear on one or both surfaces of the airfoil, which are mainly responsible for increasing sensitivity of aerodynamic performance as the Reynolds number decreases to and even below 500,000 [14–16]. For this, it is necessary to employ the active flow control methods to completely eliminate the laminar separation [17, 18]. Furthermore, it has a high probability of developing shock waves on the upper airfoil surface at transonic speed flight. These phenomena would subsequently produce the strong adverse pressure gradients and compressibility effect that greatly influence the size and position of LSB and enhance the design complexity [19, 20]. For another, the fluctuations of the flight conditions, leading to the instabilities of the sizes and locations of separation bubbles, also cause more difficulties [21].
According to these previous studies [12–21], it is badly in need of putting forward an effective optimization methodology to design the transonic NLF airfoil at LRN towards achieving direct application in the blade to enhance the propeller’s efficiency. The particle swarm optimization (PSO) algorithm has been widely used in aerodynamic optimization because of its excellent optimalsearching ability [22–24]. Khurana et al. built an adaptive mutationPSO (AMPSO) method in use for redesigning airfoil at flight envelopes encompassing lowtohigh Mach numbers [25]. Their results illustrated that the AMPSO can seek out the optimal airfoil with minimal computational resources. Tao et al. proposed an improved PSO algorithm on the basis of centroidal Voronoi tessellations (CVTs) and fitness distance ratio (FDR) [26]. The optimization of single point and robust showed that the optimized airfoil and wing both have outstanding characteristics at ondesign point. Masoumi et al. developed a hybrid algorithm combining the big bangbig crunch and PSO, which can be applied for reverse and direct optimization approaches [27]. And this proposed algorithm has proper accuracy and convergence speed that facilitates to reach the desired airfoil geometry. However, it should be mentioned that most of the above PSO algorithms would fall into a local optimum. The main reason is the evolution information for individual guidance only coming from the single group. Under such conditions, the evolutionary ability of algorithm could be declined and thereupon gives rise to an algorithm premature in local optimum [28–30]. Of particular interest is the highorder nonlinear feature of the aerodynamic design space that certainly requires the algorithm to have a strong global optimization capability. Moreover, massive computational consumption of aerodynamic characteristics needs the algorithm to have a faster convergence speed [31, 32].
With respect to the significant role of singleobjective PSO algorithms in deterministic optimization, the high robustness airfoil should be taken into consideration, in which the performances have less sensitivity in response to flight conditions’ uncertainty. Usually, the mean and variance of the characteristics were treated as the indicators for the uncertainty design. However, the mean and variance always were converted into a singleobjective optimization by a linearweight combination [33, 34]. So, the optimization results would be influenced by the treatment of weight coefficients. In this study, the mean and variance are treated as two objectives to get over the disadvantage of the weight method. In line with previous study, the multiobjective particle swarm optimization (MOPSO), employed in many multiobjective aerodynamic optimization, tried to get the Paretooptimal front [35]. Azabi et al. proposed an interactive MOPSO with the benefits of reducing computational time and avoiding the complexity of postdata analysis [36]. Huang et al. constructed an improved MOPSO with principal component analysis methodology, enhancing the convergence of Pareto front [37]. The level of nondominanted Pareto front was effectively improved, and the multiobjective requirements of a widebody aircraft were satisfied. However, the number of the Paretooptimal front will grow dramatically for multiobjective aerodynamic design. The reasonable dominant ranking and selecting for Paretooptimal sets are essential factors in enhancing the capability of the MOPSO [38].
In present work, an optimization methodology based on PSO algorithm was proposed. It aims to design a high and robust performance transonic NLF airfoil at LRN. This article is organized as follows: in Section 2.1, the PSO methodology is built up, and the algorithm performances are validated by the aerodynamic design of a supercritical airfoil; the numerical simulation and airfoil parameterization methods are given in Sections 2.2 and 2.3, respectively; in Section 3, deterministic optimization and multiobjective robust design are conducted to get an optimal airfoil to meet the design requirements; and in Section 4, several important design conclusions are summarized by the comparison of original and newly optimized airfoil.
2. Methodology
2.1. Optimization Methodology Based on the PSO Algorithm
2.1.1. Multigroups Cooperative Particle Swarm Optimization Algorithm (MCPSO)
To avoid algorithm premature, a MCPSO algorithm based on the biological systems symbiosis is put forward in this section. The populations are divided into several subgroups. Every subgroup obtains the evolution information from the group itself and the other symbiotic subgroups. The symbiotic relationships among subgroups are described by a masterslave structure. The slave subgroups and the master one separately undertake the global and local searches.
The basic principles of this algorithm are conducted as follows: the populations are divided into four subgroups, marked as the 1st, 2nd, 3rd, and 4th subgroups. The 1st subgroup is chosen as the master subgroup, which is mainly in charge of local search, and the others are specified as the slave subgroups for the global search. The 1st subgroup is to adopt the standard PSO, while the others are to adopt the improved PSO algorithms accomplished in our prior studies [39, 40]. The speed and position of the 1st subgroup is updated according to the following equation:where and are the speed and position of the particle, respectively; and are the local and global best particles of the 1st subgroup, respectively; is the global best particle of the other subgroups; and are the learning factors; and is the migration factor.
The ability to search the optimum is verified by designing the supercritical airfoil NASASC (2) 0712. The ondesign point is , and the objective is minimizing the total drag at this cruise state. The following optimized mathematical model is built up for this design problem:where is the maximum thickness of the airfoil.
In this testing aerodynamic design, the characteristics are obtained by solving the Reynoldsaveraged Navier–Stokes (RANS) equations [41]. The airfoil parameterization method adopted in this design will be displayed in Section 2.3. There are five design variables for the upper or lower surfaces, ten design variables in all for each airfoil. The cell number of the grids is 29,000.
The standard PSO is marked as GPSO. The population size of the two algorithms (GPSO and MCPSO) counts 60 particles, and the convergence condition is the maximum iteration which is set as 50. The damping strategy is utilized for the boundary treatment method [42].
Aerodynamic characteristics for the original and optimized airfoils at ondesign point are shown in Table 1. The drag coefficient of MCPSO is smallest and decreases 32.8 counts by comparing with the original one. The lifttodrag ratio increases by 24.5%.

The optimization results of GPSO and MCPSO are shown in Figure 1. The convergent curves indicate that MCPSO has less probability of falling into the local optimum than GPSO. One thing should be mentioned that the fitness of MCPSO could attain 0.13505 in the 25^{th} iteration, while such value only could be achieved by GPSO in its 40^{th} iteration. Furthermore, the GPSO falls into the local optimum at the iteration end, compared with the better fitness (0.013394) of the MCPSO. Seen from the surface pressure distributions presented in Figure 1(b), the strong shocks disappear on the upper surface of both two optimized airfoils. Notice that the suck peak of optimized surface pressure distribution is higher than that of the original one. The results indicate that the convergence speed of MCPSO is faster than that of GPSO, and the better optimal solution could be found in the design space quickly and accurately.
(a)
(b)
2.1.2. MOPSO Algorithm Based on CrowdingDistance Sorting
The MOPSO proposed in the current study is based on the crowdingdistance sorting, and the update and control of Pareto sets are as follows:(1)Assume that the maximum size of the Pareto sets named as is M, and the size of in the t iteration is . There will be new Pareto solutions in the following iteration, which puts together with into another temporary Pareto sets named as .(2)If the two solutions have the same objectives, randomly eliminate one of the duplicate solutions in the temporary Pareto sets .(3)The size of will be , when deleting the dominated solutions in the temporary Pareto sets .(4)Sort the populations of in the ascending order according to the crowing distance, marked as . The crowingdistance computation is as follows: where refers to the crowing distance of the i^{th} population; refers to the m^{th} objective function value of the i + 1^{th} individual; and are the maximum and minimum values of the m^{th} objective functions, respectively; and and are the numbers of objectives and set . Herein, the crowingdistance computation requires sorting the population corresponding to each objective function value in ascending order.(5)The sorting sets will be chosen as the new in next iteration if . Otherwise, the exact top M population members from the will be chosen as the new .
In general, the optimization results of MOPSO can be verified by generational distance (GD) and spacing metric (SP). The GD can represent the distance between the nondominated set and true Paretooptimal front, of which the formula is as follows:where n is the size of nondomination set and is the smallest distance between the i^{th} individual and true Paretooptimal front. Notice that the smaller the value of GD is, the better the convergence is.
The SP is used to calculate the distribution uniformity, and its formula is as follows:
A smaller SP denotes a better uniformity of the sets. Thus, three typical multiobjective functions labeled by ZDT1, ZDT2, and ZDT3 are to verify the performance: ZDT1: ZDT2: ZDT3:
Figure 2 shows the Paretooptimal fronts of these three functions obtained from MOPSO. The newly designed fronts are extremely consistent with the real ones. Moreover, as shown in Table 2, the GDs of MOPSO are better than those of NSGAII [43] for all three functions. As for SPs, it is hard to distinguish between the two algorithms to find which one performs better.
(a)
(b)
(c)

2.1.3. Uncertainty Quantification by Polynomial Chaos Expansion (PCE)
The main difficulties for uncertainty quantification are the excessive computational cost and unsatisfactory accuracies of the mean and variance. For instance, the widely used Monte Carlo simulation (MCS) [44] method basically requires at least 10^{4} samples to get sufficiently accurate estimate of the mean and variance, which certainly causes a much high cost of computation. In the aspect of efficient and accurate uncertainty quantification, the PCE is a promising method for accurate estimation of mean and variance of the objectives [45, 46]. A noninstrusive PCE method is adopted in this part, which is a questionanswer system and unnecessary to modify the original governing equations [47].
This formula assumes that the random variable obeys a standard normal distribution, and thus, p polynomial order of a general random process can be expressed as follows:where is the Hermite polynomial, are the vectors consisting of n Gaussian variables, , and n is the number of random variables.
Then, coefficients can be solved by the following formula:where the number of collocation points is indicated by and and are the Gaussian–Hermite quadrature coefficient and point, respectively.
According to the orthogonality of Gaussian–Hermite polynomials, the mean and variance of are presented as follows:
The accuracy and efficiency for PCE and MCS approaches are compared through the following function:where X is assumed to be normally distributed with and . The random variables are assumed as , and is expressed as .
The real mean and variance of this function would be calculated by the following equations:
Figure 3 plots the comparisons of estimated mean () and variance () between PCE and MCS methods. The sample numbers for MCS fetch two values, i.e., 100 and 500, that are utilized to appraise its estimation accuracy of the stochastic process. The results demonstrate that the mean and variance will become convergent with the increasing order P. And the estimation accuracy of mean is slightly higher than that of variance, particularly for the P value within 2–4. Additionally, the estimated and by the PCE method both fall within 0.5 percent error bands, when the polynomial order P is up to 4. Meanwhile, the estimated accuracy of agrees well with that of MCS using 500 samples. But the estimated accuracy of becomes closer to that of MCS using 500 samples with the polynomial order P increasing from 4. Thereupon, our results demonstrate that the PCE method is reasonable for the uncertain quantification with respect to its reliable accuracy and efficiency. And the same conclusion was also derived by the former study that verified the estimated results of aerodynamic characteristics by PCE [48]. Figure 3(c) shows the time cost of PCE and MCS approaches varying with the error of variances. When the error of variance is at the same level, the time cost of the MCS approach is much more than that of the PCE approach. Therefore, the PCE method is chosen to estimate the mean and variance of the objectives in present work for its efficiency.
(a)
(b)
(c)
2.2. Numerical Simulation Method for Aerodynamic Characteristics
The transition model [49, 50] is combined with the local variable of SST turbulence model to predict transition. The proposed model incorporated intermittent function with transition experience relationship. So, it can control the distribution of intermittent functions in boundary layers through the relevance functions of transition momentum thickness Reynolds number and afterward dominates turbulence through the intermittent functions. Moreover, since the momentum thickness Reynolds number is related to the strain rate, the calculation of momentum thickness can be prevented, thereby enhancing the prediction accuracy of separation transition [50].
The transition model validation is performed on NLR7301 airfoil applying the same boundary conditions adopted in the experiments. The simulation condition is . The airfoil was tested at the NLR pilot closedcircuit wind tunnel [51]. There are structured grid for the transition model solver, as shown in Figure 4. As displayed in Figure 5, the transition position (TP) and surface pressure distribution obtained from the transition model and experiment data illustrate good consistency.
(a)
(b)
(a)
(b)
(c)
(d)
Also, the validation of capturing LSB is accomplished on E387 airfoil at low speed and . The airfoil performances predicted by transition model are compared to University of Illinois UrbanaCampaign (UIUC) windtunnel data [52] in Figure 6. The results displayed in Figure 6(a) demonstrate that the liftdrag curve is in good agreement with the experimental one. As presented in Figure 6(b), both the location and size of LSB obtained from CFD are coincident with the experiment data.
(a)
(b)
Consequently, these phenomena make us to conclude that the numerical simulation technology of transition boundary layer is accurate and reliable for the aerodynamic design optimization of transonic NLF airfoil at LRN.
2.3. Airfoil Parameterization Method
The disturbing basis function, based on the class/shape function transformation (CST) [53], is introduced in airfoil parameterization. The parameterization method constructed here can not only minimize the errors caused by fitting the initial parameters but also preserve the advantages of the original CST parameterization method.
The original CST method is shown in the following equation:where and are the xcoordinates and ycoordinates of the airfoil, respectively; is the ycoordinate of the trailing edge point; is the design variables; and .
The disturbance is added to the initial design variable , and then, equation (14) can be written aswhere and are the separable ycoordinates of new and original airfoils and are the disturbances of ycoordinates, which can be denoted as
According to the aforementioned equations, the new ycoordinates can be treated as two parts, i.e., the original ycoordinates () and disturbances (). Thus, equation (16) can be denoted as the disturbing basis function based on CST.
3. Aerodynamic Design Problem
3.1. Optimization Design Definition
As mentioned above, the size and position of LSB and transition location have a great impact on aerodynamic characteristics of the transonic NLF airfoil at LRN. So, these two factors should be considered during the design process. The classical laminar bubble separation theory [54–56] indicates that the adverse pressure gradient of the upper surfaces will promote the laminar separation, and the transition will take place in free shear layer above the airfoil surface, of which the location is downstream of the SP. Consequently, the LSB would become smaller with the forward SP and RP, and the pressure drag would decrease. Conversely, the friction drag would increase by the forward TP. It is of particular interest to point out that the friction drag decreases dramatically by the backward TP. Nevertheless, an extra length of laminar flow region would result in the excessive adverse pressure gradient at the pressure recovery section, thus yielding a larger LSB on upper surface. Accordingly, the pressure drag sharply increases because of the larger LSB. Regarding to the above analysis, it should keep a balance between these two factors in the aerodynamic design of transonic NLF airfoil at LRN.
The flight condition of the propeller tip airfoil, utilized in a solar energy UAV, is transonic at LRN. As for this, the ondesign point is assumed as . For the original airfoil SD8000, its CFD results show that there is a large LSB on the upper surface and the proportion of pressure drag is 75%. Therefore, the decrease in pressure drag would correspondingly lead to a reduction in total drag. And it can be achieved by weakening the LSB. However, the friction drag would increase when the TP moves forward.
3.2. Deterministic Optimization Results and Discussion
The following optimized mathematical model is built up for the deterministic design problem:where is the lift coefficient, is the pressure drag coefficient, is the maximum thickness of airfoil, and and are the friction drag coefficients of the optimized and original airfoils, respectively. The objective of this model is to reduce the pressure drag coefficient. There are some constraints, including the growing proportion of friction drag coefficient less than 10%, the maximum thickness of airfoil greater than 8.92%, and the lift coefficient being 0.5.
There are five design variables for the upper or lower surfaces, a total of ten design variables for each airfoil. The cell number of the grids is [57], of which the distributions for the leading edge and trailing edge are demonstrated in Figures 7(a) and 7(b). With respect to the grid convergence for the original airfoil presented in Figure 7(c), when the circumferential grid number is more than 600, the variations of the total and pressure drags are within 0.0008.
(a)
(b)
(c)
In the deterministic optimization, the algorithm population in each iteration is set as 32 particles, and the maximum iteration step is set as 60. Table 3 lists the aerodynamic characteristics for both the original and optimized airfoils. It can be summarized that the total drag and pitching moment coefficients separately decrease by 48 counts and 26.6%, while the friction drag coefficient and lifttodrag ratio increase by 7.4% and 28.7% at ondesign point, respectively. The friction drag increases because of the TP moving forward, and the pressure drag reduces because of the smaller size of LSB, as displayed in Figure 8(b). From the geometric characteristics enumerated in Table 4, it can be noticed that there is a backward location of the maximum thickness on optimized airfoil. Its maximum camber is smaller than that of the original one. And the same situation happens to the leading edge radius, of which the value for optimized airfoil is also smaller.

(a)
(b)

Figure 8 displayed two airfoils and surface pressure distributions at ondesign point. Compared with the original airfoil, the optimized one with a flat upper surface has a smaller camber before the 60% chord and larger camber after the 60% chord. There is a stronger inverse pressure gradient roughly existing in 5% to 20% chord for the optimized pressure distribution, and then, this gradient becomes weaker for 50% to 75% chord. Moreover, the suction peak of the optimized airfoil obviously becomes larger. As marked in Figure 8, the SP and RP of the optimized airfoil are moved forward at the position of 15% and 60% chord, respectively. Figure 9 shows the Mach number distribution contours and the streamlines of the original and optimized airfoils at ondesign point. The separation bubble on upper surface moves forward and becomes smaller compared with that of the original one. We assume that the transition will take place in free shear layer above the airfoil surface, of which the location is the downstream of SP. The forward TP results in a slight increase of friction drag. Simultaneously, the small LSB results in a considerable decrease of pressure drag. Moreover, the height of the bubble is also smaller. It can be concluded from Figures 8 and 9 that the strong inverse pressure gradient on the leading edge could lead to the forward SP, and the early recover pressure could result in the forward RP. The forward SP and RP minimize the LSB to reduce the pressure drag. Nevertheless, the forward LSB may be conducive to the formation of a forward transition position TP, subsequently causing the increase of friction drag. The lifttodrag ratio curves of the original and optimized airfoils are shown in Figure 10. For the optimized airfoil, there is a low drag region over the range of lift coefficient from 0.3 to 0.52. It is caused by the sharp decrease in the pressure drag, as clearly presented in Figure 10.
(a)
(b)
(a)
(b)
It can be seen from Figure 10 that the drag of optimized airfoil is larger than that for the original one, when C_{L} = 0.2. By comparing the streamline distributions of the offdesign point (C_{L} = 0.2) for original airfoil in Figure 11(a), it is obviously noticed that there is a large flow separation in wake region of the optimized airfoil, which will increase the pressure drag in Figure 11(b). Furthermore, its skin friction coefficient increases in the trailing edge of the optimized airfoil in Figure 11(c), which will significantly give rise to the growth of the friction drag. Therefore, both the increases of pressure and friction drags result in increasing the total drag.
(a)
(b)
(c)
Figure 12 shows the comparison of aerodynamic characteristics between the original and optimized ones. The pressure drag of optimized airfoil shows a considerable decrease over the Mach number ranging from 0.70 to 0.76. The friction drag of optimized airfoil shows a slight increase, which meets the constraints at ondesign point. Note that the maximum increase of the friction drag is only up to 5 counts. It can be summarized from the drag divergence curves that the total drag becomes decreased significantly with the substantial decrease of the pressure drag. The optimized airfoil has the smaller drag coefficients over the range of Mach numbers from 0.70 to 0.76. The SP and RP varying with the Mach numbers are displayed in Figure 12(d). Over the Mach number ranging from 0.70 to 0.76, the LSB of the optimized airfoil is smaller than that of the original one and thus leads to a larger reduction of the pressure drag. However, when Ma > 0.75, the RP cannot be detected on the original airfoil, leading to the sharp rise of the pressure drag. Therefore, the varying trend of the pressure drag is mostly influenced by the size and position of LSB over the studied Mach number range.
(a)
(b)
(c)
(d)
Figure 13(a) shows the pressure distributions of original and optimized airfoils varying with the Mach number. It can be seen that the obvious flat pressure distribution appears around the position of separation bubble. The results also present that the strong inverse pressure gradient on the leading edge leads to a forward SP on the optimized airfoil and an inverse pressure gradient around the RP. These characteristics indicate that the laminar separation of optimized airfoil coincides with the classical LSB theory. The SP of original airfoil moves forward when the Mach number becomes larger, but the RP is difficult to be detected on the upper surface. These characteristics are consistent with the trailing edge separation bubble theory.
(a)
(b)
3.3. Multiobjective Robust Results and Discussion
Although the optimized airfoil of deterministic has an outstanding performance at ondesign point, the pressure and total drags both have the sharp increases, when Ma > 0.75. Therefore, the fluctuation of flight conditions would cause a sudden drop in the airfoil’s efficiency. A multiobjective robust design is carried out in this section aiming to reduce the performances’ sensitivity to the uncertainties of flight conditions. The PCE method is employed to obtain the mean and variance of the objectives. The four Gaussian–Hermite integration points are introduced in the nonintrusive PCE method to calculate the PCE coefficients, and then, a 3order PCE model is set up for the robust design. The population of MOPSO is set as 40, and the maximum iteration is 20.
The Pareto optimal solution is chosen to deal with the mean and variance as two objectives rather than converting them into a single objective. We assume Mach number to be normally distributed . The following optimized mathematical model is built up for the multiobjective design problem:where and are the mean and the variance of the pressure drag, respectively. and are the separate mean friction drags of the optimized and original ones.
We focus on the marked individuals of the Paretooptimal front named as optimization_2, as shown in Figure 14. The deterministic optimized airfoil named as optimization_1 also plots in Figure 14. The optimization_2 has a smaller camber before the 40% chord compared with optimization_1. However, the camber distributions for these two airfoils are aligned with each other after 40% chord.
(a)
(b)
Table 5 lists the aerodynamic force coefficients of the optimization ones at ondesign point. It can be summarized that the total and pressure drag coefficients of optimization_2 are smaller than those of optimization_1. The lifttodrag ratio of optimization_2 is better than that of optimization_1. Table 6 displays the geometric characteristics for two optimized airfoils. Compared with these of optimization_1, the positions of both the maximum thickness and camber of optimization_2 move forward slightly. And its leading edge radius becomes smaller demonstrated in Table 6.


Figure 15 shows the comparisons of aerodynamic characteristics between deterministic and robust optimizations. When Ma < 0.72, optimization_2 has the larger pressure drag coefficients than the optimization_1. However, the contradictory trend is found for them, when Ma > 0.72. With a tradeoff of the pressure drag, the varying trend of pressure drag becomes gentler over the range of Mach number from 0.70 to 0.76. The friction drag coefficients of optimization_2 are smaller than the those for optimization_1. From the drag divergence curves of the optimized airfoils, we can see that the rapid change of drag calms down through the robust design. For the performance of optimization_2, although it has a tradeoff of the pressure drag, its friction drag acquires a slighter drop. Therefore, the airfoil has a wide range of lower drag. The SP and RP varying with the Mach numbers are displayed in Figure 15(d). Two optimizations have the similar LSB length; however, the SP and RP of optimization_2 are both forward compared with those of optimization_1.
(a)
(b)
(c)
(d)
Figure 16 shows the aerodynamic characteristics of optimized airfoils at ondesign point. There is a stronger inverse pressure gradient occurring in 5% to 20% chord on optimization_2, and then, a weaker inverse pressure gradient can be found from 40% to 60% chord. Furthermore, the suction peak of the optimization_2 becomes larger than that of optimization_1. For clarity, the SPs and RPs are also marked in Figure 16. The SP and RP of optimization_2 would move forward at the positions of 11% and 58% chord, respectively. As mentioned before, the identical circumstance also happens here; that is, the strong inverse pressure gradient on the leading edge and the forward SP both lead to a forward TP. The position and length of LSB also can be seen from surface friction distributions, as shown in Figure 16. The separation bubble of optimization_2 on upper surface moves forward and becomes smaller than that of optimization_1. Figure 17 displays the aerodynamic characteristics of optimized airfoils at . There is an obvious lower drag region of optimization_2 over the range of lift coefficient from 0.1 to 0.55. Additionally, with the Mach number increasing, the suction peak becomes smaller and the pressure recovery is postponed.
(a)
(b)
(a)
(b)
The pitching moment curves of the original and optimized airfoils are shown in Figure 18. It can be seen that the pitching moments of the optimized airfoils are smaller than that of the original one. Additionally, as shown in Figures 8, 14, and 16, both the maximum camber and after load of the optimized airfoils are also smaller than those of the original one, and thereby, they all will result in the decrease in the pitching moment.
(a)
(b)
Regarding to all above discussions, the mulitobjective robust design could produce an excellent airfoil. The optimized surface pressure distribution not only keeps the low drag at ondesign point but also catches the robustness of low drag to overcome the fluctuation of flight condition.
4. Conclusions
In this paper, we carry out the aerodynamic design optimization of a propeller tip airfoil in a solar energy UAV. The flight condition of this airfoil is treated as transonic at LRN. Under such conditions, we develop an optimization methodology based on the PSO algorithm applied to the design of a transonic NLF airfoil at LRN. This new proposed methodology briefly contains two optimization design processes. They are traditional deterministic optimization at ondesign point on the basis of MCPSO and multiobjective robust design based on MOPSO. The nonintrusive PCE method is adopted to get sufficiently accurate estimations of the mean and variance. It facilitates to the fact that the optimal deterministic solution and Paretooptimal front of aerodynamic design can be captured in the space quickly and accurately. From the findings of deterministic optimizations and multiobjective robust design, we obtained some important conclusions for transonic NLF airfoil at LRN, as follows:(1)We construct a new optimization model; that is, the minimized pressure drag is typically treated as the optimization objective, while the friction drag increment is a constraint condition. According to the deterministic optimization and robust design results, this proposed model is reliable and efficacious to design a highperformance airfoil of the propeller tip.(2)The pressure drag is mostly affected by the size and position of LSB. The deterministic optimized pressure distribution can lead to a small and forward LSB. The SP and RP of the airfoil move forward by virtue of the deterministic optimized pressure distribution.(3)The multiobjective robust design can weaken the sensitivity of airfoil performance to the flight fluctuation. The proper inverse pressure gradient at the leading edge and early pressure recover show the good abilities to control the LSB.(4)The proposed optimization methodology is an efficient and reasonable approach for the aerodynamic design of transonic NLF airfoil at LRN.
Data Availability
The data used to support the findings of this study are available from the corresponding author upon request.
Conflicts of Interest
The authors declare that there are no conflicts of interest.
Acknowledgments
The authors gratefully acknowledge financial support from the China Scholarship Council. The authors are also grateful for computing support by the Supercomputing Center of Zhengzhou University. This research was sponsored by the National Science Foundation of China (NSFC) (grant number 11602226), University Science Project of Henan Province (grant number 16A130001), and China Postdoctoral Science Foundation (grant number 2019M652579).
References
 S. Morton, L. Scharber, and N. Papanikolopoulos, “Solar powered unmanned aerial vehicle for continuous flight: conceptual overview and optimization,” in Proceedings of the 2013 IEEE International Conference on Robotics and Automation, pp. 766–771, Karlsruhe, Germany, May 2013. View at: Google Scholar
 S. Leutenegger, M. Jabas, and R. Y. Siegwart, “Solar airplane conceptual design and performance estimation,” Unmanned Aerial Vehicles, vol. 61, no. 1–4, pp. 545–561, 2010. View at: Publisher Site  Google Scholar
 J. Brandt and M. Selig, “Propeller performance data at low Reynolds numbers,” in Proceedings of the 49th AIAA Aerospace Sciences Meeting including the New Horizons Forum and Aerospace Exposition, Orlando, FL, USA, 2011. View at: Google Scholar
 W. Chen and L. Bernal, “Design and performance of low Reynolds number airfoils for solarpowered flight,” in Proceedings of the 46th AIAA Aerospace Sciences Meeting and Exhibit, Reno, NV, USA, 2008. View at: Google Scholar
 R. W. Deters, G. K. A. Krishnan, and M. S. Selig, “Reynolds number effects on the performance of smallscale propellers,” in Proceedings of the 32nd AIAA Applied Aerodynamics Conference, Atlanta, GA, USA, 2014. View at: Google Scholar
 K. Wang, Z. Zhou, X. Zhu, and X. Xu, “Aerodynamic design of multipropeller/wing integration at low reynolds numbers,” Aerospace Science and Technology, vol. 84, pp. 1–17, 2019. View at: Publisher Site  Google Scholar
 B. R. Rakshith, S. M. Deshpande, R. Narasimha, and C. Praveen, “Optimal lowdrag wing planforms for tractorconfiguration propellerdriven aircraft,” Journal of Aircraft, vol. 52, no. 6, pp. 1791–1801, 2015. View at: Publisher Site  Google Scholar
 S. Zhang, H. Li, and A. A. Abbasi, “Design methodology using characteristic parameters control for low Reynolds number airfoils,” Aerospace Science and Technology, vol. 86, pp. 143–152, 2019. View at: Publisher Site  Google Scholar
 R. Ma, B. Zhong, and P. Liu, “Optimization design study of lowreynoldsnumber highlift airfoils for the highefficiency propeller of lowdynamic vehicles in stratosphere,” Science China Technological Sciences, vol. 53, no. 10, pp. 2792–2807, 2010. View at: Publisher Site  Google Scholar
 K. Wang, X. Zhu, Z. Zhou, and X. Xu, “Studying optimization design of low reynolds number airfoil using transition model,” Journal of Northwestern Polytechnical University, vol. 33, pp. 580–587, 2015. View at: Google Scholar
 F. Chen, J. Yu, and Y. Mei, “Aerodynamic design optimization for low reynolds tandem airfoil,” Proceedings of the Institution of Mechanical Engineers, Part G: Journal of Aerospace Engineering, vol. 232, no. 6, pp. 1047–1062, 2017. View at: Publisher Site  Google Scholar
 Z. Q. Zhu, X. L. Wang, Z. C. Wu, and Z. M. Chen, “Aerodynamic characteristics of small/micro unmanned aerial vehicles and their shape design,” Acta Aeronauticaet Astronautica Sinica, vol. 27, pp. 353–364, 2006. View at: Google Scholar
 F. Li and P. Bai, Aircraft Aerodynamic at Low Reynolds Number, China Aerospace Press, Beijing, China, 2017.
 R. M. Hicks and S. E. Cliff, An Evaluation of Three TwoDimensional Computational Fluid Dynamics Codes Including Low Reynolds Numbersand Transonic Mach Numbers, NASA, Silicon Valley, CA, USA, 1991.
 R. Hain, C. J. Kähler, and R. Radespiel, “Dynamics of laminar separation bubbles at lowreynoldsnumber aerofoils,” Journal of Fluid Mechanics, vol. 630, pp. 129–153, 2009. View at: Publisher Site  Google Scholar
 M. Jahanmiri, Laminar Separation Bubble: Its Structure, Dynamics and Control, Chalmers University of Technology, Gothenburg, Sweden, 2011.
 D. Kamari, M. Tadjfar, and A. Madadi, “Optimization of SD7003 airfoil performance using TBL and CBL at low reynolds numbers,” Aerospace Science and Technology, vol. 79, pp. 199–211, 2018. View at: Publisher Site  Google Scholar
 C. Chen, R. Seele, and I. Wygnanski, “Flow control on a thick airfoil using suction compared to blowing,” AIAA Journal, vol. 51, no. 6, pp. 1462–1472, 2013. View at: Publisher Site  Google Scholar
 M. Drela, “Transonic lowreynolds number airfoils,” Journal of Aircraft, vol. 29, no. 6, pp. 1106–1113, 1992. View at: Publisher Site  Google Scholar
 R. Placek and P. Ruchała, “The flow separation development analysis in subsonic and transonic flow regime of the laminar airfoil,” Transportation Research Procedia, vol. 29, pp. 323–329, 2018. View at: Publisher Site  Google Scholar
 H. Zhao, Z. Gao, Y. Gao, and C. Wang, “Effective robust design of high lift NLF airfoil under multiparameter uncertainty,” Aerospace Science and Technology, vol. 68, pp. 530–542, 2017. View at: Publisher Site  Google Scholar
 A. Banks, J. Vincent, and C. Anyakoha, “A review of particle swarm optimization. Part II: hybridisation, combinatorial, multicriteria and constrained optimization, and indicative applications,” Natural Computing, vol. 7, no. 1, pp. 109–124, 2008. View at: Publisher Site  Google Scholar
 S. Sun and J. Li, “A twoswarm cooperative particle swarms optimization,” Swarm and Evolutionary Computation, vol. 15, pp. 1–18, 2014. View at: Publisher Site  Google Scholar
 Q. Luo and D. Yi, “A coevolving framework for robust particle swarm optimization,” Applied Mathematics and Computation, vol. 199, no. 2, pp. 611–622, 2008. View at: Publisher Site  Google Scholar
 M. Khurana and K. Massey, “Swarm algorithm with adaptive mutation for airfoil aerodynamic design,” Swarm and Evolutionary Computation, vol. 20, pp. 1–13, 2015. View at: Publisher Site  Google Scholar
 J. Tao, G. Sun, X. Wang, and L. Guo, “Robust optimization for a wing at drag divergence Mach number based on an improved PSO algorithm,” Aerospace Science and Technology, vol. 92, pp. 653–667, 2019. View at: Publisher Site  Google Scholar
 H. Masoumi and F. Jalili, “Inverse and direct optimization shape of airfoil using hybrid algorithm big bangbig crunch and particle swarm optimization,” Journal of Theoretical and Applied Mechanics, vol. 57, no. 3, pp. 697–711, 2019. View at: Publisher Site  Google Scholar
 M. Kohler, M. M. B. R. Vellasco, and R. Tanscheit, “PSO+: a new particle swarm optimization algorithm for constrained problems,” Applied Soft Computing, vol. 85, Article ID 105865, 2019. View at: Publisher Site  Google Scholar
 S. Kiranyaz, J. Pulkkinen, and M. Gabbouj, “Multidimensional particle swarm optimization for dynamic environments,” in Proceedings fo the 2008 International Conference on Innovations in Information Technology, pp. 34–38, Al Lin, United Arab Emirates, December 2008. View at: Google Scholar
 A. P. Engelbrecht, “Particle swarm optimization: global best or local best?” in Proceedings of the BRICS Congress on Computational Intelligence and 11th Brazilian Congress on Computational Intelligence, pp. 124–135, Ipojuca, Brazil, September 2013. View at: Google Scholar
 H. Likeng, G. Zhenghong, and Z. Dehu, “Research on multifidelity aerodynamic optimization methods,” Chinese Journal of Aeronautics, vol. 26, no. 2, pp. 279–286, 2013. View at: Publisher Site  Google Scholar
 H. Han, L. Liang, Y. Shujin, and H. Zhifeng, “Particle swarm optimization with convergence speed controller for largescale numerical optimization,” Soft Computing, vol. 23, no. 12, pp. 4421–4437, 2018. View at: Publisher Site  Google Scholar
 M. Nemec, D. W. Zingg, and T. H. Pulliam, “Multipoint and multiobjective aerodynamic shape optimization,” AIAA Journal, vol. 42, no. 6, pp. 1057–1065, 2004. View at: Publisher Site  Google Scholar
 H. Zhao and Z. Gao, “Uncertaintybased design optimization of NLF airfoil for high altitude long endurance unmanned air vehicles,” Engineering Computations, vol. 36, no. 3, pp. 971–996, 2019. View at: Publisher Site  Google Scholar
 Y. Naranjani and J.Q. Sun, “Multiobjective optimal airfoil design for cargo aircrafts,” in Proceedings of the ASME 2016 International Mechanical Engineering Congress and Exposition, Portland, OR, USA, 2016. View at: Google Scholar
 Y. Azabi, A. Savvaris, and T. Kipouros, “The interactive design approach for aerodynamic shape design optimisation of the aegis UAV,” Aerospace, vol. 6, no. 4, 2019. View at: Publisher Site  Google Scholar
 J. Huang, Z. Zhou, Z. Gao, M. Zhang, and L. Yu, “Aerodynamic multiobjective integrated optimization based on principal component analysis,” Chinese Journal of Aeronautics, vol. 30, no. 4, pp. 1336–1348, 2017. View at: Publisher Site  Google Scholar
 A. de Campos, A. T. R. Pozo, and E. P. Duarte, “Parallel multiswarm PSO strategies for solving many objective optimization problems,” Journal of Parallel and Distributed Computing, vol. 126, pp. 13–33, 2019. View at: Publisher Site  Google Scholar
 J. Li, Z. Gao, J. Huang, and K. Zhao, “Robust design of NLF airfoils,” Chinese Journal of Aeronautics, vol. 26, no. 2, pp. 309–318, 2013. View at: Publisher Site  Google Scholar
 J. Li, Z. Gao, J. Huang, and K. Zhao, “Aerodynamic design optimization of nacelle/pylon position on an aircraft,” Chinese Journal of Aeronautics, vol. 26, no. 4, pp. 850–857, 2013. View at: Publisher Site  Google Scholar
 J. F. Wendt, Computational Fluid Dynamics, Springer, Berlin, Germany, 2009.
 J. LI, Research on the High Performance Aircraft Aerodynamic Design, Northwestern Polytechnical University, Xi’an, China, 2014.
 K. Deb, A. Pratap, S. Agarwal, and T. Meyarivan, “A fast and elitist multiobjective genetic algorithm: NSGAII,” IEEE Transactions on Evolutionary Computation, vol. 6, no. 2, pp. 182–197, 2002. View at: Publisher Site  Google Scholar
 L. Zhang, H. Wu, P. Li, Y. Hu, and C. Song, “Design, analysis, and experiment of multiring permanent magnet bearings by means of equally distributed sequences based Monte Carlo method,” Mathematical Problems in Engineering, vol. 2019, Article ID 4265698, 17 pages, 2019. View at: Publisher Site  Google Scholar
 K. Chikhaoui, D. Bitar, N. Kacem, and N. Bouhaddi, “Robustness analysis of the collective nonlinear dynamics of a periodic coupled pendulums chain,” Applied Sciences, vol. 7, no. 7, p. 684, 2017. View at: Publisher Site  Google Scholar
 H. Zhao, Z. Gao, F. Xu, Y. Zhang, and J. Huang, “An efficient adaptive forwardbackward selection method for sparse polynomial chaos expansion,” Computer Methods in Applied Mechanics and Engineering, vol. 355, pp. 456–491, 2019. View at: Publisher Site  Google Scholar
 M. Eldred and J. Burkardt, “Comparison of nonintrusive polynomial chaos and stochastic collocation methods for uncertainty quantification,” in Proceedings of the 47th AIAA Aerospace Sciences Meeting including The New Horizons Forum and Aerospace Exposition, Orlando, FL, USA, January 2009. View at: Google Scholar
 K. Zhao, Z. Gao, J. Huang, and J. Li, “Uncertainty quantification and robust design of airfoil based on polynomial chaos technique,” Chinese Journal of Theoretical and Applied Mechanics, vol. 46, pp. 10–19, 2014. View at: Google Scholar
 R. B. Langtry and F. R. Menter, “Correlationbased transition modeling for unstructured parallelized computational fluid dynamics codes,” AIAA Journal, vol. 47, no. 12, pp. 2894–2906, 2009. View at: Publisher Site  Google Scholar
 T. H. Nguyen and K. Chang, “Comparison of the pointcollocation nonintrusive polynomial (NIPC) and nonintrusive spectral projection (NISP) methods for the γRθ transition model,” Applied Sciences, vol. 9, no. 7, 2019. View at: Publisher Site  Google Scholar
 B. Ernesto and P. Rita, “Laminar to turbulent boundary layer transition investigation on a supercritical airfoil using the γ–θ transitional model,” in Proceedings of the 40th Fluid Dynamics Conference and Exhibit, Chicago, IL, USA, June 2010. View at: Google Scholar
 M. S. Selig and J. J. Guglielmo, “Highlift low reynolds number airfoil design,” Journal of Aircraft, vol. 34, no. 1, pp. 72–79, 1997. View at: Publisher Site  Google Scholar
 B. Kulfan, “Universal parametric geometry representation method,” Journal of Aircraft, vol. 45, no. 1, pp. 142–158, 2008. View at: Publisher Site  Google Scholar
 M. P. Simens and A. G. Gungor, “Influence of the transition of a laminar separation bubble on the downstream evolution of strong adverse pressure gradient turbulent boundary layers,” European Journal of MechanicsB/Fluids, vol. 67, pp. 70–78, 2018. View at: Publisher Site  Google Scholar
 S. S. Diwan and O. N. Ramesh, “Laminar separation bubbles: dynamics and control,” Sadhana, vol. 32, no. 12, pp. 103–109, 2007. View at: Publisher Site  Google Scholar
 H. P. Horton, Laminar Separation Bubbles in Two and Three Dimensional Incompressible Flow, Queen Mary College, University of London, London, UK, 1968.
 Z. Huan, Z. Gao, C. Wang, and Y. Gao, “Research on the computing grid of high speed laminar airfoil,” Chinese Journal of Applied Mechanics, vol. 2018, no. 4, pp. 351–357, 2019. View at: Google Scholar
Copyright
Copyright © 2020 Jing Li et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.