Interaction of large- and small-scale dynamos in isotropic turbulent flows from GPU-accelerated simulations

TL;DR

This study investigates how large-scale and small-scale magnetic dynamos interact in turbulent flows using GPU-accelerated simulations, revealing their dominance depending on magnetic Reynolds number and providing insights into dynamo mechanisms.

Contribution

The paper presents a systematic GPU-based simulation study of isotropic MHD dynamos, analyzing the interplay between LSD and SSD across various Reynolds numbers with new insights into their growth rates and quenching effects.

Findings

At low Rm, large-scale dynamo dominates.

At high Rm, small-scale dynamo dominates.

SSD growth rates scale logarithmically with Rm.

Abstract

Magnetohydrodynamical (MHD) dynamos emerge in many different astrophysical situations where turbulence is present, but the interaction between large-scale (LSD) and small-scale dynamos (SSD) is not fully understood. We performed a systematic study of turbulent dynamos driven by isotropic forcing in isothermal MHD with magnetic Prandtl number of unity, focusing on the exponential growth stage. Both helical and non-helical forcing was employed to separate the effects of LSD and SSD in a periodic domain. Reynolds numbers (Rm) up to $\approx 250$ were examined and multiple resolutions used for convergence checks. We ran our simulations with the Astaroth code, designed to accelerate 3D stencil computations on graphics processing units (GPUs) and to employ multiple GPUs with peer-to-peer communication. We observed a speedup of $\approx 35$ in single-node performance compared to the widely used multi-CPU MHD solver Pencil Code. We estimated the growth rates both from the averaged magnetic fields and their power spectra. At low Rm, LSD growth dominates, but at high Rm SSD appears to dominate in both helically and non-helically forced cases. Pure SSD growth rates follow a logarithmic scaling as a function of Rm. Probability density functions of the magnetic field from the growth stage exhibit SSD behaviour in helically forced cases even at intermediate Rm. We estimated mean-field turbulence transport coefficients using closures like the second-order correlation approximation (SOCA). They yield growth rates similar to the directly measured ones and provide evidence of $\alpha$ quenching. Our results are consistent with the SSD inhibiting the growth of the LSD at moderate Rm, while the dynamo growth is enhanced at higher Rm.