熵條件與唯一性

一個跳躍條件還不夠

在上一份指南裡，我們贏得了朗肯-雨貢尼奧條件：在一個純量守恆律 u_t + f(u)_x = 0 中，一個不連續面必須以速度 s = (f(u_R) - f(u_L)) / (u_R - u_L) 行進，也就是連接左右兩個狀態的割線斜率。這個單一的方程看起來似乎應該把一切都定下來——給我左值與右值，震波速度便自動跳出。但它底下藏著一樁無聲的醜聞。

醜聞在於：跳躍條件對於走錯方向的跳躍也照樣成立。想像一份資料本來就是一道平滑的斜坡——左邊低、右邊高——這種會往外攤開的流動，按物理本該只是扇形展開。然而朗肯-雨貢尼奧的代數卻欣然接受把這份資料黏成一個立在中間的尖銳不連續面，並稱之為解。平滑的扇形與立著的跳躍，兩者都在弱解的意義下滿足同一條守恆律。唯一性——每個適定問題都向我們承諾過的東西——已悄然蒸發。

兩個並排競爭的弱解

讓我們用伯格斯方程 u_t + (u^2 / 2)_x = 0 把這樁醜聞痛切地具體化，它的通量是 f(u) = u^2 / 2，特徵速度是 f'(u) = u——每個點都以自己的數值作為速度。取資料左邊 u_L = 0、右邊 u_R = 1：後面是慢的、前面是快的。因為快的數值在前頭、把慢的甩開，特徵線彼此分離，張開一道沒有任何特徵線抵達的縫隙。必須有某樣東西來填補它。

如今有兩個誠實的答案爭奪這一份資料。候選 A 是一個立著的震波：在某個跳躍的左邊保持 u = 0、右邊保持 u = 1，並讓這個跳躍以朗肯-雨貢尼奧速度 s = (f(1) - f(0)) / (1 - 0) = 1/2 行進。它恰好滿足跳躍條件，所以是一個完全合法的弱解。候選 B 則是一扇填補縫隙的連續扇形：x <= 0 時 u = 0，0 < x < t 時 u = x/t，x >= t 時 u = 1。它同樣滿足偏微分方程，而且處處沒有跳躍。同一份資料、兩個不同的解——光靠方程本身，根本無法在它們之間做出選擇。

候選 B 是一個稀疏波：它不強迫產生跳躍，而是讓那道縫隙被一扇連續的數值扇形 u = x/t 填滿，由 0 平滑地內插到 1。它正是這樣建構的，使得一條全新的特徵線穿過那塊缺失楔形的每一個點。候選 A 則是一個被跳躍條件認證、卻被物理拒絕的不連續面。我們需要一個決勝規則——而關鍵在於，這個決勝規則必須判 B 勝出、判 A 出局。

熵條件：特徵線必須撞進震波

這裡是決定一切的圖像。把震波畫成 (x,t) 平面上的一條線，並把兩側的特徵線畫成一支支小箭頭。一個不連續面唯有在特徵線從兩側流入它時，才是物理上可接受的——資訊湧進震波並在那裡被摧毀，就像車輛擠進一場交通堵塞。一個有特徵線從中向外發散、把資訊向外噴灑的跳躍，則是被禁止的「膨脹震波」：你無法憑空從一個赤裸的不連續面製造出資訊。

把那幅圖像翻成代數。震波正後方的特徵速度是 f'(u_L)、正前方是 f'(u_R)，而震波本身以速度 s 緩行。「特徵線從兩側流入」意味著後方速度大於 s、而 s 大於前方速度。對伯格斯方程這種凸通量而言，這就是拉克斯熵條件：f'(u_L) > s > f'(u_R)。對伯格斯方程它化簡為 u_L > u_R——震波唯有在後方數值大於前方數值時才可接受。

Lax entropy condition (convex flux f):    f'(u_L)  >  s  >  f'(u_R)

   characteristics behind:  speed f'(u_L)  -->|
                                              |  shock at speed s
   characteristics ahead:                  |<--  speed f'(u_R)
         (both arrows aim AT the shock line:  admissible)

Burgers (f = u^2/2, f' = u):   condition becomes   u_L > u_R

   u_L=0, u_R=1  -->  0 > 1 is FALSE  -->  Candidate A is rejected
                      (so the rarefaction fan B is the right answer)

一行寫下的拉克斯熵條件，以及它為何在 u_L=0、u_R=1 時丟掉那個膨脹震波。

現在讓我們的兩個候選通過這道關卡。對 u_L = 0、u_R = 1，拉克斯檢驗要求 0 > 1——擺明是假的——所以候選 A，那個立著的震波，是不可接受的，被丟出去。候選 B，那個稀疏波，根本沒有跳躍，存活下來。平局正如我們所願地被打破了：平滑的物理扇形勝出，虛假的震波落敗。把資料翻成 u_L = 1、u_R = 0，判決也跟著翻轉——此時 1 > 0 成立，震波是可接受的，而它才是唯一的物理答案。

規則的來歷：一絲黏性的低語

拉克斯條件不是從帽子裡變出來的。它是我們刪掉的那段物理所留下的指紋。真實的流體與真實的車流都帶有一點點擴散——一個微小的平滑項。於是考慮正則化方程 u_t + f(u)_x = epsilon u_xx，其中 epsilon > 0 很小。這個拋物型方程擁有熱方程那種平滑的力量：它完全適定、擁有平滑的唯一解、且永遠不會形成真正的不連續。那道本會變成震波的懸崖，反而被磨圓成一層薄而陡的過渡帶。

現在讓 epsilon 縮向零。那層薄薄的平滑過渡帶逐漸陡峭成一個真正的跳躍，而你抵達的極限，就是消失黏性解。深刻的定理是：這個極限是唯一的，而它恰好就是滿足熵條件的弱解——對攤開的資料是候選 B，對壓縮的資料是震波。因此熵條件並非裁判任意的判罰。它是真實物理之極限所留下的疤痕，是唯一一個身為平滑解之影子的弱解。

唯一性復原，以及該帶走什麼

把熵條件鎖上之後，這本帳終於平衡了。在成群的弱解之中，恰好有一個滿足熵條件，而那一個既存在、又對資料連續依賴。對任意維數的純量守恆律而言，這正是著名的克魯日科夫唯一性定理的內容：滿足熵條件的弱解是唯一的。我們合上了那個迴圈——它早在平滑特徵線首次相交、一個震波誕生之時便已開啟。

給定候選不連續面兩側的左右狀態，計算朗肯-雨貢尼奧速度 s = (f(u_R) - f(u_L)) / (u_R - u_L)。
計算後方與前方的特徵速度 f'(u_L) 與 f'(u_R)。
檢驗拉克斯不等式 f'(u_L) > s > f'(u_R)：特徵線是否從兩側流入這個跳躍？
若是，震波可接受，即為你的解。若否，丟棄它，改用稀疏扇形填補縫隙。

請牢牢記住兩條誠實的提醒。第一，乾淨的拉克斯形式 f'(u_L) > s > f'(u_R) 是凸通量的版本；對非凸通量 f，你需要更強的奧列尼克弦條件，而一個單一跳躍可能必須被一組震波與扇形的複合所取代。第二，那令人安心的克魯日科夫理論是純量的故事——對守恆律系統而言，存在性與唯一性遠為微妙，部分問題至今仍真正懸而未決。手握這條選擇規則，緊接著的下一份指南會把一切——熵震波與稀疏扇形——組裝成分段常數資料的完整答案：黎曼問題。