第4章-管內(nèi)氣液兩相流阻力計(jì)算

1、Chapter 4. 管內(nèi)氣液兩相流的阻力計(jì)算(Gas-liquid flow resistance calculation)西安交通大學(xué)能源與動力工程學(xué)院王樹眾 教授 2引言摩擦壓降計(jì)算加速壓降局部阻力4.14.24.34.5重位壓力降4.4 3壓力降計(jì)算是氣液兩相流研究中最重要的課題之一只有正確地進(jìn)行壓力降計(jì)算，才能使系統(tǒng)具有安全可靠和足夠的壓頭，才能為動力設(shè)備的選型以及安全經(jīng)濟(jì)運(yùn)行提供必要的依據(jù)。氣液兩相流的壓力降包括四部分，即重位壓力降、摩擦壓力降、加速壓力降和局部阻力壓力降，亦即： PT= Pg+ Pf+ Pa+ Pb式中 PT 總壓力降 Pg 重位壓力降-重力作用而引起的 Pf 摩

2、擦壓力降-摩擦阻力引起的 Pa 加速壓力降-流體速度變化而引起的 Pb局部阻力壓力降-流動方向或管截面發(fā)生變化引起. l1 -摩擦壓降計(jì)算的經(jīng)驗(yàn)方法l2 -依據(jù)流型的摩擦壓降的計(jì)算l3 -影響氣液兩相流體摩擦阻力壓力降的主要因素 l1 -摩擦壓降計(jì)算的經(jīng)驗(yàn)方法均相流模型 對單相： （4-1）將兩相看作均勻混合的介質(zhì)： （4-2）又混合相密度： （x為平均干度）（4-3） 則上式成為： （4-4）均相模型（均相模型通常計(jì)算值偏低，有時差別還相當(dāng)大）2FLVP =D22m mFVLP =D2mGL1=x1-x+2m mLFLGVLP =1x(1)D2（）蘇聯(lián)50年代鍋爐水循環(huán)計(jì)算法該式計(jì)算誤差對水

3、平均相為（20%） l1 -摩擦壓降計(jì)算的經(jīng)驗(yàn)方法均相流模型 鑒于公式4-4誤差較大1978年蘇聯(lián)修訂后的鍋爐水力計(jì)算標(biāo)準(zhǔn)狀況為： （4-5）x-為平均干度，考慮了校正系數(shù) ，其與x，p，m 有關(guān)， 可查圖求得。對于受熱管上式中x用平均值 ， 修正系數(shù)值可按照下式計(jì)算 （4-6） 出口、進(jìn)口處的干度 2m mLFLG(V )LP1x (1)D2e ei ieix -x=x -xxe xi- l1 -摩擦壓降計(jì)算的經(jīng)驗(yàn)方法均相流模型 歐美采用的均相模型計(jì)算法和蘇聯(lián)采用的計(jì)算式方法略有不同采用均相模型計(jì)算 必須用一合適的摩擦阻力系數(shù) 蘇聯(lián)水動力計(jì)算標(biāo)準(zhǔn)采用單相流體的摩擦阻力系數(shù)，且認(rèn)為流動工況已進(jìn)

4、入阻力平方區(qū)， 值已經(jīng)與Re數(shù)無關(guān)，與相對粗糙度D/K有關(guān)其計(jì)算可按下列尼古拉茲計(jì)算式21=D4 lg3.7K（）（4-7）D及K分別為管子內(nèi)直徑及管壁粗糙度歐美則都采用勃拉休斯的光滑管計(jì)算式值，其公式為：0.250.3164=Re（4-8） l1 -摩擦壓降計(jì)算的經(jīng)驗(yàn)方法均相流模型 即 （4-9）按此法計(jì)算，當(dāng)干度x=1時，不能正確反映全部為氣體流過時的數(shù)值，因?yàn)樵诖擞?jì)算法中 的函數(shù)不能轉(zhuǎn)化為氣體雷諾數(shù)的函數(shù)。為了避免這一不足，有些作者采用一個平均 的兩相動力黏度 ，其值和干度x的關(guān)系應(yīng)能滿足：當(dāng)x=0時， 當(dāng)x=1時， 的條件。歐美在計(jì)算式4-8中的Re數(shù)時又有多種算法LLLLDuRe =

5、Lf(Re )LG l1 -摩擦壓降計(jì)算的經(jīng)驗(yàn)方法均相流模型 不同作者采用的 值也不相同麥卡達(dá)姆（Mcadams）采用 希奇蒂（Cichitti）采用 杜克勒（Dukler）采用 班可夫（Bankoff）采用 如果使用麥卡達(dá)姆的式進(jìn)行計(jì)算，則 值可按下式算得： （4-10）應(yīng)用式（4-10）算得的 值比按式（4-4）算得的低，較適宜用于高質(zhì)量流速工況。GL1x(1x)GLx(1x)GLmGLx(1x)GL+ 1- （）FP1GLLGF4GGP=1+x()1x()PoFP 氣液兩相完全分開流動不考慮兩相間相互作用單位管長上的氣、液兩相的摩擦壓力降相等，且等于兩相管流的摩擦壓力降10l1 -摩擦壓

