為什麼 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 理論透過拋物柱上的局部能量不等式處理這點,但要從中提取逐點正則性,仍是困難步驟。
為什麼 3D 特別
對於標準設定下的 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 與 Navier-Stokes.