arxiv: v1 [physics.flu-dyn] 22 Nov 2018

Transcription

1 This draft was prepared using the LaTeX style file belonging to the Journal of Fluid Mechanics 1 arxiv: v1 [physics.flu-dyn] 22 Nov 2018 Reattachment vortices in hypersonic compression ramp flow: an input-output analysis Anubhav Dwivedi 1, G. S. Sidharth 1, Joseph W. Nichols 1, Graham V. Candler 1, and Mihailo R. Jovanović 2 1 Department of Aerospace Engineering and Mechanics, University of Minnesota, Minneapolis, MN 55455, USA 2 Ming Hsieh Department of Electrical Engineering, University of Southern California, Los Angeles, CA 90089, USA (Received xx; revised xx; accepted xx) We employ global input-output analysis to quantify amplification of exogenous disturbances in compressible boundary layer flows. Using the spatial structure of the dominant response to time-periodic inputs, we explain the origin of steady reattachment and Görtler-like vortices in laminar hypersonic flow over a compression ramp. Our analysis of the resulting shock/boundary layer interaction reveals that the vortices arise from selective amplification of upstream perturbations. These vortices cause heat streaks with a specific spanwise wavelength near flow reattachment and they can trigger transition to turbulence. The streak wavelength predicted by our analysis compares favorably with observations from two different hypersonic compression ramp experiments. 1. Introduction Compression corners are commonly encountered in intakes, control surfaces, and junctions. High speed flow on a compression corner is a canonical case of shock/boundary layer interaction (SBLI) (Simeonides & Haase 1995) involving flow separation and reattachment with a shock system. Even though the compression ramp geometry is homogeneous in the spanwise direction, experiments (Roghelia et al. 2017) and numerical simulations (Navarro-Martinez & Tutty 2005) show that the flow over it exhibits threedimensionality in the form of streamwise streaks near reattachment. The streaks are associated with persistent large local peaks of heat transfer; they can destabilize the boundary layer and cause transition (Simeonides & Haase 1995; Roghelia et al. 2017). Recently, Roghelia et al. (2017) and Chuvakhov et al. (2017) investigated hypersonic compression ramp flows using temperature sensitive paint (TSP) and infrared (IR) imaging. These techniques were employed to study the formation of streamwise vortices and reattachment heat flux patterns. The observed structures are often attributed to Görtler-like vortices (Inger 1977; Simeonides & Haase 1995; Chuvakhov et al. 2017). These vortices develop from upstream perturbations (Hall 1983) which are amplified by centrifugal forces that arise from concave streamline curvature near reattachment (Dwivedi et al. 2017). However, in most compression ramp studies, the comparison with theory of Görtler instability on curved walls is qualitative and does not account for the dynamics in the separation bubble. Zhuang et al. (2017) used nano-tracer planar laser scattering to visualize a Mach 3 turbulent boundary layer turning on a 25 compression ramp. They found that Görtler-like vortices not only appear after reattachment but also address for correspondence: dwive016@umn.edu

