流體動力學(xué)基礎(chǔ) 流體力學(xué)

1、第四章第四章 流體動力學(xué)基礎(chǔ)流體動力學(xué)基礎(chǔ) 第一節(jié)第一節(jié) 流體的運(yùn)動微分方程流體的運(yùn)動微分方程 第二節(jié)第二節(jié) 元流的伯努利方程元流的伯努利方程 第三節(jié)第三節(jié) 總流的伯努利方程總流的伯努利方程 第四節(jié)第四節(jié) 總流的動量方程總流的動量方程 第五節(jié)第五節(jié) 理想流體的無旋流動理想流體的無旋流動 第一節(jié)第一節(jié) 流體的運(yùn)動微分方程流體的運(yùn)動微分方程 一、理想流體運(yùn)動微分方程一、理想流體運(yùn)動微分方程 在運(yùn)動的理想流體中，取微小平行六面體(質(zhì)點(diǎn))， 正交的三個邊長dx,dy,dz,分別平行于x,y,z坐標(biāo)軸(圖4 1)。設(shè)六面體的中心點(diǎn)o，速度壓強(qiáng)，分析該微 小六面體方向的受力和運(yùn)動情況。 1. 1.表面力：

2、表面力：理想流體內(nèi)不存在切應(yīng)力只有壓強(qiáng)方 向受壓面(abcd面和abcd面)形心點(diǎn) 圖 圖4141連續(xù)性微分方程連續(xù)性微分方程 的壓強(qiáng)為： 受壓面上的壓力為： 質(zhì)量力： 由牛頓第二定律 得： ( ) -( ) + dt du dxdydz x dxpp x p M 2 1 dxpp x p N 2 1 dydzpP MM dydzpP NN dxdydzXFBx dt du x x mF dxp x p 2 1 dxp x p 2 1 dydz dxdydzX 化簡得： （41） 將加速度項(xiàng)展成歐拉法表達(dá)式 : （42） dt du z p dt du y p dt du x p z y x

3、Z Y X 1 1 1 z u zy u yx u xt u z p z u zy u yx u xt u y p z u zy u yx u xt u x p zzzz yyyy xxxx uuuZ uuuY uuuX 1 1 1 uupf t u 1 用矢量表示為： （43） 上式即理想流體運(yùn)動微分方程式，又稱歐拉運(yùn)動 微分方程式。該式是牛頓第二定律的表達(dá)式，因此是 控制理想流體運(yùn)動的基本方程式。 1755年歐拉在所著的流體運(yùn)動的基本原理中 建立了歐拉運(yùn)動微分方程式，及上一節(jié)所述的連續(xù)性 微分方程式。對于理想流體的運(yùn)動，含 有和 四個未知量，由式(330)和式(336)組成的基本方 程組，

4、滿足未知量和方程式數(shù)目一致，流動可以求解。 因此說，歐拉運(yùn)動微分方程和連續(xù)性微分方程奠定了 理想流體動力學(xué)的理論基礎(chǔ)。 zyxuuu, 二、粘性流體運(yùn)動微分方程二、粘性流體運(yùn)動微分方程 1 1粘性流體的動壓強(qiáng)粘性流體的動壓強(qiáng) 理想流體因無粘滯性，運(yùn)動時不出現(xiàn)切應(yīng)力，只 有法向應(yīng)力，即動壓強(qiáng)。用類似分析流體靜壓強(qiáng) 特性的方法，便可證明任一點(diǎn)動壓強(qiáng)的大小與作用 面的方位無關(guān)，是空間坐標(biāo)和時間變量的函數(shù)， 即 （，）。 粘性流體的應(yīng)力狀態(tài)和理想流體不同，由于粘性 作用，運(yùn)動時出現(xiàn)切應(yīng)力，使任一點(diǎn)的法向應(yīng)力的 大小與作用面的方位有關(guān)。如以應(yīng)力符號的第個 下角標(biāo)表示作用面的方位， 第二個角標(biāo)表示應(yīng)力的方

5、向，則法向應(yīng)力 進(jìn)步研究證明，任一點(diǎn)任意三個正交面上的法向應(yīng)力之和 都不變，即 據(jù)此，在粘性流體中，把某點(diǎn)三個正文面上的法向應(yīng)力的平 均值定義為該點(diǎn)的動壓強(qiáng)以p表示： （44） 如此定義，粘性流體的動壓強(qiáng)也是空間坐標(biāo)和時間變量的函 數(shù) 。 zzyyxx ppp pppppp zzyyxx tzyxpp, zzyyxx pppp 3 1 2. 2. 應(yīng)力和變形速度的關(guān)系應(yīng)力和變形速度的關(guān)系 粘性流體的應(yīng)力與變形速度有關(guān)，其中法向應(yīng)力與線變形 速度有關(guān)，切應(yīng)力則與角變形速度有關(guān) 流動中某點(diǎn)的動壓強(qiáng)是過該點(diǎn)三個相互正交平面上法向應(yīng)力 的平均值，同某一平面上的法向應(yīng)力有一定差值，稱為附加法向應(yīng) 力，以

