Fan Guo

Magnetic reconnection is a ubiquitous plasma process in which oppositely directed magnetic field lines break and rejoin, resulting in a change of the magnetic field topology. Reconnection generates magnetic islands: regions enclosed by magnetic field lines and separated by reconnection points. Proper identification of these features is important to understand particle acceleration and overall behavior of plasma. We present a contour-tree based visualization for robust and objective identification of islands and reconnection points in two-dimensional (2D) magnetic reconnection simulations. The application of this visualization to a simple simulation has revealed a physical phenomenon previously not reported, resulting in a more comprehensive understanding of magnetic reconnection.

A rotating pulsar creates a surrounding pulsar wind nebula (PWN) by steadily releasing an energetic wind into the interior of the expanding shockwave of supernova remnant or interstellar medium. At the termination shock of a PWN, the Poynting-flux-dominated relativistic striped wind is compressed. Magnetic reconnection is driven by the compression and converts magnetic energy into particle kinetic energy and accelerating particles to high energies. We carrying out particle-in-cell (PIC) simulations to study the shock structure as well as the energy conversion and particle acceleration mechanism. By analyzing particle trajectories, we find that many particles are accelerated by Fermi-type mechanism. The maximum energy for electrons and positrons can reach hundreds of TeV.

A prime question for plasma physicists is how a fraction of charged particles is accelerated to very high energy. To answer this question, physicists simulate trillions of particles with detailed dynamics and analyze their trajectories. This process requires a range of data analysis tasks with high diversity. In this paper, we present a use case of formulating various analysis tasks on terabyte-scale particle data with a novel data analysis framework called ArrayUDF. The flexibility of ArrayUDF allows it to compose a wide range of particle data operations. We also present optimization strategies to avoid frequent global reduction and to take full advantage of the data locality. Tests show that our optimization methods could accelerate these particle data analysis operations by up to 1,600 times.

Analysis of large-scale simulation output is a core element of scientific inquiry, but analysis queries may experience significant I/O overhead when the data is not structured for efficient retrieval. While in-situ processing allows for improved time-to-insight for many applications, scaling in-situ frameworks to hundreds of thousands of cores can be difficult in practice. The DeltaFS in-situ indexing is a new approach for in-situ processing of massive amounts of data to achieve efficient point and small-range queries. This paper describes the challenges and lessons learned when scaling this in-situ processing function to hundreds of thousands of cores. We propose techniques for scalable all-to-all communication that is memory and bandwidth efficient, concurrent indexing, and specialized LSM-Tree formats. Combining these techniques allows DeltaFS to control the cost of in-situ processing while maintaining 3 orders of magnitude query speedup when scaling alongside the popular VPIC particle-in-cell code to 131,072 cores.

Magnetic flux ropes (or magnetic islands) are ubiquitous space plasma structures. Recent observations suggest that they are often associated with the acceleration of charged particles, but detailed acceleration mechanisms remain unclear. In this study, we present PIC simulations studying particle acceleration due to magnetic flux ropes. We consider a simple 2D configuration of two-magnetic-island coalescence. Some electrons and protons are found to be accelerated to more than 10 times their initial kinetic energies at the end of the simulation. We use a particle tracing technique on the high-energy particles to clarify the associated acceleration mechanisms. We find that reconnection electric field and Fermi-type acceleration due to magnetic island contraction can explain the particle energy gain, which is consistent with previous simulation studies. Our results also suggest that electrons are more responsive to the island contraction mechanism compared to ions. An effective island contraction rate is derived from the simulation data. Finally we briefly discuss a statistical description of particle acceleration associated with interacting magnetic flux ropes, and how it can be connected to simulations.

In this paper we introduce the Indexed Massive Directory, a new technique for indexing data within DeltaFS. With its design as a scalable, server-less file system for HPC platforms, DeltaFS scales file system metadata performance with application scale. The Indexed Massive Directory is a novel extension to the DeltaFS data plane, enabling in-situ indexing of massive amounts of data written to a single directory simultaneously, and in an arbitrarily large number of files. We achieve this through a memory-efficient indexing mechanism for reordering and indexing data, and a log-structured storage layout to pack small writes into large log objects, all while ensuring compute node resources are used frugally. We demonstrate the efficiency of this indexing mechanism through VPIC, a widely-used simulation code that scales to trillions of particles. With DeltaFS, we modify VPIC to create a file for each particle to receive writes of that particle's output data. Dynamically indexing the directory's underlying storage allows us to achieve a 5000x speedup in single particle trajectory queries, which require reading all data for a single particle. This speedup increases with application scale while the overhead is fixed at 3% of available memory.

We study the effect of large-scale coronal magnetic field on the electron acceleration at a spherical coronal shock using a test-particle method. The coronal field is approximated by an analytical solution with a streamer-like magnetic field featured by partially open magnetic field and a current sheet at the equator atop the closed region. It shows that the closed field plays the role of a trapping agency of shock-accelerated electrons, allowing for repetitive reflection and acceleration, therefore can greatly enhance the shock-electron acceleration efficiency. It is found that, with an ad hoc pitch-angle scattering, electron injected in the open field at the shock flank can be accelerated to high energies as well. In addition, if the shock is faster or stronger, a relatively harder electron energy spectrum and a larger maximum energy can be achieved.