Source author record

J. Terradas

J. Terradas appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

33works
2topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

33 published item(s)

preprint2022arXiv

Extension and validation of the pendulum model for longitudinal solar prominence oscillations

Longitudinal oscillations in prominences are common phenomena on the Sun. These oscillations can be used to infer the geometry and intensity of the filament magnetic field. Previous theoretical studies of longitudinal oscillations made two simplifying assumptions: uniform gravity and semi-circular dips on the supporting flux tubes. However, the gravity is not uniform and realistic dips are not semi-circular. To understand the effects of including the nonuniform solar gravity on longitudinal oscillations, and explore the validity of the pendulum model with different flux-tube geometries. We first derive the equation describing the motion of the plasma along the flux tube including the effects of nonuniform gravity, yielding corrections to the original pendulum model. We also compute the full numerical solutions for the normal modes, and compare them with the new pendulum approximation. We have found that the nonuniform gravity introduces a significant modification in the pendulum model. We have also found a cut-off period, i.e. the longitudinal oscillations cannot have a period longer than 167 minutes. In addition, considering different tube geometries, the period depends almost exclusively on the radius of curvature at the bottom of the dip. We conclude that nonuniform gravity significantly modifies the pendulum model. These corrections are important for prominence seismology, because the inferred values of the radius of curvature and minimum magnetic-field strength differ substantially from those of the old model. However, we find that the corrected pendulum model is quite robust and is still valid for non-circular dips.

preprint2021arXiv

Geometrical properties of the interaction between oblique incoming coronal waves and coronal holes

Observations of coronal waves (CWs) interacting with coronal holes (CHs) show the formation of typical wave-like features, such as reflected, refracted and transmitted waves (collectively, secondary waves). In accordance with these observations, numerical evidence for the wave characteristics of CWs is given by simulations which demonstrate effects of deflection and reflection when a CW interacts with regions exhibiting a sudden density drop, such as CHs. However, secondary waves are usually weak in their signal and simulations are limited in the way the according idealisations have to be chosen. Hence, several properties of the secondary waves during a CW-CH interaction are unclear or ambiguous and might lead to misinterpretations. In this study we follow a theoretical approach and focus in particular on the geometrical properties of secondary waves caused by the interaction between oblique incoming CWs and CHs. Based on a linear theory, we derive analytical expressions for reflection and transmission coefficients, which tell us how strongly the amplitudes of the secondary waves increase and decrease with respect to the incoming wave, respectively. Additionally, we provide analytical terms for crucial incidence angles that are capable of giving information about the energy flux, the phase and the reflection properties of the secondary waves. These novel expressions provide a supplementary tool for estimating CW properties in a fast and straightforward way and, therefore, might have relevant consequences for a possible new interpretation of already studied CW-CH interaction events and the clarification of ambiguous observational data.

preprint2020arXiv

A new method for estimating global coronal wave properties from their interaction with solar coronal holes

Global coronal waves (CWs) and their interaction with coronal holes (CHs) result, among other effects, in the formation of reflected and transmitted waves. Observations of such events provide us with measurements of different CW parameters, such as phase speed and intensity amplitudes. However, several of these parameters are provided with only intermediate observational quality, other parameters, such as the phase speed of transmitted waves, can hardly be observed in general. We present a new method to estimate crucial CW parameters, such as density and phase speed of reflected as well as transmitted waves, Mach numbers and density values of the CH's interior, by using analytical expressions in combination with basic and most accessible observational measurements. The transmission and reflection coefficients are derived from linear theory and subsequently used to calculate estimations for phase speeds of incoming, reflected and transmitted waves. The obtained analytical expressions are validated by performing numerical simulations of CWs interacting with CHs. This new method enables to determine in a fast and straightforward way reliable CW and CH parameters from basic observational measurements which provides a powerful tool to better understand the observed interaction effects between CWs and CHs.

preprint2019arXiv

3D numerical simulations of oscillations in solar prominences

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. We numerically investigate the periodicity and damping of transverse and longitudinal oscillations in a 3D model of a curtain-shaped prominence. 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 (MHD) equations for a wide range of parameters, including the width and density of the prominence, and the magnetic field strength (B) of the solar corona. We studied the periodicity and attenuation of the induced oscillations. 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 density, 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.

preprint2016arXiv

The effects of magnetic-field geometry on longitudinal oscillations of solar prominences: Cross-sectional area variation for thin tubes

