白春陽,石東偉

（河南科技學院,河南新鄉(xiāng)453003）

Sobolev方程在流體穿過裂縫巖石的滲透理論、土壤中濕氣的遷移、不同介質的熱傳導等許多物理問題中有廣泛的應用[1-6].其中文獻[1]建立了特征混合元方法,并分別討論了收斂性和誤差估計.文獻[2]僅對一維Sobolev方程提出了位移有限元法.文獻[4]提出了相應的最小二乘混合元方法,該方法雖然不要求有限元空間滿足LBB穩(wěn)定性條件,但計算量大.文獻[5]提出H1-Galerkin混合有限元方法,雖不要求有限元空間滿足LBB條件,但對通量函數(shù)的光滑性要求更嚴格.文獻[6]研究了幾組簡單的低階協(xié)調元,給出了半離散格式和全離散格式的誤差估計.但是以上研究都是針對協(xié)調元進行的,而且要求剖分滿足正則條件,即,其中是單元K的最大直徑單元K的最大內切圓直徑.然而當問題定義在窄邊區(qū)域時,如果采用傳統(tǒng)的正則剖分計算量將會非常大,這時如果采用各向異性網格剖分可以用較少的自由度而得到同樣的估計結果.在這種情況下,傳統(tǒng)的引理在插值誤差分析中已不再適用,而且對于非協(xié)調元來說傳統(tǒng)邊界估計技巧也不再適用,因為此時對單元K的長邊F誤差值估計將出現(xiàn)一個因子,它可能趨于無窮,無法達到收斂結果.本文給出了與以上文獻不同的Sobloev方程的非協(xié)調混合有限元格式,在各向異性網格下進行了收斂性分析,并利用單元自身的一些特殊性質得到了超逼近結果,同時利用插值后處理技巧導出了超收斂結果.

1 單元構造

使得

由文獻[7]不難驗證上述插值具有各向異性.

2 誤差估計

相應的有限元空間為：

定義插值算子

并且逼近問題（4）的解唯一.

對非協(xié)調元來說由式（3）和式（4）可得誤差方程

3 超逼近分析

所以有

整理得

記

由Vh的構造知

則上式可化為

所以有

兩邊積分得

由Gronwall引理得

所以有

移項得

由Cauchy和ε不等式類似（8）有

把式（7）代入式（8）得

由式（7）和式（9）得

故

所以

綜合上述,估計定理結論成立.

4 結論

本文給出了一個非協(xié)調混合有限元格式,在各向異性網格下進行了收斂性分析,對Sobolev方程進行了Galerkin逼近,并且利用單元的特殊性質在不需要Ritz投影情況下得到了超逼近性,最后利用插值后處理技巧給出一類新的混合元格式的收斂性分析和超逼近結果.對求解在流體穿過裂縫巖石的滲透理論,土壤中濕氣的遷移,不同介質的熱傳導等許多物理問題中有廣泛應用的Sobolev方程提供了一定思路.

[1]趙培,陳煥貞.Sobolev方程的特征混合有限元法[J].應用數(shù)學,2003,16（4）：50-59.

[2]Nakao M T.Error estimates of a Galerkin method for some nonlinear Sobolev equations in one space dimension[J].Numerische Mathematik,1985,47（1）：139-157.

[3]Jiang Z W,Chen H Z.Error estimates for mixed finite element methods for Sobolev equation[J].Northest.Math.,2001,17（3）：301-314.

[4]Gu H M,YangDP.Least-squaresmixed finiteelementmethodsforSobolevequation[J].IndiaJ.Appl.Math.,2003,31（5）：505-517.

[5]郭玲,陳煥珍.Sobolev方程的H1-Galerkin混合有限元方法[J].系統(tǒng)科學與數(shù)學,2006（3）：301-314.

[6]鄭克龍.Sobolev方程的混合有限元方法[D].成都：四川大學,2005：1-39.

[7]Chen S C,Hai D Y,Zhao Y C.Interpolation and Quasi-Wilson element for narrow Quadrila-teral meshes[J].IMA.J.Numer.Anal.,2004,24：77-95.