高等流體力學(xué)第三講課件

第三講大雷諾數(shù)流體力學(xué)問(wèn)題的求解

本講主要內(nèi)容

（以平面流動(dòng)為例） 一、問(wèn)題概述 二、歐拉方程的求解（有勢(shì)流動(dòng)的基本理論） 三、邊界層內(nèi)運(yùn)動(dòng)的解析解法北京工業(yè)大學(xué)市政學(xué)科部——馬長(zhǎng)明高等流體(水）力學(xué)講稿第三講大雷諾數(shù)流體力學(xué)問(wèn)題的求解

一、問(wèn)題概述

1、大雷諾數(shù)近似下的歐拉方程（Euler’sEquations）

N—S方程，用相對(duì)值表示為無(wú)量綱形式，其中參考量選擇如下：

理論：在20世紀(jì)前已比較完善；成果：流場(chǎng)，升力。 問(wèn)題：阻力（達(dá)朗伯佯謬D‘Alembert’sParadox）；無(wú)滑動(dòng)條件。長(zhǎng)度：速度：壓強(qiáng)：時(shí)間：當(dāng)：時(shí)，有：歐拉方程：（Hydrodynamics）北京工業(yè)大學(xué)市政學(xué)科部——馬長(zhǎng)明高等流體(水）力學(xué)講稿第三講大雷諾數(shù)流體力學(xué)問(wèn)題的求解

2、邊界層內(nèi)流動(dòng)的量級(jí)分析—邊界層方程 （1）邊界層的提出 通過(guò)流場(chǎng)顯示結(jié)果，Prandtl（1904）提出邊界層理論。 邊界層外，流動(dòng)有勢(shì)，滿足歐拉方程； 邊界層內(nèi)，粘性運(yùn)動(dòng)，剪切層產(chǎn)生旋渦； 漩渦的粘性擴(kuò)散與對(duì)流傳輸—漩渦僅限于邊界層內(nèi)；

邊界層圖邊界層的尺度：其中：在20°C時(shí)，U=1m/s,L=1m，ν=10-6m2/s；水：Re≈106；空氣：Re≈6.7X104北京工業(yè)大學(xué)市政學(xué)科部——馬長(zhǎng)明高等流體(水）力學(xué)講稿第三講大雷諾數(shù)流體力學(xué)問(wèn)題的求解

（2）邊界層內(nèi)N—S方程的量級(jí)分析

速度：壓強(qiáng)：長(zhǎng)度—x方向：Y方向：其中：1）連續(xù)方程注意v的量級(jí)：北京工業(yè)大學(xué)市政學(xué)科部——馬長(zhǎng)明高等流體(水）力學(xué)講稿第三講大雷諾數(shù)流體力學(xué)問(wèn)題的求解

（2）邊界層內(nèi)N—S方程的量級(jí)分析速度：壓強(qiáng)：長(zhǎng)度—x方向：Y方向：其中：2）運(yùn)動(dòng)方程由此得：及即：北京工業(yè)大學(xué)市政學(xué)科部——馬長(zhǎng)明高等流體(水）力學(xué)講稿第三講大雷諾數(shù)流體力學(xué)問(wèn)題的求解

（2）邊界層內(nèi)N—S方程的量級(jí)分析速度：壓強(qiáng)：長(zhǎng)度—x方向：Y方向：2）運(yùn)動(dòng)方程，同理可得由此得：或：北京工業(yè)大學(xué)市政學(xué)科部——馬長(zhǎng)明高等流體(水）力學(xué)講稿第三講大雷諾數(shù)流體力學(xué)問(wèn)題的求解

（3）邊界層方程 按有量綱形式的方程表示

定解條件其中：在邊壁上：在邊界層界面上：北京工業(yè)大學(xué)市政學(xué)科部——馬長(zhǎng)明高等流體(水）力學(xué)講稿第三講大雷諾數(shù)流體力學(xué)問(wèn)題的求解 3、歐拉方程與邊界層方程的銜接條件在邊界層界面上有：可表示為：（銜接條件）

對(duì)平面繞流邊界靜止不動(dòng)的問(wèn)題，有v（x,0）=0，由x方向的歐拉方程，可確定邊界層壓強(qiáng)值，即：當(dāng)y→0：邊界層方程：北京工業(yè)大學(xué)市政學(xué)科部——馬長(zhǎng)明高等流體(水）力學(xué)講稿第三講大雷諾數(shù)流體力學(xué)問(wèn)題的求解 4、大雷諾數(shù)問(wèn)題解決步驟 （1）求解無(wú)邊界層的歐拉方程，得到邊界層外整個(gè)流場(chǎng)的速度分布、壓強(qiáng)分布； （2）確定物體所受的升力值，及 （3）壁面上的速度值和壓強(qiáng)值； （4）計(jì)算邊界層內(nèi)的速度分布； （5）確定邊界層厚度（可用于對(duì)歐拉計(jì)算的邊界修整）， （6）結(jié)算壁面剪應(yīng)力分布，計(jì)算壁面的阻力值。北京工業(yè)大學(xué)市政學(xué)科部——馬長(zhǎng)明高等流體(水）力學(xué)講稿第三講大雷諾數(shù)流體力學(xué)問(wèn)題的求解

二、歐拉方程的求解—?jiǎng)萘骼碚?/p>

1、基本方程與定界條件 重力場(chǎng)中的恒定的平面流動(dòng)：基本方程：在無(wú)窮遠(yuǎn)處：在物面上：當(dāng)物體靜止時(shí)：定解條件：北京工業(yè)大學(xué)市政學(xué)科部——馬長(zhǎng)明高等流體(水）力學(xué)講稿第三講大雷諾數(shù)流體力學(xué)問(wèn)題的求解 2、勢(shì)函數(shù)φ(x,y)

當(dāng)流動(dòng)無(wú)旋時(shí)，有： 上式稱(chēng)為L(zhǎng)aplace方程，為線性方程，解可疊加。 對(duì)無(wú)旋流動(dòng)的空間問(wèn)題，也存在勢(shì)函數(shù)φ(x,y,z)，滿足Laplace方程。對(duì)平面問(wèn)題：令滿足必滿足無(wú)旋條件（有勢(shì)必?zé)o旋）代入連續(xù)方程北京工業(yè)大學(xué)市政學(xué)科部——馬長(zhǎng)明高等流體(水）力學(xué)講稿第三講大雷諾數(shù)流體力學(xué)問(wèn)題的求解 3、拉格朗日與伯努利方程（積分）

（1）蘭姆——葛羅米柯（Lamb-Γpombiko）方程 （2）力勢(shì)函數(shù)與壓力勢(shì)函數(shù) 對(duì)重力場(chǎng)和不可壓流體：對(duì)歐拉方程：由矢量恒等式可得蘭姆——葛羅米柯型方程北京工業(yè)大學(xué)市政學(xué)科部——馬長(zhǎng)明高等流體(水）力學(xué)講稿第三講大雷諾數(shù)流體力學(xué)問(wèn)題的求解

（3）拉格朗日方程 將：

在恒定流條件下，有：

拉格朗日方程適用條件：運(yùn)動(dòng)無(wú)旋。

常數(shù)c適用于整個(gè)流場(chǎng)。

代入蘭姆方程：得拉格朗日方程：北京工業(yè)大學(xué)市政學(xué)科部——馬長(zhǎng)明高等流體(水）力學(xué)講稿第三講大雷諾數(shù)流體力學(xué)問(wèn)題的求解

（4）伯努利方程 在運(yùn)動(dòng)為恒定條件下，蘭姆方程為：

伯努利方程的適用條件：定常流動(dòng)。

由此可知，可將求解歐拉方程問(wèn)題轉(zhuǎn)化為：

1）由Laplace方程求出φ；

2）由勢(shì)函數(shù)求速度場(chǎng)u、v、w；

3）由拉格朗日或伯努利方程求壓強(qiáng)場(chǎng)。沿流線積分，得伯努利方程：北京工業(yè)大學(xué)市政學(xué)科部——馬長(zhǎng)明高等流體(水）力學(xué)講稿第三講大雷諾數(shù)流體力學(xué)問(wèn)題的求解 4、流函數(shù)ψ(x,y)由平面連續(xù)方程：引入ψ(x,y)滿足：當(dāng)流動(dòng)無(wú)旋時(shí)，ψ(x,y)滿足Laplce方程：注意：1）流函數(shù)僅存在于平面流動(dòng)或軸對(duì)稱(chēng)流動(dòng)問(wèn)題；

2）無(wú)論流動(dòng)是否無(wú)旋，流函數(shù)都存在。但只有滿 足無(wú)旋流動(dòng)條件，才能滿足Laplace方程。 北京工業(yè)大學(xué)市政學(xué)科部——馬長(zhǎng)明高等流體(水）力學(xué)講稿第三講大雷諾數(shù)流體力學(xué)問(wèn)題的求解 5、流函數(shù)與勢(shì)函數(shù)的一些特性 （1）流函數(shù)ψ(x,y)=const，表示流線； （2）兩流函數(shù)之差表示通過(guò)兩流線間的流量值； （3）流線與等勢(shì)線正交； （4）速度勢(shì)函數(shù)φ(x,y)不能在域內(nèi)有極大值與極小值； （5）流體質(zhì)點(diǎn)的速度的模在流域中不能達(dá)到極值； （6）壓強(qiáng)在域內(nèi)不能出現(xiàn)極小值； （7）在單連同域內(nèi)，ψ(x,y)、φ(x,y)是單值函數(shù)，再多連同域內(nèi)可以是多值函數(shù)。

其中：Q—通過(guò)多連通域封閉曲線的流量；

?！囟噙B通域封閉曲線的速度環(huán)量；

m—環(huán)繞多輛通域的次數(shù)數(shù)。北京工業(yè)大學(xué)市政學(xué)科部——馬長(zhǎng)明高等流體(水）力學(xué)講稿第三講大雷諾數(shù)流體力學(xué)問(wèn)題的求解 6、有勢(shì)流動(dòng)的主要性質(zhì) （1）凱爾文定理（Kelvin‘sTheorem） 理想正壓流體在有勢(shì)力場(chǎng)中運(yùn)動(dòng)時(shí)，連續(xù)流場(chǎng)中沿封閉流體線的速度環(huán)量不隨時(shí)間而變化。 （2）拉格朗日——渦量不生不滅定理 理想正壓流體在有勢(shì)力場(chǎng)中運(yùn)動(dòng)時(shí)，如某一時(shí)刻連續(xù)流場(chǎng)中無(wú)旋，則流場(chǎng)始終無(wú)旋。 （3）亥姆霍茨（Helmholtz）—渦管、渦線保持定理 理想正壓流體在有勢(shì)力場(chǎng)中運(yùn)動(dòng)時(shí)，滿足：

1）組成渦面的流體質(zhì)點(diǎn)永遠(yuǎn)組成渦面；

2）組成渦線的流體質(zhì)點(diǎn)永遠(yuǎn)組成渦線；

3）組成渦管的流體質(zhì)點(diǎn)永遠(yuǎn)組成渦管，且其強(qiáng)度不隨時(shí)間而變。北京工業(yè)大學(xué)市政學(xué)科部——馬長(zhǎng)明高等流體(水）力學(xué)講稿第三講大雷諾數(shù)流體力學(xué)問(wèn)題的求解 7、平面有勢(shì)流動(dòng)問(wèn)題的復(fù)勢(shì)方法—奇點(diǎn)疊加法 （1）定常平面有時(shí)流動(dòng)的Laplace方程的求解

1）求解勢(shì)函數(shù)

2）求解流函數(shù)在無(wú)窮遠(yuǎn)處：在物面上：定解條件：在無(wú)窮遠(yuǎn)處：在物面上：定解條件：北京工業(yè)大學(xué)市政學(xué)科部——馬長(zhǎng)明高等流體(水）力學(xué)講稿第三講大雷諾數(shù)流體力學(xué)問(wèn)題的求解

（2）復(fù)勢(shì)及其與有勢(shì)流動(dòng)各量間的關(guān)系

1）復(fù)勢(shì)的定義

為調(diào)和函數(shù)；且滿足：柯西—黎曼條件（Cauchy—Riemann）構(gòu)造復(fù)變函數(shù)是解析復(fù)變函數(shù)。其中：北京工業(yè)大學(xué)市政學(xué)科部——馬長(zhǎng)明高等流體(水）力學(xué)講稿第三講大雷諾數(shù)流體力學(xué)問(wèn)題的求解 2）用復(fù)勢(shì)表示有勢(shì)流動(dòng)各物理量 ①f(z)=const，表示等勢(shì)線和流線簇 ②復(fù)速度與復(fù)勢(shì)的關(guān)系 ③封閉曲線的的環(huán)量與通過(guò)的流量與復(fù)勢(shì)的關(guān)系北京工業(yè)大學(xué)市政學(xué)科部——馬長(zhǎng)明高等流體(水）力學(xué)講稿第三講大雷諾數(shù)流體力學(xué)問(wèn)題的求解 ④柱體所受流體作用力、力矩與復(fù)勢(shì)的關(guān)系作用力：作用力矩：北京工業(yè)大學(xué)市政學(xué)科部——馬長(zhǎng)明高等流體(水）力學(xué)講稿第三講大雷諾數(shù)流體力學(xué)問(wèn)題的求解

（3）復(fù)勢(shì)的確定—奇點(diǎn)疊加法

1）步驟

①研究基本復(fù)數(shù)函數(shù)所表示的流場(chǎng)； ②通過(guò)基本復(fù)勢(shì)的疊加，構(gòu)造滿足物面邊界條件的流場(chǎng)。物面是一條流線。 ③通過(guò)邊界條件（Q、Γ）確定待定常數(shù)。 ④計(jì)算流場(chǎng)速度、壓強(qiáng)分布。 ⑤計(jì)算物體所受的力。

2）基本復(fù)變函數(shù)所表示的流場(chǎng)

①均勻流f(z)=cz=c(x+iy)

其中：c是復(fù)常數(shù)。均勻流圖示北京工業(yè)大學(xué)市政學(xué)科部——馬長(zhǎng)明高等流體(水）力學(xué)講稿第三講大雷諾數(shù)流體力學(xué)問(wèn)題的求解

②源、匯流

a為實(shí)數(shù)，a>0：源；a<0：匯。

③點(diǎn)渦流

a為實(shí)數(shù)，a>0：環(huán)量逆時(shí)針轉(zhuǎn)動(dòng)；

a<0：環(huán)量順時(shí)針轉(zhuǎn)動(dòng)。源匯流圖示點(diǎn)渦流圖示柱坐標(biāo)表示的速度關(guān)系北京工業(yè)大學(xué)市政學(xué)科部——馬長(zhǎng)明高等流體(水）力學(xué)講稿第三講大雷諾數(shù)流體力學(xué)問(wèn)題的求解

④偶極子與無(wú)環(huán)量圓柱繞流 偶極子

⑤無(wú)環(huán)量圓柱繞流

偶極子圖示其中：偶極強(qiáng)度或偶極矩，方向以Q至–Q。圓柱繞流圖示其中圓柱半徑：北京工業(yè)大學(xué)市政學(xué)科部——馬長(zhǎng)明高等流體(水）力學(xué)講稿第三講大雷諾數(shù)流體力學(xué)問(wèn)題的求解

（4）復(fù)勢(shì)確定的鏡像法

1）平面鏡像

2）圓周鏡像

其中平面鏡像圖示圓周鏡像圖示a：圓周半徑，

表示在復(fù)數(shù)中除變量z外對(duì)所有復(fù)數(shù)量取共軛。北京工業(yè)大學(xué)市政學(xué)科部——馬長(zhǎng)明高等流體(水）力學(xué)講稿第三講大雷諾數(shù)流體力學(xué)問(wèn)題的求解

（5）復(fù)勢(shì)確定的保角變換法

1）保角變換——解析函數(shù)變換過(guò)程中兩線段的夾角不變，長(zhǎng)短可變。

2）變換思路

①找到一解析函數(shù)ζ=g(z)，將物理平面D中的曲線L變換到映射平面D’中的圓周L’； ②在映射平面D’中求出圓周繞流的復(fù)勢(shì)F(ζ)； ③通過(guò)z=g-1(ζ)的反函數(shù)將所求得的F(ζ)變換到物理平面中得到f(z)； ④已知復(fù)勢(shì)后其它各量可求。保角變換圖示北京工業(yè)大學(xué)市政學(xué)科部——馬長(zhǎng)明高等流體(水）力學(xué)講稿第三講大雷諾數(shù)流體力學(xué)問(wèn)題的求解

三、邊界層內(nèi)流體運(yùn)動(dòng)量的解析解法

1、邊界層（BoundaryThickness）的定義

（1）邊界層厚度（名義厚度）δ(x)

在x處的邊界層厚度定義為u(x,y)=0.99UE(x)所對(duì)應(yīng)的y值。

存在的問(wèn)題：δ(x)難以準(zhǔn)確測(cè)定。 （2）位移厚度（排擠厚度DisplacementThichness）δd(x)

由于邊界層的影響，對(duì)應(yīng)于x處通過(guò)邊界層δ(x)的流體質(zhì)量，無(wú)邊界層所需尺度為[δ(x)-δd(x)]。

即：得：北京工業(yè)大學(xué)市政學(xué)科部——馬長(zhǎng)明高等流體(水）力學(xué)講稿第三講大雷諾數(shù)流體力學(xué)問(wèn)題的求解

（3）動(dòng)量損失厚度（MomentThichness）δM(x)

由于邊界層的影響，對(duì)應(yīng)于x處，無(wú)邊界層通過(guò)尺度為[δ(x)-δd(x)]，其動(dòng)量通量與通過(guò)邊界層δ(x)的動(dòng)量通量之差成為動(dòng)量損失。定義動(dòng)量損失厚度：即動(dòng)量損失：有定義可知：位移與動(dòng)量厚度圖示北京工業(yè)大學(xué)市政學(xué)科部——馬長(zhǎng)明高等流體(水）力學(xué)講稿第三講大雷諾數(shù)流體力學(xué)問(wèn)題的求解 2、平面定常邊界層方程所給問(wèn)題是適定的。北京工業(yè)大學(xué)市政學(xué)科部——馬長(zhǎng)明高等流體(水）力學(xué)講稿第三講大雷諾數(shù)流體力學(xué)問(wèn)題的求解

3、邊界層方程的相似解

（1）相似解的概念 對(duì)邊界層速度u(x,y)存在適當(dāng)變換，使得每個(gè)橫斷面上的u(x,y)經(jīng)過(guò)變換后的分布都相同，即： 通過(guò)相似變換，將偏微分方程轉(zhuǎn)化為常微分方程。相似參數(shù)北京工業(yè)大學(xué)市政學(xué)科部——馬長(zhǎng)明高等流體(水）力學(xué)講稿第三講大雷諾數(shù)流體力學(xué)問(wèn)題的求解

（2）相似解存在條件

當(dāng)流動(dòng)無(wú)特征長(zhǎng)度（繞流邊界無(wú)限長(zhǎng)），邊界層外速度UE(x)=cxm，即為冪級(jí)數(shù)時(shí)就存在相似解。 其中α、β必須為常系數(shù)。

（3）相似解求解方法——引入流函數(shù)

引入流函數(shù)，使得對(duì)兩個(gè)變量（u、v）的方程組轉(zhuǎn)化為對(duì)一個(gè)相似參數(shù)η的常微分方程。北京工業(yè)大學(xué)市政學(xué)科部——馬長(zhǎng)明高等流體(水）力學(xué)講稿第三講大雷諾數(shù)流體力學(xué)問(wèn)題的求解

（4）邊界層相似方程

將所定義的流函數(shù)代入可得： 其中：UE’g’為對(duì)x的導(dǎo)數(shù)；f’表示對(duì)η的導(dǎo)數(shù)。 代入邊界層連續(xù)方程和運(yùn)動(dòng)方程，整理可得：邊界條件：其中：北京工業(yè)大學(xué)市政學(xué)科部——馬長(zhǎng)明高等流體(水）力學(xué)講稿第三講大雷諾數(shù)流體力學(xué)問(wèn)題的求解

（5）半無(wú)限平板恒定層流邊界層的Blasius解

1）流動(dòng)特點(diǎn)：

流動(dòng)圖示取則定解條件：方程簡(jiǎn)化為：

方程是非線性的，Blasius給出級(jí)數(shù)解，Howarth給出精度較高的數(shù)值解，見(jiàn)教材p.82表4-1所示。北京工業(yè)大學(xué)市政學(xué)科部——馬長(zhǎng)明高等流體(水）力學(xué)講稿第三講大雷諾數(shù)流體力學(xué)問(wèn)題的求解 2）數(shù)值成果分析 速度分布： 邊界層厚度： 位移厚度 動(dòng)量損失厚度

數(shù)值解成果速度分布圖北京工業(yè)大學(xué)市政學(xué)科部——馬長(zhǎng)明高等流體(水）力學(xué)講稿第三講大雷諾數(shù)流體力學(xué)問(wèn)題的求解

板面剪應(yīng)力 剪應(yīng)力系數(shù) 板長(zhǎng)為L(zhǎng)的阻力 阻力系數(shù)數(shù)值解成果北京工業(yè)大學(xué)市政學(xué)科部——馬長(zhǎng)明高等流體(水）力學(xué)講稿第三講大雷諾數(shù)流體力學(xué)問(wèn)題的求解 3）計(jì)算結(jié)果北京工業(yè)大學(xué)市政學(xué)科部——馬長(zhǎng)明高等流體(水）力學(xué)講稿第三講大雷諾數(shù)流體力學(xué)問(wèn)題的求解 4、定常平面邊界層的卡門(mén)（KarmanT.Von）動(dòng)量積分

方法特點(diǎn)： 是一種近似方法； 僅在積分意義上滿足邊界層方程； 簡(jiǎn)單，有一定精度，在工程中得到廣泛應(yīng)用。 （1）基本方程與邊界條件北京工業(yè)大學(xué)市政學(xué)科部——馬長(zhǎng)明高等流體(水）力學(xué)講稿第三講大雷諾數(shù)流體力學(xué)問(wèn)題的求解

（2）動(dòng)量積分關(guān)系 將（1-1）式乘u后與（1-2）式相加，得： 將（1-1）式乘UE后再加udUE/dx，得： 將（2-2）-（2-1）后沿邊界層積分，得：引入邊界層定義，得北京工業(yè)大學(xué)市政學(xué)科部——馬長(zhǎng)明高等流體(水）力學(xué)講稿第三講大雷諾數(shù)流體力學(xué)問(wèn)題的求解

引入定義 形狀因子 壁面摩察系數(shù) 方程特點(diǎn)：僅有一個(gè)方程，但有三個(gè)未知量δM，Hdm，cf； 當(dāng)已知邊界層速度分布后，則δM，Hdm，cf可以確定。 得卡門(mén)動(dòng)量積分方程：北京工業(yè)大學(xué)市政學(xué)科部——馬長(zhǎng)明高等流體(水）力學(xué)講稿第三講大雷諾數(shù)流體力學(xué)問(wèn)題的求解

（3）邊界層近似求解步驟

1）假設(shè)邊界層速度分布形式

2）由邊界條件確定待定系數(shù)，得到邊界層相容性方程

3）計(jì)算邊界層名義厚度

4）計(jì)算邊界層其余參量北京工業(yè)大學(xué)市政學(xué)科部——馬長(zhǎng)明高等流體(水）力學(xué)講稿第三講大雷諾數(shù)流體力學(xué)問(wèn)題的求解

（4）應(yīng)用例 對(duì)無(wú)窮遠(yuǎn)水平來(lái)流U∞=Const繞半無(wú)限平板邊界層內(nèi)流動(dòng)， 設(shè)： 確定相容性方程，并計(jì)算邊界層個(gè)參量。 解：1）由邊界（相容）條件確定a、b、c、d四個(gè)系數(shù)

北京工業(yè)大學(xué)市政學(xué)科部——馬長(zhǎng)明高等流體(水）力學(xué)講稿第三講大雷諾數(shù)流體力學(xué)問(wèn)題的求解

解得： 相容方程為

2）由動(dòng)量積分方程確定邊界層厚度 位移厚度： 動(dòng)量損失厚度： 形狀因子： 壁面剪應(yīng)力： 北京工業(yè)大學(xué)市政學(xué)科部——