Related papers: High-order exponential solver method for particle-…
Magnetic quadrupoles are essential components of particle accelerators like the Large Hadron Collider. In order to study numerically the stability of the particle beam crossing a quadrupole, a large number of particle revolutions in the…
We introduce $\texttt{A}$strophysical $\texttt{H}$ybrid-$\texttt{K}$inetic simulations with the $\texttt{flASH}$ code ($\texttt{AHKASH}$) -- a new Hybrid particle-in-cell (PIC) code developed within the framework of the multi-physics code…
The paper is concerned with overlapping domain decomposition and exponential time differencing for the diffusion equation discretized in space by cell-centered finite differences. Two localized exponential time differencing methods are…
Energy conserving particle-in-cell schemes are constructed for a class of reduced relativistic Vlasov--Maxwell equations of laser-plasma interaction. Discrete Poisson equation is also satisfied by the numerical solution. Specifically,…
The effect of electron temperature on the space-time evolution of nonlinear plasma oscillations in an inhomogeneous plasma is studied using a one-dimensional particle-in-cell (PIC) code. In contrast to the conventional wisdom, it is found…
Particle acceleration in collisionless plasma systems is a central question in astroplasma and astroparticle physics. The structure of the acceleration regions, electron-ion energy equilibration, preacceleration of particles at shocks to…
Time series alignment methods call for highly expressive, differentiable and invertible warping functions which preserve temporal topology, i.e diffeomorphisms. Diffeomorphic warping functions can be generated from the integration of…
Over the past decades, multiple gyrokinetic codes have shown to be able to simulate turbulence and associated transport in the core of Tokamak devices. However, their application to the edge and scrape-off layer (SOL) region presents…
In this work, we suggest an easy-to-code higher-order finite volume semi-discrete scheme to analyze the nonlinear behavior of the electron-plasma oscillations by solving electron fluid equations numerically. The present method employs a…
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…
We present a computationally efficient model that describes the interaction of intense, ultrashort infrared laser pulses with transparent materials in the strong ionization regime. The model is augmented with a detailed self-consistent…
Particle-in-Cell (PIC) simulations are fundamental to plasma physics but often suffer from limited scalability due to particle-grid interaction bottlenecks and particle redistribution costs. Specifically, the particle-grid interaction…
While the particle-in-cell (PIC) method is quite mature, verification and validation of both newly developed methods and individual codes has largely focused on an idiosyncratic choice of a few test cases. Many of these test cases involve…
In this study, we consider a numerical implementation of the nonlinear Rosenbluth-Trubnikov collision operator for particle simulations in plasma physics in the framework of the finite element method (FEM). The relevant particle evolution…
We present a detailed strong and weak scaling analysis of PICCANTE, an open source, massively parallel, fully-relativistic Particle-In-Cell (PIC) code. PIC codes are widely used in plasma physics and astrophysics to study the cases where…
In physics research particle accelerators are highly valued, and extraordinarily expensive, technical instruments. The high cost of particle accelerators results from the immense lengths required to accelerate particles to high energies,…
We present novel model reduction methods for rapid solution of parametrized nonlinear partial differential equations (PDEs) in real-time or many-query contexts. Our approach combines reduced basis (RB) space for rapidly convergent…
A high-order convergent numerical method for solving linear and non-linear parabolic PDEs is presented. The time-stepping is done via an explicit, singly diagonally implicit Runge-Kutta (ESDIRK) method of order 4 or 5, and for the implicit…
We present a higher-order boundary condition for atomistic simulations of dislocations that address the slow convergence of standard supercell methods. The method is based on a multipole expansion of the equilibrium displacement, combining…
Particle-in-cell methods couple mesh-based methods for the solution of continuum mechanics problems, with the ability to advect and evolve particles. They have a long history and many applications in scientific computing. However, they have…