三维位势问题边界元法中几乎奇异积分的正则化

采用一种半解析正则化算法,计算了三维位势问题边界元法中近边界点的几乎强奇异和几乎超奇异面积分.该算法适用于三角形线性等参元.对高次单元将其细分为几个三节点三角形单元即可应用该算法.由于几乎奇异性,与内点邻近的单元上的积分,采用半解析正则化积分算法计算;而远处单元的积分仍保持常规高斯积分.对三维热传导算例,计算了近边界点的温度和热流.数值结果证明了该算法的有效性和精确性.     

各向异性介质三维电磁响应模拟的Ho-GEBA算法

本文基于积分方程法研究并建立了一种模拟横向同性介质中任意各向异性异常 体三维电磁响应的高阶广义扩展Born近似(Ho-GEBA)算法. 首先利用逐次迭代技术给出积分方程的广义级数展开解, 为保证其收敛性, 引入一种各向异性条件下满足压缩映射的迭代算子. 然后利用异常体区域分解技术, 并结合扩展Born近似原理, 得到各向异性介质三维电磁响应的Ho-GEBA解. 为提高效率, 计算过程中采用并矢Green函数的解析表达式. 最后通过数值计算实例对比验证了本文算法的有效性. 关键词： 高阶广义扩展Born近似 积分方程 电磁模拟 解析Green函数     

三维静电场线性插值边界元中的解析积分方法

提出求解三维静电场的三角形线性插值边界元解析积分方法.针对含1/R和1/R²的积分项,将单元形状函数分解为常数项、含x的线性项和含y的线性项,从而将边界单元积分简化为6个基本积分组合,并导出其解析计算公式,避免了因形状函数改变而导致的重复计算.该方法不仅可以准确计算远离奇异情况下的边界元积分,而且可以准确计算一阶和二阶接近奇异积分以及一阶奇异积分.计算结果表明,在接近奇异积分和奇异积分比较突出的问题中,当数值积分方法不能给出正确结果时,用同样的边界元网格,解析积分方法可以给出正确的结果,提高了三维静电场线性插值边界元法的计算精度.     

声学时域边界元法的时间卷积计算方法研究

核函数中保留Dirac函数的原型,形成关于时间的卷积积分,是声学时域边界元法中一种稳定、有效的时间数值积分计算方法 (CQ-BEM)。然而,传统CQ-BEM中卷积积分系数的获取有计算量大、耗时长,且对不同单元需要重新计算的问题,极大地降低了CQ-BEM法计算时域声场的效率。针对传统CQ-BEM积分系数计算效率低的问题,本文利用多项式展开定理给出了待求函数泰勒系数的解析表达与数值计算方法,建立了不同单元间待求系数的转换理论,可以在一次循环迭代内完成不同单元的积分系数的计算,大幅降低了计算量,提高了CQ-BEM方法的声场计算效率。脉动球源数值算例结果表明,在相同要求下,本文方法计算时间较传统方法减少50%以上,相对误差小5个数量级以上,且计算时间随单元数的增长率仅为传统方法的2.34%。因此,本文提出的系数计算方法能够有效提高CQ-BEM方法的时域声场计算效率,拓展了CQ-BEM在大型机电设备时域声场模拟的计算规模。     

采用能量守恒和高阶Padé近似的三维水声抛物方程模型

为了充分考虑海底地形随三维空间变化的海洋环境中水平方位角耦合效应对声传播的影响,建立了一种三维柱坐标系下流体高阶抛物方程算法。该算法采用泰勒近似将二维方根算子分裂成一维方根算子,并采用分裂步进的高阶Pade近似将一维方根算子写成微分算子有理分式连乘的形式,进而应用Galerkin离散化方法来处理微分算子,最终将微分方程写成矩阵方程的形式;采用能量守恒近似来处理海底边界,以考虑复杂海底对于声传播的影响;采用交替方向隐式格式,实现了三维声场的步进计算。楔形和海底山等典型海域声场仿真计算表明,相比于已有的声场计算模型,三维柱坐标系下高阶抛物方程模型可以更加精确地计算楔形海域和海底山区域的三维声场,实现水平方位全空间声场计算。     

基于自适应单元细分法的高效高精度近奇异域积分计算

本文针对边界元法在计算薄型结构力学、裂纹扩展等物理问题时存在的积分难题,提出一种基于自适应单元细分法的高效高精度近奇异域积分计算方法,该方法基于二叉树数据结构的单元细分技术对体单元进行自适应细分,消除单元几何形状所引起的近奇异性,能直接用于计算连续核函数的近奇异域积分。针对间断核函数的近奇异域积分,在细分单元的基础上采用腔面重建算法和投影算法,重新构建源点附近的积分子单元。数值算例表明：本方法可采用较少的积分点得到准确结果,是处理近奇异域积分的一种有效方法。     

