第6章 邊界層流動(dòng)

1、 6.1 邊界層基本概念邊界層基本概念 6.1.1 邊界層流態(tài) 6.1.2 邊界層各特征厚度 6.2 二維平面邊界層流動(dòng) 6.3 二維曲面邊界層流動(dòng) 6.4 二維圓柱滑動(dòng)軸承潤(rùn)滑 6.5 圓柱和圓球繞流阻力 實(shí)際流體繞任何形狀物體的大雷諾數(shù)流動(dòng)都會(huì)在物面附近形成邊界層。圖6-1所示為空氣繞某一翼型的流動(dòng)，整個(gè)流場(chǎng)可分為邊界層、邊界層脫離翼型物面以后形成的尾流、以及邊界層和尾流以外的勢(shì)流。圖圖6-1 翼型繞流翼型繞流6.1.1 邊界層流態(tài)邊界層流態(tài) 邊界層流動(dòng)可以是層流或湍流。實(shí)際中更一般地是混合邊界層，即邊界層前緣為層流，經(jīng)過一過渡區(qū)(稱為轉(zhuǎn)捩區(qū))后轉(zhuǎn)變?yōu)橥牧?；在湍流區(qū)，緊挨物面附近還有一層流

2、底層。圖6-2所示為一均勻來流繞過平板一側(cè)所形成的邊界層流動(dòng)。圖圖6-2 平板邊界層流動(dòng)平板邊界層流動(dòng) 在湍流區(qū)，若平板表面粗糙度D大于層流底層的厚度dl，則稱之為粗糙(表面)平板；否則稱為光滑(表面)平板。當(dāng)層流區(qū)的范圍很小時(shí)，可近似地把整個(gè)邊界層看成為湍流邊界層。 為了便于判斷邊界層的流態(tài)，通常假定由層流到湍流的轉(zhuǎn)捩是在某一截面突變完成的，并稱此截面為臨界截面，它離邊界層前緣的距離稱為臨界長(zhǎng)度x*，臨界截面邊界層的厚度稱為臨界厚度d*。（圖6-2） 邊界層流態(tài)用臨界雷諾數(shù)Re*來判斷， Re*有兩種形式：Rex* = Ux*/u和 Red* = Ud*/u，對(duì)于平板繞流，Rex* = 51

3、05 3106，Red* 2800。6.1.2 邊界層各特征厚度邊界層各特征厚度 邊界層厚度邊界層厚度 邊界層理論將大雷諾數(shù)流動(dòng)的流場(chǎng)分為粘性區(qū)和無粘區(qū)兩部分，分別稱為邊界層和主流區(qū)，它們的交界面稱為邊界層(外)邊界，并人為地規(guī)定邊界層邊界上流速為主流區(qū)的99(或99.5%), 邊界層邊界到物面的距離稱為邊界層厚度d，用數(shù)學(xué)式表示即有 邊界層未脫離物面的情況下，邊界層厚度沿流程是增加的，即在迎流的前緣點(diǎn)為零，然后沿流動(dòng)方向逐漸增加，到送流的后緣點(diǎn)達(dá)到最大。 邊界層位移厚度邊界層位移厚度也稱邊界層排擠厚度。在邊界層內(nèi)，流速受到壁面的阻滯作用而減小，使通過邊界層內(nèi)的流量比理想流動(dòng)時(shí)減少，這相當(dāng)于固

4、體壁面沿其法線方向朝流場(chǎng)內(nèi)移動(dòng)了一個(gè)距離d1后理想流動(dòng)所通過的流量，這個(gè)d1就是邊界層位移厚度，如圖6-3所示。根據(jù)位移厚度d1的定義，對(duì)不可壓流動(dòng)有圖圖6-3 邊界層位移厚度邊界層位移厚度即 邊界層動(dòng)量厚度邊界層動(dòng)量厚度 與理想流動(dòng)相比，邊界層內(nèi)流速降低一方面使通過的流體質(zhì)量減少，另一方面也使通過的流體動(dòng)量減少。這種動(dòng)量減小也可以看成是相當(dāng)于將固體壁面向流場(chǎng)內(nèi)移動(dòng)了一個(gè)距離d2： 稱d2為動(dòng)量損失厚度，簡(jiǎn)稱動(dòng)量厚度。邊界層的位移厚度與動(dòng)量厚度之比稱為邊界層形狀因子： H = d1/d2。即 邊界層能量厚度邊界層能量厚度即邊界層能量損失厚度。與理想流體的流動(dòng)相比，邊界層內(nèi)流速的降低還使流體的動(dòng)

