English

Stability of mixed FEMs for non-selfadjoint indefinite second-order linear elliptic PDEs

Numerical Analysis 2022-03-10 v2 Numerical Analysis

Abstract

For a well-posed non-selfadjoint indefinite second-order linear elliptic PDE with general coefficients A,b,γ\mathbf A, \mathbf b,\gamma in LL^\infty and symmetric and uniformly positive definite coefficient matrix A\mathbf A, this paper proves that mixed finite element problems are uniquely solvable and the discrete solutions are uniformly bounded, whenever the underlying shape-regular triangulation is sufficiently fine. This applies to the Raviart-Thomas (RT) and Brezzi-Douglas-Marini (BDM) finite element families of any order and in any space dimension and leads to the best-approximation estimate in H(div)×L2H(div)\times L^2 as well as in in L2×L2L^2\times L^2 up to oscillations. This generalises earlier contributions for piecewise Lipschitz continuous coefficients to LL^\infty coefficients. The compactness argument of Schatz and Wang for the displacement-oriented problem does not apply immediately to the mixed formulation in H(div)×L2H(div)\times L^2. But it allows the uniform approximation of some L2L^2 contributions and can be combined with a recent L2L^2 best-approximation result from the medius analysis. This technique circumvents any regularity assumption and the application of a Fortin interpolation operator.

Keywords

Cite

@article{arxiv.2203.02716,
  title  = {Stability of mixed FEMs for non-selfadjoint indefinite second-order linear elliptic PDEs},
  author = {C. Carstensen and Neela Nataraj and Amiya K. Pani},
  journal= {arXiv preprint arXiv:2203.02716},
  year   = {2022}
}

Comments

Accepted for publication in Numerische Mathematik