An interpolation-based fast-multipole accelerated boundary integral equation method for the three-dimensional wave equation

Time-domain analysis of electromagnetic wave fields by boundary integral equation method

Computation of dynamic stress intensity factors by the time domain boundary integral equation method—II. Examples

Boundary integral equation method for transient analysis of 3-D cavities and inclusions

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

This paper discusses a formulation and its applications of the new Fast Multipole Method (FMM) to three-dimensional Boundary Integral Equation Method (BIEM) in elastostatic crack problems. It is shown, through numerical experiments, that the new FMM is more efficient than the original FMM.

A new hybrid finite-difference time-domain boundary integral equation (FDTD-BIE) method is proposed. The geometry under consideration is decomposed into several subregions. For some of the bounded subregions the FDTD technique is used to construct an interaction matrix, relating virtual electric and magnetic currents at the boundary of the subregion, for several frequencies at the same time. The remaining bounded and unbounded subregions are described by a boundary integral equation technique. An example is given to validate the new technique.

An enhancement of the fast time-domain boundary element method for the three-dimensional wave equation

In order to fulfill quantum information processing tasks such as quantum computation and scalable quantum networks, quantum logic gates are indispensable parts. A multiqubit logic gate can efficiently reduce the cost of physical resources and complicated operations that are necessary in the combination of fewer‐qubit (e.g., two‐qubit) logic gates to fulfill the same tasks. However, the construction of multiphoton quantum logic gates is a large challenge. Assisted by spatial and temporal degrees of freedom, a polarization multiphoton controlled one‐photon unitary gate is proposed, which can be realized with three processes including the control‐path process based on weak cross‐Kerr nonlinearities, the unitary transformation process, and the output process. In the control‐path process, the control‐path modules are applied, where the entanglement between polarization, spatial and temporal degrees of freedom is generated by beam splitters, time delayers, and nonlinear interaction in cross‐Kerr media. The scheme presented here, assisted by spatial and temporal degrees of freedom reduces the experimental complexity of the schemes assisted by only spatial degree of freedom. Moreover, it will also be a useful reference in optimally constructing multiqubit quantum logic gates with other multiple degrees of freedom in many physical systems.

One popular approach to solve the sea surface scattering and propagation in a ducting environment is the parabolic wave equation (PWE) method. An alternative method is the boundary integral equations (BIE) method. The implementation of the BIE in inhomogeneous media (ducting environments) is not straightforward, however, since the Green's function for such a medium is not usually known. In this paper, a closed-form approximation of the Green's function for a two-dimensional (2-D) ducting environment made up of a duct having a linear-square refractive index profile below a medium of constant refractive index, recently published, is used. This paper demonstrates how the BIE method can model the combined effects of surface roughness and medium inhomogeneity. Furthermore, it illustrates its capability of accurately predicting scattering in all directions including backscattering. Then, the PWE combined with the Split-Step Fourier (SSF) method is compared with this method.

In dislocation dynamics (DD) simulations, the most computationally intensive step is the evaluation of the elastic interaction forces among dislocation ensembles. Because the pair-wise interaction between dislocations is long-range, this force calculation step can be significantly accelerated by the fast multipole method (FMM). We implemented and compared four different methods in isotropic and anisotropic elastic media: one based on the Taylor series expansion (Taylor FMM), one based on the spherical harmonics expansion (Spherical FMM), one kernel-independent method based on the Chebyshev interpolation (Chebyshev FMM), and a new kernel-independent method that we call the Lagrange FMM. The Taylor FMM is an existing method, used in ParaDiS, one of the most popular DD simulation softwares. The Spherical FMM employs a more compact multipole representation than the Taylor FMM does and is thus more efficient. However, both the Taylor FMM and the Spherical FMM are difficult to derive in anisotropic elastic media because the interaction force is complex and has no closed analytical formula. The Chebyshev FMM requires only being able to evaluate the interaction between dislocations and thus can be applied easily in anisotropic elastic media. But it has a relatively large memory footprint, which limits its usage. The Lagrange FMM was designed to be a memory-efficient black-box method. Various numerical experiments are presented to demonstrate the convergence and the scalability of the four methods.

