3D numerical simulations of oscillations in solar prominences

¹ Departamet de Física, Universitat de les Illes Balears, 07122 Palma de Mallorca, Spain
e-mail: adrover@uib.es
² Institut d’Aplicacions Computacionals de Codi Comunitari (IAC 3), Universitat de les Illes Balears, 07122 Palma de Mallorca, Spain

Received: 4 October 2019
Accepted: 25 November 2019

Abstract

Context. Oscillations in solar prominences are a frequent phenomenon, and they have been the subject of many studies. A full understanding of the mechanisms that drive them and their attenuation has not been reached yet, however.

Aims. We numerically investigate the periodicity and damping of transverse and longitudinal oscillations in a 3D model of a curtain-shaped prominence.

Methods. We carried out a set of numerical simulations of vertical, transverse and longitudinal oscillations with the high-order finite-difference Pencil Code. We solved the ideal magnetohydrodynamic equations for a wide range of parameters, including the width (w_x) and density (ρ_p0) of the prominence, and the magnetic field strength (B) of the solar corona. We studied the periodicity and attenuation of the induced oscillations.

Results. We found that longitudinal oscillations can be fit with the pendulum model, whose restoring force is the field-aligned component of gravity, but other mechanisms such as pressure gradients may contribute to the movement. On the other hand, transverse oscillations are subject to magnetic forces. The analysis of the parametric survey shows, in agreement with observational studies, that the oscillation period (P) increases with the prominence width. For transverse oscillations we obtained that P increases with density and decreases with B. For longitudinal oscillations we also found that P increases with ρ_p0, but there are no variations with B. The attenuation of transverse oscillations was investigated by analysing the velocity distribution and computing the Alfvén continuum modes. We conclude that resonant absorption is the mean cause. Damping of longitudinal oscillations is due to some kind of shear numerical viscosity.

Conclusions. Our model is a good approximation of a prominence body that nearly reproduces the observed oscillations. However, more realistic simulations that include other terms such as non-adiabatic processes or partially ionised plasmas are necessary to obtain better results.