2 2 A. Dwivedi, G. S. Sidharth, J. W. Nichols, G. V. Candler and M. R. Jovanović Figure 1. Schematic of the input-output approach to compressible flow instabilities. in the separation bubble. It is thus important to understand the effect of the recirculation bubble dynamics on the formation and amplification of streamwise vortices. To include the effect of the separated flow, Sidharth et al. (2018) carried out a global stability analysis and discovered a 3D global instability in the separation bubble, which results in temperature streaks post-reattachment. The spanwise wavelength of the global instability scales with the recirculation length (Sidharth et al. 2017). This is in contrast to the spanwise wavelength observed for Görtler-like vortices (Chuvakhov et al. 2017; Roghelia et al. 2017; Navarro-Martinez & Tutty 2005), which scale with the reattachment boundary layer thickness, indicating that the global instability is not responsible for their formation. To characterize the role of external perturbations in the formation of these vortices, we consider compression ramp flows that do not exhibit 3D global instability. External perturbations are amplified as they pass through the flow field and we utilize global input-output (I/O) analysis to quantify this amplification. The I/O analysis evaluates the response (outputs) of a dynamical system to external perturbation sources (inputs). For time-periodic inputs, the transfer function maps the input forcing to output responses; see Figure 1 for an illustration. For small perturbations, the transfer function can be obtained by linearizing the compressible Navier-Stokes (NS) equations around a laminar base flow. The I/O approach has previously been employed to quantify amplification and study bypass transition mechanisms in incompressible channel flows (Jovanović 2004; Jovanović & Bamieh 2005), boundary layer flows (Brandt et al. 2011; Fosas de Pando & Schmid 2017; Cook et al. 2018; Ran et al. 2018), and jets (Jeun et al. 2016). In this paper, the input-output pair representing the dominant response of the linearized dynamics is computed to quantify the amplification of upstream perturbations in shock/boundary layer interaction. Our presentation is organized as follows. In 2, we present the linearized model and provide a brief summary of the I/O formulation. In 3, we compute the amplification in attached supersonic flat-plate boundary layers and verify our method against stateof-the-art approaches. In 4, we evaluate the dominant input-output 3D perturbation fields around a 2D laminar hypersonic base flow on a compression ramp to explain the origin of streamwise vortices near reattachment. We conclude our presentation in Formulation and background The compressible NS equations for perfect gas in conservative form are given by U t + F j x j = 0, (2.1) where F j (U) is the flux vector and U = ( ρ, ρu, E ) is the vector of conserved variables representing mass, momentum, and total energy per unit volume of the gas (Candler

3 Reattachment vortices in hypersonic compression ramp flow 3 et al. 2015). We decompose the state vector U(x, t) into a steady base component U(x) and a perturbation component U (x, t), U(x, t) = U(x)+U (x, t). The evolution of small perturbations is then governed by the linearized flow equations, t U (x, t) = A(U)U (x, t), (2.2) where A(U) represents the compressible NS operator resulting from linearization of (2.1) around the base flow U. A second order central finite volume discretization (Sidharth et al. 2018) is used to obtain the finite dimensional approximation of Eq. (2.2), d q = Aq, (2.3) dt which describes the dynamics of spatially discretized perturbation vector q. In the present work, we are interested in quantifying the amplification of exogenous disturbances in boundary layer flows (Schmid 2007). To accomplish this objective, we augment the evolution model (2.3) with external excitation sources d q = Aq + Bd, dt φ = Cq, (2.4) where d is a spatially distributed and temporally varying disturbance source (input) and φ is the quantity of interest (output). In Eq. (2.4) the matrix B specifies how the input enters into the state equation, while the matrix C extracts the output from the state q. An I/O relation is obtained by applying the Laplace transform to (2.4), φ(s) = C(sI A) 1 (q(0) + Bd(s)), (2.5) where q(0) denotes the initial condition and s is the complex number. Equation (2.5) can be used to characterize both the unforced (to initial condition) and forced (to external disturbances) responses of the flow perturbations. In boundary layer flows, the linearized flow equations are typically stable. Thus, for a time-periodic input with frequency ω, d = ˆde iωt, the steady-state output of a stable system (2.4) is given by φ = ˆφe iωt, where ˆφ = Hˆd and H is the frequency response H(ω) = C(iωI A) 1 B. (2.6) At any ω, the singular value decomposition of H(ω) can be used to quantify amplification of time-periodic inputs (Jovanović 2004; Schmid 2007; McKeon & Sharma 2010), HD = ΦΣ H = ΦΣD. (2.7) Here, ( ) denotes the complex-conjugate transpose, Φ and D are unitary matrices, and Σ is the rectangular diagonal matrix of the singular values σ i. The columns d i of the matrix D represent the input forcing directions that are mapped through the frequency response H to the corresponding columns φ i of the matrix Φ; for ˆd = d i, the output ˆφ is in the direction φ i and the amplification is determined by the corresponding singular value σ i. For a given temporal frequency ω, we use a matrix-free approach (Dwivedi et al. 2018) to compute the largest singular value σ 1 of H. Note that G = σ 1 = φ 1 E / d 1 E denotes the largest induced gain with respect to a compressible energy norm (Hanifi et al. 1996), where (d 1, φ 1 ) identify the spatial structure of the dominant I/O pair.

