Issue 
A&A
Volume 633, January 2020



Article Number  A113  
Number of page(s)  13  
Section  The Sun  
DOI  https://doi.org/10.1051/00046361/201936841  
Published online  20 January 2020 
3D numerical simulations of oscillations in solar prominences
^{1}
Departamet de Física, Universitat de les Illes Balears, 07122 Palma de Mallorca, Spain
email: adrover@uib.es
^{2}
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
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 curtainshaped prominence.
Methods. We carried out a set of numerical simulations of vertical, transverse and longitudinal oscillations with the highorder finitedifference 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 fieldaligned 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 nonadiabatic processes or partially ionised plasmas are necessary to obtain better results.
Key words: Sun: filaments, prominences / Sun: oscillations / methods: numerical
© ESO 2020
1. Introduction
Solar prominence (or filament) oscillations have been the subject of study for decades. The first detection of a winking filament dates back to 1930 (see Hyder 1966). Since then, diverse researches have aimed at establishing theoretical models that explain the examined oscillations. Observations together with the theoretical models have allowed us to expand our knowledge of solar prominences and have facilitated, for example, the estimation of physical parameters that are difficult to measure, such as the minimum magnetic field strength that supports these plasma structures, or the radius of the curvature of magnetic dips (R) (Ballester 2014; Luna et al. 2014).
Because observational techniques have improved, the oscillations are analysed in more detail so that we can classify this phenomenon as large amplitude oscillations (LAO) or small amplitude oscillations (SAO; see Oliver & Ballester 2002). The first group is characterised by velocity amplitudes of about 20 km s^{−1}, and they are associated with an energetic event, such as solar flares (Jing et al. 2003; Li & Zhang 2012), Moreton waves (Eto et al. 2002; Gilbert et al. 2008), Extreme ultraviolet Imaging Telescope (EIT) waves (Okamoto et al. 2004), or coronal shock waves (Hershaw et al. 2011; Shen et al. 2014a) (see Tripathi et al. 2009, for a more detailed review of LAOs in prominences). In general, LAOs have been considered as a rare phenomenon, but Luna et al. (2018) have catalogued 196 filament oscillations in a period of six months (January to June 2014). Ninety of these events have been classified as LAOs. This indicates that LAOs are a frequent phenomenon. The first study of LAOs was carried out by Ramsey & Smith (1966), who investigated 11 winking filaments disturbed by flares. Winking filaments are periodic motions that are intermittently observed throughout different images shifted from the H_{α} linecentre by vertical oscillations; this causes a Doppler shift. The authors determined that the examined events oscillated at their own characteristic frequency. Then, Hyder (1966) used the results of Ramsey & Smith (1966) to calculate the radial component of the filament magnetic field and the effective viscosity coefficients. He considered the movement as a damped harmonic oscillator embedded in a viscous fluid and in a uniform gravitational field, whose restoring force is the magnetic tension. Kleczek & Kuperus (1969) studied horizontal transverse oscillations as a damped free oscillator with the magnetic tension as the restoring force, and radiation of acoustic waves as the attenuation mechanism. From these assumptions, they found that the oscillation period (P) for horizontal transverse oscillations is , where 2L is the length of the prominence, B is the magnetic field strength, and ρ_{p} is the prominence density. Jing et al. (2003) later proposed different mechanisms to explain the physics of an oscillatory motion along a filament observed on 2001 October 24. Firstly, they suggested as the driving mechanism the gravitational force when a concave upwards magnetic flux tube is considered as the main configuration of the cool filament thread. The second possibility that they proposed to explain the oscillations was a pressure imbalance that drives plasma motions along field lines. Finally, they suggested as a restoring force the magnetic tension for motions perpendicular to the local magnetic field. Furthermore, Vršnak et al. (2007) detected oscillations along a filament on 2002 January 23 and proposed that the restoring force of the motion was caused by the magnetic pressure gradient along the filament axis, and Li & Zhang (2012), who studied longitudinal oscillations in a solar filament on 2012 April 7, suggested that the restoring force was the coupling between the magnetic tension and gravity.
In contrast, SAOs are events with velocity amplitudes smaller than 10 km s^{−1}, and they are not related to flare activity (Arregui et al. 2012). These periodic motions are mostly understood as the oscillatory modes of magnetohydrodynamic (MHD) waves. To comprehend the nature of the oscillatory modes, the dispersion relation of ideal MHD equations has been derived for different models. Joarder & Roberts (1992a,b, 1993) investigated in a set of three articles the normal oscillation modes of a 3D plasma slab embedded in an uniform longitudinal, transverse, and skewed magnetic field in the absence of gravity. The gravitational term was considered in Oliver et al. (1993). In a more complex configuration, Díaz et al. (2001) examined the MHD waves of a single prominence fibril surrounded by coronal medium. From these articles a wide range of oscillatory modes were obtained and the oscillatory period of the prominences was estimated.
The observational and theoretical studies of oscillations in filaments are corroborated by numerical simulations. Single 1D threads formed in magnetic flux tubes with a concave upwards shape where plasma oscillates have been studied by several authors. Luna & Karpen (2012) studied large amplitude longitudinal oscillations in 47 independent threads. They found that the main restoring force is the projected gravity along the magnetic field, and that the periodicity of longitudinal oscillations depends on the radius of the curvature of the magnetic dip as
where g is the solar gravity acceleration. This model has been called the pendulum model. Zhang et al. (2012) used 1D radiative hydrodynamic simulations that almost reproduced the periodicity of longitudinal oscillations in threads with magnetic dips observed on 2007 February 8. Zhang et al. (2013) carried out a parametric survey to investigate the behaviour of longitudinal oscillations for different filament dimensions and triggering mechanisms. They found that the swaying motion is not determined by the disturbance type and reasserted that the longitudinal oscillations follow the pendulum model. They also found that in addition to the gravity force, the pressure gradient contributes to the restoring force for short threads.
In addition to 1D simulations, 2D numerical studies of oscillations in prominences have also been investigated. Some of these works focused on studying the evolution of cold and dense plasma structures embedded in an isothermal stratified atmosphere permeated by a dipped arcade magnetic configuration. Based on this method, Terradas et al. (2013) analysed the evolution of the density distribution together with other variables such as the involved forces or the velocity field during the relaxation process of the nonequilibrium initial configuration. Then, when the prominence was in equilibrium, they studied the behaviour of induced vertical, longitudinal, and Alfvén waves solving the linearised MHD equations. For validating purposes, Luna et al. (2016) reproduced nonlinear numerical simulations of large amplitude longitudinal oscillations in a nearly equilibrium configuration. They found that the pendulum model is robust even in the range of short prominence field lines. In addition to this type of models, Kraśkiewicz et al. (2016), based on the Pikelner model, studied the dynamics of a prominence perturbed by a pressure disturbance, and Zhang et al. (2019) analysed 2D nonadiabatic MHD simulations of longitudinal oscillations.
To reproduce more realistic situations, more complex 3D MHD numerical simulations are required. In the literature, rare 3D studies have analysed prominence oscillations. Terradas et al. (2016) and Zhou et al. (2018) carried out simulations of prominences embedded in magnetic flux ropes. While Terradas et al. (2016) focused their study on the evolution of the density profile for different resolutions and magnetic twist during the relaxation process and in the development of Kelvin–Helmholtz instabilities, Zhou et al. (2018) investigated the dynamics of longitudinal, vertical, and horizontal transverse oscillations triggered by velocity disturbances. The magnetic flux rope geometry is the most widely accepted filament model because it aligns the longitudinal filament axis along the magnetic field lines. However, the magnetic arcade configuration with a curtainlike prominence have not yet been used in 3D numerical simulations of oscillations in prominences.
An important characteristic of the reported events is that in most cases, the amplitude of the oscillations attenuates. To explain the observed damping, different theoretical mechanisms have been suggested, for example thermal damping such as radiative cooling (Terradas et al. 2001) and heat conduction (Soler et al. 2007), ionneutral collisions in partially ionised plasmas (Forteza et al. 2007), mass flows (Soler et al. 2008), wave leakage (Schutgens & Tóth 1999), or resonant absorption (Goossens et al. 2002) (see the review of Arregui et al. 2018). The numerical studies reviewed in this section made an effort to clarify the nature of the attenuation. Luna & Karpen (2012) proposed that the main damping mechanism for their model is the accretion of mass onto the threads during prominence formation. Zhang et al. (2012) suggested that the nonadiabatic terms are responsible for the attenuation. They measured the importance of radiative loss and thermal conduction and concluded that in their thread model, the first term is responsible for the damping time (t_{d}). Zhang et al. (2019) studied the nonadiabatic effects in the attenuation of longitudinal oscillations in detail. They found that these processes are the primary agent that dissipates the oscillations, but wave leakage also plays an important role in dissipating the kinetic energy. Terradas et al. (2016) have shown in their model that the attenuation of the vertical transverse oscillations is produced by resonant absorption, although Kelvin–Helmholtz instabilities could also provide a mechanism to attenuate the prominences. Despite these attempts to establish attenuation mechanisms, numerical studies caution that part of the damping may be due to computational causes such as numerical viscosity or dissipation. The numerical damping complicates the study of the attenuation mechanisms, and for this reason, many studies disregard the attenuation problem. While for longitudinal oscillations the nonadiabatic processes have been analysed in many studies, the resonant absorption in 3D simulations of transverse oscillations needs a deeper study.
Longitudinal oscillations in solar prominences have been studied numerically by many authors, especially in 1D and 2D, and the pendulum model has been corroborated for different geometries. However, a parametric survey of longitudinal oscillations in 3D simulations is required. On the other hand, different models have been used to compare analytical expressions to the transverse simulated periods. For example, Terradas et al. (2013) and Zhou et al. (2018), in their simulations of vertical oscillations, used the slab and string model of Joarder & Roberts (1992b), respectively, and obtained large differences between the analytical and numerical results. Both studies obtained better correlations with the results of Díaz et al. (2001), however. This paper performs a 3D numerical study of oscillations in a curtainshaped prominence permeated by a magnetic arcade that supports the dense and cool plasma against gravity. The magnetic configuration is based on the model proposed by Kippenhahn & Schlüter (1957). Unlike other works, in which the prominence mass is suspended above the photosphere, we anchored the plasma to the bottom boundary of the computational domain, which in our model represents the photosphere. Largescale filament properties show that the cool plasma of many prominences joins the chromosphere, and possibly the photosphere, through the filament barbs (Parenti 2014). Extensive filaments with a few points of contact with the surface (Hedgerow prominences) are the most common type (Zirin 1998; Engvold et al. 2015). In addition, Shen et al. (2014a) observed transverse oscillations above the solar limb and suggested that the body oscillated as a rigid body with one end anchored on the Sun. However, this structure has not been investigated numerically so far. For these reasons, we considered that our model is a good approximation to study this rich prominence type. The aim is to compute the period and damping time for vertical and longitudinal/transverse horizontal motions triggered by velocity disturbances for different mass distributions. The primary initial configuration is described in Sect. 2 and the numerical aspects are explained in Sect. 3. The computed results are displayed in Sect. 4. Discussion and conclusions are summarised in Sect. 5.
2. Initial reference setup
To numerically study the oscillations in prominences, an isothermal stratified background atmosphere has been configured as in Terradas et al. (2015). The coronal density profile is expressed as ρ_{c} = ρ_{0}e^{−z/Λ}, where ρ_{0} is the density of the reference level located at the bottom boundary of the computational box z = 0. Unlike the works of Zhou et al. (2018) and Zhang et al. (2019), who inserted a temperature transition between the photosphere and the corona, in our model we consider that the level z = 0 represents the photosphere and that it is connected to the corona through the magnetic field (Terradas et al. 2013, 2015). is the density scale height. The velocity terms are normalised to the sound speed km s^{−1} for a coronal temperature of 1 MK. The gravity acceleration is g = 0.274 km s^{−2} and the adiabatic index γ = 5/3. The reference length has been taken as H = 10,000 km. The plasma is permeated by an archshaped magnetic field, composed of the superposition of two forcefree magnetic arcades. This quadrupolar configuration adds a magnetic dip to the structure. The depth of the dip varies with height. In Fig. 1 some magnetic field lines are plotted in the xzplane. At x = 0, the radius of the curvature takes low positive values in the bottom part of the structure and grows until z = 1.55 H, where the field line becomes flat. From this point, the dip is lost and the curvature becomes concave downwards. To avoid high values of the plasmaβ parameter (), which affects the character of the modes, the null point of the magnetic force has been located out of the numerical domain as in Luna et al. (2016). Unlike Terradas et al. (2015), the magnetic configuration has no shear and the magnetic field is invariant in the ydirection. An increase in shear produces an increase in the horizontal component of the magnetic field, leading to a higher magnetic tension. It also changes the depth and width of the dips so that the magnetic support can be different between sheared and nonsheared structures. The study of oscillations in sheared arcades is beyond the scope of this work. To obtain the previously described magnetic configuration, the magnetic vector potential is
Fig. 1. Magnetic field lines in the xzplane. In this plot B_{0} = 10 G, k_{1} = π/2L, k_{2} = 3k_{1}, L = 5 H, z_{0} = −0.2 H, and l_{1}/k_{1} = 1. The ycomponent of the magnetic field is null. The maximum initial Alfvén speed is v_{A0} = 17.7 c_{s0}. 

