A numerical investigation of quasi-static magnetoconvection with an imposed horizontal magnetic field

Abstract

Quasi-static Rayleigh-B\'enard convection with an imposed horizontal magnetic field is investigated numerically for Chandrasekhar numbers up to $Q=10^6$ with stress free boundary conditions. Both $Q$ and the Rayleigh number ($Ra$) are varied to identify the various dynamical regimes that are present in this system. We find three primary regimes: (I) a two-dimensional (2D) regime in which the axes of the convection rolls are oriented parallel to the imposed magnetic field; (II) an anisotropic three-dimensional (3D) regime; and (III) a mean flow regime characterized by a large scale horizontal flow directed transverse to the imposed magnetic field. The transition to 3D dynamics is preceded by a series of 2D transitions in which the number of convective rolls decreases as $Ra$ is increased. For sufficiently large $Q$, there is an eventual transition to two rolls just prior to the 2D/3D transition. The 2D/3D transition occurs when inertial forces become comparable to the Lorentz force, i.e. when $\sqrt{Q}/Re = O(1)$; 2D, magnetically constrained states persist when $\sqrt{Q}/Re \gtrsim O(1)$. Within the 2D regime we find heat and momentum transport scalings that are consistent with the hydrodynamic asymptotic predictions of Chini and Cox [Phys. Fluids \textbf{21}, 083603 (2009)]: the Nusselt number ($Nu$) and Reynolds number ($Re$) scale as $Nu \sim Ra^{1/3}$ and $Re \sim Ra^{2/3}$, respectively. For $Q=10^6$, we find that the scaling behavior of $Nu$ and $Re$ breaks down at large values of $Ra$ due to a sequence of bifurcations and the eventual manifestation of mean flows.