Journal of Fluid Mechanics

Papers

Instabilities in free-surface Hartmann flow at low magnetic Prandtl numbers

D. GIANNAKISa1 c1, R. ROSNERa1a2a3 and P. F. FISCHERa3

a1 Department of Physics, University of Chicago, Chicago, IL 60637, USA

a2 Department of Astronomy and Astrophysics, University of Chicago, Chicago, IL 60637, USA

a3 Argonne National Laboratory, Argonne, IL 60439, USA

Abstract

We study the linear stability of the flow of a viscous electrically conducting capillary fluid on a planar fixed plate in the presence of gravity and a uniform magnetic field, assuming that the plate is either a perfect electrical insulator or a perfect conductor. We first confirm that the Squire transformation for magnetohydrodynamics is compatible with the stress and insulating boundary conditions at the free surface but argue that unless the flow is driven at fixed Galilei and capillary numbers, respectively parameterizing gravity and surface tension, the critical mode is not necessarily two-dimensional. We then investigate numerically how a flow-normal magnetic field and the associated Hartmann steady state affect the soft and hard instability modes of free-surface flow, working in the low-magnetic-Prandtl-number regime of conducting laboratory fluids (Pm ≤ 10−4). Because it is a critical-layer instability (moderately modified by the presence of the free surface), the hard mode exhibits similar behaviour as the even unstable mode in channel Hartmann flow, in terms of both the weak influence of Pm on its neutral-stability curve and the dependence of its critical Reynolds number Rec on the Hartmann number Ha. In contrast, the structure of the soft mode's growth-rate contours in the (Re, α) plane, where α is the wavenumber, differs markedly between problems with small, but non-zero, Pm and their counterparts in the inductionless limit, Pm ↘ 0. As derived from large-wavelength approximations and confirmed numerically, the soft mode's critical Reynolds number grows exponentially with Ha in inductionless problems. However, when Pm is non-zero the Lorentz force originating from the steady-state current leads to a modification of Rec(Ha) to either a sub-linearly increasing or a decreasing function of Ha, respectively for problems with insulating or perfectly conducting walls. In insulating-wall problems we also observe pairs of counter-propagating Alfvén waves, the upstream-propagating wave undergoing an instability driven by energy transferred from the steady-state shear to both of the velocity and magnetic degrees of freedom. Movies are available with the online version of the paper.

(Received August 07 2008)

(Revised April 15 2009)

(Accepted April 18 2009)

Key words

  • Magnetohydrodynamics;
  • Interfacial flows;
  • Instability

Correspondence:

c1 Email address for correspondence: dg227@uchicago.edu

Metrics