×

On the spurious solutions in the high-order finite difference methods for eigenvalue problems. (English) Zbl 1173.74462

Summary: 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.

MSC:

74S20 Finite difference methods applied to problems in solid mechanics
76M20 Finite difference methods applied to problems in fluid mechanics
74K10 Rods (beams, columns, shafts, arches, rings, etc.)
76Q05 Hydro- and aero-acoustics
Full Text: DOI

References:

[1] Bao, G.; Wei, G. W.; Zhao, S., Local spectral time-domain method for electromagnetic wave propagation, Opt. Lett., 28, 513-515 (2003)
[2] Bermudez, A.; Duran, R.; Muschietti, M. A.; Rodriguez, R.; Solomin, J., Finite element vibration analysis of fluid-solid systems without spurious modes, SIAM J. Numer. Anal., 32, 1280-1295 (1995) · Zbl 0833.73050
[3] Boffi, D.; Brezzi, F.; Gastaldi, L., On the problem of spurious eigenvalues in the approximation of linear elliptic problems in mixed form, Math. Comput., 229, 121-140 (1999) · Zbl 0938.65126
[4] Boffi, D.; Fernandes, P.; Gastaldi, L.; Perugia, I., Computational models of electromagnetic resonators: analysis of edge element approximation, SIAM J. Numer. Anal., 36, 1264-1290 (1999) · Zbl 1025.78014
[5] Cartoixa, X.; Ting, D. Z.Y.; McGill, T. C., Numerical spurious solutions in the effective mass approximation, J. Appl. Phys., 93, 3974-3981 (2003)
[6] Chen, J. T.; Lin, S. Y.; Chen, K. H.; Chen, I. L., Mathematical analysis and numerical study of true and spurious eigenequations for free vibration of plates using real-part BEM, Comput. Mech., 34, 165-180 (2004) · Zbl 1095.74015
[7] Chen, J. T.; Lin, T. W.; Chen, K. H.; Chyuan, S. W., True and spurious eigensolutions for the problems with the mixed-type boundary conditions using BEMs, Finite Elements Anal. Des., 40, 1521-1549 (2004)
[8] Davies, E. B.; Plum, M., Spectral pollution, IMA J. Numer. Anal., 24, 417-438 (2004) · Zbl 1062.65056
[9] Foreman, B. A., Elimination of spurious solutions from eight-band \(k\)·\(p\) theory, Phys. Rev. B, 56, R12748-R12751 (1997)
[10] Fornberg, B., Calculation of weights in finite difference formulas, SIAM Rev., 40, 685-691 (1998) · Zbl 0914.65010
[11] Gutierrez, R. H.; Laura, P. A.A., Vibrations of non-uniform rings studied by means of the differential quadrature method, J. Sound Vib., 185, 507-513 (1995) · Zbl 1048.74525
[12] Jiang, B. N.; Wu, J.; Povinelli, L. A., The origin of spurious solutions in computational electromagnetics, J. Comput. Phys., 125, 104-123 (1996) · Zbl 0848.65086
[13] Pember, R. B., Numerical methods for hyperbolic conservation laws with stiff relaxation I. Spurious solutions, SIAM J. Appl. Math., 53, 1293-1330 (1993) · Zbl 0787.65062
[14] Rappaz, J.; Hubert, J. S.; Palencia, E. S.; Vassiliev, D., On spectral pollution in the finite element approximation of the elastic “membrane” shells, Numer. Math., 75, 473-500 (1997) · Zbl 0874.73066
[15] Sestieri, A.; Carcaterra, A., On the spurious solutions in complex envelope displacement analysis, J. Sound Vib., 240, 293-302 (2001)
[16] Shu, C., Differential Quadrature and its Application in Engineering (2000), Springer-Verlag: Springer-Verlag London · Zbl 0944.65107
[17] Shu, C.; Wang, C. M., Treatment of mixed and nonuniform boundary conditions in GDQ vibration analysis of rectangular plates, Engrg. Struct., 21, 125-134 (1999)
[18] Wang, Y.; Zhao, Y. B.; Wei, G. W., A note on the numerical solution of high-order differential equations, J. Comput. Appl. Math., 159, 387-398 (2003) · Zbl 1031.65087
[19] Wei, G. W., Discrete singular convolution for solution of the Fokker-Planck equations, J. Chem. Phys., 110, 8930-8942 (1999)
[20] Wei, G. W.; Zhao, Y. B.; Xiang, Y., Discrete singular convolution and its application to the analysis of plates with internal supports. I Theory and algorithm, Int. J. Numer. Methods Engrg., 55, 913-946 (2002) · Zbl 1058.74643
[21] Zhao, S.; Wei, G. W., Comparison of the discrete singular convolution and three other numerical schemes for solving Fisher’s equation, SIAM J. Sci. Comput., 25, 127-147 (2003) · Zbl 1043.35015
[22] Zhao, S.; Wei, G. W., high order FDTD methods via derivative matching for Maxwell’s equations with material interfaces, J. Comput. Phys., 200, 60-103 (2004) · Zbl 1050.78018
[23] Zhao, S.; Wei, G. W., Tensor product derivative matching for wave propagation in inhomogeneous media, Microwave Opt. Technol. Lett., 43, 69-77 (2004)
[24] S. Zhao, G.W. Wei, A new boundary implementation method for high-order central finite difference, preprint 2007.; S. Zhao, G.W. Wei, A new boundary implementation method for high-order central finite difference, preprint 2007.
[25] Zhao, S.; Wei, G. W.; Xiang, Y., DSC analysis of free-edged beams by an iteratively matched boundary method, J. Sound Vib., 284, 487-493 (2005) · Zbl 1237.74198
[26] Zhou, Y. C.; Zhao, S.; Feig, M.; Wei, G. W., High order matched interface and boundary method for elliptic equations with discontinuous coefficients and singular sources, J. Comput. Phys., 213, 1-30 (2006) · Zbl 1089.65117
This reference list is based on information provided by the publisher or from digital mathematics libraries. Its items are heuristically matched to zbMATH identifiers and may contain data conversion errors. In some cases that data have been complemented/enhanced by data from zbMATH Open. This attempts to reflect the references listed in the original paper as accurately as possible without claiming completeness or a perfect matching.