Euler 方程式與 Navier-Stokes 方程式的差別:黏性、耗散與正則性
Euler 方程式忽略粘性。Navier-Stokes 方程式包含粘性。這一個差異改變了物理、數學,以及那個百萬美元問題的面貌。
簡短回答
Euler 方程式描述一種內部摩擦為零的流體。完全沒有粘性。Navier-Stokes 方程式描述同一種流體包含粘性。
在數學上,整個差別就是一項:$\nu \Delta u$,也就是粘性擴散項。移除它,Navier-Stokes 就變成 Euler。保留它,方程式便多了一個平滑化機制;令人意想不到的是,僅僅多出這一項,就會同時改變物理與分析。
正是這一項,解釋了煙霧為何會消散、邊界層為何會沿著表面形成,以及Navier-Stokes 千禧年問題為何具有與對應 Euler 問題完全不同的性質。
在$\mathbb{R}^3$上的非壓縮 Euler 方程式為
$$\partial_t u + (u \cdot \nabla)u = -\nabla p, \qquad \nabla \cdot u = 0.$$
非壓縮 Navier-Stokes 方程式為
$$\partial_t u + (u \cdot \nabla)u = -\nabla p + \nu \Delta u, \qquad \nabla \cdot u = 0,$$
其中$\nu > 0$是運動粘度。形式上,在 Navier-Stokes 中令$\nu = 0$即可得到 Euler。但這個形式代換掩蓋了一個事實:$\nu \to 0$極限是奇異的:粘性項$\nu \Delta u$承載了系統中最高階的空間導數,而刪去它會改變 PDE 的型態,以及解所在的函數空間。
兩個方程式並列比較
以下把兩個方程式寫成標準非壓縮形式,使比較一目了然:
Euler:
$$\partial_t u + (u \cdot \nabla)u = -\nabla p$$
Navier-Stokes:
$$\partial_t u + (u \cdot \nabla)u = -\nabla p + \nu \Delta u$$
左邊相同:速度的變化率加上非線性自輸運項$(u \cdot \nabla)u$。兩者都透過$\nabla \cdot u = 0$強制非壓縮性。唯一的結構差異是粘性項$\nu \Delta u$,它位於 Navier-Stokes 的右邊。
參數$\nu$是運動粘度,即流體的一個物理常數。蜂蜜的$\nu$很大。空氣的則很小。Euler 方程式對應於$\nu = 0$:一種完全無摩擦的理想化模型,可用來近似某些遠離邊界的高 Reynolds 數流動,但在任何真實流體中都不存在。
兩個系統共享雙線性形式$(u \cdot \nabla)u$以及由無散度約束隱含決定的壓力。對動量方程式取散度並使用$\nabla \cdot u = 0$可得到壓力 Poisson 方程式
$$-\Delta p = \nabla \cdot ((u \cdot \nabla)u) = \partial_i \partial_j (u_i u_j),$$
這對兩個系統都是相同的。無論哪一種情況,壓力都是$u$的非局部泛函。
粘性項$\nu \Delta u$是作用在每個速度分量上的二階線性橢圓算子。它是拋物型正則化:它的存在使 Navier-Stokes 成為一個半線性拋物型系統,而非壓縮 Euler 則是一個具有非局部壓力耦合的一階非線性輸運方程式。這種 PDE 型態上的差異,推動了正則性理論中幾乎所有後續差別。
粘性在物理上做了什麼
粘性是流體相鄰層之間的摩擦。快速的一層旁邊有慢速的一層?粘性會在它們之間傳遞動量,抹平速度差。概念很簡單,後果卻極其重大,並以三種方式把 Navier-Stokes 與 Euler 區分開來。
- 耗散。動能轉化為熱。攪拌咖啡,然後停下來。它最終會靜止,因為粘性把運動作為熱能耗散掉。Euler 完全無法預測這點,因為方程式中沒有把動能排入熱的機制。
- 邊界層。真實流體會黏附在表面上(無滑移條件),在壁面附近產生速度快速變化的薄層。這些會造成飛機機翼上的阻力、管道中的摩擦損失,以及高速下湍流的開始。Euler 流動則滿足滑移條件,因此完全漏掉粘性壁面阻力。
- 小尺度平滑化。粘性會消除最尖銳的速度梯度。沒有它呢?就沒有任何東西阻止流動發展出無限精細、越來越尖銳的結構。這種平滑化正是使Navier-Stokes 的正則性問題成為不同於 Euler 的另一種難題的原因。
在$\mathbb{R}^3$(或週期區域)上的 Navier-Stokes 能量恆等式為
$$\frac{d}{dt}\frac{1}{2}\|u\|_{L^2}^2 = -\nu \|\nabla u\|_{L^2}^2,$$
因此動能會單調耗散。對 Euler($\nu = 0$)而言,$L^2$範數的$u$在形式上守恆。
在邊界層面,Navier-Stokes 使用無滑移條件($u|_{\partial \Omega} = 0$),而 Euler 只要求不可穿透性($u \cdot n = 0$)。當$\nu \to 0$時,這些條件之間的不匹配會產生 Prandtl 邊界層。這是一種奇異攝動現象,自 Prandtl 1904 年的論文以來,人們一直在與它搏鬥。
從物理上看,粘性就像高頻濾波器:它以速率$\nu |k|^2$阻尼 Fourier 模態,優先抑制小尺度。這種譜阻尼是湍流中 Kolmogorov 耗散尺度$\eta \sim (\nu^3 / \varepsilon)^{1/4}$背後的機制。見Reynolds 數與湍流 以了解完整的尺度變換圖像。
Euler 只是粘性為零的 Navier-Stokes 方程式嗎?
形式上?是的。令 $\nu = 0$ 就得到 Euler。但如果思考停在這裡,那就大錯特錯了。
極限 $\nu \to 0$ 是奇異的。粘性承載著方程式中的最高階導數,所以移除它並不是做一個小修正,而是會徹底改變你所面對的 PDE 類型。邊界層不會優雅地逐漸變薄;它們可能爆發成湍流。在 Navier-Stokes 方程式下原本完全光滑的解,在 Euler 方程式下可能發展出截然不同的行為。
沒錯,這兩個方程式共享數學上的 DNA。但零粘性極限是整個流體力學中最深刻的未解問題之一,不是餐巾紙上的計算。
無粘性極限 $\nu \to 0$ 是一個奇異擾動:$\nu \Delta u$ 承載系統中的最高階空間導數,所以令 $\nu = 0$ 會降低 PDE 的階數。在有邊界的區域上,這個極限與 Prandtl 邊界層展開的有效性相連,而該展開可能以驚人的方式失效(Grenier 2000, Gérard-Varet & Dormy 2010)。
在 $\mathbb{R}^3$ 或 $\mathbb{T}^3$(無邊界)上,情況會乾淨些。如果 Euler 解 $u^E$ 在 $[0,T]$ 上保持光滑,則 Navier-Stokes 方程式的解 $u^\nu$ 收斂到 $u^E$ 於 $L^2$ 當 $\nu \to 0$(Kato 1972)。比較乾淨,但仍不算乾淨:Euler 解是否能在 3D 中大域保持光滑,本身就是未解問題,而這正是整個問題的核心。
這個極限也與湍流理論正面相撞。Kolmogorov 的圖像需要 $\nu > 0$ 來定義耗散尺度。然而異常耗散,也就是當 $\nu \to 0$ 時能量耗散仍持續存在,數十年來一直是未解猜想。Onsager 猜想(如今已成定理:Isett 2018,並由 Buckmaster、De Lellis、Székelyhidi 與 Vicol 於 2019 年強化)精確刻畫了 Euler 解何時能在完全沒有粘性的情況下耗散能量。
人們什麼時候會用 Euler 而不是 Navier-Stokes 方程式?
只要粘性相較於其他作用力可以忽略,就會這麼做。這比你想像的更常見:
- 高速空氣動力學 中遠離表面的區域。遠離機翼時,氣流幾乎是無粘性的。工程師經常用 Euler 求解器處理主體流動,再在壁面附近補上邊界層修正。
- 天體物理流。 星際氣體雲、恆星內部、黑洞周圍的吸積盤。在這些尺度下,分子粘性完全無關緊要(儘管湍流有效粘性未必如此)。
- 壓縮氣體動力學。 激波。爆震。超音速飛行。主導物理的是壓力與慣性,而不是摩擦。
- 純理論。 Euler 本身就值得研究,而不只是通往 Navier-Stokes 方程式的踏腳石。它連結到 Riemannian 幾何、渦旋動力學,以及關於湍流結構本身的深層問題。
但只要摩擦、阻力或邊界行為重要(管流、靠近表面的車輛空氣動力學、血液循環、人類尺度的天氣),你就需要 Navier-Stokes 方程式。就這樣。
Euler 方程式支配高 Reynolds 數 $\mathrm{Re} = UL/\nu \gg 1$ 下的主導階行為。在這些 Reynolds 數下,粘性效應被限制在薄邊界層與內部剪切層中,而遠離壁面的主體流動可由 Euler 方程式很好地近似。
壓縮 Euler 方程式是一個具有有限傳播速度的雙曲系統,是氣體動力學的標準模型,包含激波形成與 Riemann 問題。它們不同於上面討論的非壓縮 Euler 方程式:壓縮 Euler 確實是雙曲型,而非壓縮 Euler 具有非局部壓力耦合與無限傳播速度。
在數學分析中,Euler 既作為消失粘性問題的極限對象,也是一個豐富的 PDE 系統,擁有自身的正則性理論、守恆量(螺旋度、透過微分同胚群上的 Euler-Arnold 框架得到的 Casimir 量),並與幾何力學相連。
這種差異對正則性的意義
這正是差距最重要的地方,也是真正有趣之處。
Navier-Stokes 千禧年問題 問的是一個聽起來幾乎過於簡單的問題:如果從三維中一個光滑、良態的流動開始,解會永遠保持光滑,還是可能爆發?地球上沒有人知道答案。
Euler 的同一個問題在 3D 中也仍未解。但這兩個問題的感覺完全不同:
- Navier-Stokes 方程式有粘性站在它這一邊。總是在平滑,總是在耗散能量,總是在抑制最尖銳的梯度。真正的問題是,這種平滑是否足夠強,能在線性項製造奇異性之前壓過它。
- Euler 什麼都沒有。零平滑。零耗散。非線性項可以毫無反作用力地放大速度梯度,而這是否真的會從光滑的 3D 初始資料產生有限時間奇異性,是 PDE 理論中最大的未解問題之一。
在 2D 中,兩個方程式對光滑初始資料都是大域適定的。已解決。結束了。謎團完全存在於三維,對兩個方程式都是如此,但原因根本不同。
正則性圖像:
2D: 對兩個系統而言,光滑解的大域存在性與唯一性皆已知。對具有光滑資料的 2D Euler,Wolibner(1933)在 Hölder 空間中證明了大域存在性;Yudovich(1963)對有界渦量資料建立了唯一性。對 2D Navier-Stokes 方程式,大域正則性則由 Ladyzhenskaya 不等式與渦量最大值原理推出。
3D Navier-Stokes: Leray(1934)證明了 $L^2$ 中弱解的大域存在性,但唯一性與正則性仍未解。Caffarelli、Kohn 與 Nirenberg 定理(1982)表明奇異集的一維拋物 Hausdorff 測度為零,因此若爆發發生,也會極其稀疏。粘性項提供關鍵的先驗估計 $\int_0^T \|\nabla u\|_{L^2}^2 \, dt \leq C(u_0)$,但這種 $H^1$ 控制對 3D 尺度變換而言是次臨界的,還不足以閉合 bootstrap 論證。參見 為什麼 Navier-Stokes 方程式很難 了解超臨界性缺口。
3D Euler: 對光滑資料尚無大域理論。在 Sobolev 空間 $H^s$,$s > 5/2$ 中的局部適定性是經典結果(Kato 1972,Kato and Ponce 1988)。Beale、Kato 與 Majda 判準(1984)將爆發偵測化約為 $\int_0^T \|\omega\|_{L^\infty} \, dt$:解在 $[0,T]$ 上保持光滑,若且唯若此積分有限。爆發要求渦量增長得足夠快,以至於在時間上不可積。Elgindi(2021,Annals of Mathematics)證明了對於 $C^{1,\alpha}$ 資料會形成有限時間奇異性。這是真正的突破,但仍低於光滑($C^\infty$)門檻。光滑 Euler 解在三維中是否會爆發,仍是未解問題。
粘性、湍流與串級
湍流。這正是 Euler 與 Navier-Stokes 的比較在物理上變得鮮明、幾乎可觸之處。
在湍流中,能量從大尺度(管道、機翼、風暴的尺寸)進入,並向下串級到越來越小的渦旋。這就是能量串級,也是整個物理學中最驚人的現象之一。在串級的最底端,粘性終於將動能轉化為熱。到此為止。
Euler 捕捉的是慣性區間動力學:由非線性驅動的跨尺度能量傳遞。但它沒有粘性截止。串級沒有底端。也沒有在任何確定尺度上把動能轉化為熱的機制。在無粘性極限中能量是否仍會耗散,也就是所謂的異常耗散,仍是一個深刻的未解問題。
這就是為什麼湍流建模幾乎總是使用 Navier-Stokes。Reynolds 數 $\mathrm{Re} = UL/\nu$告訴你串級有多寬:高 $\mathrm{Re}$ 意味著能量輸入與粘性耗散之間相隔許多個數量級的尺度。真實湍流存在於這種張力之中:無粘性串級把能量向下傾注,而粘性截止在最小尺度上將其摧毀。
Kolmogorov 1941 年理論預測在慣性區間中的能量譜為 $E(k) \sim \varepsilon^{2/3} k^{-5/3}$,此區域既不由大尺度外力主導,也不由粘性耗散主導。耗散尺度 $\eta = (\nu^3/\varepsilon)^{1/4}$ 設定了這個區間的底端。低於 $\eta$,粘性占上風。
Euler 方程式是形式上的 $\nu = 0$ 模型。但在湍流區域中,將 Navier-Stokes 解詮釋為收斂到 Euler 解,確實相當微妙;而能量耗散在此極限中是否持續存在(異常耗散)正是 Onsager 猜想的內容。剛性的一面,即證明對於 $u \in C^{0,\alpha}$ 且 $\alpha > 1/3$ 時沒有耗散,由 Constantin、E 與 Titi 於 1994 年證明。柔性的一面,即在低於 $C^{1/3}$ 時構造具耗散的 Euler 解,則由 Isett 於 2018 年完成,建立在 De Lellis 與 Székelyhidi 的凸積分計畫之上。
對 Navier-Stokes 而言,串級圖像直接嵌入能量平衡之中:$\varepsilon = \nu \langle |\nabla u|^2 \rangle$。未解問題十分鮮明。Navier-Stokes 解是否能保持足夠久的光滑性,使 Kolmogorov 的統計理論在數學上得到正當化?存在性與光滑性問題部分地問的正是這件事:粘性耗散是否強到足以在每一個尺度、對所有時間馴服串級。
總結:一個項,兩個不同世界
Euler 與 Navier-Stokes 的差別是一個項:$\nu \Delta u$。那個項改變了一切。
| Euler | Navier-Stokes | |
|---|---|---|
| 粘性 | 無($\nu = 0$) | 有($\nu > 0$) |
| 能量 | 守恆(形式上) | 耗散 |
| 邊界層 | 無 | 有(無滑移) |
| PDE 類型 | 一階非線性 + 非局部壓力 | 二階拋物型 + 非局部壓力 |
| 二維正則性 | 已解決 | 已解決 |
| 三維正則性 | 未解 | 未解(千禧年問題) |
Euler 不是簡化版的 Navier-Stokes。它是一個根本不同的系統,只是恰好共享了大部分結構。而這個選擇在實務上很重要:選錯模型(在粘性重要時用 Euler,或在粘性不重要時用 Navier-Stokes)可能會徹底毀掉一次模擬。完整方程式見Navier-Stokes 方程式是什麼?。關於障礙,見為什麼它很難。關於獎項,見千禧年問題。關於非壓縮與壓縮流,見非壓縮與壓縮 Navier-Stokes。
粘性項 $\nu \Delta u$ 將一個一階非線性系統(Euler)轉化為一個半線性拋物型系統(Navier-Stokes)。你得到的是:能量耗散、高階先驗估計,以及作為 Navier-Stokes 理論基礎的解析半群結構;該理論由 Leray 於 1934 年建立,並由 Fujita 與 Kato 於 1964 年延伸。
然而,這種正則化在三維中仍不足以完成大域正則性的論證。甚至遠遠不夠。Clay 千禧年問題問的正是 Navier-Stokes 中的拋物型平滑是否能對所有時間控制所有尺度上的非線性能量傳遞。對 Euler 而言,平行的問題同樣尖銳:完全沒有平滑是否保證光滑資料會發生有限時間爆發?沒有人知道。
這兩個問題都位於流體數學理論的核心;比較它們,正好揭示粘性究竟帶來了什麼,以及沒有帶來什麼。無論對哪一個系統,三維正則性問題仍是整個分析學中最困難的未解問題之一。