无黏性不可压缩流动.doc

二維勢流理論及複變函數的應用 複變數 複變數(Complex variable) Z可表示為 Z = x + iy， i = = re i? 式中，直角座標(x y)和極座標(n, ?)的關係為 r =， ? = tan?1 複變數基本的數學運算如下： Z = Z1 ? Z2 = (x1 ? x2) + i(y1 ? y2) Z = Z1Z2 = (x1x2 ? y1y2) + i(x1y2 + x2y1) = (r1e)(r2e) = r1 r2e Z + = (x + i y) + (x ? i y) = 2x Z ? = (x + i y) ? (x ? i y) = 2i y Z ? = (x + i y) ? (x ? i y) = x2 + y2 Zn = (re i?)n = r n e in? ?n Z = ?n (re i?) = ?n r + i? y Z = Z1Z2 Z1 + Z2 Z2 ?2 ? Z1 ?1 ??1 x 複變函數 複變函數(Complex function)可表示為 ? = f (Z) = f (x + i y) = ?(x, y) + i ?(x, y) 式中，?(x, y)和?(x, y)為任意實函數。 複變函數的導數可表示為 = f ?(Z) ? = + i (先?y = 0，再?x?0 ) = ?i + i (先?x = 0，再?y?0 ) 上式的結果，說明了複變函數為解析函數(analytic function)的意義，即複變函數的實數部份，?(x, y)，和虛數部份，?(x, y)，必須滿足柯希-黎曼條件(Cauchy-Riemann condition)。 即 = 和 = ? 對上兩式的微分，可進一步得到?和?的拉普拉斯(Laplace)方程式，即 + = ?2? = 0 及 + = ?2? = 0 因此，複變函數即為解析函數，它具有連續性及可微分的特性，且其實數部份和虛數部份均為拉普拉斯方程式的解，故均可稱為勢函數或是諧調函數(harmunic function)。 從上一章討論二維、不可壓縮的勢流中，分別存在流場的勢函數和流函數。若以?(x, y)和?(x, y)分別表示流場的勢函數和流函數，則知 ?2? = 0 和 ?2? = 0 以及速度分量為 u = = v = = ? 顯然地，上兩式u，v的關係，表示出滿足柯希-黎曼條件。 因此，如以二維、不可壓縮勢流的勢函數?(x, y)和流函數?(x, y)，分別代表複變函數的實數部份和虛數部份，則可定義該平面流動的複變勢函數(Complex potential)，即 W(Z) = ?(x, y) + i ?(x, y) 即其導數即為複變速度 = + i= u ? iv 而物理速度為 = ? i= u + iv 基本流場 茲以複變函數表示幾種二維、不可壓縮、無施性平面流場的複變勢W(Z)，其實數部份和虛數部份分別為流場的勢函數?和流函數?。 均勻流 設+x軸向的水平均勻流的U?勢函數為 W = ? + i? = U? Z = U? (x + i y) = U? r (cos? + i sin?) 所以， ? = U? x = U? r cos? ? = U? y = U? r sin? 若均勻流U?具有攻角?，或與+x軸向形成?角，而逆時鐘的?角為正，則 W = U? e?i? Z 其複變速度為 = U? e?i? = U? r (cos? ? i sin?) = u ? iv 由上知，具攻角?均勻流的速度分量(u, v) 為 u = U? cos? v = U? sin? 源流 假設源點在座標原點Z = 0，而其單位寬度源流的強度為K[L2T?1]，則 W =?n Z =?n(re i?) =[?nr + i ?] 故流場的勢函數和流函數，分別為 ? =?nr ? =? 若源流的原點在Z = Z0時，則 W =?n (Z ? Z0) 自由渦漩流 如渦漩流是一種無旋流動，則此種渦漩流稱為自由渦漩(free vortex)或無旋渦流(irrotetional vortex)。若流場環量的強度為? [LT?1]，順時鐘方向為正，則以原點為中心的自由渦漩流的複變勢為 W = ?i?nZ =? ? i?nr 若自由渦漩流的中心在Z0處，則 W = ?i?n(Z ? Z0) 式中