An improved non-traditional finite element formulation for solving three-dimensional elliptic interface problems

详细信息    查看全文

• 作者：Liqun Wang^a ; ^{wliqunhmily@gmail.com} ; Songming Hou^b ; ^{shou@latech.edu} ; Liwei Shi^c ; ^{sliweihmily@gmail.com}
• 关键词：Non-traditional finite element method ; Elliptic equation ; Three-Dimensional interface problems ; Jump condition ; Matrix coefficient
• 刊名：Computers & Mathematics with Applications
• 出版年：2017
• 出版时间：1 February 2017
• 年：2017
• 卷：73
• 期：3
• 页码：374-384
• 全文大小：1361 K
• 卷排序：73

文摘

Solving elliptic equations with interfaces has wide applications in engineering and science. The real world problems are mostly in three dimensions, while an efficient and accurate solver is a challenge. Some existing methods that work well in two dimensions are too complicated to be generalized to three dimensions. Although traditional finite element method using body-fitted grid is well-established, the expensive cost of mesh generation is an issue. In this paper, an efficient non-traditional finite element method with non-body-fitted grids is proposed to solve elliptic interface problems. The special cases when the interface cuts though grid points are handled carefully, rather than perturbing the cutting point away to apply the method for general case. Both Dirichlet and Neumann boundary conditions are considered. Numerical experiments show that this method is approximately second order accurate in the L∞L∞ norm and L2L2 norm for piecewise smooth solutions. The large sparse matrix for our linear system also has nice structure and properties.