Communications on Applied Mathematics and Computation ›› 2024, Vol. 6 ›› Issue (1): 311-324.doi: 10.1007/s42967-022-00248-4

• ORIGINAL PAPERS • Previous Articles     Next Articles

Numerical Investigations on the Resonance Errors of Multiscale Discontinuous Galerkin Methods for One-Dimensional Stationary Schrödinger Equation

Bo Dong1, Wei Wang2   

  1. 1. Department of Mathematics, University of Massachusetts Dartmouth, North Dartmouth, MA, 02747, USA;
    2. Department of Mathematics and Statistics, Florida International University, Miami, FL, 33199, USA
  • Received:2022-08-30 Revised:2022-12-19 Published:2024-04-16
  • Contact: Wei Wang,E-mail:weiwang1@fiu.edu;Bo Dong,E-mail:bdong@umassd.edu E-mail:weiwang1@fiu.edu;bdong@umassd.edu
  • Supported by:
    The authors consent to comply with all the Publication Ethical Standards. The research of the first author is supported by the National Science Foundation grant DMS-1818998.

Abstract: In this paper, numerical experiments are carried out to investigate the impact of penalty parameters in the numerical traces on the resonance errors of high-order multiscale discontinuous Galerkin (DG) methods (Dong et al. in J Sci Comput 66: 321-345, 2016; Dong and Wang in J Comput Appl Math 380: 1-11, 2020) for a one-dimensional stationary Schrödinger equation. Previous work showed that penalty parameters were required to be positive in error analysis, but the methods with zero penalty parameters worked fine in numerical simulations on coarse meshes. In this work, by performing extensive numerical experiments, we discover that zero penalty parameters lead to resonance errors in the multiscale DG methods, and taking positive penalty parameters can effectively reduce resonance errors and make the matrix in the global linear system have better condition numbers.

Key words: Discontinuous Galerkin (DG) method, Multiscale method, Resonance errors, One-dimensional Schrödinger equation