Open with DEXTER 
For the reference simulation we took B_{0} = 10 G. k_{1} = π/2L and k_{2} = 3k_{1} are parameters related to the lateral extension of the arcades, 2L is the width of the whole structure, l_{1} and l_{2} are a measure of vertical magnetic scale height, and z_{0} = −0.2 H is the position of the null point. The relation l_{1}/k_{1} is related to the shear of the magnetic field lines. In this work the shear is null so that l_{1}/k_{1} = 1. l_{2} is set by the constraint . The magnetic field is B = ∇ × A.
The mass deposition is artificial and instantaneous because we did not study prominence formation here (see Mackay et al. 2010, for a summary of prominence formation models). The deposition maintains the gas pressure at the same level. The density profile consists of a curtainshaped mass distribution expressed as
The term n determines the width of the prominencecorona transition region (PCTR; in this work, n = 4 by default). The width, length, and height of the curtain are w_{x} = 0.3 H, w_{y} = 2.2 H, and w_{z} = 2.6 H, respectively, and the density contrast between the corona and the core of the prominence is ρ_{p0}/ρ_{0} = 80. With these characteristics and by considering that the prominence density is 5.2 × 10^{−2} kg km^{−3}, the total mass of the filament is 1.6 × 10^{11} kg. The foot of the prominence is anchored to the base of the domain box, which represents the base of the corona, therefore the motion is subjected to this condition. The density profile is located centred on the computational box and the filament spine is aligned with the ycoordinate so the magnetic field lines thread the prominence transversely. This configuration is not the most widely used by authors, who accept that in most cases the dense plasma is almost aligned along the magnetic field lines. Observations demonstrate that the angle between the magnetic field lines and the prominence spine is quite small and has a maximum around 25° (TandbergHanssen 1995). However, we decided to use our configuration for simplicity and as a preliminary step before sheared arcades are analysed. In Fig. 2 the described density profile is represented at t = 0 with some selected magnetic field lines.
Fig. 2. Density profile together with some magnetic field lines at t = 0. In this plot n = 4, w_{x} = 0.3 H, w_{y} = 2.2 H, w_{z} = 2.6 H, and ρ_{p0} = 80 ρ_{0}. 