4 4 A. Dwivedi, G. S. Sidharth, J. W. Nichols, G. V. Candler and M. R. Jovanović Figure 2. Results of I/O analysis on (a) Mach 4.5 adiabatic boundary layer with perturbation input at x = 0.1m; (b) pressure perturbations corresponding to the principal output mode φ 1 ; and (c) comparison of spatial growth rate with LST (Ma & Zhong 2003). 3. Validation: supersonic flat plate boundary layer Before analyzing the amplification of disturbances in a hypersonic flow involving shock/boundary layer interaction, we apply I/O analysis to compute amplification in a supersonic flow over a flat plate. Our computations are verified against conventional approaches to demonstrate the agreement for canonical problems. Two amplification mechanisms are considered: two-dimensional unsteady acoustic amplification (Ma & Zhong 2003) and three-dimensional steady lift-up amplification (Zuccher et al. 2005) Two-dimensional unsteady perturbations: acoustic amplification Local spatial instabilities corresponding to acoustic perturbations dominate the transition in high speed flat plate boundary layers (Fedorov 2011). Using local spatial linear stability theory (LST) and direct numerical simulations (DNS), Ma & Zhong (2003) showed that perturbation frequencies 50 KHz < f < 183 KHz result in spatial growth due to the local instability over a part of the domain. We consider I/O analysis at f = 133 KHz and compare the region of growth with that predicted from LST (Ma & Zhong 2003); see Figure 2(a) for geometry. The base flow is computed using the finite volume compressible flow solver US3D (Candler et al. 2015) with 125 cells in the wall-normal and 1600 cells in the streamwise direction. This resolution yields grid-insensitive I/O results. As shown in Figure 2(a), we use the matrix B in Eq. (2.4) to localize the disturbance input at x 0 = 0.1 m downstream of the leading edge. This choice allows us to avoid large streamwise gradients in the base flow in the vicinity of the leading edge. The slow streamwise variation of the base flow implies that LST is approximately valid downstream of this location. Furthermore, this location is sufficiently upstream of the neutral point of the acoustic instability (x = 0.18 m). This ensures that any non-modal growth arising from the Orr-mechanism (Dwivedi et al. 2018) decays before the spatial growth rate of the local acoustic instability becomes positive. Sponge regions are used at the top and right boundaries to model non-reflecting radiation boundary conditions. We have verified independence of our results on the strength and the location of the sponge zones. The output of interest is chosen to be the perturbation field in the entire domain, i.e., φ = q. Figure 2(b) shows the spatial structure of pressure perturbation in the principal output mode φ 1. We compute the local spatial growth rate from pressure at the wall ˆp wall, α i = ( ˆp wall / x)/ˆp wall. Figure 2(c) shows that our I/O analysis correctly identifies the region of spatial instability and predicts growth rates that are close to those resulting from LST (Ma & Zhong 2003). The difference can be attributed to the fact that LST does not account for the spatially-growing nature of the base flow Three-dimensional steady perturbations: lift-up mechanism The spatially-developing boundary layer also supports significant growth of perturbations that are not related to a dominant eigenmode of the linearized dynamical generator. For example, the steady 3D streak-like perturbations, that result from the liftup mechanism (Ellingsen & Palm 1975), play an important role in transition induced by

