一個骨子裡是單一函數的流動
想像一薄層水滑過一個障礙物:每一個水平切片裡都是同一幅圖,所以整個故事是二維的。我們做出理想流體的標準理想化——流動是穩定的(圖隨時間恆定不變)、不可壓縮的(它在任何地方既不堆積也不變稀)、無旋的(丟一個小小的水車進去,它會隨波漂走卻不打轉)。在調和那一階梯你學過,恰恰是這些條件,使速度場成為某個調和函數 phi(即速度位勢)的梯度:流體順著 phi 往下流,而 phi 滿足拉普拉斯方程。所以單單一個調和函數就已主宰著這個流動。二維的魔法在於,這個實函數只不過是某個複數事物的一半。
回想調和階梯上那條連鎖。在一個漂亮的定義域上,調和的 phi 永遠有一個調和共軛 psi,這個夥伴函數使得 phi + i psi 全純。兩者由柯西-黎曼方程 u_x = v_y, u_y = -v_x 鎖在一起。本篇的頭條是:對流體流動而言,這些並非抽象的記帳——每一條柯西-黎曼方程都是一條物理定律。其中一條說流動不可壓縮;另一條說流動無旋。全純性與「漂亮的流動」,原來是同一句話,只是換了身衣裳。
複位勢:一次求導就給出速度
把速度位勢 phi 與它的共軛 psi 綁成一個全純函數。這就是複位勢 Omega(z) = phi(x, y) + i psi(x, y),整篇圍繞它而建。它的實部 phi 一如既往是速度位勢;它的虛部 psi 稱為流函數,這名字它從幾何上掙來。等值曲線 phi = 常數 是等位勢線,而等值曲線 psi = 常數 是流線——也就是一粒被流體挾帶的塵埃真正會描出的路徑。因為 phi 與 psi 是調和共軛,這兩族曲線處處互成直角相交,恰如一張被平滑彎折的方格紙上的格線。複位勢把整個流動的幾何,握在一個全純的包裹裡。
現在來到令複分析學家會心一笑的回報。速度場——兩個分量一次到位——直接從 Omega 的一次求導裡掉出來。把流速寫成 (u, v);則 Omega'(z) = u - i v。注意那個負號——導數給出的是速度的共軛,稱為複速度。所以要在任何地方求出流動的速率與方向,你只要對一個全純函數求導,再讀出答案。在 Omega'(z) = 0 之處速度消失:那是一個駐點,流體在那裡片刻靜止,宛如一塊分開溪流的岩石之鼻尖。對一個函數求一次導,便告訴你整個速度場,以及它在哪裡停頓。
Omega(z) = phi(x,y) + i psi(x,y) holomorphic complex potential phi = const equipotential lines psi = const STREAMLINES (paths the fluid follows) Cauchy-Riemann for phi, psi <==> flow is incompressible AND irrotational Omega'(z) = u - i v (complex velocity = conjugate of velocity) Omega'(z) = 0 => stagnation point
搭建流動,與共形的超能力
最簡單的複位勢就是最簡單的全純函數。取 Omega(z) = U z,其中 U 為實常數。則 phi = U x 而 psi = U y,於是流線 psi = U y 是水平直線,速度為 Omega'(z) = U,是一道以速率 U 穩定向右的流動。這就是均勻流——空空的風洞。由這顆種子與它的幾個表親,你便能組裝出一整套目錄:Omega = m log z 給出一個從原點向外噴灑流體的源(若 m < 0 則為匯);Omega = -i k log z 給出一個純渦旋,流體繞著原點打轉;而一個常數除以 z 則給出一個偶極子。因為全純函數可以相加,你便能像疊加波那樣把它們疊起來,用簡單的流動堆出精巧的流動。
現在來到讓整門學問變得強大的那一步,而它純粹是複分析。一個共形映射保角,所以它會把等位勢線與流線交織成的正交網格,從一個區域忠實地搬到另一個區域上。如果你知道某個簡單形狀上的流動——譬如均勻流,或繞圓柱的流動,這兩者你都能徒手寫出——那麼把複位勢與一個共形映射複合,就免費給了你某個複雜形狀上的流動。這正是共形映射階梯上那個移植問題的想法:在容易處求解,再把答案推過映射。流線會自動跟著過來,因為全純性在複合下保持不變。你永遠不必在那個困難的區域上重解一次拉普拉斯方程。
- 挑一個你早已知道其流動的簡單區域——最常見的是圓的外部,那裡繞圓柱的均勻流有一個整齊的封閉形式複位勢。
- 找一個共形映射,把這個簡單區域送到你真正關心的那個難搞區域上(機翼的外部、繞過一個轉角的通道)。
- 複合:困難區域上的複位勢,就是把已知位勢在逆映射處取值。全純性保證它仍描述一個有效的理想流動。
- 從移植過來的 Omega 直接讀出物理:流線取自 psi = 常數,速度取自 Omega'(z) = u - i v,駐點在該導數消失之處。
從圓到機翼:茹科夫斯基翼型
現在來為那個締造了航空的共形映射命名。茹科夫斯基映射是 w = (1/2)(z + 1/z),常省去那個一半寫作 w = z + 1/z。看它對單位圓做了什麼:令 z = e^(i theta),兩項合併為 w = cos theta,當 theta 繞行一圈時掃過 -1 到 +1 的實線段。所以這映射把單位圓壓平成一條細細的線段——一具「退化的機翼」。有趣的情形是一個通過點 z = 1、但略為偏離圓心且放大過的圓:它的像不再是線段,而是一個平滑、彎曲的翼型,前緣圓潤,後方收尖成一個銳利的尖點。
那道銳利的後緣從何而來?它是一個共形映射在單獨一點上的失靈。導數 w' = (1/2)(1 - 1/z^2) 在 z = 1 與 z = -1 處為零——這些是臨界點,是共形映射被允許打破其保角承諾的唯一地方。在臨界點處角度會被加倍,而正是這份角度加倍,把平滑的圓折出了一個角。機翼的尖點,是被映射在 z = 1 的臨界點刻意製造出來的。再來一個誠實的但書:茹科夫斯基映射在整個平面上是二對一的,因為 z 與 1/z 共用同一個像;要得到真正的共形雙射,你必須限制在單位圓的外部,而那恰好就是流體真正流動的區域。
憑空生升力:環量與庫塔-茹科夫斯基
把圓柱流動移植到機翼上後,仍有一個數是自由的:環量 Gamma,即速度繞機翼一圈的圍道積分。它量度流動繞著翼型攜帶了多少淨旋量,並由一個加進流動裡的渦旋項建成——也就是我們目錄中那個虛對數位勢。不同的 Gamma 值給出不同的、卻在數學上全都有效的繞翼流動。單憑數學並不挑出哪一個;物理必須補上那個缺失的條件。
那個缺失的條件,就是庫塔條件,是一塊偷渡進來、用以釘定那原本任意的 Gamma 的真實物理。回頭看那道由映射臨界點誕生的銳利後緣。對一般的 Gamma,理想流動的公式會預測流體以無窮的速率繞過那刀鋒般的邊緣——這在物理上荒謬。庫塔條件只是要求流動平順地離開後緣,以有限的速度乾淨地剝離,而不是繞著它纏卷。恰好有一個環量 Gamma 的值能做到這點,而這單一個要求便把它釘死。大自然挑選那個讓空氣優雅地離開機翼後方的環量。
Gamma 一旦定下,升力便以一條乾淨的公式掉出來:庫塔-茹科夫斯基定理說,單位翼展的升力為 L = rho U Gamma,其中 rho 是流體密度,U 是自由流速率,Gamma 是剛才定下的環量。升力字面上就是環量乘速率乘密度——一個圍道積分化為了血肉。對這個乾淨的故事所略去的東西要誠實:它是無黏、不可壓縮、二維流動的一種理想化。真實的空氣有黏性,還有一層緊貼機翼的薄邊界層;庫塔條件事實上是那個真正設定環量的微妙黏性效應的替身,而整幅複位勢圖像則完全忽略了阻力。它把升力預測得漂亮,卻對那阻礙飛機前進的摩擦隻字不提。
它抓住了什麼,又放掉了什麼
退一步,看清整個想法的形狀。一個物理問題——穩定、平面、無旋、無源的流動——被重新表述為一個單一的全純函數,於是前面各階梯的每件工具立刻盡數到手:調和共軛供出流線,柯西-黎曼方程化為守恆律,共形映射把容易的流動移植到困難的區域上,而一個映射的臨界點則雕出了那道使機翼成其為機翼的後緣。這同一部機器的用處遠不止空氣動力學:一模一樣的方程描述穩定熱傳(phi 是溫度,psi 是熱流線)與二維靜電學(phi 是電壓,psi 是電場線)。一個全純函數,三門物理課。
並把那些限制擺在明處,因為學生正是在那裡被燙傷的。這方法只活在二維;它需要流動無旋且無源,好讓位勢調和;茹科夫斯基升力是無黏的理想化,對阻力隻字不提;而庫塔條件是物理、不是數學,偷渡進來以挑選環量。然而在這些誠實的牆內,複位勢是整個分析學中最令人滿足的應用之一——是全純函數那抽象的剛性,伸出手把一架真實的飛機托離地面的地方。