Spectral meshless radial point interpolation (SMRPI) method to two‐dimensional fractional telegraph equation

H. Ammari In this article, an innovative technique so‐called spectral meshless radial point interpolation (SMRPI) method is proposed and, as a test problem, is applied to a classical type of two‐dimensional time‐fractional telegraph equation defined by Caputo sense for (1 < α≤2). This new methods is based on meshless methods and benefits from spectral collocation ideas, but it does not belong to traditional meshless collocation methods. The point interpolation method with the help of radial basis functions is used to construct shape functions, which play as basis functions in the frame of SMRPI method. These basis functions have Kronecker delta function property. Evaluation of high‐order derivatives is not difficult by constructing operational matrices. In SMRPI method, it does not require any kind of integration locally or globally over small quadrature domains, which is essential of the finite element method (FEM) and those meshless methods based on Galerkin weak form. Also, it is not needed to determine strict value for the shape parameter, which plays an important role in collocation method based on the radial basis functions (Kansa's method). Therefore, computational costs of SMRPI method are less expensive. Two numerical examples are presented to show that SMRPI method has reliable rates of convergence. Copyright © 2015 John Wiley & Sons, Ltd.     

Stability and convergence of spectral radial point interpolation method locally applied on two‐dimensional pseudoparabolic equation

In this article, we study a spectral meshless radial point interpolation of pseudoparabolic equations in two spatial dimensions. Shape functions, which are constructed through point interpolation method using the radial basis functions, help us to treat problem locally with the aim of high‐order convergence rate. The time derivatives are approximated by the finite difference time‐stepping method. The stability and convergence of this meshless approach are discussed and theoretically proven. Numerical results are presented to illustrate the theoretical findings. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 724–741, 2017     

A local radial basis function collocation method to solve the variable‐order time fractional diffusion equation in a two‐dimensional irregular domain

The local radial basis function (RBF) method is a promising solver for variable‐order time fractional diffusion equation (TFDE), as it overcomes the computational burden of the traditional global method. Application of the local RBF method is limited to Fickian diffusion, while real‐world diffusion is usually non‐Fickian in multiple dimensions. This article is the first to extend the application of the local RBF method to two‐dimensional, variable‐order, time fractional diffusion equation in complex shaped domains. One of the main advantages of the local RBF method is that only the nodes located in the subdomain, surrounding the local point, need to be considered when calculating the numerical solution at this point. This approach can perform well with large scale problems and can also mitigate otherwise ill‐conditioned problems. The proposed numerical approach is checked against two examples with curved boundaries and known analytical solutions. Shape parameter and subdomain node number are investigated for their influence on the accuracy of the local RBF solution. Furthermore, quantitative analysis, based on root‐mean‐square error, maximum absolute error, and maximum error of the partial derivative indicates that the local RBF method is accurate and effective in approximating the variable‐order TFDE in two‐dimensional irregular domains.     

A symmetric integrated radial basis function method for solving differential equations

In this article, integrated radial basis functions (IRBFs) are used for Hermite interpolation in the solution of differential equations, resulting in a new meshless symmetric RBF method. Both global and local approximation‐based schemes are derived. For the latter, the focus is on the construction of compact approximation stencils, where a sparse system matrix and a high‐order accuracy can be achieved together. Cartesian‐grid‐based stencils are possible for problems defined on nonrectangular domains. Furthermore, the effects of the RBF width on the solution accuracy for a given grid size are fully explored with a reasonable computational cost. The proposed schemes are numerically verified in some elliptic boundary‐value problems governed by the Poisson and convection‐diffusion equations. High levels of the solution accuracy are obtained using relatively coarse discretisations.     

Local radial point interpolation method for the fully developed magnetohydrodynamic flow

In this paper, a local radial point interpolation method (LRPIM) is presented to obtain the numerical solutions of the coupled equations in velocity and magnetic field for the fully developed magnetohydrodynamic (MHD) flow through a straight duct of rectangular section with arbitrary wall conductivity and orientation of applied magnetic field. Local weak forms are developed using weighted residual method locally for the governing equations of fully developed MHD flow. The shape functions from LRPIM possess the delta function property. Therefore, essential boundary conditions can be applied as easily as that in the finite-element method. The implementation procedure of LRPIM method is node based, and it doesn’t need any “mesh” or “element”. Computations have been carried out for different Hartmann numbers, wall conductivities and orientations of applied magnetic field.     

Solving hyperbolic conservation laws using multiquadric quasi‐interpolation

In this article, we apply the univariate multiquadric (MQ) quasi‐interpolation to solve the hyperbolic conservation laws. At first we construct the MQ quasi‐interpolation corresponding to periodic and inflow‐outflow boundary conditions respectively. Next we obtain the numerical schemes to solve the partial differential equations, by using the derivative of the quasi‐interpolation to approximate the spatial derivative of the differential equation and a low‐order explicit difference to approximate the temporal derivative of the differential equation. Then we verify our scheme for the one‐dimensional Burgers' equation (without viscosity). We can see that the numerical results are very close to the exact solution and the computational accuracy of the scheme is ??(τ), where τ is the temporal step. We can improve the accuracy by using the high‐order quasi‐interpolation. Moreover the methods can be generalized to the other equations. © 2005 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2006     

Error and stability analysis of numerical solution for the time fractional nonlinear Schrödinger equation on scattered data of general‐shaped domains

In present work, a kind of spectral meshless radial point interpolation (SMRPI) technique is applied to the time fractional nonlinear Schrödinger equation in regular and irregular domains. The applied approach is based on erudite combination of meshless methods and spectral collocation techniques. The point interpolation method with the help of radial basis functions is used to construct shape functions which play as basis functions in the frame of SMRPI. It is proved the scheme is unconditionally stable with respect to the time variable in and also convergent by the order of convergence , . In the current work, the thin plate spline are used as the basis functions and to eliminate the nonlinearity, a simple predictor‐corrector (P‐C) scheme is performed. It is shown that the SMRPI solution, as a complex function, is suitable one for the time fractional nonlinear Schrödinger equation. The results of numerical experiments are compared to analytical solutions to confirm the reliable treatment of these stable solutions. © 2016 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 33: 1043–1069, 2017     

A new multiquadric quasi‐interpolation operator with interpolation property

In this article, we discuss a class of multiquadric quasi‐interpolation operator that is primarily on the basis of Wu–Schaback's quasi‐interpolation operator and radial basis function interpolation. The proposed operator possesses the advantages of linear polynomial reproducing property, interpolation property, and high accuracy. It can be applied to construct flexible function approximation and scattered data fitting from numerical experiments. Copyright © 2013 John Wiley & Sons, Ltd.     

A steady barycentric Lagrange interpolation method for the 2D higher‐order time‐fractional telegraph equation with nonlocal boundary condition with error analysis

In this paper, we consider the numerical solution of the time‐fractional telegraph equation with a nonlocal boundary condition. A novel barycentric Lagrange interpolation collocation method is developed to solve this equation. Two difficulties have been sorted: the singularity of the integration and the higher accuracy. At the same, we put forward a steady barycentric Lagrange interpolation technique to overcome the new “Runge” phenomenon in computation. Error estimates of the barycentric Lagrange interpolation and the time‐fractional telegraph system for the present method are presented in Sobolev spaces. High convergence rates of the proposed method are obtained and are consisted with the numerical values. Especially in the time dimension, we get the error bound, for h‐refinement and for n_t‐density in the L₂ norms. The numerical results obtained show that the proposed numerical algorithm is accurate and computationally efficient for solving time‐fractional telegraph equation. Experiments demonstrate the high convergence rates of the proposed method are consisted with the theoretical values.     

Asymptotic behavior of solutions to space‐time fractional diffusion‐reaction equations

This article discusses the analyticity and the long‐time asymptotic behavior of solutions to space‐time fractional diffusion‐reaction equations in . By a Laplace transform argument, we prove that the decay rate of the solution as t→∞ is dominated by the order of the time‐fractional derivative. We consider the decay rate also in a bounded domain. Copyright © 2016 John Wiley & Sons, Ltd.     

A meshless method based on boundary integral equations and radial basis functions for biharmonic-type problems

This paper presents a meshless method, which replaces the inhomogeneous biharmonic equation by two Poisson equations in terms of an intermediate function. The solution of the Poisson equation with the intermediate function as the right-hand term may be written as a sum of a particular solution and a homogeneous solution of a Laplace equation. The intermediate function is approximated by a series of radial basis functions. Then the particular solution is obtained via employing Kansa’s method, while the homogeneous solution is approximated by using the boundary radial point interpolation method by means of boundary integral equations. Besides, the proposed meshless method, in conjunction with the analog equation method, is further developed for solving generalized biharmonic-type problems. Some numerical tests illustrate the efficiency of the method proposed.     

Finite element method for two‐dimensional time‐fractional tricomi‐type equations

In this article, we consider the finite element method (FEM) for two‐dimensional linear time‐fractional Tricomi‐type equations, which is obtained from the standard two‐dimensional linear Tricomi‐type equation by replacing the first‐order time derivative with a fractional derivative (of order α, with 1 <α< 2 ). The method is based on finite element method for space and finite difference method for time. We prove that the method is unconditionally stable, and the error estimate is presented. The comparison of the FEM results with the exact solutions is made, and numerical experiments reveal that the FEM is very effective. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2013     

A new residual posteriori error estimates of mixed finite element methods for convection‐diffusion‐reaction equations

We propose and analyze a new technique for developing residual‐based a posteriori error estimates over the stress and scalar displacement error for the lowest‐order Raviart–Thomas mixed finite element discretizations of convection‐diffusion‐reaction equations in two‐dimension space. The new technique is based on the abstract error estimates, the postprocessed approximation of the scalar displacement, and on the construction of an auxiliary problem. We consider the centered and upwind‐weighted mixed schemes, and concentrate the attention on the presence of an inhomogeneous and an anisotropic diffusion‐dispersion tensor and on a possible convection dominance. Global upper bounds can be directly computed on the base of the solution of the mixed schemes without any additional cost. Local lower bounds without any saturation assumption, hold from the case where convection or reaction are not present to convection‐ or reaction‐dominated equations, and their local efficiency depends on local or global variations in coefficients similar to Péclect number. Numerical experiments are reported to show the competitive behavior of the proposed posteriori error estimates, and to confirm the theoretical findings. © 2013 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 593–624, 2014     

Solving 2D time‐fractional diffusion equations by a pseudospectral method and Mittag‐Leffler function evaluation

Two‐dimensional time‐fractional diffusion equations with given initial condition and homogeneous Dirichlet boundary conditions in a bounded domain are considered. A semidiscrete approximation scheme based on the pseudospectral method to the time‐fractional diffusion equation leads to a system of ordinary fractional differential equations. To preserve the high accuracy of the spectral approximation, an approach based on the evaluation of the Mittag‐Leffler function on matrix arguments is used for the integration along the time variable. Some examples along with numerical experiments illustrate the effectiveness of the proposed approach. Copyright © 2016 John Wiley & Sons, Ltd.     

Exponential convergence and H‐c multiquadric collocation method for partial differential equations

The radial basis function (RBF) collocation method uses global shape functions to interpolate and collocate the approximate solution of PDEs. It is a truly meshless method as compared to some of the so‐called meshless or element‐free finite element methods. For the multiquadric and Gaussian RBFs, there are two ways to make the solution converge—either by refining the mesh size h, or by increasing the shape parameter c. While the h‐scheme requires the increase of computational cost, the c‐scheme is performed without extra effort. In this paper we establish by numerical experiment the exponential error estimate ? ～ O(λ^√c?h) where 0 < λ < 1. We also propose the use of residual error as an error indicator to optimize the selection of c. © 2003 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 19: 571–594, 2003     

A finite element–finite volume discretization of convection‐diffusion‐reaction equations with nonhomogeneous mixedboundary conditions: Error estimates

We consider a time‐dependent and a steady linear convection‐diffusion‐reaction equation whose coefficients are nonconstant. Boundary conditions are mixed (Dirichlet and Robin–Neumann) and nonhomogeneous. Both the unsteady and the steady problem are approximately solved by a combined finite element–finite volume method: the diffusion term is discretized by Crouzeix–Raviart piecewise linear finite elements on a triangular grid, and the convection term by upwind barycentric finite volumes. In the unsteady case, the implicit Euler method is used as time discretization. The ‐ and the ‐error in the unsteady case and the H¹‐error in the steady one are estimated against the data, in such a way that no parameter enters exponentially into the constants involved. © 2016Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 1591–1621, 2016     

Multiquadric quasi‐interpolation methods for solving partial differential algebraic equations

In this article, we propose two meshless collocation approaches for solving time dependent partial differential algebraic equations (PDAEs) in terms of the multiquadric quasi‐interpolation schemes. In presenting the process of the solution, the error is estimated. Furthermore, the comparisons on condition numbers of the collocation matrices using different methods and the sensitivity of the shape parameter c are given. With the use of the appropriate collocation points, the method for PDAEs with index‐2 is improved. The results show that the methods have some advantages over some known methods, such as the smaller condition numbers or more accurate solutions for PDAEs which has an modal index‐2 or an impulse solution with index‐2. Copyright © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 30: 95–119, 2014     

Finite‐time stability of fractional‐order stochastic singular systems with time delay and white noise

In this article, the finite‐time stochastic stability of fractional‐order singular systems with time delay and white noise is investigated. First the existence and uniqueness of solution for the considered system is derived using the basic fractional calculus theory. Then based on the Gronwall's approach and stochastic analysis technique, the sufficient condition for the finite‐time stability criterion is developed. Finally, a numerical example is presented to verify the obtained theory. © 2016 Wiley Periodicals, Inc. Complexity 21: 370–379, 2016     

A new stabilized finite element method for advection‐diffusion‐reaction equations

In this article, a new stabilized finite element method is proposed and analyzed for advection‐diffusion‐reaction equations. The key feature is that both the mesh‐dependent Péclet number and the mesh‐dependent Damköhler number are reasonably incorporated into the newly designed stabilization parameter. The error estimates are established, where, up to the regularity‐norm of the exact solution, the explicit‐dependence of the diffusivity, advection, reaction, and mesh size (or the dependence of the mesh‐dependent Péclet number and the mesh‐dependent Damköhler number) is revealed. Such dependence in the error bounds provides a mathematical justification on the effectiveness of the proposed method for any values of diffusivity, advection, dissipative reaction, and mesh size. Numerical results are presented to illustrate the performance of the method. © 2015 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 32: 616–645, 2016     

Use of radial basis functions for solving the second‐order parabolic equation with nonlocal boundary conditions

Nonlocal mathematical models appear in various problems of physics and engineering. In these models the integral term may appear in the boundary conditions. In this paper the problem of solving the one‐dimensional parabolic partial differential equation subject to given initial and nonlocal boundary conditions is considered. These kinds of problems have certainly been one of the fastest growing areas in various application fields. The presence of an integral term in a boundary condition can greatly complicate the application of standard numerical techniques. As a well‐known class of meshless methods, the radial basis functions are used for finding an approximation of the solution of the present problem. Numerical examples are given at the end of the paper to compare the efficiency of the radial basis functions with famous finite‐difference methods. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008