«上一篇
文章快速检索     高级检索
下一篇»
  哈尔滨工程大学学报  2018, Vol. 39 Issue (5): 870-875  DOI: 10.11990/jheu.201701006
0

引用本文  

陈凯伦, 阎昌琪, 孟兆明, 等. 套管式换热元件内流动不稳定性研究[J]. 哈尔滨工程大学学报, 2018, 39(5): 870-875. DOI: 10.11990/jheu.201701006.
CHEN Kailun, YAN Changqi, MENG Zhaoming, et al. Flow instabilities in a thimble-based heat transfer element[J]. Journal of Harbin Engineering University, 2018, 39(5): 870-875. DOI: 10.11990/jheu.201701006.

基金项目

国家自然科学基金项目(11475048)

通信作者

阎昌琪, E-mail:changqi_yan@163.com

作者简介

陈凯伦(1991-), 男, 博士研究生; 阎昌琪(1955-), 男, 教授, 博士生导师

文章历史

收稿日期:2017-01-03
网络出版日期:2018-03-29
套管式换热元件内流动不稳定性研究
陈凯伦, 阎昌琪, 孟兆明, 丁涛    
哈尔滨工程大学 核安全与仿真技术国防重点学科实验室, 黑龙江 哈尔滨 150001
摘要:针对套管式换热元件内的自然循环流动不稳定性,以单管实验装置为对象进行了不同加热功率和入口过冷度下的实验研究。分析了低压条件下系统内可能出现的不稳定现象及其形成机理,确定了流动不稳定的发生区域,并得到了用无量纲过冷度和加热功率表示的间歇泉和闪蒸不稳定边界。结果表明:在实验加热功率范围内,随入口过冷度减小,系统内相继发生间歇泉和闪蒸不稳定流动。与一般自然循环回路相比,套管式换热元件的特殊性在于其上升段流体向下降段流体传热,因此上升段流体温度降低,汽泡的冷凝现象更严重,更容易出现间歇泉不稳定,并在一定程度上推迟闪蒸不稳定的发生。闪蒸发生时流量振荡周期与流体流经上升段所需时间的比值范围为0.8~1.1。
关键词套管式换热元件    自然循环    流动不稳定性    间歇泉振荡    闪蒸不稳定    不稳定边界    
Flow instabilities in a thimble-based heat transfer element
CHEN Kailun, YAN Changqi, MENG Zhaoming, DING Tao    
Fundamental Science on Nuclear Safety and Simulation Technology Laboratory, Harbin Engineering University, Harbin 150001, China
Abstract: To understand natural circulation flow instabilities in a thimble-based heat transfer element, an experimental study was conducted on a monotube experimentation unit using different heating powers and inlet subcooling indices to analyze the unstable flow phenomena that may occur at low pressures, as well as on their driving mechanisms. The unstable region was confirmed and the geyser and the flashing instability boundary represented by dimensionless inlet subcooling and heating power were obtained. Results indicate that the geysering and flashing instability flow occurred successively, following the decrease in inlet subcooling. Compared with conventional natural circulation loops, the circulation flow of the thimble-based heat exchange element is unique in the fact that the fluid in the rising segment transfers heat to the fluid in the downcomer; therefore, the fluid temperature decreases in the rising segment and the condensation of steam bubbles becomes more serious. This effect easily causes geysering instability and delays the occurrence of flashing instability to some degree. When the flashing phenomenon occurs, the scope of the ratio between the flow oscillation period and the necessary time of fluid passing in the rising segment becomes 0.8~1.1.
Key words: thimble-based heat transfer element    natural circulation    flow instabilities    geysering oscillations    flashing instability    unstable boundary    

自然循环已被广泛应用于先进反应堆的常规热量导出和停堆后余热导出设计中[1-3]。在两相自然循环系统中,流动过程与传热过程之间存在强烈的耦合作用,当受到内部或外部扰动时,容易发生流动不稳定现象[4]。不稳定现象可能引起系统内出现强烈机械振动和热量无法及时导出的风险,因此需要对其进行全面分析和深入了解[5]

