Oblique transition in hypersonic double-wedge flow

We utilize resolvent and weakly nonlinear analyses in combination with direct numerical simulations (DNS) to identify mechanisms for oblique transition in a Mach 5 hypersonic flow over an adiabatic slender double-wedge. Even though the laminar separated flow is globally stable, resolvent analysis demonstrates significant amplification of unsteady external disturbances to the linearized flow equations. These disturbances are introduced upstream of the separation zone and they lead to the appearance of oblique waves further downstream. We demonstrate that large amplification of oblique waves arises from the growth of fluctuation shear stress due to streamline curvature of the laminar base flow in the separated shear layer. This is in contrast to the attached boundary layers, where no such mechanism exists. We also use a weakly nonlinear analysis to show that the resolvent operator associated with linearization around the laminar base flow governs the evolution of steady reattachment streaks that arise from quadratic interactions of unsteady oblique waves. These quadratic interactions generate vortical excitations in the reattaching shear layer which lead to the formation of streaks in the recirculation zone and their subsequent amplification, breakdown, and transition to turbulence downstream. Our analysis of the energy budget shows that deceleration of the base flow near reattachment is primarily responsible for amplification of steady streaks. Finally, we employ DNS to examine latter stages of transition to turbulence and demonstrate the predictive power of a weakly nonlinear input-output framework in uncovering triggering mechanisms for oblique transition in separated high-speed boundary layer flows.


Introduction
Slender double-wedges are commonly encountered in intakes, control surfaces, and junctions in high-speed supersonic and hypersonic vehicles (Dolling 2001). In this geometry, laminar boundary layer can separate at the corner because of the pressure rise that arises from deflection of the inviscid free stream. The resulting flow is characterized by separation-reattachment shocks as well as a recirculation zone and it provides a canonical setup for studying shock-wave-boundary-layer interaction (SWBLI) (Simeonides & Haase 1995). In spite of spanwise homogeneity of laminar base flows over compression corners, both experiments Roghelia et al. 2017;Dwivedi et al. 2020a) and numerical simulations (Navarro-Martinez & Tutty 2005;Dwivedi et al. 2017;Cao et al. 2021b) identify three-dimensional (3D) features in time-averaged separated flow. In particular, streamwise streaks associated with persistent local peaks of heat flux or wall temperature, that appear near reattachment, can trigger transition to turbulence downstream (Simeonides & Haase 1995;Roghelia et al. 2017).
The development of 3D flow structures in hypersonic flows was recently studied by examining the growth of small perturbations in the presence of a recirculation zone (Dwivedi 2020). For example, 2D SWBLI can become unstable inside the separation bubble when the strength of interaction increases beyond a critical value (Sidharth et al. 2017). The spanwise modulation that arises from global instability introduces streaks over compression corners (Sidharth et al. 2018) as well as oblique shocks impinging on a flat plate (Hildebrand et al. 2018) and it can trigger transition to turbulence (Cao et al. 2022). Similar 3D flow features have also been observed in hypersonic regimes where non-continuum effects are important (Sawant et al. 2022). However, recent numerical simulations and global stability analysis demonstrate that hypersonic compression corner flows can be stabilized by increasing radius of the leading edge (i.e., its bluntness) (Cao et al. 2021a) or by increasing the wall temperature (Hao et al. 2021).
Even in the absence of global instability, high-speed separated flows are highly sensitive to upstream vortical disturbances (Dwivedi et al. 2019) and small fluctuations around the laminar 2D base flow can experience significant non-modal amplification that leads to the appearance of steady reattachment streaks (Dwivedi et al. 2020b). Furthermore, recent experiments on the cone-flare configurations (Benitez et al. 2020;Butler & Laurence 2021), which represent axisymmetric counterparts of slender doublewedges, identify unsteady fluctuations in the separation zone. These fluctuations are significantly amplified in the recirculation zone and they play an important role in transition to turbulence (Butler & Laurence 2021).
In this paper, we examine amplification of unsteady fluctuations around the laminar 2D base flow in the separation/reattachment zone and investigate subsequent transition to turbulence. Free-stream disturbances (Choudhari 1996;Maslov et al. 2001) that arise from wind tunnel noise in ground experiments (Schneider 2015) or from atmospheric disturbances in free flights (Bushnell 1990;Skinner et al. 2020) can lead to the appearance of unsteady fluctuations in boundary layer flows. It is well-documented that unsteady oblique waves provide a potent mechanism for initiating transition in low-speed incompressible Rigas et al. 2021) and compressible (Chang & Malik 1994;Mayer et al. 2011) boundary layers. Even though the importance of oblique fluctuations in initiating transition in attached high-speed boundary layers has received significant attention (Ma & Zhong 2005;Sivasubramanian & Fasel 2015;Hader & Fasel 2019), their role in separated high-speed flows has not been studied. Recent experiments (Benitez et al. 2020) suggest that their amplification within the recirculation zone can trigger unsteadiness in transitional SBWLI flows. We utilize global resolvent and weakly nonlinear analyses to quantify amplification of unsteady upstream disturbances in a Mach 5 flow over a slender double-wedge and characterize their role in initiating transition to turbulence in high-speed separated boundary layers.
Resolvent analysis provides a framework for evaluating responses (outputs) of stable dynamical systems to time-periodic external disturbances (inputs) (Trefethen et al. 1993;Schmid & Henningson 2001;Schmid 2007). For time-independent globally stable base flows, the steady-state response of the linearized Navier-Stokes (NS) equations to a harmonic input with frequency ω is also harmonic with the same frequency and the frequency response operator maps the input forcing to the resulting steady-state output (Jovanović 2021). The singular value decomposition (SVD) of the frequency response characterizes amplification across frequency ω and decomposes inputs and outputs into modes whose significance is ordered by the magnitude of the corresponding singular values (Schmid 2007). In addition to providing insights into dynamics of canonical incompressible flows (Jovanović & Bamieh 2005;McKeon & Sharma 2010;Brandt et al. 2011;Sipp & Marquet 2013;Ran et al. 2019a,b), input-output analysis has also been utilized to discover mechanisms for noise generation in turbulent jets (Garnaud et al. 2013;Jeun et al. 2016;Schmidt et al. 2018), separation control on airfoils (Yeh & Taira 2019), and the appearance of reattachment streaks in hypersonic flows (Dwivedi et al. 2019).
In a Mach 5 double-wedge flow subject to unsteady disturbances, we employ resolvent analysis to demonstrate that oblique waves represent the most energetic response of the compressible linearized NS equations. We utilize the compressible energy norm (Chu 1965;Hanifi et al. 1996) to quantify energy amplification and show that unsteady upstream disturbances that are localized before flow separation induce oblique waves downstream of the double-wedge corner. Our analysis of the transport equation for the streamwise specific kinetic energy of oblique waves reveals that concave flow curvature of the separated/reattaching laminar 2D base flow is the primary source of amplification in the presence of SWBLI. We also utilize a weakly nonlinear analysis to demonstrate that quadratic interactions of oblique waves generate vortical excitations that induce reattachment streaks in the recirculation bubble. We show that the resolvent operator associated with linearization around the laminar 2D base flow governs the evolution of steady reattachment streaks and use SVD to demonstrate that the streaks are well approximated by the second output resolvent mode. Our analysis of the energy budget shows that the base flow deceleration near reattachment is primarily responsible for amplification of reattachment streaks. Finally, we conduct DNS to confirm the predictive power of our approach and provide insight into latter stages of transition to turbulence.
Recently, Rigas et al. (2021) utilized a variational framework to extend input-output analysis in the frequency domain to the nonlinear NS equations. For fundamental forcing, the disturbance that triggers transition and yields the largest skin-friction coefficient in incompressible boundary layer is given by a pair of oblique waves with temporal frequency and spanwise wavenumber which are very close to the ones identified by the resolvent analysis of the linearized NS equations (Rigas et al. 2021). While the Orrmechanism (Schmid & Henningson 2001) and the Tollmien-Schlichting linear instability (Sipp & Marquet 2013) initiate the early stages of oblique transition in the attached low-speed boundary layers (Rigas et al. 2021), even linear amplification mechanisms are poorly understood in separated compressible flows. Recent numerical simulations with inlet stochastic excitations in axisymmetric cylinder flare geometry showed that the local "first mode" instability (Mack 1984) can cause growth of oblique waves in the upstream boundary layer (i.e., before separation) and initiate transition in high-speed compressible flows with SWBLI (Lugrin et al. 2021). However, in the absence of local and global boundary layer instabilities, the role of flow separation in amplification of unsteady fluctuations and the ensuing transition has not been previously investigated. Figure 1 provides a summary of our key findings. We utilize resolvent analysis of the laminar 2D base flow with spatially localized forcing introduced in the streamwise plane immediately upstream of separation to identify the spatial structure of unsteady external disturbances that yield the most energetic response of the compressible linearized NS equations. The resulting forcing is used to trigger non-modal amplification of oblique waves in the separated shear layer and generate steady reattachment streaks, which are routinely observed in SWBLI experiments, further downstream through weakly nonlinear interactions. Interaction of streaks with oblique waves is observed after reattachment and DNS is used to demonstrate that unsteady upstream oblique disturbances can indeed trigger transition to turbulence in separated high-speed compressible flows.
Our presentation is organized as follows. In § 2, we introduce the slender doublewedge geometry along with a finite-volume compressible flow solver that we use in our computations. In § 3, we describe resolvent and weakly nonlinear analyses that we use to evaluate frequency responses of the double-wedge flow in the presence of 3D disturbances. We also utilize resolvent analysis to demonstrate large amplification of unsteady oblique disturbances to the linearized flow equations and identify the underlying physical mechanism. In § 4, we employ a weakly nonlinear analysis to demonstrate that Figure 1: Preview of key results: spatially localized unsteady upstream forcing triggers oblique waves in the separated shear layer and their quadratic interactions lead to the appearance of steady reattachment streaks further downstream. DNS validates our theoretical predictions and demonstrates the efficacy of unsteady oblique disturbances in triggering transition in globally stable separated high-speed boundary layer flows. quadratic interactions of oblique waves induce steady reattachment streaks and discuss physical mechanism responsible for their amplification in recirculation zone. In § 5, we employ DNS to validate utility of our theoretical predictions and examine latter stages of transition induced by unsteady upstream disturbances. In § 6, we analyze statistical features of the resulting transitional and turbulent boundary layers. We provide summary of our contributions and offer concluding remarks in § 7.

