An adaptive fast multipole boundary element method for three-dimensional acoustic wave problems based on the Burton–Miller formulation

The high solution costs and non-uniqueness difficulties in the boundary element method (BEM) based on the conventional boundary integral equation (CBIE) formulation are two main weaknesses in the BEM for solving exterior acoustic wave problems. To tackle these two weaknesses, an adaptive fast multipole boundary element method (FMBEM) based on the Burton–Miller formulation for 3-D acoustics is presented in this paper. In this adaptive FMBEM, the Burton–Miller formulation using a linear combination of the CBIE and hypersingular BIE (HBIE) is applied to overcome the non-uniqueness difficulties. The iterative solver generalized minimal residual (GMRES) and fast multipole method (FMM) are adopted to improve the overall computational efficiency. This adaptive FMBEM for acoustics is an extension of the adaptive FMBEM for 3-D potential problems developed by the authors recently. Several examples on large-scale acoustic radiation and scattering problems are presented in this paper which show that the developed adaptive FMBEM can be several times faster than the non-adaptive FMBEM while maintaining the accuracies of the BEM.     

A fast multipole boundary element method for solving the thin plate bending problem

A fast multipole boundary element method (BEM) for solving large-scale thin plate bending problems is presented in this paper. The method is based on the Kirchhoff thin plate bending theory and the biharmonic equation governing the deflection of the plate. First, the direct boundary integral equations and the conventional BEM for thin plate bending problems are reviewed. Second, the complex notation of the kernel functions, expansions and translations in the fast multipole BEM are presented. Finally, a few numerical examples are presented to show the accuracy and efficiency of the fast multipole BEM in solving thin plate bending problems. The bending rigidity of a perforated plate is evaluated using the developed code. It is shown that the fast multipole BEM can be applied to solve plate bending problems with good accuracy. Possible improvements in the efficiency of the method are discussed.     

Adaptive fast multipole boundary element method for three-dimensional half-space acoustic wave problems

A new adaptive fast multipole boundary element method (BEM) for solving 3-D half-space acoustic wave problems is presented in this paper. The half-space Green's function is employed explicitly in the boundary integral equation (BIE) formulation so that a tree structure of the boundary elements only for the boundaries of the real domain need to be applied, instead of using a tree structure that contains both the real domain and its mirror image. This procedure simplifies the implementation of the adaptive fast multipole BEM and reduces the CPU time and memory storage by about a half for large-scale half-space problems. An improved adaptive fast multipole BEM is presented for the half-space acoustic wave problems, based on the one developed recently for the full-space problems. This new fast multipole BEM is validated using several simple half-space models first, and then applied to model 3-D sound barriers and a large-scale windmill model with five turbines. The largest BEM model with 557470 elements was solved in about an hour on a desktop PC. The accuracy and efficiency of the BEM results clearly show the potential of the adaptive fast multipole BEM for solving large-scale half-space acoustic wave problems that are of practical significance.     

A fast algorithm for three-dimensional electrostatics analysis: fast Fourier transform on multipoles (FFTM)

In this paper, we propose a new fast algorithm for solving large problems using the boundary element method (BEM). Like the fast multipole method (FMM), the speed-up in the solution of the BEM arises from the rapid evaluations of the dense matrix–vector products required in iterative solution methods. This fast algorithm, which we refer to as fast Fourier transform on multipoles (FFTM), uses the fast Fourier transform (FFT) to rapidly evaluate the discrete convolutions in potential calculations via multipole expansions. It is demonstrated that FFTM is an accurate method, and is generally more accurate than FMM for a given order of multipole expansion (up to the second order). It is also shown that the algorithm has approximately linear growth in the computational complexity, implying that FFTM is as efficient as FMM. Copyright © 2004 John Wiley & Sons, Ltd.     

A fast multipole boundary element method for 2D multi-domain elastostatic problems based on a dual BIE formulation

