Related papers: The Energy Conserving Particle-in-Cell Method
We discuss a spectral method for the numerical solution of the Vlasov-Poisson system where the velocity space is decomposed by means of an Hermite basis. We describe a semi-implicit time discretization that extends the range of numerical…
We take issue with the conclusions in the recent publication by Savard et al. In the study, the authors implement a fully nonlinear charge- and energy-conserving implicit particle-in-cell method (ECC-IPIC), and use it to study the impact of…
A reduced-speed-of-light (RSOL) approximation is a useful technique for magnetohydrodynamic (MHD)-particle-in-cell (PIC) simulations. With a RSOL, some "in-code" speed-of-light $\tilde{c}$ is set to much lower values than the true $c$,…
We propose a new method to conserve the total energy to round-off error in grid-based codes for hydrodynamic simulations with self-gravity. A formula for the energy flux due to the work done by the the self-gravitational force is given, so…
Understanding electron energy dynamics in low-temperature plasmas such as capacitively coupled plasmas (CCPs), including energy absorption, conversion, transport, and dissipation, is essential for interpreting discharge physics and process…
A finite element discretization using a method of lines approached is proposed for approximately solving the Poisson-Nernst-Planck (PNP) equations. This discretization scheme enforces positivity of the computed solutions, corresponding to…
Discontinuous Galerkin methods are developed for solving the Vlasov-Maxwell system, methods that are designed to be systematically as accurate as one wants with provable conservation of mass and possibly total energy. Such properties in…
We introduce an extensible multi-fluid moment model in the context of collisionless magnetic reconnection. This model evolves full Maxwell equations, and simultaneously moments of the Vlasov-Maxwell equation for each species in the plasma.…
This article focuses on an energy-conservation Galerkin finite element method (FEM) for the generalized Klein-Gordon-Zakharov (KGZ) equations. This method combines the bilinear finite element method for spatial discretization with the…
In this work, we develop a new compatible finite element formulation of the thermal shallow water equations that conserves energy and mathematical entropies given by buoyancy-related quadratic tracer variances. Our approach relies on…
Using a 2-1/2-dimensional particle-in-cell (PIC) code to simulate the relativistic expansion of a magnetized collisionless plasma into a vacuum, we report a new mechanism in which the magnetic energy is efficiently converted into the…
We develop an inertial coupling method for modeling the dynamics of point-like 'blob' particles immersed in an incompressible fluid, generalizing previous work for compressible fluids. The coupling consistently includes excess (positive or…
Particle in Cell (PIC) simulations are a widely used tool for the investigation of both laser- and beam-driven plasma acceleration. It is a known issue that the beam quality can be artificially degraded by numerical Cherenkov radiation…
In this research, we introduce and investigate an approximation method that preserves the structural integrity of the non-isothermal Cahn-Hilliard-Navier-Stokes system. Our approach extends a previously proposed technique [1], which…
In this paper, we propose energy-conserving numerical schemes for the Vlasov-Amp\`{e}re (VA) systems. The VA system is a model used to describe the evolution of probability density function of charged particles under self consistent…
The effects of noise in particle-in-cell (PIC) and Vlasov simulations of the Buneman instability in unmagnetized plasmas are studied. It is found that, in the regime of low drift velocity, the linear stage of the instability in PIC…
A new methodology is developed to integrate numerically the equations of motion for classical many-body systems in molecular dynamics simulations. Its distinguishable feature is the possibility to preserve, independently on the size of the…
This paper investigates subcycling of particle orbits in variational, geometric particle-in-cell methods addressing the Vlasov--Maxwell system in magnetized plasmas. The purpose of subcycling is to allow different time steps for different…
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…
The acceleration of protons and electrons in a reconnecting current sheet (RCS) is simulated with a particle-in-cell (PIC) 2D3V code for the proton-to-electron mass ratio of 100. The electro-magnetic configuration forming the RCS…