Related papers: The Energy Conserving Particle-in-Cell Method
In this work, we develop energy-preserving iterative schemes for the (non-)linear systems arising in the Gauss integration of Poisson systems with quadratic Hamiltonian. Exploiting the relation between Gauss collocation integrators and…
This paper is concerned with the numerical integration in time of nonlinear Schr\"odinger equations using different methods preserving the energy or a discrete analog of it. The Crank-Nicolson method is a well known method of order 2 but is…
In this paper, it is shown that three-dimensional stochastic Maxwell equations with multiplicative noise are stochastic Hamiltonian partial differential equations possessing a geometric structure (i.e. stochastic mutli-symplectic…
Until recently, electromagnetic finite element PIC (EM-FEMPIC) methods that demonstrated charge conservation used explicit field solvers. It is only recently, that a series of papers developed the mathematics necessary for charge…
A simple form of the Boris solver in particle-in-cell (PIC) simulation is proposed. It employs an exact solution of the Lorentz-force part, and it is equivalent to the Boris solver with a gyrophase correction. As a favorable property for…
The study of the long time conservation for numerical methods poses interesting and challenging questions from the point of view of geometric integration. In this paper, we analyze the long time energy and magnetic moment conservations of…
The Unified Gas-Kinetic Wave-Particle (UGKWP) method is constructed for partially ionized plasma (PIP). This method possesses both multiscale and unified preserving (UP) properties. The multiscale property allows the method to capture a…
There are many interesting physical processes which involve the generation of high density plasmas in large volumes. However, when modeling these systems numerically, the large densities and volumes present a significant computational…
A new scheme for numerical integration of the 1D2V relativistic Vlasov-Maxwell system is proposed. Assuming that all particles in a cell of the phase space move with the same velocity as that of the particle located at the center of the…
We present a numerical model and a set of conservative algorithms for Non-Maxwellian plasma kinetics with inelastic collisions. These algorithms self-consistently solve for the time evolution of an isotropic electron energy distribution…
In this article we describe a unifying framework for variational electromagnetic particle schemes of spectral type, and we propose a novel spectral Particle-In-Cell (PIC) scheme that preserves a discrete Hamiltonian structure. Our work is…
We extend the asymptotic preserving and energy conserving time integrator for charged-particle motion developed in [Ricketson & Chac\'on, JCP, 2020] to include finite Larmor-radius (FLR) effects in the presence of electric-field…
In this paper, we present three-dimensional (3D) Particle-In-Cell (PIC) simulations to study the stability of 2D Bernstein-Greene-Kruskal (BGK) modes in a magnetized plasma with a finite background magnetic field. The simulations were…
We discuss the interpolation of the electric and magnetic fields within a charge-conserving Particle-In-Cell scheme. The choice of the interpolation procedure for the fields acting on a particle can be constrained by analyzing conservation…
In this paper, we develop a framework to construct energy-preserving methods for multi-components Hamiltonian systems, combining the exponential integrator and the partitioned averaged vector field method. This leads to numerical schemes…
The dynamics of two initially unmagnetized relativistic counter-streaming homogeneous ion-electron plasma beams are simulated in two dimensions using the particle-in-cell (PIC) method. It is shown that current filaments, which form due to…
The Particle-In-Cell (PIC) method has been developed by Oscar Buneman, Charles Birdsall, Roger W. Hockney, and John Dawson in the 1950s and, with the advances of computing power, has been further developed for several fields such as…
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…
The shell correction method is revisited. Contrary to the traditional Strutinsky method, the shell energy is evaluated by an averaging over the number of particles and not over the single-particle energies, which is more consistent with the…
The rapid development of AI highlights the pressing need for sustainable energy, a critical global challenge for decades. Nuclear fusion, generally seen as an ultimate solution, has been the focus of intensive research for nearly a century,…