Open Access

ARTICLE

A stabilized RBF collocation scheme for Neumann type boundary value problems

by Nicolas Ali Libre¹, Arezoo Emdadi², Edward J. Kansa³, Mohammad Rahimian², MohammadShekarchi²

1 Corresponding author, nalibre@ut.ac.ir
2 Department of civil Engineering, University of Tehran, Tehran, Iran
3 Corresponding author, ejkansa@ucdavis.edu
4 Department of Mechanical and Aeronautical Engineering, University of California, Davis, CA 95616 USA

Computer Modeling in Engineering & Sciences 2008, 24(1), 61-80. https://doi.org/10.3970/cmes.2008.024.061

Abstract

The numerical solution of partial differential equations (PDEs) with Neumann boundary 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 Neumann BC requires the approximation of the spatial 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. Increased accuracy of the spatial derivative approximation can be achieved by h-refinement reducing the spacing between discretization points or by increasing the multiquadric shape parameter, c. Increasing the MQ shape parameter is very computationally cost effective, but leads to increased ill-conditioning. We have implemented an improved version of the truncated singular value decomposition (IT-SVD) originated by Volokh and Vilnay (2000) that projects very small singular values into the null space, producing a well conditioned system of equations. To assess the proposed refinement scheme, elliptic PDEs with different boundary conditions are analyzed. Comparisons that made with analytical solution reveal superior accuracy and computational efficiency of the IT-SVD solutions.

---

Keywords

---