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

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

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

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

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

与模型自由度的平方成正比的存储量和计算量,使传统边界元无法应用到大型模型的计算。为此,发展了一种二维声学多层快速多极子边界元算法。通过二维Helmholtz核函数展开理论的简要介绍,推导了源点矩计算、源点矩转移、源点矩至本地展开转移、本地展开转移公式,并详细描述了二维声学快速多极子边界元算法的具体实现步骤。使用快速傅里叶插值进行源点矩和本地展开系数的多层传递。采用对角左预处理方法,改善边界方程的条件数,减少迭代求解次数。最后通过数值算例,验证了所发展的二维声学快速多极子算法的正确性和高效性。      

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

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

新型快速多极杂交边界点法在复合材料热传导中的应用

将杂交边界点法应用于复合材料的热传导模拟,推导一种求解复合材料的方程,该方程减少计算自由度,效率更高.将新型快速多极算法与杂交边界点法结合进行大规模计算,数值算例中对包含大量粒子的复合材料进行模拟,结果表明快速多极杂交边界点法可行,具有一定的应用前景.     

认知不确定二维声场的响应特性数值分析

针对声学参数存在认知不确定性的问题,为实现认知不确定声场声压响应的预测。提出了解决二维认知不确定声场的有限元法(Evidence Theory-based Finite Element Method,ETFEM),引入证据理论,采用焦元和基本可信度的概念来描述认知不确定参数,基于摄动法的区间分析技术,推导了认知不确定声场声压响应的标准差、期望求解公式。为验证本文方法的可行性。以认知不确定参数下的二维管道声场模型和某轿车二维声腔模型为例进行了数值计算,对比离散随机变量得到认知不确定参数的声场分析结果和相应的随机声场所得分析结果,研究表明:该方法能够有效的处理认知不确定参数下的二维声场,为工程问题中噪声预测提供可靠的分析模型。      

模糊随机参数二维声场数值计算

解决声场参数同时具有模糊性和随机性的问题,实现模糊随机声场声压响应的预测,引入了信息熵理论,利用信息熵的等效转换,将模糊随机声场转化为纯随机声场或者纯模糊声场进行求解,推导了基于摄动法的二维随机声场和模糊声场的有限元计算公式。以模糊随机参数下的二维管道声场模型和某轿车二维声腔模型为例进行了数值计算,所得结果与蒙特卡洛法(Monte Carlo Method)所预测声压变化范围基本一致,同时,转化为纯随机声场和纯模糊声场所求得声压响应变化范围也基本一致,说明了本文方法计算结果的准确性。因此本文方法能很好地应用于模糊随机参数下二维声场的预测,具有重要的工程应用价值。      

模糊随机参数二维声场数值计算

解决声场参数同时具有模糊性和随机性的问题,实现模糊随机声场声压响应的预测,引入了信息熵理论,利用信息熵的等效转换,将模糊随机声场转化为纯随机声场或者纯模糊声场进行求解,推导了基于摄动法的二维随机声场和模糊声场的有限元计算公式。以模糊随机参数下的二维管道声场模型和某轿车二维声腔模型为例进行了数值计算,所得结果与蒙特卡洛法(Monte Carlo Method)所预测声压变化范围基本一致,同时,转化为纯随机声场和纯模糊声场所求得声压响应变化范围也基本一致,说明了本文方法计算结果的准确性。因此本文方法能很好地应用于模糊随机参数下二维声场的预测,具有重要的工程应用价值。     

新型快速多极边界元法求解电荷任意分布的二维静电场

在初始快速多极边界元法(FMM)基础上提出一种适合位势问题的新型快速多极边界元格式,并用于求解静电场问题.新型算法引入对角化概念,减少了形成局部展开系数的时间,提高计算效率.最后给出数值算例,证明了新型算法的计算精度及处理大规模问题的速度优势.     

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

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

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.     

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...     

The Method of Fundamental Solutions for Solving Exterior Axisymmetric Helmholtz Problems with High Wave-Number