6、降計(jì)算的經(jīng)驗(yàn)方法分相流模型 在1944-1947年間由Lockhart和Martinelli等提出，建立在分相模型流動基礎(chǔ)上假設(shè)(PFL)G(PFL)L(PFL) 11l1 -摩擦壓降計(jì)算的經(jīng)驗(yàn)方法分相流模型 推導(dǎo)-應(yīng)用單相流動的Darcy公式(PFL)GG1DGG22GGDGG2x22G2由于GGGGGGWV AVxWGAG 所以 （1）GGx GV222FL LLL2LLLP1VG (1x)()LLD2D 2(1)同理： （2）LL(1x)GV(1) 其中， 、 分別代表混合物中的氣相和混合物中的液相FGP()LFLP()L注：在英美文獻(xiàn)中，定義摩擦阻力系數(shù)f為:摩擦切應(yīng)力 ，而 （其中p

7、為管子周長）。根據(jù)受力平衡知上式等于 而 ，顯然 。、 分別為兩相流中的氣相和液相與管壁的摩擦阻力系數(shù)（指在其流動截面AG、AL中流速時）GL20vf2dPF pdz()AdPF0pD()dzpD244D0dz4fDv22dzdpFdzDv224 f 12l1 -摩擦壓降計(jì)算的經(jīng)驗(yàn)方法分相流模型 定義分相折算系數(shù)（摩擦阻力乘積因子）?，通過它將兩相混合物的摩擦壓降梯度同氣相或液相單獨(dú)在管內(nèi)流動時的摩擦壓降梯度聯(lián)系起來）分氣相折算系數(shù)： （3） 為假定氣相單獨(dú)流過流動截面A時的摩阻。分液相折算系數(shù)： （4） 為假定液相單獨(dú)流過流動截面A時的摩阻。2FFGG0PP() / ()LLFG0P()L2

8、FFLL0PP() / ()LL(PFL)L0 13l1 -摩擦壓降計(jì)算的經(jīng)驗(yàn)方法分相流模型 其中：(PFL)G0GO1DGVSG22GODG2x22GGSGGSGGwVAVxwGAGSGGxGV（5）（因?yàn)?, 所以 ）222L SLL0FL0L0LVP1G (1x)()LD2D2（6） 14l1 -摩擦壓降計(jì)算的經(jīng)驗(yàn)方法分相流模型 將（5）、（1）代入（3）式中得： （7）（較精確的，沒有不合理假設(shè)） 都與各自的Re數(shù)有關(guān)，并按通用的勃拉修斯公式計(jì)算有： （8） （注：n值取決于流態(tài)）其中: Re=2000-105時，n=0.25， （水力光滑管）；或Re=5000-200000時，n=0

9、.2， ） （9）（光滑管區(qū)的Blasius公式為： ）21GGG2G0D1()DGG0,nGGGtGV DC ()VtRe2000,n1.0,C64tC0.314tC0.186nSGG0tGVDC ()V0.250.3164Re 15l1 -摩擦壓降計(jì)算的經(jīng)驗(yàn)方法分相流模型 又因 則（ 所以， （10）（其n存在的前提是假定氣相流過AG時的流態(tài)和其單獨(dú)流過管道A時的流態(tài)相同，因而才認(rèn)為式（8），（9）中的t，n相等）同理可得： （11）2GGSGGAD(),VVADn 522Gn 522L(1) 又定義： （12） X2稱為Lockhart-Martinelli參數(shù) 所以， （13） （14

10、）16C2(1)n52C2G2/L211C45nC2(PFL)L0/(PFL)G0 17l1 -摩擦壓降計(jì)算的經(jīng)驗(yàn)方法分相流模型 G2= C45-n+15-n2L2= C4n-5+15-n22G2L流態(tài)液相Re氣相Re紊流（液）紊流（氣）（tt）10001000紊流粘性（tv）10001000粘性紊流（vt）10001000粘性粘性（vv）10001000注: 把Re數(shù)等于1000作為層流、紊流的界限，是因?yàn)樯鲜鯮e中的速度是折算速度，因此一種相的Re數(shù)有效值會因?yàn)榱硪幌嗟拇嬖诙龃?，而?shí)際的各相雷諾數(shù)都要大。 18l1 -摩擦壓降計(jì)算的經(jīng)驗(yàn)方法分相流模型 與X的關(guān)系如圖所示LG ， 19l1

11、 -摩擦壓降計(jì)算的經(jīng)驗(yàn)方法分相流模型 在LockhartMartinelli方法計(jì)算 時，先計(jì)算出 ，求得X 值。然后再在線算圖上（或計(jì)算）求得 或者 ，最后再由 或 求得 。FPFG0FL0P / L)P /L)（、（2G2L2G2LFP該方法適用于雙組分的氣液兩相流在低壓時摩擦阻力計(jì)算，因其計(jì)算數(shù)據(jù)是建立在低壓的氣液流動基礎(chǔ)之上的。為了適用于汽水混合物的摩擦阻力壓力降的計(jì)算，Martinelli-Nelson對此方法進(jìn)行了改進(jìn)。 20l1 -摩擦壓降計(jì)算的經(jīng)驗(yàn)方法分相流模型 假設(shè)汽水兩相分開流動時都呈紊流狀態(tài)，同時利用常壓下的空氣水混合物試驗(yàn)數(shù)據(jù)和高壓汽水混合物的試驗(yàn)數(shù)據(jù)建立了 的關(guān)系曲線

