偏微分方程的數值方法

Cea 引理(Cea's lemma)

/ say-AH /

你建立了一個有限元素近似,一個自然的憂慮浮現:我算出的答案離真解有多遠?Cea 引理給出一個極為令人安心的答案。它以樸素的話說:伽遼金有限元素解本質上是所選網格與多項式空間「所能產生的最佳」近似——在一個固定常數的倍數內,你無法做得更好,而方法會自動找到接近那個最佳的東西。

以下是精確陳述。設 u 為真正的弱解、u_h 為來自子空間 V_h 的有限元素解。Cea 引理說,能量範數下的誤差滿足 norm(u - u_h) <= (C/alpha) 乘以 V_h 中所有 v_h 上 norm(u - v_h) 的最小值,其中 C 是雙線性形式的有界常數、alpha 是其強制常數。右端是最佳近似誤差——真解到該空間所含最近函數的距離。於是有限元素誤差至多是一個常數乘以不可避免的最佳可能誤差。證明的關鍵是伽遼金正交性:誤差 u - u_h 在雙線性形式下與整個子空間 V_h 正交,這是「有限元素解是真解的投影」這一事實的離散回聲。把這個擬最佳性與標準的內插估計(它說 p 次分段多項式對光滑函數近似得多好)結合,你便得到具體的收斂速率:能量範數下的誤差正比於 h^p,隨網格尺寸 h 縮小。

它為何重要:Cea 引理把拉克斯-米爾格拉姆抽象的存在保證,轉化為可用、定量的誤差界——它是有限元素法值得信任的理論理由,也是一切先驗誤差分析的基礎。兩個誠實的提醒。其一,對病態(近乎非強制)的問題,常數 C/alpha 可能很大,因此擬最佳並不意味緊。其二,它需要強制性與拉克斯-米爾格拉姆所用的同一有界性;對非強制問題(鞍點系統、混合方法),你得改用更強的 inf-sup/Babuska-Brezzi 理論。

用分段線性元素(p = 1)解 -u_xx = f,Cea 引理加內插給出正比於 h 的能量範數誤差:網格減半,誤差大致減半。改用分段二次(p = 2),同一引理預測誤差正比於 h^2——如今網格減半把誤差砍成四分之一。

有限元素誤差在一個常數倍內,就是該多項式空間所能達到的最佳。

擬最佳性假設真解夠光滑、可被良好近似;在內凹角或奇異點附近,最佳近似誤差本身就很大,於是收斂停滯——這正是適應性網格加密所要拯救的情況。

又称
Céa's lemmaquasi-optimality estimatebest-approximation estimate擬最佳性估計