No Arabic abstract
Heliospheric plasmas require multi-scale and multi-physics considerations. On one hand, MHD codes are widely used for global simulations of the solar-terrestrial environments, but do not provide the most elaborate physical description of space plasmas. Hybrid codes, on the other hand, capture important physical processes, such as electric currents and effects of finite Larmor radius, but they can be used locally only, since the limitations in available computational resources do not allow for their use throughout a global computational domain. In the present work, we present a new coupled scheme which allows to switch blocks in the block-adaptive grids from fluid MHD to hybrid simulations, without modifying the self-consistent computation of the electromagnetic fields acting on fluids (in MHD simulation) or charged ion macroparticles (in hybrid simulation). In this way, the hybrid scheme can refine the description in specified regions of interest without compromising the efficiency of the global MHD code.
Both kinetic instabilities and strong turbulence have potential to impact the behavior of space plasmas. To assess effects of these two processes we compare results from a 3 dimensional particle-in-cell (PIC) simulation of collisionless plasma turbulence against observations by the MMS spacecraft in the terrestrial magnetosheath and by the Wind spacecraft in the solar wind. The simulation develops coherent structures and anisotropic ion velocity distributions that can drive micro-instabilities. Temperature-anisotropy driven instability growth rates are compared with inverse nonlinear turbulence time scales. Large growth rates occur near coherent structures; nevertheless linear growth rates are, on average, substantially less than the corresponding nonlinear rates. This result casts some doubt on the usual basis for employing linear instability theory, and raises questions as to why the linear theory appears to work in limiting plasma excursions in anisotropy and plasma beta.
A massively parallel simulation code, called textit{dHybrid}, has been developed to perform global scale studies of space plasma interactions. This code is based on an explicit hybrid model; the numerical stability and parallel scalability of the code are studied. A stabilization method for the explicit algorithm, for regions of near zero density, is proposed. Three-dimensional hybrid simulations of the interaction of the solar wind with unmagnetized artificial objects are presented, with a focus on the expansion of a plasma cloud into the solar wind, which creates a diamagnetic cavity and drives the Interplanetary Magnetic Field out of the expansion region. The dynamics of this system can provide insights into other similar scenarios, such as the interaction of the solar wind with unmagnetized planets.
Using in situ data, accumulated in the turbulent magnetosheath by the Magnetospheric Multiscale (MMS) Mission, we report a statistical study of magnetic field curvature and discuss its role in the turbulent space plasmas. Consistent with previous simulation results, the Probability Distribution Function (PDF) of the curvature is shown to have distinct power-law tails for both high and low value limits. We find that the magnetic-field-line curvature is intermittently distributed in space. High curvature values reside near weak magnetic-field regions, while low curvature values are correlated with small magnitude of the force acting normal to the field lines. A simple statistical treatment provides an explanation for the observed curvature distribution. This novel statistical characterization of magnetic curvature in space plasma provides a starting point for assessing, in a turbulence context, the applicability and impact of particle energization processes, such as curvature drift, that rely on this fundamental quantity.
The Magnetospheric Multiscale (MMS) mission has given us unprecedented access to high cadence particle and field data of magnetic reconnection at Earths magnetopause. MMS first passed very near an X-line on 16 October 2015, the Burch event, and has since observed multiple X-line crossings. Subsequent 3D particle-in-cell (PIC) modeling efforts of and comparison with the Burch event have revealed a host of novel physical insights concerning magnetic reconnection, turbulence induced particle mixing, and secondary instabilities. In this study, we employ the Gkeyll simulation framework to study the Burch event with different classes of extended, multi-fluid magnetohydrodynamics (MHD), including models that incorporate important kinetic effects, such as the electron pressure tensor, with physics-based closure relations designed to capture linear Landau damping. Such fluid modeling approaches are able to capture different levels of kinetic physics in global simulations and are generally less costly than fully kinetic PIC. We focus on the additional physics one can capture with increasing levels of fluid closure refinement via comparison with MMS data and existing PIC simulations.
An innovative field-particle correlation technique is proposed that uses single-point measurements of the electromagnetic fields and particle velocity distribution functions to investigate the net transfer of energy from fields to particles associated with the collisionless damping of turbulent fluctuations in weakly collisional plasmas, such as the solar wind. In addition to providing a direct estimate of the local rate of energy transfer between fields and particles, it provides vital new information about the distribution of that energy transfer in velocity space. This velocity-space signature can potentially be used to identify the dominant collisionless mechanism responsible for the damping of turbulent fluctuations in the solar wind. The application of this novel field-particle correlation technique is illustrated using the simplified case of the Landau damping of Langmuir waves in an electrostatic 1D-1V Vlasov-Poisson plasma, showing that the procedure both estimates the local rate of energy transfer from the electrostatic field to the electrons and indicates the resonant nature of this interaction. Modifications of the technique to enable single-point spacecraft measurements of fields and particles to diagnose the collisionless damping of turbulent fluctuations in the solar wind are discussed, yielding a method with the potential to transform our ability to maximize the scientific return from current and upcoming spacecraft missions, such as the Magnetospheric Multiscale (MMS) and Solar Probe Plus missions.