非壓縮與壓縮 Navier-Stokes：密度、壓力與正則性

Navier-Stokes 方程式是一族系統。非壓縮流與壓縮流之間的差異並非表面上的。它會改變未知量、數學結構，以及未解問題。

物理上的區分：會改變的密度與不會改變的密度

「非壓縮與壓縮」歸根結柢就是密度。它保持常數，還是會改變？

試試看。把注射筒裝滿水，然後推動活塞。水會移動，但在日常條件下並不會明顯被壓縮。水對壓縮的抵抗非常強，因此把它視為非壓縮是一個極佳的近似。非壓縮。現在把那支注射筒裝滿空氣並封住末端。推入活塞，你會感覺到空氣退讓；同一質量的空氣在你拇指施壓下被壓縮到較小的體積中。這就是壓縮流。

在非壓縮流中，密度 $\rho$ 在整個流體中為常數，而每個微小流體團在空間中運動時都保持其體積。壓縮流則不同。密度成為一個變數，可以隨位置與時間而改變。噴射引擎周圍的空氣、爆炸中的氣體、大尺度的大氣：全都是壓縮的，也全都由密度變化所驅動。

為什麼要在意？因為這個區分會重塑 Navier-Stokes 方程式 的形式、它們所預測的內容，以及分析和求解它們的難度。

非壓縮假設把密度場 $\rho$ 設為整個流動區域中的正的常數。物理上，這表示每個物質體積元素在流映射下都保持其體積。壓縮情形則把 $\rho(x,t)$ 提升為一個完整的未知量，由其自身的演化方程式支配。

這不是同一系統的兩種記號。它們在未知量的數目、約束方程式的結構，以及壓力的性質上都不同。非壓縮系統有 $d+1$ 個純量未知量（$u$ 和 $p$ 在 $\mathbb{R}^d$ 中）；壓縮系統通常有 $d+2$ 個主要場（速度、密度，以及內能或溫度），壓力則透過狀態方程式決定。

這個區分也會改變壓力的 PDE 類型：在非壓縮模型中是橢圓型，壓力由瞬時的空間求解決定；而在具粘性的壓縮情形中，整個系統呈現混合雙曲—拋物型結構。無粘性或雙曲部分的聲學模態以有限聲速傳播，而粘性與熱傳則在數學上引入拋物型擴散。

非壓縮 Navier-Stokes 方程式

非壓縮 Navier-Stokes 方程式 描述密度為常數的流體。它們是出現在 Clay 千禧年問題 中的版本，也是本站聚焦的版本。

這個系統有兩個部分。動量方程式：

$$\partial_t u + (u \cdot \nabla)u = -\nabla p + \nu \Delta u + f$$

以及非壓縮約束：

$$\nabla \cdot u = 0$$

約束 $\nabla \cdot u = 0$ 表示速度場是 無散度：流體在任何地方既不堆積也不變稀。流入一個微小區域的量，必須以相同速率流出。這個單一條件取代了整個密度方程式。密度不變，因此不需要方程式來追蹤它。

壓力在這裡扮演特殊角色。它不是由熱力學定律（例如理想氣體定律）決定。相反地，它會在各處瞬時調整，以保持流動無散度。數學上，$p$ 解一個由約束推出的 Poisson 方程式。壓力變化以無限快的速度傳播。在非壓縮流中沒有「聲速」。

非壓縮 Navier-Stokes 系統有兩個未知場：速度 $u$ 和壓力 $p$。這種簡潔性具有欺騙性。非線性項 $(u \cdot \nabla)u$ 仍使三維中的系統極其困難。

在 $\mathbb{R}^3$ 上、運動粘性係數為 $\nu > 0$ 的非壓縮 Navier-Stokes 系統：

$$\partial_t u + (u \cdot \nabla)u = -\nabla p + \nu \Delta u + f, \qquad x \in \mathbb{R}^3,\; t > 0,$$

$$\nabla \cdot u = 0, \qquad u(x,0) = u_0(x).$$

無散度條件 $\nabla \cdot u = 0$ 是逐點約束，而不是演化方程式。它編碼了局部體積保持：流映射 $\Phi_t$ 滿足 $\det(D\Phi_t) = 1$ 對所有 $t$。

對動量方程式施加散度算子並使用非壓縮性，可得到壓力 Poisson 方程式：

$$-\Delta p = \partial_i \partial_j (u_i u_j) - \nabla \cdot f.$$

這是對每個固定時間的 橢圓型 方程式，未知量為 $p$。壓力不是獨立的熱力學變數；在標準 PDE 詮釋中，它作為 Lagrange 乘子，強制滿足無散度約束，並由速度場全域且瞬時地決定。資訊透過壓力以無限速度傳播，這是與壓縮系統不同的結構性差異，不能被輕描淡寫地掩蓋。

