Solar Phys (2018) 293:86 https://doi.org/10.1007/s11207-018-1305-6 Magnetoacoustic Waves and the Kelvin–Helmholtz Instability in a Steady Asymmetric Slab I: The Effects of Varying Density Ratios 1 1,2 M. Barbulescu · R. Erdélyi Received: 10 July 2017 / Accepted: 10 May 2018 © The Author(s) 2018 Abstract Recent observations have shown that bulk ﬂow motions in structured solar plas- mas, most evidently in coronal mass ejections (CMEs), may lead to the formation of Kelvin– Helmholtz instabilities (KHIs). Analytical models are thus essential in understanding both how the ﬂows affect the propagation of magnetohydrodynamic (MHD) waves, and what the critical ﬂow speed is for the formation of the KHI. We investigate both these aspects in a novel way: in a steady magnetic slab embedded in an asymmetric environment. The exterior of the slab is deﬁned as having different equilibrium values of the background density, pres- sure, and temperature on either side. A steady ﬂow and constant magnetic ﬁeld are present in the slab interior. Approximate solutions to the dispersion relation are obtained analytically and classiﬁed with respect to mode and speed. General solutions and the KHI thresholds are obtained numerically. It is shown that, generally, both the KHI critical value and the cut-off speeds for magnetoacoustic waves are lowered by the external asymmetry. Keywords MHD waves · Kelvin–Helmholtz instability 1. Introduction The propagation of linear magnetohydrodynamic (MHD) waves along magnetic slabs has long been a topic of study in the context of solar physics (see, e.g. Roberts, 1981). The presence of a steady ﬂow in the equilibrium state of the system affects the propagation in at least two important ways. First, perturbations may cause shearing motions in the ﬂow, which then could lead to the Kelvin–Helmholtz instability (KHI) (see Figure 1). Second, the M. Barbulescu mbarbulescu1@shefﬁeld.ac.uk R. Erdélyi robertus@shefﬁeld.ac.uk Solar Physics and Space Plasma Research Centre (SP2RC), School of Mathematics and Statistics, University of Shefﬁeld, Hounsﬁeld Road, Hicks Building, Shefﬁeld, S3 7RH, UK Department of Astronomy, Eötvös Loránd University, Pázmány P. sétány 1/A, Budapest, 1117, Hungary 86 Page 2 of 18 M. Barbulescu, R. Erdélyi Figure 1 Stages of a Kelvin–Helmholtz instability (KHI). We assume that a magnetic interface (a) separating two regions with background ﬂows in opposite directions is subject to a perturbation (b). As the system evolves in time, sufﬁciently strong ﬂows will amplify the perturbation, causing nonlinear wave steepening (c), until vortex formation occurs (d). Further evolution typically renders the system turbulent. phase speeds and the cut-off speeds of each mode of propagation are shifted proportional to the speed of the ﬂow (see, e.g. Nakariakov and Roberts, 1995). Interactions between propagating waves and ﬂows are not limited to these two instances, however. Other areas of study include negative-energy wave instabilities, if dissipative effects are taken into account (Cairns, 1979; Joarder, Nakariakov, and Roberts, 1997), or resonant ﬂow instabilities, if resonant wave excitation is considered (see Tirry et al., 1998; Taroyan and Erdélyi, 2002). More information on the above topics may be found in Taroyan and Ruderman (2011)and Ryutova (2015). The effects of steady ﬂows have been investigated in a number of different waveguide geometries and magnetic topologies. Nakariakov and Roberts (1995) studied the effect of a steady ﬂow in an inﬁnite slab of magnetised plasma in a magnetic environment. Terra- Homem, Erdélyi, and Ballai (2003) then explored the effects that a steady ﬂow has on the propagation of both linear and nonlinear waves in a straight inﬁnite cylindrical ﬂux tube. This latter work expanded on the analysis of Somasundaram, Venkatraman, and Sengottuvel (1999). For a more general approach to analysing the stability of steady MHD ﬂows, see, for example, Goedbloed (2009a,b). More recently, Soler et al. (2010) described the effects of an azimuthally dependent ﬂow on the stability of a straight ﬂux tube, while Zaqarashvili, Vörös, and Zhelyazkov (2014) investigated the stability of an incompressible, twisted cylindrical ﬂux tube, subject to a straight ﬂow, in a magnetic environment. Finally, Zaqarashvili, Zhelyazkov, and Ofman (2015) studied the stability of an incompressible, rotating, and twisted cylinder. The theo- retical results of the latter two works were applied in Kuridze et al. (2016) to determine the stability of chromospheric jets, and to estimate the growth time of the KHI. Recent observational results have reinforced the idea that plasma ﬂows are present throughout the solar atmosphere. Berger et al. (2010) and Ryutova et al. (2010) uncovered details about mass ﬂows and the formation of the KHI in solar prominences. KHI formation in the corona has also received considerable attention (see Foullon et al., 2011, 2013;Ofman and Thompson, 2011). For a recent review, see Zhelyazkov (2015). Of signiﬁcant interest are the observations by Foullon et al. (2011)ofaKHI on theﬂank of a CME. The authors interpret the system conﬁguration as consisting of three regions: the dense solar ejecta, the CME sheath, and the low-density corona, with the KHI occur- ring in the region between the ejecta sheath and the corona. A similar three-layer system is described by Möstl, Temmer, and Veronig (2013). By interpreting the CME boundary as a steady magnetic slab embedded in an asymmetric magnetic environment, the authors demonstrated that through increasing the magnetic ﬁeld strength on only one side of the slab, the ﬁeld provided a stabilising effect on that side only. This numerical study shows that Magnetoacoustic Waves and the Kelvin–Helmholtz Instability. . . Page 3 of 18 86 Figure 2 Steady magnetic slab embedded in a static asymmetric unmagnetised environment. exterior asymmetry may be an important factor when considering the physics of magnetic slabs. The three-layer system, envisioned as a slab in an asymmetric environment, has recently been studied by Allcock and Erdélyi (2017) and Zsámberger, Allcock, and Erdélyi (2018) in the context of linear wave propagation. Here, we focus on the effects that a steady ﬂow within the slab has on the propagation of magnetoacoustic waves, and on how the asymmetry affects the KHI threshold values. In Section 2, we assume that our system is governed by the ideal MHD equations, and we derive the dispersion relation for waves propagating along the slab. In Section 3, we obtain approximate solutions to the dispersion relation in the thin slab limit, and classify the modes in terms of the characteristic speeds of the system. In Section 4, we obtain general solutions to the dispersion relation and also the KHI thresholds. Finally, Section 5 summarises the results and provides context for their implications. 2. The Dispersion Relation We introduce a slab of plasma bounded by two interfaces at ±x , of density, pressure, and temperature ρ , p ,and T , respectively, and magnetic ﬁeld B = (0, 0,B ), which is subject 0 0 0 0 0 to a steady ﬂow U = (0, 0,U ). The slab is embedded in an asymmetric environment, de- 0 0 ﬁned as having density, pressure, and temperature ρ , p ,and T on the left side, and ρ , p , 1 1 1 2 2 and T on the right side, as illustrated in Figure 2. The exterior is neither subject to magnetic ﬁelds nor to ﬂows. It follows that the ﬂuid in the interior region of the slab is governed by the ideal MHD equations, while the exterior regions are described using the gas equations. We wish to obtain a governing equation describing the propagation of linear magnetoa- coustic waves along the parallel interfaces. Linearising the ideal MHD equations, subject to the previously deﬁned background conditions, allows us to write them in the form Dρ + ρ ∇· v = 0, Dt Dv 1 B ∂ b ρ =−∇ p + b B + , 0 z 0 Dt μ μ ∂z (1) Dp Dρ = c , Dt Dt Db ∂ v =−B (∇· v) + B . 0 0 Dt ∂z 86 Page 4 of 18 M. Barbulescu, R. Erdélyi Here ρ, p, b = (b ,b ,b ),and v = (v ,v ,v ) are small perturbations from the equi- x y z x y z D ∂ ∂ librium, and = + U is the material derivative. The sound speed is deﬁned as Dt ∂t ∂z c = γp /ρ . 0 0 Since we are only concerned with magnetoacoustic waves, we may disregard all depen- dence on the y -component without loss of generality. Equations 1 may, thus, be written in component form as Dv ∂ B B ∂b x 0 0 x ρ =− p + b + , 0 z Dt ∂x μ μ ∂z 0 0 Dv ∂p ρ =− , Dt ∂z Dp (2) =−c ρ ∇· v, Dt Db ∂v x x = B , Dt ∂z Db ∂v z x =−B . Dt ∂x We Fourier-decompose Equations 2 for waves propagating along the slab by assuming −i(ωt −kz) ˆ ˆ that f(r,t) = f(x)e ,where f stands for any of the small perturbations, and f is the amplitude of each perturbation. Here, ω is the angular frequency, and k is the wavenumber in the z-direction. This procedure allows us to remove all differential terms in the linearised MHD equation, except for derivatives with respect to x . Equations 2 become d B B 0 0 ˆ ˆ iρ v ˆ = p ˆ + b + ik b , 0 x z x dx μ μ 0 0 ρ v ˆ = kp, ˆ 0 z dv ˆ (3) p ˆ = c ρ −i + kv ˆ , 0 z dx b =−B kv ˆ , x 0 x dv ˆ ib = B , z 0 dx where = ω − kU is the Doppler-shifted frequency. Equations 3 may be manipulated such that, except for v ˆ , all other perturbed quantities are eliminated, leaving us with the governing equation for the velocity amplitude: 2 2 2 2 2 2 (k v − )(k c − ) 2 2 A 0 v ˆ − m v ˆ = 0,m = , (4) x 0 0 2 2 2 2 2 (c + v )(k c − ) 0 A T where the Alfvén speed v and tube speed c are deﬁned as A T 2 2 2 B c v 2 0 2 0 A v = ,c = . A T 2 2 μ ρ c + v 0 0 0 A The same scheme may be applied to the exterior layers, with the consideration that in both semi-inﬁnite layers, there are no magnetic ﬁelds or ﬂows present. The governing equations Magnetoacoustic Waves and the Kelvin–Helmholtz Instability. . . Page 5 of 18 86 for the outer layers are thus 2 2 2 v ˆ − m v ˆ = 0,m = k − , for j = 1, 2, (5) x j j where the exterior sound speeds are deﬁned as c = γp /ρ . j j We ﬁnd trapped wave solutions to Equations 4 and 5. For the solutions to Equations 5 to be realistic, they need to be evanescent (i.e. all perturbations must vanish at ±∞), meaning that m > 0 is required for j = 1, 2. This yields the general solution of Equations 4 and 5 A(cosh m x + sinh m x), x < −x , 1 1 0 v ˆ (x) = B cosh m x + C sinh m x, |x|≤ x , (6) xj 0 0 0 D(cosh m x − sinh m x), x > x , 2 2 0 where A, B , C,and D are arbitrary constants. By inspection, we establish that two wave modes are allowed to propagate under the given constraints: one that is evanescent towards the centre of the slab (for m > 0), and one that is spatially oscillatory throughout the slab (for m < 0). These modes of propagation are the so-called surface and body modes, respec- tively (see, e.g. Roberts, 1981). Equation 6 is subject to boundary conditions at the interfaces, namely, the continuity of the Lagrangian displacement, and the continuity of total pressure: v ˆ (x =−x ) v ˆ (x =−x ) x1 0 x0 0 = , v ˆ (x = x ) v ˆ (x = x ) (7) x2 0 x0 0 = , [p ] = 0, [p ] = 0, T −x T x 0 0 where the total pressure is deﬁned as iρ ω ,x< −x , ⎪ m 2 2 2 iρ (k v − ) p ˆ (x) =ˆ v (x) − , |x|≤ x , (8) T 0 xj iρ ω ⎪ 2 ⎩ ,x>x . Using Equation 6 and the associated boundary conditions 7 and 8, we obtain a system of four coupled homogeneous algebraic equations ⎛ ⎞ ⎛ ⎞ ⎛ ⎞ c − s −c ω/ s ω/ 0 A 0 1 1 0 0 ⎜ ⎟ ⎜ ⎟ ⎜ ⎟ 0 c ω/ s ω/ s − c B 0 0 0 2 2 ⎜ ⎟ ⎜ ⎟ ⎜ ⎟ = , (9) ⎝ ⎠ ⎝ ⎠ ⎝ ⎠ (c − s ) − s c 0 C 0 1 1 1 0 0 0 0 0 − s − c (c − s ) D 0 0 0 0 0 2 2 2 where, for brevity, we introduced c = cosh m x , s = sinh m x ,for j = 0, 1, 2, and j j 0 j j 0 2 2 2 iρ (k v − ) iρ ω iρ ω 0 1 2 = , = , = . 0 1 2 m m m 0 1 2 86 Page 6 of 18 M. Barbulescu, R. Erdélyi For Equation 9 to have non-trivial solutions, we require the determinant of the matrix on the left-hand side to be equal to zero. Evaluating this condition, we obtain ( s − c ω/ )( c − s ω/ ) + ( c − s ω/ )( s − c ω/ ) = 0, 0 0 1 0 0 0 2 0 0 0 1 0 0 0 2 0 (10) which, after some algebra, yields the dispersion relation for magnetoacoustic waves in a steady magnetic slab embedded in an asymmetric non-magnetic environment ρ ρ 1 ρ ρ 0 0 2 0 0 2 4 2 2 2 2 2 2 2 m ω + m m k v − − m ω k v − m + m 1 2 0 1 2 0 A A ρ ρ 2 ρ ρ 1 2 1 2 (11) × tanh(m x ) + coth(m x ) = 0. 0 0 0 0 Equation 11 is a generalisation of the dispersion relations found in Nakariakov and Roberts (1995) and Allcock and Erdélyi (2017). The dispersion relation of Allcock and Erdélyi (2017) may be immediately recovered by removing the background ﬂow, i.e. setting U = 0. On the other hand, if we retain the background ﬂow but eliminate the asymmetric density proﬁle (i.e. ρ = ρ ), we recover the dispersion relation of Nakariakov and Roberts 1 2 (1995). An interesting feature of Equation 11 is that as opposed to similar results obtained by Roberts (1981) and Nakariakov and Roberts (1995), it does not factorise into two separate equations. In these studies, due to the symmetry of the environment, the disturbances may be divided into two modes of oscillation: the sausage mode, where the two interfaces oscillate in anti-phase, and the kink mode, where they oscillate in phase. The amplitude of the velocity perturbation, v ˆ (x), is an odd function in the case of a sausage mode, and an even function in the case of a kink mode. Moreover, when factorising the dispersion relation, the equation containing tanh(m x ) corresponds to the sausage mode, and the one containing coth(m x ) 0 0 0 0 to the kink mode. It has been shown by Allcock and Erdélyi (2017) that in a slab embedded in an asymmet- ric environment, there still exist two classes of modes of oscillation analogous to Roberts (1981). However, due to density asymmetry, the amplitudes of the perturbations on either side of the slab will not be equal, meaning that the eigenfunctions are neither odd nor even. This explains why Equation 11 cannot be factorised: both tanh(m x ) and coth(m x ) are 0 0 0 0 required to describe the asymmetric solutions. Allcock and Erdélyi (2017) labelled these asymmetric modes as quasi-sausage, when the perturbations are in anti-phase, and quasi- kink, if they are in phase. 3. Mode Classiﬁcation and Analytical Solutions Information about the nature of the wave solutions may be obtained from the parameters of the dispersion relation. We have already established that in order for waves to be trapped, 2 2 the exterior parameters m and m must be positive. Modes that do not meet this condi- 1 2 tion are referred to as leaky and are excluded from the analysis in the present work. We deﬁne the phase speed as c = ω/k and deduce that for modes to be trapped, they must ph satisfy max(−c , −c )<c < min(c ,c ). It is also worth noting that the sign of the phase 1 2 ph 1 2 speed, c , determines whether modes are forward or backward propagating, a positive sign ph corresponding to the former and a negative to the latter. The parameter m offers a means of classifying the solutions obtained numerically. We have already established that surface modes satisfy the condition m > 0, while body modes 0 Magnetoacoustic Waves and the Kelvin–Helmholtz Instability. . . Page 7 of 18 86 require m < 0. We may therefore categorise all solutions of Equation 11 with respect to the 2 2 2 signs of c , m , m ,and m . ph 0 1 2 Solutions that satisfy max(c ,v )< |c − U | < min(c − U ,c − U ) are fast sur- 0 A ph 0 1 0 2 0 face or body modes, depending on the sign of m , which is determined by the ordering of the characteristic speeds. Panels (c) and (d) in Figure 3 contain forward-propagating body mode solutions in this interval. However, they are absent from panels (a) and (b) because min(c ,c )< min(c + U ,v + U ). Slow body and surface modes have phase speeds 1 2 0 0 A 0 within the interval c < |c − U | < min(c ,v ) and |c − U | <c , respectively. T ph 0 0 A ph 0 T Equation 11 is, to the best of our knowledge, insoluble analytically without the use of simplifying approximations. We thus employ the assumption that the wavelength of the propagating wave solutions is much longer than the width of the slab, i.e. that kx 1. This also implies that for surface modes, m x → 0as kx → 0, and that tanh m x ≈ m x ,and 0 0 0 0 0 0 0 −1 coth m x ≈ (m x ) . The dispersion relation, Equation 11, may then be written as 0 0 0 0 ρ ρ 0 0 2 2 4 2 2 2 m ω + m m k v − 1 2 0 A ρ ρ 1 2 (12) 1 ρ ρ 1 0 0 2 2 2 2 − m ω k v − m + m m x + = 0. 0 1 2 0 0 2 ρ ρ m x 1 2 0 0 Following Roberts (1981), we look for surface mode solutions of the form 2 2 ω = ω + kx ω + O k x . (0) 0 (1) By taking the terms of order kx in Equation 12, we ﬁnd the ﬁrst-order terms in the pertur- bation expansion, and hence obtain two solutions: one for the Doppler-shifted quasi-sausage 2 2 2 mode with ≈ k c : 2 2 2 (c − c )(c + U ) T 0 2 2 2 0 T ≈ k c 1 − 2kx , (13) T 2 2 2 1/2 2 1/2 (c −(c +U ) ) (c −(c +U ) ) 2 2 2 ρ T 0 ρ T 0 0 1 0 2 (c + v )c [ + ] 0 A ρ c ρ c 1 1 2 2 and one for its companion quasi-kink mode with ω → 0as kx → 0: 2ρ 2 2 2 2 2 ω ≈ kx k v − k U . (14) A 0 ρ + ρ 1 2 2 2 2 Roberts (1981) also found a surface sausage mode solution with ω ≈ k c , but this solution no longer exists unless a single exterior sound speed c = c = c exists. 1 2 e In order to ﬁnd body mode solutions, we must be aware that our previous assumption, that m x → 0as kx → 0, no longer holds. Instead, we must ﬁnd solutions for which m x 0 0 0 0 0 2 2 2 is non-zero and ﬁnite as kx tends to zero. We are interested in solutions with ≈ k c and m < 0. From Equation 11 we obtain two solutions, one describing the behaviour of the Doppler-shifted quasi-sausage modes, 2 2 2 2 (v − (c − U ) )(c − (c − U ) ) T 0 T 0 2 2 2 2 2 A 0 ≈ k c 1 + k x , (15) T 0 2 2 2 2 c v π j 0 A and one describing the set of Doppler-shifted quasi-kink modes, 86 Page 8 of 18 M. Barbulescu, R. Erdélyi 2 2 2 2 (v − (c − U ) )(c − (c − U ) ) T 0 T 0 2 2 2 2 2 A 0 ≈ k c 1 + k x , (16) T 0 2 2 2 2 c v π (j − ) 0 A where j is any integer. In the case of a wide slab, i.e. when the slab width is much larger than the wavelength, we demonstrate that the two interfaces that delimit the slab cease interacting. We begin by taking kx 1, which implies that for surface modes, m x 1 (Roberts, 1981). In 0 0 0 this approximation, tanh(m x ) ≈ coth(m x ) ≈ 1, which, when applied to Equation 10, 0 0 0 0 provides us with two individual dispersion relations for the two interfaces 2 2 2 2 m k v − − m ω = 0, j 0 (17) for j = 1, 2. Equations 13 to 17 may be reduced to the analogous equations in Allcock and Erdélyi (2017) by setting U = 0, andtothoseinRoberts(1981)byalsoassuming that ρ = ρ = ρ . 0 1 2 e 4. Numerical Results We now ﬁnd the general solutions to the dispersion relation, Equation 11. Since, to the best of our knowledge, these cannot be obtained analytically, we employ a numerical scheme. We ﬁrst nondimensionalise all quantities with respect to the Alfvén speed, and introduce the 2 2 2 Alfvén Mach number M = U /v , the nondimensionalised sound speeds c¯ = c /v (for A 0 A j j A 2 2 2 j = 0, 1, 2), tube speed c¯ = c /v , and phase speed c¯ = c /v = ω/kv . ph ph A A T T Dispersion diagrams displaying general solutions to Equation 11 may be found in Fig- ures 3 to 5. They illustrate the behaviour of surface and body, quasi-sausage and quasi-kink modes, under the effect of a number of different ﬂow speeds. Four types of equilibrium conditions are assumed for the slab, each of which is represented in a panel in Figures 3 to 5, respectively. Panels (a) and (b) represent the case where c <c <v and the density T 0 A inside the slab is greater than that of the exterior. Panels (c) and (d) represent the case where v <c <c and the exterior densities are greater. A T 0 In order to better visualise the differences between the symmetric and asymmetric en- vironments, we have included side-by-side phase diagrams that illustrate the change in behaviour due to the break in symmetry. Thus, in every Figure, panels (a) and (c) depict symmetric exterior proﬁles, while panels (b) and (d) represent asymmetric exterior proﬁles. The imaginary part of the solutions to Equation 11 is displayed throughout Figures 3 to 6 in order to make a distinction between stable and unstable modes. Stable modes correspond to purely real solutions, while unstable modes will have a non-zero imaginary component that will act as a growth factor since we assumed that all perturbations are proportional to −i(ωt −kz) e . Figure 3 illustrates how a background ﬂow of M = 0.4 affects the phase diagrams in all four cases. We observe that this ﬂow speed has broken the symmetry between forward- and backward-propagating solutions in all cases. Moreover, new cut-off speeds at min(c ,c ) 1 2 have been introduced by the asymmetric exterior proﬁles (panels (b) and (d)). Figure 4 displays the effects of a background ﬂow of M = 0.6 in panels (a) and (b) and M = 1.0 in panels (c) and (d). These ﬂow strengths are strong enough to cause the slow body modes, which would have been backward propagating for lesser speeds, to now become forward propagating. Different ﬂow strengths are required depending on how the characteristic speeds are ordered. We note that the behaviour is identical throughout the four Magnetoacoustic Waves and the Kelvin–Helmholtz Instability. . . Page 9 of 18 86 Figure 3 Dispersion diagrams considering an interior that is dense ((a)and (b)), and one that is evacuated ((c)and (d)), including a background ﬂow of Alfvén Mach number M = 0.4. Panels (a)and (c) illustrate the solutions obtained for symmetric exterior density proﬁles, while (b)and (d) illustrate the effects of breaking this symmetry. The asymmetric density proﬁle introduces new cut-off frequencies at min(c , c ), while the 1 2 ﬂow further breaks the symmetry by causing forward- (c¯ > 0) and backward- (c¯ < 0) propagating modes ph ph to have different phase speeds. The shaded areas represent regions for which body modes propagate. The 2 2 hatched regions contain no stable trapped solutions (m < 0or m < 0). 1 2 86 Page 10 of 18 M. Barbulescu, R. Erdélyi Figure 4 Same as Figure 3, but including background ﬂows of Alfvén Mach number M = 0.6((a)and (b)), and M = 1.0((c)and (d)). The bulk ﬂow is now strong enough to have caused the backward-propa- gating slow body modes to become forward propagating. The asymmetric density proﬁle does not affect the threshold value at which this happens. Magnetoacoustic Waves and the Kelvin–Helmholtz Instability. . . Page 11 of 18 86 Figure 5 Same as Figure 3, but including background ﬂows of Alfvén Mach number M = 0.9((a)and (b)), and M = 1.4((c)and (d)). In the symmetric case ((a)and (c)), the KHI occurs for a small interval of kx . If the exterior density proﬁle is sufﬁciently asymmetric, the sausage mode becomes KH unstable for any value of kx greater than the threshold value. 0 86 Page 12 of 18 M. Barbulescu, R. Erdélyi panels, meaning that the asymmetry in the equilibrium proﬁles does not affect the change in direction with increasing M . Figure 5 illustrates the behaviour of the system subject to a ﬂow of M = 0.9 in panels (a) and (b), and M = 1.4 in panels (c) and (d), which is strong enough for instabilities to occur. We see that in the case of symmetric equilibrium proﬁles, the instability is restricted to a short range of values of kx . However, if the exterior parameters are asymmetric, the mode that was previously unstable in only that small interval is now unstable for any value of kx greater than the instability onset value. One important point to note here is that Figures 3 – 5 reinforce the results obtained in Section 3. It is readily visible that for kx 1, there exist Doppler-shifted quasi-sausage and quasi-kink surface modes with phase speeds approximately equal to c¯ + M and 0, T A respectively. Likewise, there exist body modes with c¯ ≈±c¯ + M . The quasi-sausage ph T A mode is only present in Figure 3(c) because of the ordering of the characteristic speeds, while the rest are present throughout. In Figure 6, the phase speed has been plotted with respect to M ,for kx = 0.5and A 0 two different density ratios. The values of the density ratios and of kx were selected in order to have a clear representation of the modes in the ﬁgures. However, the chosen value of kx may also be relevant to modelling instabilities on CME ﬂanks, for example. Using the values of the wavelengths of the unstable perturbation and the width of shear layer as measured by, e.g., Foullon et al. (2011), we obtain a possible range for kx between 0.131 and 0.656. Panel (a) represents a symmetric density proﬁle, panel (b) an asymmetric one, and both satisfy c <c <v . Comparing the two panels, it is immediately apparent that T 0 A by increasing ρ , both the cut-off at c¯ and the KHI threshold are lowered. It is also worth 2 2 noting that in panel (b), modes with c > min(c ,c ) may exist as long as they are unstable ph 1 2 2 2 2 since they satisfy the condition that c > min(c ,c ) and are thus trapped. There also exist ph 1 2 no fast modes because min(c¯ , c¯ )> min(c¯ + M , 1 + M ). 1 2 0 A A Figure 7 showcases the effect of having an asymmetric density proﬁle on the KHI thresh- old value. Throughout the panels, the green and red curves (plotted for ρ = ρ = ρ and 1 2 0 ρ = ρ = 2ρ , respectively) represent the symmetric density proﬁles. In the left panel, the 1 2 0 blue curve also represents a symmetric density proﬁle, corresponding to a lower density ra- tio of ρ = ρ = 0.5ρ . This panel illustrates how, for symmetric density proﬁles, the KHI 1 2 0 threshold increases, both with increasing values of kx , but also with decreasing values of the density ratios. As suggested by Equation 17, the threshold value for a wide slab tends to that of a single interface. The middle and right panels illustrate the effect of increasing asymmetry in the density ratios. Owing to the lack of interaction between the interface when kx 1, the greater density ratio will determine the threshold value. However, if kx 1, 0 0 the densities on either side will play a role. Figure 8 compares the effects of increasing density ratios in the case of symmetric (left) and asymmetric slabs (centre). In both cases, three slab widths are considered: a thin slab (red), with kx = 0.1, an “intermediate” value of kx = 1 (green), and a wide slab (blue), 0 0 with kx = 10. In the left panel, the exterior densities are assumed to be equal (ρ = ρ = 0 1 2 ρ ), while in the centre, we only assumed that ρ /ρ = 2. Theeffect of theasymmetryis e 2 0 most intense for small kx , when there is most interaction between the interfaces. The panel on the right illustrates how the wide asymmetric slab becomes unstable when the interface corresponding to the highest density ratio becomes unstable. For ρ <ρ , the threshold 1 2 corresponds to the interface with the constant density ratio (represented by the horizontal dotted line), while for ρ >ρ , the threshold values tend to that of the interface with variable 1 2 density ratio (represented by the dot-dashed curve). Magnetoacoustic Waves and the Kelvin–Helmholtz Instability. . . Page 13 of 18 86 Figure 6 Nondimensionalised phase speed c¯ plotted with respect to the Alfvén Mach number M for ph A kx = 0.5. The shaded areas represent regions where body modes propagate. The hatched regions contain 2 2 no stable trapped solutions (m < 0or m < 0), but unstable solutions may still exist because they have both 1 2 real and imaginary components. Increasing the density on just one side of the slab decreases the KH threshold and lowers cut-off speeds. Thus, there may exist ranges of M where trapped modes, which would otherwise be able to propagate, become leaky. 86 Page 14 of 18 M. Barbulescu, R. Erdélyi Figure 7 KHI threshold values of M , calculated for values of kx from0.05to2,for symmetric andasym- A 0 metric density proﬁles. The dashed lines represent the threshold values of a single interface and correspond to the density ratios of their respective colour. Figure 8 KHI threshold values of M , calculated for symmetric (left) and asymmetric density proﬁles (centre, ρ /ρ = 2). The panel on the right compares the threshold values obtained for the wide asymmetric 2 0 slab to that of two non-interacting interfaces. The dotted horizontal line and the dot-dashed curve represent the threshold values for the interfaces with constant and variable density ratios, respectively. 5. Applications In the previous sections, we have derived the dispersion relation for a steady slab embedded in an asymmetric environment and obtained approximate and general solutions. We now wish to discuss possible applications of this model and how it compares to previous formu- lations. We primarily focus on the observations described in Foullon et al. (2011)ofaKHI at a CME ﬂank. The event observed by the Atmospheric Imaging Assembly on board the Solar Dynamics Observatory on November 3, 2010, was described as a series of Kelvin– Helmholtz vortices propagating on the ﬂank of a CME. The region including the ﬂank may Magnetoacoustic Waves and the Kelvin–Helmholtz Instability. . . Page 15 of 18 86 Figure 9 KHI detected on the ﬂank of the CME is displayed on the left.The box on the right is a schematic representation of the unstable region. For more details about the spatial and temporal evolution of this event, see Foullon et al. (2011). be interpreted as a three-layer waveguide, with the dense CME core on one side, the CME ﬂank in the middle, and the low-density solar corona on the other side, as in Figure 9.Since the core ejecta is much slower than the ﬂank on the timescale of the instability, it is reason- able to approximate it as being static. Using the parameters measured by Foullon et al. (2011) and Equation 11,wewishto estimate the densities of the CME core and ﬂank in relation to the coronal background −1 density. We begin by assuming a background Alfvén speed v = 800 km s , and sound speed c = 0.6v . The speeds of the ejecta ﬂow and of the perturbations at the interface 0 A −1 −1 were measured to be U = 833 ± 5kms and c = 417 ± 7kms , respectively. Using 0 ph these values, we calculate the Alfvén Mach number of the ﬂow, M ≈ 1.05, and the non- dimensionalised phase speed, c¯ = 0.521 ± 0.009. The wavenumber is measured as k ≈ ph 0.35 Mm, and the width of the shear layer is estimated to be 2x ≈ 2.25 ± 1.5 Mm, making kx ≈ 0.394 ± 0.263. Since kx < 1, there will be interactions between the boundaries of the 0 0 shear layer, meaning that the density asymmetry will play an important role in the formation of the KHI. Before we start our analysis, we must ﬁrst note that when the density contrast between the three regions is such that min{¯ c , c¯ } < min{¯ c + M , 1 + M }, there exist no trapped 1 2 0 A A fast modes (as is the case in Figure 6). Since we expect this to be the case, we immediately discount the fast modes. We interpret the observation as that of a slow kink mode propagat- ing along a highly asymmetric steady slab. It has been shown by Allcock and Erdélyi (2017) that for both slow and fast modes, the transverse component of the displacement is highly sensitive to the density asymmetry. The slow mode interpretation is therefore reasonable, even though one would expect little transverse displacement in the low-beta coronal plasma. The results of the numerical analysis are presented in Figure 10, where we assumed −6 density ratios of ρ /ρ = 1.7, and ρ /ρ = 10 .For M ≈ 1.05, we obtain c¯ = 0.526, 1 0 2 0 A ph which matches the observed phase speed estimate of c¯ = 0.521 +±0.009. The growth ph −1 rate of the instability, i.e. the imaginary part of ω, is calculated to be γ ≈ 0.023 s ,which −1 compares reasonably well with the observed growth rate of γ = 0.05 ± 0.03 s . We note that the choice of density ratios is signiﬁcantly more sensitive on the interface separating the core from the ﬂank. We were able to obtain values of c¯ and γ in close ph 86 Page 16 of 18 M. Barbulescu, R. Erdélyi −6 Figure 10 Slow kink mode plotted for c = 0.6v , ρ /ρ = 1.7, ρ /ρ = 10 ,and kx = 0.5. The upper 0 A 1 0 2 0 0 and lower panels contain the real and imaginary parts of the non-dimensionalised phase speed, respectively. agreement with the observations for values of ρ /ρ in the range (1.6, 1.8). On the other 1 0 −3 −6 hand, ρ /ρ may be as high as 10 , with values lower than 10 having very little further 2 0 effect. Our model is therefore in good agreement with the observations and estimates the density of the CME ejecta to be at least six orders of magnitude higher than the background −12 −3 coronal density of ≈ 10 kg m . Our interpretation is signiﬁcantly more accurate than one obtained by means of a single interface model. In such a model, one would have to assume an unrealistically low Alfvén speed in order to match the observed phase speed with a high-density contrast. Otherwise, −1 assuming a realistic Alfvén speed v = 800 km s would yield a density ratio of ρ /ρ = A 2 0 1/3 between the ﬂank and the corona, which would signiﬁcantly underestimate the density of the CME. Similarly, the high-density contrast could also not be obtained from a model of a slab in a symmetric environment. One limitation of our model is that it does not adequately explain the absence of the KHI on the inner interface, between the core and the ﬂank. It is likely that the core is permeated by a strong magnetic ﬁeld that inhibits the formation of the instability. This effect would have to be included in a more realistic interpretation. 6. Conclusion The goals of the present work are twofold: to study the effects of a steady ﬂow on the propagation of magnetoacoustic waves in a magnetic slab in an asymmetric environment, and to examine the effects of the asymmetry on the condition for occurrence of the KHI. Magnetoacoustic Waves and the Kelvin–Helmholtz Instability. . . Page 17 of 18 86 In order to accomplish this, we solved the dispersion relation, Equation 11, using analytical approximations and numerical schemes. Since our analysis is only concerned with trapped mode solutions, we ﬁrst obtained nec- essary and sufﬁcient conditions for their existence. We then classiﬁed them as surface or body, quasi-sausage or quasi-kink modes, and obtained analytical solutions using the thin slab (Equations 13 – 16), and wide slab approximations (Equation 17). Numerical solutions of the dispersion relation Equation 11, plotted in terms of the nondi- mensionalised wavenumber kx for speciﬁc values of the Alfvén Mach number, M ,are 0 A presented in Figures 3 to 5. The ﬂow causes the symmetry between forward-propagating (ω/k > 0) and backward-propagating (ω/k < 0) modes to break, causing various modes to no longer be trapped. Furthermore, it causes backward-propagating modes to become for- ward propagating after some threshold value particular to the mode. Finally, ﬂow speeds past a critical value will cause the KHI to occur. In terms of the solutions to Equation 11, this occurs when ω < 0. The imaginary part of the solution acts as the growth rate in the time evolution of the wave, causing it to steepen (see panels (b) and (c) in Figure 1). We wish to establish the qualitative effects of the asymmetry on the KHI in order to gen- eralise the results of Allcock and Erdélyi (2017) on wave propagation. The authors found that asymmetry in the density proﬁle asymmetrically modiﬁes the amplitudes of the quasi- sausage and quasi-kink modes differently. In a symmetric slab, these modes would have anti- symmetric and symmetric amplitudes about the z-axis, respectively. However, the asymmet- ric density proﬁle causes the quasi-sausage mode to increase in amplitude about the interface separating the interior from the lower density region, and decrease in amplitude about the other. The converse is true for the quasi-kink mode. Considering the above, we hypothesise that for highly asymmetric density proﬁles and for intermediate or high values of kx , the slab may become asymmetrically unstable. If this is true, a quasi-sausage wave should trigger the KHI at the boundary separating the sparser region from the interior, while the converse should be true for the quasi-kink. However, a more detailed analysis of the eigenfunctions is required in order to conﬁrm this. Highly asymmetric systems, such as the CME ﬂank in Foullon et al. (2011), are likely prone to KHIs as long as the boundaries of the slab interact. In that example, the low density of the corona stabilises the CME ﬂank, while the high-density core destabilises it, and we observe the KHI. Because this conﬁguration of CMEs is not uncommon, we suggest that the limited number of observations are not indicative of the number of instances of the KHI in these phenomena. Further study is needed in order to determine its prevalence. Applications of this model are in no way limited to CME ﬂanks, even though they received much atten- tion in this study. Any analysis of a steady conﬁguration, whether solar or magnetospheric, that may be approximated by a slab geometry, would likely beneﬁt from the inclusion of asymmetry. Acknowledgements The authors thank the anonymous referee, J.F. Mather, and M. Allcock for their valu- able suggestions and comments. The authors also thank C.J. Nelson for his help in obtaining one of the ﬁgures, and are grateful to the Science and Technology Facilities Council (STFC, grant number ST/M000826/1) UK and the Royal Society (UK) for the support received. Data in Figure 9 are courtesy of NASA/SDO and the AIA, EVE, and HMI science teams. Disclosure of Potential Conﬂicts of Interest The authors declare that they have no conﬂicts of interest. Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 Inter- national License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. 86 Page 18 of 18 M. Barbulescu, R. Erdélyi References Allcock, M., Erdélyi, R.: 2017, Magnetohydrodynamic waves in an asymmetric magnetic slab. Solar Phys. 292, 35. DOI. Berger, T.E., Slater, G., Hurlburt, N., Shine, R., Tarbell, T., Title, A., Lites, B.W., Okamoto, T.J., Ichimoto, K., Katsukawa, Y., Magara, T., Suematsu, Y., Shimizu, T.: 2010, Quiescent prominence dynamics observed with the Hinode solar optical telescope. I. Turbulent upﬂow plumes. Astrophys. J. 716, 1288. DOI. Cairns, R.A.: 1979, The role of negative energy waves in some instabilities of parallel ﬂows. J. Fluid Mech. 92(1), 1. DOI. Foullon, C., Verwichte, E., Nakariakov, V.M., Nykyri, K., Farrugia, C.J.: 2011, Magnetic Kelvin–Helmholtz instability at the Sun. Astrophys. J. Lett. 729,L8. DOI. Foullon, C., Verwichte, E., Nykyri, K., Aschwanden, M.J., Hannah, I.G.: 2013, Kelvin–Helmholtz instability of the CME reconnection outﬂow layer in the low corona. Astrophys. J. 767, 170. DOI. Goedbloed, J.P.: 2009a, New construction of the magnetohydrodynamic spectrum of stationary plasma ﬂows. I. Solution path and alternator. Phys. Plasmas 16(12), 122110. DOI. ADS. Goedbloed, J.P.: 2009b, New construction of the magnetohydrodynamic spectrum of stationary plasma ﬂows. II. Rayleigh–Taylor and Kelvin–Helmholtz instability. Phys. Plasmas 16(12), 122111. DOI. ADS. Joarder, P.S., Nakariakov, V.M., Roberts, B.: 1997, A manifestation of negative energy waves in the solar atmosphere. Solar Phys. 176, 285. DOI. Kuridze, D., Zaqarashvili, T.V., Henriques, V., Mathioudakis, M., Keenan, F.P., Hanslmeier, A.: 2016, Kelvin–Helmholtz instability in solar chromospheric jets: theory and observation. Astrophys. J. 830, 133. DOI. Möstl, U.V., Temmer, M., Veronig, A.M.: 2013, The Kelvin–Helmholtz instability at coronal mass ejection boundaries in the solar corona: observations and 2.5D MHD simulations. Astrophys. J. Lett. 766, L12. DOI. Nakariakov, V.M., Roberts, B.: 1995, Magnetosonic waves in structured atmospheres with steady ﬂows, I. Solar Phys. 159, 213. DOI Ofman, L., Thompson, B.J.: 2011, SDO/AIA observation of Kelvin–Helmholtz instability in the solar corona. Astrophys. J. Lett. 734, L11. DOI. Roberts, B.: 1981, Wave propagation in a magnetically structured atmosphere. II: Waves in a magnetic slab. Solar Phys. 69, 39. DOI. Ryutova, M.: 2015, Physics of Magnetic Flux Tubes, Astrophysics and Space Science Library, 417, Springer, Berlin. DOI. Ryutova, M., Berger, T., Frank, Z., Tarbell, T., Title, A.: 2010, Observation of plasma instabilities in quiescent prominences. Solar Phys. 267, 75. DOI. Soler, R., Terradas, J., Oliver, R., Ballester, J.L., Goossens, M.: 2010, Kelvin–Helmholtz instability in coronal magnetic ﬂux tubes due to azimuthal shear ﬂows. Astrophys. J. 712, 875. DOI. Somasundaram, K., Venkatraman, S., Sengottuvel, M.P.: 1999, Hydromagnetic surface waves along com- pressible cylindrical ﬂux tubes with steady ﬂows. Plasma Phys. Control. Fusion 41, 1421. DOI. Taroyan, Y., Erdélyi, R.: 2002, Resonant and Kelvin–Helmholtz instabilities on the magnetopause. Phys. Plasmas 9, 3121. DOI. Taroyan, Y., Ruderman, M.S.: 2011, MHD waves and instabilities in space plasma ﬂows. Space Sci. Rev. 158(2), 505. DOI. Terra-Homem, M., Erdélyi, R., Ballai, I.: 2003, Linear and non-linear MHD wave propagation in steady-state magnetic cylinders. Solar Phys. 217, 199. DOI. Tirry, W.J., Cadez, V.M., Erdelyi, R., Goossens, M.: 1998, Resonant ﬂow instability of MHD surface waves. Astron. Astrophys. 332, 786. Zaqarashvili, T.V., Vörös, Z., Zhelyazkov, I.: 2014, Kelvin–Helmholtz instability of twisted magnetic ﬂux tubes in the solar wind. Astron. Astrophys. 561, A62. DOI. Zaqarashvili, T.V., Zhelyazkov, I., Ofman, L.: 2015, Stability of rotating magnetized jets in the solar atmo- sphere. I. Kelvin–Helmholtz instability. Astrophys. J. 813, 123. DOI. Zhelyazkov, I.: 2015, On modeling the Kelvin–Helmholtz instability in solar atmosphere. J. Astrophys. As- tron. 36, 233. DOI. Zsámberger, N.K., Allcock, M., Erdélyi, R.: 2018, Magneto-acoustic waves in a magnetic slab embedded in an asymmetric magnetic environment: the effects of asymmetry. Astrophys. J. 853, 136. DOI. ADS.
Solar Physics – Springer Journals
Published: May 28, 2018
It’s your single place to instantly
discover and read the research
that matters to you.
Enjoy affordable access to
over 18 million articles from more than
15,000 peer-reviewed journals.
All for just $49/month
Query the DeepDyve database, plus search all of PubMed and Google Scholar seamlessly
Save any article or search result from DeepDyve, PubMed, and Google Scholar... all in one place.
Get unlimited, online access to over 18 million full-text articles from more than 15,000 scientific journals.
Read from thousands of the leading scholarly journals from SpringerNature, Elsevier, Wiley-Blackwell, Oxford University Press and more.
All the latest content is available, no embargo periods.
“Hi guys, I cannot tell you how much I love this resource. Incredible. I really believe you've hit the nail on the head with this site in regards to solving the research-purchase issue.”
Daniel C.
“Whoa! It’s like Spotify but for academic articles.”
@Phil_Robichaud
“I must say, @deepdyve is a fabulous solution to the independent researcher's problem of #access to #information.”
@deepthiw
“My last article couldn't be possible without the platform @deepdyve that makes journal papers cheaper.”
@JoseServera
DeepDyve Freelancer | DeepDyve Pro | |
---|---|---|
Price | FREE | $49/month |
Save searches from | ||
Create lists to | ||
Export lists, citations | ||
Read DeepDyve articles | Abstract access only | Unlimited access to over |
20 pages / month | ||
PDF Discount | 20% off | |
Read and print from thousands of top scholarly journals.
Already have an account? Log in
Bookmark this article. You can see your Bookmarks on your DeepDyve Library.
To save an article, log in first, or sign up for a DeepDyve account if you don’t already have one.
All DeepDyve websites use cookies to improve your online experience. They were placed on your computer when you launched this website. You can change your cookie settings through your browser.
ok to continue