Motion of a current-vortex sheet in the magnetic Kelvin-Helmholtz instability

In this paper, we consider the Kelvin-Helmholtz instability in the magnetohydrodynamic flow. The motion of the interface is described by a current-vortex sheet. We examine the linear stability of the current-vortex sheet model and determine the growth rate of the interface. The interface is linearly stable for M_{A}<2 where M_{A} represents the Alfvén Mach number. It is found that the interface is linearly unstable in the limit of the critical Alfvén Mach number M_{A}=2, due to resonance of eigenvalues. We perform numerical simulations for the current-vortex sheet for both regimes of M_{A}<2 and M_{A}>2. The numerical results show the stabilizing effects of the magnetic field on the evolution of the current-vortex sheet when the magnetic field is sufficiently large. For the regime M_{A}<2, the sheet oscillates both longitudinally and transversely and the transverse surface wave is pronounced for a large M_{A}. Remarkably, the interface is nonlinearly unstable for M_{A}≈2, for M_{A}<2, which may be due to the propagation of surface waves. For the regime M_{A}>2, the roll-up of the spiral is weakened and the spiral is more pinched and stretched for smaller M_{A}. A comparison of the unstable evolutions of large and small values of M_{A} shows significant differences of the magnetic field and vortex sheet strength.

Kelvin-Helmholtz instability in a compressible dust fluid flow

We report the first experimental observations of a single-mode Kelvin-Helmholtz instability in a flowing dusty plasma in which the flow is compressible in nature. The experiments are performed in an inverted [Formula: see text]-shaped dusty plasma experimental device in a DC glow discharge Argon plasma environment. A gas pulse valve is installed in the experimental chamber to initiate directional motion to a particular dust layer. The shear generated at the interface of the moving and stationary layers leads to the excitation of the Kelvin-Helmholtz instability giving rise to a vortex structure at the interface. The growth rate of the instability is seen to decrease with an increase in the gas flow velocity in the valve and the concomitant increase in the compressibility of the dust flow. The shear velocity is further increased by making the stationary layer to flow in an opposite direction. The magnitude of the vorticity is seen to become stronger while the vortex becomes smaller with such an increase of the shear velocity. A molecular dynamics simulation provides good theoretical support to the experimental findings.

Modeling Kelvin-Helmholtz Instability at the High-Latitude Boundary Layer in a Global Magnetosphere Simulation

The Kelvin-Helmholtz instability at the magnetospheric boundary plays a crucial role in solar wind-magnetosphere-ionosphere coupling, particle entry, and energization. The full extent of its impact has remained an open question due, in part, to global models without sufficient resolution to capture waves at higher latitudes. Using global magnetohydrodynamic simulations, we investigate an event when the Magnetospheric Multiscale (MMS) mission observed periodic low-frequency waves at the dawn-flank, high-latitude boundary layer. We show the layer to be unstable, even though the slow solar wind with the draped interplanetary magnetic field is seemingly unfavorable for wave generation. The simulated velocity shear at the boundary is thin ( ∼ 0.65 R E ) and requires commensurately high spatial resolution. These results, together with MMS observations, confirm for the first time in fully three-dimensional global geometry that KH waves can grow in this region and thus can be an important process for energetic particle acceleration, dynamics, and transport.

3D simulations of oxygen shell burning with and without magnetic fields

We present a first 3D magnetohydrodynamic (MHD) simulation of convective oxygen and neon shell burning in a non-rotating [Formula: see text] star shortly before core collapse to study the generation of magnetic fields in supernova progenitors. We also run a purely hydrodynamic control simulation to gauge the impact of the magnetic fields on the convective flow and on convective boundary mixing. After about 17 convective turnover times, the magnetic field is approaching saturation levels in the oxygen shell with an average field strength of [Formula: see text], and does not reach kinetic equipartition. The field remains dominated by small-to-medium scales, and the dipole field strength at the base of the oxygen shell is only [Formula: see text]. The angle-averaged diagonal components of the Maxwell stress tensor mirror those of the Reynolds stress tensor, but are about one order of magnitude smaller. The shear flow at the oxygen-neon shell interface creates relatively strong fields parallel to the convective boundary, which noticeably inhibit the turbulent entrainment of neon into the oxygen shell. The reduced ingestion of neon lowers the nuclear energy generation rate in the oxygen shell and thereby slightly slows down the convective flow. Aside from this indirect effect, we find that magnetic fields do not appreciably alter the flow inside the oxygen shell. We discuss the implications of our results for the subsequent core-collapse supernova and stress the need for longer simulations, resolution studies, and an investigation of non-ideal effects for a better understanding of magnetic fields in supernova progenitors.

Decay of Kelvin-Helmholtz Vortices at the Earth's Magnetopause Under Pure Southward IMF Conditions

