馬 力,仇性啟,王 健,鄭志偉,閻紅巧,易志勇
(1.中國石油大學(xué)(華東)化學(xué)工程學(xué)院,山東 青島 266580;2.北京石油化工工程有限公司,北京 100107)
噴霧冷卻技術(shù)廣泛應(yīng)用于石油、化工等工業(yè)領(lǐng)域[1]。在生產(chǎn)乙烯的裂解裝置中,急冷器的作用是將裂解氣從500 ℃快速冷卻到200 ℃,急冷油液滴應(yīng)該在最短的時(shí)間內(nèi)吸收熱量后完全汽化,否則未蒸發(fā)的急冷油液滴會(huì)黏附在急冷器壁面造成結(jié)焦[2]。急冷油冷卻裂解氣的實(shí)質(zhì)是液滴在高溫氣流中的蒸發(fā)問題,研究液滴在高溫氣流中的蒸發(fā)機(jī)理,將最大程度減少急冷器的結(jié)焦,提高裂解裝置的生產(chǎn)效率,并實(shí)現(xiàn)節(jié)能降耗的目的。
目前,國內(nèi)外學(xué)者對(duì)液滴蒸發(fā)進(jìn)行了大量理論研究[3-6],建立了液滴蒸發(fā)模型,分析了液滴運(yùn)動(dòng)規(guī)律,總結(jié)出液滴蒸發(fā)過程中溫度和壓力等因素對(duì)液滴溫度、直徑的影響規(guī)律。但由于實(shí)驗(yàn)條件的限制,只有少數(shù)學(xué)者對(duì)液滴蒸發(fā)進(jìn)行了實(shí)驗(yàn)研究,如蘇凌宇等[7]對(duì)負(fù)壓靜止?fàn)顟B(tài)液滴蒸發(fā)進(jìn)行了實(shí)驗(yàn)研究,但對(duì)高溫氣流中液滴蒸發(fā)的實(shí)驗(yàn)研究鮮有報(bào)道。
本工作建立了液滴蒸發(fā)實(shí)驗(yàn)裝置,以高溫氣流中的液滴為研究對(duì)象,采用懸掛液滴法研究液滴蒸發(fā)特性,探討液滴蒸發(fā)的影響因素。
液滴蒸發(fā)的實(shí)驗(yàn)裝置[8]見圖1。實(shí)驗(yàn)裝置由高速攝像儀、熱風(fēng)機(jī)、熱線風(fēng)速儀和熱電偶等主要設(shè)備組成。高速攝像儀為德國Basler公司A504KC型數(shù)字面陣攝像機(jī),拍攝的圖像以TIFF圖像格式通過圖像采集板存儲(chǔ)存在計(jì)算機(jī)中,實(shí)驗(yàn)完畢后可方便地應(yīng)用軟件進(jìn)行后續(xù)的圖像處理以獲得液滴直徑變化的信息。熱風(fēng)機(jī)能夠提供最高溫度350 ℃、最大體積流量5.5 m3/min的熱空氣;熱線風(fēng)速儀測量范圍0~50.0 m/s,精度0.01 m/s,能夠準(zhǔn)確測量氣流速度;熱電偶測量范圍0~1 000 ℃,測量精度0.1 ℃。
圖1 液滴蒸發(fā)的實(shí)驗(yàn)裝置Fig.1 Experimental device for droplet evaporation.
采用懸掛液滴法研究液滴的蒸發(fā)特性。操作過程為:由熱風(fēng)機(jī)輸送的高溫氣流,經(jīng)過直管段和擾流柵在出口位置獲得穩(wěn)定氣流;液滴在管外生成并懸掛在石英絲底部,而后轉(zhuǎn)移到高溫氣流中進(jìn)行蒸發(fā);分別利用熱電偶和熱線風(fēng)速儀測量氣流的溫度和流速,利用高速攝像儀拍照蒸發(fā)過程,拍照速率10 幀/s,每幀圖像的大小為1 280×1 024,通過后續(xù)的圖像處理軟件獲得液滴蒸發(fā)過程中的參數(shù)變化。通過調(diào)節(jié)空氣加熱器、風(fēng)機(jī)和流量調(diào)節(jié)閥可以任意設(shè)置氣流流速和溫度,以改變實(shí)驗(yàn)工況。
數(shù)據(jù)處理時(shí)選用等體積球的直徑表征液滴的蒸發(fā)過程。圖像的處理方法如圖2所示,首先用游標(biāo)卡尺測量懸掛液滴的石英絲的實(shí)際直徑,并用圖2中的圖像處理軟件測量石英絲直徑在圖像中的大小,計(jì)算圖像中的尺寸與實(shí)際尺寸的換算公式,利用換算公式計(jì)算不同時(shí)刻圖像中液滴直徑對(duì)應(yīng)的液滴實(shí)際尺寸。根據(jù)獲得的不同時(shí)刻的液滴直徑,總結(jié)液滴直徑、質(zhì)量隨時(shí)間的變化規(guī)律。
圖2 圖像處理方法Fig.2 Image processing method.
蒸發(fā)過程從液滴放入測試位置開始,到液滴全部蒸發(fā)為止,此時(shí)液滴的直徑終值為石英絲小球的直徑。高速攝像儀記錄下的液滴蒸發(fā)過程如圖3所示,在蒸發(fā)過程中液滴由球體逐漸變成紡錘體。這是因?yàn)楫?dāng)液滴進(jìn)入高溫氣流中后,除本身重力外,還受到浮力、曳力和表面張力等。初始時(shí)刻,液滴各項(xiàng)受力平衡,液滴基本呈球體;隨蒸發(fā)過程的進(jìn)行,液滴質(zhì)量減小,重力不能平衡氣流的浮力,液滴呈橢圓紡錘體直至完全蒸發(fā)。
圖3 液滴蒸發(fā)過程Fig.3 Droplet evaporation process.
液滴蒸發(fā)的曲線見圖4。由圖4可見,蒸發(fā)初始,液滴直徑(D)及直徑的平方(D2)減小的幅度較小。這是因?yàn)椋?)相比高溫氣流的溫度,液滴的初始溫度較低,液滴從高溫氣流中吸收熱量主要用于液滴升溫,蒸發(fā)量較少;2)液滴進(jìn)入高溫氣流后,液滴溫度驟然升高,液滴體積發(fā)生膨脹,當(dāng)因溫度升高使液滴直徑增大的速率大于因液滴蒸發(fā)而使直徑減小的速率時(shí),液滴直徑反而會(huì)略有增大。在蒸發(fā)較慢的實(shí)驗(yàn)工況下,液滴直徑增大這種現(xiàn)象更為明顯。液滴直徑減小幅度較小的階段為初始加熱階段,為非穩(wěn)態(tài)階段,持續(xù)時(shí)間較短。由圖4可還知,在初始加熱階段之后,液滴進(jìn)入穩(wěn)定蒸發(fā)階段,D2~t曲線基本為一直線,即dD2/dt=,為液滴平均蒸發(fā)常數(shù)。液滴的蒸發(fā)過程基本滿足D2定律。
圖4 液滴蒸發(fā)的曲線Fig.4 Droplet evaporation curves.
液滴蒸發(fā)速率(K)為單位時(shí)間內(nèi)液滴質(zhì)量的蒸發(fā)量,定義如下:
單位表面積的蒸發(fā)速率(K0)為單位時(shí)間單位表面積上液滴的蒸發(fā)量,定義如下:
液滴蒸發(fā)速率隨時(shí)間的變化見圖5。由圖5可看出,液滴蒸發(fā)初期,K和K0先減小后增大,經(jīng)歷一個(gè)短暫的非穩(wěn)態(tài)階段。在液滴置于高溫氣體環(huán)境的初始時(shí)刻,液滴表面與氣流環(huán)境間具有最大的溫度梯度和濃度梯度,即兩者間具有最大的傳熱傳質(zhì)勢差,使蒸發(fā)開始瞬間,K和K0均較大。經(jīng)歷初始時(shí)刻后K和K0迅速減小,這是因?yàn)椋?)蒸氣擴(kuò)散速度小于液滴蒸發(fā)速度,使得蒸氣在液滴周圍大量聚集,液滴與氣流環(huán)境的濃度梯度迅速減?。?)液滴在進(jìn)入高溫環(huán)境后迅速升溫,液滴與氣流溫度梯度下降;3)溫度升高液滴膨脹,液滴直徑變化更小,使K和K0迅速減小。而且液滴初始直徑越小,氣流溫度越高,熱擴(kuò)散系數(shù)越大,液滴與氣流溫度梯度和濃度梯度變化就越快,K和K0減小得越快。
圖5 液滴蒸發(fā)速率隨時(shí)間的變化Fig.5 Variation of droplet evaporation rate with t.
由圖5還可見,K和K0減小后迅速增大并趨于 穩(wěn)定。這是因?yàn)殡S著液滴蒸發(fā)的進(jìn)行,液滴與氣流溫度梯度和濃度梯度達(dá)到相對(duì)平衡狀態(tài),液滴蒸發(fā)進(jìn)入一個(gè)相對(duì)穩(wěn)定的階段,液滴直徑和液滴蒸發(fā)表面積穩(wěn)定減小,K0基本保持不變,K平穩(wěn)下降。在蒸發(fā)末期,K持續(xù)下降,但K0急劇增大。這是因?yàn)樵谡舭l(fā)末期,液滴周圍的蒸氣濃度下降,濃度梯度增大,所以蒸發(fā)推動(dòng)力增大,但此時(shí)液滴尺寸非常小,液滴表面積小,因此液滴總體蒸發(fā)速度依然下降,K很小,而K0增大。
初始加熱時(shí)間(ti)與蒸發(fā)生存時(shí)間(te)的定義見圖6。在(D/D0)2≤0.5的區(qū)域,由液滴直徑與初始直徑的比值的平方擬合出斜直線;將該斜直線延伸,與橫坐標(biāo)的交叉點(diǎn)即為te,與(D/D0)2=1.0橫線的交叉點(diǎn)即為ti。用Φi=ti/te度量液滴蒸發(fā)不穩(wěn)定度,用于表征氣流溫度和氣流速度對(duì)液滴初始加熱時(shí)間和液滴生存時(shí)間的關(guān)系。
圖6 初始加熱時(shí)間與蒸發(fā)存在時(shí)間的定義Fig.6 Definition of initial heat-up time(ti) and corrected evaporation lifetime(te).
氣流溫度和氣流速度對(duì)Φi的影響見圖7。由圖7可見,在常壓的蒸發(fā)過程中,Φi值在0~0.2之間。與文獻(xiàn)[9]報(bào)道的結(jié)果相比較,Φi值上限基本一致,Φi值下限較小,這是由實(shí)驗(yàn)方法不同造成的,文獻(xiàn)值進(jìn)行測試時(shí),液滴在蒸發(fā)實(shí)驗(yàn)開始前已經(jīng)放置于實(shí)驗(yàn)環(huán)境中,而本工作的測試中液滴放入實(shí)驗(yàn)測試位置需要一定時(shí)間,且過程中氣流與液滴間相對(duì)速度增大,液滴在拍照之前已經(jīng)開始預(yù)熱,使高速攝像儀記錄的液滴初始加熱時(shí)間縮短,導(dǎo)致Φi值下限比文獻(xiàn)值小。
由圖7還可知,當(dāng)氣流溫度小于475 K時(shí),Φi值隨氣流溫度的升高而增大,氣流溫度對(duì)液滴整個(gè)蒸發(fā)過程的影響大于對(duì)液滴初始加熱過程的影響;當(dāng)氣流溫度大于等于500 K時(shí),Φi值隨氣流溫度的升高而減小,氣流溫度對(duì)液滴整個(gè)蒸發(fā)過程的影響小于對(duì)液滴初始加熱過程的影響。所以在475~500 K內(nèi),存在一個(gè)臨界氣流溫度使Φi達(dá)到最大值,即初始加熱時(shí)間占液滴蒸發(fā)總時(shí)間的比例最大。此外,在氣流溫度低于475 K時(shí),氣流速度對(duì)Φi的影響不是很明顯;當(dāng)氣流溫度大于等于500 K時(shí),氣流速度越大,Φi值越小,氣流速度對(duì)液滴整個(gè)蒸發(fā)過程的作用更強(qiáng)。這說明氣流速度對(duì)液滴的初始加熱過程影響不大,但對(duì)液滴整個(gè)蒸發(fā)過程影響明顯,氣流速度越大,影響越明顯。
圖7 氣流溫度和氣流速度對(duì)Φi的影響Fig.7 Effects of gas flow temperature and velocity on the ratio of ti to te (Φi).
氣流溫度和氣流速度對(duì)液滴平均蒸發(fā)常數(shù)的影響見圖8。
圖8 氣流溫度和氣流速度對(duì)液滴平均蒸發(fā)常數(shù)的影響Fig.8 Effects of gas flow temperature and velocity on the average evaporation constant of droplets().
由圖8可知,隨氣流溫度升高和氣流速度增加,液滴平均蒸發(fā)常數(shù)增大。這是因?yàn)橐旱闻c氣流間的溫差越大,氣流向液滴的傳熱能力越強(qiáng),液滴蒸發(fā)加快;氣流速度越大,對(duì)流作用增強(qiáng),對(duì)流傳熱系數(shù)、質(zhì)量擴(kuò)散系數(shù)均增大,強(qiáng)化了傳熱傳質(zhì),因此液滴平均蒸發(fā)常數(shù)越大。液滴內(nèi)部的物質(zhì)分子必須到達(dá)液滴表面,暴露于氣流環(huán)境中才能蒸發(fā),其途徑有:1)液滴表面被動(dòng)的收縮;2)液滴內(nèi)部環(huán)流將內(nèi)部物質(zhì)分子輸運(yùn)到液滴表面發(fā)生蒸發(fā)。在有自然對(duì)流或強(qiáng)迫對(duì)流或兩者都存在時(shí),液滴和周圍氣流間產(chǎn)生非徑向的相對(duì)速度,氣流對(duì)表面的剪切應(yīng)力引起液滴內(nèi)部的環(huán)流。氣流速度越大,對(duì)流強(qiáng)度增大,液滴內(nèi)部的環(huán)流變得更加激烈,從而加速了液相質(zhì)量輸運(yùn)過程,縮短了液滴內(nèi)部物質(zhì)分子到達(dá)液滴表面的時(shí)間,加快了液滴蒸發(fā)??傊?,氣流溫度越高和氣流速度越大,液滴蒸發(fā)越快,液滴生存時(shí)間越短。
1)液滴蒸發(fā)過程經(jīng)歷非穩(wěn)態(tài)的初始加熱階段和穩(wěn)定的階段。在初始加熱階段,液滴的K和K0先減小再增大;在相對(duì)穩(wěn)定的階段,液滴的K平穩(wěn)減小,K0保持平穩(wěn)后增大。液滴蒸發(fā)過程基本遵循D2定律。
2)氣流溫度和氣流速度影響液滴蒸發(fā)的不穩(wěn)定度,在475~500 K內(nèi)存在一個(gè)臨界溫度使Φi值達(dá)到最大值,液滴初始加熱時(shí)間在整個(gè)蒸發(fā)時(shí)間中所占比例最大。
3)氣流溫度越高和氣流速度越大,液滴與高溫氣流間的溫差越大,對(duì)流作用增強(qiáng),液滴蒸發(fā)越快,K和K0越大。
符 號(hào) 說 明
A 液滴蒸發(fā)表面積,m2
D 液滴直徑,mm
D0液滴初始直徑,mm
K 液滴質(zhì)量蒸發(fā)速率,kg/s
K0單位表面積的蒸發(fā)速率,kg/(s·m2)
m 前一時(shí)刻液滴質(zhì)量,kg
m′ 下一時(shí)刻液滴質(zhì)量,kg
t 液滴蒸發(fā)時(shí)間,s
ti蒸發(fā)初始加熱時(shí)間,s
te蒸發(fā)生存時(shí)間,s
Δt 時(shí)間間隔,s
Φi液滴初始加熱時(shí)間與液滴生存時(shí)間的比值
[1]徐以撒,陳懋湘,趙崇敏.順丁烯二酸噴霧法脫水工藝[J].石油化工,1989,18(10):227-230.
[2]李和平,戚俊清,朱廷哲,等.乙烯生產(chǎn)裝置中急冷鍋爐結(jié)焦過程的數(shù)學(xué)模型與工業(yè)模擬[J].石油化工,2000,29(8):600-612.
[3]Spalding D B.The Combustion of Liquid Fuels[C]//Proceedings of the Fourth International Symposium in Combustion.Pittsburgh:The Combustion Institute,1953:847-864.
[4]Consolini L,Aggarwal S K,Murad S.A Molecular Dynamics Simulation of Droplet Evaporation[J].Int J Heat Mass Transfer,2003,46:3179-3188.
[5]丁繼賢,孫鳳賢,姜任秋.對(duì)流條件下環(huán)境壓力對(duì)液滴蒸發(fā)的影響研究[J].哈爾濱工程大學(xué)學(xué)報(bào),2007,28(10):1104-1108.
[6]王健,仇性啟,申玉輝.急冷油內(nèi)液滴運(yùn)動(dòng)蒸發(fā)的數(shù)值模擬[J].石油化工,2011,40(10):1074-1078.
[7]蘇凌宇,于江飛,劉衛(wèi)東.乙醇液滴蒸發(fā)過程對(duì)壓力動(dòng)態(tài)響應(yīng)的試驗(yàn)研究[J].推進(jìn)技術(shù),2009,30(3):286-291.
[8]馬力,仇性啟,王健.液滴蒸發(fā)實(shí)驗(yàn)裝置相似性設(shè)計(jì)[J].工業(yè)加熱,2012,41(4):36-39.
[9]Nomura H,Ujiie Y,Rath H J,et al.Experimental Study on High-Pressure Droplet Evaporation Using Microgravity Conditions[J].Symp(Int) Combust,1996,26(1):1267-1273.