JOVANA
Explore Library Glossary Getting Started Three Levels Fields How it works Mission
Join the mission
All guides

流體力學:尤拉方程與納維-斯托克斯方程

每一條河流、每一片機翼、每一道鋒面、每一杯被攪動的咖啡,都遵守同樣兩條定律——質量不會憑空消失,以及力等於質量乘加速度——寫成描述運動流體的偏微分方程。本文將說明尤拉方程與納維-斯托克斯方程是如何建立起來的,為何單單一個無量綱的數就決定了流動是絲滑還是混沌,以及數學究竟在何處走到了盡頭。

從 F = m a 到一整片速度場

你在整道階梯上一直在學著把偏微分方程讀作散佈於時空中的變化的記帳法。流體力學正是這份大回報:運動流體的方程不過就是物理學最古老的兩條定律——質量守恆與牛頓的 F = m a——為一個同時存在於每一點的量重新改寫而成。你不再追蹤單一顆粒子,而是追蹤一整片速度場 u(x, t):在每個位置 x 與時間 t 上,一支箭頭,告訴你流體以多快的速度、朝哪個方向滑過。密度 rho(x, t) 與壓力 p(x, t) 則作為伴隨場一同隨行。河流不是一個會移動的東西;它是一片場,逐點告訴下一瞬間該做什麼。

唯一真正全新的想法,也是一切困難的根源,是流體微團的加速度。一團流體會基於兩個理由改變它的速度:時間流逝(局部的場本身變快或變慢,這是 u_t 那一項),以及這團流體被帶往一個速度不同的區域(對流項)。跟著微團走,加速度是 u_t + (u dot grad) u。第二塊是流體在輸運它自己——速度乘上它自己的梯度——而由於 u 乘上了 u 的導數,這個方程是非線性的。這正是當一道波首次陡峭化成震波時,你遇過的那種擬線性結構:一個被它自己的值所輸運的量。流體就是把這同一種自輸運烘焙進去的守恆律。

尤拉方程:完美無摩擦的流體

和你一向的做法一樣,從理想化開始。拋掉所有內部摩擦,你就得到尤拉方程,它在 1750 年代被寫下。質量守恆給出連續方程 rho_t + div(rho u) = 0——離開一個小盒子的流體,必須表現為盒內密度的下降,這正是你替熱與交通流寫過的同一種守恆律平衡。牛頓定律給出動量方程:微團的加速度等於每單位質量所受的力,而對理想流體而言,唯一的內力是壓力由高處往低處的推擠。對不可壓縮流體(密度為常數,這是水與慢速空氣的良好模型)而言,這一對方程塌縮成一句乾淨的陳述:速度沒有源也沒有匯,而壓力會瞬間自我重排以強制這一點成立。

Incompressible Euler (no friction):

    u_t + (u . grad) u  =  - grad p / rho        (momentum: F = m a)
    div u = 0                                     (mass: nothing piles up)

Incompressible Navier-Stokes (add friction nu):

    u_t + (u . grad) u  =  - grad p / rho  +  nu * Laplacian u
    div u = 0
                       ^^^^^^^^^^^^^^^^^^^         ^^^^^^^^^^^^^^
                   nonlinear self-transport     viscous smoothing
                   (steepens, makes shocks)     (a heat-equation term)
整個故事就在這兩行上:納維-斯托克斯方程就是尤拉方程再加上一個擴散項 nu * Laplacian u。這單獨一項——喬裝過的熱算子——正是對抗非線性、並把流動平滑化的東西。

尤拉流動華麗,卻也危險。沒有摩擦,唯一對抗自我陡峭化之非線性的東西就是壓力,而對可壓縮尤拉方程——也就是氣體動力學方程——來說,這還不夠:一道平滑的流動可以一直變陡,直到它的導數在有限時間內爆掉,這正是你在非線性守恆律那裡見過的震波生成。一架超音速噴射機的弓形震波就是其中之一。而正如那裡一樣,一旦震波出現,古典解便宣告死亡,你必須轉向由熵條件所篩選出的弱解,以保住答案的唯一性。尤拉方程是誠實的第一個模型,但它沒有內建任何耗散能量的機制,因此它自身無法描述一個流體平滑地歸於靜止。

納維-斯托克斯方程:摩擦是一個隱藏的熱方程

