Next Article in Journal
Jean Cleymans, Stringy Thermal Model, Tsallis Quantum Statistics
Next Article in Special Issue
Alfvén Wave Conversion and Reflection in the Solar Chromosphere and Transition Region
Previous Article in Journal
Stochastic Variational Method for Viscous Hydrodynamics
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Existence of Purely Alvén Waves in Magnetic Flux Tubes with Arbitrary Cross-Sections

by
Michael S. Ruderman
1,2,3,* and
Nikolai S. Petrukhin
4
1
Solar Physics and Space Plasma Research Centre (SP2RC), School of Mathematics and Statistics, University of Sheffield, Hicks Building, Hounsfield Road, Sheffield S3 7RH, UK
2
Space Research Institute (IKI), Russian Academy of Sciences, Profsoyuznaya Street, 84/32, Moscow 117997, Russia
3
Moscow Center for Fundamental and Applied Mathematics, Lomonosov Moscow State University, Gubkina Street, 8, Moscow 119991, Russia
4
Higher School of Economics, Myasnitskaya Ulitsa, 20, Moscow 101000, Russia
*
Author to whom correspondence should be addressed.
Physics 2022, 4(3), 865-872; https://doi.org/10.3390/physics4030055
Submission received: 16 June 2022 / Revised: 19 July 2022 / Accepted: 20 July 2022 / Published: 29 July 2022

Abstract

