摘要
The numerical quadrature methods for dealing with the problems of singular and near-singular integrals caused by Burton-Miller method are proposed, by which the conventional and fast multipole BEMs （boundary element methods） for 3D acoustic problems based on constant elements are improved. To solve the problem of singular integrals, a Hadamard finite-part integral method is presented, which is a simplified combination of the methods proposed by Kirkup and Wolf. The problem of near-singular integrals is overcome by the simple method of polar transformation and the more complex method of PART （Projection and Angular ＆ Radial Transformation）. The effectiveness of these methods for solving the singular and near-singular problems is validated through comparing with the results computed by the analytical method and/or the commercial software LMS Virtual.Lab. In addition, the influence of the near-singular integral problem on the computational precisions is analyzed by computing the errors relative to the exact solution. The computational complexities of the conventional and fast multipole BEM are analyzed and compared through numerical computations. A large-scale acoustic scattering problem, whose degree of freedoms is about 340,000, is implemented successfully. The results show that, the near singularity is primarily introduced by the hyper-singular kernel, and has great influences on the precision of the solution. The precision of fast multipole BEM is the same as conventional BEM, but the computational complexities are much lower.
The numerical quadrature methods for dealing with the problems of singular and near-singular integrals caused by Burton-Miller method are proposed, by which the conventional and fast multipole BEMs （boundary element methods） for 3D acoustic problems based on constant elements are improved. To solve the problem of singular integrals, a Hadamard finite-part integral method is presented, which is a simplified combination of the methods proposed by Kirkup and Wolf. The problem of near-singular integrals is overcome by the simple method of polar transformation and the more complex method of PART （Projection and Angular ＆ Radial Transformation）. The effectiveness of these methods for solving the singular and near-singular problems is validated through comparing with the results computed by the analytical method and/or the commercial software LMS Virtual.Lab. In addition, the influence of the near-singular integral problem on the computational precisions is analyzed by computing the errors relative to the exact solution. The computational complexities of the conventional and fast multipole BEM are analyzed and compared through numerical computations. A large-scale acoustic scattering problem, whose degree of freedoms is about 340,000, is implemented successfully. The results show that, the near singularity is primarily introduced by the hyper-singular kernel, and has great influences on the precision of the solution. The precision of fast multipole BEM is the same as conventional BEM, but the computational complexities are much lower.
基金
supported by the National Natural Science Foundation of China(11304344,11404364)
the Project of Hubei Provincial Department of Education(D20141803)
the Natural Science Foundation of Hubei Province(2014CFB378)
the Doctoral Scientific Research Foundation of Hubei University of Automotive Technology(BK201604)
