Abstract
A finite difference method is used to do direct numerical simulation (DNS) of hypersonic unsteady flowfield under the action of freestream pulse wave. The response of the hypersonic flowfield to freestream pulse wave is studied, and the generation and evolution characteristics of the boundary layer disturbance waves are discussed. The effects of the pulse wave types on the disturbance mode in the boundary layer are investigated. Results show that the freestream disturbance waves significantly change the shock standoff distance, the distribution of flowfield parameters and the thermodynamic state of boundary layer. In the nose area, the main disturbance modes in the boundary layer are distributed near the fundamental mode. With the evolution of disturbance along with streamwise, the main disturbance modes are transformed from the dominant state of the fundamental mode to the collective leadership state of the second order and the third order harmonic frequency. The intensity of bow shock has significant effects on both the fundamental mode and the harmonic modes in each order. The strong shear structure of boundary layer under different types of freestream pulse waves reveals different stability characteristics. The effects of different types of freestream pulse waves are significant on the distribution and evolution of disturbance modes. The narrowing of frequency band and the decreasing of main disturbance mode clusters exist in the boundary layer both for fast acoustic wave, slow acoustic wave and entropy wave.
Keywords
- hypersonic flow
- boundary layer
- receptivity
- pulse wave
- numerical simulation
1. Introduction
With the progress that human achieved in the exploration of aviation and aerospace fields, the hypersonic vehicle technology has developed quickly. The design of hypersonic vehicle is a complex problem based on hypersonic aerodynamics and involving multidisciplinary and multi-domain. Problems of hypersonic aerodynamics and aerodynamic heating, structures and materials technology are the main technological difficulties. There remains much variance between hypersonic flow and supersonic flow about the problem of aerodynamics. For the hypersonic boundary layer, entropy layer and shock wave layer overlapping each other, while it is high temperature and low density flow in boundary layer, which undoubtedly makes the hypersonic boundary layer flow more complex. Therefore, hypersonic aerodynamics is a key technology in the research and development of hypersonic vehicle. Various flow disturbance problems, usually exists in the flying environment of vehicle. For instance, the explosive blast wave [1], reverse jet [2, 3], the non-uniformity flow, the instability of flight, the rough wall and so on, which would happen in flying. It can be seen disturbances are common in flow fields. It’s significantly different about the ideal hypothesis state of the steady flowfield and the flow condition existing in disturbance waves. The disturbance in the flowfield will have a significant influence on aerothermodynamics characteristics. The disturbance, whether it is strong or weak, will interfere with the flowfield, especially the shock and boundary layers. After the disturbance in the flowfield interference with shock wave and boundary layer, the disturbance wave will be induced. Then the induced wave will cause further interaction with boundary layer, and create new unstable waves. The stability characteristic of boundary layer and the laminar-turbulent transition mechanism will be significantly changed due to the induced unstable waves. Laminar-turbulent transition not only affects the aerodynamic heating of the wall of hypersonic vehicle, but directly changes the aerodynamic force. Especially, laminar-turbulent transition will greatly increases aerodynamic drag, which reduces the lift drag ratio of hypersonic vehicle significantly and increases the requirement of thermal protection.
Therefore, it is necessary to accurately predict the hypersonic unstable flowfield and the flowfield response characteristic induced disturbance waves during the design and development process of hypersonic aircraft. Considering the complexity and expensiveness of the hypersonic vehicle wind tunnel test, it is of practical significance to conduct the numerical simulation of hypersonic unsteady flowfield. In recent years, the hypersonic flowfield response induced by different disturbance waves and the influences of the disturbance wave on the stability of the boundary layer are studied by many scholars using numerical or experimental methods. Ma and Zhong [4] investigated the response of hypersonic boundary layer over a blunt cone to freestream acoustic waves at Mach 7.99. The receptivity of a flat plate boundary layer to a freestream axial vortex is discussed by Boiko [5]. Zhong [6] investigated the leading-edge receptivity to freestream disturbance waves for hypersonic flow over a parabola. The effect of wall disturbances on hypersonic flowfield, and the response of hypersonic boundary layer to wall disturbances are also widely studied [7, 8, 9, 10]. Literature [11] points out that, after the interaction between any form of freestream disturbance and the shock wave as well as the boundary layer in hypersonic flow field, three independent forms of disturbance waves, including acoustic disturbance (fast and slow acoustic disturbance), entropy wave disturbance, and vortex wave disturbance will be generated. Among these investigations, most scholars are committed to study the effects of freestream continuous disturbance or the effects of wall disturbance on the stability of boundary layer and laminar-turbulent transition [12, 13, 14]. There is still less research on the effects of freestream pulse disturbance on the stability of boundary layer, and the mechanism in this field is still not fully understood. So, this paper aims to study hypersonic flowfield and the stability characteristics of boundary layer under the freestream pulse wave. It should be mentioned that the significantly different showed in the influence mechanism of the wave types of continuous disturbance on the stability of boundary layer that had been pointed out by some scholars [6, 15, 16]. Ma and Zhong [15, 16] investigated the response of hypersonic plate boundary layer under different types of disturbance waves using direct numerical simulation and linear stability theory. It is found that acoustic disturbance has the greatest influence on the stability of boundary layer, and it is pointed out that the action mechanism of freestream entropy wave disturbance and vortex wave disturbance is similar to that of the fast acoustic wave disturbance, but different from the slow wave disturbance. Although it is significantly different between pulse disturbance wave and continuous disturbance wave in form, there is reason to believe that the response characteristics of hypersonic boundary layer are also significantly different under different freestream pulse waves. Therefore, the receptivity characteristics of hypersonic boundary layer under freestream pulse waves with different types are discussed in this paper.
In this investigation, the numerical simulations of hypersonic flowfield over a blunt wedge under the action of freestream pulse wave were conducted using a high order finite difference scheme. The response property of hypersonic flowfield under the action of freestream pulse wave analyzed, and the hypersonic boundary layer stability is investigated. The receptivity characteristic of hypersonic boundary layer under the action of different freestream pulse waves was compared, and the effects of the pulse wave types on the generation and evolution of the disturbance mode in the hypersonic boundary layer are discussed.
2. Basic equations and numerical method
2.1. Navier-Stokes equation
According to the forms of conservation equation, momentum equation and energy equation, the three basic equations of fluid governing equations can be written as a general form, that is, the two-dimensional unsteady compressive N-S equation can be expressed as:
Where the state vector
Similarly,
where
2.2. Numerical method
In this paper, a high order finite difference method is used to solve the flowfield governing equation directly. The inviscid vector flux of the Navier-Stokes equation is divided into the positive and negative convection terms using S-W method [17]. The positive and negative convection terms are discretized by the 5th order upwind WENO scheme [18]. The viscous term is discretized by the 6th order central difference scheme [19]. In order to obtain the transient information of the flowfield and reduce computation time, the 3th TVD Runge-Kutta scheme is used for time advance [20].
The spatial discretization for positive and negative convection terms can be expressed as Eqs. (7) and (8), respectively:
where
The spatial discretizations for viscous terms can be expressed as follow:
where
The TVD Runge-Kutta discretization can be expressed as follows
where Δ
3. Computing model and conditions
The model parameters and calculation conditions of flowfield calculation includes the freestream condition, model parameter, boundary conditions and meshing. For freestream conditions, the freestream temperature is equal to 69 K, and the Mach number is equal to 6. The Reynolds number, based on the nose radius, is equal to 10,000. The angle of attack is equal to 0°. Figure 1 shows the computing model and schematic diagram. Calculation mode is a blunt wedge with the wedge angle of 16°e; the nose radius
In order to explore the influence of freestream pulse wave on hypersonic flow field, the interaction process between freestream pulse slow acoustic wave and hypersonic flowfield is direct numerical simulated. In present investigation, the stable flow over a blunt wedge at Mach 6 is calculated firstly, and then the simulation of hypersonic unsteady flowfield under the freestream pulse wave is conducted by introducing freestream pulse slow acoustic wave at the upstream boundary of computing domain.
The introducing time of the freestream pulse wave is recorded as
Where
4. Results and discussion
4.1. Response of the hypersonic flowfield to freestream pulse wave
4.1.1. Response of the hypersonic flowfield
Before analyzing, it should be pointed out that the numerical simulation method and the grid independence in present work is validated in our previous investigation [21, 22, 23], which indicates the numerical method adopted in this paper is reliable. Figure 5 shows the contours of pressure under freestream pulse wave at different times.
Figure 3(a)–(d) in the figure correspond to
Figure 4 shows the contours of density under freestream pulse wave at different times. Figure 4(a)–(f) correspond to
Figure 6 shows the contours of velocity along axis
4.1.2. Response of the boundary layer
Figure 7 shows the distribution of friction factor disturbance on wall under freestream pulse wave at different times (
Figure 8 shows the comparisons of the density, temperature and pressure disturbances on wall at different times. Figure 8 indicates that, under freestream pulse wave, the variations of the density, temperature and pressure disturbances on wall are similar to that of friction factor. Zhang et al. [19] studied the evolution of continuous small disturbance in hypersonic flow using direct numerical simulation (DNS), and found that, because of the normal shock wave, the forcing disturbance in freestream is enlarged. From Figure 8, it is obtained that the amplitude of parameters disturbance on the wall is enlarged sharply relative to that of initial wave in freestream. The numerical results agree with Li’s results. It also shows that the disturbance amplitudes of pressure on the nose of blunt wedge are larger than that on non-nose, and the former are several times, even more than 10 times, larger than the latter; while the discrepancy is tiny for temperature and density.
To study the evolution of the disturbance wave in boundary layer during the interaction process, the temperature disturbance on wall along the streamwise is considered. Figure 9 shows the distribution of temperature disturbance on wall along the streamwise at
To analyze the disturbance mode distribution in the hypersonic boundary layer and its evolution characteristics along the streamwise under the action of freestream pulse wave, the fast Fourier transformation method is used to decode the time domain signal of the boundary layer disturbance. Furthermore, the signals are expanded in frequency order to make it as a frequency function, thereby converting the time domain signal of the pressure disturbance into a frequency domain signal. Figure 10 shows the results of the Fourier frequency spectral analysis of pressure disturbance in boundary layer under freestream pulse wave at different locations. Figure 10(a)–(d) are the results of the Fourier frequency spectral analysis when
Figure 11 shows the evolution of different disturbance modes along streamwise in boundary layer. Figure 11(a)–(f) are the evolution of the fundamental frequency mode, the second order harmonic frequency, the third order harmonic frequency, the fourth order harmonic frequency, the fifth order harmonic frequency and the sixth order harmonic frequency, respectively. From Figure 11(a), in the nose region, the amplitude of the fundamental frequency becomes smaller. When out of the nose region, the amplitude of the fundamental frequency gradually becomes smaller and generally remains. In Figure 11(b), it can be seen that in the nose region, the amplitude of the second order harmonic frequency mode significantly becomes smaller. Out of the nose region, the amplitude of the second order harmonic frequency mode continues to increase. From Figure 11(c), it can be seen that the amplitude of the third order harmonic frequency mode becomes smaller in the nose region, and the third order harmonic frequency mode first increases and then decreases in the downstream evolution. From Figure 11(d)–(f), it shows that the evolution of the fourth order harmonic frequency mode, the fifth order harmonic frequency mode and the sixth order harmonic frequency mode along the streamwise is similar to that of the third order harmonic frequency mode. As shown in the figure, whether it is the fundamental mode or harmonic modes, in the area of the nose boundary layer, the amplitude of the disturbance significantly becomes smaller, indicating that the intensity of the bow shock in the hypersonic flowfield becomes weak, and it has significant effects on both the fundamental mode and the harmonic modes in each order. According to the figure, it also shows that except for the fundamental mode, from the nose area to the non-nose area, the amplitude of each harmonic frequency mode all increase significantly. It is believed that this phenomenon is caused by the recompression of the flowfield, which means in this region, the effect of the shock wave intensity on the harmonic frequency modes of the boundary layer is smaller than the recompression effect, and the effect of the shock wave intensity on the fundamental mode of the boundary layer is larger than the recompression effect.
4.2. Effects of pulse wave type
In this section, to discuss the effects of the pulse wave types on the generation and evolution of the disturbance mode in the hypersonic boundary layer, the numerical simulations of hypersonic flowfield under the action of freestream pulse fast acoustic wave, slow acoustic wave and entropy wave were conducted. It should be noted that the wall temperature
Figure 13 shows the pressure disturbance in boundary layer under different freestream pulse waves. Figure 13(a) and (b) are corresponding the position
The Fourier frequency spectral analysis of pressure disturbance in boundary layer under different freestream pulse waves is shown in Figure 14. Figure 14(a)–(c) are corresponding the position
5. Conclusions
In present study, a finite difference DNS method is used to investigate the response of the hypersonic flow field to freestream pulse wave, and the generation and evolution characteristics of the disturbance mode in the hypersonic boundary layer are analyzed. The effects of the pulse wave types on the generation and evolution of the disturbance mode in the hypersonic boundary layer are investigated. The study drew the conclusions as follow:
The freestream disturbance waves significantly interact with the bow-shaped shock waves, which greatly changes the shock standoff distance and the distribution of flowfield parameters in the active region. The stronger the intensity of the shock wave is, the more intense the effect of the freestream pulse slow acoustic wave and the bow-shaped shock wave are. The freestream pulse wave’s influence on the thermodynamic state in boundary layer is significantly greater than that on the thermodynamic state the outside boundary layer. There are many obvious disturbance characteristic regions in the hypersonic boundary layer under the action of freestream pulse disturbances. There is complex interference in the boundary layer under the action of the pulse wave.
In the nose area, the main disturbance modes in the boundary layer are distributed near the fundamental frequency, and the disturbance mode component of other frequencies is relatively small. In the nose area, the disturbance mode component ratio near the fundamental frequency decreases rapidly, the high frequency disturbance harmonic mode like the second order harmonic frequency and above significantly increases. With the evolution of the disturbance waves in the boundary layer from the upstream to the downstream, the main disturbance modes in the boundary layer are transformed from the dominant state of the fundamental frequency to the collective leadership state of the second order harmonic frequency and the third order harmonic frequency. The intensity of the bow shock in the hypersonic flowfield has significant effects on both the fundamental mode and the harmonic modes in each order. From the nose area to the non-nose area, the effect of the shock wave intensity on the harmonic frequency modes in the boundary layer is less than the recompression effect, and the effect of the shock wave intensity on the fundamental mode in the boundary layer is larger than the recompression effect.
The interaction between the strong shear structure of boundary layer and different types of freestream pulse waves have different mechanisms of action. The strong shear structure of boundary layer under the action of different types of freestream pulse waves reveals different stability characteristics. The effects of different types of freestream pulse waves are significant on the distribution of modes in boundary layer and the evolution of disturbance modes along with the streamwise. The narrowing of frequency band and the decreasing of main disturbance mode clusters exist in the boundary layer both for freestream pulse fast acoustic wave, slow acoustic wave and entropy wave.
References
- 1.
Pierce HB, Manning JC. Experimental Investigation of Blast Loading on an Airfoil in Mach Number of 0.7 Airflow with Initial Angle-of-Attack Change of 20. NASA Resport, 1963, NASA-TN-D-2876. - 2.
Aso S, Hayashi K, Mizoguchi M. A study on aerodynamic heating reduction due to opposing jet in hypersonic flow. AIAA Paper. 2002; 2002 :646 - 3.
Hayashi K, Aso S, Tani Y. Numerical study on aerodynamic heating reduction by opposing jet. Memoirs of the Faculty of Engineering. 2006; 66 (1):39-54 - 4.
Ma YB, Zhong XL. Boundary-layer receptivity of Mach 7.99 Flow over a blunt cone to free-stream acoustic waves. Journal of Fluid Mechanics. 2006; 556 :55-103 - 5.
Boiko AV. Receptivity of a flat plate boundary layer to a free stream axial vortex. European Journal of Mechanics - B/Fluids. 2002; 21 :325-340 - 6.
Zhong X. Leading-edge receptivity to free stream disturbance waves for hypersonic flow over a parabola. Journal of Fluid Mechanics. 2001; 441 :315-367 - 7.
Fedorov AV, Khokhlov AP. Receptivity of hypersonic boundary layer to wall disturbances. Theoretical and Computational Fluid Dynamics. 2002; 15 :231-254 - 8.
Tumin A, Wang X, Zhong X. Numerical simulation and theoretical analysis of perturbations in hypersonic boundary layers. AIAA Journal. 2011; 49 (3):463-471 - 9.
Wang X, Zhong X, Ma Y. Response of a hypersonic boundary layer to wall blowing–suction. AIAA Journal. 2011; 49 (7):1336-1353 - 10.
Egorov IV, Fedorov AV, Sudakov VG. Direct numerical simulation of disturbances generated by periodic suction-blowing in a hypersonic boundary layer. Theoretical and Computational Fluid Dynamics. 2006; 20 (1):41-54 - 11.
Kovasznay LSG. Turbulence in supersonic flow. Journal of Aeronautical Sciences. 1953; 20 :657-682 - 12.
P. Balakumar. Receptivity of supersonic boundary layers due to acoustic disturbances over blunt cones. 37th AIAA Fluid Dynamics Conference and Exhibit; American Institute of Aeronautics and Astronautics, Miami, Florida; Reston, VA: 2007 - 13.
Kara K, Balakumar P, Kandil OA. Effects of nose bluntness on hypersonic boundary-layer receptivity and stability over cones. AIAA Journal. 2011; 49 (12):593-2606 - 14.
Xin-Liang LI, De-Xun FU, Yan-Wen MA, et al. Acoustic calculation for supersonic turbulent boundary layer flow. Chinese Physical Society. 2009; 26 (9):181-184 - 15.
Ma Y, Zhong X. Receptivity of a supersonic boundary layer over a flat plate, part 2: Receptivity to Freestream sound. Journal of Fluid Mechanics. 2003; 488 (6):79-121 - 16.
Ma Y, Zhong X. Receptivity of a supersonic boundary layer over a flat plate, part 3: Effects of different types of free-stream disturbances. Journal of Fluid Mechanics. 2005; 532 (6):63-109 - 17.
Steger J, Warming RF. Flux vector splitting of the in-viscid gasdynamic equations with application to finite difference methods. Journal of Computational Physics. 1981; 40 :263-293 - 18.
Liu XD, Osher S, Chan T. Weighted essentially non-oscillatory schemes. Journal of Computational Physics. 1994; 115 :200-212 - 19.
Zhang Y, Fu D, Ma Y, Li X. Receptivity to free-stream disturbance waves for hypersonic flow over a blunt cone. Science in China Series G: Physics, Mechanics & Astronomy. 2008; 51 (11):1682-1690 - 20.
Shu CW. Essentially Non-Oscillatory and Weighted Essentially Non-Oscillatory Schemes for Hyperbolic Conservation Laws. NASA/CR-97-206253 - 21.
Tang X, Zhu X, Hui T, et al. Receptivity characteristics of a hypersonic boundary layer under freestream slow acoustic wave with different amplitudes. European Physical Journal Applied Physics. 2017; 79 :31101-1-31101-15 - 22.
Tang X, Lv H, Meng X, Wang Z, Lv Q. Stability characteristic of hypersonic flow over a blunt wedge under freestream pulse wave. Central. European Journal of Physics. 2014; 12 (1):17-31 - 23.
Wang Z, Tang X, Lv H, Shi J. Temporal and spatial evolution characteristics of disturbance wave in a hypersonic boundary layer due to single-frequency entropy disturbance. The Scientific World Journal. 2014; 39 (2):359-371. DOI: 10.1155/2014/517242 - 24.
Xian L, Li XL, Fu DX, Ma YW. Effects of wall temperature on boundary layer stability over a blunt cone at Mach 7.99. Computers & Fluids. 2010; 39 :359-371 - 25.
Saric WS, Reed HL, Kerschen EJ. Boundary-layer receptivity to freestream disturbances. Annual Review of Fluid Mechanics. 2002; 34 :291-319