More Like this

1. Extending the Spectral Difference Method with Divergence Cleaning (SDDC) to the Hall MHD Equations

   https://doi.org/10.22191/nejcs/vol5/iss1/1  

   Hankey, Russell J.; Chen, Kuangxu; Liang, Chunlei (June 2023, Northeast Journal of Complex Systems)

   The Hall Magnetohydrodynamic (MHD) equations are an extension of the standard MHD equations that include the “Hall” term from the general Ohm’s law. The Hall term decouples ion and electron motion physically on the ion inertial length scales. Implementing the Hall MHD equations in a numerical solver allows more physical simulations for plasma dynamics on length scales less than the ion inertial scale length but greater than the electron inertial length. The present effort is an important step towards producing physically correct results to important problems, such as the Geospace Environmental Modeling (GEM) Magnetic Reconnection problem. The solver that is being modified is currently capable of solving the resistive MHD equations on unstructured grids using the spectral difference scheme which is an arbitrarily high-order method that is relatively simple to parallelize. The GEM Magnetic Reconnection problem is used to evaluate whether the Hall MHD equations have been correctly implemented in the solver using the spectral difference method with divergence cleaning (SDDC) algorithm by comparing against the reconnection rates reported in the literature. 

   more » « less
2. Extending the Spectral Difference Method with Divergence Cleaning (SDDC) to the Hall MHD Equations

   Russell J. Hankey, Kuangxu Chen (August 2023, Northeast journal of complex systems)

   The Hall Magnetohydrodynamic (MHD) equations are an extension of the standard MHD equations that include the “Hall” term from the general Ohm’s law. The Hall term decouples ion and electron motion physically on the ion inertial length scales. Implementing the Hall MHD equations in a numerical solver allows more physical simulations for plasma dynamics on length scales less than the ion inertial scale length but greater than the electron inertial length. The present effort is an important step towards producing physically correct results to important problems, such as the Geospace Environmental Modeling (GEM) Magnetic Reconnection problem. The solver that is being modified is currently capable of solving the resistive MHD equations on unstructured grids using the spectral difference scheme which is an arbitrarily high-order method that is relatively simple to parallelize. The GEM Magnetic Reconnection problem is used to evaluate whether the Hall MHD equations have been correctly implemented in the solver using the spectral difference method with divergence cleaning (SDDC) algorithm by comparing against the reconnection rates reported in the literature. 

   more » « less
3. Transition from ion-coupled to electron-only reconnection: Basic physics and implications for plasma turbulence

   Prayash Sharma Pyakurel, Michael Shay (January 2019, Physics of plasmas)

   Using 2.5 dimensional kinetic particle-in-cell (PIC) simulations, we simulate reconnection conditions appropriate for the magnetosheath and solar wind, i.e., plasma beta (ratio of gas pressure to magnetic pressure) greater than 1 and low magnetic shear (strong guide field). Changing the simulation domain size, we find that the ion response varies greatly. For reconnecting regions with scales comparable to the ion inertial length, the ions do not respond to the reconnection dynamics leading to “electron-only” reconnection with very large quasi-steady reconnection rates. Note that in these simulations the ion Larmor radius is comparable to the ion inertial length. The transition to more traditional “ion-coupled” reconnection is gradual as the reconnection domain size increases, with the ions becoming frozen-in in the exhaust when the magnetic island width in the normal direction reaches many ion inertial lengths. During this transition, the quasi-steady reconnection rate decreases until the ions are fully coupled, ultimately reaching an asymptotic value. The scaling of the ion outflow velocity with exhaust width during this electron-only to ion-coupled transition is found to be consistent with a theoretical model of a newly reconnected field line. In order to have a fully frozen-in ion exhaust with ion flows comparable to the reconnection Alfven speed, an exhaust width of at least several ion inertial lengths is needed. In turbulent systems with reconnection occurring between magnetic bubbles associated with fluctuations, using geometric arguments we estimate that fully ion-coupled reconnection requires magnetic bubble length scales of at least several tens of ion inertial lengths. 

   more » « less
4. Electron-scale current sheets and energy dissipation in 3D kinetic-scale plasma turbulence with low electron beta

   https://doi.org/10.1093/mnras/stad1931  

   Vega, Cristian; Roytershteyn, Vadim; Delzanno, Gian_Luca; Boldyrev, Stanislav (June 2023, Monthly Notices of the Royal Astronomical Society)

   ABSTRACT Three-dimensional kinetic-scale turbulence is studied numerically in the regime where electrons are strongly magnetized (the ratio of plasma species pressure to magnetic pressure is βe = 0.1 for electrons and βi = 1 for ions). Such a regime is relevant in the vicinity of the solar corona, the Earth’s magnetosheath, and other astrophysical systems. The simulations, performed using the fluid-kinetic spectral plasma solver (sps) code, demonstrate that the turbulent cascade in such regimes can reach scales smaller than the electron inertial scale, and results in the formation of electron-scale current sheets (ESCS). Statistical analysis of the geometrical properties of the detected ESCS is performed using an algorithm based on the medial axis transform. A typical half-thickness of the current sheets is found to be on the order of electron inertial length or below, while their half-length falls between the electron and ion inertial length. The pressure–strain interaction, used as a measure of energy dissipation, exhibits high intermittency, with the majority of the total energy exchange occurring in current structures occupying approximately 20 per cent of the total volume. Some of the current sheets corresponding to the largest pressure–strain interaction are found to be associated with Alfvénic electron jets and magnetic configurations typical of reconnection. These reconnection candidates represent about 1 per cent of all the current sheets identified. 

   more » « less
5. Scale filtering analysis of kinetic reconnection and its associated turbulence

   https://doi.org/10.1063/5.0185132  

   Adhikari, Subash; Yang, Yan; Matthaeus, William H; Cassak, Paul A; Parashar, Tulasi N; Shay, Michael A (February 2024, Physics of Plasmas)

   Previously, using an incompressible von Kármán–Howarth formalism, the behavior of cross-scale energy transfer in magnetic reconnection and turbulence was found to be essentially identical to each other, independent of an external magnetic (guide) field, in the inertial and energy-containing ranges [Adhikari et al., Phys. Plasmas 30, 082904 (2023)]. However, this description did not account for the energy transfer in the dissipation range for kinetic plasmas. In this Letter, we adopt a scale-filtering approach to investigate this previously unaccounted-for energy transfer channel in reconnection. Using kinetic particle-in-cell simulations of antiparallel and component reconnection, we show that the pressure–strain interaction becomes important at scales smaller than the ion inertial length, where the nonlinear energy transfer term drops off. Also, the presence of a guide field makes a significant difference in the morphology of the scale-filtered energy transfer. These results are consistent with kinetic turbulence simulations, suggesting that the pressure strain interaction is the dominant energy transfer channel between electron scales and ion scales. 

   more » « less