圆外区域求解Helmholtz方程半.PDFVIP

2006年 12月 高 等 学 校 计 算 数 学 学 报 第28卷第4期 圆外区域求解Helmholtz方程半 关晓明 林伟 (中山大学数学系，广州510275) SOIⅣINGHELMHOIJTZEQU觚 ’ION IN AN EXTERIOR CIRCULAR DOM AIN GuanXiaoming LinWei (Dept，ofMath。，ZhongshanUniversity,Guangzhou510275) Abstract Inthepresentthesiswearemainlyconcernedwiththewaveletmethod fornumericalsolutionoftheboundaryvalueproblem ofHelmholtzequationwhich raewidelyappliedinmechanicsandengineering．Asaresult，thestiffnessmatrixis circulant，symmetryornati-symmetrywherethenonzeroentrieshavefiniteexplicit expression．Sothatthecomputationalcostisextremelydecreasedna dtheaccuracy isextremelyimproved． Key words wavelet，Helmholtzequation，naturalboundray elementmethod， numericalsolution． AMS(2000)subjectclassifications35J25 中圈法分类号 O175．25 1引 言 许多科学和工程计算问题都可以归结为无界区域上的偏微分方程边值问题．而求解 椭圆方程边值问题的常用技术是有限元方法，可是对于无界区域，在用有限元方法求解时， 往往遇到困难．最简单的办法显然是直接略去区域的无界部分求解，但这样做或者导致过 低的计算精度，或者要付出很高的计算代价．边界归化，即将求解偏微分方程边值问题转 化为边界积分方程，是求解某些无界区域问题的强有力的手段．自70年代以来，有限元和 国家自然科学基金资助项目；广东省自然科学基金资助项目(021768) 收稿 日期：2004-06．13． · 290 · 关晓明等：圆外区域求解 Helmholtz方程 第4期 边界元耦合法已经逐渐成为求解无界区域的主要方法 1【,2，3】．但是由于自然边界元方法 的核函数具有强奇性，使得得到的方程组的系数矩阵不再具有稀疏性 1【2】，因而一般不能 实现快速算法．小波在积分方程数值计算中的应用，已有不少很出色的工作．在 8【】中陈 翰林和彭思龙教授应用周期拟小波求解第二类带对数核奇性的Fredholm积分方程，证明 其数值解的收敛性．在 9【】中陈仲英，Miccheli和许跃生教授应用多小波格式给出第二类 Predholm方程的Petrov-Galerkin小波多尺度数值算法，在 1【0，11】中，他们构造了不变集 上的高维小波，提出求解奇异积分方程的小波快速配鬣法和矩阵的有效的压缩策略，从而 实现快速计算．本文用Quak三角小波，得到的刚度矩阵是循环与对称或反对称的，在一 定程度实现了快速算法．80年初开始兴起了区域分解算法，这一算法把计算区域分解为 若干个子域分别求解并通过迭代获得整体近似值。由于它将问题分解，由大化小，由复杂 化简单，优越性显而易见． 基于这一思想，我们将一个平面上的无界区域用一个圆分为两个部分，其中一个部分 是有界区域，而另一部分则是以圆周为边界的无界区域．在有界区域中，可以使用有限元 的方法求解，而在无界区域中我们则可以使用自然边界元的方法来求解，这样，就可以解 决整个问题了． 许多数学物理问题，例如时间调和声波对不可穿透的障碍物的散射，海洋水下声波 的传播，电磁波的绕射与