JOVANA
Explore Library Glossary Getting Started Three Levels Fields How it works Mission
Join the mission
All guides

橢圓算子與能量法

平衡型方程——拉普拉斯算子及它更一般的表親們——全都共享一個良好的形狀:它們是橢圓型的。本篇先把這究竟意味著什麼釘牢,再揭開貫穿這整階的訣竅:別再苛求一個古典解,改成乘上一個測試函數、做分部積分,把偏微分方程化成一個關於能量的乾淨問題。光是這一步,就幾乎免費地換來了存在性。

「橢圓」究竟在命名什麼

在這道階梯很前面的地方,你曾用一個判別式的正負號來分類二階方程:拉普拉斯方程 歸為 橢圓型,熱方程為拋物型,波動方程為雙曲型。橢圓型就是那個 沒有實特徵線 的型——沒有資訊沿之傳播的特殊方向——而這正是 平衡 的標記:一個穩態的溫度、一張安定的薄膜、一個靜態的電位,每一點都同時靜靜地與它所有的鄰居協商。這一階拿起這一個型,建造它成熟的理論。第一件事,是說清楚一旦走出光禿禿的拉普拉斯算子之後,什麼才算是橢圓型。

我們在意的一般二階算子,在二階導數前面收攏了一個係數矩陣 a_ij:L u = - sum over i,j of a_ij u_(x_i x_j),外加可能的低階項。回想一下,只有 主部——最高階的那一塊——決定型別,正如判別式只看領導係數一樣。對兩個變數,橢圓性是 B^2 - A C < 0;而在任何維度,乾淨的說法是係數矩陣 (a_ij) 為 正定:它沒有平的、也沒有反向的方向。那個矩陣為正定,正是 橢圓性 的意思。

一致橢圓性與散度形式

要讓理論真正咬合,我們要求橢圓性帶有一個永不縮減到零的餘裕。一致橢圓性 要求兩個固定的數,0 < theta <= Theta,使得對每個方向 xi、每個點 x,都有 theta |xi|^2 <= sum of a_ij(x) xi_i xi_j <= Theta |xi|^2。換句話說:係數矩陣被夾在單位矩陣的兩個正倍數之間,處處同時如此。下界 theta 是真正幹活的那一個——它保證這個算子總是至少以一個固定的量往「外」推,從不鬆懈。把 theta 記牢;它會作為這一階每一個存在性證明的引擎再度現身。

我們還偏好用一種 形狀 來寫這個算子,而它的重要性超乎其外表。散度形式 把 L u 寫成 L u = - div(A grad u)——(矩陣 A 乘以 u 的梯度)的散度——而不是把二階導數一一攤開。為什麼偏好它?因為散度正是散度定理(分部積分)所為之而生、最善於處理的東西。當你一想把一個導數從 u 身上挪到一個測試函數上時,散度形式 就遞給你一個完美的、可供分部積分的導數,邊界項也自行排列整齊。非散度(「非守恆」)形式也存在,自有其理論,但我們即將遇見的能量法在散度形式裡活得最自然。

general elliptic operator:   L u = - sum_{i,j} a_ij(x) u_{x_i x_j} + (lower-order terms)

divergence form:             L u = - div( A grad u ),     A = (a_ij(x))  symmetric

uniform ellipticity:         theta |xi|^2  <=  sum_{i,j} a_ij(x) xi_i xi_j  <=  Theta |xi|^2
                             for all directions xi and all points x,   0 < theta <= Theta

plainest example:            A = identity  =>  L u = - div(grad u) = - Laplacian u   (the Poisson operator)
一般橢圓算子、它的散度形式,以及一致橢圓性的夾擠。把 A 取為單位矩陣,就回到負拉普拉斯算子。

改變一切的一步:把方程弱化

苛求一個 古典解——一個帶有貨真價實、連續的二階導數、逐點滿足 L u = f 的 u——的麻煩就在這裡。對一個粗糙的區域、一個不連續的源 f,或一個跳動的係數矩陣 A,這樣平滑的 u 也許根本不存在,然而物理上分明有個答案(一張真實的薄膜確實會安定成 某個 形狀)。補救之道不是強迫平滑,而是 要求得更少。把方程 L u = f 乘上一個在邊界上為零的平滑 測試函數 v,在區域上積分,再做一次分部積分,把一個導數從 u 身上剝下、安放到 v 上。

因為我們把 L 寫成了散度形式,那一次分部積分很乾淨:- div(A grad u) 化成 (A grad u) 與 (grad v) 內積的積分,而邊界項因 v 在邊緣上為零而消失。剩下的就是 弱形式:找一個 u,使得 (A grad u) . (grad v) 的積分等於 f v 的積分,對每一個 測試函數 v 都成立。注意這筆帳上的奇蹟——u 如今只需要 一個 導數,而非兩個,因為那第二個導數已被交付給了 v。要求 u 較少平滑性的代價,由 v 來償付,而 v 我們愛取多平滑就取多平滑。

