Inverse-average-type finite element discretizations of selfadjoint second-order elliptic problems

Peter Markowich*, Miloá A. Zlámal

*Corresponding author for this work

Research output: Contribution to journalArticlepeer-review

21 Scopus citations

Abstract

This paper is concerned with the analysis of a class of "special purpose" piecewise linear finite element discretizations of selfadjoint second-order elliptic boundary value problems. The discretization differs from standard finite element methods by inverse-average-type approximations (along element sides) of the coefficient function a(x) in the operator — div(a(x) grad u). The derivation of the discretization is based on approximating the flux density J = a grad u by constants on each element. In many cases the flux density is well behaved (moderately varying) even if a(x) and u(x) are fast varying. Discretization methods of this type have been used successfully in semiconductor device simulation for many years; however, except in the one-dimensional case, the mathematical understanding of these methods was rather limited. We analyze the stiffness matrix and prove that—under a rather mild restriction on the mesh—it is a diagonally dominant Stieltjes matrix. Most importantly, we derive an estimate which asserts that the piecewise linear interpolant of the solution u is approximated to order 1 by the finite element solution in the //'-norm. The estimate depends only on the mesh width and on derivatives of the flux density and of a possibly occurring inhomogeneity.

Original languageEnglish (US)
Pages (from-to)431-449
Number of pages19
JournalMathematics of Computation
Volume51
Issue number184
DOIs
StatePublished - Jan 1 1988

ASJC Scopus subject areas

  • Algebra and Number Theory
  • Computational Mathematics
  • Applied Mathematics

Fingerprint

Dive into the research topics of 'Inverse-average-type finite element discretizations of selfadjoint second-order elliptic problems'. Together they form a unique fingerprint.

Cite this