6、表示，它是流體微團(tuán)在法線方向上發(fā)生線變形 (伸長或縮短)引起的。 （45） 切應(yīng)力與角變形速度的關(guān)系，在簡單剪切流動中符合牛頓內(nèi) 摩擦定律 zzyy ppp xx , z u zzzz y u yyyy x u xxxx z y x pppp pppp pppp 2 2 2 dy du u 將牛頓內(nèi)摩擦定律推廣到一般空間流動,得出 （46） 3 3粘性流體運(yùn)動微分方程粘性流體運(yùn)動微分方程 采用類似于推導(dǎo)理想流體運(yùn)動微分方程式（41）的方法， 取微小平行 六面體，根據(jù)牛頓第二定律建立以應(yīng)力(包括切 應(yīng)力)表示的運(yùn)動微分方程式，并以式(45)、式(46)代人 整理，使得到粘性液體運(yùn)動微分方程： y

7、 u x u yxxy x u z u xzzx z u y u zyyz x y zx y z 用矢量表示為 式中： 拉普拉斯算子。 粘性流體運(yùn)動微分方程，又稱為納維 斯托克斯方程(簡寫為NS方程)。 z u zy u yx u xt u zz p z u zy u yx u xt u yy p z u zy u yx u xt u xx p zzzz yyyy z y xx uuuuZ uuuuY uuuuX 2 1 2 1 2 1 uuupf t u 2 1 2 2 2 2 2 2 2 zyx （48） NS方程表示作用在單位質(zhì)量流體上的質(zhì)量力、表面力(壓力 和粘性力) 的相平衡。由NS

8、方程式和連續(xù)性微分方程式組成的基 本方程組，原則上可以求解速度場和壓強(qiáng)場p,可以說粘性流體的運(yùn) 動分析，歸結(jié)為對NS方程的研究。 例41 理想流體速度場為 為常數(shù)。試求：（1）流動是否可能實(shí)現(xiàn)；（2）流線方程；（3）等 壓面方程(質(zhì)量力忽略不計(jì)) 解 (1)由連續(xù)性微分方程 滿足連續(xù)性條件，流動是可能實(shí)現(xiàn)的。 (2)由流線方程 得 ： baubxuayuzyx,0, 0 z u y u x uzyx yxu dy u dx bx dy ay dx aydybxdx 積分得流線方程 a,b同號，流線是雙曲線a,b異號,流線是圓。 (3)由歐拉運(yùn)動微分方程式，不計(jì)質(zhì)量力： 將方程組化為全微分形式:

9、 aby x u u abx y u u y xy p x yx p 1 1 )( 1 )()( 1 ydyxdxabdp ydyxdxabdy y p dx x p caybx 22 積分，得 令p=常數(shù) 即得等壓面方程 等壓面是以坐標(biāo)原點(diǎn)為中心的圓。 2 22 c yx abp cyx 22 第二節(jié)第二節(jié) 元流的伯努利方程元流的伯努利方程 一、理想流體運(yùn)動微分方程的伯努利積分一、理想流體運(yùn)動微分方程的伯努利積分 理想流體運(yùn)動微分方程式是非線性偏微分方程組，只有特 定條件下的積分，其中最為著名的是伯努利(Daniel Bernoull， 17001782，瑞士科學(xué)家)積分。 （410） z

10、u zy u yx u xz p z u zy u yx u xy p z u zy u yx u xx p zzz yyy z y x uuuZ uuuY uuuX 1 1 1 由理想流體運(yùn)動微分方程式 各式分別乘以沿流線的坐標(biāo)增量dx，dy，dz，然后相加得： 、 1.引人限定條件： 作用在流體上的質(zhì)量力只有重力：X=Y=0,Z=-g; dzdydxZdzYdyXdx z p y p x p 1 )( dz dt du dy dt du dx dt du z y x gdzZdzYdyXdx)( dt du z p dt du y p dt du x p z y x Z Y X 1 1 1

11、 .不可壓縮，恒定流：,Czyxpp, p z p y p x p ddpdzdydx 11 .恒定流流線與跡線重合：dx=uxdt，dy=uydt，dz=uzdt 則 2 222 zyx z y x uuu ddz dt du dy dt du dx dt du 相加帶入后得： dzdydxZdzYdyXdx z p y p x p 1 )( dz dt du dy dt du dx dt du z y x g up z g up z 22 2 22 2 2 11 1 Cgz u g p 2 2 Cz g u p 2 2 理想流體運(yùn)動微分方程沿流線的積分稱為伯努利積分， 由于元流的過流斷面積

12、無限小，所以沿流線的伯努利方 程就 是元流的伯努利方程。推導(dǎo)該方程引入的限定條件，就 是理想流體元流伯努利方程的應(yīng)用條件，歸納起來有：理想 流體；恒定流動；質(zhì)量力中只有重力；沿元流(流線)；不可壓 縮流體。 1. 1.物理意義式物理意義式 式子中的前兩項(xiàng) 、 和 的物理意義分 別是單位重量流體具有的比位能壓能或比勢能；單位重量 流體具有的動能。 g p z g u 2 2 g p z , 三項(xiàng)之和 是單位重量流體具有的機(jī)械能，沿 同一無流(沿同一流線)。單位重量流體的機(jī)械能守恒。伯努利 方程又稱為能量方程。 g u g p z 2 2 2. 2.流體意義流體意義 各項(xiàng)的流體力學(xué)意義為：z是位置

13、水頭， 壓強(qiáng)水頭； 兩項(xiàng)之和 是測壓管水頭， 是流速水頭，三項(xiàng)之和 稱為總水頭 g p g u 2 2 g p zHp g u g p zH 2 2 圖42水頭線 表示理想流體的恒定流動，沿同一元流(沿同一流線)各斷式 (423)則面的總水頭相等理想流體的水頭線是水平線 3. 3.幾何意義 幾何意義 各項(xiàng)的幾何意義是不同的幾何高度：z是位置高度，測壓管高 度。總結(jié)如下： p 例例42 應(yīng)用皮托(Pito，H.)管測量點(diǎn)流速 前文指出，流速水頭可直接量測，現(xiàn)以均勻管流為例加以說明。 設(shè)均勻管流，欲量測過流斷面上某點(diǎn)A的流速(圖43)。在該點(diǎn)放置 一根兩端開口，前端彎轉(zhuǎn)90的細(xì)管，使前端管口正對來

14、流方向， 另一端垂直向上，此管稱為測速管。來流在A點(diǎn)受測速管的阻滯速 度為零，動能全部轉(zhuǎn)化為壓能測速管中液面升高。 另在A點(diǎn)上游的同一流線取相距很近的o點(diǎn)，因這兩點(diǎn)相距很 近，o點(diǎn)的壓強(qiáng)p實(shí)際上等于放置測速管以前A點(diǎn)的壓強(qiáng) 應(yīng)用理想流 體元流伯努利方程： （425） （426） 0 2 2 h g p g p g u g p g u g p 2 2 圖43點(diǎn)流速的測量 式中o點(diǎn)的壓強(qiáng)水頭，由另根測壓管量測， 于是測速管和測壓管中液面的高度差,就是A 點(diǎn)的流速水頭，該點(diǎn)的流速： （427） 根據(jù)上述原理，將測速管和測風(fēng)管組合 成測量點(diǎn)流速的儀器，圖44所示，與迎流 孔(測速孔)相通的是測速管，與

15、側(cè)面順流孔 （測壓孔或環(huán)形窄縫）相通的是測壓管?？?慮到粘性流體從迎流孔至順流孔存在粘性效 應(yīng)，以及皮托管隊(duì)員流場的干擾等影響，引 用修正系數(shù)C： 02 2ghC g pp gCu 圖44 畢托管構(gòu)造 02 2gh g pp gu 式中C是修正系數(shù)數(shù)值接近于1.0，由實(shí)驗(yàn)測定。 【例4-3】 有一貯水裝置如圖（4-5）所示，水池足夠大，當(dāng)閥 門關(guān)閉時，壓強(qiáng)計(jì)讀數(shù)為2.8個大氣壓強(qiáng)。而當(dāng)將閥門全開，水從管 中流出時，壓強(qiáng)計(jì)讀數(shù)是0.6個大氣壓強(qiáng)，試求當(dāng)水管直徑d=12cm 時，通過出口的體積流量(不計(jì)流動損失)。 【解解】 當(dāng)閥門全開時列1-l、2-2截面的伯努利方程 當(dāng)閥門關(guān)閉時，根據(jù)壓強(qiáng)計(jì)的

16、讀數(shù)， 應(yīng)用流體靜力學(xué)基本方程 ， 求出值： g V g pp g p H aaa 2 6.0 00 2 2 圖45 aaa ppgHp8 . 2 OmH g p H a 2 28 9806 980608 .28 .2 所以管內(nèi)流量： sm g p HgV a /78.20 9806 980606 .0 8 .2806.92 6 .0 2 2 smVdq V /235.078.20785.0 4 3 2 2 三、粘性流體元流的伯努利方程三、粘性流體元流的伯努利方程 能量守恒原理：能量可以從一種形式轉(zhuǎn)換成另一種形式， 既不能創(chuàng)造、也不能消滅，總能量是恒定的 粘性流體流動時，單位重量流體具有的機(jī)械

17、能沿程減少， 總水頭線是沿程下降。 因此，設(shè)為粘性流體元流單位重量流體由過流斷面11 運(yùn)動至過流斷面22的機(jī)械能損失，稱為元流的水頭損失， 根據(jù)能量守恒原理，便可得到粘性流體元流的伯努利方程 水頭損失 也具有長度的量綱。 w h g up z 21 2 11 22 2 22 wg up hz 第三節(jié)第三節(jié) 總流的伯努利方程總流的伯努利方程 上一節(jié)的最后得到了粘性液體元流的伯努利方程式，為了 解決實(shí)際問題，還需要將其推廣到總流中去。 一、漸變流及其性質(zhì)一、漸變流及其性質(zhì) 漸變流： 凡質(zhì)點(diǎn)的遷移加速度(位變加速度)很小，的流動， 或者說流線近于平行直線的流動漸變流是均勻流的寬延，所以 均勻流的性質(zhì)

18、，對于漸變流都近似成立，主要是： 1漸變流的過流斷面近于平面。面上各點(diǎn)的速度方向近于 平行； 2恒定漸變流過流斷面上的動壓強(qiáng)按靜壓強(qiáng)的規(guī)律分布，即： cz p 二、總流的伯努利方程二、總流的伯努利方程 設(shè)恒定總流，過流斷面11、22為漸變流斷面，面積為 A1,A2由元流的伯努利方程： 圖47急變流和漸變流 g up z 21 2 11 22 2 22 wg up hz 單位時間通過元流兩過流斷面的能量關(guān)系 總流是由無數(shù)元流構(gòu)成的，上式對總流過流斷面積分 便得到單位時間通過總流兩過流斷面的總能量關(guān)系 分別確定三種類型的積分 第一類積分： 因所取過流斷面是漸變流斷面 dQz g up 21 2 1

19、1 dQz g up 22 2 22 dQhw 1 1 111 A p dAuz 1 2 1 112 A g u dAu 2 2 222 A p dAuz 2 2 2 222 A g u dAu Q w aQh A p udAz cz p 第二類積分： 動能校正系數(shù)， 式中 是為校正以斷面平均速度計(jì)算的動能與實(shí)際功 能的差異而引入的校正系數(shù)， 值取決于過流斷面上 的流速分布情況，分布均勻的流動 。 通常取 A g u udA 2 2 A g u udA 2 2 A g u dA 2 3 Q g v 2 2 Av dAu dA dA A A g v A g u 3 3 2 3 2 3 10. 1

20、05. 11 Q p zudAz A p 第三類積分： 為總流單位重量流體由11至22斷面的平均機(jī)械 能損失，稱總流的水頭損失 兩斷面間無分流及匯流，Q1Q2Q，并以 除 上式，得 Q w dQh w h QhdQh w Q w QhQ g v Q p z Q g v Q p z w 2 2 2 22 2 2 11 1 2gQ wg vp g vp hzz 2221 2 222 2 111 2. 2. 伯努利方程的適用條件伯努利方程的適用條件 式式(437)(437)即粘性流體總流的伯努利方程。將元流的伯努利方即粘性流體總流的伯努利方程。將元流的伯努利方 程推廣為總流的伯努利方程，引入了某些限

21、制條件，也就是程推廣為總流的伯努利方程，引入了某些限制條件，也就是 總流伯努利方程的適用條件包括：總流伯努利方程的適用條件包括： . .不可壓縮流體恒定流；不可壓縮流體恒定流； . .質(zhì)量力只有重力；質(zhì)量力只有重力； 不可壓縮流體不可壓縮流體( (以上引自粘性流體元流的伯努利方程以上引自粘性流體元流的伯努利方程) )； . .所取過流斷面為漸變流斷面；所取過流斷面為漸變流斷面； . .兩斷面間無分流和匯流；兩斷面間無分流和匯流； . .兩斷面間無能量的輸入或支出；兩斷面間無能量的輸入或支出； . .不存在相對運(yùn)動。不存在相對運(yùn)動。 3. 伯努利方程的方法步驟伯努利方程的方法步驟 . .斷面選擇

22、斷面選擇 通常選擇未知量所在的斷面和已知量最多的斷面， 它們 都必須是漸 變流斷面； . .代表點(diǎn)選擇代表點(diǎn)選擇 無壓流一般選擇自由液面，有壓流一般選在管 道中心； . .位置基準(zhǔn)面選擇位置基準(zhǔn)面選擇 習(xí)慣選擇在過各代表點(diǎn)最低者的水平面。 位置準(zhǔn)面選擇對結(jié)果無影響； . .壓強(qiáng)基準(zhǔn)面選擇壓強(qiáng)基準(zhǔn)面選擇 液體一般選取相對壓強(qiáng)；氣體一般選取絕 對壓強(qiáng)。壓強(qiáng)準(zhǔn)面選擇對結(jié)果無影響； . .列伯努利方程列伯努利方程 對于初學(xué)者，應(yīng)該分項(xiàng)列出，哪怕是零，也 應(yīng)該寫 出。但一般只用符號代替，而不代入具體數(shù)值，以 便推導(dǎo)出未知量的計(jì)算公式； . .解伯努利方程解伯努利方程 求解出題目中所要求的未知量； . .

23、給出答案給出答案 給出正確的答案 例例43 用直徑d100mm的水管從水箱引水(圖49)。水箱水面與 管道出口斷面中心的高差H4m保持恒定，水頭損失 3m水柱。 試求管道的流量。 解解 這是一道簡單的總流問題，應(yīng)用伯努利方程： 圖49管道出流 w h g vp z 21 2 11 wg vp hz 22 2 22 求解的關(guān)鍵是“三選”：選基準(zhǔn)面、計(jì)算斷面和計(jì)算點(diǎn)。為便 于計(jì)算，選通過管道出口斷面中心的水平面為基準(zhǔn)面00(圖4 9)。計(jì)算斷面應(yīng)選在漸變流斷面，并使其中一個已知量最多，另 一個含待求量。技以上原則本題選水箱水面為11斷面，計(jì)算點(diǎn) 在自由水面上、運(yùn)動參數(shù)z1=H,p1=0 (相對壓強(qiáng)

24、), v1=0 。選管道出 口斷面為22斷面，以出H斷面的中心運(yùn)動參數(shù)z2=0,p2=0, v2待求。 將各量代人總流伯努利方程： 取 得： wg v hH 2 2 2 0 . 12 smhHgvw/43. 4)(2 2 四、總流伯努利方程應(yīng)用的修正四、總流伯努利方程應(yīng)用的修正 伯努利方程是古典水動力學(xué)應(yīng)用最廣的基本方程。應(yīng)用伯努 利方程要重視方程的應(yīng)用條件，切忌不顧應(yīng)用條件，隨意套用公 式，要對實(shí)際問題做具體分析，靈活運(yùn)用。下面結(jié)合三種情況加 以討論。 1.氣體的伯努利方程 總流的伯努利方程式(436)是對不可 壓縮流體導(dǎo)出的，氣體是可壓縮流體，但 是對流速不很大(60ms)，壓強(qiáng)變化不 大

25、的統(tǒng)，如工業(yè)通風(fēng)管道、煙道等，氣流 在運(yùn)動過程中密度的變化很小，在這樣的 條件下，伯努利方程仍可用于氣流。由于 氣流的密度同外部空氣的密度是相同的數(shù)量級，在用相對壓強(qiáng)進(jìn) 行計(jì)算時，需要考慮外部大氣壓在不同高度的差值。 設(shè)恒定氣流(圖410)、氣流的密度為 外部空氣的密度為 ， 過流斷面上計(jì)算點(diǎn)的絕對壓強(qiáng) 。 列11和22斷面的伯努利方程式： a absabsPP21, 圖410恒定氣流 wg vp hz 22 2 22 g vp z 21 2 11 （438） 進(jìn)行氣流計(jì)算，通常把上式表示為壓強(qiáng)的形式 （439） 式中pw為壓強(qiáng)損失 （440） 將式(439)中的壓強(qiáng)用相對壓強(qiáng)p1,p2表示，

26、則： （441） （442） 式中 為 處的大氣壓， 為高程 處的大 壓，代人式(437)，整理得： （443） 1 21 2 2 1 11 v pz abs wabs p v pz 2 2 2 22 wwghp aabs ppp 11 1222 zzppp aaabs 12 zzp aa a p 1 z 2 z 1221 2 1 zzp a v w v pp 22 2 2 這里 稱為靜壓； 稱為動壓。 為單位體積氣體所受有效浮力， 為氣體沿 浮力方向升高的距離，乘積 為11斷面相對于22 斷面單位體積氣體的位能，稱為位壓。 式（442）就是以相對壓強(qiáng)計(jì)算的氣流伯努利方程。 當(dāng)氣流的密度和外界

27、空氣的密度相同 ，或兩計(jì)算點(diǎn)的高 度相同 時，位壓為零，式（442）化簡為： （444） 式中靜壓與動壓之和稱為全壓。 當(dāng)氣流的密度遠(yuǎn)大于外界空氣的密度( )，此時相當(dāng) 于液體總流，式(443)中 可忽略不計(jì)，認(rèn)為各點(diǎn)的當(dāng)?shù)卮髿?壓相同，式(443)化簡為： 21, p p 2 , 2 2 2 2 1 vv g a 12 zz 12 zzg a a 21 zz 21 2 1 v p w v pp 22 2 2 a a （445） 除以 ，即 （446） 由此可見，對于液體總流來說，壓強(qiáng) 不論是絕對壓強(qiáng)， 還是相對壓強(qiáng)，伯努利方程的形式不變。 2.有能量輸入或輸出有能量輸入或輸出 總流伯努利方程

