單波長溶液液膜厚度測量方法

劉暢，田井輝，丁家琦，郭曉龍，石建偉，楊薈楠，蔡小舒(上海理工大學(xué) 能源與動力工程學(xué)院，上海市動力工程多相流動與傳熱重點實驗室，上海，200093)

?

單波長溶液液膜厚度測量方法

劉暢，田井輝，丁家琦，郭曉龍，石建偉，楊薈楠，蔡小舒
(上海理工大學(xué) 能源與動力工程學(xué)院，上海市動力工程多相流動與傳熱重點實驗室，上海，200093)

摘要：提出基于朗伯?比爾定律的未知質(zhì)量分數(shù)尿素水溶液液膜厚度測量方法，選取尿素水溶液吸收率對質(zhì)量分數(shù)的變化率為0的波數(shù)位置(6 613.25Cm?1)，研制單波長半導(dǎo)體激光吸收光譜(DLAS)尿素水溶液液膜厚度測量系統(tǒng)。通過溶液液膜厚度可調(diào)(0～1mm)的標準具對該方法測量精度進行驗證。利用研制的單波長DLAS液膜測量系統(tǒng)對常溫下水平透明石英玻璃板上的尿素水溶液液膜蒸發(fā)過程進行研究。研究結(jié)果表明：尿素水溶液液膜厚度測量平均相對誤差約為1.50%，3次重復(fù)性實驗厚度標準偏差小于24.68 μm。在液膜蒸發(fā)初期，液膜厚度隨時間的變化率很大；在蒸發(fā)階段中期，液膜蒸發(fā)速率趨于穩(wěn)定；但在蒸發(fā)階段后期，劇烈的光束轉(zhuǎn)向現(xiàn)象造成了強烈的數(shù)據(jù)波動。

關(guān)鍵詞：溶液液膜；測量；厚度；吸收光譜；蒸發(fā)

液膜在熱能工程、空調(diào)制冷、航空航天、石油化工等工業(yè)中有著廣泛的應(yīng)用[1?5]，例如液膜流動(膜式流動)，因其具有流量小、溫差小、傳熱傳質(zhì)系數(shù)高、熱流密度高、氣液接觸面積大、表面更新快、結(jié)構(gòu)簡單且動力消耗小等優(yōu)點在化學(xué)工業(yè)過程中應(yīng)用廣泛[6]，如膜式蒸發(fā)、膜式氣液反應(yīng)等。但在某些工業(yè)過程中，須避免液膜的形成。例如汽車尾氣脫硝的選擇性催化還原(selectiveCatalytic reduction，SCR)系統(tǒng)[7]中汽車尾氣排放管上形成的尿素水溶液液膜[8]，其形成的三聚氰胺等混合物容易堵塞管路和影響催化反應(yīng)效率，造成安全隱患[9]，因此，應(yīng)盡可能控制尿素溶液液膜的形成。總之，對液膜厚度進行研究對優(yōu)化所涉及的各種工業(yè)過程具有重要的意義。近年來，國內(nèi)外科研人員通過數(shù)值模擬和無干擾的測量方法對溶液液膜進行了研究，如：BIRKHOLD等[9]采用數(shù)值模擬法研究了尿素水溶液液膜的厚度以及質(zhì)量分數(shù)；邱慶剛等[10?11]通過數(shù)值模擬的方法分別對水平管外降膜流動的厚度和水平管降膜蒸發(fā)液膜厚度及傳熱系數(shù)進行了研究。溶液液膜厚度的無干擾測量方法主要包括電學(xué)法和光學(xué)法。電學(xué)法測量膜厚是基于不同厚度液膜的電學(xué)性質(zhì)(例如電阻[12]、電容[13])不同，通過建立這些電學(xué)性質(zhì)與液膜厚度的關(guān)系來獲取液膜厚度信息，電學(xué)法也須對測量結(jié)果進行標定且測量精度及其測量范圍易受到液膜導(dǎo)電性、液膜形態(tài)等因素的影響[14]。而光學(xué)法測量液膜厚度則主要具有無干擾、測量結(jié)果無需標定等優(yōu)點，如：GRESZIK 等[15?16]利用激光誘導(dǎo)熒光技術(shù)和拉曼散射法獲得了透明石英片上純水水膜的二維分布信息，但測量結(jié)果須進行標定且設(shè)備昂貴；YANG等[17?19]分別利用吸收光譜法測量了液膜厚度。本文作者對不同質(zhì)量分數(shù)(10%～50%)的尿素水溶液吸收光譜特性進行研究，選取吸收率對質(zhì)量分數(shù)的變化率(dK/dc)為0的波數(shù)，建立基于朗伯?比爾定律的液膜厚度反演模型，研制單波長半導(dǎo)體激光吸收光譜(DLAS)尿素水溶液液膜厚度測量系統(tǒng)。利用液膜厚度可調(diào)的標準具對該系統(tǒng)進行驗證，通過比較液膜厚度測量值與真實值獲得該系統(tǒng)的測量誤差。在此基礎(chǔ)上，利用該系統(tǒng)對常溫下水平透明石英玻璃板上的尿素水溶液液膜蒸發(fā)過程進行研究。

