黏滯流體、雷諾數與白努利方程：從 Navier–Stokes 到理想極限

從理想到真實：黏滯性如何改寫流體圖像

當我們把一根吸管插進飲料、看見機翼把飛機托起時，背後其實藏著兩個對立的世界觀：一個是「能量守恆、滑順無摩擦」的理想流體，另一個是「層與層之間互相拖曳」的黏滯流體。普物課本常以白努利方程一筆帶過前者，卻很少把後者的數學機制攤開來。本篇要做的，正是把這兩個世界用同一組統御方程連起來，並說明「雷諾數」如何成為兩者之間的裁判。

我們的出發點是不可壓縮黏滯流體的 Navier–Stokes 方程式：

$$ \rho\left(\frac{\partial \mathbf{u}}{\partial t} + (\mathbf{u}\cdot\nabla)\mathbf{u}\right) = -\nabla p + \mu\nabla^2\mathbf{u} + \rho\mathbf{g} $$

搭配連續方程式（不可壓縮約束）$\nabla\cdot\mathbf{u}=0$。其中 $\mathbf{u}$ 是速度場、$p$ 是壓力、$\rho$ 是密度、$\mu$ 是動力黏滯係數、$\mathbf{g}$ 是體積力（如重力）。整篇文章的物理機制，幾乎都能從這條方程的不同項之間「誰主導」推導出來。

黏滯性的微觀來源與牛頓黏性定律

黏滯性的本質是動量的橫向擴散。想像流體被分成許多平行薄層，上層流得快、下層流得慢，分子的隨機熱運動會在層與層之間交換動量，使快層被拖慢、慢層被拖快。對於簡單剪切流，牛頓黏性定律寫成：

$$ \tau_{xy} = \mu\,\frac{\partial u_x}{\partial y} $$

$\tau_{xy}$ 是剪應力，$\partial u_x/\partial y$ 是速度梯度。這條線性關係定義了所謂「牛頓流體」（水、空氣、稀薄油），而血液、牙膏、玉米澱粉糊等則屬於非牛頓流體，黏度本身會隨剪切率變化。

把牛頓黏性定律代入動量守恆，並對應力張量取散度，就得到 Navier–Stokes 中的 $\mu\nabla^2\mathbf{u}$ 項。它的數學身分是一個擴散算子：黏滯性讓速度場像熱一樣往外抹平，其擴散係數正是運動黏度 $\nu = \mu/\rho$，量綱為 $\mathrm{m^2/s}$，與熱擴散係數完全同型。這也是為什麼黏滯流體與熱傳導在數學上是「同一家族」。

雷諾數：慣性與黏滯的擂台

把 Navier–Stokes 方程無因次化，是理解雷諾數的最快路徑。取特徵長度 $L$、特徵速度 $U$，令 $\mathbf{u}^*=\mathbf{u}/U$、$\nabla^*=L\nabla$、$t^*=tU/L$、$p^*=p/(\rho U^2)$，代入後（略去重力）整理得：

$$ \frac{\partial \mathbf{u}^*}{\partial t^*} + (\mathbf{u}^*\cdot\nabla^*)\mathbf{u}^* = -\nabla^* p^* + \frac{1}{Re}\nabla^{*2}\mathbf{u}^* $$

於是所有物理參數壓縮成一個無因次數——雷諾數：

$$ Re = \frac{\rho U L}{\mu} = \frac{UL}{\nu} $$

它的物理意義是「慣性項 $\rho(\mathbf{u}\cdot\nabla)\mathbf{u}$」與「黏滯項 $\mu\nabla^2\mathbf{u}$」的量級比。$Re\ll 1$（如微生物游泳、潤滑膜）時黏滯主導，流動可逆、層流穩定，方程退化為線性的 Stokes 方程；$Re\gg 1$（如河流、機翼）時慣性主導，黏滯只在貼壁的薄「邊界層」內重要，主流可近似為無黏。圓管流的臨界值約 $Re_c\approx 2300$，超過後層流逐漸失穩走向紊流。

雷諾數的另一個威力是動力相似律：兩個幾何相似的流場若 $Re$ 相同，無因次速度場就相同。這正是風洞縮尺模型能預測真實飛機氣動特性的理論根基。

白努利方程：無黏沿流線的能量帳本

白努利方程不是獨立公設，而是 Navier–Stokes 在特定極限下的積分結果。取穩態、不可壓縮、無黏（$\mu=0$）的歐拉方程，並利用向量恆等式

$$ (\mathbf{u}\cdot\nabla)\mathbf{u} = \nabla\!\left(\tfrac{1}{2}|\mathbf{u}|^2\right) - \mathbf{u}\times(\nabla\times\mathbf{u}) $$

把方程改寫後，沿著流線方向（與 $\mathbf{u}$ 平行）對渦度項取內積會消失，再沿流線積分即得：

$$ \frac{1}{2}\rho v^2 + p + \rho g h = \text{常數（沿同一流線）} $$

三項分別是動壓、靜壓與重力位能密度。關鍵限制在於它假設「無黏」與「沿流線」——一旦黏滯耗散登場（$Re$ 有限），總機械能會沿流向遞減，白努利「常數」就變成緩慢下降的量。換言之，白努利方程是 $Re\to\infty$ 理想極限的近似；現實中它只在黏滯耗散可忽略的核心流區成立，貼壁邊界層內並不適用。