28、式(437)是在兩過流斷面問除水頭損失之外， 在無能量輸入或輸出的條件下導(dǎo)出的。當(dāng)面過流斷面間有水泵、 風(fēng)機(jī)(圖411)或水輪機(jī)(圖412)等流體機(jī)械時，存在能量的輸入 或輸出。 此種情況，根據(jù)能量守恒原理，計(jì)入單位重量流體經(jīng)流體機(jī)械 獲得或失去的機(jī)械能， 1221 2 1 zzp v wg v pp 22 2 2 g g v z p 2 2 1 1 1 w p h g v z 2 2 2 2 2 21, pp 式(429)便擴(kuò)展為有能量輸入或輸出的伯努利方程式： （447） 式中：+H表示單位重量流體通過流體機(jī)械(如水泵)獲得的機(jī)械 能，對于水泵稱為水泵的揚(yáng)程； -H 表示單位重量流體給流體

29、機(jī)械(如水輪機(jī))的機(jī)械 能，又稱為水輪機(jī)的設(shè)計(jì)水頭。 H p z g v 2 1 1 2 11 w p h v z 2 2 2 2 2 圖412有能量輸出的總流 圖411有能量輸入的總流 3.兩斷面間有分流或匯流兩斷面間有分流或匯流 總流的伯努利方程式(436)，是在 兩過流斷面間無分流和匯流的條件下導(dǎo) 出的。而實(shí)際的供水供氣管道沿程多有 分流和匯流這種情況式(436)是否還 能用呢?對于兩斷面間有分流的流動(圖 413)，設(shè)想11斷面的來流，分為兩 股(以虛線劃分)分別通過22、33斷面。 對 (11斷面中的一部分)和22斷 面列伯努利方程，其間無分流： （448） 圖413沿程分流 1 1

