Fast multipole accelerated boundary element method for the Helmholtz equation in acoustic scattering problems

We apply the fast multipole method (FMM) accelerated boundary element method (BEM) for the three-dimensional (3D) Helmholtz equation, and as a result, large-scale acoustic scattering problems involving 400000 elements are solved efficiently. This is an extension of the fast multipole BEM for two-dimensional (2D) acoustic problems developed by authors recently. Some new improvements are obtained. In this new technique, the improved Burton-Miller formulation is employed to over-come non-uniqueness difficultie...     

新的对角形式快速多极边界元法求解声学Helmholtz方程

传统外部声学Helmholtz边界积分方程无法在个人计算机上求解大规模工程问题. 为了有效解决这个问题, 将快速多极方法引入到边界积分方程中, 加速系统矩阵方程组的迭代求解. 由于在边界积分方程中引入基本解的对角形式多极扩展, 新的快速多极边界元法的计算效率与传统边界元相比显著提高, 计算量和存储量减少到O(N)量级(N为问题的自由度数). 包括含有420000个自由度的大型潜艇模型数值算例验证了快速多极边界元法的准确性和高效性, 清楚表明新算法在求解大规模声学问题中的优势, 具有良好的工程应用前景.     

新的对角形式快速多极边界元法求解声学Helmholtz方程

传统外部声学Helmholtz边界积分方程无法在个人计算机上求解大规模工程问题. 为了有效解决这个问题, 将快速多极方法引入到边界积分方程中, 加速系统矩阵方程组的迭代求解. 由于在边界积分方程中引入基本解的对角形式多极扩展, 新的快速多极边界元法的计算效率与传统边界元相比显著提高, 计算量和存储量减少到O(N)量级(N为问题的自由度数). 包括含有420000个自由度的大型潜艇模型数值算例验证了快速多极边界元法的准确性和高效性, 清楚表明新算法在求解大规模声学问题中的优势,     

Transparent boundary conditions for the Helmholtz equation in some ramified domains with a fractal boundary

The paper addresses a class of boundary value problems in some self-similar ramified domains, with the Laplace or Helmholtz equations. Much stress is placed on transparent boundary conditions which allow the solutions to be computed in subdomains. A self similar finite element method is proposed and tested. It can be used for numerically computing the spectrum of the Laplace operator with Neumann boundary conditions, as well as the eigenmodes. The eigenmodes are normalized by means of a perturbation method and the spectral decomposition of a compactly supported function is carried out. Finally, a numerical method for the wave equation is addressed.     

A scaled boundary node method applied to two-dimensional crack problems

A boundary-type meshless method called the scaled boundary node method(SBNM) is developed to directly evaluate mixed mode stress intensity factors(SIFs) without extra post-processing.The SBNM combines the scaled boundary equations with the moving Kriging(MK) interpolation to retain the dimensionality advantage of the former and the meshless attribute of the latter.As a result,the SBNM requires only a set of scattered nodes on the boundary,and the displacement field is approximated by using the MK interpolation technique,which possesses the δ function property.This makes the developed method efficient and straightforward in imposing the essential boundary conditions,and no special treatment techniques are required.Besides,the SBNM works by weakening the governing differential equations in the circumferential direction and then solving the weakened equations analytically in the radial direction.Therefore,the SBNM permits an accurate representation of the singularities in the radial direction when the scaling center is located at the crack tip.Numerical examples using the SBNM for computing the SIFs are presented.Good agreements with available results in the literature are obtained.     

势问题的径向基函数——局部边界积分方程方法

将基于径向基函数构造的具有插值特性的近似函数和局部边界积分方程方法相结合，建立了求解势问题的径向基函数——局部边界积分方程方法，推导了相应离散方程.与其他边界积分方程的无网格方法相比，本文方法具有数值实现过程简单、计算量小、精度高的优点，并可直接施加边界条件.最后通过算例说明了该方法的有效性. 关键词： 径向基函数 无网格方法 局部边界积分方程 势问题     

二维Helmholtz边界超奇异积分方程解析研究