Hypersonic flow over an adiabatic slender double-wedge
Hypersonic flow over a slender double-wedge with free-stream conditions shown in figure 2 corresponds to the experiments of Yang et al. (2012). Since the enthalpy and the temperature in the flow field are low, we utilize the ideal gas law abstraction and employ the finite-volume compressible flow solver US3D (Candler et al. 2015a) to compute the solution of the compressible NS equations in conservative form, Here, F(Ψ ) := −∇ · F(Ψ ) is the dynamical generator of the compressible NS equations, F(Ψ ) is the flux vector, ∇ is the gradient, and Ψ = (ρ, ρu, E t ) is the vector of conserved variables that represent density, momentum, and total energy per unit volume of the gas. In equation (2.1) and throughout the paper, spatial coordinates are non-dimensionalized by the boundary layer thickness at separation δ 99 = 9.8 × 10 −4 m, velocity by the freestream velocity U ∞ , pressure by p ∞ , temperature by T ∞ , and time by δ 99 /U ∞ . The Reynolds number based on separation boundary layer thickness is 13.3 × 10 3 and the Mach number is 5.0.  We discretize the inviscid fluxes using the second-order-accurate modified Steger-Warming fluxes with the MUSCL limiters (Candler et al. 2015a). In previous studies, the numerical method for the computation of the base state was validated using hypersonic double-wedge and double-cone setups (Nompelis et al. 2003;Nompelis & Candler 2009). The laminar 2D base flowΨ is computed as the steady-state solution of equation (2.1), by implicit time-marching with 249 cells in the normal and 535 cells in the streamwise direction. As illustrated in Sidharth et al. (2018), this resolution is sufficient to capture separated flow and resolve the evolution of small perturbations. Figure 3 shows the contours of density gradient magnitude on the compression corner. The separation and the reattachment locations in the laminar 2D base flowΨ are marked by S and R, respectively. The 2D flow separates upstream of the corner on the first wedge, it reattaches on the second wedge, and the separated and reattaching shear layers are respectively associated with the formation of the separation and reattachment shocks. Figure 3 also shows an inset of the separation zone along with the wall-aligned coordinate system (where ξ and η denote directions that are parallel and perpendicular to the wall) and a coordinate system that is locally aligned with the streamlines of the laminar 2D flow. Both coordinate systems are used in our study of the evolution of flow fluctuations. Sidharth et al. (2018) demonstrated global linear stability of the laminar 2D base flowΨ . Recent studies of similar SWBLI configurations, such as compression ramps, revealed extreme sensitivity to upstream disturbances even in the absence of global instability (Dwivedi et al. 2019). Leading-edge roughness and free-stream disturbances provide persistent sources of external excitation and they are inevitable in realistic flows. To evaluate the role of such uncertainty in triggering early stages of transition to turbulence, we utilize input-output analysis to quantify amplification of unsteady disturbances in a hypersonic flow over slender double-wedge.
3. Input-output analysis of a high-speed double-wedge flow In this section, we employ input-output analysis to quantify amplification of small unsteady external disturbances in globally stable 2D SWBLI over a slender double-wedge and uncover physical mechanisms that trigger early stages of transition to turbulence.

