Scattered data fitting of Hermite type by a weighted meshless method

In many practical problems, it is often desirable to interpolate not only the function values but also the values of derivatives up to certain order, as in the Hermite interpolation. The Hermite interpolation method by radial basis functions is used widely for solving scattered Hermite data approximation problems. However, sometimes it makes more sense to approximate the solution by a least squares fit. This is particularly true when the data are contaminated with noise. In this paper, a weighted meshless method is presented to solve least squares problems with noise. The weighted meshless method by Gaussian radial basis functions is proposed to fit scattered Hermite data with noise in certain local regions of the problem’s domain. Existence and uniqueness of the solution is proved. This approach has one parameter which can adjust the accuracy according to the size of the noise. Another advantage of the weighted meshless method is that it can be used for problems in high dimensions with nonregular domains. The numerical experiments show that our weighted meshless method has better performance than the traditional least squares method in the case of noisy Hermite data.     

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     

Optimal control of a parabolic distributed parameter system via radial basis functions

This paper attempts to present a meshless method to find the optimal control of a parabolic distributed parameter system with a quadratic cost functional. The method is based on radial basis functions to approximate the solution of the optimal control problem using collocation method. In this regard, different applications of RBFs are used. To this end, the numerical solutions are obtained without any mesh generation into the domain of the problems. The proposed technique is easy to implement, efficient and yields accurate results. Numerical examples are included and a comparison is made with an existing result.     

Meshless method of lines for the numerical solution of generalized Kuramoto-Sivashinsky equation

In this paper, the numerical solution of the generalized Kuramoto-Sivashinsky equation is presented by meshless method of lines (MOL). In this method the spatial derivatives are approximated by radial basis functions (RBFs) giving an edge over finite difference method (FDM) and finite element method (FEM) because no mesh is required for discretization of the problem domain. Only a set of scattered nodes is required to approximate the solution. The numerical results in comparison with exact solution using different radial basis functions (RBFs) prove the efficiency and accuracy of the method.     

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.     

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.     

On the numerical solution of nonlinear Burgers’-type equations using meshless method of lines

In this paper, a meshless method of lines (MMOL) is proposed for the numerical solution of nonlinear Burgers’-type equations. This technique does not require a mesh in the problem domain, and only a set of scattered nodes provided by initial data is required for the solution of the problem using some radial basis functions (RBFs). The scheme is tested for various examples. The results obtained by this method are compared with the exact solutions and some earlier work.     

A meshless procedure for transient heat conduction in functionally graded materials

A meshless numerical procedure is developed for analyzing the transient heat conduction problem in non-homogeneous functionally graded materials. In the proposed method the time derivative of temperature is approximate by the finite difference. At each time step the original nonlinear boundary value problem is transform into a hierarchy of non-homogeneous linear problem by used the homotopy analysis method. In this method a sought solution is presented by using a finite expansion in Taylor series, which consecutive elements are solutions of series linear value problems defining differential deformations. Each of linear boundary value problems with the corresponding boundary conditions is solved by using the method of fundamental solutions and radial basis functions which are used for interpolation of the inhomogeneous term. The accuracy of the obtained approximate solution is controlled by the number of components of the Taylor series, while the convergence of the process is monitored by an additional parameter of the method. Numerical experiments demonstrate the efficiency and accuracy of the present scheme in the solution of the heat conduction problem in nonlinear functionally graded materials. (© 2013 Wiley-VCH Verlag GmbH & Co. KGaA, Weinheim)     

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 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.     

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.     

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     

Analysis of a meshless method for the time fractional diffusion-wave equation

In this paper a numerical technique is proposed for solving the time fractional diffusion-wave equation. We obtain a time discrete scheme based on finite difference formula. Then, we prove that the time discrete scheme is unconditionally stable and convergent using the energy method and the convergence order of the time discrete scheme is \(\mathcal {O}(\tau ^{3-\alpha })\). Firstly, we change the main problem based on Dirichlet boundary condition to a new problem based on Robin boundary condition and then, we consider a semi-discrete scheme with Robin boundary condition and show when \(\beta \rightarrow +\infty \) solution of the main semi-discrete problem with Dirichlet boundary condition is convergent to the solution of the new semi-discrete problem with Robin boundary condition. We consider the new semi-discrete problem with Robin boundary condition and use the meshless Galerkin method to approximate the spatial derivatives. Finally, we obtain an error bound for the new problem. We prove that convergence order of the numerical scheme based on Galekin meshless is \(\mathcal {O}(h)\). In the considered method the appeared integrals are approximated using Gauss Legendre quadrature formula. The main aim of the current paper is to obtain an error estimate for the meshless Galerkin method based on the radial basis functions. Numerical examples confirm the efficiency and accuracy of the proposed scheme.     

