On the solution of the non-local parabolic partial differential equations via radial basis functions

In this paper, the problem of solving the one-dimensional parabolic partial differential equation subject to given initial and non-local boundary conditions is considered. The approximate solution is found using the radial basis functions collocation method. There are some difficulties in computing the solution of the time dependent partial differential equations using radial basis functions. If time and space are discretized using radial basis functions, the resulted coefficient matrix will be very ill-conditioned and so the corresponding linear system cannot be solved easily. As an alternative method for solution, we can use finite-difference methods for discretization of time and radial basis functions for discretization of space. Although this method is easy to use but an accurate solution cannot be provided. In this work an efficient collocation method is proposed for solving non-local parabolic partial differential equations using radial basis functions. Numerical results are presented and are compared with some existing methods.     

Numerical solutions of KdV equation using radial basis functions

Numerical solution of the Korteweg-de Vries equation is obtained by using the meshless method based on the collocation with radial basis functions. Five standard radial basis functions are used in the method of the collocation. The results are compared for the numerical experiments of the propagation of solitons, interaction of two solitary waves and breakdown of initial conditions into a train of solitons.     

Multiscale RBF collocation for solving PDEs on spheres

In this paper, we discuss multiscale radial basis function collocation methods for solving certain elliptic partial differential equations on the unit sphere. The approximate solution is constructed in a multi-level fashion, each level using compactly supported radial basis functions of smaller scale on an increasingly fine mesh. Two variants of the collocation method are considered (sometimes called symmetric and unsymmetric, although here both are symmetric). A convergence theory is given, which builds on recent theoretical advances for multiscale approximation using compactly supported radial basis functions.     

Smooth solution of partial integro-differential equations using radial basis functions

In this work, we present the method based on radial basis functions to solve partial integro-differential equations. We focus on the parabolic type of integro-differential equations as the most common forms including the ``\emph{memory}'' of the systems. We propose to apply the collocation scheme using radial basis functions to approximate the solutions of partial integro-differential equations. Due to the presented technique, system of linear or nonlinear equations is made instead of primary problem. The method is efficient because the rate of convergence of collocation method based on radial basis functions is exponential. Some numerical examples and investigation of the experimental results show the applicability and accuracy of the method.     

Numerical solutions of the Kawahara type equations using radial basis functions

This study is carried out to investigate the numerical solutions of the Kawahara, KdV‐Kawahara, and the modified Kawahara equations by using the meshless method based on collocation with radial basis functions. Results of the meshless method with different radial basis functions are presented for the travelling wave solution of the Kawahara type equations. © 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 28: 542–553, 2012     

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.     

Error analysis of collocation method based on reproducing kernel approximation

Solving partial differential equations (PDE) with strong form collocation and nonlocal approximation functions such as orthogonal polynomials, trigonometric functions, and radial basis functions exhibits exponential convergence rates; however, it yields a full matrix and suffers from ill conditioning. In this work, we discuss a reproducing kernel collocation method, where the reproducing kernel (RK) shape functions with compact support are used as approximation functions. This approach offers algebraic convergence rate, but the method is stable like the finite element method. We provide mathematical results consisting of the optimal error estimation, upper bound of condition number, and the desirable relationship between the number of nodal points and the number of collocation points. We show that using RK shape function for collocation of strong form, the degree of polynomial basis functions has to be larger than one for convergence, which is different from the condition for weak formulation. Numerical results are also presented to validate the theoretical analysis. © 2009 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 27: 554–580, 2011     

Adaptive multiquadric collocation for boundary layer problems

An adaptive collocation method based upon radial basis functions is presented for the solution of singularly perturbed two-point boundary value problems. Using a multiquadric integral formulation, the second derivative of the solution is approximated by multiquadric radial basis functions. This approach is combined with a coordinate stretching technique. The required variable transformation is accomplished by a conformal mapping, an iterated sine-transformation. A new error indicator function accurately captures the regions of the interval with insufficient resolution. This indicator is used to adaptively add data centres and collocation points. The method resolves extremely thin layers accurately with fairly few basis functions. The proposed adaptive scheme is very robust, and reaches high accuracy even when parameters in our coordinate stretching technique are not chosen optimally. The effectiveness of our new method is demonstrated on two examples with boundary layers, and one example featuring an interior layer. It is shown in detail how the adaptive method refines the resolution.     

Adaptive multiscale methods for elliptic PDEs