Externally forced compressible NS equations
To account for the rate of change of perturbation density, momentum, and total energy, we model unsteady external disturbances d(x, t) to the compressible NS equations in the conservative form (2.1) as volumetric sources of excitation, where x := (x, y, z) is the vector of streamwise, normal, and spanwise spatial coordinates. By decomposing the flow field Ψ into the sum of the baseΨ and fluctuating ψ parts, we obtain the equation that governs the dynamics of flow fluctuations aroundΨ (x), For disturbances with small amplitude , a weakly nonlinear analysis can be utilized to examine externally forced compressible NS equations (3.3) and determine the corrections to the steady laminar 2D base flowΨ (x). Up to a second order in , the vector of flow fluctuations ψ can be represented as, where ψ (1) (x, t) satisfies the linearized flow equations Equations (3.5a) and (3.5b) are respectively obtained by neglecting O( 2 ) and O( 3 ) terms upon substitution of (3.4a) and (3.4b) to the compressible NS equations (3.3). The compressible NS operator resulting from linearization around the base flowΨ is determined by A(Ψ ) (see Candler et al. (2015b, equation (23)) and Sidharth et al. (2018, equations (A1)-(A2))) and N (2) (ψ (1) ) accounts for quadratic interactions at O( 2 ); see appendix A for details.
Several recent studies demonstrated the utility of the compressible energy norm (Chu 1965;Hanifi et al. 1996), in quantifying growth of fluctuations in high-speed boundary layer flows (Franko & Lele 2013;Sidharth et al. 2018;Quintanilha et al. 2022). This quantity is determined by the weighted L 2 norm of the vector of flow fluctuations φ := (φ 1 , φ 2 , φ 3 ) = (ρ , u , T ) in primitive variables, where · , · 2 is the standard L 2 inner product over the domain Ω, C v is the specific heat at constant volume in Ω, and is the multiplication operator determined by the pressurep, densityρ, and temperaturē T of the base flowΨ . For small amplitude disturbances, we can represent φ as where, (3.7b) As shown in appendix B, C and D are multiplication operators parameterized by the laminar 2D base flowΨ ; see equation (B 4) for their definition.
In the remainder of this section, we identify oblique waves as the most energetic responses of the linearized flow equations (3.5a) in the presence of unsteady harmonic disturbances d (1) . In § 4, we utilize equation (3.5b) to demonstrate that steady streaks can arise from quadratic interactions of unsteady oblique waves.

Amplification of exogenous disturbances to the linearized flow equations
The linearized equations (3.5a) describe the evolution of fluctuation vector ψ (1) in the presence of external disturbances d (1) and they can be written using the state-space formulation (Schmid & Henningson 2001), (3.8) Here, d (1) is a spatially distributed and temporally varying disturbance (input), ψ (1) is the state of the linearized system (which is determined by the vector of flow fluctuations in conserved variables), φ (1) is the quantity of interest (output) whose weighted L 2 norm determines energy of flow fluctuations (3.6), and A is the generator of the linearized compressible NS dynamics. The input operator B in (3.8) allows to specify spatial support of body forcing inputs and the output operator C relates the state of the linearized system ψ (1) to the vector of flow fluctuations in primitive variables φ (1) .
For the parameters shown in figure 2, the linearized system is globally stable and for a time-periodic input with frequency ω, d (1) (x, t) =d (1) (x, ω)e iωt , the steady-state output H(ω) is the frequency response operator, and R(ω) = (iωI − A) −1 is the resolvent operator associated with the linearized model (3.8). At any ω, the singular value decomposition of H(ω) can be used to quantify amplification of time-periodic inputs (Schmid 2007;Jovanović 2021), where σ i (ω) denotes the ith singular value of H(ω), · , · E is the inner product in (3.6) that induces the compressible energy norm, and d i (x, ω) and φ i (x, ω) are the left and right singular functions of H(ω) which provide orthonormal bases of the corresponding input and output spaces (with respect to · , · E ). The frequency response operator H(ω) maps the ith input mode d i (x, ω) into the response whose spatial profile is specified by the ith output mode φ i (x, ω) and the amplification is determined by the corresponding singular value σ i (ω). In other words, for Note that, at any ω, determines the largest induced gain with respect to a compressible energy norm, where (d 1 (x, ω), φ 1 (x, ω)) identify the spatial structure of the dominant input-output pair. We use a second-order central finite-volume discretization (Sidharth et al. 2018) to obtain a finite dimensional approximation of the linearized model (3.8) and employ matrix-free Arnoldi iterations (Jeun et al. 2016;Dwivedi 2020) to compute the singular values σ i (ω) of H(ω).