基于常规边界元法及超奇异边界积分方程复线性耦合的Burton-Miller方法应用于无限域声学问题的最大难点在于处理超奇异积分（二维问题）.目前,此类超奇异积分主要使用各种弱奇异/正则化方法求解,而这些弱奇异/正则化方法具有时间消耗大等弱点.基于围道积分定理,本文给出一种使用常值单元的二维Helmholtz边界超奇异积分的解析表达式.在有限部分积分意义下,所有的奇异和超奇异积分可以解析表达.数值算例表明该解析表达式是有效的.     

A moving Kriging interpolation-based boundary node method for two-dimensional potential problems

In this paper,a meshfree boundary integral equation(BIE) method,called the moving Kriging interpolationbased boundary node method(MKIBNM),is developed for solving two-dimensional potential problems.This study combines the BIE method with the moving Kriging interpolation to present a boundary-type meshfree method,and the corresponding formulae of the MKIBNM are derived.In the present method,the moving Kriging interpolation is applied instead of the traditional moving least-square approximation to overcome Kronecker’s delta property,then the boundary conditions can be imposed directly and easily.To verify the accuracy and stability of the present formulation,three selected numerical examples are presented to demonstrate the efficiency of MKIBNM numerically.     

开口/封闭薄壳体声辐射和散射的统一边界积分方程解法

推导了在二维和三维空间下开口和封闭薄壳体在任意阻抗边界条件下声辐射和散射的统一边界积分方程.相对于以前的求解方法,该方程求解声辐射和散射问题具有相同的影响矩阵,能够同时求解薄壳体气动和振动噪声的辐射和散射现象,以及分析壳体声阻抗对声波传播的影响.推导的方程可以应用于叶轮机械、管道等噪声和消声器消声性能的预测等方面.在此方程基础上,可以进一步考虑运动边界和运动介质对声辐射和散射的影响. 关键词： 薄壳体 声阻抗 积分方程 边界元方法     

Coupling of Dirichlet-to-Neumann boundary condition and finite difference methods in curvilinear coordinates for multiple scattering

The applicability of the Dirichlet-to-Neumann technique coupled with finite difference methods is enhanced by extending it to multiple scattering from obstacles of arbitrary shape. The original boundary value problem (BVP) for the multiple scattering problem is reformulated as an interface BVP. A heterogenous medium with variable physical properties in the vicinity of the obstacles is considered. A rigorous proof of the equivalence between these two problems for smooth interfaces in two and three dimensions for any finite number of obstacles is given. The problem is written in terms of generalized curvilinear coordinates inside the computational region. Then, novel elliptic grids conforming to complex geometrical configurations of several two-dimensional obstacles are constructed and approximations of the scattered field supported by them are obtained. The numerical method developed is validated by comparing the approximate and exact far-field patterns for the scattering from two circular obstacles. In this case, for a second order finite difference scheme, a second order convergence of the numerical solution to the exact solution is easily verified.     

A kernel-free boundary integral method for elliptic boundary value problems

This paper presents a class of kernel-free boundary integral (KFBI) methods for general elliptic boundary value problems (BVPs). The boundary integral equations reformulated from the BVPs are solved iteratively with the GMRES method. During the iteration, the boundary and volume integrals involving Green’s functions are approximated by structured grid-based numerical solutions, which avoids the need to know the analytical expressions of Green’s functions. The KFBI method assumes that the larger regular domain, which embeds the original complex domain, can be easily partitioned into a hierarchy of structured grids so that fast elliptic solvers such as the fast Fourier transform (FFT) based Poisson/Helmholtz solvers or those based on geometric multigrid iterations are applicable. The structured grid-based solutions are obtained with standard finite difference method (FDM) or finite element method (FEM), where the right hand side of the resulting linear system is appropriately modified at irregular grid nodes to recover the formal accuracy of the underlying numerical scheme. Numerical results demonstrating the efficiency and accuracy of the KFBI methods are presented. It is observed that the number of GMRES iterations used by the method for solving isotropic and moderately anisotropic BVPs is independent of the sizes of the grids that are employed to approximate the boundary and volume integrals. With the standard second-order FEMs and FDMs, the KFBI method shows a second-order convergence rate in accuracy for all of the tested Dirichlet/Neumann BVPs when the anisotropy of the diffusion tensor is not too strong.     

三维声学多层快速多极子边界元及其应用

