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

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

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

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

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

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

具有弱界面的柱状复合结构中轴对称声导波

利用柱状分层结构中轴对称声导波的波动方程和界面“弹簧”模型，导出了具有弱界面的双层柱状复合结构中轴对称声导波的广义色散方程，计算了刚性联接和滑移联接界面时双层复合管(腔壁为自由界面或腔内充水)和双层复合棒中轴对称声导波的色散特性，分析了弱界面的轴向劲度系数对低阶轴对称声导波的影响- 关键词：     

Navier-Stokes方程的有限元边界元耦合方法

主要研究了三维外部区域上具有Dirichlet边界条件的非定常Navier-Stokes方程的有限元边界元耦合方法,并分析了这一数值解的收敛速度.     

（2+1）维NLBQ方程和KP方程的多dromion结构和相互作用的研究

从（2+1）维双线性形式的非局域Bussinesq(NLBQ)方程和KP方程的隐线孤子解出发,可以找到与某种势所相应的各方向都指数衰减的dromion解.利用图形分析的方法,对这些dromion之间的相互作用进行了详细的研究.发现这两种模型中的dromion间的相互作用只引起位相漂移,不引起形状和速度的变化,也不引起旋转.即dromion间的相互作用是弹性的,没有能量、动量和角动量的交换. 关键词：     

(2+1)维 Zakharov-Kuznetsov 方程的精确解和孤子结构

在符号计算软件Maple的帮助下,利用改进的Riccati方程映射法得到了(2+1)维Zakharov-Kuznetsov方程(ZK)的新显式精确解. 根据得到的解,研究了ZK方程的特殊孤子结构. 关键词： 改进的Riccati方程映射法 Zakharov-Kuznetsov方程 精确解 孤子结构     

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

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

任意边界板结构的声辐射模态误差传感器的设计

获取声辐射模态伴随系数是基于声辐射模态理论进行主动结构声控制(ASAC)的重要环节。以往PVDF分布式传感器的设计难点是振速展开受边界条件的限制,其设计过程往往是针对特定边界条件展开的。本文在声辐射模态理论和两维分布式传感器的压电方程的基础上,将板表面振速分布用Legendre多项式展开,给出了两维板结构的PVDF传感器形状与边界条件无关的设计方法。这样设计得到的传感器能应用于任意边界条件和任意振速分布的两维板结构,且实时性好,拓宽了其应用范围。本文还分别以固定边界条件板及在该板中任取一小区域两种情况为例,证明了该设计方法的可行性。     

Fourier拟谱配点方法对第一类边界积分方程的应用

本文用Fourier拟谱配点方法求解有广泛应用的以对数核为主部的第一类边界积分方程,文中通过对积分算子的象征作拟谱插值来建立近似方程,利用快速Fourier变换将计算切换到频率空间进行。本文计算结果表明,用上述拟谱配点方法计算的数值精度较Galerkin配点法更为满意。     

Exact boundary condition for time-dependent wave equation based on boundary integral

An exact non-reflecting boundary conditions based on a boundary integral equation or a modified Kirchhoff-type formula is derived for exterior three-dimensional wave equations. The Kirchhoff-type non-reflecting boundary condition is originally proposed by L. Ting and M.J. Miksis [J. Acoust. Soc. Am. 80 (1986) 1825] and numerically tested by D. Givoli and D. Cohen [J. Comput. Phys. 117 (1995) 102] for a spherically symmetric problem. The computational advantage of Ting–Miksis boundary condition is that its temporal non-locality is limited to a fixed amount of past information. However, a long-time instability is exhibited in testing numerical solutions by using a standard non-dissipative finite-difference scheme. The main purpose of this work is to present a new exact boundary condition and to eliminate the long-time instability. The proposed exact boundary condition can be considered as a limit case of Ting–Miksis boundary condition when the two artificial boundaries used in their method approach each other. Our boundary condition is actually a boundary integral equation on a single artificial boundary for wave equations, which is to be solved in conjunction with the interior wave equation. The new boundary condition needs only one artificial boundary, which can be of any shape, i.e., sphere, cubic surface, etc. It keeps all merits of the original Kirchhoff boundary condition such as restricting the temporal non-locality, free of numerical evaluation of any special functions and so on. Numerical approximation to the artificial boundary condition on cubic surface is derived and three-dimensional numerical tests are carried out on the cubic computational domain.     

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

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

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

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

Vibration analysis of circular cylindrical double-shell structures under general coupling and end boundary conditions

In this paper, the free and forced vibration analysis of circular cylindrical double-shell structures under arbitrary boundary conditions is presented. This is achieved by employing the improved Fourier series method based on Hamilton’s principle. In the formulation, each displacement component of the cylindrical shells and annular plates is invariantly expanded as the superposition of a standard Fourier series with several supplementary functions introduced to remove the potential discontinuities of the original displacement and its derives at the boundaries. With the introduction of four sets of boundary springs at the coupling interfaces and end boundaries of the shell–plate combination, both elastic and rigid coupling and end boundary conditions can be easily obtained by assigning the stiffnesses of the artificial springs to certain values. The natural frequencies and mode shapes of the structures as well as frequency responses under forced vibration are obtained with the Rayleigh–Ritz procedure. The convergence of the method is validated by comparing the present results with those obtained by the finite element method. Several numerical results including natural frequencies and mode shapes are presented to demonstrate the excellent accuracy and reliability of the current method. Finally, a number of parameter studies concerning various end and coupling boundary conditions, different dimensions of shells and annular plates are also performed.     

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

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

Kinetic theory based lattice Boltzmann equation with viscous dissipation and pressure work for axisymmetric thermal flows

A lattice Boltzmann equation (LBE) for axisymmetric thermal flows is proposed. The model is derived from the kinetic theory which exhibits several features that distinguish it from other previous LBE models. First, the present thermal LBE model is derived from the continuous Boltzmann equation, which has a solid foundation and clear physical significance; Second, the model can recover the energy equation with the viscous dissipation term and work of pressure which are usually ignored by traditional methods and the existing thermal LBE models; Finally, unlike the existing thermal LBE models, no velocity and temperature gradients appear in the force terms which are easy to realize in the present model. The model is validated by thermal flow in a pipe, thermal buoyancy-driven flow, and swirling flow in vertical cylinder by rotating the top and bottom walls. It is found that the numerical results agreed excellently with analytical solution or other numerical results.     

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

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.