:
We study linear torsional Alfvén waves in a magnetic flux tube with an arbitrary cross-section. We assume that the equilibrium magnetic field is propagating in the z-direction in Cartesian coordinates x , y , and z. The tube cross-section is bounded by a smooth closed curve. Both plasma and magnetic field are homogeneous outside this curve. The magnetic field magnitude is a function of x and y, while the density is a product of two functions: one dependent on z and the other dependent on x and y. As a result, the Alfvén speed is also equal to V 0 ( x , y ) times a function of z. We define Alfvén waves as waves that do not disturb plasma density. We show that these waves can exist only when the magnetic field magnitude is a function of V 0 . When the condition of existence of Alfvén waves is satisfied, the waves are polarised in the directions tangent to the level lines of V 0 ( x , y ) and orthogonal to the equilibrium magnetic field. We found that the Alfvén wave amplitude has a specific form that depends on a particular coordinate system.

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, B . In Cartesian coordinates x , y , z , this field is in the z-direction with the magnitude dependent on x and y. It is given by B = B ( x , y ) e z , where e z is the unit vector in the z-direction. The equilibrium magnetic field and pressure are related by the condition,
p + B 2 2 μ 0 = const ,
where μ 0 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, ρ ( x , y , z ) = κ ( z ) ρ 0 ( x , y ) . Below, κ ( 0 ) = 1 is assumed. The plasma motion is governed by the set of linear ideal magnetohydrodynamic (MHD) equations:
ρ ˜ t + · ( ρ v ) = 0 , · b = 0 ,
ρ v t = p ˜ + 1 μ 0 b · B + 1 μ 0 ( B · ) b + 1 μ 0 ( b · ) B ,
b t = ( B · ) v ( v · ) B B ( · v ) ,
t p ˜ + γ p ρ ρ ˜ + ρ γ v · p ρ γ = 0 .
Here, ρ ˜ , p ˜ , and b are the perturbations of the density, pressure, and magnetic field, v 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, ρ ˜ = p ˜ = 0 in Equations (1)–(5). The assumption that the magnetic pressure is not perturbed implies that b is perpendicular to the z-axis. In addition, we assume that v is also perpendicular to the z-axis. Then, Equations (2)–(5) reduce to
· ( ρ 0 v ) = 0 , · ( B v ) = 0 , · b = 0 , b · B = 0 , v · p ρ 0 γ = 0 ,
ρ v t = B μ 0 b z , b t = B v z .
Multiplying the first equation in Equation (6) by B and subtracting the second equation in Equation (6) multiplied by ρ 0 , one obtains:
v · ( B ρ 0 ρ 0 B ) = 0 .
Using Equation (1), the last equation in Equation (6) transforms to
v · ( ρ 0 B B + μ 0 γ p ρ 0 ) = 0 .
From Equations (8), (9), and the first equation in Equation (6), it follwos that
v · B = 0 , v · ρ 0 = 0 , · v = 0 .
From the last but one equation in Equation (6) and the first equation in Equation (10), it follows that both v and b are orthogonal to B . This implies that v b , and both v and b are tangent to the level lines of B ( x , y ) . It follows from the second equation in Equation (10) that v is tangent to the level lines of ρ 0 ( x , y ) . 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 ρ 0 . 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 B = const . However, in this case, the last but one equation in Equation (6) is satisfied identically and the previously given proof that v b is not valid. However, from the second Equation (10), one can still obtain that v is tangent to the level lines of ρ 0 . Then, as it follows from Equation (7), b v . Let us emphasise that in general, we do not assume that B = const ; the only assumption is that B is a function of ρ 0 .
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 C b (see Figure 1). Outside of the domain enclosed by C b , B is a constant, while ρ can depend on z, but is independent of x and y. We define the Alfvén speed V A = B ( μ 0 ρ ) 1 / 2 , which can be rewritten as V A = κ 1 / 2 ( z ) V 0 ( x , y ) , where V 0 = B ( μ 0 ρ 0 ) 1 / 2 . V 0 is also a function of ρ 0 (or vice versa). Below, it is assumed that the level lines of V 0 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 ψ ( C ) to any level line C in such a way that ψ = 0 at the central line and ψ ( C 1 ) < ψ ( C 2 ) if the level line C 1 occurs inside the level line C 2 ; see Figure 1. If the notation ψ b ψ ( C b ) is used, then, ψ [ 0 , ψ b ] .
Now, let us introduce curvilinear orthogonal coordinates in the x y 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,
ϕ x = χ ψ y , ϕ y = χ ψ x .
The function χ is defined by the compatibility condition for the system of equations in Equation (11):
ψ x χ x + ψ y χ y + χ 2 ψ x 2 + 2 ψ y 2 = 0 .
We make χ to be non-dimensional. Then, ϕ is also non-dimensional. The ϕ -coordinate lines are defined by the condition ψ = const and coincide with the level lines of V 0 , which implies that V 0 is independent of ϕ . The ψ -coordinate lines are defined by the condition ϕ = const . They are orthogonal to the level lines of V 0 . The Jacobian of transformation from coordinates x , y to ϕ , ψ is
J = ϕ x ϕ y ψ x ψ y = χ H 2 , H 2 = ψ x 2 + ψ y 2 .
We assume that χ 0 everywhere. Since χ is a solution to Equation (12) if χ is a solution, one can always have χ > 0 . From the definition of ψ , it follows that the gradient of ψ is orthogonal to the level lines of V 0 ( x , y ) and not equal to zero anywhere except the central line, that is, H 0 . Then, one concludes that J > 0 everywhere but on the central line, and the curvilinear coordinates are right-oriented. Since v and b are tangent to the ϕ -coordinate lines that coincide with the level lines of V 0 , both v and b have only ϕ -component. In accordance with this, the new variables u and w are introduced:
v = u H e ϕ , b = w H e ϕ , e ϕ = 1 H e x ψ y e y ψ x ,
where e x and e y are the unit vectors in the x and y-direction, respectively. Using Equations (13) and (14), one obtains:
· v = χ H 2 u ϕ ,
and a similar expression for · b . Hence, the equations · v = 0 and · b = 0 reduce to the condition that u and w are independent of ϕ . Substituting expressions for v and b , given by Equation (14), in Equation (7) yields:
u t = B ρ μ 0 w z , w t = B u z .
Let us recall that ρ and, consequently, V A , 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
x = a sinh τ cosh τ cos σ a , y = a sin σ cosh τ cos σ ,
where a is a constant with the dimension of length. The coordinate origin corresponds to τ . The tube boundary C b is defined by τ = τ b , and the tube interior is defined by τ > τ b . In accordance with this, ψ = 1 / τ . Hence, ψ b = 1 / τ b . The tube radius is a ( sinh τ b ) 1 . The distance from the coordinate origin, r, is defined by
r 2 = x 2 + y 2 = a 2 1 + cosh τ 2 sinh τ + cos σ cosh τ cos σ .
Let us define the Alfvén speed in the plane z = 0 as
V 0 = V + ( V b V ) sinh τ b sinh τ 2 ,
where 0 < V < V b . One can see that V 0 is a monotonically decreasing function of τ , which implies that it gets to a minimum at the coordinate origin. In addition, V 0 = V + ( V b V ) ( r / a ) 2 for r a what corresponds to τ 1 . Hence, V 0 is a continuously differentiable function in the region bounded by C b . Using Equation (17), one obtains:
f x = 1 a f τ ( cosh τ cos σ 1 ) 1 a f σ sinh τ sin σ , f y = 1 a f τ sinh τ sin σ + 1 a f σ ( cosh τ cos σ 1 ) ,
where f is a continuously differentiable function. Using Equation (20), the formula [20],
2 f = 1 a 2 ( cosh τ cos σ ) 2 2 f σ 2 + 2 f τ 2 ,
and taking into account that ψ is independent of σ , Equation (12) reads:
τ χ d ψ d τ = 0 .
Using Equation (19), one obtains that a particular solution to this equation is χ = τ 2 . With this result and Equation (20), one obtains from Equation (11):
ϕ τ ( cosh τ cos σ 1 ) + ϕ σ sinh τ sin σ = sinh τ sin σ , ϕ τ sinh τ sin σ ϕ σ ( cosh τ cos σ 1 ) = ( cosh τ cos σ 1 ) .
From the system (23), it follows that
ϕ τ = 0 , ϕ σ = 1 .
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):
H = cosh τ cos σ a τ 2 .
Finally, e ϕ = e σ . It is interesting to compare the maximum and minimum values of the speed at the level line τ = const . Since u is independent of σ , it follows from Equations (14) and (25) that the ratio of these quantities is
| v | max | v | min = cosh τ + 1 cosh τ 1 .
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 x , y , 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 V A : that is, it can be expressed as a product of a function dependent on z and V 0 ( x , y ) . It is assumed that the level lines of V 0 ( x , y ) 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 V 0 ( x , y ) . This coordinate is zero at the central line and has the following property: If a level line C 1 is inside a level line C 2 , then this coordinate is smaller on C 1 than on C 2 .
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 V 0 . 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 V 0 . The wave amplitude varies along the level lines of V 0 . 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.