Frequency response analysis
We utilize the resolvent analysis to study amplification of harmonic disturbances with frequency ω to the linearized flow equations. In double-wedge geometry, the laminar 2D base flowΨ is a function of streamwise normal coordinates,Ψ (x, y), and owing to homogeneity in the spanwise direction, the 3D fluctuations in (3.5a) take the form, (3.14) where β = 2π/λ z is the spanwise wavenumber. Thus, in addition to ω, the frequency response operator is also parameterized by β, where A β denotes the Fourier symbol of the operator A in (3.8) obtained by replacing the spanwise differential operator ∂ z with iβ. At any pair (ω, β), H β (ω) maps the input functiond (1) of x and y into the output functionφ and SVD of H β (ω) can be used to study amplification across spatio-temporal frequencies.
We first set B = I, i.e., we introduce body forcing inputs to excite flow at every spatial location in the computational domain Ω and we choose the output operator C to examine the impact of forcing on the compressible energy norm of φ (1) in the entire Ω. The resolvent analysis is done using a resolution that yields grid-independent outputs with 545 cells in the streamwise, 249 cells in the normal direction, and numerical sponge boundary conditions near the leading edge (x = 1) and the outflow boundaries are utilized. Figure 4 shows the dependence of the input-output gain G(ω, λ z ) on the frequency ω and the wavelength λ z . There are two major amplification regions with the respective peaks at (ω = 0, λ z = 1.5) and (ω = 0.4, λ z = 3). The first peak in G identifies the largest amplification and the corresponding output is determined by reattachment streaks that result from steady vortical disturbances upstream of the recirculation zone. We observe selective amplification of disturbances with λ z ≈ 1.5 and low-pass filtering features over ω. The gain G experiences rapid decay beyond the roll-off frequency ω ≈ 0.4 and it attains its largest value at ω = 0 for λ z that scales with the reattaching shear layer thickness (Dwivedi et al. 2019). In contrast to Dwivedi et al. (2019), which focused on disturbances with ω = 0, we examine unsteady disturbances that trigger oblique waves in the reattaching shear layer, as identified by the second peak in G. This amplification region takes place in a narrow band of temporal frequencies ω over a fairly broad range of spanwise wavelengths λ z .
As demonstrated in figure 4, even when we allow disturbances to enter through the entire computational domain the largest amplification is caused by inputs that are localized upstream of the corner and the resulting response is localized downstream of the corner. The upstream disturbances are the most effective way to excite the flow because of large convection velocity of the laminar 2D base flow (Chomaz 2005;Schmid 2007) and the dominant output emerges in the separated and the reattached regions.  Figure 5: Spatial input-output analysis: (a) input is introduced at a streamwise location x in = 25 before separation and output is evaluated at x out ; (b) dependence of the inputoutput gain G xout on the streamwise location x out for streaks and oblique waves. Unsteady oblique waves with (ω = 0.4, λ z = 3) are strongly amplified throughout the separation zone.
Experimental studies of oblique transition in channel and boundary layer flows (Elofsson & Alfredsson 1998; Berlin & Henningson 1999) often utilize streamwise localized disturbances and a common criticism of the resolvent analysis is that the identified global input modes represent excitation sources that are not easy to realize experimentally. In contrast, traditional approach to the analysis of boundary layers utilizes spatially localized fluctuation sources and evaluates the streamwise growth of fluctuations as they convect downstream (Herbert 1997). However, in the presence of flow separation a parabolized approximation of the NS equations cannot be made. To evaluate amplification in different spatial regions, we restrict inputs and outputs to belong to a plane but still account for the global nature of the separated flow through the resolvent operator (iωI − A β ) −1 . As illustrated in figure 5(a), this is accomplished via a proper selection of the operators B and C in equation (3.8) by fixing the input location before flow separation, at x in = 25, and by evaluating the output at different locations downstream of the separation, x out . In this setup, G xout quantifies the largest amplification at x out of disturbances that are introduced at x in = 25. Figure 5(b) shows the dependence of G xout on x out for upstream disturbances with (ω = 0, λ z = 1.5) and (ω = 0.4, λ z = 3). The gain associated with the steady fluctuations begins to grow in the latter half of the recirculation zone, especially near the reattachment location. In contrast, unsteady perturbations with ω = 0.4 experience significant amplification throughout the separation zone. This observation suggests that the separation zone plays a critical role in the amplification of unsteady fluctuations.
3.4. Amplification of oblique waves: physical mechanism We now analyze physical mechanisms responsible for amplification of flow fluctuations within the separation zone in the presence of upstream unsteady disturbances. In particular, we examine the global response of the linearized equations to the input with (ω = 0.4, λ z = 3.0) introduced prior to separation (i.e., at x in = 25) that triggers the largest amplification in the entire domain Ω. The spatial structure of flow fluctuations is studied in the (s, n, z) coordinate system which is locally aligned with the streamlines of the laminar 2D base flowΨ ; see figure 3 for an illustration. In this frame of reference, u = (ū s , 0, 0) withū s 0, and, as discussed in Finnigan (1983) The streamwise specific kinetic energy E s := u s u s obeys the transport equation, where P, S, V, and K are the production, source, viscous, and curvature terms (see appendix C), and F is the work done by external disturbances (Dwivedi et al. 2019, appendix C). The production term P quantifies interactions of fluctuation stresses with the base flow gradients, the source term S corresponds to the perturbation component of the inviscid material derivative, the viscous term V determines dissipation of kinetic energy by viscous stresses, and K accounts for the curvature that arises from a coordinate transformation. Our computations indicate that while the dissipative viscous term V is negative throughout the domain the production term P is sign-indefinite, S and K are negligible, and F is zero downstream of the forcing plane.
Insight into physical mechanisms can be gained by the analysis of dominant production terms in (3.17) associated with the global linearized response to upstream oblique disturbances. Averaging over the time period T = 2π/ω and the spanwise wavelength λ z = 2π/β, · := (T λ z ) −1 T 0 λz 0 (·) dz dt, and neglecting the terms that do not contribute significantly to the production of the averaged streamwise specific kinetic energy E s := E s yields the following approximation to the transport equation (3.17), where R sn := u s u n denotes the averaged shear stress of the streamwise velocity fluctuations. The second term on the left-hand-side represents the production of fluctuations' energy that arises from the streamwise gradient of the base flow ∂ sūs and the term on the right-hand-side determines the production term that originates from interactions of the base flow shear ∂ nūs with the fluctuation shear stress R sn .
To understand the mechanism that facilitate the growth of E s , we now investigate the streamwise transport of R sn . In contrast to the transport equation for E s , both the production P and curvature K terms contribute significantly to the streamwise transport of R sn for fluctuations with ω = 0.4 and λ z = 3.0. As demonstrated in appendix D, omitting negligible terms leads to the following approximate transport equation for R sn , where K c and K s denote contributions that arise from the curvature normal to the streamlines and from deceleration along the streamline direction, respectively, In the (s, n, z) coordinate system,Ω = ∂ xv − ∂ yū denotes the spanwise vorticity of the base flow in the Cartesian coordinates (Finnigan 1983) and using the definition of K s , equation (3.19a) simplifies toū In summary, equations (3.18) and (3.19c) determine a coupled system of linear equations that governs the streamwise transport of E s and R sn in the separation zone for oblique fluctuations with (ω = 0.4, λ z = 3.0), (3.20) Oblique waves experience largest amplification in the separated shear layer above the recirculation bubble, i.e., in the region where the presence of flow separation leads to concave flow curvature K c < 0. Figure 6(a) shows this negative curvature along the separation streamline and figure 6(b) illustrates the physical mechanism which is absent in attached boundary layers because of negligibly small positive streamwise curvature. Concave base flow curvature (i.e., K c < 0) in the shear layer provides destabilizing effect in system (3.20) which can be understood by simplifying equation (3.20) for oblique waves. In figure 7(a), we compareū s ∂ s E s :=ū s ∂E s /∂s with the dominant production term in (3.18) to illustrate that (∂ nūs )R sn dictates the streamwise growth of E s . Furthermore, since the base shear ∂ nūs remains almost constant throughout the shear layer (cf. figure 7(b)), its streamwise derivative can be neglected, thereby leading to, As shown in figure 8(b), the concave base flow curvature (i.e., K c < 0) provides the destabilizing influence throughout the separated shear layer and a simple mechanical analogy can be used to explain amplification of oblique waves. In the regions where K c < 0 the "spring constant" 4K c (∂ nūs ) in equation (3.22) is negative and this system behaves as an inverted pendulum, which enables the spatial growth of E s .
In summary, we have utilized the resolvent analysis to identify the spatial structure of oblique fluctuations that amplify rapidly in the separation zone. Furthermore, by conducting transport analysis of the most energetic fluctuations, we have demonstrated that the resulting amplification arises from concave curvature of the laminar 2D base flow.

Nonlinear interactions of oblique waves
In § 3, we used resolvent analysis to identify oblique waves as the most energetic responses of the linearized flow equations in the presence of unsteady disturbances. Recent numerical simulations (Lugrin et al. 2021) show that, even in the presence of unsteady disturbances, the dominant response near the reattachment appears in the form of streamwise streaks. To investigate the origin of steady streaks in the presence of unsteady external disturbances we utilize a weakly nonlinear formulation based on a perturbation expansion in the amplitude of the oblique disturbances. While previous numerical studies of transition induced by oblique waves in low-speed channel (Schmid & Henningson 1992) and boundary layer (Chang & Malik 1994;Mayer et al. 2011) flows show that nonlinear interactions of oblique waves generate streaks, we focus on the origin and spatial growth of these streaks in separated high-speed compressible boundary layer flows.

