Using fully kinetic 3D simulations, the reconnection dynamics of asymmetric current sheets are examined at the Earth's magnetopause. The plasma parameters are selected to model MMS magnetopause diffusion region crossings with guide fields of 0.1, 0.4, and 1 of the reconnecting magnetosheath field. In each case, strong drift-wave fluctuations are observed in the lower-hybrid frequency range at the steep density gradient across the magnetospheric separatrix. These fluctuations give rise to cross-field electron particle transport. In addition, this turbulent mixing leads to significantly enhanced electron parallel heating in comparison to 2D simulations. We study three different methods of quantifying the anomalous dissipation produced by the drift fluctuations, based on spatial averaging, temporal averaging, and temporal averaging followed by integrating along magnetic field lines. Comparison of the different methods reveals complications in identifying and measuring the anomalous dissipation. Nevertheless, the anomalous dissipation from short wavelength drift fluctuations appears weak for each case, and the reconnection rates observed in 3D are nearly the same as in 2D models. The 3D simulations feature a number of interesting new features that are consistent with recent MMS observations, including cold beams of magnetosheath electrons that penetrate into the hotter magnetospheric inflow, the related observation of decreasing temperature in regions of increasing total density, and an effective turbulent diffusion coefficient that agrees with predictions from quasi-linear theory.

Within the resistive magnetohydrodynamic model, high-Lundquist number reconnection layers are unstable to the plasmoid instability, leading to a turbulent evolution where the reconnection rate can be independent of the underlying resistivity. However, the physical relevance of these results remains questionable for many applications. First, the reconnection electric field is often well above the runaway limit, implying that collisional resistivity is invalid. Furthermore, both theory and simulations suggest that plasmoid formation may rapidly induce a transition to kinetic scales, due to the formation of thin current sheets. Here, this problem is studied for the first time using a first-principles kinetic simulation with a Fokker-Planck collision operator in 3D. The low-β reconnecting current layer thins rapidly due to Joule heating before onset of the oblique plasmoid instability. Linear growth rates for standard (k y = 0) tearing modes agree with semi-collisional boundary layer theory, but the angular spectrum of oblique (|k y | > 0) modes is significantly narrower than predicted. In the non-linear regime, flux-ropes formed by the instability undergo complex interactions as they are advected and rotated by the reconnection outflow jets, leading to a turbulent state with stochastic magnetic field. In a manner similar to previous 2D results, super-Dreicer fields induce a transition to kinetic reconnection in thin current layers that form between flux-ropes. These results may be testable within new laboratory experiments.

No abstract

VPIC is a general purpose particle-in-cell simulation code for modeling plasma phenomena such as magnetic reconnection, fusion, solar weather, and laser-plasma interaction in three dimensions using large numbers of particles. VPIC's capacity in both fidelity and scale makes it particularly well-suited for plasma research on pre-exascale and exascale platforms. In this paper we demonstrate the unique challenges involved in preparing the VPIC code for operation at exascale, outlining important optimizations to make VPIC efficient on accelerators. Specifically, we show the work undertaken in adapting VPIC to exploit the portability-enabling framework Kokkos and highlight the enhancements to VPIC's modeling capabilities to achieve performance at exascale. We assess the achieved performance-portability trade-off through a suite of studies on nine different varieties of modern preexascale hardware. Our performance-portability study includes weakscaling runs on three of the top ten TOP500 supercomputers, as well as a comparison of low-level system performance of hardware from four different vendors.

The nonlinear saturation of crossed-beam energy transfer (CBET) for multispeckled laser beams crossing at arbitrary angles is examined using vector particle-in-cell simulations. CBET is found to saturate on fast (∼10s of picosecond) time scales involving ion trapping and excitation of oblique forward stimulated Raman scattering (FSRS). Ion trapping reduces wave damping and speckle interaction increases wave coherence length, together enhancing energy transfer; ion acoustic wave (IAW) breakup in the direction transverse to the wavenumber increases wave damping and contributes to CBET saturation. The seed beam can become unstable to oblique FSRS, which leads to beam deflection at a large angle and a frequency downshift (by the plasma frequency). FSRS saturates on fast ∼picosecond time scales by electron plasma wave self-focusing, leading to enhanced side-loss hot electrons with energy exceeding 300 keV. This may contribute to fuel preheat but FSRS can be mitigated by the presence of a density gradient. Such growth of FSRS contributes to the saturation of CBET. Scaling simulations show that CBET, as well as FSRS and hot electrons, increases with beam average intensity, beam diameter, and crossing area, but that CBET is limited by the excitation of FSRS and IAW breakups in addition to pump depletion. FSRS deflects the seed beam energy by greater than 40% of the incident beam energy and puts a few percent of the incident beam energy into hot electrons. FSRS limits the efficacy of CBET for symmetry tuning at late stages in the implosion and may account for a large portion of the “missing energy” in implosions that use gas-filled hohlraums.

No abstract

Input/Output (I/O) operations can represent a significant proportion of the run-time of parallel scientific computing applications. Although there have been several advances in file format libraries, file system design and I/O hardware, a growing divergence exists between the performance of parallel file systems and the compute clusters that they support. In this paper, we document the design and application of the RIOT I/O toolkit (RIOT) being developed at the University of Warwick with our industrial partners at the Atomic Weapons Establishment and Sandia National Laboratories. We use the toolkit to assess the performance of three industry-standard I/O benchmarks on three contrasting supercomputers, ranging from a mid-sized commodity cluster to a large-scale proprietary IBM BlueGene/P system. RIOT provides a powerful framework in which to analyse I/O and parallel file system behaviour-we demonstrate, for example, the large file locking overhead of IBM's General Parallel File System, which can consume nearly 30% of the total write time in the FLASH-IO benchmark. Through I/O trace analysis, we also assess the performance of HDF-5 in its default configuration, identifying a bottleneck created by the use of suboptimal Message Passing Interface hints. Furthermore, we investigate the performance gains attributed to the Parallel Log-structured File System (PLFS) being developed by EMC Corporation and the Los Alamos National Laboratory. Our evaluation of PLFS involves two high-performance computing systems with contrasting I/O backplanes and illustrates the varied improvements to I/O that result from the deployment of PLFS (ranging from up to 25× speed-up in I/O performance on a large I/O installation to 2× speed-up on the much smaller installation at the University of Warwick).

Conventional explicit electromagnetic particle-in-cell (PIC) algorithms do not conserve discrete energy exactly. Time-centered fully implicit PIC algorithms can conserve discrete energy exactly, but may introduce large dispersion errors in the light-wave modes. This can lead to intolerable simulation errors where accurate light propagation is needed (e.g. in laser-plasma interactions). In this study, we selectively combine the leap-frog and Crank-Nicolson methods to produce an exactly energy-and charge-conserving relativistic electromagnetic PIC algorithm. Specifically, we employ the leap-frog method for Maxwell's equations, and the Crank-Nicolson method for the particle equations. The semi-implicit algorithm admits exact global energy conservation, exact local charge conservation, and preserves the dispersion properties of the leap-frog method for the light wave. The algorithm employs a new particle pusher designed to maximize efficiency and minimize wall-clock-time impact vs. the explicit alternative. It has been implemented in a code named iVPIC, based on the Los Alamos National Laboratory VPIC code (https://github.com/losalamos/vpic).We present numerical results that demonstrate the properties of the scheme with sample test problems: relativistic two-stream instability, Weibel instability, and laser-plasma instabilities.