用VB編寫氣液兩相管流Beggs-Brill模型（課堂PPT）

1、氣液兩相管流氣液兩相管流Beggs-Brill模型模型史智慧2014年4月匯報提綱匯報提綱模型基本方程模型基本方程Beggs-Brill方法是可用于水平、垂直和任意傾斜氣液兩相管流動計算的方法，是1973年，Beggs和Brill根據(jù)在長15m，直徑1inch和1.5inch聚炳烯管中，用空氣和水進行實驗的基礎(chǔ)上提出的。氣體流量00.098 m3/s 液體流量00.0019 m3/s持液率00.87 m3/m3系統(tǒng)壓力241655 kpa(絕對壓力)壓力梯度018 kPa/m傾斜度-900+900流 型水平管流動的全部流型實驗參數(shù)范圍實驗參數(shù)范圍基本方程模型基本方程模型基本方程假設(shè)條件:假設(shè)氣

2、液混合物既未對外做功，也未受外界功，則單位質(zhì)量氣液兩相管流的壓力降消耗于位差、摩擦和加速度引起的壓力消耗?？倝毫μ荻确匠蹋簆vvDvsingzpogmmmmm12+=2模型基本方程模型基本方程(1)位差壓力梯度：消耗于混合物靜水壓頭的壓力梯度。sin1singHHgdZdpLgLL位差(2)摩擦壓力梯度：克服管壁流動阻力消耗的壓力梯度。(3)加速度壓力梯度：由于動能變化而消耗的壓力梯度。忽略液體壓縮性、考慮到氣體質(zhì)量流速變化遠(yuǎn)遠(yuǎn)小于氣體密度變化，則：vDAGDvdZdp2/22摩擦dZdppvvdZdpsg加速度Beggs - Brill 兩相水平管流型圖模型基本方程模型基本方程分離流分層流波

3、狀流環(huán)狀流間歇流團狀流段塞流分散流泡 流霧 流匯報提綱匯報提綱流型分布圖及流型判別式流型分布圖及流型判別式Beggs - Brill方法的流型分布圖圖中L1、L2、L3和L4和為四個流型區(qū)的分隔線302011316=.EL73361310=.E.L73361450=.E.L468421200092520=.E.L流型分布圖及流型判別式流型分布圖及流型判別式Beggs-Brill法流型判別條件匯報提綱匯報提綱持液率及混合物密度確定一、持液率Beggs-Brill方法計算傾斜管流時首先按水平管計算,然后進行傾斜校正。)(H)(Hll0=其中水平管持液率：( )cFrblNaEH1=0流 型a b

4、c分 離 流間 歇 流分 散 流0.98 0.4846 0.08680.845 0.5351 0.01731.065 0.5929 0.0609持液率及混合物密度確定不同E1下的傾斜校正系數(shù)圖傾斜校正系數(shù)與傾斜角、無滑脫持液率、弗洛德數(shù)及液體速度數(shù)有關(guān)。持液率及混合物密度確定傾斜校正系數(shù)的相關(guān)式為：11813183Csin( .)sin ( .) 103 . C對于垂直管：gFrfvleNNEdEC)()()(ln)1 (11其中：系數(shù)C與無滑脫持液率 、弗洛德數(shù)和液相速度數(shù)有關(guān)。流 型a b c分 離 流間 歇 流分 散 流0.98 0.4846 0.08680.845 0.5351 0.0

5、1731.065 0.5929 0.0609持液率及混合物密度確定對于過渡流型，先分別計算分離流和間歇流的持液率，之后采用內(nèi)插法確定渡流型的持液率。)()()(間歇分離LLLBHAHHALLLNBFr1232233LLNLAFr二、混合物密度利用持液率計算流動條件下混合物實際密度：)1 (lgllmHH匯報提綱匯報提綱阻力系數(shù)確定 氣液兩相流阻力系數(shù)與無滑脫氣液兩相流阻力系數(shù)的比值與持液率和無滑脫持液率(入口體積含液率)之間的關(guān)系：沿程阻力系數(shù)：1se2Re.1Re.2lg4.5332lg3.8125mmNN無滑脫阻力系數(shù)：其中：42)(ln01853. 0)(ln8725. 0ln18. 3

6、0523. 0lnyyyyS121HEy（ ）)2 . 12 . 2ln(yS當(dāng)1y1.2時，()()11111+1+=EEEEDvNglglmmRe，氣液兩相流動的雷諾數(shù)：匯報提綱匯報提綱Visual Basic 6.0簡介一、Visual Basic 6.0簡介l Visual Basic簡稱VB，是Microsoft公司1991年推出的Windows操作平臺下的可視化編程語言。l Visual Basic是在原DOS操作平臺下BASIC語言的基礎(chǔ)上發(fā)展起來的，它繼承和發(fā)展了BASIC語言的基本功能，具有簡單易學(xué)、功能強大等特點。l Visual Basic是面向?qū)ο蟮目梢暬幊坦ぞ?，具?/p>

7、快速的程序開發(fā)環(huán)境，程序界面簡單，代碼編寫量小的特點。l 一個完整的VB工程計算程序包括三部分：數(shù)據(jù)輸入、數(shù)據(jù)計算、數(shù)據(jù)輸出；不同程序間區(qū)別在于輸入數(shù)據(jù)形式、計算難易程度和輸出形式各有不同。匯報提綱匯報提綱Beggs Brill模型編寫基礎(chǔ)參數(shù)表液體密度，kg/m3氣體密度，kg/m3液體粘度，mPas氣體粘度，mPas105018.450.80.02液相體積流量，m3/s氣相體積流量，m3/s管子內(nèi)徑，m液體表面張力，N/m0.000580.010.0620.06重力加速度，m/s2管斜度，（） 平均壓力，MPa 9.81602.0 編寫B(tài)eggs-brill模型該程序有11個基本參數(shù)，這些基本參數(shù)可以由程序界面直接輸入Beggs Brill模型編寫計算壓力梯度的公式中調(diào)用了兩個函數(shù)，分別為“Dy_m函數(shù)和factor_f函數(shù)，其中Dy_m函數(shù)又調(diào)用了函數(shù)Hl_X和H l _ X 2 ； f a c t o r _ f 函 數(shù) 調(diào) 用 了 H l _ X 、 H l _ X 2 和 函 數(shù)factor_1(Nre_m)Dy_m函數(shù)factor_f函數(shù)Beggs Brill模型編寫函數(shù)Dy_m函數(shù)factor_fBeggs Brill模型編寫B(tài)eggs Brill模型編寫函數(shù)Hl_X1(x, E1, Nfr, Nvl)函數(shù)Hl_X2