The 3D MHD code GOEMHD3 for large-Reynolds-number astrophysical plasmas

TL;DR

GOEMHD3 is a highly scalable, efficient 3D MHD simulation code designed for large-Reynolds-number astrophysical plasmas, enabling high-resolution, large-domain simulations with excellent parallel performance.

Contribution

The paper introduces GOEMHD3, a new massively parallel MHD code optimized for large-scale, high-resolution astrophysical plasma simulations with modular and architecture-adaptive features.

Findings

Validated with ideal and weakly dissipative plasma tests.

Achieved high-resolution simulation of solar corona with 2048^3 grid points.

Demonstrated scalability on over 30,000 processor cores.

Abstract

The numerical simulation of turbulence and flows in almost ideal, large-Reynolds-number astrophysical plasmas motivates the implementation of almost conservative MHD computer codes. They should efficiently calculate, use highly parallelized schemes scaling well with large numbers of CPU cores, allows to obtain a high grid resolution over large simulation domains and which can easily be adapted to new computer architectures as well as to new initial and boundary conditions, allow modular extensions. The new massively parallel simulation code GOEMHD3 enables efficient and fast simulations of almost ideal, large-Reynolds-number astrophysical plasma flows, well resolved and on huge grids covering large domains. Its abilities are validated by major tests of ideal and weakly dissipative plasma phenomena. The high resolution ($2048^3$ grid points) simulation of a large part of the solar corona above an observed active region proved the excellent parallel scalability of the code using more than 30.000 processor cores.