Commun. Comput. Phys., 15 (2014), pp. 1108-1140. |
A Kernel-Free Boundary Integral Method for Variable Coefficients Elliptic PDEs Wenjun Ying ^{1*}, Wei-Cheng Wang ^{2} 1 Department of Mathematics, MOE-LSC and Institute of Natural Sciences, Shanghai Jiao Tong University, Minhang, Shanghai 200240, P.R. China.2 Department of Mathematics, National Tsing Hua University, and National Center for Theoretical Sciences, HsinChu, 300, Taiwan. Received 17 March 2013; Accepted (in revised version) 7 November 2013 Available online 21 January 2014 doi:10.4208/cicp.170313.071113s Abstract This work proposes a generalized boundary integral method for variable coefficients elliptic partial differential equations (PDEs), including both boundary value and interface problems. The method is kernel-free in the sense that there is no need to know analytical expressions for kernels of the boundary and volume integrals in the solution of boundary integral equations. Evaluation of a boundary or volume integral is replaced with interpolation of a Cartesian grid based solution, which satisfies an equivalent discrete interface problem, while the interface problem is solved by a fast solver in the Cartesian grid. The computational work involved with the generalized boundary integral method is essentially linearly proportional to the number of grid nodes in the domain. This paper gives implementation details for a second-order version of the kernel-free boundary integral method in two space dimensions and presents numerical experiments to demonstrate the efficiency and accuracy of the method for both boundary value and interface problems. The interface problems demonstrated include those with piecewise constant and large-ratio coefficients and the heterogeneous interface problem, where the elliptic PDEs on two sides of the interface are of different types. AMS subject classifications: 35J05, 65N06, 65N38Notice: Undefined variable: pac in /var/www/html/readabs.php on line 165 Key words: Elliptic partial differential equation, variable coefficients, kernel-free boundary integral method, finite difference method, geometric multigrid iteration. *Corresponding author. Email: wying@sjtu.edu.cn (W.-J. Ying), wangwc@math.nthu.edu.tw (W.-C. Wang) |