K. Murawski1 - M. Selwa 1 - L. Nocera2
1 - Institute of Physics, UMCS, ul. Radziszewskiego 10, 20-031 Lublin, Poland
2 - Theoretical Plasma Physics, IPCF - CNR, via Moruzzi 1, 56124 Pisa, Italy
Received 17 January 2005 / Accepted 8 February 2005
We consider an impulsive excitation of fast magnetosonic waves in a dense and highly magnetised curved loop that is embedded in a potential arcade. The results of our numerical simulations reveal that the period of the excited waves agrees with both the period of the fast kink mode in the arcade and with the observations by TRACE.
Key words: Sun: oscillations - Sun: corona - magnetohydrodynamics (MHD)
Oscillations of solar coronal loops and arcades have been studied for a few decades. Body MHD waves propagating throughout an arcade structure attracted much interest because the highly inhomogeneous magnetic field background of the arcade modifies the linear spectrum of the waves (Oliver et al. 1993). Also the bearings of explosive flare events on the arcade were the subject of numerical investigations. The flare-driven short period oscillations were held responsible peculiar microwave signatures from the Sun (Oliver et al. 1998).
It is, however, dense, luminous, localised loop structures that provide better opportunity for comparing theory and observation. Such structures have already been advocated as the sites for peculiar radio signatures from the Sun (e.g. Roberts et al. 1983). An even greater interest in them has recently been kindled by observations by TRACE (e.g. Nakariakov et al. 1999).
Loops with a sufficiently small aspect ratio can be treated as rectilinear ducts. Surface waves propagating along the duct's edges (aptly named ducted waves) were studied by Roberts et al. (1983, 1984) who showed that, in the assumption of cold plasma (zero ), these waves can be classified as sausage and kink waves and are akin to waves in seismology.
Zero is a good assumption for a highly magnetised solar corona. In this limit the effects of energy leakage from ducts was investigated by Murawski & Roberts (1993b). The effect of a finite is the generation of slow magnetosonic waves propagating in the duct; and these were studied by Tsiklauri et al. (2004) and by Selwa et al. (2004) who also analysed the effects of dissipation and radiative losses.
The nonlinear properties of kink and sausage modes were investigated by Murawski & Roberts (1993a), who found them prone to wave breaking. The nonlinear propagation of ducted waves was studied by Nakariakov et al. (1997), who developed a model based upon the nonlinear Schrödinger equation, and by Ruderman (1993), who used an asymptotic relaxation method.
Murawski et al. (2005) found that high-order standing fast kink mode oscillations take place within a curved loop embedded in an unrealistically high potential arcade. Waves in curved ducts were also dealt with by Smith et al. (1997), who found that high frequency kink waves are affected by curvature induced energy leakage. Van Doorsselaere et al. (2004) found that continuum damping acts on curved duct quasi-modes (virtual global modes) and drew a connection with the damping of toroidal Alfvén eigenmodes of tokamak devices.
In this short note we give an interpretation of some oscillations recently observed by TRACE (Wang & Solanki 2004) as the ducted kink oscillations of a dense, highly magnetised, curved loop of a small aspect ratio embedded in a potential arcade and subject to impulsive perturbation. The dispersion relation for linear fast waves in rectilinear ducts (Roberts 1983) provides a good approximation of the oscillation period given by the simulations. The numerous effects outlined above may explain the difference between the two periods, the most prominent effect being the inhomogeneity of the duct along its longitudinal coordinate.
Our paper is organised as follows: the physical model is described in Sect. 2; the numerical results are presented in Sect. 3; conclusions and comparison with the work by other authors are drawn in Sect. 4.
In the solar corona the sound speed attains small values in comparison to the Alfvén speed.
Therefore we describe coronal plasma by the ideal MHD equations in which gas pressure terms
are neglected in comparison to the magnetic terms. This "cold'' plasma approximation
leads to the equations
The scale parameters introduced above are shown in Table 1, and are chosen in such a way that our results agree with the observations (cf. Sect. 4).
|Figure 1: Mass density profiles within the loop according to Eq. (13). Parameters are given in Table 2.|
|Open with DEXTER|
Next we embed a loop in the arcade in such a way that its edges follow two specific magnetic field lines. We prescribe that the inner and outer field lines cross the base of the arcade at |x|=D-2d and |x|=D, respectively. D and d are given in Table 2: they uniquely specify the loop's length l, height h, and its width 2a at the top (cf. Table 2).
Table 1: Scale parameters.
Table 2: Parameters of the equilibrium.
Inside the loop we use the following density profile (cf. Fig. 1):
|Figure 2: a) Time-signatures of mass density profiles measured on the line x=0. b) The density at the loop's top in time. Parameters of the simulations are given in Tables 2 and 3.|
|Open with DEXTER|
Table 3: Parameters of the perturbation.
As this loop is perturbed by the pulse it raises up initially and subsequently exhibits oscillations. These oscillations are discernible in Fig. 2a, which shows , the profile of the mass density across the loop's apex as it develops in time. The development of mass density at the location xd=0, zd=0.9L is shown in Fig. 2. The spectra of the oscillations of Fig. 2 are displayed in Fig. 3. The peak in the wavelet spectrum agrees with the peak in the frequency spectrum and, using the values in Table 1, the main peaks in the spectra correspond to a wave period of i.e. to a frequency of . The wavelet spectrum reveals a characteristic tadpole structure which shows that the instantaneous period of oscillation decreases with time. This behaviour is akin to that of fast magnetosonic waves in coronal loops (Nakariakov et al. 1993).
|Figure 3: Power spectra of the signals of Fig. 2. a) Wavelet spectrum. b) Frequency spectrum: the continuous and dashed lines are drawn using the FFT and Maximum method respectively. Using the values in Table 1 the main peaks in the spectra correspond to a wave period of , i.e. to a frequency of .|
|Open with DEXTER|
To prove this conjecture we approximate the thin loop of Sect. 2.1 by a rectilinear slab of width 2d and extension l.
We denote the Alfvén speeds inside and outside the loop by
being the ratio of the density inside the loop to the density outside the loop.
The dispersion relation for the fast mode kink oscillation of the slab is
(Roberts et al. 1983, Eq. (2))
The most likely explanation for the mismatch is the inhomogeneity of the Alfvén speed in the arcade. Indeed the value of the Alfvén speed at the apex of the loop ( , cf. Fig. 1) is about of its value at the base (cf. Eq. (12)). If the loop is approximated by a rectilinear duct, then both and vary along the duct. Furthermore the width of the loop also varies with height. These effects were not accounted for in the derivation of the dispersion relation of Eq. (15).
We may conclude that the oscillations observed in our numerical simulations are compatible with kink oscillations of a curved magnetic arcade. A similar conclusion was reached by Murawski et al. (2005), who also report the onset of the n=3 kink oscillations, although for a different arcade model and for a finite value of the plasma .
A match for the period of the oscillations found in the simulations with the value of given by TRACE (Wang & Solanki 2004) is attained if . For an Alfvén speed of we have . According to Table 2 for the loop's height and length we have and , which agrees with the observed values (Wang & Solanki 2004).
The authors express their cordial thanks to Prof. Sami Solanki, Dr. Valery Nakariakov, and Dr. Jaume Terradas for their warm and encouraging comments. M.S.'s & K.M.'s work wasfinancially supported by the grant from the State Committee for Scientific Research Republic of Poland, with KBN grant No. 2 PO3D 016 25. LN is grateful to KBN for support during his visit to KFT-UMCS, Lublin. The magnetohydrodynamics code used in this study was authored by C. Richard DeVore and provided through the sponsorship of NASA's HPCC/ESS program. Wavelet software is disseminated by C. Torrence and G. Compo and is available at URL http://paos.colorado.edu/research/wavelets.