Author Contributions

The authors M.S.R. and N.S.P. contributed equally to this work. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Data Availability Statement

Not applicable.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Arregui, I. Wave heating of the solar atmosphere. Phil. Trans. R. Soc. A 2015, 373, 20140261. [Google Scholar] [CrossRef] [PubMed]
  2. Van Doorsselaere, T.; Srivastava, A.K.; Antolin, P.; Magyar, N.; Farahani, S.V.; Tian, H.; Kolotkov, D.; Ofman, L.; Guo, M.; Arregui, I.; et al. Coronal heating by MHD waves. Space Sci. Rev. 2020, 216, 140. [Google Scholar] [CrossRef]
  3. Tomczyk, S.; McIntosh, S.W.; Keil, S.L.; Judge, P.G.; Schad, T.; Seeley, D.H.; Edmondson, J. Alfvén waves in the solar corona. Science 2007, 317, 1192–1196. [Google Scholar] [CrossRef] [PubMed]
  4. De Pontieu, B.; McIntosh, S.W.; Carlsson, M.; Hansteen, V.H.; Tarbell, T.D.; Schrijver, C.J.; Title, A.M.; Shine, R.A.; Tsuneta, S.; Katsukawa, Y.; et al. Chromospheric Alfvénic waves strong enough to power the solar wind. Science 2007, 318, 1574–1577. [Google Scholar] [CrossRef] [PubMed]
  5. Erdélyi, R.; Fedun, V. Are there Alfvén waves in the solar atmosphere? Science 2007, 318, 1572–1574. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Van Doorsselaere, T.; Nakariakov, V.M.; Verwichte, E. Detection of waves in the solar corona: Kink or Alfvén? Astrophys. J. Lett. 2008, 676, L73–L75. [Google Scholar] [CrossRef] [Green Version]
  7. Jess, D.B.; Mathioudakis, M.; Erdélyi, R.; Crockett, P.J.; Keenan, F.P.; Christian, D.J. Alfvén waves in the lower solar atmosphere. Science 2009, 323, 1582–1585. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. De Pontieu, B.; Carlsson, M.; van der Voort, L.H.M.R.; Rutten, R.J.; Hansteen, V.H.; Watanabe, H. Ubiquitous torsional motions in type II spicules. Astrophys. J. Lett. 2012, 752, L12. [Google Scholar] [CrossRef] [Green Version]
  9. Liu, J.; Nelson, C.J.; Snow, B.; Wang, Y.; Erdélyi, R. Evidence of ubiquitous Alfvén pulses transporting energy from the photosphere to the upper chromosphere. Nature Com. 2019, 10, 3504. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  10. Aschwanden, M.J.; Wang, T.J. Torsional Alfvénic oscillations discovered in the magnetic free energy during solar flares. Astrophys. J. 2020, 891, 99. [Google Scholar] [CrossRef] [Green Version]
  11. Kohutova, P.; Verwichte, E.; Froment, C. First direct observation of a torsional Alfvén oscillation at coronal heights. Astron. Astrophys. 2020, 633, L6. [Google Scholar] [CrossRef] [Green Version]
  12. Stangalini, M.; Erdélyi, R.; Boocock, C.; Tsiklauri, D.; Nelson, C.J.; Del Moro, D.; Berrilli, F.; Korsós, M.B. Torsional oscillations within a magnetic pore in the solar photosphere. Nat. Astron. 2021, 5, 691–696. [Google Scholar] [CrossRef]
  13. Alfvén, H. Existence of electromagnetic-hydrodynamic waves. Nature 1942, 150, 405–406. [Google Scholar] [CrossRef]
  14. Spruit, H.C. Propagation speed and acoustic damping of waves in magnetic flux tubes. Solar Phys. 1982, 75, 3–17. [Google Scholar] [CrossRef]
  15. Ruderman, M.S.; Berghmans, D.; Goossens, M.; Poedts, S. Direct excitation of resonant torsional Alfvén waves by footpoint motions. Astron. Astrophys. 1997, 320, 305–318. Available online: https://adsabs.harvard.edu/full/1997A%26A...320..305R (accessed on 21 July 2022).
  16. Ruderman, M.S. Coronal loop heating by torsional Alfvén waves directly driven by footpoint motions: Harmonic driving versus stochastic driving. Astrophys. J. 1999, 521, 851–858. [Google Scholar] [CrossRef] [Green Version]
  17. Zaqarashvili, T.V. Observation of coronal loop torsional oscillation. Astron. Astrophys. 2003, 399, L15–L18. [Google Scholar] [CrossRef] [Green Version]
  18. Soler, R.; Terradas, J.; Oliver, R.; Ballester, J.L. Propagation of torsional Alfvén waves from the photosphere to the corona: Reflection, transmission, and heating in expanding flux tubes. Astrophys. J. 2017, 840, 20. [Google Scholar] [CrossRef]
  19. Roberts, B. MHD Waves in the Solar Atmosphere; Cambridge University Press: Cambridge, UK, 2019. [Google Scholar] [CrossRef] [Green Version]
  20. Korn, G.A.; Korn, T.M. Mathematical Handbook for Scientists and Engineers; Dover Publications, Inc.: Mineola, NY, USA, 1961; Available online: https://webéducation.com/mathematical-handbook-for-scientists-and-engineers-definitions-theorems-and-formulas-for-reference-and-review-pdf/ (accessed on 21 July 2022).
