Related papers: Multiple Boris integrators for particle-in-cell si…
We introduce a particle-based simulation method for granular material in interactive frame rates. We divide the simulation into two decoupled steps. In the first step, a relatively small number of particles is accurately simulated with a…
Modeling multi-scale collisionless magnetized processes constitutes an important numerical challenge. By treating electrons as a fluid and ions kinetically, the so-called hybrid Particle-In-Cell (PIC) codes represent a promising…
The identification of pre-acceleration mechanisms for cosmic ray ions in supernova remnant shocks is an important problem in astrophysics. Recent particle-in-cell (PIC) shock simulations have shown that inclusion of the full electron…
A basic leapfrog integrator and its energy-preserving and variational / symplectic variants are proposed and studied for the numerical integration of the equations of motion of relativistic charged particles in an electromagnetic field. The…
In this article we introduce a novel coupled algorithm for massively parallel direct numerical simulations of electrophoresis in microfluidic flows. This multiphysics algorithm employs an Eulerian description of fluid and ions, combined…
Particulate flows have been largely studied under the simplifying assumptions of one-way coupling regime where the disperse phase do not react-back on the carrier fluid. In the context of turbulent flows, many non trivial phenomena such as…
We provide an algorithm for evolving general spin-$s$ Gross-Pitaevskii / non-linear Schr\"odinger systems carrying a variety of interactions, where the $2s+1$ components of the `spinor' field represent the different spin-multiplicity…
Use of explicit integration methods for power electronic circuits with ideal switch models significantly improves simulation speed. The PLECS package [1] has effectively used this idea; however, the implementation details involved in PLECS…
Symplectic N-body integrators are widely used to study problems in celestial mechanics. The most popular algorithms are of 2nd and 4th order, requiring 2 and 6 substeps per timestep, respectively. The number of substeps increases rapidly…
Exposing a molecule to intense light pulses may bring this molecule to a nonstationary quantum state, thus launching correlated dynamics of electronic and nuclear subsystems. Although much had been achieved in the understanding of…
We present several implementations of the Metropolis method, an adaptive Monte Carlo algorithm, which allow for the calculation of multi-dimensional integrals over arbitrary on-shell four-momentum phase space. The Metropolis technique…
In a recent paper Lima, Panario and Wang have provided a new method to multiply polynomials in Chebyshev basis which aims at reducing the total number of multiplication when polynomials have small degree. Their idea is to use Karatsuba's…
Computer simulation methods, such as Monte Carlo or Molecular Dynamics, are very powerful computational techniques that provide detailed and essentially exact information on classical many-body problems. With the advent of ab-initio…
This paper presents a novel 3-D full electromagnetic particle-in-cell (PIC) code called JefiPIC, which uses Jefimenko's equations as the electromagnetic (EM) field solver through a full-space integration method. Leveraging the power of…
We present particle-in-cell simulations of one dimensional relativistic electromagnetic shocks in a uniform magnetic field, for a range of magnetic field strengths, plasma temperatures and numerical initial conditions. We show that the…
Particle-in-cell with Monte Carlo collisions (PIC/MCC) is a fully kinetic, particle based numerical simulation method with increasing popularity in the field of low temperature gas discharge physics. Already in its simplest form…
The accurate prediction of occurrence and strength of kinetic instabilities in plasmas remains a significant challenge in nuclear fusion research. To accurately capture the plasmas dynamics one is required to solve the Vlasov equation for…
Integrability conditions for systems of bosons or fermions with seniority conserving hamiltonians are derived. The conditions are shown to be invariant under a large class of transformations of the interaction matrix elements. Previously…
We construct numerical integrators for Hamiltonian problems that may advantageously replace the standard Verlet time-stepper within Hybrid Monte Carlo and related simulations. Past attempts have often aimed at boosting the order of accuracy…
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…