Open with DEXTER 
3. Numerical aspects
The ideal MHD equations we solved are
The set of governing equations are expressed in terms of the entropy rather than the more usual thermal energy or pressure variables. c_{p} is the heat capacity at constant pressure and was assumed equal to one. Φ_{grav} is the gravity potential and j is the electric current. The induction equation is solved in terms of the magnetic vector potential. D/Dt = ∂/∂t + u · ∇ is the advective derivative. The rest of characters have their normal meaning.
The code we used to solve the MHD equations is the Pencil Code^{1}. It is a publicly available model that uses sixthorder finitedifference schemes. Highorder derivative schemes, such as the one used in the Pencil Code, reduce the numerical dissipation, but to obtain consistency in numerical solutions, the code needs small amounts of diffusion to damp out the modes near the Nyquist frequency (Brandenburg 2003). For this reason we introduced a simplified secondorder hyperviscosity term (see Haugen & Brandenburg 2004) and a shock viscosity term at the equation of motion. We did not use artificial viscosity in the other equations. The hyperviscosity is proportional to a diffusion coefficient. This coefficient should be as small as possible but sufficient to reduce the wiggles in the results. In addition, to avoid wiggles, we used fifthorder upwind derivatives for the advection terms u · ∇ ln ρ and u · ∇s. The time step is a thirdorder Runge–Kutta scheme.
The Pencil Code uses three layers of ghost points to implement boundary conditions. In this work, closed boundary conditions were applied (see Terradas et al. 2016). This means that linetying conditions were imposed at all the boundaries of the computational box. A linetying boundary condition sets the three components of the velocity to 0, the normal component of the magnetic field is kept constant, and the density and entropy variables have their spatial derivatives be null. This condition imposed at the bottom boundary is crucial to mimic the purely reflecting conditions of the photosphere and to obtain the magnetic support. Because the Pencil Code uses the magnetic vector potential instead of the magnetic field, to fix the magnetic component perpendicular to the boundary, we applied to the three components of A the condition of antisymmetry relative to the boundary value that vanishes the second derivative of A. The numerical domain consists of a box of 180 × 144 × 90 mesh points. The dimensions of the box are 100 Mm in the xdirection, 80 Mm in the ydirection, and 50 Mm in the vertical component. In this way, we imposed an equidistant grid of 556 km.
4. Results
4.1. Relaxation process
As was mentioned in Sect. 2, the mass deposition is instantaneous. The prominence plasma breaks the initial equilibrium, and then starts to oscillate vertically. Initially, the forcefree magnetic condition implies equilibrium between magnetic pressure and magnetic tension, and gas pressure force counteracts the gravitational force, so that the entire system is stable. However, when the prominence is introduced, the system loses its balance and the gravity force pushes the structure downwards. Immediately, the plasma deforms the field lines, increasing the depth of the magnetic dip. In panel a of Fig. 3 we show the final state of the field lines (blue solid curves) compared with the initial state (black dashed curves). This distortion implies an increase in the zcomponent of the magnetic tension (and to a lesser extent, the magnetic pressure) in such a way that it counteracts the movement, restoring the prominence upwards. After an evolution of approximately 80 min, the Lorentz force almost balances the gravity force.
Fig. 3. Density distribution of the final equilibrium state at t = 90 min. Top panel: density isocontours and magnetic field lines in the plane y = 0. Bottom panel: density isocontours in the plane x = 0. Dashed lines correspond to the initial reference setup, and solid lines show the final equilibrium state at t = 90 min. 

Open with DEXTER 
The total mass of the system remains practically constant, although Fig. 4 shows that the maximum density contrast between the reference level and the plasma slab increases. This involves the redistribution of the density profile. The density distribution of the final equilibrium state at t = 90 min is shown in the top and bottom panels of Fig. 3. In panel a we observe that the prominence widens slightly in the xdirection near its base, where the plasma gathers. To compensate for the accumulation of mass at the base, the filament narrows at the top and reduces its height. In panel b we find that in the ydirection, the redistribution of mass is more noticeable, especially in the inner layers of the prominence, where as a result of the accumulation of plasma, the isocontour ρ = 80 ρ_{0} expands. The evolution of the zcomponent of the centre of mass is represented in Fig. 5. To suppress interactions with the background, the centroid was averaged from a delimited box that only included the prominence body. The green dotted line corresponds to the reference simulation, and, initially, the centre of mass is located approximately at z = 1.1 H, however, it rapidly drops. After two cycles, the prominence oscillates around a new equilibrium height located at z = 1.01 H. This fall agrees with the mass redistribution.
Fig. 4. Time evolution of density, v_{z}, and magnetic field lines for the relaxation process. v_{z} is represented in the yzplane at x = 0 as a 2D slice. The initial configuration for this case is shown in Fig. 2. 

Open with DEXTER 
Fig. 5. Time evolution of the zcomponent of the centre of mass for n = 2 (red), n = 4 (green), and n = 6 (blue). The density average was computed for a delimited region (x ∈ [−0.5 H, 0.5 H], y ∈ [−3 H, 3 H], and z ∈ [0, 3.5 H]). 

Open with DEXTER 
The evolution of the vertical velocity at the point x = 0, y = 0, z = 1.3 H is represented in Fig. 6 (dotted line). At the beginning of the simulation, the system is completely at rest, but the prominence body immediately drops, increasing the zcomponent of the velocity (v_{z}), which rapidly reaches its maximum value. Approximately at t = 3.5 min, the speed reaches positive values, which indicates that the movement begins to ascend towards its original position. This periodic motion is repeated approximately 16 times for about 90 min. Figure 6 shows that the amplitude of the oscillations decreases with time. The observations indicate that the movement describes an attenuated oscillation. We can fit an exponentially decayed sine curve with respect to time, expressed as
where f_{0} is the signal amplitude, ω the frequency, and t_{d} is the damping time. For the reference simulation we obtained a period of P = 2π/ω = 5.6 min, a damping time of t_{d} = 25.1 min, and a velocity amplitude of f_{0} = −0.082 c_{s0}. We represent the fitted curve in Fig. 6 as a solid line.
As in Terradas et al. (2015), in Fig. 4 we represent the evolution of v_{z} in a yzslice at x = 0. As we showed in Fig. 6, the top panel of Fig. 4 shows that at the beginning of the simulation v_{z} takes a negative value at the prominence body, describing a downwards movement. Nevertheless, the global evolution is more complex, and the system also develops vertical shear motions at the lateral edges of the prominence. In the first steps of the simulation, the system evolves a single wide lateral strip, but in the bottom panel of Fig. 4 we observe that at t = 8.4 min the number of bands with opposite speed has grown and their width has narrowed. This is clearer in Fig. 7, where we plot the time evolution of v_{z} as a function of ycoordinate at x = 0 and z = 1.3 H. At t = 2.4 min (solid line) and t = 8.4 min (dashed line), we observe the same pattern as Fig. 4, and at t = 14.4 min (dotted line), we obtain that these shear motions remain localised at the PCTR and that the typical spatial scales have decreased with time. This velocity pattern is generated as a result of resonant absorption, and it is the origin of part of the attenuation. Terradas et al. (2015) also found this positive/negative velocity pattern in v_{z} (see Fig. 3 of their work). They associated this behaviour with the process of mode conversion and phase mixing. Our results agree with the resonant absorption process because the energy conversion takes place at the inhomogeneous layer of the prominence, that is, at the sides of the structure in the transverse direction with respect to the field lines and the direction of motion; in this case, the yedges. This process is also characterised by the decrease with time of the typical spatial scales of the shear stripes.
Fig. 6. v_{z} as a function of time for the relaxation process at a fixed point (0, 0, and 1.3 H). Dots are the values obtained from the numerical simulation, and the continuous grey curve is the fitting result based on Eq. (8). 

Open with DEXTER 
Fig. 7. v_{z} as a function of y (x = 0, z = 1.3 H) at t = 2.4 min (solid line), t = 8.4 min (dashed line), and t = 14.4 min (dotted line). 

