Next Article in Journal
Research on Game Theory of Air Traffic Management Cyber Physical System Security
Next Article in Special Issue
Aerodynamic Optimization of Transonic Rotor Using Radial Basis Function Based Deformation and Data-Driven Differential Evolution Optimizer
Previous Article in Journal
Speech GAU: A Single Head Attention for Mandarin Speech Recognition for Air Traffic Control
Previous Article in Special Issue
A Multi-Objective Factorial Design Methodology for Aerodynamic Off-Takes and Ducts
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Aerodynamic Performance Uncertainty Analysis and Optimization of a Conventional Axisymmetric Vehicle Based on Parallel Polynomial Chaos Expansions

1
School of Astronautics, Beihang University, Beijing 100191, China
2
Key Laboratory of Spacecraft Design Optimization and Dynamic Simulation Technologies, Ministry of Education, Beijing 100091, China
3
Aircraft and Propulsion Laboratory, Ningbo Institute of Technology, Beihang University, Ningbo 315832, China
*
Author to whom correspondence should be addressed.
Aerospace 2022, 9(8), 396; https://doi.org/10.3390/aerospace9080396
Submission received: 20 May 2022 / Revised: 10 July 2022 / Accepted: 16 July 2022 / Published: 22 July 2022
(This article belongs to the Special Issue Aerodynamic Shape Optimization for Aerospace Engineering Applications)

Abstract

:
In this study, the aerodynamic uncertainty analysis and optimization of a conventional axisymmetric vehicle with an aerodynamic configuration were investigated. The prediction precision of the typical aerodynamic performance estimating methods, namely, engineering estimation and numerical simulation, was compared using the wind tunnel test data of the vehicle. Then, using a modified missile data compendium (DATCOM) software, a high-efficiency and high-precision method was developed, which was applied to analyze and characterize the aerodynamic parameters of the axisymmetric vehicle. To enhance the robustness and reliability of aerodynamic performance, an uncertainty-based design optimization (UDO) framework was established. The design space was scaled by parameter sensitivity analysis, and improved computational efficiency was achieved by developing parallel polynomial chaos expansions (PCEs). The optimized results show that the modified method exhibits high accuracy in predicting aerodynamic performance. For the same constraints, the results of the deterministic design optimization (DDO) showed that compared with the initial scheme, the probability of the controllability-to-stability ratio satisfying the constraint decreased from 98.8% to 72.4%, and this value increased to 99.9% in the case of UDO. Compared with the results of the initial scheme and DDO, UDO achieved a considerable reduction in mean values and standard deviation of aerodynamic performances, which can ensure a higher probability of constraints meeting the design requirements, thereby, realizing a reliable and robust design.

1. Introduction

In the coming decades, physical modeling and highly accurate predictions of aerodynamic performance will be the focus of vehicle research. It can be seen from current research that aerodynamic data prediction of vehicles is always based on deterministic design, which does not consider the errors and uncertainties of physical models, accidental uncertainties caused by natural variables, and errors or uncertainties from test data. The precision of long-range strike weapons and the design of new axisymmetric vehicles require robust and reliable aerodynamic data. Therefore, research should focus on potential uncertainties in aerodynamic prediction and optimization, such as the dimensional deviations between the physical model and the actual vehicle, errors in the numerical model, actual flight conditions, etc.
Several studies have investigated various uncertainty factors in aerodynamic experiments and numerical calculations. Burner et al. [1] supplemented uncertainty analyses with assessments of experimental errors, and the utilization of experimental methods with modern designs, which helped explore the influence of various uncertainties, such as sensors and state variables on the measurement accuracy of model pitch angle in wind tunnel tests. Some academics [2,3,4,5] focused on the influence of shape deviation on aerodynamic performance. Some recent studies have shown that deterministic designs are easily influenced by uncertainties in certain geometric variables or environmental parameters [6]. Moreover, there are many uncertain factors, such as the uncertainty of model assumptions, when using a mathematical model to represent a physical process. Some details of the source and classification of uncertainties in computational fluid dynamics (CFD) simulations [7,8] have been discussed. Zhang et al. [9] focused on quantifying blunt body surface thermal uncertainty, along with flow parameters, and jet geometry in the presence of free-stream perturbations. Andrea et al. explored the impact of epistemic uncertainties and CFD computational errors on aerodynamic performance [10]. These studies, however, have been limited to vehicle components or two-dimensional airfoils. More evidence is needed to determine whether local alterations considerably influence the main aerodynamic performance of vehicles. When we compare airfoil or component design with vehicle shape design, the number of design variables, and the discrepancy of feature dimensions increase considerably. The changes in variables can provide crucial insights for optimizing the acquisition of more high-fidelity data and faster prediction of data in shape design optimization. Thus, it is crucial to analyze the uncertainty of aerodynamic performance at the vehicle system level and provide clear answers.
In recent years, several scholars have conducted uncertainty analyses using methods such as the sensitivity method, moment method, Monte Carlo simulation (MCS), and polynomial chaos expansions (PCEs) [11]. The sensitivity and moment methods are suitable for a system with minuscule parameter uncertainty and an approximately linear model. The MCS method is a statistical method that can present an accurate analysis of the uncertainty of a system only after drawing large numbers of samples. The PCE method has been widely used as it can quantify and characterize uncertainty through polynomial coefficients. However, its drawback is that the chain reaction between the calculated amounts increases considerably, and the dimensions are overinflated, which is known as a “dimensional disaster”. Therefore, developing a reasonable dimensionality reduction method is an urgent problem.
The conventional (wing-tail shape) axisymmetric aerodynamic configuration is widely adopted in air vehicles, especially for supersonic air-to-air and surface-to-air missiles. This configuration can provide greater maneuvering ability while maintaining flight speed. In this study, we proposed an uncertainty-based aerodynamic shape optimization framework for a target drone—with a conventional axisymmetric aerodynamic configuration—to investigate the aerodynamic characteristics of vehicles with such configurations. Analysis and quantification of all uncertainties of the vehicle system were studied based on random theory; moreover, the influences of uncertainties on the aerodynamic performance of conventional axisymmetric aerodynamic vehicles were analyzed. The result of uncertainty-based optimization was compared with deterministic optimization and the initial scheme.
The rest of the paper is organized as follows: In Section 2, the aerodynamic profile parameters of the vehicle and flight conditions are introduced; in Section 3, the precision of the two prediction methods of aerodynamic data is compared with the wind tunnel test data, based on which a high-efficiency and high-precision method is developed, and the potential uncertainty factors are analyzed and classified sequentially; in Section 4, the parallel PCE is described and tested; in Section 5, a mathematical model of deterministic design optimization (DDO) and uncertainty-based design optimization (UDO) are introduced; in Section 6, the results of DDO and UDO are presented, followed by the comparative analysis results.

2. Description

2.1. Shape and Parameters

The vehicle examined in this study was a target drone (Figure 1), which has the typical characteristics of a conventional axisymmetric aerodynamic configuration. It has a Haack series-shaped nose cone, a cylindrical body, four trapezoidal wings, and four trapezoidal tail surfaces. The wing profile (airfoil) is hexagonal, and the tail profile (airfoil) is rhomboidal. Two rectangular cable covers run through both sides of the projectile body. The main wing is installed near the mass center of the vehicle, and the tail is installed after the main wing. There is a convex part on the belly of the vehicle, which results in the fuselage being an imperfect rotating body. The sliders are regarded as part of the fuselage without design optimization.
We used 30 geometric shape parameters and design variables to establish the axisymmetric vehicle model. Numerical variations in the vehicle’s geometric parameters directly affected the aerodynamic characteristics of the vehicle, as they changed the geometric shape of the vehicle. All shape variables of this vehicle are listed and illustrated in Table 1.

2.2. Flight Conditions

