1、工程流體力學,第八章 粘性流體繞物體的流動,第八章 粘性流體繞物體的流動,實際流動都是有粘流動，目前對粘性流動研究方法主要有：1、基于N-S方程的紊 流模擬2、流體實驗,流動分類,根據工程的實際情況，流動可分為： 內流和外流。 內流 ：如右上圖。 外流： 如右下圖。,本章的主要內容,本章主要討論繞流問題，即外流問題。首先將介紹粘性流體的運動微分方程，然后將給出邊界層的概念及其控制方程

2、，最后針對繞流流動現象的一些具體問題進行了討論。 ◆空間流動三維問題，N—S方程及其求解 ◆擾流阻力及其計算 ◆附面層的問題,,,第一節(jié) 不可壓縮粘性流體的運動微分方程,以流體微元為分析對象，流體的運動方程可寫為如下的矢量形式： 這里 ： 是流體微團的加速度，微分符號： 稱為物質導數或隨體導數，它表示流體微團的某性質 時間的變化率。,,,,

3、(8-1),(8-2),(8-3),一、微元體的受力分析和運動微分方程的推導,如圖所示，控制體的各邊長分別為dx,dy,dz，微元體的體積為： （ 8－4） 作用在微元體上的質量力為

4、 ，其可用 三個分量 表示為： （8－5）這里： （8－6）如果的三個分量是 ，則：

5、 （8－7）,,,,,,★作用在微元體上的表面力,將微元體六個面上的應力分別投影到三個坐標方向上如圖,◇作用于微元體個面上的x軸方向的應力,把作用于控制體上x方向的力疊加起來，得到作用在微元體上的表面力在x方向的分量為：,,◇作用于微元體個面上的Y、Z軸方向的應力,同理，表面力在y方向的分量為：表面力在z方向的分量為：,,,★作用在微元體上的表面力,如果用 ,

6、和 表示單位體積的表面力，則： （ 8－8）,,,,,★作用在微元體上的表面力,將上式和式（8－7）代入式（8－1）則得： （8－9） 這就是微分形式的運動方

7、程。,二、本構方程,本構方程是確立應力和應變率之間關系的方程式。斯托克斯通過將牛頓內摩擦定律推廣到了粘性流體的任意流動中，建立了牛頓流體的本構方程： （8－10） 上式也稱為廣義牛頓定律,,三、納維－斯托克斯方程（簡稱N－S方程）,將式（8－10）代入式（8－9）可得：

8、

9、 （ 8－11） 上式稱納維－斯托克斯（Naver-Stokes）方程，是粘性流體運動微

10、分方程的又一種形式。,,對于不可壓流體，其連續(xù)方程為：對于不可壓縮粘性流體，粘性體膨脹應力為零，其運動方程為： （8－12）,,,三、納維－斯托克斯方程（簡稱N－S方程）,●并考慮到拉普拉斯算子： 不可壓縮粘性流體的運動方

11、程還可寫為： （8－13）,,,三、納維－斯托克斯方程（簡稱N－S方程）,如果質量力只有重力作用，用 代表重力加速度，不可壓縮粘性流體的運動方程的矢量形式為：

12、 (8-14) 右端第一項表示單位質量的質量力；第二項代表作用于單位質量流體的壓強梯度力；第三項代表黏性變形應力。,,,三、納維－斯托克斯方程（簡稱N－S方程）,對理想流動，認為流體無粘性， ，這時運動方程簡化為歐拉方程：

13、 （8－15） 或矢量形式 （8－16）,,,,三、納維－斯托克斯方程（簡稱N－S方程）,●當流體靜止不動時， ，則運動方程簡化為：

14、 （8－17）,,,三、納維－斯托克斯方程（簡稱N－S方程）,第二節(jié) 蠕動流動,蠕動流動：雷諾數很低的流動。特點：流動的尺度和流動的速度均很小如：熱電廠鍋爐爐膛氣流中繞煤粉顆粒、 油滴等的流動；滑動軸承間隙中的