Open with DEXTER 
The studies of the resonant absorption as a mechanism for the damping of filament oscillations are usually based on simple models of magnetic flux tubes (see e.g. Arregui et al. 2008; Soler et al. 2009). In these works, the damping time of the fast kink waves caused by the resonant absorption process depends on the width of PCTR. To study the dependence of the attenuation on the PCTR, the relaxation process was simulated by modifying the density profile but keeping the total mass constant, so that for n = 2, ρ_{p0} = 144 ρ_{0}, and for n = 6, ρ_{p0} = 63.7 ρ_{0}. The prominence with the term n = 6 has the steepest profile, and that of n = 2 is the least steep, so that when n increases, the density slope is stronger and leads to thinner transitions between the core and the corona. In Fig. 5 the evolution of the centre of mass is plotted for the three cases. The evolutions of the three curves are similar each other, but we clearly see that for n = 6 (blue curve), the evolution has a softer attenuation and for n = 2 (red curve) we obtain the strongest damping. In agreement with Arregui et al. (2008), we obtain that the wider the PCTR, the stronger the attenuation. We are interested in estimating the ratio of the damping time to the period. From the fitted curve of Fig. 6 we computed t_{d} from the v_{z} series at a localised point. However, in complex configurations we can detect different attenuation times at distinct points. For this reason, in Fig. 5 we calculated the time evolution of the spatially averaged density. For the relaxation process, to calculate the damping time of the three signals of Fig. 5, we firstly applied an empirical mode decomposition (see Terradas et al. 2004) to decompose the signal and select the first intrinsic mode function. In this way, we subtracted the trend of the centre of mass drop. After this, we obtained a damping per period of t_{d}/P = 1.9 for n = 2, t_{d}/P = 3.7 for n = 4, and t_{d}/P = 5.9 for n = 6.
To extend the study of the relaxation process, we studied its evolution for a range of values of the filament width, the density contrast, and the magnetic field strength. The results are shown in Fig. 8. First we investigated the dependence of the period on w_{x}. We varied the width of the prominence from 0.25 H to 0.5 H, and for the other parameters, we used the reference values. When the width is doubled, the total mass is multiplied by a factor of two. To compute the periodicity of the oscillations, we fitted the results for the different time series of v_{z} with Eq. (8) following the same procedure as used in the reference process. The curtain with the widest width, and therefore with the greatest total mass, shows a higher initial zvelocity amplitude and a greater fall of the centroid into the new equilibrium point. Panel a of Fig. 8 indicates that the oscillation period varies with the prominence width. We obtained that the wider the prominence, the longer the period.
Fig. 8. Scatter plots of the oscillation period during the relaxation process as a function of w_{x} (panel a), ρ_{p0} (panel b), and B_{0} (panel c). The values of P were obtained from the fitted curve based on Eq. (8) for the evolution of v_{z} at a fixed point (0, 0, and 1.3 H). 

Open with DEXTER 
We also studied the dependence of the periodicity on prominence mass. We varied ρ_{p0} and kept the reference values of the other variables. Now the total mass of the parametric survey ranges from 7.9 × 10^{10} kg to 1.8 × 10^{11} kg. The highest mass prominence suffers the greatest drop into the new equilibrium point and has the largest zvelocity amplitude. In panel b of Fig. 8 we show an increase in period with total mass. In panels a and b we observe an almost lineal dependence of the period on the other two parameters. The parameter values were chosen so that the curtain with a width of 0.25 H and 0.35 H would have the same total mass as the curtain with ρ_{p0} = 66.67 ρ_{0} and ρ_{p0} = 93.33 ρ_{0}, respectively. Figure 8 shows that configurations with the same total mass have practically the same oscillation period.
The last parameter we modified is the magnetic field strength. We varied B_{0} so that the maximum Alfvén speed ranges from 13.3 c_{s0} to 22.1 c_{s0}. Terradas et al. (2015) analysed the dependence on β of the relaxation process for a suspended prominence. They found that for low β, the magnetic field of their configuration is able to keep the prominence suspended. However, for low magnetic field strength, the dense plasma falls down and finally forms a hedgerow or curtainshaped prominence. These downward flows due to Rayleigh–Taylor instabilities complicate the dynamics of the system. Nevertheless, we avoided these flows by anchoring the density enhancement to the computational box. We obtained that for a weak magnetic field, the magnetic tension can hardly maintain the structure so that the material oscillates with a larger amplitude and at a lower equilibrium point. Panel c of Fig. 8 shows that the period decreases with the magnetic field strength. These results agree with Terradas et al. (2013), who found in their numerical study that for vertical oscillations, the period increases when the total mass is increased and decreases for weaker magnetic field strength. They also found that the period decreases slightly with the total length of the magnetic field line that crosses the prominence centre.
We analysed the attenuation of the vertical oscillations for the different cases. From the v_{z} series we obtained that the damping per period t_{d}/P does not vary significantly with width and ranges from 4.4 to 5.2. However, we found that the attenuation depends strongly on ρ_{p0} and B_{0}. t_{d}/P decreases with ρ_{p0} from 7.9 to 3.7 and increases with B_{0} from 2.2 to 7.6. When instead of analysing the v_{z} series, we study the centroid, as has been previously explained, no important changes in periodicity or attenuation are observed. This means, as in Terradas et al. (2008), that the global mode is dominant everywhere inside the prominence and that the damping time is basically the same everywhere.
4.2. Longitudinal oscillations
When the dense plasma enters an almost stationary state, we introduce in the prominence body a velocity perturbation in the xdirection to trigger longitudinal oscillations. The perturbation is
We refer to longitudinal oscillations as the movement along the magnetic field lines. The spatial distribution of the velocity disturbance was defined in such a way that it produces a global motion of the prominence but maintains the filament foot fixed to the photosphere. The disturbance has a maximum speed v_{x0} = 0.05 c_{s0} located at the centre of the prominence at a height of z_{v0} = 1.3 H, and its shape fits the density profile (w_{vx} = 0.3 H, w_{vy} = 2.2 H, and w_{vz} = 1.1 H). In the top panel of Fig. 9 the evolution of the xcomponent of the velocity (v_{x}) at the xzslice passing through y = 0 at the new equilibrium instant t = 0 is plotted with some selected density contours and some magnetic field lines.
Fig. 9. Time evolution of v_{x} at the xzplane passing through y = 0 for the evolution of longitudinal oscillations. The 12 blue curves correspond to some selected magnetic field lines. The orange isocontours represent the density profile. The initial parameters are the same as those in Fig. 2. 

Open with DEXTER 
Initially, the prominence moves in the direction given by the disturbance. The displacement reaches a different maximum amplitude depending on the perturbation profile. The greatest displacement occurs in the centre of the structure where we located the maximum speed. After the oscillation reaches the maximum displacement, the backward movement restores the motion at a different instant for each height. The simulation shows that the restoring movement occurs earlier in the lower part of the structure than in the upper part. This fact causes the motion to describe a serpentine movement that can be seen in the middle and bottom panels of Fig. 9. Moreover, in the middle and bottom panel of Fig. 9 we observe that v_{x} displays different velocity stripes of opposite sign aligned with the magnetic field lines that cross the prominence body. These positive/negative bands indicate that the body oscillates out of phase. On the whole, the centroid oscillates with a period of P = 37 min, a damping time of t_{d} = 47.6 min, and a displacement of f_{0} = 0.093 H.
To study the dependence of the period on height and the driving mechanism of the longitudinal oscillations, we compared the results with the pendulum model. Based on Luna et al. (2016), we selected 11 magnetic field lines located at the central xzplane of our numerical domain. The position of the centre of mass along each magnetic line (s_{cm}) was calculated. The left foot position of the selected field lines ranges from x = −4 H to x = −3.65 H with an increase of 0.035 H. The height of the selected lines at the centre of the prominence varies approximately from z = 0.3 H to z = 1.4 H for the reference simulation. At these heights all the selected field lines are concave upwards. We calculated the oscillation periods from the time evolution of s_{cm} of each selected field line and represent them in Fig. 10 as a function of the radius of the field line curvature. To calculate R we ignored its changes during the evolution of the oscillations, and we averaged the different values of the radius along a distance of 0.25 H from the prominence centre. The three panels of Fig. 10 show that the numerical results for the reference simulation (filled circles line) agree with the theoretical model expressed by Eq. (1) (solid line).
Fig. 10. Scatter plots of the oscillation periods for longitudinal oscillations as a function of the radius of the field line curvature. The top panel represents the different cases as a function of w_{x}, the middle panel as a function of ρ_{p0}, and the bottom panel as a function of B_{0}. The filled circles correspond to the reference simulation. The solid line represents Eq. (1). 