The flight trajectory of the vehicle includes the take-off, climb phase, cruise phase, and dive phase (Figure 2). The design altitude of the cruise phase in this ballistic scheme is 20 km, and the design velocity is 2.1 Mach number. As a target drone, maximum flight time should be designed for the aircraft. Under the condition that the propulsion system carries a fixed mass of propellant, the time of the cruise phase can be improved by reducing propellant consumption per unit time, which can be achieved by reducing the drag in the cruise phase.
The main optimization objective of this study was to realize the design of a minimum drag in the cruise phase. There were also other implicit design conditions for this state, including the ability to provide enough lift for retaining flight altitude, flight stability, and controllability. The stability of a vehicle refers to its ability to resist the influence of disturbance, which is characterized by its longitudinal static stability. Controllability is the response characteristic of a projectile body to tail deflection, and in this study, the controllability-to-stability ratio reflects the influence of tail deflection on attitude change. Considering the mutual restrictions between static stability and maneuverability, the constraints of longitudinal static stability and the controllability-to-stability ratio were specified.

3. Aerodynamic Performance Prediction Method and System Uncertainty Factors

3.1. Methods for Estimating Aerodynamic Performance Parameters

From conceptual design to detailed design, aerodynamic characteristics considerably influence the selection of subsequent configuration parameters, such as lift force and drag force. Three prediction methods are used for aerodynamic numerical characteristics, namely, engineering estimation, numerical simulation, and wind tunnel tests. In this study, the aerodynamic performance of the vehicle was examined using wind tunnel tests under two conditions, including free-stream Mach numbers of Ma = 1.75 and 2.53. The comparison of these two sets of aerodynamic performance helps determine the errors in the aerodynamic data.
(1)
Engineering estimation
Engineering estimation is realized using missile data compendium (DATCOM) software, which adopts component combination and modular methods; this enables accurate and efficient prediction of various aerodynamic data of vehicles with traditional aerodynamic configurations [12]. DATCOM also has good adaptability and high estimation accuracy for the conventional axisymmetric configuration, and it can be used to describe various types of airframe surface projections. In DATCOM, the designer defines one vehicle shape by describing its form and the dimensional values of each part (the parameters given in Table 1) and sets the flight conditions such as flight height and Reynold’s number. For a vehicle with a specific shape, all aerodynamic performance data are obtained with an overall calculation time < 2 s, indicating the efficiency and convenience of the engineering estimation method.
(2)
Numerical simulation
Numerical simulations based on CFD are carried out using ANSYS Fluent software [13]. The 3D compressible RANS equations are solved using the shear-stress transport (SST) renormalization group K-ω turbulence model. The SST K-ω model is adopted as it has been widely used in the numerical calculation of multiple flow problems, such as airfoil boundary layer [14], supersonic vehicle design [15], supersonic flow [16], and especially in aerodynamic shape optimization [17]. The second-order upwind scheme is used for the discretization of the convective terms in all transport equations.
Based on the objective of this research, that is, to determine the aerodynamic values of the vehicle and calculate its aerodynamic coefficients, it is important to simplify the calculation model and only retain the features of the main components during parametric modeling in CFD. A fluent mesh was used to generate an unstructured grid in the computing domain. Grid independence was verified before further simulation to exclude the influence of the grid on numerical results (Figure 3). The numerical simulation focuses on the whole turbulent layer; therefore, there is one important relation y + = Δ y ρ u τ / μ , where Δ y is the distance of the first grid point off, and ρ , u τ , μ are the density, friction velocity, and molecular viscosity, respectively. For simulating complicated flowfields, ρ , u τ , μ have different values in different flowfield regions; as a result, Δ y satisfying y + < 1 will also have different conditions [18]. ρ , u τ , μ depend on the standard atmospheric characteristics of the flight state; when y+ is 1, the first layer grid height of the boundary layer can be calculated as Δ y = 1 . 83957 × 10 6   m , which will be used as the reference value for the mesh near the vehicle surface. Grid quality directly affects the solving accuracy of numerical solutions. To ensure better calculation accuracy of models with different grid elements, the height of the first-layer wall grid corresponding to different grids is given in Table 2. At the completion of the numerical simulation, the y+ of the wall surface is counted and displayed as a percentage, in which the proportion of y+ in the interval from 0 to 1 is the focus of attention.
Models with different degrees of grid resolution were established, and the mean pressure coefficients of the vehicle were calculated with different grid numbers. When there were 8.56 × 106 unstructured grids in the computing domain, the numerical results obtained using CFD were highly accurate, and all y+ were in the range of 0–1, and the maximum y+ value was approximately 0.83. As the number of grids increased, all aerodynamic coefficients did not change substantially. The entire computing domain was set as a sphere (Figure 4), and the grid model schematic diagram of the vehicle body is as shown in Figure 5.
(3)
Wind tunnel
Aerodynamic tests of the vehicle were carried out in an FD12 wind tunnel built by the China Academy of Aerospace Aerodynamics. The FD12 wind tunnel test equipment had a cross-section of 1.2 m × 1.2 m and a Mach number range of Ma = 0.4–4. It belongs to the subsonic, transonic, and supersonic wind tunnels, also known as a three-sonic wind tunnel. A 1:6 physics scale model (Figure 6) was used in the wind tunnel tests using air as the test medium. The aerodynamic parameters of the vehicle were measured using wind tunnel tests at an altitude of 10 km under two conditions, i.e., free-stream Mach numbers of Ma = 1.75 and 2.53. For these tests, the flow conditions were set as follows: angle of attack (AOA) α = −2°, 0°, 2°, 4°, 6°, 8°, 11°, 14°, 17°, and 20°, and Reynolds number of Re = 8.82 × 107 and 1.27 × 108. The wind tunnel test adopted a support sting, as shown in Figure 7.
(4)
Comparative analysis
DATCOM, CFD simulations, and wind tunnel tests determined the axial force coefficient (CA), normal force coefficient (CN), and pressure core coefficient (Xp) of the axisymmetric vehicle, as shown in Figure 8 and Figure 9.
All aerodynamic coefficient data are shown in Figure 8 and Figure 9. Aerodynamic coefficient deviations due to the utilization of two computing methods were calculated. For the different Reynolds numbers, the accuracy of CFD was higher than that of DATCOM. In particular, aerodynamic coefficients and their tendency with changes in AOA were close to the wind tunnel data, based on the calculated CA and CN using CFD. The DATCOM projections were appropriate in some (e.g., the center of pressure coefficient), but not in all areas. The empirical database that underpins these forecasts of DATCOM is for a conventional vehicle without a vortical flow. Related studies have shown that Euler simulation cannot predict the flowfield correctly in similar cases, and detached eddy simulations would be a better option [19]. Due to the limitations of the turbulence model, vortex characteristics were not adequately captured. This is attributed to the appearance and disappearance of coherent eddy currents near the fuselage, which disturbs the accuracy of the pitching moment and normal force to some extent. With an increase in AOA, the calculated values of the pressure center coefficient using CFD, and wind tunnel tests show obvious deviations.
The vehicle creates a vortex near the fuselage at a Mach number of Ma = 1.75, such as a wing-tip vortex and a body-shedding vortex. With an increase in AOA, the circulation of the vortex and the vorticity increase considerably. The eddy motion causes energy loss and local pressure on the vehicle surface. This explains the differences between the values of CA predicted by DATCOM using CFD and the wind tunnel tests. This phenomenon also occurs in Ma = 2.53. The appearance of vortical flow is the primary distinction between the three prediction methods; however, this aspect is not considered in DATCOM.
(5)
Modified engineering prediction method
Based on the objective of this study to conduct uncertainty analyses on aerodynamic performance, the computational efficiency of DATCOM matched the research requirements, and the prediction accuracy of the center of pressure coefficient was satisfactory. Therefore, DATCOM was selected to obtain the aerodynamic data of the vehicle. However, there were some errors in the aerodynamic coefficient of DATCOM due to the vortex. In that case, the model for the prediction of aerodynamic data should be corrected based on the DATCOM software. In the subsequent verification of model accuracy, CFD was used as a standard to test the predicted model.
To facilitate optimization, a numerically modified engineering prediction method for force coefficient was developed using the wind tunnel test data. Considering the influence of vortical flow on the aerodynamic coefficient, CA was corrected based on the relative error between the DATCOM data and wind tunnel test data. First, a set of statistical errors Y e can be obtained based on known test data.
Y e ( M a , α ) = C A t e s t ( M a , α ) C A D A T C O M ( M a , α ) C A D A T C O M ( M a , α )
where Ma = 1.75 and 2.53, α = −2°, 0°, 2°, 4°, 6°, 8°, 11°, 14°, 17°, and 20°. Therefore, when Mach number and AOA are satisfied, 1.75 M a 2.53 , and 2 α 20 , respectively, the relative errors of CA  y ˜ C A ( M a , α ) in different states can be obtained by interpolation.
y ˜ C A ( M a , α ) = Y e ( M a , α 1 ) + ( Y e ( M a , α 2 ) Y e ( M a , α 1 ) ) × ( α α 1 ) / ( α 2 α 1 ) Y e ( M a 1 , α ) + ( Y e ( M a 2 , α ) Y e ( M a 1 , α ) ) × ( M a M a 1 ) / ( M a 2 M a 1 )
According to the relative error, the CA predicted by DATCOM can be corrected as:
C A c o D A T C O M = C A D A T C O M × ( 1 + y ˜ C A )
However, CN is corrected using the derivative of CN based on AOA, since the derivative of the CN calculated using wind tunnel data and DATCOM differs considerably with increase in AOA. First, the derivative of CN based on AOA was calculated, as given below:
Δ C N ( M a , α ) = ( d C N t e s t ( M a , α ) d α d C N D A T C O M ( M a , α ) d α ) / d C N D A T C O M ( M a , α ) d α
Similarly, the relative error correction function of CN can be obtained by interpolation through data statistics with wind tunnel data, which is written as y ˜ c n Δ C N ( M a , α ) .
C N c o D A T C O M , α = C N D A T C O M , α 1 × ( 1 + d C N D A T C O M d α × ( 1 + y ˜ c n ) )
The values of CA and CN are calculated based on the modified model, as shown in Figure 10.
After the aerodynamic coefficients obtained using DATCOM were corrected, the aerodynamic coefficient values obtained by three different calculations were compared at Ma = 2.1. Experimental data were not available for this stage since the vehicle in Ma = 2.1 did not undergo wind tunnel testing. For further comparison, CFD was selected as the criteria as it performed well for both Mach numbers, Ma = 1.75 and 2.53. The CA trend predicted by the corrected DATCOM considerably improved, and it was consistent with the trend predicted by CFD. Since the CFD model also has certain errors, some value deviations were allowed. The calculation methods of uncertainty caused by the error of corrected DATCOM are presented in Section 3.2.