未知量是 $u : \mathbb{R}^3 \times [0,T) \to \mathbb{R}^3$ 和 $p : \mathbb{R}^3 \times [0,T) \to \mathbb{R}$。在 Clay 表述（Fefferman, 2000）中，$u_0 \in C^\infty(\mathbb{R}^3)$ 是無散度的，而問題是 $u$ 是否保持在 $C^\infty(\mathbb{R}^3 \times [0,\infty))$ 中並具有有界能量。

壓縮 Navier-Stokes 方程式

壓縮 Navier-Stokes 方程式 支配密度會變化的流動。更大的系統。更多未知量。更多方程式。

你仍然有一個動量方程式，但現在密度 $\rho$ 會明確出現：

$$\partial_t (\rho u) + \nabla \cdot (\rho u \otimes u) = -\nabla p + \nabla \cdot \tau + \rho f$$

約束 $\nabla \cdot u = 0$ 不見了。取而代之的是一個 連續方程式，用來追蹤密度如何演化：

$$\partial_t \rho + \nabla \cdot (\rho u) = 0$$

這表示質量守恆：密度會改變，是因為流動壓縮或膨脹流體團。

系統還需要一個能量方程式和一個 狀態方程式，以及像 $p = \rho R T$（理想氣體定律）這樣把壓力與密度、溫度聯繫起來的熱力學關係。壓力不再只是被動地施加約束。它的聲學部分以聲速傳播，而完整的黏性系統還包含擴散效應。

壓縮系統對於高速空氣動力學、天體物理氣體動力學、燃燒，以及任何密度變化很重要的流動都是不可或缺的。但它與非壓縮方程式是截然不同的數學物件。未知量更多、方程式更多，PDE 結構也完全不同。

壓縮 Navier-Stokes 系統耦合了速度 $u(x,t)$、密度 $\rho(x,t)$、壓力 $p(x,t)$，以及比內能 $e(x,t)$（或溫度 $\theta$）。以守恆形式寫為：

連續方程式： $$\partial_t \rho + \nabla \cdot (\rho u) = 0$$

動量： $$\partial_t (\rho u) + \nabla \cdot (\rho u \otimes u) + \nabla p = \nabla \cdot \tau + \rho f$$

能量： $$\partial_t (\rho E) + \nabla \cdot ((\rho E + p)u) = \nabla \cdot (\tau \cdot u) + \nabla \cdot (\kappa \nabla \theta) + \rho f \cdot u$$

其中 $E = e + \tfrac{1}{2}|u|^2$ 是總比能，$\tau$ 是粘性應力張量（對 Newtonian 流體，$\tau = \mu(\nabla u + \nabla u^T) + \lambda (\nabla \cdot u)I$，其中體積粘性係數為 $\lambda$），而 $\kappa$ 是熱導率。

封閉系統需要狀態方程，例如 $p = (\gamma - 1)\rho e$，對於絕熱指數為 $\gamma$ 的理想氣體。

在壓縮系統的無粘性或雙曲部分中，聲學擾動以有限速度傳播，其特徵聲速為 $c = \sqrt{\partial p / \partial \rho |_s}$。完整的黏性 Navier-Stokes-Fourier 系統則是混合雙曲—拋物型，因此粘性與熱擴散在數學上會引入無限傳播。壓縮方程式支援震波、稀疏波與接觸不連續；這些在非壓縮流中沒有對應物。

Mach 數：何時必須考慮可壓縮性？

何時必須考慮可壓縮性？由一個數字決定：Mach 數。

$$\text{Ma} = \frac{|u|}{c}$$

$|u|$ 是流速。$c$ 是聲速。它們的比值告訴你，流動相對於壓力擾動在介質中傳播的速度有多快；而這個比較決定了你是否可以安全地忽略密度變化，或是密度變化會主導物理。

當 $\text{Ma} < 0.3$ 時，密度變化通常不到大約 5%。這時非壓縮方程式通常可行。室內的空氣、管中的水、建築物周圍的風，都是低 Mach 數流動：壓力擾動傳播得比流動本身快得多，因此密度幾乎不動。

高於 $\text{Ma} \approx 0.3$ 時，可壓縮性開始變得顯著；而在約 $\text{Ma} \approx 1$ 時，會進入跨音速區域，局部超音速區塊出現並形成震波。戰鬥機。火箭噴嘴。重返大氣層的太空船。

這不是二元開關。大多數日常流體流動，以及 Clay 千禧年問題，都穩穩地位於低 Mach 數區域，在此非壓縮方程式適用。

Mach 數 $\text{Ma} = |u|/c$ 參數化了可壓縮性的重要程度，其中 $c = \sqrt{\partial p / \partial \rho |_s}$ 是等熵聲速。形式上，非壓縮方程式作為壓縮系統的低 Mach 極限而出現。

