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-refmement 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.
N. A. Libre et al., "A Stabilized RBF Collocation Scheme for Neumann Type Boundary Value Problems," CMES - Computer Modeling in Engineering and Sciences, vol. 24, no. 1, pp. 61-80, Tech Science Press, Jan 2008.
Civil, Architectural and Environmental Engineering
Keywords and Phrases
Boundary value problem; Improved truncated singular value decomposition; Meshfree; Multiquadric; Neumann condition; RBF collocation; Stability; Strong form
International Standard Serial Number (ISSN)
Article - Journal
© 2008 Tech Science Press, All rights reserved.
01 Jan 2008