Navier-Stokes 精確解：Poiseuille 流、Couette 流與 Stokes 擴散

Navier-Stokes 方程式以非線性著稱。那麼，怎麼可能有人能精確求解它們？

對稱性。這就是全部訣竅。

當流動的幾何足夠簡單（一直管、兩塊平板、無限平面）時，速度只能指向一個方向，並沿一個或兩個座標變化。在許多這類對稱設定中，非線性項要麼消失，要麼大幅簡化，使方程式退化為線性 PDE。在穩態情形下，通常只剩下一個可以用紙筆求解的 ODE。

直覺是這樣的。Navier-Stokes 方程式描述所有可能的流體運動。但如果你把流體限制在非常有秩序的情境中，沒有旋渦、沒有混沌，一切都朝同一方向前進，那麼方程式的大部分複雜性就變得無關緊要。Navier-Stokes 的困難之處在於流體推動自身運動的回饋迴路。在這些對稱流中，沒有什麼可供它推擠。非線性自平流項就直接消失了。

這些精確解並非奇聞軼事。它們是流體力學教育的基礎、數值程式的基準，也是理解真實流動何時以及如何失控的起點。

Navier-Stokes 方程式的一個精確解是速度場 $u$ 和壓力 $p$ 滿足

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

以閉形式給出，不依賴數值近似。

關鍵機制是非線性平流項 $(u \cdot \nabla)u$ 的消失或簡化。對於 平行流（即 $u = (u(y,t),\, 0,\, 0)$ 在直角座標中，或 $u = (0,\, 0,\, u(r,t))$ 在圓柱座標中）速度場自動無散度，而平流項恆等消失，因為速度在流向方向上沒有分量，也不沿流向方向變化。動量方程因此化為線性 PDE；若為穩態流，則化為 ODE。

當流動也是 穩態（$\partial_t u = 0$）時，剩下的是

$$\nu \Delta u = \nabla p,$$

一個 Poisson 方程，其解就是粘性流的經典精確解：Poiseuille 流、Couette 流以及它們的相關流動。

對於非穩態平行流（$\partial_t u \neq 0$ 但平流仍然消失），方程式變成擴散方程 $\partial_t u = \nu\, \partial_{yy} u$，由此得到 Stokes 的第一與第二問題。

Poiseuille 流（也稱為 Hagen-Poiseuille 流）是 Navier-Stokes 方程式最重要的精確解，也是多數工程師最先接觸到的那一個。

想像水穩定地流過一根又長又直的管子。兩端之間的恆定壓差驅動流動。管壁不動，因此接觸管壁的流體被固定在零速度（這就是無滑移條件）。離管壁越遠，流體速度越快。中心處移動最快。

速度剖面？是一條 拋物線。在管壁為零。在中心最大。中間是一條光滑曲線。把管子切開看橫截面：它像一個倒扣的碗。

關鍵的定量事實是：總流量隨管半徑的 四次方 成比例。四次方。半徑加倍，流量不是加倍，甚至也不是變成四倍；而是變成十六倍。這就是 Hagen-Poiseuille 定律，也解釋了為什麼即使動脈出現一點點狹窄，也可能阻斷血液供應。

假設： Poiseuille 流假設流體是非壓縮且 Newtonian（粘度為常數），流動為穩態且充分發展（不是仍在入口處加速），並且流動為層流。光滑。有序。實際上，管流在 Reynolds 數約為 2,300 時會轉變為湍流；這是一個經驗觀察，尚無人能從底層理論推導出來。

考慮半徑為 $R$ 的圓管中，穩態、充分發展、軸對稱的流動，由均勻壓力梯度 $dp/dx < 0$ in the axial direction $x$. In cylindrical coordinates $(r, \theta, x)$ the velocity field has the form $u = (0,\, 0,\, u(r))$.

非壓縮條件 $\nabla \cdot u = 0$ 自動滿足。平流項消失，因為 $u$ 不依賴於 $x$。軸向動量方程化為

$$0 = -\frac{dp}{dx} + \mu\left(\frac{d^2 u}{dr^2} + \frac{1}{r}\frac{du}{dr}\right),$$

或等價地

$$\frac{1}{r}\frac{d}{dr}\left(r\frac{du}{dr}\right) = \frac{1}{\mu}\frac{dp}{dx}.$$

