第3章 低速平面位流ppt課件

1、1/653 31 1 理想不可壓縮流體平面位流的基本方程理想不可壓縮流體平面位流的基本方程3 32 2 幾種簡單的二維位流幾種簡單的二維位流3 32 21 1 直勻流直勻流3 32 22 2 點源點源3 32 23 3 偶極子偶極子3 32 24 4 點渦點渦3 33 3 一些簡單的流動迭加舉例一些簡單的流動迭加舉例3 33 31 1 直勻流加點源直勻流加點源3 33 32 2 直勻流加偶極子直勻流加偶極子3 33 33 3 直勻流加偶極子加點渦直勻流加偶極子加點渦3 34 4 二維對稱物體繞流的數值解二維對稱物體繞流的數值解2/65 本章討論怎樣求解不可壓理想流體無旋運動的規(guī)本章討論怎樣求解

2、不可壓理想流體無旋運動的規(guī)律。律。在理想不可壓條件下歐拉方程和連續(xù)方程包括四個在理想不可壓條件下歐拉方程和連續(xù)方程包括四個方程和四個未知函數方程和四個未知函數u,v,w,p），理論上是可解），理論上是可解的的由于飛行器的外形都比較復雜，要在滿足如此復雜由于飛行器的外形都比較復雜，要在滿足如此復雜的邊界條件下求該偏微分方程組的解析解是非常的邊界條件下求該偏微分方程組的解析解是非常困難的，原因在于方程包含非線性項，而且方程困難的，原因在于方程包含非線性項，而且方程中速度與壓強相互耦合，需要一并求出中速度與壓強相互耦合，需要一并求出3/65人們發(fā)現(xiàn)在無旋條件下問題可以得到大大簡化人們發(fā)現(xiàn)在無旋條件下

3、問題可以得到大大簡化，尤其是可以將速度和壓強分開求解，這是因，尤其是可以將速度和壓強分開求解，這是因為無旋條件可使關于速度位的方程化為線性方為無旋條件可使關于速度位的方程化為線性方程，從而便于單獨求得速度位即求出速度，而程，從而便于單獨求得速度位即求出速度，而壓強可利用伯努利方程求解壓強可利用伯努利方程求解本章的思路是，先針對理想不可壓無旋流求得本章的思路是，先針對理想不可壓無旋流求得一些典型的速度位基本解，將這些基本解進行一些典型的速度位基本解，將這些基本解進行疊加得到滿足非常簡單邊界條疊加得到滿足非常簡單邊界條 件的流動。對復件的流動。對復雜外形的繞流，介紹用基本解進行疊加的數值雜外形的繞

4、流，介紹用基本解進行疊加的數值解法大意解法大意4/65有無旋條件，就有位函數 存在，并且位函數與速度分量之間滿足：平面流動的連續(xù)方程是：結合兩式，得平面不可壓位流必須滿足的方程：該方程稱為拉普拉斯方程，是個只與速度有關的線性方程，給定適當邊界條件方程是容易求解的。ux02222yxvy0yvxu5/65對于二維不可壓縮流動，微分形式的質量方程可以寫為：對于二維不可壓縮流動，微分形式的質量方程可以寫為： 0yvxuyvxu數學上這是使數學上這是使 成為某個函數成為某個函數 的全微分的的全微分的充要條件充要條件 ，即，即 udyvdxudyvdxdvx其中：uydyydxx或：或：6/65代入無旋

5、條件：也滿足拉普拉斯方程：這也是只與速度有關的線性方程，給定邊條容易求解。位函數與流函數的關系稱為柯西黎曼條件：xyyx,02222yxyuyv7/65 拉普拉斯方程可用算子 2 表為 20。它是個線性方程，可以用疊加原理求復合的解。 所謂疊加原理是說如果有 分別滿足拉普拉斯方程，則這些函數的線性組合也必滿足拉普拉斯方程： 此外，由于速度分量與位函數之間的關系是線性的因此也滿足疊加原理： 而壓強與速度間關系為非線性故不滿足疊加原理nnnnuauaxaxaxu.111112,.,n nnaa.118/65 邊界條件是在流場邊界上規(guī)定的條件，邊界通邊界條件是在流場邊界上規(guī)定的條件，邊界通常分為內邊