15、流 動等等。,一、蠕動流動的微分方程,對于定常流動，忽略慣性力和質量力，在直角坐標系下，可把納維爾――斯托克斯方程（8－14）組簡化成 ： （8－18）,,一、蠕動流動的微分方程,●如果流動是不可壓縮流體，則連續(xù)性方程為：

16、 （8－19） 將式（8－18）依次求 、 、 ，然后相加，并結合連續(xù)性方程，即得： 即蠕動流動的壓力場滿足拉普拉

17、斯方程。,,,,,,(8- 20),二、繞球的蠕動流動,對如圖所示的無窮遠來流以速度 均勻平行流沿 軸繞半徑為 的靜止圓球流動，得速度與壓 強分布為： （8－21）,,,,,二、繞球的蠕動流動,式中 為無窮遠處來流的壓力。 圓球以

18、很小的速度在靜止流體中作等速運動時，在流場中通過x軸的平面上的流譜如圖所示。,,二、繞球的蠕動流動,在圓球的前后兩駐點A和B處的壓強是壓強的最高點和最低點，分別為：在前駐點A（ ＝180° ） （8－22） 在后駐點B（ ＝0°）：

19、 （8－23）而切應力的最大值，發(fā)生在C（ ＝90°）為： （8－24） 等于

20、A、B點處的壓強與無窮遠處的壓強之差的絕對值。,,,,,,二、繞球的蠕動流動,球面上的壓強和剪切應力也可根據速度分布公式算出，為： （8-25） 對上述兩式積分，可分別得到作用在球面上的壓強和切應力的合力。將這兩個合力在流動方

21、向的分量相加，可得到流體作用在圓球上的阻力為： （8-26） 這就是圓球的斯托克斯阻力公式。式中d=2 為圓球的直徑。,,,,第三節(jié) 邊界層的概念,邊界層：物體壁面附近存在大的速度梯度的薄層。 我們可以用如圖

22、所示的繞平板的流動情況說明邊界層的概念。,,★邊界層的定義,粘性流體繞流物體時，由于粘性的作用，在物體的表面附近，存在一速度急劇變化的薄層——邊界層。 例如：來流 的流體繞流平板時，在平板表面形成邊界層。,,,,,★邊界層的定義,在平板的前部邊界層呈層流狀態(tài)，隨著流程的增加，邊界層的厚度也在增加，層流變?yōu)椴环€(wěn)定狀態(tài)，流體的質點運動變得不規(guī)則，最終發(fā)展為紊流，這一變化發(fā)生在一段很短的長度范圍，稱之為轉捩區(qū)，轉類區(qū)的開始點

23、稱為轉捩點。轉類區(qū)下游邊界層內的流動為紊流狀態(tài)。在轉捩區(qū)和紊流區(qū)的壁面附近，由于流體的質點的隨機脈動受到平板壁面的限制，因此在靠近壁面的更薄的區(qū)域內，流動仍保持為層流狀態(tài)，稱為層流底層和粘性底層。,◆邊界層的特點,邊界層內速度梯度很大，旋渦強度大，有旋流動慣性力和粘性具有相同的數量級，同時考慮。邊界層外部速度梯度很小，可以作為理想流體的勢流處理。邊界層厚度隨 的增大而增大，隨 的增大而減小。由于邊界層很薄，因而可以近似

24、認為，邊界層任一截面上各點壓強相等。,◆邊界層的分類,按流動狀態(tài)，可分為層流邊界層和紊流邊界層?！衽袆e準則——雷諾準則： 平板上的臨界雷諾數 = ~ ●邊界層的構成： 1.層流邊界層，當 較小時，邊界層內全為層流，稱為層流邊界層。 2.混合邊界層：除前部起始部分有一小片層流區(qū)，其余大部分為紊流區(qū)，稱為混合邊界層。,◆邊界層的厚度,兩個

25、流動區(qū)域之間并沒有明顯的分界線。邊界層的厚度：通常，取壁面到沿壁面外法線上速度達到勢流區(qū)速度的99％處的距離作為邊界層的厚度，以δ表示，這一厚度也稱邊界層的名義厚度。邊界層的厚度取決于慣性和粘性作用之間的關系，即取決于雷諾數的大小。雷諾數越大，邊界層就越?。环粗?，隨著粘性作用的增長，邊界層就變厚。沿著流動方向由繞流物體的前緣點開始，邊界層逐漸變厚。,第四節(jié) 平面層流邊界層的微分方程,在這一節(jié)里，將利用邊界層流動的特點如流體的粘度大小

