论文标题

来自gpu加速模拟的各向同性湍流中的大型和小尺度的发电机的相互作用

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

论文作者

Väisälä, Miikka S., Pekkilä, Johannes, Käpylä, Maarit J., Rheinhardt, Matthias, Shang, Hsien, Krasnopolsky, Ruben

论文摘要

在存在湍流的许多不同天体情况下出现磁流失动力学(MHD)发电机,但是大规模(LSD)和小规模的Dynmos(SSD)之间的相互作用尚未完全了解。我们对磁力MHD中的各向同性强迫驱动的湍流发电机进行了系统研究,并具有磁性统一数量的磁性统一数,重点是指数级生长阶段。螺旋强迫和非螺旋强迫都被用来分开LSD和SSD在周期域中的影响。检查了雷诺数字(RM)最多$ \ $ \ 250美元,并用于收敛检查的多种决议。我们使用Astaroth代码进行了模拟,该代码旨在加速图形处理单元(GPU)的3D模具计算,并使用与对等通信的多个GPU。与广泛使用的多CPU MHD MHD求解器铅笔代码相比,我们观察到了$ \ $ \ 35 $的速度约为35美元。我们估计了平均磁场及其功率光谱的增长率。在低RM处,LSD的生长占主导地位,但在高RM下,SSD在螺旋和非螺旋强迫病例中似乎都占主导地位。纯SSD的增长率遵循对数缩放率随RM的函数。磁场的概率密度功能从生长阶段发出,即使在中间RM处,在螺旋强迫的情况下也表现出SSD行为。我们使用二阶相关近似(SOCA)等封闭来估计平均场湍流传输系数。它们产生的增长率类似于直接测量的增长率,并提供了$α$淬火的证据。我们的结果与SSD一致,SSD抑制LSD在中等RM处的生长,而Dynamo的生长在较高的RM处增强。

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 $α$ 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.

扫码加入交流群

加入微信交流群

微信交流群二维码

扫码加入学术交流群,获取更多资源