Abstract

Flow field reconfiguration and drag reduction have been investigated by numerically solving Navier-Stokes equations, with the 2nd order dispersion controlled dissipative scheme adopted for the convective term discretion and the 2nd order central difference scheme for viscous term. Spikes and sideward jets nearby the spike noses were set up in the front of the blunt body to reconfigure the flow field and reduce drag during hypersonic flight. The spikes and the sideward jets led to the formation of the conical shock, the reattached shock, and the circumfluence, instead of a bow shock. According to the quantitative analysis, the reconfiguration of the flow field played a dominant role in drag reduction. Basing on the locations of the reattached shock and its interaction with the conical shock, a combined parameter was proposed for a primary criterion to design the spiked bodies with sideward jets for the purpose of drag reduction.

1. Introduction

The concept of a spiked blunt body was first proposed by Bogdonoff [1]. However, investigations of flow separation in front of blunt bodies at supersonic speeds have been made since the early 1950s [25]. Research has shown that spikes can construct a separation zone over blunt bodies to significantly lower the aerodynamic heat rate and pressure distribution, which was valuable for thermal protection and drag reduction. Such a simple and efficient configuration has made it appealing. In the 1960s, much research was conducted on separation characteristics and corresponding flow instability. Maull investigated unsteady flow field properties caused by spike length and shape on blunt bodies, concluding that flow oscillation was due to two factors, separation induced by shock wave and flow reattachment [6]. Wood investigated spiked cone cylinders flying at Mach 10 and defined five different flow patterns and corresponding scopes using a combination of spike lengths and cone angles [7]. Hahn introduced mass injection at different spiked blunt body positions flying at Mach 3.3, including the tip, base, and the reattachment region [8]. Reding et al. considered unsteady aerodynamics for a spiked drag reduction device based on structure deflection coupled with thermal expansion caused by aerodynamic heating [9]. Calarese and Hankey experimentally and numerically studied shockwave oscillation modes on a spiked body at Mach 3, observing both symmetrical and asymmetrical oscillations that demonstrated and confirmed the existence of standing rotational waves and their corresponding effects [10]. A study by Hutt and Howe first considered the preliminary optimization design of vehicle configuration, in which the effects of spikes on bodies with different cross sections were compared by measuring force coefficients and the pressure center in a supersonic wind tunnel [11].

Alongside the rapid expansion of computer technology beginning in the 1990s, many studies on numerical simulations and solutions have allowed for a remarkable level of understanding with regards to spiked blunt bodies. Cooperation between experiments and numerical simulationss validated many of these studies. Mehta studied the aerodynamic heat flux and pressure distribution over spiked bodies at Mach 6.8 and the relationship between them [12]. Milićev et al. studied the influence of spike shape on aerodynamic coefficients at supersonic flow past blunt-nosed bodies, including drag, lift, pitching moment, and pressure center [13]. Panaras and Drikakis investigated spiked blunt bodies in supersonic and hypersonic flow and showed that unsteady flow was characterized by periodic radial inflation and conical separation bubble collapse [14]. They also pointed out that there was a certain spike length, beyond which flow is stable. Menezes et al. studied spiked blunt body drag in a hypersonic shock tunnel and showed that drag reduction significantly decreased as attack angle increased. In some critical cases, the spiked blunt body changed into a resistance element [15, 16]. Morgenstern showed that an asymmetric flow field in a three-dimensional simulation resulted in unsteady fluctuation during flight with attack angle, leading to axial pulsation and normal forces, which resulted in an extreme disadvantage for flight and controlling [17]. Thurman [18] and Schulein [19] achieved some positive results by changing aerodynamic characteristics utilizing a rotating spike, which required a special system to control and rotate the rod and was too complicated for practical application. Liu and Jiang gave some experimental investigations on pressure distribution and flow field structure of the spiked blunt body with unexpanded sideward jets [20]. Ma et al. furtherly investigated the flow unsteadiness in the circumfluence formed by the spiked bodies [21]. Eghlima and Mansour developed a method for drag reduction by the combination of a spike and a counter jet on blunt bodies [22]. Bibi et al. numerically investigated the use of the active opposing jet concept in combination with geometric variations of the opposing jet nozzle to alleviate high wave drag formation [23]. Huang et al. numerically studied the influence of the cavity configuration on the drag and heat flux reduction mechanism of a blunt body [24]. Han set up the new nonablation concept utilizing adaptive drag reduction and thermal protection system (NADTPS) for hypersonic vehicles basing on a spiked blunt body with sideward jet nearby the spike nose [25]. The NADTPS uses aerodynamic heating on the blunt body and pressure changes coolant into vapor, thus generating a sideward jet and convective cooling at the spike nose, preventing the spike from being destroyed. This sideward jet can lower shock wave interaction on the blunt body shoulder, thus optimizing aerodynamic characteristics and efficiently demonstrating experiments in a hypersonic tunnel with an attack angle.

