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 "Consistent" Quasidiffusion Method for Solving Particle Transport Problems

Stabilizing the Consistent Quasidiffusion Method with Linear Prolongation

On the Accuracy of Transport Spatial Discretization Methods for Diffusive Regimes

An improved method for solving quasilinear convection diffusion problems on a coarse mesh

Qualitatively accurate spectral schemes for advection and transport

The Extended Diamond Difference - Constant Nodal Method with Decoupled Cell Iteration Scheme in Two-Dimensional Discrete Ordinate Transport Problems

Numerical Study of Lattice Boltzmann Methods for a Convection–diffusion Equation Coupled with Navier–Stokes Equations

A high order conservative semi-Lagrangian discontinuous Galerkin method for two-dimensional transport simulations

A Characteristics-Mixed Finite Element Method for Advection-Dominated Transport Problems

A consistency-driven particle-advection formulation for weakly-compressible smoothed particle hydrodynamics

Consistent Second Moment Methods with Scalable Linear Solvers for Radiation Transport

A Discontinuous Galerkin Finite Element Method Based Axial SN for the 2D/1D Transport Method

Particle Simulation of Fractional Diffusion Equations

A high-order generalized differential quadrature method with lattice Boltzmann flux solver for simulating incompressible flows

A finite volume method preserving the invariant region property for the quasimonotone reaction-diffusion systems

Efficient finite element method for 2D singularly perturbed parabolic convection diffusion problems with discontinuous source term

Upscaling transport in heterogeneous media featuring local-scale dispersion: flow channeling, macro-retardation and parameter prediction

A conservative Eulerian finite element method for transport and diffusion in moving domains

Delta-form-based method of solving high order spatial discretization schemes for neutron transport

A Less Time-Consuming Upwind Compact Difference Method with Adjusted Dissipation Property for Solving the Unsteady Incompressible Navier-Stokes Equations