由於 $dp/dx$ 為常數，這是一個關於 $r$ 的 ODE。配合邊界條件 $u(R) = 0$（無滑移）和 $du/dr|_{r=0} = 0$（對稱性）積分兩次，得到 拋物線速度剖面：

$$u(r) = \frac{1}{4\mu}\left(-\frac{dp}{dx}\right)(R^2 - r^2).$$

中心線處的最大速度為 $u_{\max} = \frac{R^2}{4\mu}\left(-\frac{dp}{dx}\right)$，平均速度為 $\bar{u} = u_{\max}/2$。

對橫截面積分得到 Hagen-Poiseuille 定律，給出體積流率：

$$Q = \frac{\pi R^4 \Delta p}{8 \mu L},$$

其中 $\Delta p > 0$ 是長度為 $L$ 的管段上的壓降。

這個 $R^4$ 依賴性是其決定性特徵：半徑的微小變化會造成流量的大幅變化。在血流動力學中，這就是血流對動脈狹窄如此敏感的原因。

適用性： 此解適用於非壓縮、Newtonian、穩態、充分發展的層流。實驗上，在 $\text{Re} = \bar{u} D / \nu \approx 2{,}300$ 時會發生向湍流的轉變，其中 $D = 2R$ 是管徑。這個閾值是一項經驗觀察；沒有定理能從 Navier-Stokes 方程式預測它。關於 Reynolds 數與層流—湍流轉變的討論，請見 Reynolds 數、湍流，以及為什麼小尺度很重要.

Couette 流 是流體被夾在兩塊平行平板之間、其中一塊平板運動而另一塊保持靜止時的精確解。它是流體力學中最簡單的精確解之一。

想像一疊紙牌平放在桌上。把最上面的牌往側邊拖動，下面的牌也會跟著移動，只是每張都比上面那張移動得少一些。就是這樣。速度從底板的零線性變化到頂板的速度，除此之外沒有更多複雜之處。

一條直線速度剖面。底板：靜止。頂板：運動。中間的一切只是線性插值，不需要壓力梯度，不需要複雜設定，運動純粹由移動邊界拖曳流體而驅動。

當你還沿通道施加壓力梯度時，事情就更有趣了，因為這時是在同一個間隙中結合剪切驅動與壓力驅動的流動。速度剖面會變形成疊加在線性剖面上的拋物線，有時稱為 平面 Poiseuille-Couette 流，而且取決於壓力梯度相對於平板速度的強度，甚至可能在一側壁面附近出現回流。

把移動平板完全拿掉，讓兩側壁面都保持靜止，並讓壓力差完成所有工作。這就是 平面 Poiseuille 流，也就是管流的平板類比。拋物線形。在中間最快。兩側壁面為零。

考慮兩塊相距間隙 $h$ 的無限平行平板之間的穩態流。令 $y$ 為垂直於平板的座標，底板位於 $y = 0$，頂板位於 $y = h$。流動是平行的：$u = (u(y),\, 0,\, 0)$.

簡單 Couette 流。 頂板以速度 $U$ 沿 $x$ 方向運動；底板靜止；$dp/dx = 0$。動量方程式化簡為 $d^2u/dy^2 = 0$，得到

$$u(y) = U\frac{y}{h}.$$

這是 Navier-Stokes 方程式最簡單的非平凡精確解：一個完全由邊界條件驅動的線性速度剖面。

平面 Poiseuille 流。 兩塊平板都靜止；一個常數壓力梯度 $dp/dx < 0$ drives the flow. The momentum equation becomes

$$\mu \frac{d^2 u}{dy^2} = \frac{dp}{dx},$$

且 $u(0) = u(h) = 0$。解為

$$u(y) = \frac{1}{2\mu}\left(-\frac{dp}{dx}\right)y(h - y),$$

一個關於中心線 $y = h/2$ 對稱的拋物線剖面。

一般 Couette-Poiseuille 流。 將移動頂板與壓力梯度結合，得到疊加

$$u(y) = U\frac{y}{h} + \frac{1}{2\mu}\left(-\frac{dp}{dx}\right)y(h - y).$$

第一項是剪切驅動的貢獻；第二項是壓力驅動的貢獻。取決於 $dp/dx$ 相對於 $\mu U / h^2$ 的符號與大小，剖面可以是單調的、有內部極大值，甚至在一側壁面附近出現回流。

