脉动热管换热器传热性能的研究.pdf

脉动热管换热器传热性能的研究倪摇 兵1,2,郝玉刚1,吴秀杰1,徐建炎1(1.中国恩菲工程技术有限公司,北京 100038;2.大连理工大学,辽宁 大连 116024)摘摇 要 通过脉动热管的理论基础及实验,研究了脉动热管的启动特性、准稳态运行特性以及换热工质为含不凝汽的蒸汽时,脉动热管换热器温度的变化规律。同时以传统重力热管换热器的计算为依据,结合脉动热管自身特点,针对工程常见的烟气进行设计计算。通过实验数据分析,确定了脉动热管换热器在实验条件下,设备响应较快,启动性能较好。通过设计计算,提供了脉动热管换热器的结构尺寸。关键词 脉动热管换热器;运行特性;余热回收;设计计算中图分类号 TK121摇 摇 摇 文献标志码 B摇 摇 摇 文章编号 1003-8884(2023)04-0041-07DOI:10.19611/11鄄2919/tg.2023.04.008收稿日期 2023-05-16第一作者 倪兵(1984),吉林长春人,博士,主要从事余热回收及废水处理相关的设备研发工作。引用格式 倪兵,郝玉刚,吴秀杰,等.脉动热管换热器传热性能的研究J.有色设备,2023,37(4):41-47.0摇 引言余热利用对节约能源、增加产量、提高产品质量和降低生产成本等方面起着越来越大的作用,已成为工业生产中不可分割的组成部分。但由于生产方法、生产工艺、生产设备以及原料、条件的不同和工艺上千变万化的需要,给余热回收利用带来很多困难。脉动热管换热器作为一种先进的换热装置,已被应用于电子冷却、空调余热回收、太阳能集热器等领域,并且可实现对废水及废气进行余热回收,是一种高效相变传热元件。在相同温差下,管内流体以汽化潜热的方式传递热量,并且脉动热管适用于不同加热方式和加热位置。将脉动热管技术应用于余热回收水处理装置中,不仅能降低废水及废气的温度,还能将热量回收利用,达到节能减排的目的1-2。20 世纪 40 年代初,美国 Gaugler 就已经提出热管这个理论。1963 年第一个热管问世,至此热管才进入现实生活中,并不断受到人们的重视,成为一种具有高传热效率的元件。热管的理论阶段一直持续到 20 世纪 70 年代后才进入应用阶段,但仍然存在技术不成熟、成本较高的问题。所以只有航天、核能等一些高端领域范围得以应用3。到80 年代后,由于技术不断成熟,使成本降低,热管技术才逐渐进入不同的生产领域。影响脉动热管传热性能的因素有很多,但其工作原理无法很明确的被解释,运行机理非常复杂,涉及到的因素较多,带有很大的随机性。自从脉动热管的概念提出以来,国内外研究者进行较多的研究。已有文献显示,研究人员研究了各种工作液体的脉动热管4-8,从他们的研究结果中可知,直径并不是影响脉动热管传热性能的唯一参数。脉动热管的管长、弯头数、每段长度(加热段、绝热段和冷凝段)、工作液体的热物理特性、加热/冷却的范围和模式9-13等都会影响脉动热管的传热性能。在余热回收领域中,Rittidech 等14将脉动热管应用到空气预热器中,用于干燥尾气回收,通过实验可知,当热空气温度从 60 益升高到 80 益,热传递效率逐渐上升;工作液体由水换成 R123,热传递速率升高。Khandekar15研究了脉动热管换热器在余热处理中的应用。通过实验显示,脉动热管换热器可以在较高热流密度条件下应用,但在实际应用中,需要更系统和详细的研究来建立该系统。Meena等16针对干燥循环系统设计了闭合环路脉动热管空气预热器,用于降低尾气湿度并回收其热量。脉动热管的内径为 2 mm,蒸发及冷凝段长为 0郾 19 m,14有色设备摇 2023 年第 4 期绝热段长为 0郾 08 m,总长度为 3郾 58 m。实验过程中,热空气的流速为 0郾 5 3 m/s,温度为50 80 益,相对湿度为 100%,热 管 内 填 充 的 工 作 液 体 为R134a,充液率为 50%。实验结果显示,随热空气温度上升,传热速率增大,相对湿度由 89%100%降至 54%72%。经过对脉动热管传热特性的不断研究,我国脉动热管技术已经应用到加热炉余热再生利用系统 辽河油田 10 吨注汽锅炉烟气余热回收系统,该系统用于加热注汽炉的进水,可以将注汽炉进水温度提高 20 益 左右,排出烟气温度降低100 益,燃料节约 5%,在采用脉动热管进行余热回收后,每台注汽炉每年节约燃料费用可达到 30 万元以上2。本文对脉动热管换热器进行运行特性的研究和设计计算。结果可应用于余热回收及水处理装置中,为化工、石化、电力、冶金等重点耗能行业推广提供设计方法。1摇 实验说明1郾 1摇 实验系统根据实验原理建立脉动热管测试实验系统,其示意图如图 1 所示。实验系统由脉动热管、蒸汽发生器、空气加热器、空气压缩机、数据采集仪、流量计、阀门、管道及管件等组成,实验图片如图 2 所示。实验冷热水箱均选用聚丙烯材质,壁厚为 10 mm,水箱尺寸长 伊 宽 伊 高为 800 mm 伊80 mm 伊300 mm。图 1摇 实验系统示意图摇图 2摇 实验系统图片摇摇 摇 脉动热管换热器内部安装一组铜水闭合回路脉动热管。本实验中选用铜作为管壳材料。管内封装的工作液体为去离子水,充液率为 50%,绝对压力为1 伊10-2Pa。管外径为 4 mm,为保证强度,管壁厚0郾 5 mm。脉动热管分为加热段、绝热段和冷凝段三部分,其长度分别为250 mm、200 mm 和150 mm。系统安装了蒸汽发生器,为了得到预期的热风进口温度,在热风管的空气入口前,安装了电加热器。两个气体流量计分别连接空气和蒸汽管上,用以监测其流量。蒸汽与经空气加热器加热过的空气混合后进入热管换热器的热段,与冷段的冷却水通过脉动热管进行热交换,冷却水的流量由转子流量计进行测量,混合气体中不凝汽的含量可通过调节蒸汽和空气的流量来控制。1郾 2摇 实验方法实验前确定脉动热管热负荷范围,在负荷范围内实验,确定启动功率和启动时间。在不同输入功率时,观察壁面温度曲线的跳跃现象,当壁面温度发生急剧跳跃现象时,系统启动并监测对应启动温度、24有色设备摇 2023 年第 4 期启动功率、启动时间等。启动后,运行约 10 min,使热管达到稳定运行状态,稳定后进行实验数据监测,时间约 10 30 min。开启蒸汽发生器,通过蒸汽流量计观察流量大小,调节蒸汽发生器出口阀门,使流量达到所需大小。运行稳定后开始测量各测点的温度和压力。启用空气压缩机和加热器,通过空气流量计观察流量大小,调节空气压缩机使流量达到所需大小。经加热器加热后的热空气和蒸汽混合后进入换热器热段。1郾 3摇 结果与讨论对于实际工程项目应用过程,目前主要通过启动性能和传热性能两方面来对脉动热管性能进行评价。通常将脉动热管从开始加热到进入稳定振荡运行定义为启动过程,反映的是热管对工作环境温度变化做出反应的过程。传热性能作为热管在一定条件下传热能力的体现,是热管制作研究的重要目标。当脉动热管壁面温度曲线第一次出现波峰、波谷振荡就说明脉动热管开始启动。但在低热负荷下,有时工质开始流动,但并未发生明显相变过程,这样在温度曲线上不会出现波峰和波谷的变化。因此,在温度曲线光滑上升阶段,如果 dt/d子 的值出现明显改变时,同样可以认为脉动热管开始启动。脉动热管进入稳态运行前的阶段包括脉动热管的启动、间歇震荡和准稳态震荡阶段。在脉动热管启动后,脉动热管要经历升温过程,所以在进入稳态运行后,不能完全避免间歇震荡。图 3摇 脉动热管壁面温度变化趋势在分析启动特性和准稳态运行特性时,选择启动时温度跳跃幅度相对较大,稳定性相对较差的实验组数据进行分析。图 3 是冷水流量为 80 L/h,蒸汽流量为 4郾 4 m3/h 时,脉动热管壁面温度的变化趋势。当系统运行后,脉动热管壁面温度随时间 子 的增长而升高。当 子1达到 40 s 后,脉动热管壁面温度发生震荡。在此之前,脉动热管并未启动,内部工质处于静止状态,对比图中相同时间,热段、绝热段和冷凝段的壁面温度,温度沿脉动热管的轴向略有升高,这是通过铜管导热作用使热量进行传递。当时间达到 40 s 后,蒸发段管内液体工质吸收热量发生相变,产生气泡,气泡被输送到冷凝段被冷却,释放汽化潜热,此时脉动热管启动。受压差作用,脉动热管内部工质发生剧烈震荡,从而造成启动时温度跳跃。间歇性振荡是脉动热管启动后出现的反复“停止-启动冶的现象。当 子1分别达到 78 s、218 s 和 410s 后,脉动热管震荡停止,温度呈平滑变化后再次震荡。随着实验的进行,脉动热管内部工质不断积累热量,直到工质能连续震荡。间歇震荡会影响脉动热管传热过程的稳定性。当脉动热管倾角为 90 益或热负荷较大时会减少间歇震荡,本组实验输入热负荷较小,所以实验过程出现间歇震荡。在实际工程应用中热负荷相对较大,因此,可以减少间歇震荡对传热过程稳定性的影响。当时间达到 478 s 后,脉动热管进入准稳态震荡阶段。在此条件下,脉动热管的启动时间为 40 s,启动时温度跳跃幅度(启动时最高温度和最低温度的差值)为 4 益,跳跃幅度较小,间歇震荡时间为 7min,在此阶段出现了温度急剧跳跃的现象,冷凝段壁面温差达到 30 益。但在此阶段停留时间较短,所以设备响应较快,启动性能较好。图 4摇 不同含气率时流体和壁面的平均温度变化趋势在热源工质中含有硫化物时,为防止低温腐蚀,将工质的温度控制高于露点 20 益左右,因此换热器壁面温度和流体平均温度的温差直接影响了换热器34倪摇 兵等:脉动热管换热器传热性能的研究的应用。图 4 为不同含气率时流体和壁面温度的变化趋势。随着含气率的增大,壁面平均温度和流体平均温度的温差逐渐升高。当含气率高于 20%后壁面平均温度和流体平均温度的温差高于 20 益,当含气率高于 40%后其温差趋于稳定,为 26 益。2摇 脉动热管设计计算2郾 1摇 脉动热管的设计原理以冷源侧水和热源侧烟气换热为例。1)换热量计算热端放出热量 Qa计算式见式(1)。Qa=籽1q1Cv1(t忆1-t1)(1)冷端吸收热量 Qb计算式见式(2)。Qb=(1+灼)Qa(2)定性温度计算式见式(3)。tm1=t1+t忆12(3)式中:籽1和 Cv1分别为热端定性温度 tm1时,烟气的密度(kg/m3)和比热(kJ/(kg K);q1为烟气流量,kg/s;t1和 t忆1分别为烟气进出口温度,益;灼 为散热器散热损失,取 5%10%2)热端换热系数计算烟气侧雷诺数 Re1计算式见式(4)。Re1=u1d0淄1(4)烟气侧努谢尔特数 Nu1计算式见式(5)。Nu1=CRem1(5)对流换热系数 h1计算式见式(6)。h1=Nu姿1d0(6)式中:淄1和 姿1分别为热端定性温度 tm1时,烟气的黏度(m2/s)和导热系数(W/(mK);d0为脉动热管的直径,m;u1为烟气侧的迎风速度,m/s。3)冷端换热系数计算冷端定性温度计算式见式(7)。tm2=t2+t忆22(7)水侧雷诺数 Re2计算式见式(8)。Re2=u2d0淄2(8)烟气侧努谢尔特数 Nu2计算式见式(9)。Nu2=CRem2(9)对流换热系数 h2计算式见式(10)。h2=Nu姿2d0(10)式中:t2和 t忆2分别为水侧进出口温度,益;淄2和 姿2分别为热端定性温度 tm2时,水的黏度(m2/s)和导热系数(W/(m K);u2为水侧的迎风速度,m/s。4)换热系数计算管外流体与管外壁的换热量 Q1计算式见式(11)。Q1=h0A(tf-two)(11)管外壁到管内壁的到热量 Q2计算式见式(12)。Q2=2仔姿wl(two-twi)ln(d0/di)(12)管内壁与管内工质的换热量 Q3计算式见式(13)(15)。Q3=hiA(tw,i-tn)(13)Q1=Q2=Q3=KA(tf-tn)(14)K=11h0+d0ln(d0/di)2姿w+1hi(15)式中:h0和 hi分别为管内、外流体的换热系数,W/(m2 K);姿w为脉动热管导热系数,W/(mK);di为脉动热管内径,m。管内流体的换热系数 hi较大,所以忽略1/hi,传热公式简化为式(16)。K=11h0+d0ln(d0/di)2姿w(16)5)换热温差计算脉动热管壁面的温度 tw计算式见式(17)。tw=h0A1tm1+h1A2tm2h0A1+h1A2(17)6)换热面积计算根据公式 Q=KA驻t 计算换热面积。热端换热面积 A1计算式见式(18)。A1=Q1K1驻t1(18)冷端换热面积 A2计算式见式(19)。A2=Q2K2驻t2(19)总换热面积 A 计算式见式(20)。A=A1+A2(20)7)脉动热管结构计算44有色设备摇 2023 年第 4 期热端迎风面积 Ay1计算式见式(21)。Ay1=q1u1(21)冷端迎风面积 Ay2计算式见式(22)。Ay2=q2u2(22)式中:ql和 q2分为冷热端流体的流量 kg/s。迎风面宽 E 计算式见式(23)。E=AL(23)迎风面上管排数 B 计算式见式(24)。B=EPl(24)式中:L 为热管长度;Pl为脉动热管迎风面上管间距。8)脉动热管尺寸计算0郾 7滓g(籽l-籽g)臆D臆2滓g(籽l-籽g)(25)2郾 2摇 脉动热管的设计结果2郾 2郾 1摇 冷源无相变本案例工质:冷源侧为水,热源侧为烟气,初始参数如表 1 所示。表 1摇 初始参数参数热端冷端工质烟气水入口温度200/益60/益出口温度120/益80/益流量65 000/Nm3/h 摇 摇 脉动热管的初始设计参数如表 2 所示,脉动热管换热器的初始计算参数如表 3 所示,脉动热管设计计算结果如表 4 所示。根据计算,管脉动热管换热器尺寸长 伊 宽 伊 高为 1郾 73 m 伊0郾 93 m 伊1郾 7 mm。表 2摇 脉动热管初始参数参数直径充液率壁厚Pl/d0迎风距离工作流体值4/mm50%0郾 6/mm28/mm去离子水表 3摇 脉动热管换热器初始参数参数 热端长 冷端长 绝热段长热端迎风风速冷端迎风风速值1 m0郾 2 m0郾 2 m6 m/s0郾 5 m/s表 4摇 脉动热管换热器设计计算参数器件参数蒸发端冷凝端平均温差73郾 7/益11郾 2/益换热器迎风面积2郾 25/m20郾 03/m2换热系数 216/W/(m2 K)5 670/W/(m2 K)脉动热管总长23 768/m5 645/m面积309/m273郾 5/m22郾 2郾 2摇 冷源有相变在低温多效蒸发中,低温蒸汽可以直接做为热源。因此,本案例将脉动热管换热器冷源室设置于真空条件下,按冷端出口为蒸汽进行设计计算。工质:冷源为水,热源为烟气,初始参数如表 5所示。表 5摇 初始参数参数热端冷端工质烟气水/汽入口温度200/益60/益出口温度120/益80/益压力0郾 1/MPa0郾 047/MPa流量65 000/Nm3/h 摇 摇 脉动热管的初始设计参数如表 6 所示,脉动热管换热器的初始计算参数如表 7 所示,脉动热管设计计算结果如表 8 所示。根据计算,管脉动热管换热器尺寸长 伊 宽 伊 高为 1郾 73 m 伊0郾 58 m 伊1郾 7 m。表 6摇 脉动热管初始参数参数直径充液率壁厚Pl/d0迎风距离工作流体值4/mm50%0郾 6/mm28/mm去离子水表 7摇 脉动热管换热器初始参数参数 热端长 冷端长 绝热段长热端迎风风速冷端迎风风速值1/m0郾 2/m0郾 2/m6/m/s0郾 5/m/s3摇 结论本文建立了脉动热管测试实验台,分析了脉动热管的启动特性和准稳态运行特性,并分析了54倪摇 兵等:脉动热管换热器传热性能的研究摇 摇表 8摇 脉动热管换热器设计计算参数器件参数蒸发端冷凝端平均温差73郾 5/益11郾 5/益换热器迎风面积3/m20郾 01/m2换热系数 242/W/(m2 K)9 353/W/(m2 K)脉动热管总长21 254/m58/m面积277/m20郾 8/m2在换热工质为水和蒸汽含不凝汽时脉动热管换热器的温度。根据脉动热管的理论基础,通过已有文献的数据及相关项目的参数,对脉动热管换热器进行了结构尺寸的设计计算,主要结论如下:1)在实验条件下,脉动热管在启动过程中可能会产生较大的跳跃温度,间歇震荡明显,但间歇震荡时间较短,设备响应较快,启动性能较好,在实际工程中选用脉动热管换热器回收余热,流体平均温度和壁面平均温度的温差较小,有效减小腐蚀对换热设备的影响。2)脉动热管换热器在设计计算条件下,烟气侧的换热系数最高可达 242 W/(m2K),相比传统换热设备,脉动热管换热器的换热系数可以达到普通换热设备的 2 倍以上。3)以传统重力热管为基础,针对热源侧为烟气的工程实例设计脉动热管换热器,计算了脉动热管及换热器内部结构尺寸。参考文献1 柴本银.振荡流热管换热器换热性能及其在干燥系统中的应用研究D.济南:山东大学,2010.2 韩雅芳.脉动热管及其换热器传热特性研究D.北京:北京工业大学,2013.3 马永昌,张宪峰.热管技术的原理、应用与发展J.变频器世界,2009(7):70-75.4 Burban G,Ayel V,Alexandre A,et al.Experimental in鄄vestigation of a pulsating heat pipe for hybrid vehicle appli鄄cationsJ.Applied Thermal Engineering,2013,50(1):94-103.5 Liu S,Li J,Dong X,et al.Experimental study of flow pat鄄terns and improved configurations for pulsating heat pipesJ.Journal of Thermal Science,2007,16(1):56-62.6 Moghari R M,Mujumdar A S,Shariat M,et al.Investiga鄄tion effect of nanoparticle mean diameter on mixed convec鄄tion Al2O3鄄water nanofluid flow in an annulus by two phasemixture modelJ.International Communications in Heatand Mass Transfer,2013,49:25-35.7 Saha M,Feroz C M,Ahmed F,et al.Thermal perform鄄ance of an open loop closed end pulsating heat pipeJ.Heat and Mass Transfer,2012,48(2):259-265.8 Tseng C Y,Yang K S,Chien K H,et al.Investigation ofthe performance of pulsating heat pipe subject to uniform/alternating tube diametersJ.Experimental Thermal andFluid Science,2014,54:85-92.9 Yang H,Jiang S,Shen W,et al.Certificateless ProvableGroup Shared Data Possession with Comprehensive PrivacyPreservation for Cloud StorageJ.Future Internet,2018,10(6):49-60.10 Nikolayev V S.Effect of tube heat conduction on the sin鄄gle branch pulsating heat pipe start鄄upJ.InternationalJournal of Heat and Mass Transfer,2016,95:477-487.11 Pachghare P R,Mahalle A M.Thermo鄄hydrodynamics ofclosed loop pulsating heat pipe:an experimental studyJ.Journal of Mechanical Science and Technology,2014,28(8):3387-3394.12 Charoensawan P,Terdtoon P.Thermal performance ofhorizontal closed鄄loop oscillating heat pipesJ.AppliedThermal Engineering,2008,28(5):460-466.13 Pastukhov V G,Maydanik Y F.Development of a pulsa鄄ting heat pipe with a directional circulation of a workingfluidJ.Applied Thermal Engineering,2016,109:155-161.14 Rittidech S T P,Tantakon P,et al.effect of inclinationangles,evaporator section lengths and working fluid prop鄄erties on heat transfer characteristics of a closed鄄end oscil鄄lation heat pipeC.6th Internationl Heat Pipe Symposi鄄um,Chinag Mai,Thailand,2000.15 S.Khandekar.Pulsating heat pipe based heat exchangersC.The 21st International Symposium on TransportPhenomena,Kaohsiung City,Taiwan,2010.16 Meena P,Rittidech S,Poomsa鄄Ad N.Closed鄄loop oscil鄄lating heat鄄pipe with check valves(CLOHP/CVs)air鄄preheater for reducing relative humidity in drying systemsJ.Applied Energy,2007,84(4):363-373.64有色设备摇 2023 年第 4 期Study on Heat Transfer Performance of Pulsating Heat Pipe Heat ExchangerNI Bing,HAO Yu鄄gang,WU Xiu鄄jie,XU Jian鄄yanAbstract:The start鄄up characteristics,quasi鄄steady operation characteristics,and the temperature varia鄄tion regular of pulsating heat pipe heat exchanger when working fluid containing non condensable steam insteam were studied and experimented on.The pulsating heat pipe heat exchanger were designed and cal鄄culated using flue gas as the working fluid.Pulsating heat pipe heat exchangers with flue gas as workingfluid were design and calculation,based on the calculation of traditional gravity heat pipe heat exchangersand combined with the characteristics of pulsating heat pipes themselves.Experimental data by analysis,it is found that under experimental conditions,the pulsating heat pipe heat exchanger has a fast responseand good start鄄up performance.The design dimensions of pulsating heat pipes are provided by design cal鄄culations.Key words:pulsating heat pipe heat exchanger;starting characteristics;design calculation蕚蕚蕚蕚蕚蕚蕚蕚蕚蕚蕚蕚蕚蕚蕚蕚蕚蕚蕚蕚蕚蕚蕚蕚蕚蕚蕚蕚蕚蕚蕚蕚蕚蕚蕚蕚蕚蕚蕚蕚蕚蕚蕚蕚蕚蕚蕚蒉(上接第 19 页)Automated Repair Solution for the Full Process ofAnode Guide Rod AssemblyWU Dao鄄lin,LI Meng,WU Zu鄄huai,ZHANG Ye,ZHOU BingAbstract:The anode guide rod component is one of the important components in the production of elec鄄trolytic aluminum.At present,the anode guide rod components are mostly repaired manually.This meth鄄od has problems such as high labor intensity,low repair efficiency,and unstable quality,which can leadto detachment,cracking,or reduced lifespan of the anode guide rod components during the productionstage.This article introduces a full process automated repair solution for anode guide rod assembly.Thisplan covers narrow gap full section welding process and equipment,anode steel claw repair technologyand equipment based on resistance thermal pressure welding,and automation technology and equipmentfor the entire process of aluminum guide rod repair.The application of this solution can repair more than40 sets of components per day,extend the service life by more than half,reduce the labor intensity of op鄄erators,and save more than 15 kW h of electricity per ton of aluminum in the electrolytic aluminum in鄄dustry.Key words:electrolytic aluminum;anode guide rod;narrow gap full section welding;resistance thermalpressure welding蒉74倪摇 兵等:脉动热管换热器传热性能的研究