vix.ing · top · new · best · stats · spec

Interaction of Large- and Small-scale Dynamos in Isotropic Turbulent Flows from GPU-accelerated Simulations

2020/12/16 by Miikka S. Väisälä, Johannes Pekkilä, M. J. Käpylä +5
Biochemistry, Genetics and Molecular Biology · Physics and Astronomy · #Astro and Planetary Science #Computational physics #Direct numerical simulation #Dynamo #Geomagnetism and Paleomagnetism Studies #Homogeneous isotropic turbulence #Magnetic Prandtl number #Magnetic Reynolds number #Magnetic field #Magnetohydrodynamics #Mechanics #Physics #Quantum mechanics #Reynolds number #Solar and Space Plasma Dynamics #Turbulence #astro-ph.SR #physics.comp-ph #physics.flu-dyn

paper · pdf · doi:10.3847/1538-4357/abceca

22 pages, 23 figures, 2 tables, Accepted for publication in the Astrophysical Journal

arxiv created 2020/12/16 · openalex created_date 2020/12/21 · openalex publication_date 2021/02/01 · arxiv updated 2021/02/10 · openalex updated_date 2026/08/06

Abstract

Abstract Magnetohydrodynamical (MHD) dynamos emerge in many different astrophysical situations where turbulence is present, but the interaction between large-scale dynamos (LSDs) and small-scale dynamos (SSDs) 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 nonhelical forcing was employed to separate the effects of LSD and SSD in a periodic domain. Reynolds numbers ( ) up to ≈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 ≈35 in single-node performance compared to the widely used multi-CPU MHD solver Pencil Code . We estimated the growth rates from both the averaged magnetic fields and their power spectra. At low LSD growth dominates, but at high SSD appears to dominate in both helically and nonhelically forced cases. Pure SSD growth rates follow a logarithmic scaling as a function of . Probability density functions of the magnetic field from the growth stage exhibit SSD behavior in helically forced cases even at intermediate . 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 α quenching. Our results are consistent with the SSD inhibiting the growth of the LSD at moderate , while the dynamo growth is enhanced at higher .

Citations