# Modeling Kelvin-Helmholtz Instability in Soft X-Ray Solar Jets.

1. IntroductionJets are considered to be ubiquitous confined plasma ejecta in the solar atmosphere. They have been extensively observed in the solar atmosphere in various wavebands, such as Ha [1, 2], Ca ii H [3-5], EUV [6], and soft X-ray [7] in order to understand their multitemperature characteristics. X-ray jets were discovered by the Soft X-Ray Telescope (SXT) on board Yohkoh [8], as transient X-ray energy release and enhancement with apparent collimated ballistic motions of the plasma associated with the flares in X-ray bright points, emerging flux regions, or active regions (for details, see Shibata et al. [7]). As it has been pointed out by Shimojo et al. [9], jets from X-ray bright points in active regions most likely appear at the western edge of preceding sunspots and exhibit a recurrent plasma propulsion in the solar atmosphere. These X-rayjets are confined plasma dynamics with typical morphological properties, for example, (1-40) x 104 km length, and the width of 5 x [10.sup.3] -[10.sup.5] km. Such jets possess apparent velocities of 10-1000 km [s.sup.-1] and lifetime of 100-16,000 s [9]. The electron densities of the X-ray jets are of the orders of (0.7-4) x [10.sup.9] [cm.sup.-3]. Their temperatures lie in the range of 3-8 MK with an average temperature of 5.6 MK [10].

In terms of spatial location jets can be classified as polar jets [11] and active region jets [12]. A study of polar jet parameters based on Hinode XRT observations was carried out by Savcheva et al. [13] who showed that jets preferably occur inside the polar coronal holes. Culhane et al. [11] have found from Hinode's Extreme-ultraviolet Imaging Spectrometer (EIS) 40" slot observations of a polar coronal hole that jet temperature ranges from 0.4 to 5.0 MK. The jet velocities had typical values that are mostly less than the Sun's escape velocity (618 km [s.sup.-1]); therefore, in consequence most of the jets fall back in the lower solar atmosphere after their triggering. Using the XTR on Hinode, Cirtain et al. [14] conclude that X-ray jets in polar coronal holes have two distinct velocities: one near the Alfven speed (~800km [s.sup.-1]) and another near the sound speed (200 km [s.sup.-1]). Moreover, they were the first to give an evidence for the propagation of Alfven waves in solar X-ray jets. Kim et al. [15] presented the morphological and kinematic characteristics of three small-scale X-ray/EUV jets simultaneously observed by the Hinode XRT and the Transition Region and Coronal Explorer (TRACE). While observing the coronal jets, for two different wavelength bands, they obtain matching characteristics for their projected speed (90-310 km [s.sup.-1]), lifetime (100-2000 s), and size (1.1-5 x [10.sup.5] km). Chifor at al. [12] have reported 2007 January 15/16 observations of a recurring jet situated on the west side of NOAA active region 10938. A strong blue-shifted component and an indication of a weak red-shifted component at the base of the jet were observed around [T.sub.e] = 1.6 MK in these jets. The upflow velocities were observed exceeding up to 150 km [s.sup.-1]. These jets were seen over a range of temperatures between 0.25 and 2.5 MK, while their estimated electron densities lie above [10.sup.11] [cm.sup.-3] for the high-velocity upflow components. Yang et al. [16] presented simultaneous observations of three recurring jets in EUV and soft X-ray (SXR), which occurred in an active region on 2007 June 5. On comparing their morphological and kinematic properties, the authors have found that EUV and SXR jets had similar onset locations, directions, size, and terminal velocities. The three observed jets were having maximum Doppler velocities ranging from 25 to 121 km [s.sup.-1] in the Fe xii [lambda]195 line and from 115 to 232 km [s.sup.-1] in the He ii [lambda]256 line. Extensive multi-instrument observations obtained simultaneously with the SUMER spectrometer on board the Solar and Heliospheric Observatory (SoHO), with EIS and XRT on board Hinode, and with the Extreme-ultraviolet imagers (EUVI) of the SunEarth Connection Coronal and Heliospheric Investigation (SECCHI) instrument suite on board the Ahead and Behind STEREO spacecrafts were performed by Madjarska [17]. The dynamic process of X-ray jet formation and evolution has been derived in great detail. In particular, for the first time there was found spectroscopically a temperature of 12 MK (Fe xxiii 263.76 A) and density of 4 x [10.sup.10] [cm.sup.-3] in the quiet Sun. The author has clearly identified two types of upflows in which the first one was the collimated upflow along the open magnetic fields, and the second was the formation of a plasma cloud from the expelled bright point small-scale loops. Chandrashekhar et al. [18] studied the dynamics of two jets seen in a polar coronal hole with a combination of EIS andXRT/Hinode data. They found no evidence of helical motions in these events but detected a significant shift of the jet position in a direction normal to the jet axis, with a drift velocity of about 27 and 7 km [s.sup.-1], respectively.

The launch of the Solar Dynamics Observatory (SDO) [19] with the Atmospheric Imaging Assembly (AIA) [20, 21] opens a new page in observing the solar jets. Moschou et al. [22] have reported high cadence observations of solar coronal jets observed in the extreme-ultraviolet (EUV) 304 A using Atmospheric Imaging Assembly (AIA) instrument on board SDO. They registered, in fact, coronal hole jets, with speeds of 94 to 760 km [s.sup.-1] and lifetimes of the order of several tens of minutes. A detailed description of the dynamical behavior of a jet in an on-disk coronal hole observed with AIA/SDO was presented by Chandrashekhar et al. [23]. Their study reveals new evidence of plasma flows prior to the jet's initiation along the small-scale loops at the base of the jet. The authors have also found further evidence that flows along the jet consisting of multiple, quasi-periodic small-scale plasma ejections. In addition, spectroscopic analysis estimates temperature as Log 5.89 [+ or -] 0.08 K and electron densities as Log 8.75 [+ or -] 0.05 [cm.sup.-3] in the observed jet. Measured properties of the registered transverse wave have provided evidence that strong damping of the wave occurred as it propagates along the jet with speeds of ~110km [s.sup.-1]. Using the magnetoseismological inversion, observed plasma, and wave parameters, the jet's magnetic field is estimated as B = 1.21 [+ or -] 0.2 G. Recently, Sterling et al. [24] have reported high-resolution X-ray and extreme-ultraviolet observations of 20 randomly selected X-ray jets that form in coronal holes at the solar polar caps. In each jet, converse to the widely accepted emerging magnetic flux model, a miniature version of the filament eruptions that initiated coronal mass ejections drove the jet producing reconnection process. Formation of a rotating jet during the filament eruption on 2013 April 10-11 on the base of multiwavelength and multiviewpoint observations with STEREO/SECCHI/EUVI and SDO/AIA was reported by Filippov et al. [25]. The confined eruption of the filament within a null-point topology, which is also known as an Eiffel tower magnetic field configuration, forms a twisted jet after magnetic reconnection near the null point. The sign of the helicity in the jet is observed the same as that of the sign of the helicity in the filament. It is noteworthy that the untwisting motion of the reconnected magnetic field lines gives rise to the accelerating plasma along the jet.

It is well established that the magnetic reconnection at different heights in the solar atmosphere plays a key role in triggering jet-like events. The magnetic reconnection between open and closed fields (standard reconnection scenario) is one of the well-known processes of the jet's occurrence [15, 26]. The jets emerging by this kind of mechanism are known as standard jets [3, 7, 27, 28]. Some other observational and simulation studies showed that the reconnection at the magnetic null in a fan-spine magnetic topology can also trigger jet-like events [29-32]. Eruptions of small arches, filaments, and flux ropes from within this type of magnetic field configurations can be responsible for reconnection and jets' occurring. These types of jets are known as blowout jets [28, 33, 34]. Moore et al. [35] used the full-disk He ii 304 [Angstrom] movies from the Atmospheric Imaging Assembly on SDO to study the cool (T ~ [10.sup.5] K) component of X-ray jets observed in polar coronal holes by XRT. The AIA 304 [Angstrom] movies revealed that most polar X-ray jets spin as they erupt. The authors examined 54 X-ray jets that were found in polar coronal holes in XRT movies sporadically taken during the first year of continuous operation of AIA (2010 May through 2011 April). These 54 jets were big and bright enough in the XRT images to be categorized as a standardjet or as a blowout jet. From the X-ray movies, 19 of the 54 jets appeared like standard jets, 32 appeared as blowout jets, and three were ambiguous and were not falling in any category. Moore et al. [36] have studied 14 large-scale solar coronal jets observed in Sun's pole. In EUV movies from the SDO/AIA, each jet was very similar to most X-ray and EUV jets erupting in coronal holes. However, each was exceptional in that it went higher than most of the standard coronal jets. They were detected in the outer corona beyond 2.2[R.sub.[dot encircle]] in images as observed from the Solar and Heliospheric Observatory/Large Angle Spectroscopic Coronagraph (LASCO/C2 coronagraph [37]).

Schmieder et al. [38] proposed a new model between the standard and blowout models, where magnetic reconnection occurs in the bald patches around some twisted field lines, one of whose foot points is open. Pariat et al. [39] included the magnetic field inclination and photospheric field distribution and performed another 3D numerical MHD model for the two different types of jet events: standard and blowout jets. We note also that stereoscopic studies (multiple points of view) have been carried out by the EUVI/SECCHI imagers on board the twin STEREO spacecraft to estimate the expected speed, motion, and morphology of polar coronal jets [40]. To sum up, the flux emergence [41,42] and flux cancellation [43, 44] are the two main triggering processes that are known to be responsible for jets' occurrence. In few observational and primarily numerical studies the wave-induced reconnection has also been suggested as a cause for the onset of jet-like events [42, 45-47].

We consider magnetically structured X-ray solar jets as moving cylindrical magnetic flux tubes that support the excitation/propagation of various kind of magnetohydrodynamic (MHD) oscillations and waves. While in static solar atmospheric plasma the propagating MHD modes are stable, the axial motion of the flux tubes engenders a velocity jump at the tube surface which can trigger a Kelvin-Helmholtz (KH) instability. The KH instability arises at the interface of two fluid layers that move with different speeds (see, e.g., Chandrasekhar [48])--then a strong velocity shear arises near the interface between these two fluids forming a vortex sheet. This vortex sheet becomes unstable to the spiral-like perturbations at small spatial scales [49]. In cylindrical geometry, when a magnetic flux tube is axially moving, such a vortex sheet is evolved near tube's boundary and it may become unstable against KH instability provided that the tube axial velocity exceeds a critical value [50]. Further on, in the nonlinear stage of KH instability, this vortex sheet causes the conversion of the directed flow energy into turbulent energy making an energy cascade at smaller spatial scales [51].