12、。定義 ，全液相折算系數(shù)兩相壓降與假設(shè)汽水混合物全部為液相時的摩擦壓力降（以總的質(zhì)量流向相等為原則）之比 為假設(shè)兩相混合物全部為液相時的摩擦阻力系數(shù)2FL00Pf( x,p)P20oLL GPD 222L LFFLLL2LLLLVLG (1x)P(P )D2D2(1)o2FL00Pf( x,p)P 21l1 -摩擦壓降計(jì)算的經(jīng)驗(yàn)方法分相流模型 與x，P的關(guān)系如下2L0圖A圖B 22l1 -摩擦壓降計(jì)算的經(jīng)驗(yàn)方法分相流模型 按圖A根據(jù)干度x及壓力P查出 ，再乘上算出的 值，即可得出不受熱管中的 值。對于進(jìn)口處干度x=0，出口處x=xe的受熱管，可按出口干度xe及P值在圖B中查出自x=0到x=xe

13、的 平均值 ，然后按同法求得 值。一般認(rèn)為此方法適用于低質(zhì)量流速工況。2L00PFP2L0FP 23l1 -摩擦壓降計(jì)算的經(jīng)驗(yàn)方法分相流模型 奇斯霍姆對兩相流動摩擦阻力壓降進(jìn)行了許多研究，得出了 與X值的關(guān)系為或者 其中 式中C為系數(shù)，可按下式確定： 式中 系數(shù) C2系數(shù) 、 氣體及液體的比容，m3/kg. 系數(shù)C是壓力P和干度x以及質(zhì)量流速的函數(shù)2L2FL2L0PC11PXXG21CCC2X2(PFL)L0/(PFL)G00.50.50.5GLGL2GLGC(C)()()()GL 24l1 -摩擦壓降計(jì)算的經(jīng)驗(yàn)方法分相流模型 Chisholm推薦C值如下LiquidGasSubscriptC

14、TurbulentViscousTurbulentViscousTurbulentTurbulentViscousViscousttvttvvv2012105對壓力P3MPa的汽水混合物，C值的計(jì)算A:當(dāng)質(zhì)量流速 時對光滑管： C值按P31頁公式計(jì)算，此時 。 對粗糙管： 此時 。NoImage2u* 2000kg/ ms2u*0.75,Cu2u*1500kg / ms2u*1,Cu 25l1 -摩擦壓降計(jì)算的經(jīng)驗(yàn)方法分相流模型 B:當(dāng)質(zhì)量光流速 時 對光滑管（ 和粗糙管 式中 ， ， 值可按下法計(jì)算； C值仍然按照P31公式計(jì)算。對于粗糙管 ；對光滑管 。對于粗糙管n=0，對于光滑管n=0.

15、25.由上所述可見，奇斯霍姆計(jì)算法的計(jì)算過程較麻煩。但是計(jì)算結(jié)果和試驗(yàn)偏差較小。uu*2u2000kg / ms)2( u1500kg / ms)2L2C11XXC22C1C11 / 1TTTT0.50.5GLLGvvC()()vv( 2 n )/2n 20.5LLGGxvT()()()1xv10.752u*Cu 26l1 -摩擦壓降計(jì)算的經(jīng)驗(yàn)方法分相流模型 Friedel把一個有25000個點(diǎn)的數(shù)據(jù)庫用他本人的關(guān)系式以及當(dāng)時存在的其他公式進(jìn)行了比較，得出全液相折算系數(shù)2L0FF00.0450.0353.24FH(P /L) / (P /L)BFrWe其中 分別為在總質(zhì)量流速相同的情況下，具有

16、氣相或液相物性的單相流體在管內(nèi)流動下的摩擦阻力系數(shù)，即全氣相、全液相時的摩擦阻力系數(shù)） 適用于垂直向上流動與水平流動（應(yīng)用于垂直向下流動的關(guān)系式稍有不同，其對單組分流動的標(biāo)準(zhǔn)偏差約為30%，對雙組分流動標(biāo)準(zhǔn)偏差約為4050%，可見誤差仍然相當(dāng)大）其中假設(shè)兩相混合物的總質(zhì)量全部以液體流過時的壓降22LoGGoLfB(1x )xfoGoLf, f0.760.24Fx(1x)0.910.190.7GGLGLLH()()(1)222mmGG DFr,WegD 1mGLx1x() 27l1 -摩擦壓降計(jì)算的經(jīng)驗(yàn)方法分相流模型 u當(dāng) 時，應(yīng)用Friedel關(guān)系式u當(dāng) 時，應(yīng)用Chisholm關(guān)系式u當(dāng) 時

17、，應(yīng)采用Martinelli關(guān)系式在“傳熱和流動服務(wù)中心”的專利數(shù)據(jù)庫基礎(chǔ)上得到的計(jì)算結(jié)果對以上幾個分相模型公式的應(yīng)用建議：（通過試驗(yàn)比較給出的建議）LG/1000LG/1000,G100LG/1000,G100 28l1 -摩擦壓降計(jì)算的經(jīng)驗(yàn)方法我國摩擦壓力降的計(jì)算方法作者等應(yīng)用能譜理論和混合長度理論相結(jié)合的方法導(dǎo)得一種計(jì)算管內(nèi)氣液兩相流動的摩擦阻力壓力降的方法設(shè)有一絕熱的均勻氣液兩相紊流沿水平管道流動。氣液兩相流沿管道中心線方向的速度為u，在垂直中心線方向的速度為v。令u為速度u的脈動分量，v為速度v的脈動分量，為密度的脈動分量。則可得到氣液兩相流的摩擦切應(yīng)力如下： （4-1）式中 ， 均