Open with DEXTER 
In order to perform a parametric study, we introduced a v_{x} disturbance in each of the different stationary states obtained in Sect. 4.1, and analysed the periodicity of longitudinal oscillations as a function of w_{x}, ρ_{p0}, and B_{0}. In Fig. 10 we show the scatter plots of the oscillation periods as a function of R for the different events. Panels a, b, and c correspond to simulations of different w_{x}, ρ_{p0} and B_{0}, respectively. As we showed in Sect. 4.1, the structures with a thinner density distribution, a lower density contrast, and a greater magnetic field strength suffer a smaller fall of their centre of mass during the relaxation process. Consequently, the prominence mass deforms the magnetic field lines less strongly, so that the radius of each field line is longer. As a result, in Fig. 10 we observe that the dotted lines for w_{x} = 0.25 H (panel a), ρ_{p0} = 40 ρ_{0} (panel b), and B_{0} = 12.5 G (panel c) reach longer radii for field lines with the same footpoint. In panel a of Fig. 10 we show that the period increases with the prominence width. As a result, the computed results become different from those of the pendulum model. This result agrees with Zhang et al. (2013). They reported a numerical study of a parametric survey of longitudinal oscillations and found that P increases slightly with the width of the thread. The authors proposed that the dependence of the period on the prominence width arises because in shorter threads the gas pressure gradient increases, contributing to the restoring force and shortening the periodicity. In panel b of Fig. 10 we also show that the oscillation period increases with the density contrast of the structure. In agreement with this result, Terradas et al. (2013) found that higher mass prominences have longer periods than lowmass prominences. They proposed that the cause of the increase in periodicity is the decrement of the sound speed when the mass is increased. In contrast, in panel c of Fig. 10 the dispersion in the scatter plot is small, and this means that there is no dependence of the period on the magnetic field. In a lowbeta regime, longitudinal oscillations are associated with the slow oscillatory modes whose characteristic speed is essentially c_{s0}. The period therefore does not change significantly with B_{0}. In panels a and b of Fig. 10 we also show that the differences obtained with respect to the reference simulation and therefore with the pendulum model are greater for larger radii. These results indicate that in addition to the curvature radius, the prominence width and density contrast also modify the period of the longitudinal oscillations. This means that the pendulum model is only a first approximation of longitudinal oscillations. This result agrees with Luna et al. (2012). They determined that for larger radii the pressure force contributes to the restoring force so that the pressuredriven term introduces a correction into the pendulum model. In addition, the discrepancies with the pendulum model could be due to the calculation of the radii of the magnetic field line curvature. R is not uniform along the streamline. The distance along the field lines we used to calculate the averaged radii can modify the results considerably. Furthermore, in this research, the variations in time that R suffers due to the residual oscillations of the relaxation process were ignored.
To estimate the damping time per period, we computed these two parameters from the evolution of the xcomponent of the centroid. We obtained that t_{d}/P decreases slightly with ρ_{p0} and B_{0}, from 1.7 to 1.2 and from 1.8 to 1.1, respectively. t_{d}/P does not vary significantly with w_{x}, whose values range between 1.2 and 1.5 (these results do not change if we calculate the damping from the v_{x} series). The values of t_{d}/P are much more lower for longitudinal oscillations than for the relaxation process. This significant damping in contrast to the transverse motion was also observed in the simulations of Luna et al. (2016), who associated the damping mechanism for longitudinal motions with a numerical viscosity that cancels the motion by stress. The variations of t_{d}/P along with ρ_{p0} and B_{0} could be understood by the correlation between the attenuation and the oscillation amplitude. For the same speed disturbance, the higher the density, the greater the inserted kinetic energy, and therefore the greater the oscillation amplitude. From the simulations we also obtained a slight increase in amplitude with magnetic force strength and so the variation of the attenuation, but the amplitude does not change with respect to w_{x}. Another mechanism that can explain the attenuation of longitudinal oscillations is the wave leakage. Zhang et al. (2019) found in a 2D nonadiabatic filament simulations that when the ratio of gravity to Lorentz force is close to unity, longitudinal oscillations deform the magnetic field, generating the subsequent transverse waves, which propagate away from the filament body. In our model we did not find significant deformation of the field lines and the system does not develop waves perpendicular to the magnetic lines (not shown here). This means that no wave leakage is observed in the simulations. This is probably due to the low β regime considered in this work.
4.3. Horizontal transverse oscillations
To trigger horizontal transverse oscillations, we induced a velocity disturbance pointing in the ydirection. The spatial distribution of the velocity perturbation is the same as that introduced in Sect. 4.2 for longitudinal oscillations, expressed as Eq. (9). In this case, the maximum speed located at the centre of the prominence at a height of z = 1.3 H is v_{y0} = 0.2 c_{s0}. In the top panel of Fig. 11 we show a yzslice of the ycomponent of the velocity (v_{y}) passing through the central plane at t = 0 together with some contours of the density profile. The plasma motion is somewhat complex. Initially, the inserted disturbance pushes the curtain, which is anchored at the base, tilting it towards the imposed direction. The displacement amplitude is greater in the internal parts of the filament where the velocity disturbance is stronger. Immediately, the plasma of the corona reacts to the movement in the opposite direction and tends to fill the displaced plasma. In addition, it is observed that the recovery movement starts at distinct times for different parts of the curtain. Regarding v_{y}, it is observed that during the first steps of the simulation the movement spreads as much throughout the prominence body as throughout the corona. As time progresses, the oscillations are more localised at the prominence and its PCTR but attenuate at the corona. In the middle and bottom panel of Fig. 11 we show that v_{y} develops different stripes of opposite speed. These bands are localised in the prominence body but also in the background. The distribution of the velocity traces horizontal bands at the corona, but these stripes bend when they approach the prominence and take the shape of the density contour. These clear shear motions indicate differences in periodicity with height. The numerical results show that from the evolution of the ycomponent of the centre of mass, the oscillation period is P = 9.41 min, the damping time is t_{d} = 11.9 min, and the displacement amplitude is f_{0} = 0.14 H. As we expected, P is shorter than that of longitudinal oscillations because in a lowβ regime slow modes have longer periods. We obtained that for transverse oscillations P is longer than that of the relaxation process, however. As was mentioned, P varies with height, therefore we plot in Fig. 12 the oscillation period as a function of z. We computed P by fitting the temporary evolution of v_{y} for different heights with Eq. (8) in which the sine was changed by a cosine. We obtained that the periodicity decreases with height. This result was also found by Zhou et al. (2018).
Fig. 11. Evolution of v_{y} at the yzplane passing through x = 0 for horizontal transverse oscillations with the prominence density (orange isocontours). The initial parameters are the same as those in Fig. 2. 

Open with DEXTER 
Fig. 12. Scatter plot of the period as a function of height for horizontal transverse oscillations. The initial parameters are the same as those in Fig. 2. 

Open with DEXTER 
We again analysed the behaviour of horizontal transverse oscillations from the different events studied in Sect. 4.1. In Fig. 13 the scatter plots of P are represented as a function of w_{x}, ρ_{p0} and B_{0}. Panels a and b show that the periodicity increases with the prominence width and with the initial density contrast. However, panel c indicates that P decreases with the magnetic field strength. This behaviour coincides with the results observed in Fig. 8 for the relaxation process. In both cases the decrease with the strength of the magnetic field is strong and the increase with density and prominence width is weaker. The fast mode of transverse oscillations is essentially a magnetic wave that is driven by magnetic forces; for this reason, we expected to find a stronger dependence of the period on magnetic field strength.
Fig. 13. Same as Fig. 6, but for horizontal transverse oscillations. 