The KH instability studying in various solar jets over the past decade arose from the fact that KH vortices were observed in solar prominences [52-54], in Sweet-Parker current sheets [55], in a coronal streamer [56], and in coronal mass ejections [57-61]. All these observations stimulated the modeling of KH instability in moving twisted magnetic flux tubes in nonmagnetic environment [62], in magnetic tubes of partially ionized plasma [63], in spicules [64-66], in photospheric tubes [67], in high-temperature and cool surges [68, 69], in dark mottles [70], at the boundary of rising coronal mass ejections [60, 61, 71, 72], in rotating, tornado-like magnetized jets [49], and in a chromospheric jet (fast disappearance of rapid red-shifted and blue-shifted excursions alongside a larger scale H[alpha] jet) [73]. A more extensive review on modeling the KH instability in solar atmosphere jets the reader canbe seen in Zhelyazkov [74] and references therein.

The first modeling of KH instability in X-ray jets was carried out by Vasheghani Farahani et al. [75] who explored transverse wave propagation along the detected by Cirtain et al. [14] coronal hole soft X-ray jets. Vasheghani Farahani et al. analyzed analytically, in the limit of thin magnetic flux tube, the dispersion relation of the kink MHD mode and have obtained that this mode is unstable against the KH instability when the critical jet velocity is equal to 4.47[v.sub.A] = 3576 km [s.sup.-1] ([v.sub.A] = 800 km [s.sup.-1] is the Alfven speed inside the jet). Numerical solving of the same dispersion relation when considering the jet and its environment as cold magnetized plasmas, carried out by Zhelyazkov [76, 77], yielded a little bit lower critical flow speed for the instability onset; namely, 4.31[v.sub.A] = 3448 km [s.sup.-1]. The lowest critical jet speed of 4.025[v.sub.A] = 3220 km [s.sup.-1] was derived by numerically solving the wave dispersion relation without any approximations, that is, treating both media as compressible plasmas. But even the latter critical jet speed is still too high for the KH instability to be detected/observed in coronal hole soft X-ray jets. The reason for obtaining such high critical speeds is the circumstance that Vasheghani Farahani et al. [75] assumed an electron number density of the order of [10.sup.8] [cm.sup.-3] and magnetic field strength of 10 G.

Here, we study the propagation of kink and higher MHD modes in standard active region soft X-ray jets, notably jets #8, #11, and #16 of Shimojo and Shibata's set of sixteen observed flares and jets [10], and have shown that with one order higher electron densities, ~[10.sup.9] [cm.sup.-3], and moderate magnetic field, ~7 G, MHD modes in high-speed jets, like these ones, can become unstable against the KH instability at accessible jets speeds, except for jet #16 which requires a higher flow velocity. In the next section we list the basic physical parameters of jet #11 and the topology of magnetic fields inside and outside the moving flux tube modeling jet and also derive the wave dispersion equation for both untwisted and twisted tubes. The physical parameters of other two jets (#8 and #16) will be provided en route in the next section. The numerical solving MHD wave dispersion relations and the discussion of the conditions under which the KH instability can develop in such moving structures are presented in Section 3. The last section summarizes the main results obtained in this article and outlooks our future studies of KH instability in more complex (rotating) solar atmosphere jets.

2. Geometry, Magnetic Field Topology, and MHD Wave Dispersion Relations

We consider the soft X-ray jet as a straight cylinder of radius a and density [[rho].sub.i] embedded in a uniform field environment with density [[rho].sub.e]. We study the propagation of MHD waves in two magnetic configurations, notably in untwisted and twisted flux tubes. For an untwisted tube magnetic fields in both media are homogeneous and directed along the z-axis of our cylindrical coordinate system (r,[phi],z): [B.sub.i] = (0, 0 , [B.sub.i]) and [B.sub.e] = (0, 0, [B.sub.e]), respectively (see Figure 1).

The magnetic field inside the twisted tube is helicoid, [B.sub.i] = (0, [B.sub.i[phi]](r), [B.sub.iz](r)), while outside the tube the magnetic field is uniform and directed along the tube axis, [B.sub.e] = (0,0, [B.sub.e]). Note that we assume a magnetic fields' equilibrium with uniform twist for which the magnetic field inside the tube is [B.sub.i] = (0, Ar, [B.sub.iz]), where A and [B.sub.iz] are constant. The parameter that characterizes the uniform magnetic field twist is the ratio [B.sub.i[phi]](a)/[B.sub.iz] [equivalent to] [epsilon], that is, [epsilon] = Aa/[B.sub.iz]. Our frame of reference for studying the wave propagation in the jet is attached to the surrounding magnetoplasma--thus [v.sub.0] = (0,0, [v.sub.0]) represents the relative j et velocity, if there is any flow in the environment. The jump of the tangential velocity at the tube boundary then initiates the magnetic KH instability onset when the jump exceeds a critical value.

Before dealing with governing MHD equations, it is necessary to specify what kind of plasma each medium is (the moving tube and its environment). As seen from the Events List in [10], the electron density in jet #11 is [n.sub.jet] = [n.sub.i] = 2.9 x [10.sup.9] [cm.sup.-3], the temperature is [T.sub.jet] = [T.sub.i] = 5.5-6.4 MK, and jet speed (in their notation) is [V.sub.jet] = 437 km [s.sup.-1]. Our choice for environment magnetic field, electron density, and temperature is [B.sub.e] = 6.7 G, [n.sub.e] = 2.6 x [10.sup.9] [cm.sub.-3], and [T.sub.e] = 2.0 MK, respectively. With [T.sub.i] = 5.5 MK, the total pressure balance equation (equality of the sum of thermal and magnetic pressures in both media), that is,

[p.sub.i] + [B.sup.2]/2[mu] = [p.sub.e] + [B.sup.2.sub.e]/2[mu], (1)

where [mu] is the magnetic permeability of vacuum, yields the following basic sound and Alfven speeds in the jet and its environment: [c.sub.si] = 275 km [s.sup.-1] and [v.sub.Ai] = 111 km [s.sup.-1] and [c.sub.se] = 166km [s.sup.-1] and [v.sub.Ae] = 286km [s.sup.-1], respectively. Accordingly, the plasma betas in both media are [[beta].sub.i] = 7.369 and [[beta].sub.e] = 0.403--this implies that, in principle, one can treat the jet as incompressible plasma and its surrounding medium as cool magnetized plasma. When studying the MHD wave propagation in the untwisted magnetic flux tube we shall use two approaches, namely, of compressible plasmas in both media and the simplified limit of incompressible and cool plasmas--a similarity of dispersion curves patterns, obtained from corresponding dispersion relations, will eventually justify the usage of the second approach in particular in the case of twisted tube. There are two important input parameters in the modeling KH instability in moving magnetic flux tubes, which are the density contrast, [[rho].sub.e]/ [[rho].sub.i] = = 0.896, and the ratio of the axial external and internal magnetic fields, [B.sub.e]/[B.sub.i] = b = 2.44. For the twisted tube the second parameter has the form [b.sub.twist] = [B.sub.e]/[B.sub.iz].

In a system of cylindrical coordinates, the equilibrium physical variables (density, fluid velocity, and pressure) are functions of the radial coordinate r only. Then, their perturbations can be Fourier-analyzed putting them proportional to exp[i(-[omega]t + m[phi] [k.sub.z]z)], where [omega] is the angular wave frequency (that, in general, can be a complex quantity), m is the mode number (a positive or negative integer), and [k.sub.z] is the axial wave number. We can eliminate all except two of the perturbations (perturbation [p.sub.tot] of the total (thermal + magnetic) pressure and the radial component [[xi].sub.r] of the Lagrangian displacement [xi]) to get the following governing equations [78]:

D d/dr = (r[[xi].sub.r]) = [C.sub.1]r[[xi].sub.r] - [C.sub.2]r[p.sub.tot]. D d[p.tot]/dr = [C.sub.3][[xi].sub.r] - [C.sub.1][p.sub.tot]. (2)

The coefficients D, [C.sub.1], [C.sub.2], and [C.sub.3] are functions of the equilibrium variables [[rho].sub.0], [B.sub.0], and [v.sub.0] and of the Doppler-shifted frequency [ohm] = [omega] - k x [v.sub.0] and have the following forms:

[mathematical expression not reproducible], (3)

where

[mathematical expression not reproducible]. (4)

Here [[omega].sub.A] is the Alfven frequency and [[omega].sub.c] is the cusp frequency; the other notation is standard.

Eliminating [[xi].sub.r] from (2), one obtains the well-known second-order ordinary differential equation [79-81]

