搜索
您的当前位置:首页正文

基于CFD数值模拟的换热器外导流筒结构选型

来源:爱够旅游网
—116—

      

石 油 机 械

CHINAPETROLEUMMACHINERY

2017年 第45卷 第2期

  

◀炼油与化工机械▶

基于CFD数值模拟的换热器外导流筒结构选型

王庆锋 郝 帅 李 凯 李 中

(北京化工大学机电工程学院)

摘要:管壳式换热器壳程流场均匀分布能够改善传热效果,外导流筒的合理选型对流场均匀性有重要影响。然而,外导流筒结构多样,其选型设计多依靠经验,缺乏指导性方法。鉴于此,以筒体直径2200mm、入口接管外径720mm的竖直管降膜式蒸发器换热段为研究对象,设计了外导流筒,建立了3种结构的三维实体模型。基于CFD数值模拟的方法,以平均速度和流场均匀性为衡量指标,研究了3种导流筒结构壳程纵截面和入口截面的流体分布效果,证明了外导流筒合理选型确实能够改善换热器壳程流场均匀性。定性结合定量分析得出:变截面圆弧过渡式外导流筒平均速度相对合理,流场均匀性最佳,是较适宜的导流筒结构。研究结果对外导流筒的选型设计具有重要意义。

关键词:管壳式换热器;外导流筒;CFD;数值模拟;选型

中图分类号:TE965 文献标识码:A doi:10􀆰16082/j􀆰cnki􀆰issn􀆰1001-4578􀆰2017􀆰02􀆰026

CFD⁃basedSelectionofExternalFlowDistributor

StructureforHeatExchanger

(CollegeofMechanicalandElectricalEngineering,BeijingUniversityofChemicalTechnology)

WangQingfeng HaoShuai LiKai LiZhong

improvetheheattransfereffect.Therationalselectionoftheexternalflowdistributorhasimportantinfluenceontheuniformityoftheflowfield.However,theexternalflowdistributorhasvariousstructures,anditsdesignismoredependentonexperienceandlackofguidancemethod.Inviewofthis,three⁃dimensionalsolidmodelswiththreekindsofstructureshavebeenestablished,andtheexternalflowdistributorhasbeendesignedfortheheatexchangeof720mm.BasedonCFDnumericalsimulationmethod,andtakingtheaveragevelocityandtheuniformityofflowofthethreekindsofexternalflowdistributorhavebeenstudied.Itisprovedthattherationalselectionoftheexter⁃transitionexternalflowdistributorisrelativelyreasonableandtheuniformityoftheflowfieldisthebest,thuspres⁃entingasuitableflowdistributorstructure.Thestudyresultsareofgreatsignificanceforthetypeselectionofexter⁃nalflowdistributor􀆰lection

Keywords:shell⁃and⁃tubeheatexchanger;externalflowdistributor;CFD;numericalsimulation;typese⁃

Abstract:Theuniformdistributionofflowfieldintheshell⁃sideoftheshell⁃and⁃tubeheatexchangercould

sectionoftheverticaltubefalling⁃filmevaporatorwithacylinderdiameterof2200mmandtheinletpipediameterfieldastheindexes,thefluiddistributionsalongthelongitudinalsectionandtheinletcrosssectionoftheshellpassnalflowdistributorcouldimprovetheflowfielduniformityoftheshellpassofheatexchanger.Thecombinationofthequantitativeandqualitativeanalysisshowedthattheaveragevelocityofthevariable⁃cross⁃sectioncircular⁃arc

国石油科技创新基金项目“油气长输管道系统泄漏风险评估方法及应急救援策略”(2015D-5006-0606);北京市重点学科建设项目“城市管网安全监测预警及事故防治”(21510028008)。

∗基金项目:国家高技术研究发展计划项目“基于物联网的石化工业风险识别与安全应急系统研发与示范应用”(2014AA041806);中

Copyright©博看网 www.bookan.com.cn. All Rights Reserved.2017年 第45卷 第2期王庆锋等:基于CFD数值模拟的换热器外导流筒结构选型      —117—

0 引 言