This paper investigated the spiked blunt bodies with and without sideward jets flying at hypersonic speed, using numerical simulations. The incoming flow Mach number changed from 5 to 10, and the spike length ratio over the blunt body diameter changed from 0.5 to 1.0. The flow structure reconfiguration, especially the shock wave structure, was analyzed by flow field visualization. Aerodynamic force characteristics were calculated and compared with those of blunt bodies to analyze drag reduction effects and the corresponding mechanisms.

2. Governing Equations and Numerical Methods

The two-dimensional axisymmetric compressible Navier-Stokes equation set was applied as the governing equations, with the thermally perfect gas model and ideal gas state equation [26]. The flow medium can be considered as a single component with diffusion neglected. Under the Jacobian transformation between physical space and computational space, the governing equation is and each vector component transformation can be written as follows: where ξ and η stand for computational space, and x and y stand for physical space. J stands for the Jacobian determinant:

The detailed formula for each vector in (2) can be written as follows:

The viscosity coefficient μ was calculated with Sutherland formula, and the heat conduction coefficient k can be calculated via the Prandtl number definition derived as Pr = k/μcp, where cp denotes the specific heat at constant pressure and Pr = 0.72.

Jiang et al. analyzed the dispersion terms in the modified equation and theoretically derived the dispersion conditions for nonoscillatory shock capturing schemes, on which a series of dispersion controlled dissipative (DCD) schemes were set up [27, 28]. The 2nd order DCD scheme was adopted to discrete the convective terms, and the Steger-Warming method [29] was applied for flux vector splitting. The 2nd order central difference scheme [26] was used for viscous term discretization. The explicit method was used for time integration, with the CFL condition for time-step calculation and scheme stability requirements. The local time step method [30] was introduced to accelerate the simulation convergence process using different time steps at grid points.

The configuration and dimension of the spiked blunt body are shown in Figure 1(a). The diameter D of the spherical head cylinder was chosen as the characteristic length, and the diameter of the spike was set to d/D = 0.1. The L/D denotes spike length and was set from 0.5 to 1.0 with an interval of 0.1. Rt stands for transition part radius at the spike root and was set as Rt/D = 0.2. The sideward jet position was set at the range from (0, r) to (d/10, r) nearby the spike nose in the Cartesian coordinate system. The physical model and computational domain are shown in Figure 1(b). The RL method [31] was applied for grid generation by solving an elliptic partial differential equation set. The spiked blunt body model and the corresponding grid system are shown in Figure 1(b), in which orthogonality is satisfied in each grid system and refined by algebraic method near the wall boundary.

Four different grid resolutions have been applied for primary investigation, including 101 × 401, 101 × 801, 201 × 801, and 401 × 1601. Much finer structure can be captured by the resolution of both 201 × 801 and 401 × 1601, especially the shock waves, spline line, and flow structure interactions. And the results basing on the resolution of 201 × 801 can match the result of the finer resolution, which can get to a convergent solution much faster. Therefore, the computational grids were set as 201 grid points in a vertical direction to the body surface and 801 in the direction along the surface. The boundary condition was set as a constant temperature wall with Tw = 300.0 K. The sideward jets were set to be sonic jets, where the static temperature and density were equal to the local flow condition, and the speed equal to the corresponding sound speed. And the sideward jet was set as a slit with the width of d/10 near by the spike nose from (0, r) to (d/10, r) in the Cartesian coordinate system in Figure 1.

The numerical solution validation was carried out by comparing the numerical wave structure and experimental flow field as shown in Figure 2. The conditions were set as M = 6.0 and L/D = 1.0. Numerical and experimental results are in agreement, indicating the main flow field structures, including conical shock, reattached shock, slip line, shear layer, and separation zone. Aerodynamic force validation calculations utilized flow simulation around a sphere flying at hypersonic speed with the diameter of 0.1 m as listed in Table 1. The drag coefficient scope was from 0.8713 to 0.8910 in the simulations, close to the 0.87 to 0.90 scope identified in the experiments [32].