1　理論方法

尿素水溶液的光譜吸收率與溶液質(zhì)量分數(shù)相關(guān)，圖1所示為通過傅里葉變換紅外光譜儀獲得的質(zhì)量分數(shù)為10%～50%的尿素水溶液在近紅外區(qū)域(6 000～8000Cm?1)的吸收率K(v,c)[20]。由圖1可見：隨著尿素水溶液質(zhì)量分數(shù)的增加，吸收峰向著波數(shù)減小的方向偏移。通過求取吸收率對質(zhì)量分數(shù)的導(dǎo)數(shù)(dK/dc)可知：當波數(shù)為6 613.25Cm?1時，吸收率隨質(zhì)量分數(shù)的變化率為0(這里變化率近似為0)，即質(zhì)量分數(shù)變化對吸收率無影響。因此，本文選取該波長，建立不同質(zhì)量分數(shù)下的尿素水溶液液膜厚度反演模型。即基于朗伯?比爾定律，假設(shè)波長為v時，液膜厚度L可表示為

式中：It和 I0分別為透射光強和入射光強；K 為尿素水溶液在6 613.25Cm?1時的吸收率，該值為8.03Cm?1。

圖1　質(zhì)量分數(shù)為10%～50%的尿素水溶液的FT?IR譜Fig.1　FT?IR spectra of urea-water solutions for mass fractions of10%?50%

2 實驗研究

2.1實驗驗證

本文利用液膜厚度可調(diào)(0～1mm)的標準具對單波長尿素水溶液液膜測量方法的精度進行驗證。標準具內(nèi)含有2塊平行的石英玻璃板，通過調(diào)節(jié)石英玻璃板之間的距離可以改變?nèi)芤阂耗さ暮穸取嶒炑b置圖如圖2所示。利用激光控制器將分布反饋式半導(dǎo)體激光器的溫度和電流控制在恒定值，以確保輸出激光波數(shù)為 6613.25Cm?1。激光通過準直器聚焦后穿過標準具內(nèi)的尿素水溶液液膜，之后被銦砷化鎵探測器(Thorlabs,PDA10CS?EC)接收。實驗數(shù)據(jù)的采集(采樣頻率為1kHz)和處理過程均在LabVIEW中進行。

圖2實驗裝置示意圖Fig.2Schematic drawing of experimental setup

分別對室溫下質(zhì)量分數(shù)為10%～50%的尿素水溶液液膜進行 3次重復(fù)測量，其測量結(jié)果如圖 3所示。由圖3可以看出：隨著石英玻璃板間距的增加，不同質(zhì)量分數(shù)的尿素水溶液液膜厚度的測量值與真實值的相對偏差均減小，尿素水溶液液膜厚度測量平均相對誤差約為1.50%，3 次重復(fù)性實驗厚度標準偏差小于24.68 μm。