5、能通量減少。類似于動(dòng)量厚度，可以定義不可壓流動(dòng)的邊界層能量厚度d3：以上定義式表示邊界層實(shí)際的流量具有的理想流動(dòng)動(dòng)能與實(shí)際流動(dòng)動(dòng)能之差。容易證明，在邊界層任一截面，恒有：d d1 d3 d2。即 6.1 邊界層基本概念 6.2 二維平面邊界層流動(dòng)二維平面邊界層流動(dòng) 6.2.1 微分方程及其精確解 6.2.2 積分方程及其近似解 6.3 二維曲面邊界層流動(dòng) 6.4 二維圓柱滑動(dòng)軸承潤(rùn)滑 6.5 圓柱和圓球繞流阻力 二維平面不可壓邊界層流動(dòng)是最簡(jiǎn)單的一類粘性流動(dòng)，即便如此也只有極少數(shù)情況能通過邊界層微分方程求得精確解，大多數(shù)情況只能通過邊界層積分方程求近似解。6.2.1 微分方程及其精確解微分方程

6、及其精確解 微分方程微分方程 在直角坐標(biāo)系，定常、不可壓、不計(jì)重力的二維流動(dòng)N-S方程為 根據(jù)小粘度二維平面邊界層流動(dòng)的特點(diǎn)d L以及uy ux對(duì) N-S方程中各變量和參數(shù)作數(shù)量級(jí)估計(jì)，有 量級(jí)1的量： dx, dx2; ux, dux, d2ux; p, dp; r 量級(jí)e 1的量：dy; uy, duy, d2uy 量級(jí)e 2的量： dy2; u依照以上量級(jí)對(duì)N-S方程進(jìn)行簡(jiǎn)化分析，可得微微分分方方程程及及其其精精確確解解 以上就是二維平面邊界層流動(dòng)的微分方程，由普朗特在1904年首次提出。雖然普朗特邊界層微分方程相對(duì)N-S方程大為簡(jiǎn)化，但仍然是非線性的，只能對(duì)特殊情況下的某些層流邊界層求

7、得精確解。 求解邊界層微分方程時(shí)，首先要得到邊界層外部勢(shì)流的速度，使壓強(qiáng)p成為已知量，這樣未知量只有ux和uy，由邊界層微分方程x分式和連續(xù)方程一起構(gòu)成封閉的求解系。注意：普朗特邊界層微分方程不適用于d /x 1條件得不到滿足的邊界層前緣部分，該部分對(duì)應(yīng)的雷諾數(shù)范圍一般為Rex 25。微微分分方方程程及及其其精精確確解解 微分方程的精確解微分方程的精確解 應(yīng)用邊界層微分方程解決粘性流動(dòng)問題的一個(gè)最簡(jiǎn)單的例子，是流體繞順流放置平板的層流邊界層流動(dòng)，即均勻來流繞過沿平行于流動(dòng)方向放置的一塊薄平板(其厚度假設(shè)為零)并在平板一側(cè)附近所產(chǎn)生的流動(dòng)。微微分分方方程程及及其其精精確確解解圖圖6-4 平板層流

8、邊界層平板層流邊界層 微分方程的精確解微分方程的精確解 如圖6-4所示，取平板前緣為直角坐標(biāo)系的原點(diǎn)，則平板前方未受擾動(dòng)的均勻來流速度U與平板平行。由伯努利方程知，在繞平板流動(dòng)的勢(shì)流部分，U = U、dp/dx = 0；而由邊界層微分方程知，在邊界層中壓強(qiáng)沿y方向是均勻分布的，即邊界層內(nèi)任一點(diǎn)處的壓強(qiáng)都與同x坐標(biāo)處邊界層外勢(shì)流的壓強(qiáng)相等。微微分分方方程程及及其其精精確確解解圖圖6-4 平板層流邊界層平板層流邊界層 在邊界層微分方程和連續(xù)方程中引入流函數(shù)y，則由流函數(shù)定義有: y/x = -uy，y/y = ux，連續(xù)方程ux/x + uy/y = 0自動(dòng)滿足，邊界層微分方程成為微微分分方方程程

