On the spurious solutions in the high-order finite difference methods for eigenvalue problems

In this paper, the origins of spurious solutions occurring in the high-order finite difference methods are studied. Based on a uniform mesh, spurious modes are found in the high-order one-sided finite difference discretizations of many eigenvalue problems. Spurious modes are classified as spectral pollution and non-spectral pollution. The latter can be partially avoided by mesh refinement, while the former persists when the mesh is refined. Through numerical studies of some prototype eigenvalue problems, such as those of the Helmholtz and beam equations, we show that perfect central differentiation schemes do not produce any spurious modes. Nevertheless, high-order central difference schemes encounter difficulties in implementing complex boundary conditions. We further show that the central difference schemes will produce spurious modes as well, if asymmetric approximation is involved in boundary treatments. In general, central difference schemes are less likely to produce spurious modes than one-sided difference schemes.

[1]  Bo-nan Jiang,et al.  The Origin of Spurious Solutions in Computational Electromagnetics , 1996 .

[2]  D. Boffi,et al.  Computational Models of Electromagnetic Resonators: Analysis of Edge Element Approximation , 1999 .

[3]  C. Shu.  Differential Quadrature and Its Application in Engineering , 2000 .

[4]  B. Fornberg.  CALCULATION OF WEIGHTS IN FINITE DIFFERENCE FORMULAS∗ , 1998 .

[5]  B. A. Foreman,et al.  Elimination of spurious solutions from eight-band k.p theory , 1997 .

[6]  Shiang-Woei Chyuan,et al.  True and spurious eigensolutions for the problems with the mixed-type boundary conditions using BEMs , 2004 .

[7]  E B Davies,et al.  Spectral Pollution , 2002 .

[8]  Mathematical analysis and numerical study of true and spurious eigenequations for free vibration of plates using real-part BEM , 2004 .

[9]  Bengt Fornberg,et al.  Classroom Note: Calculation of Weights in Finite Difference Formulas , 1998, SIAM Rev..

[10]  Yang Xiang,et al.  Discrete singular convolution and its application to the analysis of plates with internal supports. Part 1: Theory and algorithm , 2002 .

[11]  Shan Zhao,et al.  High-order FDTD methods via derivative matching for Maxwell's equations with material interfaces , 2004 .

[12]  A. Carcaterra,et al.  ON THE SPURIOUS SOLUTIONS IN COMPLEX ENVELOPE DISPLACEMENT ANALYSIS , 2001 .

[13]  Shan Zhao,et al.  Tensor product derivative matching for wave propagation in inhomogeneous media , 2004 .

[14]  Shan Zhao,et al.  Local spectral time-domain method for electromagnetic wave propagation. , 2003, Optics letters.

[15]  Shan Zhao,et al.  High order matched interface and boundary method for elliptic equations with discontinuous coefficients and singular sources , 2006, J. Comput. Phys..

[16]  A. Bermúdez,et al.  Finite element vibration analysis of fluid-solid systems without spurious modes , 1995 .

[17]  J. Rappaz,et al.  On spectral pollution in the finite element approximation of thin elastic “membrane” shells , 1997 .

[18]  P.A.A. Laura,et al.  Vibrations of non-uniform rings studied by means of the differential quadrature method , 1995 .

[19]  Daniele Boffi,et al.  On the problem of spurious eigenvalues in the approximation of linear elliptic problems in mixed form , 2000, Math. Comput..

[20]  Shan Zhao,et al.  Comparison of the Discrete Singular Convolution and Three Other Numerical Schemes for Solving Fisher's Equation , 2003, SIAM J. Sci. Comput..

[21]  Chang Shu,et al.  Treatment of mixed and nonuniform boundary conditions in GDQ vibration analysis of rectangular plates , 1999 .

[22]  Guo-Wei Wei,et al.  Discrete singular convolution for the solution of the Fokker–Planck equation , 1999 .

[23]  T. C. Mcgill,et al.  Numerical spurious solutions in the effective mass approximation , 2003 .

[24]  Richard B. Pember,et al.  Numerical Methods for Hyperbolic Conservation Laws With Stiff Relaxation I. Spurious Solutions , 1993, SIAM J. Appl. Math..

[25]  Guo-Wei Wei,et al.  A note on the numerical solution of high-order differential equations , 2003 .

[26]  Y. Xiang,et al.  DSC analysis of free-edged beams by an iteratively matched boundary method , 2005 .