Aritomi等[6]发现自然循环系统在低压条件下的启动过程中存在三种不稳定现象,即间歇泉不稳定流动、自然循环不稳定流动和密度波振荡,并对前两者的形成机制进行推测。Jiang等[7]通过实验总结了不同压力下自然循环系统内的不稳定现象,其中低压和中压条件下发生间歇泉、闪蒸和高频流量振荡,高压条件下发生低质量含汽率的密度波振荡。Prasad等[8]已对出现在两相自然循环系统中的不稳定性进行了广泛研究。目前对自然循环流动不稳定性的研究多集中于结构相对简单的回路[9-10]。套管式换热元件是熔盐堆非能动余热排出系统可能采用的方案之一。Beall等[11]提出将双层套管式换热元件应用于美国熔盐堆实验(MSRE)中,并验证了这种设计的可行性。套管式换热元件的优点在于采用了嵌套结构,使得回路结构紧凑,节省空间,而且双层隔离结构可以有效降低冷却剂与高温熔盐接触的风险;其缺点在于结构较为复杂,设计加工有一定难度。与一般自然循环回路相比,套管式回路中的流动和传热情况更加复杂,目前关于这种套管系统中流动不稳定性的研究还十分不足。由于流动不稳定性受系统几何结构、压力、热流密度等参数共同影响,有必要对套管式换热元件中的自然循环流动不稳定性进行研究,为熔盐堆非能动余热排出系统的设计和运行提供依据。基于此,本文以去离子水为工作介质,对自然循环条件下套管式换热元件内的不稳定特性和发生条件进行了实验研究。

1 流动不稳定性实验装置

实验装置主要由套管式换热元件、汽包、冷凝水箱和电加热装置组成,如图 1所示。冷却水由入口进入中心管后向下流动,到达底部后向上流动进入加热段,在环隙内水被加热沸腾,这样下降段中的单相流体和上升段中的两相流体产生一定的密度差,驱动自然循环持续带走热量。换热元件产生的蒸汽进入汽包后沿着蒸汽管线流向冷凝水箱,冷凝后在重力作用下又返回至汽包中。套管式换热元件包括中心管(12.7 mm×1 mm)和外层套管(25.4 mm×1.5 mm)两部分。套管的加热段长度为1.4 m,上升段长度为1.7 m。所述套管式换热元件的几何参数选取基于文献[12-13]中对熔盐堆非能动余热排出系统中套管式换热元件的初步设计。加热段上部装有一段长20 cm的玻璃视镜,用来观察环隙内的流动情况。加热段为外层套管的一部分,由不锈钢制成,加热段上下端与电源之间采用铜排连接,铜排电阻远小于加热段电阻,电源为低压直流电源,实验功率范围1~4 kW。整个实验装置包覆了保温层,实验工况下估算的散热量占电功率的0.2%左右,可认为电功率即为加热功率。

Download:
图 1 实验装置示意图 Fig. 1 Schematic diagram of experimental apparatus

由于自然循环回路驱动压头较小,为了测量循环流量时不引入过大的阻力,在汽包和中心管入口之间安装了电磁流量计(精度等级0.2)。在中心管入口、上升段出口、加热段出口和套管底部分别布置了直径1 mm的K型铠装热电偶测量当地流体温度。加热段出口和套管底部的压力通过精度等级0.1的压力传感器测得。在加热段外壁面沿轴向等距布置了7个直径0.5 mm的K型铠装热电偶,间距为167 mm。实验使用的所有热电偶均经过专门标定,测量误差小于0.7 ℃。温度、压力、流量信号都通过solartron数据采集系统实时记录并储存。所用热电偶的热响应时间小于0.05 s,电磁流量计和压力传感器的响应时间小于0.1 s,满足动态过程测量要求。

实验采用去离子水作为工质,所有实验均在大气压下进行。实验前将回路加热至饱和状态并维持一段时间,以去除不凝气体,再冷却至室温。实验时,采集不同加热功率和入口过冷度下系统内各参数变化。

2 实验结果与分析 2.1 典型流动特性变化过程

在不同入口过冷度下,换热元件内流体与加热面的换热形式发生改变,并伴随着显著不同的自然循环流动特性。图 2表示了加热功率为1.5 kW时自然循环流量的时序变化行为,可分为四个阶段:高入口过冷度下的稳定流动(a);间歇泉引起的不稳定流动(b);闪蒸引起的不稳定流动(c);最终进入稳定两相流动(d)。值得注意的是,实验中未观察到单相自然循环,这是由于套管式换热元件的流道较窄,回路阻力较大,单相时冷热流体密度差无法产生足够的驱动力。随着入口过冷度减小,系统先后经历两个流动不稳定区域,其流量振荡幅度和频率存在很大区别。

Download:
图 2 自然循环流量和入口过冷度的时序变化行为 Fig. 2 Time evolution of natural circulation flow rate and inlet subcooling
2.2 间歇泉不稳定流动