3.2. Source of Aerodynamic Deviations

The uncertainty analysis of aerodynamic characteristics should quantitatively evaluate the input and model uncertainties of the modified engineering prediction method. The possible errors of the input parameters can be roughly divided into the following categories: the errors and uncertainties of physical modeling and random aerodynamic interference during flight conditions. The model uncertainties are mainly attributed to the errors in the engineering prediction method.
(1)
Input uncertainties
Compared with the real final product of the vehicle, there are many errors and deviations with the geometric shape parameters in one physical vehicle model due to the processes of manufacturing, installation, and positioning [20]. All these uncertainties are characterized and represented based on design code or engineering experience as shown in Table 1. There are also some random aerodynamic disturbances during flight conditions, such as the disturbance caused by vehicle propulsion system, GPS positioning error, and disturbance due to strong wind weather to vehicle trajectory. Hosder et al. [21] studied the uncertainty analysis of aerodynamic characteristics considering the uncertainty of flight conditions. Loeven et al. [22] conducted uncertainty analysis of subsonic aerodynamic characteristics based on the uncertainty of free-flow velocity. Simon [23] and Chassaing [24] analyzed the uncertainty of the airfoil aerodynamic load distribution considering the uncertainty of Mach number and AOA. Resmini et al. [25] implemented subsonic aerodynamic characteristic analysis of NACA0015 considering the uncertainty of flight conditions and geometric shapes. Considering different parameters as uncertain variables, the results of a robust aerodynamic design differ considerably [26]. These studies show that deviations in flight conditions due to flight disturbances must also be considered; the uncertainty factors are quantitatively characterized, as shown in Table 3. Here, the upper and lower boundaries of the uncertainty range of a flight state are determined, so the uncertainty factors can be represented by the probability theory. They are also regarded as the uniformly distributed random variables.
(2)
Model uncertainties
The errors of the aerodynamic prediction methods also cause uncertainties, which are attributed to cognitive uncertainties. The impact of model uncertainty on data correctness was investigated by Kim et al. [27]. Jens et al. [28] considered the modeling bias and uncertainties in the prediction model. Thus, in this study, model uncertainty was determined according to the relative deviation between the predicted value and the standard value. There are some new uncertainty expression methods, such as probability boxes, Dempster–Shafer theory for representing cognitive uncertainty [29]. Herein, Tchebycheff’s inequality [30] was used to calculate the boundary of random variable variability. The maximum relative error U of independent aerodynamic performance was selected as the model uncertainty, expressed as follows:
U C A = C A D A T C O M C A S T D C A D A T C O M max U C N = C N D A T C O M C N S T D C N D A T C O M max U X c p = X c p D A T C O M X c p S T D X c p D A T C O M max
where the subscripts of “DATCOM” and “STD” indicate that the aerodynamic parameters are calculated using the modified DATCOM and wind tunnel, respectively. The error bar of drag coefficient (Cd) is provided by the maximum value between CA and CN.
U C d = M A X C A ( 1 + U C A ) cos ( a l p h a ) + C N ( 1 + U C N ) sin ( a l p h a ) C A cos ( a l p h a ) + C N sin ( a l p h a )
It is challenging to calculate η using CFD, and η is highly related to Xcp and CN; therefore, its model uncertainty is expressed as a function of U C N and U X c p , as follows:
U η _ u p p e r = ( 1 + U X c p ) ( 1 + U C N ) 1 U η _ l o w e r = = 1 ( 1 U X c p ) ( 1 U C N )
The effect of model uncertainties is utilized to adjust the upper and lower bounds of the output response during the computation of the system response. The upper bound of the output response is established once the mean and standard deviation of the output response have been determined. According to U η _ u p p e r , the numerical difference between the upper limit of the output response and the mean value is recalculated. The distance between the lower end of the range and the mean is changed by U η _ l o w e r . With the same sort of probability distribution, a new mean and standard deviation may be inferred when the upper and lower limits are modified.

4. Uncertainty Analysis Method

Uncertainty analysis is used to calculate the output uncertainties of a system’s performance using the input uncertainties, model uncertainties, and the system’s mathematical model. Herein, the corrected DATCOM was selected as the system’s mathematical model to calculate the aerodynamic performance parameters because it is faster than CFD and more accurate than DATCOM. Both the input uncertainties and model uncertainties of this method are illustrated as probability uncertainties in Section 3.2. There are several commonly used probability uncertainty analysis methods, such as MCS, Taylor series approximation, and PCEs. Compared to MCS and the Taylor series approximation, the advantage of PCE decomposition lies in that it can obtain a system’s output uncertainties with high precision under limited sample points. Therefore, we used PCE in this study.

4.1. Polynomial Chaos Decomposition