定量小範例：文氏管測流量

考慮水平文氏管，入口截面積 $A_1 = 50\ \mathrm{cm^2}$、喉部 $A_2 = 20\ \mathrm{cm^2}$，水（$\rho=1000\ \mathrm{kg/m^3}$）流過，測得入口與喉部壓差 $\Delta p = p_1 - p_2 = 3000\ \mathrm{Pa}$。求體積流量 $Q$。

步驟一：連續方程 $A_1 v_1 = A_2 v_2 = Q$，故 $v_1 = Q/A_1$、$v_2 = Q/A_2$。

步驟二：水平管 $h_1=h_2$，白努利方程化為 $\tfrac12\rho v_1^2 + p_1 = \tfrac12\rho v_2^2 + p_2$，即

$$ \Delta p = \tfrac{1}{2}\rho\left(v_2^2 - v_1^2\right) = \tfrac{1}{2}\rho Q^2\left(\frac{1}{A_2^2} - \frac{1}{A_1^2}\right). $$

步驟三：解出

$$ Q = \sqrt{\frac{2\,\Delta p}{\rho\left(\dfrac{1}{A_2^2} - \dfrac{1}{A_1^2}\right)}}. $$

代入 $A_1 = 5\times10^{-3}\ \mathrm{m^2}$、$A_2 = 2\times10^{-3}\ \mathrm{m^2}$：

$$ \frac{1}{A_2^2} - \frac{1}{A_1^2} = 2.5\times10^5 - 4.0\times10^4 = 2.1\times10^5\ \mathrm{m^{-4}}. $$

$$ Q = \sqrt{\frac{2\times 3000}{1000\times 2.1\times10^5}} = \sqrt{2.857\times10^{-5}} \approx 5.35\times10^{-3}\ \mathrm{m^3/s}. $$

步驟四（自洽檢驗）：喉部速度 $v_2 = Q/A_2 \approx 2.67\ \mathrm{m/s}$。若喉部直徑約 $5\ \mathrm{cm}$、水 $\nu\approx10^{-6}\ \mathrm{m^2/s}$，則 $Re = v_2 d/\nu \approx 1.3\times10^5 \gg 1$，慣性遠大於黏滯，白努利的無黏假設在此核心流區成立，計算自洽。實務上仍需乘上排放係數 $C_d\approx0.98$ 修正薄邊界層造成的微小損失。

深入探討（研究所視角）

在更高的理論層次上，流體力學可被嵌入變分與場論架構。對理想流體，Clebsch 表述把速度場寫成 $\mathbf{u} = \nabla\phi + \alpha\nabla\beta$，使歐拉方程從一個拉格朗日密度 $\mathcal{L} = -\rho(\partial_t\phi + \alpha\,\partial_t\beta + \tfrac12|\mathbf{u}|^2 + e(\rho))$ 經最小作用量原理導出。此時白努利「常數」其實是時間平移對稱性下的守恆量，與 Noether 定理直接相連；而 Kelvin 環流定理 $\frac{d}{dt}\oint_C \mathbf{u}\cdot d\mathbf{l}=0$ 則對應到粒子重標籤對稱（particle relabelling symmetry）所保護的拓撲不變量——這正是渦度動力學與螺旋度（helicity）守恆的源頭。

哈密頓表述則揭示理想流體是一個無窮維、具非正則 Poisson 結構的可積系統，其 Casimir 不變量（如螺旋度 $\int \mathbf{u}\cdot\boldsymbol\omega\,dV$）對應渦線的拓撲纏繞，無法被任何理想演化改變。Arnold 在 1960 年代更指出，理想流體運動等價於體積保持微分同胚群 $\mathrm{SDiff}$ 上的測地線流，把流體不穩定性重新詮釋為這個無窮維流形上的負曲率——這是現代幾何流體力學的奠基洞見。

黏滯項的加入則打破上述守恆結構，引入耗散與不可逆性。當 $Re$ 趨於極大，能量沿尺度串級的圖像由 Kolmogorov 1941 理論刻畫：在慣性次區間，能量譜遵循 $E(k)\sim \varepsilon^{2/3}k^{-5/3}$，其中 $\varepsilon$ 是能量耗散率、$k$ 是波數，能量從大渦無耗散地傳到 Kolmogorov 微尺度 $\eta=(\nu^3/\varepsilon)^{1/4}$ 才被黏滯抹除。這套標度律至今仍是紊流研究的核心，而三維 Navier–Stokes 解的整體光滑性與唯一性，更是 Clay 數學研究所七大千禧年問題之一，尚未解決。

在跨主題連結上，Navier–Stokes 的 $\nu\nabla^2\mathbf{u}$ 與熱傳導方程同型，使「黏滯—熱擴散」可用 Prandtl 數 $Pr=\nu/\alpha$ 統一描述，把本文與「熱」主題群直接縫合；而流體的波動（聲波、表面重力波、毛細波）則是把可壓縮性或自由邊界重新放回方程後的線性化結果，與「波動」主題群相通。從一條 Navier–Stokes 方程出發，理想與真實、層流與紊流、守恆與耗散，最終都收束於同一張數學版圖。