聪明文档网

聪明文档网

最新最全的文档下载
当前位置: 首页> 安全验证

安全验证

时间:2023-11-20 23:28:20    下载该word文档
界面传质中Rayleigh对流的定量分析
陈曼;赵嵩;曾爱武;于海路
【摘要】SchlierensystemwasbuilttovisualizethephenomenonofRayleighconvectionperpendiculartothegas-liquidinterfaceduringtheabsorptionofcarbondioxideintopureethanolinthespecialmasstransferdevice.Astheabsorptionprogressed,thestabilityofliquidlayerweakenedanddisturbancesaggravatedtheinstability,whichledtoturbulenceintheinterfacialvicinityandwentontodevelopplumeconvectiondowntothebulkliquidwiththegrowthandintegrationofconvectivecells.Quantitativemethodwasappliedtoobtaintheconcentrationcontoursoftheliquidphase.ConcentrationdistributionandvariationofinstantaneousmasstransfercoefficientrepresentedtheonsetanddevelopmentofRayleighconvectionanditsreinforcementeffectonthemasstransferprocess.InterfacialconcentrationdistributionandcriticalRayleighnumberexplainedthemechanisminwhichthenonuniformmasstransfercontributedtothetriggerofturbulence.Theexperimentalresultsshowedthatthedownwardhighconcentrationflowspromotedtheexchangeoftheliquidbetweeninterfacialvicinityandliquidbulk,showingthatRayleighconvectioncanenhancethemasstransferprocess.%通过纹影光路观察了特定气液传质装置中乙醇吸收CO2过程所引发的Rayleigh对流在垂直界面方向上的发展过程。随着溶质吸收的进行,液层的流体稳定性变弱,扰动加剧气液界面失稳并发生湍动,进而发展为羽状流并逐步向液相主体发展,在此过程中伴随着对流胞的融合与增长。液层的浓度分布可通过对相应液层纹影图像进行定量分
析获得。液层浓度分布和瞬时传质系数变化表征了Rayleigh对流的引发与发展及其对传质过程的强化效果,界面浓度分布及临界Rayleigh数解释了非均匀传质对湍动的引发机理。羽状流将高浓度液体快速带入主体,加速了近界面液层与主体液层的混合,增强了气液传质。【期刊名称】《化工学报》【年(,期】2016(067011【总页数】8(P4566-4573
【关键词】Rayleigh对流;纹影;吸收;界面;传质;定量分析【作者】陈曼;赵嵩;曾爱武;于海路
【作者单位】天津大学化工学院,化学工程联合国家重点实验室,天津300072;天津大学化工学院,化学工程联合国家重点实验室,天津300072;天津大学化工学院,化学工程联合国家重点实验室,天津300072;天津大学化工学院,化学工程联合国家重点实验室,天津300072【正文语种】【中图分类】TQ028
Rayleigh对流是相际传热与传质过程中因相界面处流体密度变化所产生的界面湍动现象[1-2]。由于温度或浓度差异导致的近界面处流体密度差异将引发流体的不稳定性,密度梯度产生的推动力将驱动界面处高密度流体向下运动到主体,增强了相际传递过程。
纹影法[3]、投影法[4]以及干涉法[5]作为非介入式方法被广泛应用于Rayleigh流的实验研究中。在现象观察实验中,研究人员主要对垂直于界面方向的对流发展
进行了研究。Okhotsimskii[6]运用纹影技术在多种有机溶剂吸收CO2的过程中观测到向液相主体发展的羽状对流。Arendt[7]在对水吸收CO2过程的纹影观测中观察到类似现象,并计算了相应过程的传质速率,结果表明自然对流的产生增强了气液传质。利用放大投影法,沙勇等[4]在垂直于界面方向上观察到了丙酮-水溶液解吸传质过程中的羽状对流结构;此外,王勇等[8]Sun[9]在平行于界面方向上观察了有机溶剂在水平层流运动过程中吸收CO2的传质过程,得到了滚筒状和多边形的Rayleigh对流结构。以上研究主要是对自然对流现象的定性研究。为获得自然对流的更多信息,数值模拟(FEM[10]FVM[11-12]LBM[13-15]获取界面对流相应的浓度场与速度场已有广泛应用。在实验方面,Guo[5]运用激光干涉技术获取了气液传质过程中的近界面浓度分布;Liu[16]与于海路等[17]利用定量纹影分析了气液传质中Marangoni对流,获取了界面浓度与表面张力分布。
本实验以乙醇吸收CO2的过程作为研究对象,单独研究Rayleigh对流。在获得界面对流纹影图像的同时对其进行定量分析,得到相应图片的液层浓度分布,并以此分析了Rayleigh对流的引发与发展过程,更为深入地研究其对于传质的增强机理。
1.1实验试剂
实验用的分析纯(纯度≥99.7%)乙醇由天津康科德有限公司提供;N2CO2(纯度≥99.99%)由天津市北方气体公司提供。1.2实验装置与操作
本文所用纹影系统(图1)及实验装置(图2)与于海路等[17]所用相同,实验流程略有变化。调节纹影光路使平行光路垂直通过传质装置中的竖直液层,得到界面对流的主视图。为获取较清晰的图片,CCD的帧数调整至1/s,图片像素为1236×1628,其中的像素格间距所对应的实际尺寸为0.0384mm

实验前,将纯乙醇进行超声除气的预处理,用注射泵将其注入到装置内部,液层深度约为60mm,厚度为2mm。为避免传质过程中液相乙醇的挥发,实验所用的高纯CO2在通过装有乙醇的鼓泡装置饱和后,由装置上端的进气口进入。在此选取气体流量为100ml·min-1,以减少吹扫气对气液界面造成的扰动。气液传质实验在常温25℃、常压101.3kPa下进行,并利用PID温度控制仪保证吹扫气与液层温度相同,以消除因温度变化引发的自然对流。
本实验中,由于传质装置的空间限制,腔室内的液层宽度均一,平行光线通过时,液层中折射率梯度的存在是光线发生偏折的唯一因素。液相的折射率随CO2的浓度升高而降低[17],向下偏折的光线被刀口遮挡使得相应的区域变暗,纹影图片中的亮暗分布直接反映了液层各处的折射率梯度分布。
当平行光线垂直通过较窄液层时,因折射率存在导致的偏折角度[3]
式中,(xy为近似二维液层的空间坐标;nl为此处液层的折射率;L为液层厚度。光线穿过液层进入空气后,因介质的改变将再次发生偏折,相应的偏折角度可表示
式中,ε'y为光线穿过液层到达空气中的偏折角度;na为外部环境中的折射率,近似为1。在面光源纹影系统中,传质纹影图像各点灰度(If)与相应背景灰度(Ik存在如下关系
式中,ΔI为传质纹影图像与背景的灰度差值;ak为背景图像对应的未被刀口遮挡的光源像宽度;f2为第二面凹面镜焦距。联立式(1~式(3,可得到灰度变化与液层折射率梯度间的关系
通过传质纹影图像的灰度变化分析可得到各点处的折射率梯度,进而对式(4进行积分运算,可获得对应的折射率分布[3,18]式中,积分起点的折射率n0已知,相应区域未受传质影响,位于分析区域下方。通过折射率与浓度的关系曲线的转换,可将求取的折射率分布转换为该液层的浓度
分布,25℃时乙醇-CO2体系的折射率变化关系[17]式中,c为乙醇中的CO2浓度。
1所示为实验体系物性[6,19],其中,ΔρsatΔσsat分别为饱和CO2的乙醇溶液与纯乙醇的密度差与界面张力差;csat为饱和CO2的乙醇溶液浓度。在乙醇吸收CO2的过程中,随着CO2浓度的增加,近界面处的液体密度增加表面张力减小,故该传质过程为Rayleigh不稳定,Marangoni稳定[2,6],便于单一地考Rayleigh对流。
3为吸收尚未开始时的液层纹影图像,此时液层图像亮暗一致,无对流产生。其中的矩形区域为选取的定量分析区域,气液界面处的黑色条纹是液相与光学玻璃间的接触角所致。
3.1乙醇吸收CO2过程中Rayleigh对流定性观察
通过不同时刻下的纹影图片对比可知,该吸收过程可分为湍动发生前后两个阶段。初始阶段,选取局部区域(图4)进行分析讨论。湍动未发生前(09s),暗条纹自气液接触开始出现在前面提及的黑色条纹下方,并逐渐向主体扩散。当到达一定深度(CO2的渗透距离),原本分布较为均匀的条纹出现变化(10s),其厚度变得不再一致,局部出现向下突出的趋势,其相邻区域相对变窄。突出的区域将以羽状流的形式向主体发展,至此,湍动开始。
湍动发生后,选取全部定量分析区域(图5)进行分析讨论。临近界面的区域内颜色相比于湍动发生前略微变亮,可以观察到沿气液界面分布着许多黑色点状区域,这些黑点在吸收过程中不断产生,并随着界面处液体的横向运动汇聚到羽状流中。此外,通过观察系列纹影图片可以看到,向下发展的羽状流会不断产生,并且会出现流股间的合并与发展,其相应的特征尺寸及深度不断增加。由于高密度液体向下的流动,会挤压其周围的液体,使得主体液体向上运动,并因此导致羽状流的弯曲变形,最终液相主体的纹影图像中对流胞不断合并交错。纹影图像直观反映出
Rayleigh对流加强了不同浓度液体的混合,加速了传质过程。3.2乙醇吸收CO2过程的液层浓度分布
对纹影图片进行定量分析,可获得相应时刻下的液层浓度分布,结果如图6所示。上界为气液界面,横坐标为图3分析区域的液层宽度,纵坐标为液层深度,不同图片中的各颜色代表不同的溶质浓度。初始阶段是以分子扩散为主导的传质过程,CO2以较为均匀的分布逐步向液层扩散,此时液层溶质浓度较低,且无对流出现。传质开始到10s的过程中,界面局部吸收速率不均而产生的随机浓度扰动,吹扫气进出系统而产生的随机外力扰动,以及装置边壁处对流体流动的固有限制使得液层失稳,产生的局部高浓度区域将在重力的驱动下产生向下运动趋势。此后界面湍动开始,高浓度液体以羽状流的形式向下发展。由图6(f可知,羽状流浓度分布由上到下逐渐减小,中间区域的液相浓度要比其边缘区域更大,且伴随着羽状流向下发展的过程中,主体液层的浓度逐渐增加。3.3界面液层浓度分析
在以往的研究中,界面浓度往往默认为气液平衡浓度,尚未有实验精确测量近界面浓度的真实变化过程。Baumann[20]Ma[21]分别通过实验和理论说明了气液界面传质阻力的存在,即在气液传质过程中,其界面并不能瞬间到达其平衡状态,真实的界面浓度要小于平衡浓度。
如图7所示,传质初期,界面浓度随气液接触时间增加而增加,且其横向的浓度分布较为一致。随着渗透过程的进行,传质体系的流体力学稳定性遭到破坏,界面处的扰动导致了非均匀传质的出现,使得界面浓度分布不再均一,高浓度区域因较大的密度差异将引发对流,向下发展的羽状对流胞将牵引近界层两侧液体,使得该区域传质得到加强。
实验体系中,乙醇溶液的表面张力随CO2浓度升高而降低。由界面浓度横向分布可知,湍动区域的浓度较高,表面张力较低不利于两侧液体向羽状流汇聚,但两股
羽状流间的区域因略低的浓度而具有较高的表面张力,所形成的表面张力差异将促进羽状对流胞的融合。3.4主体液层浓度分析
8显示的是不同时刻下主体液层各深度平均浓度的变化。在此以cavg/cI的变化曲线进行表征。在湍动发生前,截面平均浓度随液层深度迅速下降,其相应液层深度仅为界面下12mm。这说明在此液层中存有较大的传质阻力。湍动后,随Rayleigh对流的发展,相应的浓度下降区域变大,曲线的下降趋势变缓,说明相应的传质阻力变小。此外,如图8所示,在湍动后的近界面处出现了拐点,对应图9(b定量分析结果,在临近界面的液层中浓度分布不再是单一下降趋势。在前面的纹影图片(图5)中可知近界面处出现略亮的区域,即正折射率梯度,如图9(a所示,出现负的浓度梯度。形成这种现象的原因是部分被羽状流挤压向上运动的主体低浓度液体到达了此处,部分区域形成了负的浓度梯度。不同浓度液体在此混合,这也是Rayleigh对流能加强气液传质的一个主要原因。3.5临界Rayleigh数的计算
通常以特征数Ra来表征Rayleigh对流强度,其在传质Rayleigh对流中的表达式
理论临界Ra和理论临界引发时间tc反映不同物系的特性。Sparrow[22]通过对自然对流过程研究,得到不同温度边界下的理论临界RaTan[23]类比传热,基于渗透理论推导出传质中的最大Ra表达式,如式(8所示,对应渗透深度为最大渗透深度。Fu[14]使用LBM法模拟Rayleigh对流的引发和发展,当系统的浓度梯度即密度梯度达到某一临界点,并伴有适当扰动时,才会有对流产生。而此浓度梯度的临界点可以用临界Ra表示,并基于前人研究成果,提出不同Rac下的临界对流引发时间计算公式,即式(9
式中,,其中,假设流体满足Boussinesq近似,t为时间,cI值一般为标准状态
下反映物性特征的界面饱和浓度csat。而对于本文研究的乙醇-CO2体系中的Rayleigh对流,临界Ra1028[22],临界对流时间tc51.8s[15,24]本文将纹影图片与背景图片进行比较,通过液层相应灰度变化确定临界时刻最大渗透深度约为1.35mm,并以此作为特征尺寸,界面与主体的密度差值通过∆ρ=∆ρsatcI/csat来计算。
观察界面各处,可以发现在羽状流产生的区域存在着如图10所示的浓度分布变形,传质较快的区域有较大的Ra,并在高于临界数值后发生湍动。本文所计算得到平均临界Ra843(小于理论数值约18%),这与实际的界面浓度cI小于理论平衡浓度csat有关,如图7所示,虽有cI增大的趋势,但直至对流发生(10s)时,界面浓度始终低于饱和浓度。实验中获得临界对流时间为10s,远小于理论值51.8s,是因为实验值包含了吸收过程中扰动导致界面浓度分布不均匀的影响,而(9计算得到的临界Ra没有考虑此影响,所以实验值与理论值相差较远。3.6Rayleigh对流对传质的影响
对液层各处溶质含量的积分可获得整体液层的溶质总量
式中,m为乙醇吸收CO2的质量;Vcell为每个像素点所占用的液体体积。如前文所述,由于界面传质阻力的存在,界面浓度低于气液平衡浓度,真实的液相传质推动力[21]为界面液相侧的浓度cI与主体浓度cavg间的差值:cI−cavg,并非csat−cavg,则瞬时传质通量及液相传质系数的计算公式如下
Rayleigh对流发生前,如图6(a(b所示,界面附近浓度垂直方向递减,水平方向均一。如图11所示,液相传质系数随着溶质渗透而逐渐减小,趋势符合菲克第二定律,在此阶段主要是分子扩散。10s后湍动发生,传质系数迅速增大。与此同时,如图7所示,随着界面浓度的增大,传质通量变大,使得吸收总量快速增加。如图6(e(f所示,伴随着羽状对流在界面不同位点相继产生,并不断在主体中发展耗散,图11中液相瞬时传质系数产生无规则的起伏变化。在此阶段以
对流传质为主导,辅以分子扩散过程。总体来讲,湍动后传质速率大大增加,Rayleigh对流较为明显地加强了传质过程。
在获取清晰纹影图像的同时利用定量分析获得浓度场,进而获得瞬时传质系数随时间变化,更为深入地探讨了Rayleigh对流的形成与发展过程。湍动发生前,乙醇CO2的吸收过程以分子扩散为主,界面附近浓度垂直分布,水平方向均一,瞬时传质系数逐渐降低;湍动发生后,失稳界面出现向下运动的羽状流股,瞬时传质系数猛增,伴随流股的融合与耗散以及界面处新流股的出现,瞬时传质系数出现起伏变化。在界面传质阻力的影响下,气液界面在传质开始后并未瞬间到达平衡,而是界面浓度随时间增大,且扰动使得浓度分布不再均一,进而导致对流引发,故理论临界对流引发时间和由此计算的临界Ra与理论值存在差异。在吸收液层的纹影图像中观测到近界层处的较亮区域,对此进行定量分析,发现部分被羽状流挤压向上运动的低浓度主体液体与近界层高浓度流体在此混合形成负的浓度梯度(拐点),更新了近界层流体,从而加强气液传质。
[1]RAYLEIGHL.Onconvectioncurrentsinahorizontallayeroffluid,whenthehighertemperatureisontheunderside[J].TheLondon,Edinburgh,andDublinPhilosophicalMagazineandJournalofScience,1916,32(192:529-546.[2]余国琮,袁希钢.化工计算传质学导论[M].天津:天津大学出版社,2011:263-327.YUGC,YUANXG.IntroductiontoComputationalMassTransfer[M].Tianjin:TianjinUniversityPress,2011:263-327.[3]PANIGRAHIPK,MURALIDHARK.SchlierenandShadowgraphMethodsinHeatandMassTransfer[M].NewYork:Springer,2012:23-46.[4]沙勇,李樟云,林芬芬,.气液传质界面湍动现象投影观察[J].化工学报,2010,61(4:844-847.SHAY,LIZY,LINFF,etal.Shadowgraphobservation
oninterfacialturbulencephenomenaingas-liquidmasstransfer[J].CIESCJournal,2010,61(4:844-847.[5]GUOY,YUANXG,ZENGAW,etal.Measurementofliquidconcentrationfieldsnearinterfacewithcocurrentgas-liquidflowabsorptionusingholographicinterferometry[J].Chin.J.Chem.Eng.,2006,14(6:747-753.[6]OKHOTSIMSKIIA,HOZAWAM.Schlierenvisualizationofnaturalconvectioninbinarygas-liquidsystems[J].Chem.Eng.Sci.,1998,53(14:2547-2573.[7]ARENDTB,DITTMARD,EGGERSR.InteractionofinterfacialconvectionandmasstransfereffectsinthesystemCO2-water[J].Int.J.HeatMassTransfer,2004,47(17/18:3649-3657.[8]王勇,张泽廷.气液传质界面湍动现象的实验及分析[J].北京化工大学学报,2002,29(2:11-14.WANGY,ZHANGZT.Opticalobservationandanalysisoftheinterfacialturbulenceinmasstransferprocessesofgas-liquidsystems[J].JournalofBUCT,2002,29(2:11-14.[9]SUNZF,YUKT,WANGSY,etal.Absorptionanddesorptionofcarbondioxideintoandfromorganicsolvents:effectofRayleighandMarangoniinstability[J].Ind.Eng.Chem.Res.,2002,41(7:1905-1913.[10]SHAY,CHENGH,YUYH.Numericalanalysisofthegas-liquidabsorptionprocessaccompaniedbyRayleighconvection[J].Chin.J.Chem.Eng.,2002,10(5:539-544.[11]FARAJZADEHR,SALIMIH,ZITHAPLJ,etal.Numericalsimulationofdensity-drivennaturalconvectioninporousmediawithapplicationfor
CO2injectionprojects[J].Int.J.HeatMassTransfer,2007,50(25/26:5054-5064.
[12]SHAY,CHENH,YINY,etal.CharacteristicsoftheMarangoniconvectioninducedininitialquiescentwater[J].Ind.Eng.Chem.Res.,2010,49(18:8770-8777.[13]CHENSY,FUB,YUANXG,etal.LatticeBoltzmannforsimulationofsolutalinterfacialconvectioningas-liquidsystem[J].Ind.Eng.Chem.Res.,2012,51(33:10955-10967.[14]FUB,LIUBT,YUANXG,etal.ModelingofRayleighconvectioningas-liquidinterfacialmasstransferusinglatticeBoltzmannmethod[J].Chem.Eng.Res.Des.,2013,91(3:437-447.[15]GUOK,LIUC,CHENS,etal.SpatialscaleeffectsonRayleighconvectionandinterfacialmasstransfercharacteristicsinCO2absorption[J].Chem.Eng.Technol.,2015,38(1:23-32.[16]LIUCX,ZENGAW,YUANXG,etal.Experimentalstudyonmasstransferneargas-liquidinterfacethroughquantitativeSchlierenmethod[J].Chem.Eng.Res.Des.,2008,86(2:201-207.[17]于海路,曾爱武.气液传质过程中Marangoni对流的观测与定量分析[J].工学报,2014,65(10:3760-3768.YUHL,ZENGAW.VisualizationandquantitativeanalysisforMarangoniconvectioninprocessofgas-liquidmasstransfer[J].CIESCJournal,2014,65(10:3760-3768.[18]ALVAREZ-HERRERAC,MORENO-HERNÁNDEZD,BARRIENTOS-GARCÍAB,etal.TemperaturemeasurementofairconvectionusingaSchlierensystem[J].Optics&LaserTechnology,2009,41(3:233-240.
[19]TAKAHASHIM,KOBAYASHIY,TAKEUCHIH.Diffusioncoefficientsandsolubilitiesofcarbondioxideinbinarymixedsolvents[J].J.Chem.Eng.Data,1982,27(3:328-331[20]BAUMANNKH,MÜHLFRIEDELK.Masstransferandconcentrationprofilesnearphaseboundaries[J].Therm.Sci.,2010,40(5:425-436.[21]MAYG,YUGC,LIHZ.Noteonthemechanismofinterfacialmasstransferofabsorptionprocesses[J].Int.J.HeatMassTransfer,2005,48(16:3454-3460.[22]SPARROWEM,GOLDSTEINRJ,JONSSONVK.Thermalstabilityinahorizontalfluidlayer:effectofboundaryconditionsandnon-lineartemperatureprofile[J].J.FluidMech.,1964,18(4:522-523.[23]TANKK,THORPERB.Gasdiffusionintoviscousandnon-Newtonianliquids[J].Chem.Eng.Sci.,1992,47(13/14:3565-3572.[24]TANKK,THORPERB.Theonsetofconvectioninducedbybuoyancyduringgasdiffusionindeepfluids[J].Chem.Eng.Sci.,1999,54(19:4179-4187.2016-05-18收到初稿,2016-08-03收到修改稿。联系人:曾爱武。

免费下载 Word文档免费下载: 安全验证

  • 29.8

    ¥45 每天只需1.0元
    1个月 推荐
  • 9.9

    ¥15
    1天
  • 59.8

    ¥90
    3个月

选择支付方式

  • 微信付款
郑重提醒:支付后,系统自动为您完成注册

请使用微信扫码支付(元)

订单号:
支付后,系统自动为您完成注册
遇到问题请联系 在线客服

常用手机号:
用于找回密码
图片验证码:
看不清?点击更换
短信验证码:
新密码:
 
绑定后可用手机号登录
请不要关闭本页面,支付完成后请点击【支付完成】按钮
遇到问题请联系 在线客服