30、 g v g p z 2 1 1 2 1 2 1 2 22 2 2 w h g v g p z 因所取11斷面為漸變流斷面。面上各點(diǎn)的勢能相等，則： （449） 如11斷面流速分布較為均勻，則： （450） 故 （451） 近似成立。同理可得： （452） 由以上分析，對于實(shí)際I程中沿程分流的總流，當(dāng)所取過流斷面為漸 變流斷面，斷面上流速分布較為均勻，并計(jì)人相應(yīng)斷面之間的水頭 損失。 g P Z g P Z 1 11 g v g P Z g v g P Z 22 2 1 1 1 2 11 1 21 2 2 2 1 2 1 22 21 w P h g v g Z g v g p Z 31 2 3

31、 3 3 2 1 1 1 22 wh g v g P Z g v g P Z g v g v 22 2 1 2 1 第四節(jié)第四節(jié) 總流的動量方程總流的動量方程 總流的動量方程是繼連續(xù)性方程式、伯努利方程式(436)之后 的第三個積分形式基本方程，它們在流體力學(xué)及水力學(xué)中習(xí)慣地 被稱為三大方程，下面由動量原理，推導(dǎo)總流的動量方程。 一、總流的動量方程一、總流的動量方程 設(shè)恒定總流，取過流斷面、為漸變流斷面，面積 為以過流斷面及總流的例表面圍成的空間為控制體(圖314)?？?制體內(nèi)的流體，經(jīng)dt時間，由運(yùn)動到位置。 在流過控制體的總流內(nèi)，任取元流12,斷面面積dA1,dA2，點(diǎn) 流速為 ，dt時間