According to the Wiener–Askey polynomial chaotic expansion theory [31], in the Askey scheme, each subset of the orthogonal polynomials has a different weighting function in its orthogonality polynomial. The uncertain variables of different distribution types correspond to different PCEs, and they converge in the sense of a norm. For example, Hermite polynomials are associated with the Gaussian distribution, while Legendre polynomials are related to the uniform distribution.
The function Y = F ( X ) can be defined as a physical model under deterministic mapping. In this study, Y = { y 1 , y 2 , , y m } T R m , m 1 is the vector of quantities of interest provided by the model, and X = { x 1 , x 2 , x n } T R n , n 1 is used to describe the random vector of the input variables. To reduce the calculated amount, the random response Y = F ( X ) can be orthogonally expanded in the chaotic polynomial space, and the polynomial can be truncated. The finite-term polynomial whose order is not greater than P is reserved as follows:
Y = F ( X ) F p ( X ) = 0 α p a α ψ α ( X ) + ε
where α = i = 1 n α i , ε refers to the truncation error.
According to the Gameron–Martin theorem [32], the Fourier–Hermite series converges to any L2 Function in the L2 sense. The PC coefficients [33] are estimated using the regression approach expressed as follows:
a ^ = A r g min a R p i = 1 N ( F ( x ( i ) ) α T ψ ( x ( i ) ) ) 2
In this study, variance was used as a standard to evaluate the fitting degree of the PCE model to aerodynamic characteristics. Building an aerodynamic PCE model can be defined as a regression analysis, which is the investigation of the functional relationship between system uncertainty variables and aerodynamic characteristics. In the PCE model, the statistical moments of the response PC expansion can be analytically derived from their coefficients. The mean and variance are expressed as follows:
μ Y = a 0 σ Y 2 = 0 < | α | p a α 2

4.2. Parallel PCEs

When PCE is taken as the uncertainty analysis method, the efficiency of the design optimization based on aerodynamic uncertainty decreases substantially with the increase in the number of design variables [34]. A properly designed program can be much more efficient in parallel computing than in a serial one. Open multi-processing (OpenMP) [35] is a common programming specification in high-performance computing. The shared memory can be accessed by all OpenMP threads, and this programming approach is mainly used in multi-core shared memory scenarios. Therefore, parallel PCE is applied herein based on OpenMP, which is introduced into the derivation and aggregation programming model. When a loop occurs in a program, it can derive multiple threads to form a thread group solution. When the program loop ends, the thread group returns to the main program thread aggregation. The parallel PCE method developed by OpenMP can effectively speed up the uncertainty analysis under each set of design variables, and the efficiency of solving the mean and variance of response values does not decrease when the data of sampling points increase greatly due to system requirements.
Uncertainty analysis of aerodynamic characteristics was carried out using a test sample, which revealed the relationship between the number of computer threads and each sample’s test time, as shown in Figure 11. When the program is run serially, the time of uncertainty analysis spent on a single point is 183.011 s. With the application of parallel execution and an increased number of computer threads, the time decreases gradually. When the number of computer threads is 20, the test time of a sample point is 30 s. The calculation efficiency of the same model is improved by 86.3% with parallel computing. Statistical timing was performed for single CFD evaluation, single DATCOM calculation, PCE model, and parallel PCE model, and the time comparisons are shown in Table 4. All PCE models are calculated with corrected DACTOM here. All the above calculations were performed on the same computer, with an AMD EPYC 7742 64-Core Processor 2.25 GHz CPU.

4.3. Dimension Descending Method Based on Parameter Sensitivity Analysis

Sensitivity analysis [36] is defined as a method of assigning model output changes to different input sources; it is used to study the influence of a single input variable on the model. In this study, the Sobol exponent [37] was used for sensitivity analysis. Before using PCE, this sensitivity analysis was introduced to reduce the dimension of global variables by screening out the variables with a weak influence degree.
Any model can be viewed as a function Y = f ( X ) , where X refers to a vector of d uncertain model inputs X 1 , X 2 , , X d , and Y is defined as an univariate model output. Furthermore, it is assumed that the inputs are independently and uniformly distributed in the unit hypercube, i.e., X i [ 0 , 1 ] i = 1 , 2 , , d . f ( X ) may be decomposed in the following way:
Y = f 0 + i = 1 d f i ( X i ) + i < j d f i j ( X i , X j ) + + f 1 , 2 , , d ( X 1 , X 2 , , X d )
where f 0 refers to a constant, f i stands for a function of X i , and f i j represents a function of X i , X j , etc.
Further, it is assumed that f ( X ) is square-integrable, and variance expression can be described as:
V a r ( V ) = i = 1 d V i + i < j d V i j + + V 12 d
where partial variance is expressed as follows:
V i = V a r X i ( E X ~ i ( Y | X i ) ) V i , j = V a r X i j ( E X ~ i j ( Y | X i , X j ) ) V i V j
where V i represents the influence of the corresponding single-variable on the output, X ~ i indicates the set of all variables except X i , and V i 1 , i 2 , , i s represents the influence of the interactions among input variables on the output.
S X i = V a r ( Y ) V a r X ~ i ( E X i ( Y | X ~ i ) ) V a r ( Y )
The sensitivity index of the Sobol index method includes the main effect index and the total effect index [38]. The total effect index reflects the influence of the main effect of the variable and the cross-effects of other variables on the variance of the function, as shown in Formula (15). According to the vehicle’s aerodynamic model considered herein, the uniform sampling design method was adopted to calculate the aerodynamic performance under the influence of different design variables.

5. Results and Discussion

5.1. Design Optimization

5.1.1. Mathematical Model of Optimized Design

The design state is supersonic flight (Ma = 2; altitude = 20 km; Re = 2.7 × 107), and the uncertainties of system parameters, such as Ma, altitude, and shape parameters, such as CHORD and SPAN, were considered. All variables are shown in Table 1 and Table 2. Since the total length and the cross-section diameters of the fuselage considerably influence the propulsion system and structural system, which cannot be changed, nine related variables were excluded from the design optimization and set as constants. Thus, only 21 design variables were involved in the optimization process, as shown in Table 1, and their value ranges are shown in Table 5.
The optimization objective at trim AOA is to lower the drag’s mean and variance and to turn the multi-objective optimization issue into a single-objective optimization problem by weighting parameter k for optimization computation. As for the constraint settings, the static stability margin was 2.5–5%, and the controllability-to-stability ratio range was 0.5 η 1.5 , where η refers to the ratio of elevator trim angle to trim AOA. In UDO, all constraint conditions are transformed to the possibility of design value to meet requirements, i.e., >0.95. The deterministic optimization is described by Formula (16).
min F [ C d m e a n ( x 1 , x 2 , , x n ) , C d s i g m a ( x 1 , x 2 , , x n ) ] α = α trim s . t . 0.025 X c p ( x 1 , x 2 , , x n ) 0.05 0.5 η ( x 1 , x 2 , , x n ) 1.5 x L x x U
The uncertainty-based optimization model is described using Formula (17):
min F [ C d m e a n ( x 1 , x 2 , , x n ) , C d s i g m a ( x 1 , x 2 , , x n ) ] α = α trim s . t . P ( 0.025 X c p ( x 1 , x 2 , , x n ) 0.05 ) 0.95 P ( 0.5 η ( x 1 , x 2 , , x n ) 1.5 ) 0.95 x L x x U
In this study, corrected DATCOM was used to evaluate the aerodynamic performance of the vehicle. During the calculation in each round of optimization, different vehicle shapes can be described with each set of special design variable values. The influence of uncertain factors on a system input can also be represented by numerical values, and they are all processed into random numbers, allowing the design values to vary equally within a small range. Simple calculations were used to determine the trim of AOA and its corresponding aerodynamic performance. Under the condition of trim AOA, this vehicle generates more lift than its own gravity, as depicted in Figure 12.

5.1.2. Optimization Algorithm

