单波长溶液液膜厚度测量方法

刘畅，田井辉，丁家琦，郭晓龙，石建伟，杨荟楠，蔡小舒(上海理工大学 能源与动力工程学院，上海市动力工程多相流动与传热重点实验室，上海，200093)



单波长溶液液膜厚度测量方法

刘畅，田井辉，丁家琦，郭晓龙，石建伟，杨荟楠，蔡小舒
(上海理工大学 能源与动力工程学院，上海市动力工程多相流动与传热重点实验室，上海，200093)

摘要：提出基于朗伯−比尔定律的未知质量分数尿素水溶液液膜厚度测量方法，选取尿素水溶液吸收率对质量分数的变化率为0的波数位置(6 613.25Cm−1)，研制单波长半导体激光吸收光谱(DLAS)尿素水溶液液膜厚度测量系统。通过溶液液膜厚度可调(0～1mm)的标准具对该方法测量精度进行验证。利用研制的单波长DLAS液膜测量系统对常温下水平透明石英玻璃板上的尿素水溶液液膜蒸发过程进行研究。研究结果表明：尿素水溶液液膜厚度测量平均相对误差约为1.50%，3次重复性实验厚度标准偏差小于24.68 μm。在液膜蒸发初期，液膜厚度随时间的变化率很大；在蒸发阶段中期，液膜蒸发速率趋于稳定；但在蒸发阶段后期，剧烈的光束转向现象造成了强烈的数据波动。

关键词：溶液液膜；测量；厚度；吸收光谱；蒸发

液膜在热能工程、空调制冷、航空航天、石油化工等工业中有着广泛的应用[1−5]，例如液膜流动(膜式流动)，因其具有流量小、温差小、传热传质系数高、热流密度高、气液接触面积大、表面更新快、结构简单且动力消耗小等优点在化学工业过程中应用广泛[6]，如膜式蒸发、膜式气液反应等。但在某些工业过程中，须避免液膜的形成。例如汽车尾气脱硝的选择性催化还原(selectiveCatalytic reduction，SCR)系统[7]中汽车尾气排放管上形成的尿素水溶液液膜[8]，其形成的三聚氰胺等混合物容易堵塞管路和影响催化反应效率，造成安全隐患[9]，因此，应尽可能控制尿素溶液液膜的形成。总之，对液膜厚度进行研究对优化所涉及的各种工业过程具有重要的意义。近年来，国内外科研人员通过数值模拟和无干扰的测量方法对溶液液膜进行了研究，如：BIRKHOLD等[9]采用数值模拟法研究了尿素水溶液液膜的厚度以及质量分数；邱庆刚等[10−11]通过数值模拟的方法分别对水平管外降膜流动的厚度和水平管降膜蒸发液膜厚度及传热系数进行了研究。溶液液膜厚度的无干扰测量方法主要包括电学法和光学法。电学法测量膜厚是基于不同厚度液膜的电学性质(例如电阻[12]、电容[13])不同，通过建立这些电学性质与液膜厚度的关系来获取液膜厚度信息，电学法也须对测量结果进行标定且测量精度及其测量范围易受到液膜导电性、液膜形态等因素的影响[14]。而光学法测量液膜厚度则主要具有无干扰、测量结果无需标定等优点，如：GRESZIK 等[15−16]利用激光诱导荧光技术和拉曼散射法获得了透明石英片上纯水水膜的二维分布信息，但测量结果须进行标定且设备昂贵；YANG等[17−19]分别利用吸收光谱法测量了液膜厚度。本文作者对不同质量分数(10%～50%)的尿素水溶液吸收光谱特性进行研究，选取吸收率对质量分数的变化率(dK/dc)为0的波数，建立基于朗伯−比尔定律的液膜厚度反演模型，研制单波长半导体激光吸收光谱(DLAS)尿素水溶液液膜厚度测量系统。利用液膜厚度可调的标准具对该系统进行验证，通过比较液膜厚度测量值与真实值获得该系统的测量误差。在此基础上，利用该系统对常温下水平透明石英玻璃板上的尿素水溶液液膜蒸发过程进行研究。

1　理论方法

尿素水溶液的光谱吸收率与溶液质量分数相关，图1所示为通过傅里叶变换红外光谱仪获得的质量分数为10%～50%的尿素水溶液在近红外区域(6 000～8000Cm−1)的吸收率K(v,c)[20]。由图1可见：随着尿素水溶液质量分数的增加，吸收峰向着波数减小的方向偏移。通过求取吸收率对质量分数的导数(dK/dc)可知：当波数为6 613.25Cm−1时，吸收率随质量分数的变化率为0(这里变化率近似为0)，即质量分数变化对吸收率无影响。因此，本文选取该波长，建立不同质量分数下的尿素水溶液液膜厚度反演模型。即基于朗伯−比尔定律，假设波长为v时，液膜厚度L可表示为

式中：It和 I0分别为透射光强和入射光强；K 为尿素水溶液在6 613.25Cm−1时的吸收率，该值为8.03Cm−1。