真實的流體會黏住它自己。移動快的流層會拖著慢的一起走;這種內部摩擦就是黏滯性,把它加進尤拉方程,就得到納維-斯托克斯方程。新增的力正比於 nu * Laplacian u,其中 nu 是運動黏滯係數。仔細看那一項,你會認出一位老朋友:那個 Laplacian 正是熱方程 u_t = k u_xx 的算子。黏滯性讓速度像熱一樣擴散——快與慢的區域彼此向中間平均,陡峭的速度梯度被抹開。所以納維-斯托克斯方程是兩股你早已熟稔之力的拉鋸:試圖把流動陡峭化成震波的非線性對流,以及試圖把它平滑回去的擴散。

與其追蹤流體的速度,不如追蹤它的自旋。對速度取旋度,你就得到渦度,一個量度每個微小微團旋轉多少的場——就是你看到船槳劃過留下尾流時的那股漩渦。把納維-斯托克斯方程改寫成用渦度表達,會讓壓力完全消失,並把真正的戲劇張力暴露出來:在三維中,對流項能像拉麥芽糖一樣拉伸一根渦管,而拉伸一根旋轉的管子會使它轉得更快(角動量守恆,就像花式溜冰選手收回手臂一樣)。這種渦管拉伸正是把能量往越來越小的尺度上灌注的引擎,而它是三維所獨有的——它在二維中根本不存在。就是這一個差別,是二維溫馴、三維狂野的核心所在。

一個數決定一切:雷諾數

你如何能事先知道一道流動會是玻璃般光滑、還是翻騰的一團亂?驚人的是,它歸結於單單一個無量綱的比值。把非線性對流項的強度與黏滯平滑項相比較,你就得到雷諾數 Re = U L / nu,由典型速度 U、典型尺度 L 與黏滯係數 nu 組成。它是想搗亂的慣性想讓一切平靜下來的摩擦之間的比值。當 Re 很小,擴散獲勝:流動是層流,井然有序,一層滑過一層。當 Re 很大,非線性獲勝,流動就翻覆進入湍流

同一個數解釋了整座行為的動物園。滴落的蜂蜜(極小的 Re)由斯托克斯方程描述——把非線性項整個拿掉,你就剩下一個線性問題 nu * Laplacian u = grad p,可時間反演、溫和,是游動細菌與潤滑油膜的領域。一條河流或一片飛機機翼(巨大的 Re)則位於另一個極端。介於兩者之間,雷諾數掌控著邊界層的厚度——邊界層是固體壁面附近那一薄片,流體在其中被迫從全速降到表面上的零速。那層邊界層正是阻力誕生之處、流動分離之處,也是工程成敗的所在。

湍流,以及那道始終敞開的前線

在高雷諾數下,流動碎裂成各種尺寸的渦旋,湍流就是其結果:不規則、旋轉、細節上不可預測,卻又貫穿著統計上的規律。值得隨身攜帶的圖像是理查森的能量串級——大渦旋碎成較小的渦旋,較小的又碎成更小的,把能量沿著這條鏈往下傳,直到在最微小的尺度上,黏滯性終於咬住,把運動轉化為熱。前面提過、三維獨有的渦管拉伸,正是驅動能量沿這道階梯往下走的機制。湍流不是疊加在流動之上的雜訊;它是同一組納維-斯托克斯方程貨真價實、完全非線性的解,而且沒有任何封閉公式能描述它。

這裡是這個梯級必須說出的、誠實而令人謙卑的真相。我們擁有納維-斯托克斯方程已近兩個世紀,我們每天用數值方法求解它來設計機翼、預報天氣,然而對於三維不可壓縮的情形,沒有人知道平滑的解是否會永遠保持平滑。有可能渦管拉伸從一組完美良好的初始資料出發,在有限時間內把速度梯度逼到單一點上的無窮——一場有限時間爆破——而正反兩面都不存在證明。這就是納維-斯托克斯存在性與光滑性問題,是克雷數學研究所七大千禧年大獎難題之一,百萬美元的獎金至今仍無人領取。

既然二維已被完全解決,為什麼這還如此之難?因為在二維中,渦度只是被輸運著、有界,所以黏滯性的平滑化總是獲勝;而在三維中,渦管拉伸原則上能跑贏擴散。已知的工具——索伯列夫空間中的能量估計,正是你先前建立的那套機器——能給出全域的解(勒雷,1934 年),卻無法保證那些弱解是光滑的或唯一的。一個僅僅存在的弱解,與一個行為良好的古典解之間的鴻溝,正是這道千禧年難題要你去填補的。這是這個梯級完美的收尾音符:全物理學中實用上最成功的方程,其核心至今仍藏著一個未解的問題。