A stabilized RBF collocation scheme for Neumann type boundary value problems

The numerical solutionof partial dif- ferential equations (PDEs) with Neumann bound- ary conditions (BCs) resulted from strong form collocation scheme are typically much poorer in accuracy compared to those with pure Dirichlet BCs. In this paper, we show numerically that the reason of the reduced accuracy is that Neu- mann BC requires the approximation of the spa- tial derivatives at Neumann boundaries which are significantly less accurate than approximation of main function. Therefore, we utilize boundary treatment schemes that based upon increasing the accuracy of spatial derivatives at boundaries. In- creased accuracy of the spatial derivative approx- imation can be achieved by h-refinement reduc- ing the spacing between discretization points or byincreasingthemultiquadricshape parameter, c. Increasing the MQ shape parameter is very com- putationally cost effective, but leads to increased ill-conditioning. We have implemented an im- proved version of the truncated singular value de- composition (IT-SVD) originated by Volokh and Vilnay (2000) that projects very small singular values into the null space, producing a well con- ditioned system of equations. To assess the pro- posed refinement scheme, elliptic PDEs with dif- ferent boundary conditions are analyzed. Com- parisons that made with analytical solution reveal superior accuracy and computationalefficiency of the IT-SVD solutions.

[1]  S. Timoshenko,et al.  Theory of elasticity , 1975 .

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

[3]  J. Z. Zhu,et al.  The finite element method , 1977 .

[4]  C. Micchelli.  Interpolation of scattered data: Distance matrices and conditionally positive definite functions , 1986 .

[5]  W. Madych,et al.  Multivariate interpolation and condi-tionally positive definite functions , 1988 .

[6]  R. L. Hardy.  Theory and applications of the multiquadric-biharmonic method : 20 years of discovery 1968-1988 , 1990 .

[7]  M. Buhmann.  Multivariate interpolation in odd-dimensional euclidean spaces using multiquadrics , 1990 .

[8]  W. R. Madych,et al.  Miscellaneous error bounds for multiquadric and related interpolators , 1992 .

[9]  Robert Schaback,et al.  Error estimates and condition numbers for radial basis function interpolation , 1995, Adv. Comput. Math..

[10]  Xiong Zhang,et al.  Meshless methods based on collocation with radial basis functions , 2000 .

[11]  E. Kansa,et al.  Circumventing the ill-conditioning problem with multiquadric radial basis functions: Applications to elliptic partial differential equations , 2000 .

[12]  K. Yu. Volokh,et al.  Pin-pointing solution of ill-conditioned square systems of linear equations , 2000, Appl. Math. Lett..

[13]  E. Kansa,et al.  Improved multiquadric method for elliptic partial differential equations via PDE collocation on the boundary , 2002 .

[14]  T. Driscoll,et al.  Interpolation in the limit of increasingly flat radial basis functions , 2002 .

[15]  YuanTong Gu,et al.  A meshfree method: meshfree weak–strong (MWS) form method, for 2-D solids , 2003 .

[16]  N Mai Duy,et al.  APPROXIMATION OF FUNCTION AND ITS DERIVATIVES USING RADIAL BASIS FUNCTION NETWORKS , 2003 .

[17]  B. Fornberg,et al.  Some observations regarding interpolants in the limit of flat radial basis functions , 2003 .

[18]  E. Kansa,et al.  Exponential convergence and H‐c multiquadric collocation method for partial differential equations , 2003 .

[19]  A. U.S.,et al.  Stable Computation of Multiquadric Interpolants for All Values of the Shape Parameter , 2003 .

[20]  Leevan Ling,et al.  Preconditioning for radial basis functions with domain decomposition methods , 2004, Math. Comput. Model..

[21]  Holger Wendland,et al.  Scattered Data Approximation: Conditionally positive definite functions , 2004 .

[22]  C. S. Chen,et al.  A mesh free approach using radial basis functions and parallel domain decomposition for solving three‐dimensional diffusion equations , 2004 .

[23]  D. L. Young,et al.  Solution of Maxwell's Equations Using the MQ Method , 2005 .

[24]  B. Šarler.  A Radial Basis Function Collocation Approach in Computational Fluid Dynamics , 2005 .

[25]  D. A. Shirobokov,et al.  Using radial basis functions in a ``finite difference mode'' , 2005 .

[26]  Renato Natal Jorge,et al.  Static deformations and vibration analysis of composite and sandwich plates using a layerwise theory and multiquadrics discretizations , 2005 .

[27]  J. Wertz,et al.  The role of the multiquadric shape parameters in solving elliptic partial differential equations , 2006, Comput. Math. Appl..

[28]  T. Tran-Cong,et al.  Computation of Laminated Composite Plates using Integrated Radial Basis Function Networks , 2007 .

[29]  Hojatollah Adibi,et al.  Numerical solution for biharmonic equation using multilevel radial basis functions and domain decomposition methods , 2007, Appl. Math. Comput..

[30]  Y. Hon,et al.  Geometrically Nonlinear Analysis of Reissner-Mindlin Plate by Meshless Computation , 2007 .

[31]  A. Cheng,et al.  Error estimate, optimal shape factor, and high precision computation of multiquadric collocation method , 2007 .

[32]  Jiun-Shyan Chen,et al.  Weighted radial basis collocation method for boundary value problems , 2007 .

[33]  Solving partial differential equations with point collocation and one-dimensional integrated interpolation schemes , 2007 .

[34]  E. Kansa,et al.  Stable PDE Solution Methods for Large Multiquadric Shape Parameters , 2008 .