×

Implicit boundary conditions in partial differential equations discretizations: identifying spurious modes and model reduction. (English) Zbl 07811330

Summary: We revisit the problem of spurious modes that are sometimes encountered in partial differential equations discretizations. It is generally suspected that one of the causes for spurious modes is due to how boundary conditions are treated, and we use this as the starting point of our investigations. By regarding boundary conditions as algebraic constraints on a differential equation, we point out that any differential equation with homogeneous boundary conditions also admits a typically infinite number of hidden or implicit boundary conditions. In most discretization schemes, these additional implicit boundary conditions are violated, and we argue that this is what leads to the emergence of spurious modes. These observations motivate two definitions of the quality of computed eigenvalues based on violations of derivatives of boundary conditions on the one hand, and on the Grassmann distance between subspaces associated with computed eigenspaces on the other. Both of these tests are based on a standardized treatment of boundary conditions and do not require a priori knowledge of eigenvalue locations. The effectiveness of these tests is demonstrated on several examples known to have spurious modes. In addition, these quality tests show that in most problems, about half of the computed spectrum of a differential operator is of low quality. The tests also specifically identify the low accuracy modes, which can then be projected out as a type of model reduction scheme.

MSC:

76Exx Hydrodynamic stability
65Lxx Numerical methods for ordinary differential equations
76Mxx Basic methods in fluid mechanics

References:

[1] Bourne, D., Hydrodynamic stability, the Chebyshev tau method and spurious eigenvalues. Contin. Mech. Thermodyn., 6, 571-579 (2003) · Zbl 1063.76074
[2] Gardner, D. R.; Trogdon, S. A.; Douglass, R. W., A modified tau spectral method that eliminates spurious eigenvalues (1989), Tech. Rep. 1 · Zbl 0661.65084
[3] Huang, W.; Sloan, D. M., The pseudospectral method for solving differential eigenvalue problems. J. Comput. Phys., 2, 399-409 (1994) · Zbl 0799.65091
[4] Lindsay, K. A.; Ogden, R. R., A practical implementation of spectral methods resistant to the generation of spurious eigenvalues. Int. J. Numer. Methods Fluids, 1277-1294 (1992) · Zbl 0762.76082
[5] McFadden, G. B.; Murray, B. T.; Boisvert, R. F., Elimination of spurious eigenvalues in the Chebyshev tau spectral method. J. Comput. Phys., 1, 228-239 (1990) · Zbl 0717.65063
[6] Zebib, A., Removal of spurious modes encountered in solving stability problems by spectral methods. J. Comput. Phys., 2, 521-525 (1987) · Zbl 0612.65049
[7] Xu, Y. S., Origin and elimination of spurious modes in the solution of field eigenvalue problems by the method of moments, 465-468
[8] Straughan, B.; Walker, D. W., Two very accurate and efficient methods for computing eigenvalues and eigenfunctions in porous convection problems. J. Comput. Phys., 128-141 (1996) · Zbl 0858.76064
[9] Fornberg, B., An improved pseudospectral method for initial-boundary value problems. J. Comput. Phys., 381-397 (1990) · Zbl 0717.65073
[10] Danilov, S.; Kutsenko, A., On the geometric origin of spurious waves in finite-volume discretizations of shallow water equations on triangular meshes. J. Comput. Phys. (2019) · Zbl 1453.65245
[11] Schroeder, W.; Wolff, I., The origin of spurious modes in numerical solutions of electromagnetic field eigenvalue problems. IEEE Trans. Microw. Theory Tech., 4, 644-653 (1994)
[12] Manning, M. L.; Bamieh, B.; Carlson, J. M., Descriptor approach for eliminating spurious eigenvalues in hydrodynamic equations, 1-13 (2007)
[13] Phillips, T. N.; Roberts, G. W., The treatment of spurious pressure modes in spectral incompressible flow calculations. J. Comput. Phys., 1, 150-164 (1993) · Zbl 0766.76059
[14] Bewley, T. R.; Liu, S., Optimal and robust control and estimation of linear paths to transition. J. Fluid Mech., 305-349 (1998) · Zbl 0924.76028
[15] Orszag, S. A., Accurate solution of the Orr-Sommerfeld stability equation. J. Fluid Mech., 4, 689-703 (1971) · Zbl 0237.76027
[16] Canuto, C.; Quarteroni, A., On the boundary treatment in spectral methods for hyperbolic systems. J. Comput. Phys., 1, 100-110 (1987) · Zbl 0626.65123
[17] Hespanha, J. P., Linear Systems Theory (2018), Princeton University Press · Zbl 1381.93001
[18] Canuto, C.; Hussaini, M. Y.; Quarteroni, A.; Zang, T. A., Spectral Methods: Evolution to Complex Geometries and Applications to Fluid Dynamics (2007), Springer Science & Business Media · Zbl 1121.76001
[19] Bayliss, A.; Class, A.; Matkowsky, B. J., Roundoff error in computing derivatives using the Chebyshev differentiation matrix. J. Comput. Phys., 2, 380-383 (1995) · Zbl 0826.65014
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.