Kink oscillations of cooling coronal loops with variable crosssection
^{1} Solar Physics and Space Plasma Research Centre (SP 2RC), University of Sheffield, Hicks Building, Hounsfield Road, Sheffield, S3 7RH, UK
email: m.s.ruderman@sheffield.ac.uk
^{2} Space Research Institute (IKI) Russian Academy of Sciences, 117342 Moscow, Russia
^{3} ITMO University, Kronverkskii ave 49, 197101 SaintPetersburg, Russia
^{4} Debrecen Heliophysical Observatory (DHO), Research Centre for Astronomy and Earth Sciences, Hungarian Academy of Sciences, 4010 Debrecen, PO Box 30, Hungary
^{5} Key Laboratory of Solar Activity, National Astronomical Observatories, Chinese Academy of Sciences, 100012 Beijing, PR China
Received: 29 November 2016
Accepted: 6 March 2017
We study kink waves and oscillations in a thin expanding magnetic tube in the presence of flow. The tube consists of a core region and a thin transitional region at the tube boundary. In this region the plasma density monotonically decreases from its value in the core region to the value outside the tube. Both the plasma density and velocity of background flow vary along the tube and in time. Using the multiscale expansions we derive the system of two equations describing the kink oscillations. When there is no transitional layer the oscillations are described by the first of these two equations. We use this equation to study the effect of plasma density variation with time on kink oscillations of an expanding tube with a sharp boundary. We assume that the characteristic time of the density variation is much greater than the characteristic time of kink oscillations. Then we use the WentzelKramerBrillouin (WKB) method to derive the expression for the adiabatic invariant, which is the quantity that is conserved when the plasma density varies. The general theoretical results are applied to the kink oscillations of coronal magnetic loops. We consider an expanding loop with the halfcircle shape and assume that the plasma temperature inside a loop decays exponentially with time. We numerically calculated the dependences of the fundamental mode frequency, the ratio of frequencies of the first overtone and fundamental mode, and the oscillation amplitude on time. We obtained that the oscillation frequency and amplitude increase and the frequency ratio decreases due to cooling. The amplitude increase is stronger for loops with a greater expansion factor. This effect is also more pronounced for higher loops. However, it is fairly moderate even for loops that are quite high.
Key words: magnetohydrodynamics (MHD) / plasmas / waves / Sun: corona / Sun: oscillations
© ESO, 2017
1. Introduction
After transverse oscillations of coronal magnetic loops were first observed by the Transition Region and Coronal Explorer (TRACE) spacecraft and reported by Aschwanden et al. (1999) and Nakariakov et al. (1999), they received ample attention from solar physicists. These oscillations have been interpreted as fast standing kink body waves in cylindrical magnetic flux tubes. In the earlier studies dealing with the interpretation of transverse oscillations, the simplest model of a coronal loop was used: a straight homogeneous magnetic flux tube (e.g. Ryutov & Ryutova 1976; Edwin & Roberts 1983). Later, more complex and realistic models of coronal loops were developed. For a review of theory of transverse coronal loop oscillations see e.g. Ruderman & Erdélyi (2009).
One particularly important property of coronal loops is the plasma density variation along the loop. Andries et al. (2005b) showed that this effect leads to the reduction of the ratio of the fundamental mode to the first overtone frequency below the canonical value of 2 of a homogeneous loop. On the basis of this result, Andries et al. (2005a) suggested the method of estimation of the atmospheric scale height in the solar corona, using the frequency ratio (see also McEwan et al. 2008, or for a review Andries et al. 2009). Dymova & Ruderman (2005) showed that kink oscillations of thin magnetic tubes with the density varying along the tube can be described by the wave equation with the phase speed equal to the local kink speed. Although typically the coronal loop expansion is relatively small, the ratio of the loop crosssection radii at the apex and at the footpoints can still reach a value of 1.5 (Watko & Klimchuk 2000; Klimchuk et al. 2000). It is surprising that in the chromosphere the expansion of vertical magnetic flux tubes can reach a value of up to a few hundred (e.g. Tsuneta et al. 2008). Ruderman et al. (2008) and Verth & Erdélyi (2008) generalised the derivation by Dymova & Ruderman (2005) for the case where the density and crosssection radius both vary along the loop. Using the wave equation that they derived, they showed that the two effects act in the opposite direction. While the density variation causes the reduction of the frequency ratio, the tube expansion causes its increase. Interestingly, these inhomogeneity effects show similar behaviour for sausage oscillations (LunaCardozo et al. 2012).
Observations by the Solar and Heliospheric Observatory (SoHO; Brekke et al. 1997; Winebarger et al. 2002), TRACE (Winebarger et al. 2001) and Hinode (Chae et al. 2008; Ofman & Wang 2008; Terradas et al. 2008) showed that flows are ubiquitously present in active region loops. Hence, it is important to study how flows in coronal loops modify the properties of kink oscillations. Gruszecki et al. (2008) have studied numerically the kink oscillations of coronal loops in the presence of fieldaligned flow using the slab model. Terradas et al. (2010) studied the effect of flow on the period and damping rate of propagating kink oscillations in a magnetic tube homogeneous in the longitudinal direction. Erdélyi et al. (2014) studied how a “jet”, i.e. a flow pulse, would affect the eigenoscillations as a jet propagates along the waveguide from one end to another. Ruderman (2010) derived the equation describing the kink oscillations of a magnetic tube with flow in the thin tube approximation. He assumed that the tube crosssection radius is constant; however, the density and flow velocity can vary both along the tube and in time. Using this equation he studied the effect of flow on the eigenfrequencies of kink oscillations.
It is very often observed that oscillating coronal loops are in a highly dynamic state. In particular, they can cool quickly with a characteristic cooling time of the order of a few periods of the kink oscillation (e.g. Aschwanden & Terradas 2008; Aschwanden & Schrijver 2011). This cooling can strongly affect the kink oscillations. The kink oscillations of cooling coronal loops have been studied by Morton & Erdélyi (2009, 2010). In particular, Morton & Erdélyi (2009) found that cooling causes the decrease in the period of kink oscillations. Ruderman (2011b) used the equation derived by Ruderman (2010) to study the kink oscillations of a coronal loop with the density varying along the loop and in time. In the approximation of slow temporal variation he derived the expression for an adiabatic invariant, which is a quantity conserved during the density evolution. Using this invariant he showed that cooling causes the amplification of kink oscillations. Similar results were obtained for longitudinal oscillations (see AlGhafri & Erdélyi 2013; AlGhafri et al. 2014).
Ruderman (2010) assumed that the density can vary along the loop, but it does not vary in the transverse direction. Ruderman (2011a) generalised this study to also include the density variation in the transverse direction from a larger value inside the loop to a smaller value in the surrounding plasma. He derived the system of equations describing the kink oscillations of a loop with such a density transition and studied the kink oscillations of cooling loops damped due to resonant absorption. Recently the kink oscillations of radiatively cooling coronal loops have been studied numerically by Magyar et al. (2015).
In this article, we generalise the studies by Dymova & Ruderman (2005), Ruderman et al. (2008), and Ruderman (2010, 2011a,b), and derive the system of equations describing the kink oscillations of an expanding loop with the density varying both along and across the loop, as well as in time, and in the presence of background flow. We then use this system to study the effect of cooling on kink oscillations of an expanding magnetic loop.
The paper is organised as follows. In the next section we describe the equilibrium state and present the governing equations. In Sect. 3 we transform the linearised MHD equations. In particular, we introduce the plasma displacement. In Sect. 4 we derive the governing equation for the kink oscillations in the thin tube approximation. In Sect. 5 we present the general analysis of the eigenvalue problem describing the kink oscillations in the presence of a stationary background flow. In Sect. 6 we derive the expression for the adiabatic invariant for the kink oscillations of an expanding magnetic tube with the slowly varying density. We then use this invariant to study the effect of cooling on the coronal loop kink oscillations. Section 7 contains the summary of the obtained results and our conclusions.
2. Equilibrium state and governing equations
We model a coronal loop as a straight, thin, and expanding magnetic tube with a circular crosssection. The tube can be either finite or semiinfinite, so that we can consider both standing and propagating waves. The tube consists of a core and a transition region where the density decreases from a higher value inside the tube to the lower value representing the surrounding plasma. In cylindrical coordinates r, φ, z with the zaxis coinciding with the tube axis, the plasma density is defined by (1)where l is a constant determining the thickness of a transitional layer, ρ(t,r,z) is a continuous function, and ρ_{t}(t,r,z) is a monotonically decreasing function of r. We explain below why we split the tube into the core and transition layer. The equilibrium magnetic field is B = (B_{r}(r,z),0,B_{z}(r,z)). The divergencefree condition for the magnetic field is written as (2)It follows from this equation that B can be expressed in terms of the magnetic flux function ψ as (3)We assume that the boundaries of the transitional layer are magnetic surfaces. This implies that the equations r = R(z)(1−l/ 2) and r = R(z)(1 + l/ 2) can be written as ψ = ψ_{i} = const. and ψ = ψ_{e} = const., respectively, where i and e refer to internal and external values.
In what follows we use the cold plasma approximation. Then, to satisfy the equilibrium condition, the magnetic field must be potential. This implies that B must satisfy the equation (4)which can be rewritten in terms of the magnetic flux function as (5)There is a timedependent background flow with the velocity U = (U_{r}(t,r,z),0,U_{z}(t,r,z)). The flow velocity is assumed to be parallel to the equilibrium magnetic field, U ∥ B. The plasma density and velocity are related by the mass conservation equation (6)The perturbations of the magnetic field and plasma velocity, b = (b_{r},b_{φ},b_{z}) and u = (u_{r},u_{φ},u_{z}), are described by the linearised MHD equations in the cold plasma approximation,
where μ_{0} is the magnetic permeability of free space.
3. Transformation of the governing equations
To derive the governing equations for kink oscillations we follow Ruderman et al. (2008). However now there are two features that make the derivation more complex. The first is the presence of the transitional layer between the tube core and the external plasma. The crucial step in the derivation of the governing equation for kink oscillations of a tube with a sharp boundary is to use two conditions at the tube boundary: the continuity of the plasma displacement in the direction normal to the boundary and the magnetic pressure perturbation. Here we introduce the jumps of these quantities across the transitional layer. After that the derivation of the governing equation is the same as in the case of a sharp boundary. The second feature is the presence of flow. In our derivation we need a closed system of equations for the plasma displacement and magnetic pressure perturbation. To obtain this system the perturbations of the plasma velocity and magnetic field need to be expressed in terms of the plasma displacement. When there is no background flow these expressions are trivial. It is still possible to obtain these expressions even when there is background flow. However, its presence makes them quite cumbersome. The aim of this section is to derive these expressions and then to obtain a closed system of equations for the plasma displacement and magnetic pressure perturbation. Hence, the main result of this section is the system of Eqs. (38)–(42).
Below we use the following identities (e.g. Korn & Korn 1961): (10)where F and G are arbitrary vectorfunctions.
We consider a Lagrangian particle with the initial position defined by vector a. A current position of this particle is x(t,a), so x(0,a) = a. The trajectory of this particle is given by x = x(t,a). In the unperturbed flow it is x = x_{0}(t,a). We introduce the plasma displacement ξ = x−x_{0}. In what follows, we consider ξ as a function of Eulerian coordinates, meaning that it is a function of t and x. The particle velocity is given by (11)where the subscripts a and x indicate that the derivative with respect to t is calculated at constant a and constant x, respectively. Now, ∂x_{0}/∂t = U(t,x_{0}). In the linear approximation we can substitute U(t,x_{0}) for in the last term on the righthand side of Eq. (11). Now, the Eulerian perturbation of the velocity is . Then, using the approximate relation U(t,x) ≈ U(t,x_{0}) + (ξ·∇)U we obtain from Eq. (11) in the linear approximation (12)where here and below the time derivative is calculated at constant x. Using the first identity in Eq. (10) we transform this expression to (13)Now we assume that the velocity U is potential, ∇ × U = 0. Then, again using the first identity in Eq. (10), we obtain (14)where u_{∥} = b_{0}·u and b_{0} = B/B is the unit vector in the magnetic field direction. Substituting this result in Eq. (7) and taking the scalar product of the obtained equation with b_{0} yields (15)Assuming that u_{∥} = 0 for t = 0 we obtain from this equation that u_{∥} = 0 for t> 0. Now, taking the scalar product of Eq. (13) with b_{0} and using the second identity in Eq. (10) yields (16)where ξ_{∥} = b_{0}·ξ and ξ_{1} = ξ−ξ_{∥}b_{0}.
With the aid of the third identity in Eq. (10) we obtain (17)where (18)Using Eq. (12) we transform Eq. (8) to (19)With the aid of Eq. (17) we reduce this equation to (20)Then, using the identity (21)we finally arrive at (22)It follows from this equation that b−f + B(∇·ξ) = 0 for t> 0 if it is zero at t = 0, which is assumed below. Hence, we obtain (23)Using the third identity in Eq. (10) and recalling that the magnetic field is divergentfree we can rewrite this expression as (24)We now introduce the magnetic pressure perturbation P = B·b/μ_{0}. Taking the scalar product of Eq. (24) with B, using the second identity in Eq. (10), and recalling that the equilibrium magnetic field is potential, we easily obtain (25)Substituting Eqs. (14) and (24) in Eq. (7), and taking into account that u_{∥} = 0 yields (26)Equations (13), (16), and (26) constitute the system of equations for ξ.
Now we introduce the components of the displacement and velocity that are perpendicular to the magnetic field lines and in the plains containing the zaxis, (27)Then, recalling that B is solenoidal and potential, we transform Eqs. (16) and (25) to where w = Bξ_{⊥}. To obtain the expression for u in terms of ξ we use Eq. (13). First, we take the scalar product of this equation with the vector (b_{0z},0,−b_{0r}) and use the third identity in Eq. (10) to obtain the expression for u_{⊥}. Then we take the φcomponent of Eq. (13) to obtain the expression for u_{φ}. As a result we arrive at
Now we take the scalar product of Eq. (26) with the vector (b_{0z},0,−b_{0r}). Using the third identity in Eq. (10) yields (32)Then we take the φcomponent of Eq. (26), and again use the third identity in Eq. (10) to obtain (33)Using Eq. (29) we reduce Eqs. (32) and (33) to (34)(35)We see that Eq. (28) is detached from the other equations and only serves to determine ξ_{∥}. Hence we do not use it below. Equations (29)–(31), (34), and (35) constitute the system of equations for ξ_{⊥}, ξ_{φ}, u_{⊥}, u_{φ}, and P.
Now, following Ruderman et al. (2008), we use the magnetic flux function ψ as an independent variable instead of r. With the aid of the chain rule we obtain the following relations: (36)where f is an arbitrary function, and the subscripts r and ψ indicate that the derivative is taken at constant r and constant ψ, respectively. In particular, taking f = r in these relations, we obtain (37)Using Eqs. (36) and (37) we transform Eqs. (29)–(31), (34), and (35) to
The system of Eqs. (38)–(42) is used in the next section to derive the governing equations for the kink oscillations of a thin and expanding magnetic flux tube. Obviously, u_{⊥} and u_{φ} can be eliminated using Eqs. (39) and (40). Hence, Eqs. (38)–(42) can be considered as a system of equations for w = Bξ_{⊥}, ξ_{φ}, and P.
4. Derivation of the governing equations
We now assume that the magnetic tube is thin, , where ϵ ≪ 1 and L is the characteristic scale of variation of all equilibrium quantities, like the tube radius and magnetic field magnitude, along the tube. In the case of standing waves we can take L equal to the tube length. We also assume that the characteristic scale of variation of perturbations along the tube is much larger than the tube radius. However, we do not make any assumption about the relation between L and the characteristic scale of perturbation variation. In particular, they can be of the same order. In accordance with the assumption that the tube is thin we introduce the stretching variable along the tube defined as Z = ϵz. The characteristic time of the problem is L/V_{∗}, where V_{∗} is the characteristic Alfvén speed. Its ratio to the Alfvénic time determined with respect to the tube radius, R_{∗}/V_{∗}, where R_{∗} is the characteristic value of the tube radius, is of the order of ϵ^{1}. Accordingly, we introduce the stretching time T = ϵt. The characteristic scale of variation of perturbations in the radial direction is R_{∗}. However, we assume that the characteristic scale of variation of ρ, B, and U outside of the transitional layer, i.e. in the core of the tube (r<R(z)(1−l/ 2)) and outside of the tube (r>R(z)(1 + l/ 2)), is L. While the characteristic scale of variation of B is L even in the transitional layer, ρ and U can vary in this layer on the scale lR_{∗}. The ratio of the radial and azimuthal component of the equilibrium magnetic field is of the order of ϵ, and the same estimate is valid for the radial and azimuthal component of the equilibrium flow. In accordance with these estimates, we introduce the scaled radial components of the magnetic field and bulk flow velocity, namely and .
From now on the derivation of governing equation is quite similar to that in Ruderman et al. (2008). It is based on using asymptotic expansions with respect to ϵ. The main result of this section is Eq. (75) with the righthand side given by Eq. (76).
Below we impose the restriction B_{z}> 0. Since the tube is thin, in the vicinity of the magnetic tube we can approximate ψ by the first term of its expansion in the Taylor series. Since the tube axis is a magnetic field line itself, ψ = const. at r = 0. Without loss of generality, we can take ψ = 0 at r = 0. Then it follows from Eq. (3) that the first term of the Taylor expansion of ψ with respect to r is proportional to r^{2} and we can use the approximate expression (43)Now it follows from Eq. (3) that the equilibrium magnetic field is given by (44)where the prime indicates the derivative. For the magnetic field magnitude we obtain (45)It follows from the magnetic flux conservation that (46)For the flow speed we obtain (47)It follows from the symmetry that ∂U/∂r = 0 at r = 0. Then, using Taylor’s theorem, we obtain in the core region (r<R(z)(1−l/ 2)) (48)where 0 <r_{∗}<R(z)(1−l/ 2). Outside of the tube in its vicinity (R(z)(1 + l/ 2) <r ≪ L) we have (49)where R(z)(1 + l/ 2) <r_{∗} ≪ L. Expressions similar to Eqs. (48) and (49) are also valid for ρ outside of the transitional layer and for B everywhere not far from the tube (r ≪ L). Hence, in the leading order approximation with respect to ϵ, we can consider U and ρ as quantities independent of r outside of the transitional layer, while B can be considered as a quantity independent of r everywhere in the region r ≪ L.
We now integrate Eq. (6) over the area of the tube core crosssection, that is over a circle of radius R(z)(1−l/ 2). Below we assume that l ≪ 1, meaning that one can neglect l in comparison with unity. Then we obtain (50)Since the velocity is parallel to the magnetic field, we have with the aid of Eqs. (44) and (46) (51)Substituting this result in Eq. (50) and recalling that B_{z} ≈ B, we eventually arrive at (52)This equation is valid in the core of the tube.
In the new variables we transform Eqs. (38)–(42) to
where . Below we consider ξ_{⊥} and ξ_{φ} as quantities of the order of unity. Then it follows from Eqs. (54) and (55) that and . With the aid of these estimates we obtain from Eqs. (56) and (57) that . In accordance with these estimates we put , , and P = ϵ^{2}B^{2}Q. Then, assuming that perturbations of all quantities are proportional to e^{iφ}, we reduce the system of Eqs. (53)–(57) to
Eliminating , , and ξ_{φ} and only keeping the leading terms with respect to ϵ we obtain the system of equations for w and Q in the leading order approximation: (63)(64)where W = rw = rBξ_{⊥}. Using Eq. (43) and recalling that, outside of the transitional layer, we can consider ρ, B, and U as quantities independent of r and, consequently, on ψ, we reduce Eqs. (63) and (64) to
We differentiate Eq. (66) with respect to ψ and then eliminate W from the obtained equation and Eq. (65). This yields (67)The solution to this equation must be regular at ψ = 0 and decay as ψ → ∞. Hence, we obtain (68)where Q_{i}(T,Z) and Q_{e}(T,Z) are arbitrary functions. Recall that ψ = ψ_{i} and ψ = ψ_{e}, where ψ_{i} and ψ_{e} are constants, are the equations of the boundaries of the transitional layer. It follows from Eq. (43) that (69)Equations (66) and (68) imply that inside the tube core region. Then, recalling the definition of W and using Eq. (43) we conclude that ξ_{⊥} is independent of r inside the tube core region, which implies that, in the thin tube approximation, the tube oscillates as a solid.
Substituting Eq. (68) in Eq. (66) we arrive at 10pt
where W_{i} and W_{e} are calculated at ψ = ψ_{i} and ψ = ψ_{e}, respectively. We introduce the new variable (72)We also introduce the jumps of this new variable and the magnetic pressure perturbation across the transitional layer, (73)We have the estimates , . Then, with the aid of Eq. (68), we obtain (74)Now we multiply Eq. (70) by ρ_{i}, Eq. (71) by ρ_{e}, and add the results; use Eqs. (46), (73), and (74); and return to the original nonscaled independent variables. As a result, we arrive at We note that we only kept terms of the order of l when deriving Eq. (76), while we neglected terms of higher orders with respect to l. When there is no transitional layer (l = 0) we have ℒ = 0, and Eq. (75) governs the kink oscillations of a thin expanding magnetic tube with a background fieldaligned flow. However, in the presence of the transitional layer the system of Eqs. (75) and (76) is not closed. To close it one needs to express δη and δP in terms of η. We do not do this here. These expressions are obtained when we use Eqs. (75) and (76) to study the damping of kink oscillations of expanding magnetic flux tubes due to resonant absorption.
When there is no flow (U = 0) and no transitional region (l = 0), Eq. (75) reduces to the equation describing the kink oscillations of a thin expanding magnetic tube derived by Ruderman et al. (2008). On the other hand, when there is no expansion (R = const.), Eqs. (75) and (76) reduce to Eqs. (23) and (24) in Ruderman (2011a) describing resonantly damped kink oscillations of a thin magnetic tube with nonstationary density and flow. We note that the expression for ℒ given by Eq. (76) with R = const. is slightly different from that given by Equations (24) in Ruderman (2011a). However, it is straightforward to show using Eq. (75) that the difference between Rℒ with ℒ given by Eq. (76) and ℒ given by Eqs. (24) in Ruderman (2011a) is of the order of l^{2}.
It is expedient to compare the results of this section with that obtained by Andries & Cally (2011). They considered waves in static magnetic tubes with sharp boundaries. They assumed that the characteristic scale of the tube crosssection variation along the tube is much larger than the characteristic transverse scale of the tube; however, they did not assume that the perturbation wavelength is much larger than the characteristic transverse scale. In particular, they managed to obtain a single equation describing kink oscillations. However, this equation is quite complicated because it contains functions of differential operators. It is quite difficult to work with them. We see that the thin tube approximation enormously simplifies the study of kink oscillations. Even when the tube is nonstationary and there is background flow, kink oscillations are described by a relatively simple differential Eq. (75).
5. Eigenvalue problem in the presence of stationary flow
In what follows we consider kink oscillations of magnetic tubes with the sharp boundary and take l = 0. We assume that there is no flow outside the tube, U_{e} = 0. We also assume that both the density and background flow velocity are independent of time. Then Eq. (52) reduces to (77)where we have dropped the subscript “i” at the internal flow speed U. Using Eq. (45) we rewrite Eq. (46) in the approximate form (78)Since we have taken l = 0 it follows that ℒ = 0. Below we study standing waves and impose the boundary conditions (79)where L is the tube length. We look for eigenfrequencies of kink oscillations, so we take η ∝ e^{− iωt}. Then Eq. (75) with ℒ = 0 reduces to (80)This equation, together with the boundary conditions Eq. (79), constitutes an eigenvalue problem. The aim of this section is to study general properties of this problem. Below we assume that μ_{0}ρ_{i}U^{2}< 2B^{2}, that is , and this condition is satisfied for all z ∈ [−L/ 2,L/ 2]. Then we make the variable substitution (81)Now, multiplying the obtained equation by R^{4}, we reduce Eq. (80) to (82)where (83)Equation (82) with the boundary conditions Eq. (79) constitutes the classical SturmLiouville problem with the eigenvalue ω^{2}. It follows from the general theory of the SturmLiouville problem (e.g. Coddington & Levinson 1955) that the eigenvalues are real and constitute a monotonically increasing unbounded sequence. The first eigenvalue is the square of the fundamental frequency, and the corresponding eigenfunction has no nodes in (− L/ 2,L/ 2). All other eigenvalues are the squares of frequencies of overtones. The eigenfunction corresponding to the nth overtone has n−1 nodes in (− L/ 2,L/ 2).
It is obvious that we can always take q to be real. Multiplying Eq. (82) by q, integrating the obtained equation, and using Eq. (79) we obtain (84)Since W(z) > 0, it follows from this equation that ω^{2}> 0, i.e. that all the eigenfrequencies are real. This implies that the inequality is a sufficient condition for the stability of the tube with respect to long kink perturbations. We emphasise that this condition is sufficient but not necessary. On the other hand, this analysis does not prove that the tube is stable if the condition is satisfied. While it guaranties that the tube is stable with respect to long kink perturbations, it can be unstable with respect to other types of perturbations.
6. Kink oscillations of coronal loops with slowly varying density
The analysis in this section is motivated by the observations of kink oscillations of cooling coronal loops (e.g. Aschwanden & Terradas 2008; Aschwanden & Schrijver 2011). Below we consider kink oscillations of a magnetic tube with slowly varying density. We assume that the boundary condition Eq. (79) is satisfied. We also assume that the characteristic time of the density variation, t_{ch}, is much greater than the characteristic time of oscillations. We note that this is not a very restrictive assumption. It is usually believed that the characteristic time of oscillations is the oscillation period. However, this is not true. The characteristic time of oscillations is 1 /ω, where ω is the oscillation frequency. This time is about 6 times smaller than the oscillation period. Hence, the condition that t_{ch} should be much larger than the characteristic time of oscillations is satisfied even when t_{ch} is of the order of the oscillation period. We also assume that the loop has a sharp boundary, which is l = 0 and, consequently, ℒ = 0.
6.1. Derivation of adiabatic invariant
We denote the ratio of t_{ch} to the characteristic time of oscillations as ν^{1}, where ν ≪ 1. Hence, the characteristic time of oscillations is νt_{ch}. On the other hand, it is also of the order of the loop length divided by the characteristic kink speed, which is equal to , where ρ_{ch} is the characteristic density. As a result, we have
This estimate inspired us to introduce the scaled magnetic field . Then we rewrite Eq. (75) with ℒ = 0 as (85)Now we use the WentzelKramerBrillouin (WKB) method (e.g. Bender & Orszag 1978) and look for the solution to this equation in the form (86)Then we expand S in the series (87)Substituting Eq. (86) in Eq. (85), using Eq. (87), collecting terms of the order of ν^{2}, and taking into account the boundary condition Eq. (79) we obtain the SturmLiouville problem (88)where (89)This order approximation is called the approximation of geometrical optics (e.g. Bender & Orszag 1978). The SturmLiouville problem Eq. (88) coincides with a problem obtained by Dymova & Ruderman (2005) to describe the kink oscillations of magnetic tubes whose density varies along the tube, and by Ruderman et al. (2008) to describe the kink oscillations of magnetic tubes whose density and crosssection radius vary along the tube. Below we assume that Ω^{2} is an eigenvalue and S_{0} the corresponding eigenfunction. In accordance with the general theory (e.g. Coddington & Levinson 1955), Ω^{2} is real. It is easy to prove that Ω^{2}> 0. Obviously we can always assume that S_{0} is a real function.
Now we collect the terms of the order of ν^{1} to obtain (90)This is the approximation of the physical optics. S_{1} must satisfy the boundary conditions (91)The homogeneous counterpart of the boundary value problem constituted by Eq. (90) and the boundary conditions Eq. (91) has a nontrivial solution S_{1} = S_{0}. This implies that the boundary value problem determining S_{1} has a solution only if the righthand side of Eq. (90) satisfies the compatibility condition, which is the condition that it should be orthogonal to S_{0}. To obtain this condition we multiply Eq. (90) by S_{0} and then integrate with respect to z from −L/ 2 to L/ 2. After some algebra we obtain (92)Using the integration by parts, Eq. (52), and the relation yields (93)Substituting Eq. (93) in Eq. (92) and using Eq. (89) we eventually obtain (94)where (95)The lefthand side of Eq. (94) is the adiabatic invariant. Equation (94) states that this invariant is conserved. It is worth noting that the adiabatic invariant is the same as that obtained by Ruderman (2011b) for oscillations of magnetic tubes with the constant crosssection radius. Hence, the tube expansion only affects the time evolution of kink oscillations of magnetic tubes with the varying density through C_{k} and S_{0} that both depend on R(z).
6.2. Effect of cooling on the kink oscillations of coronal magnetic loops
In this subsection we consider the kink oscillations of cooling coronal loops. The main cause of the loop cooling is radiation. Its intensity is proportional to the plasma density squared. This implies that the energy deposition that is sufficient to cover the energy losses in the rarefied plasma outside the loop may be too small to cover the energy losses inside the loop. Hence, it seems to be a viable assumption that while the plasma inside the loop is cooling, the temperature of the plasma outside the loop remains constant. Following Aschwanden & Terradas (2008) and Morton & Erdélyi (2010), we approximate the temperature evolution inside the loop by an exponentially decaying function, (96)where we have assumed that the cooling starts at the initial time t = 0. We also assume that the plasma temperature inside and outside the loop was the same at t = 0, implying that the temperature of the external plasma remains equal to T_{0}.
Below we consider loops with a halfcircle shape. We neglect the effect of the loop curvature on kink oscillations. Hence, the loop shape only determines the density variation along the loop. In the external plasma it is given by the barometric formula (97)where ρ_{f} is the plasma density at the footpoints inside the loop (assumed to be constant), and ζ the ratio of the plasma densities inside and outside the loop at t = 0, (98)k_{B} is the Boltzmann constant, g the gravity acceleration, and m the mean mass per particle equal to one half of the proton mass for a protonelectron plasma. The cooling causes the plasma flow inside the loop. This flow affects the plasma density meaning that it is not described by the barometric formula. However, Ruderman (2011b) showed that for typical coronal conditions and observed cooling times the effect of this flow on the plasma density is fairly week. Hence, the barometric formula provides a sufficiently good approximation for the plasma density inside the loop. It reads (99)where H(t) is defined by Eq. (98) with T(t) substituted for T_{0}.
To describe the crosssection radius variation along the loop we use the same expression as in Ruderman et al. (2008), (100)where R_{f} is the crosssection radius at the footpoints, l_{c} an arbitrary positive constant with the dimension of length, and λ = R(0) /R_{f} the expansion factor. It is shown in Ruderman et al. (2008) that in order to have the zcomponent of the magnetic field positive everywhere in the region  z  ≤ L/ 2, the expansion factor must satisfy the inequality λ<λ_{m}, where (101)The parameter λ_{m} is a monotonically increasing function of L/l_{c}, λ_{m} → 1 as L/l_{c} → 0, and λ_{m} → 1.87 as L/l_{c} → ∞. As we have already pointed out, the coronal loop expansion factor typically does not exceed 1.5. Hence, by properly choosing the ratio L/l_{p} we can cover the whole range of the expansion factor variation.
We introduce the dimensionless variables and parameters (102)where the kink speed at the footpoints is defined by (103)Then Eq. (88) reduces to (104)When deriving this equation we dropped the prime at S_{0} and used the relation BΛ^{2} = B_{f}. The function Λ(Z) is defined by Eq. (100). The eigenfunction corresponding to the fundamental mode is even, while the eigenfunction corresponding to the first overtone is odd. This enables us to solve Eq. (104) not on the interval [− 1,1], but on the interval [0,1] The solution corresponding to the fundamental mode must satisfy the boundary conditions (105)and the solution corresponding to the first overtone must satisfy the boundary conditions (106)Since the eigenfunction S_{0} is accurately determined up to the multiplication by an arbitrary function of τ, we can always fix its value at one particular point. Let X(τ,Z) be an eigenfunction corresponding to the fundamental mode that satisfies the condition X(τ,0) = 1. Then the general solution to the eigenvalue problem corresponding to the fundamental mode is S_{0}(τ,Z) = A(τ)X(τ,Z), where A(τ) is the oscillation amplitude at Z = 0. Then Eq. (94) reduces to (107)Equations (104) and (107) have been solved numerically for ζ = 3 and l_{c} = L/ 6. In this case λ_{m} ≈ 1.67. In Figs. 1–3 the dependences of the dimensionless fundamental mode frequency, the ratio of frequencies of the first overtone and fundamental mode, and the oscillation amplitude at Z = 0 on τ are shown for various values of κ and λ.
Fig. 1 Dependence of the dimensionless fundamental mode frequency on τ. The upper, middle, and lower panels correspond to κ = 0.5, κ = 1, and κ = 2. The dotted, solid, dashed, and dashdotted lines correspond to λ = 1,1.15,1.3, and 1.5, respectively. 