At the Earth's low-latitude magnetopause, clear signatures of the Kelvin-Helmholtz (KH) waves have been frequently observed during periods of the northward interplanetary magnetic field (IMF), whereas these signatures have been much less frequently observed during the southward IMF. Here, we performed the first 3-D fully kinetic simulation of the magnetopause KH instability under the southward IMF condition. The simulation demonstrates that fast magnetic reconnection is induced at multiple locations along the vortex edge in an early nonlinear growth phase of the instability. The reconnection outflow jets significantly disrupt the flow of the nonlinear KH vortex, while the disrupted turbulent flow strongly bends and twists the reconnected field lines. The resulting coupling of the complex field and flow patterns within the magnetopause boundary layer leads to a quick decay of the vortex structure, which may explain the difference in the observation probability of KH waves between northward and southward IMF conditions.

Energy Transport by Kelvin-Helmholtz Instability at the Magnetopause

By means of the formation of vortices in the nonlinear phase, the Kelvin Helmholtz instability is able to redistribute the flux of energy of the solar wind that flows parallel to the magnetopause. The energy transport associated with the Kelvin Helmholtz instability contributes significantly to the magnetosphere and magnetosheath dynamics, in particular at the flanks of the magnetopause where the presence of a magnetic field perpendicular to the velocity flow does not inhibit the instability development. By means of a 2D two-fluid simulation code, the behavior of the Kelvin Helmholtz instability is investigated in the presence of typical conditions observed at the magnetopause. In particular, the energy penetration in the magnetosphere is studied as a function of an important parameter such as the solar wind velocity. The influence of the density jump at the magnetopause is also discussed.

Strongly localized magnetic reconnection by the super-Alfvénic shear flow

We demonstrate that the dragging of the magnetic field by the super-Alfvénic shear flows out of the reconnection plane can strongly localize the reconnection x-line in collisionless pair plasmas, reversing the current direction at the x-line. Reconnection with this new morphology, which is impossible in resistive-magnetohydrodynamics, is enabled by the particle inertia. Surprisingly, the quasi-steady reconnection rate remains of order 0.1 even though the aspect ratio of the local x-line geometry is larger than unity, which completely excludes the role of tearing physics. We explain this by examining the transport of the reconnected magnetic flux and the opening angle ma de by the upstream magnetic field, concluding that the reconnection rate is still limited by the constraint imposed at the inflow region. Based on these findings, we propose that this often observed fast rate value of order 0.1 itself, in general, is an upper bound value determined by the upstream constraint, independent of the localization mechanism and dissipation therein.

Observing Kelvin-Helmholtz instability in solar blowout jet

Kelvin-Helmholtz instability (KHI) is a basic physical process in fluids and magnetized plasmas, with applications successfully modelling e.g. exponentially growing instabilities observed at magnetospheric and heliospheric boundaries, in the solar or Earth's atmosphere and within astrophysical jets. Here, we report the discovery of the KHI in solar blowout jets and analyse the detailed evolution by employing high-resolution data from the Interface Region Imaging Spectrograph (IRIS) satellite launched in 2013. The particular jet we focus on is rooted in the surrounding penumbra of the main negative polarity sunspot of Active Region 12365, where the main body of the jet is a super-penumbral structure. At its maximum, the jet has a length of 90 Mm, a width of 19.7 Mm, and its density is about 40 times higher than its surroundings. During the evolution of the jet, a cavity appears near the base of the jet, and bi-directional flows originated from the top and bottom of the cavity start to develop, indicating that magnetic reconnection takes place around the cavity. Two upward flows pass along the left boundary of the jet successively. Next, KHI develops due to a strong velocity shear (∼204 km s-1) between these two flows, and subsequently the smooth left boundary exhibits a sawtooth pattern, evidencing the onset of the instability.

Turbulent mass transfer caused by vortex induced reconnection in collisionless magnetospheric plasmas

Magnetic reconnection is believed to be the main driver to transport solar wind into the Earth's magnetosphere when the magnetopause features a large magnetic shear. However, even when the magnetic shear is too small for spontaneous reconnection, the Kelvin-Helmholtz instability driven by a super-Alfvénic velocity shear is expected to facilitate the transport. Although previous kinetic simulations have demonstrated that the non-linear vortex flows from the Kelvin-Helmholtz instability gives rise to vortex-induced reconnection and resulting plasma transport, the system sizes of these simulations were too small to allow the reconnection to evolve much beyond the electron scale as recently observed by the Magnetospheric Multiscale (MMS) spacecraft. Here, based on a large-scale kinetic simulation and its comparison with MMS observations, we show for the first time that ion-scale jets from vortex-induced reconnection rapidly decay through self-generated turbulence, leading to a mass transfer rate nearly one order higher than previous expectations for the Kelvin-Helmholtz instability.