26、、速度與溫度梯度大和邊界層的厚度與物體的特征長度相比為一小量等對N-S方程進行簡化從而導出層流邊界層微分方程。在簡化過程中，假定流動為二維不可壓定常流，不考慮質量力，則流動的控制方程N-S方程為： （8-27）,,第四節(jié) 平面層流邊界層的微分方程,將上述方程組無量綱化。為此考慮如圖所示的一半無窮繞流平板，假定無

27、窮遠來流 的速度 ，流動繞過平板時在平板附近形成邊界層，其厚度為 ，平板前緣至某點的距離為 。取 和 為特征量，可定義如下 的無量綱量： / / / /

28、 /（ ）,,,,,,,,,,,,,,,,,,,,,,,,代入方程組（8－27），整理后得： （8-28）式中雷諾數,,,第四節(jié) 平面層流邊界層的微分方程,與 相比較是很小的 ，即 << 或 /

29、<< 1，同時注意到， 與 、 與 、 與 具有同一數量級，于是 、 、 和 的量級均為1，并可以得到： ～1 ～1 ～ 1 ～ 為了估計其他各量的數量級，由連續(xù)性方程可得： ＝ ～1,,,,,,,,,,,,,,,,,,,,

30、,,第四節(jié) 平面層流邊界層的微分方程,第四節(jié) 平面層流邊界層的微分方程,因此 ~ ，于是又得到： ~ ~ ~ 1 ~ 通過分析方程組（8－28）各項的數量級，方程組（8－28）中第二式中各慣性項可以忽略掉 ，同時可以略 去 、 、 。于是在方程組（8－

31、28）的粘性 項中只剩第一式中的一項 。,,,,,,,,,,,,如果僅保留數量級為1的項，而將數量級比1小的各項全部略去，再恢復到有量綱的形式，便可以得到層流邊界層的微分方程組為： （8-29） 沿邊界層上緣由伯努利可知：

32、 常數 上式對 求導，得：,,,,,第四節(jié) 平面層流邊界層的微分方程,這樣，層流邊界層的微分方程又可寫為： （8-30） 方程組（8－30） 即為在物體壁面為平面的假設下得到的邊界層微分方程 。,,第四節(jié) 平面層流

33、邊界層的微分方程,第五節(jié) 邊界層的動量積分關系式,邊界層的動量積分方程是對邊界層內流動的再簡化。其推導過程有兩種方法：一種是沿邊界層厚度方向積分邊界層的方程組，一種是在邊界層內直 接應用動量守恒原理。下面的推導采用第二種方法。,◆邊界層動量積分方程的推導,如圖所示為不可壓縮流體的定常二維邊界層流動 ，設物體表面型線的曲率很小。 取一個單位厚度的微小控制體，它的投影面ABDC 。用動量定理來建立該控制體內的流體在單位時間內沿x方向

34、的動量變化和外力之間的關系。,,◆邊界層動量積分方程的推導,設壁面上的摩擦應力為 根據邊界層的控制方程組，邊界層內的壓強僅近似地依賴于 而與 無關，設AB面上的壓強為 ，DC上的壓強為 控制面AC為邊界層的外邊界 其外部為理想流體的勢流 ，只有與之垂直的壓力 ，設AC上的壓強為A，C兩點壓強的平均值 。作用在控制體上的表面力沿方向的合力為：,◆邊界層動量積分方程的推導,式中為邊界層外邊界AC與方

35、向的夾角，由幾何關系可知： ，上式經整理并略去高階小量，得：單位時間內沿方向經過AB流入控制體的質量和動量分別為：經過CD面流出的質量和動量分別為：定常流動條件下，可知從控制面AC流入控制體中的流量為：由此引起流入的動量為：,◆邊界層動量積分方程的推導,式中V為邊界層外邊界上的速度。這樣，可得單位時間內該控制體內沿x方向的動量 變化為 根據動量定理， ，則可得邊界層的動量積分方程

