1. Introduction
One of the most intriguing problems in solar physics is the problem of coronal heating. It is still not completely clear what mechanisms are responsible for the increase of temperature in the solar atmosphere from about 5000 K in the photosphere to more than one million Kelvin in the corona. One popular mechanism of coronal heating is wave heating [
1,
2]. Any theory of wave heating has first to explain how to deposit the wave energy from the lower part of the solar atmosphere to the corona and second how to dissipate it there. The Alfvén waves are rather promising candidates for transporting wave energy to the corona.
The observation of Alfvén waves in the solar atmosphere is a rather complicated task because these waves do not perturb the plasma density. Still, at present, there are quite convincing observational results proving that Alfvén waves are ubiquitously present in the solar atmosphere. Probably, the first results of observation of Alfvén waves were reported in [
3,
4]. However, these results were questioned in [
5,
6]; there, it was argued that the results, reported in [
3,
4], are the observational signatures of not Alfvén waves but kink waves in magnetic flux tubes. Later, more robust observational evidence of Alfvén wave presence in the solar atmosphere was reported in [
7]. After that, other observations of Alfvén waves were presented. In [
8], observations of torsional Alfvén waves in type II spicules were reported. In [
9], the observations showing that Alfvénic pulses, propagating from the photosphere to the upper chromosphere are generated by swirls in the photosphere, were described. In [
10], the results of observations of torsional Alfvén oscillations in solar flares were presented. The first direct observational evidence of a fully resolved torsional Alfvén oscillation in a large-scale structure occurring at coronal heights was reported in [
11]. In [
12], the detection of torsional oscillations observed in a magnetic pore in the photosphere is described.
The theory of Alfvén waves was initiated in [
13], where the existence of the new type of waves was predicted. Then, it became evident that the solar atmosphere is highly inhomogeneous, Alfvén waves were studied in a magnetically structured atmosphere. One of the most popular structures modeling real magnetic structures in the solar atmosphere is an isolated magnetic tube. An axisymmetric magnetic flux tube can support torsional Alfvén waves. These waves have been studied in various publications, see, e.g., [
14,
15,
16,
17,
18]; see also the book by Roberts [
19]. In particular, in [
17], a method for observing torsional Alfvén waves was suggested. It was also shown that the observation of torsional Alfvén waves can be used for determining the Alfvén speed in the solar corona.
The torsional Alfvén waves were first studied in a straight homogeneous magnetic flux tube. Later, this study was extended to an expanding magnetic flux tube with the plasma density varying both across and along the tube. In all these studies, axisymmetric equilibria were considered. However, real flux tubes in the solar atmosphere could be not axisymmetric. Then, a question arises as to whether purely Alfvén waves can exist in such magnetic tubes. We aim to address this problem.
The paper is organised as follows.
Section 2 describes the equilibrium configuration and presents the governing equations. In
Section 3, the propagation of Alfvén waves in a magnetic flux tube with an arbitrary cross-section is considered. In
Section 4, one particular example of Alfvén waves in a non-axisymmetric equilibrium is studied.
Section 5 contains the summary of the results and conclusions.
2. Equilibrium and Governing Equations
Here, we consider a static equilibrium with the unidirectional magnetic field,
. In Cartesian coordinates
, this field is in the
z-direction with the magnitude dependent on
x and
y. It is given by
, where
is the unit vector in the
z-direction. The equilibrium magnetic field and pressure are related by the condition,
where
is the magnetic permeability of free space, and the gravity is neglected. This equation implies that
p also depends on
x and
y. What concerns equilibrium density,
, it can depend on all three Cartesian coordinates: that is, it can vary along the tube. However, in this case, the density must have the factorised form,
. Below,
is assumed. The plasma motion is governed by the set of linear ideal magnetohydrodynamic (MHD) equations:
Here, , , and are the perturbations of the density, pressure, and magnetic field, is the velocity, and is the ratio of specific heats.
3. Alfvén Waves in Ideal MHD
We define Alfvén waves as the waves that do not perturb the density, plasma pressure, and magnetic pressure. Hence,
in Equations (
1)–(
5). The assumption that the magnetic pressure is not perturbed implies that
is perpendicular to the
z-axis. In addition, we assume that
is also perpendicular to the
z-axis. Then, Equations (
2)–(
5) reduce to
Multiplying the first equation in Equation (
6) by
B and subtracting the second equation in Equation (
6) multiplied by
, one obtains:
Using Equation (
1), the last equation in Equation (
6) transforms to
From Equations (
8), (
9), and the first equation in Equation (
6), it follwos that
From the last but one equation in Equation (
6) and the first equation in Equation (
10), it follows that both
and
are orthogonal to
. This implies that
, and both
and
are tangent to the level lines of
. It follows from the second equation in Equation (
10) that
is tangent to the level lines of
. Then, one concludes that Alfvén waves can exist in the equilibrium as it is considered in this paper only when
B is a function of
. This condition is well satisfied in the solar corona where the magnetic pressure strongly dominates the plasma pressure. Consequently, from Equation (
1), it follows that
. However, in this case, the last but one equation in Equation (
6) is satisfied identically and the previously given proof that
is not valid. However, from the second Equation (
10), one can still obtain that
is tangent to the level lines of
. Then, as it follows from Equation (
7),
. Let us emphasise that in general, we do not assume that
; the only assumption is that
B is a function of
.
We see that the level lines of the equilibrium density play an important role in our analysis. Below, the density is assumed only to vary in a magnetic tube bounded by a simple smooth closed curve
(see
Figure 1). Outside of the domain enclosed by
,
B is a constant, while
can depend on
z, but is independent of
x and
y. We define the Alfvén speed
, which can be rewritten as
, where
.
is also a function of
(or vice versa). Below, it is assumed that the level lines of
constitute a nested set of simple smooth closed curves. This means that any level line is a smooth curve without intersection, and any two level lines do not intersect. The only exception is the central line. It can be either a dot or a not-closed line without intersection. When it represents a point, then we put the coordinate origin at this point; otherwise, the coordinate origin is put at one of the points along the central line. We prescribe a number
to any level line
in such a way that
at the central line and
if the level line
occurs inside the level line
; see
Figure 1. If the notation
is used, then,
.
Now, let us introduce curvilinear orthogonal coordinates in the
plane as follows. One of these coordinates is
. The second coordinate is
, defined by the condition that the system of curvilinear coordinates
and
is orthogonal,
The function
is defined by the compatibility condition for the system of equations in Equation (
11):
We make
to be non-dimensional. Then,
is also non-dimensional. The
-coordinate lines are defined by the condition
and coincide with the level lines of
, which implies that
is independent of
. The
-coordinate lines are defined by the condition
. They are orthogonal to the level lines of
. The Jacobian of transformation from coordinates
to
is
We assume that
everywhere. Since
is a solution to Equation (
12) if
is a solution, one can always have
. From the definition of
, it follows that the gradient of
is orthogonal to the level lines of
and not equal to zero anywhere except the central line, that is,
. Then, one concludes that
everywhere but on the central line, and the curvilinear coordinates are right-oriented. Since
and
are tangent to the
-coordinate lines that coincide with the level lines of
, both
and
have only
-component. In accordance with this, the new variables
u and
w are introduced:
where
and
are the unit vectors in the
x and
y-direction, respectively. Using Equations (
13) and (
14), one obtains:
and a similar expression for
. Hence, the equations
and
reduce to the condition that
u and
w are independent of
. Substituting expressions for
and
, given by Equation (
14), in Equation (
7) yields:
Let us recall that
and, consequently,
,
u, and
w can depend on
z. Equations (
16) are quite similar to those in an axisymmetric equilibrium. The only difference is that in an axisymmetric equilibrium,
u and
w are the
-components of the velocity and magnetic field perturbation, while in a non-axisymmetric equilibrium, the relation between the velocity and magnetic field perturbation and
u and
w is more involved and given by Equation (
14). Finally, to note is that the Alfvén waves studied here do not disturb the tube boundary.
4. Example: Tube with Circular Cross-Section and Non-Axisymmetric Alfvén Speed Distribution
In order to make more clear the general results, obtained in
Section 3, it is instructive to consider a particular example going step by step through the calculations made above. In this example, Alfvén waves are studied in a magnetic flux tube with a circular cross-section and non-axisymmetric Alfvén speed distribution. We use bi-cylindrical coordinates
and
(e.g., [
20]) (see
Figure 2). Cartesian coordinates are expressed in terms of curvilinear coordinates as
where
a is a constant with the dimension of length. The coordinate origin corresponds to
. The tube boundary
is defined by
, and the tube interior is defined by
. In accordance with this,
. Hence,
. The tube radius is
. The distance from the coordinate origin,
r, is defined by
Let us define the Alfvén speed in the plane
as
where
. One can see that
is a monotonically decreasing function of
, which implies that it gets to a minimum at the coordinate origin. In addition,
for
what corresponds to
. Hence,
is a continuously differentiable function in the region bounded by
. Using Equation (
17), one obtains:
where
f is a continuously differentiable function. Using Equation (
20), the formula [
20],
and taking into account that
is independent of
, Equation (
12) reads:
Using Equation (
19), one obtains that a particular solution to this equation is
. With this result and Equation (
20), one obtains from Equation (
11):
From the system (
23), it follows that
These equations imply that one can take
. Although we found the expressions for
and
, it is more convenient to use
and
in this particular example. Using Equations (
13) and (
20):
Finally,
. It is interesting to compare the maximum and minimum values of the speed at the level line
. Since
u is independent of
, it follows from Equations (
14) and (
25) that the ratio of these quantities is
One can see that this ratio increases with the distance from the coordinate origin and can be large far from it, what is for sufficiently small values of .
5. Summary and Conclusions
In this paper, torsional Alfvén waves were studied in a magnetic flux tube with an arbitrary cross-section. We assumed that the magnetic field is propagating in the z-direction in Cartesian coordinates and z. The tube cross-section is bounded by a smooth closed curve. Both the plasma pressure and the magnetic field are homogeneous outside this curve, while the plasma density still can vary in the z-direction. Inside the tube, all equilibrium quantities can depend on x and y. The density can also depend on z; however, it can be factorised and expressed as a product of two functions: one dependent on x and y and the other dependent on z. Hence, the same is valid for the Alfvén speed : that is, it can be expressed as a product of a function dependent on z and . It is assumed that the level lines of are smooth closed curves and constitute a nested set of curves. The only exception is the central curve that is either a point or a not closed line without self-intersections. We introduced an orthogonal curvilinear coordinate system with one of the two curvilinear coordinates being constant on the level lines of . This coordinate is zero at the central line and has the following property: If a level line is inside a level line , then this coordinate is smaller on than on .
Here, Alfvén waves were defined as waves that do not disturb plasma density. We showed that these waves can exist only when the magnetic field magnitude is a function of . We consider this as the main result of the paper. When the condition of existence of Alfvén waves is satisfied, the latter are polarised in the directions been orthogonal to the magnetic field and tangent to the level lines of . The wave amplitude varies along the level lines of . As an example, torsional Alfvén waves in a magnetic tube with a circular cross-section but non-axisymmetric distribution of the Alfvén speed were studied.