New Geometric Immersed Interface Multigrid Solvers
New Geometric Immersed Interface Multigrid Solvers
复制标题
DOI:
10.1137/s1064827503421707
复制
发表时间:
2004-05
期刊:
影响因子:
--
通讯作者:
L. Adams;T. Chartier
中科院分区:
文献类型:
--
作者:
L. Adams;T. Chartier
In [L. Adams and Z. Li, SIAM J. Sci. Comput., 24 (2002), pp. 463--479], a multigrid method was designed specifically for interface problems that have been discretized using the methods described therein, and in [Z. Li and K. Ito, SIAM J. Sci. Comput., 23 (2001), pp. 339--361] for elliptic interface problems using the maximum principle preserving schemes. In this paper, we improve on this method by giving a new interpolator for grid points near the immersed interface and a new restrictor that guarantees the coarse-grid matrices are M-matrices. We compare this new restrictor to injection and the transpose of interpolation. We show that the number of V-cycles is constant as the mesh size decreases and increases only slightly as the ratio of the discontinuous problem coefficient grows at the interface only when this new restrictor is used.