×

A spatially adaptive iterative method for a class of nonlinear operator eigenproblems. (English) Zbl 1295.65109

Summary: 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 E. Jarlebring et al. [Numer. Math. 122, No. 1, 169–195 (2012; Zbl 1256.65043)]. 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.

MSC:

65N25 Numerical methods for eigenvalue problems for boundary value problems involving PDEs
35P15 Estimates of eigenvalues in context of PDEs
35P30 Nonlinear eigenvalue problems and nonlinear spectral theory for PDEs

Citations:

Zbl 1256.65043