Fast multipole cell-based domain integration method for treatment of volume potentials in 3D elasticity problems

Abstract

Purpose Domain integrals, known as volume potentials in 3D elasticity problems, exist in many boundary-type methods, such as the boundary element method (BEM) for inhomogeneous partial differential equations. The purpose of this paper is to develop an accurate and reliable technique to effectively evaluate the volume potentials in 3D elasticity problems. Design/methodology/approach An adaptive background cell-based domain integration method is proposed for treatment of volume potentials in 3D elasticity problems. The background cells are constructed from the information of the boundary elements based on an oct-tree structure, and the domain integrals are evaluated over the cells rather than volume elements. The cells that contain the boundary elements can be subdivided into smaller sub-cells adaptively according to the sizes and levels of the boundary elements. The fast multipole method (FMM) is further applied in the proposed method to reduce the time complexity of large-scale computation. Findings The method is a boundary-only discretization method, and it can be applied in the BEM easily. Much computational time is saved by coupling with the FMM. Numerical examples demonstrate the accuracy and efficiency of the proposed method.. Originality/value Boundary elements are used to create adaptive background cells, and domain integrals are evaluated over the cells rather than volume elements. Large-scale computation is made possible by coupling with the FMM.