Navier-Stokes 方程式是什麼:流體方程的意義與數學形式
從簡單直觀到數學形式,清楚介紹描述流體運動的偏微分方程式
Navier-Stokes 方程式是什麼
Navier-Stokes 方程式 是一組偏微分方程式,用來描述水和空氣等粘性流體的運動。
它們用於描述管中的水、飛機機翼周圍的空氣、動脈中的血液,以及無數其他流動。
從高層次來看,它們說的是:流體會因壓力、粘性以及作用於其上的任何力而改變其運動。壓力推動流體四處流動,粘性會抹平運動中的尖銳差異,而重力等外力可以驅動流動。
這些方程式不只是物理口號。它們是流體力學、工程學與計算模擬中很大一部分的工作語言。
Navier-Stokes 方程式是粘性 Newtonian 流體的動量平衡方程式。在非壓縮設定中,它們透過一個非線性 PDE 系統耦合速度場 $u(x,t)$ 與壓力 $p(x,t)$。
嚴格來說,Navier-Stokes 是一個方程式系統,而不是單一方程式。它們將動量守恆與本構律一起建模;該本構律指出粘性應力與對稱速度梯度成正比。非壓縮性則加入局部體積保持的約束。
對於 Clay 千禧年問題以及本網站大部分內容,相關設定是 三維非壓縮 系統。
用簡單語言理解 Navier-Stokes 方程式
在最簡單且常見的形式中,方程式如下:
$$\partial_t u + (u \cdot \nabla)u = -\nabla p + \nu \Delta u + f$$
$$\nabla \cdot u = 0$$
其中:
- $u$ 是流體的速度
- $p$ 是壓力
- $\nu$ 是運動粘度
- $f$ 是任何外力,例如重力
左邊描述速度如何隨時間改變,以及流體如何輸送自身的運動。右邊包含推動並平滑流動的各種力。
在 $\mathbb{R}^3$ 上的非壓縮 Navier-Stokes 系統為
$$\partial_t u + (u \cdot \nabla)u = -\nabla p + \nu \Delta u + f,$$
$$\nabla \cdot u = 0, \qquad u(x,0)=u_0(x).$$
各項有標準解釋:
- $\partial_t u$:局部時間變化
- $(u \cdot \nabla)u$:非線性平流,意即流體輸送自身的速度
- $-\nabla p$:壓力力
- $\nu \Delta u$:粘性擴散
- $f$:外部強迫
- $\nabla \cdot u = 0$:非壓縮約束
這是本網站關於千禧年問題、困難之處與證明策略的各頁所使用的形式。
Navier-Stokes 方程式從何而來
這些方程式來自一個簡單想法:將 Newton 第二定律應用於一小團流體。該流體微團的質量乘以加速度,必須等於作用於其上的總力。
對於粘性流體,這些力主要來自壓力與內部摩擦。當你在流體中每一點仔細寫下這個平衡,就會得到 Navier-Stokes 方程式。
因此這些方程式並非任意而來。它們是連續介質力學版本的力等於質量乘以加速度。完整的逐步推導請見Navier-Stokes 方程式是如何推導出來的。
推導從連續介質的動量守恆開始。先在物質體積上寫出線動量平衡,再將該恆等式局部化以得到 PDE。
對於 Newtonian 非壓縮流體,Cauchy 應力張量具有以下形式
$$T = -pI + 2\mu D(u), \qquad D(u)=\frac{1}{2}(\nabla u + \nabla u^T),$$
其中 $\mu$ 是動力粘度。將此本構律代入動量方程式並除以密度,即得到熟悉的非壓縮系統,其運動粘度為 $\nu = \mu/\rho$。
在標準的常密度非壓縮設定中,連續方程式化為 $\nabla \cdot u = 0$。關於從動量平衡與 Newtonian 本構律出發的完整推導,請見 推導。
為什麼 Navier-Stokes 方程式很困難
困難之處在於非線性項 $(u \cdot \nabla)u$。流體不只是回應外力;它也會推動自身四處運動。這種回饋使湍流與看似混沌的運動成為可能。
在二維空間中,方程式的表現要好得多。在三維中,我們仍不知道每個光滑初始流動是否都會永遠保持光滑。
這就是為什麼這些方程式的知名度遠超工程領域:它們直接導向 Navier-Stokes 千禧年問題。
它們的用途
Navier-Stokes 方程式每天都被用於科學與工程。典型應用包括:
- 機翼與車輛周圍的氣流
- 天氣與氣候模型
- 海洋環流
- 工業流體輸送
- 血流與其他生物輸送問題
在實務上,人們通常以數值方式求解這些方程式的近似,並常加入額外的建模假設。這種實務上的成功,正是剩餘數學問題如此引人注目的原因之一。
應用工作通常在特定狀態下使用 Navier-Stokes 方程式或相關模型的數值近似:非壓縮流、壓縮流、湍流閉合、邊界層近似,以及降階模型。
直接數值模擬、大渦模擬與 Reynolds 平均閉合都可追溯到同一個連續介質 PDE 框架,但它們並不能消除三維中的基礎正則性問題。
實務有效性與理論不完整之間的這種分裂,正是此主題如此引人入勝的部分原因。
接下來讀什麼
如果你的主要問題是這個問題是否已被解決,請先讀 Navier-Stokes 問題解決了嗎?。
如果你想了解其廣泛的數學意義,請繼續讀 千禧年問題。
如果你想理解 Navier-Stokes 方程式與無粘性 Euler 方程式有何不同,以及粘性為何重要,請讀 Euler 與 Navier-Stokes。
如果你想了解湍流與小尺度背後的物理直覺,請讀 Reynolds 數、湍流,以及小尺度為何重要。
如果你想了解主要障礙,請前往 為何困難。
本網站上自然的下一步:
- Navier-Stokes 問題解決了嗎? 了解狀態問題,以及弱存在性與大域光滑理論之間的區別
- 千禧年問題 了解 Clay 的精確表述
- Euler 與 Navier-Stokes 了解粘性的作用,以及無粘性與粘性系統之間的關鍵差異
- 非壓縮與壓縮 Navier-Stokes 方程式 了解常密度與變密度表述之間的差異
- Reynolds 數、湍流,以及小尺度為何重要 了解多尺度流動背後的狀態層級直覺
- 為何困難 了解尺度變換、超臨界性與爆發情境