scholarly journals A Fourier Continuation Method for the Solution of Elliptic Eigenvalue Problems in General Domains

2015 ◽  
Vol 2015 ◽  
pp. 1-15 ◽  
Author(s):  
Oscar P. Bruno ◽  
Timothy Elling ◽  
Ayon Sen

We present a new computational method for the solution of elliptic eigenvalue problems with variable coefficients in general two-dimensional domains. The proposed approach is based on use of the novel Fourier continuation method (which enables fast and highly accurate Fourier approximation of nonperiodic functions in equispaced grids without the limitations arising from the Gibbs phenomenon) in conjunction with an overlapping patch domain decomposition strategy and Arnoldi iteration. A variety of examples demonstrate the versatility, accuracy, and generality of the proposed methodology.

Author(s):  
Anthony M.J Davis ◽  
Stefan G Llewellyn Smith

Motivated by problems involving diffusion through small gaps, we revisit two-dimensional eigenvalue problems with localized perturbations to Neumann boundary conditions. We recover the known result that the gravest eigenvalue is O (|ln  ϵ | −1 ), where ϵ is the ratio of the size of the hole to the length-scale of the domain, and provide a simple and constructive approach for summing the inverse logarithm terms and obtaining further corrections. Comparisons with numerical solutions obtained for special geometries, both for the Dirichlet ‘patch problem’ where the perturbation to the boundary consists of a different boundary condition and for the gap problem, confirm that this approach is a simple way of obtaining an accurate value for the gravest eigenvalue and hence the long-term outcome of the underlying diffusion problem.


Sign in / Sign up

Export Citation Format

Share Document