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

TL;DR

This paper develops a theoretical model describing the transition in heat transport behavior between buoyancy-driven and Lorentz-force-driven regimes in quasistatic magnetoconvection, supported by numerical simulations and literature data.

Contribution

It introduces a unified scaling relation linking heat transport exponents in both regimes and predicts the transition onset based on magnetic and buoyancy forces.

Findings

Derived a relation between scaling exponents in different regimes.

Validated the model with direct numerical simulations.

Identified the transition scaling with magnetic and buoyancy parameters.

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.