Related papers: Code-Verification Techniques for Particle-in-Cell …
A second-order accurate divergence-conserving hybrid particle-in-cell code Maximus has been developed for microscopic modeling of collisionless plasmas. The main specifics of the code include a constrained transport algorithm for exact…
A number of problems arise when long-range forces, such as those governed by Bessel functions, are used in particle-particle simulations. If a simple cut-off for the interaction is used, the system may find an equilibrium configuration at…
We analyze here in some detail, the derivation of the Particle and Monte Carlo methods of plasma simulation, such as Particle in Cell (PIC), Monte Carlo (MC) and Particle in Cell / Monte Carlo (PIC/MC) from formal manipulation of transport…
In this paper we develop a direct simulation Monte Carlo (DSMC) method for simulating highly nonequilibrium dynamics of nearly degenerate ultra-cold gases. We show that our method can simulate the high-energy collision of two thermal clouds…
Molecular Dynamics simulations can help scientists to gather valuable insights for physical processes on an atomic scale. This work explores various techniques for SIMD vectorization to improve the pairwise force calculation between…
The present study addresses the challenge of enhancing computational efficiency without compromising accuracy in numerical simulations of vacuum gas dynamics using the direct simulation Monte Carlo (DSMC) method. A technique termed "fixed…
We report an accessible and robust tool for evaluating the effects of Coulomb collisions on a test particle in a plasma that obeys Maxwell-J\"uttner statistics. The implementation is based on the Beliaev-Budker collision integral which…
Computer simulation is an important tool for scientific progress, especially when lab experiments are either extremely costly and difficult or lack the required resolution. However, all of the simulation methods come with limitations. In…
Rigid bodies, made of smaller composite beads, are commonly used to simulate anisotropic particles with molecular dynamics or Monte Carlo methods. To accurately represent the particle shape and to obtain smooth and realistic effective pair…
We formulate and test a hybrid fluid-Monte Carlo scheme for the treatment of elastic collisions in gases and plasmas. While our primary focus and demonstrations of applicability are for moderately collisional plasmas, as described by the…
A simulation package employing a Particle in Cell (PIC) method is developed to study the high current beam transport and the dynamics of plasmas. This package includes subroutines those are suited for various planned projects at University…
In the present work, we present a novel numerical algorithm to couple the Direct Simulation Monte Carlo method (DSMC) for the solution of the Boltzmann equation with a finite volume like method for the solution of the Euler equations.…
Large-scale plasma simulations are critical for designing and developing next-generation fusion energy devices and modeling industrial plasmas. BIT1 is a massively parallel Particle-in-Cell code designed for specifically studying plasma…
A treatment of direct simulation Monte Carlo method (DSMC) as a Markov process with a master equation is given and the corresponding master equation is derived. A hierarchy of equations for the reduced probability distributions is derived…
We present the standard electromagnetic Particle-in-Cell method, starting from the discrete approximation of derivatives on a uniform grid. The application to second-order, centered, finite-difference discretization of the equations of…
Simulation studies are used to evaluate and compare the properties of statistical methods in controlled experimental settings. In most cases, performing a simulation study requires knowledge of the true value of the parameter, or estimand,…
In a world made of atoms, the computer simulation of molecular systems, such as proteins in water, plays an enormous role in science. Software packages that perform these computations have been developed for decades. In molecular…
We present a new, for plasma physics, highly efficient multilevel Monte Carlo numerical method for simulating Coulomb collisions. The method separates and optimally minimizes the finite-timestep and finite-sampling errors inherent in the…
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…
While boundary plasmas in present-day tokamaks generally fall in a fluid regime, neutral species near the boundary often require kinetic models due to long mean-free-paths compared to characteristic spatial scales in the region. Monte-Carlo…