Streaks generated by oblique waves: a weakly nonlinear analysis
In the presence of a small external disturbance, − (x, y)e −2iωt e 2iβz . (4.2) Here,Ψ (x, y) is the 2D laminar base flow,d ± andψ (1) ± are the principal oblique input and state modes resulting from the linearized analysis of § 3, whereasψ Here, A 2β is the Fourier symbol of the dynamical generator in the linearized statespace model (3.8) andd (2) ± ) is the forcing term that arises from quadratic interactions of O( ) oblique waves with the spanwise wavenumber β; see appendix A for details. Thus, the resolvent operator associated with (3.8) evaluated at (ω = 0, 2β) maps the nonlinear modulationd (2) 0 ( · , · ) (x, y). (4.4) To investigate the emergence of streaks from unsteady disturbances, we introduce forcing inputs with (ω = ±0.4, λ z = 3) and examine a weakly nonlinear evolution of the resulting oblique waves. These forcing inputs are introduced at the upstream plane x in = 25 and their spatial structure is identified using the resolvent analysis of § 3 to generate the most energetic response at the reattachment (i.e., at x out = 60). Figure 9(a) illustrates the setup in which disturbances corresponding to a pair of oblique input modesd (1) ± with (ω = ±0.4, λ z = 3) are introduced at x in = 25. The resulting response of the linearized dynamics consists of oblique waves with opposite phase velocities, leading to a checkerboard wave pattern in the spanwise direction.  (2) 0 in that region. Figure 10(b) shows the wall-normal profiles of the forcing term to the mass, momentum, and temperature equations in (4.3) before reattachment at x = 58. We observe the strongest contribution of the forcing  to the wall-normal and spanwise components of the momentum equations, thereby demonstrating its vortical nature. Figure 10(c) illustrates the spatial structure of the forcing near reattachment in the (z, η) plane. The forcing termd (2) 0 which forms counterrotating vortices in the separated shear layer is 90 • out of phase relative to the induced streak response u s . In contrast to the dominant vortical forcing resulting from the linearized analysis, the vortical source term that arises from weakly nonlinear interactions of oblique waves primarily lies downstream of recirculation zone.
We utilize DNS to verify the predictions of weakly nonlinear analysis. In DNS, the input oblique modesd (1) ± (with ω = ±0.4, λ z = 3, and reference amplitude A 0 = 1.0) are introduced in the plane x in = 25. Additional details about the grid resolution and numerical implementation in our DNS study are provided in § 5. Figure 11 shows the spatial evolution of steady streamwise velocity fluctuations u s along the base flow separation streamline that are triggered by unsteady oblique disturbances of different amplitudes. As shown in figure 11(a), steady streaks undergo similar spatial growth across the range of forcing amplitudes. Figure 11(b) illustrates that the magnitude of streamwise velocity fluctuations collapses when scaled with A 2 0 . This demonstrates an excellent agreement between streak profiles resulting from DNS and a weakly nonlinear analysis (that leads to equation (4.3)). Deviations from predictions of the weakly nonlinear analysis are only observed for the largest amplitude considered and they are manifested by saturation of streaks in post-reattachment region.
Collapse of spatial profiles that characterize amplification of streaks irrespective of the amplitude of oblique disturbances, which differ by O(10 4 ), shows that streaks generated via quadratic interactions of oblique waves undergo linear amplification in the separation zone. This demonstrates predictive power of the weakly nonlinear analysis across the range of forcing amplitudes. In what follows, we utilize the input-output modes obtained from the resolvent analysis to characterize evolution of O( 2 ) steady streamwise streaks and uncover corresponding amplification mechanisms.