6、界和外邊界。對飛行器或物體而言，內常分為內邊界和外邊界。對飛行器或物體而言，內邊界即飛行器或物體表面，外邊界為無窮遠。邊界即飛行器或物體表面，外邊界為無窮遠。數學上滿足拉氏方程的函數稱為調和函數。故要找數學上滿足拉氏方程的函數稱為調和函數。故要找一代表具體的定常不可壓理想位流運動，就是要找一代表具體的定常不可壓理想位流運動，就是要找一個能符合具體流動邊界條件的調和函數，求出位一個能符合具體流動邊界條件的調和函數，求出位函數或流函數之后，即可求出速度分布，然后用伯函數或流函數之后，即可求出速度分布，然后用伯努利方程求解壓強分布。努利方程求解壓強分布。9/65按照在邊界上所給條件是針對位函數自身還

7、是位函數的按照在邊界上所給條件是針對位函數自身還是位函數的法向導數，邊界條件分為三種類型：法向導數，邊界條件分為三種類型：（1第一邊值問題狄利希特問題）：給出邊界上位函第一邊值問題狄利希特問題）：給出邊界上位函數自身值數自身值（2第二邊值問題諾曼問題）：給出邊界上位函數的第二邊值問題諾曼問題）：給出邊界上位函數的法向導數值法向導數值（3第三邊值問題龐卡萊問題）：給出部分邊界上位第三邊值問題龐卡萊問題）：給出部分邊界上位函數自身值，部分邊界上位函數的法向導數值函數自身值，部分邊界上位函數的法向導數值氣動問題大多數屬于第二邊值問題氣動問題大多數屬于第二邊值問題10/65將坐標系與飛行器或物體固連，

8、則外邊界在遠離物體處，將坐標系與飛行器或物體固連，則外邊界在遠離物體處，速度為速度為 V ，內邊界是物體表面，不允許流體穿過或表，內邊界是物體表面，不允許流體穿過或表面法向速度為零面法向速度為零外邊界外邊界內邊界內邊界 n為物面法向為物面法向可以證明，拉普拉斯方程的解若在給定邊界上能滿足上可以證明，拉普拉斯方程的解若在給定邊界上能滿足上述條件，則解是唯一的述條件，則解是唯一的求不可壓理想無旋流繞物體的流動問題就轉化為求解拉求不可壓理想無旋流繞物體的流動問題就轉化為求解拉普拉斯方程的滿足給定邊條的特解這一數學問題普拉斯方程的滿足給定邊條的特解這一數學問題Vx0zy0nV11/65 速度位函數由無

9、旋條件定義，位函數值可以差任意常數而不影響流動。 (2) 速度位函數沿著某一方向的偏導數等于該方向的速度分量，速度位函數沿著流線方向增加。 (3) 對于理想不可壓縮無旋流動，速度位函數滿足拉普拉斯方程，是調和函數，滿足解的線性迭加原理。12/65(4) 速度位函數相等的點連成的線稱為等位線，速度方向速度位函數相等的點連成的線稱為等位線，速度方向垂直于等位線。垂直于等位線。(5) 連接任意兩點的速度線積分等于該兩點的速度位函數連接任意兩點的速度線積分等于該兩點的速度位函數之差。速度線積分與路徑無關，僅決定于兩點的位置。對之差。速度線積分與路徑無關，僅決定于兩點的位置。對封閉曲線，速度環(huán)量為零。封

10、閉曲線，速度環(huán)量為零。13/65(1) 流函數由平面不可壓縮連續(xù)條件定義，流函數值可以差任意常數而不影響流動。(2) 等流函數線是流線。即等流函數線的切線方向等流函數線是流線。即等流函數線的切線方向與速度矢量方向重合。與速度矢量方向重合。(3) 對于理想不可壓縮無旋流動，流函數滿足拉普對于理想不可壓縮無旋流動，流函數滿足拉普拉斯方程，是調和函數，解也滿足疊加原理。拉斯方程，是調和函數，解也滿足疊加原理。14/65(5) 平面內任兩點流函數的差等于通過此兩點連線的平面內任兩點流函數的差等于通過此兩點連線的流量。流量。(4) 等流函數線與等位線正交。等流函數線與等位線正交。1KKK0:,K0:,2

