Exponential convergence through linear finite element discretization of stratified subdomains

详细信息    查看全文

• 作者：Murthy N. Guddati^a ; ^{murthy.guddati@ncsu.edu" class="auth_mail" title="E-mail the corresponding author} ; Vladimir Druskin^b ; Ali Vaziri Astaneh^a
• 关键词：Optimal grids ; Padé ; approximants ; Dirichlet-to-Neumann maps ; Rational approximation ; Spectral element methods
• 刊名：Journal of Computational Physics
• 出版年：2016
• 出版时间：1 October 2016
• 年：2016
• 卷：322
• 期：Complete
• 页码：429-447
• 全文大小：789 K

文摘

Motivated by problems where the response is needed at select localized regions in a large computational domain, we devise a novel finite element discretization that results in exponential convergence at pre-selected points. The key features of the discretization are (a) use of midpoint integration to evaluate the contribution matrices, and (b) an unconventional mapping of the mesh into complex space. Named complex-length finite element method (CFEM), the technique is linked to Padé approximants that provide exponential convergence of the Dirichlet-to-Neumann maps and thus the solution at specified points in the domain. Exponential convergence facilitates drastic reduction in the number of elements. This, combined with sparse computation associated with linear finite elements, results in significant reduction in the computational cost. The paper presents the basic ideas of the method as well as illustration of its effectiveness for a variety of problems involving Laplace, Helmholtz and elastodynamics equations.