A&A 366, 306-310 (2001)
DOI: 10.1051/0004-6361:20000027
K. Murawski^{1} - V. M. Nakariakov^{2} - E. N. Pelinovsky^{3}
1 - Technical University of Lublin,
Department of Environmental Physics, ul. Nadbystrzycka 40,
20-618 Lublin, Poland
2 -
Physics Department, University of Warwick, Coventry CV4 7AL, UK
3 -
Institute of Applied Physics, 46 Ul'anova Street 603600,
Nizhny Novgorod, Russia
Received 9 June 2000 / Accepted 3 November 2000
Abstract
The propagation of fast magnetoacoustic waves in a randomly structured solar
corona is considered in the linear and cold plasma limits. The random field
is assumed to be static and associated with plasma density inhomogeneities
only. A transcendental dispersion relation for the fast magnetoacoustic waves
which propagate perpendicularly to the magnetic field is derived in the
weak random field approximation.
It is shown analytically that the fast magnetosonic waves experience
acceleration, attenuation, and dispersion in comparison to
the homogeneous case. These analytical findings are essentially
confirmed by numerical simulations for a wide-spectrum pulse, except that
the waves were found decelerated.
It is concluded that
the coronal Moreton waves can be applied to MHD seismology of the solar
corona.
Key words: oscillations MHD waves - methods: analytical - Sun: corona
MHD waves play an important role in the dynamics and formation of the solar corona as they are believed to be associated with the long-standing problem of the plasma heating and the solar wind acceleration. The EUV imaging telescopes onboard SOHO (EIT) and TRACE spacecrafts made it possible to directly observe the MHD waves in the corona. In particular, compressive longitudinal waves were detected in the polar plumes (Deforest & Gurman 1998; Ofman et al. 1999). These waves were interpreted as slow magnetoacoustic waves (Ofman et al. 1999; Ofman et al. 2000). Similar waves were detected in long coronal loops (Berghmans & Clette 1999).
Recently, Thompson et al. (1999) reported observations of the so-called coronal Moreton waves (or flare waves). These waves are seen as perturbations of the EUV emission intensity, propagating predominantly radially from a flare site. Most probably, these waves are fast magnetoacoustic waves because they are compressive and they propagate essentially isotropically in the low- coronal plasma. They are believed to be a manifestation of the waves which were discovered in the chromosphere in the H line (1961). These waves provide a unique tool for MHD seismology of the solar corona. Indeed, measuring parameters of these waves, e.g. the evolution of the wave amplitude and speed as a function of the distance from the epicenter, one can build up a map of the local Alfvén speed. However, application of this method requires a well-developed theory of MHD wave propagation in the corona.
Theoretical analysis of MHD waves in the corona shows that a key element of the coronal wave theory is interaction of the waves with coronal structures (e.g., Roberts 1991; Roberts & Ulmschneider 1997) which affect the MHD waves by a number of different ways: they can guide the waves due to refraction and reflection, resonantly absorb the wave energy when the wave phase speed coincides with local Alfvén or tube speeds, make the waves dispersive, lead to linear coupling of the waves and enhance efficiency of their nonlinear coupling. These and other aspects of the interaction of MHD waves with coronal structures are well-understood in the case of a deterministic structuring.
Recent interest in the effect of a random field on coronal MHD waves is arising in connection with the significant progress in the spatial and temporal resolution of coronal observations. For a wide band of large-scale wave motions, the coronal plasma is a randomly structured medium. Indeed, the typical scale of the inhomogeneity across the magnetic field is determined by the cross-section diameter of coronal loops, which randomly varies around a few Mm. This fact can dramatically affect the MHD wave dynamics. In particular, a random field leads to the attenuation of wave amplitude, dispersion, and to enhanced coupling of MHD modes. Valley (1971, 1974) derived dispersion relations for MHD waves propagating through plasma with a randomly structured magnetic field. He found that random irregularities in the magnetic field strength lead to attenuation of the waves. However, his choice of an analytically inconvenient correlation function did not allow him to estimate the effect quantitatively (see the discussion at the end of his Sect. V). Lou & Rosner (1986) reconsidered Valley's results for damping of incompressible MHD (pure Alfvén) waves in plasma with time-dependent fluctuations of the magnetic and velocity fields. Li & Zweibel (1987) developed the previous studies by considering time-dependent random irregularities. They found that as a consequence of the random field, an efficient conversion of Alfvén waves into compressive magnetoacoustic waves occurs.
This paper is organized as follows. The physical model is described in Sect. 2. The dispersion relation is derived and analyzed in Sect. 3. Numerical simulations for the wave equation are described in Sect. 4. Possible implications of these studies for the coronal Moreton waves are discussed in Sect. 5.
We consider a fast magnetoacoustic wave that propagates perpendicularly to the ambient magnetic field. As a consequence of this we can assume that the coronal plasma is cold. The fast magnetosonic waves are described then by cold ideal MHD equations, viz.
(3) |
The solar coronal plasma consists of randomly distributed coronal loops. Consequently, for the fast magnetoacoustic waves, the corona is a random environment with a characteristic spatial scale (the correlation length l_{z}) of a few Mm (the coronal loop cross-section diameter) and a density ratio (the density ratio inside and outside the loop).
In what follows we assume one-dimensional motions with
and V_{x}=V_{y}=0 and consider the
equilibrium state, in which the magnetic field
points in the
x-direction and is uniform. The equilibrium quantities are given by
(5) |
The standard linearization procedure of Eqs. (1-4) leads
to the following wave equation:
(7) |
Following Howe (1971), we use the expansion
(8) |
From the dispersion relation (12) it follows that the dependence of the cyclic frequency on the wave number k differs from the coherent dispersion relation, . A random field leads to the appearance of dispersion. This dispersion is described by the real part of . As a consequence of scattering by turbulent flow, the energy of the fast magnetoacoustic wave is partially transformed into the turbulent field. This phenomenon is associated with the imaginary part of the frequency, .
We introduce dimensionless variables
(14) |
In principle, one can
consider solutions of Eq. (16) in some limiting
cases, applying asymptotic expansions for the plasma dispersion function
(Fried & Conte 1961). However, in the general case, determination of
the dependency of the wave frequency
on the wave number K
requires numerical methods.
Figure 1: The difference between the real parts of the fast magnetoacoustic wave frequency for the random () and homogeneous () media (top panel) and imaginary part of (bottom panel) as functions of the wave number K. Different curves correspond to different variances: (solid), (broken), (dash-dotted), and (dash-dot-dot-dot). The cyclic frequency and the wavenumber K are shown in units of l_{z}/V_{A0} and l_{z}, respectively | |
Open with DEXTER |
Figure 1 (top panel) displays the difference between the real part of that is evaluated from Eq. (16) and , as a function of the wave number K. A more intensive random field (with a higher value of ) increases more significantly the frequency of the fast magnetoacoustic wave with respect to the homogeneous medium case. This increase is also higher for higher values of K. For low values of and K the influence of the random field is negligible. In the case of and the real part of the cyclic frequency departs significantly from for K>4 and K>2, respectively.
Figure 1 (bottom panel) illustrates the imaginary part of the frequency . It is noteworthy that always attains negative values which correspond to the wave damping. A higher damping occurs for higher values of which represents a stronger random field. This result confirms the qualitative conclusions of Valley (1971). It is interesting that the damping decrement has a minimum value for a certain value of the wave number K, prescribed by the variance . For a higher , the minimum of , which corresponds to the strongest damping, is shifted towards lower K.
To illustrate the theory developed above, we perform numerical simulations
of Eq. (6).
Attention is paid to evolution of a wide-spectrum pulse.
As different spectral components of the pulse propagate with different phase and
group speeds, due to a dispersion brought by the random field, the pulse shape evolves in time. Our aim is to determine a typical
profile of the small-amplitude pulse which interacts with randomly distributed
coronal structures. A typical profile of the Alfvén speed is shown in Fig.
2. The average Alfvén speed is unity.
Figure 2: A typical profile of the local Alfvén speed | |
Open with DEXTER |
The numerical simulations are performed with the use of the existing
and well-tested solver to the problem
(LeVeque 1997).
The code utilizes the Godunov-type method which
yields accurate results near steep profiles and complex flows.
The computational domain has a typical dimension
and is divided
into a grid of at least 20000 cells. The boundaries of the simulation region
are set at z=0 and z=200l_{z}. At the left and right boundaries free-streaming
conditions
are applied. Tests were performed with larger domains and a finer numerical
grid to ensure grid convergence and to be sure that the accuracy of the results was
not being compromised by the grid size.
Initially, at t=0, an impulse in the perturbed velocity V_{z} is launched
at the initial position z = 5l_{z}
to excite waves which propagate through the simulation region.
Figure 3: The profiles of the fast magnetoacoustic pulses at . The solid line corresponds to , dotted to and dashed to . The pulse amplitude is normalized to the initial amplitude at t=0 | |
Open with DEXTER |
Figure 4: The delay (top) and attenuation (bottom) of the fast magnetoacoustic pulse in randomly structured medium with respect to the non-structured case for different values of the parameter . The crosses correspond to particular realizations results. The ensemble averaged data is joined by the dashed lines | |
Open with DEXTER |
Figure 3 shows typical profiles of the initially Gaussian pulse at the time . In the case of almost uniform medium (the variance of the random field is negligibly small, ), the pulse approximately keeps its Gaussian form and initial amplitude (solid curve). In a randomly structured medium, the pulse is attenuated and decelerated (dotted and dashed curves). The tails appearing behind the pulses are typical manifestations of the effect of dispersion.
The deceleration and attenuation of the pulse are obviously proportional to the
variance of the random field. Figure 4 (top) shows the relative
delay of the pulse, defined as follows:
The attenuation of the pulse due to a random field is shown in Fig. 4 (bottom). The relative amplitude of the pulses is shown as a function of the variance at . The amplitude is normalized by the pulse amplitude in the case of a homogeneous medium . We see that a random field leads to very efficient attenuation of the fast magnetoacoustic pulse amplitude, about 50% for .
We considered a simple 1D model of a perpendicular fast magnetoacoustic wave interaction with randomly distributed coronal structures. Our findings can be summarized as follows: a) even a relatively weak random density field, with variance less than 0.5, leads to significant (up to 50%) attenuation of the fast wave amplitude. This effect is stronger for wavelengths shorter than the correlation length of the random field. b) Numerical results reveal that fast waves are decelerated by random field. However, this effect is less pronounced than attenuation. For , the deceleration can be about 25% for some Fourier components. A total deceleration of the wide-spectrum pulse is found to be for . c) A fast magnetoacoustic pulse experiences a dispersion that is introduced by the random field.
Our results are in qualitative agreement with the numerical findings by Shapiro & Kneib (1993) who have shown that a random pulse experiences broadening, decaying, and its spectrum is shifted towards lower frequencies as well as by Karweit & Blanc-Benon (1995) who have observed retardation of random waves.
These results suggest that the theory of interaction of MHD waves with randomly distributed coronal structures can be applied to seismology of the corona. Indeed, observing the damping and attenuation of coronal waves and applying the theory developed above, we obtain an additional tool for the determination of physical parameters of the solar corona, such as the coronal loop filling factor (associated with the correlation length l_{z} and the variance ).
There are several physical effects which have been left out of the model: in particular, coupling of the waves with slow and Alfvén modes and refraction. These neglected mechanisms can also lead to deceleration of the fast waves and the amplitude attenuation. However, consideration of these effects has to be at least 2D and so requires a much more developed technique. Also, any analytical treatment of such problems does not seem to be possible. From this point of view, this work can be considered as a first step toward the construction of the theory for wave interaction with a randomly structured corona.
One interesting application of the coronal MHD wave theory is the explanation of the coronal Moreton wave speed determined observationally as about 200-600 kms^{-1} (Thomson et al. 1999). As is discussed above, these waves are believed to be fast magnetoacoustic waves. In the low-coronal plasma, the fast magnetoacoustic wave speed has to be between the Alfvén speed V_{A0} and , depending upon the angle between the wave vector and the local magnetic field. Here, c is the sound speed. So, the wave speed has to be slightly higher than the Alfvén speed. But, according to accepted estimations, the coronal Alfvén speed is over 700-800 kms^{-1}, which is supported by recent observational findings (Nakariakov et al. 1999). As a consequence of that the coronal Moreton waves propagate with speeds significantly lower than the expected speed. This discrepancy is to be explained.
The numerical simulations presented above suggest that the fast magnetoacoustic waves are indeed decelerated in a randomly structured medium. These simulations provide a qualitative explanation of the Moreton wave speed. An extension of the simulations which allow the quantitative comparison with the observations will have to incorporate 2D effects such as a random orientation of the magnetic field along the path of the wave.
Acknowledgements
The authors express their sincere thanks to Dr. Luigi Nocera and to Prof. Bernard Roberts for the stimulating discussions. K. M.'s work was financially supported by the KBN grant No. 2 PO3D 017 17. The work of K. M. and V. M. N. was in the frame of the Nato grant PST.CLG.976135, while V. M. N. and E. N. P. acknowledge the INTAS grant 97 31931. A part of the numerical calculations have been performed on a SUN workstation at the Institute of Physics, UMCS Lublin.