NDLI: Unstructured Additive Schwarz--Conjugate Gradient Method for Elliptic Problems with Highly Discontinuous Coefficients

Content Provider	Society for Industrial and Applied Mathematics (SIAM)
Author	Graham, I. G. Hagger, M. J.
Copyright Year	1999
Abstract	This paper concerns the iterative solution of symmetric elliptic problems with piecewise constant coefficients in two or three space dimensions discretized by linear finite element methods on unstructured triangular or tetrahedral meshes. The effect of the discontinuous coefficients is studied by first postulating that there are d fixed regions of the domain where the coefficient takes constant positive values a = (a1 , . . . , ad )and then considering certain (positive) coefficient sequences {a(m) } in which some of these values approach 0 or $\infty$ as $m \rightarrow \infty$. We consider the performance of additive Schwarz domain decomposition preconditioners constructed from local solves on automatically generated subdomains together with a global solve on some coarser grid. Assuming no relationship between the regions on which the coefficient function is constant and either the subdomains or the coarse grid, we show that the preconditioned conjugate gradient method converges (in both the energy and the Euclidean norms) with a number of iterations which grows only logarithmically in the size of the maximum jump ${\cal J}^{(m)} := \max\{a_k^{(m)}/a_l^{(m)} \ : \ k, l = 1, \ldots , d\} $, as $m\rightarrow \infty$. The result is obtained by a careful analysis of the preconditioned matrix. It cannot be obtained by the usual procedure of estimating condition numbers: a simple example is given in which the condition number of both the original and the preconditioned stiffness matrices degrade linearly in ${\cal J}^{(m)}$. Recent results of Chan, Smith, and Zou have shown that by using this preconditioner together with the conjugate gradient method, the number of iterations can be bounded independently of the mesh diameter provided the subdomains have overlap commensurate with the size of the coarse mesh. Our results now show that this method is also highly resilient to discontinuous coefficients, even if no attention is paid to the coefficient discontinuity in the construction of the solver.
Starting Page	2041
Ending Page	2066
Page Count	26
File Format	PDF
ISSN	10648275
DOI	10.1137/S1064827596305593
e-ISSN	10957197
Issue Number	6
Volume Number	20
Language	English
Publisher	Society for Industrial and Applied Mathematics
Publisher Date	2006-07-25
Access Restriction	Subscribed
Subject Keyword	discontinuous coefficients Iterative methods for linear systems domain decomposition overlapping subdomains diagonal scaling additive Schwarz conjugate gradient method Finite elements, Rayleigh-Ritz and Galerkin methods, finite methods
Content Type	Text
Resource Type	Article
Subject	Applied Mathematics Computational Mathematics

Sl.	Authority	Responsibilities	Communication Details
1	Ministry of Education (GoI), Department of Higher Education	Sanctioning Authority	https://www.education.gov.in/ict-initiatives
2	Indian Institute of Technology Kharagpur	Host Institute of the Project: The host institute of the project is responsible for providing infrastructure support and hosting the project	https://www.iitkgp.ac.in
3	National Digital Library of India Office, Indian Institute of Technology Kharagpur	The administrative and infrastructural headquarters of the project	Dr. B. Sutradhar bsutra@ndl.gov.in
4	Project PI / Joint PI	Principal Investigator and Joint Principal Investigators of the project	Dr. B. Sutradhar bsutra@ndl.gov.in Prof. Saswat Chakrabarti will be added soon
5	Website/Portal (Helpdesk)	Queries regarding NDLI and its services	support@ndl.gov.in
6	Contents and Copyright Issues	Queries related to content curation and copyright issues	content@ndl.gov.in
7	National Digital Library of India Club (NDLI Club)	Queries related to NDLI Club formation, support, user awareness program, seminar/symposium, collaboration, social media, promotion, and outreach	clubsupport@ndl.gov.in
8	Digital Preservation Centre (DPC)	Assistance with digitizing and archiving copyright-free printed books	dpc@ndl.gov.in
9	IDR Setup or Support	Queries related to establishment and support of Institutional Digital Repository (IDR) and IDR workshops	idr@ndl.gov.in

Two-Level Overlapping Schwarz Algorithms for a Staggered Discontinuous Galerkin Method

Schwarz domain decomposition preconditioners for discontinuous Galerkin approximations of elliptic problems: non-overlapping case

Additive Average Schwarz Method for Elliptic Mortar Finite Element Problems with Highly Heterogeneous Coefficients

TwoLevel Schwarz Algorithms with Overlapping Subregions for Mortar Finite Elements

Multiplicative Schwarz Methods for Discontinuous Galerkin Approximations of Elliptic Problems

Domain Decomposition for Less Regular Subdomains: Overlapping Schwarz in Two Dimensions

An Alternative Coarse Space for Irregular Subdomains and an Overlapping Schwarz Algorithm for Scalar Elliptic Problems in the Plane

An Optimal Domain Decomposition Preconditioner for the Finite Element Solution of Linear Elasticity Problems

An Overlapping Schwarz Algorithm for Raviart--Thomas Vector Fields with Discontinuous Coefficients

Unstructured Additive Schwarz--Conjugate Gradient Method for Elliptic Problems with Highly Discontinuous Coefficients

Similar Documents

Two-Level Overlapping Schwarz Algorithms for a Staggered Discontinuous Galerkin Method

Schwarz domain decomposition preconditioners for discontinuous Galerkin approximations of elliptic problems: non-overlapping case

Additive Average Schwarz Method for Elliptic Mortar Finite Element Problems with Highly Heterogeneous Coefficients

TwoLevel Schwarz Algorithms with Overlapping Subregions for Mortar Finite Elements

Multiplicative Schwarz Methods for Discontinuous Galerkin Approximations of Elliptic Problems

Domain Decomposition for Less Regular Subdomains: Overlapping Schwarz in Two Dimensions

An Alternative Coarse Space for Irregular Subdomains and an Overlapping Schwarz Algorithm for Scalar Elliptic Problems in the Plane

An Optimal Domain Decomposition Preconditioner for the Finite Element Solution of Linear Elasticity Problems

An Overlapping Schwarz Algorithm for Raviart--Thomas Vector Fields with Discontinuous Coefficients

Unstructured Additive Schwarz--Conjugate Gradient Method for Elliptic Problems with Highly Discontinuous Coefficients