全部 标题 作者
关键词 摘要

OALib Journal期刊
ISSN: 2333-9721
费用:99美元

查看量下载量

相关文章

更多...
Physics  2012 

An efficient solver for volumetric scattering based on fast spherical harmonics transforms

Full-Text   Cite this paper   Add to My Lib

Abstract:

The Helmholtz equation arises in the study of electromagnetic radiation, optics, acoustics, etc. In spherical coordinates, its general solution can be written as a spherical harmonic series which satisfies the radiation condition at infinity, ensuring that the wave is outgoing. The boundary condition at infinity is hard to enforce with a finite element method since a suitable approximation needs to be made within reasonable distance from scatterers. Luckily, the Helmholtz equation can be represented as a Lippmann-Schwinger integral equation which removes the necessity of the boundary approximations and its Green's function can be expanded as a spherical harmonic series which leads to our numerical scheme based on spherical harmonic polynomial transform. In this paper, we present an efficient solver for the Helmholtz equation which costs $O(N\log N)$ operations, where $N$ is the number of the discretization points. We use the fast spherical harmonic transforms which are originally developed in \cite{suda}. The convergence order of the method is tied to the global regularity of the solution. At the lower end, it is second order accurate for discontinuous material properties. The order increases with increasing regularity leading to spectral convergence for globally smooth solutions.

Full-Text

Contact Us

service@oalib.com

QQ:3279437679

WhatsApp +8615387084133