TY - JOUR
T1 - Solving eigenvalue problems on curved surfaces using the Closest Point Method
AU - Macdonald, Colin B.
AU - Brandman, Jeremy
AU - Ruuth, Steven J.
N1 - KAUST Repository Item: Exported on 2020-10-01
Acknowledged KAUST grant number(s): KUK-C1-013-04
Acknowledgements: The work of this author was supported by an NSERC postdoctoral fellowship, NSF grant No. CCF-0321917, and by Award No. KUK-C1-013-04 made by King Abdullah University of Science and Technology (KAUST).The work of this author was supported by an NSF Mathematical Sciences Postdoctoral Research Fellowship.The work of this author was partially supported by a Grant from NSERC Canada.
This publication acknowledges KAUST support, but has no KAUST affiliated authors.
PY - 2011/6
Y1 - 2011/6
N2 - Eigenvalue problems are fundamental to mathematics and science. We present a simple algorithm for determining eigenvalues and eigenfunctions of the Laplace-Beltrami operator on rather general curved surfaces. Our algorithm, which is based on the Closest Point Method, relies on an embedding of the surface in a higher-dimensional space, where standard Cartesian finite difference and interpolation schemes can be easily applied. We show that there is a one-to-one correspondence between a problem defined in the embedding space and the original surface problem. For open surfaces, we present a simple way to impose Dirichlet and Neumann boundary conditions while maintaining second-order accuracy. Convergence studies and a series of examples demonstrate the effectiveness and generality of our approach. © 2011 Elsevier Inc.
AB - Eigenvalue problems are fundamental to mathematics and science. We present a simple algorithm for determining eigenvalues and eigenfunctions of the Laplace-Beltrami operator on rather general curved surfaces. Our algorithm, which is based on the Closest Point Method, relies on an embedding of the surface in a higher-dimensional space, where standard Cartesian finite difference and interpolation schemes can be easily applied. We show that there is a one-to-one correspondence between a problem defined in the embedding space and the original surface problem. For open surfaces, we present a simple way to impose Dirichlet and Neumann boundary conditions while maintaining second-order accuracy. Convergence studies and a series of examples demonstrate the effectiveness and generality of our approach. © 2011 Elsevier Inc.
UR - http://hdl.handle.net/10754/599671
UR - https://linkinghub.elsevier.com/retrieve/pii/S0021999111003858
UR - http://www.scopus.com/inward/record.url?scp=81155160907&partnerID=8YFLogxK
U2 - 10.1016/j.jcp.2011.06.021
DO - 10.1016/j.jcp.2011.06.021
M3 - Article
SN - 0021-9991
VL - 230
SP - 7944
EP - 7956
JO - Journal of Computational Physics
JF - Journal of Computational Physics
IS - 22
ER -