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

TL;DR

This study numerically explores quasi-static magnetoconvection under an imposed horizontal magnetic field, identifying distinct dynamical regimes and scaling laws, and analyzing transitions influenced by magnetic and inertial forces.

Contribution

It provides a detailed numerical analysis of the regimes and transitions in horizontal magnetic field magnetoconvection, extending understanding of flow behaviors and scaling laws in this context.

Findings

Identification of three primary regimes: 2D, anisotropic 3D, and mean flow.

Transition to 3D occurs when inertial forces balance Lorentz forces, with specific scaling conditions.

Heat and momentum transport follow predicted scalings in the 2D regime, with breakdown at high Rayleigh numbers.

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.