3. Results and Discussion

3.1. Typical Flow Field Structure

Former investigations on spiked blunt bodies usually focused on the formation and stability of the separation zone. However, this paper considers that drag reduction mainly depends on wave drag reduction due to flow field reconfiguration, in which the bow shock is changed into a conical shock interacting with a reattached shock. In this section, the shock wave structures around the spiked blunt bodies are illustrated using 6 different Mach numbers and 6 different spike lengths, including those with spiked blunt bodies and sideward sonic jets near the spike nose, as shown in Figures 314.

The typical shock wave structures of pure spiked blunt bodies with defined labels are shown in Figure 15. Density gradient in the y direction was used to describe a Schlieren photograph and reconfigured flow field. The spike nose became the new stagnation with a little bow shock in front and conical shock around. High-speed flow after the conical shock traveled to the blunt body shoulder, leading to reattachment. Hence, the reattached shock was formed on the blunt body shoulder, resulting in a zone with high pressure. The inverse pressure gradient caused some of the gas flow to the spike nose, in the opposite direction of incoming flow. Thus, a circumfluence was generated around the spike and the blunt body nose. In the circumfluence, the gas flow velocity was lower than that after the conical shock, generating a shear layer. The shear flow reduced some of the fluid in the circumfluence, while the reattachment fed back some fluid so that the circumfluence remained in a dynamic equilibrium state. If such an equilibrium state was broken, it causes the circumfluence and shock wave structures to vibrate. The relationship between the conical shock and reattached shock might be the real mechanism for the vibration or collapse of the circumfluence, while the former investigations took the circumfluence vibration as the unsteadiness mechanism. Due to the circumfluence, a separation shock will be generated near the spike nose before the shear layer. The conical shock interacts with the reattached shock on the blunt body shoulder, forming a slip line to match the transition between flow across two shock waves, which is exactly a shear layer due to viscosity.

The reattached shock and its interaction with conical shock were pushed away from the central line when spike length increased or the Mach number was lowered, as shown in Figures 38. The conical shock angle increased when the sideward sonic jet was applied near the spike nose, which pushed the interaction between conical shock and reattached shock away from the central line, as shown in the 36 figures in Figures 914. The sonic sideward jet affected the flow structure locally by keeping the whole macrostructure similar to that around the pure spiked blunt bodies. The detailed flow structure of the spiked blunt body with sonic sideward jet is shown in Figure 16. The sideward jet is set near the spike nose with equal velocity to local sound speed. The numerical Schlieren photograph shows that the sideward jet flows in the shear layer exchanging with the flow after conical shock and circumfluence, rather than directly flowing into the circumfluence. Due to the sideward jet, the boundary layer near the spike nose was forced to rise and flow direction was forced to swerve, putting the separation shock position forward and resulting in secondary vortex structure formation in the wake of the separation shock. The relative motion between the flow after separation shock and the sideward jet was transformed into one part of the shear layer due to the forward motion from the separation shock, causing the shear layer starting point to be drawn forward. The separation shock angle was increased by the sideward jet compared to the pure spike body flow field. Therefore, the angle of conical shock increased and the interaction between conical shock and reattachment shock was pushed away from the axial line. This can significantly affect the shockwave interaction on the shoulders, which is considered as a typical element of flow field reconfiguration.

3.2. Shock Wave Interaction Characteristics

Spikes reconfigured the flow field around the blunt bodies, and the interaction between shock waves was the most typical feature during the reconfiguration. The spikes set in front of the blunt bodies transformed the bow shock into a conical shock. Such a transformation must be considered as an optimization of aerodynamic characteristics, despite shock wave interactions generated on the blunt body shoulder. The main flow field configurations were similar for the pure spiked blunt bodies and spiked blunt bodies with sideward jets, including the separation shock, the conical shock, the reattached shock, the shear layer, and the circumfluence. The geometrical position of the shear layer was located near the inner tangent line connecting the spike nose and the blunt body shoulder. The flow after the conical shock was supersonic or even hypersonic, which experienced a deflection on the blunt body shoulder and generated the reattached shock. The conical shock also interacted with the reattached shock, results in the formation of a slip line, which is exactly a shear layer due to viscosity.