18、為統(tǒng)計(jì)平均值；y為從管子內(nèi)壁算起的距離。由混合長度理論可知： （4-2） （4-3）式中 ， 分別為速度分布及密度分布的混合長度。du= uvu vdy u vvdldyuduuvldy 29l1 -摩擦壓降計(jì)算的經(jīng)驗(yàn)方法我國摩擦壓力降的計(jì)算方法 根據(jù)能譜理論，前述兩種速度脈動分量之間的關(guān)系為： （4-4） 式中 k紊流流動的外尺寸； 無窮大 譜密度譜密度可用下式表示： （4-5）式中s任意值； b常數(shù)； 總能量； v運(yùn)動黏度。將式（4-5）代入 （4-4），可得 （4-6）9/45 1/4(1 s )/ s7(P)b( v )16.9x22k2uv( P)dP3( P )6229/45 1/

19、4(1 s )/ s7k2Akuvb( v )16.9xdx36 30l1 -摩擦壓降計(jì)算的經(jīng)驗(yàn)方法我國摩擦壓力降的計(jì)算方法式中令 。將公式（4-6）代入（4-1），并令 ， 可得 （4-7）速度u和密度之間的關(guān)系可用下式表明： （4-8） 上式中的A和B均為常數(shù)。應(yīng)用下式邊界條件：y=0處（管壁上） ， 處（管子中心）， ，并將上述條件代入式（4-8），可以得到下列密度和速度之間的關(guān)系式： （4-9） 上式中 及 分別表示液體密度及管子中心處氣液兩相流密度。9/45 1/4(1 s )/ s2Ab( v )16.93ulll62Akdd(l u)6dydyAuBL,u0;yRcc,uuc c

20、LLccuu(uu )Lc 31l1 -摩擦壓降計(jì)算的經(jīng)驗(yàn)方法我國摩擦壓力降的計(jì)算方法令 （4-10） 則上式（4-9）可改寫成為下列形式： （4-11） （4-12）設(shè)兩相流動混合長度和單相流動的相同，則可應(yīng)用下式計(jì)算混合長度 （4-13） 式中K1及K2均為無因次混合長度常數(shù)。聯(lián)立解式（4-7）、（4-11）、（4-12）及（4-13），并令 （4-14） 可得 （4-15） 式中 并 令 Lc cLcu (1)u()c c2Lcudud1dy() dy 12ylK y1exp()KLc1c cLaKu6222222Akadyy1()()1exp() ()()6dyKadK1dyLLLLu

21、yy/ ();uL22LLKK/ ();Luu/ 32l1 -摩擦壓降計(jì)算的經(jīng)驗(yàn)方法我國摩擦壓力降的計(jì)算方法進(jìn)一步推導(dǎo)可得 （4-16）式（4-16）的邊界條件為 。應(yīng)用這一邊界條件，且在a，K1， 和 給定時，通過逐步積分可得出 和y之間的函數(shù)關(guān)系。對于管子半徑為R的圓管，平均混合物物密度 可由下式求得： （4-17） 式（4-17）也可以改寫成 （4-18） 式中 （4-19）式中 ReL按總流量和液體物性計(jì)算的雷諾數(shù)22222212dyyad(1)a() 1exp() ()KdyKdy1,y02Kum1m0yy2(1) d()RRRm022( Ry )dyR1L01aKRRe21 33l

22、1 -摩擦壓降計(jì)算的經(jīng)驗(yàn)方法我國摩擦壓力降的計(jì)算方法 （4-20） 質(zhì)量流速 （4-21）速度分布 （4-22） 將式 （4-22）及前述的 等計(jì)算式代入式（4-21），可得 （4-23）或改寫成 （4-24）當(dāng)總流量為液體時，摩擦切應(yīng)力可按下式計(jì)算： （4-25） 兩相摩擦阻力壓力降 和單相液體摩擦阻力壓力降 之比可寫作為 （4-26） 或者 （4-27）當(dāng)常數(shù)K1， 和 已知時，應(yīng)用式（4-16）、（4-18） （4-19）和式（4-27）可得出兩相摩擦阻力壓力降 的解。萊維（Levy）的研究表明：K1=0.4， =26， 。因而根據(jù)上述4個方程式可解得 和 的值如下圖所示LLGDRe10

23、yyG2(1) ud()RR1(1)u(K a)2y ,K ,uLLm1/G(1)K a 22210L2m2K a(u/ 2)(1)20LLLu142 FP0P222F1m0LLP8K a(1)P222L1mPF8K a (1)Po2KFP2K1LF0P /P m 34l1 -摩擦壓降計(jì)算的經(jīng)驗(yàn)方法我國摩擦壓力降的計(jì)算方法 35l1 -摩擦壓降計(jì)算的經(jīng)驗(yàn)方法我國摩擦壓力降的計(jì)算方法若知 可按 值在圖上查得 值。 可按照多列查爾（Dolezal）計(jì)算法 算得， 為單相液體的摩擦阻力系數(shù)，可按式（18）計(jì)算,可求出 。 ，其中 按林宗虎氣液兩相流和沸騰傳熱書中公式（4-28）計(jì)算，因而 （4-28

24、）式中的截面含氣率 值可按照按林宗虎氣液兩相流和沸騰傳熱書中第四章中式（4-4）計(jì)算.mLF0P /P LRe2LoLPouD2LFPmmL/mLGmL(1) 36l1 -摩擦壓降計(jì)算的經(jīng)驗(yàn)方法我國摩擦壓力降的計(jì)算方法陳立勛計(jì)算采用的型式和蘇聯(lián)61年鍋爐水力計(jì)算標(biāo)準(zhǔn)方法建議使用的公式一樣； （4-29）只是校正系數(shù)的值計(jì)算方法有所不同。陳立勛以林宗虎氣液兩相流和沸騰傳熱書中圖3-2的試驗(yàn)曲線為基礎(chǔ)得出的這些試驗(yàn)曲線校正系數(shù)的近似數(shù)學(xué)表達(dá)式如下：當(dāng) 1500kg/m2s時， 值按照下式計(jì)算 （4-30） 1500kg/m2s時， 值按照下式計(jì)算 （4-31）式中x蒸汽干度。按照此法計(jì)算結(jié)果與前述