一個對所有測試函數都滿足這個積分恆等式的 u,稱為 弱解,而那個恆等式本身就是 弱形式。u 的自然居所不再是二次可微函數的空間,而是你上一階建起的 索伯列夫空間——只帶有 一個 L^2 中的 弱導數 的函數,並且(為了把零邊界值編碼進去)取其子空間 空間 H^1 sub zero。一切都已安排妥當,使弱問題被擺在一個希爾伯特空間裡。這絕非偶然;它正是為了那句妙語而布的局。

能量:雙線性形式,以及為何強制性致勝

把弱形式的左邊看成一台機器,它吃進 兩個 函數 u 與 v,吐回一個數:B(u, v) = (A grad u) . (grad v) 的積分。它對每一個槽各自都是線性的,所以它是一個 雙線性形式。整個弱問題如今讀來簡潔得驚人:找 H^1 sub zero 中的 u,使得對所有 v 都有 B(u, v) = (f, v),其中 (f, v) 不過是 f v 的積分。我們已把一個偏微分方程熬煮成希爾伯特空間上的一個雙線性恆等式——而這個問題抽象到足以一般地作答。

現在一致橢圓性付出它的紅利。把 同一個 函數餵進兩個槽:B(u, u) = (A grad u) . (grad u) 的積分。下界 theta |xi|^2 說這個被積函數至少是 theta |grad u|^2,所以 B(u, u) >= theta 乘以 |grad u|^2 的積分。這個雙線性形式餵進它自己時,被一個貨真價實的 能量 從下方界住——而這個性質有個名字,強制性:B(u, u) >= alpha ||u||^2,其中 alpha > 0 是個固定常數。(把「|grad u|^2 的積分」變成對 u 完整範數的界,用到了 龐加萊不等式,那是在 H^1 sub zero 中工作、邊界被釘住所賜的禮物。)強制性 恰恰就是「能量永不鬆懈」這句話——是 theta 永不趨零的抽象迴響。

拉克斯-米爾格拉姆:幾乎免費的存在性

我們在希爾伯特空間 H^1 sub zero 上,為我們的雙線性形式 B 集齊了恰好三件事實:它是雙線性的、它是 有界的(上界 Theta 讓 B(u, v) 不超過某常數乘以 ||u|| ||v||)、它是 強制的(B(u, u) >= alpha ||u||^2)。拉克斯-米爾格拉姆定理 說這三件事實就是 你所需要的全部:對任何合理的右端 f,都存在唯一一個空間中的 u,使得對每個 v 都有 B(u, v) = (f, v)。弱解的存在性與唯一性,一舉交付,沒有 u 的公式、沒有巧妙的構造——只憑雙線性、有界性與強制性。

  1. 從 L u = f 出發,其中 L 在散度形式下一致橢圓,且帶零邊界條件。
  2. 乘上一個在邊界上為零的測試函數 v,在區域上積分,再做一次分部積分,落入弱形式 B(u, v) = (f, v)。
  3. 把問題擺進希爾伯特空間 H^1 sub zero,使 u 只需要一個弱導數。
  4. 由橢圓性上界 Theta 驗證有界性,再由下界 theta 加上龐加萊不等式驗證強制性。
  5. 援引拉克斯-米爾格拉姆:唯一的弱解 u 存在。完成——不需要任何公式。

值得細細品味這有多便宜。古典做法會去逐區域獵捕一個顯式的 u;能量法 則只驗兩個不等式,便讓 拉克斯-米爾格拉姆 把存在性整批交來。當算子 對稱 時(A 對稱、沒有一階項),同一枚硬幣還有一張更漂亮的面:弱解恰好就是某個能量積分的極小化元——這就是 狄利克雷原理,它正是通往本階第五篇變分法的橋。

你買到了什麼,又還欠著什麼

對這筆交易要誠實。拉克斯-米爾格拉姆給的是一個 解——一個帶一個導數的索伯列夫空間成員——而非古典的、二次可微的解。它還沒告訴你這個解是平滑的。把平滑性贖回來是另一套工程,橢圓正則性,它是第三篇的核心:它證明對良好的係數與資料,弱解其實暗地裡要多平滑有多平滑,所以那份弱只是一筆暫時的借款,會如數歸還。在那之前,「解」誠實地就指「弱解」。

方法本身也有一道邊界。強制性可能失效——能量碗可能變平——比方當一個低階項把算子傾斜,或在某些特殊頻率上、齊次問題 L u = 0 帶有非平凡解之時。那時拉克斯-米爾格拉姆便不再直接適用,存在性化為一個由 弗雷德霍姆二擇一 所支配的精微的非此即彼:要嘛 L u = f 永遠可解,要嘛齊次方程有非零解、擋住了某些 f。那個二分法,以及標出強制性究竟在何處崩潰的特徵值,正是第二篇與第四篇的題目。

並且,把型別認誠實。這套能量法之所以屬於 橢圓型 問題,正因為 B(u, u) 是一個正定的能量——就是當初讓這個型成為橢圓型的那同一份平衡性格。試著把類似的論證搬到波動方程,那個「能量」便不再是單一正負號;試著把它逆著時間搬到熱方程,它便爆掉。橢圓性在這裡不是裝飾;它是承重牆,而一致橢圓性正是讓這道牆在區域裡任何地方都不致裂開的東西。