表面电荷法平面元电荷系数的半解析解法

针对在线性及高次电荷估计下,表面电荷法中全解析法对平面元电荷系数求解和实现较复杂的问题,提出一种半解析法.将系数积分由局部坐标系变换到整体坐标系下分离参数的二重积分,利用内层积分存在解析解的特点,将二维问题降为一维,方便数值积分计算.对于对数奇异积分问题采用辅助函数法消除.通过算例与全解析法计算精度进行比较,结果表明,在一定的单元划分下,可达到很高精度,具有一定可行性.     

三维随机粗糙面上导体目标散射的解析-数值混合算法

提出三维导体目标与导体粗糙面复合散射的解析-数值混合迭代算法,推导出三维目标与粗糙面的耦合积分方程,以及粗糙面散射的Kirchhoff近似(KA)计算式.粗糙面的KA解析计算大大降低了粗糙面求解的复杂度,与目标矩量法的混合迭代保证了计算结果的精度,使得三维体-面目标复合散射计算变得可行.由于体-面两者的高阶耦合作用明显减小,保证了该混合迭代算法的收敛性.与镜像Green函数方法的比较表明该混合算法的有效性,并讨论了粗糙面长度选择对计算结果的影响.结合Monte-Carlo方法,数值分析了理想导体Gauss 关键词： 复合散射 Kirchhoff近似 共轭梯度法 互耦迭代     

三维随机粗糙面上导体目标散射的解析-数值混合算法

提出三维导体目标与导体粗糙面复合散射的解析-数值混合迭代算法，推导出三维目标与粗糙面的耦合积分方程，以及粗糙面散射的Kirchhoff近似(KA)计算式.粗糙面的KA解析计算大大降低了粗糙面求解的复杂度，与目标矩量法的混合迭代保证了计算结果的精度，使得三维体-面目标复合散射计算变得可行.由于体-面两者的高阶耦合作用明显减小，保证了该混合迭代算法的收敛性.与镜像Green函数方法的比较表明该混合算法的有效性，并讨论了粗糙面长度选择对计算结果的影响.结合Monte-Carlo方法，数值分析了理想导体Gauss     

基于等几何分析的边界元法求解Helmholtz问题

将基于一类局部双变量B样条函数的等几何分析方法和Burton-Miller方法相结合,分析三维Helmholtz问题.对于某些从二维参数域映射到三维空间具有奇异点的参数曲面,该方法可以有效地避免奇异点处大量奇异与近奇异积分的计算.数值算例表明该方法具有较好的计算精度和计算效率.复杂问题的分析表明,该方法具有良好的工程应用前景.     

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

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

A fast wave superposition spectral method with complex radius vector combined with two-dimensional fast Fourier transform algorithm for acoustic radiation of axisymmetric bodies

Based on the two-dimensional fast Fourier transform (2D FFT) algorithm, a wave superposition spectral method with complex radius vector has been proposed to efficiently analyze the acoustic radiation from an axisymmetric body. First, the complex Fourier series are used along both circumferential and meridian directions, to expand the integral kernel function and unknown source strength density distributed function. Then, by means of the rectangular integral formula, the radiation sound pressure is described in the form of two-dimensional discrete Fourier transform and generalized through 2D FFT algorithm. Finally, several numerical examples are performed to verify the accuracy and efficiency of the present method. Comparing with the other existing analysis ways, the present method has three different characteristics: (1) there is no singular integral in the numerical computation; (2) the unique solution can be given for all eigen wavenumbers owing to the application of the virtual boundary technology with complex radius vector; and (3) the computational efficiency is improved remarkably because all Fourier terms are calculated simultaneously through 2D FFT algorithm.     

Two-dimensional coupled mode equation for grating waveguide excitation by a focused beam

The problem of grating coupling of a focused incident beam under non-normal incidence into a slab waveguide is given a complete three-dimensional (3D) solution. The diffracted field is expressed as the Fourier integral of a regular part and of a singular part resulting from the existence of the coupled guided mode. A suitable expression of the field in the neighborhood of the pole and a rigorous definition of the modal field lead to a generalized coupled mode equation relating the incident field and the two-dimensional (2D) modal field propagating in the plane of the slab waveguide. The phenomenological parameters involved in the coupled wave equation: the propagation constant, the radiation coefficient as well as the modal field shape are derived from the exact treatment of plane wave diffraction in the same structure. The solution of the complete coupling problem is given in the particular case of a Gaussian incident beam, and of a high index step-index waveguide.     

Surface acoustic wave distribution and acousto-optic interaction in proton exchanged LiNbO3 waveguides