圖3　標準具驗證實驗結(jié)果Fig.3　Validation results measured inCalibration tool

2.2水平石英玻璃板上液膜蒸發(fā)研究

利用研制的單波長尿素水溶液液膜測量系統(tǒng)對常溫下水平透明石英玻璃板上的初始質(zhì)量濃度為20%的尿素水溶液液膜蒸發(fā)過程進行研究。在實驗數(shù)據(jù)開始采集前，將透明石英玻璃板上無液膜時的激光光強記錄為入射光強。在開始采集實驗數(shù)據(jù)后，利用注射器噴射尿素水溶液，從而在透明石英玻璃板上形成液膜。圖4所示為測得的液膜厚度隨時間的變化趨勢。由圖4 可知：由于噴射過程中(數(shù)據(jù)采集后 0～11.88 s)劇烈的光束轉(zhuǎn)向現(xiàn)象導(dǎo)致穿過液膜后的激光光束不能被接收光纖接收，從而造成了強烈的數(shù)據(jù)波動；在液膜蒸發(fā)初期(數(shù)據(jù)采集后11.88～606.82 s)，液膜厚度隨時間的變化率很大，此時的液膜厚度減小可能是由液膜蒸發(fā)和液膜向外擴散共同引起；在蒸發(fā)階段中期(數(shù)據(jù)采集后606.82～5 754.40 s)，液膜蒸發(fā)速率基本趨于穩(wěn)定(0.104 μm/s)，測量厚度的浮動可能是不均勻蒸發(fā)引起液膜重組所致；在蒸發(fā)階段后期(數(shù)據(jù)采集后5 754.40～6000.00 s)，液膜蒸發(fā)速率較大(約0.520 μm/s)。大約在數(shù)據(jù)采集后6 000.00 s，激光測量區(qū)域附近尿素已基本結(jié)晶析出，因此，劇烈的光束轉(zhuǎn)向現(xiàn)象再次造成強烈的測量厚度波動，這種劇烈的光束轉(zhuǎn)向現(xiàn)象可以在進一步工作中通過在探測器前放置積分球來消除。

圖4水平透明石英板上尿素水溶液液厚度隨時間的變化Fig.4Liquid film thickness variation of urea?water solutions on transparent quartz plate

3 結(jié)論

1)提出了基于朗伯?比爾定律的不同質(zhì)量分數(shù)的尿素水溶液液膜厚度測量方法，選取尿素水溶液吸收率對質(zhì)量分數(shù)的變化率為0的波數(shù)位置(6 613.25Cm?1)，研制了單波長半導(dǎo)體激光吸收光譜(DLAS)尿素水溶液液膜厚度測量系統(tǒng)。

2)利用液膜厚度可調(diào)(0～1mm)的標準具對該系統(tǒng)進行驗證可知，隨著標準具內(nèi)液膜厚度的增加，不同質(zhì)量分數(shù)的尿素水溶液液膜厚度的測量值與真實值的相對偏差均減小，測量平均相對誤差約為1.50%。

3)劇烈的光束轉(zhuǎn)向造成了液膜形成和蒸發(fā)后期強烈的數(shù)據(jù)波動，但是這種劇烈的光束轉(zhuǎn)向現(xiàn)象可以在進一步工作中通過在探測器前放置積分球來消除。

4)不同溫度下尿素水溶液在近紅外區(qū)域的吸收光譜會出現(xiàn)偏移，下一步工作中將對本文提出新方法進行優(yōu)化，以實現(xiàn)非常溫下的液膜多參數(shù)研究。

參考文獻：

