為何我們還需要更多工具
到目前為止,你能用三種方式求解熱方程:在區間上用分離變數法、在直線上用熱核,以及對平滑化、無限速度與不可逆性的定性解讀。但這每一種都倚賴把解寫出來——一個傅立葉級數、一個高斯卷積。仍有兩個誠實的缺口。第一:你能否不必求解、只憑資料,就對一個解說出某些確定的事?第二:當棒子內部有一個加熱器在運轉,使方程帶有一個源項、那些乾淨的公式不再適用時,會發生什麼事?
本指南把兩個缺口都補上。最大值原理回答第一個:它只用初始與邊界資料就把解處處界定住,立即見效,完全不需要解本身的公式。杜阿梅原理回答第二個:它把任意源視為一連串微小的初始脈衝,重複使用你已信賴的無源求解器,從而構造出對源的回應。兩者合起來,正是整個梯級在定性與構造上的壓軸。
最大值原理:尖峰住在拋物邊界上
把一根金屬棒加熱,固定它的兩端,然後觀察。內部沒有任何一點會自發地變得比一開始更熱、或比兩端維持的溫度更熱——你會看到的最高溫,早已存在於初始分布或邊界之上。擴散只會抹平,它絕不會憑空製造出新的尖峰。拋物型最大值原理對熱方程講的正是這件事:在一個時空區域上,u 的最大值(與最小值)在它的邊界上達到,而非在內部。
微妙而關鍵的一點在於「哪一部分」邊界算數。對一個在區間上、跨越從 t = 0 到 t = T 時間窗的問題,相關的集合——稱為拋物邊界——是底部(初始切片 t = 0)連同側面(兩端,對所有時刻),但「特別不」包含頂部(最終切片 t = T)。原理回溯到初始資料與側面資料;它永遠不會前進到最終時刻。這一個不對稱——過去重要、未來不重要——正是擴散的不可逆性以純粹幾何的面貌現身。
t = T (final time) <-- NOT part of the parabolic boundary
+-----------------+
| |
x = 0 | interior: | x = L <-- sides ARE part of it
(side) | no new max | (side)
| forms here |
+-----------------+
t = 0 (initial slice) <-- bottom IS part of it
max of u = max over { bottom (t=0) + sides (x=0, x=L) }
EXCLUDES the top edge t = T為何它成立,以及它立刻為你換來什麼
證明是一個簡短而漂亮的把戲,值得看一次。假設 u 滿足 u_t = k u_xx 且在某點有一個嚴格的內部最大值。在空間上的內部最大值處,分布向下彎,所以 u_xx <= 0;而在時間上的最大值處,u 不再上升,所以 u_t >= 0。但方程要求 u_t = k u_xx 且 k > 0,這會迫使一個非負量等於一個非正量——只有兩者皆為零才可能,於是排除了嚴格的尖峰。要解決臨界情形,改看 v = u - epsilon t,其中 epsilon > 0 極小:則 v_t - k v_xx = -epsilon < 0,於是 v 根本不可能有內部最大值,它的最大值落在拋物邊界上,再讓 epsilon 趨於零,就把結論交還給 u。
從這一個原理,三項獎賞一次掉出。唯一性:兩個初始與邊界資料相同的解,其差 w 滿足熱方程且拋物邊界資料為零,於是它的最大與最小皆為零,故 w = 0——熱方程的初邊值問題恰有唯一一個解。穩定性:把資料移動至多 delta,解就處處移動至多 delta,這正是使問題適定的連續依賴。還有免費的界限:整個內部都被拋物邊界上的最大與最小值框住,這是一個你在求解之前就握有的真正先驗估計。
最有用的重新包裝是比較原理:較大的資料給出較大的解,逐點如此。若一根棒子在每一點的起始與邊界都比另一根熱,它就處處、永遠保持較熱——把最大值原理套用在兩個解的差上即可。一個乾淨的特例是正性保持:非負的初始與邊界資料迫使解非負,這正是為什麼溫度(或機率密度、或濃度)在擴散下絕不會變負。常識,終於得證。
杜阿梅原理:在加熱器開著時組裝解
現在來補第二個缺口。直到此刻,棒子都沒有內部熱源——只有一個散開的初始分布。但假設內部開了一個加熱器,隨時間持續加入熱量,於是方程變成非齊次的 u_t = k u_xx + f(x,t),其中 f 是源。高斯卷積與裸傅立葉級數是為無源方程打造的;它們並不直接涵蓋 f。杜阿梅原理就是那座優雅的橋,而它完全建立在熱方程豐沛擁有的那一個性質上:線性,也因此是疊加原理。
這個想法是一幅生動的圖像。把運轉中的源想成不是一次連續的推動,而是一條不間斷的微小瞬間脈衝之流。在某瞬間 s 附近一個短暫間隔內注入的熱,表現得像一個在時刻 s 開啟的全新初始條件 f(x,s) ds——而從那一刻起,它就只是自行向前擴散,在剩下的時間 t - s 裡服從無源方程。每個脈衝都忘了自己曾是源,只是像普通的初始資料那樣散開。由疊加,完整的受迫解就是所有那些獨立擴散脈衝之和(積分)。
- 在某瞬間 s 把源凍結:讀出它在那裡送出的熱 f(x,s) ds,把那一片當作在時刻 s 開啟的全新初始條件。
- 讓那個凍結的脈衝在無源熱方程下,對剩餘的時間 t - s 自由擴散——用的正是你已經精通的求解器(熱核,或特徵函數級數)。
- 對從 0 到 t 的每個點火時刻 s 把貢獻加總。用 S(tau) 表示對時間 tau 的無源演化,則受迫解為 u(x,t) = integral from 0 to t of S(t - s) f(x,s) ds。
- 若初始資料不為零,再用疊加加上單獨演化的初始條件——杜阿梅處理源,普通演化處理起點,兩者直接相加。
用核寫清楚,對加熱器的回應是 u(x,t) = integral from 0 to t of [ integral of G(x - y, t - s) f(y,s) dy ] ds,其中 G 是熱核。這正是你在 ODE 學過的參數變異公式在 PDE 上的對應,而內層的物件 G(x - y, t - s) 是時間方向的因果格林函數:它只活在 t - s > 0 之處,編碼了「現在的溫度只能依賴過去加入的熱、絕不依賴未來」。
誠實的限度,以及兩個原理如何契合
兩件工具都附帶誠實的細則。杜阿梅原理徹頭徹尾是一個線性陳述:它之所以奏效,全靠對個別脈衝的回應能夠疊加。對一個非線性源——比如 f 依賴於 u 自身,如同反應擴散模型——乾淨的求和就崩解了,你必須把杜阿梅公式變成一個積分方程、朝一個不動點迭代,而非一筆讀出答案。它也悄悄假設無源流本身是適定的,這正是為什麼你能把它向前時間跑、卻絕不能向後。
最大值原理也有它自己的限度,而它們正是那條把擴散與波分開的界線。它根本上是一個拋物與橢圓型的現象,生於擴散那種平滑、保序的本性。波動方程完全不服從最大值原理:它的解在資料的上下自由振盪,一個內部的波峰能升得比邊界上任何東西都高,比較原理徹底失效。若你哪天發現自己對波也期待「較大的資料給出較大的解」,停下——那直覺屬於擴散,不屬於波。