Related papers: Numerical Method for Space-Time Fractional Diffusi…
This paper develops solutions of fractional Fokker-Planck equations describing subdiffusion of probability densities of stochastic dynamical systems driven by non-Gaussian L\'evy processes, with space-time-dependent drift, diffusion and…
Solving the Fokker-Planck equation for high-dimensional complex dynamical systems remains a pivotal yet challenging task due to the intractability of analytical solutions and the limitations of traditional numerical methods. In this work,…
The aim of this paper is to introduce a new Monte Carlo method based on importance sampling techniques for the simulation of stochastic differential equations. The main idea is to combine random walk on squares or rectangles methods with…
The fractional material derivative appears as the fractional operator that governs the dynamics of the scaling limits of L\'evy walks - a stochastic process that originates from the famous continuous-time random walks. It is usually defined…
In this paper, we consider a numerical method for the multi-term Caputo-Fabrizio time-fractional diffusion equations (with orders $\alpha_i\in(0,1)$, $i=1,2,\cdots,n$). The proposed method employs a fast finite difference scheme to…
The first passage time (FPT) problem is ubiquitous in many applications. In finance, we often have to deal with stochastic processes with jump-diffusion, so that the FTP problem is reducible to a stochastic differential equation with…
We study a time--space nonlocal diffusion equation driven by additive time--space white noise, where the time derivative is the Caputo derivative of order $\alpha\in(0,2)$. The model couples local diffusion with a nonlocal convolution…
We explain how the invariant subspace method can be extended to a scalar and coupled system of time-space fractional partial differential equations. The effectiveness and applicability of the method have been illustrated through time-space…
In this paper, we develop a Monte Carlo method for solving PDEs involving an integral fractional Laplacian (IFL) in multiple dimensions. We first construct a new Feynman-Kac representation based on the Green function for the fractional…
Hamiltonian Monte Carlo (HMC) sampling methods provide a mechanism for defining distant proposals with high acceptance probabilities in a Metropolis-Hastings framework, enabling more efficient exploration of the state space than standard…
Accurate prediction of rarefied gas flows is important for space vehicle design, particularly in rarefied regimes where the Navier-Stokes equations are no more valid. While the direct simulation Monte Carlo (DSMC) method acts as a numerical…
Stochastic collocation methods for approximating the solution of partial differential equations with random input data (e.g., coefficients and forcing terms) suffer from the curse of dimensionality whereby increases in the stochastic…
Recently, a class of efficient spectral Monte-Carlo methods was developed in \cite{Feng2025ExponentiallyAS} for solving fractional Poisson equations. These methods fully consider the low regularity of the solution near boundaries and…
The use of reaction-diffusion models rests on the key assumption that the underlying diffusive process is Gaussian. However, a growing number of studies have pointed out the prevalence of anomalous diffusion, and there is a need to…
A physical-mathematical approach to anomalous diffusion may be based on fractional diffusion equations and related random walk models. The fundamental solutions of these equations can be interpreted as probability densities evolving in time…
The basic problem in equilibrium statistical mechanics is to compute phase space average, in which Monte Carlo method plays a very important role. We begin with a review of nonlocal algorithms for Markov chain Monte Carlo simulation in…
In this paper, we propose a dynamically low-dimensional approximation method to solve a class of time-dependent multiscale stochastic diffusion equations. A dynamically bi-orthogonal (DyBO) method was developed to explore low-dimensional…
Solving multiscale diffusion problems is often computationally expensive due to the spatial and temporal discretization challenges arising from high-contrast coefficients. To address this issue, a partially explicit temporal splitting…
Over the past few decades, there has been substantial interest in evolution equations that involving a fractional-order derivative of order $\alpha\in(0,1)$ in time, due to their many successful applications in engineering, physics, biology…
High-dimensional fractional reaction-diffusion equations have numerous applications in the fields of biology, chemistry, and physics, and exhibit a range of rich phenomena. While classical algorithms have an exponential complexity in the…