Related papers: A second-order particle Fokker-Planck-Master metho…
We describe a hybrid Direct Simulation Monte Carlo (DSMC) code for simultaneously solving the collisional Boltzmann equation for gas and the collisionless Boltzmann equation for stars and dark matter for problems important to galaxy…
Models and methods that are able to accurately and efficiently predict the flows of low-speed rarefied gases are in high demand, due to the increasing ability to manufacture devices at micro and nano scales. One such model and method is a…
A high order, deterministic direct numerical method is proposed for the nonrelativistic $2D_{\bf x} \times 3D_{\bf v}$ Vlasov-Maxwell system, coupled with Fokker-Planck-Landau type operators. Such a system is devoted to the modelling of…
In this paper, we propose and analyze a second order accurate (in both time and space) numerical scheme for the Poisson-Nernst-Planck-Navier-Stokes system, which describes the ion electro-diffusion in fluids. In particular, the…
Solving the Fokker-Planck equation for high-dimensional complex turbulent dynamical systems is an important and practical issue. However, most traditional methods suffer from the curse of dimensionality and have difficulties in capturing…
Second-order Volume-preserving algorithms (VPAs) for simulating charged particle motion in electromagnetic fields have been generalized to a rotating angle formulation by using the matrix decomposition methods. Based on this method, the…
A double-distribution-function based lattice Boltzmann method (DDF-LBM) is proposed for the simulation of polyatomic gases in the supersonic regime. The model relies on an extended equilibrium state that is constructed to reproduce the…
Particle flow filters solve Bayesian inference problems by smoothly transforming a set of particles into samples from the posterior distribution. Particles move in state space under the flow of an McKean-Vlasov-Ito process. This work…
The unified gas-kinetic scheme (UGKS) provides a framework for simulating multiscale transport with the updates of both gas distribution function and macroscopic flow variables on the cell size and time step scales. The multiscale dynamics…
Photomultiplier tubes (PMTs) are widely deployed at neutrino and dark matter experiments for photon counting. When multiple photons hit a PMT consecutively, their photo-electron (PE) pulses pile up to hinder the precise measurements of the…
We propose a method to efficiently determine the optimal coarse-grained force field in mesoscopic stochastic simulations of Newtonian fluid and polymer melt systems modeled by dissipative particle dynamics (DPD) and energy conserving…
This paper proposes and validates two new particle regularization techniques for the Smoothed Particle Hydrodynamics (SPH) numerical method to improve its stability and accuracy for free surface flow simulations. We introduce a general form…
We present a new formulation of the central moment lattice Boltzmann (LB) method based on a continuous Fokker-Planck (FP) kinetic model, originally proposed for stochastic diffusive-drift processes (e.g., Brownian dynamics), by adapting it…
We propose an explicit particle method for the Vlasov-Fokker-Planck equation that conserves energy at the fully discrete level. The method features two key components: a deterministic and conservative particle discretization for the…
Validity of fluid models breaks down for non-thermal or weakly collisional plasmas which often occur e.g. in the solar wind. In these regimes one has to resort to modelling through the first-principle Vlasov-Maxwell system, but its…
The Direct Simulation Monte Carlo (DSMC) method was widely used to simulate low density gas flows with large Knudsen numbers. However, DSMC encounters limitations in the regime of lower Knudsen numbers (Kn<0.1). In such cases, approaches…
Uniform error estimates of a bi-fidelity method for a kinetic-fluid coupled model with random initial inputs in the fine particle regime are proved in this paper. Such a model is a system coupling the incompressible Navier-Stokes equations…
Coulomb collisions in particle simulations for weakly coupled plasmas are modeled by the Landau-Fokker-Planck equation, which is typically solved by Monte-Carlo (MC) methods. One of the main disadvantages of MC is the timestep accuracy…
Stochastic reaction-diffusion models are now a popular tool for studying physical systems in which both the explicit diffusion of molecules and noise in the chemical reaction process play important roles. The Smoluchowski diffusion-limited…
Improving the efficiency of the direct simulation Monte Carlo (DSMC) method has become increasingly urgent with the rapid development of space exploration. To address this issue, the direct intermittent general synthetic iteration (DIG)…