机械求积法及其外推算法求解Helmholtz非线性边界积分方程.pdfVIP

andMechanics Mathematics 应用数学和力学，第32卷第12期 Applied 2011年12月15日出版 V01．32，No．12，Dec．15，2011 1000—0887 文章编号：1000—0887(2011)12—1405—10 ④应用数学和力学编委会，1SSN 机械求积法及其外推算法求解Helmholtz 非线性边界积分方程+ 程 攀1，黄 晋2，王柱2’3 (1．重庆交通大学理学院，重庆400074； 2．电子科技大学数学科学院，成都611731； 3．维吉尼亚综合技术大学数学系，布莱克斯堡VA24061，美国) 摘要： 当Helmholtz微分方程转化为非线性边界积分方程后，可以利用机械求积法求得近似解， 此方法具有较高的收敛精度阶D(h3)和较低的计算复杂度．构造机械求积法时，一个非线性方程 系统通过离散非线性积分方程得到．此外，每个矩阵元素的值都不需要计算任何奇异积分．根据渐 法，收敛精度阶可以提高到0(h5)．为了求解非线性方程组，利用Ostrowski不动点定理研究了 Newton的解的收敛性．几个算例从数值上说明了本算法的有效性． 关键词：Helmholtz方程；机械求积法； Newton迭代法；非线性边界条件 中图分类号：024；039文献标志码：A DOI：10．3879／j．issn．1000—0887．2011．12．002 引 言 调和波在绕过障碍物时的传播和发散本质上可以转换为一个两维的问题，通常可以通过 一个Helmhohz方程以及它对应的边界条件来描述这个物理过程．本文主要讨论具有下述非线 性边界条件的Helmhohz方程： J△反～Ot2藐=0， 在力内， (1) 在，上， LOFt／On=一G(五)+．厂(工)=一g(x，冠)+八工)， 其中，ot是一个描述传播系数的正实数，门cR2是一个具有光滑边界，的有界、单连通区域， a／(On)是边界厂上的单位外法向导数，以及g(x，面)，“X)是已知函数． kinson和Chandler[7 3利用求积方法得到具有非线性边界条件的Laplace方程的高精度解． ·收稿日期： 2010．11一i8；修订日期：2011—11—02 基金项目： 作者简介： 程攀(1976一)，男，重庆人，讲师，博士(联系人．E—mail：cheng．pass@sina．corn)． 1405 1406 程 攀 黄 晋 王 柱 D(h3)的收敛速度． 由转化为边界积分方程的Green公式旧引，Helmholtz方程可等价于 J，∈n· (2) Jr Jr nlI 面(y)=Jr|j}+(J，，工)藐(z)d气一Jr^+(y，x)嚣d％ 令Y趋于边界，上的点，并代人方程(1)的边界条件，得到边界，上的微分方程 h+(Y，X)G(巍)ds，= ÷五(y)一fk+(Y，x)u(x)ds，一l 厶 Ji? Jfl