TY - JOUR

T1 - Adaptive Finite Element Methods for Elliptic Problems with Discontinuous Coefficients

AU - Bonito, Andrea

AU - DeVore, Ronald A.

AU - Nochetto, Ricardo H.

N1 - KAUST Repository Item: Exported on 2020-10-01
Acknowledged KAUST grant number(s): KUS-C1-016-04
Acknowledgements: Department of Mathematics, Texas A&M University, College Station, TX 77843 (bonito@math.tamu.edu, rdevore@math.tamu.edu). The first author was partially supported by NSF grant DMS-1254618 and ONR grant N000141110712. The second author was partially supported by ONR grants N00014-12-1-0561 and N00014-11-1-0712, NSF grant DMS-12227151, and award KUS-C1-016-04 made by King Abdullah University of Science and Technology (KAUST).Department of Mathematics and Institute for Physical Science and Technology, University of Maryland, College Park, MD 20742 (rhn@math.umd.edu). This author was partially supported by NSF grant DMS-1109325.
This publication acknowledges KAUST support, but has no KAUST affiliated authors.

PY - 2013/1

Y1 - 2013/1

N2 - Elliptic PDEs with discontinuous diffusion coefficients occur in application domains such as diffusions through porous media, electromagnetic field propagation on heterogeneous media, and diffusion processes on rough surfaces. The standard approach to numerically treating such problems using finite element methods is to assume that the discontinuities lie on the boundaries of the cells in the initial triangulation. However, this does not match applications where discontinuities occur on curves, surfaces, or manifolds, and could even be unknown beforehand. One of the obstacles to treating such discontinuity problems is that the usual perturbation theory for elliptic PDEs assumes bounds for the distortion of the coefficients in the L∞ norm and this in turn requires that the discontinuities are matched exactly when the coefficients are approximated. We present a new approach based on distortion of the coefficients in an Lq norm with q < ∞ which therefore does not require the exact matching of the discontinuities. We then use this new distortion theory to formulate new adaptive finite element methods (AFEMs) for such discontinuity problems. We show that such AFEMs are optimal in the sense of distortion versus number of computations, and report insightful numerical results supporting our analysis. © 2013 Societ y for Industrial and Applied Mathematics.

AB - Elliptic PDEs with discontinuous diffusion coefficients occur in application domains such as diffusions through porous media, electromagnetic field propagation on heterogeneous media, and diffusion processes on rough surfaces. The standard approach to numerically treating such problems using finite element methods is to assume that the discontinuities lie on the boundaries of the cells in the initial triangulation. However, this does not match applications where discontinuities occur on curves, surfaces, or manifolds, and could even be unknown beforehand. One of the obstacles to treating such discontinuity problems is that the usual perturbation theory for elliptic PDEs assumes bounds for the distortion of the coefficients in the L∞ norm and this in turn requires that the discontinuities are matched exactly when the coefficients are approximated. We present a new approach based on distortion of the coefficients in an Lq norm with q < ∞ which therefore does not require the exact matching of the discontinuities. We then use this new distortion theory to formulate new adaptive finite element methods (AFEMs) for such discontinuity problems. We show that such AFEMs are optimal in the sense of distortion versus number of computations, and report insightful numerical results supporting our analysis. © 2013 Societ y for Industrial and Applied Mathematics.

UR - http://hdl.handle.net/10754/597455

UR - http://epubs.siam.org/doi/10.1137/130905757

UR - http://www.scopus.com/inward/record.url?scp=84892946078&partnerID=8YFLogxK

U2 - 10.1137/130905757

DO - 10.1137/130905757

M3 - Article

SN - 0036-1429

VL - 51

SP - 3106

EP - 3134

JO - SIAM Journal on Numerical Analysis

JF - SIAM Journal on Numerical Analysis

IS - 6

ER -