Theory and observations of upward field-aligned currents at the magnetopause boundary layer

The dependence of the upward field-aligned current density (J_‖) at the dayside magnetopause boundary layer is well described by a simple analytic model based on a velocity shear generator. A previous observational survey confirmed that the scaling properties predicted by the analytical model are applicable between 11 and 17 MLT. We utilize the analytic model to predict field-aligned currents using solar wind and ionospheric parameters and compare with direct observations. The calculated and observed parallel currents are in excellent agreement, suggesting that the model may be useful to infer boundary layer structures. However, near noon, where velocity shear is small, the kinetic pressure gradients and thermal currents, which are not included in the model, could make a small but significant contribution to J_‖. Excluding data from noon, our least squares fit returns log(J_{‖,max_cal}) = (0.96 ± 0.04) log(J_{‖_obs}) + (0.03 ± 0.01) where J_{‖,max_cal} = calculated J_‖,max and J_{‖_obs} = observed J_‖.

Kelvin-Helmholtz versus Hall magnetoshear instability in astrophysical flows

We study the stability of shear flows in a fully ionized plasma. Kelvin-Helmholtz is a well-known macroscopic and ideal shear-driven instability. In sufficiently low-density plasmas, also the microscopic Hall magnetoshear instability can take place. We performed three-dimensional simulations of the Hall-magnetohydrodynamic equations where these two instabilities are present, and carried out a comparative study. We find that when the shear flow is so intense that its vorticity surpasses the ion-cyclotron frequency of the plasma, the Hall magnetoshear instability is not only non-negligible, but it actually displays growth rates larger than those of the Kelvin-Helmholtz instability.

Role of the Kelvin-Helmholtz instability in the evolution of magnetized relativistic sheared plasma flows

We explore, via analytical and numerical methods, the Kelvin-Helmholtz (KH) instability in relativistic magnetized plasmas, with applications to astrophysical jets. We solve the single-fluid relativistic magnetohydrodynamic (RMHD) equations in conservative form using a scheme which is fourth order in space and time. To recover the primitive RMHD variables, we use a highly accurate, rapidly convergent algorithm which improves upon such schemes as the Newton-Raphson method. Although the exact RMHD equations are marginally stable, numerical discretization renders them unstable. We include numerical viscosity to restore numerical stability. In relativistic flows, diffusion can lead to a mathematical anomaly associated with frame transformations. However, in our KH studies, we remain in the rest frame of the system, and therefore do not encounter this anomaly. We use a two-dimensional slab geometry with periodic boundary conditions in both directions. The initial unperturbed velocity peaks along the central axis and vanishes asymptotically at the transverse boundaries. Remaining unperturbed quantities are uniform, with a flow-aligned unperturbed magnetic field. The early evolution in the nonlinear regime corresponds to the formation of counter-rotating vortices, connected by filaments, which persist in the absence of a magnetic field. A magnetic field inhibits the vortices through a series of stages, namely, field amplification, vortex disruption, turbulent breakdown, and an approach to a flow-aligned equilibrium configuration. Similar stages have been discussed in MHD literature. We examine how and to what extent these stages manifest in RMHD for a set of representative field strengths. To characterize field strength, we define a relativistic extension of the Alfvénic Mach number M(A). We observe close complementarity between flow and magnetic field behavior. Weaker fields exhibit more vortex rotation, magnetic reconnection, jet broadening, and intermediate turbulence. Sufficiently strong fields (M(A)<6) completely suppress vortex formation. Maximum jet deceleration, and viscous dissipation, occur for intermediate vortex-disruptive fields, while electromagnetic energy is maximized for the strongest fields which allow vortex formation. Highly relativistic flows destabilize the system, supporting modes with near-maximum growth at smaller wavelengths than the shear width of the velocity. This helps to explain early numerical breakdown of highly relativistic simulations using numerical viscosity, a long-standing problem. While magnetic fields generally stabilize the system, we have identified many features of the complex and turbulent reorganization that occur for sufficiently weak fields in RMHD flows, and have described the transition from disruptive to stabilizing fields at M(A)≈6. Our results are qualitatively similar to observations of numerous jets, including M87, whose knots may exhibit vortex-like behavior. Furthermore, in both the linear and nonlinear analyses, we have successfully unified the HD, MHD, RHD, and RMHD regimes.

[Not Available]