管壳式换热器外导流筒起到减轻流体介质对换热管冲蚀、均匀分布入口介质、强化传热性能和膨导流结构的导流性能。杜明等[4]阐述了管壳式换热器外导流筒结构的设计依据,介绍了相关国际标准及专利内容。张立文[5]以某乙烯装置中环形外导流筒换热器为例,介绍了环形外导流筒结构尺寸的设计。郭雪华[6]详细阐述了换热器外导流筒的结构设计原则,并深入分析了各结构尺寸的计算方法和相互制约关系。古新等[7-8]研究了夹套式变截面导流筒换热器外导流筒的优化设计方法。王庆锋等[9]采用数值模拟方法研究了管壳式换热器传热效率的影响因素。吴金星等[10]研究了纵向流管壳式换热器的强化传热方法,指出进、出口处设置导流筒可使传热效率大大提高。壳程流体分配的不均匀会引起换热器效率下降[11],李平等[12]基于AN⁃SYSFLOTRANCFD分析了高压加热器蒸汽冷却段经过导流筒后,壳程流体流量的均匀性对换热效率的影响。然而,外导流筒结构多样,其选型设计目前依然缺乏指导性方法。

笔者以筒体公称直径为2200mm、入口接管公称直径为700mm的竖直管降膜式蒸发器换热段为研究对象,设计了外导流筒,建立了3种结构模型,基于CFD数值模拟,以壳程纵截面和入口截面的平均速度及流场均匀性为衡量指标,选出了较胀节作用[1]。张哲等[2-3]研究了板翅式换热器不同

适宜的外导流筒结构。

1 外导流筒几何尺寸的确定

竖直管降膜式蒸发器换热段筒体内径为2200mm,入口接管ø720mm×8,负压饱和蒸汽流速75导流筒结构如图1所示。根据标准及相关文献[6,13]设计计算,得出基本尺寸如表1所示。

m/s,换热管规格为ø45mm×3,管数1100根,管长6000mm,正三角形布管,单管程。开槽式外

图1 开槽式外导流筒结构

Fig􀆰1 Structureoftheslottedtypeexternalflowdistributor

表1 外导流筒基本尺寸

d/mm70426K1L1/mm10981􀆰1A0/mm2L2/mm150A1/mm2L3/mm3242818H/mm188428do/mmb/mm462720H/doK2Table1 GeneraldimensionsoftheexternalflowdistributorB/mm11390􀆰262􀆰0Do/mm1112D/mm2200Di/mm2600n/个1δ1/mm13θ/(°)100δ2/mm389256L4/mmA2/mm277851  表中:d为开孔直径;K1为系数,与A0、A1

倾角。

有关;A0为接管流通面积;A1为导流筒环隙流通面积;H为导流筒高度;do为接管外径;B为外导δ1为导流筒的公称厚度;δ2为环板的厚度;L1为管板到第1块折流板的距离;L2为管板与环板外侧的距离;L3为环板内侧到缺口的距离;L4为折流板到缺口的距离;b为内衬筒导流缺口宽度;K2为系数,与A2、A1有关;A2为壳程入口流通面流筒的宽度;Do为筒体外径;Di为导流筒内径;

2 模型建立与网格划分

2􀆰1 模型简化假设

(1)流场内流体物理性质稳定不发生相变;(3)不考虑换热管和折流板;(4)不考虑传热;

(2)流场内流体为各向同性的均匀、连续性

介质;

积;D为筒体内径;n为导流缺口数目;θ为入口

Copyright©博看网 www.bookan.com.cn. All Rights Reserved.—118—      石 油 机 械2017年 第45卷 第2期

(5)不考虑重力作用;2􀆰2 建立几何模型

(6)流场内流体运动为稳态流动。

外导流筒结构概括起来有4种[14-15]:无圆弧

1695mm,宽度462mm,缺口位于进口正后方。

离Bo=1192mm,导流筒缺口角度100°,弦长

2)、有圆弧过渡的导流筒(结构3)和夹套式变截面导流筒(结构4)。无圆弧过渡的导流筒受力情况不好,故不考虑采用;夹套式变截面导流筒会增大管板直径,增加材料费用,也不考虑采用。于是选择结构2和结构3进行对比,结构3可改为变截面式,将其作为结构5。这3种结构的几何模型分别如图2、图3和图4所示。

图4 结构5几何模型Fig􀆰4 Geometricmodelofstructure5

