Differential quadrature method for space-fractional diffusion equations on 2D irregular domains

X. G. Zhu,Z. B. Yuan,F. Liu,Y. F. Nie
DOI: https://doi.org/10.48550/arXiv.1703.10308
2018-01-02
Abstract:In mathematical physics, the space-fractional diffusion equations are of particular interest in the studies of physical phenomena modelled by Lévy processes, which are sometimes called super-diffusion equations. In this article, we develop the differential quadrature (DQ) methods for solving the 2D space-fractional diffusion equations on irregular domains. The methods in presence reduce the original equation into a set of ordinary differential equations (ODEs) by introducing valid DQ formulations to fractional directional derivatives based on the functional values at scattered nodal points on problem domain. The required weighted coefficients are calculated by using radial basis functions (RBFs) as trial functions, and the resultant ODEs are discretized by the Crank-Nicolson scheme. The main advantages of our methods lie in their flexibility and applicability to arbitrary domains. A series of illustrated examples are finally provided to support these points.
Numerical Analysis
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to solve the spatial fractional - order diffusion equation on two - dimensional irregular domains. Specifically, the authors developed the differential quadrature (DQ) method to handle such equations. These methods can transform the original spatial fractional - order partial differential equations into a set of ordinary differential equations (ODEs), thereby simplifying the solution process of the problem. The following is a detailed description of this problem: ### Problem Background In mathematical physics, the spatial fractional - order diffusion equation is an important tool for studying physical phenomena modeled by Lévy processes, and is sometimes also called the super - diffusion equation. Such equations perform well in describing physical processes with memory effects and historical dependencies, such as anomalous dispersion in complex heterogeneous aquifers and fractal geometries. ### Research Objectives The main objective of the paper is to develop a differential quadrature (DQ) method for the spatial fractional - order diffusion equation on two - dimensional irregular domains. These methods transform the original equation into a set of ordinary differential equations (ODEs) by introducing reasonable DQ formulas based on the function values at discrete nodes within the problem domain. Specifically, the paper addresses the following aspects: 1. **Applicability**: Most of the existing numerical algorithms are only applicable to one - dimensional or rectangular domain problems, while the method proposed in this paper can be applied to two - dimensional irregular domains of arbitrary shapes. 2. **Flexibility**: The proposed method uses radial basis functions (RBFs) as trial functions to calculate the required weighting coefficients and adopts the Crank - Nicolson scheme to discretize the resulting system of ordinary differential equations, thereby improving the flexibility and applicability of the method. 3. **Efficiency**: By using the differential quadrature method, high - precision numerical solutions can be achieved with fewer nodes, while avoiding the complexity and computational burden brought by traditional mesh generation methods. ### Mathematical Model The paper considers the spatial fractional - order diffusion equation in the following form: \[ \frac{\partial u(x,y,t)}{\partial t}-\kappa\int_{0}^{2\pi}D_{\theta}^{\alpha}u(x,y,t)P(\theta)d\theta = f(x,y,t),\quad(x,y;t)\in\Omega\times(0,T] \] where $\Omega$ is a bounded region, $\partial\Omega$ is its boundary, $\kappa$ is a non - negative diffusion parameter, $P(\theta)d\theta$ is a probability measure on the unit disk, and $f(x,y,t)$ is a preset source function. $D_{\theta}^{\alpha}u(x,y,t)$ represents the spatial fractional - order directional derivative in the Caputo sense, defined as: \[ D_{\theta}^{\alpha}u(x,y,t)=I_{\theta}^{2-\alpha}D_{\theta}^{2}u(x,y,t),\quad0\leq\theta < 2\pi \] ### Method Innovation Points The main innovation points of the paper are: - Proposing a differential quadrature method for the spatial fractional - order diffusion equation applicable to two - dimensional irregular domains; - Using radial basis functions (RBFs) to calculate weighting coefficients, improving the flexibility and applicability of the method; - Discretizing the system of ordinary differential equations through the Crank - Nicolson scheme to ensure the stability of time stepping. In summary, this paper aims to provide a flexible and efficient numerical method for solving the spatial fractional - order diffusion equation on two - dimensional irregular domains to address the limitations of existing methods in complex domain problems.