Subscribe to the gold package and get unlimited access to Shamra Academy
Register a new userPegasus: A New Hybrid-Kinetic Particle-in-Cell Code for Astrophysical Plasma Dynamics
234
0
0.0
(
0
)
Authors
Matthew W. Kunz
Ask ChatGPT about the research
No Arabic abstract
We describe Pegasus, a new hybrid-kinetic particle-in-cell code tailored for the study of astrophysical plasma dynamics. The code incorporates an energy-conserving particle integrator into a stable, second-order--accurate, three-stage predictor-predictor-corrector integration algorithm. The constrained transport method is used to enforce the divergence-free constraint on the magnetic field. A delta-f scheme is included to facilitate a reduced-noise study of systems in which only small departures from an initial distribution function are anticipated. The effects of rotation and shear are implemented through the shearing-sheet formalism with orbital advection. These algorithms are embedded within an architecture similar to that used in the popular astrophysical magnetohydrodynamics code Athena, one that is modular, well-documented, easy to use, and efficiently parallelized for use on thousands of processors. We present a series of tests in one, two, and three spatial dimensions that demonstrate the fidelity and versatility of the code.
rate research
Read More
We present the Photon-Plasma code, a modern high order charge conserving particle-in-cell code for simulating relativistic plasmas. The code is using a high order implicit field solver and a novel high order charge conserving interpolation scheme for particle-to-cell interpolation and charge deposition. It includes powerful diagnostics tools with on-the-fly particle tracking, synthetic spectra integration, 2D volume slicing, and a new method to correctly account for radiative cooling in the simulations. A robust technique for imposing (time-dependent) particle and field fluxes on the boundaries is also presented. Using a hybrid OpenMP and MPI approach the code scales efficiently from 8 to more than 250.000 cores with almost linear weak scaling on a range of architectures. The code is tested with the classical benchmarks particle heating, cold beam instability, and two-stream instability. We also present particle-in-cell simulations of the Kelvin-Helmholtz instability, and new results on radiative collisionless shocks.
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.
SMILEI is a collaborative, open-source, object-oriented (C++) particle-in-cell code. To benefit from the latest advances in high-performance computing (HPC), SMILEI is co-developed by both physicists and HPC experts. The codes structures, capabilities, parallelization strategy and performances are discussed. Additional modules (e.g. to treat ionization or collisions), benchmarks and physics highlights are also presented. Multi-purpose and evolutive, SMILEI is applied today to a wide range of physics studies, from relativistic laser-plasma interaction to astrophysical plasmas.
Charge exchange X-ray emission provides unique insights into the interactions between cold and hot astrophysical plasmas. Besides its own profound science, this emission is also technically crucial to all observations in the X-ray band, since charge exchange with the solar wind often contributes a significant foreground component that contaminates the signal of interest. By approximating the cross sections resolved to $n$ and $l$ atomic subshells, and carrying out complete radiative cascade calculation, we create a new spectral code to evaluate the charge exchange emission in the X-ray band. Comparing to collisional thermal emission, charge exchange radiation exhibits enhanced lines from large-$n$ shells to the ground, as well as large forbidden-to-resonance ratios of triplet transitions. Our new model successfully reproduces an observed high-quality spectrum of comet C/2000 WM1 (LINEAR), which emits purely by charge exchange between solar wind ions and cometary neutrals. It demonstrates that a proper charge exchange model will allow us to probe remotely the ion properties, including charge state, dynamics, and composition, at the interface between the cold and hot plasmas.
We perform two-dimensional particle-in-cell simulations of reconnection in magnetically dominated electron-positron plasmas subject to strong Compton cooling. We vary the magnetization $sigmagg1$, defined as the ratio of magnetic tension to plasma inertia, and the strength of cooling losses. Magnetic reconnection under such conditions can operate in magnetically dominated coronae around accreting black holes, which produce hard X-rays through Comptonization of seed soft photons. We find that the particle energy spectrum is dominated by a peak at mildly relativistic energies, which results from bulk motions of cooled plasmoids. The peak has a quasi-Maxwellian shape with an effective temperature of $sim 100$ keV, which depends only weakly on the flow magnetization and the strength of radiative cooling. The mean bulk energy of the reconnected plasma is roughly independent of $sigma$, whereas the variance is larger for higher magnetizations. The spectra also display a high-energy tail, which receives $sim 25$% of the dissipated reconnection power for $sigma=10$ and $sim 40$% for $sigma=40$. We complement our particle-in-cell studies with a Monte-Carlo simulation of the transfer of seed soft photons through the reconnection layer, and find the escaping X-ray spectrum. The simulation demonstrates that Comptonization is dominated by the bulk motions in the chain of Compton-cooled plasmoids and, for $sigmasim 10$, yields a spectrum consistent with the typical hard state of accreting black holes.
Log in to be able to interact and post comments
comments
Fetching comments
Sign in to be able to follow your search criteria