Though computational molecular dynamics is an effective tool for nano-scale phenomenon analysis, computational costs associated with its computer simulation are extremely high. There are two major computational steps associated with computer simulation of dynamics of molecular structures. They are calculation of interatomic forces and formation and solution of the equations of motion. Currently, these two computational steps are treated separately in most commonly-used methods. For example, Fast Multipole Method (FMM) and Cell Multipole Method (CMM) have been used for calculation of interatomic forces, and Cartesian Coordinate Method (CCM) and Internal Coordinate Molecular Dynamics Method (ICMD) have been created for the formation and solution of equations of motion of an atomistic molecular system. In this paper, a new procedure is presented through a proper integration between multibody molecular algorithms (MMA) and fast multipole methods to improve computational efficiency for computer simulation of the dynamical behaviors of multibody molecular structures in polymers and biopolymers. For the computational costs associated with interatomic forces, a fast multipole method is used to calculate the interatomic forces due to the potentials. For the computational costs associated with formation and solution of equations of motion, a multibody molecular algorithm developed by the author in his previous work will be utilized to integrate with fast multipole methods. The algorithm significantly improves computational efficiency when comparing with its counterpart procedures. The fast multipole method begins by scaling all atoms into a box with coordinate ranges to ensure numerical stability of subsequent operations. The parent box is then divided into half in the direction of each Cartesian axis and each child box is then subdivided to form a computational family tree. The flow of calculations is carried out along the tree structure with five passes. The fast multipole method has been improved and modified to achieve better effectiveness and higher efficiency since it was created. The multibody molecular algorithm starts with numbering subsets, forming bond graph, and developing three computing passes along the tree structure of an atomistic molecular system. Computing data flows in the fast multipole method and the multibody molecular algorithm will properly line up with the parent-child recursive relationship along the configuration of the tree structure due to linear recursive natures of both fast multipole method and multibody molecular algorithm. Then the time spent on the recursive simulation passes in the fast multipole method for computing forces may overlap with the time spent on the three recursive computational passes in the multibody molecular algorithm for forming and solving equations of motion.

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

Fast multipole singular boundary method for Stokes flow problems

Applying the concept of minimum potential energy and the variational method proposed in this paper, one can derive the governing equation and transversality conditions for the critical slip surface of a cohesive land slope described by simplified Janbu's model where both horizontal and vertical inter-slice forces are neglected. The governing equation, transversality conditions and boundary conditions were solved by the boundary integral equation method, which is a one dimensional BIEM, so that the critical slip surface and its associated minimal factor of safety can be determined effectively. By comparison of the results gotten from the boundary integral equation method and other numerical methods, it can be concluded, that, for some simplified cases, by using the boundary integral equation method on slope stability analysis of a cohesive land slope, a more reasonable result can be obtained.

A hybrid formulation of the 2.5D elastodynamic scattering problem combining the finite-element method and boundary integral equation method is presented and validated. The formulation of the 2.5D boundary integral equation method that is used was presented in detail by Papageorgiou and Pei (1998) and is an extension of the discrete wavenumber boundary integral equation method originally proposed by Kawase (1988) for 2D scattering problems. Modeling of the wave field in the domain of the scatterer is based on the variational principle of virtual displacements, and discretization of the domain is accomplished using the finite-element method.The formulation may be used to study the wave field in models of sedimentary deposits (e.g., valleys) or topography (e.g., canyons or ridges) with a 2D variation in structure but obliquely incident plane waves. The hybrid method exploits the versatility of the finite element method for modeling the scatterer and the effectiveness of the boundary integral equation method for taking care of the radiation condition in the half-space. The advantage of the 2.5D formulation is that it provides the means for calculations of 3D wave fields in scattering problems by requiring a storage comparable to that of the corresponding 2D calculations.