[1]陳世昌,馬建平,張先明,等.豎直降液膜流動在反應(yīng)工程中的應(yīng)用[J].化工進展,2014,33(10): 2528?2534.CHEN Shichang,MA Jianping,ZHANG Xianming,et al.Advances in application of falling liquid film for reaction engineering[J].Chemical Industry and Engineering Progress,2014,33(10): 2528?2534.

[2]SUN Y Z,SONG X F,JIN M M,et al.Gas-liquid reactiveCrystallization of lithiumCarbonate by a falling filmColumn[J].Industrial EngineeringChemistry Research,2013,52(49):17598?17606.

[3]PARKC D,NOSOKO T.Three-dimensional wave dynamics on a falling film and associated mass transfer[J].AIChE Journal,2003,49(11): 2715?2727.

[4]YEONG K K,GAVRIILIDIS A,ZAPF R,et al.Catalyst preparation and deactivation issues for nitrobenzene hydrogenation in a microstructure falling film reactor[J].Catalysis Today,2003,81(4): 641?651.

[5]齊復(fù)東,賈樹本,馬義偉.電站凝汽設(shè)備和冷卻系統(tǒng)[M].北京: 水利電力出版社,1990: 28?32.QI Fudong,JIA Shuben,MA Yiwei.Power plantCondenser device andCooling system[M].Beijing: Water Resources and Electric Power Press,1990: 28?32.

[6]戴干策.化學(xué)工程基礎(chǔ): 流體流動、 傳熱及傳質(zhì)[M].北京: 中國石化出版社,1991: 25?32.DAI Gance.Chemical engineering foundation: fluid flow,heat transfer and mass transfer[M].Beijing:China Petrochemical Press,1991: 25?32.

[7]SHIAEDH S.Development of multi-component evaporation models and 3D modeling of NOx-SCR reduction system[D].Toulouse: University of Toulouse.Institute National Polytechnique de Toulouse,2011: 83?84.

[8]GIESHOFF J,PFEIFER M,SCHAEFER-SINDLINGER A,et al.Advanced urea SCRCatalysts for automotive applications[C]//SAE Technical Papers.Detroit,United States: SAE WorldCongress,2001: 0514-1?0514-10.

[9]BIRKHOLD F,MEINGAST U,WASSERMANN P,et al.Analysis of the injection of urea-water-solution for automotive SCR DeNOx-systems: modeling of two-phase flow and spray/wall-interaction[C]//SAE Technical Papers.Detroit,United States: SAE WorldCongress,2006: 0643-1?0643-13.

[10]邱慶剛,呂多.水平管外降膜流動液膜厚度數(shù)值模擬[J].大連理工大學(xué)學(xué)報,2013,53(3): 359?363.QIU Qinggang,Lü Duo.Numerical simulation of falling film thickness around a horizontal tube[J].Journal of Dalian University of Technology,2013,53(3): 359?363.

[11]王海濤,陳曄.水平管降膜蒸發(fā)液膜厚度及傳熱系數(shù)研究[J].海洋技術(shù),2013,32(4): 81?84.WANG Haitao,CHEN Ye.Study on film thickness and evaporation heat transferCoefficient of horizontal tube falling film evaporator[J].Ocean Technology,2013,32(4): 81?84.

[12]崔潔,陳雪莉,王清立,等.電導(dǎo)法測量新型旋風(fēng)分離器內(nèi)液膜的分布規(guī)律[J].化工學(xué)報,2009,60(6):1487?1493.CUI Jie,CHEN Xueli,WANG Qingli,et al.Double-parallelConductance probe for measuring thickness of liquid film in new-typeCyclone separator[J].Journal ofChemical Industry and Engineering,2009,60(6):1487?1493.

[13]張鋒,耿皎,馬少玲,等.利用電容法測量降膜厚度[J].計量技術(shù),2005,5(6): 34?35.ZHANG Feng,GENG Jiao,MA Shaoling,et al.Measurement of the falling film thickness byCapacitance method[J].Measurement Technique,2005,5(6): 34?35.