5 Reattachment vortices in hypersonic compression ramp flow 5 Figure 3. (a) Optimal I/O amplification for steady perturbations; (b) streamwise velocity perturbation for different spanwise wavenumbers β along the base-flow streamlines; and (c) contours of streamwise vorticity (d 1) and isosurfaces of streamwise velocity perturbation (φ 1 ) along the B/L (shown using base flow streamwise velocity). The dashed curve in (a) indicates the spanwise wavenumber from spatial transient growth calculations (Zuccher et al. 2005).. distributed surface roughness (Reshotko 2001). Zuccher et al. (2005) used the linearized boundary layer (BL) equations to compute spatial transient growth and analyze this mechanism. For the verification purposes, we compare the spanwise wavelength of the maximally amplified streaks resulting from the linearized BL equations and the I/O analysis. We specifically consider the conditions in Zuccher et al. (2005) corresponding to a boundary layer on a 2D flat plate in a supersonic free-stream. A grid with 250 cells in the wall-normal and 600 cells in the streamwise direction is used to compute the base flow and conduct I/O analysis. The input is localized to a plane at streamwise location x/l = 0.3, where L denotes the plate length. As in the previous subsection, the output is the perturbation field in the entire domain. Due to homogeneity in the spanwise direction, 3D perturbations are of the form q = q(x, y)e i(βz ωt), where β = 2π/λ z is the spanwise wavenumber. To capture the steady lift-up mechanism, we conduct the I/O analysis for ω = 0. In Figure 3(a), we illustrate the β-dependence of the gain G resulting from the I/O analysis and the value of β at which the largest spatial transient growth takes place (Zuccher et al. 2005). We attribute a slight difference in the location of maxima to different base flow profiles; while we use a numerically computed 2D base flow, Zuccher et al. (2005) used an analytical self-similar base flow profile. The input d 1 (shown in Figure 3(c)) consists of streamwise vortical perturbations and the output φ 1 consists of a rapid development of streamwise velocity streaks. The algebraic nature of the growth (as opposed to exponential) is illustrated in Figure 3(b) for different β. As expected, a large initial transient growth is followed by eventual downstream decay. The above results show that I/O analysis correctly captures the physical mechanisms responsible for amplification in canonical supersonic flows. As we demonstrate in the next section, this analysis can also be readily applied to study the early stages of transition in complex hypersonic compression ramp flow with shock/boundary layer interaction. 4. Input-output analysis of hypersonic compression ramp flow Streamwise streaks in wall temperature are often observed in compression ramp experiments. Although they are typically attributed to the centrifugal amplification that arises from concave flow curvature near reattachment (Navarro-Martinez & Tutty 2005; Chuvakhov et al. 2017), quantifying amplification in the presence of a recirculation bubble is an open challenge. Herein, we employ the I/O framework to study the amplification of spanwise periodic upstream disturbances in the compression ramp flow and explain the origin of heat streaks at reattachment. Recently, Roghelia et al. (2017) and Chuvakhov et al. (2017) reported multiple hypersonic compression ramp experiments in two different facilities with matched freestream

6 6 A. Dwivedi, G. S. Sidharth, J. W. Nichols, G. V. Candler and M. R. Jovanović Re L p T U ρ Pa 55 K 1190 m/s Kg/m Pa 55 K 1188 m/s Kg/m 3 Table 1. Free-stream conditions for experiments reported in Chuvakhov et al. (2017) and Roghelia et al. (2017), respectively. Reynolds number Re L is based on plate length L. Mach and Reynolds numbers. Temperature sensitive paint (TSP) and infrared thermography measurements of reattachment heat-flux wall patterns revealed quantitatively similar streaks. The effects of freestream Reynolds number and leading edge radius on the spanwise wavelength λ z of the streaks were also reported. In Inger (1977); Navarro-Martinez & Tutty (2005), scaling of λ z with the incoming boundary layer was analyzed. Our objective is to identify the streak wavelength λ z that is selected by the shock/boundary layer interaction. We consider the experiments performed in the UT-1M Ludwig tube (Chuvakhov et al. 2017) at Mach 8. The geometry consists of an L = 50 mm isothermal flat plate with a sharp leading edge and wall temperature T w = 293 K, followed by an inclined ramp with 15 ; see Figure 4(a). Table 1 summarizes the two free-stream conditions that we considered in our study. Figure 4(b) provides comparison of the experimental schlieren image with the 2D base-flow density gradient magnitude field that is computed using US3D (Candler et al. 2015). Apart from the fact that the experimental schlieren contains the side-wall shock, the overall agreement is excellent. Figure 4(c) uses variation of the Stanton number St along the surface (Schlichting & Gersten 2016) to compare the heat fluxes resulting from experiment and simulations with different grid resolutions. Except near separation and post-reattachment regions, the 2D flow captures the heat flux trends correctly. The observed discrepancy can be attributed to the presence of 3D flow structures in the experiments. Since the flow is globally stable with respect to 3D perturbations (Sidharth et al. 2018), we conjecture that these flow structures arise from non-modal amplification of 3D perturbations around the 2D base flow. To verify our hypothesis, we employ global I/O analysis to quantify the amplification of exogenous disturbances and understand mechanisms that trigger transition in a hypersonic compression ramp flow. Because of the persistence of heat streaks in the experimental TSP images, we confine our analysis to amplification of inputs with zero temporal frequency (i.e., we set ω = 0). Our I/O analysis of steady disturbances reveals the importance of incoming streamwise vortical perturbations as inputs which lead to reattachment streaks as the response. This observation holds for all values of λ z that we consider in our study. We also note that responses to unsteady low-frequency inputs (with ω 0) are nearly identical to the steady case (Dwivedi et al. 2018) and are omited for brevity. The I/O analysis is conducted on a grid with 412 cells in the streamwise and 249 cells in the wall-normal direction (labeled as G3 in Figure 4(c)). Via the proper selection of the matrix B in Eq. (2.4), we restrict the inputs to the domain prior to separation (i.e., x/l < 0.5) and choose the perturbation field in the entire domain as the output, φ = q. Numerical sponge boundary conditions are applied near the leading edge (x/l < 0.02) and the outflow (x/l > 1.6). Figure 5(a) shows the dependence of the gain G on λ z. For both Reynolds numbers, the amplification curve achieves its maximum for a particular value of λ z. This indicates that SBLI selectively amplifies upstream perturbations with a specific spanwise wavelength. For both free-stream conditions, the experimental estimates of λ z resulting from the TSP

