摘要

This article proposes an algebraic multigrid (AMG) approach to solve linear systems arising from applications where strong discontinuities are modeled by the extended finite element method. The application of AMG methods promises optimal scalability for solving large linear systems. However, the straightforward (or black-box') use of existing AMG techniques for extended finite element method problems is often problematic. In this paper, we highlight the reasons for this behavior and propose a relatively simple adaptation that allows one to leverage existing AMG software mostly unchanged. Numerical tests demonstrate that optimal iterative convergence rates can be attained that are comparable with AMG convergence rates associated with linear systems for standard finite element approximations without discontinuities. Published 2012. This article is a US Government work and is in the public domain in the USA.

  • 出版日期2013-4-20