[14]WANG R,LEE B A,LEE J S,et al.Analytical estimation of liquid film thickness in two-phase annular flow using electrical resistance measurement[J].Applied Mathematical Modelling,2012,36(7): 2833?2840.

[15]GRESZIK D,YANG H,DREIER T,et al.Laser-based diagnostics for the measurement of liquid water film thickness[J].Applied Optics,2011,50(4): A60?A67.

[16]GRESZIK D,YANG H,DREIER T,et al.Measurements of water film thickness by laser-induced fluorescence and Raman imaging[J].Applied Physics B: Laser and Optics,2011,102(1):123?132.

[17]YANG H,GRESZIK D,WLOKAS I,et al.Tunable diode laser absorption sensor for the simultaneous measurement of water film thickness,liquid- and vapor-phase temperature[J].Applied Physics B: Laser and Optics,2011,104(7): 21?27.

[18]YANG H,GRESZIK D,DREIER T,et al.Simultaneous measurement of liquid water film thickness and vapor temperature using near-infrared tunable diode laser spectroscopy[J].Applied Physics B: Laser and Optics,2010,99(3): 385?390.

[19]SALAZER R P,MARSHALL E.Time average local thickness measurement in falling liquid film flow[J].International Journal of Multiphase Flow,1978,4(4): 408?412.

[20]YANG Huinan,GUO Xiaolong,SU Mingxun,et al.Novel method for simultaneous measurement of film thickness and mass fraction of urea-water solution[J].Chinese Optics Letters,2014,12(12):123102-1?123102-4.

(編輯 劉錦偉)

Film thickness measurement of aqueous solution based on single diode laser absorption spectroscopy

LIUChang,TIAN Jinghui,DING Jiaqi,GUO Xiaolong,SHI Jianwei,YANG Huinan,CAI Xiaoshu
(Shanghai Key Laboratory of Multiphase Flow and Heat Transfer in Power Engineering,School of Energy and Power Engineering,University of Shanghai for Science and Technology,Shanghai 200093,China)

Abstract:A novel method was proposed to determine liquid film thickness of urea-water solutions at different mass fractions.The wavenumber position at6613.25Cm?1wasChosen,where the derivative of the absorptionCoefficient of urea-water solutions with respect to massConcentration was zero.A single diode laser absorption spectroscopy(DLAS)sensor was developed to obtain liquid film thickness of urea-water solutions.The measurement accuracy of this method was validated with aCalibration tool with which liquid film thicknessCan be adjusted from0to1mm.Furthermore,the evaporation process of liquid film of urea-water solutions on transparent quartz plate at room temperature was also studied.The results show that average relative error between the measured thickness and known liquid film thickness is about1.50%，and the standard deviation of film thickness is less than 24.68 μm for three repeated experiments.The evaporation rate of film thickness is large at the beginning of film evaporation process,and it keepsConstant in the middle of the film evaporation process.However,at the end of the evaporation process,strong fluctuation of the film thickness occurs because of the serious beam steering effect.

Key words:liquid film of aqueous solution? measurement? thickness? absorption spectroscopy? evaporation

中圖分類號：TK31

文獻標志碼：A

文章編號：1672?7207(2016)01?0286?04

DOI:10.11817/j.issn.1672-7207.2016.01.039

收稿日期：2015?01?31；修回日期：2015?03?20

基金項目(Foundation item)：國家自然科學(xué)基金資助項目(51306123)；博士點基金聯(lián)合資助課題(20133120120008)(Project(51306123)supported by the National Natural Science Foundation ofChina? Project(20133120120008)supported by the Joint Specialized Research Found for the Doctoral Program of Higher Education)

通信作者：楊薈楠，博士，副教授，從事激光光譜測量方法與應(yīng)用的研究；E-mail: yanghuinan@usst.edu.cn