Related papers: Novel Kinetic 3D MHD Algorithm for High Performanc…
Under the low density and high temperature conditions typical of heliospheric plasmas, the macroscopic evolution of the heliosphere is strongly affected by the kinetic plasma physics governing fundamental microphysical mechanisms. Kinetic…
An algorithm is proposed for finding numerical solutions of a kinetic equation that describes an infinite system of point articles placed in $\mathbb{R}^d (d \geq 1)$. The particles perform random jumps with pair wise repulsion, in the…
To simulate plasma phenomena, large-scale computational resources have been employed in developing high-precision and high-resolution plasma simulations. One of the main obstacles in plasma simulations is the requirement of computational…
We study numerical methods for dissipative particle dynamics (DPD), which is a system of stochastic differential equations and a popular stochastic momentum-conserving thermostat for simulating complex hydrodynamic behavior at mesoscales.…
We describe a novel splitting approach to numerical relativistic magnetohydrodynamics (RMHD) designed to expand its applicability to the domain of ultra-high magnetisation (high-$\sigma$). In this approach, the electromagnetic field is…
The paper describes an explicit multi-dimensional numerical scheme for Special Relativistic Two-Fluid Magnetohydrodynamics of electron-positron plasma and a suit of test problems. The scheme utilizes Cartesian grid and the third order WENO…
Astrophysical plasmas in relativistic spacetimes, such as black hole accretion flows, are often weakly collisional and require kinetic modeling to capture non-local transport and particle acceleration. However, the extreme scale separation…
Space plasma simulations have seen an increase in the use of magnetohydrodynamic (MHD) with embedded Particle-in-Cell (PIC) models. This combined MHD-EPIC algorithm simulates some regions of interest using the kinetic PIC method while…
Most plasmas are only partially ionized. To better understand the dynamics of these plasmas, the behaviors of a mixture of neutral species and plasma in ideal magnetohydrodynamic states are investigated. The current approach is about the…
We propose a hybrid moment method for the multi-scale kinetic equations in the framework of the regularized moment method [7]. In this method, the fourth order moment system is chosen as the governing equations in the fluid region. When…
Symplectic schemes are powerful methods for numerically integrating Hamiltonian systems, and their long-term accuracy and fidelity have been proved both theoretically and numerically. However direct applications of standard symplectic…
We present novel numerical schemes for ideal magnetohydrodynamic (MHD) simulations aimed at enhancing stability against numerical shock instability and improving the accuracy of low-speed flows in multidimensions. Stringent benchmark tests…
We present an algorithm specifically tailored for solving kinetic equations onto GPUs. The efficiency of the algorithm is demonstrated by solving the one-dimensional shock wave structure problem and a two-dimensional low Mach number driven…
Numerical methods for solving the ideal magnetohydrodynamic (MHD) equations in more than one space dimension must either confront the challenge of controlling errors in the discrete divergence of the magnetic field, or else be faced with…
Numerical simulations of turbulent flows are well known to pose extreme computational challenges due to the huge number of dynamical degrees of freedom required to correctly describe the complex multi-scale statistical correlations of the…
This article proposes a new statistical numerical method to address gas kinetics problems obeying the Boltzmann equation. This method is inspired from some Monte-Carlo algorithms used in linear transport physics, where virtual particles are…
Studying the propagation of uncertainties in a nonlinear dynamical system usually involves generating a set of samples in the stochastic parameter space and then repeated simulations with different sampled parameters. The main difficulty…
We present a new algorithm which is named the Dynamical Functional Particle Method, DFPM. It is based on the idea of formulating a finite dimensional damped dynamical system whose stationary points are the solution to the original…
For quantum systems with competing potentials, the conventional perturbation theory often yields an asymptotic series and the subsequent numerical outcome becomes uncertain. To tackle such kind of problems, we develop a general solution…
Quasiparticle dynamics in relativistic plasmas associated with hot, weakly-coupled gauge theories (such as QCD at asymptotically high temperature $T$) can be described by an effective kinetic theory, valid on sufficiently large time and…