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

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

三维声场边界元法几乎奇异积分问题分析

以三维声场问题为例,提出一种准确计算高阶单元几乎奇异积分的半解析算法.首先分析高阶单元几何特征,构造近似几何量,然后应用扣除法,将奇异积分核函数分解为规则核函数与近似几何量表达的奇异核函数.规则核函数积分采用常规Gauss数值积分计算,奇异核函数积分采用半解析算法计算.给出三维声场内问题和外问题经典算例,计算了近边界点的声压,结果证明本文半解析算法的有效性和准确性.     

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

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

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

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

热应力边界元法中几乎超奇异积分的计算

通过分部积分变换将热弹性力学应力边界积分方程中的超奇异积分转化为强奇异积分,然后与另一个强奇异积分求和,得到仅含几乎强奇异的热应力自然边界积分方程.再对其中的几乎强奇异积分施以正则化,消除了热弹性力学边界元法中的几乎奇异积分,可以准确计算出热弹性力学问题中近边界内点的热应力.算例证明了方法的有效性.     

应用有限元和边界元法计算方坯软接触结晶器的电磁场

给出了用有限元和边界元相结合的方法计算方坯软接触结晶器内钢液电磁场分布的全过程,并对4面体单元基础上的Whitney边元素,H-Φ方程及边界积分方程的离散方法作了重点解释.采用有限元和边界元相结合的方法来计算电磁场的分布可以大大减少计算工作量和计算时间.自行开发了三维电磁场计算程序,将数值模拟结果与物理实际进行了比较.     

两种构造Birkhoff表示的新方法

给出两种构造一阶系统Birkhoff表示的新方法,可以从微分方程直接计算得到Birkhoff函数B和Birkhoff函数组R_μ. 举例说明所得结果的应用. 关键词： 分析力学 Birkhoff方程 Birkhoff表示 一阶微分方程     

标量波传播问题的双互易精细积分法

为避免使用计算多种特征频率下的声场响应,采用双互易方法将边界积分方程中时间二次导数项的域积分转化为边界积分.首先,将计算场点配置在边界上并考虑边界条件,可以获得由内部节点上声压量线性表示的边界节点上的物理量;其次,将计算场点配置于域内离散节点上,将所得边界积分方程组中关于边界物理量用内部节点的声压量线性表示,获得关于声压量的二阶常微分方程组;第三,引入声压变化速度作为未知量,将二阶常微分方程组转化为一阶常微分方程组;最后,采用精细积分法精确求解常微分方程组.数值算例验证了双互易精细积分法的正确性和稳定性.     

子波变换在声辐射和声散射数值解中的应用

通过将边界变量用于波展开,获得了求解二维及三维轴对称声辐射和声散射的边界积分方程的子波谱方法,既可求解Ｄｉｒｉｃｈｌｅｔ、Ｎｅｕｍａｎｎ问题,也可求解混合过值问题;它能处理任意边界条件的油对称体．用三维元方法解决了三维轴对称问题边界于波谱方法的奇异积分。给出了二维问题奇异积分的近似积分公式。给出了子波谱方法的系数计算方法,它与传统的边界元系数计算方法相似,易于计算机程序实现,能处理复杂的边界几何形状。该方法的优点是可以获得稀疏的系数矩阵。算例表明：该方法收敛较快,精度高。     

Burton-Miller改进型边界方程的多频计算方法

提出了一种采用Burton-Miller改进型边界积分方程进行多频计算的方法。将Burton-Miller方程中的高奇异积分转化为弱奇异积分形式,获得Burton-Miller改进型边界积分方程;将方程中格林函数进行Taylor级数展开,并把波数从方程中分离出来,从而使随波数变化的计算矩阵表示为波数的矩阵级数形式。数值分析表明,本方法不仅保证了解在全波数范围内的唯一性,并且计算频率点数较多时可以节约大量时间,提高计算效率。     

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

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

声学Burton-Miller方程边界元法GPU并行计算

