自适应边界元法的后验误差估计.pdf

文幸约号;1006-7329(2000)06一0016-04 自适应边界元法的后验误差估计 ‘ 金朝高 (重庆建筑大学 应用科学与技术系，重庆 400045) 奋 摘要:蓄充在Sobolev空问的框架下.对一般的葬于方祖的Galerkin遥近给出了后验误差 J 尹 . ‘ ， 枯计的结果。然后。对以有限平面为屏蔽物的声敬封问题 (共数学摸型是三维Helmholtz r 声 方桩以有限平面为边界的Neumann问理)在三角剑分下给出了其 自适应边界元解法的后 脸诀里佑计的其体表达式。 关 . 调;Sobolev空问;自迷应边界无法，后脸误里枯计 中日分组号:0e4- 文献标识码:A 一飞 在Galerkin愈义下计算边界带有奇点的边界积分方程的数值解时，由于临近奇点处未知量的 变化是急剧的，需要对临近奇点的边界单元进行细分。自适应边界元法就是处理此类问题的一种有 效方法。这种方法要求对精确解与用边界元法求得的近似解之间的误差的大小和分布情况都有准 确的掌握并可以计算。通常，从数值解得到关于局部误差的误差指标，以及用某种范数表示的整体 误差的估计I，而且这个估计t仅由箱入数据和数值解就可计算并能用局部误差指标表示，这就是 所谓后脸误差估计。利用后验误差估计，若发现临近奇点的单元上的误差指标不能满足设定的精度 要求，在计算时自动细分这些单元，如此反复，直至局部误差指标及整体误差估计量同时满足设定 的精度要求。可以说，实现后验误差估计是使用自适应边界元法进行计算的先决条件。 近年来，C.Carstensen,E.P.Stephan,M.Feistaur,G.Hsiao等人就这一课题发表了一些研究成 果.本文从不同的视角进行研究。首先，在Sobolev空间的框架下，对一般算子方程的Galerkin通近 给出了后验误差估计的结果。其次，对以有限平面为屏蔽物的声散射问题的数学模型— 三维 Helmltolt:方程以有限平面为边界的Neuman。间压，进行了较为详尽的讨论。我们用双层位势将该 间题划归为第一类Fredholm边界积分方程并用Galerkin通近求数值解。这里既有方程本身的高阶 奇异性，有限平面的各边也都具有奇异性.我们给出了可供数值计算的变分方程，并在三角剖分下， 给出了后脸误差估计的具体表达式。 1 一般葬子方程Galerkin退近的后验误差估计 设X,Y是两个Sobolev空间，A是自X至Y的一个同构映射(即双连续的线性双射)。于是存 在正常教C使得 RUG(CIIAully,V。 〔X ) 1 ( 且算子方程 Au二 了 ) 2 ( 对任愈fEY有唯一解 “EX, 设Y’是Y的对偶空间，《·，·)是YXY上的对偶配对。对应于算子方程((2)的变分方程是 v,Au=v,f， VvEY` (3) .收枯B期,2000-05-01 作者曰介.金朗衡(1943-),男 四川三台人.教授，主要从事边界元及金触工怪学的研究