Poiseuille 流與 Couette 流都是穩態的。它們不隨時間改變。Stokes 問題 是最簡單的非穩態精確解，並揭示了一件美麗的事：粘性使動量在流體中擴散，就像熱在固體中擴散一樣。

Stokes 第一問題（也稱為 Rayleigh 問題）：想像一大片靜止的流體停留在一塊平板上方。在時間零點，平板突然開始以恆定速度側向滑動。緊貼平板的流體會立刻被拖著走，但較遠處的流體需要一段時間才會察覺。光滑的邊界層從平板向外成長，隨時間推移而變厚。

平板上方任意高度的速度取決於該高度與特徵擴散長度 $\sqrt{\nu t}$ 的比值，其中 $\nu$ 是粘性係數，而 $t$ 是經過的時間。流體越粘？運動向上傳播得越快。

Stokes 第二問題：相同的設定，但現在平板不是以恆定速度運動，而是作正弦式來回振盪。振盪只會穿透到流體中的有限距離。再往上，流體幾乎感覺不到。運動振幅隨高度呈指數衰減，形成薄的振盪邊界層。這就是振盪邊界層背後的機制：振盪的平板使附近流體運動，但擾動會隨著離平板距離增加而指數衰減。

兩個 Stokes 問題都涉及位於 $y > 0$ 處平板上方的半無限流體（$y = 0$）。流動是平行的：$u = (u(y,t),\, 0,\, 0)$。平流項消失，控制方程式是一維擴散方程式：

$$\frac{\partial u}{\partial t} = \nu \frac{\partial^2 u}{\partial y^2}.$$

Stokes 第一問題（Rayleigh 問題）。 初始條件：$u(y,0) = 0$ 對於 $y > 0$。邊界條件：$u(0,t) = U$ 對於 $t > 0$。遠場：$u \to 0$ 當 $y \to \infty$。

相似變數 $\eta = y / (2\sqrt{\nu t})$ 將 PDE 化簡為 ODE。解為

$$u(y,t) = U\,\operatorname{erfc}\!\left(\frac{y}{2\sqrt{\nu t}}\right),$$

其中 $\operatorname{erfc}$ 是互補誤差函數。邊界層厚度按 $\delta \sim \sqrt{\nu t}$ 成長，這是擴散展布的標誌。

Stokes 第二問題。 平板振盪：$u(0,t) = U\cos(\omega t)$。尋找形如 $u(y,t) = \operatorname{Re}[\hat{u}(y)\, e^{i\omega t}]$ 的解，得到

$$\hat{u}(y) = U\exp\!\left(-(1+i)\frac{y}{\delta_s}\right), \qquad \delta_s = \sqrt{\frac{2\nu}{\omega}},$$

因此

$$u(y,t) = U\exp\!\left(-\frac{y}{\delta_s}\right)\cos\!\left(\omega t - \frac{y}{\delta_s}\right).$$

振幅隨穿透深度呈指數衰減 $\delta_s$，且相位線性滯後：這是一個橫波，其能量完全由粘性耗散。頻率越高，穿透越淺。

Poiseuille、Couette 與 Stokes 流受到最多關注。不過，它們絕不是唯一的精確解。遠遠不是。

• Taylor-Green 渦旋：二維中一種衰減的旋轉渦旋圖樣，具有真正的渦量結構。測試過 CFD 程式碼嗎？你很可能已經拿它跑過基準測試。這是大家最先想到的標準測試案例，而且幾十年來一直如此。
• Jeffery-Hamel 流：楔形通道中的流動，通道可以收斂或發散。它捕捉流體如何加速進入變窄的間隙，或在擴張的通道中減速。
• Hiemenz 駐點流：流體正面撞上一面平牆，在表面減速到零並向側方偏轉。風吹向建築物。噴流撞擊平板。

對稱性。這些解每一個都利用特定的幾何對稱性，使方程式變得可處理；它們在專門情境中很重要，但對日常入門流體力學而言，Poiseuille 與 Couette 仍然承擔了主要工作。

除了平行流之外，還有若干其他精確解族利用特定對稱性或自相似結構：

