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

Title: Two multigrid methods for three-dimensional problems with discontinuous and anisotropic coefficients

Journal Article · · SIAM J. Sci. Stat. Comput.; (United States)
DOI:https://doi.org/10.1137/0908059· OSTI ID:6134243

The subject of this paper is two multigrid methods for the numerical solution of -del . (D(x,y,z)delU(x,y,z)) +sigma(x,y,z)U(x,y,z)=F(x,y,z) in a bounded region ..cap omega.. or R/sup 3/, where D = (D/sup 1/, D/sup 2/, D/sup 3/), D/sup i/ is positive, i - 1,2,3, and D/sup i/, sigma, and F are allowed to be discontinuous across internal boundaries of ..cap omega... The first method differs in two important ways from previously considered multigrid methods for such problems; first, if standard coarsening is used, then anisotropic problems require plane relaxation to obtain a good smoothing factor; for problems with discontinuous coefficients, this was previously done with an ICCG method; here it is done with a two-dimensional multigrid method. Second, a slightly different form of interpolation is used which improves performance for nearly singular problems. Previously considered multigrid methods for problems with discontinuous coefficients have defined the coarse grid operators by Galerkin coarsening; in the second method, such coarsening is also used but with auxiliary intermediate grids obtained by semicoarsening successively in each of the three independent variables; this artifice simplifies coding considerably and results in comparable convergence.

Research Organization:
Theoretical Div., Los Alamos National Lab., Los Alamos, NM 87545
DOE Contract Number:
W-7405-ENG-36
OSTI ID:
6134243
Journal Information:
SIAM J. Sci. Stat. Comput.; (United States), Vol. 8:5
Country of Publication:
United States
Language:
English