An Augmented Two-Scale Finite Element Method for Eigenvalue Problems

Abstract

In this paper, an augmented two-scale finite element method is proposed for a class of linear and nonlinear eigenvalue problems on tensor-product domains. Through a correction step, the augmented two-scale finite element solution is obtained by solving an eigenvalue problem on a low-dimensional augmented subspace. Theoretical analysis and numerical experiments show that the augmented two-scale finite element solution achieves the same order of accuracy as the standard finite element solution on a fine grid, but the computational cost required by the former solution is much lower than that demanded by the latter. The augmented two-scale finite element method also improves the approximation accuracy of eigenfunctions in the $L^2(\varOmega )$ norm compared with the two-scale finite element method.

Keywords

Two-scale / Finite element / Augmented subspace method / Eigenvalue problem / Partial differential equation / 65N15 / 65N25 / 65N30 / 65N50

References

[1] Babuska I, Osborn JE. Finite element-Galerkin approximation of the eigenvalues and eigenvectors of selfadjoint problems. Math. Comput.. 1989, 52(186): 275-297

[2] Babuska, I., Osborn, J.E.: Eigenvalue problems. In: Handbook of Numerical Analysis, vol. II, pp. 641–787. North-Holland, Amsterdam (1991)

[3] Bao, W.: The nonlinear Schrödinger equation and applications in Bose-Einstein condensation and plasma physics. Master Review, Lecture Note Series, vol. 9. IMS, NUS (2007)

[4] Bao W, Du Q. Computing the ground state solution of Bose-Einstein condensates by a normalized gradient flow. SIAM J. Sci. Comput.. 2004, 25: 1674-1697

[5] Cancès E, Chakir R, He L, Maday Y. Two-grid methods for a class of nonlinear elliptic eigenvalue problems. IMA J. Numer. Anal.. 2018, 38: 605-645

[6] Cancès E, Chakir R, Maday Y. Numerical analysis of nonlinear eigenvalue problems. J. Sci. Comput.. 2010, 45: 90-117

[7] Cancès E, Chakir R, Maday Y. Numerical analysis of the planewave discretization of some orbital-free and Kohn-Sham models. Math. Model. Numer. Anal.. 2012, 46: 341-388

[8] Cancès E, Defranceschi M, Kutzelnigg W, Le Bris C, Maday Y. Ciarlet PG, Le Bris C. Computational quantum chemistry: a primer. Handbook of Numerical Analysis, Volume X: Special Volume: Computational Chemistry. 2003North-Holland, Amsterdam: 3270

[9] Chen H, Gong X, He L, Yang Z, Zhou A. Numerical analysis of finite dimensional approximations of Kohn-Sham models. Adv. Comput. Math.. 2013, 38: 225-256

[10] Chen H, Gong X, Zhou A. Numerical approximations of a nonlinear eigenvalue problem and applications to a density functional model. Math. Methods Appl. Sci.. 2010, 33: 1723-1742

[11] Chen H, He L, Zhou A. Finite element approximations of nonlinear eigenvalue problems in quantum physics. Comput. Methods Appl. Mech. Eng.. 2011, 200: 1846-1865

[12] Chen H, Liu F, Zhou A. A two-scale higher-order finite element discretization for Schrödinger equation. J. Comput. Math.. 2009, 27: 315-337

[13] Ciarlet PG. The Finite Element Method for Elliptic Problems. 1978, Amsterdam, North-Holland

[14] Dai X, Zhou A. Three-scale finite element discretizations for quantum eigenvalue problems. SIAM J. Numer. Anal.. 2008, 46: 295-324

[15] Dauge M. Elliptic boundary value problems on corner domains. Lecture Notes in Mathematics. 1988, Berlin, Springer1341

[16] Gao X, Liu F, Zhou A. Three-scale finite element eigenvalue discretizations. BIT. 2008, 48(3): 533-562