9、及及其其精精確確解解圖圖6-4 平板層流邊界層平板層流邊界層 因?yàn)閐 f2 f3 )上式第一個(gè)(即一階)漸近解就是勢(shì)流解，即f1 = h + bb 為積分常數(shù)。令f(h)的二階漸近解為f = f1 + f2并代入原常微分方程2f+ff =0積分，得微微分分方方程程及及其其精精確確解解g 為另一積分常數(shù)。 類似還可得三階漸近解f = f1 + f2 + f3甚至更高階漸近解，本問題中僅考慮到二階。 級(jí)數(shù)解由邊界層靠近壁面向外求解，漸進(jìn)解則由邊界層外的勢(shì)流向內(nèi)求解，兩種解在邊界層內(nèi)某一點(diǎn)必須匹配，即兩種解在這一點(diǎn)的f、f 、f 值都相等，由此得到A2 = 0.332, b = 1.72, g =

10、0.231整個(gè)流動(dòng)問題得解。繼布拉休斯之后，其他學(xué)者也對(duì)二維平面層流邊界層流動(dòng)即方程2f + ff = 0進(jìn)行了數(shù)值求解，其中霍華斯在1938年得到的結(jié)果對(duì)照實(shí)驗(yàn)具有更好的準(zhǔn)確度。 根據(jù)霍華斯的結(jié)果，在h = 5.0處ux /U = u /U = f = 0.99155，將它作為邊界層邊界，通過積分可得平板層流邊界層各特征量如下微微分分方方程程及及其其精精確確解解微微分分方方程程及及其其精精確確解解 邊界層厚度： 邊界層位移厚度： 邊界層動(dòng)量厚度： 壁面切應(yīng)力系數(shù)： 摩擦阻力系數(shù)：t0為壁面切應(yīng)力、FDf為整個(gè)平板受到的力，即 以上結(jié)果得到試驗(yàn)的證實(shí)。圖6-5表示順流放置平板層流邊界層的布拉休

11、斯精確解，以及據(jù)此繪制的邊界層厚度的沿程變化和流速分布。微微分分方方程程及及其其精精確確解解 圖圖6-5 順流放置平板層流邊界層流動(dòng)順流放置平板層流邊界層流動(dòng) 對(duì)于非順流放置平板的繞流流動(dòng)，理論指出，只要?jiǎng)萘髁魉賃與x坐標(biāo)(沿平板表面)成冪指數(shù)關(guān)系：U = C x m (C為常數(shù)、m為有理數(shù))，邊界層微分方程微微分分方方程程及及其其精精確確解解就存在相似性解，這時(shí)流速u(x, y)的分布具有這樣的性質(zhì)：如果把任意斷面x上的流速分布圖形u-y的u和y坐標(biāo)分別用有關(guān)尺度因子變換為量綱一的坐標(biāo)u0和y0 ，則在任何x斷面上u0-y0的分布圖形都相同。 順流放置平板繞流的精確解只是邊界層微分方程相似性

12、解中的一個(gè)特例，對(duì)應(yīng)于m = 0。6.2.2 積分方程及其近似解積分方程及其近似解 積分方程積分方程 對(duì)定常不可壓二維平面邊界層流動(dòng)，取控制體1221進(jìn)行分析，如圖6-8所示。在截面1-1和2-2上，流體參數(shù)分別為 圖圖6-8 平板邊界層流動(dòng)平板邊界層流動(dòng) 控制體流體在x方向受到的總作用力為整理并忽略高階小量后，簡(jiǎn)化為 通過控制面進(jìn)入和離開控制體的流體在x方向的動(dòng)量分別為1-1 截面：2-2 截面：1-2 截面：積積分分方方程程及及其其近近似似解解 將以上4個(gè)式子代入動(dòng)量方程x分式，就得在上式中代入以下邊界條件并整理得由于ux u，上式兩個(gè)積分項(xiàng)分別為位移厚度和動(dòng)量厚度，所以邊界層動(dòng)量積分方程