In this paper, we investigate the method of fundamental solutions (MFS) for solving exterior Helmholtz problems with high wave-number in axisymmetric domains. Since the coefficient matrix in the linear system resulting from the MFS approximation has a block circulant structure, it can be solved by the matrix decomposition algorithm and fast Fourier transform for the fast computation of large-scale problems and meanwhile saving computer memory space. Several numerical examples are provided to demonstrate its applicability and efficacy in two and three dimensional domains.     

Solution of Two-Dimensional Stokes Flow Problems Using Improved Singular Boundary Method

In this paper, an improved singular boundary method (SBM), viewed as one kind of modified method of fundamental solution (MFS), is firstly applied for the numerical analysis of two-dimensional (2D) Stokes flow problems. The key issue of the SBM is the determination of the origin intensity factor used to remove the singularity of the fundamental solution and its derivatives. The new contribution of this study is that the origin intensity factors for the velocity, traction and pressure are derived, and based on that, the SBM formulations for 2D Stokes flow problems are presented. Several examples are provided to verify the correctness and robustness of the presented method. The numerical results clearly demonstrate the potentials of the present SBM for solving 2D Stokes flow problems.     

The method of fundamental solutions for 2D and 3D Stokes problems

A numerical scheme based on the method of fundamental solutions (MFS) is proposed for the solution of 2D and 3D Stokes equations. The fundamental solutions of the Stokes equations, Stokeslets, are adopted as the sources to obtain flow field solutions. The present method is validated through other numerical schemes for lid-driven flows in a square cavity and a cubic cavity. Test results obtained for a rectangular cavity with wave-shaped bottom indicate that the MFS is computationally efficient than the finite element method (FEM) in dealing with irregular shaped domain. The paper also discusses the effects of number of source points and their locations on the numerical accuracy.     

A multilayer method of fundamental solutions for Stokes flow problems

The method of fundamental solutions (MFS) is a meshless method for the solution of boundary value problems and has recently been proposed as a simple and efficient method for the solution of Stokes flow problems. The MFS approximates the solution by an expansion of fundamental solutions whose singularities are located outside the flow domain. Typically, the source points (i.e. the singularities of the fundamental solutions) are confined to a smooth source layer embracing the flow domain. This monolayer implementation of the MFS (monolayer MFS) depends strongly on the location of the user-defined source points: On the one hand, increasing the distance of the source points from the boundary tends to increase the convergence rate. On the other hand, this may limit the achievable accuracy. This often results in an unfavorable compromise between the convergence rate and the achievable accuracy of the MFS. The idea behind the present work is that a multilayer implementation of the MFS (multilayer MFS) can improve the robustness of the MFS by efficiently resolving different scales of the solution by source layers at different distances from the boundary. We propose a block greedy-QR algorithm (BGQRa) which exploits this property in a multilevel fashion. The proposed multilayer MFS is much more robust than the monolayer MFS and can compute Stokes flows on general two- and three-dimensional domains. It converges rapidly and yields high levels of accuracy by combining the properties of distant and close source points. The block algorithm alleviates the overhead of multiple source layers and allows the multilayer MFS to outperform the monolayer MFS.     

Singular meshless method using double layer potentials for exterior acoustics

Time-harmonic exterior acoustic problems are solved by using a singular meshless method in this paper. It is well known that the source points cannot be located on the real boundary, when the method of fundamental solutions (MFS) is used due to the singularity of the adopted kernel functions. Hence, if the source points are right on the boundary the diagonal terms of the influence matrices cannot be derived. Herein we present an approach to obtain the diagonal terms of the influence matrices of the MFS for the numerical treatment of exterior acoustics. By using the regularization technique to regularize the singularity and hypersingularity of the proposed kernel functions, the source points can be located on the real boundary and therefore the diagonal terms of influence matrices are determined. We also maintain the prominent features of the MFS, that it is free from mesh, singularity, and numerical integration. The normal derivative of the fundamental solution of the Helmholtz equation is composed of a two-point function, which is one of the radial basis functions. The solution of the problem is expressed in terms of a double-layer potential representation on the physical boundary based on the potential theory. The solutions of three selected examples are used to compare with the results of the exact solution, conventional MFS, boundary element method, and Dirichlet-to-Neumann finite element method. Good numerical performance is demonstrated by close agreement with other solutions.