流體力學(xué)物體繞流流動

流體力學(xué)物體繞流流動第一頁，共六十一頁，編輯于2023年，星期日2023/6/11例1：空氣運動粘度大Re數(shù)流動是常見現(xiàn)象.設(shè)汽車?yán)?：水運動粘度設(shè)船第二頁，共六十一頁，編輯于2023年，星期日2023/6/12用N-S方程可以得到小雷諾數(shù)流動條件下的近似解，工程上涉及到大雷諾數(shù)流動，要尋求新的近似方法。若采用歐拉方程,同時在固體壁面上采用滑移條件(而不是無滑移粘附條件),這就是理想流體的模型。在理想流體模型的范圍內(nèi),算出的物體表面的壓力分布,在流動不分離或在接近尾緣處有小分離區(qū)的情況下與實測結(jié)果比較符合。但無法解決阻力問題。在實際流體繞流固體時，固體邊界上的流速為0，在固體邊界的外法線方向上的流體速度從0迅速增大，在邊界附近的流區(qū)存在相當(dāng)大的速度梯度，在這個流區(qū)內(nèi)粘性作用不能忽略，邊界附近的流區(qū)稱為邊界層（或附面層），邊界層外流區(qū)，粘性作用可以忽略，當(dāng)作理想流體來處理。

第三頁，共六十一頁，編輯于2023年，星期日2023/6/13如圖，平板前方均勻來流的速度v∞，從平板前緣開始形成邊界層，其厚度沿流增加。在邊界層外緣附近流速漸近于當(dāng)?shù)赝饬魉俣取ＵJ(rèn)為邊界層厚度是沿表面法線方向從到的一段距離。邊界層定義：繞流物體表面上一層厚度很小且其中的流動具有很大法向速度梯度的流動區(qū)域。第四頁，共六十一頁，編輯于2023年，星期日2023/6/14整個流場可以明顯地分成性質(zhì)很不相同的兩個區(qū)域：(1)緊貼物面非常薄的一層區(qū)域稱為邊界層。在該區(qū)域內(nèi),速度分量ux沿物面的法向變化非常迅速,它比沿切向的變化高一個數(shù)量級。即甚大。雖然在大數(shù)情況,很小,但因很大,故粘性應(yīng)力仍然可以達(dá)到很高的數(shù)值。2)邊界層外的整個流動區(qū)域稱為外部流動區(qū)域。在該區(qū)域內(nèi),很小,因此粘性應(yīng)力在大Re數(shù)情況下的確比慣性力小得多,可以將粘性力全部略去,因而把流體近似地看成是理想的。對于均勻來流繞過物體的流動而言,在整個外部流動區(qū)域中不僅可把流體視為理想的,而且可視為運動是無旋的。

注意：對于平板繞流，邊界層外緣，對于彎曲固壁，邊界層外緣。

邊界層的外邊界線與流線不重合，外流區(qū)域中的流體質(zhì)點可以連續(xù)地穿過邊界層的外緣進(jìn)入邊界層內(nèi)。第五頁，共六十一頁，編輯于2023年，星期日2023/6/15一、邊界層特點普朗特理論：邊界層內(nèi)慣性力與粘性力量級相等。邊界層很薄當(dāng)邊界層厚度增長10.1.2邊界層基本概念第六頁，共六十一頁，編輯于2023年，星期日2023/6/16二、邊界層厚度邊界層內(nèi)流態(tài)實驗測量表明邊界層內(nèi)層流態(tài)向湍流態(tài)轉(zhuǎn)捩的雷諾數(shù)為1名義厚度δ定義為速度達(dá)到外流速度99%的厚度。對平板層流邊界層第七頁，共六十一頁，編輯于2023年，星期日2023/6/172.排擠厚度（位移厚度）將由于不滑移條件造成的質(zhì)量虧損折算成無粘性流體的流量相應(yīng)的厚度δd

。又稱為質(zhì)量流量虧損厚度uouo動量損失厚度δm將由于不滑移條件造成的動量流量虧損折算成無粘性流體的動量流量相應(yīng)的厚度δm

。動量損失厚度<排擠厚度第八頁，共六十一頁，編輯于2023年，星期日2023/6/1810.2平板邊界層流動10.2.1普郎特邊界層方程第九頁，共六十一頁，編輯于2023年，星期日2023/6/19今以來流速度為特征速度,平板長度為特征長度。引進(jìn)如下無量綱量

