skip to main content
OSTI.GOV title logo U.S. Department of Energy
Office of Scientific and Technical Information

Title: Multilevel Methods for Elliptic Problems with Highly Varying Coefficients on Nonaligned Coarse Grids

Abstract

We generalize the analysis of classical multigrid and two-level overlapping Schwarz methods for 2nd order elliptic boundary value problems to problems with large discontinuities in the coefficients that are not resolved by the coarse grids or the subdomain partition. The theoretical results provide a recipe for designing hierarchies of standard piecewise linear coarse spaces such that the multigrid convergence rate and the condition number of the Schwarz preconditioned system do not depend on the coefficient variation or on any mesh parameters. One assumption we have to make is that the coarse grids are sufficiently fine in the vicinity of cross points or where regions with large diffusion coefficients are separated by a narrow region where the coefficient is small. We do not need to align them with possible discontinuities in the coefficients. The proofs make use of novel stable splittings based on weighted quasi-interpolants and weighted Poincaré-type inequalities. Finally, numerical experiments are included that illustrate the sharpness of the theoretical bounds and the necessity of the technical assumptions.

Authors:
 [1];  [2];  [3]
  1. Univ. of Bath (United Kingdom). Dept. of Mathematical Sciences
  2. Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
  3. Pennsylvania State Univ., University Park, PA (United States). Dept. of Mathematics
Publication Date:
Research Org.:
Lawrence Livermore National Lab. (LLNL), Livermore, CA (United States)
Sponsoring Org.:
USDOE
OSTI Identifier:
1226946
Report Number(s):
LLNL-JRNL-451252
Journal ID: ISSN 0036-1429
DOE Contract Number:  
AC52-07NA27344
Resource Type:
Journal Article
Journal Name:
SIAM Journal on Numerical Analysis
Additional Journal Information:
Journal Volume: 50; Journal Issue: 3; Journal ID: ISSN 0036-1429
Publisher:
Society for Industrial and Applied Mathematics
Country of Publication:
United States
Language:
English
Subject:
97 MATHEMATICS, COMPUTING, AND INFORMATION SCIENCE; coarse spaces; multigrid; overlapping Schwarz method; large coefficient jumpts

Citation Formats

Scheichl, Robert, Vassilevski, Panayot S., and Zikatanov, Ludmil T. Multilevel Methods for Elliptic Problems with Highly Varying Coefficients on Nonaligned Coarse Grids. United States: N. p., 2012. Web. doi:10.1137/100805248.
Scheichl, Robert, Vassilevski, Panayot S., & Zikatanov, Ludmil T. Multilevel Methods for Elliptic Problems with Highly Varying Coefficients on Nonaligned Coarse Grids. United States. doi:10.1137/100805248.
Scheichl, Robert, Vassilevski, Panayot S., and Zikatanov, Ludmil T. Thu . "Multilevel Methods for Elliptic Problems with Highly Varying Coefficients on Nonaligned Coarse Grids". United States. doi:10.1137/100805248. https://www.osti.gov/servlets/purl/1226946.
@article{osti_1226946,
title = {Multilevel Methods for Elliptic Problems with Highly Varying Coefficients on Nonaligned Coarse Grids},
author = {Scheichl, Robert and Vassilevski, Panayot S. and Zikatanov, Ludmil T.},
abstractNote = {We generalize the analysis of classical multigrid and two-level overlapping Schwarz methods for 2nd order elliptic boundary value problems to problems with large discontinuities in the coefficients that are not resolved by the coarse grids or the subdomain partition. The theoretical results provide a recipe for designing hierarchies of standard piecewise linear coarse spaces such that the multigrid convergence rate and the condition number of the Schwarz preconditioned system do not depend on the coefficient variation or on any mesh parameters. One assumption we have to make is that the coarse grids are sufficiently fine in the vicinity of cross points or where regions with large diffusion coefficients are separated by a narrow region where the coefficient is small. We do not need to align them with possible discontinuities in the coefficients. The proofs make use of novel stable splittings based on weighted quasi-interpolants and weighted Poincaré-type inequalities. Finally, numerical experiments are included that illustrate the sharpness of the theoretical bounds and the necessity of the technical assumptions.},
doi = {10.1137/100805248},
journal = {SIAM Journal on Numerical Analysis},
issn = {0036-1429},
number = 3,
volume = 50,
place = {United States},
year = {2012},
month = {6}
}