Nonlinear Long-Wavelength Torsional Alfvén Waves
Total Page:16
File Type:pdf, Size:1020Kb
A&A 526, A80 (2011) Astronomy DOI: 10.1051/0004-6361/201016063 & c ESO 2010 Astrophysics Nonlinear long-wavelength torsional Alfvén waves S. Vasheghani Farahani1,V.M.Nakariakov1,2, T. Van Doorsselaere3, and E. Verwichte1 1 Centre for Fusion, Space and Astrophysics, Physics Department, University of Warwick, Coventry CV4 7AL, UK e-mail: [email protected] 2 Central Astronomical Observatory of the Russian Academy of Sciences at Pulkovo, 196140 St Petersburg, Russia 3 Centrum voor Plasma-Astrofysica, Mathematics Department, Celestijnenlaan 200B bus 2400, 3001 Leuven, Belgium Received 3 November 2010 / Accepted 23 November 2010 ABSTRACT Aims. We investigate the nonlinear phenomena accompanying long-wavelength torsional waves in solar and stellar coronae. Methods. The second order thin flux-tube approximation is used to determine perturbations of a straight untwisted and non-rotating magnetic flux-tube, nonlinearly induced by long-wavelength axisymmetric magnetohydrodynamic waves of small, but finite ampli- tude. Results. Propagating torsional waves induce compressible perturbations oscillating with double the frequency of the torsional waves. In contrast with plane shear Alfvén waves, the amplitude of compressible perturbations is independent of the plasma-β and is pro- portional to the torsional wave amplitude squared. Standing torsional waves induce compressible perturbations of two kinds, that grow with the characteristic time inversely proportional to the sound speed, and that oscillate at double the frequency of the inducing torsional wave. The growing density perturbation saturates at the level, inversely proportional to the sound speed. Key words. magnetohydrodynamics (MHD) – waves – Sun: corona – Sun: activity – magnetic fields 1. Introduction The compressible flows induced by nonlinear Alfvén waves have been intensively studied in the context of solar wind accel- Alfvén waves are often considered as the primary candidate for eration (e.g. Ofman & Davila 1998). A modified Cohen-Kulsrud the acceleration of solar (e.g. Cranmer 2009) and stellar (e.g. equation describing spherical Alfvén waves in a stratified atmo- Charbonneau & MacGregor 1995; Suzuki 2007) winds and coro- sphere with a radial magnetic field was derived in Nakariakov nal heating (e.g. Ofman 2005). Alfvén waves are also considered et al. (2000b) in application to coronal holes. Analytical results in the context of the collimation and confinement of astrophys- are consistent with numerical MHD modelling performed in ical jets (Bisnovatyi-Kogan 2007), and in core-collapse super- Torkelsson & Boynton (1998). More advanced one-dimensional nova explosions (Suzuki et al. 2008). However, observational models for nonlinear Alfvén waves in solar and stellar atmo- ´ spheres in open-field regions include upflows, super-radial mag- evidence of Alfven waves in astrophysical plasmas still remains ff indirect, e.g. as a possible interpretation of non-thermal broad- netic field geometry, non-adiabatic e ects and various dissipa- ening of coronal emission lines (e.g. Banerjee et al. 2009)in tion mechanisms (Suzuki 2004; Suzuki & Inutsuka 2005; Suzuki coronal holes. 2007, 2008). Secularly growing compressible perturbations, in- duced by standing Alfvén waves, have been considered in the One of the key ingredients of the theoretical modelling of context of coronal loop hydrostatics by Litwin & Rosner (1998), Alfvén waves in solar and stellar atmospheres is the concept and in connection with magnetospheric field-line resonances by of the nonlinear cascade. It is needed to explain the transfer Tikhonchuk et al. (1995). Secular compressible effects associ- of wave energy from the low-frequency injection range (that is ated with travelling in opposite directions Alfvén waves were around 10 mHz in the solar atmosphere) to the high- found in Verwichte et al. (1999). All those studies were carried frequencydissipation range. In one-dimensional models, the out in terms of a shear Alfvén wave, one-dimensional model. nonlinear cascade is connected with nonlinear generation of One-dimensional models mentioned above are based upon higher harmonics, and hence steepening of the waves, caus- the assumption that the waves are plane. For long-period Alfvén ing the onset of nonlinear dissipation or of non-MHD dissipa- waves, with periods of the order of typical time scales of lower- tive processes. In a uniform plasma, this process is analytically atmospheric dynamics, this condition is not fulfilled. For exam- described by a weakly-nonlinear evolutionary equation, known ple, for a period of 10 min and an Alfvén speed of 1 Mm/s, the as the Cohen–Kulsrud equation (Cohen & Kulsrud 1974). The longitudinal wavelength is 600 Mm. In a plane wave, the trans- cubically-nonlinear term in this equation accounts for the non- verse wavelength should be much larger than the longitudinal linear self-interaction of linearly or elliptically polarised, plane wavelength. Hence, for the generation of a plane Alfvén wave of Alfvén waves via the wave-induced perturbation of the local a 10-min period, the wave driver should be of the size exceeding Alfvén speed. These perturbations are often referred to as the the solar diameter. Also there should be no transverse structur- nonlinearly-induced compressible motions in Alfvén waves. In ing of the plasma in the Alfvén speed, otherwise the wavefront contrast with the parallel magnetoacoustic waves (e.g., slow is distorted, and the wave becomes non-planar (e.g. Botha et al. waves in the low-β plasma of the solar corona), these perturba- 2000). Thus, the study of the initial stage of the nonlinear cas- tions exist even in the zero-β regime. Circularly polarised plane cade in the corona requires consideration of non-planar Alfvén Alfvén waves are not subject to this effect. waves. Article published by EDP Sciences A80, page 1 of 4 A&A 526, A80 (2011) In non-uniform plasma structures, Alfvén waves are situ- force connected with the azimuthal rotation of the plasma, the ated on magnetic surfaces. In the field-aligned structured coro- magnetic tension force caused by the magnetic field curvature, nal plasmas, Alfvén waves can be present in the form of tor- and the ponderomotive force that is connected with the longitu- sional modes (see discussion in Van Doorsselaere et al. 2008a,b). dinal gradients of the magnetic pressure perturbation in the tor- Torsional waves are intensively studied theoretically in the con- sional wave. The first two forces are absent from the plane wave text of coronal heating (e.g. Antolin et al. 2008; Copil et al. 2008; theory of Alfvén waves and appear because of plasma structur- Antolin & Shibata 2010), coronal seismology (Zaqarashvili & ing. These forces can modify the flux-tube cross-sectional area, Murawski 2007; Verth et al. 2010) and particle acceleration in hence inducing compressible plasma motions in the longitudi- solar flares (Fletcher & Hudson 2008). There is some indirect nal and radial directions. In the second-order thin flux-tube ap- evidence of torsional standing modes (Zaqarashvili 2003)and proximation, these effects are taken into account in the trans- propagating waves (Banerjee et al. 2009) in spectroscopic data, verse force-balance equation. The ponderomotive force causes and also in microwave emission (Tapping 1983; Grechnev et al. the nonlinear self-interaction of Alfvén waves in the Cohen- 2003). The aim of this article is to study compressible pertur- Kulsrud equation formalism. bations induced by long-wavelength weakly-nonlinear torsional We consider a weakly-nonlinear torsional wave and restrict waves, which are essentially non-plane. our attention to the linear terms of the compressible variables, described by the equations 2. Compressible flows induced by torsional waves 1 A ρ ∂V A B ∂2B p + B B − 0 0 − 0 z0 z π z0 z π ∂ π2 ∂ 2 We are interested in torsional waves of wavelength much longer 4 2 t 16 z 2 2 than the diameter of the wave-guiding magnetic flux tube with A0 J A0 ρ0 Ω = pext + − , (4) cylindrical coordinates (r,ϕ,z). In an untwisted and non-rotating T 8π2 2π tube with the equilibrium magnetic field Bz0, mass density ρ0, and circular cross-section area A , linear torsional waves are 0 ∂ ∂ ∂ twisting azimuthal motions vϕ accompanied by the appearance u p 1 2 J ρ0 + = − JR , (5) of the azimuthal component of the magnetic field Bϕ.Atthe ∂ t ∂ z 4π ∂z axis of the flux-tube, both quantities vanish, and hence can- where pext is the total pressure in the external medium, and not be described by the first order thin magnetic flux theory √ T = /π of Roberts & Webb (1978), while they appear in the second- R A0 is the flux-tube radius. The nonlinear terms associ- order thin flux-tube approximation of Zhugzhda (1996). The ated with the torsional wave are on the right hand side. In Eq. (5) approximation allows one to describe axisymmetric (m = 0, the term responsible for the ponderomotive force appears after where m is the azimuthal wave number) flows of plasma, in- accounting for the higher-order terms in the thin flux-tube ex- cluding sausage, longitudinal and torsional modes (Zhugzhda & pansion. Nakariakov 1999; Vasheghani Farahani et al. 2010) and their in- The compressible variables are expressed through the lon- teraction. Linear torsional waves are governed by the equations gitudinal component of the magnetic field perturbation Bz, with for the quantities Ω=vϕ/r and J = Bϕ/r, which in the thin the use of the linear expressions flux tube approximation correspond to the vorticity and electric ∂ρ ∂ u current density, respectively, + ρ + 2ρ V = 0, (6) ∂ t 0 ∂ z 0 ∂ Ω ∂ − Bz0 J = , ∂ πρ ∂ 0 (1) t 4 0 z ∂ B z + 2B V = 0, (7) ∂ t z0 ∂ J ∂ Ω − B = 0, (2) ∂ t z0 ∂ z ∂ p ∂ρ − C2 = 0, (8) where r is the radial coordinate. Also it is worth mentioning that ∂ t s ∂ t in the magnetic cylinder, a torsional wave could exist on any cylindrical shell (magnetic surface), having an arbitrary depen- where Cs is the sound speed.