Data-driven transient growth analysis

The transient growth of disturbances made possible by the non-normality of the linearized Navier-Stokes equations plays an important role in bypass transition for many shear flows. Transient growth is typically quantified by the maximum energy growth among all possible initial disturbances, which is given by the largest squared singular value of the matrix exponential of the linearized Navier-Stokes operator. In this paper, we propose a data-driven approach to studying transient growth wherein we calculate optimal initial conditions, the resulting responses, and the corresponding energy growth directly from flow data. Mathematically, this is accomplished by optimizing the growth over linear combinations of input and output data pairs. We also introduce a regularization to mitigate the sensitivity to noisy measurements and unwanted nonlinearity. The data-driven method simplifies and broadens the application of transient growth analysis -- it removes the burden of writing a new code or linearizing an existing one, alleviates the computational expense for large problems, eliminates the challenge of obtaining a well-posed spatial propagator for spatial growth analyses, and enables the direct application of transient growth analysis to experimental data. We validate the data-driven method using a linearized Ginzburg-Landau model problem corrupted by process and measurement noise and obtain good agreement between the data-driven and the standard operator-based results. We then apply the method to study the spatial transient growth of disturbances in a transitional boundary layer using data from the Johns Hopkins Turbulence Database. Our method successfully identifies the optimal output response and provides plausible estimates of the transient spatial energy growth at various spanwise wavenumbers.