Abstract:This work explores different particle-based approaches to the simulation of one-dimensional fractional subdiffusion equations in unbounded domains. We rely on smooth particle approximations, and consider four methods for estimating the fractional diffusion term. The first method is based on direct differentiation of the particle representation, it follows the Riesz definition of the fractional derivative and results in a non-conservative scheme. The other three methods follow the particle strength exchange (PSE) methodology and are by construction conservative, in the sense that the total particle strength is time invariant. The first PSE algorithm is based on using direct differentiation to estimate the fractional diffusion flux, and exploiting the resulting estimates in an integral representation of the divergence operator. Meanwhile, the second one relies on the regularized Riesz representation of the fractional diffusion term to derive a suitable interaction formula acting directly on the particle representation of the diffusing field. A third PSE construction is considered that exploits the Green's function of the fractional diffusion equation. The performance of all four approaches is assessed for the case of a one-dimensional diffusion equation with constant diffusivity. This enables us to take advantage of known analytical solutions, and consequently conduct a detailed analysis of the performance of the methods. This includes a quantitative study of the various sources of error, namely filtering, quadrature, domain truncation, and time integration, as well as a space and time self-convergence analysis. These analyses are conducted for different values of the order of the fractional derivatives, and computational experiences are used to gain insight that can be used for generalization of the present constructions.

A Practical Guide to Deterministic Particle Methods

On a practical implementation of particle methods

Particle Methods for Dispersive Equations

A Particle Method for the Semiclassical Limit of the Schrödinger Equation and the Vlasov--Poisson Equations

Particle Simulation of Fractional Diffusion Equations

Particle method for the numerical simulation of the path-dependent McKean-Vlasov equation

A particle method for non-local advection–selection–mutation equations

Discrete Particle Simulation of Particulate Systems: Theoretical Developments

Integration of the EPDiff Equation by Particle Methods

An efficient stochastic particle method for high-dimensional nonlinear PDEs

Implicit Temporal Discretization and Exact Energy Conservation for Particle Methods Applied to the Poisson–Boltzmann Equation

Efficient methods for particle-resolved direct numerical simulation

A particle method for non-local advection-selection-mutation equations

A Multilevel Monte Carlo Asymptotic-Preserving Particle Method for Kinetic Equations in the Diffusion Limit

Particle representations for stochastic partial differential equations with boundary conditions

Generalized Energy-Conserving Dissipative Particle Dynamics with Mass Transfer. Part 2: Applications and Demonstrations

Numerical Solution of Stochastic Differential Equations: Diffusion and Jump-Diffusion Processes

ADVANCES OF NUMERICAL METHODS FOR GENERAL DYNAMIC EQUATION DESCRIBING DISPERSED SYSTEM

Kinetic Equation for Particle Transport in Turbulent Flows

Generalised dissipative particle dynamics with energy conservation: density- and temperature-dependent potentials.

A practical guide to stochastic simulations of reaction-diffusion processes