A new fast multipole formulation for the hypersingular BIE (HBIE) for 2D elasticity is presented in this paper based on a complex-variable representation of the kernels, similar to the formulation developed earlier for the conventional BIE (CBIE). A dual BIE formulation using a linear combination of the developed CBIE and HBIE is applied to analyze multi-domain problems with thin inclusions or open cracks. Two pre-conditioners for the fast multipole boundary element method (BEM) are devised and their effectiveness and efficiencies in solving large-scale problems are discussed. Several numerical examples are presented to study the accuracy and efficiency of the developed fast multipole BEM using the dual BIE formulation. The numerical results clearly demonstrate the potentials of the fast multipole BEM for solving large-scale 2D multi-domain elasticity problems. The method can be applied to study composite materials, functionally-graded materials, and micro-electro-mechanical-systems with coupled fields, all of which often involve thin shapes or thin inclusions.     

Analytical integration of the moments in the diagonal form fast multipole boundary element method for 3-D acoustic wave problems

A diagonal form fast multipole boundary element method (BEM) is presented in this paper for solving 3-D acoustic wave problems based on the Burton-Miller boundary integral equation (BIE) formulation. Analytical expressions of the moments in the diagonal fast multipole BEM are derived for constant elements, which are shown to be more accurate, stable and efficient than those using direct numerical integration. Numerical examples show that using the analytical moments can reduce the CPU time by a lot as compared with that using the direct numerical integration. The percentage of CPU time reduction largely depends on the proportion of the time used for moments calculation to the overall solution time. Several examples are studied to investigate the effectiveness and efficiency of the developed diagonal fast multipole BEM as compared with earlier p³ fast multipole method BEM, including a scattering problem of a dolphin modeled with 404,422 boundary elements and a radiation problem of a train wheel track modeled with 257,972 elements. These realistic, large-scale BEM models clearly demonstrate the effectiveness, efficiency and potential of the developed diagonal form fast multipole BEM for solving large-scale acoustic wave problems.     

On the preconditioners for fast multipole boundary element methods for 2D multi-domain elastostatics

Fast multipole method (FMM) has been successfully applied to accelerate the numerical solvers of boundary element method (BEM). However, the coefficient matrix implicitly formed by using FMM is sometimes ill-conditioned in cases when mixed boundary conditions exist, resulting in poor rate of convergence for iteration. So preconditioning is a critical part in the development of efficient FMM solver for BEM. In this paper, preconditioners based on sparse approximate inverse type are used for fast multipole BEM to deal with 2D elastostatics. Several sparsity patterns of the preconditioner are considered for single- and multi-domain problems, especially for 2D elastic body with large number of inclusions or cracks. Algorithms and cost analysis of preconditioning under different prescribed sparsity patterns are discussed. GMRES is used as the iterative solver. Numerical results show this type of preconditioner achieves satisfactory rate of convergence for fast multipole BEM and performs well for problems of fairly large sizes.     

A fast multipole accelerated method of fundamental solutions for potential problems

The fast multipole method (FMM) is a very effective way to accelerate the numerical solutions of the methods based on Green's functions or fundamental solutions. Combined with the FMM, the boundary element method (BEM) can now solve large-scale problems with several million unknowns on a desktop computer. The method of fundamental solutions (MFS), also called superposition or source method and based on the fundamental solutions but without using integrals, has been studied for several decades along with the BEM. The MFS is a boundary meshless method in nature and offers more flexibility in modeling of a problem. It also avoids the singularity of the kernel by placing the source at some auxiliary points off the problem domain. However, like the traditional BEM, the conventional MFS also requires O(N²) operations to compute the system of equations and another O(N³) operations to solve the system using direct solvers, with N being the number of unknowns. Combining the FMM and MFS can potentially reduce the operations in formation and solution of the MFS system, as well as the memory requirement, all to O(N). This paper is an attempt in this direction. The FMM formulations for the MFS is presented for 2D potential problem. Issues in implementation of the FMM for the MFS are discussed. Numerical examples with up to 200,000 DOF's are solved successfully on a Pentium IV PC using the developed FMM MFS code. These results clearly demonstrate the efficiency, accuracy and potentials of the fast multipole accelerated MFS.     