32、，元流動量的增量 （453） （454） 21, uu 2121 KKKd dtt KK 22212111 KK 1122 KKKd 2221KK 2111 KK dt時間，總流動量的增量，因?yàn)檫^流斷面為漸變流斷面，各點(diǎn)的 流速平行，按平行矢量和的法則，定義為方向的基 本單位向量，為方向的基本單位向量 （455） 對于不可壓縮液體，并引入校正系數(shù)，以斷面 平均流速v代替點(diǎn)流速 積分得： （456） 式中 是為校正以斷面平均速度計(jì)算的動量與實(shí)際動量的差 異而引入的校正系數(shù)，稱為流速分布不均勻動量校正系數(shù)： （457） 2i2u 1i1u Kd 2 2222 2 iudtdAu A 1 1111

33、1 iudtdAu A Kd 21 2 222 vAvdt 1 111 vAvdt 1122 vvdtQ dtF Av dAu A 2 2 值取決于過流斷面上的速度分布，速度分布較均勻的流動， 1.021. 05，通常取 1.0 由動量原理，質(zhì)點(diǎn)系動員的增量等于作用于該質(zhì)點(diǎn)系上的外 力的沖量： (458) 投影式： (459) 式(458)、式(459)就是恒定總流的動量方程。方程表 明，作用于控制體內(nèi)流體上的外力，等于單位時間控制體流出動 量與流人動量之差。綜合推導(dǎo)式(447)規(guī)定的條件，總流動量方 程的應(yīng)用條件有：恒定流；過流斷面為漸變流斷面，不可壓縮流 體。 dtF 1122 vvdtQ

34、 F 1122 vvQ zzz yyy xxx vvQF vvQF vvQF 1122 1122 1122 二、動量方程應(yīng)用舉例二、動量方程應(yīng)用舉例 【例例39】 水平放置在混凝土支座上的變直徑彎管，彎管兩端 與等直徑管相連接處的斷面11上壓力表讀數(shù)p1=17.6104Pa ， 管中流量qv=0.1m3/s，若直徑d1=300，d2=200，轉(zhuǎn)角=60， 如圖414所示。求水對彎管作用力F的大小。 【解解】 水流經(jīng)彎管，動量 發(fā)生變化，必然產(chǎn)生作用力F。而 F與管壁對水的反作用力R平衡。 管道水平放置在xoy面上，將R分 解成Rx和Ry兩個分力。 取管道進(jìn)、出兩個截面和管內(nèi)壁 為控制面，如圖所

35、示，坐標(biāo)按圖示方向設(shè)置。 圖414 .根據(jù)連續(xù)性方程可求得： .列管道進(jìn)出口的伯努利方程 ，則： .所取控制體受力分析，進(jìn)、出口控制面上得總壓力： sm d q v v /42. 1 3 . 0 41 . 0 4 2 2 1 1 sm d q v v /18. 3 2 . 0 41 . 0 4 2 2 2 2 g v g p g v g p 22 2 22 2 11 2 2 2 2 112 vvpp218. 342. 11000106 .17 223 Pa 3 102 .17 43.123 .0 4 106 .17 23 111 ApP 40.52 .0 4 106 .17 23 222 Ap

36、P （kN） （kN） 壁面對控制體內(nèi)水的反力Rx、Ry，其方向先假定如圖(414)所示。 .寫出動量方程 選定坐標(biāo)系后，凡是作用力（包括其分力）與坐標(biāo)軸方向一致 的，在方程中取正值；反之，為負(fù)值。 沿x軸方向 沿y軸方向 coscos 1 221 vvqRPP Vx coscos 122 1 PPvvqR Vx 568. 060cos43.1240. 560cos42. 118. 31 . 0 (KN) sin0sin 11 vqRP Vx sinsin 11 vqPR Vy 88.1060sin42. 11 . 060sin43.12 (KN) 管壁對水的反作用力 水流對彎管的作用力F與R

37、大小相等，方向相反?？偭鲃恿糠匠?是動量原理的總流表達(dá)式，方程給出了總流動量變化與作用力之 間的關(guān)系。根據(jù)這一特點(diǎn)，求總流與邊界面之間的相互作用力問 題，以及因水頭損失難以確定運(yùn)用伯努利方程受到限制的問題， 適于用動量方程求解。 三、動量矩方程三、動量矩方程 上面對動量定理的推導(dǎo)過程中所用之方法、步驟，對動量矩定 理也完全適用，而所得結(jié)果與動量定理完全相似，只要在以上的 相應(yīng)式個，將動量換成動量短就成為動量矩定理；這里不作重復(fù) 的推演。 恒定流動的動量矩定理為： 89.1088.10568. 0 2 2 22 YX RRR (KN) 上式表明，在流出面上的流出動量 矩與流入面上的流入動量矩之差

