Navier-Stokes 方程式的弱解、強解與光滑解:正則性的差距

千禧年獎問的是光滑解。對任意初始資料,我們能證明大域存在的只有弱解。這道落差就是整個問題。

什麼是弱解?

情況是這樣。千禧年獎懸賞一百萬美元,要求解決三維 Navier-Stokes 是否永遠都有能持續存在的光滑解,或者是否可能發生爆發。光滑的意思是速度場表現得完全良好:沒有突然跳躍、沒有無限大的速度、沒有數學失效的點。但在近一個世紀的努力中,任何人所證明的最佳存在性結果,只保證了某種較弱的東西。這些稱為弱解

那麼什麼是弱解?它不是近似。也不是「差不多正確」。它是方程式的精確解,只是遵守較寬鬆的規則。一般的(「古典」)解要求速度足夠光滑,使你能在每一個點計算它的變化率。弱解則跳過這個要求。你不是逐點檢查方程式,而是在空間區域上「平均地」檢查它們。

這裡有個類比。古典解就像一個學生,解每道考題都一步一步寫出完整過程。弱解則像一個學生,無法向你展示中間步驟,但對你可能提出的每個問題,其最終答案都可證明是正確的。你看不到他們怎麼算,但答案總是驗算無誤。

為什麼要接受這種東西?因為有時方程式對古典解而言太狂野了。流體可能發展出某些區域,在那裡速度變化劇烈到你根本無法計算變化率。數學失效了。弱解讓你能在古典解放棄的地方繼續前進。它們是在情況變得粗糙時,讓方程式仍然存活的安全網。

問題在於:弱解可能不唯一。從完全相同的流動出發,你可能得到多個弱解,而且沒有人能告訴你哪一個才是「真正答案」。這是個問題,因為物理上流體應該做一件特定的事,而不是好幾件事。而且弱解也可能不光滑。光滑性正是千禧年獎所要求的,也是沒有人能證明的。

Navier-Stokes 方程式的弱解,將逐點的 PDE 替換為分佈意義下的表述。考慮$\mathbb{R}^3 \times (0,T)$上的非壓縮系統:

$$\partial_t u + (u \cdot \nabla)u = -\nabla p + \nu \Delta u, \qquad \nabla \cdot u = 0.$$

一個無散度向量場$u \in L^2_{\mathrm{loc}}(\mathbb{R}^3 \times [0,T))$若對每個光滑、緊支撐、無散度的測試函數$\varphi$:

$$\int_0^T \!\int_{\mathbb{R}^3} \bigl[-u \cdot \partial_t \varphi - (u \otimes u) : \nabla \varphi - \nu \, u \cdot \Delta \varphi \bigr] \, dx \, dt = \int_{\mathbb{R}^3} u_0 \cdot \varphi(x,0) \, dx.$$

壓力從這個表述中消失,因為$\varphi$是無散度的。所有導數都透過分部積分從$u$移到$\varphi$上。重點是:$u$不需要在古典意義下可微。它只需要有足夠可積性,使這些積分收斂。

弱表述不是近似。逐點滿足 PDE 的古典解,也會對每個測試函數滿足弱表述(往另一個方向分部積分即可)。反之不成立:弱解未必光滑到能逐點滿足 PDE。

這點很重要,因為:(1)對$L^2$初始資料,弱解大域存在(Leray,1934),而三維中的古典解目前只知道局部時間存在;(2)一般而言,弱解是否唯一仍未知;(3)弱解是否總是光滑的問題,對適當初始資料而言等價於Clay 千禧年問題

Leray 與第一個存在性證明(1934)

1934 年,Jean Leray 做了一件至今仍定義著這個領域的事。在一篇 73 頁的論文中,他證明了三維 Navier-Stokes 方程式的弱解對所有時間都存在,並且可從任何合理的初始流動出發。任何一個。只要起始速度不是具有無限能量或在物理上毫無意義,Leray 就保證你會得到一個永遠持續的解。這是首次有人為三維方程式證明大域存在性結果;九十多年後,它仍然是我們擁有的最強無條件存在性定理。

他的策略很巧妙。真正的方程式太難直接求解,因為流體速度會回饋到自身(這就是非線性)。所以 Leray 稍微模糊化方程式,就像替影像加上一個很小的高斯濾波器。模糊化後的方程式溫順到足以求解。接著他把模糊程度調向零,並證明這些解不會分崩離析。它們會收斂到某個在弱意義下滿足原始、未模糊化方程式的東西。

