A global spectral-Galerkin investigation of a Rayleigh-Taylor instability in plasma using an MHD-Boussinesq model

This paper presents a new efficient algorithm based on the spectral-Galerkin numerical approximations complemented by a magnetohydrodynamics-Boussinesq model and a new solver for studying the development of a Rayleigh-Taylor (RT) instability. We use the Shenfun computational framework in the Cartesian coordinates, which gives the spectral order and accuracy for the considered model based on the magnetohydrodynamics equations and the Boussinesq conjecture. The numerical simulations were conducted for each two- and three-dimensional case, both with and without an external static magnetic field. The validity of the numerical results was examined by comparing the calculated squared L-2-norm of the density parameter with the linear stability analysis. We also examined the effects of a uniform tangential magnetic field on the onset and growth of an RT instability at different magnetic field strengths. The analysis of the effectiveness of the presented method suggests that it can be modified for further research on two-component plasma.(c) 2023 Author(s). All article content, except where otherwise noted, is licensed under a Creative Commons Attribution (CC BY) license(http://creativecommons.org/licenses/by/4.0/)

Automated summary

This paper presents a new efficient algorithm based on spectral-Galerkin numerical approximations and a magnetohydrodynamics-Boussinesq model for studying the development of a Rayleigh-Taylor instability. The validity of the algorithm was confirmed through numerical simulations and comparisons with linear stability analysis. The effects of an external magnetic field on the instability were also investigated.