On the subdomain‐Galerkin/least squares method for 2‐ and 3‐D mixed elliptic problems with reaction terms

Abstract

In this article we apply the subdomain-Galerkin/least squares method, which is first proposed by Chang and Gunzburger for first-order elliptic systems without reaction terms in the plane, to solve second-order non-selfadjoint elliptic problems in two-and three-dimensional bounded domains with triangular or tetrahedral regular triangulations. This method can be viewed as a combination of a direct cell vertex finite volume discretization step and an algebraic least-squares minimization step in which the pressure is approximated by piecewise linear elements and the flux by the lowest order Raviart-Thomas space. This combined approach has the advantages of both finite volume and least-squares methods. Among other things, the combined method is not subject to the Ladyzhenskaya-Babuska-Brezzi condition, and the resulting linear system is symmetric and positive definite. An optimal error estimate in the H 1 (⍀) ϫ H(div; ⍀) norm is derived. An equivalent residual-type a posteriori error estimator is also given.