Open with DEXTER 
The restoring force for the horizontal transverse oscillations is the magnetic tension force instead of gravity or pressure gradient. The different driving mechanisms for longitudinal and transverse oscillations reveal the different periodicity of the two oscillation modes.
While for vertical oscillations we obtained that t_{d}/P varies with ρ_{p0} and B_{0}, the damping in horizontal transverse oscillations does not change considerably with the considered parameters. When we compute t_{d}/P from series of the ycomponent of the centroid, all events show values that range from 1.25 to 1.34. These values are higher (from 2.79 to 2.92) when the calculation is made from the series of v_{y} at the point x = 0, y = 0, z = 1.3 H. This means that we have different damping at different parts of the prominence. As was explained in Sect. 4.1, the attenuation mechanism of the transverse modes is resonant absorption. For horizontal oscillations the absorption would occur in the upper part of the curtain, in its PCTR. However, as we showed in Fig. 11, the shear motions produced as a consequence of the attenuation are displaced towards the inside of the prominence. This behaviour is studied in more detail in Sect. 4.4. The attenuation of transverse oscillations also depends on the PCTR. As we saw for the vertical oscillations, we obtained that the smoother the PCTR, the stronger the attenuation. We obtained that for n = 2, t_{d}/P = 0.7, n = 4 t_{d}/P = 1.27 and n = 6 t_{d}/P = 1.87, respectively.
4.4. Continuum modes
We have computed the Alfvén continuum modes of the MHD equations based on the expressions found in Goossens et al. (1985). They derived the continuous spectrum of a 2D magnetostatic equilibrium. Their study is invariant with respect to y, therefore we considered our structure as 2D slices in the xzplane. This is the main approximation because our prominence model is in 3D. The expression for the modes of the Alfvén continuum is written in Terradas et al. (2013, see their Eq. (15)) in terms of the distance along the magnetic field lines. According to this equation, we can numerically compute the Alfvén frequency for different field lines by solving the eigenfrequency problem. To solve the equation, we need to extract the coordinates of each line and the variations of the magnetic field strength along the field line for the equilibrium state. We selected diverse magnetic field lines whose footpoints range from −4.315 H to −3.65 H in the xdirection and from 0 to 3.5 H in the ydirection. The Alfvén spectrum is represented in Fig. 14 in terms of the oscillatory period as a function of the position of each line at the central plane. The computation of the local Alfvén frequency is useful for understanding the resonant absorption process because, essentially, the resonant surface is the one where this frequency coincides with the frequency of the global mode inferred from the timedependent simulations. In Fig. 14 we show the contour of the global mode for the relaxation process (dashed black line) and for the horizontal transverse oscillation (dashed green line). For vertical oscillations, the resonant surface is more external than in the horizontal case. Terradas et al. (2008) found that the resonant surface can be located inside the structure, as we have obtained for the horizontal case. In Sects. 4.1 and 4.3, according to the velocity distribution of Figs. 4 and 11, we have estimated where the energy conversion takes place: at the yedges of the prominence body for the relaxation case and at the sheared bands inside the prominence for the horizontal transverse oscillations. However, the location of the resonant absorption process is clearer when we study the kinetic energy of the system. The top panel of Fig. 15 shows that the resonance surface of the vertical motion (dashed black curve) matches in space the location where the vertical kinetic energy increases due to the energy conversion. Only one snapshot of E_{kin, z} at t = 8.4 min is represented, but the maximum of E_{kin, z} is always localised around this point until when the oscillations attenuate. In the bottom panel of Fig. 15 we show a snapshot at t = 32 min of the ycontribution of the kinetic energy for horizontal transverse oscillations. Again, the resonance layer almost coincides with the maximum energy of the system. The exact behaviour of 3D Alfvén resonances has been studied in Wright & Elsden (2016) and Elsden & Wright (2017). They found that in 3D an infinite number of possible resonant solutions exist within a resonant zone, where the energy is accumulated, on which a dominant contour stands out from different ridges. It is beyond the scope of the present work to calculate the 3D Alfvén resonances, and it is enough to mention that the 2D approach provides reasonable results.
Fig. 14. Alfvén continuum spectrum in terms of the periodicity as a function of the position of each field lines at x = 0. Dashed black and green curves are the isocontours corresponding to the periods of the vertical and transverse oscillations, respectively, inferred from the timedependent simulations. 

Open with DEXTER 
Fig. 15. Top panel: snapshot of the vertical kinetic energy and the density profile at t = 8.4 min with the vertical resonant surface (dashed black curve) obtained in Fig. 14. The orange density contours are the same as those in the bottom panel of Fig. 4. Bottom panel: snapshot of the ycomponent of the kinetic energy and the density profile at t = 32 min with the horizontal resonant surface (dashed green curve) obtained in Fig. 14. The orange density contours are the same as those in the bottom panel of Fig. 11. 