Reasonable and feasible optimization algorithms should be selected for the UDO framework; reliable convergence should be ensured, and the efficiency of the convergence should be improved. The multi-island genetic algorithm (MIGA) [39], as a parallel genetic algorithm, is based on the traditional genetic algorithm through the addition of sectionalization. In this case, the individual population was divided into several subgroups, and the subgroups migrated to different islands, with each island executing the genetic algorithm. Due to the screening of heredity and variation, the algorithm can be used to determine the location of the optimal solution.
To prevent large computation costs in uncertain optimization, the number of stages of external circulation should be minimized. The parameters related to the MIGA algorithm are briefly tested to ensure that the optimization algorithm achieves global convergence in both deterministic and uncertain optimization under the same number of iterations, as well as the optimization information and global search ability in multiple populations. In the optimization process in this study, the group size, generation, and the number of islands were set as 10, 100, and 10, respectively. Its iteration curve is shown in Figure 13. The MIGA method successfully considers both the convergence speed and the approximation accuracy of the ideal solution under this parameter setting; even when the number of cycles increases, the approximation accuracy of the optimal solution no longer improves substantially. Consequently, the parameter selection mentioned above is thought to be appropriate for the optimization problem stated in this study.

5.1.3. Sensitivity Analysis of Uncertain Parameters

In this optimization problem, there are 35 uncertain parameters, indicating a considerably low computation efficiency when PCE is used. It is necessary to scale the design space of uncertain parameters for the improvement of the efficiency of PCE used in uncertainty analysis. This strategy can considerably reduce calculation costs [40]. Therefore, the influence of all design variables on aerodynamic performance, which is calculated by the corrected DATCOM, is sorted by the sensitivity analysis. Therefore, the influence of all design variables on aerodynamic performance is sorted by the sensitivity analysis. The feasible solutions of DDO are used as samples in the sensitivity analysis. The main effect indexes of CA, CN, Xcp, and η are shown in Figure 14, which were obtained with the Sobol sensitivity analysis of all uncertainty parameters.
The total effect index shows how aerodynamic performance is influenced by the cross-effects between a single variable and other variables. The variables with a larger total effect index indicate that their linear effect and the interaction effect between them and other variables are larger than others, that is, their change has a greater influence on the change in aerodynamic performance. Specific values of the sensitivity analyses are shown in Figure 14. For the sensitivity analysis of CN, the total effect for AOA was the largest, followed by flight velocity Mach and nose diameter DNOSE. The semi-span length of the wing and the thickness-to-chord ratio of the wing root also substantially influence the value of CN. During the actual flight, pressure drag is the main reason for CA value variation, and the impact of shape parameters, such as the leading-edge radius and the thickness-to-chord ratio of the pressure on the surface of the vehicle. In the meantime, flight velocity determines the dynamic pressure; therefore, it is closely related to CA and CN. The static stability margin and operating stability ratio represent the stability and maneuverability of the aircraft, respectively, and their values are closely related to the position of the vehicle’s mass center. In this case, they were considerably influenced by Xc. The uncertainties from different sources and with a considerable effect on the vehicle’s aerodynamic performance should be considered in uncertainty characterization and design optimization. Therefore, the aforementioned six uncertain parameters are selected in the UDO.

5.1.4. Uncertainty Analysis Method

At this point, the PCE is established by sampling in the uncertain space, which can characterize the system’s performance by the mean and variance. In PCE, the number of terms M raises its own polynomials with is both the total degree P and the number of input variables N.
M = ( N + P P )
In the meantime, the number of input variables N and the total degree in the polynomial P are adjusted, and the fitting quality of PCE to the system uncertainty is judged according to R 2 and R a d j 2 [33]. In this study, we selected input variables N = 6 and total degree P = 5. According to the sequencing results of the above sensitivity analysis, the six input variables with the greatest impact on each aerodynamic performance were selected to form the respective PCEs for CA, CN, Xcp, and η . Moreover, the test of PCE efficiency and accuracy uses the feasible solutions of DDO as the sample.
MCS is considered a well-established tool for the determination of the impact of parameter variations [41]. It is also considered suitable for investigating uncertainty in complex systems, which is used as the standard in this study for comparison with PCE. As shown in Figure 15, the convergence of the Cd mean is derived using PCE, and a comparison with MCS with an increasing number of sampling points. N = 500 is sufficient for the PCE model to achieve high computing accuracy. MCS is used to calculate the dependability index based on N = 10 5 number of model evaluations, which is used as the reference solution to the problem considered in this study.

5.2. Comparison and Validation of Different Results

All optimization results are presented in Table 6. The comparison between the deterministic optimization results and the original design configuration is shown in Figure 16. In the deterministic optimization scheme, the position of the vehicle’s main wing moves forward, the span length increases, the sweep angle of the tail increases considerably, and the ratio of the tail root increases slightly; these findings are different from those of the original design scheme. In the uncertain optimization result scheme, the position of the main wing moves forward, the three variables of the root chord tip and length of the chord increase, and the position of the tail moves forward; these findings are different from those of the original design scheme.
According to the calculation results shown in Table 7, the error between the aerodynamic coefficient of optimization results and CFD simulation results is small, indicating that the improved DATCOM method has high calculation accuracy. Compared with MCS, the error in PCE results is small, indicating that PCE is also characterized by high precision.
Compared with the original scheme, both optimization results reduce the drag coefficient in horizontal flight by decreasing AOA, indicating that the lift-to-drag ratio of the aircraft in the benchmark scheme is not optimal, and that the lift-to-drag ratio and drag reduction can be achieved by increasing the wingspan to a certain extent through optimization.
The standard deviation of Cd of the uncertain optimization results is smaller than that of the deterministic scheme, and the probability of constraint satisfaction is greater than 99%. The uncertainty optimization improved the P s t a t i c = 87.2% to 99.8% of the deterministic scheme, showing that uncertainty-based optimization can effectively improve the robustness and reliability of the scheme.
The main difference between the two configurations and the original scheme lies in trim AOA, and the flowfield in the head cone region of the CFD calculation results has changed, as shown in Figure 17. The deterministic optimization scheme shows large areas with low pressure at the tail of the main wing and waist. This phenomenon does not exist in uncertainty-based optimization, which may mainly account for the fact that the controllability-to-stability ratio does not satisfy the constraint.

6. Conclusions

In this study, a system-level uncertainty analysis was conducted, and an optimization method based on the parallel PCE method was developed. The aerodynamic performance of a target drone in a typical conventional axisymmetric aerodynamic configuration was optimized and analyzed. The main conclusions were drawn as follows:
(1)
The corrected DATCOM showed high accuracy in aerodynamic prediction of conventional axisymmetric aerodynamic configuration. Compared with traditional MCS, PCE also has the characteristics of high precision based on the selection of partial variables, showing the effectiveness of the dimensionality reduction method.
(2)
The parallel PCE uncertainty analysis method can be used to effectively characterize the uncertainty of the system. Compared with the general PCE, the cost time is shortened from 183.011 s to 29.287 s, thereby, effectively reducing the cost of uncertainty analysis.
(3)
With the same optimizer, uncertainty-based optimization can reduce a vehicle’s drag coefficient in the cruise phase, and it can effectively improve the robustness and reliability of the system’s performance.
(4)
The results of sensitivity analysis showed that the uncertainties of AOA and vehicle mass center have considerable influence on the aerodynamic performance of conventional axisymmetric vehicles.

Author Contributions

