Kelvin-Helmholtz instability and induced magnetic reconnection at the Earth’s magnetopause: 3D simulation based on satellite data

Author(s):  
Matteo Faganello ◽  
Manuela Sisti ◽  
Francesco Califano ◽  
Benoit Lavraud

Abstract A 3D two-fluid simulation, using plasma parameters as measured by MMS on September 8th 2015, shows the nonlinear development of the Kelvin-Helmholtz instability at the Earth’s magnetopause. It shows an extremely rich dynamics, including the development of a complex magnetic topology, vortex merging and secondary instabilities. Vortex induced and mid-latitude magnetic reconnection coexist and produce an asymmetric distribution of magnetic reconnection events. Off-equator reconnection exhibits a predominance of events in the southern hemisphere during the early nonlinear phase, as observed by satellites at the dayside magnetopause. The late nonlinear phase shows the development of vortex pairing for all latitudes while secondary Kelvin-Helmholtz instability develops only in the northern hemisphere leading to an enhancement of the occurrence of off-equator reconnection there. Since vortices move tailward while evolving, this suggests that reconnection events in the northern hemisphere should dominate at the nightside magnetopause.

Fluids ◽  
2019 ◽  
Vol 4 (4) ◽  
pp. 189 ◽  
Author(s):  
Francesco Palermo

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.


2011 ◽  
Vol 29 (6) ◽  
pp. 1169-1178 ◽  
Author(s):  
F. Palermo ◽  
M. Faganello ◽  
F. Califano ◽  
F. Pegoraro

Abstract. The nonlinear behaviour of the Kelvin-Helmholtz instability is investigated with a two-fluid simulation code in both sub-magnetosonic and super-magnetosonic regimes in a two-dimensional configuration chosen so as to represent typical conditions observed at the Earth's magnetopause flanks. It is shown that in super-magnetosonic regimes the plasma density inside the vortices produced by the development of the Kelvin-Helmholtz instability is approximately uniform, making the plasma inside the vortices effectively stable against the onset of secondary instabilities. However, the relative motion of the vortices relative to the plasma flow can cause the formation of shock structures. It is shown that in the region where the shocks are attached to the vortex boundaries the plasma conditions change rapidly and develop large gradients that allow for the onset of secondary instabilities not observed in sub-magnetosonic regimes.


2021 ◽  
Author(s):  
Yann Pfau-Kempf ◽  
Minna Palmroth ◽  
Andreas Johlander ◽  
Lucile Turc ◽  
Markku Alho ◽  
...  

<p>Dayside magnetic reconnection at the magnetopause, which is a major driver of space weather, is studied for the first time in a three-dimensional (3D) realistic setup using the Vlasiator hybrid-Vlasov kinetic model. A noon–midnight meridional plane simulation is extended in the dawn–dusk direction to cover 7 Earth radii. The southward interplanetary magnetic field causes magnetic reconnection to occur at the subsolar magnetopause. Perturbations arising from kinetic instabilities in the magnetosheath appear to modulate the reconnection. Its characteristics are consistent with multiple, bursty, and patchy magnetopause reconnection. It is shown that the kinetic behavior of the plasma, as simulated by the model, has consequences on the applicability of methods such as the four-field junction to identify and analyse magnetic reconnection in 3D kinetic simulations.</p>


2015 ◽  
Vol 33 (5) ◽  
pp. 525-530 ◽  
Author(s):  
T. Ishida ◽  
Y. Ogawa ◽  
A. Kadokura ◽  
K. Hosokawa ◽  
Y. Otsuka

Abstract. Ionospheric blobs are localized plasma density enhancements, which are mainly produced by the transportation process of plasma. To understand the deformation process of a blob, observations of plasma parameters with good spatial–temporal resolution are desirable. Thus, we conducted the European Incoherent Scatter radar observations with high-speed meridional scans (60–80 s) during October and December 2013, and observed the temporal evolution of a blob during a substorm on 4 December 2013. This paper is the first report of direct observations of blob deformation during a substorm. The blob deformation arose from an enhanced plasma flow shear during the substorm expansion phase, and then the blob split into two smaller-scale blobs, whose scale sizes were more than ~100 km in latitude. Our analysis indicates that the Kelvin–Helmholtz instability and dissociative recombination could have deformed the blob structure.