25、奇斯霍姆計(jì)算式的計(jì)算結(jié)果相近，但是計(jì)算過程要簡便得多。2m mLFLG(V )LP1x(1)D2LGLG1500 x(1x)(1)u11x()1LGLG1500 x(1x)(1)u11(1x)()1uu 37l1 -摩擦壓降計(jì)算的經(jīng)驗(yàn)方法我國摩擦壓力降的計(jì)算方法我國電站鍋爐水動力計(jì)算方法中計(jì)算汽水混合物兩相摩擦阻力壓力降的計(jì)算式的計(jì)算型式和式（4-29）型式相同，但是校正系數(shù)值的計(jì)算方法主要采用（4-30）和式（4-31）的型式。根據(jù)林宗虎的研究，如果將式（4-30）和式（4-31）所適用的質(zhì)量流速分界值 改寫成 ，并將式（4-30）和式（4-31）中的1500改寫成1000，則計(jì)算所得的摩擦

26、阻力壓力降值和其他實(shí)驗(yàn)值以及運(yùn)行數(shù)據(jù)吻合的更好。壓力為4.5-10.5MPa，質(zhì)量流速為500-2700Kg/m2s，熱流密度為0-0.57MW/m2,進(jìn)口干度為0-0.81的試驗(yàn)條件所進(jìn)行的水平管中汽水混合物摩擦阻力壓力降試驗(yàn)證實(shí)了這一點(diǎn)。2u1500kg / ms2u1000kg / ms 38l1 -摩擦壓降計(jì)算的經(jīng)驗(yàn)方法我國摩擦壓力降的計(jì)算方法 在我國電站鍋爐水動力計(jì)算方法中，汽水混合物在水平、垂直以及傾斜管中流動時，摩擦阻力壓力降按下式計(jì)算： （4-32） 式中 -汽-水混合物的質(zhì)量流速，kg/m2s，也可用 代替，因?yàn)樵诘戎睆街惫苤懈鹘孛娴馁|(zhì)量流速相等；X管內(nèi)平均質(zhì)量含氣率； 摩擦

27、阻力系數(shù) 2m mLFLG(u )LP1x(1)D2m muL 0u 39l1 -摩擦壓降計(jì)算的經(jīng)驗(yàn)方法我國摩擦壓力降的計(jì)算方法摩擦壓降校正系數(shù)按以下方法計(jì)算： （1） ；（2）時， （4-33）（3） 時， （4-34）2u1000kg / ms12u1000kg / msLGLG1000 x(1x)(1)u11x(1)2u1000kg / msLGLG1000 x(1x)(1)u11(1x)(1) 40l1 -摩擦壓降計(jì)算的經(jīng)驗(yàn)方法我國摩擦壓力降的計(jì)算方法n值也可以按照線算圖查得。n對于不受熱管， 值按照林宗虎氣液兩相流和沸騰傳熱書中圖3-12查得；n對于受熱管 值按林一書上圖3-13查得

28、。n對于管子入口處工質(zhì)已經(jīng)是汽水混合物的受熱管，當(dāng)出口干度和進(jìn)口干度之差xe-xi 0.1時， 值按下式求得： （4-35） 式中 及 分別按xe及xi 在林宗虎氣液兩相流和沸騰傳熱書中圖3-13查出。e ei ieixx=xxei 41l1 -摩擦壓降計(jì)算的經(jīng)驗(yàn)方法我國摩擦壓力降的計(jì)算方法n當(dāng) xe-xi 0.1時， 值按下式計(jì)算： （4-35）式中， 值為按 查出的 值， 及 則分別按xe及xi在林宗虎氣液兩相流和沸騰傳熱書中圖3-12中查出。ei1=46( xexi ) / 2ei 42l2 -依據(jù)流型的摩擦壓降的計(jì)算對氣液兩相流摩擦阻力壓力降的計(jì)算至今已做了不少研究工作，提出了許多純經(jīng)

29、驗(yàn)的或半理論半經(jīng)驗(yàn)的計(jì)算式或計(jì)算方法。但大多數(shù)都存在著適用范圍較窄和考慮影響因素不夠全面的缺點(diǎn)。對基于某種流型的具體特點(diǎn)，建立在理論分析基礎(chǔ)上的摩擦阻力壓力降計(jì)算方法，也曾進(jìn)行了一些研究。Mandhane曾經(jīng)收集了近萬個氣液兩相流體在水平管中流動的試驗(yàn)數(shù)據(jù)，并和16個摩擦阻力壓力降計(jì)算式的結(jié)果進(jìn)行了比較，表明根據(jù)管中流型應(yīng)用相應(yīng)計(jì)算式的計(jì)算結(jié)果最精確。本節(jié)就將根據(jù)不同流型的具體特點(diǎn)，分別建立相應(yīng)的摩擦阻力壓力降計(jì)算式。 在泡狀流中氣相彌散在連續(xù)的液相中，這時摩擦阻力的產(chǎn)生主要是由以速度v運(yùn)動的混合物中的液相造成的，因而泡狀流的摩擦阻力壓力降可表示為： 式中 fM為兩相混合的摩擦阻力系數(shù)，M為兩

