No Arabic abstract
In mathematical physics, the space-fractional diffusion equations are of particular interest in the studies of physical phenomena modelled by L{e}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.
This article studies a direct numerical approach for fractional advection-diffusion equations (ADEs). Using a set of cubic trigonometric B-splines as test functions, a differential quadrature (DQ) method is firstly proposed for the 1D and 2D time-fractional ADEs of order $(0,1]$. The weighted coefficients are determined, and with them, the original equation is transformed into a group of general ordinary differential equations (ODEs), which are discretized by an effective difference scheme or Runge-Kutta method. The stability is investigated under a mild theoretical condition. Secondly, based on a set of cubic B-splines, we develop a new Crank-Nicolson type DQ method for the 2D space-fractional ADEs without advection. The DQ approximations to fractional derivatives are introduced and the values of the fractional derivatives of B-splines are computed by deriving explicit formulas. The presented DQ methods are evaluated on five benchmark problems and the concrete simulations of the unsteady propagation of solitons and Gaussian pulse. In comparison with the existing algorithms in the open literature, numerical results finally illustrate the validity and accuracy.
We propose and test the first Reduced Radial Basis Function Method (R$^2$BFM) for solving parametric partial differential equations on irregular domains. The two major ingredients are a stable Radial Basis Function (RBF) solver that has an optimized set of centers chosen through a reduced-basis-type greedy algorithm, and a collocation-based model reduction approach that systematically generates a reduced-order approximation whose dimension is orders of magnitude smaller than the total number of RBF centers. The resulting algorithm is efficient and accurate as demonstrated through two- and three-dimensional test problems.
This article aims to develop a direct numerical approach to solve the space-fractional partial differential equations (PDEs) based on a new differential quadrature (DQ) technique. The fractional derivatives are approximated by the weighted linear combinations of the function values at discrete grid points on problem domain with the weights calculated via using three types of radial basis functions (RBFs) as test functions. The method in presence is robust, straight forward to apply, and highly accurate under the condition that the shape parameters of RBFs are well chosen. Numerical tests are provided to illustrate its validity and capability.
We propose some multigrid methods for solving the algebraic systems resulting from finite element approximations of space fractional partial differential equations (SFPDEs). It is shown that our multigrid methods are optimal, which means the convergence rates of the methods are independent of the mesh size and mesh level. Moreover, our theoretical analysis and convergence results do not require regularity assumptions of the model problems. Numerical results are given to support our theoretical findings.
In this work, we present a collocation method based on the Legendre wavelet combined with the Gauss--Jacobi quadrature formula for solving a class of fractional delay-type integro-differential equations. The problem is considered with either initial or boundary conditions and the fractional derivative is described in the Caputo sense. First, an approximation of the unknown solution is considered in terms of the Legendre wavelet basis functions. Then, we substitute this approximation and its derivatives into the considered equation. The Caputo derivative of the unknown function is approximated using the Gauss--Jacobi quadrature formula. By collocating the obtained residual at the well-known shifted Chebyshev points, we get a system of nonlinear algebraic equations. In order to obtain a continuous solution, some conditions are added to the resulting system. Some error bounds are given for the Legendre wavelet approximation of an arbitrary function. Finally, some examples are included to show the efficiency and accuracy of this new technique.