7 Reattachment vortices in hypersonic compression ramp flow 7 Figure 4. (a) Flow geometry and 2D steady streamwise velocity at Re L = ; (b) comparison to experimental schlieren; and (c) dependence of Stanton number (St) on x/l. Curves G1-G4 denote computational grids at varying resolution, with G1 ( ), G2 ( ), G3 ( ), and G4 ( ); the first entry denotes number of streamwise and second wall normal cells in the grid. images in Figure 5(b) agree well with the predictions of I/O analysis. Both approaches identify λ z /δ R 1.8 for maximally amplified streaks, where δ R is the reattachment boundary layer thickness. This value is also consistent with previous studies (Inger 1977; Navarro-Martinez & Tutty 2005). We finally note that, in the presence of random upstream perturbations, λ z corresponding to the peak in the amplification curve is predicted to be the dominant wavelength of the reattachment streaks (Nichols 2018). To gain insight into the structure of the most amplified flow perturbations, we examine the spatial evolution of the streamwise velocity and helicity perturbations for three different values of λ z = 1, 3, 10 mm. These respectively identify the flow structures with small, dominant, and large spanwise wavelengths. Figure 6(a) shows the streamwise velocity perturbation magnitude u for different values of λ z along the 2D steady flow streamline through the boundary layer edge. After initial amplification post-separation (S), the amplification curve exhibits two distinct profiles in the separated shear layer. For small wavelengths, λ z /δ R 1, u saturates in the separated shear layer, persists over the length of the bubble, and undergoes rapid amplification near the reattachment (R). In contrast, for λ z = 10 mm, we observe a steady amplification in the bubble followed by a weak amplification near reattachment. For the wavelength that corresponds to the largest amplification, λ z /δ R = 1.8, the large λ z trends are followed until reattachment and significant amplification is achieved in the separated flow. Near reattachment, u starts to follow the profile corresponding to flow perturbations with small λ z. These amplification trend are further illustrated in Figure 6(b) by plotting the corresponding perturbation helicity, u ω, for different values of λ z. Here, u is the base flow velocity vector, ω is the perturbation vorticity vector, and the perturbation helicity determines the projection of ω on the base flow streamline. We note that a contribution from the separation bubble increases with increase in λ z. Thus, the selected wavelength represents incoming vortical perturbations that grow post-separation and feed on the centrifugal effects at reattachment. In what follows, we refer to these spatially amplifying perturbations as reattachment vortices. We now demonstrate that this post-separation amplification in the reattachment vortices is associated with the growth of perturbations in the recirculation bubble. We compare and contrast this mechanism with Görtler vortices where the growth exclusively arises from centrifugal effects. To quantify the difference, we repeat the I/O analysis in the absence of perturbation development in the recirculation bubble for higher value of the Reynolds number. By introducing a numerical sponge zone within the bubble,

