Unifying heat transport model for the transition between buoyancy-dominated and Lorentz-force-dominated regimes in quasistatic magnetoconvection

Abstract

In magnetoconvection, the flow of electromagnetically conductive fluid is driven by a combination of buoyancy forces, which create the fluid motion due to thermal expansion and contraction, and Lorentz forces, which distort the convective flow structure in the presence of a magnetic field. The differences in the global flow structures in the buoyancy-dominated and Lorentz-force-dominated regimes lead to different heat transport properties in these regimes, reflected in distinct dimensionless scaling relations of the global heat flux (Nusselt number $\textrm{Nu}$) versus the strength of buoyancy (Rayleigh number $\textrm{Ra}$) and electromagnetic forces (Hartmann number $\textrm{Ha}$). Here, we propose a theoretical model for the transition between these two regimes for the case of a quasistatic vertical magnetic field applied to a convective fluid layer confined between two isothermal, a lower warmer and an upper colder, horizontal surfaces. The model suggests that the scaling exponents $\gamma$ in the buoyancy-dominated regime, $\textrm{Nu}\sim\textrm{Ra}^\gamma$, and $\xi$ in the Lorentz-force-dominated regime, $\textrm{Nu}\sim(\textrm{Ha}^{-2}\textrm{Ra})^\xi$, are related as $\xi=\gamma/(1-2\gamma)$, and the onset of the transition scales with $\textrm{Ha}^{-1/\gamma}\textrm{Ra}$. These theoretical results are supported by our Direct Numerical Simulations for $10\leq \textrm{Ha}\leq2000$, Prandtl number $\textrm{Pr}=0.025$ and $\textrm{Ra}$ up to $10^9$ and data from the literature.