Related papers: Direct simulation Monte Carlo schemes for Coulomb …
We simulate TIP3P water using a constrained Monte Carlo algorithm to generate electrostatic interactions eliminating the need to sum over long ranged Coulomb interactions. We study discretization errors when interpolating charges using…
We investigate the stochastic Direct Simulation Monte Carlo method (DSMC) for numerically solving the collision-term in heavy-ion transport theories of the Boltzmann-Uehling-Uhlenbeck (BUU) type. The first major modification we consider is…
In [C. Mouhot and L. Pareschi, "Fast algorithms for computing the Boltzmann collision operator," Math. Comp., to appear; C. Mouhot and L. Pareschi, C. R. Math. Acad. Sci. Paris, 339 (2004), pp. 71-76], fast deterministic algorithms based on…
Large-angle Coulomb collisions lead to an avalanching generation of runaway electrons in a plasma. We present the first fully conservative large-angle collision operator, derived from the relativistic Boltzmann operator. The relation to…
We analyze and compare the computational complexity of different simulation strategies for Monte Carlo in the setting of classically scaled population processes. This allows a range of widely used competing strategies to be judged…
This work explores the capability of simulating complex fluid flows by directly solving the Boltzmann equation. Due to the high-dimensionality of the governing equation, the substantial computational cost of solving the Boltzmann equation…
A jump-diffusion process along with a particle scheme is devised as an accurate and efficient particle solution to the Boltzmann equation. The proposed process (hereafter Gamma-Boltzmann model) is devised to match the evolution of all…
This is a review of recent developments in Monte Carlo methods in the field of ultra cold gases. For bosonic atoms in an optical lattice we discuss path integral Monte Carlo simulations with worm updates and show the excellent agreement…
We present a Direct Simulation Monte Carlo (DSMC) method for the spatially homogeneous multispecies Landau-Fokker-Planck equation. The scheme is derived from a first-order approximation of the multispecies Boltzmann operator in the grazing…
An accurate algorithm is proposed to improve the prediction of a particle in collision with a moving wall within the direct simulation Monte Carlo (DSMC) framework for the simulation of unsteady rarefied flows. This algorithm is able to…
We study discrete-time simulation schemes for stochastic Volterra equations, namely the Euler and Milstein schemes, and the corresponding Multi-Level Monte-Carlo method. By using and adapting some results from Zhang [22], together with the…
We combine ab initio path integral Monte Carlo (PIMC) simulations with fixed ion configurations from density functional theory molecular dynamics (DFT-MD) simulations to solve the electronic problem for hydrogen under warm dense matter…
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…
We present a combined theoretical and experimental study of primary and post-collision mechanisms involved when colliding low energy multiply charged ions with van der Waals dimers. The collision dynamics is investigated using a classical…
We develop a novel Monte Carlo strategy for the simulation of the Boltzmann-BGK model with both low-collisional and high-collisional regimes present. The presented solution to maintain accuracy in low-collisional regimes and remove…
We present the formulation of a kinetic mapping scheme between the Direct Simulation Monte Carlo (DSMC) and the Lattice Boltzmann Method (LBM) which is at the basis of the hybrid model used to couple the two methods in view of efficiently…
The linearized model collision operator for multiple species plasmas given by H. Sugama, T.-H. Watanabe, and M. Nunami [Phys.\ Plasmas {\bf 16}, 112503 (2009)] is improved to be properly applicable up to the highly collisional regime. The…
The construction of discrete velocity models or numerical methods for the Boltzmann equation, may lead to the necessity of computing the collision operator as a sum over lattice points. The collision operator involves an integral over a…
We propose and investigate a new multi-level Monte Carlo scheme for numerical solutions of the kinetic Boltzmann equation for neutral species in edge plasmas. In particular, this method explicitly exploits a key structural property of…
Recent progress in treating the dynamically screened nature of the Coulomb interaction in strongly correlated lattice models and materials is reviewed with a focus on computational schemes based on the dynamical mean field approximation. We…