随着入口过冷度不断减小,自然循环流量逐渐变得不稳定,当入口过冷度达到31 ℃时,进入间歇泉不稳定阶段。典型的间歇泉现象如图 3所示,系统内发生间歇沸腾现象,自然循环流量呈明显的周期性波动(见图 3(a))。Paniagua等[14]对间歇泉现象的发展过程进行了总结:过冷流体进入竖直通道的加热段产生蒸汽,蒸汽向上移动时由于静压头减小体积逐渐增加,形成大汽泡;蒸汽继续向上进入上升段被过冷液体冷凝,过冷液体重新充满通道,沸腾结束;这个过程周期性重复形成流动不稳定。在本实验中,由于套管式换热元件结构的特殊性,加热段出口流体在上升过程中会受到下降段过冷流体的冷却作用(见图 3(b)),温度降低7~14 ℃。这进一步加强了上升段内汽泡的冷凝,使间歇泉更易发生。

Download:
图 3 间歇泉不稳定引起的温度和流量变化 Fig. 3 Geysering induced variations of temperature and flow rate

图 3(c)中可以发现,入口流体通过下降段时会受到上升段与加热段内流体的加热作用,温升达24 ℃,这相当于对入口流体进行了预热。Prasad等[8]认为入口过冷度达到一定阈值是间歇泉不稳定发生的必要条件。在一般的自然循环系统中,下降段不发生能量交换。而在套管式换热元件中,从下降段入口进入的流体虽然具有较大的过冷度,由于上升段与加热段的预热作用,到达加热段入口时其过冷度已显著减小,这会促进间歇泉提前发生。

2.3 闪蒸不稳定流动

当入口过冷度低于18 ℃时,观察到自然循环流量出现典型的自持振荡,振荡周期约12 s,如图 4所示。此时系统内其他参数如加热段壁面温度、加热段出口温度、底部压力等也出现周期性振荡。这种不稳定现象是由闪蒸引起的。

Download:
图 4 闪蒸不稳定引起的流量、压力、温度变化 Fig. 4 Flashing instability induced variations of flow rate, pressure and temperature

图 4Tw1~Tw4分别代表加热壁面外侧自下至上布置的4个温度测点,A点代表加热段壁面温度达到极小值的时刻,此时由于自然循环流量低,不能及时带走传递到加热面的热量,壁面温度持续上升。由于壁面过热度不断增加,汽泡开始在壁面上核化并生长。汽泡在B时刻脱离壁面进入主流,上升段与下降段中的流体密度差产生驱动力,自然循环流量开始增加。在向上流动过程中,饱和流体由于重位压降的减小形成过热,在C时刻发生闪蒸现象。在D时刻加热段出口处的流体也开始闪蒸。闪蒸极大地增加了截面含汽率,因此循环流量快速升高,强化的换热效果使壁面温度持续下降。由于流速较快,流体通过加热段的时间较短,造成截面含汽率的下降。自然循环流量回到较低的水平,加热段壁面温度再一次开始上升,形成新的循环。

结合图 4(a)(b)可以发现,自然循环流量与加热段底部压力存在反相的关系。这是由于上升段某位置发生闪蒸时,会引起气泡核化区域下方的压力降低,底部压力降低又促进了闪蒸的发生,使得闪蒸不断向下发展,流量显著增加,底部压力不断下降,因此自然循环流量的增加与加热段底部压力的降低几乎无时间差。当C时刻(图 4(c))上升段某点开始闪蒸时,加热段出口还未发生闪蒸现象(见图 4(d)),直到约1 s后的D时刻才开始闪蒸,这也印证了闪蒸向下方发展的推论。图 4(d)中通过加热段出口压力查得当地饱和温度。

Jiang等[15]指出加热段出口温度高于上升段出口饱和温度是闪蒸发生的条件,但在套管式换热元件中,尽管加热段出口温度达到101 ℃,高于上升段出口饱和温度,上升段内仍未发生闪蒸现象(见图 3)。这是由于上升过程中流体受到下降段的冷却,无法达到饱和状态。这种冷却作用在一定程度上推迟了闪蒸的发生。

2.4 两相稳定流动

在大部分入口过冷度下,流动是不稳定的,只有当入口过冷度减小到2 ℃时,自然循环进入两相稳定流动阶段,如图 5所示。此时流体在加热段中大部分区域发生饱和沸腾,加热段出口处截面含汽率已经达到较高的水平。虽然在上升段中仍会发生闪蒸,但其对截面含汽率的影响微弱,驱动力不会发生明显改变,因此流动不稳定现象消失。

Download:
图 5 两相稳定流动 Fig. 5 Two-phase stable flow
2.5 密度波不稳定的判定方法

