有限元素法與有限體積法

Céa 引理(Céa's lemma)

/ say-AH /

伽遼金法在能量範數下找到子空間所能提供的「最佳」近似——但若你在意的是一般的精度,「能量意義下的最佳」是個奇怪的尺度。Céa 引理是那座令人安心的橋樑:它證明有限元素解本質上「好到」幾乎與所選形狀函數所能做到的最佳一樣,至多差一個已知的常數倍。它把抽象的伽遼金正交性性質,轉換成一個關於精度的具體承諾。

精確地說,對一個其雙線性形式 a 有界(常數 M)且強制(coercive,常數 alpha)的問題,Céa 引理斷言伽遼金解 u_h 的誤差滿足 ||u - u_h|| <= (M / alpha) * inf over v_h in V_h of ||u - v_h||。仔細讀右邊:那個下確界是 u 在子空間 V_h 中「最佳可能」近似的誤差——任何形狀函數組合所能最接近的程度。所以 Céa 說:你算出的 u_h 落在最佳值的 M/alpha 倍固定因子之內。這稱為「擬最佳性(quasi-optimality)」。實務上的後果是,有限元素的收斂速率化約成一個純粹的「逼近」問題:網格尺寸 h 上的 p 次分段多項式能把 u 逼近得多好?插值理論對光滑的 u 回答以 O(h^p)——而 Céa 引理立刻把這個速率轉移到實際的有限元素誤差上。

Céa 引理是入門有限元素誤差理論的基石,但它的誠實藏在常數 M/alpha 裡。對對稱強制的問題,這個因子不大,方法確實近乎最佳。但對強制性很弱的問題——對流主導的傳輸,其中 alpha 極小——常數 M/alpha 可能巨大,那個界雖仍為真,卻變得幾乎無用:伽遼金解可能離真解很遠(正是促成迎風與穩定化方法的虛假振盪)。所以 Céa 引理保證的是擬最佳性,而非最佳性,而那道落差恰恰是困難又有趣的問題所在之處。

對 -u'' = f,用線性元素逼近一個光滑解,插值理論給出能量範數下的最佳近似誤差 O(h)。Céa 引理於是保證實際的有限元素誤差也是 O(h),並透過對偶(Aubin-Nitsche)論證在 L2 範數下為 O(h^2)。方法的收斂速率完全來自帽函數對 u 的逼近能力。

有限元素誤差 <= (M/alpha) x 最佳可能逼近誤差——收斂化約為逼近理論。

Céa 引理證明的是擬最佳性,而非最佳性:當問題勉強強制(強對流)時,常數 M/alpha 可能爆增,使界雖真卻鬆,且單純的伽遼金解會振盪。那種情形正是穩定化與迎風方法的用武之地。

又稱
Cea's lemmaquasi-optimality estimatebest-approximation property擬最佳性最佳逼近性質