Resolvent mode representation of O( 2 ) streaks
As described in § 3.2, the left and right singular function of the frequency response operator provide orthonormal bases of input and output spaces that can be used to study responses of the double-wedge flow to external excitations. In particular, O( 2 ) steady streaks resulting from weakly nonlinear interactions of O( ) unsteady oblique waves (cf. equation (4.4)) can be represented using SVD of the frequency response operator associated with the linearized system (3.8) at ω = 0 and the spanwise wavenumber 2β, (4.5) Here, a i : ± ), quantifies the contribution of the ith output mode φ i of H 2β (0) to O( 2 ) steady streaksφ (2) 0 , σ i is the ith singular value, (d i , φ i ) are the corresponding input-output modes of H 2β (0), and for λ z = 1.5 the inner product ·, · E is carried over the entire flow domain in (x, y). Figure 12(a) shows 15 largest singular values of the resolvent for the linearized system (3.8) with (ω = 0, λ z = 1.5). Even though the principal singular value σ 1 is an order    Figure 13(a) shows the wall-normal profiles of the streamwise velocity component u associated withφ (2) 0 and the first two output modes (φ 1 , φ 2 ) of the resolvent. We observe striking similarity between O( 2 ) steady streaks and φ 2 in the post-reattachment region, at x = 65. Similarly, figure 13(b) compares the wall-normal shapes of the corresponding input modes d 1 and d 2 with the forcingd the second output mode of the resolvent associated with the linearized equations at (ω = 0, λ z = 1.5). To gain insight into the amplification mechanism that generates O( 2 ) streaks, we examine the dominant terms in the energy transport equation for the output mode φ 2 . Similar to the analysis in § 3.4, we utilize the (s, n, z) coordinate system which is locally aligned with the streamlines of the base flow (ū s , 0, 0). The transport of the spanwise averaged specific kinetic energy of streamwise velocity fluctuations E s = u s u s and fluctuation shear stress R sn = u s u n associated with the second output mode is approximately governed by equation (3.20) in § 3.4. Figure 14(a) shows the streamwise evolution of the wall-normal integral of E s for the mode φ 2 associated with O( 2 ) streaks with (ω = 0, λ z = 1.5). In contrast to the oblique waves (cf. figure 8(a)), which experience monotonic amplification throughout the separation shear layer above the recirculation bubble, we observe a non-monotonic xdependence of E s for the streaks within the bubble. To gain physical insight, we evaluate the terms of the coefficient matrix in equation (3.20) for steady streaks. Figure 14(b) shows that both streamwise deceleration (i.e., ∂ sūs < 0) and shear ∂ nū contribute to amplification of O( 2 ) streaks in different regions of the separated laminar flow. Towards reattachment, the base flow curvature K c is positive inside the recirculation bubble (cf. figure 14(c)) and the streamwise deceleration term (i.e., ∂ sūs < 0) is primarily responsible for the energy E s amplification. In this region, the transport equation for E s can be approximated by,ū This is in concert with the analysis of the contribution of the first (i.e., most amplified) output mode φ 1 to the amplification of streaks in a compression ramp flow (Dwivedi et al.

2019)
, which also revealed dominance of streamwise deceleration near reattachment. On the other hand, in the post-reattachment region, the shear term ∂ nūs < 0 dominates and the coupled system of equations (3.20) for E s and R sn simplifies to the following second order equation ( Thus, in the post reattachment region, the concave streamline curvature of the laminar 2D base flow (i.e. K c < 0) is primarily responsible for amplification of O( 2 ) streaks.

Direct numerical simulations of streak breakdown
Weakly nonlinear analysis demonstrates that small unsteady disturbances induce steady streaks in a hypersonic flow over a double-wedge. These streaks result from quadratic interactions of oblique waves and they undergo rapid amplification near reattachment. Downstream of reattachment, the streaks grow large enough to modify the time-averaged 2D flow and we utilize DNS to examine breakdown of streaks. We also report instantaneous and statistical properties of the flow as it transitions to turbulence.

Numerical setup
To study the onset of transition, we extend computational domain in the streamwise direction downstream of reattachment. We place the inflow boundary at x = 20 and, at this location, we interpolate the inflow profile that results from 2D base flow computations. To avoid spurious numerical errors we utilize a non-reflecting numerical sponge near the outflow. The wall is assumed to be adiabatic and periodic boundary conditions in the spanwise direction are applied.
Computational domain of size 80×13×9 in the streamwise, wall-normal, and spanwise directions is discretized using 900×249×384 grid points (i.e., 86 million cells). The grid is constructed to ensure uniform spacing in the streamwise and spanwise directions and, in the normal direction, the mesh is stretched to ensure y + < 0.22 at the wall. In appendix E, we provide a grid convergence study and compare numerical resolution that we use with recent simulations of supersonic and hypersonic transitional and turbulent flows. Our simulations utilize low-dissipation sixth-order spatially accurate kinetic-energy-consistent (KEC) fluxes (Subbareddy & Candler 2009) for spatial discretization. The low-dissipation KEC fluxes were previously employed in high-fidelity simulations of transitional and turbulent hypersonic boundary layers (Subbareddy et al. 2014;Subbareddy & Candler 2012). The time integration is carried out using the explicit third-order Runge-Kutta scheme with CFL number 0.5.

Secondary instability and breakdown
To simulate breakdown of streaks we excite the laminar 2D flow using the oblique input modes with (ω = ±0.4, λ z = 3). The disturbance amplitude is set to A ob = 2.50 × 10 2 A 0 and is determined based on the results reported in figure 11. With this amplitude, fluctuations grow linearly in the recirculation region but saturate nonlinearly post-reattachment (i.e., beyond x = 60).
In figure 15, we show time evolution of streaks in the plane close to the wall, at η = 0.25, for three different values of x. These plots demonstrate that, in the presence of unsteady oblique disturbances, O( 2 ) steady streaks undergo spanwise motion close to reattachment. The identified spanwise oscillations become stronger as we progress downstream and their period corresponds to the time period of oblique wave inputs. We note the presence of a "sinuous-subharmonic" motion, where two adjacent low-speed  streaks oscillate out of phase, and observe that the amplitude of spanwise oscillations almost doubles as we move from x = 55 to x = 65.
In figure 16, we illustrate the effect of streak oscillations on the mean flow by visualizing the root-mean-square of temporal fluctuations in the streamwise velocity. At different streamwise locations, fluctuations are plotted against ∂ū/∂z, which characterizes steady spanwise variations of the mean flow. Initially, unsteadiness is restricted to the oblique wave pair, which is located further away from the wall relative to streaks. However, at x = 65, unsteadiness is observed in the region of the largest spanwise shear because of the presence of streaks. This "locking-in" of u rms with streak oscillations identifies late stages of the streak evolution just before smaller-scales (associated with higher frequencies) set in.
Amplification of streaks and their unsteadiness induce rapid steepening of spanwise and wall-normal mean flow gradients, thereby leading to the emergence of inflection points in the mean (time-and spanwise-averaged) flow. Figure 17(a) shows the resulting inflectional mean flow profile and figure 17(b) shows the location of unsteady fluctuations with respect to the streaks. As we move from x = 65 to x = 70, we note the appearance of higher spanwise wavenumbers in the streaks as well as in the unsteady fluctuations. As discussed by Hall & Horseman (1991); Yu & Liu (1994), inflectional points in the mean velocity serve as an indicator of its susceptibility to the growth of broadband fluctuations. Amplification of high-frequency harmonics is also observed in the temporal spectra of the fluctuations. Therefore, at x = 70, flow is well within fully nonlinear stages of transition. We also note strong spatial correlation of unsteady fluctuations with the spanwise shear associated with the streaks, even at this nonlinear stage.
In figure 18, we plot Q-isosurfaces of the instantaneous flow-field. These visualize vortical structures that arise from sinuous-subharmonic motion of the streaks, prior to breakdown of the flow. We note the appearance of staggered lambda vortices, similar to the ones identified in transition induced by oblique waves in incompressible boundary layers ). These vortical structures are also sometimes referred to as X-vortices, for reasons illustrated in the (η, z) slice array plot in figure 18. Initially, at x ≈ 62, we observe that the interaction with the oblique waves causes the roll-up of the low-speed streak as they come together due to sinuous-subharmonic motion. Further downstream, as the flow structures associated with the X-vortices spread apart, we see that wall normal jets of low-speed flow form mushroom structures. As these lowspeed regions oscillate further, they interact to generate fluctuations with higher spanwise wavenumbers. Finally, at x ≈ 70, the shear introduced by the upward jets of low-speed streaks becomes significant enough to cause a local Kelvin-Helmholtz instability that induces streamwise rollers (Reddy et al. 1998). At this stage, the laminar boundary layer flow has started to break down to turbulence.

Transition to turbulence
Our DNS study demonstrates that amplification of steady streaks, that result from quadratic interactions of oblique waves, leads to a formation of a 3D inflectional boundary layer profile. During this process, fluctuations with multiple spatial and temporal scales develop and, despite periodic nature of upstream forcing, the flow downstream of reattachment becomes turbulent. In this section, we utilize wall and boundary layer statistics to illustrate the onset of turbulence on the second wedge.

Wall statistics: skin friction
Boundary layer transition is characterized by a rapid increase in wall friction. Figure 19(a) utilizes instantaneous wall-shear distribution to identify transitional and turbulent regions. The wall-shear stress experiences sinuous-subharmonic modulation downstream of the reattachment (at x = 60, i.e., Re x = 8.2 × 10 5 ) and finer spanwise scales emerge after x = 80 (i.e., Re x ≈ 11×10 5 ). As shown in figure 19(b), this location coincides with the highest value of time-averaged wall shear. Even though significant attenuation of spanwise variations of the wall-shear stress associated with the initial streaks takes place by x = 80, nonlinear interactions within the transition zone lead to the appearance of new streaks further downstream.
The spatial extent of transition zone is visualized in figure 19(c) by showing streamwise development of the skin friction coefficient, Here, τ * w is the dimensional shear stress at the wall, ρ * e and u * e are the dimensional density and streamwise velocity at the boundary layer edge, L z is the spanwise extent of the computational domain, and T = 4L x /u * ∞ . The values of C f in laminar flows over a 12 degree wedge and over the double-wedge are also plotted for comparison. The skin friction coefficient drops because of flow separation but it grows again after reattachment. Near and downstream of reattachment, we observe a significant difference between laminar and turbulent values of C f . After Re x = 11.5 × 10 5 , when C f starts to decay with x, the wallfriction coefficient is about seven times larger than its laminar counterpart. A comparison with the Van-Driest turbulent correlation for standard compressible boundary layers demonstrates that the flow on the second wedge is approaching fully developed turbulent values towards the end of the computational domain. In addition to skin friction, transition also has a significant impact on wall temperature. Additional details about the mean temperature and the wall statistics are included in appendix F.

