Three-dimensional global instability analysis for high-speed boundary layer flow

Herein, a numerical approach for the three-dimensional global instability analysis of high-speed boundary layer is developed. The approach is based on the Krylov–Schur method for large-scale eigenproblems. A shift-and-invert spectral transformation is employed to accelerate the convergence for eigenvalues of interest. Iterative linear solver is adopted in the transformation to reduce memory usage by avoiding the fill-in problem of direct solver and the explicit construction of matrices. A Robin condition based on BiGlobal instability analysis is imposed at the inflow boundary, and a non-reflecting condition derived from inviscid flow is used for the far-field boundary. The approach is validated in a Blasius boundary layer. According to a comparison between the performances of several algorithms in this case, the conjugate gradient squared method with no preconditioning is chosen as the linear solver. The results provided by three-dimensional global instability analysis consist with those from linear stability theory and BiGlobal instability analysis. A further demonstration of our approach is provided for a three-dimensional roughness-disturbed boundary layer. The fluctuations extracted from direct numerical simulation results illustrate that the symmetric component is the dominant unstable mode in roughness wake, as it grows twice as fast as the antisymmetric component. Two eigenvalues are obtained from the three-dimensional global instability analysis, and these represent the symmetric and antisymmetric unstable modes. The wavelengths and growth rates of both modes delivered by instability analysis agree well with the direct numerical simulation results, and the profiles of disturbance amplitude match the features of the varicose and sinuous modes.