A pFFT Accelerated Linear Strength BEM Potential Solver
D.J. Willis, J.K. White and J. Peraire
Massachusetts Institute of Technology, US
Keywords: pFFT, linear, potential, BEM, Galerkin
A linear strength, Galerkin Boundary Element Meth- od (BEM) for the solution of the three dimensional, direct potential boundary integral equation is presented. The method incorporates node based linear shape functions of the single and double layers on flat triangular elements. The BEM solution is accelerated using a precorrected Fast Fourier Transform algorithm (pFFT). Due to the extended compact support of the linear basis, there exist several approaches for implementing a linear strength pFFT. In this paper, two approaches are discussed and results are presented for the simpler of the two implementations. The work presented in this paper is applied to potential flow problems. Results are presented for flow solutions around spheres and aircraft wings. The results of the sphere simulations are compared with analytical solutions, while the solutions for the wings are compared with 2-Dimensional results. The results indicate accurate solutions of the potential flow around 3-Dimensional bodies. The linear basis shows improved accuracy when compared with the constant basis approach; however, the error of the linear BEM solution converges at a similar rate to the constant panels. This is due to the domination of the surface discretization error, which converges in the first order for planar element representations of curved surfaces.
Nanotech 2004 Conference Technical Program Abstract