An advanced meshless method for time fractional diffusion equation

Recently, because of the new developments in sustainable engineering and renewable energy, which are usually governed by a series of fractional partial differential equations (FPDEs), the numerical modelling and simulation for fractional calculus are attracting more and more attention from researchers. The current dominant numerical method for modeling FPDE is Finite Difference Method (FDM), which is based on a pre-defined grid leading to inherited issues or shortcomings including difficulty in simulation of problems with the complex problem domain and in using irregularly distributed nodes. Because of its distinguished advantages, the meshless method has good potential in simulation of FPDEs. This paper aims to develop an implicit meshless collocation technique for FPDE. The discrete system of FPDEs is obtained by using the meshless shape functions and the meshless collocation formulation. The stability and convergence of this meshless approach are investigated theoretically and numerically. The numerical examples with regular and irregular nodal distributions are used to validate and investigate accuracy and efficiency of the newly developed meshless formulation. It is concluded that the present meshless formulation is very effective for the modeling and simulation of fractional partial differential equations

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

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.Comment: 25 pages, 25 figures, 7 table

MESHLESS METHODS FOR SOLVING REACTION-DIFFUSION PROBLEMS-A BRIEF REVIEW

Reaction-diffusion equations represent many important and critical applications in engineering and science. Numerical techniques play an important role for solving such equations accurately and efficiently. This paper presents a brief review of meshless methods for solving general diffusion equations, including reaction-diffusion systems

An alternating direction implicit spectral method for solving two dimensional multi-term time fractional mixed diffusion and diffusion-wave equations

In this paper, we consider the initial boundary value problem of the two dimensional multi-term time fractional mixed diffusion and diffusion-wave equations. An alternating direction implicit (ADI) spectral method is developed based on Legendre spectral approximation in space and finite difference discretization in time. Numerical stability and convergence of the schemes are proved, the optimal error is $O(N^{-r}+\tau^2)$, where $N, \tau, r$ are the polynomial degree, time step size and the regularity of the exact solution, respectively. We also consider the non-smooth solution case by adding some correction terms. Numerical experiments are presented to confirm our theoretical analysis. These techniques can be used to model diffusion and transport of viscoelastic non-Newtonian fluids

Fast Method of Particular Solutions for Solving Partial Differential Equations

Method of particular solutions (MPS) has been implemented in many science and engineering problems but obtaining the closed-form particular solutions, the selection of the good shape parameter for various radial basis functions (RBFs) and simulation of the large-scale problems are some of the challenges which need to overcome. In this dissertation, we have used several techniques to overcome such challenges. The closed-form particular solutions for the Matérn and Gaussian RBFs were not known yet. With the help of the symbolic computational tools, we have derived the closed-form particular solutions of the Matérn and Gaussian RBFs for the Laplace and biharmonic operators in 2D and 3D. These derived particular solutions play an important role in solving inhomogeneous problems using MPS and boundary methods such as boundary element methods or boundary meshless methods. In this dissertation, to select the good shape parameter, various existing variable shape parameter strategies and some well-known global optimization algorithms have also been applied. These good shape parameters provide high accurate solutions in many RBFs collocation methods. Fast method of particular solutions (FMPS) has been developed for the simulation of the large-scale problems. FMPS is based on the global version of the MPS. In this method, partial differential equations are discretized by the usual MPS and the determination of the unknown coefficients is accelerated using a fast technique. Numerical results confirm the efficiency of the proposed technique for the PDEs with a large number of computational points in both two and three dimensions. We have also solved the time fractional diffusion equations by using MPS and FMPS

A RBF-based differential quadrature method for solving two-dimensional variable-order time fractional advection-diffusion equation

The file attached to this record is the author's final peer reviewed version. The Publisher's final version can be found by following the DOI link.Numerical simulation technique of two-dimensional variable-order time fractional advection-diffusion equation is developed in this paper using radial basis function-based differential quadrature method (RBF-DQ). To the best of the authors’ knowledge, this is the first application of this method to variable-order time fractional advection-diffusion equations. For the general case of irregular geometries, the meshless local form of RBF-DQ is used and the multiquadric type of radial basis functions is selected for the computations. This approach allows one to define a reconstruction of the local radial basis functions to treat accurately both the Dirichlet and Neumann boundary conditions on the irregular boundaries. The method is validated by the well documented test examples involving variable-order fractional modelling of air pollution. The numerical results demonstrate that the proposed method provides accurate solutions fortwo-dimensional variable-order time fractional advection-diffusion equations