Skip to main content
SpringerLink
Log in

On the numerical treatment of the eigenparameter dependent boundary conditions

  • Original Paper
  • Published:
Numerical Algorithms Aims and scope Submit manuscript

Abstract

In this paper, we consider the numerical treatment of singular eigenvalue problems supplied with eigenparameter dependent boundary conditions using spectral methods. On the one hand, such boundary conditions hinder the construction of test and trial space functions which could incorporate them and thus providing well-conditioned Galerkin discretization matrices. On the other hand, they can generate surprising behavior of the eigenvectors hardly detected by analytic methods. These singular problems are often indirectly approximated by regular ones. We argue that spectral collocation as well as tau method offer remedies for the first two issues and provide direct and efficient treatment to such problems. On a finite domain, we consider the so-called Petterson-König’s rod eigenvalue problem and on the half line, we take into account the Charney’s baroclinic stability problem and the Fourier eigenvalue problem. One boundary condition in these problems depends on the eigenparameter and additionally, this also could depend on some physical parameters. The Chebyshev collocation based on both, square and rectangular differentiation and a Chebyshev tau method are used to discretize the first problem. All these schemes cast the problems into singular algebraic generalized eigenvalue ones which are solved by the QZ and/or Arnoldi algorithms as well as by some target oriented Jacobi-Davidson methods. Thus, the spurious eigenvalues are completely eliminated. The accuracy of square Chebyshev collocation is roughly estimated and its order of approximation with respect to the eigenvalue of interest is determined. For the problems defined on the half line, we make use of the Laguerre-Gauss-Radau collocation. The method proved to be reliable, accurate, and stable with respect to the order of approximation and the scaling parameter.

This is a preview of subscription content, log in via an institution to check access.

Access this article

Log in via an institution

Price excludes VAT (USA)
Tax calculation will be finalised during checkout.

Instant access to the full article PDF.

Institutional subscriptions

Similar content being viewed by others