36、為： （8-51） 上式也稱為卡門動量積分關系式。該式是針對邊界層流動在二維定常流動條件下導出的，并沒有涉及邊界層的流態(tài)，所以其對層流和紊流邊界層都能適用。,◆積分方程的求解,實際上可以把 、 和 看作已知數，而未知數只有 、 和 三個。

37、再補充兩個關系式： 一、沿邊界層厚度的速度分布 = (y) 二、切向應力與邊界層厚度的關系式 一般在應用邊界層的動量積分關系式（8－51）來求解邊界層問題時，邊界層內的速度分布是按照已有的經驗來假定的。假定的 愈接近實際，則所得到的結果愈正確。所以選擇邊界層內的速度分布函數 是求解邊界層問題的重要關鍵。,,,,,,,,,,,,第六

38、節(jié)邊界層的位移厚度和動量損失厚度,邊界層的厚度 ，表示粘性影響的范圍。 位移厚度 動量損失厚度根據伯努力方程可知：又由于：帶入（8-51）得 或 （8-52）,◆邊界層厚度計算式的推導,因此在邊

39、界層內由于粘性影響使體積流量的減小量 ，即上式中第一項積分。 位移厚度或排擠厚度 可表示成： （8-53）同理動量損失厚度 可表示為：

40、 （8-54） 將 和 代入式（8－51）, 得

41、 （8-55）,,◆邊界層厚度計算式的推導,式（8-55）是另一種形式的平面不可壓縮粘性流體邊界層動量積分關系式 。 、 和 都是未知數，它們決定于邊界層內速度的分布規(guī)律。 將式（8－55）化為無因次形式，統(tǒng)除以 ，得 （8－5

