无穷凹角区域各向异性问题的自然边界元法.pdfVIP

第32卷第5期 应用数学学报 V01．32No．5 2009年9月 ACTAMATHEMATICAEAPPLICATAESINICA Sep．，2009 无穷凹角区域各向异性问题的 自然边界元法木 陈亚军 杜其奎 (南京师范大学数学科学学院，南京210097) (E-mail：duqikui@njnu．edu．cn) 摘要本文研究无穷凹角区域上一类各向异性问题的自然边界元法．利用自然边界归化原理，获得该问 题的Poisson积分公式和自然积分方程，给出了自然积分方程的数值方法，以及逼近解的收敛性和误差估 计，最后给出r数值例子，以示方法的可行性和有效性． 关键词无穷凹角区域；各向异性问题；自然边界元法 MR(2000)主题分类65N30 中图分类0242．21 1 引言 科学与工程计算中常常遇到无界区域问题，数值求解这类问题一直是人们关注的 热点．由于区域的无界性，这给数值计算带来了一定的困难．虽然经典的边界元方法可 用于求解无界区域问题，但在实际处理时，涉及到大量奇异积分的数值计算，计算量较 大．20世纪70年代末由冯康教授和余德浩教授首创并发展起来的自然边界元方法(可 参见【1—4】)，与经典的边界元方法相比具有独特的优点：易实现，数值稳定性好，与有限 元基于同一变分原理，可与有限元自然直接地耦合．用有限元与自然边界元耦合法求解 凹角型外问题时，人们通常选取圆弧作人工边界，但对具有“凹角长条型”内边界的外 问题，以圆弧作人工边界显然并非最佳选择，它将导致一些不必要的计算，甚至不能获 得令人满意的数值结果．可以预测，用一个接近凹角长条型区域边界形状的人工边界 (如椭圆弧)可能会更好．在此基础上也可研究一类各向异性常系数椭圆型方程的外边 值问题，得到圆弧或椭圆弧上精确的人工边界条件，即自然积分方程，并发展相应的耦 合算法和区域分解算法． 本文2008年7月9日收到．2009年3月26日收到修改稿． ’国家自然科学基金资助项目． 应 用 数 学 学 报 32卷 设Q为具有角度u的无穷凹角区域，0u≤2丌，区域的边界为互不重叠的三部 件．问题的描述如下： o，(z，们∈‘2 (1) n耄+6笔-0(删)∈Q n丽+6研2 U u=0，(z，Y)∈FoF。， (2) 【z，ⅣJ∈1， (3) 。n。关+‰瓦Ou：g，(训)∈r， on2瓦+溉Ⅳ瓦2， 其中a，b均为正常数，他为r上关于Q的单位外法线方向．为确保问题(1)一(3)的可解 性，还需恰当的无穷远边界条件． 本文在坐标变换的基础上，利用自然边界归化原理，获得该问题的Poisson积分公 式和自然积分方程，给出了自然积分方程的数值方法，以及逼近解的收敛性和误差估 计，最后给出了数值例子，以示方法的可行性和有效性． 2 自然边界!J|化 2．1凹角扇形外区域上问题(1)-(3)的自然边界归化 设ba0，区域Q，边界r，ro和L描述如下： Q=_((r，口)frR，0伊u}，F=((r，矽)fr=R，0护Ⅳ，， Fo={(r，口)lrR，p=o)，F。=．[(r，口)IrR，口=u)， 其中(n护)为极坐标．它与直角坐标(z，Y)有如下关系： X=rCOSp，