Optimal variable shape parameter for multiquadric based RBF-FD method

In this follow up paper to our previous study in Bayona et al. (2011) [2], we present a new technique to compute the solution of PDEs with the multiquadric based RBF finite difference method (RBF-FD) using an optimal node dependent variable value of the shape parameter. This optimal value is chosen so that, to leading order, the local approximation error of the RBF-FD formulas is zero. In our previous paper (Bayona et al., 2011) [2] we considered the case of an optimal (constant) value of the shape parameter for all the nodes. Our new results show that, if one allows the shape parameter to be different at each grid point of the domain, one may obtain very significant accuracy improvements with a simple and inexpensive numerical technique. We analyze the same examples studied in Bayona et al. (2011) [2], both with structured and unstructured grids, and compare our new results with those obtained previously. We also find that, if there are a significant number of nodes for which no optimal value of the shape parameter exists, then the improvement in accuracy deteriorates significantly. In those cases, we use generalized multiquadrics as RBFs and choose the exponent of the multiquadric at each node to assure the existence of an optimal variable shape parameter.

[1]  Gregory E. Fasshauer,et al.  Analysis of natural frequencies of composite plates by an RBF-pseudospectral method , 2007 .

[2]  Guirong Liu,et al.  A point interpolation meshless method based on radial basis functions , 2002 .

[3]  Manuel Kindelan,et al.  Optimal constant shape parameter for multiquadric based RBF-FD method , 2011, J. Comput. Phys..

[4]  Benjamin Seibold,et al.  Minimal positive stencils in meshfree finite difference methods for the Poisson equation , 2008, 0802.2674.

[5]  Robert Vertnik,et al.  Meshfree explicit local radial basis function collocation method for diffusion problems , 2006, Comput. Math. Appl..

[6]  C. Shu,et al.  Multiquadric Finite Difference (MQ-FD) Method and its Application , 2009 .

[7]  Y. V. S. S. Sanyasiraju,et al.  Local radial basis function based gridfree scheme for unsteady incompressible viscous flows , 2008, J. Comput. Phys..

[8]  Bengt Fornberg,et al.  Scattered node compact finite difference-type formulas generated from radial basis functions , 2006, J. Comput. Phys..

[9]  B. Fornberg,et al.  Radial basis function interpolation: numerical and analytical developments , 2003 .

[10]  S. Chantasiriwan.  Investigation of the use of radial basis functions in local collocation method for solving diffusion problems , 2004 .

[11]  Renato Natal Jorge,et al.  Analysis of Functionally Graded Plates by a Robust Meshless Method , 2007 .

[12]  Oleg Davydov,et al.  Adaptive meshless centres and RBF stencils for Poisson equation , 2011, J. Comput. Phys..

[13]  António J.M. Ferreira,et al.  A local radial basis functions—Finite differences technique for the analysis of composite plates , 2011 .

[14]  Manuel Kindelan,et al.  RBF-FD formulas and convergence properties , 2010, J. Comput. Phys..

[15]  E. Kansa.  MULTIQUADRICS--A SCATTERED DATA APPROXIMATION SCHEME WITH APPLICATIONS TO COMPUTATIONAL FLUID-DYNAMICS-- II SOLUTIONS TO PARABOLIC, HYPERBOLIC AND ELLIPTIC PARTIAL DIFFERENTIAL EQUATIONS , 1990 .

[16]  H. Ding,et al.  Error estimates of local multiquadric‐based differential quadrature (LMQDQ) method through numerical experiments , 2005 .

[17]  Manuel Kindelan,et al.  Use of singularity capturing functions in the solution of problems with discontinuous boundary conditions , 2009 .

[18]  B. Fornberg,et al.  A numerical study of some radial basis function based solution methods for elliptic PDEs , 2003 .

[19]  David Stevens,et al.  The use of PDE centres in the local RBF Hermitian method for 3D convective-diffusion problems , 2009, J. Comput. Phys..

[20]  S. C. Fan,et al.  Local multiquadric approximation for solving boundary value problems , 2003 .

[21]  Bengt Fornberg,et al.  Stabilization of RBF-generated finite difference methods for convective PDEs , 2011, J. Comput. Phys..

[22]  J. Halton.  On the efficiency of certain quasi-random sequences of points in evaluating multi-dimensional integrals , 1960 .

[23]  E. Kansa.  Multiquadrics—A scattered data approximation scheme with applications to computational fluid-dynamics—I surface approximations and partial derivative estimates , 1990 .

[24]  A. I. Tolstykh,et al.  On using radial basis functions in a “finite difference mode” with applications to elasticity problems , 2003 .

[25]  R. L. Hardy.  Multiquadric equations of topography and other irregular surfaces , 1971 .

[26]  C. D. Gelatt,et al.  Optimization by Simulated Annealing , 1983, Science.

[27]  C. Shu,et al.  Local radial basis function-based differential quadrature method and its application to solve two-dimensional incompressible Navier–Stokes equations , 2003 .

[28]  Y. Sanyasiraju,et al.  Local RBF‐FD solutions for steady convection–diffusion problems , 2007 .