無量綱的基本方程組：

第十頁，共六十一頁，編輯于2023年，星期日2023/6/110在大數(shù)情況下邊界層內(nèi)的流動有如下兩個主要性質(zhì)：(1)邊界層的厚度δ較物體的特征長度L小得多,即是一個小量。(2)在邊界層內(nèi)粘性力和慣性力同數(shù)量級。對“量級”的兩點說明：(1)估計量級必須有個標(biāo)準(zhǔn)。在邊界層問題中,我們?nèi)ˇ?作為估計量級的標(biāo)準(zhǔn)。并采用符號Ο,例如Ο（δ*）表示和δ*同量級。(2)所謂量級不是指該物理量或幾何量的具體數(shù)值,而是指該量在整個區(qū)域內(nèi)相對于標(biāo)準(zhǔn)小參數(shù)(在邊界層問題中,即相對于δ*)而言的平均水平。所以允許Ο(1)的量(即與1同量級的量)甚至比1大一個量級的量在個別點上或局部小區(qū)域內(nèi)取較小的數(shù)值,甚至等于零,重要的是它的平均水平是與1同量級,或比1大一個量級就行了。

第十一頁，共六十一頁，編輯于2023年，星期日2023/6/1111.及其各階導(dǎo)數(shù)的量級

2.及其各階導(dǎo)數(shù)的量級

由連續(xù)方程得

積分得第十二頁，共六十一頁，編輯于2023年，星期日2023/6/1123.及的量級

壓力梯度是起調(diào)節(jié)作用的被動的力。它們的量級由方程中其它類型力中的最大量級所決定

第十三頁，共六十一頁，編輯于2023年，星期日2023/6/113第十四頁，共六十一頁，編輯于2023年，星期日2023/6/114無量綱形式的普朗特層流邊界層方程

第十五頁，共六十一頁，編輯于2023年，星期日2023/6/115有量綱形式的普朗特層流邊界層方程為：

邊界條件

①在物面上,②在邊界層外邊界上,

第十六頁，共六十一頁，編輯于2023年，星期日2023/6/116內(nèi)蒙古工業(yè)大學(xué)工程流體力學(xué)電子課件10.2.2平板層流邊界層的精確解因為平板沒有厚度,當(dāng)理想流體沿平板方向流過平板時,平板對流動沒有擾動,因此外流的速度場是均勻的且等于常數(shù)。根據(jù)伯努利方程,壓力也均勻

常數(shù),

第十七頁，共六十一頁，編輯于2023年，星期日2023/6/117普朗特邊界層方程式

（10-10）邊界條件是：x≥0，

第十八頁，共六十一頁，編輯于2023年，星期日2023/6/118說明：①第三式表明邊界層內(nèi)y方向壓強梯度為零，說明外部壓強可穿透邊界層直接作用在平板上。外部壓強由勢流決定②第二式右邊得到簡化（x方向二階偏導(dǎo)數(shù)消失），有利于數(shù)值計算。利用該方程就可計算壁切應(yīng)力和流動阻力，具有里程碑式意義。第十九頁，共六十一頁，編輯于2023年，星期日2023/6/119無量綱化方程和邊界條件為：

第二十頁，共六十一頁，編輯于2023年，星期日2023/6/120返回有量綱形式為

由于討論的是半無限長平板上邊界層內(nèi)的流動,不存在特征長度,因此在最后的解中不應(yīng)該出現(xiàn)特征長度,即和應(yīng)該與無關(guān)。這就要求自變量以下列組合第二十一頁，共六十一頁，編輯于2023年，星期日2023/6/121第二十二頁，共六十一頁，編輯于2023年，星期日2023/6/122這種原有兩個自變量的偏微分方程組,若其解只依賴于一個組合變量,能使偏微分方程變?yōu)槌Ｎ⒎址匠?則稱此方程式具有相似性解。ξ（η）所滿足的方程式是一個非線性的三階常微分方程,形式雖然簡單,但卻無法找出封閉形式的解析解來。

布拉修斯于1908年用級數(shù)銜接法求出此問題的近似解,而后托普費爾(Topfer)、哥德斯坦(Goldstein)、霍華斯(Howarth)、哈脫利(Hartree)等人分別用數(shù)值方法求出了此問題的解。

