A fast multipole boundary element method for acoustics in viscothermal fluids

A fast multilevel multipole algorithm (FMM) is derived for the Helmholtz equation and adopted to the symmetric Galerkin boundary element method (BEM) for acoustics.The FMM allows to evaluate a matrix-vector product of the BEM with a computational cost of O(N log 2 N ), thus leading to a significant reduction of computation time and memory requirements compared to standard BEM formulations.This allows the simulation of large scale acoustic models.A coupling algorithm based on Lagrange multipliers is proposed for the simulation of structure-acoustic field interaction.Finite plate elements are coupled to a Galerkin boundary element formulation of the acoustic domain.The interface pressure is interpolated as a Lagrange multiplier, thus, allowing the coupling of non-matching grids.The resulting saddle point problem is solved by an approximate Uzawa-type scheme in which the matrix-vector products of the boundary element operators are evaluated efficiently by the fast multipole boundary element method.The algorithm is demonstrated on the example of a cavity backed elastic panel.

A modified dual-level fast multipole boundary element method based on the Burton–Miller formulation for large-scale three-dimensional sound field analysis

A numerical solution for the diffraction of linear waves around a vertical cylinder is presented in this paper using the fast multipole boundary element method (FMBEM). The fast multipole method (FMM) is employed to accelerate the traditional boundary element method. The accuracy and efficiency of the FMBEM for 2D Helmholtz equation is demonstrated, and it’s proved that this approach is suitable for solving large-scaled problems. In addition, numerical results of water dynamics coefficients, which are related to the diffraction problems of linear waves around a cylinder, are supplied in this paper.

The complex variable fast multipole boundary element method for the analysis of strongly inhomogeneous media

Finite Element Method (FEM1) is pervasively used in most of 3D product design analysis, in which Computer Aided Design (CAD) models need to be converted in to mesh models first and then enriched with some material features and boundary conditions data, etc. The interaction between CAD models and FEM models is intensive. Boundary Element Method (BEM) has been expected to be advantageous in large-scale problems in recent years owing to its reduction of the dimensionality and its reduced complexity in mesh generation. However, the BEM application has so far been limited to relatively small problems due to the memory and computational complexity for matrix buildup are O(N2). The fast multipole BEM (FMBEM) combined with BEM and fast multipole method (FMM) can overcome the defect of the traditional BEM, and provides an effective method to solve the large-scale problem. Combining GPU parallel computing with FMBEM can further improve its efficiency significantly. Based on the three-dimensional elastic mechanics problems, the parallelisms of the multipole moment (ME), multipole moment to multipole moment (M2M) translation, multipole moment to local expansion (M2L) translation, local expansion to local expansion (L2L) translation and near-field direct calculation were analyzed respectively according to the characteristics of the FMM, and the parallel strategies under CUDA were presented in this paper. Three main major parts are included herein: (1) FMBEM theory in 3D elastostatics, (2) the parallel FMBEM algorithm using CUDA, and (3) comparison the GPU parallel FMBEM with BEM, FEM and FMBEM respectively by engineering examples. Numerical example results show the 3D elastostatics GPU FMBEM not only can speed up the boundary element calculation process, but also save memory which can be effective to solve the large-scale engineering problems.

Fast multipole boundary element method using the binary tree structure with tight bounds: application to a calculation of an electrostatic force for the manipulation of a metal micro particle ☆

It is suitable to solve the acoustic problems by using the fast multipole boundary element method (FMBEM), since the FMBEM can accelerate the matrix-vector multiplication dramatically by reducing the CPU time and memory of conventional boundary element method to O(N log2N) and O(N) respectively. We propose a 3D acoustic FMBEM based on Burton-Miller formulation in this paper. A new adaptive algorithm is applied to the diagonal form FMBEM, and a new proposed analytical moment formulation is used in the moment computation. Both of them further improve the efficiency of FMBEM. Acoustic scattering of soft sphere at resonant frequency is investigated to validate the accuracy of solution using Burton-Miller formulation. Comparisons of solution to the multi-radiating spheres problem with the one solved by Bapat's program demonstrate the accuracy and the efficiency of our algorithm in solving large-scale acoustic problems. In the end, we use our algorithm to analyze the inner sound filed of a car and dolphin acoustic scattering.

Eddy current nondestructive evaluation (NDE) of airframe structures involves the detection of electromagnetic field irregularities due to non-conducting inhomogeneities in an electrically conducting material. Usually, the eddy current NDE problem can be formulated by the boundary integral equations (BIE) and discretized into matrix equations by the method of moments (MoM) or the boundary element method (BEM). The fast multipole method (FMM) is a well-established and effective method for accelerating numerical solutions of the matrix equations. Accelerated by the FMM, the BIE method can now solve large-scale electromagnetic wave propagation and diffusion problems. The traditional BIE method 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; in contrast, the BIE method accelerated by the two-level FMM can potentially reduce the operations and memory requirement to O(N). Moreover, several approaches have been proposed for the field calculation in the presence of flaws in three dimensional NDE; however, seldom work has been done in applying efficient methods to seek rapid solution in eddy current NDE simulation. As elaborated in the dissertation, we introduce a fast multipole BIE method for two-dimensional diffusive scalar problem and an efficient BIE method for three-dimensional eddy current NDE. Firstly, we work with the two-dimensional Helmholtz equation with a complex wave number for non-trivial boundary geometry. We describe the FMM acceleration procedure of the BIE method and its features briefly, explaining that the FMM is not only efficient in meshing complicated geometries, accurate for solving singular fields or

An adaptive dual-information FMBEM for 3D elasticity and its GPU implementation

