Related papers: Multiple Boris integrators for particle-in-cell si…
We extend the recently developed entropic and conservative variance reduction framework [M. Sadr, N. G. Hadjiconstantinou, A variance-reduced direct Monte Carlo simulation method for solving the Boltzmann equation over a wide range of…
An expansion procedure using third kind Chebyshev polynomials as base functions is suggested for solving second type Volterra integral equations with logarithmic kernels. The algorithm's convergence is studied and some illustrative examples…
The electromagnetic behavior of plasmonic structures can be predicted after discretizing and solving a linear system of equations, derived from a continuous surface integral equation (SIE) and the appropriate boundary conditions, using a…
In low-temperature high-density plasmas quantum effects of the electrons are becoming increasingly important. This requires the development of new theoretical and computational tools. Quantum Monte Carlo methods are among the most…
We show how to improve the molecular dynamics step of Hybrid Monte Carlo, both by tuning the integrator using Poisson brackets measurements and by the use of force gradient integrators. We present results for moderate lattice sizes.
We propose new local error estimators for splitting and composition methods. They are based on the construction of lower order schemes obtained at each step as a linear combination of the intermediate stages of the integrator, so that the…
A method of solving the Schr\"{o}dinger equation based on the use of constant particle-particle interaction potential surfaces (IPS) is proposed. The many-body wave function is presented in a configuration interaction form, with…
Self-consistent multi-particle simulation plays an important role in studying beam-beam effects and space charge effects in high-intensity beams. The Poisson equation has to be solved at each time-step based on the particle density…
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…
The particle-in-cell (PIC) method is successfully used to study magnetized plasmas. However, this requires large computational costs and limits simulations to short physical run-times and often to setups in less than three spatial…
In this paper, energy-preserving methods are formulated and studied for solving charged-particle dynamics. We first formulate the scheme of energy-preserving methods and analyze its basic properties including algebraic order and symmetry.…
Immersed boundary-lattice Boltzmann method (IB-LBM) has been widely used for simulation of particle-laden flows recently. However, it was limited to small-scale simulations with no more than O(103) particles. Here, we expand IB-LBM for…
We present a collocated-grid framework for Direct Numerical Simulations of polydisperse particles submerged in a viscous fluid. The fluid-particle forces are coupled with the Immersed Boundary Method (IBM) while the particle-particle forces…
First-principles particle-in-cell (PIC) simulation is a powerful tool for understanding plasma behavior, but this power often comes at great computational expense. Artificially reducing the ion/electron mass ratio is a time-honored practice…
In this article, we introduce a novel dimensionality reduction formulation for the Poisson's equation in the Vlasov-Poisson system that yields a reduced-order particle-in-cell scheme. This scheme allows a remarkable reduction in the…
In particle-in-cell simulations, excessive or even unfeasible computational demands can be caused by the growth of the number of particles in the course of prolific ionization or cascaded pair production due to the effects of quantum…
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…
Trapped Bosons exhibit fundamental physical phenomena and are potentially useful for quantum technologies. We present a method for simulating Bosons using path integral molecular dynamics. A main challenge for simulations is including all…
In this paper, we extend several time reversible numerical integrators to solve the Lorentz force equations from second order accuracy to higher order accuracy for relativistic charged particle tracking in electromagnetic fields. A fourth…
We propose a novel time stepping method for linear poroelasticity by extending a recent iterative decoupling approach to the second-order case. This results in a two-step scheme with an inner iteration and a relaxation step. We prove…