Figure 1. Sketch of magnetic tube cross-section. The thick curve C b shows the tube boundary. The plasma density and magnetic field magnitude are constant beyond this curve. The thin curves C 1 and C 2 are the level lines of V 0 . See text for details.
Figure 1. Sketch of magnetic tube cross-section. The thick curve C b shows the tube boundary. The plasma density and magnetic field magnitude are constant beyond this curve. The thin curves C 1 and C 2 are the level lines of V 0 . See text for details.
Physics 04 00055 g001
Figure 2. Curvilinear coordinates σ and τ . The thick curve shows the tube boundary C b . The non-concentric circles are the σ -coordinate lines defined by the equations τ = const . The non-circular lines are the τ -coordinate lines defined by the equations σ = const .
Figure 2. Curvilinear coordinates σ and τ . The thick curve shows the tube boundary C b . The non-concentric circles are the σ -coordinate lines defined by the equations τ = const . The non-circular lines are the τ -coordinate lines defined by the equations σ = const .
Physics 04 00055 g002
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ruderman, M.S.; Petrukhin, N.S. Existence of Purely Alvén Waves in Magnetic Flux Tubes with Arbitrary Cross-Sections. Physics 2022, 4, 865-872. https://doi.org/10.3390/physics4030055

AMA Style

Ruderman MS, Petrukhin NS. Existence of Purely Alvén Waves in Magnetic Flux Tubes with Arbitrary Cross-Sections. Physics. 2022; 4(3):865-872. https://doi.org/10.3390/physics4030055

Chicago/Turabian Style

Ruderman, Michael S., and Nikolai S. Petrukhin. 2022. "Existence of Purely Alvén Waves in Magnetic Flux Tubes with Arbitrary Cross-Sections" Physics 4, no. 3: 865-872. https://doi.org/10.3390/physics4030055

APA Style

Ruderman, M. S., & Petrukhin, N. S. (2022). Existence of Purely Alvén Waves in Magnetic Flux Tubes with Arbitrary Cross-Sections. Physics, 4(3), 865-872. https://doi.org/10.3390/physics4030055

Article Metrics

Back to TopTop