密度波不稳定性在两相自然循环系统中具有极易影响自然循环传热和流动的特性。目前广泛认为闪蒸是密度波不稳定性的一种[16-17]图 6给出了本实验中闪蒸引起的流量振荡周期与流体通过上升段所需时间的关系,随着加热功率和入口过冷度等参数改变,振荡周期发生明显变化,但两者的比值范围几乎都保持在0.8~1.1。Furuya等[18]将流量振荡周期为流体流经上升段所需时间的1.5~2倍作为密度波不稳定发生的判定标准。由于密度波不稳定的成因是沸腾管道内循环流量、蒸汽产生速率、压降的相互反馈,这些反馈会受到回路结构尺寸、热流密度等参数影响,因此不同系统内振荡周期与上升时间的比值可能发生改变。本实验得到的比值范围可用于判别套管式换热元件内的密度波不稳定性。

Download:
图 6 闪蒸不稳定周期与流体流经上升段所需时间的关系 Fig. 6 Relationship between time required to pass through riser and oscillation period
2.6 不稳定边界

图 7所示为套管式换热元件内自然循环流动不稳定边界。其中Nsub表示过冷度数,其计算公式如下

Download:
图 7 不稳定边界 Fig. 7 Instability boundary
$ {N_{{\rm{sub}}}} = \frac{{{h_{\rm{f}}}-{h_{{\rm{in}}}}}}{{{h_{{\rm{fg}}}}}}\left( {\frac{{{\rho _{\rm{f}}}}}{{{\rho _{\rm{g}}}}}-1} \right) $ (1)

式中:hf为换热元件入口饱和水比焓,kJ/kg;hin为换热元件入口水比焓,kJ/kg;hfg为换热元件入口压力下的汽化潜热,kJ/kg;ρf为换热元件入口饱和水密度,kg/m3ρg为换热元件入口饱和蒸汽密度,kg/m3

在不同加热功率下,随着入口过冷度的减小流动都会进入不稳定状态。系统内存在两个明显的不稳定边界,即高过冷度下的不稳定边界和低过冷度下的不稳定边界。在高过冷度和低加热功率下,流动是趋于稳定的,但随着加热功率增加,系统内很容易产生高过冷度下的间歇泉不稳定,这是上升段和加热段对下降段的预热作用以及加热段内强烈的过冷沸腾共同引起的。闪蒸不稳定区域也随着加热功率增加不断扩大。当Nsub为35左右时,系统始终处于闪蒸不稳定状态,且在本实验范围内与加热功率大小无关。进一步降低入口过冷度,闪蒸对驱动力的影响逐渐减小,系统最终再次4进入两相稳定流动区域。在低过冷度条件下,提高加热功率可以增大两相稳定自然循环流动的范围,使系统更加稳定。

流动不稳定边界对系统实际运行有重要指导意义。若在较高入口过冷度下启动系统,间歇泉和闪蒸不稳定将无法避免。若将入口流体提前加热至饱和温度附近,则系统直接进入两相稳定流动状态。

3 结论

1) 随着入口过冷度减小,系统内相继发生间歇泉和闪蒸不稳定流动。与一般自然循环回路相比,套管式换热元件的上升段流体向下降段流体传热,因此上升段流体温度更低,汽泡的冷凝现象更严重,更容易出现间歇泉不稳定性。且上升段越长,这种间歇泉现象越强烈。

2) 闪蒸发生时流量振荡周期与流体流经上升段所需时间的比值范围为0.8~1.1,这与常规自然循环系统有所区别。

3) 得到了实验范围内不同入口过冷度和加热功率下的自然循环不稳定边界。若从高入口过冷度下启动系统,由于水温不断升高,间歇泉和闪蒸不稳定将无法避免。当入口过冷度很小时,提高加热功率可使两相稳定流动区域扩大。

本文中得到的不稳定边界适用于所选取的套管式换热元件。而套管几何参数改变对自然循环不稳定边界的影响,还有待于进一步研究。