• Taylor-Green 渦旋。 在二維中，$u = (A\cos(ax)\sin(by),\, -B\sin(ax)\cos(by))$ 且 $aA = bB$（非壓縮）以及指數型時間衰減 $e^{-\nu(a^2+b^2)t}$。這是完整二維 Navier-Stokes 方程式的一個精確解，包括非線性項在內（結果它是一個梯度，並被吸收到壓力中）。它是 DNS 程式碼的標準驗證案例。
• Jeffery-Hamel 流。 穩態、純徑向流 $u_r(r,\theta)$ 位於半角為 $\alpha$ 的楔形區域中。流函數 $\psi = \nu f(\theta)$ 滿足一個關於 $\theta$ 的三階非線性 ODE。收斂與發散通道皆存在解，且在較高 Reynolds 數下具有豐富的分岔結構。
• Hiemenz 駐點流。 描述流動撞擊平面牆的二維相似解。流函數 $\psi = \sqrt{a\nu}\, x\, f(\eta)$，$\eta = y\sqrt{a/\nu}$，將 Navier-Stokes 化約為 Hiemenz ODE：$f''' + ff'' - f'^2 + 1 = 0$ 且 $f(0)=f'(0)=0$，$f'(\infty)=1$。

我們可以在上述所有情形中精確求解 Navier-Stokes 方程式。那為什麼仍然有一個百萬美元的公開問題？

技巧。每一個解都依賴一個技巧。幾何被選得如此精巧，以至於方程式最困難的部分（非線性項）要麼完全消失，要麼化約為可處理的東西；問題之所以可解，正是因為它已經被抽掉了所有使 Navier-Stokes 困難的成分。管流？一維。Couette？一條線。Taylor-Green 渦旋把它的非線性藏在壓力裡。

這個 千禧年獎問題問的是一般三維流動。沒有對稱性技巧。沒有簡化幾何。光滑的無散初始資料，所有尺度上的完整非線性相互作用。在那個情境中，沒有人證明解總是保持光滑，也沒有人證明它們會爆發。我們確實不知道。

所以這些精確解告訴了我們一些事，但遠遠不夠。它們證明，當你給方程式強大的對稱性可依靠時，方程式可以產生明確的光滑解。那個百萬美元問題是：在標準三維表述中，對任意光滑無散初始資料，光滑性是否總是成立，或者在湍流的全部狂暴之中，某處是否會發生災難性且不可逆的錯誤。

上面討論的每個精確解，都是透過消除或平凡化非線性對流項 $(u \cdot \nabla)u$ 來取得可處理性。在平行流中它恆等為零。在 Taylor-Green 渦旋中，它是被吸收到壓力裡的梯度。在相似解中，它透過變數變換化約為較低維的問題。

Clay 千禧年獎問題關注的是初值問題，針對具有任意光滑、快速衰減初始資料的三維非壓縮 Navier-Stokes 方程式；這正是上述簡化全都不適用的區域。

問題是以下方程式的解是否

$$\partial_t u + (u \cdot \nabla)u = -\nabla p + \nu \Delta u, \quad \nabla \cdot u = 0, \quad u(x,0)=u_0(x)$$

在 $\mathbb{R}^3$ 上對所有時間都保持在 $C^\infty(\mathbb{R}^3 \times [0,\infty))$ 中，給定光滑、無散、快速衰減的初始資料。精確解無法回答這點，因為它們棲居在解空間的子空間中，在那裡完整的非線性動力學從未真正啟動。

簡言之：精確解在高度對稱的 regime 中提供明確的光滑解。公開問題是適定性是否能延伸到不受限制的三維流動。精確表述請見 Navier-Stokes 存在性與光滑性。

若要理解方程式本身以及每一項的意義，請從 什麼是 Navier-Stokes 方程式？

若要了解這些方程式如何從第一原理建立起來，請閱讀 Navier-Stokes 方程式是如何推導出來的。

若要理解像 Poiseuille 與 Couette 這樣的層流何時會崩解為湍流，請閱讀 Reynolds 數、湍流，以及為什麼小尺度很重要。

若要理解精確解無法回答的百萬美元問題，請閱讀 千禧年獎問題：存在性與光滑性。

本站推薦接續頁面：

• 什麼是 Navier-Stokes 方程式？（控制系統及其各項）
• Navier-Stokes 方程式如何推導（從動量平衡與 Newtonian 本構律到 PDE）
• Reynolds 數、湍流，以及為何小尺度很重要（從此處描述的層流狀態過渡到完全湍流）
• 千禧年獎問題：存在性與光滑性（精確解無法解決的開放問題之精確表述）