但以下是 Leray 沒有證明的事。唯一性。他的方法產生至少一個弱解,但從同一個流動出發也可能有其他弱解。他無法排除這點。他也沒有證明光滑性。他的解具有有限能量,並滿足能量不等式:摩擦可以耗散能量,但能量不能憑空自發出現。僅此而已。沒有更多。

Leray 本人懷疑奇異性可能形成。他勾勒了可能的樣貌:流體朝某一點坍縮,越來越快,把所有能量集中到越來越小的區域,就像一個渦旋以無限速度收縮成一點。1996 年,Nečas、Růžička 與 Šverák 證明這種精確的自相似坍縮不可能發生。Leray 對潛在爆發形狀的猜測是錯的。至於爆發是否會以任何形式發生?沒有人知道。

1951 年,Eberhard Hopf 將 Leray 的構造推廣到有界容器中的流體(不只是整個無限空間),所得的類別後來稱為Leray-Hopf 弱解:滿足能量不等式的弱解。這是標準概念。當研究者沒有進一步限定而說「弱解」時,他們幾乎總是指這個。

還有一點。即使在 Leray-Hopf 弱解之中,也有一個更挑剔的子類,稱為適合弱解。這些解不只在大域上滿足能量不等式(總能量不增加),也在局部上滿足:能量不能暗中在流體的一角堆積,同時從另一角耗散。Caffarelli、Kohn 與 Nirenberg(CKN)在 1982 年證明其著名的部分正則性結果,正是針對這個較小類別。不要混淆兩者:CKN 適用於適合弱解,而不是所有 Leray-Hopf 解。

Leray 1934 年的論文建立了以下結果:對任意無散度的$u_0 \in L^2(\mathbb{R}^3)$,至少存在一個弱解$u$滿足$\mathbb{R}^3 \times (0,\infty)$上的 Navier-Stokes 方程式:

  • $u \in L^\infty(0,\infty; L^2(\mathbb{R}^3)) \cap L^2(0,\infty; \dot{H}^1(\mathbb{R}^3))$
  • 能量不等式:$\|u(t)\|_{L^2}^2 + 2\nu \int_0^t \|\nabla u(s)\|_{L^2}^2 \, ds \leq \|u_0\|_{L^2}^2$對幾乎每個$t > 0$

此構造透過磨光化進行。將非線性項$(u \cdot \nabla)u$替換為$(u_\varepsilon \cdot \nabla)u$,其中$u_\varepsilon = J_\varepsilon * u$是空間磨光化。正則化系統具有大域光滑解(磨光化消除了最糟的非線性交互作用)。Leray 為正則化解取得一致能量界,接著抽取弱收斂子序列。極限滿足弱表述與能量不等式。

Leray 沒有建立唯一性。緊性論證給出至少一個聚點的存在;不同子序列可能收斂到不同極限。三維中 Leray-Hopf 弱解的唯一性至今仍是開放問題。

Hopf(1951)將此構造改編到有界區域$\Omega \subset \mathbb{R}^3$並搭配 Dirichlet 邊界條件,使用 Galerkin 近似(投影到有限維子空間)而非磨光化。所得類別,也就是滿足能量不等式的弱解,帶有兩人的名字:Leray-Hopf 弱解

Caffarelli-Kohn-Nirenberg 定理(1982)關注一個更受限制的類別:適合弱解,它們另外滿足如下形式的局部能量不等式

$$\partial_t \left(\frac{|u|^2}{2}\right) + \nabla \cdot \left(\left(\frac{|u|^2}{2} + p\right)u\right) + \nu |\nabla u|^2 \leq \nu \Delta \left(\frac{|u|^2}{2}\right)$$

在分佈意義下。CKN 證明,對任何適合弱解,時空中奇異集的一維拋物 Hausdorff 測度為零。這表示奇異性若存在,會極其稀疏(它們具有零一維拋物 Hausdorff 測度)。但該定理並未說明奇異性是否真的發生,而且它只適用於適合弱解,不適用於所有 Leray-Hopf 解。