[17] Gong X, Shen L, Zhou A. Finite element approximations for Schrödinger equations with applications to electronic structure computations. J. Comput. Math.. 2008, 26: 1-14

[18] Hou P, Liu F. Two-scale finite element discretizations for nonlinear eigenvalue problems in quantum physics. Adv. Comput. Math.. 2021, 47: 59

[19] Hu G, Xie H, Xu F. A multilevel correction adaptive finite element method for Kohn-Sham equation. J. Comput. Phys.. 2018, 355: 436-449

[20] Jia S, Xie H, Xie M, Xu F. A full multigrid method for nonlinear eigenvalue problems. Sci. China Math.. 2016, 59: 2037-2048

[21] Lieb EH. Thomas-Fermi and related theories of atoms and molecules. Rev. Mod. Phys.. 1981, 53: 603-641

[22] Lin Q, Xie H. A multi-level correction scheme for eigenvalue problems. Math. Comput.. 2015, 84(291): 71-88

[23] Lin Q, Yan N, Zhou A. A sparse finite element method with high accuracy. Part I. Numer. Math.. 2001, 88(4): 731-742

[24] Liu F, Stynes M, Zhou A. Postprocessed two-scale finite element discretizations, part I. SIAM J. Numer. Anal.. 2011, 49: 1947-1971

[25] Liu F, Zhou A. Two-scale finite element discretizations for partial differential equations. J. Comput. Math.. 2006, 24: 373-392

[26] Liu F, Zhou A. Localizations and parallelizations for two-scale finite element discretizations. Commun. Pure Appl. Anal.. 2007, 6(3): 757-773

[27] Liu F, Zhou A. Two-scale Boolean Galerkin discretizations for Fredholm integral equations of the second kind. SIAM J. Numer. Anal.. 2007, 45: 296-312

[28] Liu F, Zhu J. Two-scale sparse finite element approximations. Sci. China Math.. 2016, 59(4): 789-808

[29] Lyu T, Shih T, Liem C. Splitting Extrapolation and Combination Technique: a New Technology for Solving Multidimensional Problems in Parallel. 1998, Beijing, Science Press(in Chinese)

[30] Martin RM. Electronic Structure: Basic Theory and Practical Methods. 2020, Cambridge, Cambridge University Press

[31] Pflaum C, Zhou A. Error analysis of the combination technique. Numer. Math.. 1999, 84: 327-350

[32] Pousin J, Rappaz J. Consistency, stability, a priori and a posteriori errors for Petrov-Galerkin methods applied to nonlinear problems. Numer. Math.. 1994, 69: 213-231

[33] Wang YA, Carter EA. Schwartz SD. Orbital-free kinetic-energy density functional theory. Theoretical Methods in Condensed Phase Chemistry. 2000, Dordrecht, Kluwer: 117184

[34] Xie, H.: An augmented subspace method and its applications. J. Numer. Methods Comput. Appl. 41(3), 23 (2020) (in Chinese)

[35] Xie H, Xie M. A multigrid method for the ground state solution of Bose-Einstein condensates. J. Comput. Phys.. 2016, 19(3): 648-662

[36] Xu J, Zhou A. A two-grid discretization scheme for eigenvalue problems. Math. Comput.. 2001, 70: 17-25

[37] Xu Y, Zhou A. Fast Boolean approximation methods for solving integral equations in high dimensions. J. Integral Equ. Appl.. 2004, 16: 83-110

[38] Zhou A. An analysis of finite-dimensional approximations for the ground state solution of Bose-Einstein condensates. Nonlinearity. 2004, 17: 541-550

[39] Zhou A. Finite dimensional approximations for the electronic ground state solution of a molecular system. Math. Methods Appl. Sci.. 2007, 30: 429-447

Funding

Key Research and Development Program of Sichuan Province(2019YFA0709600 and 2019YFA0709601)

RIGHTS & PERMISSIONS

Shanghai University