第二十三頁，共六十一頁，編輯于2023年，星期日2023/6/123由數(shù)值解繪制的無量綱速度廓線與尼古拉茲實驗測量結(jié)果吻合。布拉修斯解邊界層厚度排擠厚度動量損失厚度第二十四頁，共六十一頁，編輯于2023年，星期日2023/6/124布拉修斯的精確解有兩個重要的前提條件：邊界層內(nèi)的流動是層流；沿平板的壓力梯度為零。工程中常遇到的問題要復(fù)雜些，如邊界層是湍流，沿壁面的壓力梯度不為零等等，這就促使人們?nèi)で蠼平夥?。在近似解法中．?yīng)用最廣泛的是馮·卡門的邊界層近似積分法。

10.2.3

邊界層動量方程

第二十五頁，共六十一頁，編輯于2023年，星期日2023/6/125

流體繞流中作用在物體上的力可以分為垂直于來流方向的升力和平行于來流方向的阻力，繞流阻力可以分成摩擦阻力與形狀阻力，都與邊界層有關(guān)。繞流阻力作用表現(xiàn)在于邊界層內(nèi)流速的降低，引起動量的變化。通過建立邊界層的動量方程來研究摩擦阻力。沿物體的曲面取x軸，沿物體表面法線取y軸，在物體表面取邊界層微元段ABCD，把它放大，x軸便成為直線，線段BD長為dx，AC為邊界層外邊界，AB、CD垂直于物體表面。第二十六頁，共六十一頁，編輯于2023年，星期日2023/6/126假設(shè)：①

不計質(zhì)量力②

流動為定常流動③

dx無限小，BD、AC可看成直線由動量方程由控制面AB沿x方向流入動量（1）由控制面CD沿x方向流出動量

由控制面AC沿x方向流入動量

第二十七頁，共六十一頁，編輯于2023年，星期日2023/6/127因為，所以邊界層內(nèi)邊界就是物體表面，其流速為0，其壓強等于邊界層外邊界的壓強，即沿物體表面的法線y方向壓強不變，p與y無關(guān)，可用全微分代替偏微分，上式可寫作（5）將（2）、（3）、（4）、（5）代入（1）得到（6）方程（6）就是邊界層積分方程，由馮·卡門首先推導(dǎo)出來的，稱作卡門動量積分方程。第二十八頁，共六十一頁，編輯于2023年，星期日2023/6/128討論：（1）如果dp/dx=0，um為常數(shù)（2）動量損失厚度第二十九頁，共六十一頁，編輯于2023年，星期日2023/6/12910.2.4平板層流邊界層的近似解繞平板的穩(wěn)態(tài)不可壓縮流動：①x方向的壓力梯度dp／dx取邊界層外邊界上外流(即勢流)的只值；②對于平扳來說，邊界層外邊界上的速度等于來流速度uo；③假定一個合理的速度分布．通常認(rèn)為沿y方向的速度分布曲線在任意x處類似，即認(rèn)為u是y／δ的函數(shù)；④根據(jù)牛頓剪切定律，用

代替τ。

第三十頁，共六十一頁，編輯于2023年，星期日2023/6/1301.速度剖面族的選取 我們選取如下三次多項式其中依據(jù)條件,可得：

由條件,可得由條件,可得設(shè)定平板上為層流邊界層，首先補充邊界層流速分布關(guān)系式，第三十一頁，共六十一頁，編輯于2023年，星期日2023/6/131由條件,可得

選定的速度剖面為

2.單參數(shù)δ（x）的確定第三十二頁，共六十一頁，編輯于2023年，星期日2023/6/132第三十三頁，共六十一頁，編輯于2023年，星期日2023/6/133積分得（）3壁面摩擦阻力第三十四頁，共六十一頁，編輯于2023年，星期日2023/6/134局部阻力系數(shù)布拉修斯的精確解則對長為寬為且兩邊浸沒在流體中的平板所受到的總摩擦阻力為：第三十五頁，共六十一頁，編輯于2023年，星期日2023/6/135總阻力系數(shù)為（）10.2.5平板湍流邊界層的近似解為了方便,在工程上往往采用冪次公式作為近似速度剖面族。例如采用經(jīng)驗公式(以下表示時均值的上標(biāo)“一”略去,并用代替)：

一、湍流邊界層內(nèi)速度剖面的選取