过渡的导流筒(结构1)、开槽式导流筒(结构

图3中外导流筒高度H=188mm,环板外壁距离Bo=1192mm,内筒高端间隙取L=118mm。

图4中外导流筒高度H=188mm,环板外壁距离Bo=1192mm,内筒高端间隙取L=118mm,入2􀆰3 网格无关性验证口倾角θ=24°。

导流筒模型结构复杂,模拟中采用非结构网

格。关于网格数量,笔者进行了网格无关性验证。

图2 结构2几何模型Fig􀆰2 Geometricmodelofstructure2

具体表述如下:对于相同模型,通过调整参数加密网格进行模拟计算,然后对比模拟结果,直至与上次结果几乎一致,则认为上次的网格数量能够满足计算要求。

选择结构5作为验证模型,控制最小单元尺寸为12mm,最大单元尺寸分别为38、30和24mm,网格总数分别为145万、267万和496万,依次增长接近1倍。对比这3种情况下壳程纵截面和入口截面速度云图,结果分别如图5和图6所示。267万的主流区速度分布差异较大,而网格总数为267万与496万的主流区速度分布差异很小,几乎一致。这说明总网格单元数为267万左右时,计算结果几乎不再变化。

由图5和图6可以看出:网格总数为145万与

图3 结构3几何模型

图2中外导流筒高度H=188mm,环板外壁距

Fig􀆰3 Geometricmodelofstructure3

图5 壳程纵截面主流区速度云图

Fig􀆰5 Velocitydistributionofthemainflowareaalongthelongitudinalsectionofshellpass

Copyright©博看网 www.bookan.com.cn. All Rights Reserved.2017年 第45卷 第2期王庆锋等:基于CFD数值模拟的换热器外导流筒结构选型      —119—

图6 壳程入口截面速度云图

267万和496万的模型中任意取出2点进行速度对比,结果如表2所示。

表2 2种网格模型任意2点速度对比Table2 2⁃pointsvelocitycontrastofthe

twokindsofgridmodelv1/(m·s-1)27􀆰657􀆰7坐标点/mQ1(0􀆰5,0􀆰5,0􀆰5)Q2(0􀆰5,1􀆰0,0􀆰3)v2/(m·s-1)27􀆰556􀆰8f/%  为了进一步体现一致性,笔者在网格总数为

Fig􀆰6 Velocitydistributionoftheinletsectionofshellpass

用压力基求解器,对于高速流动问题,常采用密度

基求解器[16]。空气动力学中,马赫数大于0􀆰3时,需考虑密度变化的影响,视为可压缩流体;反之,视为不可压缩流体[17]。据此,对模拟中的介质负压饱和蒸汽进行判断[18]。

蒸汽的气体常数:

Rg=R/Mc=

kRgT

(1)(2)(3)

0􀆰361􀆰56  当地声速:  马赫数:

v2分别为267万网格模型和496万网格模型对应速度;f为同一点不同网格模型的速度差值百分率。

由表2可看出:对于点Q1和Q2,2种网格模

表中:Q1、Q2分别为几何模型上2坐标点;v1、Ma=v/c

式中:Rg为蒸汽气体常数,J/(kg·K);R为摩kg/mol;c为当地声速,m/s;k为比热比;T为温度,K;Ma为马赫数;v为蒸汽流速,m/s。根据式(1)~式(3)进行计算,结果如表3所示。其中Ma=0􀆰174<0􀆰3,故该介质视为不可压3􀆰2 计算模型

缩流体,则选用压力基求解器。

模拟中着重于流场分析,入口雷诺数为

尔气体常数,J/(mol·K);M为蒸汽摩尔质量,

型速度差异分别为0􀆰36%和1􀆰56%,此差异在工程允许范围内。此外,笔者还针对267万左右的网格进行了稀疏探索,将网格总数稀疏到225万,发现网格总数225万与网格总数267万的结果差异很小。因此,模拟中网格总数选择267万左右。

3 Fluent设置与计算

3􀆰1 求解类型的选择

对于低速流动密度变化不大的流场计算,多选

46200,流动为湍流,不考虑相变和传热,只需设置湍流模型。选择标准κ⁃ε湍流模型,该模型具有较高的稳定性、经济性和计算精度[17]。

表3 马赫数计算结果

R/(J·mol-1·K-1)M/(kg·mol-1)0􀆰018kv/(m·s-1)75Table3 Machnumbercalculationtable1􀆰1354623􀆰3 边界条件8􀆰3145Rg/(J·kg-1·K-1)354􀆰35T/Kc/(m·s-1)4310􀆰174Ma3􀆰3􀆰1 入口边界条件设置

模拟着重关注壳程入口的流场特性,实际问题中出口充分发展,所以将自由流边界设置为出口边3􀆰3􀆰3 壁面边界条件设置3􀆰4 收敛标准界条件。

模拟中负压饱和蒸汽流动视为不可压缩流动,

所以选择速度入口作为入口边界条件。根据进口流量和管径计算得出速度为75m/s,压力-50kPa,3􀆰3􀆰2 出口边界条件设置

湍流强度4􀆰18%,水力直径0􀆰704m。

无滑移静止壁面,保持默认设置。

模拟中残差标准保持默认设置。收敛判定标准

Copyright©博看网 www.bookan.com.cn. All Rights Reserved.—120—      石 油 机 械2017年 第45卷 第2期

如下:残差曲线下降到设置标准,随步数增加稳定0􀆰001,则认为收敛。

不变,则认为收敛;进出口质量流量差值小于

3􀆰5 计算及结果分析

5几种外导流筒结构进行模拟计算。计算结果如图方法,使问题描述地更加直观清晰。

mm、外导流筒内外筒间隙H=188mm、环板外壁距离B0=1192mm不变,对结构2、结构3和结构

控制筒体内径2200mm、入口接管外径720

7和图8所示。结果分析采用定性与定量相结合的

图7 壳程纵截面速度云图

Fig􀆰7 Velocitydistributionalongthelongitudinalsectionofshellpass

3􀆰5􀆰1 定性分析

图8 壳程入口截面速度云图

Fig􀆰8 Velocitydistributionattheinletsectionoftheshellpass

图7所示。按照结构2设计,饱和蒸汽进入壳程后,沿竖直方向整体速度梯度大,局部区域速度高,流体分布均匀性差;与结构2相比,结构3整体速度梯度较好,无高速区,但局部区域梯度明5的整体速度梯度很小,无高速区,流体分布均匀性更佳。

3种导流筒结构壳程入口截面的速度分布效果显,流体分布均匀性次之;相比前2种结构,结构

3种导流筒结构壳程纵截面的速度分布效果如

范围内即可;标准差越小,表示流场均匀性越好。对模拟结果进行数理统计,如表4和表5所示。

表4为结构2、3、5壳程纵截面和入口截面的速度平均值。根据工程经验,常压下,壳程一般气体的流速范围为2~15m/s;-40kPa时,蒸汽的最大允许流速为45m/s[19-20]。由表4可知,按照结构2设计,蒸汽进入壳程后,入口截面和纵截面平均流速都接近最大允许速度,阻力消耗较大,压降较大;相比结构2,结构3纵截面平均流速较为合理,但入口截面流速依然较大;而结构5壳程入口截面和纵截面平均流速均相对合理。

表4 速度平均值

 分析截面壳程入口截面壳程纵截面如图8所示。按照结构2设计,饱和蒸汽进入壳程后,在入口截面上整体速度梯度大,进口区域速度较高,流体分布均匀性差;与结构2相比,结构3整体速度梯度较好,无高速区,局部区域梯度明5的整体速度梯度均匀,无高速区,流体分布3􀆰5􀆰2 定量分析

均匀性更佳。

显,流体分布均匀性次之;相比前两种结构,结构

m·s-1

结构522􀆰2210􀆰02Table4 Theaveragevelocity结构242􀆰5745􀆰55结构331􀆰1918􀆰08m·s-1

表5为结构2、3、5壳程纵截面和入口截面的2标准差较大,流场均匀性较差,结构3和结构5速度标准差。由表可知,对于壳程入口截面,结构标准差接近,均远小于结构2,流场均匀性较好;

采用速度平均值p和速度标准差S来定量地描

述壳程蒸汽速度和流场均匀性。其中,速度在合理

Copyright©博看网 www.bookan.com.cn. All Rights Reserved.2017年 第45卷 第2期王庆锋等:基于CFD数值模拟的换热器外导流筒结构选型      —121—

对于壳程纵截面,结构2标准差较大,流场均匀性5标准差更小,流场均匀性更好。

较差,结构3标准差较小,流场均匀性较好,结构

表5 速度标准差

结构227􀆰9433􀆰02结构315􀆰4711􀆰99结构516􀆰928􀆰38[7] 古新,朱培纳,刘敏珊,等.变截面导流筒换热器[8] 夏立荣,岳希明,刘敏珊.换热器夹套式变截面导

流筒数值模拟与结构优化[C]∥全国换热器学术会议.2007:46-50.

[9] 王庆锋,庞鑫,赵双.管壳式换热器传热效率影响[10] 吴金星,魏新利,董其伍,等.纵向流管壳式换热[11] 刘敏珊,孙爱芳,董其伍.紧凑型板壳式换热器导[12] 李平,蒋丹,何克强.高压加热器蒸汽冷却段导流[13] 中华人民共和国国家质量监督检验检疫总局,中

国国家标准化管理委员会.压力容器:GB150—2011[S].北京:中国标准出版社,2011.

[14] 毛希澜.化工设备设计全书之换热器设计[M].

(5):762-765.32-35.(4):49-52.(10):102-107.63(12):3839-3846.

进[J].石油化工设备技术,2007,28(4):1-3.

入口流场均化性能数值仿真[J].化工学报,2012,

 分析截面Table5 Thestandarddeviationofthevelocity壳程入口截面壳程纵截面因素及数值模拟分析[J].石油机械,2015,43

4 结 论

大,流场均匀性差;与结构2相比,结构3整体速度梯度较好,但局部区域梯度明显,流场均匀性较好;相比前2种结构,结构5的整体速度梯度更均匀,流体分布均匀性更佳。

(2)结构2入口截面和纵截面平均流速较大;(1)结构2进口区域速度高,壳程速度梯度

器强化传热研究与发展[J].石油机械,2002,30

流结构优化设计[J].压力容器,2006,23(8):

筒内流动的数值模拟[J].动力工程,2007,27

相比结构2,结构3纵截面平均流速较小,但入口截面平均流速依然较大;相比前2种结构,结构5壳程入口截面和纵截面平均流速均较小且更合理。

(3)相比结构3和结构5,结构2标准差较

大,流场均匀性较差;结构3和结构5入口截面标准差接近,纵截面标准差结构5小于结构3,故结构5流场均匀性更好。

(4)综合来看,结构5流速合理,流场均匀

参 考 文 献

[1] 兰州石油机械研究所.换热器[M].2版.北京:[2] 张哲,厉彦忠,田津津.板翅式换热器导流片结构[3] 焦安军,厉彦忠,张瑞,等.板翅式换热器不同结[4] 杜明,陈建文,惠宁,等.管壳式换热器壳侧外导[5] 张立文.换热器环形外导流筒的设计及制造要点[6] 郭雪华,段瑞.换热器外导流筒的结构设计分析与改

[J].石油化工设备技术,2005,26(6):15-16.2014(10):23-26.2001,35(11):1113-1117.-1314.

中国石化出版社,2013.

[15] 史美中,王中铮.热交换器原理与设计[M].4[16] 胡坤,李振北.ANSYSICEMCFD工程实例详解[17] 于勇.FLUENT入门与进阶教程[M].北京:北京[18] 沈维道,童均耕.工程热力学[M].4版.北京:[19] 朱聘冠.换热器原理及计算[M].北京:清华大[20] 孙兰义,马占华,王志刚,等.换热器工艺设计

[M].北京:中国石化出版社,2015:118.学出版社,1987:253-254.高等教育出版社,2012:76.理工大学出版社,2008:40-44.

[M].北京:人民邮电出版社,2014:385.版.南京:东南大学出版社,2009:52-53.

上海:上海科学技术出版社,1988:185-186.

性更好,为笔者研究问题更适宜的导流筒结构。

的数值模拟[J].化工学报,2002,53(12):1311

构导流片导流性能的研究[J].西安交通大学学报,

  第一作者简介:王庆锋,副教授,生于1972年,2011年毕业于北京化工大学化工机械专业,获博士学位,研究方向为设备设计与风险评估管理。地址:(100029)北京市朝阳区。

通信作者:郝 帅,E⁃mail:1010489131@qq.com。

 

收稿日期:2016-08-21

(本文编辑 刘 峰)

流筒结构设计与工程应用[J].石油和化工设备,

Copyright©博看网 www.bookan.com.cn. All Rights Reserved.

因篇幅问题不能全部显示,请点此查看更多更全内容

Top