為什麼 Navier-Stokes 問題困難：非線性、壓力與湍流

許多人最早在物理中接觸到的方程式是線性的：輸入加倍，反應也加倍。Navier-Stokes 方程式不是這樣。

Navier-Stokes 方程式？ 非線性。流體的速度會影響它自身的變化率，這表示流體會推動自身。想像一下，要預測一群人會往哪裡走，而每個人的移動都取決於周圍所有人在做什麼，且那些人在做什麼又取決於他們周圍所有人，如此永遠向外盤旋擴散。這就是你正在面對的情況。

罪魁禍首是自交互作用項 $(u \cdot \nabla)u$。它會產生回饋迴路，使微小擾動放大成巨大擾動，也正因如此，流體湍流才會如此極端複雜（更多內容見子問題）。

對流非線性項 $(u \cdot \nabla)u$ 是根本障礙。在渦度形式 $\omega = \nabla \times u$ 中，方程式變為

$$\partial_t \omega + (u \cdot \nabla)\omega = (\omega \cdot \nabla)u + \nu \Delta \omega$$

渦旋拉伸項 $(\omega \cdot \nabla)u$ 沒有固定符號。它可以無界地放大渦度。在 2D 中此項消失（因為 $\omega$ 是垂直於流動的純量），這就是為什麼2D 大域正則性已知（Ladyzhenskaya，1969）。在 3D 中，渦旋拉伸是有限時間爆發的主要候選機制。

關鍵在於：非線性項對 二次 依賴於 $u$。$H^1$ 能量估計給出 $\|\nabla u\|_{L^2}$，但要控制 $(u \cdot \nabla)u$ 在 $L^2$ 中，通常需要更強的資訊，例如 $u \in L^\infty$ 連同 $\nabla u \in L^2$，或等價的臨界尺度控制；僅靠能量類並不能提供這一點。

Navier-Stokes 方程式具有一種尺度變換對稱性。把一個解放大觀察，按正確比例讓一切變得更小且更快，你就會得到另一個完全有效的解。這種對稱性在數學上很自然，但在分析上很危險。

為什麼？我們唯一能可靠控制的量是流體的總能量，而它位於完全錯誤的尺度上，只告訴我們整體圖像，卻對爆發實際形成的微觀尺度上正在發生的事毫無說明。

這就像監測一座城市的總用電量來偵測一粒火花。有用嗎？當然。有足夠細緻嗎？差得遠。這個落差就是整個問題。

在自然的尺度變換 $u_\lambda(x,t) = \lambda u(\lambda x, \lambda^2 t)$ 下，臨界 Sobolev 空間是 $\dot{H}^{1/2}(\mathbb{R}^3)$（等價地 $L^3$）。若某個受控範數在放大觀察時變小，則相對於那種控制，方程式就是超臨界的：這個界低於看見集中現象所需的尺度。

能量不等式給出對 $u \in L^\infty_t L^2_x \cap L^2_t \dot{H}^1_x$ 的控制。就能量層級而言，3D 問題是超臨界的：

$$\|u_\lambda\|_{L^2} = \lambda^{-1/2} \|u\|_{L^2}, \quad \|u_\lambda\|_{L^2_t \dot{H}^1_x} = \lambda^{-1/2} \|u\|_{L^2_t \dot{H}^1_x}$$

當 $\lambda \to \infty$（放大觀察）時，這些界相對於不變量 $\|u\|_{L^3}$ 與 $\|u\|_{\dot{H}^{1/2}}$ 會變弱。非線性項與粘性項在尺度上是一起變換的；問題在於已知的先驗估計位於臨界尺度之下。從能量類推進到臨界範數，是核心困難。

觀察一條河流。流體運動變得混沌。湍流。大渦分裂成較小的渦，較小的渦又分裂成更小的渦，一路級串到微觀尺度，最後由粘性把一切抹平。

Kolmogorov 在 1941 年描述了這種能量級串，而 Navier-Stokes 方程式漂亮地捕捉了它。但真正讓人夜不能寐的是：如果能量集中到越來越小的區域，其速度快過粘性所能耗散的速度，會怎樣？那就是爆發。

它真的會發生嗎？還是粘性總會勝出？這就是懸而未決的問題，句點。關於從 Reynolds 數到這幅小尺度圖像的物理橋樑，請見Reynolds 數、湍流，以及為何小尺度很重要。

Kolmogorov 的 K41 理論預測能量譜 $E(k) \sim \varepsilon^{2/3} k^{-5/3}$ 在慣性範圍 $k_f \ll k \ll k_\eta$ 中，其中 $k_\eta \sim (\varepsilon/\nu^3)^{1/4}$ 是 Kolmogorov 耗散波數。能量通量在各尺度之間保持恆定。這就是乾淨的圖像。