第三十六頁，共六十一頁，編輯于2023年，星期日2023/6/136當(dāng)在5×105到107范圍內(nèi),隨著ReL增加,1/n下降

如取在湍流外邊界處（*）第三十七頁，共六十一頁，編輯于2023年，星期日2023/6/137式中δ（x）是依賴于x的未知函數(shù)δ=f（x）。為了完全確定速度剖面,還需要確定δ（x）。

為了確定δ（x）,需要應(yīng)用式

二、單參數(shù)δ（x）的確定

由（*）式得又因第三十八頁，共六十一頁，編輯于2023年，星期日2023/6/138所以將以上代入卡門積分動量方程式，得通過采用近似方法,認(rèn)為湍流邊界層從平板前緣就已形成,即認(rèn)為：第三十九頁，共六十一頁，編輯于2023年，星期日2023/6/139 三、摩擦阻力

局部阻力系數(shù)長為L,寬為b的平板受到的總阻力系數(shù)為：

同實驗結(jié)果比較表明,若把上式中的0.072改為0.074,即

適用范圍第四十頁，共六十一頁，編輯于2023年，星期日2023/6/140當(dāng)時,式（10-51）就不太準(zhǔn)確了。需要采用對數(shù)速度剖面。在的范圍內(nèi),史里希丁用對數(shù)速度分布與積分關(guān)系式聯(lián)合求解得出的摩擦阻力系數(shù)公式為：它表明阻力和來流速度的1.8次方成正比,而層流時阻力和來流速度的1.5次方成正比,故湍流邊界層的摩阻比層流邊界層的摩阻大。邊界層內(nèi)流動形態(tài)轉(zhuǎn)變的典型情況表示在圖中。：

第四十一頁，共六十一頁，編輯于2023年，星期日2023/6/141四、平板混合邊界層的近似計算前面假定整個平板上是層流或紊流邊界層，實際上，當(dāng)Re增大到一定數(shù)值時，平板長度達(dá)到一定長度，即

L＞xer時，平板前部是層流邊界層，后部是紊流邊界層，中間有一過渡段，這種邊界層稱為混合邊界層。計算時引入假設(shè)：（1）層流邊界層轉(zhuǎn)變?yōu)槲闪鬟吔鐚邮窃谔幫蝗话l(fā)生，無過渡段；（2）混合邊界層的紊流邊界層可以看作是從平板的首端開始的紊流邊界層的一部分。普朗特建議：在邊界層轉(zhuǎn)捩位置以前采用層流的摩擦阻力系數(shù),在其后采用的摩擦阻力系數(shù),于是混合邊界層的總的阻力系數(shù)為第四十二頁，共六十一頁，編輯于2023年，星期日2023/6/142其中第四十三頁，共六十一頁，編輯于2023年，星期日2023/6/143近似計算方法第四十四頁，共六十一頁，編輯于2023年，星期日2023/6/14410.3邊界層分離與壓差阻力

10.3.1邊界層分離現(xiàn)象

流體繞過非線型鈍頭物體時，較早脫離物體表面，在物體后部形成較寬闊的尾流區(qū)，在邊界層內(nèi)，流體質(zhì)點在某些情況下向邊界層外流動的現(xiàn)象稱為邊界層從固體分離。

以圓柱繞流為例，虛線為邊界層外邊界。第四十五頁，共六十一頁，編輯于2023年，星期日2023/6/145注意：C點的位置，這是由于在加速減壓和減速增壓的過程中，還存在克服流動阻力所消耗的能量損失

由伯努利方程知，愈靠近圓柱，流速越小，壓強越大，在貼近圓柱面A處流速為0，壓強最大，A點稱為駐點。由于液體不可壓縮，繼續(xù)流來的液體質(zhì)點在駐點的壓強的作用下，將壓能轉(zhuǎn)化為動能，從而改變流向，沿圓柱面兩側(cè)繼續(xù)向前流動。由于圓柱面的阻滯作用，在表面產(chǎn)生邊界層，從A點經(jīng)1/4圓周到B點之前，柱面向外凸出，流線趨于密集，邊界層內(nèi)流體處在加速減壓情況，，這時由于壓能減小部分還能夠補償動能增加和由于克服流動阻力而消耗的能量損失，因此此時B點處邊界層內(nèi)流體質(zhì)點速度不為0。第四十六頁，共六十一頁，編輯于2023年，星期日2023/6/146

