A fast multipole boundary element method for 2D viscoelastic problems

A fast multipole formulation for 2D linear viscoelastic problems is presented in this paper by incorporating the elastic–viscoelastic correspondence principle. Systems of multipole expansion equations are formed and solved analytically in Laplace transform domain. Three commonly used viscoelastic models are introduced to characterize the time-dependent behavior of the materials. Since the transformed multipole formulations are identical to those for the 2D elastic problems, it is quite easy to implement the 2D viscoelastic fast multipole boundary element method. Besides, all the integrals are evaluated analytically, leading to highly accurate results and fast convergence of the numerical scheme. Several numerical examples, including planar viscoelastic composites with single inclusion or randomly distributed multi-inclusions, as well as the problem of a crack in a pressured viscoelastic plane, are presented. The results are verified by comparison with the developed analytical solutions to illustrate the accuracy and efficiency of the approach.     

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 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 multipole hybrid boundary node method for composite materials

This article presents a multi-domain fast multipole hybrid boundary node method for composite materials in 3D elasticity. The hybrid boundary node method (hybrid BNM) is a meshless method which only requires nodes constructed on the surface of a domain. The method is applied to 3D simulation of composite materials by a multi-domain solver and accelerated by the fast multipole method (FMM) in this paper. The preconditioned GMRES is employed to solve the final system equation and precondition techniques are discussed. The matrix–vector multiplication in each iteration is divided into smaller scale ones at the sub-domain level and then accelerated by FMM within individual sub-domains. The computed matrix–vector products at the sub-domain level are then combined according to the continuity conditions on the interfaces. The algorithm is implemented on a computer code written in C + +. Numerical results show that the technique is accurate and efficient.     

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.     

A wideband fast multipole boundary element method for three dimensional acoustic shape sensitivity analysis based on direct differentiation method

This paper presents a wideband fast multipole boundary element approach for three dimensional acoustic shape sensitivity analysis. The Burton-Miller method is adopted to tackle the fictitious eigenfrequency problem associated with the conventional boundary integral equation method in solving exterior acoustic wave problems. The sensitivity boundary integral equations are obtained by the direct differentiation method, and the concept of material derivative is used in the derivation. The iterative solver generalized minimal residual method (GMRES) and the wideband fast multipole method are employed to improve the overall computational efficiency. Several numerical examples are given to demonstrate the accuracy and efficiency of the present method.     

A multipole expansion-based boundary element method for axisymmetric potential problem

The multipole expansion is an approximation technique used to evaluate the potential field due to sources located in the far field. Based on the multipole expansion, we describe a new technique to calculate the far potential field due to ring sources which are encountered in the boundary element method (BEM) formulation of axisymmetric problems. As the sources in the near field are processed by the slower conventional BEM, it is important to maximize the amount of multipole calculations taking advantage of both interior and exterior multipole expansions. Numerical results are presented for an axisymmetric potential test problem with Neumann and Dirichlet boundary conditions. The complexity of the proposed method remains O(N²), which is equal to that of the conventional BEM. However, the proposed technique coupled with an iterative solver speeds up the solution procedure. The technique is significantly advantageous when medium and large numbers of elements are present in the domain.     

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 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.     

A low-frequency fast multipole boundary element method based on analytical integration of the hypersingular integral for 3D acoustic problems

A low-frequency fast multipole boundary element method (FMBEM) for 3D acoustic problems is proposed in this paper. The FMBEM adopts the explicit integration of the hypersingular integral in the dual boundary integral equation (BIE) formulation which was developed recently by Matsumoto, Zheng et al. for boundary discretization with constant element. This explicit integration formulation is analytical in nature and cancels out the divergent terms in the limit process. But two types of regular line integrals remain which are usually evaluated numerically using Gaussian quadrature. For these two types of regular line integrals, an accurate and efficient analytical method to evaluate them is developed in the present paper that does not use the Gaussian quadrature. In addition, the numerical instability of the low-frequency FMBEM using the rotation, coaxial translation and rotation back (RCR) decomposing algorithm for higher frequency acoustic problems is reported in this paper. Numerical examples are presented to validate the FMBEM based on the analytical integration of the hypersingular integral. The diagonal form moment which has analytical expression is applied in the upward pass. The improved low-frequency FMBEM delivers an algorithm with efficiency between the low-frequency FMBEM based on the RCR and the diagonal form FMBEM, and can be used for acoustic problems analysis of higher frequency.     

A fast elasto-plastic formulation with hierarchical matrices and the boundary element method

Boundary element methods offer some advantages for the simulation of tunnel excavation since the radiation condition is implicitly fulfilled and only the excavation and ground surfaces have to be discretized. Hence, large meshes and mesh truncation, as required in the finite element method, are avoided. Recently, capabilities for efficiently dealing with inelastic behavior and ground support have been developed, paving the way for the use of the method to simulate tunneling. However, for large scale three dimensional problems one drawback of the boundary element method becomes prominent: the computational effort increases quadratically with the problem size. To reduce the computational effort several fast methods have been proposed. Here a fast boundary element solution procedure for small strain elasto-plasticity based on a collocation scheme and hierarchical matrices is presented. The method allows the solution of problems with the computational effort and sparse storage increasing almost linearly with respect to the problem size.     

An adaptive fast multipole boundary face method with higher order elements for acoustic problems in three-dimension

An adaptive fast multipole boundary face method using higher order elements based on the well-known Burton-Miller equation is presented in this paper for solving the large-scale three-dimensional exterior acoustic wave problems. The fast multipole boundary face method is referred to as FMBFM. In the FMBFM, the boundary integration and field variables approximation are both performed in the parametric space of each boundary face exactly the same as the B-rep data structure in standard solid modeling packages. In this FMBFM, higher order elements are employed to improve the computational accuracy and efficiency, and an adaptive tree structure is constructed to improve the efficiency of the FMBFM. Numerical examples for large-scale acoustic radiation and scattering problems in this paper demonstrated the accuracy, efficiency and validity of the adaptive FMBFM. Comparison study showed that the FMBFM with high order elements out-performs the FMBFM with constant elements respect to accuracy and CPU time at the same number of the nodes. In addition, the CAD models, even with complicated geometry, are directly converted into the FMBFM models, thus the present method provides a new way toward automatic simulation.     

Acoustic problems analysis of 3D solid with small holes by fast multipole boundary face method

In this paper, an adaptive fast multipole boundary face method is introduced to implement acoustic problems analysis of 3D solids with open-end small tubular shaped holes. The fast multipole boundary face method is referred as FMBFM. These holes are modeled by proposed tube elements. The hole is open-end and intersected with the outer surface of the body. The discretization of the surface with circular inclusions is achieved by applying several special triangular elements or quadrilateral elements. In the FMBFM, the boundary integration and field variables approximation are both performed in the parametric space of each boundary face exactly the same as the B-rep data structure in standard solid modeling packages. Numerical examples for acoustic radiation in this paper demonstrated the accuracy, efficiency and validity of this method.