Legolas 2.0: Improvements and extensions to an MHD spectroscopic framework

Abstract

We report on recent extensions and improvements to the Legolas code, which is an open-source, finite element-based numerical framework to solve the linearised (magneto)hydrodynamic equations for a three-dimensional force- and thermally balanced state with a nontrivial one-dimensional variation. The standard Fourier modes imposed give rise to a complex, generalised non-Hermitian eigenvalue problem which is solved to quantify all linear wave modes of the given system in either Cartesian or cylindrical geometries. The framework now supports subsystems of the eight linearised MHD equations, allowing for pure hydrodynamic setups, only one-dimensional density/temperature/velocity variations, or the option to treat specific closure relations. We discuss optimisations to the internal datastructure and eigenvalue solvers, showing a considerable performance increase in both execution time and memory usage. Additionally the code now has the capability to fully visualise eigenfunctions associated with given wave modes in multiple dimensions, which we apply to standard Kelvin-Helmholtz and Rayleigh-Taylor instabilities in hydrodynamics, thereby providing convincing links between linear stability analysis and the onset of non-linear phenomena.

arXiv commentAccepted for CPC

References in corpus (12)

• Enzo: An Adaptive Mesh Refinement Code for Astrophysics
• The Black Hole Accretion Code
• MPI-AMRVAC for Solar and Astrophysics
• MPI-AMRVAC 2.0 for Solar and Astrophysical Applications
• Computation of multi-region relaxed magnetohydrodynamic equilibria
• Constrained transport and adaptive mesh refinement in the Black Hole Accretion Code
• MPI-AMRVAC 3.0: updates to an open-source simulation framework
• Linear MHD stability studies with the STARWALL code
• Legolas: a modern tool for magnetohydrodynamic spectroscopy
• Magnetohydrodynamic spectroscopy of a non-adiabatic solar atmosphere
• The Super-Alfvénic Rotational Instability in accretion disks about black holes
• The effect of linear background rotational flows on magnetoacoustic modes of a photospheric magnetic flux tube

Cited by in corpus (5)

• A spectroscopic investigation of thermal instability for cylindrical equilibria with background flow
• Parametric Survey of Nonaxisymmetric Accretion Disk Instabilities: Magnetorotational Instability to Super-Alfvénic Rotational Instability
• Modification of the resistive tearing instability with Joule heating by shear flow
• Neural network classification of eigenmodes in the magnetohydrodynamic spectroscopy code Legolas
• Thermal instability in coronal loops: linking eigenvalue spectra to time-dependent evolution