38、等于 外力矩之和。 常見的流體機(jī)械中，離心式水泵、 風(fēng)機(jī)都是將其機(jī)械能轉(zhuǎn)換為流體的動 能和壓能的。水輪機(jī)則是利用流體的 動能使葉片機(jī)械轉(zhuǎn)動向外輸出功率， 其工作原理都是相同的。 圖415表示水輪機(jī)葉輪的兩個葉 片所形成的槽道，流體自葉輪外徑 的圓周面流入槽道， 經(jīng)葉輪內(nèi) 徑的 圓周面流出槽道，進(jìn)入葉輪 中心區(qū)域的導(dǎo)管沿軸向流出；葉輪葉 片就是在流體流動時獲得力矩而轉(zhuǎn)動 向外作功的。 iin A n A FrdAVVrdAVVr INou (460) 圖415 1r 2r 假定葉片數(shù)目足夠多，則葉片間的槽道可近似為一元流動，各 截面上的速度是均勻的。還假定葉輪作等角速 的旋轉(zhuǎn)，則葉 輪個流場雖為

39、不定常，但葉輪中的總體動量矩不隨時間變化，可 適用定常的動量矩公式，下面我們來導(dǎo)出水輪機(jī)(也稱渦輪機(jī))的動 量矩公式。 先選取控制面：半徑的 進(jìn)口圓周團(tuán)和半徑 的出口圓周團(tuán) 之間的流體表面，其中包括各葉片與流體的接觸面； 現(xiàn)在分析控制面上的運(yùn)動情況及受力情況。設(shè)流體以相對速度 經(jīng)半徑 的圓周團(tuán)流入葉片槽道，由于半徑 的圓周速度即牽 連速度 ，則流體流入槽道的絕對速度為 (461) 設(shè)絕對速度為 與圓周切向夾角為 則其徑向分量 和周 向分量 的大小分別為： (462) 1r 2 r rv 1 r 1 r 1 1rVe 111erVVV 111sinVVn 111cosVVt 1V1 1Vn 1V

40、t (463) 同理，流體在流出半徑 圓周面上的相對速度 ，牽連速 度 ，則絕對速度為 (464) 設(shè)絕對速度為 與圓周切向夾角為 ，則其徑向分量 和周 向分量 的大小分別為： (465) (466) 在流量為Q的情況下，流出控制面的動量矩為其切向動量 與半徑 的乘積，即： (467) 同理，流入控制團(tuán)的動量矩為其切向動量與半徑之乘積，即 ： (468) 假定無粘性力作用，則控制面中的兩圓周面上的壓力合力不產(chǎn) 生力矩，只有葉片對流體的作用力矩。 2r 2Vr 22rVe 222erVVV 2V2 2Vn 2Vt 222sinVVn 222cosVVt 2tQV 2r 22222cosrQVrQ

41、Vt 11111cosrQVrQVt 則根據(jù)動量矩定理，(464)式減(465)式等于外力矩： (469) 根據(jù)作用反作用原理，葉片上獲得流體所給的作用力矩力 （470） 這就是歐拉渦輪方程式，是渦輪機(jī)械的基本方程式。葉輪所獲得 的功率為 （471） 當(dāng)流出葉片槽道的絕對速度 的方向取半徑方向，即 時， 則葉輪獲得的力矩公式變?yōu)?（472） 相應(yīng)地，葉輪所獲得的功率公式為 （473） 11221112220coscosrVtrVtQrVrVQM 22112221110coscosrVtrVtQrVrVQM 22112211110coscosVeVtVeVtQVeVVeVQMp 2V 90 11

42、1110cosrQVrQVMt 11111coseteVQVVQVP 第五節(jié)第五節(jié) 理想流體的無旋流動理想流體的無旋流動 在第三章中，在微團(tuán)運(yùn)動分析的基礎(chǔ)上，見流體的運(yùn)動分為有旋流動和無旋 流動。理論研究證明只有不可壓縮理想流體，運(yùn)動初始無旋。嚴(yán)格地說，粘性流 體的運(yùn)動都是有旋流動，但在實(shí)際流動中，多有粘性的影響很小，從靜止轉(zhuǎn)入流 動（初始無旋）的情況，諸如通風(fēng)車間，用吸風(fēng)裝置抽氣，工作區(qū)內(nèi)形成的氣流; 水庫中的靜水，因閘門開啟形成的閘孔出流或堰流；以及空氣或水繞物體流動 時，在邊界層外面，廣闊區(qū)域的流動等，都可視為無旋流動。 一、勢函數(shù)一、勢函數(shù)： 根據(jù)曲線積分定理，無旋流的條件式（550）

43、是表達(dá)式 成為某一函數(shù)的全微分的必要和充分條件 （474） （475） 得： ， ， （ 476） zzyyxxdududu dzudyudxud zyx dzdydxd zyx xx u yy u zz u gradu 函數(shù) 仿照應(yīng)力場勢函數(shù)，靜電場勢函數(shù)的定義，稱為 速度勢函數(shù)。由此得出，無旋流是有速勢的流動，簡稱勢流; 反之，有速勢的流動是無旋流，兩者含義相同。 將式（4-37）不可壓縮流體連續(xù)性微分方程 ： （478） 即： （479） 式中 拉普拉斯算子 式（478）是著名的拉普拉斯方程，滿足拉普拉斯方程的函 數(shù)是調(diào)和函數(shù)。所以，調(diào)和函數(shù)的一切性質(zhì)，也是速度勢函數(shù) 擁有的性質(zhì)。 ),