過B點之后，流線逐漸疏散，邊界層內(nèi)流體處于減速增壓的情況，動能轉(zhuǎn)化成壓能，同時也用以克服流動阻力而消耗的能量。在C點處邊界層內(nèi)流體質(zhì)點速度下降為0。流體質(zhì)點在C點停滯下來，形成新的停滯點，繼續(xù)流來的流體質(zhì)點將脫離原來的流線，沿另一流線CE流去，從而使邊界層脫離了圓柱面，這樣就形成了邊界層的分離現(xiàn)象，C點為分離點。分離點的位置與繞流物的形狀、粗糙程度、流動的Re數(shù)和來流與物體的相對方向有關(guān)。邊界層分離后，邊界層和圓柱面之間，由于分離點下游壓強大，從而使流體發(fā)生反向回流，形成旋渦區(qū)。第四十七頁，共六十一頁，編輯于2023年，星期日2023/6/147邊界層分離：邊界層脫離壁面

邊界層分離2.分離的原因—

粘性圓柱后部：貓眼1.分離現(xiàn)象在順壓梯度區(qū)（BC)：流體加速在逆壓梯度區(qū)（CE)：CS段減速S點停止SE段倒流。3.分離的條件—

逆壓梯度4.分離的實際發(fā)生—

微團(tuán)滯止和倒流第四十八頁，共六十一頁，編輯于2023年，星期日2023/6/148分離實例從靜止開始邊界層發(fā)展情況擴(kuò)張管（上壁有抽吸）

邊界層分離第四十九頁，共六十一頁，編輯于2023年，星期日2023/6/14910.3.2壓差阻力CD=CDf+CDp形狀阻力（壓差阻力）：粘性流體繞流時，在物體表面上所作用的壓力的合力在流動方向上的投影。對非流線型物體，是由于邊界層的分離，在物體尾部形成旋渦，旋渦區(qū)的壓強較物體前部低，在流動方向上產(chǎn)生了壓強差，形成了作用于物體上的阻力，稱為壓差阻力。壓差阻力主要取決于物體的形狀。第五十頁，共六十一頁，編輯于2023年，星期日2023/6/150一摩擦阻力是由于流體的粘性引起的，當(dāng)流體繞流物體時，在表面上形成了邊界層，邊界層內(nèi)速度梯度大，粘性的牽制作用使物體受到阻力。阻力發(fā)生在運動物體表面上。

摩擦阻力特點阻力系數(shù)強烈地依賴于雷諾數(shù)；2)對相同雷諾數(shù)，層流態(tài)的阻力明顯低于湍流態(tài)；3)對湍流邊界層，光滑壁面的阻力最小，粗糙度增加使阻力系數(shù)增大；4)摩擦阻力與壁面面積成正比。第五十一頁，共六十一頁，編輯于2023年，星期日2023/6/151二壓差阻力

與邊界層的分離現(xiàn)象密切相關(guān)。當(dāng)流體流過一個圓頭尖尾的回轉(zhuǎn)體時，在物體前端形成減速區(qū)，在前端頂點A形成駐點，流體壓強隨流速變化而變化，在駐點處最大，離開駐點，壓強逐漸減小，從B點處開始變成負(fù)值，過最大速度點C后，流速減小，壓強上升，壓強又變成正值。第五十二頁，共六十一頁，編輯于2023年，星期日2023/6/152壓強分布如實線所示，虛線理想壓強分布。從圖中可以看出，前端的正壓強產(chǎn)生一個向后的水平合力，后端的正壓強產(chǎn)生一個向前的水平合力，中段壓強為負(fù)值，產(chǎn)生吸力，其前半部合成一向前的水平力，后半部合成一向后的水平力，這兩者數(shù)值相差不大，幾乎相互抵消。因此，物體所受的水平合力取決于前端正壓強造成的向后的較大的力與后端正壓強造成的向前的較小的力，相互抵消后，還剩下向后的反物體前進(jìn)的力，即壓差阻力。物體形狀→后部逆壓梯度→壓強分布→壓強合力用實驗方法確定形狀阻力→阻力曲線第五十三頁，共六十一頁，編輯于2023年，星期日2023/6/153內(nèi)蒙古工業(yè)大學(xué)工程流體力學(xué)電子課件第五十四頁，共六十一頁，編輯于2023年，星期日2023/6/154內(nèi)蒙古工業(yè)