8 8 A. Dwivedi, G. S. Sidharth, J. W. Nichols, G. V. Candler and M. R. Jovanović Figure 5. (a) Amplification map for steady inputs with ω = 0; and (b) comparison of experiments and dominant output mode at reattachment. Dashed lines denote low and solid lines denote high Reynolds number cases. The red streaks denote regions of increased heating. An inset of T isosurface corresponding to φ 1 is shown in (a). we ensure that the perturbation growth is dominated by centrifugal effects. In this case, the I/O analysis reveals streamwise Görtler vortices as dominant outputs and the amplification map G(λ z ) is shown in Figure 7(a). We see that eliminating the role of the bubble causes the spanwise wavelength that corresponds to the largest amplification to decrease from λ z /δ R = 1.8 to λ z /δ R = 1.2. The difference between the Görtler and the reattachment vortices is illustrated by plotting the corresponding perturbation helicity. While the amplification in the separation bubble for the reattachment vortices is caused by streamwise vortical motions in the separation zone, Görtler vortices only consist of counter-rotating streamwise vortices that travel parallel to a surface (akin to an artificial wall with the curvature properties of the separation streamline). Figure 7(b) shows the evolution with streamwise coordinate of u and the Görtler number G = (U l/ν) l/r close to the separation streamline, where, l and R respectively denote the local viscous length scale and radius of curvature. The spatial growth of Görtler vortices is compared with that of reattachment vortices. Both experience strong amplification in regions of high G but the difference is attributed to the amplification of reattachment vortices in the separated shear layer. We note that, for small values of λ z, the I/O analysis of the shock/boundary layer interaction reveals Görtler vortices as the dominant response. Their negligible spatial growth in the separated shear layer was observed in Figure 6(b) in the presence of perturbation dynamics in the bubble. This is consistent with Figure 7(a), which shows almost identical gains for small λ z. In summary, we have quantified the amplification of 3D upstream perturbations in a 2D base flow with shock/boundary layer interaction. The I/O analysis demonstrates that vortical perturbations with a specific spanwise wavelength are most amplified. In addition to centrifugal effects near reattachment, perturbation evolution in the recirculation bubble also contributes to their amplification. 5. Concluding remarks We have employed an input-output analysis to investigate amplification of disturbances in compressible boundary layer flows. Our approach utilizes global linearized dynamics to propagate flow perturbations and analyze the spatial structure of the dominant response. We have verified that I/O approach captures both acoustic and vortical spatial growth mechanisms on a supersonic flat plate boundary layer without any a priori knowledge of the perturbation form. In an effort to explain steady heat streaks near reattachment, we have also examined the origin of Görtler-like vortices in shock/boundary layer interaction on Mach 8 flow

9 Reattachment vortices in hypersonic compression ramp flow 9 Figure 6. (a) Development of streamwise velocity perturbation magnitude for different spanwise wavelengths λ z along the separation streamline. Perturbation vorticity sketch in the (y, z)-plane at x/l = 0.75 for large and small values of λ z are also shown. Separation streamline is abbreviated as sep.sl. (b) Perturbation helicity in the (x, y)-plane shows larger contribution from separation bubble for larger λ z. The bold black line denotes the separation streamline. Figure 7. (a) Amplification map for steady inputs plotted with insets of perturbation helicity in the presence and absence of perturbation dynamics in the recirculation bubble; and (b) development of streamwise velocity perturbation for optimal λ z plotted with Görtler number G. over 15 compression ramp. In spite of global stability, the I/O analysis predicts large amplification of incoming streamwise vortical perturbations with a specific spanwise lengthscale. We have demonstrated that, in addition to centrifugal effects, perturbations in the separation bubble significantly contribute to amplification of streamwise reattachment vortices and that our predictions agree well with two recent experimental observations. The I/O approach provides a useful framework to quantify the spatial evolution of external perturbations in shock/boundary layer interactions. Improved understanding of amplification mechanisms can provide important insight and pave the way for the development of predictive transition models in complex high speed flows. Acknowledgements We acknowledge support from the Air Force Office of Scientific Research (FA and FA ) and the Office of Naval Research (N ). REFERENCES Brandt, L., Sipp, D., Pralits, J. O. & Marquet, O Effect of base-flow variation in noise amplifiers: the flat-plate boundary layer. J. of Fluid Mech. 687, Candler, G. V., Johnson, H. B., Nompelis, I., Gidzak, V. M., Subbareddy, P. K. & Barnhardt, M Development of the US3D code for advanced compressible and reacting flow simulations. In 53rd AIAA Aerospace Sciences Meeting. AIAA

