Vyacheslav Olshevsky

Intermittent energy dissipation by turbulent reconnection

Magnetic reconnection—the process responsible for many explosive phenomena in both nature and laboratory—is efficient at dissipating magnetic energy into particle energy. To date, exactly how this dissipation happens remains unclear, owing to the scarcity of multipoint measurements of the “diffusion region” at the sub-ion scale. Here we report such a measurement by Cluster—four spacecraft with separation of 1/5 ion scale. We discover numerous current filaments and magnetic nulls inside the diffusion region of magnetic reconnection, with the strongest currents appearing at spiral nulls (O-lines) and the separatrices. Inside each current filament, kinetic-scale turbulence is significantly increased and the energy dissipation, E′ ⋅ j, is 100 times larger than the typical value. At the jet reversal point, where radial nulls (X-lines) are detected, the current, turbulence, and energy dissipations are surprisingly small. All these features clearly demonstrate that energy dissipation in magnetic reconnection occurs at O-lines but not X-lines.

How to find magnetic nulls and reconstruct field topology with MMS data

In this study, we apply a new method—the first-order Taylor expansion (FOTE)—to find magnetic nulls and reconstruct magnetic field topology, in order to use it with the data from the forthcoming MMS mission. We compare this method with the previously used Poincare index (PI), and find that they are generally consistent, except that the PI method can only find a null inside the spacecraft (SC) tetrahedron, while the FOTE method can find a null both inside and outside the tetrahedron and also deduce its drift velocity. In addition, the FOTE method can (1) avoid limitations of the PI method such as data resolution, instrument uncertainty (Bz offset), and SC separation; (2) identify 3-D null types (A, B, As, and Bs) and determine whether these types can degenerate into 2-D (X and O); (3) reconstruct the magnetic field topology. We quantitatively test the accuracy of FOTE in positioning magnetic nulls and reconstructing field topology by using the data from 3-D kinetic simulations. The influences of SC separation (0.05~1 di) and null-SC distance (0~1 di) on the accuracy are both considered. We find that (1) for an isolated null, the method is accurate when the SC separation is smaller than 1 di, and the null-SC distance is smaller than 0.25~0.5 di; (2) for a null pair, the accuracy is same as in the isolated-null situation, except at the separator line, where the field is nonlinear. We define a parameter ξ ≡ |( λ1 + λ2 + λ3 )|/|λ|max in terms of the eigenvalues (λi) of the null to quantify the quality of our method—the smaller this parameter the better the results. Comparing to the previously used parameter (η≡|∇ ⋅ B|/|∇ × B|), ξ is more relevant for null identification. Using the new method, we reconstruct the magnetic field topology around a radial-type null and a spiral-type null, and find that the topologies are well consistent with those predicted in theory. We therefore suggest using this method to find magnetic nulls and reconstruct field topology with four-point measurements, particularly from Cluster and the forthcoming MMS mission. For the MMS mission, this null-finding algorithm can be used to trigger its burst-mode measurements.

THEORETICAL MODELING OF PROPAGATION OF MAGNETOACOUSTIC WAVES IN MAGNETIC REGIONS BELOW SUNSPOTS

We use two-dimensional numerical simulations and eikonal approximation to study properties of magnetohydrodynamic (MHD) waves traveling below the solar surface through the magnetic structure of sunspots. We consider a series of magnetostatic models of sunspots of different magnetic field strengths, from 10 Mm below the photosphere to the low chromosphere. The purpose of these studies is to quantify the effect of the magnetic field on local helioseismology measurements by modeling waves excited by subphotospheric sources. Time-distance propagation diagrams and wave travel times are calculated for models of various field strengths and compared to the nonmagnetic case. The results clearly indicate that the observed time-distance helioseismology signals in sunspot regions correspond to fast MHD waves. The slow MHD waves form a distinctly different pattern in the time-distance diagram, which has not been detected in observations. The numerical results are in good agreement with the solution in the short-wavelength (eikonal) approximation, providing its validation. The frequency dependence of the travel times is in good qualitative agreement with observations.

Energetics of kinetic reconnection in a three-dimensional null point cluster

We perform three-dimensional particle-in-cell simulations of magnetic reconnection with multiple magnetic null points. Magnetic field energy conversion into kinetic energy is about five times higher than in traditional Harris sheet configuration. More than 85% of initial magnetic field energy is transferred to particle energy during 25 reversed ion cyclofrequencies. Magnetic reconnection in the cluster of null points evolves in three phases. During the first phase, ion beams are excited, then give part of their energy back to the magnetic field in the second phase. In the third phase, magnetic reconnection occurs in many small patches around the current channels formed along the stripes of a low magnetic field. Magnetic reconnection in null points essentially presents three-dimensional features, with no two-dimensional symmetries or current sheets.

Electron and Ion Dynamics of the Solar Wind Interaction with a Weakly Outgassing Comet

Using a 3D fully kinetic approach, we disentangle and explain the ion and electron dynamics of the solar wind interaction with a weakly outgassing comet. We show that, to first order, the dynamical interaction is representative of a four-fluid coupled system. We self-consistently simulate and identify the origin of the warm and suprathermal electron distributions observed by ESAs Rosetta mission to comet 67P/Churyumov-Gerasimenko and conclude that a detailed kinetic treatment of the electron dynamics is critical to fully capture the complex physics of mass-loading plasmas.