Conceptualization, H.Z. and D.X.; methodology, X.P. and H.Z.; software, X.P. and M.X.; validation, X.P., H.Z. and M.X.; formal analysis, W.W.; investigation, H.Z.; resources, G.C.; data curation, D.X.; writing—original draft preparation, X.P.; writing—review and editing, D.X. and H.Z.; visualization, X.P. and M.X.; supervision, W.W.; project administration, G.C. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by Hao Zhu grant number 20180151001. And the APC was funded by Aeronautical Science Foundation of China.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Burner, A.; Liu, T.; DeLoach, R. Uncertainty of Videogrammetric Techniques used for Aerodynamic Testing. In Proceedings of the 22nd AIAA Aerodynamic Measurement Technology and Ground Testing Conference, St. Louis, MI, USA, 24–26 June 2002. [Google Scholar] [CrossRef] [Green Version]
  2. Salehi, S.; Raisee, M.; Cervantes, M.J.; Nourbakhsh, A. An efficient multifidelity ℓ1-minimization method for sparse polynomial chaos. Comput. Methods Appl. Mech. Eng. 2018, 334, 183–207. [Google Scholar] [CrossRef]
  3. Yu, X.; Li, M.; An, G.; Liu, B. A Coupled Effect Model of Two-Position Local Geometric Deviations on Subsonic Blade Aerodynamic Performance. Appl. Sci. 2020, 10, 8976. [Google Scholar] [CrossRef]
  4. Shi, W.; Chen, P.; Li, X.; Jing, R.; Jiang, H. Uncertainty quantification of the effects of squealer tip geometry deviation on aerothermal performance. Proc. Inst. Mech. Eng. Part A J. Power Energy 2020, 234, 1026–1038. [Google Scholar] [CrossRef]
  5. Fusi, F.; Congedo, P.M.; Guardone, A.; Quaranta, G. Shape optimization under uncertainty of morphing airfoils. Acta Mech. 2017, 229, 1229–1250. [Google Scholar] [CrossRef]
  6. Ryan, K.M. Robust Multi-Objective Optimization of Hypersonic Vehicles under Asymmetric Roughness-Induced Boundary-Layer Transition. Ph.D. Dissertation, University of Maryland, College Park, MD, USA, 2014. [Google Scholar]
  7. Pelletier, D.; Turgeon, E.; Lacasse, D.; Borggaard, J. Ccedil Adaptivity, Sensitivity, and Uncertainty: Toward Standards of Good Practice in Computational Fluid Dynamics. AIAA J. 2003, 41, 1925–1933. [Google Scholar] [CrossRef]
  8. Luckring, J.; Hemsch, M.; Morrison, J. Uncertainty in computational aerodynamics. In 41st Aerospace Sciences Meeting and Exhibit; AIAA: Reno, NV, USA, 2003; p. 409. [Google Scholar] [CrossRef] [Green Version]
  9. Zhang, W.; Wang, S.; Wang, Q.; Ju, S.; Yan, C. Reverse jet parameters study on aerodynamic thermal uncertainty of a blunt body. Aerosp. Sci. Technol. 2020, 107, 106260. [Google Scholar] [CrossRef]
  10. Da Ronch, A.; Drofelnik, J.; van Rooij, M.P.; Kok, J.C.; Panzeri, M.; Voß, A. Aerodynamic and aeroelastic uncertainty quantification of NATO STO AVT-251 unmanned combat aerial vehicle. Aerosp. Sci. Technol. 2019, 91, 627–639. [Google Scholar] [CrossRef]
  11. Walters, R.W.; Huyse, L. Uncertainty Analysis for Fluid Mechanics with Applications. NASA/CR-2001-210648. Available online: https://www.cs.odu.edu/~mln/ltrs-pdfs/icase-2002-1.pdf (accessed on 1 July 2022).
  12. Williams, J.E.; Vukelich, S.R. The USAF Stability and Control Digital DATCOM; Users Manual [R]; Usaf Stability & Control Digital Datcom: St. Louis, MN, USA, 1979; Volume I. [Google Scholar]
  13. Rybdylova, O.; Al Qubeissi, M.; Braun, M.; Crua, C.; Manin, J.; Pickett, L.M.; De Sercey, G.; Sazhina, E.M.; Sazhin, S.S.; Heikal, M. A model for droplet heating and its implementation into ANSYS Fluent. Int. Commun. Heat Mass Transf. 2016, 76, 265–270. [Google Scholar] [CrossRef] [Green Version]
  14. Kavvadias, I.S.; Papoutsis-Kiachagias, E.M. The continuous adjoint approach to the k–ω SST turbulence model with applications in shape optimization. Eng. Optim. 2015, 47, 1523–1542. [Google Scholar] [CrossRef]
  15. Bramantya, M.A.; Nararya, J.P. Effects caused by the implementation of the area rule theory to the fuselage of a surveillance unmanned aerial vehicle using computational fluid dynamics methods. In AIP Conference Proceedings; AIP Publishing LLC: Yogyakarta, Indonesia, 2020. [Google Scholar] [CrossRef]
  16. Georgiadis, N.; Yoder, D.; Towne, C.; Engblom VBhagwandin, W.; Lankford GPower, D.; Nelson, C. Wind-US Code Physical Modeling Improvements to Complement Hypersonic Testing and Evaluation. In Proceedings of the 47th AIAA Aerospace Sciences Meeting including The New Horizons Forum and Aerospace Exposition, Orlando, FL, USA, 5–8 January 2009; p. 193. [Google Scholar]
  17. Şumnu, A.; Güzelbey, İ.H.; Öğücü, O. Aerodynamic Shape Optimization of a Missile Using a Multiobjective Genetic Algorithm. Int. J. Aerosp. Eng. 2020, 2020, 1528435. [Google Scholar] [CrossRef]
  18. Gao, Z.; Jiang, C.; Lee, C. Improvement and application of wall function boundary condition for high-speed compressible flows. Sci. China Technol. Sci. 2013, 56, 2501–2515. [Google Scholar] [CrossRef]
  19. Ghoreyshi, M.; Badcock, K.J.; Da Ronch, A.; Marques, S.; Swift, A.; Ames, N. Framework for Establishing Limits of Tabular Aerodynamic Models for Flight Dynamics Analysis. J. Aircr. 2011, 48, 42–55. [Google Scholar] [CrossRef]
  20. Ju, Y.; Liu, Y.; Jiang, W.; Zhang, C. Aerodynamic analysis and design optimization of a centrifugal compressor impeller considering realistic manufacturing uncertainties. Aerosp. Sci. Technol. 2021, 115, 106787. [Google Scholar] [CrossRef]
  21. Hosder, S.; Walters, R.W.; Balch, M. Point-Collocation Nonintrusive Polynomial Chaos Method for Stochastic Computational Fluid Dynamics. AIAA J. 2010, 48, 2721–2730. [Google Scholar] [CrossRef]
  22. Loeven, G.J.A.; Witteveen, J.A.S.; Bijl, H. Probabilistic collocation: An efficient non-intrusive approach for arbitrarily distributed parametric uncertainties. In 45th AIAA Aerospace Sciences Meeting and Exhibit; AIAA: Reno, NV, USA, 2007; Volume 317. [Google Scholar] [CrossRef]
  23. Simon, F.; Guillen, P.; Sagaut, P.; Lucor, D. A gPC-based approach to uncertain transonic aerodynamic. Comput. Methods Appl. Mech. Eng. 2010, 199, 1091–1099. [Google Scholar] [CrossRef]
  24. Chassaing, J.C.; Lucor, D. Stochastic investigation of flows about airfoils at transonic speeds. AIAA J. 2010, 48, 938–950. [Google Scholar] [CrossRef]
  25. Resmini, A.; Peter, J.; Lucor, D. Sparse grids-based stochastic approximations with applications to aerodynamics sensitivity analysis. Int. J. Numer. Methods Eng. 2015, 106, 32–57. [Google Scholar] [CrossRef]
  26. Du, X.; Leifsson, L. Optimum aerodynamic shape design under uncertainty by utility theory and metamodeling. Aerosp. Sci. Technol. 2019, 95, 105464. [Google Scholar] [CrossRef]
  27. Kim, J.-T.; Stubbs, N. Model-Uncertainty Impact and Damage-Detection Accuracy in Plate Girder. J. Struct. Eng. 1995, 121, 1409–1417. [Google Scholar] [CrossRef]
  28. Refsgaard, J.C.; van der Sluijs, J.P.; Brown, J.; van der Keur, P. A framework for dealing with uncertainty due to model structure error. Adv. Water Resour. 2006, 29, 1586–1597. [Google Scholar] [CrossRef] [Green Version]
  29. Zhang, J.; Yin, J.; Wang, R. Basic Framework and Main Methods of Uncertainty Quantification. Math. Probl. Eng. 2020, 2020, 6068203. [Google Scholar] [CrossRef]
  30. Oberguggenberger, M.; Fellin, W. Reliability bounds through random sets: Non-parametric methods and geotechnical applications. Comput. Struct. 2008, 86, 1093–1101. [Google Scholar] [CrossRef]
  31. Xiu, D.; Karniadakis, G.E. The Wiener-Askey Polynomial Chaos for Stochastic Differential Equations. SIAM J. Sci. Comput. 2003, 24, 619–644. [Google Scholar] [CrossRef]
  32. Cameron, R.H.; Martin, W.T. The orthognoal development of nonlinear functionals in series of Fourier-Herimte functionals. Ann. Math. 1947, 48, 386–392. [Google Scholar] [CrossRef]
  33. Blatman, G.; Sudret, B. An adaptive algorithm to build up sparse polynomial chaos expansions for stochastic finite element analysis. Probabilistic Eng. Mech. 2010, 25, 183–197. [Google Scholar] [CrossRef]
  34. Zhang, W.; Wang, Q.; Zeng, F.; Yan, C. An adaptive Sequential Enhanced PCE approach and its application in aerodynamic uncertainty quantification. Aerosp. Sci. Technol. 2021, 117, 106911. [Google Scholar] [CrossRef]
  35. Chapman, B.; Jost, G.; Van Der Pas, R. Using OpenMP: Portable Shared Memory Parallel Programming (Scientific and Engineering Computation); MIT Press: Cambridge, MA, USA; London, UK, 2007. [Google Scholar]
  36. Saltelli, A.; Annoni, P. Sensitivity Analysis; Wiley. 2009. Available online: https://www.wiley.com/en-us/Sensitivity+Analysis-p-9780470743829 (accessed on 19 May 2022).
  37. Sobol, I.M. Sensitivity estimates for nonlinear mathematical models. Math. Model. Comput. Exp. 1993, 1, 407–414. [Google Scholar]
  38. Homma, T.; Saltelli, A. Importance measures in global sensitivity analysis of nonlinear models. Reliab. Eng. Syst. Saf. 1996, 52, 1–17. [Google Scholar] [CrossRef]
  39. Alba, E.; Luna, F.; Nebro, A.J.; Troya, J.M. A survey of parallel distributed genetic algorithms. Complexity 1999, 4, 31–52. [Google Scholar] [CrossRef]
  40. Dai, J.; Liu, P.; Qu, Q.; Li, L.; Niu, T. Aerodynamic optimization of high-lift devices using a 2D-to-3D optimization method based on deep reinforcement learning and transfer learning—ScienceDirect. Aerosp. Sci. Technol. 2022, 121, 107348. [Google Scholar] [CrossRef]
  41. Argentini, T.; Pagani, A.; Rocchi, D.; Zasso, A. Monte Carlo analysis of total damping and flutter speed of a long span bridge: Effects of structural and aerodynamic uncertainties. J. Wind Eng. Ind. Aerodyn. 2014, 128, 90–104. [Google Scholar] [CrossRef]