基于GPU,对声学Burton-Miller积分方程的边界元解法进行并行计算.提出并行计算格式和程序实现方法,以及Burton-Miller方程中各类奇异(包括强奇异、超奇异)积分的GPU计算和局部修正方法.典型算例结果表明,在特征频率处可获得正确的解,具有较高精度,可在普通个人计算机上快速完成自由度超过2×10⁵的声学边界元分析.为计算声学及相关工程领域的中、大规模声场分析问题提供一种快速、高效、简便的数值计算工具.     

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

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

变系数瞬态热传导问题边界元格式的特征正交分解降阶方法

提出了一种基于边界元法求解变系数瞬态热传导问题的特征正交分解(POD)降阶方法,重组并推导出变系数瞬态热传导问题适合降阶的边界元离散积分方程,建立了变系数瞬态热传导问题边界元格式的POD降阶模型,并用常数边界条件下建立的瞬态热传导问题的POD降阶模态,对光滑时变边界条件瞬态热传导问题进行降阶分析.首先,对一个变系数瞬态热传导问题,建立其边界域积分方程,并将域积分转换成边界积分;其次,离散并重组积分方程,获得可用于降阶分析的矩阵形式的时间微分方程组;最后,用POD模态矩阵对该时间微分方程组进行降阶处理,建立降阶模型并对其求解.数值算例验证了本文方法的正确性和有效性.研究表明:1)常数边界条件下建立的低阶POD模态矩阵,能够用来准确预测复杂光滑时变边界条件下的温度场结果;2)低阶模型的建立,解决了边界元法中采用时间差分推进技术求解大型时间微分方程组时求解速度慢、算法稳定性差的问题.     

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-type singular boundary method for acoustic radiation and scattering

This paper proposes the singular boundary method (SBM) in conjunction with Burton and Miller?s formulation for acoustic radiation and scattering. The SBM is a strong-form collocation boundary discretization technique using the singular fundamental solutions, which is mathematically simple, easy-to-program, meshless and introduces the concept of source intensity factors (SIFs) to eliminate the singularities of the fundamental solutions. Therefore, it avoids singular numerical integrals in the boundary element method (BEM) and circumvents the troublesome placement of the fictitious boundary in the method of fundamental solutions (MFS). In the present method, we derive the SIFs of exterior Helmholtz equation by means of the SIFs of exterior Laplace equation owing to the same order of singularities between the Laplace and Helmholtz fundamental solutions. In conjunction with the Burton–Miller formulation, the SBM enhances the quality of the solution, particularly in the vicinity of the corresponding interior eigenfrequencies. Numerical illustrations demonstrate efficiency and accuracy of the present scheme on some benchmark examples under 2D and 3D unbounded domains in comparison with the analytical solutions, the boundary element solutions and Dirichlet-to-Neumann finite element solutions.     

Shape design sensitivity analysis for the radiated noise from the thin-body

Many industrial applications generally use thin-body structures in their design. To calculate the radiated noise from vibrated structure including thin bodies, the conventional boundary element method (BEM) using the Helmholtz integral equation is not an effective resolution. Thus, many researchers have studied to resolve the thin-body problem in various physical fields. No major study in the design sensitivity analysis (DSA) fields for thin-body acoustics, however, has been reported.A continuum-based shape DSA method is presented for the radiated noise from the thin-body. The normal derivative integral equation is employed as an analysis formulation. And, for the acoustic shape design sensitivity formulation, the equation is differentiated directly by using material derivative concept. To solve the normal derivative integral equation, the normal velocities on the surface should be calculated. In the acoustic shape sensitivity formulation, not only the normal velocities on the surface are required but also derivative coefficients of the normal velocities (structural shape design sensitivity) are also required as the input. Hence, the shape design sensitivity of structural velocities on the surface, with respect to the shape change, should be calculated. In this research, the structural shape design sensitivities are also obtained by using a continuum approach. And both a modified interpolation function and the Cauchy principle value are used to regularize the singularities generated from the acoustic shape design sensitivity formulation.A simple annular disk is considered as a numerical example to validate the accuracy and efficiency of the shape design sensitivity equations derived in this research. The commercial BEM code, SYSNOISE, is utilized to confirm the results of the developed in-house code based on a normal derivative integral equation. To validate the calculated design sensitivity results, central finite difference method (FDM) is employed. The error between FDM and the analytical result are less than 3%. This comparison demonstrates that the proposed design sensitivities of the radiated pressure are very accurate.