Leray 本人考慮過如下形式的自相似爆發可能性$u(x,t) = (T-t)^{-1/2} U(x / (T-t)^{1/2})$。Nečas、Růžička 與 Šverák(1996)證明,對$L^3(\mathbb{R}^3)$中的解,不存在這種自相似爆發;Tsai(1998)則在相應假設下排除了某些漸近自相似爆發情境。潛在奇異性的形狀若有,仍然未知。

強解與正則性

弱解大域存在。但它們可能不唯一,也可能不光滑。我們能做得更好嗎?

可以,但只能暫時。強解是升級版:它有足夠正則性,讓方程式幾乎處處成立,而不只是「平均地」成立。光滑解或古典解則更強:方程式在每一點都成立。對三維中的光滑初始資料,強解會在短時間內存在。多短?這取決於起始流動有多狂野。平靜、溫和的流動能得到較長的保證。劇烈、湍流般的起始條件?可能只有微秒。

而且沒有人能證明這些強解最終不會爆發。

1962 年,James Serrin 證明了一種類似升級規則的結果。它是這樣的:如果一個弱解碰巧保持足夠良好(不太大、不把能量集中到越來越小的區域),那麼它其實從頭到尾都是光滑的。你可以把它升級。而根據稱為弱-強唯一性的原理,它也是具有這些起始條件的唯一弱解。一個解,光滑且唯一,結案。但如果你無法驗證該解保持溫順?那就什麼也沒有。你卡住了。

這是一個條件性結果。若解不太狂野,則它表現得完全良好。整個困難就在於證明那個「若」。

二維,能量估計強到足以讓每個弱解自動通過 Serrin 的檢驗。完成。這就是為什麼二維情形已經解決。在三維中,這些估計距離所需條件只差一點點,而補上這個缺口就是整個問題的核心。

研究者也找到了其他條件式檢驗,每一個都是不同的進攻角度:「只要證明解具備這一個特定性質,我就免費給你光滑性。」若能無條件證明其中任何一個,就會解決千禧年問題。沒有人辦到。關於人們嘗試過的不同證明策略,有一整頁專門介紹。

Navier-Stokes 方程式的強解,是指具有足夠正則性,使得 PDE 逐點(幾乎處處)成立,且非線性項$(u \cdot \nabla)u$可作為函數良好定義,而不只是分佈。通常這表示$u \in L^\infty(0,T; H^1(\mathbb{R}^3)) \cap L^2(0,T; H^2(\mathbb{R}^3))$。一個相關但不同的框架是 Fujita-Kato(1964),它在臨界空間中構造局部 mild 解。

對足夠正則的 Sobolev 初始資料,強解的局部存在性已相當成熟。在臨界空間中,Fujita-Kato(1964)框架透過不動點論證構造局部 mild 解:

$$u(t) = e^{\nu t \Delta} u_0 - \int_0^t e^{\nu(t-s)\Delta} \mathbb{P} \nabla \cdot (u \otimes u)(s) \, ds,$$

其中$\mathbb{P}$是投影到無散度向量場上的 Leray 投影。這個積分方程式在適當的函數空間中,依壓縮映射原理有唯一的局部解。對臨界空間中的小資料($L^3$,$\dot{H}^{1/2}$,$BMO^{-1}$),解是大域的。

問題在於大資料強解是否能對所有時間持續存在。關鍵的條件式結果是 Serrin(1962):若一個 Leray-Hopf 弱解滿足$u \in L^p(0,T; L^q(\mathbb{R}^3))$且

$$\frac{2}{p} + \frac{3}{q} \leq 1, \qquad q > 3,$$

則$u$在$(0,T] \times \mathbb{R}^3$上是光滑的,並且是具有給定初始資料的唯一 Leray-Hopf 弱解。這些稱為Prodi-Serrin 條件(Prodi 1959 證明了一個相關結果)。

端點$q = 3$($p = \infty$)由 Escauriaza、Seregin 與 Šverák(2003)解決:若$u \in L^\infty(0,T; L^3(\mathbb{R}^3))$,則$u$不會在時間$T$爆發。這是 Prodi-Serrin 尺度的端點,也是目前已知最尖銳的延拓判準之一。

