三维变形管内降膜蒸发特性.pdf

流 体 机 械2023 年 7 月76 第 51 卷第 7 期 收稿日期：2022-05-17 修稿日期：2023-02-27基金项目：国家重点研发计划项目（2022YFB4202005）；广州市重点研发计划项目（202206010124）；吉林省与中国科学院科技合作高 技术产业化专项资金项目（2022SYHZ0027）；中科院创新研究院-榆林学院联合基金项目doi：10.3969/j.issn.1005-0329.2023.07.011应用技术三维变形管内降膜蒸发特性陈杭生1，2，3，朱冬生1，2，陈二雄1，2，林成迪4，胡广涛5（1.中国科学院 广州能源研究所，广州 510640；2.中国科学院 可再生能源重点实验室，广州 510640；3.中国科学院大学，北京 100049；4.三峡新能源阳江发电有限公司，广东阳江 529532；5.榆林学院，陕西榆林 719000）摘 要：为了提高降膜蒸发器的传热性能，将三维变形管用于MVR 降膜蒸发器中，利用FLUENT 软件探究了管内液膜的分布特点，分析了长、短轴比对三维变形管的传热性能的影响。结果表明：相比圆管，三维变形管内液膜的厚度会逐渐减小，流体速度场与温度场的协同性相对来说也更好；随着长、短轴比的增大，三维变形管内液膜的厚度会逐渐减小，液膜速度更为缓慢，而且进一步提高了速度场与温度场的协同性，还增强了液膜的扰动性，最终努塞尔数从 0.249 提高到 0.276。将三维变形管应用于MVR 蒸发器中，具有良好的换热效果，可节省换热材料。关键词：三维变形管；降膜蒸发；液膜厚度；长、短轴比中图分类号：TH12 文献标志码：A Characteristics of falling film evaporation in three-dimensional deformed tubeCHENHangsheng1，2，3，ZHUDongsheng1，2，CHENErxiong1，2，LINChengdi4，HUGuangtao5（1.GuangzhouInstituteofEnergyConversion，ChineseAcademyofSciences，Guangzhou 510640，China；2.KeyLaboratoryofRenewableEnergy，ChineseAcademyofSciences，Guangzhou 510640，China；3.UniversityofChineseAcademyofSciences，Beijing 100049，China；4.ThreeGorgesRenewablesYangjiangPowerGenerationCo.，Ltd.，Yangjiang 529532，China；5.YulinUniversity，Yulin 719000，China）Abstract：Inordertoimprovetheheattransferperformanceofthefallingfilmevaporator，athree-dimensionaldeformedtubewasusedintheMVRfallingfilmevaporator.ThedistributioncharacteristicsoftheliquidfilminthetubewereexploredbyFLUENT，andtheinfluenceofmajor-minoraxisratioontheheattransferperformanceofthethree-dimensionaldeformedtubewasanalyzed.Theresearchresultsshowthatthethicknessoftheliquidfilminthethree-dimensionaldeformedtubewillgraduallydecreasecomparedwiththeroundtube，andthesynergybetweenthefluidvelocityfieldandthetemperaturefieldisrelativelybetter；withtheincreaseinthemajor-minoraxisratio，thethicknessoftheliquidfilminthethree-dimensionaldeformedtubewillgraduallydecrease.Thethicknessoftheliquidfilmwillgraduallydecrease，thespeedoftheliquidfilmisslower，andthesynergybetweenthevelocityfieldandthetemperaturefieldisfurtherimproved，andtheturbulenceoftheliquidfilmisenhanced.Finally，theNusseltnumberisincreasedfrom0.249to0.276.Applyingthree-dimensionaldeformedtubeinMVRevaporatorhasgoodheatexchangeeffectandcansaveheatexchangematerial.Key words：three-dimensionaldeformedtube；fallingfilmevaporation；liquidfilmthickness；major-minoraxisratio0 引言降膜蒸发器是工业上常见的换热设备，具有成膜均匀、蒸发速率稳定、处理量大、耗能低等优点，因此广泛应用于化学工业、制冷制热以及海水淡化等方面1。近几十年来，众多学者在实验和理论推导及数值模拟方面对降膜蒸发过程做了大量的研究工作。77PARKEN 等2通过实验研究，分析水在 2种不同管径光滑管下的降膜蒸发情况，并且给出蒸发和沸腾条件下的传热关联式。2003年，FEDDAOUI 等3通过求解液膜与气膜耦合的控制方程，研究表明流动液膜的对流传热是壁面排热的主要机制。郝丽等4数值模拟结果表明，竖直降膜蒸发器加热管具有液膜薄，均匀性好的特点，管表面容易形成柱状流，液膜厚度为 1.0mm。BIGHAM 等5采用两相流法研究降膜蒸发现象，结果表明平均换热系数随着入口质量流量及温差的增大而增大。WANG 等6使用FLUENT 软件对竖直管内R113 在不同结构条件下的气液两相逆流流动特性进行了数值研究，研究表明液膜波动随着管长的增加而增强，当管长大于 12m 时，不能完全湿润壁面。谢迎春等7通过数值模拟得出，竖直管内降膜过程中逆向气流扰动会改变稳定段在管内的分布位置与范围。张天娇等8通过试验研究了布液器参数对降膜流动的影响，为避免有“干斑”出现，布液器环隙间距不应小于0.5mm。本文将对液膜在不同长、短轴比的三维变形管内蒸发特性进行数值模拟，探究三维变形管内液膜流动以及分布特点，并对比分析圆管和三维变形管内的场协同性，为此类降膜蒸发器的设计、优化提供依据。1 三维变形管三维变形管结构如图 1（a）所示，其中结构参数A 为长轴，B 为短轴，P 为扭矩。三维变形管由基础圆管压扁后扭曲而成，沿管长方向呈现出一种螺旋形状，液体在变形管内以螺旋方式流动，流动过程中会在垂直于主流的方向上产生二次流，加强液体的扰动，有利于破坏边界层，提高液体的传热效率9。图 1（b）为三维变形管自支撑示意，经钢带捆扎后的三维变形管管束在最大变径凸点处相互接触，可形成自支撑结构，减少振动10-11。（a）结构示意 （b）自支撑示意图 1 三维变形管Fig.1 Three-dimensionaldeformedtube2 模型与计算2.1 几何模型本文对三维变形管进行降膜蒸发的过程采用插件型布膜器使液体成膜。布膜器和管子之间存在一定的间隙，从而使液体能够沿着间隙向下流动铺展成膜，建模区域如图 2 所示。图 2 建模区域示意Fig.2 Schematicdiagramofmodelingarea利用SolidWorks 软件建立相应的圆管模型。在本文研究中，所有换热管的长度都为 620mm，布膜间隙为 1.2mm。采用圆管的直径为 25mm，三维变形管一共有 4 种模型：螺距P 为 300mm，长、短轴比（即A/B）分别为 1.16，1.30，1.47，1.67。其具体的几何模型的结构尺寸见表 1。表 1 换热管几何模型结构尺寸Tab.1 Geometricmodelstructuresizeofheattransfertube型号换热管尺寸/mm管型ID=25，L=620圆管IIP=300，L=620，A/B=1.16三维变形管IIIP=300，L=620，A/B=1.30三维变形管IVP=300，L=620，A/B=1.47三维变形管VP=300，L=620，A/B=1.67三维变形管在对换热管几何模型进行建立时忽略换热管和布膜器的壁厚，以其中一个三维变形管为例展示其几何模型，如图 3 所示。图 3 三维变形管几何模型示意Fig.3 Schematicdiagramofgeometricmodelofthree-dimensionaldeformedtube陈杭生，等：三维变形管内降膜蒸发特性78FLUID MACHINERYVol.51，No.7，20232.2 数值模型及边界条件采用二维双精度求解器进行求解计算。选用VOF 模型和RNG-模型，并采用标准壁面函数法对近壁区进行处理。数值求解方法选用分离解法，时间离散格式采用非稳态，压力-速度耦合采用PISO 算法进行求解。此外，压力项采用PRESTO 算法，动量、能量以及湍流参量的求解均采用二阶迎风格式，气液界面的追踪采用精度较高的Geo-Reconstruction 界面重构技术。本文模拟过程中采用的蒸发介质水。入口边界条件采用速度入口，并且定义了相应的速度以及入口温度；出口边界条件选用压力出口，并且定义了相应的回流压力和回流温度；壁面条件给定了恒定温度并将壁面设置为无滑移壁面。2.3 网格划分及无关性验证采用ICEM 软件进行网格划分，选用精度较高的六面体结构化网格，图 4 示出V 号三维变形管的局部网格示意。为了保证近壁面处计算精度以及气液界面的准确捕捉，对近壁面液膜流动区域进行加密，以保证计算结果的精度。图 4 结构化网格局部示意Fig.4 SchematicdiagramofPartialstructuredmesh液膜平均厚度随网格数量变化如图 5 所示。图 5 液膜平均厚度随网格数量变化Fig.5 Changeofaveragethicknessoftheliquidfilmwiththenumberofgrids为了满足计算结果对网格数量的要求，进行网格无关性验证。对近壁面处的网格进行不同程度的加密，得到网格数量为 50 万140 万。由图 5可知，当网格数量达到 115 万时，增加网格数量对冷态模拟（无传热）情况下的液膜平均厚度影响不大，在兼顾计算精度和时间成本的前提下，选用网格数量为 115 万进行计算。2.4 模拟结果的正确性为了确认数值模拟的正确性，对饱和水在圆管内进行降膜蒸发模拟，与前人的试验结果进行对比分析。CHUN 等12降膜蒸发关联式被广泛使用，因此将模拟结果所得的传热系数与其经验公式值进行对比，其结果如图 6 所示。模拟值比经验值大，主要是因为模拟过程对模型进行了简化，避免了介质纯度、管壁热阻、测量仪器等因素的影响，但两者的最大误差仅为 5%，说明了模拟过程中选用方法的可靠性和准确性。图 6 努塞尔数Nu 模拟值与经验值对比Fig.6 ComparisonbetweennumericalvalueandempiricalvalueofNu3 结果与分析3.1 模拟结果初步分析针对I 号圆管和V 号三维变形管的几何模型，流体速度入口为 1.2m/s，流体进口温度设置为372.15K，壁面温度为 374K。3.1.1 管内速度分布从图 7 中管内速度分布可看出，I 号圆管内流体在前半段分布较为均匀，速度沿管子中心呈同心圆状，流动较为平稳；后半段流体速度分布产生了一定的波动；在V 号三维变形管内，流速呈现扭曲状，主要是因为三维变形管的螺旋管型使得79管内流体在沿轴向运动时会同时受到离心力的作用，导致流体沿着管子扭曲方向发生偏转，越靠近长短轴交界处，所受到的扰动就越激烈。图 7 管内速度分布Fig.7 Distributionchartofvelocityintube3.1.2 管内液膜分布图 8 示出沿管长方向不同截面的管内液膜分布，在I 号圆管中不同截面液膜的分布基本一致，在换热管前半段分布较为均匀，而后半段会产生一定的波动。在V 号三维变形管内，不同截面液膜的分布规律基本一致，顺着换热管扭曲的方向，液膜会聚集在长轴与短轴之间，形成较厚的液膜，而在相反的方向所形成的液膜厚度就较薄。图 8 管内液膜分布Fig.8 Distributionchartofliquidfilmintube3.1.3 液膜厚度分布取Y=0 剖切面进行分析，图 9 示出液膜厚度沿管长方向变化。图 9 液膜厚度沿管长方向变化Fig.9 Changeofliquidfilmthicknessalongthelengthofthetube在圆管内前半段液膜流动较平稳，而后半段液膜会产生激烈的波动，是因为后半段的液膜开始蒸发，所产生的蒸汽会对液膜的分布有所影响。在三维变形管内，因扭矩的存在而使得液膜的分布具有周期性的变化规律，且呈现相应的“波峰”和“波谷”。在三维变形管内沿管长方向的液膜平均厚度为 1.18mm，而在圆管内平均值为1.36mm，前者比后者减小了 15.3%，说明三维变形管内液膜的传热热阻会比圆管的小，更有利于热量的传递。3.1.4 液膜表面速度分布取Y=0 剖切面进行分析，图 10 示出液膜表面速度沿管长方向变化。在圆管内，前半段液膜表面速度分布较为平稳，而在后半段会出现一定的波动。在三维变形管内液膜表面速度分布因扭矩的存在而呈现周期性的变化规律，分布更为平稳。三维变形管内的液膜表面速度平均值为 1.04m/s，而圆管内的液膜表面速度平均值为1.16m/s，前者比后者减小了 10.3%，说明在三维变形管内液膜所受到的阻力比圆管的大，因此流动得更为缓慢。图 10 液膜表面速度沿管长方向变化Fig.10 Changeofsurfacevelocityofliquidfilmalongthelengthofthetube3.1.5 液膜表面温度分布取Y=0 剖切面进行分析，图 11 示出液膜表面温度沿管长方向变化。在圆管内，液膜表面温度沿着管长的方向在逐渐地增大，且在后半段会出现微小的波动；沿着管长方向温度的变化率在逐渐地降低。在三维变形管内，液膜表面温度总体上沿着管长的方向呈现逐渐增长的趋势，但在中间管段因为螺旋流动出现小幅度的下降，然而温度依然处于高位。在出口处三维变形管内液膜表面温度比圆管的高，说明在三维变形管内液膜温度比圆管变化得更快。陈杭生，等：三维变形管内降膜蒸发特性80FLUID MACHINERYVol.51，No.7，2023图 11 液膜表面温度沿管长方向变化Fig.11 Changeofsurfacetemperatureofliquidfilmalongthelengthofthetube3.1.6 液膜速度与温度梯度协同角分布换热管的传热性能与速度和温度梯度之间的协同角有关，协同角的表达式为13：=()arccos?uTuT式中，u为速度矢量；T?为温度梯度矢量。减薄热边界层以及增加流体中的扰动其实质都是要减小速度与温度梯度间的夹角。取Y=0 剖切面进行分析，并在该剖切面上取液膜的等速线u=0.4m/s，从图 12 中可以看出，在圆管内前半段速度与温度梯度协同角的分布较为平稳，基本没有发生太大的变化，而在后半段会出现一定的波动。在三维变形管内速度与温度梯度协同角波动得非常激烈，变化幅度相对圆管较大。总体来看，三维变形管内的速度与温度梯度协同角会比圆管的小，说明在三维变形管内速度场与温度场的协同性比圆管的好，更有利于热量的传递。图 12 液膜速度与温度梯度协同角沿管长方向变化Fig.12 Changeofsynergisticanglebetweenliquidfilmvelocityandtemperaturegradientalongthelengthofthetube3.2 长、短轴比对传热性能的影响针对IIV 号三维变形管的几何模型，流体速度入口为 1.2m/s，流体进口温度设置为 372.15K，壁面温度为 374K。3.2.1 长、短轴比对液膜厚度的影响取Y=0 剖切面进行分析，并对液膜厚度取平均值，图 13 示出了不同长、短轴比下液膜厚度的变化情况。随着长、短轴比的增大。液膜厚度在逐渐地减小，液膜的传热热阻就会越小，越有利于热量的传递。长、短轴比从 1.16 增加到 1.67 时，液膜厚度从 1.48mm 减小到了 1.24mm，降低了 16.2%。图 13 液膜厚度随长、短轴比变化Fig.13 Changeofliquidfilmthicknesswithmajor-minoraxisratio3.2.2 长、短轴比对液膜表面速度的影响取Y=0 剖切面进行分析，并对液膜表面速度取平均值，图 14 示出了不同长、短轴比下液膜表面速度的变化情况。随着长、短轴比的增大。液膜表面速度在逐渐地减小，说明长、短轴比越大，液膜所受到的流动阻力就会越大，从而流体的流动就越缓慢。长、短轴比从 1.16 增加到 1.67 时，液膜表面速度从 1.28m/s 减小到了 1.13m/s，降低了 11.7%。图 14 液膜表面速度随长、短轴比变化Fig.14 Changeofsurfacevelocityofliquidfilmwithmajor-minoraxisratio813.2.3 长、短轴比对液膜速度与温度梯度协同角的影响取Y=0 剖切面上液膜等速线u=0.4m/s 进行分析，图 15 示出该等速线上液膜速度与温度梯度协同角平均值随长、短轴比变化。随着长、短轴比的增大，液膜速度与温度梯度协同角变得越来越小，说明长、短轴比越大，速度场与温度场的协同性就越好，越有利于热量的传递。图 15 液膜速度与温度梯度协同角随长短轴比变化Fig.15 Changeofsynergisticanglebetweenliquidfilmvelocityandtemperaturegradientwithmajor-minoraxisratio3.2.4 长、短轴比对努塞尔数的影响图 16 示出不同长、短轴比下努塞尔数的变化情况，当长、短轴比从 1.16 变化到 1.67 时，努塞尔数从 0.249 提高到 0.276，增加 10.8%。由前面分析知，长、短轴比的增大不仅有利于减薄液膜厚度，而且提高了速度场与温度场的协同性，还增强液膜的扰动性，这 3 方面都有利于提高液膜的换热性能，因此液膜努塞尔数会随着长、短轴比的增大而增大。图 16 努塞尔数Nu 随长、短轴比变化Fig.16 ChangechartofNuwithmajor-minoraxisratio4 应用实例某纤维素生产企业的废水零排放项目，采用MVR 强制循环蒸发结晶系统，蒸发器采用钛管TA2 作为换热元件，壳体材质为 316L 不锈钢。三维变形管MVR 蒸发器优化设计参数与普通圆管蒸发器的设计参数对比见表 2，通过对比可以发现，在总换热量同为 10753.5kW 的情况下，采用三维变形管MVR 蒸发器较普通圆管蒸发器设计工况计算综合传热系数提高了 35.8%，换热面积减少 26.7%，体积缩小 35.7%。通过现场标定证实已达到设计要求，表明采用三维变形管MVR 蒸发器可以节省换热材料和运输成本，具有良好的强化传热效果。表 2 性能参数对比Tab.2 Comparisonchartofperformanceparameters参数普通圆管蒸发器三维变形管蒸发器废水流量/（m3 h-1）38503850废水进口温度/97.8597.85废水出口温度/100100管内压力/kPa150150管内压降/kPa30.437.2蒸汽流量/（103kg h-1）17.3617.36蒸汽温度/110110壳程压力/kPa143143总换热量/kW10753.510753.5换热面积/m214821086总传热系数/（W m-2 K-1）8751188裕量/（%）2019.7流程1/11/1换热管尺寸/（mmmmmm）321.512000321.59000A=38，B=22，P=300换热管数目/根12281200筒体尺寸/（mmmm）2100140001900110005 结论（1）相对于圆管，三维变形管内液膜的流动为螺旋流动，速度更为缓慢，且具有一定的周期性的变化规律。在三维变形管内液膜的厚度会比圆管减少 15.3%，温度场与速度场的协同性比圆管的好，有利于提升管内的传热效率。陈杭生，等：三维变形管内降膜蒸发特性82FLUID MACHINERYVol.51，No.7，2023（2）在长、短轴比为 1.161.67 的范围内，随着长、短轴比的增大，三维变形管内液膜的厚度会逐渐地减小，速度更为缓慢，温度场与速度场的协同性越好，最终努塞尔数会从 0.249 提高到0.276。（3）将三维变形管应用于MVR 蒸发器中，具有良好的换热效果，可节省换热材料。参考文献：1 ZHAOCY，JIWT，JINPH，etal.HeattransfercorrelationofthefallingfilmevaporationonasinglehorizontalsmoothtubeJ.AppliedThermalEngineering，2016，103：177-186.2 PARKENWH，FLETCHERLS，HANJC，etal.HeattransferthroughfallingfilmevaporationandboilingonhorizontaltubesJ.JournalofHeatTransfer，1990，112，3（3）：744-750.3 FEDDAOUIM，MIRA，BELAHMIDIE.CocurrentturbulentmixedconvectionheatandmasstransferinfallingfilmofwaterinsideaverticalheatedtubeJ.InternationalJournalofHeatandMassTransfer，2003，46（18）：3497-3509.4 郝丽，肖晓明，张吕鸿，等.竖直管降膜蒸发器流动与传热过程的数值模拟 J.化学工业与工程，2014，31（2）：60-65.HAOL，XIAOXM，ZHANGLH，etal.NumericalsimulationoftheflowandheattransferprocessofverticaltubefallingfilmevaporatorJ.ChemicalIndustryandEngineering，2014，31（2）：60-65.5 BIGHAMS，KOUHIKAMALIR，ABADIS.Two-phaseflownumericalsimulationandexperimentalverificationoffallingfilmevaporationonahorizontaltubebundleJ.Desalination&WaterTreatment，2014，55（8）：1-14.6 WANGW，LIAOY，YANY，etal.NumericalstudyonfallingfilmflowingcharacteristicsofR113insideverticaltubeunderdifferentstructuralconditionsJ.ChineseJournalofChemicalEngineering，2020，28（1）：23-32.7 谢迎春，刘军，徐畅，等.竖直降膜蒸发冷凝管内液膜流动特性研究 J.流体机械，2022，50（6）：39-46.XIE Y C，LIU J，XU C，et al.Study on flowcharacteristicsofliquidfilminverticalfallingfilmevaporativecondensertubeJ.FluidMachinery，2022，50（6）：39-46.8 张天娇，马良栋，张吉礼，等.布液器入口参数对竖直管外降膜流动特性的影响J.流体机械，2020，48（9）：1-5.ZHANGTJ，MALD，ZHANGJL，etal.EffectofinletparametersofliquiddistributoronfallingfilmflowcharacteristicsoutsideverticaltubeJ.FluidMachinery，2020，48（9）：1-5.9 TODDL.Somecommentsonsteady，laminar-flowthroughtwistedpipesJ.JournalofEngineeringMathematics，1977，11（1）：29-48.10 TANXH，ZHUDS，ZHOUGY，etal.ExperimentalandnumericalstudyofconvectiveheattransferandfluidflowintwistedovaltubesJ.InternationalJournalofHeatandMassTransfer，2012，55（17-18）：4701-4710.11 莫逊，朱冬生，张洁娜.扭曲椭圆管在MVR 系统降膜蒸发器上的应用研究J.化学工程，2016，44（9）：24-28.MOX，ZHUDS，ZHANGJN.PracticalresearchontwistedellipticaltubeinthefallingfilmevaporatorofMVRsystemJ.ChemicalEngineering（China）.2016，44（9）：24-28.12 CHUNKR，SEBANRA.HeattransfertoevaporatingliquidfilmsJ.J.HeatTransfer，1971，96：391-396.13 杨世铭，陶文铨.传热学M.4 版.北京：高等教育出版社，2016.YANGSM，TAOWQ.HeattransferM.4thedition.Beijing：HigherEducationPress，2016.作者简介：陈杭生（1990），男，硕士研究生，主要从事余热回收利用、强化传热及节能环保研究，E-mail：。通信作者：朱冬生（1964），男，教授，博士，博士生导师，主要从事余热回收利用、强化传热及节能环保研究，通信地址：510640 广东省广州市天河区能源路 2 号中国科学院广州能源研究所，E-mail：。本文引用格式：陈杭生，朱冬生，陈二雄，等.三维变形管内降膜蒸发特性 J.流体机械，2023，51（7）：76-82.CHENHS，ZHUDS，CHENEX，etal.Characteristicsoffallingfilmevaporationinthree-dimensionaldeformedtubeJ.FluidMachinery，2023，51（7）：76-82.