Open with DEXTER 
5. Discussion and conclusions
A 3D numerical study of transverse and longitudinal oscillations in filaments has been presented. A curtainshaped prominence is anchored to the base of the corona so that the system does not develop downwards flows due to Rayleigh–Taylor instability, which complicates the analysis of the results of the oscillatory motion. In addition to the fixed foot of the structure, the behaviour of the system is subject to closed boundary conditions, meaning that we have imposed linetying conditions at all sides of the computational box. This perfectly reflecting condition applied in the bottom plane is essential to support the structure. To study the attenuation mechanism for transverse oscillations, we analysed the velocity distribution of the motion, the dependency of the attenuation on the PCTR, and the continuum Alfvén modes. In agreement with our simulations, we proposed that the main damping mechanism is the resonant absorption. As the theoretical models predict, we obtained that the attenuation is stronger for wider PCTR. From the analysis of the Alfvén continuum modes, we obtained that the location of the resonant surface coincides with the concentration of the kinetic energy, confirming that the resonant absorption is the main mechanism for transverse oscillations. However, we obtained that t_{d}/P ranges from 2.2 to 7.9 for the relaxation process of different scenarios, values that in principle are higher than expected from the observations. This means that in addition to the resonant absorption, another type of damping mechanism must be considered in our model. On the other hand, for transverse and longitudinal horizontal oscillations we obtained a stronger attenuation. One reason that can explain this significant damping for longitudinal oscillations could be changes in periodicity with height. The plasma of two adjacent magnetic field lines with different curvature radius oscillates at a different frequency, so that some shear may brake the movement. This type of sheared viscosity mechanism can also contribute to the attenuation of horizontal transverse oscillations. The shear motions inside the prominence due to the resonant absorption can brake the oscillation and increase the attenuation. Regardless, our model is relatively simple to study the observed damping, and other mechanisms such as nonadiabatic processes, mass flows, or partial ionisation effects must be considered.
Regarding the periodicity, we reproduced oscillations with periods that range approximately from 4.2 min to 7.3 min for vertical oscillations, from 28.2 min to 42.3 min for longitudinal oscillations, and from 7.4 min to 12.0 min for horizontal transverse oscillations. These results are somewhat shorter than the observations. In Tripathi et al. (2009) the periods associated with vertical motions range from 15 min to 29 min and in Shen et al. (2014b) from 11 min to 22 min. In Luna et al. (2018) the periods for horizontal transverse oscillations range from 30 min to 80 min and for longitudinal events from 30 min to 110 min. We recall that in the field of observations, in contrast to this work, longitudinal oscillations are motions along the filament axis. Mackay et al. (2010) reported that the magnetic field rotates with respect to the filament axis by an angle in the range from 15 to 30°. Then, Luna et al. (2018) in their catalogue considered as longitudinal oscillations those whose velocity vector forms an angle with the filament axis smaller than 40°. This suggests that a certain error is made in classifying the oscillation polarisation. Despite the short results obtained in P, for longitudinal oscillations we performed a reference simulation that agrees with the pendulum model. However, for values of w_{x} and ρ_{po} different from those of the reference simulation, the results move away from the theoretical expression. This means that another mechanism may drive the longitudinal motions, such as pressure force. For example, we obtained that the periodicity of the oscillations increases with the prominence width independently of the polarisation direction. This is an interesting result because it can be compared with diverse observed events. Luna et al. (2018) assembled a valuable catalogue from which we can extract much information, especially from their statistics section. In agreement with our results, in Fig. 21 (panel d) of Luna et al. (2018) they show that P tends to increase with width for quiescent prominences. In addition, the catalogue shows that high values of t_{d}/P are reported, which means that our results are not as different as the observed results, just as we had thought. Moreover, strong damping of about t_{d}/P ∼ 1.25 as the values obtained for longitudinal oscillations is the rule. As we found for horizontal transverse and longitudinal oscillations, Luna et al. (2018) did not find a correlation between t_{d}/P and the filament dimensions.
Another interesting result that we can compare with observations is the variation of P with height. Hershaw et al. (2011) studied the dependence of the periodicity and the attenuation on height of two successive trains of transverse oscillations observed on 2005 July 30. They found that the periods range between 90 and 110 min and increase slightly with height. Conversely, Shen et al. (2014a) did not find differences of the period with the altitude in a transverse oscillation observed on 9 August 2012, so that they suggested that the prominence oscillated as a rigid body. These behaviours differ from our results and from those of Zhou et al. (2018). We obtained that the period increases with height for longitudinal oscillations and decreases for transverse ones. These discrepancies suggest that a deeper analysis of the variations of P with height is necessary. For example, Hershaw et al. (2011) reported long periods, and if we take our simulations into account together with the long periods observed in their study, we can suggest that the oscillations are longitudinal instead of transverse. In any case, knowing the oscillation polarisation with respect to the magnetic field is essential for understanding the nature of the observations.
We have exhaustively analysed vertical oscillations for different scenarios. The origin of these oscillations is the relaxation process of our prominence model. For the stationary state, we induced a v_{z} perturbation in the reference simulation to trigger vertical motions. The results (not shown in this work) reveal basically the same behaviour as the relaxation process. According to this, we could consider the vertical oscillations of the relaxation process as a winking filament if the velocity amplitude of the movement were sufficiently large. We did not take the nature of the disturbance into account here. For this reason, in a future work we will study the dependence of the oscillations varying the characteristics of the perturbation. We are interested in simulating oscillations triggered by a distant perturbation or by pressure pulses and changing the energy of the wave.
Acknowledgments
We acknowledge the support from grant AYA201785465P (MINECO/AEI/FEDER, UE) and the Conselleria d’Innovació, Recerca i Turisme del Govern Balear to IAC^{3}. We are grateful to the ISSI Team led by Manuel Luna “LargeAmplitude Oscillations as a Probe of Quiescent and Erupting Solar Prominences” for inviting us to be part of the Team. A. A. acknowledges the Spanish “Ministerio de Economía, Industria y Competitividad” for the ’Ayuda para contratos predoctorales’ grant BES2015075040.
References
 Arregui, I., Oliver, R., & Ballester, J. L. 2012, Liv. Rev. Sol. Phys., 9, 2 [Google Scholar]
 Arregui, I., Oliver, R., & Ballester, J. L. 2018, Liv. Rev. Sol. Phys., 15, 3 [Google Scholar]
 Arregui, I., Terradas, J., Oliver, R., & Ballester, J. L. 2008, ApJ, 682, L141 [NASA ADS] [CrossRef] [Google Scholar]
 Ballester, J. L. 2014, in Nature of Prominences and their Role in Space Weather, eds. B. Schmieder, J. M. Malherbe, & S. T. Wu, IAU Symp., 300, 30 [NASA ADS] [Google Scholar]
 Brandenburg, A. 2003, in Computational Aspects of Astrophysical MHD and Turbulence, eds. A. FerrizMas, & M. Núñez, 269 [Google Scholar]
 Díaz, A. J., Oliver, R., Erdélyi, R., & Ballester, J. L. 2001, A&A, 379, 1083 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Elsden, T., & Wright, A. N. 2017, J. Geophys. Res. (Space Physics), 122, 3247 [NASA ADS] [CrossRef] [Google Scholar]
 Engvold, O. 2015, in Solar Prominences, eds. J. C. Vial, & O. Engvold, Astrophys. Space Sci. Lib., 415, 31 [NASA ADS] [CrossRef] [Google Scholar]
 Eto, S., Isobe, H., Narukage, N., et al. 2002, PASJ, 54, 481 [NASA ADS] [CrossRef] [Google Scholar]
 Forteza, P., Oliver, R., Ballester, J. L., & Khodachenko, M. L. 2007, A&A, 461, 731 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Gilbert, H. R., Daou, A. G., Young, D., Tripathi, D., & Alexander, D. 2008, ApJ, 685, 629 [NASA ADS] [CrossRef] [Google Scholar]
 Goossens, M., Andries, J., & Aschwanden, M. J. 2002, in Solar Variability: From Core to Outer Frontiers, ed. A. Wilson, ESA Spec. Publ., 506, 629 [NASA ADS] [Google Scholar]
 Goossens, M., Poedts, S., & Hermans, D. 1985, Sol. Phys., 102, 51 [NASA ADS] [CrossRef] [Google Scholar]
 Haugen, N. E. L., & Brandenburg, A. 2004, Phys. Rev. E, 70, 026405 [NASA ADS] [CrossRef] [Google Scholar]
 Hershaw, J., Foullon, C., Nakariakov, V. M., & Verwichte, E. 2011, A&A, 531, A53 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Hyder, C. L. 1966, Z. Astrophys., 63, 78 [NASA ADS] [Google Scholar]
 Jing, J., Lee, J., Spirock, T. J., et al. 2003, ApJ, 584, L103 [NASA ADS] [CrossRef] [Google Scholar]
 Joarder, P. S., & Roberts, B. 1992a, A&A, 256, 264 [NASA ADS] [Google Scholar]
 Joarder, P. S., & Roberts, B. 1992b, A&A, 261, 625 [NASA ADS] [Google Scholar]
 Joarder, P. S., & Roberts, B. 1993, A&A, 277, 225 [NASA ADS] [Google Scholar]
 Kippenhahn, R., & Schlüter, A. 1957, Z. Astrophys., 43, 36 [Google Scholar]
 Kleczek, J., & Kuperus, M. 1969, Sol. Phys., 6, 72 [NASA ADS] [CrossRef] [Google Scholar]
 Kraśkiewicz, J., Murawski, K., Solov’ev, A., & Srivastava, A. K. 2016, Sol. Phys., 291, 429 [NASA ADS] [CrossRef] [Google Scholar]
 Li, T., & Zhang, J. 2012, ApJ, 760, L10 [NASA ADS] [CrossRef] [Google Scholar]
 Luna, M., & Karpen, J. 2012, ApJ, 750, L1 [NASA ADS] [CrossRef] [Google Scholar]
 Luna, M., Díaz, A. J., & Karpen, J. 2012, ApJ, 757, 98 [NASA ADS] [CrossRef] [Google Scholar]
 Luna, M., Knizhnik, K., Muglach, K., et al. 2014, ApJ, 785, 79 [NASA ADS] [CrossRef] [Google Scholar]
 Luna, M., Terradas, J., Khomenko, E., Collados, M., & de Vicente, A. 2016, ApJ, 817, 157 [NASA ADS] [CrossRef] [Google Scholar]
 Luna, M., Karpen, J., Ballester, J. L., et al. 2018, ApJS, 236, 35 [NASA ADS] [CrossRef] [Google Scholar]
 Mackay, D. H., Karpen, J. T., Ballester, J. L., Schmieder, B., & Aulanier, G. 2010, Space Sci. Rev., 151, 333 [NASA ADS] [CrossRef] [Google Scholar]
 Okamoto, T. J., Nakai, H., Keiyama, A., et al. 2004, ApJ, 608, 1124 [NASA ADS] [CrossRef] [Google Scholar]
 Oliver, R., & Ballester, J. L. 2002, Sol. Phys., 206, 45 [NASA ADS] [CrossRef] [Google Scholar]
 Oliver, R., Ballester, J. L., Hood, A. W., & Priest, E. R. 1993, ApJ, 409, 809 [NASA ADS] [CrossRef] [Google Scholar]
 Parenti, S. 2014, Liv. Rev. Sol. Phys., 11, 1 [Google Scholar]
 Ramsey, H. E., & Smith, S. F. 1966, AJ, 71, 197 [NASA ADS] [CrossRef] [Google Scholar]
 Schutgens, N. A. J., & Tóth, G. 1999, A&A, 345, 1038 [NASA ADS] [Google Scholar]
 Shen, Y., Liu, Y. D., Chen, P. F., & Ichimoto, K. 2014a, ApJ, 795, 130 [NASA ADS] [CrossRef] [Google Scholar]
 Shen, Y., Ichimoto, K., Ishii, T. T., et al. 2014b, ApJ, 786, 151 [NASA ADS] [CrossRef] [Google Scholar]
 Soler, R., Oliver, R., & Ballester, J. L. 2007, A&A, 471, 1023 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Soler, R., Oliver, R., & Ballester, J. L. 2008, ApJ, 684, 725 [NASA ADS] [CrossRef] [Google Scholar]
 Soler, R., Oliver, R., Ballester, J. L., & Goossens, M. 2009, ApJ, 695, L166 [NASA ADS] [CrossRef] [Google Scholar]
 TandbergHanssen, E. 1995, The Nature of Solar Prominences (Netherlands: Springer Science & Business Media), 199 [Google Scholar]
 Terradas, J., Oliver, R., & Ballester, J. L. 2001, A&A, 378, 635 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Terradas, J., Oliver, R., & Ballester, J. L. 2004, ApJ, 614, 435 [NASA ADS] [CrossRef] [Google Scholar]
 Terradas, J., Arregui, I., Oliver, R., et al. 2008, ApJ, 679, 1611 [NASA ADS] [CrossRef] [Google Scholar]
 Terradas, J., Soler, R., Díaz, A. J., Oliver, R., & Ballester, J. L. 2013, ApJ, 778, 49 [NASA ADS] [CrossRef] [Google Scholar]
 Terradas, J., Soler, R., Luna, M., Oliver, R., & Ballester, J. L. 2015, ApJ, 799, 94 [NASA ADS] [CrossRef] [Google Scholar]
 Terradas, J., Soler, R., Luna, M., et al. 2016, ApJ, 820, 125 [NASA ADS] [CrossRef] [Google Scholar]
 Tripathi, D., Isobe, H., & Jain, R. 2009, Space Sci. Rev., 149, 283 [NASA ADS] [CrossRef] [Google Scholar]
 Vršnak, B., Veronig, A. M., Thalmann, J. K., & Žic, T. 2007, A&A, 471, 295 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Wright, A. N., & Elsden, T. 2016, ApJ, 833, 230 [NASA ADS] [CrossRef] [Google Scholar]
 Zhang, Q. M., Chen, P. F., Xia, C., & Keppens, R. 2012, A&A, 542, A52 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Zhang, Q. M., Chen, P. F., Xia, C., Keppens, R., & Ji, H. S. 2013, A&A, 554, A124 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Zhou, Y.H., Xia, C., Keppens, R., Fang, C., & Chen, P. F. 2018, ApJ, 856, 179 [NASA ADS] [CrossRef] [Google Scholar]
 Zhang, L. Y., Fang, C., & Chen, P. F. 2019, ApJ, 884, 74 [NASA ADS] [CrossRef] [Google Scholar]
 Zirin, H. 1998, The Astrophysics of the Sun (Cambridge: Cambridge University Press) [Google Scholar]