Towards a turbulent compressible boundary layer
To illustrate the onset of turbulence, we also examine first-and second-order statistics in the latter stages of transition. At a given streamwise location, we report statistics in terms of the inner coordinate η + which is obtained by non-dimensionalizing the wall-normal coordinate with the local viscous length δ ν . Figure 20(a) shows the mean streamwise velocity u + non-dimensionalized by the local friction velocity u τ . The mean profile is obtained by averaging in time over 2L x /u ∞ and in the spanwise direction over the extent of the computational domain L z . In the presence of density variation along the compressible boundary layer, we utilize the Van-Driest transformation, to compare the transitional velocity profiles with the 'log-law' observed in incompressible turbulent boundary layers (White & Majdalani 2006), where · denotes averaging over time and spanwise direction, and ρ w is the mean density at the wall. A fully developed turbulent boundary layer is characterized by the presence of the viscous sublayer close to the wall (η + < 10), where u + c = η + , and further away from the wall the mean velocity obeys the log-law, where κ = 0.41 and C = 5.2. Figure 20(a) shows that the mean velocity in the boundary layer approaches the fully turbulent profile at Re x = 1.2 × 10 6 . Upstream of this location, within the transition zone, the boundary layer has a significantly greater momentum and most of it lies away from the wall. Furthermore, closer to the wall, the slope of the streamwise velocity profile is larger than the slope obtained using viscous sublayer profile of a fully turbulent flow. This observation is consistent with the overshoot in skin-friction coefficient within the transition zone; cf. figure 19(b).
In addition to the mean velocity, the mean temperature profile plays an important role in heat transfer and material response analysis of hypersonic flows. Walz's modified Crocco-Busemann relation, is commonly used to describe the relation between temperature and velocity in a zero pressure gradient turbulent boundary layer. Here, u e and T e denote the mean boundary layer edge velocity and temperature, respectively, T r is the mean recovery temperature and, since the wall is adiabatic, the mean wall temperature is determined by  We next examine spatial development of flow fluctuations by evaluating second-order statistics and comparing them with those observed in a Mach 5 fully developed turbulent flat plate boundary layer (Duan et al. 2011). Figures 21(a)-(c) show the streamwise variation of the density-weighted root-mean-square (RMS) values of the streamwise, wallnormal, and spanwise velocity fluctuation components. In the transition zone, we observe large values of fluctuations away from the wall in all three plots. Further downstream, at Re x = 1.2 × 10 6 , the RMS values closer to the wall are in agreement with fully turbulent values. However, away from the wall, the RMS values of u and w deviate from the flat plate profiles. We attribute this deviation to the presence of 3D oblique waves and streaks that persist downstream because of continuous upstream forcing in our setup.
To illustrate the broadband nature of velocity fluctuations in the later stages of transition, in figure 21(d) we evaluate the spanwise wavenumber dependence of the individual contributions of velocity fluctuations to the 1D energy spectrum at Re x = 1.2 × 10 6 and η + = 40. At this location, the energy spectrum of streamwise velocity fluctuations clearly exhibits an inertial and a dissipative subrange features, indicating that the flow is approaching a fully turbulent stage.
Our DNS also provides data for evaluating non-dimensional parameters which can be utilized for low-complexity modeling of high-speed compressible turbulent flows using time-averaged NS equations (Hirsch 2007). In particular, we examine the spatial evolution of fluctuation Mach number, M t , fluctuation Prandtl number, P r t , and Huang's modified strong Reynolds analogy parameter (Huang et al. 1995), ψ sra , Here, u i u i = u u + v v + w w , γ and R denote the specific heat ratio and the gas constant, respectively, and T 0 is the stagnation temperature. In a fully developed turbulent boundary layer, ψ sra provides a measure of correlation between velocity and temperature fluctuations (Huang et al. 1995). Figure 22 shows that towards the end of the computational domain, at Re x = 1.2×10 6 , these parameters are close to those observed in canonical hypersonic boundary layers, i.e., M t ≈ 0.2 − 0.3, P r t ≈ 0.9, and ψ sra ≈ 1 (Pirozzoli & Bernardini 2011). However, upstream of the breakdown region, at Re x = 1.0 × 10 6 , there is a significant deviation compared to these canonical values. Here, M t becomes as high as 0.45 which suggests that from canonical values in the outer region of boundary layer are also observed in ψ sra . We conjecture that, in the presence of persistent upstream excitations, the resulting unsteady fluctuations are primarily responsible for these discrepancies.