The main flow field characteristics can be distinguished, including the shock waves, the corresponding interactions, and the circumfluence. Two nondimensionalized parameters have been defined for the main structure description in the flow field, denoted as y1/D and y2/D, as shown in Figure 17(a). The first parameter y1/D stands for the reattached shock position, which is equal to the distance between the reattachment position and the central line. The second parameter y2/D stands for the position of the interaction between conical shock and reattached shock, which is also equal to the distance to the central line. Furthermore, combined with the spike lengths, y1/D can also denote the shape or volume of the circumfluence generated by the flow field reconfiguration.

In this paper, 36 combinations of the incoming flow Mach numbers and the spike lengths for pure spiked bodies were taken into account so did the 36 combinations of the spiked bodies with sonic sideward jets. The data of (y1/D, y2/D) in the 72 cases were tallied and plotted, as shown in Figure 17(b), with y2/D for the horizontal coordinate axis and y1/D for the vertical coordinate axis. The scope of the 36 points for spiked bodies with sonic sideward jets was generally higher than that for pure spiked bodies, in both y1/D and y2/D direction of the coordinate system. It demonstrates that the sonic sideward jets led both the reattached shock and the conical shock interactions to be pushed away from the central line, affecting the aerodynamic force of the entire system. A third order polynomial can be used to fit the 72 discrete points, which generally describes a position’s variation tendency of the reattached shock and its interaction. Compared with the reattached shock, the interaction between reattached shock and conical shock was more sensitive to the conditions, such as the spike lengths, the income flow Mach numbers, or the sideward jets.

The tendency of the 3rd order polynomial cure shows that there might exist a limited value of y1/D, which maintained the reattached shock position as the spike length increased. Based on the numerical results, the limit value was about 0.43, located at the angle about 60 degrees to the central line. Different from the reattached shock, the interaction between the conical shock and the reattached shock moved away from the central line as the spike length increased so does the slip line resulted by the interaction.

3.3. Aerodynamic Drag

The spikes reconfigure the flow field around the blunt body, changing the bow shock into three main aspects: the conical shock, the reattached shock, and the circumfluence. The flow field reconfiguration greatly reduced the shock wave drag by eliminating the strong bow shock. The reattached shock and its interaction with the conical shock will be pushed further away from the central line as the spike length increased, significantly lowering the strength and changing the aerodynamic characteristics. The drag coefficients of the pure spiked blunt bodies with 6 different spike lengths L/D flying at 6 different Mach numbers are shown in Figure 18(a), where L/D denotes the horizontal coordinate axis. The drag coefficient Cd decreased as L/D increased from 0.5 to 1.0, with the drag reduction of 30–40%. Both the conical shock and the reattached shock are pushed to the central line as Mach number increases, which leads to an increase in drag. However, the existence of the circumfluence will prevent the conical shock and the reattached shock from being close to the central line and forming a balance. If such a balance was broken, the shocks and the circumfluence become unstable and start vibrating, during which the circumfluence seems to collapse [7]. Hence, the effect on drag by Mach number will be much less than that by the spike length in these 36 ordinary cases.

Because the wave drag plays a dominant role, the drag coefficient for a sphere flying at hypersonic speed is close to that of a sphere cylinder, which is 0.87 to 0.90 experimentally [32] and 0.87 to 0.89 in numerical results. Comparing with the pure blunt sphere cylinders, the spikes reduce wave drag as much as 35–65%, which has been confirmed by other studies [1517]. Bushnell pointed out that the effective load can be increased by 5–10% if the aerodynamic drag got discounted by as much as 1%, which is likely why so many studies have been conducted regarding spiked blunt bodies [33].

The drag coefficients of the spiked blunt bodies flying at hypersonic speed with sonic sideward jets near the spike nose are shown in Figure 18(b). These include the spike lengths and Mach numbers corresponding to the 36 cases of the pure spiked blunt bodies. The effect on Cd caused by the spike length L/D is similar to the cases without jets. The sonic sideward jet will push the conical shock and the reattached shock away from the central line, thus lowering the reattached shock strength and its interaction with conical shock. Therefore, the sonic sideward jet lowers spiked body drag more than the pure spiked bodies. Comparing with the 36 pure spiked blunt body cases, the sideward jets decreased the drag at a range of about 5–20%. It was especially true for cases with short spikes. Thus, drag reduction caused by the spikes and the sonic sideward jets will be as much as 45–70% when compared with the blunt body drag. The effect on drag coefficients by the Mach number due to the sonic sideward jets is much lesser than those without the jets. Though Mach number varies, the data points are concentrated at the same spike length. Thus, drag reduction might increase based on the jet variety and strength, which requires further researches.

3.4. Reattached Shock and Circumfluence Criterion