参考文献
[1]
SCHULZ T L. Westinghouse AP1000 advanced passive plant[J]. Nuclear engineering and design, 2006, 236(14-16): 1547-1557. DOI:10.1016/j.nucengdes.2006.03.049 (0)
[2]
STOSIC Z V, BRETTSCHUH W, STOLL U. Boiling water reactor with innovative safety concept:The Generation Ⅲ+ SWR-1000[J]. Nuclear engineering and design, 2008, 238(8): 1863-1901. DOI:10.1016/j.nucengdes.2007.12.014 (0)
[3]
HINDS D, MASLAK C. Next-generation nuclear energy:The ESBWR[J]. Nuclear news, 2006, 49(1): 35-40. (0)
[4]
RUSPINI L C, MARCEL C P, CLAUSSE A. Two-phase flow instabilities:A review[J]. International journal of heat and mass transfer, 2014, 71: 521-548. DOI:10.1016/j.ijheatmasstransfer.2013.12.047 (0)
[5]
MANERA A, PRASSER H M, VAN D H T, et al. A comparison of void-fraction measurements during flashing-induced instabilities obtained with a wire-mesh sensor and a gamma-transmission set-up[C]//ICMF-2001, New Orleans, May, 2001. (0)
[6]
ARITOMI M, CHIANG J H, NAKAHASHI T, et al. Fundamental study on thermo-hydraulics during start-up in natural circulation boiling water reactors, (I) thermo-hydraulic instabilities[J]. Journal of nuclear science and technology, 1992, 29(7): 631-641. DOI:10.1080/18811248.1992.9731576 (0)
[7]
JIANG S Y, WU X X, ZHANG Y J. Experimental study of two-phase flow oscillation in natural circulation[J]. Nuclear science and engineering, 2000, 135(2): 177-189. DOI:10.13182/NSE00-A2133 (0)
[8]
PRASAD G V D, PANDEY M, KALRA M S. Review of research on flow instabilities in natural circulation boiling systems[J]. Progress in nuclear energy, 2007, 49(6): 429-451. DOI:10.1016/j.pnucene.2007.06.002 (0)
[9]
侯晓凡, 孙中宁, 孙秋南, 等. 开式自然循环系统中间歇喷泉流动的实验研究[J]. 哈尔滨工程大学学报, 2016, 37(4): 556-561.
HOU Xiaofan, SUN Zhongning, SUN Qiunan, et al. Experimental study on geysering in an open natural circulation system[J]. Journal of Harbin Engineering University, 2016, 37(4): 556-561. (0)
[10]
戚展飞, 佟立丽, 曹学武. 加热系统间歇泉流动现象实验研究[J]. 核动力工程, 2011, 33(2): 76-80.
QI Zhanfei, TONG Lili, CAO Xuewu. Preliminary experimental study on geysering phenomena in heating system[J]. Nuclear power engineering, 2011, 33(2): 76-80. (0)
[11]
BEALL S E, HAUBENREICH P N, LINDAUER R B, et al. MSRE design and operations report. Part Ⅴ. Reactor safety analysis report[R]. United States: Oak Ridge National Lab, Tenn, 1964. https://www.osti.gov/biblio/4034157 (0)
[12]
CHEN Kailun, YAN Changqi, MENG Zhaoming, et al. Experimental analysis on passive residual heat removal in molten salt reactor using single cooling thimble test system[J]. Energy, 2016, 112: 1049-1059. DOI:10.1016/j.energy.2016.07.004 (0)
[13]
SUN Lu, SUN Licheng, YAN Chengqi, et al. Conceptual design and analysis of a passive residual heat removal system for a 10MW molten salt reactor experiment[J]. Progress in nuclear energy, 2014, 70: 149-158. DOI:10.1016/j.pnucene.2013.09.013 (0)
[14]
PANIAGUA J, ROHATGI U S, PRASAD V. Modeling of thermal hydraulic instabilities in single heated channel loop during startup transients[J]. Nuclear engineering and design, 1999, 193(1/2): 207-226. (0)
[15]
JIANG S Y, YAO M S, BO J H, et al. Experimental simulation study on start-up of the 5 MW nuclear heating reactor[J]. Nuclear engineering and design, 1995, 158(1): 111-123. DOI:10.1016/0029-5493(95)01020-I (0)
[16]
VIJAYAN P K, NAYAK A K, SAHA D, et al. Effect of Loop Diameter on the steady state and stability behaviour of single-phase and two-phase natural circulation loops[J]. Science and technology of nuclear installations, 2008, 2008: 672-704. (0)
[17]
MANERA A, VAN DER HAGEN T H J J. Stability of natural-circulation-cooled boiling water reactors during startup:experimental results[J]. Nuclear technology, 2003, 143(1): 77-88. DOI:10.13182/NT03-A3399 (0)
[18]
FURUYA M, INADA F, VAN DER HAGEN T H J J. Flashing-induced density wave oscillations in a natural circulation BWR-mechanism of instability and stability map[J]. Nuclear engineering and design, 2005, 235(15): 1557-1569. DOI:10.1016/j.nucengdes.2005.01.006 (0)