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

第33卷第4期 应用数学学报 V01．33No．4 2010年7月 ACTAMATHEMATICAEAPPLICATAESINICA July,2010 无穷凹角区域各向异性问题的 自然边界元与有限元耦合法木 陈亚军 (上海海事大学高等技术学院，上海200136) (南京师范大学数学科学学院，南京210046) (E-mail：chenyajun@shmtu．edu．an) 杜其奎 (南京师范大学数学科学学院，南京210046) 摘要本文研究无穷凹角区域上一类各向异性问题的自然边界元与有限元耦合法．利用自然边界归化原 理，获得圆弧或椭圆弧人工边界上的自然积分方程，给出了耦合的变分形式及其数值方法，以及逼近解的收 敛性和误差估计，最后给出了数值例子，以示方法的可行性和有效性． 关键词无穷凹角区域；各向异性问题；自然边界归化；耦合法 MR(2000)主题分类65N30 中图分类0242．21 1 引言 科学与工程计算中常常遇到无界区域问题，数值求解这类问题一直是人们关注的 热点．由于区域的无界性，这给数值计算带来了一定的困难．虽然经典的边界元方法可 用于求解无界区域问题，但在实际处理时，涉及到大量奇异积分的数值计算，计算量 较大． 20世纪70年代末由冯康教授和余德浩教授首创并发展起来的自然边界元方法 [i-4]，与经典的边界元方法相比具有独特的优点：易实现，数值稳定性好，与有限元基 于同一变分原理，可与有限元自然直接地耦合． 由于方法本身的限制，此前采用有限元与自然边界元耦合法求解凹角型内边界外 问题时，人们通常选取圆弧作人工边界【5】’但对具有“凹角长条型”内边界的外问题， 以圆弧作人工边界显然并非最佳选择，它将导致一些不必要的计算，甚至不能获得令人 本文2009年6月23日收到．2010年5月23日收到修改稿 +国家自然科学基金资助项目． 4期 陈亚军，杜其奎：无穷凹角区域各向异性问题的自然边界元与有限元耦合法 759 满意的数值结果．另外，对某些问题，以圆弧作为人工边界较难获得Poisson积分公式 和自然积分方程【6J． 鉴于此，本文在坐标变换及凹角椭圆外区域上调和问题的自然边界归化的基础上， 提出了求解无穷凹角区域上一类各向异性常系数椭圆方程的自然边界元与有限元耦合 法，给出了耦合的变分形式及其数值方法，以及逼近解的收敛性和误差估计，最后给出 了数值例子，以示方法的可行性和有效性． 2 各向异性问题的自然边界归化 设i2为具有角度u的无穷凹角区域，0u≤2r，区域的边界琊Z分为互不重叠的 三部分r，ro和L．考虑如下椭圆型微分方程的Dirichlet边值问题 f。孬02u+6酽02u扎Q内， u (1) {u：0， ro L上， 【u：“。， r上 和Dirichlet—Neumann混合边值问题 f。礤02u¨象_o，Q内， u (2) {u：o， ro L上， L 9， r上，I’上， 【。礼。五O—u+觇可．O磊u：9， o礼z瓦+％瓦2 其中n满足恰当的无穷远边界条件，a，b均为正常数，佗为r上关于Q的单位外法 线方向． 今设b