A fast multipole boundary element method based on the improved Burton–Miller formulation for three-dimensional acoustic problems

A fast multipole boundary element method (FMBEM) based on the improved Burton–Miller formulation is presented in this paper for solving large-scale three-dimensional (3D) acoustic problems. Some improvements can be made for the developed FMBEM. In order to overcome the non-unique problems of the conventional BEM, the FMBEM employs the improved Burton–Miller formulation developed by the authors recently to solve the exterior acoustic problems for all wave numbers. The improved Burton–Miller formulation contains only weakly singular integrals, and avoids the numerical difficulties associated to the evaluation of the hypersingular integral, it leads to the numerical implementations more efficient and straightforward. In this study, the fast multipole method (FMM) and the preconditioned generalized minimum residual method (GMRES) iterative solver are applied to solve system matrix equation. The block diagonal preconditioner needs no extra memory and no extra CPU time in each matrix–vector product. Thus, the overall computational efficiency of the developed FMBEM is further improved. Numerical examples clearly demonstrate the accuracy, efficiency and applicability of the FMBEM based on improved Burton–Miller formulation for large-scale acoustic problems.     

A fast multipole boundary element method for solving two-dimensional thermoelasticity problems

A fast multipole boundary element method (BEM) for solving general uncoupled steady-state thermoelasticity problems in two dimensions is presented in this paper. The fast multipole BEM is developed to handle the thermal term in the thermoelasticity boundary integral equation involving temperature and heat flux distributions on the boundary of the problem domain. Fast multipole expansions, local expansions and related translations for the thermal term are derived using complex variables. Several numerical examples are presented to show the accuracy and effectiveness of the developed fast multipole BEM in calculating the displacement and stress fields for 2-D elastic bodies under various thermal loads, including thin structure domains that are difficult to mesh using the finite element method (FEM). The BEM results using constant elements are found to be accurate compared with the analytical solutions, and the accuracy of the BEM results is found to be comparable to that of the FEM with linear elements. In addition, the BEM offers the ease of use in generating the mesh for a thin structure domain or a domain with complicated geometry, such as a perforated plate with randomly distributed holes for which the FEM fails to provide an adequate mesh. These results clearly demonstrate the potential of the developed fast multipole BEM for solving 2-D thermoelasticity problems.     

A dual BIE approach for large‐scale modelling of 3‐D electrostatic problems with the fast multipole boundary element method

A dual boundary integral equation (BIE) formulation is presented for the analysis of general 3‐D electrostatic problems, especially those involving thin structures. This dual BIE formulation uses a linear combination of the conventional BIE and hypersingular BIE on the entire boundary of a problem domain. Similar to crack problems in elasticity, the conventional BIE degenerates when the field outside a thin body is investigated, such as the electrostatic field around a thin conducting plate. The dual BIE formulation, however, does not degenerate in such cases. Most importantly, the dual BIE is found to have better conditioning for the equations using the boundary element method (BEM) compared with the conventional BIE, even for domains with regular shapes. Thus the dual BIE is well suited for implementation with the fast multipole BEM. The fast multipole BEM for the dual BIE formulation is developed based on an adaptive fast multiple approach for the conventional BIE. Several examples are studied with the fast multipole BEM code, including finite and infinite domain problems, bulky and thin plate structures, and simplified comb‐drive models having more than 440 thin beams with the total number of equations above 1.45 million and solved on a PC. The numerical results clearly demonstrate that the dual BIE is very effective in solving general 3‐D electrostatic problems, as well as special cases involving thin perfect conducting structures, and that the adaptive fast multipole BEM with the dual BIE formulation is very efficient and promising in solving large‐scale electrostatic problems. Copyright © 2007 John Wiley & Sons, Ltd.     

