Charles Explorer logo
🇬🇧

3D Finite Element MHD code for simulations in Solar atmosphere

Publication at Faculty of Mathematics and Physics |
2019

Abstract

We present an innovative numerical code designed to study multiscale phenomena in the solar atmosphere such as magnetic reconnection, solar flares, quiescent and eruptive prominences as well as their oscillations. The three-dimensional code is based on Discontinuous Galerkin Finite Element Method and is fully adaptive.

It is also designed to handle shocks well and can guarantee zero divergence of magnetic field naturally. Therefore, it should be well suited for simulations that would fill the gap between particle processes and observed large-scale structures.

This makes the code relevant for comparing with both solar radio spectra and potentially also images. We show several tests of the code including those connected to the Solar atmosphere such as perturbed long current sheets and 3D Titov-Demoulin equilibrium.

We also compare some of the results to an older code based on Finite Differences.