The efficiency of acoustooptic (AO) interaction in YZ-cut proton exchanged (PE) LiNbO₃ waveguides is theoretically analysed by determining the overlap between the optical and acoustic field distributions. The present analysis takes into account the perturbed SAW field distribution due to the presence of the PE layer on the LiNbO₃ substrate determined by the rigorous layered medium approach. The overlap is found to be significant upto very high acoustic frequencies of the order of 5 GHz, whereas in the earlier analysis by vonHelmolt and Schaffer [6] for diffused waveguides, it was shown that the overlap integral rolls down to nearly zero at this high frequency range.     

Multilevel fast multipole algorithm for acoustic wave scattering by truncated ground with trenches

The multilevel fast multipole algorithm (MLFMA) is extended to solve for acoustic wave scattering by very large objects with three-dimensional arbitrary shapes. Although the fast multipole method as the prototype of MLFMA was introduced to acoustics early, it has not been used to study acoustic problems with millions of unknowns. In this work, the MLFMA is applied to analyze the acoustic behavior for very large truncated ground with many trenches in order to investigate the approach for mitigating gun blast noise at proving grounds. The implementation of the MLFMA is based on the Nystrom method to create matrix equations for the acoustic boundary integral equation. As the Nystrom method has a simpler mechanism in the generation of far-interaction terms, which MLFMA acts on, the resulting scheme is more efficient than those based on the method of moments and the boundary element method (BEM). For near-interaction terms, the singular or near-singular integrals are evaluated using a robust technique, which differs from that in BEM. Due to the enhanced efficiency, the MLFMA can rapidly solve acoustic wave scattering problems with more than two million unknowns on workstations without involving parallel algorithms. Numerical examples are used to demonstrate the performance of the MLFMA with report of consumed CPU time and memory usage.     

Numerical quadrature for singular and near-singular integrals of boundary element method and its applications in large-scale acoustic problems

The numerical quadrature methods for dealing with the problems of singular and near-singular integrals caused by Burton-Miller method are proposed,by which the conventional and fast multipole BEMs(boundary element methods) for 3D acoustic problems based on constant elements are improved.To solve the problem of singular integrals,a Hadamard finite-part integral method is presented,which is a simplified combination of the methods proposed by Kirkup and Wolf.The problem of near-singular integrals is overcome by the simple method of polar transformation and the more complex method of PART(Projection and Angular Radial Transformation).The effectiveness of these methods for solving the singular and near-singular problems is validated through comparing with the results computed by the analytical method and/or the commercial software LMS Virtual.Lab.In addition,the influence of the near-singular integral problem on the computational precisions is analyzed by computing the errors relative to the exact solution.The computational complexities of the conventional and fast multipole BEM are analyzed and compared through numerical computations.A large-scale acoustic scattering problem,whose degree of freedoms is about 340,000,is implemented successfully.The results show that,the near singularity is primarily introduced by the hyper-singular kernel,and has great influences on the precision of the solution.The precision of fast multipole BEM is the same as conventional BEM,but the computational complexities are much lower.     

边界元奇异与近奇异数值积分方法及其应用于大规模声学问题

提出了综合处理Burton-Miller方法所导致的奇异积分与近奇异积分问题的数值求积方法,以此改进了基于常量元素的常规边界元和低频快速多极边界元方法。对于奇异积分问题,利用Hadamard有限积分方法进行解决;对于近奇异积分问题,则采用极坐标变换法和PART方法(Projection and Angular&;Radial Transformation)进行克服。与解析解和LMS Virtual.Lab商业软件的结果比较验证了方法的正确性,并对比分析了奇异积分与近奇异积分对计算精度的影响。采用低频快速多极子方法以加速常规边界元法的计算效率,计算分析了计算复杂度,并成功实现了34万自由度大规模问题的计算。结果表明,近奇异积分问题主要由超奇异核函数引起,对计算精度的影响不容忽略;快速多极边界元法的精度与常规边界元法一致,但计算复杂度要远低于后者。      

High-order quadratures for the solution of scattering problems in two dimensions

We construct an iterative algorithm for the solution of forward scattering problems in two dimensions. The scheme is based on the combination of high-order quadrature formulae, fast application of integral operators in Lippmann–Schwinger equations, and the stabilized bi-conjugate gradient method (BI-CGSTAB). While the FFT-based fast application of integral operators and the BI-CGSTAB for the solution of linear systems are fairly standard, a large part of this paper is devoted to constructing a class of high-order quadrature formulae applicable to a wide range of singular functions in two and three dimensions; these are used to obtain rapidly convergent discretizations of Lippmann–Schwinger equations. The performance of the algorithm is illustrated with several numerical examples.