MAGNETIC NULL POINTS IN KINETIC SIMULATIONS OF SPACE PLASMAS

We present a systematic attempt to study magnetic null points and the associated magnetic energy conversion in kinetic Particle-in-Cell simulations of various plasma configurations. We address three-dimensional simulations performed with the semi-implicit kinetic electromagnetic code iPic3D in different setups: variations of a Harris current sheet, dipolar and quadrupolar magnetospheres interacting with the solar wind; and a relaxing turbulent configuration with multiple null points. Spiral nulls are more likely created in space plasmas: in all our simulations except lunar magnetic anomaly and quadrupolar mini-magnetosphere the number of spiral nulls prevails over the number of radial nulls by a factor of 3-9. We show that often magnetic nulls do not indicate the regions of intensive energy dissipation. Energy dissipation events caused by topological bifurcations at radial nulls are rather rare and short-lived. The so-called X-lines formed by the radial nulls in the Harris current sheet and lunar magnetic anomaly simulations are rather stable and do not exhibit any energy dissipation. Energy dissipation is more powerful in the vicinity of spiral nulls enclosed by magnetic flux ropes with strong currents at their axes (their cross-sections resemble 2D magnetic islands). These null lines reminiscent of Z-pinches efficiently dissipate magnetic energy due to secondary instabilities such as the two-stream or kinking instability, accompanied by changes in magnetic topology. Current enhancements accompanied by spiral nulls may signal magnetic energy conversion sites in the observational data.

Progress towards physics-based space weather forecasting with exascale computing

Space weather is a rapidly growing field of science which studies processes occurring in the area of space between the Sun and the Earth. The development of space weather forecasting capabilities i ...

Role of Z-pinches in magnetic reconnection in space plasmas

A widely accepted scenario of magnetic reconnection in collisionless space plasmas is the breakage of magnetic field lines in X-points. In laboratory, reconnection is commonly studied in pinches, current channels embedded into twisted magnetic fields. No model of magnetic reconnection in space plasmas considers both null-points and pinches as peers. We have performed a particle-in-cell simulation of magnetic reconnection in a three-dimensional configuration where null-points are present initially, and Z-pinches are formed during the simulation along the lines of spiral null-points. The non-spiral null-points are more stable than spiral ones, and no substantial energy dissipation is associated with them. On the contrary, turbulent magnetic reconnection in the pinches causes the magnetic energy to decay at a rate of ~1.5% per ion gyro period. Dissipation in similar structures is a likely scenario in space plasmas with large fraction of spiral null-points.

A new Particle-in-Cell method for modeling magnetized fluids

Abstract We present a new Particle-in-Cell method for plasma simulations. This is based on the original algorithm of FLIP-MHD, which uses a Lagrangian formulation of the macroscopic equations. A finite-difference approximation of the equations of motion is solved on a fixed (non-moving) grid, while convection of the quantities is modeled with the support of Lagrangian particles. Interpolation with first-order b-splines is used to project the conserved quantities from particles to the grid and back. In this work, we introduce two modifications of the original scheme. A particle volume evolution procedure is adopted to reduce the computational error, based on the formulation used in the Material Point Method for computational mechanics. The additional step introduces little to none computational diffusion and turns out to efficiently suppress the so-called ringing instability, allowing the use of explicit time differencing. Furthermore, we eliminate the need for a Poisson solver in the magnetic field computation with the use of a vector potential in place of the particles’ magnetic moment. The vector potential evolution is modeled with a moving grid and interpolated to the fixed grid points at each time step to obtain a solenoidal magnetic field. The method is tested with a number of standard hydrodynamic and magnetohydrodynamic tests to show the efficiency of the new approach. The results show good agreement with the reference solutions and rather fast time and space convergence.

Slurm: Fluid particle-in-cell code for plasma modeling

Abstract With the approach of exascale computing era, particle-based models are becoming the focus of research due to their excellent scalability. We present a new code, Slurm, which implements the classic particle-in-cell algorithm for modeling magnetized fluids and plasmas. It features particle volume evolution which damps the numerical finite grid instability, and allows modeling of key physical instabilities such as Kelvin–Helmholtz and Rayleigh–Taylor. The magnetic field in Slurm is handled via the electromagnetic vector potential carried by particles. Numerical diffusion of the magnetic flux is extremely low, and the solenoidality of the magnetic field is preserved to machine precision. A double-linked list is used to carry particles, thus implementation of open boundary conditions is simple and efficient. The code is written in C++ with OpenMP multi-threading, and has no external dependencies except for Boost. It is easy to install and use on multi-core desktop computers as well as on large shared-memory machines. Slurm is an ideal tool for its primary goal, modeling of space weather events in the heliosphere. This article walks the reader through the physical model, the algorithm, and all important details of implementation. Ideally, after finishing this paper, the reader should be able to either use Slurm for solving the desired problem, or create a new fluid PIC code.