13、為積積分分方方程程及及其其近近似似解解 上式由卡門在1921年根據(jù)動(dòng)量定理首次導(dǎo)出，故又稱為卡門動(dòng)量積分方程，其邊界條件為 邊界層動(dòng)量積分方程對(duì)層流和湍流都適用，對(duì)于順流放置平板的邊界層流動(dòng)則簡(jiǎn)化為 邊界層動(dòng)量積分方程還可由邊界層微分方程在邊界層內(nèi)對(duì)y進(jìn)行積分獲得。此外，用流速u乘以邊界層微分方程中的每一項(xiàng)并對(duì)y進(jìn)行積分，還可得到邊界層能量積分方程。積積分分方方程程及及其其近近似似解解 積分方程的近似解積分方程的近似解 邊界層動(dòng)量積分方程中包含壁面切應(yīng)力t0，邊界層位移厚度d1和動(dòng)量厚度d2三個(gè)未知量；由d1和d2的定義式以及壁面邊界條件還可以補(bǔ)充三個(gè)方程，但又出現(xiàn)另外兩個(gè)未知量(流速u和邊界

14、層厚度d)，因此邊界層動(dòng)量積分方程在數(shù)學(xué)上是不封閉的，只宜采用近似方法求解。通常的做法是，首先假定某種速度分布，據(jù)此算得d1(d)、d2 (d)和t0(d)，然后將它們代入邊界層動(dòng)量積分方程，最后通過積分求得邊界層厚度d及阻力系數(shù)CDf等特征量。 上述做法的特點(diǎn):只在物面及邊界層外緣滿足邊界層微分方程；假設(shè)的邊界層內(nèi)流速分布與實(shí)際不一定吻合。積積分分方方程程及及其其近近似似解解 1) 二維平面層流邊界層的近似解二維平面層流邊界層的近似解 設(shè)邊界層內(nèi)流速分布u/U = sin(py/2d)，則有積積分分方方程程及及其其近近似似解解將以上的d2和t0代入邊界層動(dòng)量積分方程，得由上式解得邊界層厚度d

15、，并計(jì)算其他特征厚度和系數(shù)得微微分分方方程程及及其其精精確確解解 邊界層厚度： 邊界層位移厚度： 邊界層動(dòng)量厚度： 壁面切應(yīng)力系數(shù)： 摩擦阻力系數(shù)： 假設(shè)不同的邊界層流速分布，得到的邊界層 各特征厚度和系數(shù)也不相同。 2) 二維平面湍流邊界層的近似解二維平面湍流邊界層的近似解 應(yīng)用水力光滑圓管湍流的實(shí)驗(yàn)成果，可設(shè)光滑壁面平板湍流邊界層流速分布、壁面切應(yīng)力分別為積積分分方方程程及及其其近近似似解解就有積分并利用平板前緣點(diǎn)條件x = 0: d = 0，得光滑平板湍流邊界層各特征厚度和系數(shù)如下將d2和t0代入邊界層動(dòng)量積分方程，得 邊界層厚度： 邊界層位移厚度： 邊界層動(dòng)量厚度： 壁面切應(yīng)力系數(shù)：

16、摩擦阻力系數(shù)：根據(jù)試驗(yàn)數(shù)據(jù)，上面的摩擦阻力系數(shù)應(yīng)修正為 積積分分方方程程及及其其近近似似解解 在實(shí)際中，靠近平板前緣總有一部分是層流邊界層，因此摩擦阻力系數(shù)計(jì)算式須作進(jìn)一步修正。如圖6-9所示，假定層流向湍流的轉(zhuǎn)捩在某一斷面突然發(fā)生并完成，這樣整個(gè)平板的阻力就只需將轉(zhuǎn)捩斷面之前的那部分湍流阻力代之以層流阻力、其余部分湍流阻力則保持不變。積積分分方方程程及及其其近近似似解解 圖圖6-9 平板混合邊界層平板混合邊界層 轉(zhuǎn)捩斷面前湍流阻力與層流阻力之差為積積分分方方程程及及其其近近似似解解相應(yīng)摩擦阻力系數(shù)之差為所以光滑平板湍流邊界層的實(shí)際阻力系數(shù)為式中的A值與臨界雷諾數(shù)Re*的對(duì)應(yīng)關(guān)系如下表 6.1