二維中,能量界$u \in L^\infty_t L^2_x \cap L^2_t H^1_x$結合 Ladyzhenskaya 不等式$\|f\|_{L^4}^2 \leq C\|f\|_{L^2}\|\nabla f\|_{L^2}$(二維特有)給出$u \in L^4_t L^4_x$,它滿足 Serrin 條件$2/4 + 2/4 = 1$(在二維版本中,$2/p + 2/q \leq 1$)。在三維中,能量界透過 Sobolev 嵌入給出$u \in L^{10/3}_t L^{10/3}_x$,它滿足$2/(10/3) + 3/(10/3) = 3/2 > 1$。Serrin 條件失敗的差距,正好對應於超臨界性的缺口。關於這個結構性障礙,參見為什麼 Navier-Stokes 如此困難以了解更多。

光滑解與千禧年問題

光滑解是黃金標準。速度場在所有地方、所有時間都表現完好。沒有突兀跳躍。沒有無限速度。無論你放大到多細,解都一直保持良好。

Clay 千禧年獎問題由 Charles Fefferman 於 2000 年表述,提出了一個可以寫在索引卡上的問題。從任意光滑且物理上合理、充滿三維空間的速度場開始。Navier-Stokes 方程式是否總會產生一個永遠存在的光滑解,還是你能找到某個初始流動,使得解最終爆發?

任一答案都值一百萬美元。

目前的狀況是這樣。沒有人證明三維中光滑解總是大域存在,也沒有人構造出爆發。自 Leray 1934 年的論文以來,我們一直卡在中間;這是整個數學中最困難的未解問題之一,已經超過九十年,而我們仍不知道答案會落在哪一邊。

短時間呢?沒問題。對光滑初始資料,方程式確實會在某段時間內產生光滑解。流體開始運動,數學運作正常,一切乾淨俐落。但之後會發生什麼?解會永遠保持光滑,還是會碰到某個點,速度一路飆到無窮大?

如果它確實保持光滑,就會發生一件好事。那個光滑解也會自動滿足弱解的放寬規則,因此它是弱解。而由弱-強唯一性可知,具有相同初始條件的其他弱解不可能存在。因此若有人證明大域光滑性,整個層級就會塌縮:弱解、強解、光滑解全都變成同一件事,也就是一個在所有時間都表現完好、單一且唯一的解。這正是這個問題如此吸引人也如此困難的原因。我們能證明存在的東西(弱解)與我們想要的東西(光滑解)之間的落差,正是一百萬美元問題的內容。

Clay 千禧年問題(Fefferman 2000)問的是:對任意無散度的$u_0 \in C^\infty(\mathbb{R}^3)$滿足$|\partial^\alpha u_0(x)| \leq C_{\alpha K}(1 + |x|)^{-K}$對所有$\alpha, K$,是否存在$u \in C^\infty(\mathbb{R}^3 \times [0,\infty))$以及$p \in C^\infty(\mathbb{R}^3 \times [0,\infty))$滿足 Navier-Stokes 方程式,且對所有 $t \geq 0$ 都有 $\int_{\mathbb{R}^3} |u(x,t)|^2 \, dx < C$?

另一種選項(同樣有獎):找到$u_0$屬於上述類別,使得這樣的光滑解不存在。

已知結果:

  • 局部存在性: 對$u_0 \in H^s(\mathbb{R}^3)$,$s \geq 1/2$,在$[0,T^*)$上存在一個唯一的局部 mild/強解,其中$T^* > 0$取決於$\|u_0\|_{H^s}$,且它對每個正時間$t > 0$都是光滑的。若 $T^* < \infty$,則當 $t \to T^*$ 時有 $\|u(t)\|_{H^s} \to \infty$。
  • 小資料大域存在性: 若$\|u_0\|_{\dot{H}^{1/2}}$(或$\|u_0\|_{L^3}$,或$\|u_0\|_{BMO^{-1}}$) 小於某個普適常數,則解是大域且光滑的。主要參考文獻:$\dot{H}^{1/2}$(Fujita-Kato 1964)、$BMO^{-1}$(Koch-Tataru 2001),以及 $L^3$ 中的對應結果(Kato 1984)。
  • 弱-強唯一性: 若強解存在於 $[0,T]$,則每個具有相同初始資料的 Leray-Hopf 弱解在 $[0,T]$ 上都與它一致。這由 Serrin(1962)建立,並由後續工作加以精煉。這意味著證明正則性也會解決 Leray-Hopf 類中的唯一性。