Open with DEXTER 
Fig. 2 Dependence of the ratio of frequencies of the first overtone and fundamental mode on τ. The upper, middle, and lower panels correspond to κ = 0.5, κ = 1, and κ = 2. The dotted, solid, dashed, and dashdotted lines correspond to λ = 1,1.15,1.3, and 1.5, respectively. 

Open with DEXTER 
Fig. 3 Dependence of the dimensionless amplitude of the fundamental mode at Z = 0 on τ. The upper, middle, and lower panels correspond to κ = 0.5, κ = 1, and κ = 2. The dotted, solid, dashed, and dashdotted lines correspond to λ = 1,1.15,1.3, and 1.5, respectively. 

Open with DEXTER 
We see in Fig. 1 that the oscillation frequency grows for all values of the expansion factor and for all values of the loop height. The same result was obtained by Morton & Erdélyi (2009, 2010) and Ruderman (2011b) for nonexpanding magnetic loops. This result has a simple physical explanation. When the loop is cooling the plasma flow with the velocity directed to the loop footpoints appears. The plasma is evacuated from the loop meaning that the loop mass decreases. At the same time the restoring force created by the magnetic tension remains the same. As a result the oscillation frequency increases.
Figure 2 shows that the ratio of frequencies of the first overtone and fundamental harmonic decreases with time. Again the same result was obtained by Morton & Erdélyi (2009, 2010) and Ruderman (2011b) for nonexpanding magnetic loops. And again there is a simple physical explanation for this result. Owing to cooling, the plasma density at the loop decreases everywhere except for the loop footpoints. In particular, it decreases at the loop apex. As a result the kink speed at the loop apex increases, while it remains the same at the loop footpoints. Ruderman et al. (2016) showed that the frequency ratio is a monotonically decreasing function of the ratio of the kink speed at the loop apex and footpoints.
The most interesting result of this section is shown in Fig. 3. We see that cooling enhances the oscillation amplitude. Previously the same result was obtained by Ruderman (2011b) for nonexpanding loops. We also see that the loop expansion makes the oscillation amplification stronger. This effect becomes more pronounced when the loop height increases. However, even for loops whose height is equal to two times the initial atmospheric scaleheight, the effect is fairly weak.
7. Summary and conclusions
In this article, we studied kink waves and oscillations of thin expanding magnetic tubes in the presence of flow. Both the plasma density and flow velocity were allowed to vary along the tube and depend on time. These physical quantities can also vary in the radial direction. However, it is assumed that the characteristic scale of their variation in the radial direction in the core region inside the tube and in the region outside the tube is the same as that along the loop. They can only vary on short spatial scales inside a thin transitional layer at the tube boundary.
Using the multiscale expansion we derived the system of two equations (see Eqs. (75) and (76)) governing the kink waves and oscillations. This is the main result of the article. We emphasise that this system describes both standing and propagating waves. When there is no transitional layer the oscillations are only governed by the first equation. On the other hand, in the presence of the transitional layer the system of two equations is not closed because it contains the jumps of the plasma displacement and magnetic pressure across the transitional layer. To close this system these jumps must be expressed in terms of the tube displacement.
We used the derived equations to study the effect of the plasma density variation with time. We assumed that the characteristic time of the density variation is much greater than the characteristic time of oscillations. Then we used the WKB method to derive the expression for the adiabatic invariant, which is the quantity that is conserved while the plasma density changes. The general theory was applied to kink oscillations of cooling coronal loops. We considered expanding loops with a halfcircle shape and assumed that the plasma temperature inside the loop decays exponentially. Then we numerically calculated the time dependence of the fundamental mode frequency, the ratio of frequencies of the first overtone and fundamental mode, and the oscillation amplitude.
The numerical results can be summarised as follows. Cooling causes an increase in the oscillation frequency and a decrease in the ratio of frequencies of the first overtone and fundamental mode. This is the generalisation of the results previously obtained by Morton & Erdélyi (2009, 2010) and Ruderman (2011b) for nonexpanding loops. Cooling also causes the enhancement of the oscillation amplitude. Previously, this result was obtained by Ruderman (2011b) for nonexpanding loops. The amplitude enhancement due to cooling becomes stronger when the loop expansion increases. This effect is more pronounced for higher loops. However, it is very moderate even for loops where the height is equal to the two times the initial atmospheric scale height.
Acknowledgments
M.S.R. and R.E. acknowledge the support by the STFC grant. R.E. acknowledges the support received from the Chinese Academy of Sciences President’s International Fellowship Initiative, Grant No. 2016VMA045. R.E. also acknowledges the excellent hospitality of the National Astronomical Observatory of China (NAOC), Chinese Academy of Sciences (CAS), where a part of this work was carried out.
References
 AlGhafri, K. S., & Erdélyi, R. 2013, Sol. Phys., 283, 413 [NASA ADS] [CrossRef] [Google Scholar]
 AlGhafri, K. S., Ruderman, M. S., Williamson, A., & Erdélyi, R. 2014, ApJ, 786, A36 [NASA ADS] [CrossRef] [Google Scholar]
 Andries, J., & Cally, P. S. 2011, ApJ, 743, 164 [NASA ADS] [CrossRef] [Google Scholar]
 Andries, J., Arregui, I., & Goossens, M. 2005a, ApJ, 624, L57 [NASA ADS] [CrossRef] [Google Scholar]
 Andries, J., Goossens, M., Hollweg, J. V., Arregui, I., & Van Doorsselaere, T. 2005b, A&A, 430, 1109 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Andries, J., Van Doorsselaere, T., Roberts, B., et al. 2009, Space Sci. Rev., 149, 3 [NASA ADS] [CrossRef] [Google Scholar]
 Aschwanden, M. J., & Schrijver, C. J. 2011, ApJ, 736, A102 [NASA ADS] [CrossRef] [Google Scholar]
 Aschwanden, M. J., & Terradas, J. 2008, ApJ, 686, L127 [NASA ADS] [CrossRef] [Google Scholar]
 Aschwanden, M. J., Fletcher, L., Schrijver, C. J., & Alexander, D. 1999, ApJ, 520, 880 [NASA ADS] [CrossRef] [Google Scholar]
 Bender, C. M., & Orszag, S. A. 1978, Advanced Mathematical Methods for Scientists and Engineers (New York: McGrawHill) [Google Scholar]
 Brekke, P., KjeldsethMoe, O., & Harrison, R. A. 1997, Sol. Phys., 175, 511 [NASA ADS] [CrossRef] [Google Scholar]
 Chae, J., Ahn, K., Lim, E. K., Choe, G. S., & Sakurai, T. 2008, ApJ, 689, L73 [NASA ADS] [CrossRef] [Google Scholar]
 Coddington, E. A., & Levinson, N. 1955, Theory of Ordinary Differential Equations (New York: McGrawHill) [Google Scholar]
 Dymova, M. V., & Ruderman, M. S. 2005, Sol. Phys., 229, 79 [NASA ADS] [CrossRef] [Google Scholar]
 Edwin, P. M., & Roberts, B. 1983, Sol. Phys., 88, 179 [NASA ADS] [CrossRef] [Google Scholar]
 Erdélyi, R., Hague, A., & Nelson, C. J. 2014, Sol. Phys., 289, 167 [NASA ADS] [CrossRef] [Google Scholar]
 Gruszecki, M., Murawski, K., & Ofman, L. 2008, A&A, 488, 757 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Klimchuk, J. A., Antiochos, S. K., Norton, D., & Watko, J. A. 2000, AAS, SPD meeting, 31 [Google Scholar]
 Korn, G., & Korn, T. 1961, Mathematical Handbook for Scientists and Engineers (New York: McGrawHill) [Google Scholar]
 LunaCardozo, M., Verth, G., & Erdélyi, R. 2012, ApJ, 748, 110 [NASA ADS] [CrossRef] [Google Scholar]
 Magyar, N., Van Doorsselaere, T., & Marcu, A. 2015, A&A, 582, A117 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 McEwan, M. P., Díaz, A. J., & Roberts, B. 2008, A&A, 481, 819 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Morton, R. J., & Erdélyi, R. 2009, ApJ, 707, 750 [NASA ADS] [CrossRef] [Google Scholar]
 Morton, R. J., & Erdélyi, R. 2010, A&A, 519, A43 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Nakariakov, V. M., Ofman, L., Deluca, E. E., Roberts, B., & Davila, J. M. 1999, Science, 285, 862 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Ofman, L., & Wang, T. J. 2008, A&A, 482, L9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Ruderman, M. S. 2010, Sol. Phys., 267, 377 [NASA ADS] [CrossRef] [Google Scholar]
 Ruderman, M. S. 2011a, A&A, 534, A78 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Ruderman, M. S. 2011b, Sol. Phys., 271, 41 [NASA ADS] [CrossRef] [Google Scholar]
 Ruderman, M. S., & Erdélyi, R. 2009, Space Sci. Rev., 149, 199 [NASA ADS] [CrossRef] [Google Scholar]
 Ruderman, M. S., Verth, G., & Erdélyi, R. 2008, ApJ, 686, 694 [NASA ADS] [CrossRef] [Google Scholar]
 Ruderman, M. S., Petrukhin, N. S., & Pelinovsky, E. 2016, Sol. Phys., 291, 1143 [NASA ADS] [CrossRef] [Google Scholar]
 Ryutov, D. D., & Ryutova, M. P. 1976, Sov. Phys. – JETP, 43, 491 [NASA ADS] [Google Scholar]
 Terradas, J., Arregui, I., Oliver, R., & Ballester, J. L. 2008, ApJ, 678, L153 [NASA ADS] [CrossRef] [Google Scholar]
 Terradas, J., Goossens, M., & Ballai, I. 2010, A&A, 515, A46 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Tsuneta, S., Ichimoto, K., Katsukawa, Y., et al. 2008, ApJ, 688, 1374 [NASA ADS] [CrossRef] [Google Scholar]
 Verth, G., & Erdélyi, R. 2008, A&A, 486, 1015 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Watko, J. A., & Klimchuk, J. A. 2000, Sol. Phys., 193, 77 [NASA ADS] [CrossRef] [Google Scholar]
 Winebarger, A. R., DeLuca, E. E., & Golub, L. 2001, ApJ, 553, L81 [NASA ADS] [CrossRef] [Google Scholar]
 Winebarger, A. R., Warren, H., van Ballegooijen, A., DeLuca, E. E., & Golub, L. 2002, ApJ, 567, L89 [NASA ADS] [CrossRef] [Google Scholar]
All Figures
Fig. 1 Dependence of the dimensionless fundamental mode frequency on τ. The upper, middle, and lower panels correspond to κ = 0.5, κ = 1, and κ = 2. The dotted, solid, dashed, and dashdotted lines correspond to λ = 1,1.15,1.3, and 1.5, respectively. 

Open with DEXTER  
In the text 
Fig. 2 Dependence of the ratio of frequencies of the first overtone and fundamental mode on τ. The upper, middle, and lower panels correspond to κ = 0.5, κ = 1, and κ = 2. The dotted, solid, dashed, and dashdotted lines correspond to λ = 1,1.15,1.3, and 1.5, respectively. 

Open with DEXTER  
In the text 
Fig. 3 Dependence of the dimensionless amplitude of the fundamental mode at Z = 0 on τ. The upper, middle, and lower panels correspond to κ = 0.5, κ = 1, and κ = 2. The dotted, solid, dashed, and dashdotted lines correspond to λ = 1,1.15,1.3, and 1.5, respectively. 

Open with DEXTER  
In the text 