以 $\text{Ma}^2$ 的冪次所作的漸近展開（見 Klainerman & Majda，1981、1982；Schochet，1986）顯示，當 $\text{Ma} \to 0$ 且初始資料適當時，壓縮解會收斂到非壓縮解。在具備良好準備資料的標準低 Mach 無量綱化中，人們常將壓力寫成一個幾乎空間均勻的熱力學背景，加上一個較小的動態修正；後者在極限中施加非壓縮約束。

這是一個奇異極限：聲速 $c \to \infty$，而壓縮系統的雙曲特性退化為非壓縮流的橢圓壓力方程式。聲學模態變得無限快，並與渦旋動力學解耦。

$\text{Ma} < 0.3$ 這個區間是一條工程上的經驗法則。它反映了經驗觀察：在此範圍內，相對密度變化 $\delta\rho / \rho \sim \text{Ma}^2 / 2$ 會維持在約 5% 以下。數學上的正當性來自低 Mach 極限的收斂定理，而該定理需要良好準備的初始資料與相容的邊界條件。

為何千禧年問題討論的是非壓縮情形

Clay 千禧年問題提出了一個精確問題：給定 $\mathbb{R}^3$ 上一個光滑、無散度的初始速度，非壓縮 Navier-Stokes 系統是否總是產生一個對所有時間都存在的光滑解？

為什麼特別是非壓縮？有三個理由。

第一，它已經夠困難了。 自 Leray 於 1934 年的奠基工作以來，三維非壓縮方程式的大域正則性一直未能被證明。加入可變密度、熱力學與震波，只會讓問題困難得多，而不是更容易處理。

第二，困難是純粹的流體力學。 非壓縮系統隔離出核心數學挑戰：非線性平流 $(u \cdot \nabla)u$ 與粘性耗散 $\nu \Delta u$ 之間的競爭，且沒有熱力學或聲學上的複雜性。它是提出正則性問題最乾淨的場域。

第三，物理很清楚。 非壓縮方程式描述最常見的日常流動。它們是否能從光滑資料產生奇異性，是關於古典流體力學數學一致性的根本問題。

壓縮系統有其自身深刻的未解問題（大資料下大域解的存在性、震波的形成與相互作用），但那些是具有不同結構的不同問題。Clay 獎聚焦於非壓縮情形，因為那正是 Fefferman 為三維 Navier-Stokes 所表述的特定正則性問題。

Clay 官方表述（Fefferman，2000）指定的是 非壓縮 系統，定義在 $\mathbb{R}^3$ 上：

$$\partial_t u + (u \cdot \nabla)u = -\nabla p + \nu \Delta u, \qquad \nabla \cdot u = 0, \qquad u|_{t=0} = u_0,$$

其中 $u_0 \in C^\infty(\mathbb{R}^3)$ 無散度，而問題是是否 $u \in C^\infty(\mathbb{R}^3 \times [0,\infty))$ 且 $\int_{\mathbb{R}^3} |u(x,t)|^2\,dx$ 對所有 $t \geq 0$ 都有界。

選擇非壓縮系統有其數學動機。關鍵未解問題，也就是 Leray-Hopf 弱解（其大域存在，但可能不唯一或不光滑）與古典光滑解（其局部存在，但可能爆發）之間的落差，是三維非壓縮方程式所特有的。在二維中，光滑非壓縮 Navier-Stokes 解的大域正則性已知成立；三維情形仍然未解。

壓縮系統引入了性質上不同的困難：震波形成（即使對於 Euler 方程式 具有光滑資料）、真空狀態（$\rho \to 0$），以及渦度與聲學模態之間的耦合。這些都是重要的開放問題，但它們在結構上不同於非壓縮正則性問題。

非壓縮問題將能量超臨界非線性與粘性耗散之間的競爭單獨凸顯出來。自然的能量估計給出 $u \in L^\infty_t L^2_x \cap L^2_t \dot{H}^1_x$，在 3D 中比尺度臨界空間 $L^\infty_t \dot{H}^{1/2}_x$ 少了半個導數。彌合這個差距，或證明它無法被彌合，正是千禧年問題的核心。

接下來讀什麼

從這裡開始。想把非壓縮系統中的每一項都拆開，從頭解釋每個部分的物理意義與數學角色嗎？ 什麼是 Navier-Stokes 方程式？

這個系統從何而來？ Navier-Stokes 方程式的推導。

去掉粘性就得到 Euler 方程式；它早了整整一個世紀，在紙面上看起來更簡單，但在某些方面從數學上甚至更難理解，因為你失去了擴散項的平滑化效果。 Euler 與 Navier-Stokes。

獎題。 Navier-Stokes 存在性與光滑性問題。

從這裡出發的建議路徑：