二维Helmholtz方程的插值型边界无单元法.PDFVIP

　 ， ３９ ４ Ａｐｐｌｉｅｄ Ｍａｔｈｅｍａｔｉｃｓ ａｎｄ Ｍｅｃｈａｎｉｃｓ　 　 ２０１８年４月１ 日出版 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 　 Ｖｏｌ．３９，Ｎｏ．４，Ａｐｒ．１，２０１８ 文章编号：１０００⁃０８８７（２０１８）０４⁃０４７０⁃１５ ⓒ 应用数学和力学编委会，ＩＳＳＮ １０００⁃０８８７ 二维 Ｈｅｌｍｈｏｌｔｚ方程的插值型边界无单元法∗ 陈林冲，　 李小林 （重庆师范大学 数学科学学院，重庆４０１３３１） 摘要：　 针对二维 Ｈｅｌｍｈｏｌｔｚ 方程的内外边值问题，提出了插值型边界无单元法（ｉｎｔｅｒｐｏｌａｔｉｎｇ ｂｏｕｎｄａｒｙ ｅｌｅｍｅｎｔ⁃ｆｒｅｅ ｍｅｔｈｏｄ）．在间接位势理论的基础上，利用Ｌａｐｌａｃｅ方程基本解的特性，建立了 求解Ｈｅｌｍｈｏｌｔｚ方程Ｎｅｕｍａｎｎ边值内外问题的正则化形式，有效消除了强奇异积分的计算．其次通 过引入全局距离展开成局部距离的幂级数，详细推导了距离函数的导数和法向导数差值的极限表 达式．最后给出了４个插值型边界无单元法的数值算例，表明了该方法可取得较高的可行性和有 效性． 关　 键　 词：　 内外边值问题；　 插值型边界无单元法；　 正则化；　 强奇异积分；　 幂级数 中图分类号：　 Ｏ２４２．２　 　 　 文献标志码：　 Ａ ＤＯＩ：１０．２１６５６／ １０００⁃０８８７．３８０２０２ 引　 　 言 Ｈｅｌｍｈｏｌｔｚ方程在工程技术、电磁场理论、散射理论、力学等较多领域有着广泛的应用，研 ［１］ 究其数值解不仅有广泛的实际意义，也有着重要的理论价值 ．传统的数值计算方法如有限差 ［２］ 分法、有限元法、边界元法 等，必须要建立与插值节点拓扑相关的网格．在处理大变形以及动 态裂纹扩展等问题时，网格可能会发生畸变，此时必须要进行网格重构，增加了计算时间，降低 了精度．２０世纪后期以来，产生了一种不需要与节点相关联网格的新型数值计算方法———无 ［３⁃４］ 网格方法，并且近年来得到了较大的发展 ． ［５］ 最早的无网格方法是 １９７７年Ｌｕｃｙ提出的光滑粒子法 ，主要用于研究无边界的天体问 ［６］ 题和流体问题，但稳定性不是很好．１９８１年，Ｌａｎｃａｓｔｅｒ等建立了移动最小二乘法 ，采用完备 多项式作为基函数，利用误差的加权平方和构造泛函，得到的逼近函数精度高、光滑性好，函数 本身和其导函数都连续．２０世纪９０ 年代，Ｎａｙｒｏｌｅｓ等将移动最小二乘法引入微分方程边值问 ［７］ 题的数值求解中 ．后来Ｍｕｋｈｅｒｊｅｅ等对移动最小二乘法进行了改进，提出了求解势问题和弹 ［８］ 性力学问题的边界点法（ｂｏｕｎｄａｒｙ ｎｏｄｅ ｍｅｔｈｏｄ） ．边界点法是将移动最小二乘法导出的插值 公式与边界积分方程结合起来，通过边界积分方程的边值条件，建立边界上的线性方程组．但 ［９］ 是该方法容易形成病态方程组，不易施加边界条件．Ｗａｎｇ等提出了插值型边界无单元法 ，将 改进的插值型移动最小二乘近似和边界积分方程相结合，这一方法的主要优点是形函数满足 ∗ 收稿日期：　 ２０１７⁃０７⁃２０；修订日期：　 ２０１７⁃１１⁃２７ 基金项目：　 国家自然科学基金（面上项目）；重庆市