Solar prominences are subject to both field-aligned (longitudinal) and transverse oscillatory motions, as evidenced by an increasing number of observations. Large-amplitude longitudinal motions provide valuable information on the geometry of the filament-channel magnetic structure that supports the cool prominence plasma against gravity. Our pendulum model, in which the restoring force is the gravity projected along the dipped field lines of the magnetic structure, best explains these oscillations. However, several factors can influence the longitudinal oscillations, potentially invalidating the pendulum model. The aim of this work is to study the influence of large-scale variations in the magnetic field strength along the field lines, i.e., variations of the cross-sectional area along the flux tubes supporting prominence threads. We studied the normal modes of several flux tube configurations, using linear perturbation analysis, to assess the influence of different geometrical parameters on the oscillation properties. We found that the influence of the symmetric and asymmetric expansion factors on longitudinal oscillations is small.}{We conclude that the longitudinal oscillations are not significantly influenced by variations of the cross-section of the flux tubes, validating the pendulum model in this context.

preprint2015arXiv

Apparent cross-field superslow propagation of magnetohydrodynamic waves in solar plasmas

In this paper we show that the phase mixing of continuum Alfvén waves and/or continuum slow waves in magnetic structures of the solar atmosphere as, e.g., coronal arcades, can create the illusion of wave propagation across the magnetic field. This phenomenon could be erroneously interpreted as fast magnetosonic waves. The cross-field propagation due to phase mixing of continuum waves is apparent because there is no real propagation of energy across the magnetic surfaces. We investigate the continuous Alfvén and slow spectra in 2D Cartesian equilibrium models with a purely poloidal magnetic field. We show that apparent superslow propagation across the magnetic surfaces in solar coronal structures is a consequence of the existence of continuum Alfvén waves and continuum slow waves that naturally live on those structures and phase mix as time evolves. The apparent cross-field phase velocity is related to the spatial variation of the local Alfvén/slow frequency across the magnetic surfaces and is slower than the Alfvén/sound velocities for typical coronal conditions. Understanding the nature of the apparent cross-field propagation is important for the correct analysis of numerical simulations and the correct interpretation of observations.

preprint2015arXiv

Dynamics of coronal rain and descending plasma blobs in solar prominences: II. partially ionized case

Coronal rain clumps and prominence knots are dense condensations with chromospheric to transition region temperatures that fall down in the much hotter corona. Their typical speeds are in the range 30--150~km~s$^{-1}$ and of the order of 10--30~km~s$^{-1}$, respectively, i.e., they are considerably smaller than free fall velocities. These cold blobs contain a mixture of ionized and neutral material that must be dynamically coupled in order to fall together, as observed. We investigate this coupling by means of hydrodynamic simulations in which the coupling arises from the friction between ions and neutrals. The numerical simulations presented here are an extension of those of \citet{oliver2014} to the partially ionized case. We find that, although the relative drift speed between the two species is smaller than 1~m~s$^{-1}$ at the blob center, it is sufficient to produce the forces required to strongly couple charged particles and neutrals. The ionization degree has no discernible effect on the main results of our previous work for a fully ionized plasma: the condensation has an initial acceleration phase followed by a period with roughly constant velocity and, in addition, the maximum descending speed is clearly correlated with the ratio of initial blob to environment density.

preprint2015arXiv

On the robustness of the pendulum model for large-amplitude longitudinal oscillations in prominences

Large-amplitude longitudinal oscillations (LALOs) in prominences are spectacular manifestations of the solar activity. In such events nearby energetic disturbances induce periodic motions on filaments with displacements comparable to the size of the filaments themselves and with velocities larger than 20 km/s. The pendulum model, in which the gravity projected along a rigid magnetic field is the restoring force, was proposed to explain these events. However, it can be objected that in a realistic situation where the magnetic field reacts to the mass motion of the heavy prominence, the simplified pendulum model could be no longer valid. We have performed non-linear time-dependent numerical simulations of LALOs considering a dipped magnetic field line structure. In this work we demonstrate that for even relatively weak magnetic fields the pendulum model works very well. We therefore validate the pendulum model and show its robustness, with important implications for prominence seismology purposes. With this model it is possible to infer the geometry of the dipped field lines that support the prominence.

preprint2015arXiv

On the support of neutrals against gravity in solar prominences

Cool and dense prominences found in the solar atmosphere are known to be partially ionized because of their relative low temperature. In this Letter, we address the long-standing problem of how the neutral component of the plasma in prominences is supported against gravity. Using the multiple fluid approach we solve the time-dependent equations in two dimensions considering the frictional coupling between the neutral and ionized components of the magnetized plasma representative of a solar prominence embedded in a hot coronal environment. We demonstrate that given an initial density enhancement in the two fluids, representing the body of the prominence, the system is able to relax in the vicinity of magnetic dips to a stationary state in which both neutrals and ionized species are dynamically suspended above the photosphere. Two different coupling processes are considered in this study, collisions between ions and neutrals and charge exchange interactions. We find that for realistic conditions ions are essentially static while neutrals have a very small downflow velocity. The coupling between ions and neutrals is so strong at the prominence body that the behavior is similar to that of a single fluid with an effective density equal to the sum of the ion and neutral species. We also find that the charge exchange mechanism is about three times more efficient sustaining neutrals than elastic scattering of ions with neutrals.

preprint2015arXiv

Propagation and dispersion of sausage wave trains in magnetic flux tubes

A localized perturbation of a magnetic flux tube produces a pair of wave trains that propagate in opposite directions along the tube. These wave packets disperse as they propagate, where the extent of dispersion depends on the physical properties of the magnetic structure, on the length of the initial excitation, and on its nature (e.g., transverse or axisymmetric). In Oliver et al. (2014) we considered a transverse initial perturbation, whereas the temporal evolution of an axisymmetric one is examined here. In both papers we use a method based on Fourier integrals to solve the initial value problem. Previous studies on wave propagation in magnetic wave guides have emphasized that the wave train dispersion is influenced by the particular dependence of the group velocity on the longitudinal wavenumber. Here we also find that long initial perturbations result in low amplitude wave packets and that large values of the magnetic tube to environment density ratio yield longer wave trains. To test the detectability of propagating transverse or axisymmetric wave packets in magnetic tubes of the solar atmosphere (e.g., coronal loops, spicules, or prominence threads) a forward modelling of the perturbations must be carried out. This is left for a future work.

preprint2015arXiv

Solar prominences embedded in flux ropes: morphological features and dynamics from 3D MHD simulations

The temporal evolution of a solar prominence inserted in a three-dimensional magnetic flux rope is investigated numerically. Using the model of Titov Demoulin (1999) under the regime of weak twist, the cold and dense prominence counteracts gravity by modifying the initially force-free magnetic configuration. In some cases a quasi-stationary situation is achieved after the relaxation phase, characterized by the excitation of standing vertical oscillations. These oscillations show a strong attenuation with time produced by the mechanism of continuum damping due to the inhomogeneous transition between the prominence and solar corona. The characteristic period of the vertical oscillations does not depend strongly on the twist of the flux rope. Nonlinearity is the responsible for triggering the Kelvin-Helmholtz instability associated to the vertical oscillations and that eventually produces horizontal structures. Contrary to other configurations in which the longitudinal axis of the prominence is permeated by a perpendicular magnetic field, like in unsheared arcades, the orientation of the prominence along the flux rope axis prevents the development of Rayleigh-Taylor instabilities and therefore the appearance of vertical structuring along this axis.

preprint2014arXiv

Morphology and dynamics of solar prominences from 3D MHD simulations

In this paper we present a numerical study of the time evolution of solar prominences embedded in sheared magnetic arcades. The prominence is represented by a density enhancement in a background stratified atmosphere and is connected to the photosphere through the magnetic field. By solving the ideal magnetohydrodynamic (MHD) equations in three dimensions we study the dynamics for a range of parameters representative of real prominences. Depending on the parameters considered, we find prominences that are suspended above the photosphere, i.e., detached prominences, but also configurations resembling curtain or hedgerow prominences whose material continuously connects to the photosphere. The plasma$-β$ is an important parameter that determines the shape of the structure. In many cases magnetic Rayleigh-Taylor (MRT) instabilities and oscillatory phenomena develop. Fingers and plumes are generated, affecting the whole prominence body and producing vertical structures in an essentially horizontal magnetic field. However, magnetic shear is able to reduce or even to suppress this instability.

preprint2014arXiv

Propagation and dispersion of transverse wave trains in magnetic flux tubes

The dispersion of small amplitude, impulsively excited wave trains propagating along a magnetic flux tube is investigated. The initial disturbance is a localized transverse displacement of the tube that excites a fast kink wave packet. The spatial and temporal evolution of the perturbed variables (density, plasma displacement, velocity, ...) is given by an analytical expression containing an integral that is computed numerically. We find that the dispersion of fast kink wave trains is more important for shorter initial disturbances (i.e. more concentrated in the longitudinal direction) and for larger density ratios (i.e. for larger contrasts of the tube density with respect to the environment density). This type of excitation generates a wave train whose signature at a fixed position along a coronal loop is a short event (duration ~ 20 s) in which the velocity and density oscillate very rapidly with typical periods of the order of a few seconds. The oscillatory period is not constant but gradually declines during the course of this event. Peak values of the velocity are of the order of 10 km/s and are accompanied by maximum density variations of the order of 10-15% the unperturbed loop density.

preprint2014arXiv

Rayleigh-Taylor instabilities with sheared magnetic fields

Magnetic Rayleigh-Taylor (MRT) instabilities may play a relevant role in many astrophysical problems. In this work the effect of magnetic shear on the growth rate of the MRT instability is investigated. The eigenmodes of an interface and a slab model under the presence of gravity are analytically calculated assuming that the orientation of the magnetic field changes in the equilibrium, i.e., there is magnetic shear. We solve the linearised magnetohydrodynamic (MHD) equations in the incompressible regime.We find that the growth rate is bounded under the presence of magnetic shear. We have derived simple analytical expressions for the maximum growth rate, corresponding to the most unstable mode of the system. These expressions provide the explicit dependence of the growth rate on the various equilibrium parameters. For small angles the growth time is linearly proportional to the shear angle, and in this regime the single interface problem and the slab problem tend to the same result. On the contrary, in the limit of large angles and for the interface problem the growth time is essentially independent of the shear angle. In this regime we have also been able to calculate an approximate expression for the growth time for the slab configuration. Magnetic shear can have a strong effect on the growth rates of the instability. As an application of the results found in this paper we have indirectly determined the shear angle in solar prominence threads using their lifetimes and the estimation of the Alfvén speed of the structure.

preprint2013arXiv

Alfven waves in a partially ionized two-fluid plasma

Alfvén waves are a particular class of magnetohydrodynamic waves relevant in many astrophysical and laboratory plasmas. In partially ionized plasmas the dynamics of Alfvén waves is affected by the interaction between ionized and neutral species. Here we study Alfvén waves in a partially ionized plasma from the theoretical point of view using the two-fluid description. We consider that the plasma is composed of an ion-electron fluid and a neutral fluid, which interact by means of particle collisions. To keep our investigation as general as possible we take the neutral-ion collision frequency and the ionization degree as free parameters. First, we perform a normal mode analysis. We find the modification due to neutral-ion collisions of the wave frequencies and study the temporal and spatial attenuation of the waves. In addition, we discuss the presence of cut-off values of the wavelength that constrain the existence of oscillatory standing waves in weakly ionized plasmas. Later, we go beyond the normal mode approach and solve the initial-value problem in order to study the time-dependent evolution of the wave perturbations in the two fluids. An application to Alfvén waves in the low solar atmospheric plasma is performed and the implication of partial ionization for the energy flux is discussed.

preprint2013arXiv

Dynamics of coronal rain and descending plasma blobs in solar prominences: I. Fully ionised case

Observations of active regions and limb prominences often show cold, dense blobs descending with an acceleration smaller than that of free fall. The dynamics of these condensations falling in the solar corona is investigated in this paper using a simple fully ionised plasma model. We find that the presence of a heavy condensation gives rise to a dynamical rearrangement of the coronal pressure that results in the formation of a large pressure gradient that opposes gravity. Eventually this pressure gradient becomes so large that the blob acceleration vanishes or even points upwards. Then, the blob descent is characterised by an initial acceleration phase followed by an essentially constant velocity phase. These two stages can be identified in published time-distance diagrams of coronal rain events. Both the duration of the first stage and the velocity attained by the blob increase for larger values of the ratio of blob to coronal density, for larger blob mass, and for smaller coronal temperature. Dense blobs are characterised by a detectable density growth (up to 60% in our calculations) and by a steepening of the density in their lower part, that could lead to the formation of a shock. They also emit sound waves that could be detected as small intensity changes with periods of the order of 100 s and lasting between a few and about ten periods. Finally, the curvature of the falling path is only relevant when a very dense blob falls along inclined magnetic field lines.

preprint2013arXiv

MHD waves in two-dimensional prominences embedded in coronal arcades

Solar prominence models used so far in the analysis of MHD waves in such structures are quite elementary. In this work, we calculate numerically magnetohydrostatic models in two-dimensional configurations under the presence of gravity. Our interest is in models that connect the magnetic field to the photosphere and include an overlying arcade. The method used here is based on a relaxation process and requires solving the time-dependent nonlinear ideal MHD equations. Once a prominence model is obtained, we investigate the properties of MHD waves superimposed on the structure. We concentrate on motions purely two-dimensional neglecting propagation in the ignorable direction. We demonstrate how by using different numerical tools we can determine the period of oscillation of stable waves. We find that vertical oscillations, linked to fast MHD waves, are always stable and have periods in the 4-10 min range. Longitudinal oscillations, related to slow magnetoacoustic-gravity waves, have longer periods in the range of 28-40 min. These longitudinal oscillations are strongly influenced by the gravity force and become unstable for short magnetic arcades.

preprint2012arXiv

Analytic approximate seismology of propagating MHD waves in the solar corona

Observations show that propagating magnetohydrodynamic (MHD) waves are ubiquitous in the solar atmosphere. The technique of MHD seismology uses the wave observations combined with MHD wave theory to indirectly infer physical parameters of the solar atmospheric plasma and magnetic field. Here we present an analytical seismological inversion scheme for propagating MHD waves. This scheme uses in a consistent manner the observational information on wavelengths and damping lengths, along with observed values of periods or phase velocities, and is based on approximate asymptotic expressions for the theoretical values of wavelengths and damping lengths. The applicability of the inversion scheme is discussed and an example is given.

preprint2012arXiv

Prominence seismology

Given the difficulty in directly determining prominence physical parameters from observations, prominence seismology stands as an alternative method to probe the nature of these structures. We show recent examples of the application of magnetohydrodynamic (MHD) seismology techniques to infer physical parameters in prominence plasmas. They are based on the application of inversion techniques using observed periods, damping times, and plasma flow speeds of prominence thread oscillations. The contribution of Hinode to the subject has been of central importance. We show an example based on data obtained with Hinode's Solar Optical Telescope. Observations show an active region limb prominence, composed by a myriad of thin horizontal threads that flow following a path parallel to the photosphere and display synchronous vertical oscillations. The coexistence of waves and flows can be firmly established. By making use of an interpretation based on transverse MHD kink oscillations, a seismological analysis of this event is performed. It is shown that the combination of high quality Hinode observations and proper theoretical models allows flows and waves to become two useful characteristics for our understanding of the nature of solar prominences.

preprint2012arXiv

Surface Alfven waves in solar flux tubes

Magnetohydrodynamic (MHD) waves are ubiquitous in the solar atmosphere. Alfven waves and magneto-sonic waves are particular classes of MHD waves. These wave modes are clearly different and have pure properties in uniform plasmas of infinite extent only. Due to plasma non-uniformity MHD waves have mixed properties and cannot be classified as pure Alfven or magneto-sonic waves. However, vorticity is a quantity unequivocally related to Alfven waves as compression is for magneto-sonic waves. Here, we investigate MHD waves superimposed on a one-dimensional non-uniform straight cylinder with constant magnetic field. For a piecewise constant density profile we find that the fundamental radial modes of the non-axisymmetric waves have the same properties as surface Alfven waves at a true discontinuity in density. Contrary to the classic Alfveen waves in a uniform plasma of infinite extent, vorticity is zero everywhere except at the cylinder boundary. If the discontinuity in density is replaced with a continuous variation of density, vorticity is spread out over the whole interval with non-uniform density. The fundamental radial modes of the non-axisymmetric waves do not need compression to exist unlike the radial overtones. In thin magnetic cylinders the fundamental radial modes of the non-axisymmetric waves with phase velocities between the internal and the external Alfven velocities can be considered as surface Alfven waves. On the contrary, the radial overtones can be related to fast-like magneto-sonic modes.

preprint2012arXiv

The role of Rayleigh-Taylor instabilities in filament threads

Many solar filaments and prominences show short-lived horizontal threads lying parallel to the photosphere. In this work the possible link between Rayleigh-Taylor instabilities and thread lifetimes is investigated. This is done by calculating the eigenmodes of a thread modelled as a Cartesian slab under the presence of gravity. An analytical dispersion relation is derived using the incompressible assumption for the magnetohydrodynamic (MHD) perturbations. The system allows a mode that is always stable, independently of the value of the Alfvén speed in the thread. The character of this mode varies from being localised at the upper interface of the slab when the magnetic field is weak, to having a global nature and resembling the transverse kink mode when the magnetic field is strong. On the contrary, the slab model permits another mode that is unstable and localised at the lower interface when the magnetic field is weak. The growth rates of this mode can be very short, of the order of minutes for typical thread conditions. This Rayleigh-Taylor unstable mode becomes stable when the magnetic field is increased, and in the limit of strong magnetic field it is essentially a sausage magnetic mode. The gravity force might have a strong effect on the modes of oscillation of threads, depending on the value of the Alfvén speed. In the case of threads in quiescent filaments, where the Alfvén speed is presumably low, very short lifetimes are expected according to the slab model. In active region prominences, the stabilising effect of the magnetic tension might be enough to suppress the Rayleigh-Taylor instability for a wide range of wavelengths.

preprint2012arXiv

Transverse kink oscillations in the presence of twist

Magnetic twist is thought to play an important role in coronal loops. The effects of magnetic twist on stable magnetohydrodynamic (MHD) waves is poorly understood because they are seldom studied for relevant cases. The goal of this work is to study the fingerprints of magnetic twist on stable transverse kink oscillations. We numerically calculated the eigenmodes of propagating and standing MHD waves for a model of a loop with magnetic twist. The azimuthal component of the magnetic field was assumed to be small in comparison to the longitudinal component. We did not consider resonantly damped modes or kink instabilities in our analysis. For a nonconstant twist the frequencies of the MHD wave modes are split, which has important consequences for standing waves. This is different from the degenerated situation for equilibrium models with constant twist, which are characterised by an azimuthal component of the magnetic field that linearly increases with the radial coordinate. In the presence of twist standing kink solutions are characterised by a change in polarisation of the transverse displacement along the tube. For weak twist, and in the thin tube approximation, the frequency of standing modes is unaltered and the tube oscillates at the kink speed of the corresponding straight tube. The change in polarisation is linearly proportional to the degree of twist. This has implications with regard to observations of kink modes, since the detection of this variation in polarisation can be used as an indirect method to estimate the twist in oscillating loops.

preprint2012arXiv

Wave Leakage and Resonant Absorption in a Loop Embedded in a Coronal Arcade

We investigate the temporal evolution of impulsively generated perturbations in a potential coronal arcade with an embedded loop. As the initial configuration we consider a coronal loop, represented by a density enhancement, which is unbounded in the ignorable direction of the arcade. The linearized time-dependent magnetohydrodynamic equations have been numerically solved in field-aligned coordinates and the time evolution of the initial perturbations has been studied in the zero-beta approximation. For propagation constrained to the plane of the arcade, the considered initial perturbations do not excite trapped modes of the system. This weakness of the model is overcome by the inclusion of wave propagation in the ignorable direction. The inclusion of perpendicular propagation produces two main results. First, damping by wave leakage is less efficient because the loop is able to act as a wave trap of vertical oscillations. Second, the consideration of an inhomogeneous corona enables the resonant damping of vertical oscillations and the energy transfer from the interior of the loop to the external coronal medium.

preprint2011arXiv

Seismology of transversely oscillating coronal loops with siphon flows

There are ubiquitous flows observed in the solar atmosphere of sub-Alfvénic speeds, however after flaring and coronal mass ejection events flows can become Alfvénic. In this Letter, we derive an expression for the standing kink mode frequency due to siphon flow in coronal loops, valid for both low and high speed regimes. It is found that siphon flow introduces a linear spatially dependent phase shift along coronal loops and asymmetric eigenfunctions. We demonstrate how this theory can be used to determine the kink and flow speed of oscillating coronal loops with reference to an observational case study. It is shown that the presence of siphon flow can cause the underestimation of magnetic field strength in coronal loops using the traditional seismological methods.

preprint2010arXiv

On the linear coupling between fast and slow MHD waves due to line-tying effects

Oscillations in coronal loops are usually interpreted in terms of uncoupled magnetohydrodynamic (MHD) waves. Examples of these waves are standing transverse motions, interpreted as the kink MHD modes, and propagating slow modes, commonly reported at the loop footpoints. Here we study a simple system in which fast and slow MHD waves are coupled. The goal is to understand the fingerprints of the coupling when boundary conditions are imposed in the model. The reflection problem of a fast and slow MHD wave interacting with a rigid boundary, representing the line-tying effect of the photosphere, is analytically investigated. Both propagating and standing waves are analysed and the time-dependent problem of the excitation of these waves is considered. An obliquely incident fast MHD wave on the photosphere inevitably generates a slow mode. The frequency of the generated slow mode at the photosphere is exactly the same as the frequency of the incident fast MHD mode, but its wavelength is much smaller, assuming that the sound speed is smaller than the Alfvén speed. The main signatures of the generated slow wave are density fluctuations at the loop footpoints. We have derived a simple formula that relates the velocity amplitude of the transverse standing mode with the density enhancements at the footpoints due to the driven slow modes. Using these results it is shown that there are possible evidences in the observations of the coupling between these two modes.

preprint2010arXiv

Selective spatial damping of propagating kink waves due to resonant absorption

There is observational evidence of propagating kink waves driven by photospheric motions. These disturbances, interpreted as kink magnetohydrodynamic (MHD) waves are attenuated as they propagate upwards in the solar corona. In this paper we show that resonant absorption provides a simple explanation to the spatial damping of these waves. Kink MHD waves are studied using a cylindrical model of solar magnetic flux tubes which includes a non-uniform layer at the tube boundary. Assuming that the frequency is real and the longitudinal wavenumber complex, the damping length and damping per wavelength produced by resonant absorption are analytically calculated. The damping length of propagating kink waves due resonant absorption is a monotonically decreasing function of frequency. For kink waves with low frequencies the damping length is exactly inversely proportional to frequency and we denote this as the TGV relation. When moving to high frequencies the TGV relation continues to be an exceptionally good approximation of the actual dependency of the damping length on frequency. This dependency means that resonant absorption is selective as it favours low frequency waves and can efficiently remove high frequency waves from a broad band spectrum of kink waves. It is selective as the damping length is inversely proportional to frequency so that the damping becomes more severe with increasing frequency. This means that radial inhomogeneity can cause solar waveguides to be a natural low-pass filter for broadband disturbances. Hence kink wave trains travelling along, e.g., coronal loops, will have a greater proportion of the high frequency components dissipated lower down in the atmosphere. This could have important consequences with respect to the spatial distribution of wave heating in the solar atmosphere.

preprint2010arXiv

Three-Dimensional Propagation of Magnetohydrodynamic Waves in Solar Coronal Arcades

We numerically investigate the excitation and temporal evolution of oscillations in a two-dimensional coronal arcade by including the three-dimensional propagation of perturbations. The time evolution of impulsively generated perturbations is studied by solving the linear, ideal magnetohydrodynamic (MHD) equations in the zero-beta approximation. As we neglect gas pressure the slow mode is absent and therefore only coupled MHD fast and Alfven modes remain. Two types of numerical experiments are performed. First, the resonant wave energy transfer between a fast normal mode of the system and local Alfven waves is analyzed. It is seen how, because of resonant coupling, the fast wave with global character transfers its energy to Alfvenic oscillations localized around a particular magnetic surface within the arcade, thus producing the damping of the initial fast MHD mode. Second, the time evolution of a localized impulsive excitation, trying to mimic a nearby coronal disturbance, is considered. In this case, the generated fast wavefront leaves its energy on several magnetic surfaces within the arcade. The system is therefore able to trap energy in the form of Alfvenic oscillations, even in the absence of a density enhancement such as that of a coronal loop. These local oscillations are subsequently phase-mixed to smaller spatial scales. The amount of wave energy trapped by the system via wave energy conversion strongly depends on the wavelength of perturbations in the perpendicular direction, but is almost independent from the ratio of the magnetic to density scale heights.

preprint2010arXiv

Transverse oscillations of a multi-stranded loop

We investigate the transverse oscillations of a line-tied multi-stranded coronal loop composed of several parallel cylindrical strands. First, the collective fast normal modes of the loop are found with the T-matrix theory. There is a huge quantity of normal modes with very different frequencies and a complex structure of the associated magnetic pressure perturbation and velocity field. The modes can be classified as bottom, middle, and top according to their frequencies and spatial structure. Second, the temporal evolution of the velocity and magnetic pressure perturbation after an initial disturbance are analyzed. We find complex motions of the strands. The frequency analysis reveals that these motions are a combination of low and high frequency modes. The complexity of the strand motions produces a strong modulation of the whole tube movement. We conclude that the presumed internal fine structure of a loop influences its transverse oscillations and so its transverse dynamics cannot be properly described by those of an equivalent monolithic loop.

preprint2009arXiv

Excitation of standing kink oscillations in coronal loops

In this work we review the efforts that have been done to study the excitation of the standing fast kink body mode in coronal loops. We mainly focus on the time-dependent problem, which is appropriate to describe flare or CME induced kink oscillations. The analytical and numerical studies in slab and cylindrical loop geometries are reviewed. We discuss the results from very simple one-dimensional models to more realistic (but still simple) loop configurations. We emphasise how the results of the initial value problem complement the eigenmode calculations. The possible damping mechanisms of the kink oscillations are also discussed.

preprint2009arXiv

Kelvin-Helmholtz instability in coronal magnetic flux tubes due to azimuthal shear flows

Transverse oscillations of coronal loops are often observed and have been theoretically interpreted as kink magnetohydrodynamic (MHD) modes. Numerical simulations by Terradas et al. (2008, ApJ 687, L115) suggest that shear flows generated at the loop boundary during kink oscillations could give rise to a Kelvin-Helmholtz instability (KHI). Here, we investigate the linear stage of the KHI in a cylindrical magnetic flux tube in the presence of azimuthal shear motions. We consider the basic, linearized MHD equations in the beta = 0 approximation, and apply them to a straight and homogeneous cylindrical flux tube model embedded in a coronal environment. Azimuthal shear flows with a sharp jump of the velocity at the cylinder boundary are included in the model. We obtain an analytical expression for the dispersion relation of the unstable MHD modes supported by the configuration, and compute analytical approximations of the critical velocity shear and the KHI growth rate in the thin tube limit. A parametric study of the KHI growth rates is performed by numerically solving the full dispersion relation. We find that fluting-like modes can develop a KHI in time-scales comparable to the period of kink oscillations of the flux tube. The KHI growth rates increase with the value of the azimuthal wavenumber and decrease with the longitudinal wavenumber. However, the presence of a small azimuthal component of the magnetic field can suppress the KHI. Azimuthal motions related to kink oscillations of untwisted coronal loops may trigger a KHI, but this phenomenon has not been observed to date. We propose that the azimuthal component of the magnetic field is responsible for suppressing the KHI in a stable coronal loop. The required twist is small enough to prevent the development of the pinch instability.

preprint2009arXiv

On the nature of kink MHD waves in magnetic flux tubes

Magnetohydrodynamic (MHD) waves are often reported in the solar atmosphere and usually classified as slow, fast, or Alfvén. The possibility that these waves have mixed properties is often ignored. The goal of this work is to study and determine the nature of MHD kink waves. This is done by calculating the frequency, the damping rate and the eigenfunctions of MHD kink waves for three widely different MHD waves cases: a compressible pressure-less plasma, an incompressible plasma and a compressible plasma with non-zero plasma pressure which allows for MHD radiation. In all three cases the frequency and the damping rate are for practical purposes the same as they differ at most by terms proportional to $(k_z R)^2$. In the magnetic flux tube the kink waves are in all three cases, to a high degree of accuracy incompressible waves with negligible pressure perturbations and with mainly horizontal motions. The main restoring force of kink waves in the magnetised flux tube is the magnetic tension force. The total pressure gradient force cannot be neglected except when the frequency of the kink wave is equal or slightly differs from the local Alfvén frequency, i.e. in the resonant layer. Kink waves are very robust and do not care about the details of the MHD wave environment. The adjective fast is not the correct adjective to characterise kink waves. If an adjective is to be used it should be Alfvénic. However, it is better to realize that kink waves have mixed properties and cannot be put in one single box.

preprint2009arXiv

The effect of longitudinal flow on resonantly damped kink oscillations

The most promising mechanism acting towards damping the kink oscillations of coronal loops is resonant absorption. In this context most of previous studies neglected the effect of the obvious equilibrium flow along magnetic field lines. The flows are in general sub-Alfvénic and hence comparatively slow. Here we investigate the effect of an equilibrium flow on the resonant absorption of linear kink MHD waves in a cylindrical magnetic flux tube with the aim of determining the changes in the frequency of the forward and backward propagating waves and in the modification of the damping times due to the flow. A loop model with both the density and the longitudinal flow changing in the radial direction is considered. We use the thin tube thin boundary (TTTB) approximation in order to calculate the damping rates. The full resistive eigenvalue problem is also solved without assuming the TTTB approximation. Using the small ratio of flow and Alfvén speeds we derive simple analytical expressions to the damping rate. The analytical expressions are in good agreement with the resistive eigenmode calculations. Under typical coronal conditions the effect of the flow on the damped kink oscillations is small when the characteristic scale of the density layer is similar or smaller than the characteristic width of the velocity layer. However, in the opposite situation the damping rates can be significantly altered, specially for the backward propagating wave which is undamped while the forward wave is overdamped.

preprint2008arXiv

Transverse oscillations of systems of coronal loops

We study the collective kinklike normal modes of a system of several cylindrical loops using the T-matrix theory. Loops that have similar kink frequencies oscillate collectively with a frequency which is slightly different from that of the individual kink mode. On the other hand, if the kink frequency of a loop is different from that of the others, it oscillates individually with its own frequency. Since the individual kink frequency depends on the loop density but not on its radius for typical 1 MK coronal loops, a coupling between kink oscillations of neighboring loops take place when they have similar densities. The relevance of these results in the interpretation of the oscillations studied by \citet{schrijver2000} and \citet{verwichte2004}, in which transverse collective loop oscillations seem to be detected, is discussed. In the first case, two loops oscillating in antiphase are observed; interpreting this motion as a collective kink mode suggests that their densities are roughly equal. In the second case, there are almost three groups of tubes that oscillate with similar periods and therefore their dynamics can be collective, which again seems to indicate that the loops of each group share a similar density. All the other loops seem to oscillate individually and their densities can be different from the rest.