Figure 1. Axisymmetric normal vehicle.
Figure 1. Axisymmetric normal vehicle.
Aerospace 09 00396 g001
Figure 2. Project ballistic trajectory of a vehicle.
Figure 2. Project ballistic trajectory of a vehicle.
Aerospace 09 00396 g002
Figure 3. Grid independence verification of the axisymmetric vehicle. (a) Axial force coefficient, (b) normal force coefficient, (c) pitching moment coefficient.
Figure 3. Grid independence verification of the axisymmetric vehicle. (a) Axial force coefficient, (b) normal force coefficient, (c) pitching moment coefficient.
Aerospace 09 00396 g003
Figure 4. CFD unstructured mesh of the axisymmetric vehicle.
Figure 4. CFD unstructured mesh of the axisymmetric vehicle.
Aerospace 09 00396 g004
Figure 5. CFD unstructured mesh of the vehicle surface.
Figure 5. CFD unstructured mesh of the vehicle surface.
Aerospace 09 00396 g005
Figure 6. Model photographs.
Figure 6. Model photographs.
Aerospace 09 00396 g006
Figure 7. Support sting in wind tunnel test.
Figure 7. Support sting in wind tunnel test.
Aerospace 09 00396 g007
Figure 8. Aerodynamic coefficient at Ma= 1.75. (a) CA, (b) CN, and (c) Xp.
Figure 8. Aerodynamic coefficient at Ma= 1.75. (a) CA, (b) CN, and (c) Xp.
Aerospace 09 00396 g008
Figure 9. Aerodynamic coefficient at Ma = 2.53. (a) CA, (b) CN, (c) Xp.
Figure 9. Aerodynamic coefficient at Ma = 2.53. (a) CA, (b) CN, (c) Xp.
Aerospace 09 00396 g009
Figure 10. Aerodynamic coefficient at Ma = 2.1. (a) CA and (b) CN.
Figure 10. Aerodynamic coefficient at Ma = 2.1. (a) CA and (b) CN.
Aerospace 09 00396 g010
Figure 11. Different number of threads and their corresponding elapsed time.
Figure 11. Different number of threads and their corresponding elapsed time.
Aerospace 09 00396 g011
Figure 12. Uncertainty-based optimization process with parallel PCE and corrected DATCOM.
Figure 12. Uncertainty-based optimization process with parallel PCE and corrected DATCOM.
Aerospace 09 00396 g012
Figure 13. Change curve of the objective function with the number of iterations.
Figure 13. Change curve of the objective function with the number of iterations.
Aerospace 09 00396 g013
Figure 14. Total effect index of each variable.
Figure 14. Total effect index of each variable.
Aerospace 09 00396 g014aAerospace 09 00396 g014b
Figure 15. Relationship between sampling points and Cd mean.
Figure 15. Relationship between sampling points and Cd mean.
Aerospace 09 00396 g015
Figure 16. Three aerodynamic shape schemes of vehicle.
Figure 16. Three aerodynamic shape schemes of vehicle.
Aerospace 09 00396 g016
Figure 17. Comparison between flowfield cloud images and the three schemes considered herein. (a) Initial configuration design (α = 10°), (b) deterministic optimization (α = 8°), (c) and uncertainty-based optimization (α = 8°).
Figure 17. Comparison between flowfield cloud images and the three schemes considered herein. (a) Initial configuration design (α = 10°), (b) deterministic optimization (α = 8°), (c) and uncertainty-based optimization (α = 8°).
Aerospace 09 00396 g017aAerospace 09 00396 g017b
Table 1. Parametric modeling of vehicle.
Table 1. Parametric modeling of vehicle.
Design VariableUnitPhysical SignificancePerturbation VariableRelative DeviationOptimization VariableProbability Distribution
LNOSEmNose lengthx1±0.001nouniform
LCENTRmFuselage lengthx2±0.001nouniform
DNOSEmNose diameter at basex3±0.001nouniform
DCENTRmFuselage diameterx4±0.001nouniform
BNOSEmNose bluntness radius or radius of truncationx5±0.001nouniform
Xc---Mass center of vehiclex6±0.01nouniform
XLE1mDistance from missile nose to chord leading edge at wing root chord stationx7±0.01yesuniform
SSPAN11mSemi-span length of wingx8±0.001yesuniform
CHORD11mPanel chord length at wing root chord locationx9±0.001yesuniform
CHORD12mPanel chord length at wing tip chord locationx10±0.001yesuniform
LER11mLeading edge radius at wing root chord stationx11±0.001yesuniform
LER12mLeading edge radius at wing tip chord stationx12±0.001yesuniform
ZUPPER11---Thickness-to-chord ratio of upper surface at wing root chord stationx13±0.001yesuniform
ZUPPER12---Thickness-to-chord ratio of upper surface at wing tip chord stationx14±0.001yesuniform
LMAXU11---Fraction of chord from leading edge to maximum thickness of upper surface at root chord stationx15±0.001yesuniform
LMAXU12---Fraction of chord from leading edge to maximum thickness of upper surface at tip chord stationx16±0.001yesuniform
LFLATU11---Fraction of chord of constant thickness section of upper surface at root chord stationx17±0.001yesuniform
LFLATU12---Fraction of chord of constant thickness section of upper surface at tip chord stationx18±0.001yesuniform
XLE21mDistance from missile nose to chord leading edge at tail root chord stationx19±0.01yesuniform
XLE22mDistance from missile nose to chord leading edge at tail root chord stationx20±0.01yesuniform
SSPAN21mSemi-span length of tailx21±0.001yesuniform
CHORD21mPanel chord length at tail root chord locationx22±0.001yesuniform
CHORD22mPanel chord length at tail tip chord locationx23±0.001yesuniform
LER21mLeading edge radius at tail root chord stationx24±0.001yesuniform
LER22mLeading edge radius at tail tip chord stationx25±0.001yesuniform
ZUPPER21---Thickness-to-chord ratio of upper surface at tail root chord stationx26±0.002yesuniform
ZUPPER22---Thickness-to-chord ratio of upper surface at tail tip chord stationx27±0.002yesuniform
THETA1degsetting angle of wingsx28±1nouniform
THETA2degsetting angle of tailx29±1nouniform
Xhinge---Hinge axis positionx30±0.01nouniform
Table 2. Aerodynamic coefficients under different grid numbers.
Table 2. Aerodynamic coefficients under different grid numbers.
Aerodynamic Coefficients
Number of Elements Δ y CACNCM0 < y+ ≤ 1
1,583,7950.00030.61633.41052.095792.13%
2,334,7770.00020.60373.41542.110499.16%
3,981,4900.00020.60193.42372.117199.73%
4,947,0500.00020.59683.41292.103699.82%
6,350,8020.000080.58763.42092.109499.90%
7,123,4090.000080.57933.42152.116999.96%
8,562,3530.000080.57663.42122.1173100%
9,690,5070.000080.57673.42132.1174100%
Table 3. Range and category of flight condition uncertain variables.
Table 3. Range and category of flight condition uncertain variables.
Design VariableUnitPhysical SignificancePerturbation VariableDeviation Range
MACH---Mach numberx31±0.02
ALTMAltitudex32±50
BETAdegAngle of Sideslipx33±1
DELTAdegAngle of tail reflectionx34±1
ALPHAdegAngle of attackx35±0.5
Table 4. Statistical time under different aerodynamic performance evaluation methods.
Table 4. Statistical time under different aerodynamic performance evaluation methods.
CFD EvaluationDATCOM EvaluationPCE ModelParallel PCE Model
Statistical time86,605.43 s0.352 s183.011 s29.474 s
Table 5. Optimization range of design variables.
Table 5. Optimization range of design variables.
Design VariablesRangeDesign VariablesRange
XLE1[3.2, 3.6]XLE22[5.57, 5.8]
LER11[0.0001, 0.001]LER21[0.001, 0.008]
LER12[0.001, 0.008]LER22[0.001, 0.008]
CHORD11[1.0, 1.5]CHORD21[0.2, 0.4]
CHORD12[0.1, 0.4]CHORD22[0.05, 0.2]
SSPAN1[0.4, 0.6]SSPAN2[0.2, 0.4]
LMAXU11[0.1, 0.3]ZUPPER11[0.01, 0.04]
LMAXU12[0.1, 0.3]ZUPPER12[0.01, 0.04]
LFLATU11[0.5, 0.6]ZUPPER21[0.05, 0.1]
LFLATU12[0.5, 0.6]ZUPPER22[0.05, 0.1]
XLE21[5.5, 5.7]
Table 6. Deterministic optimization and uncertainty-based optimization results of vehicle.
Table 6. Deterministic optimization and uncertainty-based optimization results of vehicle.
VariablesRangeInitialDeterministic OptimizationUncertainty-Based Optimization
XLE1[3.2, 3.7]3.373.3113.302
SSPAN1[0.4, 0.7]0.450.5430.569
CHORD11[1.0, 1.5]1.21.1091.076
CHORD12[0.1, 0.4]0.20.2620.235
LER11[0.001, 0.008]0.0010.0010.001
LER12[0.001, 0.008]0.0010.0010.001
ZUPPER11[0.01, 0.04]0.020.010.011
ZUPPER12[0.01, 0.04]0.0250.010.011
XLE21[5.5, 5.7]5.6795.5075.514
XLE22[5.57, 5.8]5.7545.7455.787
SSPAN2[0.2, 0.4]0.250.2530.217
CHORD21[0.2, 0.4]0.250.2030.213
CHORD22[0.05, 0.2]0.10.0540.063
LER21[0.001, 0.008]0.0010.0010.001
LER22[0.001, 0.008]0.0010.0010.001
ZUPPER21[0.05, 0.1]0.0640.0510.052
ZUPPER22[0.05, 0.1]0.070.0510.054
LMAXU11[0.1, 0.3]0.20.2610.117
LMAXU12[0.1, 0.3]0.20.2050.146
LFLATU11[0.5, 0.6]0.60.5180.534
LFLATU12[0.5, 0.6]0.60.5010.537
Table 7. PCE and MCS of different optimization.
Table 7. PCE and MCS of different optimization.
Uncertainties Calculated by PCEUncertainties Recalculated by MCS
Initial SchemeDDO ResultUDO ResultInitial SchemeDDO ResultUDO Result
α t r i m 10881088
C d C F D 1.3251.18211.18891.3251.18211.1889
C d m e a n 1.320261.19051.19431.33411.17031.1716
C d s i g m a 0.05540.06250.06130.04530.07050.0621
s t a t i c m e a n 0.04510.03960.03390.04440.03980.0345
s t a t i c s i g m a 0.00150.00170.00170.00140.00190.0016
η m e a n 0.69410.51930.51810.70430.51490.4981
η s i g m a 0.02480.02420.03380.02110.02710.0319
p s t a t i c 0.9820.9760.9990.9920.9990.999
p η 0.9880.7240.9990.9960.8720.998
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Peng, X.; Zhu, H.; Xu, D.; Xiao, M.; Wang, W.; Cai, G. Aerodynamic Performance Uncertainty Analysis and Optimization of a Conventional Axisymmetric Vehicle Based on Parallel Polynomial Chaos Expansions. Aerospace 2022, 9, 396. https://doi.org/10.3390/aerospace9080396

AMA Style

Peng X, Zhu H, Xu D, Xiao M, Wang W, Cai G. Aerodynamic Performance Uncertainty Analysis and Optimization of a Conventional Axisymmetric Vehicle Based on Parallel Polynomial Chaos Expansions. Aerospace. 2022; 9(8):396. https://doi.org/10.3390/aerospace9080396

Chicago/Turabian Style

Peng, Xun, Hao Zhu, Dajun Xu, Mingyang Xiao, Weizong Wang, and Guobiao Cai. 2022. "Aerodynamic Performance Uncertainty Analysis and Optimization of a Conventional Axisymmetric Vehicle Based on Parallel Polynomial Chaos Expansions" Aerospace 9, no. 8: 396. https://doi.org/10.3390/aerospace9080396

APA Style

Peng, X., Zhu, H., Xu, D., Xiao, M., Wang, W., & Cai, G. (2022). Aerodynamic Performance Uncertainty Analysis and Optimization of a Conventional Axisymmetric Vehicle Based on Parallel Polynomial Chaos Expansions. Aerospace, 9(8), 396. https://doi.org/10.3390/aerospace9080396

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