Stability of pulsatile quasi-two-dimensional duct flows under a transverse magnetic field

TL;DR

This study numerically investigates how pulsatile wall oscillations can significantly reduce flow instability thresholds in magnetohydrodynamic duct flows under transverse magnetic fields, relevant for fusion reactor cooling systems.

Contribution

It demonstrates that pulsatile wall oscillations can greatly lower the critical Reynolds number in MHD duct flows, especially at high magnetic field strengths, with potential applications in fusion reactor cooling.

Findings

Up to 70% reduction in critical Reynolds number at low magnetic fields.

Over 90% reduction achieved at higher magnetic fields.

Nonlinear simulations show no turbulence despite linear instability.

Abstract

This manuscript has been accepted for publication in Physical Review Fluids, see https://journals.aps.org/prfluids/accepted/53075Se8O0b1b109b1cc0061b280aaa122f0f92dc. The stability of a pulsatile quasi-two-dimensional duct flow was numerically investigated. Flow was driven, in concert, by a constant pressure gradient and by the synchronous oscillation of the lateral walls. This prototypical setup serves to aid understanding of unsteady magnetohydrodynamic flows in liquid metal coolant ducts subjected to transverse magnetic fields, motivated by the conditions expected in magnetic confinement fusion reactors. A wide range of wall oscillation frequencies and amplitudes were simulated. Focus was placed on the driving pulsation optimized for the greatest reduction in the critical Reynolds number, for a range of friction parameters $H$ (proportional to magnetic field strength). An almost $70$% reduction in the critical Reynolds number, relative to that for the steady base flow, was obtained toward the hydrodynamic limit ($H=10^{-7}$), while just over a $90$% reduction was obtained by $H=10$. For all oscillation amplitudes, increasing $H$ consistently led to an increasing percentage reduction in the critical Reynolds number. This is a promising result, given fusion relevant conditions of $H \geq 10^4$. These reductions were obtained by selecting a frequency that both ensures prominent inflection points, and a growth in perturbation energy in phase with deceleration of the base flow. Nonlinear simulations at the optimized frequency and amplitude still satisfied the no net growth condition at the greatly reduced critical Reynolds numbers. However, although the linear mode undergoes a symmetry breaking process, turbulence was not triggered. Nonlinear base flow modulation also arrested the linear decay of the perturbation, with exponential growth not observed at supercritical Reynolds numbers.