All Figures
Fig. 1. Magnetic field lines in the xzplane. In this plot B_{0} = 10 G, k_{1} = π/2L, k_{2} = 3k_{1}, L = 5 H, z_{0} = −0.2 H, and l_{1}/k_{1} = 1. The ycomponent of the magnetic field is null. The maximum initial Alfvén speed is v_{A0} = 17.7 c_{s0}. 

Open with DEXTER  
In the text 
Fig. 2. Density profile together with some magnetic field lines at t = 0. In this plot n = 4, w_{x} = 0.3 H, w_{y} = 2.2 H, w_{z} = 2.6 H, and ρ_{p0} = 80 ρ_{0}. 

Open with DEXTER  
In the text 
Fig. 3. Density distribution of the final equilibrium state at t = 90 min. Top panel: density isocontours and magnetic field lines in the plane y = 0. Bottom panel: density isocontours in the plane x = 0. Dashed lines correspond to the initial reference setup, and solid lines show the final equilibrium state at t = 90 min. 

Open with DEXTER  
In the text 
Fig. 4. Time evolution of density, v_{z}, and magnetic field lines for the relaxation process. v_{z} is represented in the yzplane at x = 0 as a 2D slice. The initial configuration for this case is shown in Fig. 2. 

Open with DEXTER  
In the text 
Fig. 5. Time evolution of the zcomponent of the centre of mass for n = 2 (red), n = 4 (green), and n = 6 (blue). The density average was computed for a delimited region (x ∈ [−0.5 H, 0.5 H], y ∈ [−3 H, 3 H], and z ∈ [0, 3.5 H]). 

Open with DEXTER  
In the text 
Fig. 6. v_{z} as a function of time for the relaxation process at a fixed point (0, 0, and 1.3 H). Dots are the values obtained from the numerical simulation, and the continuous grey curve is the fitting result based on Eq. (8). 

Open with DEXTER  
In the text 
Fig. 7. v_{z} as a function of y (x = 0, z = 1.3 H) at t = 2.4 min (solid line), t = 8.4 min (dashed line), and t = 14.4 min (dotted line). 

Open with DEXTER  
In the text 
Fig. 8. Scatter plots of the oscillation period during the relaxation process as a function of w_{x} (panel a), ρ_{p0} (panel b), and B_{0} (panel c). The values of P were obtained from the fitted curve based on Eq. (8) for the evolution of v_{z} at a fixed point (0, 0, and 1.3 H). 

Open with DEXTER  
In the text 
Fig. 9. Time evolution of v_{x} at the xzplane passing through y = 0 for the evolution of longitudinal oscillations. The 12 blue curves correspond to some selected magnetic field lines. The orange isocontours represent the density profile. The initial parameters are the same as those in Fig. 2. 

Open with DEXTER  
In the text 
Fig. 10. Scatter plots of the oscillation periods for longitudinal oscillations as a function of the radius of the field line curvature. The top panel represents the different cases as a function of w_{x}, the middle panel as a function of ρ_{p0}, and the bottom panel as a function of B_{0}. The filled circles correspond to the reference simulation. The solid line represents Eq. (1). 

Open with DEXTER  
In the text 
Fig. 11. Evolution of v_{y} at the yzplane passing through x = 0 for horizontal transverse oscillations with the prominence density (orange isocontours). The initial parameters are the same as those in Fig. 2. 

Open with DEXTER  
In the text 
Fig. 12. Scatter plot of the period as a function of height for horizontal transverse oscillations. The initial parameters are the same as those in Fig. 2. 

Open with DEXTER  
In the text 
Fig. 13. Same as Fig. 6, but for horizontal transverse oscillations. 

Open with DEXTER  
In the text 
Fig. 14. Alfvén continuum spectrum in terms of the periodicity as a function of the position of each field lines at x = 0. Dashed black and green curves are the isocontours corresponding to the periods of the vertical and transverse oscillations, respectively, inferred from the timedependent simulations. 

Open with DEXTER  
In the text 
Fig. 15. Top panel: snapshot of the vertical kinetic energy and the density profile at t = 8.4 min with the vertical resonant surface (dashed black curve) obtained in Fig. 14. The orange density contours are the same as those in the bottom panel of Fig. 4. Bottom panel: snapshot of the ycomponent of the kinetic energy and the density profile at t = 32 min with the horizontal resonant surface (dashed green curve) obtained in Fig. 14. The orange density contours are the same as those in the bottom panel of Fig. 11. 

Open with DEXTER  
In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.