一個時間用罄的解
到目前為止,這一階把一個反射操練得很熟:在非線性偏微分方程裡你不能把解相加,所以疊加原理那份老安慰沒了,而全新的行為也就被放了進來。上一篇指南展示了它的一張臉——一個反應擴散方程把自己組織成一道行波鋒面或一個圖靈圖案。這篇指南要展示相反的那張臉,那張暴烈的臉。有時候一個完美平滑、完全無辜的初始狀態會無界地增長,並在某個有限的時刻抵達無窮大。過了那一瞬,就根本沒有解可談了。我們稱這為有限時間爆破。
認識它最乾淨的地方,是半線性熱方程 u_t = Laplacian u + u^p,其中 p > 1,初始溫度非負。把它讀成一場拔河。拉普拉斯算子就是你早從熱方程認得的擴散:它把熱攤開、把尖峰壓平,獨自運作的話會讓一切平靜下來。而 u^p 這一項是一個源項——一個製造更多 u 的反應,而且因為指數 p 大於 1,u 越大它製造得越快。這篇指南的整齣戲,就是哪一邊會贏:是平息一切的擴散,還是自我放大的反應。
預言這場災難的那個小小常微分方程
為了看清「有限時間」真的可以是「有限時間」,先把擴散丟開一會兒,只留下反應。設想溫度在空間中處處相同,於是 Laplacian u = 0、什麼都不擴散。那麼 u 不再依賴位置,偏微分方程就塌縮成一個你能徒手解出的常微分方程:y'(t) = y^p。這是沒有任何冷卻的純反應——對熱來說是最壞的情況,對災難來說是最好的情況——而且它小到可以一行就積出來。
y'(t) = y^p , y(0) = a > 0 , p > 1 Separate: y^(-p) dy = dt Integrate: y^(1-p) / (1 - p) = t + C Use y(0)=a: C = a^(1-p) / (1 - p) Solve: y(t) = [ a^(1-p) - (p-1) t ]^(-1/(p-1)) Denominator hits 0 when t = T* = a^(1-p) / (p - 1) => y(t) -> +infinity as t -> T* (BLOW-UP at a finite time T*) Contrast p = 1: y' = y => y = a e^t (finite for every finite t)
好好盯著爆破時間 T* = a^(1-p) / (p - 1)。三個誠實的教訓直接掉出來。第一,爆破是真的,而且時刻是明確的——y 的分母在 T* 處歸零,y 衝向無窮,沒有任何含糊。第二,起點越熱爆破越早:較大的初始值 a 讓 a^(1-p) 變小(因為 1 - p < 0),於是 T* 縮短。第三,這公式在 T* 之後毫無意義;解並不會「再降回來」——它停止存在,這正是為什麼我們必須談一個最大存在時間,並在那裡打住。這個小小的常微分方程,就是整個現象的心跳。
把它變嚴格:用一個爆破的常微分方程把偏微分方程困在下方
懷疑論者應該立刻反對:你刪掉擴散才得到那個常微分方程,但真實的擴散仍在運轉,而擴散會冷卻。拉普拉斯算子難道不會洩掉足夠的熱,把解救回來嗎?對一個夠大的熱團來說,答案是不會,而證明這點的工具,就是你在拋物理論裡見過的比較原理。它的口號簡單而有力:對這些方程,若一個解一開始在另一個解之下、並遵守同一種不等式,它就會永遠待在下方。起點的大小順序,會在時間中保持。我們要用它,把真正的偏微分方程解夾在一個我們已知會爆炸的量之上。
- 挑一個候選的下方屏障。取上一節那個空間上為常數的常微分方程解 y(t),起始值 a 設為初始資料的最小值。它精確地滿足 y' = y^p,沒有擴散項。
- 驗證它是一個下解。把 y 代入偏微分方程算子:y_t - Laplacian y - y^p = y^p - 0 - y^p = 0,而因為 y 在空間中為常數,擴散項消失。所以 y 以等號滿足方程——它是 u_t = Laplacian u + u^p 一個如假包換的下解。
- 把起始資料排序。在 t = 0 我們把 a 選成真實初始溫度所取到的最小值,於是在空間中處處 u(x,0) >= y(0)。真正的解一開始就坐在屏障之上。
- 套用比較原理。因為 u 一開始就在下解 y 之上(或相等)、而兩者遵守同一個反應擴散不等式,所以只要 u 存在,就有 u(x,t) >= y(t)。真正的解被釘在屏障之上。
- 把結論讀出來。y(t) 在有限時間 T* 爆破,而 u 坐在它上面,所以 u 必定不晚於 T* 爆破。擴散沒能救下這個解;它只不過無法把 u 壓到一個正在爆炸的地板之下。
這就是整個證明,而它的優雅值得停下來品味。我們從未解出偏微分方程;我們只蓋了一個徒手就能解的卑微對象,並證明真正的解逃不到它底下去。這正是先驗估計的精神——在尚未、且不必真正建構出解之前,就抽出一個關於它的硬事實(這裡是:它必定在 T* 之前死去)。也請留意那份不對稱:這個論證需要初始熱團夠熱,熱到連它最低處的反應都贏。又小又涼的資料可能反而擴散到一無所有。你是爆破還是淡去,可以取決於你起始時東西的大小。
臨界指數:擴散何時能贏?
比較論證顯示大資料會爆破,但它留下一個更尖銳的開放問題:對小資料而言,擴散究竟能不能贏、讓解永遠活下去?由藤田(Fujita)找到的驚人答案是:這取決於單一一個數——一個臨界指數。對整個 n 維空間上的 u_t = Laplacian u + u^p,分界線是 p* = 1 + 2/n。把指數 p 拿來和這道門檻相比,越過它時判決就乾淨地翻轉。
用文字說,這道門檻講的是這件事。對 1 < p <= p*——次臨界情形——無論資料多小、多溫和,每一個非平凡的非負解都會在有限時間內爆破。對 p > p*——超臨界——夠小的資料能全域存在、活到永遠,而大資料仍可能爆破。把維度代進去,這個數就很具體:一維時 p* = 3(n = 1)、平面上 p* = 2(n = 2)、三維時 p* = 5/3(n = 3)。請留意,當 n 增大,門檻會朝 1 下降,所以較高的維度讓全域存活變得更容易——擴散有更多方向可以攤開。
好好體會那個次臨界區間有多麼怪異。當 1 < p <= p* 時,根本沒有任何安全的初始狀態:把資料取得要多小、多溫和都行,解依舊能找到辦法在有限時間內爆破。在低維度裡,擴散就是太弱了,弱到連一個微弱的反應製造熱的速度都排不掉。一旦越過臨界指數,天平就傾斜了:此時擴散在較高維度有更大的空間攤開,能夠贏過小資料,於是存在全時間的解——雖然一個夠大的熱團仍會被上面的比較論證炸掉。臨界指數,正是這兩個世界之間那道精確的刀鋒。
關於臨界情形本身,有一個誠實的但書。在這裡,邊界 p = p* 歸在爆破的那一側——恰在門檻上,連小資料都仍會爆破——但那個端點很微妙,是在藤田原始工作之後才被確定下來的。在整個非線性偏微分方程裡都成立的更廣教訓是:臨界指數很少是一道客氣的切點:邊界情形通常需要它自己一套細緻的論證,而且往往表現得像你最意想不到的那一側。當你遇到這樣一道門檻時,請把等號的情形當成一個獨立的問題,而不是一條腳註。
爆破實際長什麼樣——以及我們仍不知道的事
把這場爆炸畫在空間中、而不只是在時間中,會很有幫助。解往往不是處處同時爆破;它會集中。一根尖峰在某個區域上長起來——一個爆破點——而場的其餘部分保持有界,並且在那一點附近,輪廓常常把自己組織成一個自相似的形狀:同一個隆起,按剩餘時間 T* - t 的某些次方重新縮放高度與寬度,一面朝奇異時刻變尖、一面自我複製。這是非線性偏微分方程一再出現的把戲——在奇異點附近,用恰當的縮放放大進去,一個乾淨而普適的輪廓就會從混沌中浮現。
退一步留意:我們如今已見過一個平滑解能在有限時間內崩潰的兩種截然不同的方式,把它們分清楚是值得的。在這座階梯前頭的一階守恆律裡,值 u 維持有界,但特徵線交叉、斜率趨於無窮——那是梯度爆破,是震波的誕生,靠過渡到弱解來收拾。而這裡,在半線性熱方程裡,是值 u 本身衝向無窮。這是真正不同的兩場災難:一個是解能以較弱的意義存活下來的一堵陡牆,另一個則是一個徹底逃出數線的振幅。
用一個關於知識邊緣的誠實註腳作結。對這個乾淨的模型方程,爆破的故事大致已被理解,但你一旦把它弄複雜——不同的反應項、方程組、其他幾何——要分辨全域存在與爆破,很快就變得困難,而對許多重要的非線性偏微分方程而言,這仍然是真正開放的。頭號例子就在隔壁一階:三維那維-斯托克斯方程的平滑解究竟會不會在有限時間內爆破、還是活到永遠,是一個著名的未解難題。你方才學到的那套簡單的「常微分方程加比較」機制,至今仍是專家們第一個會抓起來的工具,正因為這個一般問題並沒有一般的答案。