17、 邊界層基本概念 6.2 二維平面邊界層流動(dòng) 6.3 二維曲面邊界層流動(dòng)二維曲面邊界層流動(dòng) 6.3.1 邊界層方程 6.3.2 邊界層分離 6.3.3 層流邊界層的卡門-波爾毫森解法 6.3.4 湍流邊界層的海特近似解法 6.4 二維圓柱滑動(dòng)軸承潤(rùn)滑 6.5 圓柱和圓球繞流阻力 流體繞曲面物體流動(dòng)時(shí)，邊界層外勢(shì)流流速將隨曲面曲率的變化而改變，壓強(qiáng)也隨之變化。 分析彎曲壁面附近的邊界層流動(dòng)通常采用隨體坐標(biāo)系，或邊界層坐標(biāo)系。這是一種特殊的正交曲線坐標(biāo)系，它以壁面前駐點(diǎn)O為原點(diǎn)、以沿壁面指向下游為x坐標(biāo)、自壁面算起沿壁面外法線為y坐標(biāo)，圖6-11和圖6-12分別表示二維曲面和軸對(duì)稱曲面的隨體坐標(biāo)系

18、。 邊界層內(nèi)任一點(diǎn)的坐標(biāo)為x = OP0，y = P0P；若Q為P的鄰點(diǎn)并且PQ = ds(圖6-11)，則ds在過P點(diǎn)的x和y坐標(biāo)上的投影分別為式中，h1、h2分別為坐標(biāo)x, y的拉梅系數(shù)。以R(x)表示曲面在P0點(diǎn)的曲率半徑，df表示點(diǎn)P0和Q0處曲率半徑間的夾角，則有 若為軸對(duì)稱曲面邊界層，則R(x)是子午面內(nèi)壁輪廓線的曲率半徑，r(x)為邊界層內(nèi)任意點(diǎn)到對(duì)稱軸的距離、即回轉(zhuǎn)半徑，r0(x)為軸對(duì)稱曲面上任意點(diǎn)的回轉(zhuǎn)半徑。圖圖6-11 二維曲面隨體坐標(biāo)二維曲面隨體坐標(biāo)圖圖6-12 軸對(duì)稱曲面隨體坐標(biāo)軸對(duì)稱曲面隨體坐標(biāo) 6.3.1 邊界層方程邊界層方程 微分方程微分方程 對(duì)于定常不可壓二維

19、曲面和軸對(duì)稱曲面邊界層流動(dòng)，采用類似于二維平面邊界層流動(dòng)的量級(jí)分析方法，可得相應(yīng)的微分方程，即式中k = 0(二維曲面)或 k = 1(軸對(duì)稱曲面)。 二維曲面和軸對(duì)稱曲面邊界層微分方程的邊界條件為y = 0：u = 0；y = ：u = U邊界層外部勢(shì)流的速度U由理想流體繞同一物面流動(dòng)的歐拉方程解確定，然后利用伯努利方程得到壓強(qiáng)梯度。 和二維平面邊界層微分方程類似，二維曲面和軸對(duì)稱曲面邊界層微分方程必須滿足限制條件uy / ux 1，/L 1。除此外，還要求R (x) L、r0 (x) L，即曲面的曲率半徑、回轉(zhuǎn)半徑與流動(dòng)方向的曲面總長(zhǎng)度為相同量級(jí)。 積分方程積分方程 對(duì)于二維任意形狀物體的

20、繞流，采用邊界層坐標(biāo)后，只要uy/ux 1的條件得以滿足，就仍可使用卡門動(dòng)量積分方程： 二維曲面或軸對(duì)稱曲面邊界層外勢(shì)流的流速和壓強(qiáng)不再是常數(shù)，導(dǎo)致在逆壓梯度(壓強(qiáng)沿流動(dòng)方向增大)的地方有可能發(fā)生邊界層分離，這時(shí)邊界層內(nèi)的流體在曲面的某些部位脫離曲面，使這部分曲面不再起“導(dǎo)流”作用，引起受粘性影響的流場(chǎng)范圍和流動(dòng)阻力迅速增大。6.3.2 邊界層分離邊界層分離 分離現(xiàn)象分離現(xiàn)象 實(shí)際流體在繞曲面流動(dòng)途中，邊界層內(nèi)流體有可能在外部勢(shì)流區(qū)逆壓梯度的作用下從曲面某個(gè)部位開始脫離曲面，使部分曲面不再起導(dǎo)流作用，這種現(xiàn)象稱為邊界層分離。邊界層從壁面分離后，如果外部勢(shì)流區(qū)的壓強(qiáng)梯度改善為順壓梯度或零梯度，則

