Open Access

ARTICLE

Fast BEM Solvers for 3D Poisson-Type Equations

by Xuefei He¹, Kian-Meng Lim¹, Siak-Piang Lim¹

1 Singapore-MIT Alliance, E4-04-10, 4 Engineering Drive 3, Singapore 117576.
2 Department of Mechanical Engineering, National University of Singapore, Singapore 119260.
3 Corresponding author. Email: limkm@nus.edu.sg.

Computer Modeling in Engineering & Sciences 2008, 35(1), 21-48. https://doi.org/10.3970/cmes.2008.035.021

Abstract

The boundary element method (BEM) is known to have the advantage of reducing the dimension of problem by discretizing only the boundary of the domain. But it becomes less attractive for solving Poisson-type equations, due to the need to evaluate the domain integral which is computationally expensive. In this paper, we present the extension of a recently developed fast algorithm for Laplace equation, based on fast Fourier transform on multipoles (FFTM), to solve large scale 3D Poisson-type equations. We combined the Laplace solver with two fast methods for handling the domain integral based on fast Fourier transform (FFT). The first method uses the FFT on multipoles to accelerate the domain integral, while the second method solves the domain integral as a particular solution using FFT. The particular solution method is found to be faster and more accurate, and it is extended to solve non-linear Poisson-type equations. The algorithm is shown to be efficient when it is used in the inner loop of the iterative solver for the non-linear equations.

---

Keywords

---