10 10 A. Dwivedi, G. S. Sidharth, J. W. Nichols, G. V. Candler and M. R. Jovanović Chuvakhov, P. V., Borovoy, V. Y., Egorov, I. V., Radchenko, V. N., Olivier, H. & Roghelia, A Effect of small bluntness on formation of Görtler vortices in a supersonic compression corner flow. J. Appl. Mech. Tech. Phys. 58 (6), Cook, D. A., Thome, J., Brock, J. M., Nichols, J. W. & Candler, G. V Understanding effects of nose-cone bluntness on hypersonic boundary layer transition using input-output analysis. In 2018 AIAA Aerospace Sciences Meeting. AIAA Dwivedi, A., Nichols, J. W., Jovanović, M. R. & Candler, G. V Optimal spatial growth of streaks in oblique shock/boundary layer interaction. In 8th AIAA Theoretical Fluid Mechanics Conference. AIAA Dwivedi, A., Sidharth, G. S., Candler, G. V., Nichols, J. W. & Jovanović, M. R Input-output analysis of shock boundary layer interaction. In 2018 Fluid Dynamics Conference. AIAA Ellingsen, T. & Palm, E Stability of linear flow. Phys. Fluids 18 (4), Fedorov, A Transition and stability of high-speed boundary layers. Annu. Rev. Fluid Mech. 43, Hall, P The linear development of Görtler vortices in growing boundary layers. J. Fluid Mech. 130, Hanifi, A., Schmid, P. J. & Henningson, D. S Transient growth in compressible boundary layer flow. Phys. Fluids 8 (3), Inger, G. R Three-dimensional heat-and mass-transfer effects across high-speed reattaching flows. AIAA J. 15 (3), Jeun, J., Nichols, J. W. & Jovanović, M. R Input-output analysis of high-speed axisymmetric isothermal jet noise. Phys. Fluids 28 (4), Jovanović, M. R Modeling, analysis, and control of spatially distributed systems. PhD thesis, University of California, Santa Barbara. Jovanović, M. R. & Bamieh, B Componentwise energy amplification in channel flows. J. Fluid Mech. 534, Ma, Y. & Zhong, X Receptivity of a supersonic boundary layer over a flat plate. Part 1: Wave structures and interactions. J. Fluid Mech. 488, McKeon, B. J. & Sharma, A. S A critical-layer framework for turbulent pipe flow. J. Fluid Mech. 658, Navarro-Martinez, S. & Tutty, O. R Numerical simulation of Görtler vortices in hypersonic compression ramps. Comput. Fluids 34 (2), Nichols, J. W Input/output analysis of f and s mode synchronization in hypersonic boundary layers. Phys. Rev. Fluids Submitted. Fosas de Pando, M. & Schmid, P. J Optimal frequency-response sensitivity of compressible flow over roughness elements. J. Turbul. 18 (4), Ran, W., Zare, A., Hack, M. J. P. & Jovanović, M. R Stochastic receptivity analysis of boundary layer flow. Phys. Rev. Fluids Submitted; also arxiv: Reshotko, E Transient growth: a factor in bypass transition. Phys. Fluids 13 (5), Roghelia, A., Olivier, H., Egorov, I. & Chuvakhov, P Experimental investigation of Görtler vortices in hypersonic ramp flows. Exp. Fluids 58 (10), 139. Schlichting, H. & Gersten, K Boundary-layer theory. Springer. Schmid, P. J Nonmodal stability theory. Annu. Rev. Fluid Mech. 39, Sidharth, G. S., Dwivedi, A., Candler, G. V. & Nichols, J. W Global linear stability analysis of high speed flows on compression ramps. In 47th AIAA Fluid Dynamics Conference. AIAA Sidharth, G. S., Dwivedi, A., Candler, G. V. & Nichols, J. W Onset of threedimensionality in supersonic flow over a slender double wedge. Phys. Rev. Fluids 3 (9), Simeonides, G. & Haase, W Experimental and computational investigations of hypersonic flow about compression ramps. J. Fluid Mech. 283, Zhuang, Y., Tan, H., Liu, Y., Zhang, Y. & Ling, Y High resolution visualization of Görtler-like vortices in supersonic compression ramp flow. J. Vis. 20 (3), Zuccher, S., Tumin, A. & Reshotko, E Optimal disturbances in compressible boundary layers-complete energy norm analysis. In 4th AIAA Theoretical Fluid Mechanics Meeting. AIAA