A new fast multipole boundary element method for solving large‐scale two‐dimensional elastostatic problems

A new fast multipole boundary element method (BEM) is presented in this paper for large‐scale analysis of two‐dimensional (2‐D) elastostatic problems based on the direct boundary integral equation (BIE) formulation. In this new formulation, the fundamental solution for 2‐D elasticity is written in a complex form using the two complex potential functions in 2‐D elasticity. In this way, the multipole and local expansions for 2‐D elasticity BIE are directly linked to those for 2‐D potential problems. Furthermore, their translations (moment to moment, moment to local, and local to local) turn out to be exactly the same as those in the 2‐D potential case. This formulation is thus very compact and more efficient than other fast multipole approaches for 2‐D elastostatic problems using Taylor series expansions of the fundamental solution in its original form. Several numerical examples are presented to study the accuracy and efficiency of the developed fast multipole BEM formulation and code. BEM models with more than one million equations have been solved successfully on a laptop computer. These results clearly demonstrate the potential of the developed fast multipole BEM for solving large‐scale 2‐D elastostatic problems. Copyright © 2005 John Wiley & Sons, Ltd.     

An <i>O(N)</i> Fast Multipole Hybrid Boundary Node Method for 3D Elasticity

The Hybrid boundary node method (Hybrid BNM) is a boundary type meshless method which based on the modified variational principle and the Moving Least Squares (MLS) approximation. Like the boundary element method (BEM), it has a dense and unsymmetrical system matrix and needs to be speeded up while solving large scale problems. This paper combines the fast multipole method (FMM) with Hybrid BNM for solving 3D elasticity problems. The formulations of the fast multipole Hybrid boundary node method (FM-HBNM) which based on spherical harmonic series are given. The computational cost is estimated and an O(N) algorithm is obtained. The algorithm is implemented on a computer code written in C++. Numerical results demonstrate the accuracy and efficiency of the proposed technique.     

激光干涉仪反射镜三维温度场的快速多极边界元分析

温差造成光刻机激光干涉仪反射镜热变形,从而影响光刻机的精度。该文将基于二次单元的快速边界元法用于激光干涉仪反射镜的大规模温度场模拟。不连续单元的引入可以有效处理角点问题;新型快速多极算法用于边界元法的加速求解。建立统一的二次单元多极展开格式以处理混合边界。数值算例分析了快速多极边界元法的计算精度和效率,并和常规算法比较;使用该算法对激光干涉仪反射镜进行了大规模温度场计算,并和有限元法比较。结果表明:基于二次单元的快速多极边界元法可以高精度求解大规模三维传热问题。     

快速多极与常规边界元法机群并行计算的比较

以三维弹性力学问题为例,对快速多极与常规边界元法机群并行计算进行了比较。其中常规边界元法求解方程采用高斯消去法,通过调用标准并行求解函数库ScaLAPACK实现;快速多极边界元法并行计算程序采用ANSIC++语言、调用MPI并行通信库自行编写。两种程序均运行于同一机群并行环境。数值算例表明,在同样的机群条件下,采用快速多极边界元法可使解题规模有数量级的提高,计算速度明显高于常规边界元法,并行效率也优于常规边界元法。     

Application of fast multipole Galerkin boundary integral equation method to elastostatic crack problems in 3D

Fast multipole method (FMM) has been developed as a technique to reduce the computational cost and memory requirements in solving large‐scale problems. This paper discusses an application of FMM to three‐dimensional boundary integral equation method for elastostatic crack problems. The boundary integral equation for many crack problems is discretized with FMM and Galerkin's method. The resulting algebraic equation is solved with generalized minimum residual method (GMRES). The numerical results show that FMM is more efficient than conventional methods when the number of unknowns is more than about 1200 and, therefore, can be useful in large‐scale analyses of fracture mechanics. Copyright © 2001 John Wiley & Sons, Ltd.     