2016 ◽  
Vol 34 (3) ◽  
pp. 357-367 ◽  
Author(s):  
Takuma Nakamura ◽  
Rumi Nakamura ◽  
Hiroshi Haseagwa

Abstract. Spatial dimensions of the detailed structures of the electron diffusion region in anti-parallel magnetic reconnection were analyzed based on two-dimensional fully kinetic particle-in-cell simulations. The electron diffusion region in this study is defined as the region where the positive reconnection electric field is sustained by the electron inertial and non-gyrotropic pressure components. Past kinetic studies demonstrated that the dimensions of the whole electron diffusion region and the inner non-gyrotropic region are scaled by the electron inertial length de and the width of the electron meandering motion, respectively. In this study, we successfully obtained more precise scalings of the dimensions of these two regions than the previous studies by performing simulations with sufficiently small grid spacing (1∕16–1∕8 de) and a sufficient number of particles (800 particles cell−1 on average) under different conditions changing the ion-to-electron mass ratio, the background density and the electron βe (temperature). The obtained scalings are adequately supported by some theories considering spatial variations of field and plasma parameters within the diffusion region. In the reconnection inflow direction, the dimensions of both regions are proportional to de based on the background density. Both dimensions also depend on βe based on the background values, but the dependence in the inner region ( ∼ 0.375th power) is larger than the whole region (0.125th power) reflecting the orbits of meandering and accelerated electrons within the inner region. In the outflow direction, almost only the non-gyrotropic component sustains the positive reconnection electric field. The dimension of this single-scale diffusion region is proportional to the ion-electron hybrid inertial length (dide)1∕2 based on the background density and weakly depends on the background βe with the 0.25th power. These firm scalings allow us to predict observable dimensions in real space which are indeed in reasonable agreement with past in situ spacecraft observations in the Earth's magnetotail and have important implications for future observations with higher resolutions such as the NASA Magnetospheric Multiscale (MMS) mission.


2009 ◽  
Vol 27 (2) ◽  
pp. 895-903 ◽  
Author(s):  
D. G. Sibeck

Abstract. We present an analytical model for the magnetic field perturbations associated with flux transfer events (FTEs) on the dayside magnetopause as a function of the shear between the magnetosheath and magnetospheric magnetic fields and the ratio of their strengths. We assume that the events are produced by component reconnection along subsolar reconnection lines with tilts that depend upon the orientation of the interplanetary magnetic field (IMF), and show that the amplitudes of the perturbations generated during southward IMF greatly exceed those during northward IMF. As a result, even if the distributions of magnetic reconnection burst durations/event dimensions are identical during periods of northward and southward IMF orientation, events occurring for southward IMF orientations must predominate in surveys of dayside events. Two factors may restore the balance between events occurring for northward and southward IMF orientations on the flanks of the magnetosphere. Events generated on the dayside magnetopause during periods of southward IMF move poleward, while those generated during periods of northward IMF slip dawnward or duskward towards the flanks. Due to differing event and magnetospheric magnetic field orientations, events that produce weak signatures on the dayside magnetopause during intervals of northward IMF orientation may produce strong signatures on the flanks.


Author(s):  
S. A. Fuselier ◽  
S. Haaland ◽  
P. Tenfjord ◽  
G. Paschmann ◽  
S. Toledo‐Redondo ◽  
...  

Author(s):  
Rouhollah Fatehi ◽  
Mostafa Safdari Shadloo ◽  
Mehrdad T Manzari

Instability of the interface between two immiscible fluids representing the so-called Kelvin–Helmholtz instability problem is studied using smoothed particle hydrodynamics method. Interfacial tension is included, and the fluids are assumed to be inviscid. The time evolution of interfaces is obtained for two low Richardson numbers [Formula: see text] and [Formula: see text] while Bond number varies between zero and infinity. This study focuses on the effect of Bond and Richardson numbers on secondary instability of a two-dimensional shear layer. A brief theoretical discussion is given concerning the linear early time regime followed by numerical investigation of the growth of secondary waves on the main billow. Results show that for [Formula: see text], at all Bond numbers, secondary instabilities start in the early times after a perturbation is imposed, but they grow only for Bond numbers greater than 1. For [Formula: see text], however, secondary instabilities appear only at Bond numbers greater than 10. Finally, based on numerical simulations and using an energy budget analysis involving interfacial potential energy, a quantitative measure is given for the intensity of secondary instabilities using interfacial surface area.


Sign in / Sign up

Export Citation Format

Share Document