A posteriori error estimates for discontinuous Galerkin methods using non-polynomial basis functions. Part II: Eigenvalue problems
ESAIM: Mathematical Modelling and Numerical Analysis , Tome 51 (2017) no. 5, pp. 1733-1753

Voir la notice de l'article provenant de la source Numdam

We present the first systematic work for deriving a posteriori error estimates for general non-polynomial basis functions in an interior penalty discontinuous Galerkin (DG) formulation for solving eigenvalue problems associated with second order linear operators. Eigenvalue problems of such types play important roles in scientific and engineering applications, particularly in theoretical chemistry, solid state physics and material science. Based on the framework developed in [L. Lin and B. Stamm, ESAIM: M2AN 50 (2016) 1193–1222] for second order PDEs, we develop residual type upper and lower bound error estimates for measuring the a posteriori error for eigenvalue problems. The main merit of our method is that the method is parameter-free, in the sense that all but one solution-dependent constants appearing in the upper and lower bound estimates are explicitly computable by solving local and independent eigenvalue problems, and the only non-computable constant can be reasonably approximated by a computable one without affecting the overall effectiveness of the estimates in practice. Compared to the PDE case, we find that a posteriori error estimators for eigenvalue problems must neglect certain terms, which involves explicitly the exact eigenvalues or eigenfunctions that are not accessible in numerical simulations. We define such terms carefully, and justify numerically that the neglected terms are indeed numerically high order terms compared to the computable estimators. Numerical results for a variety of problems in 1D and 2D demonstrate that both the upper bound and lower bound are effective for measuring the error of eigenvalues and eigenfunctions in the symmetric DG formulation. Our numerical results also demonstrate the sub-optimal convergence properties of eigenvalues when the non-symmetric DG formulation is used, while in such case the upper and lower bound estimators are still effective for measuring the error of eigenfunctions.

DOI : 10.1051/m2an/2016081
Classification : 65J10, 65N15, 65N30
Keywords: Discontinuous Galerkin method, a posteriori error estimation, non-polynomial basis functions, eigenvalue problems

Lin, Lin 1 ; Stamm, Benjamin 2

1 Department of Mathematics, University of California, Berkeley and Computational Research Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA.
2 Center for Computational Engineering, Mathematics Department, RWTH Aachen University, Schinkelstr. 2, 52062 Aachen, Germany and Computational Biomedicine, Institute for Advanced Simulation IAS-5 and Institute of Neuroscience and Medicine INM-9, Forschungszentrum Jülich, Germany
@article{M2AN_2017__51_5_1733_0,
     author = {Lin, Lin and Stamm, Benjamin},
     title = {A posteriori error estimates for discontinuous {Galerkin} methods using non-polynomial basis functions. {Part} {II:} {Eigenvalue} problems},
     journal = {ESAIM: Mathematical Modelling and Numerical Analysis },
     pages = {1733--1753},
     publisher = {EDP-Sciences},
     volume = {51},
     number = {5},
     year = {2017},
     doi = {10.1051/m2an/2016081},
     mrnumber = {3731547},
     zbl = {1384.65080},
     language = {en},
     url = {http://geodesic.mathdoc.fr/articles/10.1051/m2an/2016081/}
}
TY  - JOUR
AU  - Lin, Lin
AU  - Stamm, Benjamin
TI  - A posteriori error estimates for discontinuous Galerkin methods using non-polynomial basis functions. Part II: Eigenvalue problems
JO  - ESAIM: Mathematical Modelling and Numerical Analysis 
PY  - 2017
SP  - 1733
EP  - 1753
VL  - 51
IS  - 5
PB  - EDP-Sciences
UR  - http://geodesic.mathdoc.fr/articles/10.1051/m2an/2016081/
DO  - 10.1051/m2an/2016081
LA  - en
ID  - M2AN_2017__51_5_1733_0
ER  - 
%0 Journal Article
%A Lin, Lin
%A Stamm, Benjamin
%T A posteriori error estimates for discontinuous Galerkin methods using non-polynomial basis functions. Part II: Eigenvalue problems
%J ESAIM: Mathematical Modelling and Numerical Analysis 
%D 2017
%P 1733-1753
%V 51
%N 5
%I EDP-Sciences
%U http://geodesic.mathdoc.fr/articles/10.1051/m2an/2016081/
%R 10.1051/m2an/2016081
%G en
%F M2AN_2017__51_5_1733_0
Lin, Lin; Stamm, Benjamin. A posteriori error estimates for discontinuous Galerkin methods using non-polynomial basis functions. Part II: Eigenvalue problems. ESAIM: Mathematical Modelling and Numerical Analysis , Tome 51 (2017) no. 5, pp. 1733-1753. doi: 10.1051/m2an/2016081

Cité par Sources :