30、相混合物密度，G為混合物的質(zhì)量流速. 22fMMddzDfPGSLSGLGGVVl2 -依據(jù)流型的摩擦壓降的計(jì)算 Beattie & Whalley通過與大量的水平或垂直管中的壓降數(shù)據(jù)的分析比較認(rèn)為，對泡狀流如fM、M按以下方法計(jì)算，將獲得較好的結(jié)果： 式中 /D為管壁相對糙度，ReM為兩相混合物的雷諾數(shù)。 16/2100ee13.48 42( /) 9.35/()2100logee10fRRMMMDfRRMMMfM 當(dāng) 當(dāng)eMMDGRl2 -依據(jù)流型的摩擦壓降的計(jì)算 混合物密度M按均相模型計(jì)算如下; 兩相混合物的動力粘度M可由下式計(jì)算; 11MgLxx11MGLxxl2 -依據(jù)流型的

31、摩擦壓降的計(jì)算 間歇流中，每一個流動單元的摩擦阻力壓降由兩部分組成：液彈段摩擦阻力壓力降和Taylor氣泡段摩擦阻力壓力降，即： 2222SFMTsMSGGfGdDffVVLLPDl2 -依據(jù)流型的摩擦壓降的計(jì)算 式中fs，fG分別為液彈的摩擦阻力系數(shù)和Taylor氣泡段的摩擦阻力系數(shù)DG為Taylor氣泡的當(dāng)量直徑fs，fG可分別根據(jù)ReMS，由Chen式求得 其中， 液彈段混合物密度 液彈段混合物粘度MSMSMMSVReGGTGGVDRe()GLSLSLMS1(1)LSLSGLMSl2 -依據(jù)流型的摩擦壓降的計(jì)算 Taylor氣泡的當(dāng)量直徑DG對垂直管為： DGD2LL為Taylor氣泡段

32、液膜平均厚度， 因此()LFLD112211(11)GLFDDDLFl2 -依據(jù)流型的摩擦壓降的計(jì)算 對水平或微傾斜管;SG,Si可由LF按以下方法求得;由于于是 可由該式中求得值則無量綱液位高度iGGGSSAD4p2sinLF2sinLFp22cos1Lhl2 -依據(jù)流型的摩擦壓降的計(jì)算 在式中，間歇流的平均摩擦阻力壓力降梯度為：在以下計(jì)算中，所有的間歇流特性參數(shù)，如LS、LF 、LS、LU、VT均等可由第四章中相應(yīng)的間歇流特行參數(shù)的封閉方程組中求得。1cos (21)GLDSh21(21)GLDSh(1)GLFAA2222MTfsMSSGGFUGUddzDffVVPLLLDLl2 -依據(jù)流

33、型的摩擦壓降的計(jì)算 傾角為的氣芯和液膜的動量方程為： sin0GGiiGdpdzSAAsin0GLBBiiLdpdzSSAAl2 -依據(jù)流型的摩擦壓降的計(jì)算以上兩式相加消去iSi，并整理得： 又因AL/A=1-且M= G+(1+ )L,于是上式可寫成: ()() sinGBBLLGAAdpgdzAAAS()sinBBMdpgdzAS 總的壓力降由兩部分組成組成，即摩擦阻力壓力降（右式第一項(xiàng)）和重位壓力降（右邊第二項(xiàng)），加速壓力降在上面的模型中被忽略。于是環(huán)狀流的摩擦阻力壓力降為： B為液膜與管壁的摩擦切應(yīng)力，SB為液膜濕周,A為管道流通截面積。()BBdpdzASl2 -依據(jù)流型的摩擦壓降的計(jì)

34、算 將以上三式代入（5-30）中得： 22222LSLLLBLFLFLVVffBSDp24ADp22244()22(1)SLSLLFLLFLLdpdzDDffVV 液膜摩擦阻力系數(shù)式中DL為液膜的當(dāng)量直徑，由第（3.3.3）節(jié)中的分析知：于是 當(dāng)ReLF20002500kg/m2s）采用均相模型較為適宜。（兩相湍動度大，兩相分布更為均勻）流動方向 l3 -影響氣液兩相流體摩擦阻力壓力降的主要因素 研究較少，通常管壁粗糙物突出貼壁液膜較多時，兩相流體流動沿粗糙管流動時的摩擦阻力系數(shù)約比沿光滑管流動時大一倍。但當(dāng)液膜能蓋住管壁粗糙度的凸尖時，則光滑管和粗糙管的摩擦阻力壓力降相同。 不同的流動型態(tài)造

35、成兩相流動規(guī)律的不同，因此根據(jù)不同的流動形態(tài)的具體特點(diǎn)建立相應(yīng)的數(shù)學(xué)物理模型，從而求得相應(yīng)的摩擦阻力壓力降，在本質(zhì)上能接近真實(shí)的流動狀態(tài)。6.管內(nèi)粗糙度 D l3 -影響氣液兩相流體摩擦阻力壓力降的主要因素該方法也是目前正在廣泛研究的課題之一。無非是 等的影響，歸結(jié)起來很大程度上都是由于流型的變化，而使得摩擦阻力發(fā)生變化。統(tǒng)一的經(jīng)驗(yàn)的方法是需要綜合考慮各種因素。顯然是力不從心。建立依據(jù)流型的壓力降計(jì)算模型和準(zhǔn)確確立管內(nèi)的流型的流型的轉(zhuǎn)換界限是緊密相聯(lián)的。v DP xD、 、目前這方法研究主要有：Bendiksen（Slug Flow），陳宣政（垂直上升管），王樹眾(各種布置狀況下)、Kokal

