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

Mechanical Quadrature Methods and Extrapolation for Solving Nonlinear Boundary Integral Equations of Helmholtz Equation

作　　者：程攀[1] 黄晋[2] 王柱[2,3] 

机构地区：[1]重庆交通大学理学院,重庆400074 [2]电子科技大学数学科学院,成都611731 [3]维吉尼亚综合技术大学数学系

出　　处：《应用数学和力学》2011年第12期1405-1414,共10页Applied Mathematics and Mechanics

基　　金：国家自然科学基金资助项目(10871034);重庆市自然科学基金资助项目(CSTC20-10BB8270);AFOSR(FA9550-08-1-0136);NSF(OCE-0620464)的部分资助项目

摘　　要：当Helmholtz微分方程转化为非线性边界积分方程后,可以利用机械求积法求得近似解,此方法具有较高的收敛精度阶O(h3)和较低的计算复杂度.构造机械求积法时,一个非线性方程系统通过离散非线性积分方程得到.此外,每个矩阵元素的值都不需要计算任何奇异积分.根据渐近紧理论和Stepleman定理,整个系统的稳定性和收敛性得到了证明.利用h3-Richardson外推算法,收敛精度阶可以提高到O(h5).为了求解非线性方程组,利用Ostrowski不动点定理研究了Newton的解的收敛性.几个算例从数值上说明了本算法的有效性.Mechanical quadrature methods（MQMs） for solving nonlinear boundary integral equations of Helmholtz equation,which possessed high accuracy order O（h3） and low computing complexities,were presented.Moreover,the mechanical quadrature methods were simple without computing any singular integration.A nonlinear system was constructed by discretizing the nonlinear boundary integral equations.The stability and convergence of the system were proved based on asymptotical compact theory and Stepleman theorem.Using the h3-Richardson extrapolation algorithms（EAs）,the accuracy order to O（h5） was improved.For solving the nonlinear system,Newton iteration was discussed extensively by Ostrowski fixed point theorem.The efficiency of the algorithms was illustrated by numerical examples.

关 键 词：HELMHOLTZ方程 机械求积法 NEWTON迭代法 非线性边界条件 

分 类 号：O24[理学—计算数学] O39[理学—数学]