A numerical algorithm to reduce ill-conditioning in meshless methods for the Helmholtz equation

Some meshless methods have been applied to the numerical solution of boundary value problems involving the Helmholtz equation. In this work, we focus on the method of fundamental solutions and the plane waves method. It is well known that these methods can be highly accurate assuming smoothness of the domains and the boundary data. However, the matrices involved are often ill-conditioned and the effect of this ill-conditioning may drastically reduce the accuracy. In this work, we propose a numerical algorithm to reduce the ill-conditioning in both methods. The idea is to perform a suitable change of basis. This allows to obtain new basis functions that span exactly the same space as the original meshless method, but are much better conditioned. In the case of circular domains, this technique allows to obtain errors close to machine precision, with condition numbers of order O(1), independently of the number of basis functions in the expansion.     

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     

Solving differential equations with radial basis functions: multilevel methods and smoothing

Some of the meshless radial basis function methods used for the numerical solution of partial differential equations are reviewed. In particular, the differences between globally and locally supported methods are discussed, and for locally supported methods the important role of smoothing within a multilevel framework is demonstrated. A possible connection between multigrid finite elements and multilevel radial basis function methods with smoothing is explored. Various numerical examples are also provided throughout the paper. This revised version was published online in June 2006 with corrections to the Cover Date.     

Radial basis functions method for solving of a non-local boundary value problem with Neumann’s boundary conditions

In this paper, the problem of solving the two-dimensional diffusion equation subject to a non-local condition involving a double integral in a rectangular region is considered. The solution of this type of problems are complicated. Therefore, a simple meshless method using the radial basis functions is constructed for the non-local boundary value problem with Neumann’s boundary conditions. Numerical examples are included to demonstrate the reliability and efficiency of this method. Also N_e and root mean square errors are obtained to show the convergence of the method.     

A dynamic meshless method for the least squares problem with some noisy subdomains

The interpolation method by radial basis functions is used widely for solving scattered data approximation. However, sometimes it makes more sense to approximate the solution by least squares fit. This is especially true when the data are contaminated with noise. A meshfree method namely, meshless dynamic weighted least squares (MDWLS) method, is presented in this paper to solve least squares problems with noise. The MDWLS method by Gaussian radial basis function is proposed to fit scattered data with some noisy areas in the problem’s domain. Existence and uniqueness of a solution is proved. This method has one parameter which can adjusts the accuracy according to the size of noises. Another advantage of the developed method is that it can be applied to problems with nonregular geometrical domains. The new approach is applied for some problems in two dimensions and the obtained results confirm the accuracy and efficiency of the proposed method. The numerical experiments illustrate that our MDWLS method has better performance than the traditional least squares method in case of noisy data.     

Corrective meshless particle formulations for time domain Maxwell''s equations

In this paper a meshless approximation of electromagnetic (EM) field functions and relative differential operators based on particle formulation is proposed. The idea is to obtain numerical solutions for EM problems by passing up the mesh generation usually required to compute derivatives, and by employing a set of particles arbitrarily placed in the problem domain. The meshless Smoothed Particle Hydrodynamics method has been reformulated for solving the time domain Maxwell's curl equations. The consistency of the discretized model is investigated and improvements in the approximation are obtained by modifying the numerical process. Corrective algorithms preserving meshless consistency are presented and successfully used. Test problems, dealing with even and uneven particles distribution, are simulated to validate the proposed methodology, also by introducing a comparison with analytical solution.     

Convergence order estimates of meshless collocation methods using radial basis functions

We study meshless collocation methods using radial basis functions to approximate regular solutions of systems of equations with linear differential or integral operators. Our method can be interpreted as one of the emerging meshless methods, cf. T. Belytschko et al. (1996). Its range of application is not confined to elliptic problems. However, the application to the boundary value problem for an elliptic operator, connected with an integral equation, is given as an example. Although the method has been used for special cases for about ten years, cf. E.J. Kansa (1990), there are no error bounds known. We put the main emphasis on detailed proofs of such error bounds, following the general outline described in C. Franke and R. Schaback (preprint). This revised version was published online in June 2006 with corrections to the Cover Date.