A&A 395, 669-675 (2002)
Q. Y. Luo - F. S. Wei - X. S. Feng
Key Laboratory for Space Weather, CAS, POB 8701, Beijing 100080, PR China
Received 26 March 2002 /Accepted 4 July 2002
The turbulent excitation of torsional Alfvén waves is presented on the basis of Lighthill-Stein theory. Due to the special properties of torsional modes, we can apply the theory to the inhomogeneous magnetic flux tube embedded in the solar photosphere, to evaluate the wave energy generated from the turbulent source in the convection zone. Torsional wave spectrum distribution along the r coordinate is achieved. In order to estimate the net wave energy transported to the chromosphere, we investigate the wave propagation and dissipation in the photospheric flux tube by simple phase mixing due to the inhomogeneity. Results for several cases with different model parameters are compared. Conclusions are drawn on the significant role of the energy carried by torsional Alfvén waves in the chromospheric and coronal heating, and the wave spectrum may be helpful in understanding the problem of footpoint azimuthal motion of coronal loops.
Key words: waves - Sun: photosphere - Sun: magnetic fields
Magnetic flux tubes are ubiquitous in space and in solar plasma, from small-scale structures like supergranulation, networks, coronal loops, and sunspots, to large-scale structures like prominences and coronal holes, all of which can be analyzed as tube topology. With more observations available from spacecrafts such as SOHO, TRACE, and ground stations, we get more data to understand the importance of magnetic flux tubes in the activities of the sun, e.g., explosive and transient phenomena. Hence, magnetic flux tubes are an essential subject in current space and solar physics. The dynamic motion of magnetic flux tubes was investigated in the convection zone and chromosphere in the classic paper by Spruit (1981), and the tube wave modes were discussed in Spruit (1982). More recently, Roberts & Ulmschneider (1997) reviewed extensively the research on the various dynamical processes and wave effects in magnetic flux tubes.
The mechanism of coronal heating has been an open question for many years, and the particular effects of magnetohydrodynamic (MHD) waves in the magnetic flux tubes are believed to play roles here. Many efforts have been made (De Groof et al. 2002; Sakai et al. 2000; Berghmans & De Bruyne 1995; Berghmans & Tirry 1997). The original work on the wave generation from turbulence, which concerns us most, was by Lighthill (1952), in which he studied the excitation of acoustic waves from fluid turbulence. Stein (1967) developed this wave excitation theory, and applied it to MHD wave generation. Furthermore, Collins (1989a,b, 1992) explored the general property of flux tube wave modes in the case of the period velocity field. Lighthill-Stein theory was also extended by Musielak et al. (1989), and his results highlighted the enhancement of MHD wave generation in magnetic flux tubes. More applications of the turbulent excitation theory had been applied to research on sunspots (Lee 1993), late-type dwarfs and subgiants (Musielak et al. 2000).
In the AC heating mechanism of chromosphere and corona, Alfvén waves are the most hopeful candidates to carry energy from the convection zone to the upper atmosphere. Therefore, Alfvén waves are the key element in the problem of MHD wave propagation in the solar atmosphere (Hollweg 1991). The dissipation mechanism of these waves have also been studied with many theories such as resonance absorption, shock dissipation, turbulent dissipation, nonlinear evolution, etc. After being first proposed by Heyvaerts & Priest (1983), phase mixing was accepted as one of the most important dissipation mechanisms of Alfvén waves, and was applied to explain the heating of coronal and solar winds. Ruderman et al. (1999) investigated the phase mixing effect in the two dimensionally axisymmetric magnetic flux tube with steady flow, and provided a possible description of the wave energy supply. More research is reviewed in Narain & Ulmschneider (1990, 1996).
In this paper, we first explore the problem of incompressible linear torsional Alfvén waves generated in the inhomogeneous magnetic flux tube, where the wave source is the turbulent motion located in the convection zone. To tackle the problem, we describe the whole flux tube as being divided into a bunch of "elementary tubes'' (see Sect. 2), where the plasma is uniform. In the process of wave excitation, the interference among the waves in different elementary tubes can be ignored under the assumption of a strong turbulent source model and the independence of the torsional mode. Then we can obtain the torsional wave energy flux in the elementary tube using Lighthill-Stein theory, and the waves are generated from the local turbulent source by the resonance excitation (unlike the work of Musielak, who studied the longitudinal and transverse modes). Due to the lack of knowledge about turbulence features in the convection zone, similar to previous work, the Kolmogorov turbulent spectrum is taken as the description of our source model. To estimate the net wave energy flux that can transfer through the photosphere, we evaluate the wave energy dissipation by taking account of the simple phase mixing effect. Finally, the magnitude of the torsional wave energy, which is available in the heating of chromosphere and corona, is given in our results.
In the next section, we give a brief introduction to the wave modes in magnetic flux tubes, and describe the inhomogeneous tube model in detail. In Sect. 3, we analyze the torsional Alfvén wave energy flux in the local elementary flux tube. Then in Sect. 4, to estimate the dissipation of the torsional waves, we employ the Heyvaerts and Priest result about phase mixing. Section 5 presents our results and relevant discussion and Sect. 6 our conclusions.
The basic modes of oscillation in magnetic flux tubes are determined by the property of the tube magnetic configuration. In terms of the geometrical topology, we distinguish them into sausage modes, kink modes, and fluting modes; in terms of the fluctuation pattern, they may be longitudinal modes, transverse modes, and torsional modes. For more details, see Roberts & Ulmschneider (1997).
Based on the Lighthill-Stein theory, the MHD wave energy flux generated from turbulent motion was studied for different modes: longitudinal modes in Musielak et al. (1989, 1995), transverse modes in Musielak & Ulmschneider (2001). Here, what we are concerned with is the dynamical process of excitation and dissipation of the linear incompressible torsional Alfvén mode in flux tubes.
We consider an individual untwisted magnetic flux tube, vertically
oriented, with a circular section and initial internal axial
(=B0(r)ez), and in
equilibrium with the surrounding medium. The magnetic field
in the tube declines with radial coordinate raccording to
|Figure 1: The elementary tube, as showed with the shaded area. R0 is the radius of the inhomogeneous magnetic flux tube, Bthe magnetic field declining with increasing radial coordinate, and the velocity perturbation.|
|Open with DEXTER|
In our problem, we neglect the solar atmosphere stratification due to gravitational force, and concentrate on the dynamical process of torsional modes, which is not so realistic but simplifies the problem considerably.
During the derivation of wave excitation, we assume that the
resistivity and viscosity of the magnetofluid can be ignored on
small spatial scales. Also, similar to the procedure of
Nakariakov et al. (1997) but in cylindrical coordinates, the MHD
equations of the
component for the inhomogeneous flux
tube are obtained,
Next, we will employ the MHD Eqs. (3), (4) in the local dynamical wave excitation. As showed in Fig. 1, an elementary portion can be separated from the inhomogeneous magnetic flux tube, and the variation of field parameters in this zone can be neglected. Here we decouple the torsional Alfvén wave from the magnetosonic waves by the constraint . And obviously, if we apply Eqs. (3), (4) to the local elementary portion, shown by the shaded area in Fig. 1, the component of perturbation can be replaced by the transverse x component in the "elementary tube'', as defined in this paper. In other words, the torsional modes in the whole magnetic flux tube can be considered as the sum of local waves generated from the local turbulence source. Furthermore, as the premise of the Lighthill-Stein theory, the strong turbulence pressure is assumed to be overwhelming in the nonlinear terms (see Musielak et al. 1995), i.e., the right hand sides of Eqs. (3), (4). In this sense, we can omit the effect from the interference of the waves among different elementary tubes during the process of wave excitation. Therefore, with these descriptions above, we assume that the efficiency of wave excitation can only be directly determined by the local hydromagnetic parameters and the turbulence strength. In what follows, we shall derive the generation of torsional Alfvén waves in an individual elementary tube first.
In this section, we deduce the energy flux of torsional Alfvén waves generated from the local turbulent source.
From the discussion in Sect. 2, if the
perturbations in Eqs. (3), (4) are replaced by the
x components in the elementary tube, we get the wave equation
The expression of wave energy flux in the elementary tube can be
derived from the general definition of MHD energy flux,
In the local elementary tube, we can Fourier-transform Eq. (5)
space (hereafter, we explicitly drop the
subscript z in kz),
During the above process of resonance excitation, we neglected the wave propagating antiparallel to the z direction. Then, we have the energy flux in Eq. (13), except for the expression of the source term, which is given in the next subsection.
Similar to Stein (1967), we take time and space average of the
terms, and substitute the arguments in Eq. (11) with
where we suppose there are no correlations between
longitudinal and transverse turbulence, and the fourth order
correlation can be reduced to the product of two second order
Thus, from our knowledge about the turbulence power spectrum,
which here satisfies the well-known Kolmogorov spectrum, we can
describe the one dimensional longitudinal and transverse
turbulence spectra in Eq. (18) as
With these expressions about the turbulent source term in Sect. 3.3, we can get the final expression of monochromatic wave energy
flux from Eq. (13) as
Then, we get the wave energy flux spectrum in the local elementary tube. Nevertheless, as mentioned above, not the total wave energy flux generated from the turbulence motion in the convection zone can propagate through the photosphere. If we want to estimate the torsional wave energy available in the heating of chromosphere and corona, we have to take the effect of phase mixing into account, which is the dominant dissipation mechanism of the torsional mode in our problem.
Propagation of MHD waves in the magnetofluid may produce some particular instabilities and dissipation owing to the inhomogeneity. Phase mixing is a kind of wave refraction in nature. When the phase velocity of the Alfvén wave depends on the spatial coordinates in the directions perpendicular to equilibrium magnetic field, large gradients of the wave amplitude in the perpendicular directions may be built up. As a consequence, it strongly increases the dissipation of Alfvén wave energy.
The phase mixing effect is a particular property of torsional Alfvén waves, and here we directly cite the result of phase mixing for propagating waves in Heyvaerts & Priest (1983) with only a slight modification.
The parameters involved are:
Finally, we can integrate Eq. (22) with respect to over the frequency region lower than to obtain the net wave energy flux, which is transferred up into chromosphere. In the same way, if we integrate it over the frequency region higher than , we then get the torsional Alfvén wave energy dissipated in the photosphere. The result of our evaluation is presented in the next section.
To evaluate the quantitative results of energy flux and the phase mixing effect, we give the necessary model parameters first. During the calculation, the integration of Eq. (22) with source depth is simplified to a factor z. We choose this value with the same order as the pressure scale height, here , and suppose the deeper sources (> ) with larger turbulent scale length are negligible in our problem. In our model, the mixing length coefficient is 1 (see Musielak & Ulmschneider 2001), and the pressure scale height is km. We choose the length of the magnetic flux tube as 500 km, which is the magnitude of the height of the photosphere, to estimate the possible energy carried by torsional waves up to the chromosphere after phase mixing.
The solar observation shows that there exists the horizontal velocity, 3 km s-1 at the surface of photosphere (Muller 1994), and even faster motion interacting with the magnetic flux tubes, which all can be sources to excite MHD waves. The rms values of these velocities are the turbulent velocities in our problem. Huang et al. (1995) suggested a rational range of turbulent velocity as km s-1. We choose 1.5 km s-1 in our estimation.
For comparison, we take some parameters similar to those of Musielak. The flux tube temperature 5770 K, external pressure dyn/cm2, and magnetic field at the tube center Gauss, all of which are general magnitudes for photospheric magnetic flux tubes. Using them in the formulas in previous sections, such as Eqs. (22) and (28), we can evaluate the results of wave energy flux and dissipation of torsional Alfvén waves.
Using Eq. (22) and the model parameters given in Sect. 5.1, we can obtain the torsional Alfvén wave energy flux spectrum given in Fig. 2 (where ), which shows that the energy flux density increases with radial coordinate r in the same tube. For example, the maximum of the wave energy flux increases from Hz to , while the respective r moves from the tube center to .
|Figure 2: Wave energy flux spectra in different r, denoted by with the magnetic radial decline coefficient and radius . The wave energy flux is the total torsional Alfvén wave energy flux generated from turbulent source.|
|Open with DEXTER|
Another aspect is that the profile of the spectrum tends to move towards the lower frequency direction when r increases. From Fig. 2, when increases from 0 to 0.8, the corresponding frequency, where the energy flux approaches its maximum, changes slightly from 7.08 Hz to 5.62 Hz. However, the respective phase mixing damping frequency , given by Eq. (28), decreases rapidly. When is 0, 0.4, and 0.8, the is 41.32, 3.64, and 0.43, respectively. Therefore, the strength of phase mixing at the exterior is much stronger than that at the interior of the flux tube.
The problem of how the fluid motion energy, carried by MHD waves from the convection zone, propagates through magnetic flux tubes up into the upper solar atmosphere has attracted solar physicists for many years. Musielak et al. (1989) investigated the interaction between thin magnetic flux tubes and the ambient turbulent medium, and worked out the longitudinal wave energy flux of about the order of 107 erg cm-2 s-1. Applying the Lighthill-Stein theory in the sunspot model, Lee (1993) evaluated the longitudinal and Alfvén wave energy flux to be about 108 erg cm-2 s-1. A larger energy flux given in Muller (1994) is about 1010 erg cm-2 s-1, but the difference may come from the different initial conditions. Based on the Spruit (1981) magnetic flux tube equations, Musielak & Ulmschneider (2001) estimated the transverse modes energy flux to be about the order of 108 erg cm-2 s-1. Progress was also made by numerical simulations, such as ny Huang et al. (1995), which simulated the nonlinear wave excitation response to the transverse perturbation at different heights relative to the magnetic flux tubes in the solar atmosphere, and obtained an energy flux density to be about the order of 109 erg cm-2 s-1.
Unlike these works mentioned above, our emphasis here is on the torsional modes generated from an inhomogeneous magnetic flux tube, and for a tube R0=500 km, we gain the average net energy flux density :
|Figure 3: Radial distribution of wave energy flux density in flux tubes with different radii, when . For , the energy flux density decreases near the edge of tube, due to the stronger enhancement of the phase mixing effect, which is far beyond the increase wave energy excitation.|
|Open with DEXTER|
|Figure 4: Radial distribution of wave energy flux density with different , when radius . For the stronger inhomogeneous magnetic field where , the energy flux density fades greatly owing to the stronger effects of phase mixing.|
|Open with DEXTER|
Figures 3 and 4 show that the exterior region of the tube, where r/R0>0.5, contributes largely to the total wave energy flux. The wave energy flux density increases rapidly with increasing coordinate r, which is due to the relaxation of strong magnetic field confinement in the tube center. It is similar to the stiffness effect on the flux tube wave excitation (see Musielak & Ulmschneider 2001).
As mentioned above, we neglect the stratification effect of the solar atmosphere, which will add a cutoff frequency into the wave equations (see Lee 1993), . The MHD waves with frequencies lower than can never transfer energy up to the chromosphere. Thus, the cutoff frequency may block the transfer of a part of the energy, which should be included in further research. Moreover, we note that the energy fluxes given here are the net wave energy fluxes, and those parts of energy dissipated by phase mixing are given below.
|Figure 5: Radial distribution of wave energy dissipation from phase mixing to different radii, , , and , when .|
|Open with DEXTER|
|Figure 6: Radial distribution of wave energy dissipation from phase mixing to different magnetic radial decline coefficient , when radius .|
|Open with DEXTER|
As a special property of torsional Alfvén waves, phase mixing should be taken into account when describing their propagation through an inhomogeneous magnetic flux tube. Here, we evaluate only the simple phase mixing effect, and obtain its heating efficiency for the flux tube.
For the flux tube with a magnetic radial decline coefficient , the average phase mixing dissipation is:
From the results given above, we find strong dependence of phase mixing on , which show the inhomogeneity of the tube. Comparing Fig. 4 with Fig. 6, when , phase mixing can be ignored; but when , for a tube with R0=500 km, the energy dissipation of phase mixing is overwhelming, and drags the net wave energy flux down fast when .
According to our results, thin flux tubes (here, R0=125 km) are affected more strongly by phase mixing, which can even deplete the wave energy of torsional modes. In the tubes with the same magnetic radial decline coefficient , e. g., for a flux tube with R0=2000 km, the dissipation is only about of the total wave energy generated from a turbulent source, while this value is about 14% for a tube with R0=500 km, and greater than for a tube with R0=125 km. But just like the consideration before, thin magnetic flux tubes are more prone to be distorted under external disturbance, so that its torsional modes couple with other modes, and the strong heating effect perhaps deviates the magnetofluid parameters greatly from their initial value in our models. Therefore, more careful consideration is necessary for this problem.
In summary, the results and main points of our work are listed as follows:
This work was supported by the National Natural Science Foundation of China under grants 49390450 & 499925412, and Chinese National Basic Research Project (Grant No. NkBRSF G20000984).