This paper describes an application of the fast multipole boundary element method (FMBEM) in simulating the viscoelastic behavior of asphalt concrete (AC). FMBEM has some merits including easy meshing of complicated geometries, accuracy in solving singular fields, high calculation efficiency, and low data storage requirements. These merits lead to many applications in multiple time-step and multi-inclusion problems, such as prediction of AC creep behavior. In this paper, a fast multipole formulation and an algorithm for two-dimensional linear viscoelastic problems are introduced first. To consider the bonding between the asphalt matrix and aggregates, a viscoelastic interface model is used to simulate the interfacial imperfection. A numerical model of asphalt concrete is then constructed based on image processing techniques, and the model is capable of taking into account the effects of the size and content of the air voids. The creep compliance of AC is predicted by the developed FMBEM arithmetic and subsequently compared with the experimental results. It has been shown that the developed method can describe and simulate the correct creep evolutionary trend of AC. DOI: 10.1061/(ASCE)MT.1943- 5533.0000603. © 2013 American Society of Civil Engineers. CE Database subject headings: Asphalts; Concrete; Micromechanics; Viscoelasticity; Imaging techniques; Boundary element method. Author keywords: Asphalt concrete; Micromechanics; Fast multipole boundary element method; Viscoelastic properties; Image processing technique.

Abstract: In this paper, a half space adaptive fast multipole boundary facemethod (FMBFM) is presented for solving the three-dimensional half space ex-terior acoustic problems. In the presented method, the Burton-Miller equationbased on the conventional boundary integral equation (CBIE) and its hyper-singularboundary integral equation (HBIE) is used to deal with the fictitious eigenfrequen-cies problem. The half space Green’s function is employed, thus the tree structurein the fast multipole method can be used only for the real domain. The higher orderelements and an adaptive tree structure are used to improve the efficiency of theFMBFM. This half space adaptive FMBFM for half space acoustic problems is anextension of the adaptive FMBFM for full space acoustic problems developed bythe authors. Numerical examples for half space acoustic problems in this paperdemonstrate the efficiency and validity of this method.Keywords: fast multipole boundary face method, Burton-Miller equation, acous-tic problems, half space, modified hyper-singular boundary integral equation1 IntroductionAcoustic problems in half space are a major subject in engineering. Many researchworks have been done on this aspect (Seybert et al., 1988; Secgin et al., 2010; Caoet al., 2010). The half space Green’s function (Seybert et al., 1988) can be appliedby adding the source solution at the image point to the original full-space Green’sfunction for dealing with half space problems. Sound source localization and iden-tification of a refrigerator were accomplished by using the half space boundaryelement method (BEM) (Secgin et al., 2010). Using the half space Green’s func-tion, the discretization of the rigid infinite plane is removed and only the boundaries

The computation of the acoustic scattering from objects in the ocean is useful to analyze the characteristics of the targets, which can be numerically modelled by acoustic boundary element method (BEM). The BEM has a defect of fictitious frequency problem, Burton-Miller method is studied to overcome the singular problems. In this paper, the Hadamard finite-part integral method based on collocation method is studied to overcome the singularity. It is found that, the near singularity is also quite significant and deteriorates largely the accuracy of BEM. A method based on polar transformation is proposed in this paper to overcome such a defect. Once the problems of singularity and near-singularity are solved, the conventional BEM model is established to serve as a foundation for fast BEM.It is well known that conventional BEM yields a full matrix, which is of O(N <sup xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">2</sup> ) memory complexity. Furthermore, the time complexity could be up to O(N <sup xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">3</sup> ), if Gaussian elimination method is adopted to solve the linear equation of BEM. Many fast methods are proposed to solve the problems, such as the GMRES iterative solver, fast multipole method and H-matrix (Hierarchical Matrices) method. H-matrix method is applied to overcome the faults of conventional BEM, by using the low-rank approximations to compress blocks of the matrix, which is a pure algebraic method. The ACA (Adaptive Cross Approximation) method is proposed to compute the outer-product form of low-rank matrices, instead of SVD (Singular Value Decomposition), which consumes high CPU time and is of no usability. ACA-BEM is used to predict the acoustic scattering from objects in the ocean, especially the large-scale problems. The directivities of a rigid sphere and a rigid ellipsoid are computed by ACA BEM, the time consumed is very short. For a problem of more than 20,000 DOFs, it can be solved within 2 minutes. Therefore, ACA BEM can be used to solve large-scale acoustic scattering problems very efficiently.

Boundary element method is one of the most powerful tools to solve partial differential equations. However, it takes a lot of computational time when applied to a problem with a large number of unknowns. In this paper, we adopted the binary tree structure with tight bounds and the downward pass for that structure for the fast multipole method to reduce the number of M2L translations. The fast multipole boundary element method using the binary tree structure with tight bounds is developed and the analyses of an electrostatic force acting on a matal micro particle above a complex structure to confirm the effectiveness of this method.

—By employing the analysis of storage and computation time in fast multi-pole algorithm, this article introduces a new fast multi-pole algorithm to solve the problem of large computation time in multi-pole transfer process of the original fast multi-pole boundary element method. The new method uses the index expansion form to complete the transfer process of multi-pole moments at the expense of a small amount of storage. The obtained computational efficiency is significantly improved. Finally, the two-dimensional plate capacitor model is used to verify the precision of the authors’ method, and the three-dimensional spherical capacitor model is used to verify the superiority of the new fast multi-pole algorithm in solving multi-freedom electric field distribution problems. The electric field distribution of a 110 kV substation switching field model is calculated by using the new fast multi-pole algorithm. Compared with the calculated and measured values, the error is less than 10% upon optimization of the model; therefore, this method is better able to solve large-scale, multi degree-of-freedom problems.

A fast multipole boundary element method for 2D viscoelastic problems