正則性問題問的是更陰暗的事：這個級聯能否退化？耗散尺度能否$k_\eta^{-1}$在有限時間內縮小到零，且$\|\nabla u\|_{L^2} \to \infty$ 同時總能量仍保持有限？

Onsager 的耗散異常 猜想（1949）在某種意義上說是的：在粘性趨於零的極限中$\nu \to 0$，能量耗散仍然存在，而且弱解 可以使 Euler 方程式耗散能量。這在 Hölder 指數低於$1/3$ 的情況下已獲確認（Isett，2018；Buckmaster 等，2018）。但這對 Navier-Stokes 正則性意味著什麼，仍然真正不清楚。關於層級性的直觀理解，請見Reynolds 數、湍流，以及為什麼小尺度很重要。

壓力在 Navier-Stokes 方程式中很奇特。它根本不是獨立變數；速度透過一個單一約束完全決定它：流體是非壓縮的，所以不能被擠壓。

這使得壓力成為非局部。在非壓縮模型中，改變某一區域的速度會影響全域的壓力場，因為壓力由該時刻整個速度場決定。

對分析而言，這是毀滅性的。你無法只研究單一點發生的事，而不同時考慮整個流體。局部推理？別想了。

非壓縮約束$\nabla \cdot u = 0$ 透過 Poisson 方程決定壓力

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

所以$p = (-\Delta)^{-1} \partial_i \partial_j (u_i u_j)$，其中涉及 Riesz 變換（奇異積分算子）。壓力是速度的非局部函數，而這種非局部性正是逐點估計或空間局部估計的核心障礙。

標準的最大值原理論證在這裡失效：即使粘性項$\nu \Delta u$ 是耗散性的，壓力梯度$-\nabla p$ 仍可能把遠處區域的能量集中起來。Caffarelli-Kohn-Nirenberg 理論透過拋物柱上的局部能量不等式處理這點，但要從中提取逐點正則性，仍是困難步驟。

對於標準設定下的 2D 非壓縮 Navier-Stokes 方程式，已知存在大域光滑解；這在包括 Ladyzhenskaya（1969）在內的經典工作中已建立。見為什麼 2D 較容易。

三維呢？一切都崩解，而原因歸結於一個機制：渦旋伸展。在 2D 中，渦旋可以旋轉和合併，但不能伸展。在 3D 中，流體可以抓住渦管，把它們拉得越來越細、越來越細、越來越細，可能把每一點能量都集中到一條無限細的絲狀結構中。

這種伸展能否在有限時間內失控到無限大，還是粘性總會介入？這就是那個百萬美元問題。字面意義上。

這個二分非常尖銳。

2D： 渦量$\omega$ 是滿足$\partial_t \omega + u \cdot \nabla \omega = \nu \Delta \omega$ 的純量。最大值原理給出$\|\omega(t)\|_{L^\infty} \leq \|\omega_0\|_{L^\infty}$，BKM 蘊含大域正則性，而渦旋伸展項$(\omega \cdot \nabla)u$ 恆等於零。

3D： 渦量$\omega \in \mathbb{R}^3$ 滿足$\partial_t \omega + (u \cdot \nabla)\omega = (\omega \cdot \nabla)u + \nu \Delta \omega$，其中伸展項$(\omega \cdot \nabla)u$ 可以超線性地放大$|\omega|$（形式上$\sim |\omega|^2$ 經由 Biot-Savart）。沒有可用的最大值原理。

渦量平方積分$\|\omega\|_{L^2}^2$ 滿足

$$\frac{d}{dt}\|\omega\|_{L^2}^2 \leq C\|\omega\|_{L^2}^2 \|\nabla u\|_{L^\infty} - 2\nu \|\nabla \omega\|_{L^2}^2$$

控制$\|\nabla u\|_{L^\infty}$ 需要$\omega \in L^\infty$，而這又需要控制$\|\nabla u\|_{L^\infty}$。循環。現有技術尚未打破這個迴圈。

本文是問題的一部分。

這些障礙使數學家轉而將問題分解為子問題，並為每個子問題發展專門方法。

關於為什麼粘性項有幫助但仍不夠的背景，請見Euler vs. Navier-Stokes。

本文是問題.

正則性問題可分解為較易處理的部分；見子問題。為克服這些障礙而建立的分析工具，見方法。為何粘性項有所幫助，卻仍無法彌合缺口？Euler 與 Navier-Stokes.