44、(zyx x u x y u y z u z 0 2 2 2 2 2 2 zyx 0 2 2 2 2 2 2 2 2 zyx 由以上分析可知，不可壓縮流體無旋流動的問題，歸結(jié)為在 給定的邊界條件下，求解拉普拉斯方程，一旦求得速度勢 ， 就可由式（476）求得流速 ，解得壓強(qiáng)，問題得 到解決。 二、流函數(shù)二、流函數(shù) 對于平面運(yùn)動，有連續(xù)性微分方程 ，移項(xiàng)得 根據(jù)曲線積分定理，前式是表達(dá)式 成為某一函數(shù) 的全微分的必要和充分條件 （480） 比較 （481） 得 （482） ),(zyxuuuu x u x y u y 0 x u x y u y dxudyu yx yx, dxudyud yx

45、dydxd yx xx u yy u 函數(shù) 稱為流函數(shù)。由流函數(shù)的引出條件可知，凡是不可 壓縮流體的平面的流動，連續(xù)性微分方程成立，不論無旋流動 或有旋流動，都存在流函數(shù)，而只有無旋流動才有流速勢，可 見流函數(shù)比流速勢更具有普遍意義。 1.流函數(shù)具有以下性質(zhì)： 流函數(shù)的等值線是流線 證明： 流函數(shù)值相等 ，由式得流函數(shù)等值線方 程 則 上式即平面流動的流線方程，故 流函數(shù)的等值線是流線，給流線以不同值，便得到族流線給流 函數(shù)以不同值，便得到流線族。 .兩條流線的流函數(shù)的差值，等于通過該兩流線間的單 寬流量： yx, 0,dc 0dxudyu yx yx u dy u dx （483） 這一性質(zhì)

46、也可表述為：平面流動中，通過任一曲線的單寬 流量，等于該曲線兩端流函數(shù)的差值。 .平面無旋流動的等流函數(shù)線（流線）與等勢線正交。 證明：對于平面無旋流動，同時存 在流速勢函數(shù)和流函數(shù)，由等流函數(shù)線方程 某一點(diǎn)的斜率 由等勢線方程 dlynuxnudludq yxn ,cos,cos dluu dl dx ydl dy x dxudyu yx d 12 2 1 2 1 ddqq q q 0dxudyud yx x y u u dx dy m 1 0dyudxud yx 圖416流函數(shù) 同一點(diǎn)等勢線斜率 （484） 等流函數(shù)線與等勢線正交，故等勢線也就是過流斷面線。 .平面無旋流動，流函數(shù)是調(diào)和函

47、數(shù)。 證明：因?yàn)槠矫鏌o旋流動 則 得 帶入上式，得 （485） x y u u dx dy m 2 1 21 y x x y u u u u mm 0 2 1 y u x u x y z x u y u yx , 0 y u x u x y 0 2 2 2 2 yx 即： 平面無旋流動的流函數(shù)滿足拉普拉斯方程，是調(diào)和函數(shù)。 式中 拉樸拉斯算子 （486） 式即柯西黎曼條件。滿足拉普拉斯方程和柯西黎曼條 件，是一對共軛調(diào)和函數(shù)。 三、幾種常見的基本平面勢流：三、幾種常見的基本平面勢流： 拉普拉斯方程在復(fù)雜的邊界條件下，雖然難以求解，一 些簡單的平面勢流，其流速勢和流函數(shù)卻不難求得。研究這些 簡單

48、的平面勢流的意義在于通過簡單勢流的疊加，往往能組合 成符合某些給定邊界條件的復(fù)雜流場。 0 2 2 2 2 2 2 yx xy yx 1均勻直線流均勻直線流 均勻直線流是流場中各點(diǎn)速度大小相等，方向相同的流動，是 一種最簡單的平面勢。速度場 ， ; 速度勢 （487） （488） 若均勻直線流流速平行于軸 （489） 若均勻直線流流速平行于軸 （490） au x bu y dyudxu yx byax dxudyu yx bxay ayaxuy, 0 bxbyux, 0 圖417均勻直線流 2源流源流 如圖418所示，在平面勢流中，源流就是流體從潭點(diǎn)均勻地向各個 方向出流的流動。組成這種流型的線，就是源點(diǎn)所在平面勢流中i面 上，從源點(diǎn)0出發(fā)的一族射線。 速度場 速度 流函數(shù) 等勢線方程 等勢線是以o點(diǎn)為圓心的同心圓。 流線方程 ，流線是由o點(diǎn)引出的射線以直角坐標(biāo) 系表示。 r q r u 2 0 u rd udrur rdr q r q ln 22 drurdur 22 q r qr d crc, 圖418平面源流 cc, （491） （492） 3匯流匯流 流體從四周沿經(jīng)向均勻的流入的流動稱之為匯流。流入?yún)R 點(diǎn)單位厚度流量稱為匯流強(qiáng)度-q。匯流的速度勢和流函數(shù)的表 達(dá)式與源流的相同，符號相反