石 油 机 械
CHINAPETROLEUMMACHINERY
2017年 第45卷 第2期
◀炼油与化工机械▶
∗
基于CFD数值模拟的换热器外导流筒结构选型
王庆锋 郝 帅 李 凯 李 中
(北京化工大学机电工程学院)
摘要:管壳式换热器壳程流场均匀分布能够改善传热效果,外导流筒的合理选型对流场均匀性有重要影响。然而,外导流筒结构多样,其选型设计多依靠经验,缺乏指导性方法。鉴于此,以筒体直径2200mm、入口接管外径720mm的竖直管降膜式蒸发器换热段为研究对象,设计了外导流筒,建立了3种结构的三维实体模型。基于CFD数值模拟的方法,以平均速度和流场均匀性为衡量指标,研究了3种导流筒结构壳程纵截面和入口截面的流体分布效果,证明了外导流筒合理选型确实能够改善换热器壳程流场均匀性。定性结合定量分析得出:变截面圆弧过渡式外导流筒平均速度相对合理,流场均匀性最佳,是较适宜的导流筒结构。研究结果对外导流筒的选型设计具有重要意义。
关键词:管壳式换热器;外导流筒;CFD;数值模拟;选型
中图分类号:TE965 文献标识码:A doi:1016082/jcnkiissn1001-4578201702026
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⁃nalflowdistributorlection
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 开槽式外导流筒结构
Fig1 Structureoftheslottedtypeexternalflowdistributor
表1 外导流筒基本尺寸
d/mm70426K1L1/mm109811A0/mm2L2/mm150A1/mm2L3/mm3242818H/mm188428do/mmb/mm462720H/doK2Table1 GeneraldimensionsoftheexternalflowdistributorB/mm113902620Do/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 模型建立与网格划分
21 模型简化假设
(1)流场内流体物理性质稳定不发生相变;(3)不考虑换热管和折流板;(4)不考虑传热;
(2)流场内流体为各向同性的均匀、连续性
介质;
积;D为筒体内径;n为导流缺口数目;θ为入口
Copyright©博看网 www.bookan.com.cn. All Rights Reserved.—118— 石 油 机 械2017年 第45卷 第2期
(5)不考虑重力作用;22 建立几何模型
(6)流场内流体运动为稳态流动。
外导流筒结构概括起来有4种[14-15]:无圆弧
1695mm,宽度462mm,缺口位于进口正后方。
离Bo=1192mm,导流筒缺口角度100°,弦长
2)、有圆弧过渡的导流筒(结构3)和夹套式变截面导流筒(结构4)。无圆弧过渡的导流筒受力情况不好,故不考虑采用;夹套式变截面导流筒会增大管板直径,增加材料费用,也不考虑采用。于是选择结构2和结构3进行对比,结构3可改为变截面式,将其作为结构5。这3种结构的几何模型分别如图2、图3和图4所示。
图4 结构5几何模型Fig4 Geometricmodelofstructure5
过渡的导流筒(结构1)、开槽式导流筒(结构
图3中外导流筒高度H=188mm,环板外壁距离Bo=1192mm,内筒高端间隙取L=118mm。
图4中外导流筒高度H=188mm,环板外壁距离Bo=1192mm,内筒高端间隙取L=118mm,入23 网格无关性验证口倾角θ=24°。
导流筒模型结构复杂,模拟中采用非结构网
格。关于网格数量,笔者进行了网格无关性验证。
图2 结构2几何模型Fig2 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,环板外壁距
Fig3 Geometricmodelofstructure3
图5 壳程纵截面主流区速度云图
Fig5 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)276577坐标点/mQ1(05,05,05)Q2(05,10,03)v2/(m·s-1)275568f/% 为了进一步体现一致性,笔者在网格总数为
Fig6 Velocitydistributionoftheinletsectionofshellpass
用压力基求解器,对于高速流动问题,常采用密度
基求解器[16]。空气动力学中,马赫数大于03时,需考虑密度变化的影响,视为可压缩流体;反之,视为不可压缩流体[17]。据此,对模拟中的介质负压饱和蒸汽进行判断[18]。
蒸汽的气体常数:
Rg=R/Mc=
kRgT
(1)(2)(3)
036156 当地声速: 马赫数:
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=0174<03,故该介质视为不可压32 计算模型
缩流体,则选用压力基求解器。
模拟中着重于流场分析,入口雷诺数为
尔气体常数,J/(mol·K);M为蒸汽摩尔质量,
型速度差异分别为036%和156%,此差异在工程允许范围内。此外,笔者还针对267万左右的网格进行了稀疏探索,将网格总数稀疏到225万,发现网格总数225万与网格总数267万的结果差异很小。因此,模拟中网格总数选择267万左右。
3 Fluent设置与计算
31 求解类型的选择
对于低速流动密度变化不大的流场计算,多选
46200,流动为湍流,不考虑相变和传热,只需设置湍流模型。选择标准κ⁃ε湍流模型,该模型具有较高的稳定性、经济性和计算精度[17]。
表3 马赫数计算结果
R/(J·mol-1·K-1)M/(kg·mol-1)0018kv/(m·s-1)75Table3 Machnumbercalculationtable113546233 边界条件83145Rg/(J·kg-1·K-1)35435T/Kc/(m·s-1)4310174Ma331 入口边界条件设置
模拟着重关注壳程入口的流场特性,实际问题中出口充分发展,所以将自由流边界设置为出口边333 壁面边界条件设置34 收敛标准界条件。
模拟中负压饱和蒸汽流动视为不可压缩流动,
所以选择速度入口作为入口边界条件。根据进口流量和管径计算得出速度为75m/s,压力-50kPa,332 出口边界条件设置
湍流强度418%,水力直径0704m。
无滑移静止壁面,保持默认设置。
模拟中残差标准保持默认设置。收敛判定标准
Copyright©博看网 www.bookan.com.cn. All Rights Reserved.—120— 石 油 机 械2017年 第45卷 第2期
如下:残差曲线下降到设置标准,随步数增加稳定0001,则认为收敛。
不变,则认为收敛;进出口质量流量差值小于
35 计算及结果分析
5几种外导流筒结构进行模拟计算。计算结果如图方法,使问题描述地更加直观清晰。
mm、外导流筒内外筒间隙H=188mm、环板外壁距离B0=1192mm不变,对结构2、结构3和结构
控制筒体内径2200mm、入口接管外径720
7和图8所示。结果分析采用定性与定量相结合的
图7 壳程纵截面速度云图
Fig7 Velocitydistributionalongthelongitudinalsectionofshellpass
351 定性分析
图8 壳程入口截面速度云图
Fig8 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的整体速度梯度均匀,无高速区,流体分布352 定量分析
均匀性更佳。
显,流体分布均匀性次之;相比前两种结构,结构
m·s-1
结构522221002Table4 Theaveragevelocity结构242574555结构331191808m·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 速度标准差
结构227943302结构315471199结构51692838[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.
因篇幅问题不能全部显示,请点此查看更多更全内容