快速多极子边界元算法可以加速矩阵和向量乘法运算, 将传统边界元算法的计算量和内存占用量分别降为O(N log²N)和O(N), 适用于大型声学模型模拟计算. 本文发展了一种基于Burton-Miller方程的三维多层声学快速多极子边界元算法. 将新的自适应树状算法应用到对角形式的快速多极子边界元算法, 并使用最新提出的解析式源点矩计算公式, 进一步提高了快速多极子边界元的计算效率. 绝对软球体在内部共振频率处的散射声场计算, 验证了所发展算法在共振频率处求解的正确性. 与Bapat所提供的程序在多脉动球体辐射声场计算精度的比较, 验证了算法及程序在大型模型声学计算中的准确性, 同时显示了其求解的高效性. 最后, 将该算法用于车内声场及水下声学探测的分析计算.     

三维Helmholtz方程外问题的自然积分方程及其数值解

用文[2,3]提出的自然边界归化方法来处理三维Helmholtz方程的外边值问题。在简要介绍如何用球谐展开的方法得到Helmholtz问题在外球域上的自然积分方程后,给出求解该自然积分方程的一种数值方法及相应的数值算例。     

Application of the extended traction boundary element-free method to the fracture of two-dimensional infinite magnetoelectroelastic solid

A novel extended traction boundary element-free method is proposed to analyze the crack problems of two-dimensional infinite magnetoelectroelastic solid.An extended traction boundary integral equation only involving Cauchy singularity is firstly derived.Then,the extended dislocation densities on the crack surface are expressed as the combination of a characteristic term and unknown weight functions,and the radial point interpolation method is adopted to approximate the unknown weight functions.The numerical...     

导电平板上任意缝隙填充各向异性介质时TM波散射和传输特性分析

简单介绍了用于分析任意口径问题的通用频域方法——广义网络原理,并利用边界积分法和广义网络原理分析导电平板上任意缝隙填充各向异性介质时TM波的散射及传输特性.由于缝隙填充各向异性介质的情形尚未见公开文献报道,作为验证,将本方法退化计算各向同性介质填充时缝隙的散射和传输特性,并与文献结果进行比较.最后,给出了缝隙填充各向异性介质时的算例.     

On the indefinite Helmholtz equation: Complex stretched absorbing boundary layers,iterative analysis,and preconditioning

This paper studies and analyzes a preconditioned Krylov solver for Helmholtz problems that are formulated with absorbing boundary layers based on complex coordinate stretching. The preconditioner problem is a Helmholtz problem where not only the coordinates in the absorbing layer have an imaginary part, but also the coordinates in the interior region. This results into a preconditioner problem that is invertible with a multigrid cycle. We give a numerical analysis based on the eigenvalues and evaluate the performance with several numerical experiments. The method is an alternative to the complex shifted Laplacian and it gives a comparable performance for the studied model problems.     

Some new exact solutions to the Burgers--Fisher equation and generalized Burgers--Fisher equation

Some new exact solutions of the Burgers--Fisher equation and generalized Burgers--Fisher equation have been obtained by using the first integral method. These solutions include exponential function solutions, singular solitary wave solutions and some more complex solutions whose figures are given in the article. The result shows that the first integral method is one of the most effective approaches to obtain the solutions of the nonlinear partial differential equations.     

Simulating regular wave propagation over a submerged bar by boundary element method model and Boussinesq equation model

Numerical models based on the boundary element method and Boussinesq equation are used to simulate the wave transform over a submerged bar for regular waves.In the boundary-element-method model the linear element is used,and the integrals are computed by analytical formulas.The Boussinesq-equation model is the well-known FUNWAVE from the University of Delaware.We compare the numerical free surface displacements with the laboratory data on both gentle slope and steep slope,and find that both the models simulate the wave transform well.We further compute the agreement indexes between the numerical result and laboratory data,and the results support that the boundary-element-method model has a stable good performance,which is due to the fact that its governing equation has no restriction on nonlinearity and dispersion as compared with Boussinesq equation.     

模糊和散射对闪光照相边界检测的影响

借助蒙特卡罗方法模拟闪光照相过程,研究了模糊和散射因素对边界检测的影响。结果表明：模糊与边界退化量的变化关系依赖于相邻飞层的材料、密度和尺寸;相同客体结构下,模糊越大,边界退化量也越大;客体结构不同,相同模糊下的边界退化量也不同;散直比越大,散射对边界检测结果的影响也越大,当散直比大于50时,散射引起的边界检测误差超过1个像素。