21、邊界層可重回壁面附近，稱之為邊界層重新附著；例如流體繞順流放置平板上的一個(gè)階梯流動(dòng)時(shí)，邊界層在階梯后緣發(fā)生分離，在經(jīng)過階梯后一段距離將重新回附在平板附近。邊界層分離將導(dǎo)致受粘性影響的流場(chǎng)范圍和流動(dòng)阻力迅速增大，在實(shí)際中通常需要避免邊界層發(fā)生分離。 圖6-13表示實(shí)際的二維圓柱繞流流場(chǎng)。如果整個(gè)流場(chǎng)均為無粘勢(shì)流，則流體從圓柱的前緣D至頂點(diǎn)E是加速的、從頂點(diǎn)E至后緣F是減速的。由伯努利方程知，在DE流段壓強(qiáng)沿流動(dòng)方向逐漸減小(dp/dx0)，稱為逆壓梯度；在圓柱后緣點(diǎn)F壓強(qiáng)恢復(fù)到前緣點(diǎn)D的數(shù)值，即恢復(fù)到駐點(diǎn)壓強(qiáng)：rU2/2 (表壓)。 實(shí)際流動(dòng)中，圓柱附近為邊界層，其中流體因粘性作用而損耗能量，導(dǎo)

22、致在DE流段壓能的降低一部分轉(zhuǎn)化為動(dòng)能，其余則克服粘性阻力而消耗掉；在EF流段，流體動(dòng)能的降低一部分轉(zhuǎn)化為壓能，其余用于克服粘性阻力。圖圖6-13 二維圓柱繞流流場(chǎng)示意圖二維圓柱繞流流場(chǎng)示意圖 在圓柱后緣點(diǎn)F，壓強(qiáng)不能恢復(fù)到前緣點(diǎn)D的數(shù)值，而是在EF流段的某點(diǎn)S處，物面附近的流體動(dòng)能被消耗怡盡、流速降為零； 在S點(diǎn)的下游，外部勢(shì)流的壓強(qiáng)較高，導(dǎo)致流體在逆壓梯度的作用下發(fā)生回流，將邊界層內(nèi)的來流擠向主流而使邊界層脫離壁面、造成分離。S點(diǎn)稱為邊界層分離點(diǎn)，在分離點(diǎn)下游形成受粘性影響的回流和尾流區(qū)，其間滿布了大大小小的旋渦，造成較大的能量損失。尾流中壓強(qiáng)比無粘流動(dòng)時(shí)低，因此鈍形物體繞流形成的壓差阻力

23、遠(yuǎn)大于細(xì)長(zhǎng)的流線形物體。 飛機(jī)機(jī)翼是典型的流線型物體，其尾部逆壓梯度很小，使得分離點(diǎn)很靠近尾部而減小阻力。 邊界層分離前后流速分布邊界層分離前后流速分布 在分離點(diǎn)上游，dp/dx0，2u/y20、壁面附近產(chǎn)生流速為負(fù)值的回流區(qū)，且在y=0處有u/y0；在分離點(diǎn)S，壁面上(y=0)有u/y=0，實(shí)際中常根據(jù)這一條件(即速度的法向?qū)?shù)在物面為零)來確定分離點(diǎn)位置。分離點(diǎn)處流線與物面形成的角度a與雷諾數(shù)有關(guān)。 分離點(diǎn)下游的粘性流場(chǎng)范圍迅速增大，破壞了邊界層方程的限制條件：d 12時(shí)，出現(xiàn)u/U 1，顯然這是不允許的。因此，L的界限為-12 L12。 在速度分布多項(xiàng)式中引入形狀因子L并未增加新的變量，