The Kelvin-Helmholtz instability gained scientific attention after observations at Venus by the spacecraft Pioneer Venus Orbiter gave rise to speculations that the instability contributes to the loss of planetary ions through the formation of plasma clouds. Since then, a handful of studies were devoted to the Kelvin-Helmholtz instability at the ionopause and its implications for Venus. The aim of this study is to investigate the stability of the two instability-relevant boundary layers around Venus: the induced magnetopause and the ionopause. We solve the 2D magnetohydrodynamic equations with the total variation diminishing Lax-Friedrichs algorithm and perform simulation runs with different initial conditions representing the situation at the boundary layers around Venus. Our results show that the Kelvin-Helmholtz instability does not seem to be able to reach its nonlinear vortex phase at the ionopause due to the very effective stabilizing effect of a large density jump across this boundary layer. This seems also to be true for the induced magnetopause for low solar activity. During high solar activity, however, there could occur conditions at the induced magnetopause which are in favour of the nonlinear evolution of the instability. For this situation, we estimated roughly a growth rate for planetary oxygen ions of about 7.6 × 10(25) s(-1), which should be regarded as an upper limit for loss due to the Kelvin-Helmholtz instability.

Kinetic effects on the Kelvin-Helmholtz instability in ion-to-magnetohydrodynamic scale transverse velocity shear layers: Particle simulations

Ion-to-magnetohydrodynamic scale physics of the transverse velocity shear layer and associated Kelvin-Helmholtz instability (KHI) in a homogeneous, collisionless plasma are investigated by means of full particle simulations. The shear layer is broadened to reach a kinetic equilibrium when its initial thickness is close to the gyrodiameter of ions crossing the layer, namely, of ion-kinetic scale. The broadened thickness is larger in B⋅Ω<0 case than in B⋅Ω>0 case, where Ω is the vorticity at the layer. This is because the convective electric field, which points out of (into) the layer for B⋅Ω<0 (B⋅Ω>0), extends (reduces) the gyrodiameters. Since the kinetic equilibrium is established before the KHI onset, the KHI growth rate depends on the broadened thickness. In the saturation phase of the KHI, the ion vortex flow is strengthened (weakened) for B⋅Ω<0 (B⋅Ω>0), due to ion centrifugal drift along the rotational plasma flow. In ion inertial scale vortices, this drift effect is crucial in altering the ion vortex size. These results indicate that the KHI at Mercury-like ion-scale magnetospheric boundaries could show clear dawn-dusk asymmetries in both its linear and nonlinear growth.

Magnetic effects on the coalescence of Kelvin-Helmholtz vortices

We simulate the coalescence process of MHD-scale Kelvin-Helmholtz vortices with the electron inertial effects taken into account. Reconnection of highly stretched magnetic field lines within a rolled-up vortex destroys the vortex itself and the coalescence process, which is well known in ordinary fluid dynamics, is seen to be inhibited. When the magnetic field is initially antiparallel across the shear layer, on the other hand, multiple vortices are seen to coalesce continuously because another type of magnetic reconnection prevents the vortex decay. This type of reconnection at the hyperbolic point also changes the field line connectivity and thus leads to large-scale plasma mixing across the shear layer.

Competing mechanisms of plasma transport in inhomogeneous configurations with velocity shear: the solar-wind interaction with earth's magnetosphere

Two-dimensional simulations of the Kelvin-Helmholtz instability in an inhomogeneous compressible plasma with a density gradient show that, in a transverse magnetic field configuration, the vortex pairing process and the Rayleigh-Taylor secondary instability compete during the nonlinear evolution of the vortices. Two different regimes exist depending on the value of the density jump across the velocity shear layer. These regimes have different physical signatures that can be crucial for the interpretation of satellite data of the interaction of the solar wind with the magnetospheric plasma.

Transport of solar wind into Earth's magnetosphere through rolled-up Kelvin-Helmholtz vortices

Establishing the mechanisms by which the solar wind enters Earth's magnetosphere is one of the biggest goals of magnetospheric physics, as it forms the basis of space weather phenomena such as magnetic storms and aurorae. It is generally believed that magnetic reconnection is the dominant process, especially during southward solar-wind magnetic field conditions when the solar-wind and geomagnetic fields are antiparallel at the low-latitude magnetopause. But the plasma content in the outer magnetosphere increases during northward solar-wind magnetic field conditions, contrary to expectation if reconnection is dominant. Here we show that during northward solar-wind magnetic field conditions-in the absence of active reconnection at low latitudes-there is a solar-wind transport mechanism associated with the nonlinear phase of the Kelvin-Helmholtz instability. This can supply plasma sources for various space weather phenomena.

Decay of MHD-scale Kelvin-Helmholtz vortices mediated by parasitic electron dynamics

We have simulated nonlinear development of MHD-scale Kelvin-Helmholtz (KH) vortices by a two-dimensional two-fluid system including finite electron inertial effects. In the presence of moderate density jump across a shear layer, in striking contrast to MHD results, MHD KH vortices are found to decay by the time one eddy turnover is completed. The decay is mediated by smaller vortices that appear within the parent vortex and stays effective even when the shear layer width is made larger. It is shown that the smaller vortices are basically of MHD nature while the seeding for these is achieved by the electron inertial effect. Application of the results to the magnetotail boundary layer is discussed.