36、 和Stanislav（1989）水平管等。 l作業(yè)直徑D=5.08cm管子，P=180bar，進(jìn)口流量M=2.14kg/s,進(jìn)口為飽和水， 粗糙管 ，出口干度 ，管長100m，求 。 分別用MN法、Chisholm方法（經(jīng)驗(yàn)的C公式）、蘇聯(lián)78年計(jì)算標(biāo)準(zhǔn)、我國水動力計(jì)算方法。（參考：陳之航、曹柏林、趙在三，“氣液雙相流流動和傳熱”，機(jī)械工業(yè)出版社）=0.002Dex =0.1825FP 在實(shí)際工程應(yīng)用中，在多相管流的壓降梯度中，加速壓降所占的份額很低。某些研究者導(dǎo)出多相管流壓降關(guān)系式時，考慮了加速壓降，有的則忽略加速壓降。本章對加速壓降的計(jì)算和特點(diǎn)進(jìn)行了討論。 氣液兩相流體在管路中流動時，加

37、速壓力降的產(chǎn)生通常由以下兩個原因造成：（1）由于加熱、冷卻以及壓力變化，使兩相混合物的組成或流速發(fā)生變化；（2）由于管路流通截面積發(fā)生變化，從而引起兩相流速發(fā)生變化。當(dāng)兩相混合物在等截面直管中流動時，第二部分的加速壓力降就消失了。 在通常的氣液兩相流操作速度范圍內(nèi)，由于兩相流速較低，加速壓力降常忽略不計(jì)。當(dāng)兩相混合物流速較高或兩相混合物由于壓力、溫度等的變化而凝析或氣化，造成兩相流速發(fā)生較大變化時，就可能會導(dǎo)致較大的加速壓力降。這時，必須計(jì)及加速壓力降的影響。對加速壓力降的計(jì)算基本上有兩種方法:1、基于氣液兩相完全分開流動的分相模型進(jìn)行；2、根據(jù)均相模型進(jìn)行推導(dǎo)。由于后者不考慮兩相間的相對滑動

38、，因而誤差較大。 本節(jié)將采用分相模型計(jì)算方法來導(dǎo)得氣液兩相流的加速壓力降。 由兩相流的動量方程知穩(wěn)定流動的氣液兩相流的加速壓力降為: 式中，G為兩相混合物的質(zhì)量流速，流動過程中質(zhì)量守恒。2221(1)(1)aLGddxxAdzA dzPG 222222(1)(1)(1)(1)LGLGxxdGxxdAGdzAdz 在上式中:第一項(xiàng)為由于兩相混合物組成（x）或壓力變化（P）而引起的壓力降;第二項(xiàng)為由于管路流通截面積發(fā)生變化而造成的加速壓力降。若在等截面直管中，氣液兩相混合物從位置Z1流到Z2，則上式中第二項(xiàng)為零，于是加速度壓力降梯度為： 122222211222211(1)(1)(1)(1)()(

39、1)(1)LGaLGLGxxddPxxxxGGdzdz 管路內(nèi)氣液兩相流的重位壓力降為 其中，為流動方向和水平線的夾角，向上流動時為正。a為相應(yīng)流型下的截面含氣率。在間歇流中，由于每個流動單元總是由氣泡段和液彈段兩部分組成，環(huán)繞Taylor氣泡的是貼壁液膜，而Taylor氣泡內(nèi)壓力沿長度方向基本上是恒定的，因而Taylor氣泡造成的重位壓力降可以忽略不計(jì)，這樣段塞流的重位壓力梯度可表示為：(1)singGLdPgdz L( sin)( sin)singSFLFMSUUSFLFMSLUUdggdzgPLLLLLLLL 式中，MS為液彈中兩相混合物的密度，可按下式計(jì)算： 試驗(yàn)表明，重位壓力降是管內(nèi)

40、氣液兩相流豎直流動時總壓降的重要組成部分，因此準(zhǔn)確計(jì)算重位壓力降，對提高管路內(nèi)兩相流壓力降的計(jì)算精度具有重要意義。而重位壓力降計(jì)算的關(guān)鍵就在于對截面含氣率（或持液率）的準(zhǔn)確計(jì)算。(1)LSLSMSLG 目前，對氣液兩相流流過局部阻力件時，局部阻力壓力降的研究還很不夠，還主要是依據(jù)試驗(yàn)確定的經(jīng)驗(yàn)的方法。通常通過人工查表求得各種經(jīng)驗(yàn)系數(shù)，然后再進(jìn)行計(jì)算。下面重點(diǎn)介紹彎頭、節(jié)流式測量元件、突擴(kuò)、突縮接頭、三通、閥門的局部阻力壓力降。 汽液兩相流通過彎頭時的阻力可分為兩部分：在彎頭內(nèi)部產(chǎn)生的，是由于兩相流體通過彎頭時出現(xiàn)二次流和流場變化而引起的；由于汽液兩相滑動比發(fā)生了變化，在下游又要恢復(fù)到直管內(nèi)的滑