[d.sup.2][p.sub.tot]/[dr.sup.2] + [[C.sup.3]/rD] [d/dr] (rD/[C.sub.3])]d[p.sub.tot]/dr + [[C.sup.3]/rD] [d/dr] (r[C.sub.1/[C.sub.3]) + 1/[D.sup.2]([C.sub.2][C.sub.3] - [C.sup.2.sub.1)][p.sub.tot] = 0. (5)

By means of the solutions to (5) in both media, one can find the corresponding expressions for [[xi].sub.r], and after merging these solutions, together with those for [p.sub.tot], through appropriate boundary conditions at the interface r = a, one can derive the dispersion relation of the normal modes propagating in the moving magnetic flux tube.

2.1. Dispersion Relation of MHD Modes in an Untwisted Flux Tube. In an untwisted magnetic flux tube, the coefficient [C.sub.1] = 0, while [C.sub.3] = [[rho].sub.0]D([[OMEGA].sup.2] - [[omega].sup.2.sub.A])--then (5) takes the form

[d.sup.2][p.sub.tot]/[dr.sup.2] + [1/r] [d[p.sub.tot]/dr] - ([m.sup.2.sub.0] + [m.sup.2]/[r.sup.2]) [p.sub.tot] = 0, (6)

where

[m.sup.2.sub.0] = ([[OMEGA].sup.2] - [k.sup.2.sub.z][c.sup.2.sub.s])/([[OMEGA].sup.2] - [k.sup.2.sub.z][v.sup.2.sub.A])/([c.sup.2.sub.s] + [v.sup.2.sub.A])([[OMEGA].sup.2] - [[omega].sup.2.sub.c]. (7)

The cusp frequency, [[omega].sup.c], is usually expressed via the so-called tube speed, [c.sub.T], notably [[omega].sub.c] = [k.sub.z][c.sub.T], where [82]

[C.sub.T] = [C.sub.s][V.sub.A]/[square root of [c.sup.2.sub.s] + [v.sup.2.sub.A] (8)

The solutions for [p.sub.tot] can be written in terms of modified Bessel functions: [I.sub.m]([m.sub.0i]r) inside the jet and [K.sub.m]([m.sub.0e]r) in its surrounding plasma. We note that wave attenuation coefficients, [m.sub.oi] and [m.sub.oe], in both media are calculated from (7) with replacing the sound and Alfveen speeds with the corresponding values for each medium. Recall that, in evaluating [m.sub.oe], the wave frequency is not Doppler-shifted--it is simply [omega]. By expressing the Lagrangian displacements [[xi].sub.ir] and [[xi].sub.er] in both media via the derivatives of corresponding Bessel functions and by applying the boundary conditions for continuity of the pressure perturbation [p.sub.tot] and [[xi].sub.r] across the interface, r = a, one obtains the dispersion relation of normal MHD modes propagating in a flowing compressible jet surrounded by a static compressible plasma [65, 83, 84]

[mathematical expression not reproducible], (9)

Due to the flowing plasma, the wave frequency is Doppler-shifted inside the jet. We recall that for the kink mode (m = 1) one defines the so-called kink speed [82]

[c.sub.k] = [([[rho].sub.i][v.sup.2.sub.Ai] + [[rho].sub.e][v.sup.2.sub.Ae]/[[rho].sub.i] + [[rho].sub.e]).sup.1/2] = [(1 + [B.sup.2.sub.e]/[B.sup.2.sub.i]/1 + [[rho].sub.e]/ [[rho].sub.i]).sup.1/2] [v.sub.A]i, (10)

which, as seen, is independent of sound speeds and characterizes the propagation of transverse perturbations. We will show that notably the kink mode can become unstable against KH instability.

When the jet is considered as incompressible plasma and its environment as a cool one, the MHD wave dispersion relation (9) keeps its form, but the two attenuation coefficients [m.sub.oie] become much simpler, namely,

[m.sub.oi] = [k.sub.z], [m.sub.oe] = ([k.sup.2.sub.z][v.sup.2.sub.Ae] - [[omega].sup.2]/[v.sub.Ae] (11)

respectively.

2.2. Dispersion Relation of MHD Normal Modes in a Twisted Flux Tube. Inside the tube (r [less than or equal to] a), where [B.sub.i[phi]] = Ar, the quantities [f.sub.B] and [[omega].sub.Ai] take the forms

[f.sub.B] = mA + [k.sub.z][B.sub.iZ], [[omega].sub.Ai] = mA + [k.sub.z][B.sub.iz]/[square root of [mu][[rho].sub.i] (12)

respectively. For incompressible plasma, we redefine (without the loss of generality) the coefficients D, [C.sub.1], [C.sub.2], and [C.sub.3] by dividing them by [C.sub.4] to obtain

[mathematical expression not reproducible]. (13)

Radial displacement is expressed through the total pressure perturbation as

[[xi].sub.r] = [D/[C.sub.3]] d[p.sub.tot]/dr + [C.sub.1]/[C.sub.3] [p.sub.tot] (14)

The solution to this equation obviously depends upon the magnetic field and density profile.

With aforementioned coefficients D, [C.sub.1], [C.sub.2], and [C.sub.3], evaluated for the jet's medium, (5) reduces to the modified Bessel equation

[[d.sup.2]/[dr.sup.2] + [1/4] [d/dr] - ([m.sup.2.sub.0i] + [m.sup.2]/[r.sup.2])][p.sub.tot] = 0, (15)

where

[m.sup.2.sub.01]; [k.sup.2.sub.z] [1 - 4[A.sub.2][[omega].sup.2.sub.Ai]/ [mu][[rho].sub.i][([[OMEGA].sup.2] - [[omega].sup.2.sub.Ai]).sup.2]]. (16)

The solution to (15) bounded at the tube axis is

[p.sub.tot] (r [less than or equal to] a) = [[alpha].sub.i][I.sub.m] ([m.sub.0i]r), (17)

where [I.sub.m] is the modified Bessel function of order m and [[alpha].sub.i] is a constant. Lagrangian displacement [[xi].sub.ir], by using (14) can be written as

[mathematical expression not reproducible] (18)

where the prime sign means a differentiation with respect the Bessel function argument.

For the cool environment with a straight-line magnetic field [B.sub.ez] - [B.sub.e] and homogeneous density [[rho].sub.e], the [C.sub.1-3] and D coefficients take the form

[mathematical expression not reproducible]. (19)

The total pressure perturbation outside the tube obeys the same Bessel equation as (15), but [m.sup.2.sub.0i] is replaced by

[m.sup.2.sub.0e] = [k.sup.2.sub.z] (l - [[omega].sup.2]/[[omega].sup.2.sub.Ae]), (20)

which coincides with the attenuation coefficient in the cool environment of an untwisted flux tube. The solution bounded at infinity now is

[p.sub.tot] (r > a) = [[alpha].sub.e][K.sub.m] ([m.sub.0e]r), (21)

where [K.sub.m] is the modified Bessel function of order m and [[alpha].sub.e] is a constant.

In this case, the Lagrangian displacement can be written as

[[xi].sub.er] = [[[alpha].sub.e]/r] [m.sub.0e]r [K'.sub.m]([m.sub.0e]r)/[[rho].sub.e]([[omega].sup.2] - [[omega].sup.2.sub.Ae]), (22)

and the Alfven frequency is simplified to

[[omega].sub.Ae] = [k.sub.z][B.sub.ez]/[square root of ([mu][[rho].sub.e]) = [k.sub.z][v.sub.Ae], (23)

Here, [v.sub.Ae] - [B.sub.e]/[square root of ([mu][[rho].sub.e]) is the Alfven speed in the surrounding magnetized plasma.

The boundary conditions which merge the solutions of the Lagrangian displacement and total pressure perturbation inside and outside the twisted magnetic flux tube have the forms [85]

[[xi].sub.ir]|[sub.r=a] = [[xi].sub.er]|[sub.r=a], [p.sub.toti] - [B.sup.2.sub.i[phi]]/[mu]a[[xi].sub.ir]|[sub.r=a] = [p.sub.tote]|[sub.r=a], (24)

where total pressure perturbations [p.sub.toti] and [p.sub.tote] are given by (17) and (21), respectively. With the help of these boundary conditions we derive the dispersion relation of the normal MHD modes propagating along a twisted magnetic flux tube with axial mass flow [v.sub.0]

[mathematical expression not reproducible], (25)

where, as we already mentioned, [OMEGA] = [omega] - k x [v.sub.0] is the Doppler-shifted wave frequency in the moving medium,

[F.sub.m]([m.sub.0i]a)[I'.sub.m]([m.sub.0i]a)/[I.sub.m]([m.sub.0i]a)/[K.sub.m] ([m.sub.0e]a). (26)

This dispersion equation is similar to the dispersion equation of normal MHD modes in a twisted flux tube surrounded by incompressible plasma [65]--the only difference is that, there in (10), [K.sub.e] = [m.sub.0e] - [k.sub.z].

3. Numerical Solutions and Results

Firstly we shall study the dispersion characteristics of the kink (m - 1) mode in untwisted moving magnetic flux tube (in two approaches, notably (i) considering the jet and its surrounding plasma as compressible media and (ii) treating the jet as an incompressible medium while its environment is assumed to be cool plasma) and, later on, explore the same thing for the kink (m - 1) and higher (m > 1) MHD modes propagating in a twisted moving flux tube. Two dispersion relations (9) and (25) are transcendent equations in which the wave frequency, [omega], is a complex quantity: Re([omega]) + Im([omega]), while the axial wave number, [k.sub.z], is a real variable. The appearance of KH instability (i.e., Im([omega]) > 0) is determined primarily by the jet velocity and in searching for a critical/threshold value of it, we will gradually change velocity magnitude, [v.sub.0], from zero to that critical value (and beyond). Our numerical task is to solve the dispersion relation in complex variables, obtaining the real and imaginary parts of the wave frequency, or as is usually done, of the wave phase velocity [v.sub.ph] = [omega]/[k.sub.z], as functions of the axial wave number, [k.sub.z], at various magnitudes of the velocity shear between the soft X-ray jet and its environment, [v.sub.0].

For numerical solving the wave dispersion relations, it is practically to normalize all variables and we do that normalizing the wave length, [lambda] = 2[pi]/[k.sub.z], to the tube radius, a, that implies a dimensionless wave number [k.sub.z]a. All speeds are normalized with respect to the Alfven speed inside the jet, [v.sub.Ai]. For normalizing the Alfven speed in the ambient coronal plasma, [v.sub.Ae], we need the density contrast, [eta], and the ratio of the magnetic fields b = [B.sub.e]/[B.sub.i], to find [v.sub.Ae]/[v.sub.Ai] = b/ [square root of [eta]. The normalization of sound speeds in both media requires the specification of the reduced plasma betas, [[??].sub.ie] = [c.sup.2.sub.si,e]/[v.sup.2.sub.Ai,e]. In the dimensionless computations, the flow speed, [v.sub.0], will be presented by the Alfven Mach number [M.sub.A] = [v.sub.0]/[v.sub.Ai].

3.1. Kelvin-Helmholtz Instability in Untwisted Flux Tubes. Among the various MHD wave spectra that exist in a static magnetic field flux tube of compressible plasma, surrounded by compressible medium, the most interesting for us is the kink-speed wave whose speed for jet #11, according to (10), is equal to

[c.sub.k] = [square root of 1 + [b.sub.2]/1 + [eta] [v.sub.Ai] = 212.7 km [s.sup.-1] (27)

or in dimensionless form [c.sub.k]/[v.sub.Ai] = 1.9166.

The normal modes propagating in a static homogeneously magnetized flux tube can be pure surface waves, pseudo-surface (body) waves, or leaky waves (see Cally [86]). The type of the wave crucially depends on the ordering of the basic speeds in both media (the flux tube and its surrounding plasma), more specifically of sound and Alfven speeds as well as corresponding tube speeds. In our case the ordering is

[V.sub.Ai] < [c.sub.se] < [c.sub.si] > [v.sub.Ae] (28)

which with [c.sub.Ti] = 103km [s.sup.-1] and [c.sub.Te] = 143.6km [s.sup.-1] after normalizing them with respect to the external sound speed, [c.sub.se], yields (in Cally's notation)

[mathematical expression not reproducible]. (29)

Since A < C and A < [C.sub.Te] along with A < [C.sub.Te] < C, according to Cally's classification, the kink mode in a rest flux tube must be pure surface mode type [S.sup.-.sub.+] (for detail, see Cally [86]). For an S- type wave V should lie between A and C: A [less than or equal to] V [less than or equal to] C, and also V < [C.sub.Te]. In our normalization, the normalized wave velocity should be bracketed between 1 and [square root of [[??].sub.i] = 2.478, that is, 1 [less than or equal to] [omega] /([k.sub.z][v.sub.Ai]) [less than or equal to] 2.478. The typical wave dimensionless wave velocity [omega]/([k.sub.z][v.sub.Ai]) is the normalized kink speed [C.sub.k]/[v.sub.Ai] = 1.9166, which, as expected, lies between 1 and 2.478. Concerning the relation between normalized wave phase velocity and external tube speed, in our case we have the opposite inequality; that is, V > [C.sub.Te], that in Cally's normalization reads as 1.2805 > 0.8649 and in ours as 1.9166 > 1.2937. The reason for that unexpected change in the mutual relations between two velocities is the relatively bigger value of parameter b (=2.44) that yields a larger magnitude of [C.sub.k], than, for example, in the case when b is close to 1--then the inequality V < [C.sub.Te] is satisfied (see, e.g., Zhelyazkov et al. [87]).

A reasonable question is how the flow will change the dispersion characteristics of the kink (m = 1) mode. Calculations show that the flow shifts upwards the kink-speed dispersion curve and splits it into two separate curves [65]. (A similar duplication is observed for the tube-speed [v.sub.Ti] dispersion curves, too.) The evolution of the pair of kink-speed dispersion curves can be seen in Figure 2(a)--the input parameters for solving the dispersion equation (9) of the kink mode (m = 1) traveling in a moving flux tube of compressible plasma surrounded by compressible coronal medium are [eta] = 0.896, b = 2.44, [[??].sub.i] = 6.141, and [[??].sub.e] = 0.336. During calculations the Alfven Mach number, [M.sub.A], was varied from zero (static plasma) to values at which we obtained unstable solutions. Let us first note that at [M.sup.A] = 0 we get the kink-speed dispersion curve which at a very small dimensionless wave number [k.sub.z]a (=0.005) yields the value of the normalized wave phase velocity equal to 1.9168--very close to the previously estimated magnitude of 1.9166. This observation implies that our code for solving the wave dispersion relation is correct. For small Alfven Mach numbers the pair of kink-speed modes travel with velocities [M.sub.A] [bar.+] [C.sub.k]/[v.sub.Ai] [65] (in that case, their dispersion curves go practically parallel). At higher MA, however, the behavior of each curve of the pair [M.sub.A] [bar.+] [C.sub.k]/[v.sub.Ai] turns out to be completely different. As seen from Figure 2(a), for [M.sub.A] [greater than or equal to] 3.785 both kink curves break and merge forming a family of semiclosed dispersion curves. A further increasing in [M.sub.A] leads to a separation of those curves in opposite directions. This behavior of the kink (m = 1) mode dispersion curves signals us that we are in a range of the Re([v.sub.ph]/[v.sub.Ai])-[k.sub.z]a-plane where one can expect the occurrence of KH instability. A prediction of the Alfven Mach number [M.sub.A] at which the instability will occur can be found from the inequality [88]

[absolute value of m][M.sup.2.sub.A] > (1 + 1/[eta])([absolute value of m][b.sup.2] +1) (30)

that gives for the kink (m = 1) mode [M.sub.A] > 3.839. In our case the KH instability starts at that [M.sub.A], for which the left-hand side semiclosed curve disappears--this happens at [M.sup.cr.sub.A] = 3.865. The growth rates of unstable kink modes are plotted in Figure 2(b). The red curves in all diagrams denote the marginal dispersion/growth rate curves: for values of the Alfven Mach number smaller than [M.sup.cr.sub.A] the kink mode is stable; otherwise it becomes unstable and the instability is of the KH type. With [M.sup.cr.sub.A] = 3.865 the kink mode will be unstable when the velocity of the moving flux tube is higher than 429 km [s.sup.-1]--a speed, which is below than the observationally measured jet speed of 437 km [s.sup.-1]. We would like to underline that all the stable kink modes in the untwisted ([epsilon] = 0) moving tube modeling jet #11 are pure surface modes while the unstable ones are not--the latter becomes partly surface and partly leaky modes (their external attenuation coefficients, [m.sub.0e]s, are complex quantities with positive imaginary parts). This circumstance means that wave energy is radiated outward in the surrounding medium, which allows us to claim the KH instability plays a dual role: once in its nonlinear stage the instability can trigger wave turbulence and simultaneously the propagating KH-mode is radiating its energy outside.

When numerically solving (9) for the kink (m =1) mode treating the jet's plasma as incompressible medium and its environment as cool plasma, we obtain dispersion curves' and growth rates' patterns very similar to those shown in Figure 2 (see Figure 3). Computations show, as expected, that the stable kink mode is a nonleaky surface mode, while the unstable one possesses a real internal attenuation coefficient, [m.sub.0i] = [k.sub.z], not changed by the instability, but the external one, [m.sub.0e], becomes a complex quantity with positive real and negative imaginary parts. Now the threshold Alfven Mach number is a little bit bigger and yields a critical jet speed for the instability onset of [congruent to] 442km [s.sup.-1], being with 5 km [s.sup.-1] higher than the observationally measured 437 km [s.sup.-1]. We note that, as a rule, jet's incompressible plasma approximation yields slightly higher threshold Alfven Mach numbers than the model of compressible media. This ascertainment shows how sensitive to jet's and its environment's media treatment the occurrence of KH instability of the kink mode in our jet is--one can expect instability onset at an accessible jet speed only if both media are considered as compressible magnetized plasmas.

The basic physical parameters of jet #8, namely [n.sub.i] = 4.1 x [10.sup.9] [cm.sup.-3] and [T.sub.i] = 5.2 MK, for a low density contrast of [eta] = 0.963 assuming as before that the temperature of surrounding plasma is [T.sub.e] = 2.0 MK and the background magnetic field, [B.sub.e], is equal to 7 G, yield the following sound and Alfven speeds: [c.sub.se] = 166 km [s.sup.-1], [c.sub.si] = 267.5 km [s.sup.-1], [v.sub.Ae] = 242.8 km [s.sup.-1], and [v.sub.Ai] = 52.22 km [s.sup.-1]. Accordingly, plasma betas in both media are [beta] = 31.5 and [[beta].sub.e] = 0.56, while the magnetic fields ratio b = 4.563 defines [B.sub.i] = 1.53 G. The speed ordering, [v.sub.A]i < [c.sub.se] < [v.sub.Ae] < [c.sub.si], according to Cally's [86] classification, tells us that in a rest magnetic flux tube the propagating wave must be a pure surface mode of type [S.sup.-.sub.+]. The numerical computations confirm this, as well as reproduce the normalized kink speed of 3.3342 within three places after the decimal point. With the flow inclusion, the behavior of the pair kink-speed dispersion curves in the region of the instability onset, as seen from Figure 4(a), is completely different; more specifically, the lower semiclosed kink-speed dispersion curves at [M.sub.A] = 6.625 and 6.65 (the orange curve) correspond to pseudosurface (body) waves, while the higher kink-speed dispersion curves are associated with surface waves. Note that at [M.sub.A] = 6.668 the surface wave dispersion curves brakes in two parts: the right-hand one merges with the lower kink-speed curve at [k.sub.z]a = 0.2337 (see the blue line) while its left-hand side forms a narrow semiclosed dispersion curve of pure surface mode. The threshold Alfveen Mach number is equal to 6.6704 (the prediction one according to (30) is 6.67) which implies that the critical flow speed for KH instability onset is equal to 348.3km [s.sup.-1]--a value far below the jet speed of 385km [s.sup.-1]. This circumstance guarantees that even if one considers the jet medium as incompressible plasma and its environment as cool plasma we will get an accessible jet's speed for instability onset.

In a similar way, calculating the sound and Alfven speeds for jet #16 and its environment (with [n.sub.i] = 1.0 x [10.sub.9] [cm.sup.-3], [T.sub.i] = 7.4 MK, [eta] = 0.95, [T.sub.e] = 2.0 MK, and Be = 6.7 G) we get [c.sub.se] = 166km [s.sup.-1], [c.sub.si] [congruent to] 319km [s.sup.-1], [v.sub.Ae] [congruent to] 474km [s.sup.-1], [v.sub.Ai] = 350km [s.sup.-1], and plasma betas [[beta].sub.i] [congruent to] 1.0 and [[beta].sub.e] = 0.14. The magnetic fields ratio, b, is equal to 1.32 giving [[beta].sub.i] 5.1 G. In this case, the speed ordering is [c.sub.se] < [c.sub.si] < [v.sub.Ai] < [v.sub.Ae] and the kink (m = 1) mode propagating in a rest magnetic flux tube is a pseudosurface (body) wave of type [B.sup.-.sub.+] (see Table I in [86]). The numerical calculations confirm this mode type and yield a normalized value of the kink speed very close to the predicted one of 1.1858. Now the pattern of stable kink-speed dispersion curves is more complicated (see Figure 5(a))--while the regular lower kink-speed curves correspond to pseudosurface (body) waves, the higher kink-speed dispersion curves have the form of narrow and long semiclosed up to [M.sub.A] = 2.3 loops and are in fact bulk waves: both attenuation coefficients are purely imaginary numbers. As expected, the dispersion curves of unstable kink mode possess complex attenuation coefficients (with positive real and imaginary parts for [m.sub.0i] and positive real and negative imaginary part for [m.sub.0e]). The threshold Alfven Mach number being equal to 2.359305 determines a critical flow velocity of 825.8km [s.sup.-1] for KH instability onset which is much higher than the jet #16 speed of 532 km [s.sup.-1]. This consideration shows that even a relatively high-speed observed soft X-ray jet cannot become unstable--the reason for that in this case is the low density of the jet's plasma.

3.2. Kelvin-Helmholtz Instability in Twisted Flux Tubes. A natural question that immediately raises is how the twist of the internal magnetic field, [B.sub.i], will change the conditions for instability occurrence and the shape of the dispersion curves. When we model the soft X-ray jet as a moving twisted magnetic flux tubes, there are two types of instabilities that can develop in the jet: namely, kink instability due to the twist of the magnetic field and KH instability owing to the tangential discontinuity of plasma velocity at the tube boundary. According to Dungey and Loughhead [89], the kink instability will occur provided [B.sub.i[phi]] (a) > 2[B.sub.iz], or in our notation, as [epsilon] > 2. Further on, in order to prevent the developing of kink instability, we will consider weakly twisted flux tubes, [epsilon] < 1, and, thus, only KH instability can take place in such a jet configuration. When numerically solving dispersion equation (25) we need the normalization of the local Alfven frequencies [omega.sub.Aie]. The normalization of these frequencies is performed by multiplying each of them by the tube radius, a, and dividing by the Alfven speed [v.sub.Ai] = [B.sub.iz/[square root of [mu][[rho].sub.i] to get

a[omega.sub.Ai]/[V.sub.Ai] = m Aa/[B.sub.iz] + [k.sub.z]a = m[epsilon] + [k.sub.z]a, a[omega.sub.Ae]/[V.sub.Ai] = [k.sub.z]a [B.sub.e]/[B.sub.iz]/[square root of [eta]] = [k.sub.z]a [b.sub.twist]/[square root of [eta]], (31)

where [b.sub.twist] = [B.sub.e]/[B.sub.iz] = b[square root of (1 + [[epsilon].sup.2])]. For small values of [epsilon], we will take [b.sub.twist] [congruent to] b.

We begin our calculations for the kink (m =1) mode in the moving twisted tube (modeling jet #11) with [eta] = 0.896, [b.sub.twist] [congruent to] b = 2.44, and [epsilon] = 0.025. The results of numerical task are presented in Figure 6, from which one sees that the kink mode dispersion and growth rate curves are very similar to those shown in Figure 3. The critical flow velocity for emerging KH instability now is [v.sup.cr.sub.0] = 442.6 km [s.sup.-1] calculated by using the "reduced" Alfven speed [V.sub.Ai]/[square root of (1 + [[epsilon].sup.2]) = 110.96km [s.sup.-1]. KH instability of the (m = -1) mode will occur if the jet speed exceeds 442 km [s.sup.-1], of which value is beyond the observationally derived jet speed of 437km [s.sup.-1]. We note that while the wave attenuation coefficient of the m = 1 mode propagating in untwisted moving magnetic flux tube of incompressible plasma surrounded by a cool medium is not changed by the flow, in the twisted flux tube (in the same approximations) that attenuation coefficient becomes a complex number with positive imaginary part, like the attenuation coefficient in the cool environment.

According to the instability criterion (30), the occurrence of KH instability of higher MHD modes would require lower threshold Alfven Mach numbers (and accordingly lower critical jet speeds) than those of the kink (m = [+ or -]1) mode. Our computations show, however, that, with [epsilon] = 0.025, the threshold Alfveen Mach numbers are generally still high to initiate the KH instability. Their magnitudes for the m = 2, 3, and 4 modes are 3.965, 3.943, and 3.926, respectively. The corresponding critical velocities for instability onset are accordingly equal to 440.0, 437.5, and 435.6km [s.sup.-1]. As seen, only the m = 4 MHD mode can become unstable against KH instability. A clutch of unstable dispersion curves and normalized growth rates are shown in Figure 7. Here we observe a new phenomenon: the KH instability starts at some critical [k.sub.za]-number along with the corresponding threshold Alfven Mach number. That critical dimensionless wave number for the m = 4 MHD mode (look at Figure 7) is equal to 0.528. If we assume that jet #11 has a width [DELTA]l = 5x[l0.sup.3] km, then the critical wavelength for a KH instability emergence is [[lambda].sup.m=4.sub.cr] [congruent to] 29.8 Mm.

A substantial decrease in the threshold Alfven Mach number in moving twisted magnetic flux tubes can be obtained by increasing the magnetic field twist parameter, say taking, for instance, [epsilon] = 0.4. The dispersion curves and normalized wave growth rates of the kink (m = [+ or -]1) mode are plotted in Figures 8 and 9. One is immediately seeing the rather complicated forms of both the dispersion curves and dimensionless growth rates of the m = -1 kink mode--such a complication was absent at [epsilon] = 0.025. We note that computations were performed at [b.sub.twist] = b/[square root of (1 + [[epsilon].sup.2]], and the reference "reduced" Alfven speed inside the jet now is [v.sub.Ai]/[square root of (1 + [[epsilon].sup.2]] = 103 km [s.sup.-1]. With this Alfven speed the critical velocities for the instability occurrence are equal to 416km [s.sup.-1] (for the m =1 mode) and [congruent to] 408 km [s.sup.-1] for the m = -1 mode, respectively.

In solving dispersion equation (25) for the m = 2 MHD mode, we start with a little bit bigger Alfven Mach number than that predicted by instability criterion (30) and have obtained three different kinds of dispersion curves and normalized growth rates (see Figure 10), notably curves with distinctive change of the normalized wave phase velocity in a relatively narrow [k.sub.z]a-interval (the orange and blue curves in Figure 10(a)), two piece-wise dispersion curves (the purple and green curves in the same plot), and one almost linear dispersion curve (the red one which is in fact the marginal dispersion curve obtained for the threshold Alfven Mach number equal to 4.062). Not less interesting are the dimensionless growth rate curves corresponding to these three kinds of dispersion curves. For the threshold Alfveen Mach number of 4.062 we obtained one instability window and instability starts at the critical dimensionless wave number of 1.944, or equivalently at [[lambda].sup.m=4.sub.cr] [congruent to] 8.0 Mm. For the slightly superthreshold Alfven Mach numbers of 4.075 and 4.1 we got two different instability windows (see the green and purple curves in Figure 10(b)) and a further increase in Ma leads to merging of those two instability windows--see the blue and orange curves in the same plot. The critical flow velocity at which KH instability arises is equal to [congruent to] 418 km [s.sup.-1].

We have skipped the complicated dispersion and growth rate curves pictured in Figure 10 for the m = 3 and m = 4 MHD modes--we have calculated only a few curves for Alfven Mach numbers close to the corresponding threshold ones (see Figures 11 and 12). One can observe that for the magnetic field twist parameter [epsilon] = 0.4 the critical dimensionless wave numbers are shifted far on the right--their values for both modes are equal to 3.619 and 4.512, respectively, that yield [[lambda].sup.m=3.sub.cr] [congruent to] 4.3 Mm and [[lambda].sup.m=4.sub.cr] [congruent to] 3.5 Mm. The critical jet speeds at which KH instability emerges are equal correspondingly to [congruent to] 422 and 424 km [s.sup.-1]. We note that all the threshold Alfveen Mach numbers of higher MHD modes (m [greater than or equal to] 2) traveling on the moving twisted magnetic flux tube are larger than the predicted ones, but while for [epsilon] = 0.025 they (Alfven Mach numbers) are decreasing with increasing the mode number, at [epsilon] = 0.4 we have just the opposite ordering. We would like to notice that finding the marginal dispersion and growth rate curves of the higher modes at [epsilon] = 0.4 turns out to be a laborious computational task.

To finish our survey on KH instability of MHD modes in moving twisted magnetic flux tubes, we will briefly consider how a weak twist of the internal magnetic field, [epsilon] = 0.025, will change the critical jet speeds for the occurrence of KH instability in jet #8. We are not going to graphically present the dispersion curves and growth rates of unstable modes because simply there is nothing special in their shape--in the incompressible plasma jet approximation and cool environment at small magnetic field twist the curves look, more or less, similar. Recall that at density contrast of 0.963 and magnetic fields ratio b = 4.56, the threshold Alfven Mach number for the kink (m = 1) mode is equal to 6.9047, which for the reference Alfven speed of 52.22 km [s.sup.-1] yields a critical flow velocity of 360.6 km [s.sup.-1] that is lower than the jet speed of 385km [s.sup.-1]. The KH instability of m = 2, 3, and 4 MHD modes should occur at flow velocities of [congruent to] 259,357, and 355.5 km [s.sup.-1] and start at critical dimensionless wave numbers [([k.sub.z]a).sub.cr] equal, respectively, to 0.054, 0.119, and 0.205.

4. Summary and Conclusion

In this paper, we have explored the conditions under which MHD waves with various mode numbers (m = [+ or -] 1, 2, 3, and 4) propagating along untwisted and twisted magnetic flux tubes (representing three observed X-ray solar jets) can become unstable against the KH instability. We have used two models of a soft X-ray jet, namely, considering the jet and its surrounding coronal medium as compressible magnetized plasmas and another simplified model representing the jet as an incompressible medium and its environment as cool coronal plasma. A comparison of the dispersion curves and normalized wave growth rates of the kink (m = 1) mode calculated on the base of these two models shows that one obtains similar dispersion curves' and growth rates' patterns (see Figures 2 and 3). This circumstance justifies the usage of the second, simplified, model in studying the propagation characteristics of MHD modes in moving twisted magnetic flux tubes. We must, however, immediately underline that the onset of KH instability in incompressible magnetic flux tube always requires slightly greater threshold Alfveen Mach numbers and correspondingly higher critical jet speeds.

Our numerical computations show the following:

(i) The flow does not change the nature of propagating stable MHD modes being pure surface waves in a rest magnetic flux tube, but unstable ones due to high enough jet speed become partly surface and partly leaky waves. That is why it is more appropriate to term them generalized surface waves.

(ii) The KH instability of the kink mode (m = [+ or -] 1) starts at Alfven Mach numbers very close to the predicted ones. Numerically found threshold MAs for the higher MHD modes are, in fact, greater than their initial evaluations. Accessible flow velocities for instability onset in jet #11 have been derived for moving untwisted magnetic flux tube of compressible plasma (429 km [s.sup.-1]) and for the m = 4 MHD mode propagating in weakly twisted magnetic flux tube ([epsilon] = 0.025) of incompressible plasma surrounded by cool coronal medium (435.6 km [s.sup.-1]). A nonnegligible decrease of the critical jet speed for all considered modes running in twisted moving flux tube is obtained when the magnetic field twist parameter is e = 0.4--then the critical velocities of the m = [+ or -] 1, 2, 3, and 4 modes are equal to 416, [congruent to] 408, 418, [congruent to] 422, and 424km [s.sup.-1], respectively. The instability onset of all considered MHD modes in jet #8 occurs at flow velocities, in general, lower than the jet speed of 385 km [s.sup.-1]--their magnitudes for the kink (m = 1) mode running on untwisted or slightly twisted magnetic flux tube (with magnetic field twist parameter [epsilon] = 0.025) are equal to 348.3 and 360.6km [s.sup.-1], respectively. The corresponding critical flow velocities for the higher (m = 2, 3, and 4) MHD modes turn out to lie between aforementioned speeds, namely, being equal to [congruent to] 359, 35,7 and 355.5km [s.sup.-1]. Jet #16 possessing the highest flow speed of 532km [s.sup.-1] in Shimojo and Shibata's Events List [10] is, however, stable against the KH instability--the critical flow velocity at which the kink (m = 1) mode would become unstable propagating in untwisted magnetic flux tube is equal to [congruent to] 826 km [s.sup.-1]--much higher than the jet speed. The critical velocities in twisted tubes are even larger, equal to 852.4, 848.3, 844.8, and 842.6km [s.sup.-1] for the m = 1-4 MHD modes, respectively.

(iii) The instability onset for a given mode m is very sensitive to the main input parameters: density contrast, [eta], magnetic field ratio, b, the magnetic field twist [epsilon], and especially the background magnetic field Be. If we assume that, for jet #11 [B.sub.e] = 10 G, Alfven speeds in the two media have new values, notably [v.sub.Ai] = 320.38km [s.sup.-1] and [v.sub.Ae] = 427.5km [s.sup.-1]. With these new Alfveen speeds the input parameters for solving wave dispersion relations (9) and (25) have changed and more specifically [[??].sub.i] = 0.7376, [[??].sub.e] = 0.1506, and b = 1.2636. Now the critical jet velocity for the instability onset dramatically increases--it becomes 751 km [s.sup.-1] for the kink mode in untwisted tube and 774 km [s.sup.-1] in weakly twisted ([epsilon] = 0.025) flux tube.

(iv) If we accept the Paraschivs et al. [90] suggestion that plasma density in the region crossed by coronal jets has to be the sum [n.sub.cor] + [n.sub.jet] or, on average, two times higher than the plasma in the surrounding corona [n.sub.cor]; then the total pressure balance equation, again for jet #11, can be satisfied at [B.sub.e] = 10 G and we have the following set of input parameters: [eta] = 0.437, [v.sub.Ai] = 105.7km [s.sup.-1], [v.sub.Ae] = 427.5km [s.sup.-1], [[??].sub.i] = 6.7782, [[??].sub.e] = 0.1506, and b = 3.781. With these parameters the critical flow velocities for KH instability occurrence reduce to 560 km [s.sup.-1] for the kink mode in untwisted flux tube and to 575km [s.sup.-1] for the same mode in weakly twisted ([epsilon] = 0.025) flux tube. Note that within this approach one obtains critical jet speeds that are approximately with 200 km [s.sup.-1] less than those calculated for the same background magnetic field of 10 G. Here one raises the big question: "which approach, the standard one or that of Paraschiv's et al. [90], will be appropriate for studying/modeling the KH instability in solar X-ray (and other) jets?" The perfect answer to this question will be obtained when we will try to model a really observationally detected KH instability in soft X-ray jets like that of Foullon et al. [57].

It is intriguing to see what critical jet speeds we will obtain for a numerically modeled solar X-jet--a case in point is the study of Miyagoshi and Yokoyama [91]. These authors have presented MHD numerical simulations of solar X-ray jets based on magnetic reconnection model that includes chromospheric evaporation. Peculiar to their study is that total pressure balance equation excludes the magnetic pressure inside the jet saying that it is very weak; that is (in their notation),

[P.sub.jet] = [P.sub.cor] + [B.sub.2]/2[mu] (32)

where B [equivalent to] [B.sub.e] is the background (coronal) magnetic field. Supposing that [n.sub.cor] = 109 [cm.sup.-3], [T.sub.cor] = 106 K, B = 10 G, and [n.sub.jet] = 4.5 x 109 [cm.sup.-3], they obtain from their total balance equation the temperature of the jet, [T.sub.jet] = 6.7 x [10.sup.6] K--indeed a reasonable value. However, to satisfy the standard total pressure balance (1) (assuming a very small, but finite value for the magnetic field in the jet, [B.sub.i]), with the aforementioned values for both the temperature and jet density, we were forced to double the coronal density, that is, to take [n.sub.cor] = [n.sub.e] = 2.0 x [10.sup.9] [cm.sup.-3].In such a case the sound and Alfven speeds in the jet and its environment are csi s 304km [s.sup.-1], [c.sub.se] = 117km [s.sup.-1], [v.sub.Ai] = 47.67km [s.sup.-1], and [v.sub.Ae] = 487.5 km [s.sup.-1]. The ordering of sound and Alfven speeds is the same as for the #11 X-ray jet in the Shimojo and Shibata paper [10] that implies the propagation of pure surface stable modes. The magnetic field inside the jet is [B.sub.i] [congruent to] 1.47 G--hence the magnetic fields ratio is b = 6.817. Thus the input parameters for solving (9) and (25) are [eta] = 0.444, [[??].sub.i] = 40.5872, and [[??].sub.e] = 0.0695 and, as we already found, b = 6.817. Note that this relatively high value of b would suspect a rather big value of the threshold Alfveen Mach number--the instability criterion yields [M.sub.A] > 12.42. Numerically found threshold [M.sub.A] for initiating KH instability of the kink (m = 1) mode in moving untwisted magnetic flux tube of compressible plasma is 12.435 that yields a critical velocity of 607km [s.sup.-1]. For the case of weakly twisted flux tube ([epsilon] = 0.025) in the approximation incompressible jet and cool plasma environment we obtained for the same kink mode a critical velocity of 620 km [s.sup.-1]. The critical speeds for the higher (m = 2-4) modes are accordingly of 615, 609, and 606km [s.sup.-1].

The three examples of coronal X-jets embedded in a background magnetic field of 10 G show that KH instability of basically the kink (m = 1) mode can occur if jets' speeds lie in the range of 500-700 km [s.sup.-1]. Such velocities are generally accessible for high-speed jets; lower critical speeds between 400 and 500km [s.sup.-1] can be observed/detected at moderate external magnetic fields in the range of 6-7 G. This requirement is in agreement with Pucci et al. [34] evaluation that the magnetic field inside a standard X-ray jet would be equal to 2.8 G, while for a blowout jet that value should be around 4.5 G. To be honest, we can say that the most of soft X-ray solar jets with speeds below 400-450 km [s.sup.-1] are stable against the KH instability--this instability can develop only in relatively dense high-speed jets. Note, however, that each soft X-ray jet is a unique event--it requires a separate careful exploration of the instability conditions. Arising KH instability in turn might trigger wave turbulence considered as an effective mechanism for coronal heating.

Our simplified model of studying the possibilities for emerging of KH instability in fast soft X-ray jets has to be improved by assuming a radial density gradient or some magnetic field and flow velocity shears. A challenge also is the modeling of the same instability in rotating and oscillatory swaying coronal hole X-ray jets.

http://dx.doi.org/10.1155/2017/2626495

Competing Interests

The authors declare that they have no competing interests.

Acknowledgments

This work was supported by the Bulgarian Science Fund and the Department of Science & Technology, Government of India Fund, under Indo-Bulgarian Bilateral Project DNTS/INDIA 01/7, /Int/Bulgaria/P-2/12. The authors are indebted to Snezhana Yordanova for drawing one figure.

References

[1] J. M. Beckers, "Solar spicules," Annual Review of Astronomy and Astrophysics, vol. 10, no. 1, pp. 73-100, 1972.

[2] J. R. Roy, "The magnetic properties of solar surges," Solar Physics, vol. 28, no. 1, pp. 95-114,1973.

[3] K. Shibata, T. Nakamura, T. Matsumoto et al., "Chromospheric anemone jets as evidence of ubiquitous reconnection," Science, vol. 318, no. 5856, pp. 1591-1594, 2007.

[4] C. Chifor, H. Isobe, E. Mason et al., "Magnetic flux cancellation associated with a recurring solar jet observed with Hinode, RHESSI, and STEREO/EUVI," Astronomy & Astrophysics, vol. 491, no. 1, pp. 279-288, 2008.

[5] N. Nishizuka, M. Shimizu, T. Nakamura et al., "Giant chromospheric anemone jet observed with Hinode and comparison with magnetohydrodynamic simulations: Evidence of propagating Alfven waves and magnetic reconnection," The Astrophysical Journal Letters, vol. 683, no. 1, pp. L83-L86, 2008.

[6] D. Alexander and L. Fletcher, "High-resolution observations of plasma jets in the solar corona," Solar Physics, vol. 190, no. 1-2, pp. 167-184, 1999.

[7] K. Shibata, Y. Ishido, L. Acton et al., "Observations of X-ray jets with the Yohkoh soft X-ray telescope," Publications of the Astronomical Society of Japan, vol. 44, no. 5, pp. L173-L179, 1992.

[8] S. Tsuneta, L. Acton, M. Bruner et al., "The soft X-ray telescope for the SOLAR-A mission," Solar Physics, vol. 136, no. 1, pp. 37-67, 1991.

[9] M. Shimojo, S. Hashimoto, K. Shibata, T. Hirayama, H. S. Hudson, and L. W. Acton, "Statistical study of solar X-ray jets observed with the Yohkoh soft X-ray telescope," Publications of the Astronomical Society of Japan, vol. 48, no. 1, pp. 123-136, 1996.

[10] M. Shimojo and K. Shibata, "Physical parameters of solar X-ray jets," The Astrophysical Journal, vol. 542, no. 2, pp. 1100-1108, 2000.

[11] L. Culhane, L. K. Harra, D. Baker et al., "Hinode EUV study of jets in the Sun's south polar corona," Publications of the Astronomical Society of Japan, vol. 59, supplement 3, pp. S751-S756, 2007

[12] C. Chifor, P. R. Young, H. Isobe et al., "An active region jet observed with Hinode," Astronomy & Astrophysics, vol. 481, no. 1, pp. L57-L60, 2008.

[13] A. Savcheva, J. Cirtain, E. E. DeLuca et al., "A study of polar jet parameters based on Hinode XRT observations," Publications of the Astronomical Society of Japan, vol. 59, supplement 3, pp. S771-S778, 2007.

[14] J. W. Cirtain, L. Golub, L. Lundquist et al., "Evidence for Alfven waves in solar X-ray jets," Science, vol. 318, no. 5856, pp. 1580-1582, 2007.

[15] Y.-H. Kim, Y.-J. Moon, Y.-D. Parket al., "Small-scaleX-ray/EUV jets seen in Hinode XRT and TRACE," Publications of the Astronomical Society of Japan, vol. 59, supplement 3, pp. S763-S769, 2007.

[16] L.-H. Yang, Y.-C. Jiang, J.-Y. Yang, Y. Bi, R.-S. Zheng, and J.-C. Hong, "Observations ofEUV and soft X-ray recurring jets in an active region," Research in Astronomy and Astrophysics, vol. 11, no. 10, pp. 1229-1242, 2011.

[17] M. S. Madjarska, "Dynamics and plasma properties of an X-ray jet from SUMER, EIS, XRT, and EUVI A & B simultaneous observations," Astronomy & Astrophysics, vol. 526, no. 2, article A19, 2011.

[18] K. Chandrashekhar, A. Bemporad, D. Banerjee, G. R. Gupta, and L. Teriaca, "Characteristics of polar coronal hole jets," Astronomy & Astrophysics, vol. 561, article A104, 2014.

[19] W. D. Pesnell, B. J. Thompson, and P C. Chamberlin, "The Solar Dynamics Observatory (SDO)," Solar Physics, vol. 275, no. 1, pp. 3-15, 2012.

[20] J. R. Lemen, A. M. Title, D. J. Akin et al., "The Atmospheric Imaging Assembly (AIA) on the Solar Dynamics Observatory (SDO)," Solar Physics, vol. 275, no. 1, pp. 17-40, 2012.

[21] P Boerner, C. Edwards, J. Lemen et al., "Initial calibration of the Atmospheric Imaging Assembly (AIA) on the Solar Dynamics Observatory (SDO)," Solar Physics, vol. 275, no. 1, pp. 41-66, 2012.

[22] S. P Moschou, K. Tsinganos, A. Vourlidas, and V. Archontis, "SDO observations of solar jets," Solar Physics, vol. 284, no. 2, pp. 427-438, 2013.

[23] K. Chandrashekhar, R. J. Morton, D. Banerjee, and G. R. Gupta, "The dynamical behaviour of a jet in an on-disk coronal hole observed with AIA/SDO," Astronomy & Astrophysics, vol. 562, article A98, 2014.

[24] A. C. Sterling, R. L. Moore, D. A. Falconer, and M. Adams, "Small-scale filament eruptions as the driver of X-ray jets in solar coronal holes," Nature, vol. 523, no. 7561, pp. 437-440, 2015.

[25] B. Filippov, A. K. Srivastava, B. N. Dwivedi et al., "Formation of a rotating jet during the filament eruption on 2013 April 10-11," Monthly Notices of the Royal Astronomical Society, vol. 451, no. 1, pp. 1117-1129, 2015.

[26] M. Shimojo, K. Shibata, T. Yokoyama, and K. Hori, "One-dimensional and pseudo-two-dimensional hydrodynamic simulations of solar X-ray jets," Astrophysical Journal, vol. 550, no. 2, pp. 1051-1063, 2001.

[27] S. Kamio, H. Hara, T. Watanabe et al., "Velocity structure of jets in a coronal hole," Publications of the Astronomical Society of Japan, vol. 59, supplement 3, pp. S757-S762, 2007

[28] R. L. Moore, J. W. Cirtain, A. C. Sterling, and D. A. Falconer, "Dichotomy of solar coronal jets: Standard jets and blowout jets," The Astrophysical Journal, vol. 720, no. 1, pp. 757-770, 2010.

[29] E. Pariat, S. K. Antiochos, and C. R. DeVore, "A model for solar polar jets," The Astrophysical Journal, vol. 691, no. 1, pp. 61-74, 2009.

[30] B. Filippov, L. Golub, and S. Koutchmy, "X-ray jet dynamics in a polar coronal hole region," Solar Physics, vol. 254, no. 2, pp. 259-269, 2009.

[31] E. Pariat, S. K. Antiochos, and C. R. Devore, "Three-dimensional modeling of quasi-homologous solar jets," Astrophysical Journal, vol. 714, no. 2, pp. 1762-1778, 2010.

[32] F. Jiang, J. Zhang, and S. Yang, "Interaction between an emerging flux region and a pre-existing fan-spine dome observed by IRIS and SDO," Publications of the Astronomical Society of Japan, vol. 67, no. 4, article 78, 2015.

[33] A. C. Sterling, L. K. Harra, and R. L. Moore, "Fibrillar chromospheric spicule-like counterparts to an extreme-ultraviolet and soft X-ray blowout coronal jet," The Astrophysical Journal, vol. 722, no. 2, pp. 1644-1653, 2010.

[34] S. Pucci, G. Poletto, A. C. Sterling, and M. Romoli, "Physical parameters of standard and blowout jets," Astrophysical Journal, vol. 776, article 16, 2013.

[35] R. L. Moore, A. C. Sterling, D. A. Falconer, and D. Robe, "The cool component and the dichotomy, lateral expansion, and axial rotation of solar X-ray jets," The Astrophysical Journal, vol. 769, no. 2, article 134, 2013.

[36] R. L. Moore, A. C. Sterling, and D. A. Falconer, "Magnetic untwisting in solar jets that go into the outer corona in polar coronal holes," Astrophysical Journal, vol. 806, no. 1, article 11, 2015.

[37] G. E. Brueckner, R. A. Howard, M. J. Koomen et al., "The Large Angle Spectroscopic Coronagraph (LASCO)," Solar Physics, vol. 162, no. 1, pp. 357-402, 1995.

[38] B. Schmieder, Y. Guo, F. Moreno-Insertis et al., "Twisting solar coronal jet launched at the boundary of an active region," Astronomy and Astrophysics, vol. 559, article A1, 2013.

[39] E. Pariat, K. Dalmasse, C. R. DeVore, S. K. Antiochos, and J. T. Karpen, "Model for straight and helical solar jets: I. Parametric studies of the magnetic field geometry," Astronomy and Astrophysics, vol. 573, article 130, 2015.

[40] S. Patsourakos, E. Pariat, A. Vourlidas, S. K. Antiochos, and J. P. Wuelser, "STEREO SECCHI stereoscopic observations constraining the initiation of polar coronal jets," The Astrophysical Journal, vol. 680, no. 1, pp. L73-L76, 2008.

[41] C. Gontikakis, V. Archontis, and K. Tsinganos, "Observations and 3D MHD simulations of a solar active region jet," Astronomy & Astrophysics, vol. 506, no. 3, pp. L45-L48, 2009.

[42] R. Chandra, G. R. Gupta, S. Mulay, and D. Tripathi, "Sunspot waves and triggering of homologous active region jets," Monthly Notices of the Royal Astronomical Society, vol. 446, no. 4, pp. 3741-3748, 2015.

[43] J. Chae, H. Wang, C.-Y. Lee, P R. Goode, and U. Schuhle, "Photospheric magnetic field changes associated with transition region explosive events," Astrophysical Journal, vol. 497, no. 2, pp. L109-L112, 1998.

[44] L. R. B. Rubio and C. Beck, "Magnetic flux cancellation in the moat of sunspots: results from simultaneous vector spectropolarimetry in the visible and the infrared," Astrophysical Journal, vol. 626, no. 2, pp. L125-L128, 2005.

[45] L. Heggland, B. De Pontieu, and V. H. Hansteen, "Observational signatures of simulated reconnection events in the solar chromosphere and transition region," The Astrophysical Journal, vol. 702, no. 1, pp. 1-18, 2009.

[46] D. E. Innes, R. H. Cameron, and S. K. Solanki, "EUV jets, type III radio bursts and sunspot waves investigated using SDO/AIA observations," Astronomy and Astrophysics, vol. 531, article L13, 2011.

[47] A. K. Srivastava and K. Murawski, "Observations of a pulse-driven cool polar jet by SDO/AIA," Astronomy & Astrophysics, vol. 534, article A62, 2011.

[48] S. Chandrasekhar, Hydrodynamic and Hydromagnetic Stability, The International Series of Monographs on Physics, Clarendon Press, Oxford, UK, 1961.

[49] T. V Zaqarashvili, I. Zhelyazkov, and L. Ofman, "Stability of rotating magnetized jets in the solar atmosphere. I. Kelvin-Helmholtz instability," Astrophysical Journal, vol. 813, no. 2, article 123, 2015.

[50] D. Ryu, T. W. Jones, and A. Frank, "The magnetohydrodynamic Kelvin-Helmholtz instability: A three-dimensional study of nonlinear evolution," The Astrophysical Journal, vol. 545, no. 1, pp. 475-493, 2000.

[51] S. A. Maslowe, "Shear flow instabilities and transition," in Hydromagnetic Instabilities and the Transition to Turbulence, H. L. Swinney and J. P. Gollub, Eds., pp. 181-228, Springer, Berlin, Germany, 1985.

[52] T. E. Berger, G. Slater, N. Hurlburt et al., "Quiescent prominence dynamics observed with the Hinode solar optical telescope. I. Turbulent upflow plumes," The Astrophysical Journal, vol. 716, no. 2, pp. 1288-1307, 2010.

[53] M. Ryutova, T. Berger, Z. Frank, T. Tarbell, and A. Title, "Observation of plasma instabilities in quiescent prominences," Solar Physics, vol. 267, no. 1, pp. 75-94, 2010.

[54] D. Martinez-Gomez, R. Soler, and J. Terradas, "Onset of the Kelvin-Helmholtz instability in partially ionized magnetic flux tubes," Astronomy and Astrophysics, vol. 578, article A104, 2015.

[55] N. F. Loureiro, A. A. Schekochihin, and D. A. Uzdensky, "Plasmoid and Kelvin-Helmholtz instabilities in Sweet-Parker current sheets," Physical Review E: Statistical, Nonlinear, and Soft Matter Physics, vol. 87, no. 1, article 013102, 2013.

[56] L. Feng, B. Inhester, and W. Q. Gan, "Kelvin-Helmholtz instability of a coronal streamer," Astrophysical Journal, vol. 774, no. 2, article 141, 2013.

[57] C. Foullon, E. Verwichte, V. M. Nakariakov, K. Nykyri, and C. J. Farrugia, "Magnetic Kelvin-Helmholtz instability at the Sun," Astrophysical Journal Letters, vol. 729, no. 1, article L8, 2011.

[58] L. Ofman and B. J. Thompson, "SDO/AIA observation of Kelvin-Helmholtz instability in the solar corona," The Astrophysical Journal Letters, vol. 734, no. 1, article L11, 2011.

[59] C. Foullon, E. Verwichte, K. Nykyri, M. J. Aschwanden, and I. G. Hannah, "Kelvin-Helmholtz instability of the CME reconnection outflow layer in the low corona," The Astrophysical Journal, vol. 767, no. 2, article 170, 2013.

[60] U. V. Mostl, M. Temmer, and A. M. Veronig, "The Kelvin-Helmholtz instability at coronal mass ejection boundaries in the solar corona: Observations and 2.5D MHD simulations," Astrophysical Journal Letters, vol. 766, no. 1, article 12, 2013.

[61] K. Nykyri and C. Foullon, "First magnetic seismology of the CME reconnection outflow layer in the low corona with 2.5-D MHD simulations of the Kelvin-Helmholtz instability," Geophysical Research Letters, vol. 40, no. 16, pp. 4154-4159, 2013.

[62] T. V. Zaqarashvili, A. J. Diaz, R. Oliver, and J. L. Ballester, "Instability of twisted magnetic tubes with axial mass flows," Astronomy & Astrophysics, vol. 516, article A84, 2010.

[63] R. Soler, A. J. Diaz, J. L. Ballester, and M. Goossens, "Kelvin-Helmholtz instability in partially ionized compressible plasmas," The Astrophysical Journal, vol. 749, no. 2, article 163, 2012.

[64] T. V. Zaqarashvili, "Solar spicules: recent challenges in observations and theory," in Proceedings of the 3rd School and Workshop on Space Plasma Physics, I. Zhelyazkov and T. Mishonov, Eds., vol. 1356 of AIP Conference Proceedings, pp. 106-116, AIP Publishing LLC, Melville, NY, USA, 2011.

[65] I. Zhelyazkov, "Magnetohydrodynamic waves and their stability status in solar spicules," Astronomy & Astrophysics, vol. 537, article A124, 2012.

[66] A. Ajabshirizadeh, H. Ebadi, R. E. Vekalati, and K. Molaverdikhani, "The possibility of Kelvin-Helmholtz instability in solar spicules," Astrophysics and Space Science, vol. 357, no. 1, article 33, 2015.

[67] I. Zhelyazkov and T. V. Zaqarashvili, "Kelvin-Helmholtz instability of kink waves in photospheric twisted flux tubes," Astronomy & Astrophysics, vol. 547, article A14, 2012.

[68] I. Zhelyazkov, R. Chandra, A. K. Srivastava, and T. Mishonov, "Kelvin-Helmholtz instability of magnetohydrodynamic waves propagating on solar surges," Astrophysics and Space Science, vol. 356, no. 2, pp. 231-240, 2015.

[69] I. Zhelyazkov, T. V. Zaqarashvili, R. Chandra, A. K. Srivastava, and T. Mishonov, "Kelvin-Helmholtz instability in solar cool surges," Advances in Space Research, vol. 56, no. 12, pp. 2727-2737, 2015.

[70] I. Zhelyazkov, R. Chandra, and A. K. Srivastava, "Can magnetohydrodynamic waves traveling on solar dark mottles become unstable?" Bulgarian Journal of Physics, vol. 42, no. 1, pp. 68-87, 2015.

[71] I. Zhelyazkov, T. V. Zaqarashvili, and R. Chandra, "Kelvin-Helmholtz instability in coronal mass ejecta in the lower corona," Astronomy and Astrophysics, vol. 574, 2015.

[72] I. Zhelyazkov, R. Chandra, and A. K. Srivastava, "Kelvin-Helmholtz instability in coronal mass ejections and solar surges," in Proceedings of the 5th School and Workshop on Space Plasma Physics, I. Zhelyazkov and T. Mishonov, Eds., vol. 1714 of AIP Conference Proceedings, AIP Publishing LLC, Melville, NY, USA, article 030005, 2016.

[73] D. Kuridze, V Henriques, M. Mathioudakis et al., "The dynamics of rapid redshifted and blueshifted excursions in the solar Ha line," The Astrophysical Journal, vol. 802, no. 1, article 26, 2015.

[74] I. Zhelyazkov, "On modeling the Kelvin-Helmholtz instability in solar atmosphere," Journal of Astrophysics and Astronomy, vol. 36, no. 1, pp. 233-254, 2015.

[75] S. Vasheghani Farahani, T. Van Doorsselaere, E. Verwichte, and V. M. Nakariakov, "Propagating transverse waves in soft X-ray coronal jets," Astronomy & Astrophysics, vol. 498, no. 2, pp. L29-L32, 2009.

[76] I. Zhelyazkov, "Review of the magnetohydrodynamic waves and their stability in solar spicules and X-ray jets," in Topics in Magnethydrodynamics, L. Zheng, Ed., chapter 6, InTech, Rieka, Croatia, 2012.

[77] I. Zhelyazkov, "Kelvin-Helmholtz instability of kink waves in photospheric, chromospheric, and X-ray solar jets," in Proceedings of the 4th School and Workshop on Space Plasma Physics, I. Zhelyazkov and T. Mishonov, Eds., vol. 1551, pp. 150-164, AIP Publishing LLC, Melville, NY, USA, 2013.

[78] M. Goossens, J. V. Hollweg, and T. Sakurai, "Resonant behaviour of MHD waves on magnetic flux tubes. III. Effect of equilibrium flow," Solar Physics, vol. 138, no. 2, pp. 233-255, 1992.

[79] K. Hain and R. Lust, "Zur Stabilitat zylindersymmetrischer Plasmakonfigurationen mit Volumenstromen," Zeitschrift fur Naturforschung, vol. 13, pp. 936-940, 1958.

[80] J. P Goedbloed, "Stabilization of magnetohydrodynamic instabilities by force-free magnetic fields: I. Plane plasma layer," Physica, vol. 53, no. 3, pp. 412-444, 1971.

[81] T. Sakurai, M. Goossens, and J. V. Hollweg, "Resonant behaviour of MHD waves on magnetic flux tubes. I. Connection formulae at the resonant surfaces," Solar Physics, vol. 133, no. 2, pp. 227-245, 1991.

[82] P M. Edwin and B. Roberts, "Wave propagation in a magnetic cylinder," Solar Physics, vol. 88, no. 1-2, pp. 179-191, 1983.

[83] M. Terra-Homem, R. Erdelyi, and I. Ballai, "Linear and nonlinear MHD wave propagation in steady-state magnetic cylinders," Solar Physics, vol. 217, no. 2, pp. 199-223, 2003.

[84] V. M. Nakariakov, "MHD oscillations in solar and stellar coronae: current results and perspectives," Advances in Space Research, vol. 39, no. 12, pp. 1804-1813, 2007

[85] K. Bennett, B. Roberts, and U. Narain, "Waves in twisted magnetic flux tubes," Solar Physics, vol. 185, no. 1, pp. 41-59, 1999.

[86] P S. Cally, "Leaky and non-leaky oscillations in magnetic flux tubes," Solar Physics, vol. 103, no. 2, pp. 277-298, 1986.

[87] I. Zhelyazkov, R. Chandra, and A. K. Srivastava, "Kelvin-Helmholtz instability in an active region jet observed with Hinode," Astrophysics and Space Science, vol. 361, no. 2, article 51, 2016.

[88] T. V. Zaqarashvili, Z. Voros, and I. Zhelyazkov, "Kelvin-Helmholtz instability of twisted magnetic flux tubes in the solar wind," Astronomy & Astrophysics, vol. 561, article A62, 2014.

[89] J. W. Dungey and R. E. Loughhead, "Twisted magnetic fields in conducting fluids," Australian Journal of Physics, vol. 7, no. 1, pp. 5-13, 1954.

[90] A. R. Paraschiv, A. Bemporad, and A. C. Sterling, "Physical properties of solar polar jets: A statistical study with Hinode XRT data," Astronomy & Astrophysics, vol. 579, article A96, 2015.

[91] T. Miyagoshi and T. Yokoyama, "Magnetohydrodynamic numerical simulations of solar X-ray jets based on the magnetic reconnection model that includes chromospheric evaporation," The Astrophysical Journal, vol. 593, no. 2, pp. L133-L136, 2003.

Ivan Zhelyazkov, (1) Ramesh Chandra, (2) and Abhishek K. Srivastava (3)

(1) Faculty of Physics, Sofia University, 1164 Sofia, Bulgaria

(2) Department of Physics, Kumaun University, Nainital 263001, India

(3) Department of Physics, Indian Institute of Technology, Banaras Hindu University, Varanasi 221005, India

Correspondence should be addressed to Ivan Zhelyazkov; izh@phys.uni-sofia.bg

Received 8 August 2016; Revised 16 October 2016; Accepted 1 December 2016; Published 8 February 2017

Academic Editor: Valery Nakariakov

Caption: Figure 1: Equilibrium magnetic fields of a soft X-ray solar jet in an untwisted flux tube (a) and in a weakly twisted flux tube (b).

Caption: Figure 2: (a) Dispersion curves of stable and unstable kink (m = 1) MHD mode propagating in a moving untwisted magnetic flux tube of compressible plasma (modeling jet #11) at [eta] = 0.896 and b = 2.44. Unstable dispersion curves located above the middle of the plot have been calculated for four values of the Alfven Mach number [M.sub.A] = 3.865, 3.89, 3.915, and 3.94. (b) The normalized growth rates of the unstable mode for the same values of [M.sub.A]. Red curves in both plots correspond to the onset of KH instability.

Caption: Figure 3: (a) Dispersion curves of stable and unstable kink (m = 1) MHD mode propagating in a moving untwisted magnetic flux tube of incompressible plasma surrounded by cool medium at the same input parameters as in Figure 2. The dispersion curves of unstable modes located above the middle of the plot have been calculated for [M.sub.A] = 3.985, 4.01, 4.035, and 4.06. Alfven Mach numbers associated with the nonlabeled black, green, and purple curves are equal to 3.796,3.815, and 3.83, respectively, while those of far right green and purple semiclosed curves have magnitudes of 3.925 and 3.965. (b) Similar plots as in Figure 2.

Caption: Figure 4: (a) Dispersion curves of stable and unstable kink (m = 1) MHD mode propagating in a moving untwisted magnetic flux tube of compressible plasma (modeling jet #8) at [eta] = 0.963 and b = 4.56. Unstable mode's dispersion curves are pictured in the middle of the plot for [M.sub.A] = 6.6704, 6.68, 6.69, and 6.7. (b) The corresponding normalized growth rates of unstable modes for the same values of [M.sub.A].

Caption: Figure 5: (a) Dispersion curves of stable and unstable kink (m = 1) MHD mode propagating in a moving untwisted magnetic flux tube of compressible plasma (modeling jet #16) at [eta] = 0.95 and b = 1.32. Dispersion curves of unstable modes located in the middle of the plot have been calculated for [M.sub.A] = 2.359305, 2.37, 2.38, and 2.39. (b) Plots of corresponding normalized growth rates.

Caption: Figure 6: (a) Dispersion curves of unstable kink (m = 1) MHD mode propagating in a moving twisted flux tube of incompressible plasma (modeling jet #11) at [eta] = 0.896 and [epsilon] = 0.025 and for [M.sub.A] = 3.989, 4.0 (green curve), 4.011 (purple curve), and 4.022. (b) The normalized growth rates for the same values of [M.sub.A].

Caption: Figure 7: (a) Dispersion curves of unstable (m = 4) MHD mode propagating in a moving twisted flux tube of incompressible plasma at the same parameters as in Figure 6 for [M.sub.A] = 3.926, 3.934 (green curve), 3.942 (purple curve), and 3.95. (b) The normalized growth rates for the same values of [M.sub.A].

Caption: Figure 8: (a) Dispersion curves of unstable (m = 1) MHD mode propagating in a moving twisted flux tube of incompressible plasma (modeling jet #11) at [eta] = 0.896 and e = 0.4 and for [M.sub.A] = 4.04, 4.07 (green curve), 4.1 (purple curve), and 4.13. (b) The normalized growth rates for the same values of [M.sub.A].

Caption: Figure 9: (a) Dispersion curves of unstable (m = -1) MHD mode propagating in a moving twisted flux tube of incompressible plasma (modeling jet #11) at the same parameters as in Figure 8 and for [M.sub.A] = 3.9565, 3.98 (green curve), 4.0 (purple curve), and 4.02. (b) The normalized growth rates for the same values of [M.sub.A].

Caption: Figure 10: (a) Dispersion curves of unstable (m = 2) MHD mode propagating in a moving twisted flux tube of incompressible plasma (modeling jet #11) at the same parameters as in Figure 8 and for [M.sub.A] = 4.2,4.15, 4.1, 4.075, and 4.062. (b) The normalized growth rates for the same values of [M.sub.A].

Caption: Figure 11: (a) Dispersion curves of unstable (m = 3) MHD mode propagating in a moving twisted flux tube of incompressible plasma (modeling jet #11) at the same parameters as in Figure 8 and for [M.sub.A] = 4.093, 4.1, 4.109, and 4.125. (b) The normalized growth rates for the same values of [M.sub.A].

Caption: Figure 12: (a) Dispersion curves of unstable (m = 4) MHD mode propagating in a moving twisted flux tube of incompressible plasma (modeling jet #11) at the same parameters as in Figure 8 and for [M.sub.A] = 4.12,4.13,4.14, and 4.15. (b) The normalized growth rates for the same values of [M.sub.A].

Printer friendly Cite/link Email Feedback | |

Author: | Zhelyazkov, Ivan; Chandra, Ramesh; Srivastava, Abhishek K. |
---|---|

Publication: | Advances in Astronomy |

Article Type: | Report |

Date: | Jan 1, 2017 |

Words: | 14272 |

Previous Article: | Capability of the HAWC Gamma-Ray Observatory for the Indirect Detection of Ultrahigh-Energy Neutrinos. |

Next Article: | Astronomical Tasks for Tests of X-Ray Optics in VZLUSAT-1 Nanosatellite. |

Topics: |