二维Poisson方程的Legendre Tau方法的误差估计
2011-01-31沈婷婷马和平
沈婷婷, 马和平
(上海大学理学院,上海200444)
谱方法作为数值求解偏微分方程的有效工具之一,近年来得到了广泛的应用.根据选取检验函数的不同,谱方法可分为Galerkin方法、tau方法和配置法.关于Galerkin方法和配置法的一些理论分析和数值计算结果可参见文献[1-3].本研究主要讨论tau方法在二维问题中的收敛性态.

文献[12]虽然给出了Legendre tau方法求解一维二阶微分方程的L2模的最优误差估计,但对于高维情况下的收敛结果却没有具体讨论.由于tau方法在高维情况下仅有次优的误差估计,因此,本研究对于高维情况下 tau方法的收敛性态更感兴趣.
本研究考虑二维Poisson方程的Dirichlet问题,即




本研究主要目的是证明Legendre tau方法对于求解二维Poisson方程Dirichlet问题具有H1模和L2模的最优误差估计.文献[12]取类似(1-x2)-1uN作为检验函数,其中uN∈(I),I=(-1,1).本研究将该方法应用到二维的情形.
1 基本引理
对任意的u∈H2(Ω),v∈L2(Ω),定义双线性形式为

引理1 对任意的u∈H2(Ω),v∈H1(Ω),存在常数C,使得

证明 由文献[2]中的结论,得到

另一方面,有

因此,

如果考虑对u,v∈H1(Ω)定义a(u,v)(弱形式),则同样可以得到以上结果.
下面引入2个正交投影算子.

引理2[1]如果 u∈Hr(Ω),且0≤l≤2≤r,则有


引理3[1]如果u∈Hr(Ω)∩H10(Ω),且0≤l≤1≤r,则有

2 收敛性分析
由式(1)和(2),可以得到如下误差方程:


假设ω-1,-1U∈H2,设u*=ω1,1PN-22(ω-1,-1U),e=uN-u*,有∶=ω-1,-1e∈N-2,因此,


并且,

因此,由式(4),得

根据引理2,可得

另一方面,

由三角不等式和Poincaré不等式,可以得到以下定理.

式中,C为依赖于‖ω-1,-1U‖r的正常数.
下面利用对偶技巧来估计‖U-uN‖.

式中,C为依赖于‖ω-1,-1U‖r的正常数.
证明 考虑如下问题:对于g∈L2(Ω),令φ= φ(g),满足

根据文献[1],可知式(5)有唯一解,并且φ满足


因此,

由定理1和……