This paper deals with the numerical solution of elliptic partial differential equations with an adaptive multiscale algorithm, which uses symmetric collocation and radial basis functions. (© 2017 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

Dispersion and stability properties of radial basis collocation method for elastodynamics

Strong form collocation with radial basis approximation, called the radial basis collocation method (RBCM), is introduced for the numerical solution of elastodynamics. In this work, the proper weights for the boundary collocation equations to achieve the optimal convergence in elastodynamics are first derived. The von Neumann method is then introduced to investigate the dispersion characteristics of the semidiscrete RBCM equation. Very small dispersion error (< 1%) in RBCM can be achieved compared to linear and quadratic finite elements. The stability conditions of the RBCM spatial discretization in conjunction with the central difference temporal discretization are also derived. We show that the shape parameter of the radial basis functions not only has strong influence on the dispersion errors, it also has profound influence on temporal stability conditions in the case of lumped mass. Further, our stability analysis shows that, in general, a larger critical time step can be used in RBCM with central difference temporal discretization than that for finite elements with the same temporal discretization. Our analysis also suggests that although RBCM with lumped mass allows a much larger critical time step than that of RBCM with consistent mass, the later offers considerably better accuracy and should be considered in the transient analysis. © 2012 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 2013     

Solitary wave solutions of the MRLW equation using radial basis functions

In this study, traveling wave solutions of the modified regularized long wave (MRLW) equation are simulated by using the meshless method based on collocation with well‐known radial basis functions. The method is tested for three test problems which are single solitary wave motion, interaction of two solitary waves and interaction of three solitary waves. Invariant values for all test problems are calculated, also L₂, L_∞ norms and values of the absolute error for single solitary wave motion are calculated. Numerical results by using the meshless method with different radial basis functions are presented. Figures of wave motions for all test problems are shown. Altogether, meshless methods with radial basis functions solve the MRLW equation very satisfactorily.© 2010 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq 28: 235–247, 2012     

Numerical solution of nonlinear Volterra-Fredholm-Hammerstein integral equations via collocation method based on radial basis functions

A numerical technique based on the spectral method is presented for the solution of nonlinear Volterra-Fredholm-Hammerstein integral equations. This method is a combination of collocation method and radial basis functions (RBFs) with the differentiation process (DRBF), using zeros of the shifted Legendre polynomial as the collocation points. Different applications of RBFs are used for this purpose. The integral involved in the formulation of the problems are approximated based on Legendre-Gauss-Lobatto integration rule. The results of numerical experiments are compared with the analytical solution in illustrative examples to confirm the accuracy and efficiency of the presented scheme.     

An Adaptive Greedy Algorithm for Solving Large RBF Collocation Problems

The solution of operator equations with radial basis functions by collocation in scattered points leads to large linear systems which often are nonsparse and ill-conditioned. But one can try to use only a subset of the data for the actual collocation, leaving the rest of the data points for error checking. This amounts to finding sparse approximate solutions of general linear systems arising from collocation. This contribution proposes an adaptive greedy method with proven (but slow) linear convergence to the full solution of the collocation equations. The collocation matrix need not be stored, and the progress of the method can be controlled by a variety of parameters. Some numerical examples are given.     

Strongly elliptic pseudodifferential equations on the sphere with radial basis functions

Spherical radial basis functions are used to define approximate solutions to strongly elliptic pseudodifferential equations on the unit sphere. These equations arise from geodesy. The approximate solutions are found by the Galerkin and collocation methods. A salient feature of the paper is a unified theory for error analysis of both approximation methods.     

Numerical solution of three-dimensional problems of filtration consolidation with regard for the influence of technogenic factors by the method of radial basis functions

We use a meshless method to find an approximate solution of the problem that describes a mathematical model of the filtration consolidation process in a three-dimensional domain. It is based on the collocation method using radial basis functions. The performed numerical experiments testify to the efficiency of the proposed approach.     

A local meshless method for time fractional nonlinear diffusion wave equation

Numerical Algorithms - We present a radial basis function-based local collocation method for solving time fractional nonlinear diffusion wave equation.The main beauty of the local collocation...     

Meshless method with ridge basis functions

Meshless collocation methods for the numerical solutions of PDEs are increasingly adopted due to their advantages including efficiency and flexibility, and radial basis functions are popularly employed to represent the solutions of PDEs. Motivated by the advantages of ridge basis function representation of a given function, such as the connection to neural network, fast convergence as the number of terms is increased, better approximation effects and various applications in engineering problems, a meshless method is developed based on the collocation method and ridge basis function interpolation. This method is a truly meshless technique without mesh discretization: it neither needs the computation of integrals, nor requires a partition of the region and its boundary. Moreover, the method is applied to elliptic equations to examine its appropriateness, numerical results are compared to that obtained from other (meshless) methods, and influence factors of accuracy for numerical solutions are analyzed.     

Radial basis collocation method and quasi‐Newton iteration for nonlinear elliptic problems

This work presents a radial basis collocation method combined with the quasi‐Newton iteration method for solving semilinear elliptic partial differential equations. The main result in this study is that there exists an exponential convergence rate in the radial basis collocation discretization and a superlinear convergence rate in the quasi‐Newton iteration of the nonlinear partial differential equations. In this work, the numerical error associated with the employed quadrature rule is considered. It is shown that the errors in Sobolev norms for linear elliptic partial differential equations using radial basis collocation method are bounded by the truncation error of the RBF. The combined errors due to radial basis approximation, quadrature rules, and quasi‐Newton and Newton iterations are also presented. This result can be extended to finite element or finite difference method combined with any iteration methods discussed in this work. The numerical example demonstrates a good agreement between numerical results and analytical predictions. The numerical results also show that although the convergence rate of order 1.62 of the quasi‐Newton iteration scheme is slightly slower than rate of order 2 in the Newton iteration scheme, the former is more stable and less sensitive to the initial guess. © 2007 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2008     

A meshless method for numerical solution of a linear hyperbolic equation with variable coefficients in two space dimensions

A meshless method is proposed for the numerical solution of the two space dimensional linear hyperbolic equation subject to appropriate initial and Dirichlet boundary conditions. The new developed scheme uses collocation points and approximates the solution employing thin plate splines radial basis functions. Numerical results are obtained for various cases involving variable, singular and constant coefficients, and are compared with analytical solutions to confirm the good accuracy of the presented scheme. © 2008 Wiley Periodicals, Inc. Numer Methods Partial Differential Eq, 2009     

A meshless method for numerical solution of the one-dimensional wave equation with an integral condition using radial basis functions

The hyperbolic partial differential equation with an integral condition arises in many physical phenomena. In this paper, we propose a numerical scheme to solve the one-dimensional hyperbolic equation that combines classical and integral boundary conditions using collocation points and approximating the solution using radial basis functions (RBFs). The results of numerical experiments are presented, and are compared with analytical solution and finite difference method to confirm the validity and applicability of the presented scheme.