图1　质量分数为10%～50%的尿素水溶液的FT−IR谱Fig.1　FT−IR spectra of urea-water solutions for mass fractions of10%−50%

2 实验研究

2.1实验验证

本文利用液膜厚度可调(0～1mm)的标准具对单波长尿素水溶液液膜测量方法的精度进行验证。标准具内含有2块平行的石英玻璃板，通过调节石英玻璃板之间的距离可以改变溶液液膜的厚度。实验装置图如图2所示。利用激光控制器将分布反馈式半导体激光器的温度和电流控制在恒定值，以确保输出激光波数为 6613.25Cm−1。激光通过准直器聚焦后穿过标准具内的尿素水溶液液膜，之后被铟砷化镓探测器(Thorlabs,PDA10CS−EC)接收。实验数据的采集(采样频率为1kHz)和处理过程均在LabVIEW中进行。

图2实验装置示意图Fig.2Schematic drawing of experimental setup

分别对室温下质量分数为10%～50%的尿素水溶液液膜进行 3次重复测量，其测量结果如图 3所示。由图3可以看出：随着石英玻璃板间距的增加，不同质量分数的尿素水溶液液膜厚度的测量值与真实值的相对偏差均减小，尿素水溶液液膜厚度测量平均相对误差约为1.50%，3 次重复性实验厚度标准偏差小于24.68 μm。

图3　标准具验证实验结果Fig.3　Validation results measured inCalibration tool

2.2水平石英玻璃板上液膜蒸发研究

利用研制的单波长尿素水溶液液膜测量系统对常温下水平透明石英玻璃板上的初始质量浓度为20%的尿素水溶液液膜蒸发过程进行研究。在实验数据开始采集前，将透明石英玻璃板上无液膜时的激光光强记录为入射光强。在开始采集实验数据后，利用注射器喷射尿素水溶液，从而在透明石英玻璃板上形成液膜。图4所示为测得的液膜厚度随时间的变化趋势。由图4 可知：由于喷射过程中(数据采集后 0～11.88 s)剧烈的光束转向现象导致穿过液膜后的激光光束不能被接收光纤接收，从而造成了强烈的数据波动；在液膜蒸发初期(数据采集后11.88～606.82 s)，液膜厚度随时间的变化率很大，此时的液膜厚度减小可能是由液膜蒸发和液膜向外扩散共同引起；在蒸发阶段中期(数据采集后606.82～5 754.40 s)，液膜蒸发速率基本趋于稳定(0.104 μm/s)，测量厚度的浮动可能是不均匀蒸发引起液膜重组所致；在蒸发阶段后期(数据采集后5 754.40～6000.00 s)，液膜蒸发速率较大(约0.520 μm/s)。大约在数据采集后6 000.00 s，激光测量区域附近尿素已基本结晶析出，因此，剧烈的光束转向现象再次造成强烈的测量厚度波动，这种剧烈的光束转向现象可以在进一步工作中通过在探测器前放置积分球来消除。

图4水平透明石英板上尿素水溶液液厚度随时间的变化Fig.4Liquid film thickness variation of urea−water solutions on transparent quartz plate

3 结论

1)提出了基于朗伯−比尔定律的不同质量分数的尿素水溶液液膜厚度测量方法，选取尿素水溶液吸收率对质量分数的变化率为0的波数位置(6 613.25Cm−1)，研制了单波长半导体激光吸收光谱(DLAS)尿素水溶液液膜厚度测量系统。

2)利用液膜厚度可调(0～1mm)的标准具对该系统进行验证可知，随着标准具内液膜厚度的增加，不同质量分数的尿素水溶液液膜厚度的测量值与真实值的相对偏差均减小，测量平均相对误差约为1.50%。

3)剧烈的光束转向造成了液膜形成和蒸发后期强烈的数据波动，但是这种剧烈的光束转向现象可以在进一步工作中通过在探测器前放置积分球来消除。

4)不同温度下尿素水溶液在近红外区域的吸收光谱会出现偏移，下一步工作中将对本文提出新方法进行优化，以实现非常温下的液膜多参数研究。

参考文献：

[1]陈世昌,马建平,张先明,等.竖直降液膜流动在反应工程中的应用[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]齐复东,贾树本,马义伟.电站凝汽设备和冷却系统[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]戴干策.化学工程基础: 流体流动、 传热及传质[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]邱庆刚,吕多.水平管外降膜流动液膜厚度数值模拟[J].大连理工大学学报,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]王海涛,陈晔.水平管降膜蒸发液膜厚度及传热系数研究[J].海洋技术,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]崔洁,陈雪莉,王清立,等.电导法测量新型旋风分离器内液膜的分布规律[J].化工学报,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].计量技术,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)：国家自然科学基金资助项目(51306123)；博士点基金联合资助课题(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)

通信作者：杨荟楠，博士，副教授，从事激光光谱测量方法与应用的研究；E-mail: yanghuinan@usst.edu.cn