Dual BIE approaches for modeling electrostatic MEMS problems with thin beams and accelerated by the fast multipole method

Three boundary integral equation (BIE) formulations are investigated for the analysis of electrostatic fields exterior to thin-beam structures as found in some micro-electro-mechanical systems (MEMS). The three BIE formulations are: (1) the regular BIE using only the single-layer potential; (2) the dual BIE (a) using the regular BIE on one surface of a beam and the gradient BIE on the other surface; and (3) the dual BIE (b) using a linear combination of the regular BIE and gradient BIE on all the surfaces of the beam. Similar to crack problems in elasticity, the regular BIE degenerates when the beam thickness tends to zero, while the two dual BIE formulations do not degenerate. Most importantly, the dual BIE (b) is found to be well conditioned for all the values of the beam thickness, and thus well suited for implementation with the fast multipole BEM. The fast multipole BEM for both the regular BIE and the dual BIE (b) formulations are developed and tested on a simplified comb-drive model. The numerical results clearly show that the dual BIEs are very effective in solving MEMS problems with thin beams and the fast multipole BEM with the dual BIE (b) formulation is very efficient in solving large-scale MEMS models.     

Application of new fast multipole boundary integral equation method to crack problems in 3D

Fast multipole method (FMM) has been developed as a technique to reduce the computational cost and memory requirements in solving large scale problems. This paper discusses an application of the new version of FMM to three-dimensional boundary integral equation method (BIEM) for crack problems for the Laplace equation. The boundary integral equation is discretised with collocation method. The resulting algebraic equation is solved with generalised minimum residual method (GMRES). The numerical results show that the new version of FMM is more efficient than the original FMM.     

A fast multipole boundary element method for 3D multi-domain acoustic scattering problems based on the Burton-Miller formulation

A fast multipole boundary element method (FMBEM) for 3D multi-domain acoustic scattering problems based on the Burton-Miller formulation is presented in this paper. A multi-tree structure is designed for the multi-domain FMBEM. It results in mismatch of leaves and well separate cells definition in different domains and complicates the implementation of the algorithm, especially for preconditioning. A preconditioner based on boundary blocks is devised for the multi-domain FMBEM and its efficiency in reducing the number of iterations in solving large-scale multi-domain scattering problems is demonstrated. In addition to the analytical moment, another method, based on the anti-symmetry of the moment kernel, is developed to reduce the moment computation further by a factor of two. Frequency sweep analysis of a penetrable sphere shows that the multi-domain FMBEM based on the Burton-Miller formulation can overcome the non-unique solution problem at the fictitious eigenfrequencies. Several other numerical examples are presented to demonstrate the accuracy and efficiency of the developed multi-domain FMBEM for acoustic problems. In spite of the high cost of memory and CPU time for the multi-tree structure in the multi-domain FMBEM, a large BEM model studied with a PC has 0.3 million elements corresponding to 0.6 million unknowns, which clearly shows the potential of the developed FMBEM in solving large-scale multi-domain acoustics problems.     

基于有限元与宽频快速多极边界元的二维流固耦合声场分析

采用有限元/快速多极边界元法进行水下弹性结构的辐射和散射声场分析。Burton-Miller法用于解决传统单Helmholtz边界积分方程在求解外边界值问题时出现的非唯一解的问题。该文采用GMRES和快速多极算法加速求解系统方程。针对传统快速算法在高频处效率低和对角式快速算法在低频处不稳定这一问题,该文通过结合这两种快速算法形成宽频快速算法来克服。同时该文通过观察不同参数条件设置下,宽频快速多极法得到的数值结果在计算精度和计算时间上的变化,得到最优的参数组合值。最后通过数值算例验证该文算法的正确性和有效性。