此層級向上塌縮:光滑 $\Rightarrow$ 強 $\Rightarrow$ 弱,而弱-強唯一性意味著若光滑解存在,它就是唯一的 Leray-Hopf 弱解。因此 千禧年問題 等價於問:所有具有光滑初始資料的 Leray-Hopf 弱解本身是否都是光滑的?「弱解存在」(Leray 1934)與「光滑解存在」(未解)之間的缺口,正是這道懸賞問題。

為什麼這個區分重要

如果弱解存在並描述了流體,為什麼還要在意光滑性?

有三個理由。

第一,唯一性。物理要求一個答案。給我流體的初始狀態,我應該能精確地告訴你它接下來會怎麼做,而不是「這裡有幾種可能,隨你挑」。但弱解並不保證這一點。多個弱解可能從同一個初始流動中產生,卻無法判定真實流體會遵循哪一個。方程式會變成菜單,而不是食譜。那不是物理。

第二,數值可靠性。許多重要的流體模擬都基於 Navier-Stokes 方程式或與其密切相關的模型:天氣預報、空氣動力學、動脈中的血流等等。把網格加細,模擬應該收斂到真正的答案。沒有光滑性與唯一性的保證呢?沒有定理說這在每個三維情境中都確實會發生。模擬能用;但我們無法完全解釋為什麼。

第三,極端物理。如果奇異性可以形成,那就是自然在向我們傳遞訊息。Navier-Stokes 方程式,這個我們描述流體運動的最佳模型,將會內建一個到期日:在某個極端尺度下,模型本身停止適用,而方程式是在告訴我們:「你需要新的物理。」

這不是技術細節。它是貫穿一切的斷層線。一邊是存在性(Leray,1934,已完成)。另一邊是光滑性(未解,一百萬美元)。為什麼跨越如此困難? 三維中的能量估計剛好差了所需的一點點;九十年的努力、數百篇論文、整個職涯投入其中,仍沒有人補上這個缺口。

目前正在追求的每一種 證明策略 都是在嘗試彌合這道鴻溝。證明弱解是光滑的。或者證明它們不是。兩三個字而已。

三維非壓縮 Navier-Stokes 方程式的解層級如下:

$$\text{smooth} \subset \text{strong} \subset \text{Leray-Hopf weak} \subset \text{distributional weak}$$

各層級已知情況:

類別大域存在性唯一性正則性
分佈意義弱解是(無能量控制)可以非常粗糙
Leray-Hopf 弱解是(Leray 1934,Hopf 1951)未解$L^\infty_t L^2_x \cap L^2_t \dot{H}^1_x$
適合弱解未解CKN:奇異集合具有 $\mathcal{P}^1 = 0$
強解/溫和解未解(小資料時為是)是(在存在區間內)$C^\infty$ 若存在
光滑($C^\infty$)未解(= Clay 千禧年問題)是(弱-強唯一性)$C^\infty$ 依定義

Clay 千禧年問題位於大域 Leray-Hopf 弱解存在性與大域光滑正則性之間的缺口。核心困難在於:能量不等式給出 $u \in L^\infty_t L^2_x \cap L^2_t \dot{H}^1_x$,這比臨界尺度變換低了半個導數 $\dot{H}^{1/2}$。彌合這個超臨界缺口,等價於證明大域正則性。

弱-強唯一性意味著,若光滑解大域存在,則此層級會塌縮。更精確地說:若對給定的 $u_0 \in C^\infty$ 有一個光滑解 $u$ 存在於 $[0,T]$,則每個具有相同資料的 Leray-Hopf 弱解都等於 $u$ 在 $[0,T]$ 上。因此大域光滑性 $\Rightarrow$ Leray-Hopf 類中的大域唯一性。

反過來說,若對光滑可容許初始資料的 Leray-Hopf 弱解出現非唯一性,則會意味著該資料類的大域光滑解不可能持續存在(因為光滑解會迫使唯一性)。近年關於凸積分的工作(建立在 De Lellis-Székelyhidi 對 Euler 方程式的工作之上,並由 Buckmaster-Vicol 2019 延伸,用以構造低於 Leray-Hopf 正則性的 Navier-Stokes 方程式非唯一弱解)顯示,分佈意義弱解可以高度非唯一。這種非唯一性是否延伸到 Leray-Hopf 類,是一個重大未解問題,並且對 Clay 千禧年問題有直接影響。

關於目前的 證明策略 關於如何攻克這個缺口,以及它為何在結構上如此難以處理,請參見 為什麼 Navier-Stokes 如此困難.