The transient response of optimal linear perturbations of liquid metal flow under a strong axial magnetic field in an electrically insulated rectangular duct is considered. The focus is on the subcritical regime, below the onset of von Kármán vortex shedding, to determine the role of optimal disturbances in developing wake instabilities. In this configuration, the flow is quasi-two-dimensional and can be solved over a two-dimensional domain. Parameter ranges considered include Reynolds numbers \documentclass[12pt]{minimal}\begin{document}$50 \le \mbox{\textit {Re}}\lesssim 2100$\end{document}50≤Re≲2100, modified Hartmann numbers \documentclass[12pt]{minimal}\begin{document}$50 \le {\mbox{\textit {Ha}}^\star }\lesssim 500$\end{document}50≤Ha★≲500, and blockage ratios 0.1 ⩽ β ⩽ 0.4. In some instances, the optimal disturbances are found to generate energy growth of greater than four orders of magnitude. Variation in the wake recirculation length in the steady flow regime is determined as a function of Reynolds number, Hartman number, and blockage ratio, and a universal expression is proposed. For all β, the energy amplification of the disturbances is found to decrease significantly with increasing Hartmann number and the peak growth shifts towards smaller times. The optimal initial disturbances are consistently located in the vicinity of the boundary layer separation from the cylinder, and the structure of these disturbances is consistent for all Hartmann numbers and blockage ratios considered in this study. The time evolution of the optimal perturbations is presented, and is shown to correspond to sinuous oscillations of the shear layer downstream of the wake recirculation. The critical Reynolds number for the onset of growth at different Hartmann numbers and blockage ratios is determined. It is found that it increases rapidly with increasing Hartmann number and blockage ratio. For all β, the peak energy amplification grows exponentially with \documentclass[12pt]{minimal}\begin{document}$\mbox{\textit {Re}}$\end{document}Re at low and high Hartmann numbers. Direct numerical simulation in which the inflow is perturbed by a random white noise confirms the predictions arising from the transient growth analysis: that is, the perturbation excites and feeds energy into the global mode.
Read full abstract