基本解方法与边界节点法数值计算Helmholtz方程的比较研究.pdfVIP

基本解方法与边界节点法数值计算 Helmholtz方程的比较研究 * 姜欣荣 陈文 河海大学工程力学系，工程与科学数值模拟软件中心，江苏 南京 210098 摘要：基本解方法和边界节点法是基于径向基函数的两种重要无网格边界离散 数值技术。针对Helmholtz方程，本文比较研究了这两种数值方法在不同计算区 域上的计算精度、插值矩阵对称性和病态性、计算成本。结果表明数值计算精度 依赖于数值方法中的插值矩阵的秩。另外，基本解方法可以调整配置点的位置使 得插值矩阵满秩。边界节点法不一定能得到满秩矩阵，但可以得到对称的插值矩 阵，同时计算成本相对较低。 关键词：基本解方法，边界节点法，径向基函数，Helmholtz方程 [中图分类号] O302 [文献标识码] A A Comparative Study of Numerical Solution of Helmholtz Equation by Method of Fundamental Solution and Boundary Knot Method Abstract : Both based on the radial basis function (RBF), method of fundamental solution (MFS) and boundary knot method (BKM) are two important meshless boundary discrete numerical technologies. For the Helmholtz equation, this paper made a comparative study of these two numerical methods in different regions by calculation accuracy, symmetry, pathological of interpolation matrix, and computing cost. The results show that numerical accuracy relates with the interpolation matrix rank. Particularly, the MFS can adjust the location of the source points and collocation points to get nonsingular interpolation matrix, meanwhile the BKM can make a symmetric interpolation matrix, and take lower cost. Keywords: method of fundamental solution, boundary knot method, radial basis function, Helmholtz equation *基金项目：国家自然科学基金面上项目，项目批准号 作者简介：姜欣荣（1985- ），硕士研究生，jxrhhu@hhu.edu.cn 陈文（1967- ），博士生导师，chenwen@hhu.edu.cn 1 第 637 页 1 引言 近十几年来，无网格方法数值求解偏微分方程由于克服了传统的基于网格类 的网格依赖性，减少了因网格畸变而引起的计算困难，引起了广泛的关注[1-4] 。 根据近似方案和离散方案的不同，无网格方法可分为多种，各种方法的详细原理 [