11、12211故：斜率，可得由斜率，可得由vuvdyudxCuvudyvdxCABBABABAddxxdyydsnVQjdsdxidsdyn , jxiyj viuV)(xyABdsnVo15/65 位函數位函數 和流函數和流函數 之間滿足柯西之間滿足柯西-黎曼條件：黎曼條件： 速度分量與位函數和流函數之間的關系是：速度分量與位函數和流函數之間的關系是：rrrrxyyx 標：坐極笛卡兒坐標：rrVrrVxyvyxur , , 標：坐極笛卡兒坐標：16/65直勻流是一種速度不變的最簡單的平行流動。其流速為流動是無旋的，由速度位全微分積分可得位函數：又可求出流函數： 流線與等位線是正交的如圖 bxay

12、byax bdyadxdcbxay cbyaxau bv 17/65常用的是這樣的直勻流，它與 x 軸平行，從左面遠方流來，流速為 。此時VyVxV18/65 點源是從流場上某一點有一定的流量向四面八方流開去的一種流動。源可以有正負。負源又名匯是一種與正源流向相反的向心流動。如果把源放在坐標原點上，那末這流動便只有 Vr，而沒有 V 。xy位于原點的點源實驗演示的點源19/65設半徑為 r 處的流速是 Vr ，那末這個源的總流量是流量是常數，故流速 Vr 與半徑成反比 rQVr2rrvQ2rVx、y 向的速度可分別寫為向的速度可分別寫為代入速度與位函數關系代入速度與位函數關系 可積分求位函數。

13、可積分求位函數。cosrVu 2222sinyxxQryrQVvryvxu,rxrQ2222yxxQ20/65比較簡便的是利用極座標下位函數與速度的關系：rVrVr,由由 位函數由上式積分得：位函數由上式積分得：rQVrr2)ln(4ln222yxQrQ（注：等位線（注：等位線C 是一系列同心圓）是一系列同心圓）21/65流函數由流函數由積分得：積分得：rrxyarctgQQ22（注：流線（注：流線c1 即即c2 是一系列射線）是一系列射線）此外注意上式中此外注意上式中的值域為的值域為-2,2,但反正切函但反正切函數的值域為數的值域為-/2,/2，故兩種表達有一定區(qū)別，故兩種表達有一定區(qū)別。r

14、QVr2xy22/65如果源的位置不在坐標原點，而在 A,處，那么22)()(ln2yxQxyarctgQ2相應的速度分量為：2222()2()()()2()()xyQxvxxyQyvyxyuv除奇點處速度無定義之外，流場其他區(qū)域都是是無旋的。23/65點源加1/2強度點匯點源加等強度點匯偶極子24/65. p 等強度的一個源和一個匯，放在x軸線上，源放在（-h，0處，匯放在0，0處。從源出來的流量都進入匯，流動情況如圖:ln)(ln22222yxyhxQ)(221Qhxyarctg1xyarctg2 其中其中1 、2 分別是點分別是點 P 與源和匯的連線與正與源和匯的連線與正 x 的夾角的夾

15、角 應用疊加原理，位函數和流函數如下應用疊加原理，位函數和流函數如下25/65現(xiàn)在我們考慮一種極限情況，當 h0 但同時 Q 增大，使 保持不變的極限情況。這時位函數變成顯然等位線=C是一系列圓心在 x 軸上的圓，且都過原點。2222202ln4),(limyxhxhyxQyxhMQh2)1ln(0(,24220limxxxyxhxQh時當22yxxM除奇點處速度無定義之外，流場其他區(qū)域都是是無旋的。26/65求流函數：上述位函數可寫為:rMcos22yMxy 利用極座標下流函數與位函數的關系：rr對積分得：rMsin即：顯然流線=C是一些圓心在 y 軸上的圓，且均過原點。2cosrM27/6

16、5兩個分速的表達式是合速要注意偶極子有軸線方向，上述布于 x 軸上的正負源形成的偶極子其軸線在x方向，對于指向正 x 方向的偶極子，上述位函數、流函數和速度分布都要改變符號。2222222cos)()(rMyxxyMxvxu22222sin)()2(rMyxxyMyvyv222rMvuV28/65如果偶極子軸線和 x 軸成角，正向指向第三象限如圖所示，在 xy 坐標系中的位函數及流函數可寫為：yxxysincossincos,xyyyxx根據二坐標系的旋轉變換關系：2,2,2,2,yxyMyxxM29/65代入上述位函數和流函數表達，并注意到坐標旋轉時向徑不變：x2+y2 = x2+y2 ，得

17、到在 (x,y) 坐標系中的偶極子：如果偶極子位于，），軸線和 x軸成角，正向指向第三象限，那么 22()cos()sin()()xyMxy22()cos()sin()()yxMxy 2222sincossincosyxxyMyxyxMyxxy30/65實際旋渦包含有旋的渦核和渦核實際旋渦包含有旋的渦核和渦核外的被誘導的無旋流場。外的被誘導的無旋流場。rV rVk / rr0p實際旋渦的渦核內為有旋流實際旋渦的渦核內為有旋流渦核外為無旋流白色粉末顯示實際旋渦的周向誘導速白色粉末顯示實際旋渦的周向誘導速度隨半徑增大而減小度隨半徑增大而減小渦核誘導流場31/65 點渦可以看成實際旋渦的渦核直徑趨于

18、零時的一種極限情況，除渦所在一點外，整個平面流場是無旋的，流體被點渦誘導繞點渦作圓周運動，流線是一些同心圓，流速只有周向速度 ，而沒有徑向速度 。繞點渦的環(huán)量是個確定的常數，例如繞半徑為 r 的圓環(huán)作環(huán)量計算，有：式中的 是個常數稱為點渦的強度，反時針方向為正。從而周向速度與離開中心點的距離 r 成反比：VrV這與無限長渦線產生的誘導速度一致。這與無限長渦線產生的誘導速度一致。rV)2(rVrV232/65由幾何條件可立刻寫出由幾何條件可立刻寫出 u 、 v 分量：分量：sinVu2222cosyxxrxrVvxyuvV位函數可由上式代入位函數可由上式代入 等后積分求出，但方便等后積分求出，但

19、方便的還是利用極座標關系：的還是利用極座標關系：uxrVr2積分后得：積分后得：xyarctg22顯然等位線顯然等位線=C是是一系列射線一系列射線ryr2222yxy33/65求流函數可由極座標下流函數與位函數的柯西黎曼關系：求流函數可由極座標下流函數與位函數的柯西黎曼關系：rr積分得：積分得：)ln(4ln222yxr顯然流線顯然流線 = C 是一系列同心圓，可見點渦與點源的位函數是一系列同心圓，可見點渦與點源的位函數與流函數只是對調了一下上述負號只是代表渦轉向）。與流函數只是對調了一下上述負號只是代表渦轉向）。如果點渦的位置不在原點，而在如果點渦的位置不在原點，而在，），則點渦的位函），則

20、點渦的位函數和流函數的式子分別是：數和流函數的式子分別是： 2yarctgx22ln2xy rV234/65 事實上沿任意形狀的圍線計算環(huán)量，值都是 ，只要這個圍線把點渦包圍在內。但不包含點渦在內的圍線，其環(huán)量卻是等于零的。點渦是實際旋渦的一種數學近似。點渦的速度在半徑點渦是實際旋渦的一種數學近似。點渦的速度在半徑 r0 時將使時將使 V 勢必使壓強勢必使壓強 p ，這是不現(xiàn)實的，這時粘性必然要起作，這是不現(xiàn)實的，這時粘性必然要起作用，因此實際的旋渦存在一個渦核，核內流體用，因此實際的旋渦存在一個渦核，核內流體 V與半徑成正比為有與半徑成正比為有旋流，核外為無旋流。實際渦核尺寸與粘性和渦強弱有

21、關，一般不大旋流，核外為無旋流。實際渦核尺寸與粘性和渦強弱有關，一般不大，故數學上抽象為一個點，形成點渦模型。，故數學上抽象為一個點，形成點渦模型。)(2222222233112211IJGICDEFCDABIJGNEFCDABrrrrrrrrrr35/65bxaybyax rQln2點源：2Q22yxxM偶極子：22yMxy 直勻流：直勻流：2點渦：rln2xy基本解位函數、流函數小結：基本解位函數、流函數小結：ab36/65 用基本的位流解可疊加成不同的流動，例如用二點源可形成如下的模擬壁面的流動：實驗演示二等強度點源模擬固璧實驗演示二等強度點源模擬固璧用固壁代替右端點源用固壁代替右端點源

22、37/65實驗演示的直均流加點源用曲線固壁代替點源用直均流加點源可模擬如下半無限壁面的流動：38/65 在一個平行于 x 軸由左向右流去的直勻流里，加一個強度為Q的源會產生如圖的流動 把坐標原點放在源所在的地方，迭加得到的位函數是：22ln4ln2),(yxQxvrQxvyxVV39/65在 x 軸上有一個合速度為零的點稱為駐點A，令 即得駐點 xA 坐標為：0AAvu0Ay 兩個分速是兩個分速是此處速度為零是因為點源速度恰好與直勻流速度相互抵消。222yxyQyvyv222yxxQvxvxuVVQxA2該速度分布的特點之一是該速度分布的特點之一是 x時，時，uV，v0。40/65 我們可以把

23、外部流動看作是在直勻流中放了一個BAB那樣形狀的物體所造成的流動，反過來也可認為繞該物體的流動可以用直勻流加點源來構造。 該半無限體在+x無限遠處，其寬度y向尺寸趨向一個漸近值D。過駐點過駐點A的流線的流線BAB是一條特殊是一條特殊的流線，把流場劃分成為兩部分的流線，把流場劃分成為兩部分。外面的是直勻流繞此圍墻的流。外面的是直勻流繞此圍墻的流動，里面的是源流在此圍墻限制動，里面的是源流在此圍墻限制之內的流動。之內的流動。41/65 流線BAB的形狀可以根據流函數=c 畫出來，也可以從流量關系推算出來。由流函數表達： CQyV2由駐點坐標y=0,=) 定常數c，得 cQ / 2 ，從而得流線BA

24、B的方程為：)(2VQy用直角坐標表達，注意到反正切的值域為-/2,/2：)32()41 (像限、像限、xyarctgxyarctg42/65 )41 ()(2象限、xyarctgVQy該流線與 y 軸交于 處，當VQ4VQyx2時，即流線在無窮遠處趨于寬度為 的直線。 VQyD2)32()(2象限、xyarctgVQy從物理上這個結果很好理解，從源流出的流量只能限制在圍線中，由速度分布知：0,vVux時，而源的流量為而源的流量為Q，以速度，以速度 V 流過時將占據寬度流過時將占據寬度 D=Q / V 43/65另一方面，流線另一方面，流線BAB的方程：的方程：)(2VQy可寫為：)(2sin

25、QrV左邊是直勻流左邊是直勻流 V 流過高流過高 y =rsin的寬度的流量，右邊則的寬度的流量，右邊則是從中心角為是從中心角為 （中流出的流量，二者相抵消，從而中流出的流量，二者相抵消，從而得流線方程的極座標表達為：得流線方程的極座標表達為：)(sin2Dr44/65 通常將壓強表為無量綱的壓強系數，其定義是當地靜壓減去來流靜壓再除以來流的動壓頭這樣得到的結果與來流參數具體值 p 、V 無關，具有通用性）：221VppCp 流場上的壓強可以用伯努利公式表達出來：流場上的壓強可以用伯努利公式表達出來： 得到表面壓強系數的表達為：得到表面壓強系數的表達為：)(2)2(222vuVpp面面2222

26、211VVVvuCp45/65將速度分布和表面流線幾何關系代入上式得到表面將速度分布和表面流線幾何關系代入上式得到表面壓強系數的結果為：壓強系數的結果為：Cp 沿沿 x 軸分布的曲線特點如圖：軸分布的曲線特點如圖：面2221VvuCp2sin2sinpC46/65直均流加變強度點源直均流加等強度點源、點匯用分布的點源、點匯構造物面直均流加偶極子實驗演示的直均流加點源和點匯的其他例子實驗演示的直均流加點源和點匯的其他例子47/65 設直勻流 平行于 x 軸，由左向右流。再把一個軸線指向負 x 的偶極子放在坐標原點處。這時，將產生如圖繞圓的流動： 流函數是：流函數是：22yxyMyV 流動的位函數

27、是：22yxxMxVV48/65圓的半徑可從駐點A的坐標定出，令：02422rMxrMVx解得：從而位函數和流函數分別寫為：sin)(2rarVcos)()(222rarVxraxV22/aVMr49/65 =0 是一條特殊的流線，這時是一條特殊的流線，這時 sin=0 ，即，即 或或 ，這就是這就是 x 軸線，還有圓表面：軸線，還有圓表面：r =a。00兩個分速的式子是：兩個分速的式子是：2sin)2cos1 (2222raVyvraVxu用在用在 的圓上時，有：的圓上時，有：ar 2sin)2cos1 (VvVu50/65將上述速度分布代入壓強系數可得：將上述速度分布代入壓強系數可得：22

28、2221sin411VVVppCp該壓強系數的分布特點如圖：該壓強系數的分布特點如圖：繞圓流動在表面上只有周向速度，沒有徑向速度：繞圓流動在表面上只有周向速度，沒有徑向速度：sin2022VvuVvvr可見在可見在/2 處速度達到最大為處速度達到最大為 2V 。51/65 達朗培爾疑題達朗培爾DAlembert，18世紀法國著名數學家提出，在理想不可壓流中，任何一個封閉物體的繞流，其阻力都是零。這個結論不符合事實。這個矛盾多少耽誤了一點流體力學的發(fā)展，那時人們以為用無粘的位流去處理實際流動是沒有什么價值的。后來才知道，這樣撇開粘性來處理問題，是一種很有價值的合乎邏輯的抽后來才知道，這樣撇開粘性

29、來處理問題，是一種很有價值的合乎邏輯的抽象，它能使我們把影響流動的各種因素分開來看清楚。譬如，早期由經驗象，它能使我們把影響流動的各種因素分開來看清楚。譬如，早期由經驗得出來的良好翼型，最大的升阻比不過是幾十比一，后來在位流理論指導得出來的良好翼型，最大的升阻比不過是幾十比一，后來在位流理論指導下，設計出來的翼型的最大升阻比竟達三百比一。這就是無粘抽象的指導下，設計出來的翼型的最大升阻比竟達三百比一。這就是無粘抽象的指導意義意義 。粘性流體繞圓柱的流動顯示實驗粘性流體繞圓柱的流動顯示實驗粘性流體繞圓柱的數值模擬粘性流體繞圓柱的數值模擬52/65xy 在直勻流加偶極子的流動之上再在圓心處加一個強

30、度為（ ）的點渦順時針轉為負），將形成如下圖的流動 這時位函數和流函數分別是：這時位函數和流函數分別是：21),(2xravyxVryravyxln21),(2V53/65在極坐標下，兩個分速是： 仍是一條流線。在這個圓上：ar 可見由于引入環(huán)量可見由于引入環(huán)量，在，在/2 處的最大速度將大于處的最大速度將大于 2V 。cos122ravrvrrravrv2sin1122VV0rv 2sin2vva V54/65224sssyaxVy或寫出駐點的直角坐標表達：或寫出駐點的直角坐標表達： ayss1sin0V 駐點的位置現(xiàn)在不在駐點的位置現(xiàn)在不在=和和=0處了，其位置可從處了，其位置可從 定出來

31、：定出來：aVs4sinxy55/65 s在第三和第四象限內，前后駐點對 y 軸是對稱的。這個角度離開和0 的多少決定于環(huán)量 對 4aV 之比值； 越大，駐點越往下移。 當點渦強度變大到 = 4aV 時，s = /2 ，二個駐點在/2處重合。 當點渦強度進一步增大使 4aV 時，駐點將離開圓柱表面，且位于圓柱之下。aVs4sinxy56/65下圖給出幾種不同點渦強度下駐點位置圖畫：顯然，有環(huán)量的繞圓流動其左右仍是對稱的，但上下已不對稱了，因此在垂直于來流的 y 方向合力就不會為零。垂直于來流方向的空氣動力分力稱為升力，可以通過沿圓柱表面壓強積分利用伯努利方程將壓強表為速度分布后積分求得），或者

32、利用動量方程求出合力。57/65下面從動量定理出發(fā)計算繞圓柱的有環(huán)量流動的升力。以原點為中心，畫一個半徑為 r1 很大的控制面 S，整個控制面還包括圓的表面 S1 以及連接 S 和 S1 的兩條割線(第二類控制體) 。注意這兩條割線上的壓力和動量注意這兩條割線上的壓力和動量進出都對消了。進出都對消了。S1 上的壓力積分是物體所受的合上的壓力積分是物體所受的合力。受力情況左右對稱，不會有力。受力情況左右對稱，不會有X 方方向合力。僅計算向合力。僅計算 Y 方向合力方向合力 L 即可。即可。設徹體力略去不計、流動定常，根據動量方程圓柱所受到設徹體力略去不計、流動定常，根據動量方程圓柱所受到的升力的

33、升力 L 可表為：可表為：dSVvdSynpLSsn)()(),cos(58/65第一個積分中的 p 按伯努利公式用速度來表達，結果得： 在在 r1 大圓上，大圓上， ，sin),cos(yndrdS12222112sin2dVvrdprLr第二個積分得： 2122/2/2212121sin1ravdravVV212212122121121cos2cos2sin1222ravdrravravrVVV59/65結果與 r1大小無關，總之合力 L 等于來流的密度乘速度 V 再乘以環(huán)量 。方向等于把直勻流的指向逆著環(huán)流轉/2，稱為升力，該結果稱為庫塔-儒可夫斯基升力定理。VraraVL2122121

34、121所以: VL考慮到速度、環(huán)量和升力之間的向量關考慮到速度、環(huán)量和升力之間的向量關系，升力定理可寫為：系，升力定理可寫為：VL 只要是封閉物體，代表其作用的正負源強度總和必須等只要是封閉物體，代表其作用的正負源強度總和必須等于零，在遠離物體的地方其作用和一個偶極子沒有什么區(qū)別于零，在遠離物體的地方其作用和一個偶極子沒有什么區(qū)別,說明物形對升力沒有直接的關系，關鍵在于必須有繞物體的說明物形對升力沒有直接的關系，關鍵在于必須有繞物體的環(huán)量存在。有了環(huán)量又有一個直勻流，便有了升力。環(huán)量存在。有了環(huán)量又有一個直勻流，便有了升力。60/65 環(huán)量之所以能產生一個 Y 向的合力，也可以從圓柱體上的壓力

35、分布直接看到。其中有環(huán)量和無環(huán)量繞流情況作了對比。無環(huán)量時，上半圓由至0上的壓力分布和下半圓由至2上的壓力分布對稱，結果是合力為零。有環(huán)量時，上半圓上的負壓遠有環(huán)量時，上半圓上的負壓遠遠超過下半圓上的負壓，所以遠超過下半圓上的負壓，所以有一個向上的合力，即升力。有一個向上的合力，即升力。這個力的來源主要靠上半圓上這個力的來源主要靠上半圓上的吸力。的吸力。61/65足球中的弧線球現(xiàn)象就是環(huán)量產生升力的例子之一：足球中的弧線球現(xiàn)象就是環(huán)量產生升力的例子之一：62/65機翼的特殊形狀使它不用旋轉就能產生環(huán)量，上部流機翼的特殊形狀使它不用旋轉就能產生環(huán)量，上部流速加快形成吸力，下部流速減慢形成壓力。速

36、加快形成吸力，下部流速減慢形成壓力。63/65 下面用解二維對稱物體繞流的例子來說明奇點疊加數值解法的應用。無迎角的對稱物體沒有升力，根據上述分析和演示，提示我們把直勻流和分布的偶極子或總強度為零的分布的點源和點匯，無環(huán)量疊加起來，得到組合流動對稱封閉物體繞流。 設直勻流速度為 V ，在 x 軸上(a,b) 范圍內，連續(xù)分布單位長度內強度設為 mx的偶極子。稱為偶極子密度。 該組合流動對任一空間點 p(x,y) 處的流函數為： dyxymyvba22V64/65 對這種無升力物體的外形可以用零流線來表示，改變不對這種無升力物體的外形可以用零流線來表示，改變不同的偶極子密度分布，可以獲得不同形狀

37、的封閉物體。同的偶極子密度分布，可以獲得不同形狀的封閉物體。 由流函數與速度的關系確定速度分布，由速度與壓強的由流函數與速度的關系確定速度分布，由速度與壓強的關系即伯努利方程確定壓強分布關系即伯努利方程確定壓強分布 。 對于實際問題，往往是給定物體的外形來確定其流動的對于實際問題，往往是給定物體的外形來確定其流動的特性。特性。 待求方程是一個積分方程，求它的解是比較困難的，但待求方程是一個積分方程，求它的解是比較困難的，但是隨著計算機技術的發(fā)展，可以用數值方法比較迅速地是隨著計算機技術的發(fā)展，可以用數值方法比較迅速地獲得這種方程的有一定準確度的數值解。獲得這種方程的有一定準確度的數值解。65/

38、65l 首先，我們把偶極子分布區(qū)域分成等寬度的 n 段，設每段的寬度為，段數 n 可根據計算機容量及結果的準確度要求而確定。l 某一定點 P(x,y) 處流函數為：l l 式中 為第 j 段的中點離原點的距離； 為第 j段內偶極子密度的平均值； 表示第 j 段內偶極子的強度。jjmjm221yxymyvjjnjV66/65l用物面邊界條件來確定待求的偶極子密度現(xiàn)在即物面為零流線，滿足0），對于給定物體外形上的 n 個已知點xi，yi），就可以得到一個對未知函數的 n 元一次聯(lián)立代數方程組： l其中 Cij 為影響系數，表示 處的單位偶極子密度對物體表面某點 Pi (xi，yi) 處的流函數的貢

39、獻。22ijiiijyxyCni,2,1jnjjijiimCyv10V67/65l 展開上式，即l 利用解一次方程組的各種計算方法，求解上面方程組，確定偶極子密度 mj。 l 一旦解得所給定物體外形的偶極子密度分布，則可確定一旦解得所給定物體外形的偶極子密度分布，則可確定流場內任意點處的流函數，此后可由流函數與速度的關流場內任意點處的流函數，此后可由流函數與速度的關系及伯努利方程，確定流場內各點處的速度及壓強值。系及伯努利方程，確定流場內各點處的速度及壓強值。11212111yvmCmCmCnn22222121yvmCmCmCnnnnnnnnyvmCmCmC2211VVV68/65 在上述過程中，我們實際上是把第在上述過程中，我們實際上是把第 j 段中分布的偶極子用段中分布的偶極子用集中在該段中點處的等強度的偶極子來代替了。顯然，集中在該段中點處的等強度的偶極子來代替了。顯然，如果分段數量較多，這種近似表示才有一定的準確性。如果分段數量較多，這種近似表示才有一定的準確性。理論上，當段數理論上，當段數 n 趨于無限大時，偶極子密度分布的數趨于無限大時，偶極子密度分布的數值結果趨近于精確解。值結果趨近于精確解。 在實際應用時，由于計算機容量和計算機機時的限制，在實際應用時，由于計