-
Understanding Large-Scale Plasma Simulation Challenges for Fusion Energy on Supercomputers
Authors:
Jeremy J. Williams,
Ashish Bhole,
Dylan Kierans,
Matthias Hoelzl,
Ihor Holod,
Weikang Tang,
David Tskhakaya,
Stefan Costea,
Leon Kos,
Ales Podolnik,
Jakub Hromadka,
JOREK Team,
Erwin Laure,
Stefano Markidis
Abstract:
Understanding plasma instabilities is essential for achieving sustainable fusion energy, with large-scale plasma simulations playing a crucial role in both the design and development of next-generation fusion energy devices and the modelling of industrial plasmas. To achieve sustainable fusion energy, it is essential to accurately model and predict plasma behavior under extreme conditions, requiri…
▽ More
Understanding plasma instabilities is essential for achieving sustainable fusion energy, with large-scale plasma simulations playing a crucial role in both the design and development of next-generation fusion energy devices and the modelling of industrial plasmas. To achieve sustainable fusion energy, it is essential to accurately model and predict plasma behavior under extreme conditions, requiring sophisticated simulation codes capable of capturing the complex interaction between plasma dynamics, magnetic fields, and material surfaces. In this work, we conduct a comprehensive HPC analysis of two prominent plasma simulation codes, BIT1 and JOREK, to advance understanding of plasma behavior in fusion energy applications. Our focus is on evaluating JOREK's computational efficiency and scalability for simulating non-linear MHD phenomena in tokamak fusion devices. The motivation behind this work stems from the urgent need to advance our understanding of plasma instabilities in magnetically confined fusion devices. Enhancing JOREK's performance on supercomputers improves fusion plasma code predictability, enabling more accurate modelling and faster optimization of fusion designs, thereby contributing to sustainable fusion energy. In prior studies, we analysed BIT1, a massively parallel Particle-in-Cell (PIC) code for studying plasma-material interactions in fusion devices. Our investigations into BIT1's computational requirements and scalability on advanced supercomputing architectures yielded valuable insights. Through detailed profiling and performance analysis, we have identified the primary bottlenecks and implemented optimization strategies, significantly enhancing parallel performance. This previous work serves as a foundation for our present endeavours.
△ Less
Submitted 29 June, 2024;
originally announced July 2024.
-
Understanding the Impact of openPMD on BIT1, a Particle-in-Cell Monte Carlo Code, through Instrumentation, Monitoring, and In-Situ Analysis
Authors:
Jeremy J. Williams,
Stefan Costea,
Allen D. Malony,
David Tskhakaya,
Leon Kos,
Ales Podolnik,
Jakub Hromadka,
Kevin Huck,
Erwin Laure,
Stefano Markidis
Abstract:
Particle-in-Cell Monte Carlo simulations on large-scale systems play a fundamental role in understanding the complexities of plasma dynamics in fusion devices. Efficient handling and analysis of vast datasets are essential for advancing these simulations. Previously, we addressed this challenge by integrating openPMD with BIT1, a Particle-in-Cell Monte Carlo code, streamlining data streaming and s…
▽ More
Particle-in-Cell Monte Carlo simulations on large-scale systems play a fundamental role in understanding the complexities of plasma dynamics in fusion devices. Efficient handling and analysis of vast datasets are essential for advancing these simulations. Previously, we addressed this challenge by integrating openPMD with BIT1, a Particle-in-Cell Monte Carlo code, streamlining data streaming and storage. This integration not only enhanced data management but also improved write throughput and storage efficiency. In this work, we delve deeper into the impact of BIT1 openPMD BP4 instrumentation, monitoring, and in-situ analysis. Utilizing cutting-edge profiling and monitoring tools such as gprof, CrayPat, Cray Apprentice2, IPM, and Darshan, we dissect BIT1's performance post-integration, shedding light on computation, communication, and I/O operations. Fine-grained instrumentation offers insights into BIT1's runtime behavior, while immediate monitoring aids in understanding system dynamics and resource utilization patterns, facilitating proactive performance optimization. Advanced visualization techniques further enrich our understanding, enabling the optimization of BIT1 simulation workflows aimed at controlling plasma-material interfaces with improved data analysis and visualization at every checkpoint without causing any interruption to the simulation.
△ Less
Submitted 27 June, 2024;
originally announced June 2024.
-
Experience and Analysis of Scalable High-Fidelity Computational Fluid Dynamics on Modular Supercomputing Architectures
Authors:
Martin Karp,
Estela Suarez,
Jan H. Meinke,
Måns I. Andersson,
Philipp Schlatter,
Stefano Markidis,
Niclas Jansson
Abstract:
The never-ending computational demand from simulations of turbulence makes computational fluid dynamics (CFD) a prime application use case for current and future exascale systems. High-order finite element methods, such as the spectral element method, have been gaining traction as they offer high performance on both multicore CPUs and modern GPU-based accelerators. In this work, we assess how high…
▽ More
The never-ending computational demand from simulations of turbulence makes computational fluid dynamics (CFD) a prime application use case for current and future exascale systems. High-order finite element methods, such as the spectral element method, have been gaining traction as they offer high performance on both multicore CPUs and modern GPU-based accelerators. In this work, we assess how high-fidelity CFD using the spectral element method can exploit the modular supercomputing architecture at scale through domain partitioning, where the computational domain is split between a Booster module powered by GPUs and a Cluster module with conventional CPU nodes. We investigate several different flow cases and computer systems based on the modular supercomputing architecture (MSA). We observe that for our simulations, the communication overhead and load balancing issues incurred by incorporating different computing architectures are seldom worthwhile, especially when I/O is also considered, but when the simulation at hand requires more than the combined global memory on the GPUs, utilizing additional CPUs to increase the available memory can be fruitful. We support our results with a simple performance model to assess when running across modules might be beneficial. As MSA is becoming more widespread and efforts to increase system utilization are growing more important our results give insight into when and how a monolithic application can utilize and spread out to more than one module and obtain a faster time to solution.
△ Less
Submitted 9 May, 2024;
originally announced May 2024.
-
Optimizing BIT1, a Particle-in-Cell Monte Carlo Code, with OpenMP/OpenACC and GPU Acceleration
Authors:
Jeremy J. Williams,
Felix Liu,
David Tskhakaya,
Stefan Costea,
Ales Podolnik,
Stefano Markidis
Abstract:
On the path toward developing the first fusion energy devices, plasma simulations have become indispensable tools for supporting the design and development of fusion machines. Among these critical simulation tools, BIT1 is an advanced Particle-in-Cell code with Monte Carlo collisions, specifically designed for modeling plasma-material interaction and, in particular, analyzing the power load distri…
▽ More
On the path toward developing the first fusion energy devices, plasma simulations have become indispensable tools for supporting the design and development of fusion machines. Among these critical simulation tools, BIT1 is an advanced Particle-in-Cell code with Monte Carlo collisions, specifically designed for modeling plasma-material interaction and, in particular, analyzing the power load distribution on tokamak divertors. The current implementation of BIT1 relies exclusively on MPI for parallel communication and lacks support for GPUs. In this work, we address these limitations by designing and implementing a hybrid, shared-memory version of BIT1 capable of utilizing GPUs. For shared-memory parallelization, we rely on OpenMP and OpenACC, using a task-based approach to mitigate load-imbalance issues in the particle mover. On an HPE Cray EX computing node, we observe an initial performance improvement of approximately 42%, with scalable performance showing an enhancement of about 38% when using 8 MPI ranks. Still relying on OpenMP and OpenACC, we introduce the first version of BIT1 capable of using GPUs. We investigate two different data movement strategies: unified memory and explicit data movement. Overall, we report BIT1 data transfer findings during each PIC cycle. Among BIT1 GPU implementations, we demonstrate performance improvement through concurrent GPU utilization, especially when MPI ranks are assigned to dedicated GPUs. Finally, we analyze the performance of the first BIT1 GPU porting with the NVIDIA Nsight tools to further our understanding of BIT1 computational efficiency for large-scale plasma simulations, capable of exploiting current supercomputer infrastructures.
△ Less
Submitted 15 April, 2024;
originally announced April 2024.
-
Physics-based adaptivity of a spectral method for the Vlasov-Poisson equations based on the asymmetrically-weighted Hermite expansion in velocity space
Authors:
Cecilia Pagliantini,
Gian Luca Delzanno,
Stefano Markidis
Abstract:
We propose a spectral method for the 1D-1V Vlasov-Poisson system where the discretization in velocity space is based on asymmetrically-weighted Hermite functions, dynamically adapted via a scaling $α$ and shifting $u$ of the velocity variable. Specifically, at each time instant an adaptivity criterion selects new values of $α$ and $u$ based on the numerical solution of the discrete Vlasov-Poisson…
▽ More
We propose a spectral method for the 1D-1V Vlasov-Poisson system where the discretization in velocity space is based on asymmetrically-weighted Hermite functions, dynamically adapted via a scaling $α$ and shifting $u$ of the velocity variable. Specifically, at each time instant an adaptivity criterion selects new values of $α$ and $u$ based on the numerical solution of the discrete Vlasov-Poisson system obtained at that time step. Once the new values of the Hermite parameters $α$ and $u$ are fixed, the Hermite expansion is updated and the discrete system is further evolved for the next time step. The procedure is applied iteratively over the desired temporal interval. The key aspects of the adaptive algorithm are: the map between approximation spaces associated with different values of the Hermite parameters that preserves total mass, momentum and energy; and the adaptivity criterion to update $α$ and $u$ based on physics considerations relating the Hermite parameters to the average velocity and temperature of each plasma species. For the discretization of the spatial coordinate, we rely on Fourier functions and use the implicit midpoint rule for time stepping. The resulting numerical method possesses intrinsically the property of fluid-kinetic coupling, where the low-order terms of the expansion are akin to the fluid moments of a macroscopic description of the plasma, while kinetic physics is retained by adding more spectral terms. Moreover, the scheme features conservation of total mass, momentum and energy associated in the discrete, for periodic boundary conditions. A set of numerical experiments confirms that the adaptive method outperforms the non-adaptive one in terms of accuracy and stability of the numerical solution.
△ Less
Submitted 31 May, 2023; v1 submitted 30 August, 2022;
originally announced August 2022.
-
Breaking Down the Parallel Performance of GROMACS, a High-Performance Molecular Dynamics Software
Authors:
Måns I. Andersson,
N. Arul Murugan,
Artur Podobas,
Stefano Markidis
Abstract:
GROMACS is one of the most widely used HPC software packages using the Molecular Dynamics (MD) simulation technique. In this work, we quantify GROMACS parallel performance using different configurations, HPC systems, and FFT libraries (FFTW, Intel MKL FFT, and FFT PACK). We break down the cost of each GROMACS computational phase and identify non-scalable stages, such as MPI communication during th…
▽ More
GROMACS is one of the most widely used HPC software packages using the Molecular Dynamics (MD) simulation technique. In this work, we quantify GROMACS parallel performance using different configurations, HPC systems, and FFT libraries (FFTW, Intel MKL FFT, and FFT PACK). We break down the cost of each GROMACS computational phase and identify non-scalable stages, such as MPI communication during the 3D FFT computation when using a large number of processes. We show that the Particle-Mesh Ewald phase and the 3D FFT calculation significantly impact the GROMACS performance. Finally, we discuss performance opportunities with a particular interest in developing GROMACS for the FFT calculations.
△ Less
Submitted 29 August, 2022;
originally announced August 2022.
-
Large-Scale Direct Numerical Simulations of Turbulence Using GPUs and Modern Fortran
Authors:
Martin Karp,
Daniele Massaro,
Niclas Jansson,
Alistair Hart,
Jacob Wahlgren,
Philipp Schlatter,
Stefano Markidis
Abstract:
We present our approach to making direct numerical simulations of turbulence with applications in sustainable shipping. We use modern Fortran and the spectral element method to leverage and scale on supercomputers powered by the Nvidia A100 and the recent AMD Instinct MI250X GPUs, while still providing support for user software developed in Fortran. We demonstrate the efficiency of our approach by…
▽ More
We present our approach to making direct numerical simulations of turbulence with applications in sustainable shipping. We use modern Fortran and the spectral element method to leverage and scale on supercomputers powered by the Nvidia A100 and the recent AMD Instinct MI250X GPUs, while still providing support for user software developed in Fortran. We demonstrate the efficiency of our approach by performing the world's first direct numerical simulation of the flow around a Flettner rotor at Re=30'000 and its interaction with a turbulent boundary layer. We present one of the first performance comparisons between the AMD Instinct MI250X and Nvidia A100 GPUs for scalable computational fluid dynamics. Our results show that one MI250X offers performance on par with two A100 GPUs and has a similar power efficiency.
△ Less
Submitted 23 June, 2022;
originally announced July 2022.
-
A Deep Learning-Based Particle-in-Cell Method for Plasma Simulations
Authors:
Xavier Aguilar,
Stefano Markidis
Abstract:
We design and develop a new Particle-in-Cell (PIC) method for plasma simulations using Deep-Learning (DL) to calculate the electric field from the electron phase space. We train a Multilayer Perceptron (MLP) and a Convolutional Neural Network (CNN) to solve the two-stream instability test. We verify that the DL-based MLP PIC method produces the correct results using the two-stream instability: the…
▽ More
We design and develop a new Particle-in-Cell (PIC) method for plasma simulations using Deep-Learning (DL) to calculate the electric field from the electron phase space. We train a Multilayer Perceptron (MLP) and a Convolutional Neural Network (CNN) to solve the two-stream instability test. We verify that the DL-based MLP PIC method produces the correct results using the two-stream instability: the DL-based PIC provides the expected growth rate of the two-stream instability. The DL-based PIC does not conserve the total energy and momentum. However, the DL-based PIC method is stable against the cold-beam instability, affecting traditional PIC methods. This work shows that integrating DL technologies into traditional computational methods is a viable approach for developing next-generation PIC algorithms.
△ Less
Submitted 5 July, 2021;
originally announced July 2021.
-
The Old and the New: Can Physics-Informed Deep-Learning Replace Traditional Linear Solvers?
Authors:
Stefano Markidis
Abstract:
Physics-Informed Neural Networks (PINN) are neural networks encoding the problem governing equations, such as Partial Differential Equations (PDE), as a part of the neural network. PINNs have emerged as a new essential tool to solve various challenging problems, including computing linear systems arising from PDEs, a task for which several traditional methods exist. In this work, we focus first on…
▽ More
Physics-Informed Neural Networks (PINN) are neural networks encoding the problem governing equations, such as Partial Differential Equations (PDE), as a part of the neural network. PINNs have emerged as a new essential tool to solve various challenging problems, including computing linear systems arising from PDEs, a task for which several traditional methods exist. In this work, we focus first on evaluating the potential of PINNs as linear solvers in the case of the Poisson equation, an omnipresent equation in scientific computing. We characterize PINN linear solvers in terms of accuracy and performance under different network configurations (depth, activation functions, input data set distribution). We highlight the critical role of transfer learning. Our results show that low-frequency components of the solution converge quickly as an effect of the F-principle. In contrast, an accurate solution of the high frequencies requires an exceedingly long time. To address this limitation, we propose integrating PINNs into traditional linear solvers. We show that this integration leads to the development of new solvers whose performance is on par with other high-performance solvers, such as PETSc conjugate gradient linear solvers, in terms of performance and accuracy. Overall, while the accuracy and computational performance are still a limiting factor for the direct use of PINN linear solvers, hybrid strategies combining old traditional linear solver approaches with new emerging deep-learning techniques are among the most promising methods for developing a new class of linear solvers.
△ Less
Submitted 5 July, 2021; v1 submitted 11 March, 2021;
originally announced March 2021.
-
Kinetic Modeling of Magnetospheres
Authors:
Stefano Markidis,
Vyacheslav Olshevsky,
Gabor Toth,
Yuxi Chen,
Ivy Peng,
Giovanni Lapenta,
Tamas Gombosi
Abstract:
This paper presents the state of the art of kinetic modeling techniques for simulating plasma kinetic dynamics in magnetospheres. We describe the critical numerical techniques for enabling large-scale kinetic simulations of magnetospheres: parameter scaling, implicit Particle-in-Cell schemes, and fluid-kinetic coupling. We show an application of these techniques to study particle acceleration and…
▽ More
This paper presents the state of the art of kinetic modeling techniques for simulating plasma kinetic dynamics in magnetospheres. We describe the critical numerical techniques for enabling large-scale kinetic simulations of magnetospheres: parameter scaling, implicit Particle-in-Cell schemes, and fluid-kinetic coupling. We show an application of these techniques to study particle acceleration and heating in asymmetric magnetic reconnection in the Ganymede magnetosphere.
△ Less
Submitted 11 December, 2020;
originally announced December 2020.
-
Automatic Particle Trajectory Classification in Plasma Simulations
Authors:
Stefano Markidis,
Ivy Peng,
Artur Podobas,
Itthinat Jongsuebchoke,
Gabriel Bengtsson,
Pawel Herman
Abstract:
Numerical simulations of plasma flows are crucial for advancing our understanding of microscopic processes that drive the global plasma dynamics in fusion devices, space, and astrophysical systems. Identifying and classifying particle trajectories allows us to determine specific on-going acceleration mechanisms, shedding light on essential plasma processes.
Our overall goal is to provide a gener…
▽ More
Numerical simulations of plasma flows are crucial for advancing our understanding of microscopic processes that drive the global plasma dynamics in fusion devices, space, and astrophysical systems. Identifying and classifying particle trajectories allows us to determine specific on-going acceleration mechanisms, shedding light on essential plasma processes.
Our overall goal is to provide a general workflow for exploring particle trajectory space and automatically classifying particle trajectories from plasma simulations in an unsupervised manner. We combine pre-processing techniques, such as Fast Fourier Transform (FFT), with Machine Learning methods, such as Principal Component Analysis (PCA), k-means clustering algorithms, and silhouette analysis. We demonstrate our workflow by classifying electron trajectories during magnetic reconnection problem. Our method successfully recovers existing results from previous literature without a priori knowledge of the underlying system.
Our workflow can be applied to analyzing particle trajectories in different phenomena, from magnetic reconnection, shocks to magnetospheric flows. The workflow has no dependence on any physics model and can identify particle trajectories and acceleration mechanisms that were not detected before.
△ Less
Submitted 11 October, 2020;
originally announced October 2020.
-
Magnetohydrodynamic with embedded particle-in-cell simulation of the Geospace Environment Modeling dayside kinetic processes challenge event
Authors:
Yuxi Chen,
Gabor Toth,
Heli Hietala,
Sarah Vines,
Ying Zou,
Yukitoshi Nishimura,
Marcos Silveira,
Zhifang Guo,
Yu Lin,
Stefano Markidis
Abstract:
We use the MHD with embedded particle-in-cell model (MHD-EPIC) to study the Geospace Environment Modeling (GEM) dayside kinetic processes challenge event at 01:50-03:00 UT on 2015-11-18, when the magnetosphere was driven by a steady southward IMF. In the MHD-EPIC simulation, the dayside magnetopause is covered by a PIC code so that the dayside reconnection is properly handled. We compare the magne…
▽ More
We use the MHD with embedded particle-in-cell model (MHD-EPIC) to study the Geospace Environment Modeling (GEM) dayside kinetic processes challenge event at 01:50-03:00 UT on 2015-11-18, when the magnetosphere was driven by a steady southward IMF. In the MHD-EPIC simulation, the dayside magnetopause is covered by a PIC code so that the dayside reconnection is properly handled. We compare the magnetic fields and the plasma profiles of the magnetopause crossing with the MMS3 spacecraft observations. Most variables match the observations well in the magnetosphere, in the magnetosheath, and also during the current sheet crossing. The MHD-EPIC simulation produces flux ropes, and we demonstrate that some magnetic field and plasma features observed by the MMS3 spacecraft can be reproduced by a flux rope crossing event. We use an algorithm to automatically identify the reconnection sites from the simulation results. It turns out that there are usually multiple X-lines at the magnetopause. By tracing the locations of the X-lines, we find the typical moving speed of the X-line endpoints is about 70~km/s, which is higher than but still comparable with the ground-based observations.
△ Less
Submitted 13 January, 2020;
originally announced January 2020.
-
Automated classification of plasma regions using 3D particle energy distributions
Authors:
Vyacheslav Olshevsky,
Yuri V. Khotyaintsev,
Ahmad Lalti,
Andrey Divin,
Gian Luca Delzanno,
Sven Anderzen,
Pawel Herman,
Steven W. D. Chien,
Levon Avanov,
Andrew P. Dimmock,
Stefano Markidis
Abstract:
We investigate the properties of the ion sky maps produced by the Dual Ion Spectrometers (DIS) from the Fast Plasma Investigation (FPI). We have trained a convolutional neural network classifier to predict four regions crossed by the MMS on the dayside magnetosphere: solar wind, ion foreshock, magnetosheath, and magnetopause using solely DIS spectrograms. The accuracy of the classifier is >98%. We…
▽ More
We investigate the properties of the ion sky maps produced by the Dual Ion Spectrometers (DIS) from the Fast Plasma Investigation (FPI). We have trained a convolutional neural network classifier to predict four regions crossed by the MMS on the dayside magnetosphere: solar wind, ion foreshock, magnetosheath, and magnetopause using solely DIS spectrograms. The accuracy of the classifier is >98%. We use the classifier to detect mixed plasma regions, in particular to find the bow shock regions. A similar approach can be used to identify the magnetopause crossings and reveal regions prone to magnetic reconnection. Data processing through the trained classifier is fast and efficient and thus can be used for classification for the whole MMS database.
△ Less
Submitted 21 September, 2021; v1 submitted 15 August, 2019;
originally announced August 2019.
-
Studying dawn-dusk asymmetries of Mercury's magnetotail using MHD-EPIC simulations
Authors:
Yuxi Chen,
Gabor Toth,
Xianzhe Jia,
James Slavin,
Weijie Sun,
Stefano Markidis,
Tamas Gombosi,
Jim Raines
Abstract:
MESSENGER has observed a lot of dawn-dusk asymmetries in Mercury's magnetotail, such as the asymmetries of the cross-tail current sheet thickness and the occurrence of flux ropes, dipolarization events and energetic electron injections. In order to obtain a global pictures of Mercury's magnetotail dynamics and the relationship between these asymmetries, we perform global simulations with the magne…
▽ More
MESSENGER has observed a lot of dawn-dusk asymmetries in Mercury's magnetotail, such as the asymmetries of the cross-tail current sheet thickness and the occurrence of flux ropes, dipolarization events and energetic electron injections. In order to obtain a global pictures of Mercury's magnetotail dynamics and the relationship between these asymmetries, we perform global simulations with the magnetohydrodynamics with embedded particle-in-cell (MHD-EPIC) model, where Mercury's magnetotail region is covered by a PIC code. Our simulations show that the dawnside current sheet is thicker, the plasma density is larger, and the electron pressure is higher than the duskside. Under a strong IMF driver, the simulated reconnection sites prefer the dawnside. We also found the dipolarization events and the planetward electron jets are moving dawnward while they are moving towards the planet, so that almost all dipolarization events and high-speed plasma flows concentrate in the dawn sector. The simulation results are consistent with MESSENGER observations.
△ Less
Submitted 14 April, 2019;
originally announced April 2019.
-
Particle-in-Cell Simulations of Plasma Dynamics in Cometary Environment
Authors:
Chaitanya Prasad Sishtla,
Vyacheslav Olshevsky,
Steven W. D. Chien,
Stefano Markidis,
Erwin Laure
Abstract:
We perform and analyze global Particle-in-Cell (PIC) simulations of the interaction between solar wind and an outgassing comet with the goal of studying the plasma kinetic dynamics of a cometary environment. To achieve this, we design and implement a new numerical method in the iPIC3D code to model outgassing from the comet: new plasma particles are ejected from the comet "surface" at each computa…
▽ More
We perform and analyze global Particle-in-Cell (PIC) simulations of the interaction between solar wind and an outgassing comet with the goal of studying the plasma kinetic dynamics of a cometary environment. To achieve this, we design and implement a new numerical method in the iPIC3D code to model outgassing from the comet: new plasma particles are ejected from the comet "surface" at each computational cycle. Our simulations show that a bow shock is formed as a result of the interaction between solar wind and outgassed particles. The analysis of distribution functions for the PIC simulations shows that at the bow shock part of the incoming solar wind, ions are reflected while electrons are heated. This work attempts to reveal kinetic effects in the atmosphere of an outgassing comet using a fully kinetic Particle-in-Cell model.
△ Less
Submitted 28 January, 2019;
originally announced January 2019.
-
PolyPIC: the Polymorphic-Particle-in-Cell Method for Fluid-Kinetic Coupling
Authors:
Stefano Markidis,
Vyacheslav Olshevsky,
Chaitanya Prasad Sishtla,
Steven Wei-der Chien,
Erwin Laure,
Giovanni Lapenta
Abstract:
Particle-in-Cell (PIC) methods are widely used computational tools for fluid and kinetic plasma modeling. While both the fluid and kinetic PIC approaches have been successfully used to target either kinetic or fluid simulations, little was done to combine fluid and kinetic particles under the same PIC framework. This work addresses this issue by proposing a new PIC method, PolyPIC, that uses polym…
▽ More
Particle-in-Cell (PIC) methods are widely used computational tools for fluid and kinetic plasma modeling. While both the fluid and kinetic PIC approaches have been successfully used to target either kinetic or fluid simulations, little was done to combine fluid and kinetic particles under the same PIC framework. This work addresses this issue by proposing a new PIC method, PolyPIC, that uses polymorphic computational particles. In this numerical scheme, particles can be either kinetic or fluid, and fluid particles can become kinetic when necessary, e.g. particles undergoing a strong acceleration. We design and implement the PolyPIC method, and test it against the Landau damping of Langmuir and ion acoustic waves, two stream instability and sheath formation. We unify the fluid and kinetic PIC methods under one common framework comprising both fluid and kinetic particles, providing a tool for adaptive fluid-kinetic coupling in plasma simulations.
△ Less
Submitted 13 July, 2018;
originally announced July 2018.
-
Global three-dimensional simulation of Earth's dayside reconnection using a two-way coupled magnetohydrodynamics with embedded particle-in-cell model: initial results
Authors:
Yuxi Chen,
Gabor Toth,
Paul Cassak,
Xianzhe Jia,
Tamas I. Gombosi,
James A. Slavin,
Stefano Markidis,
Ivy Bo Peng,
Vania K. Jordanova
Abstract:
We perform a three-dimensional (3D) global simulation of Earth's magnetosphere with kinetic reconnection physics to study the flux transfer events (FTEs) and dayside magnetic reconnection with the recently developed magnetohydrodynamics with embedded particle-in-cell model (MHD-EPIC). During the one-hour long simulation, the FTEs are generated quasi-periodically near the subsolar point and move to…
▽ More
We perform a three-dimensional (3D) global simulation of Earth's magnetosphere with kinetic reconnection physics to study the flux transfer events (FTEs) and dayside magnetic reconnection with the recently developed magnetohydrodynamics with embedded particle-in-cell model (MHD-EPIC). During the one-hour long simulation, the FTEs are generated quasi-periodically near the subsolar point and move toward the poles. We find the magnetic field signature of FTEs at their early formation stage is similar to a `crater FTE', which is characterized by a magnetic field strength dip at the FTE center. After the FTE core field grows to a significant value, it becomes an FTE with typical flux rope structure. When an FTE moves across the cusp, reconnection between the FTE field lines and the cusp field lines can dissipate the FTE. The kinetic features are also captured by our model. A crescent electron phase space distribution is found near the reconnection site. A similar distribution is found for ions at the location where the Larmor electric field appears. The lower hybrid drift instability (LHDI) along the current sheet direction also arises at the interface of magnetosheath and magnetosphere plasma. The LHDI electric field is about 8 mV/m and its dominant wavelength relative to the electron gyroradius agrees reasonably with MMS observations.
△ Less
Submitted 12 April, 2017;
originally announced April 2017.
-
On the ions acceleration via collisionless magnetic reconnection in laboratory plasmas
Authors:
Emanuele Cazzola,
Davide Curreli,
Stefano Markidis,
Giovanni Lapenta
Abstract:
This work presents an analysis of the ion outflow from magnetic reconnection throughout fully kinetic simulations with typical laboratory plasmas values. A symmetric initial configuration for the density and magnetic field is considered across the current sheet. After analyzing the behavior of a set of nine simulations with a reduced mass ratio and with a permuted value of three initial electron t…
▽ More
This work presents an analysis of the ion outflow from magnetic reconnection throughout fully kinetic simulations with typical laboratory plasmas values. A symmetric initial configuration for the density and magnetic field is considered across the current sheet. After analyzing the behavior of a set of nine simulations with a reduced mass ratio and with a permuted value of three initial electron temperature and magnetic field intensity, the best ion acceleration scenario is further studied with a realistic mass ratio in terms of the ion dynamics and energy budget. Interestingly, a series of shock waves structures are observed in the outflow, resembling the shock discontinuities found in recent magnetohydrodynamic (MHD) simulations. An analysis of the ion outflow at several distances from the reconnection point is presented, in light of possible laboratory applications. The analysis suggests that magnetic reconnection could be used as a tool for plasma acceleration, with applications ranging from electric propulsion to production of ion thermal beams.
△ Less
Submitted 28 October, 2016;
originally announced October 2016.
-
On the Electron Agyrotropy during Rapid Asymmetric Magnetic Island Coalescence in Presence of a Guide Field
Authors:
Emanuele Cazzola,
Maria Elena Innocenti,
Martin V. Goldman,
David L. Newman,
Stefano Markidis,
Giovanni Lapenta
Abstract:
We present an analysis of the properties of the electron velocity distribution during island coalescence in asymmetric reconnection with and without guide field. In a previous study, three main domains were identified, in the case without guide field, as X-, D- and M-regions featuring different reconnection evolutions {Cazzola et al. 2015). These regions are also identified here in the case with g…
▽ More
We present an analysis of the properties of the electron velocity distribution during island coalescence in asymmetric reconnection with and without guide field. In a previous study, three main domains were identified, in the case without guide field, as X-, D- and M-regions featuring different reconnection evolutions {Cazzola et al. 2015). These regions are also identified here in the case with guide field. We study the departure from isotropic and gyrotropic behavior by means of different robust detection algorithms proposed in the literature. While in the case without guide field these metrics show an overall agreement, when the guide field is present a discrepancy in the agyrotropy within some relevant regions is observed, such as at the separatrices and inside magnetic islands. Moreover, in light of the new observations from the Multiscale MagnetoSpheric mission, an analysis of the electron velocity phase-space in these domains is presented.
△ Less
Submitted 2 August, 2016;
originally announced August 2016.
-
Energy exchanges in reconnection outflows
Authors:
Giovanni Lapenta,
Martin V. Goldman,
Davd L. Newman,
Stefano Markidis
Abstract:
Reconnection outflows are highly energetic directed flows that interact with the ambient plasma or with flows from other reconnection regions. Under these conditions the flow becomes highly unstable and chaotic, as any flow jets interacting with a medium. We report here massively parallel simulations of the two cases of interaction between outflow jets and between a single outflow with an ambient…
▽ More
Reconnection outflows are highly energetic directed flows that interact with the ambient plasma or with flows from other reconnection regions. Under these conditions the flow becomes highly unstable and chaotic, as any flow jets interacting with a medium. We report here massively parallel simulations of the two cases of interaction between outflow jets and between a single outflow with an ambient plasma. We find in both case the development of a chaotic magnetic field, subject to secondary reconnection events that further complicate the topology of the field lines. The focus of the present analysis is on the energy balance. We compute each energy channel (electromagnetic, bulk, thermal, for each species) and find where the most energy is exchanged and in what form. The main finding is that the largest energy exchange is not at the reconnection site proper but in the regions where the outflowing jets are destabilized.
△ Less
Submitted 2 July, 2016;
originally announced July 2016.
-
Magnetic null points in kinetic simulations of space plasmas
Authors:
Vyacheslav Olshevsky,
Jan Deca,
Andrey Divin,
Ivy Bo Peng,
Stefano Markidis,
Maria Elena Innocenti,
Emanuele Cazzola,
Giovanni Lapenta
Abstract:
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 interactin…
▽ More
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.
△ Less
Submitted 7 December, 2015;
originally announced December 2015.
-
On the electron dynamics during island coalescence in asymmetric magnetic reconnection
Authors:
Emanuele Cazzola,
Maria Elena Innocenti,
Stefano Markidis,
Martin V. Goldman,
David L. Newman,
Giovanni Lapenta
Abstract:
We present an analysis of the electron dynamics during rapid island merging in asymmetric magnetic reconnection. We consider a doubly periodic system with two asymmetric transitions. The upper layer is an asymmetric Harris sheet initially perturbed to promote a single reconnection site. The lower layer is a tangential discontinuity that promotes the formation of many X-points, separated by rapidly…
▽ More
We present an analysis of the electron dynamics during rapid island merging in asymmetric magnetic reconnection. We consider a doubly periodic system with two asymmetric transitions. The upper layer is an asymmetric Harris sheet initially perturbed to promote a single reconnection site. The lower layer is a tangential discontinuity that promotes the formation of many X-points, separated by rapidly merging islands. Across both layers the magnetic field and the density have a strong jump, but the pressure is held constant. Our analysis focuses on the consequences of electron energization during island coalescence. We focus first on the parallel and perpendicular components of the electron temperature to establish the presence of possible anisotropies and non-gyrotropies. Thanks to the direct comparison between the two different layers simulated, we can distinguish three main types of behavior characteristic of three different regions of interest. The first type represents the regions where traditional asymmetric reconnections take place without involving island merging. The second type of regions instead show reconnection events between two merging islands. Finally, the third regions identifies the regions between two diverging island and where typical signature of reconnection is not observed. Electrons in these latter regions additionally show a flat-top distribution resulting from the saturation of a two-stream instability generated by the two interacting electron beams from the two nearest reconnection points. Finally, the analysis of agyrotropy shows the presence of a distinct double structure laying all over the lower side facing the higher magnetic field region. The distinguishing features found for the three types of regions investigated provide clear indicators to the recently launched MMS NASA mission for investigating magnetopause reconnection involving multiple islands.
△ Less
Submitted 18 November, 2015;
originally announced November 2015.
-
Where should MMS look for electron diffusion regions?
Authors:
G. Lapenta,
M. Goldman,
D. Newman,
S. Markidis
Abstract:
A great possible achievement for the MMS mission would be crossing electron diffusion regions (EDR). EDR are regions in proximity of reconnection sites where electrons decouple from field lines, breaking the frozen in condition. Decades of research on reconnection have produced a widely shared map of where EDRs are. We expect reconnection to take place around a so called x-point formed by the inte…
▽ More
A great possible achievement for the MMS mission would be crossing electron diffusion regions (EDR). EDR are regions in proximity of reconnection sites where electrons decouple from field lines, breaking the frozen in condition. Decades of research on reconnection have produced a widely shared map of where EDRs are. We expect reconnection to take place around a so called x-point formed by the intersection of the separatrices dividing inflowing from outflowing plasma. The EDR forms around this x-point as a small electron scale box nested inside a larger ion diffusion region. But this point of view is based on a 2D mentality. We have recently proposed that once the problem is considered in full 3D, secondary reconnection events can form [Lapenta et al., Nature Physics, 11, 690, 2015] in the outflow regions even far downstream from the primary reconnection site. We revisit here this new idea confirming that even using additional indicators of reconnection and even considering longer periods and wider distances the conclusion remains true: secondary reconnection sites form downstream of a reconnection outflow causing a sort of chain reaction of cascading reconnection sites. If we are right, MMS will have an interesting journey even when not crossing necessarily the primary site. The chances are greatly increased that even if missing a primary site during an orbit, MMS could stumble instead on one of these secondary sites.
△ Less
Submitted 13 October, 2015;
originally announced October 2015.
-
Energetics of kinetic reconnection in a three-dimensional null points cluster
Authors:
Vyacheslav Olshevsky,
Giovanni Lapenta,
Stefano Markidis
Abstract:
We performed three-dimensional Particle-in-Cell simulations of magnetic reconnection with multiple magnetic null points. Magnetic field energy conversion into kinetic energy was about five times higher than in traditional Harris sheet configuration. More than 85% of initial magnetic field energy was transferred to particle energy during 25 reversed ion cyclofrequencies. Magnetic reconnection in th…
▽ More
We performed three-dimensional Particle-in-Cell simulations of magnetic reconnection with multiple magnetic null points. Magnetic field energy conversion into kinetic energy was about five times higher than in traditional Harris sheet configuration. More than 85% of initial magnetic field energy was transferred to particle energy during 25 reversed ion cyclofrequencies. Magnetic reconnection in the cluster of null points evolved in three phases. During the first phase, ion beams were excited, that then gave part of their energy back to 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 low magnetic field. Magnetic reconnection in null points presents essentially three-dimensional features, with no two dimensional symmetries or current sheets.
△ Less
Submitted 26 September, 2015;
originally announced September 2015.
-
Role of Z-pinches in magnetic reconnection in space plasmas
Authors:
Vyacheslav Olshevsky,
Giovanni Lapenta,
Stefano Markidis,
Andrey Divin
Abstract:
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…
▽ More
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.
△ Less
Submitted 26 September, 2015;
originally announced September 2015.
-
Energy dissipation in magnetic null points at kinetic scales
Authors:
Vyacheslav Olshevsky,
Andrey Divin,
Elin Eriksson,
Stefano Markidis,
Giovanni Lapenta
Abstract:
We use kinetic particle-in-cell and magnetohydrodynamic simulations supported by an observational dataset to investigate magnetic reconnection in clusters of null points in space plasma. The magnetic configuration under investigation is driven by fast adiabatic flux rope compression that dissipates almost half of the initial magnetic field energy. In this phase powerful currents are excited produc…
▽ More
We use kinetic particle-in-cell and magnetohydrodynamic simulations supported by an observational dataset to investigate magnetic reconnection in clusters of null points in space plasma. The magnetic configuration under investigation is driven by fast adiabatic flux rope compression that dissipates almost half of the initial magnetic field energy. In this phase powerful currents are excited producing secondary instabilities, and the system is brought into a state of `intermittent turbulence' within a few ion gyro-periods. Reconnection events are distributed all over the simulation domain and energy dissipation is rather volume-filling. Numerous spiral null points interconnected via their spines form null lines embedded into magnetic flux ropes; null point pairs demonstrate the signatures of torsional spine reconnection. However, energy dissipation mainly happens in the shear layers formed by adjacent flux ropes with oppositely directed currents. In these regions radial null pairs are spontaneously emerging and vanishing, associated with electron streams and small-scale current sheets. The number of spiral nulls in the simulation outweighs the number of radial nulls by a factor of 5\---10, in accordance with Cluster observations in the Earth's magnetosheath. Twisted magnetic fields with embedded spiral null points might indicate the regions of major energy dissipation for future space missions such as Magnetospheric Multiscale Mission (MMS).
△ Less
Submitted 26 September, 2015;
originally announced September 2015.
-
Kinetic plasma turbulence during the nonlinear stage of the Kelvin-Helmholtz instability
Authors:
Koen Kemel,
Pierre Henri,
Giovanni Lapenta,
Francesco Califano,
Stefano Markidis
Abstract:
Using a full kinetic, implicit particle-in-cell code, iPiC3D, we studied the properties of plasma kinetic turbulence, such as would be found at the interface between the solar wind and the Earth magnetosphere at low latitude during northwards periods. In this case, in the presence of a magnetic field B oriented mostly perpendicular to the velocity shear, turbulence is fed by the disruption of a Ke…
▽ More
Using a full kinetic, implicit particle-in-cell code, iPiC3D, we studied the properties of plasma kinetic turbulence, such as would be found at the interface between the solar wind and the Earth magnetosphere at low latitude during northwards periods. In this case, in the presence of a magnetic field B oriented mostly perpendicular to the velocity shear, turbulence is fed by the disruption of a Kelvin-Helmholtz vortex chain via secondary instabilities, vortex pairing and non-linear interactions.
We found that the magnetic energy spectral cascade between ion and electron inertial scales, $d_i$ and $d_e$, is in agreement with satellite observations and other previous numerical simulations; however, in our case the spectrum ends with a peak beyond $d_e$ due to the occurrence of the lower hybrid drift instability. The electric energy spectrum is influenced by effects of secondary instabilities: anomalous resistivity, fed by the development of the lower hybrid drift instability, steepens the spectral decay and, depending on the alignment or anti-alignment of B and the shear vorticity, peaks due to ion-Bernstein waves may dominate the spectrum around $d_i$. These waves are generated by counter-streaming flow structures, through flux freezing also responsible for reconnection of the in-plane component of the magnetic field, which then generates electron pressure anisotropy and flattening of the field-aligned component of the electron distribution function.
△ Less
Submitted 3 December, 2014;
originally announced December 2014.
-
Signatures of Secondary Collisionless Magnetic Reconnection Driven by Kink Instability of a Flux Rope
Authors:
S. Markidis,
G. Lapenta,
G. L. Delzanno,
P. Henri,
M. V. Goldman,
D. L. Newman,
T. Intrator,
E. Laure
Abstract:
The kinetic features of secondary magnetic reconnection in a single flux rope undergoing internal kink instability are studied by means of three-dimensional Particle-in-Cell simulations. Several signatures of secondary magnetic reconnection are identified in the plane perpendicular to the flux rope: a quadrupolar electron and ion density structure and a bipolar Hall magnetic field develop in proxi…
▽ More
The kinetic features of secondary magnetic reconnection in a single flux rope undergoing internal kink instability are studied by means of three-dimensional Particle-in-Cell simulations. Several signatures of secondary magnetic reconnection are identified in the plane perpendicular to the flux rope: a quadrupolar electron and ion density structure and a bipolar Hall magnetic field develop in proximity of the reconnection region. The most intense electric fields form perpendicularly to the local magnetic field, and a reconnection electric field is identified in the plane perpendicular to the flux rope. An electron current develops along the reconnection line in the opposite direction of the electron current supporting the flux rope magnetic field structure. Along the reconnection line, several bipolar structures of the electric field parallel to the magnetic field occur making the magnetic reconnection region turbulent. The reported signatures of secondary magnetic reconnection can help to localize magnetic reconnection events in space, astrophysical and fusion plasmas.
△ Less
Submitted 5 August, 2014;
originally announced August 2014.
-
Separatrices: the crux of reconnection
Authors:
Giovanni Lapenta,
Stefano Markidis,
Andrey Divin,
David Newman,
Martin Goldman
Abstract:
Reconnection is one of the key processes in astrophysical and laboratory plasmas: it is the opposite of a dynamo. Looking at energy, a dynamo transforms kinetic energy in magnetic energy while reconnection takes magnetic energy and returns is to its kinetic form. Most plasma processes at their core involve first storing magnetic energy accumulated over time and then releasing it suddenly. We focus…
▽ More
Reconnection is one of the key processes in astrophysical and laboratory plasmas: it is the opposite of a dynamo. Looking at energy, a dynamo transforms kinetic energy in magnetic energy while reconnection takes magnetic energy and returns is to its kinetic form. Most plasma processes at their core involve first storing magnetic energy accumulated over time and then releasing it suddenly. We focus here on this release. A key concept in analysing reconnection is that of the separatrix, a surface (line in 2D) that separates the fresh unperturbed plasma embedded in magnetic field lines not yet reconnected with the hotter exhaust embedded in reconnected field lines. In kinetic physics, the separatrices become a layer where many key processes develop. We present here new results relative to the processes at the separatrices that regulate the plasma flow, the energisation of the species, the electromagnetic fields and the instabilities developing at the separatrices.
△ Less
Submitted 24 June, 2014;
originally announced June 2014.
-
Electromagnetic Energy Conversion in Downstream Fronts from 3D Kinetic Reconnection
Authors:
Giovanni Lapenta,
Martin Goldman,
David Newman,
Stefano Markidis,
Andrey Divin
Abstract:
The electromagnetic energy equation is analyzed term by term in a 3D simulation of kinetic reconnection previously reported by \citet{vapirev2013formation}. The evolution presents the usual 2D-like topological structures caused by an initial perturbation independent of the third dimension. However, downstream of the reconnection site, where the jetting plasma encounters the yet unperturbed pre-exi…
▽ More
The electromagnetic energy equation is analyzed term by term in a 3D simulation of kinetic reconnection previously reported by \citet{vapirev2013formation}. The evolution presents the usual 2D-like topological structures caused by an initial perturbation independent of the third dimension. However, downstream of the reconnection site, where the jetting plasma encounters the yet unperturbed pre-existing plasma, a downstream front (DF) is formed and made unstable by the strong density gradient and the unfavorable local acceleration field. The energy exchange between plasma and fields is most intense at the instability, reaching several $pW/m^3$, alternating between load (energy going from fields to particles) and generator (energy going from particles to fields) regions. Energy exchange is instead purely that of a load at the reconnection site itself in a region focused around the x-line and elongated along the separatrix surfaces.
Poynting fluxes are generated at all energy exchange regions and travel away from the reconnection site transporting an energy signal of the order of about $\mathbf S \approx 10^{-3} W/m^2$.
△ Less
Submitted 1 February, 2014;
originally announced February 2014.
-
Nonlinear evolution of the magnetized Kelvin-Helmholtz instability: from fluid to kinetic modeling
Authors:
P. Henri,
S. S. Cerri,
F. Califano,
F. Pegoraro,
C. Rossi,
M. Faganello,
O. Šebek,
P. M. Trávníček,
P. Hellinger,
J. T. Frederiksen,
Å. Nordlund,
S. Markidis,
R. Keppens,
G. Lapenta
Abstract:
The nonlinear evolution of collisionless plasmas is typically a multi-scale process where the energy is injected at large, fluid scales and dissipated at small, kinetic scales. Accurately modelling the global evolution requires to take into account the main micro-scale physical processes of interest. This is why comparison of different plasma models is today an imperative task aiming at understand…
▽ More
The nonlinear evolution of collisionless plasmas is typically a multi-scale process where the energy is injected at large, fluid scales and dissipated at small, kinetic scales. Accurately modelling the global evolution requires to take into account the main micro-scale physical processes of interest. This is why comparison of different plasma models is today an imperative task aiming at understanding cross-scale processes in plasmas. We report here the first comparative study of the evolution of a magnetized shear flow, through a variety of different plasma models by using magnetohydrodynamic, Hall-MHD, two-fluid, hybrid kinetic and full kinetic codes. Kinetic relaxation effects are discussed to emphasize the need for kinetic equilibriums to study the dynamics of collisionless plasmas in non trivial configurations. Discrepancies between models are studied both in the linear and in the nonlinear regime of the magnetized Kelvin-Helmholtz instability, to highlight the effects of small scale processes on the nonlinear evolution of collisionless plasmas. We illustrate how the evolution of a magnetized shear flow depends on the relative orientation of the fluid vorticity with respect to the magnetic field direction during the linear evolution when kinetic effects are taken into account. Even if we found that small scale processes differ between the different models, we show that the feedback from small, kinetic scales to large, fluid scales is negligable in the nonlinear regime. This study show that the kinetic modeling validates the use of a fluid approach at large scales, which encourages the development and use of fluid codes to study the nonlinear evolution of magnetized fluid flows, even in the colisionless regime.
△ Less
Submitted 29 October, 2013;
originally announced October 2013.
-
Acceleration of solar wind ions to 1 MeV by electromagnetic structures upstream of the Earth's bow shock
Authors:
K. Stasiewicz,
S. Markidis,
B. Eliasson,
M. Strumik,
M. Yamauchi
Abstract:
We present measurements from the ESA/NASA Cluster mission that show in situ acceleration of ions to energies of 1 MeV outside the bow shock. The observed heating can be associated with the presence of electromagnetic structures with strong spatial gradients of the electric field that lead to ion gyro-phase breaking and to the onset of chaos in ion trajectories. It results in rapid, stochastic acce…
▽ More
We present measurements from the ESA/NASA Cluster mission that show in situ acceleration of ions to energies of 1 MeV outside the bow shock. The observed heating can be associated with the presence of electromagnetic structures with strong spatial gradients of the electric field that lead to ion gyro-phase breaking and to the onset of chaos in ion trajectories. It results in rapid, stochastic acceleration of ions in the direction perpendicular to the ambient magnetic field. The electric potential of the structures can be compared to a field of moguls on a ski slope, capable of accelerating and ejecting the fast running skiers out of piste. This mechanism may represent the universal mechanism for perpendicular acceleration and heating of ions in the magnetosphere, the solar corona and in astrophysical plasmas. This is also a basic mechanism that can limit steepening of nonlinear electromagnetic structures at shocks and foreshocks in collisionless plasmas.
△ Less
Submitted 5 July, 2013;
originally announced July 2013.
-
The Fluid-Kinetic Particle-in-Cell Solver for Plasma Simulations
Authors:
Stefano Markidis,
Pierre Henri,
Giovanni Lapenta,
Kjell Ronnmark,
Maria Hamrin,
Zakaria Meliani,
Erwin Laure
Abstract:
A new method that solves concurrently the multi-fluid and Maxwell's equations has been developed for plasma simulations. By calculating the stress tensor in the multi-fluid momentum equation by means of computational particles moving in a self-consistent electromagnetic field, the kinetic effects are retained while solving the multi-fluid equations. The Maxwell's and multi-fluid equations are disc…
▽ More
A new method that solves concurrently the multi-fluid and Maxwell's equations has been developed for plasma simulations. By calculating the stress tensor in the multi-fluid momentum equation by means of computational particles moving in a self-consistent electromagnetic field, the kinetic effects are retained while solving the multi-fluid equations. The Maxwell's and multi-fluid equations are discretized implicitly in time enabling kinetic simulations over time scales typical of the fluid simulations. The fluid-kinetic Particle-in-Cell solver has been implemented in a three-dimensional electromagnetic code, and tested against the ion cyclotron resonance and magnetic reconnection problems. The new method is a promising approach for coupling fluid and kinetic methods in a unified framework.
△ Less
Submitted 5 June, 2013;
originally announced June 2013.
-
Kinetic Simulations of Plasmoid Chain Dynamics
Authors:
Stefano Markidis,
Pierre Henri,
Giovanni Lapenta,
Andrey Divin,
Martin Goldman,
David Newman,
Erwin Laure
Abstract:
The dynamics of a plasmoid chain is studied with three dimensional Particle-in-Cell simulations. The evolution of the system with and without a uniform guide field, whose strength is 1/3 the asymptotic magnetic field, is investigated. The plasmoid chain forms by spontaneous magnetic reconnection: the tearing instability rapidly disrupts the initial current sheet generating several small-scale plas…
▽ More
The dynamics of a plasmoid chain is studied with three dimensional Particle-in-Cell simulations. The evolution of the system with and without a uniform guide field, whose strength is 1/3 the asymptotic magnetic field, is investigated. The plasmoid chain forms by spontaneous magnetic reconnection: the tearing instability rapidly disrupts the initial current sheet generating several small-scale plasmoids, that rapidly grow in size coalescing and kinking. The plasmoid kink is mainly driven by the coalescence process. It is found that the presence of guide field strongly influences the evolution of the plasmoid chain. Without a guide field, a main reconnection site dominates and smaller reconnection regions are included in larger ones, leading to an hierarchical structure of the plasmoid-dominated current sheet. On the contrary in presence of a guide field, plasmoids have approximately the same size and the hierarchical structure does not emerge, a strong core magnetic field develops in the center of the plasmoid in the direction of the existing guide field, and bump-on-tail instability, leading to the formation of electron holes, is detected in proximity of the plasmoids.
△ Less
Submitted 5 June, 2013;
originally announced June 2013.
-
Rethinking Electrostatic Solvers in Particle Simulations for the Exascale Era
Authors:
Stefano Markidis,
Giovanni Lapenta,
Rossen Apostolov,
Erwin Laure
Abstract:
In preparation to the exascale era, an alternative approach to calculate the electrostatic forces in Particle Mesh (PM) methods is proposed. While the traditional techniques are based on the calculation of the electrostatic potential by solving the Poisson equation, in the new approach the electric field is calculated by solving the Ampere's law. When the Ampere's law is discretized explicitly in…
▽ More
In preparation to the exascale era, an alternative approach to calculate the electrostatic forces in Particle Mesh (PM) methods is proposed. While the traditional techniques are based on the calculation of the electrostatic potential by solving the Poisson equation, in the new approach the electric field is calculated by solving the Ampere's law. When the Ampere's law is discretized explicitly in time, the electric field values on the mesh are simply updated from the previous values. In this way, the electrostatic solver becomes an embarrassingly parallel problem, making the algorithm extremely scalable and suitable for exascale computing platforms. An implementation of a one dimensional PM code is presented to show that the proposed method produces correct results, and it is a very promising algorithm for exascale PM simulations.
△ Less
Submitted 28 May, 2012; v1 submitted 10 May, 2012;
originally announced May 2012.
-
Three dimensional density cavities in guide field collisionless magnetic reconnection
Authors:
Stefano Markidis,
Giovanni Lapenta,
Andrey Divin,
Martin V. Goldman,
D. Newman,
Laila Andersson
Abstract:
Particle-in-Cell simulations of collisionless magnetic reconnection with a guide field reveal for the first time the three dimensional features of the low density regions along the magnetic reconnection separatrices, the so-called "cavities". It is found that structures with further lower density develop within the cavities. Because their appearance is similar to the rib shape, these formations ar…
▽ More
Particle-in-Cell simulations of collisionless magnetic reconnection with a guide field reveal for the first time the three dimensional features of the low density regions along the magnetic reconnection separatrices, the so-called "cavities". It is found that structures with further lower density develop within the cavities. Because their appearance is similar to the rib shape, these formations are here called "low density ribs". Their location remains approximately fixed in time and their density progressively decreases, as electron currents along the cavities evacuate them. They develop along the magnetic field lines and are supported by a strong perpendicular electric field that oscillates in space. In addition, bipolar parallel electric field structures form as isolated spheres between the cavities and the outflow plasma, along the direction of the low density ribs and of magnetic field lines.
△ Less
Submitted 10 March, 2012;
originally announced March 2012.
-
Collisionless magnetic reconnection in a plasmoid chain
Authors:
Stefano Markidis,
Pierre Henri,
Giovanni Lapenta,
Andrey Divin,
Martin V. Goldman,
David Newman,
Stefan Eriksson
Abstract:
The kinetic features of plasmoid chain formation and evolution are investigated by two dimensional Particle-in-Cell simulations. Magnetic reconnection is initiated in multiple X points by the tearing instability. Plasmoids form and grow in size by continuously coalescing. Each chain plasmoid exhibits a strong out-of plane core magnetic field and an out-of-plane electron current that drives the coa…
▽ More
The kinetic features of plasmoid chain formation and evolution are investigated by two dimensional Particle-in-Cell simulations. Magnetic reconnection is initiated in multiple X points by the tearing instability. Plasmoids form and grow in size by continuously coalescing. Each chain plasmoid exhibits a strong out-of plane core magnetic field and an out-of-plane electron current that drives the coalescing process. The disappearance of the X points in the coalescence process are due to anti-reconnection, a magnetic reconnection where the plasma inflow and outflow are reversed with respect to the original reconnection flow pattern. Anti-reconnection is characterized by the Hall magnetic field quadrupole signature. Two new kinetic features, not reported by previous studies of plasmoid chain evolution, are here revealed. First, intense electric fields develop in-plane normally to the separatrices and drive the ion dynamics in the plasmoids. Second, several bipolar electric field structures are localized in proximity of the plasmoid chain. The analysis of the electron distribution function and phase space reveals the presence of counter-streaming electron beams, unstable to the two stream instability, and phase space electron holes along the reconnection separatrices.
△ Less
Submitted 13 February, 2012;
originally announced February 2012.
-
A Multi Level Multi Domain Method for Particle In Cell Plasma Simulations
Authors:
M. E. Innocenti,
G. Lapenta,
S. Markidis,
A. Beck,
A. Vapirev
Abstract:
A novel adaptive technique for electromagnetic Particle In Cell (PIC) plasma simulations is presented here. Two main issues are identified in designing adaptive techniques for PIC simulation: first, the choice of the size of the particle shape function in progressively refined grids, with the need to avoid the exertion of self-forces on particles, and, second, the necessity to comply with the stri…
▽ More
A novel adaptive technique for electromagnetic Particle In Cell (PIC) plasma simulations is presented here. Two main issues are identified in designing adaptive techniques for PIC simulation: first, the choice of the size of the particle shape function in progressively refined grids, with the need to avoid the exertion of self-forces on particles, and, second, the necessity to comply with the strict stability constraints of the explicit PIC algorithm. The adaptive implementation presented responds to these demands with the introduction of a Multi Level Multi Domain (MLMD) system (where a cloud of self-similar domains is fully simulated with both fields and particles) and the use of an Implicit Moment PIC method as baseline algorithm for the adaptive evolution. Information is exchanged between the levels with the projection of the field information from the refined to the coarser levels and the interpolation of the boundary conditions for the refined levels from the coarser level fields. Particles are bound to their level of origin and are prevented from transitioning to coarser levels, but are repopulated at the refined grid boundaries with a splitting technique. The presented algorithm is tested against a series of simulation challenges.
△ Less
Submitted 30 January, 2012;
originally announced January 2012.
-
Bipolar Electric Field Signatures of Reconnection Separatrices for a Hydrogen Plasma at Realistic Guide Fields
Authors:
G. Lapenta,
S. Markidis,
A. Divin,
M. Goldman,
D. Newman
Abstract:
In preparation for the MMS mission we ask the question: how common are bipolar signatures linked to the presence of electron holes along separatrices emanating from reconnection regions? To answer this question, we conduct massively parallel simulations for realistic conditions and for the hydrogen mass ratio in boxes larger than considered in similar previous studies.
The magnetic field configu…
▽ More
In preparation for the MMS mission we ask the question: how common are bipolar signatures linked to the presence of electron holes along separatrices emanating from reconnection regions? To answer this question, we conduct massively parallel simulations for realistic conditions and for the hydrogen mass ratio in boxes larger than considered in similar previous studies.
The magnetic field configuration includes both a field reversal and a out of plane guide field, as typical of many space situations. The guide field is varied in strength from low values (typical of the Earth magnetotail) to high values comparable to the in plane reconnecting field (as in the magnetopause). In all cases, along the separatrices a strong electron flow is observed, sufficient to lead to the onset of streaming instabilities and to form bipolar parallel electric field signatures. The presence of bipolar structures at all guide fields allows the control of the MMS mission to consider the presence of bipolar signatures as a general flag of the presence of a nearby reconnection site both in the nightside and in the dayside of the magnetosphere.
△ Less
Submitted 11 August, 2011;
originally announced August 2011.
-
The Energy Conserving Particle-in-Cell Method
Authors:
Stefano Markidis,
Giovanni Lapenta
Abstract:
A new Particle-in-Cell (PIC) method, that conserves energy exactly, is presented. The particle equations of motion and the Maxwell's equations are differenced implicitly in time by the midpoint rule and solved concurrently by a Jacobian-free Newton Krylov (JFNK) solver. Several tests show that the finite grid instability is eliminated in energy conserving PIC simulations, and the method correctly…
▽ More
A new Particle-in-Cell (PIC) method, that conserves energy exactly, is presented. The particle equations of motion and the Maxwell's equations are differenced implicitly in time by the midpoint rule and solved concurrently by a Jacobian-free Newton Krylov (JFNK) solver. Several tests show that the finite grid instability is eliminated in energy conserving PIC simulations, and the method correctly describes the two-stream and Weibel instabilities, conserving exactly the total energy. The computational time of the energy conserving PIC method increases linearly with the number of particles, and it is rather insensitive to the number of grid points and time step. The kinetic enslavement technique can be effectively used to reduce the problem matrix size and the number of JFNK solver iterations.
△ Less
Submitted 9 August, 2011;
originally announced August 2011.
-
Development and performance analysis of a UPC Particle-in-Cell code
Authors:
Stefano Markidis,
Giovanni Lapenta
Abstract:
The development and the implementation of a Particle-in-Cell code written in the Unified Parallel C (UPC) language for plasma simulations with application to astrophysics and fusion nuclear energy machines are presented. A simple one dimensional electrostatic Particle-in-Cell code has been developed first to investigate the implementation details in the UPC language, and second to study the UPC pe…
▽ More
The development and the implementation of a Particle-in-Cell code written in the Unified Parallel C (UPC) language for plasma simulations with application to astrophysics and fusion nuclear energy machines are presented. A simple one dimensional electrostatic Particle-in-Cell code has been developed first to investigate the implementation details in the UPC language, and second to study the UPC performance on parallel computers. The initial simulations of plasmas with the UPC Particle-in-Cell code and a study of parallel speed-up of the UPC code up to 128 cores are shown.
△ Less
Submitted 9 August, 2011;
originally announced August 2011.
-
Kinetic simulations of magnetic reconnection in presence of a background O+ population
Authors:
Stefano Markidis,
Giovanni Lapenta,
Lapo Bettarini,
Martin V. Goldman,
David Newman,
Laila Andersson
Abstract:
Particle-in-Cell simulations of magnetic reconnection with an H+ current sheet and a mixed background plasma of H+ and O+ ions are completed using physical mass ratios. Four main results are shown. First, the O+ presence slightly decreases the reconnection rate and the magnetic reconnection evolution depends mainly on the lighter H+ ion species in the presented simulations. Second, the Hall magnet…
▽ More
Particle-in-Cell simulations of magnetic reconnection with an H+ current sheet and a mixed background plasma of H+ and O+ ions are completed using physical mass ratios. Four main results are shown. First, the O+ presence slightly decreases the reconnection rate and the magnetic reconnection evolution depends mainly on the lighter H+ ion species in the presented simulations. Second, the Hall magnetic field is characterized by a two-scale structure in presence of O+ ions: it reaches sharp peak values in a small area in proximity of the neutral line, and then decreases slowly over a large region. Third, the two background species initially separate in the outflow region because H+ and O+ ions are accelerated by different mechanisms occurring on different time scales and with different strengths. Fourth, the effect of a guide field on the O+ dynamics is studied: the O+ presence does not change the reconnected flux and all the characteristic features of guide field magnetic reconnection are still present. Moreover, the guide field introduces an O+ circulation pattern between separatrices that enhances high O+ density areas and depletes low O+ density regions in proximity of the reconnection fronts. The importance and the validity of these results are finally discussed.
△ Less
Submitted 9 August, 2011;
originally announced August 2011.
-
Jet deflection by very weak guide fields during magnetic reconnection
Authors:
Martin Goldman,
Giovanni Lapenta,
David Newman,
Stefan Markidis,
Haihong Che
Abstract:
Simulations of antiparallel reconnection have shown collimated electron jets outflowing from the x-point, and associated highly elongated "outer electron diffusion regions." New PIC simulations with ion/electron mass ratios as large as 1836 show that jets are deflected towards the magnetic separatrix by out-of-plane guide fields, Bg, as small as 0.05 times the asymptotic reconnecting field, B0. Th…
▽ More
Simulations of antiparallel reconnection have shown collimated electron jets outflowing from the x-point, and associated highly elongated "outer electron diffusion regions." New PIC simulations with ion/electron mass ratios as large as 1836 show that jets are deflected towards the magnetic separatrix by out-of-plane guide fields, Bg, as small as 0.05 times the asymptotic reconnecting field, B0. The outer electron diffusion region is distorted and broken up, but the diffusion rate is unchanged. These results are interpreted in terms of electron dynamics and are compared to recent measurements of reconnection jets in the magnetosheath.
△ Less
Submitted 30 October, 2010;
originally announced November 2010.