References

  1. Boyd, J.P.: Orthogonal rational functions on a semi-infinite interval. J. Comp. Phys. 70, 63–88 (1987)

    Article  MathSciNet  MATH  Google Scholar 

  2. Boyd, J.P.: Chebyshev domain truncation is inferior to fourier domain truncation for solving problems on an infinite interval. J. Sci. Comput. 3, 109–120 (1988)

    Article  MathSciNet  MATH  Google Scholar 

  3. Boyd, J.P., Rangan, C., Bucksbaum, P.H.: Pseudospectral methods on a semi-infinite interval with application to the hydrogen atom: a comparison of the mapped Fourier sine method with Laguerre series and rational Chebyshev expansions. J. Comp. Phys. 188, 56–74 (2003)

    Article  MATH  Google Scholar 

  4. Boyd, J.P.: Five themes in Chebyshev spectral methods applied to the regularized Charney eigenproblem: extra numerical boundary conditions, a boundary-layer-resolving change of coordinate, parametrzing a curve which is singularat at an endpoint, extending the tau method to log-and-polynomials and finding the roots of a polynomial-and-log approximation. Comput. Math. Appl. 71, 1277–1241 (2016)

    Article  Google Scholar 

  5. Branscome, L.E.: The charney baroclinic stability problem: approximate solutions and modal structures. J. Atmos. Sci., 1393–1409 (1983)

  6. Driscoll, T.A., Hale, N., Trefethen, L.N.: Chebfun Guide. Pafnuty Publications, Oxford (2014)

    Google Scholar 

  7. Driscoll, T.A., Hale, N.: Rectangular spectral collocation. IMA J. Numer. Anal. doi:10.1093/imanum/dru062

  8. Gheorghiu, C.I., Pop, I.S.: A modified Chevyshev-Tau method for a hydrodynamic stability problem. Approximation and Optimization Proceedings of the International Conference on Approximation and Optimization (Romania)-ICAOR Cluj-Napoca, vol. II, pp 119–126 (1996)

  9. Gheorghiu, C.I., Rommes, J.: Application of the Jacobi-Davidson method to accurate analysis of singular linear hydrodynamic stability problems. Int. J. Numer. Method Fl. (2012)

  10. Gheorghiu, C.I.: Spectral methods for non-standard eigenvalue problems. Fluid and Structural Mechanics and Beyond. Springer Cham Heidelberg New York Dondrecht London (2014)

  11. Gheorghiu, C.I.: Stable spectral collocation solutions to a class of Benjamin Bona Mahony initial value problems. Appl. Math. Comput. 273, 1090–1099 (2016)

    MathSciNet  Google Scholar 

  12. Gheorghiu, C.I.: Spectral collocation solutions to systems of boundary layer type. Numer. Algor. doi:10.1007/s11075-015-0083-6

  13. Gottlieb, D., Orszag, S.A.: Numerical analysis of spectral methods: theory and applications. SIAM, Philadelphia (1977)

    Book  MATH  Google Scholar 

  14. Goussis, D.A., Pearlstein, A.J.: Removal of infinite eigenvalues in the generalized matrix eigenvalue problems. J. Comput. Phys. 84, 242–246 (1998)

    Article  MathSciNet  MATH  Google Scholar 

  15. Maohzu, Z., Sun, J., Zettl, A.: The spectrum of singular Sturm-Liouville problems with eigenparameter dependent boundary conditions and its approximation. Results. Math. 63, 1311—1330 (2013)

    MathSciNet  MATH  Google Scholar 

  16. van Noorden, T., Rommes, J.: Computing a partial generalized real Schur form using the Jacobi-Davidson method. Numer. Linear Algebra Appl. 14, 197–215 (2007)

    Article  MathSciNet  MATH  Google Scholar 

  17. Rommes, J.: Arnoldi and Jacobi-Davidson methods for generalized eigenvalue problems A x = λ B x with B singular. Math. of Comput. 77, 995–1015 (2008)

    Article  MathSciNet  MATH  Google Scholar 

  18. Shen, J., Tang, T., Wang, L.-L.: Spectral methods. Algorithms, Analysis and Applications Springer Heidelberg Dordrecht London New York (2011)

  19. Trefethen, L.N.: Computation of pseudospectra. Acta Numer. 8, 247–295 (1999)

    Article  MathSciNet  MATH  Google Scholar 

  20. Tretter, Ch.: Nonselfadjoint spectral problems for linear pencils N−λ P of ordinary differential operators with λ-linear boundary conditions: completeness results. Integr. Equat. Oper. Th. 26, 222–248 (1996)

    Article  MathSciNet  MATH  Google Scholar 

  21. Tretter, C.h.: A linearization for a class of λ− nonlinear boundary eigenvalue problems. J. Math. Analysis Appl. 247, 331–355 (2000)

    Article  MathSciNet  MATH  Google Scholar 

  22. Tretter, Ch.: Boundary eigenvalue problems for differential equations N η = λ P η with λ-polynomial boundary conditions. Integr. J. Diff. Eqns. 170, 408–471 (2001)

    Article  MathSciNet  Google Scholar 

  23. Xu, K., Hale, N.: Explicit construction of rectangular differential matrices. IMA J. Numer. Anal. doi:10.1093/imanum/drv013

  24. Xu, K., Hale, N.: The Chebyshev points of the first kind. Appl. Numer. Math. 102, 17–30 (2016)

    Article  MathSciNet  MATH  Google Scholar 

  25. Wang, C.M., Wang, C.Y., Reddy, J.N.: Exact solutions for buckling of structural members. CRC Press (2005)

  26. Weideman, J.A.C., Reddy, S.C.: A MATLAB differentiation matrix suite. ACM T. Math. Software 26, 465–519 (2000)

    Article  MathSciNet  Google Scholar 

  27. von Winckel, G.: Fast Chebyshev Transform, MathWorksⓇ. File Exchanges (2005)

Download references

Acknowledgements

The author is grateful to both anonymous referees for their careful reading of the paper and suggestions which resulted in considerable improvement. He is also fairly indebted to Prof. John Boyd for illuminating discussions on the Charney’s problem during our visit to NAG in Oxford with the occasion of Chebyshev Day, 2016. With the same opportunity, Prof. Nick Trefethen provided us hints on rectangular differentiation.

Author information

Authors and Affiliations

  1. Romanian Academy, “T. Popoviciu” Institute of Numerical Analysis, Cluj-Napoca, Romania

    Călin-Ioan Gheorghiu

Authors
  1. Călin-Ioan Gheorghiu
    View author publications

    You can also search for this author in PubMed Google Scholar

Corresponding author

Correspondence to Călin-Ioan Gheorghiu.

Rights and permissions

Reprints and permissions

About this article

Check for updates. Verify currency and authenticity via CrossMark

Cite this article

Gheorghiu, CI. On the numerical treatment of the eigenparameter dependent boundary conditions. Numer Algor 77, 77–93 (2018). https://doi.org/10.1007/s11075-017-0305-1

Download citation

  • Received:

  • Accepted:

  • Published:

  • Issue Date:

  • DOI: https://doi.org/10.1007/s11075-017-0305-1

Keywords

Search

Navigation