1. Introduction
The laminar-to-turbulent transition in hypersonic flows is vital for aerodynamic and aerothermal vehicle design because it is accompanied by a significant increase in skin friction and heat transfer. Nevertheless, accurately predicting the transition position remains exceedingly challenging, as this process is markedly influenced by external factors such as noise and turbulence levels and wall roughness. Morkovin [
1] systematically summarized the potential existence of multiple boundary layer development transition pathways, with the natural transition process being the most probable in actual flight conditions. Generally, the natural transition process can be divided into four stages: receptivity, linear instability, nonlinear instability, and turbulence.
As the first transition stage, receptivity is the process by which external disturbances excite mode disturbances within the boundary layer. These external disturbances can arise from freestream or localized surface imperfections. Receptivity can be classified into natural and forced receptivity based on different external disturbance types. Exploring the receptivity process is fundamentally directed at elucidating the “scale conversion” relationship between external and mode disturbances within the boundary layer.
For natural receptivity, external disturbances can be acoustic, entropy, and vorticity waves. Goldstein [
2] computed receptivity coefficients for acoustic waves at a zero incidence angle for a flat-plate boundary layer, while Kerschen [
3] investigated the receptivity of vorticity waves to mode disturbances within the boundary layer. Hammerton’s [
4,
5] works extended and applied the receptivity theory from flat plates to parabolic leading edges. Notably, these studies primarily centered around subsonic flow. However, the boundary layer’s strong compressibility introduces new receptivity mechanisms for hypersonic flows. The best known is the mode synchronization theory proposed by Fedorov [
6], suggesting two disturbance mode synchronizations, near the leading edge and at a location downstream, which completely excite the most unstable Mack mode. This theory has been confirmed by a direct numerical simulation conducted by Zhong [
7].
Unlike natural receptivity, forced receptivity involves directly exciting instability modes through wall surface actions like periodic blowing–suction, local temperature changes, and wall roughness. Theoretical studies have found that even very small localized imperfections can trigger significant local receptivity, so this mechanism is important. The “finite Reynolds number theory” developed by Goldstein [
8,
9,
10] is well-known in forced receptivity research, and relevant studies can be carried out using this theory whether the wall disturbance excitation is isolated or distributed.
In the hypersonic boundary layer context, Ruban [
11] studied steady longitudinal vortex excitation using spectral methods, while Duan [
12] explored forced receptivity over flat-plate boundary layers using direct numerical simulation. Beyond the “finite Reynolds number theory” and direct numerical simulations, a practical compromise involves solving adjoint equations. This method offers computational efficiency and mathematical rigor, making it accessible to a wider audience. Based on this approach, Tempelmann [
13] examined crossflow mode receptivity on a subsonic swept wing driven by roughness. Wang [
14] investigated plasma actuation’s receptivity at the wall, and Xi [
15] explored the stability and receptivity of a swept blunt body.
In addition, the friction between the vehicle surface and the air generates elevated temperatures as freestream Mach numbers increase, reaching into the range of thousands or even tens of thousands of degrees within the boundary layer. This region experiences various thermal and chemical processes, causing significant real gas effects, and the traditional assumption of a calorically perfect gas (CPG) fails. To capture the flow’s attributes with precision, researchers have employed a variety of gas models, including the thermally perfect gas (TPG) model accounting for vibrational energy excitation, the chemical nonequilibrium (CNE) gas model, the thermochemical nonequilibrium (TCNE) gas model, and the thermochemical equilibrium (TCE) gas model, assuming both thermal and chemical equilibrium states are achieved [
16].
Over the past two decades, numerous studies have focused on the real gas effect implications on flow stability. Malik [
17,
18] investigated the impact of dissociation effects on flow stability using a self-similar solution of a flat-plate boundary layer for a thermochemical equilibrium gas. They discovered that the growth rate of the second mode was higher with the chemical equilibrium model, and the range of unstable frequencies shifted toward lower frequencies. Stuckert [
19,
20] conducted a similar study on the boundary layer using a chemical nonequilibrium model considering finite chemical reaction rates. Their findings suggested that heat absorption through chemical reactions expands the supersonic region within the boundary layer, reducing the second mode’s frequency.
Both thermochemical and various thermodynamic and transport models can significantly impact flow stability [
21,
22,
23,
24,
25]. Miró Miró [
25] systematically compared the effects of these models on flow stability. In addition to the commonly observed second mode, supersonic modes in high-temperature boundary layers have garnered considerable attention [
26,
27,
28]. In calorically perfect gas cases, supersonic modes are limited to extremely cold wall conditions. However, considering real gas effects appears to facilitate the occurrence of supersonic modes. Research indicates that significant nose bluntness can elevate the instability of the supersonic mode even beyond the second mode, making it the dominant mode [
28]. While most prior studies focused on the linear growth stage of disturbances, recent research has considered the nonlinear stage, exploring the impact of real gas effects using methods like direct numerical simulation (DNS) [
29] and nonlinear parabolized stability equations (NPSE) [
30,
31].
There has been extensive research on real gas effects on the growth stage of disturbances. However, few studies have focused on the impact of a real gas on the receptivity stage. To accurately predict the transition position, it is crucial to consider not only the growth stage but also the initial amplitude of the disturbance, which is an essential factor influencing the transition position. In the case of blunt body swept flows, the dominant unstable mode is the crossflow mode, which can be categorized into two classical types: traveling and stationary modes, based on their frequencies. The former is associated with the interaction of freestream disturbances and wall roughness, while the latter is only related to wall roughness. These make the receptivity coefficient of a stationary mode significantly larger than that of a traveling mode. As a result, the stationary mode is typically more critical than the traveling mode because it has a significantly larger receptivity amplitude.
This research on the effect of the real gas on the receptivity of hypersonic stationary crossflow modes to roughness is the primary focus. This paper is structured as follows:
Section 2 introduces the governing equations and methodologies for solving for the base flow, stability, and receptivity.
Section 3 examines the influence of real gas effects and the Mach number on the stability and receptivity of stationary crossflow modes and provides predictions for transition positions. Finally, the conclusions from this study are presented in
Section 4.
3. Results and Discussion
3.1. Baseflow and Linear Stability Analysis
In this paper, the maximum Mach number is 20, but the maximum streamwise velocity component u is limited to about 4200 m/s due to the sweep angle. We calculated the Ma = 20 case using two air reaction models, 5-species and 11-species. It was found that the ionization reaction has not yet occurred at this time, and the flow fields obtained under the two models largely overlap. Therefore, the 5-species air reaction model was chosen for subsequent calculations, which is appropriate for the cases examined in this paper.
Figure 3 shows the temperature distributions around the nose region at different Mach numbers, highlighting the significant impact of the chemical nonequilibrium effect on the flow field characteristics. The endothermic chemical reaction significantly reduces the temperature near the stagnation point and decreases the shock stand-off distance.
Figure 4 shows the oxygen/nitrogen mass fraction distribution around the nose. The series of chemical reactions in the air have not yet occurred when the temperature is below 2500 K. It can be seen that when Ma = 10, most of the flow field remains chemically frozen. The results of the chemical nonequilibrium model exhibit only slight deviations from the calorically perfect gas model in this case. As the Mach number increases, chemical reactions occur, with oxygen dissociating initially. At Ma = 20, oxygen undergoes nearly complete dissociation, while nitrogen begins to dissociate.
Figure 5 illustrates the temperature and velocity profiles in the boundary layer for various Mach numbers. It is apparent that the chemical nonequilibrium effect reduces both the temperature and velocity within the boundary layer, leading to modifications in the flow properties at the boundary layer’s edge. The variation in flow characteristics at the boundary layer’s edge can be attributed to changes in the shock wave intensity. As shown in
Figure 6a, the shock angle for the CNE model in the downstream region significantly deviates from that of the CPG model. The post-shock parameters also exhibit variations because of the differing shockwave intensities between the two models.
Figure 6b shows the boundary layer momentum thickness, indicating that there is no noticeable difference in the momentum thicknesses for different gas models when the Mach number is small. However, as the Mach number increases, the momentum thickness for the CNE model decreases noticeably. Some previous studies also considered real gas effects on boundary layer flow, but the wall conditions were mostly adiabatic. These studies discovered that chemical reactions contribute to a decreased wall temperature, exhibiting an influence similar to wall cooling. This study finds that chemical reactions cooled the flow within the boundary layer for isothermal walls, decreasing its thickness despite the constant wall temperature.
The unstable growth rate distributions for various Mach numbers and gas models are shown in
Figure 7, including the second mode and stationary and traveling crossflow modes. The unstable frequency range and growth rate of the second mode appear to escalate, while the crossflow mode stabilizes with an increasing Mach number. As shown in
Figure 7b–e and c–f, the chemical nonequilibrium effect further destabilizes the second mode and crossflow modes. The growth rate peaks for the CNE model are notably larger than those of the CPG model. When the Mach number reaches 15, the peak growth rate of the second mode significantly surpasses that of the crossflow mode.
Figure 8 shows the N-factor envelope of the second and stationary crossflow modes for the CNE model at Ma = 20. While the second mode exhibits a higher peak growth rate at this point, its shorter integration path results in a smaller N-factor. By contrast, the stationary crossflow mode’s longer growth path allows it to dominate downstream transitions.
Figure 9 illustrates the N-factor of stationary crossflow modes at various Mach numbers.
Figure 9a,b shows the N-factor envelope curves for the CPG and CNE models, respectively, while
Figure 9c presents the N-factor at
x = 2 m. As shown in
Figure 9c, the N-factor decreases with increasing Mach number in the calorically perfect gas model. In other words, increasing the Mach number makes the stationary crossflow modes more stable, referred to as the Mach number effect.
Chemical nonequilibrium effects become non-negligible with increasing Mach number. The N-factor for the CNE model is significantly larger than that of the CPG model at the same Mach number, indicating that the chemical nonequilibrium effect further destabilizes the stationary crossflow modes. For the CNE model, a non-monotonic N-factor trend is observed as the Mach number increases. Initially, the N-factor decreases as the Mach number increases from 10 to 12. At this stage, the Mach number is relatively low, chemical nonequilibrium effects are weak, and the Mach number effects are dominant.
As the Mach number increases to 17, the influence of chemical nonequilibrium effects becomes substantial, resulting in an increase in the N-factor with rising Mach number. However, as the Mach number further increases to 20, the N-factor declines and approaches the value observed at Ma = 15. This phenomenon could be attributed to the Mach number effect’s enhanced stabilizing impact on the stationary crossflow modes despite intensifying chemical nonequilibrium effects as the Mach number continues to rise. At this point, the Mach number effect regains dominance. In essence, the Mach number effect and the chemical nonequilibrium effect exhibit opposing influences on the stability of the stationary crossflow modes. These two effects compete, ultimately resulting in a non-monotonic variation of the N-factor in the CNE model.
3.2. Receptivity
We assume a Gaussian roughness distribution in the streamwise direction with a periodic distribution in the spanwise direction,
where
h is the roughness height,
a is the shape factor,
x0 is the center position of the roughness, and
β is the spanwise wavenumber. The specific parameters associated with these distributions are presented in
Table 2, where
δ is the scale of roughness in the streamwise direction.
The typical roughness height is in the micron range. In this paper, the roughness height is about 1% of the boundary layer thickness, permitting the use of boundary compensation methods to address the small-scale roughness. Specifically, the roughness surface (
y =
ε) retains the no-slip and isothermal boundary conditions, while a first-order Taylor expansion is implemented at the original wall (
y = 0), resulting in equivalent boundary conditions at the wall,
Figure 10 illustrates the stationary crossflow mode receptivity coefficients obtained using ALST for both gas models. With the CPG model, the receptivity coefficient increases with the Mach number. Conversely, the CNE model demonstrates a rapid decline in the receptivity coefficient as the Mach number increases. The receptivity coefficient variation with Mach number at
β = 60 m
−1 is depicted in
Figure 10c. The receptivity coefficients from the CNE model are significantly smaller than those from the CPG model at the same Mach number, indicating that the chemical nonequilibrium effect significantly reduces the receptivity of the stationary crossflow modes.
Figure 10c also presents the receptivity coefficients for both roughness distributions, showing that the receptivity coefficient increases when roughness is closer to the upstream.
It is worth noting that when Ma = 10–12, the receptivity coefficient of the CNE model is larger than that of the CPG model. However, based on the preceding analysis of the base flow, it is evident that the boundary layer is mostly chemically frozen at low Mach numbers. Therefore, the observed increase in the receptivity coefficient for the CNE model is not attributable to chemical nonequilibrium effects. In addition to chemical reactions, the two gas models differ because the CPG model does not consider vibrational energy excitation, and their transport coefficients vary. However, exactly which effect leads to an increase in the stationary crossflow modes’ receptivity coefficients at low Mach numbers requires further study.
Figure 11 shows the amplitude evolutions of stationary crossflow modes with different spanwise wavenumbers excited by roughness elements, and the amplitude envelope is also provided. Although the chemical nonequilibrium effect significantly reduces the receptivity of the stationary crossflow modes, it also results in a larger growth rate. At a Mach number of 20, the initial disturbance amplitude is low for the CNE model. However, upon exiting the excitation section, the disturbance amplitude rapidly increases, surpassing that of the CPG model, as shown in
Figure 11b.
Figure 12 shows the amplitude envelope for different Mach numbers. The trend of the amplitude envelope is similar to the N-factor envelope. For the CPG model, the amplitude gradually declines as the Mach number increases, whereas the CNE model exhibits a non-monotonic amplitude trend with increasing Mach number
Figure 13 illustrates the transition positions based on A = 5% and N = 8 predictions, with the latter not accounting for the effect of receptivity. Both methods show similar trends. The receptivity difference influences the transition position to a certain degree, but the more critical factors are the effects of Mach number and gas model on the stability of stationary crossflow modes. Chemical nonequilibrium effects appear for the blunt body swept flow dominated by stationary crossflow modes as the Mach number increases, working in conjunction with the Mach number effect to cause two transition position reversals.
At lower Mach numbers, the transition experiences a delay as the Mach number increases, followed by an advancement and subsequent delay in the transition position as the Mach number continues to rise. The CPG model predicts a transition position close to that of the CNE model at lower Mach numbers. However, when the Mach number exceeds 12, the CPG model predicts a transition position farther downstream, deviating significantly from the CNE model. This finding implies that it is crucial to account for the influence of chemical nonequilibrium effects when dealing with high-Mach-number conditions. Neglecting these effects may lead to erroneous predictions of the transition position.
4. Conclusions
This paper investigates the effects of the gas model and Mach number on the stability and receptivity of stationary crossflow modes using linear stability analysis (LST) and its adjoint (ALST). It was found that chemical nonequilibrium effects significantly alter the shock angle and, consequently, the flow quantities at the boundary layer’s edge. As the Mach number increases, the second mode gradually becomes unstable while the crossflow mode stability improves. The growth rate of the second mode exceeds that of the crossflow mode after the Mach number exceeds a particular critical value. However, the stationary crossflow modes always dominate the transition because of their longer growth paths. Chemical nonequilibrium effects further destabilize the second and crossflow modes. The Mach number and chemical nonequilibrium effects have opposite effects on the stability of stationary crossflow modes. These two effects alternately dominate with increasing Mach number.
The stationary crossflow modes are excited by distributed roughness. The results indicate that the receptivity coefficients of stationary crossflow modes gradually increase with increasing Mach number, while the chemical nonequilibrium effect causes a significant decrease in the receptivity coefficients of stationary crossflow modes.
The combined influence of the Mach number and the gas model on receptivity and stability was examined by comparing disturbance amplitude envelopes. It was found that the difference in receptivity only affects the transition position to a certain extent. The effects of Mach number and gas model on the stability of the constant crossflow modes are more important. In a blunt body swept flow dominated by stationary crossflow modes, chemical nonequilibrium and Mach number effects combine to result in two reversals in the transition position. This study highlights the significant influence of real gas effects on the receptivity and stability of stationary crossflow modes. Specifically, it is necessary to consider the influence of chemical nonequilibrium effects at high Mach numbers to prevent incorrect predictions of the transition position.