42、6） 或式中H＝ 。計算曲面邊界層時，用上式較為方便。,第七節(jié) 平板邊界層流動的近似計算,平板層流邊界層的近似計算 對于式（8－51），如果邊界層外部的壓強梯度為零，方程變?yōu)椋?(8-57） 假定平板非常薄，對流動沒有

43、影響。邊界層外層流動： 則上式可變?yōu)椋?(8-58） 兩個補充關系式:一、馮卡門假定，二、牛頓內摩擦定律。平板紊流邊界層的近似計算 采用將邊界層內的速度分布與圓管內充分發(fā)展紊流的速度分布規(guī) 律進行類比的方法。,◆平板層流邊界層的近似計算,選擇

44、一三次項式速度分布： (8-59) 根據下列邊界條件來確定待定系數 和 . (1)在平板壁面上的速度為零，即在 處 (2)在邊界層外邊界上的速度等于來流速度 ,即在 處

45、 ， (3)在邊界層外邊界上，摩擦切應力 為零，即在 處 ， (4)由于在平板壁面上的速度為零，即 ，由方程組（8－50）的第一式得,◆平板層流邊界層的近似計算,速度分布的四個系數可確定為： 于是，層流邊界層中速度的分布規(guī)律為

46、

47、 (8-60) 第二個補充關系式：利用牛頓內摩擦定律和式（8－60）得出

48、 (8-61)式中為動力粘性系數。將速度分布方程（8－60）帶入方程（8－61）并積分得：分離變量，并積分得： (8-62),◆平板層流邊界層的近似計算,式中為 運動粘性系數，為基于長度的雷諾數 。

49、合并方程（8－62）和(8－61)得到： (8-63)如果表面摩擦系數 為： (8-

50、64) 那么 ，為： (8-65) 根據動量損失厚度的定義式（8－54），并考慮式（8－62），可得動量損失厚度為：

51、 (8-66) 同理，位移厚度為: (8-67)上述計算

52、結果是依賴于所假設的速度分布規(guī)律的，不同階次的速度分布，可以得出不同的結果。表8.1 給出幾種不同的情況。,表8.1不同階次的速度分布所得結果比較,,二、平板紊流邊界層的近似計算,如前所述由于流動的混參以及速度和壓力的波動，紊流邊界層的速度分布都采用一些模型假定。普朗特建議，當邊界層雷諾數 時，邊界層內的速度分布可采用 次方規(guī)律，即：

53、 （8-68） 該式不能直接應用于邊界層的內邊界。通常認為粘性底層內的速度分布為線形分布。 雷諾數取 時的摩擦阻力系數為： 當時 普朗特和施利希廷 （ H. Schlichting）采用對數速度分布，得到如下的半經驗公式：,層流與紊流邊界層的

54、近似計算公式匯總,平板的層流邊界層和紊流邊界層的重大差別有：紊流邊界層內沿平板壁面發(fā)向截面上的速度比層流邊界層的速度增加得快 沿平板壁面紊流邊界層的厚度比層流邊界層的厚度增加得快 在其它條件相同的情況下，平板壁面上的切向應力 沿著壁面的減小在紊流邊界層中要比層流邊界層減小得慢。在同一 下，紊流邊界層得摩擦阻力系數比層流邊界層的大得多 實際情況下，邊界層是層流和紊流同時存在的混合邊界層,,,`,層流與紊流邊界層的近

55、似計算公式匯總,,第八節(jié) 邊界層的分離與卡門渦街,一、邊界層的分離,以如圖所示的圓柱繞流為例 在勢流流動中流體質點從D到E是加速的，為順壓強梯度;從E到F則是減速的, 為逆壓強梯度流體質點由D到E過程，由于流體壓能向動能的轉變，不發(fā)生邊界層分離E到F 段動能只存在損耗，速度減小很快，在S點處出現粘滯 ，由于壓力的升高產生回流導致邊界層分離，并形成尾渦。如圖為邊界層分離示意圖。,,從以上的分析中可得如下結論：粘性流體在壓

56、力降低區(qū)內流動（加速流動），決不會出現邊界層的分離，只有在壓力升高區(qū)內流動（減速流動），才有可能出現分離，形成漩渦。尤其是在主流減速足夠大的情況下，邊界層的分離就一定會發(fā)生。,圖8－14 邊界層分離示意圖,,一、邊界層的分離,二、卡 門 渦 街,上圖表示 不同雷諾數條件下繞圓柱的流動圖譜 討論圓柱繞流問題：隨著雷諾數的增大邊界層首先出現分離，分離點并不斷的前移，當雷諾數大到一定程度時，會形成兩列幾乎穩(wěn)定的、非對稱性的、交替脫落的、

57、旋轉方向相反的旋渦，并隨主流向下游運動，這就是卡門渦街 ，如右圖。卡門對渦街進行運動分析得出了阻力、渦釋放頻率以及斯特羅哈數的經驗公式?？ㄩT渦街會產生共振，危害很大；也可應用于流量測量。,,,第九節(jié) 物體的阻力與減阻,物體繞流時會受到升力和阻力的作用。物體阻力包括摩擦阻力和壓差阻力。摩擦阻力與物體表面積大小有關，壓差阻力與物體的形狀有關系。物體的阻力目前都是用實驗測得。,,激波阻力,◆物體阻力的減小辦法,減小摩擦阻力：可以

58、使層流邊界層盡可能的長，即層紊流轉變點盡可能向后推移，計算合理的最小壓力點的位置。在航空工業(yè)上采用一種“層流型”的翼型 ，便是將最小壓力點向后移動來減阻，并要求翼型表面的光滑程度。減小壓差阻力：使用翼型使得后面的“尾渦區(qū)”盡可能小。也就是使邊界層的分離點盡可能向后推移 。例如采用流線性物體就可以達到這樣的目的。工程上習慣用無因次的阻力系數 來代替阻力 (8-85）,,,,,

59、,,◆物體阻力的大小與雷諾數的關系,按相似定律可知，對于不同的不可壓縮流體中的幾何相似的物體，如果雷諾數相同，則它們的阻力系數也相同 在不可壓縮粘性流體中，對于與來流方向具有相同方位角的幾何相似體，其阻力系數：,,繞圓柱流動的阻力系數與雷諾數的關系,典型物體的阻力系數,圓 柱,半 管,半 管,方 柱,平 板,橢 柱,橢 柱,球,半 球,半 球,方 塊,方 塊,矩 形 板