A hybrid coupling interface method for elliptic complex interface problems |
| |
Authors: | I-Liang Chern Yu-Chen Shu |
| |
Institution: | 1. Taida Institute for Mathematical Sciences, Taipei, Taiwan;2. National Center for Theoretical Science at Taipei, Taiwan
National Taiwan University, Taipei, Taiwan
Division of Mechanics, Research Center for Applied Sciences, Adademia Sinica, Taiwan |
| |
Abstract: | We propose a coupling interface method (CIM) under Cartesian grid for solving elliptic complex interface problems in arbitrary d dimensions, where the coefficients, the source terms, and the solutions may be discontinuous or singular across the interfaces. It consists of a first-order version (CIM1) and a second-order version (CIM2). In one dimension, this finite difference method at a grid point adjacent to the interface is derived based on piecewise linear (CIM1) or quadratic (CIM2) approximation of the solution and two jump conditions. The method is extended to high dimensions through a dimensionby-dimension approach. To connect information from each dimension, a coupled equation for the principal derivatives is derived through the jump conditions in each coordinate direction. For CIM2, one-side interpolation for cross derivatives is need. This coupling approach reduces number of grid point in the finite difference stencil. The hybrid method uses CIM1 or CIM2 adaptly for complex interface. Numerical tests demonstrate that CIM1 and CIM2 are respectively first order and second order in the maximal norm with less error as compared with other methods. In addition, the hybrid CIM can solve complex interface problems in two and three dimensions. (© 2008 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim) |
| |
Keywords: | |
|
|