Premium
Energy‐minimizing coarse spaces for two‐level Schwarz methods for multiscale PDEs
Author(s) -
Van lent Jan,
Scheichl Robert,
Graham Ivan G.
Publication year - 2009
Publication title -
numerical linear algebra with applications
Language(s) - English
Resource type - Journals
SCImago Journal Rank - 1.02
H-Index - 53
eISSN - 1099-1506
pISSN - 1070-5325
DOI - 10.1002/nla.641
Subject(s) - preconditioner , mathematics , basis function , basis (linear algebra) , partial differential equation , additive schwarz method , simple (philosophy) , schwarz alternating method , finite element method , partition (number theory) , space (punctuation) , mathematical optimization , multigrid method , domain decomposition methods , elliptic partial differential equation , function space , computer science , mathematical analysis , geometry , iterative method , philosophy , physics , epistemology , combinatorics , thermodynamics , operating system
Two‐level overlapping Schwarz methods for elliptic partial differential equations combine local solves on overlapping domains with a global solve of a coarse approximation of the original problem. To obtain robust methods for equations with highly varying coefficients, it is important to carefully choose the coarse approximation. Recent theoretical results by the authors have shown that bases for such robust coarse spaces should be constructed such that the energy of the basis functions is minimized. We give a simple derivation of a method that finds such a minimum energy basis using one local solve per coarse space basis function and one global solve to enforce a partition of unity constraint. Although this global solve may seem prohibitively expensive, we demonstrate that a one‐level overlapping Schwarz method is an effective and scalable preconditioner and we show that such a preconditioner can be implemented efficiently using the Sherman–Morrison–Woodbury formula. The result is an elegant, scalable, algebraic method for constructing a robust coarse space given only the supports of the coarse space basis functions. Numerical experiments on a simple two‐dimensional model problem with a variety of binary and multiscale coefficients confirm this. Numerical experiments also show that, when used in a two‐level preconditioner, the energy‐minimizing coarse space gives better results than other coarse space constructions, such as the multiscale finite element approach. Copyright © 2009 John Wiley & Sons, Ltd.