Uncovering the Varieties of Three-dimensional Hall-MHD Turbulence

TL;DR

This study uses detailed numerical simulations to explore how the statistical properties of 3D Hall-MHD turbulence depend on magnetic Prandtl number and reveals the interplay between different length scales and wave phenomena.

Contribution

It provides new insights into the spectral scaling laws and the dependence on magnetic Prandtl number in 3D Hall-MHD turbulence through comprehensive DNS and theoretical analysis.

Findings

Two inertial regions with different power-law scalings were identified.

Spectral exponents depend on magnetic Prandtl number, with specific values for different regimes.

Theoretical relations between kinetic and magnetic energy spectra were confirmed by DNS results.

Abstract

We carry out extensive pseudospectral direct numerical simulations (DNSs) of decaying three-dimensional (3D) Hall magnetohydrodynamics (3D HMHD) plasma turbulence at three magnetic Prandtl numbers $Pr_{m}=0.1$, $1.0$ and $10.0$. Our DNSs have been designed to uncover the dependence of the statistical properties of 3D HMHD turbulence on $Pr_m$ and to bring out the subtle interplay between three lengths, the kinetic and magnetic dissipation length scales $\eta_u$, and $\eta_b$ and the ion-inertial scale $d_i$, below which we see the manifestations of the Hall term. This interplay, qualitatively apparent from isosurface plots of the moduli of the vorticity and the current density, is exposed clearly by the kinetic-energy and magnetic-energy spectra, $E_u(k)$ and $E_b(k)$, respectively. We find two different inertial regions, In the first inertial region $k<k_{i}\sim1/d_i$, both the kinetic-energy and magnetic-energy spectra, $E_u(k)$ and $E_b(k)$, respectively, display power-law regions with an exponent that is consistent with Kolmogorov-type $-5/3$ scaling, at all values of $Pr_m$. In the second inertial region $k > k_{i}$, the scaling of $E_b(k)$ depends upon $Pr_M$: At $Pr_{m}=0.1$, the spectral-scaling exponent is $-17/3$, but for $Pr_{m}=1$ and $10$ this exponent is $-11/3$. We then show theoretically that $E_u(k) \sim k^2 E_b(k)$ for $Pr_m \ll 1$ and $E_b(k) \sim k^2 E_u(k)$ for $Pr_m \gg 1$; our DNS results are consistent with our theoretical predictions. We examine, furthermore, left- and right-polarised fluctuations of the fields that lead, respectively, to the dominance of ion-cyclotron or whistler waves.