Driving solar coronal MHD simulations on high-performance computers

DOI10.1080/03091929.2019.1643849zbMATH Open1482.76141arXiv1908.08557OpenAlexW2965429044WikidataQ114640051 ScholiaQ114640051MaRDI QIDQ5065705FDOQ5065705

Authors: Philippe-A. Bourdin 

Publication date: 22 March 2022

Published in: Geophysical & Astrophysical Fluid Dynamics (Search for Journal in Brave)

Abstract: The quality of today's research is often tightly limited to the available computing power and scalability of codes to many processors. For example, tackling the problem of heating the solar corona requires a most realistic description of the plasma dynamics and the magnetic field. Numerically solving such a magneto-hydrodynamical (MHD) description of a small active region (AR) on the Sun requires millions of computation hours on current high-performance computing (HPC) hardware. The aim of this work is to describe methods for an efficient parallelization of boundary conditions and data input/output (IO) strategies that allow for a better scaling towards thousands of processors (CPUs). The Pencil Code is tested before and after optimization to compare the performance and scalability of a coronal MHD model above an AR. We present a novel boundary condition for non-vertical magnetic fields in the photosphere, where we approach the realistic pressure increase below the photosphere. With that, magnetic flux bundles become narrower with depth and the flux density increases accordingly. The scalability is improved by more than one order of magnitude through the HPC-friendly boundary conditions and IO strategies. This work describes also the necessary nudging methods to drive the MHD model with observed magnetic fields from the Sun's photosphere. In addition, we present the upper and lower atmospheric boundary conditions (photospheric and towards the outer corona), including swamp layers to diminish perturbations before they reach the boundaries. Altogether, these methods enable more realistic 3D MHD simulations than previous models regarding the coronal heating problem above an AR -- simply because of the ability to use a large amount of CPUs efficiently in parallel.

Full work available at URL: https://arxiv.org/abs/1908.08557

Recommendations

• Using high performance Fortran for magnetohydrodynamic simulations
• scientific article; zbMATH DE number 4053284
• Non-Fourier description of heat flux evolution in 3D MHD simulations of the solar corona
• scientific article; zbMATH DE number 1968601
• A parallel implementation of an MHD code for the simulation of mechanically driven, turbulent dynamos in spherical geometry

zbMATH Keywords

coronamagneto-hydrodynamicshigh-performance computingastrophysicsSun

Mathematics Subject Classification ID

Magnetohydrodynamics and electrohydrodynamics (76W05) Hydrodynamic and hydromagnetic problems in astronomy and astrophysics (85A30)

Cites Work

• Transport Phenomena in a Completely Ionized Gas

Cited In (3)

• Introduction
• High-order shock-capturing methods for modeling dynamics of the solar atmosphere
• Title not available (Why is that?)

Uses Software

• MURaM

This page was built for publication: Driving solar coronal MHD simulations on high-performance computers

Report a bug (only for logged in users!)Click here to report a bug for this page (MaRDI item Q5065705)