24、L中的未知量d最終由邊界層積分方程的求解而確定。 將速度分布多項(xiàng)式代入邊界層位移厚度、動(dòng)量厚度、壁面切應(yīng)力定義式后進(jìn)行積分或運(yùn)算，同時(shí)對(duì)速度分布形狀因子L進(jìn)行微分，最后應(yīng)用邊界層動(dòng)量積分方程，就得到二維曲面繞流層流邊界層的求解方程組，即 卡卡門門 波波爾爾豪豪森森解解法法卡卡門門 波波爾爾豪豪森森解解法法 初始截面條件為x = 0：d = 0, L= 0(尖前緣)；或x = 0：U = 0, L = 7.052(鈍前緣)。 求解時(shí)，首先解得二維曲面物體的勢(shì)流速度分布U(x)及其一階導(dǎo)數(shù)U和二階導(dǎo)數(shù)U，然后對(duì)速度分布因子L的微分方程進(jìn)行數(shù)值求解得到L(x)，再由L的定義式解得d、由L代入速度分布

25、多項(xiàng)式解得u/U，最后由簡(jiǎn)化后的邊界層動(dòng)量積分方程、邊界層的位移厚度式和動(dòng)量厚度式解得d1、d2和t0，并由限制條件t0 = 0 (或L = -12)得到邊界層分離點(diǎn)位置xs，至此，整個(gè)層流邊界層問題得解?？ㄩT門 波波爾爾豪豪森森解解法法6.3.4 湍流邊界層的海特近似解法湍流邊界層的海特近似解法 流體繞二維曲面湍流邊界層的近似解法屬半經(jīng)驗(yàn)方法，涉及的經(jīng)驗(yàn)公式多且在不斷完善。 海特卷吸法以邊界層動(dòng)量積分方程和卷吸積分式為基本方程。因這2個(gè)基本方程是不封閉的，故采用壁面切應(yīng)力經(jīng)驗(yàn)公式Cf = Cf (H, Red2)作為補(bǔ)充方程，同時(shí)還進(jìn)一步假定一個(gè)卷吸速度分布 ，其中的 也由經(jīng)驗(yàn)公式給出。

26、以上基本方程和經(jīng)驗(yàn)公式一起構(gòu)成數(shù)學(xué)上封閉的方程組，在給定初值后就能進(jìn)行數(shù)值求解。 卷吸積分關(guān)系式卷吸積分關(guān)系式 卷吸速度UE是單位時(shí)間內(nèi)通過邊界層邊界的單位面積從外部勢(shì)流進(jìn)入邊界層的體積流量，如圖6-15所示，圖中虛線部分為控制體。由質(zhì)量守恒易得 上式稱為卷吸(流量)積分關(guān)系式。海特應(yīng)用的卷吸速度分布以及湍流邊界層形狀參數(shù)的經(jīng)驗(yàn)公式分別為 圖圖6-15 湍流邊界層的卷吸速度湍流邊界層的卷吸速度 當(dāng)?shù)刈枇ο禂?shù)經(jīng)驗(yàn)公式當(dāng)?shù)刈枇ο禂?shù)經(jīng)驗(yàn)公式 海特建議采用的當(dāng)?shù)刈枇ο禂?shù)經(jīng)驗(yàn)公式為 求解方法求解方法 將邊界層動(dòng)量積分方程和卷吸關(guān)系式改寫成然后進(jìn)行數(shù)值求解。具體步驟為，先求解邊界層外部勢(shì)流得到U(x)，再由給出的初始截面(通常為轉(zhuǎn)捩截面)上的d2和H值開始進(jìn)行數(shù)值積分(例如采用常用的龍格-庫塔積分方法進(jìn)行)，求解 ，直至Cf = 0、得到分離點(diǎn)位置為止。實(shí)際計(jì)算表明，分離點(diǎn)出現(xiàn)的位置對(duì)應(yīng)的H值位于1.8和2.8之間。 6.1 邊界層基本概念 6.2 二維平面邊界層流動(dòng) 6.3 二維曲面邊界層流動(dòng) 6.4 二維圓柱滑動(dòng)軸承潤(rùn)滑二維圓柱滑動(dòng)軸承潤(rùn)滑 6.4.1 雷諾潤(rùn)滑方程 6.4.2 二維圓柱滑動(dòng)軸承潤(rùn)滑 6.5 圓柱和圓球繞流阻力 直接接觸的固體壁面作相對(duì)運(yùn)動(dòng)時(shí)，因壁面粗糙部分的碰撞及不同壁面材料間的粘合而產(chǎn)生較大摩擦阻力，若不采取潤(rùn)滑措施減小摩擦將