Related papers: The Magnetohydrodynamic-Particle-In-Cell Module in…
Magnetic fields play an important role in astrophysics on a wide variety of scales, ranging from the Sun and compact objects to galaxies and galaxy clusters. Here we discuss a novel implementation of ideal magnetohydrodynamics (MHD) in the…
The Pencil Code is a highly modular physics-oriented simulation code that can be adapted to a wide range of applications. It is primarily designed to solve partial differential equations (PDEs) of compressible hydrodynamics and has lots of…
Context. Relativistic jets in active galactic nuclei are known for their exceptional energy output, and imaging the synthetic synchrotron emission of numerical jet simulations is essential for a comparison with observed jet polarization…
Recent developments in compact object astrophysics, especially the discovery of merging neutron stars by LIGO, the imaging of the black hole in M87 by the Event Horizon Telescope (EHT) and high precision astrometry of the Galactic Center at…
In this era of precision cosmology, a detailed physical understanding on the evolution of cosmic baryons is required. Cosmic magnetic fields, though still poorly understood, may represent an important component in the global cosmic energy…
Recent observations in X-ray and Gamma-ray suggest that the emission region of the pulsar magnetosphere can be multifold. In particular, the open-close boundary of the magnetic field, so-called the Y-point, can be the new candidate place…
This paper introduces a novel formulation of the Particle-In-Cell (PIC) method for the modeling of relativistic plasmas, which leverages the ability of the Pseudo-Spectral Analytical Time-Domain solver (PSATD) to handle arbitrary time…
This paper introduces the Sheffield Magnetohydrodynamics Algorithm Using GPUs (SMAUG+), an advanced numerical code for solving magnetohydrodynamic (MHD) problems, using multi-GPU systems. Multi-GPU systems facilitate the development of…
Ion cyclotron resonance energization (ICRE) such as ion cyclotron resonance heating (ICRH) is widely applied to magnetic confinement fusion and high-power electric propulsion. Since ICRE involves cyclotron resonance processes, a kinetic…
A hybrid spectral/finite-element code is developed to numerically solve the resistive finite-pressure magnetohydrodynamic equilibria without the necessity of postulating nested magnetic flux surfaces in the non-axisymmetric toroidal…
We present an implementation of Smoothed Particle Hydrodynamics (SPH) in an adaptive-mesh PPPM algorithm. The code evolves a mixture of purely gravitational particles and gas particles. The code retains the desirable properties of previous…
The particle-in-cell approach has proven effective at modeling neutron star and black hole magnetospheres from first principles, but global simulations are plagued with an unrealistically small separation between the scales where…
We present the first plasma simulations obtained with the code dHybridR, a hybrid particle-in-cell code with fluid electrons and both thermal and energetic ions that retain relativistic dynamics. dHybridR is constructed to study…
Solving the problem of fast eruptive events in magnetically dominated astrophysical plasmas requires the use of particularly well adapted numerical tools. Indeed, the central mechanism based on magnetic reconnection is determined by a…
We present a new code designed to solve the equations of classical ideal magneto-hydrodynamics (MHD) in three dimensions, submitted to a constant gravitational field. The purpose of the code centers on the analysis of solar phenomena within…
Our fully kinetic, implicit Particle-in-Cell (PIC) simulations of global magnetospheres on up to 32,768 of El Capitan's AMD Instinct MI300A Accelerated Processing Units (APUs) represent an unprecedented computational capability that…
Relativistic macroscopic plasma dynamics can be described by general-relativistic magnetohydrodynamics. In many high-energy astrophysical settings, such as the interior dynamics of magnetized stars, the ideal GRMHD approximation, in which…
For simulating laser-plasma interactions, we developed a parallel, multi-dimensional, fully relativistically particle-in-cell (PIC) code, named YUNIC. The core algorithm is introduced, including field solver, particle pusher, field…
This work describes a new 1D hybrid approach for modeling atmospheric pressure discharges featuring complex chemistry. In this approach electrons are described fully kinetically using Particle-In-Cell/Monte-Carlo (PIC/MCC) scheme, whereas…
Particle-in-Cell (PIC) simulation codes have wide applicability to first-principles modeling of multidimensional nonlinear plasma phenomena, including wake-field accelerators. This review addresses both finite difference and pseudo-spectral…