The main characteristics of reconfigured flow fields have been analyzed. The reconfiguration using spikes and sonic sideward jets leads to the formation of the conical shock, the reattached shock, and the circumfluence. The roles played by each factor and the relative importance of these factors have been discussed in former sections. The aerodynamic force characteristics associated with these factors have also been analyzed. To simplify flow field description, two nondimensional geometry parameters were defined as y1/D and y2/D, standing for the distance from the reattached shock and its interaction with the conical shock to the central line, respectively. Edney investigated the interaction between oblique shocks and bow shocks and classified six different types of flow fields according to interaction position and flow structure characteristics [34]. Though these shock interactions in this paper are different, flow characteristic classification can still be taken as a reference for analysis.

The interactions between the conical shock and the reattached shock in all cases, including the pure spiked bodies and the spiked bodies with sideward jets, were similar to type V or type VI interaction determined by different combinations of the spike length and the Mach number, as shown in Figures 314. The supersonic or hypersonic flow after the conical shock reattaches to the blunt body shoulder and leads to the formation of the reattached shock. The circumfluence comes from this flow reattachment due to the pressure gradient generated. The circumfluence feeds back after the conical shock, affects the flow reattachment, and lowers the strength. Therefore, a dynamical balance or couple between the reattached shock and the circumfluence was created. And both the reattached shock and the circumfluence were so stable that there was no shock interaction as strong as type IV existing in the 72 cases. However, no single parameter, including Ma, Re, L/D, y1/D, or y2/D could characterize this balance. Hence, the reattached shock and the circumfluence are preferred to be focused, based on former analysis. These two factors might be responsible for flow field reconfiguration. The balance can be quantitatively described as the distance from the reattached shock to the central line and the circumfluence volume. Therefore, a new parameter can be set up by combining these two factors, denoted as y12L/D3.

The drag coefficients in the 72 cases of the pure spiked bodies and the spiked bodies with sonic sideward jets were plotted together, as shown in Figure 19. These results varied based on different horizontal coordinates. Generally, the drag coefficients decreased as the horizontal coordinate increased, remaining nearly constant when the three parameters exceeded some critical value. However, Figures 19(a) and 19(b) indicate obvious disparity among data points, suggesting that neither y1/D nor y2/D can describe the flow field reconfiguration mechanism by the spikes and the sideward jets. With y12L/D3 taken as the horizontal coordinate, the distribution of data points was significantly changed and more concentrated in the two drawn lines, as shown in Figure 19(c). Thus, the parameter y12L/D3 might be a simple and suitable description for the two main factors during flow field reconfiguration. In addition, drag coefficient tendency changes and becomes more gently when y12L/D3 is greater than 0.14. The reattached shock and the circumfluence criterion for drag reduction could be used by determining the critical value of y12L/D3. Therefore, y12L/D3 ≈ 0.14 can be considered as a criterion for the design optimization regarding the spike length and the strength of sideward jets.

4. Conclusions

The 72 flow field reconfiguration cases were numerically investigated in this paper, including 36 pure spiked blunt bodies and 36 spiked bodies with sideward jets. Main flow characteristic structures and corresponding effects on aerodynamic drag have been analyzed. The spikes generated a conical shock in front of the blunt body, and the flow after conical shock reattached to the shoulder to form a reattached shock. The pressure gradient caused by the reattachment resulted in the circumfluence formation, which reversely effected the conical shock and the reattachment. Thus, a dynamical balance can be set up among the conical shock, the reattached shock, and the circumfluence, playing a dominant role in aerodynamic drag reduction. The variation tendency depends on parameter combination and might be determined in utilizing complex mechanisms, despite the force variation tendency being drawn over a single variable, such as the spike length, the position of reattached shock, and the shock/shock interaction. A criterion was primarily proposed basing on main flow structure geometric parameters y12L/D3, explicitly standing for the effect both of the reattached shock and the circumfluence. Since the connection between y1/D and y2/D can be set up utilizing a third order polynomial, the combined parameter contains the effect and strength of the conical shock and the shock interaction on the shoulder. The drag coefficient variation can be divided into two stages, and the drag reduction might be neglected when y12L/D3 is greater than 0.14 for both the pure spiked blunt bodies and the spiked bodies with sonic sideward jets.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

Acknowledgments

This study was supported by the National Natural Science Foundation of China (Grant nos. 90916028, 11472281, and 11532014) and the National Key Research and Development Program (Grant no. 2016YFA0401201).