Concluding remarks
Axisymmetric cone-flare experiments (Benitez et al. 2020;Butler & Laurence 2021) identified unsteady fluctuations in the separation zone and it is believed that these play an important role in initiating hypersonic flow transition. As demonstrated in figure 23, we observe strong qualitative similarity between spatial structures of fluctuations observed using schlieren measurements in Mach 6 reattaching flow on axisymmetric coneflare (Butler & Laurence 2021) and the dominant oblique density fluctuations that we identify via resolvent analysis of separated flow over a slender double-wedge. Inspired by these observations, we have examined transition mechanisms in a Mach 5 hypersonic flow over a slender double-wedge subject to unsteady upstream disturbances.
To investigate the early stages of transition, we employ resolvent analysis to evaluate responses of the laminar 2D base flow to exogenous time-periodic inputs. This allows us to identify prevailing spatio-temporal scales, the spatial structure of disturbances that most effectively excite the double-wedge flow, as well as the spatial structure of the resulting responses and the underlying physical mechanisms. In the presence of flow separation, our analysis shows that two types of disturbances are strongly amplified by the linearized dynamics: steady streamwise vortices and unsteady oblique waves. While amplification of steady upstream vortical disturbances has been studied in Dwivedi et al. (2019), oblique waves that amplify within separation-reattachment zone have not been investigated. Recently, Lugrin et al. (2021) examined the growth of unsteady perturbations that arise from the "first mode" instability of an axisymmetric boundary layer over cylinder-flare geometry. In the presence of stochastic disturbances in the inlet of the computational domain, DNS was used to demonstrate that oblique waves that emerge in the upstream boundary layer (i.e., before separation) can undergo nonlinear interactions similar to those observed in attached compressible boundary layers (Chang & Malik 1994;Mayer et al. 2011) and cause transition in separated high-speed flows (Lugrin et al. 2021). In contrast, we show that unsteady disturbances that are localized upstream of the corner trigger oblique waves downstream of the corner even in the absence of local or global boundary layer instabilities. These oblique waves experience significant amplification within separation-reattachment zone and their role in initiating transition in flows with SWBLI has not been studied before. By carrying out resolvent analysis of the linearized flow equations subject to disturbances that are introduced in a plane immediately upstream of separation, we identify the physical mechanism responsible for non-modal amplification of oblique waves in presence of flow separation. The subsequent nonlinear interaction of dominant unsteady oblique fluctuations is examined using a weakly nonlinear analysis to demonstrate the emergence of steady reattachment streaks inside the recirculation bubble. DNS confirms the predictions of our analysis and provides a detailed characterization of transition initiated by time-periodic upstream oblique disturbances.
We next briefly summarize our main contributions: (i) Amplification of oblique waves by base flow curvature. We analyze fluctuations' kinetic energy in a streamline-aligned orthogonal coordinate system to identify physical mechanisms responsible for amplification of oblique waves in the separation zone. Large energy amplification arises from the growth of the fluctuation shear stress due to streamline curvature in the separated shear layer. This is in contrast to the attached boundary layers, where no such mechanism exists. To compare separated and attached boundary layers, we also conduct resolvent analysis of the flow over a wedge that does not contain the compression corner (this wedge is identical to the first wedge in the doublewedge configuration analyzed in the paper). Figure 24 demonstrates that the presence of a recirculation zone in the double-wedge flow significantly increases amplification relative to the single-wedge flow. The amplification profiles of the fluctuation shear stresses differ in these two cases because of fundamentally different physical mechanisms.
(ii) Steady reattachment streaks through base flow deceleration. We utilize a weakly nonlinear analysis to show that the resolvent operator associated with the linearized dynamics governs the evolution of steady streaks that arise from quadratic interactions of unsteady oblique waves. Vortical excitations in the reattaching shear layer generated by these interactions lead to the formation of streaks in the recirculation zone and their subsequent amplification downstream. Additionally, we use SVD of the resolvent operator to demonstrate that secondary reattachment streaks are well approximated by the second output resolvent mode. Similar to the most amplified steady output in compression ramp flow (Dwivedi et al. 2019), our analysis of the energy budget shows that deceleration of the laminar base flow near reattachment is responsible for amplification of reattachment streaks associated with this sub-optimal mode.
(iii) Transition to turbulence. We use DNS to examine nonlinear stages of the evolution of flow fluctuations. In the presence of strong upstream oblique disturbances, steady streaks saturate after reattachment and experience sinuous sub-harmonic oscillations. The resulting 3D boundary layer breaks down further downstream and the observed flow structures are similar to those in other canonical configurations: nonlinear interactions of streaks with oblique waves lead to the development of staggered patterns of lambda vortices, followed by a rapid emergence of higher harmonics in fluctuations and multiple inflection points in the mean velocity profile before breakdown to turbulence (e.g., see Hall & Horseman 1991;Yu & Liu 1994;Reddy et al. 1998). As the flow transitions to turbulence, the wall friction increases rapidly before settling to the values predicted by turbulent correlations. Within the transitional zone, non-dimensional parameters that are critical for modeling temperature and compressibility effects in high-speed turbulent flows (Hirsch 2007), e.g., the fluctuation Prandtl and Mach numbers, significantly differ from their fully developed turbulent values. Post-breakdown, the boundary layer develops mean and fluctuation statistics that agree well with observations made in attached hypersonic turbulent boundary layers, thereby demonstrating the efficacy of unsteady oblique waves in triggering transition in separated high-speed boundary layer flows.
Unsteady disturbances in hypersonic boundary layers can arise from free-stream turbulence in wind tunnel experiments (Schneider 2001(Schneider , 2015, interactions of unsteady free-stream disturbances with surface roughness (Wu 2001;Gonzalez & Wu 2019), atmospheric particulates associated with ice-clouds and volcanic dust (Turco 1992;Chuvakhov et al. 2019), and atmospheric turbulence (Bushnell 1990). Novel physical mechanisms that we identify are unique to high-speed boundary layers with a separation-reattachment zone. We expect that insights about transition mechanisms that we provide using a combination of resolvent and weakly nonlinear analyses with DNS will motivate a systematic evaluation of specific disturbance environments that appear in wind tunnels or free flights and guide the development of low-complexity models for fast and accurate prediction of transition in hypersonic flows under realistic in-flight conditions. As shown in § 4.1,d 0,ρv = (γ − 1) 0,ρw = iβρ (γ − 1) 0,Et = −iβρ w (1) + (1) where we drop the caret notation on the right-hand side of equation (A 1) for brevity.
Here, C v = R/(γ − 1), where R is the gas constant, γ = 1.4 is the ratio of the specific heat capacities, and µ(x, y) and λ(x, y) are the coefficients of viscosity and bulk viscosity, respectively. We utilize Sutherland's law for computing viscosity and assume λ = −2µ/3.

Appendix B. Relation between the state Ψ and the output Φ
Herein, we utilize a weakly nonlinear expansion to establish the relation between the components of the state vector Ψ = (Ψ 1 , Ψ 2 , Ψ 3 ) := (ρ, ρu, E t ) in conserved variables and the components of the vector Φ = (Φ 1 , Φ 2 , Φ 3 ) := (ρ, u, T ) in primitive variables. Here, E t is the total energy per unit volume of the gas, C v is the specific heat at constant volume, and |u| 2 := u T u. Within the weakly nonlinear framework, we can decompose Ψ and Φ into the sums of base and fluctuating components, and utilize the following relations between the components of Ψ and Φ, Duan (1) Mayer (2) Pirozzoli (3) Franko (4)

Appendix F. Distribution of wall temperature
In addition to the skin friction, the distribution of wall temperature provides insights into the thermal effects encountered in compressible boundary layer flows along the transition zone. In flows with adiabatic walls, there is no heat transfer to the wall and viscous dissipation near the wall converts kinetic into internal energy. This leads to high temperatures near the wall as well as in the associated thermal boundary layer and 3D patterns in the wall temperature are caused by temperature transport within the boundary layer by flow fluctuations.
Figures 26(a) and 26(b) illustrate instantaneous and mean wall temperatures in the transition zone. In contrast to the instantaneous skin friction, where the streaks determine the spanwise modulation prior to flow transition, the instantaneous wall temperature contains strong imprints of unsteady oblique waves. The role of the streaks and higher spanwise harmonics becomes apparent when we examine the time-averaged wall temperature.
In Figure 26(c), we illustrate the mean wall-temperature along the double-wedge. Even though temperature variations are not significant under present conditions, its analysis can be informative for different free-stream conditions. Comparison with the double wedge laminar solution shows that the wall temperature is higher in the 3D flow field immediately after the separation point and that it rises rapidly post-reattachment. In the transition zone, we observe an overshoot before reduction to its turbulent value.