A SPATIALLY ADAPTIVE ITERATIVE METHOD FOR A CLASS OF NONLINEAR OPERATOR EIGENPROBLEMS
2014 (English)In: Electronic Transactions on Numerical Analysis, ISSN 1068-9613, Vol. 41, 21-41 p.Article in journal (Refereed) Published
We present a new algorithm for the iterative solution of nonlinear operator eigenvalue problems arising from partial differential equations (PDEs). This algorithm combines automatic spatial resolution of linear operators with the infinite Arnoldi method for nonlinear matrix eigenproblems proposed by Jarlebring et al. [Numer. Math., 122 (2012), pp. 169-195]. The iterates in this infinite Arnoldi method are functions, and each iteration requires the solution of an inhomogeneous differential equation. This formulation is independent of the spatial representation of the functions, which allows us to employ a dynamic representation with an accuracy of about the level of machine precision at each iteration similar to what is done in the Chebfun system with its chebop functionality although our function representation is entirely based on coefficients instead of function values. Our approach also allows nonlinearity in the boundary conditions of the PDE. The algorithm is illustrated with several examples, e.g., the study of eigenvalues of a vibrating string with delayed boundary feedback control.
Place, publisher, year, edition, pages
2014. Vol. 41, 21-41 p.
Arnoldi's method, nonlinear eigenvalue problems, partial differential equations, Krylov subspaces, delay-differential equations, Chebyshev polynomials
IdentifiersURN: urn:nbn:se:kth:diva-161594ISI: 000348498600003ScopusID: 2-s2.0-84909994515OAI: oai:DiVA.org:kth-161594DiVA: diva2:798019
QC 201503252015-03-252015-03-132015-03-25Bibliographically approved