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
, where
is the distance of the first grid point off, and
are the density, friction velocity, and molecular viscosity, respectively. For simulating complicated flowfields,
have different values in different flowfield regions; as a result,
satisfying
will also have different conditions [
18].
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
, 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 × 10
6 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 × 10
7 and 1.27 × 10
8. 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
can be obtained based on known test data.
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,
, and
, respectively, the relative errors of
CA in different states can be obtained by interpolation.
According to the relative error, the
CA predicted by
DATCOM can be corrected as:
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:
Similarly, the relative error correction function of
CN can be obtained by interpolation through data statistics with wind tunnel data, which is written as
.
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:
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.
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
and
, as follows:
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 , 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 . With the same sort of probability distribution, a new mean and standard deviation may be inferred when the upper and lower limits are modified.