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

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.

Keywords

Discontinuous Galerkin (DG) method / Multiscale method / Resonance errors / One-dimensional Schrödinger equation

References

[1.] Aarnes J, Heimsund B-O. . Multiscale discontinuous Galerkin methods for elliptic problems with multiple scales. Multiscale Methods in Science and Engineering. Lecture Notes in Computational Science and Engineering, 2005 Berlin Springer 1-20

[2.] Arnold A, Ben Abdallah N, Negulescu C. WKB-based schemes for the oscillatory 1D Schrödinger equation in the semiclassical limit. SIAM J. Numer. Anal., 2011, 49: 1436-1460,

[3.] Ben Abdallah N, Mouis M, Negulescu C. An accelerated algorithm for 2D simulations of the quantum ballistic transport in nanoscale MOSFETs. J. Comput. Phys., 2007, 225: 74-99,

[4.] Ben Abdallah N, Pinaud O. Multiscale simulation of transport in an open quantum system: resonances and WKB interpolation. J. Comput. Phys., 2006, 213: 288-310,

[5.] Cockburn B, Dong B. An analysis of the minimal dissipation local discontinuous Galerkin method for convection-diffusion problems. J. Sci. Comput., 2007, 32: 233-262,

[6.] Dong B, Shu C-W, Wang W. A new multiscale discontinuous Galerkin method for the one-dimensional stationary Schrödinger equation. J. Sci. Comput., 2016, 66: 321-345,

[7.] Dong B, Wang W. High-order multiscale discontinuous Galerkin methods for the one-dimensional stationary Schrödinger equation. J. Comput. Appl. Math., 2020, 380: 1-11,

[8.] Dong, B., Wang, W.: A high-order multiscale discontinuous Galerkin method for two-dimensional Schrödinger equation in quantum transport. J. Comput. Appl. Math. 418, 114701 (2023)

[9.] Dong, B., Wang, W.: A new multiscale discontinuous Galerkin method for a class of second-order equations with oscillatory solutions in two-dimensional space. In: Spectral and High Order Methods for Partial Differential Equations ICOSAHOM 2020+1. Lecture Notes in Computational Science and Engineering. Springer Nature Switzerland AG (accepted)

[10.] Efendiev YR, Hou TY. . Multiscale Finite Element Methods: Theory and Applications, 2009 New York Springer

[11.] Gabard G. Discontinuous Galerkin methods with plane waves for time-harmonic problems. J. Comput. Phys., 2007, 225: 1961-1984,

[12.] Gittelson C, Hiptmair R, Perugia I. Plane wave discontinuous Galerkin methods: analysis of the h-version. ESAIM: M2AN Math. Model. Numer. Anal., 2009, 43: 297-331,

[13.] Guo L, Xu Y. Local discontinuous Galerkin methods for the 2D simulation of quantum transport phenomena on quantum directional coupler.. Commun. Comput. Phys., 2014, 15: 1012-1028,

[14.] Negulescu C. Numerical analysis of a multiscale finite element scheme for the resolution of the stationary Schrödinger equation. Numer. Math., 2008, 108: 625-652,

[15.] Negulescu C, Ben Abdallah N, Polizzi E, Mouis M. Simulation schemes in 2D nanoscale MOSFETs: a WKB based method. J. Comput. Electron., 2004, 3: 397-400,

[16.] Polizzi E, Ben Abdallah N. Subband decomposition approach for the simulation of quantum electron transport in nanostructures. J. Comput. Phys., 2005, 202: 150-180,

[17.] Wang W, Guzmán J, Shu C-W. The multiscale discontinuous Galerkin method for solving a class of second order elliptic problems with rough coefficients. Int. J. Numer. Anal. Model., 2011, 8: 28-47

[18.] Wang W, Shu C-W. The WKB local discontinuous Galerkin method for the simulation of Schrödinger equation in a resonant tunneling diode. J. Sci. Comput., 2009, 40: 360-374,

[19.] Yuan L, Shu C-W. Discontinuous Galerkin method based on non-polynomial approximation spaces. J. Comput. Phys., 2006, 218: 295-323,

[20.] Yuan L, Shu C-W. Discontinuous Galerkin method for a class of elliptic multi-scale problems. Int. J. Numer. Methods Fluids, 2008, 56: 1017-1032,

[21.] Zhang Y, Wang W, Guzmán J, Shu C-W. Multi-scale discontinuous Galerkin method for solving elliptic problems with curvilinear unidirectional rough coefficients. J. Sci. Comput., 2014, 61: 42-60,

Funding

National Science Foundation(DMS-1818998)