41、動比所引起的動量變化。 Chisholm提出了一種半理論半經(jīng)驗(yàn)的方法來計(jì)算流過彎頭的局部阻力壓力降。 對第一部分壓力降按均相模型有：PBL0為假設(shè)兩相流體全部為液體時，流過彎頭的摩擦壓力降：式中, 單相流體摩擦阻力系數(shù) L/D彎頭當(dāng)量長度 G兩相混合物的質(zhì)量流速，kg/(m2s)1011LBBLGPPx()22022BLLLLVGLPDD 對單相流體的摩擦阻力系數(shù)，Churchill推薦在整個流動范圍內(nèi)（在全Re數(shù)范圍內(nèi)）可按下式計(jì)算： 其中 B=(37530/Re)16()1/12123/2818eRAB()()160.912.45ln7 /0.27/eARD 彎頭當(dāng)量長度LD，對于 90彎

42、頭可根據(jù)彎頭相對彎曲半徑RD求得。通過數(shù)據(jù)擬合 LD和 RD的關(guān)系如下： 對第二部分壓力降按分相模型計(jì)算： 式中MFL0為假定兩相混合物全部為液體時單位面積的動量MFL0 =G2/L()()8765432/1.249634/5.886433/0.117155(/)1.28073(/)8.36847(/)33.3258(/)79.1471(/)102.565(/)67.7579LDERDERDRDRDRDRDRDRD201 () (1)LBLGPMFB xx 系數(shù)B反映了流過彎頭過程中，出口處和進(jìn)口處的滑動比的變化于是兩相流體流過彎頭時的局部阻力為：經(jīng)整理得 11.12/BSR D12BBBPP

43、P201()(1)LBBLGPPBxx 式中系數(shù)B對90彎頭有： 當(dāng)彎曲角度小于90o時，取BB90當(dāng)彎曲角度大于90o時，其中 單相液體流過90度彎頭的阻力系數(shù) 單相液體時，彎曲角度為 的彎頭阻力系數(shù) ()9090()2.212/BLDBR D ()90()90()11BLDBLDBB ()90BLD()90BLDLD()BLD()0()BLDkB 其中 0 彎頭的原始阻力系數(shù)，與其形狀和相對彎曲半徑有關(guān)； k考慮管壁粗糙度的影響系數(shù)； (k0)對緩轉(zhuǎn)彎與焊接彎頭可按下式給出（經(jīng)數(shù)據(jù)擬合）：B是與彎曲角度有關(guān)的修正系數(shù)，經(jīng)數(shù)據(jù)擬合可按下式計(jì)算：式中的單位為度 ()3203210.9596(

44、/)29.6667( /)27.8393( /)9.38418 0.00106273( /)0.0157747( /)0.0981203( /)kR DR DR DR DR DR D321.9547378.8757950.01743150.00047619BEE 應(yīng)該注意的是，上式僅用于質(zhì)量含氣率x0.04的情形。因?yàn)樵趚0.04時彎頭的局部阻力計(jì)算如下方法:式中, PL0假定兩相全部為液體時，流過彎頭的局部阻力壓力降,2(1)3BLLOGPxxP()22LOBLDLGP -單相液體流過彎曲角度為 的彎頭時的局部阻力系數(shù)；相對彎曲半徑R/D的修正系數(shù)，可由下式查出經(jīng)多項(xiàng)式數(shù)據(jù)擬合可按下式計(jì)算(

45、)BLDR/D3456700.71.31.82.020/30.0785714(/)1.29571(/) 3.19714/3R DR DR DR Da當(dāng)當(dāng) 孔板、文丘利管和噴嘴是節(jié)流式測量元件，可以測流量和干度。因?yàn)樗鼈兘Y(jié)構(gòu)簡單、無轉(zhuǎn)動部件、使用簡便、運(yùn)行可靠且有足夠多的精度，所以主要用作兩相測量設(shè)備。它們的工作原理相同，下面以孔板為例。1）對單相流體： (kg/s) Y 流體膨脹系數(shù)，若認(rèn)為流體不可壓縮，則 = 1 孔板的膨脹系數(shù) 孔口直徑和管子內(nèi)徑之比， A 孔口截面積， c 流出系數(shù)，discharge coefficient，其值與取壓方式（角接取壓、法蘭取壓、徑距取壓等）、Re數(shù)等有關(guān)

46、 孔板前的流體密度PAcYW214 阻力損失 孔板前后壓力差 動能變化造成的壓力變化定義：流量系數(shù) 對不可壓縮流體且忽略孔板熱膨脹時， P41cPAW2 2）對兩相流體，流過孔板時的壓力降（流過孔板的流動很復(fù)雜，單純的均相模型和分相模型誤差很大）。自50年代以來，發(fā)表了不少計(jì)算其壓降的計(jì)算式，如Murdock、James、Bizon、Collins、林、Chisholm、Smith和均相、分相模型公式等，但大多數(shù)誤差較大，較好的計(jì)算公式有： 均相模型 分相模型 Murdock計(jì)算式 Chisholm計(jì)算式 林宗虎計(jì)算式 下面重點(diǎn)介紹Chisholm計(jì)算式和林宗虎計(jì)算式。 假定兩相流經(jīng)孔板時作分享流動，應(yīng)用動量守恒方程可得出 兩相流中液相單獨(dú)流過孔板時的壓力降，可根據(jù)其流量 由單相的孔板公式求得 兩相混合物流過孔板時的壓力降 Chisholm（1977）建議：系數(shù) 211ZZcPPLOTPLOPLWPAWLOL2TPP()LGGOLOxxPPZ2221/時，當(dāng)時，當(dāng)1Z1Z2/12/14/14/1mGGmLGGLc 其中，c和壓力及滑動比有關(guān) 該計(jì)算方法，當(dāng) x0.1， P=17MPa時，誤差為20%左右，其計(jì)算值和Murdock、James的試驗(yàn)不符合，但和其他研究