留言板

尊敬的读者、作者、审稿人, 关于本刊的投稿、审稿、编辑和出版的任何问题, 您可以本页添加留言。我们将尽快给您答复。谢谢您的支持!

姓名
邮箱
手机号码
标题
留言内容
验证码

竖直上升管内超临界CO2异常传热机理研究

胡娟娟 李增光 高一民 杨建 刘浪

胡娟娟, 李增光, 高一民, 等. 竖直上升管内超临界CO2异常传热机理研究[J]. 中国舰船研究, 2020, 15(X): 1–11 doi: 10.19693/j.issn.1673-3185.01921
引用本文: 胡娟娟, 李增光, 高一民, 等. 竖直上升管内超临界CO2异常传热机理研究[J]. 中国舰船研究, 2020, 15(X): 1–11 doi: 10.19693/j.issn.1673-3185.01921
HU J J, LI Z G, GAO Y M, et al. Investigation into abnormal heat transfer mechanism of supercritical carbon dioxide in vertically upward Tube[J]. Chinese Journal of Ship Research, 2020, 15(X): 1–11 doi: 10.19693/j.issn.1673-3185.01921
Citation: HU J J, LI Z G, GAO Y M, et al. Investigation into abnormal heat transfer mechanism of supercritical carbon dioxide in vertically upward Tube[J]. Chinese Journal of Ship Research, 2020, 15(X): 1–11 doi: 10.19693/j.issn.1673-3185.01921

竖直上升管内超临界CO2异常传热机理研究

doi: 10.19693/j.issn.1673-3185.01921
详细信息
    作者简介:

    胡娟娟,女,1996年生,硕士,助理工程师。研究方向:换热器强化换热,流固耦合换热仿真技术。E-mail:330686773@qq.com

    通讯作者:

    胡娟娟

  • 中图分类号: U664.15

Investigation into abnormal heat transfer mechanism of supercritical carbon dioxide in vertically upward Tube

  • 摘要:   目的  小通道内的换热已广泛应用于核反应堆与航天发动机冷却等方面。为保证该系统内热力装置安全、高效运行,须对超临界CO2的传热规律与异常传热机理进行系统、深入的研究。  方法  通过CFD数值模拟,对不同质量流速下内径为2 mm的竖直上升管中超临界CO2的换热随热流密度与压力的变化规律展开研究。并对传热强化与传热恶化工况下不同截面处径向物性参数的变化规律进行研究,深入分析异常传热机理。  结果  研究结果表明,质量流速为500~1 000 kg·m−2·s−1时,随着热流密度的升高,超临界CO2由传热强化转变为传热恶化,且在传热恶化工况下,随着热流密度的升高,超临界CO2传热恢复区的换热系数的峰值有所降低,直至消失;随着运行压力的升高,超临界CO2的异常传热程度均被削弱。通过机理分析发现,在传热强化工况下,超临界CO2的换热在拟临界点附近达到最强,主要原因在于管截面上大比热容区流体的份额显著增加;在传热恶化工况下,壁温峰值对应的截面处的速度呈“M”型分布,拐点处速度梯度降为0,湍动能也降到最低,引起传热能力变弱。  结论  研究结果对于改善小通道内超临界CO2的传热性能具有一定的指导意义。
  • 图  1  物理模型与边界条件设置

    Figure  1.  Physical model and the setting of boundary condition

    图  2  网格划分

    Figure  2.  Mesh generation

    图  3  网格无关性验证

    Figure  3.  Mesh sensitivity validation

    图  4  模型验证

    Figure  4.  Model validation

    图  5  G=500 kg·m−2·s−1时超临界CO2的换热特性随热流密度的变化

    Figure  5.  The variation of heat transfer characteristics of S-CO2 with heat flux under the condition of G=500 kg·m−2·s−1

    图  6  G=750 kg·m−2·s−1时超临界CO2的换热特性随热流密度的变化

    Figure  6.  The variation of heat transfer characteristics of S-CO2 with heat flux under the condition of G=750 kg·m−2·s−1

    图  7  G=1 000 kg·m−2·s−1时超临界CO2的换热特性随热流密度的变化

    Figure  7.  The variation of heat transfer characteristics of S-CO2 with heat flux under the condition of G=1000 kg·m−2·s−1

    图  8  G=1 000 kg·m−2·s−1q=100 kW·m−2时超临界CO2的换热特性随压力的变化

    Figure  8.  The variation of heat transfer characteristics of S-CO2 with pressure under the condition of G=1000 kg·m−2·s−1 and q=100 kW·m−2

    图  9  G=1 000 kg·m−2·s−1q=200 kW·m−2时超临界CO2的换热特性随压力的变化

    Figure  9.  The variation of heat transfer characteristics of S-CO2 with pressure under the condition of G=1 000 kg·m−2·s−1 and q=200 kW·m−2

    图  10  G=500 kg·m−2·s−1q=100 kW·m−2时超临界CO2的换热特性随压力的变化

    Figure  10.  The variation of heat transfer characteristics of S-CO2 with pressure under the condition of G=500 kg·m−2·s−1 and q=100 kW·m−2

    图  11  传热强化工况下壁温与对流换热系数随流体焓值的变化曲线

    Figure  11.  The variation of wall temperature and convective heat transfer coefficient with fluid enthalpy under the condition of heat transfer enhancement

    图  12  传热强化工况下的径向物性分布

    Figure  12.  Radial physical properties distribution under the condition of heat transfer enhancement

    图  13  传热恶化工况下壁温与对流换热系数随流体焓值的变化曲线

    Figure  13.  The variation of wall temperature and convective heat transfer coefficient with fluid enthalpy under the condition of heat transfer deterioration

    图  14  传热恶化工况下的径向物性分布

    Figure  14.  Radial physical properties distribution under the condition of heat transfer deterioration

    表  1  各截面的平均比热容

    Table  1.   Average specific heat in each section

    截面编号ib/(kJ·kg−1)cp,ave/(kJ·kg−1·K−1)
    12322.66
    22756.98
    330013.42
    432524.99
    535026.05
    637512.64
    下载: 导出CSV
  • [1] 严新平, 王佳伟, 孙玉伟, 等. 船舶余热利用sCO2布雷顿循环发电技术综述[J]. 中国机械工程, 2019, 30(8): 939–946. doi: 10.3969/j.issn.1004-132X.2019.08.009

    YAN X P, WANG J W, SUN Y W, et al. Review on sCO2 Brayton cycle power generation technology based on ship waste heat recovery utilization[J]. China Mechanical Engineering, 2019, 30(8): 939–946 (in Chinese). doi: 10.3969/j.issn.1004-132X.2019.08.009
    [2] 吴攀, 高春天, 单建强. 超临界二氧化碳布雷顿循环在核能领域的应用[J]. 现代应用物理, 2019, 10(3): 031202.

    WU P, GAO C T, SHAN J Q. Application of supercritical carbon dioxide Brayton cycle in nuclear engineering[J]. Modern Applied Physics, 2019, 10(3): 031202 (in Chinese).
    [3] 芦树平, 刘怡雯, 林原胜, 等. 基于CEFR的超临界二氧化碳循环系统特性分析[J]. 节能, 2018, 37(7): 21–26.

    LU S P, LIU Y W, LIN Y S, et al. Parameters optimization of supercritical carbon dioxide Brayton cycle based on CEFR[J]. Energy Conservation, 2018, 37(7): 21–26 (in Chinese).
    [4] 张振兴, 武林林. 应用于钠冷快堆的超临界二氧化碳布雷顿循环系统[J]. 硅谷, 2014, 7(10): 192–193.

    ZHANG Z X, WU L L. Application of supercritical carbon dioxide Brayton cycle in sodium-cooled fast reactor[J]. Echnology Forum, 2014, 7(10): 192–193 (in Chinese).
    [5] 程源源, 李涛, 张晓明, 等. 超临界流体技术在工业领域的应用研究进展[J]. 河南化工, 2018, 35(5): 6–11.

    CHENG Y Y, LI T, ZHANG X M, et al. Application research progress of supercritical fluid technology in industrial field[J]. Henan Chemical Industry, 2018, 35(5): 6–11 (in Chinese).
    [6] KIM D E, KIM M H. Experimental study of the effects of flow acceleration and buoyancy on heat transfer in a supercritical fluid flow in a circular tube[J]. Nuclear Engineering and Design, 2010, 240(10): 3336–3349. doi: 10.1016/j.nucengdes.2010.07.002
    [7] KIM D E, KIM M H. Experimental investigation of heat transfer in vertical upward and downward supercritical CO2 flow in a circular tube[J]. International Journal of Heat and Fluid Flow, 2011, 32(1): 176–191. doi: 10.1016/j.ijheatfluidflow.2010.09.001
    [8] KIM D E, KIM M H. Two layer heat transfer model for supercritical fluid flow in a vertical tube[J]. The Journal of Supercritical Fluids, 2011, 58(1): 15–25. doi: 10.1016/j.supflu.2011.04.014
    [9] ZAHLAN H, GROENEVELD D, TAVOULARIS S. Measurements of convective heat transfer to vertical upward flows of CO2 in circular tubes at near-critical and supercritical pressures[J]. Nuclear Engineering and Design, 2015, 289: 92–107. doi: 10.1016/j.nucengdes.2015.04.013
    [10] 张丽娜, 刘敏珊, 董其伍. 超临界二氧化碳细管内流动与换热特性分析研究[J]. 郑州大学学报(工学版), 2009, 30(4): 69–72.

    ZHANG L N, LIU M S, DONG Q W. Analysis of flow and heat transfer characteristics of supercritical CO2 in circular tube[J]. Journal of Zhengzhou University (Engineering Science), 2009, 30(4): 69–72 (in Chinese).
    [11] HE S, KIM W S, JACKSON J D. A computational study of convective heat transfer to carbon dioxide at a pressure just above the critical value[J]. Applied Thermal Engineering, 2008, 28(13): 1662–1675. doi: 10.1016/j.applthermaleng.2007.11.001
    [12] BAE J H, YOO J Y, CHOI H. Direct numerical simulation of turbulent supercritical flows with heat transfer[J]. Physics of Fluids, 2005, 17(10): 105104. doi: 10.1063/1.2047588
    [13] 瞿家港. 寒冷地区电厂主厂房自平衡供暖数值模拟[D]. 哈尔滨: 哈尔滨工业大学, 2018.

    QU J G. Numerical study on self-balancing heating of main power builfingin cold regions[D]. Harbin: Harbin Institute of Technology, 2018 (in Chinese).
    [14] LI Z H, WU Y X, TANG G L, et al. Comparison between heat transfer to supercritical water in a smooth tube and in an internally ribbed tube[J]. International Journal of Heat and Mass Transfer, 2015, 84: 529–541. doi: 10.1016/j.ijheatmasstransfer.2015.01.047
    [15] LIU X X, XU X X, LIU C, et al. Numerical study of the effect of buoyancy force and centrifugal force on heat transfer characteristics of supercritical CO2 in helically coiled tube at various inclination angles[J]. Applied Thermal Engineering, 2017, 116: 500–515. doi: 10.1016/j.applthermaleng.2017.01.103
    [16] VAN ABEL MARK E N, ANDERSON H, CORRADINI M L. Numerical investigation of pressure drop and local heat transfer of supercritical CO2 in printed circuit heat exchangers[C]//Proceedings of Supercritical CO2 Power Cycle Symposium. Boulder, Colorado: [s.n.], 2011.
    [17] MENTER F R. Two-equation eddy-viscosity turbulence models for engineering applications[J]. AIAA Journal, 1994, 32(8): 1598–1605. doi: 10.2514/3.12149
    [18] KOO G W, LEE S M, KIM K Y. Shape optimization of inlet part of a printed circuit heat exchanger using surrogate modeling[J]. Applied Thermal Engineering, 2014, 72(1): 90–96. doi: 10.1016/j.applthermaleng.2013.12.009
    [19] NATHAN K. An experimental study of heat transfer deterioration at supercritical pressures[D]. Ottawa: University of Ottawa, 2017.
    [20] SONG J H, KIM H Y, KIM H, et al. Heat transfer characteristics of a supercritical fluid flow in a vertical pipe[J]. The Journal of Supercritical Fluids, 2008, 44(2): 164–171. doi: 10.1016/j.supflu.2007.11.013
    [21] 张鸣远, 景思睿, 李国君. 高等工程流体力学[M]. 西安: 西安交通大学出版社, 2007.

    ZHANG M Y, JING S R, LI G J. Advanced Engineering Fluid Mechanics[M]. Xi'an: Xi'an Jiaotong University Press, 2007 (in Chinese).
  • [1] 刘浪, 李增光, 杨建, 何文强, 王青青.  超临界水在套管换热器中的流动传热试验研究 . 中国舰船研究, doi: 10.19693/j.issn.1673-3185.01809
    [2] 张涛, 蔡薇, 张天奇, 丁芳芳, 龚鑫鸿.  基于多学科融合的邮轮总布置设计机理探究(一) . 中国舰船研究, doi: 10.19693/j.issn.1673-3185.01805
    [3] 薛庆雨, 武珅, 顾湘男, 翟树成.  冰阻塞诱导的连体涡空泡及其发生机理分析 . 中国舰船研究, doi: 10.19693/j.issn.1673-3185.01687
    [4] 朱军, 朱韬, 王智宇, 夏齐强, 黄昆仑, 葛义军.  船舶参数横摇数值计算与力学机理分析 . 中国舰船研究, doi: 10.19693/j.issn.1673-3185.01618
    [5] 吴崇建, 雷智洋, 吴有生.  结构噪声核心价值与理论逻辑解读第三部分:WPA法与机理分析 . 中国舰船研究, doi: 10.19693/j.issn.1673-3185.01408
    [6] 计晨, 李素云.  石墨烯增强铝基SiC复合材料抗侵彻机理试验与数值仿真 . 中国舰船研究, doi: 10.19693/j.issn.1673-3185.01147
    [7] 赵延杰, 郝轶, 刘建湖, 张攀.  水下接触爆炸作用下泡沫夹芯板耗能机理研究 . 中国舰船研究, doi: 10.19693/j.issn.1673-3185.01168
    [8] 李瑞彪, 徐荣武, 崔立林, 余文晶.  双层圆柱壳异常噪声源定位试验研究与应用 . 中国舰船研究, doi: 10.3969/j.issn.1673-3185.2017.04.022
    [9] 华宏星, 俞强.  船舶艉部激励耦合振动噪声机理研究进展与展望 . 中国舰船研究, doi: 10.3969/j.issn.1673-3185.2017.04.002
    [10] 刘俊峰, 胡志强.  船舶碰撞机理三维解析法实现及恢复系数研究 . 中国舰船研究, doi: 10.3969/j.issn.1673-3185.2017.02.011
    [11] 杨建, 牛茂升, 郭栋, 竺晓程, 杜朝辉.  某超超临界汽轮机高压进汽腔整圈斜置静叶的数值仿真 . 中国舰船研究, doi: 10.3969/j.issn.1673-3185.2016.03.017
    [12] 王春旭, 吴崇建, 陈乐佳, 邱昌林, 熊济时.  流致噪声机理及预报方法研究综述 . 中国舰船研究, doi: 10.3969/j.issn.1673-3185.2016.01.008
    [13] 郑福明, 王畅.  摇摆运动对过冷沸腾传热特性影响的机理分析 . 中国舰船研究, doi: 10.3969/j.issn.1673-3185.2015.01.018
    [14] 张彦魁, 毕大强, 刘同和.  船舶中压电网高阻接地方式机理研究 . 中国舰船研究, doi: 10.3969/j.issn.1673-3185.2014.02.016
    [15] 马俊.  一种基于频域信息的机械系统异常振动识别方法 . 中国舰船研究, doi: 10.3969/j.issn.1673-3185.2013.03.012
    [16] 杨成, 黄贤俊, 刘培国.  基于压控导电结构的能量选择表面防护机理和仿真实验研究 . 中国舰船研究, doi: 10.3969/j.issn.1673-3185.2012.03.019
    [17] 李跃华, 徐筱欣.  舰船空气净化材料吸附CO2的分子模拟研究 . 中国舰船研究, doi: 10.3969/j.issn.1673-3185.2012.01.021
    [18] 郭君, 朱枫, 姚熊亮, 周维星.  被撞船舷侧典型结构单元的吸能机理分析 . 中国舰船研究, doi: 10.3969/j.issn.1673-3185.2011.01.003
    [19] 袁东红, 华锁宝, 顾则红, 潘政广, 鲁飞.  往复泵振动和噪声机理分析及减振降噪措施 . 中国舰船研究, doi: 10.3969/j.issn.1673-3185.2009.05.016
    [20] 宋东安, 易学勤, 温定娥.  金属挡板抑制辐射干扰的机理研究 . 中国舰船研究, doi: 10.3969/j.issn.1673-3185.2007.01.010
  • 加载中
图(14) / 表 (1)
计量
  • 文章访问数:  21
  • HTML全文浏览量:  11
  • PDF下载量:  1
  • 被引次数: 0
出版历程
  • 收稿日期:  2020-04-07
  • 修回日期:  2020-06-27
  • 网络出版日期:  2020-12-10

竖直上升管内超临界CO2异常传热机理研究

doi: 10.19693/j.issn.1673-3185.01921
    作者简介:

    胡娟娟,女,1996年生,硕士,助理工程师。研究方向:换热器强化换热,流固耦合换热仿真技术。E-mail:330686773@qq.com

    通讯作者: 胡娟娟
  • 中图分类号: U664.15

摘要:   目的  小通道内的换热已广泛应用于核反应堆与航天发动机冷却等方面。为保证该系统内热力装置安全、高效运行,须对超临界CO2的传热规律与异常传热机理进行系统、深入的研究。  方法  通过CFD数值模拟,对不同质量流速下内径为2 mm的竖直上升管中超临界CO2的换热随热流密度与压力的变化规律展开研究。并对传热强化与传热恶化工况下不同截面处径向物性参数的变化规律进行研究,深入分析异常传热机理。  结果  研究结果表明,质量流速为500~1 000 kg·m−2·s−1时,随着热流密度的升高,超临界CO2由传热强化转变为传热恶化,且在传热恶化工况下,随着热流密度的升高,超临界CO2传热恢复区的换热系数的峰值有所降低,直至消失;随着运行压力的升高,超临界CO2的异常传热程度均被削弱。通过机理分析发现,在传热强化工况下,超临界CO2的换热在拟临界点附近达到最强,主要原因在于管截面上大比热容区流体的份额显著增加;在传热恶化工况下,壁温峰值对应的截面处的速度呈“M”型分布,拐点处速度梯度降为0,湍动能也降到最低,引起传热能力变弱。  结论  研究结果对于改善小通道内超临界CO2的传热性能具有一定的指导意义。

English Abstract

胡娟娟, 李增光, 高一民, 等. 竖直上升管内超临界CO2异常传热机理研究[J]. 中国舰船研究, 2020, 15(X): 1–11 doi: 10.19693/j.issn.1673-3185.01921
引用本文: 胡娟娟, 李增光, 高一民, 等. 竖直上升管内超临界CO2异常传热机理研究[J]. 中国舰船研究, 2020, 15(X): 1–11 doi: 10.19693/j.issn.1673-3185.01921
HU J J, LI Z G, GAO Y M, et al. Investigation into abnormal heat transfer mechanism of supercritical carbon dioxide in vertically upward Tube[J]. Chinese Journal of Ship Research, 2020, 15(X): 1–11 doi: 10.19693/j.issn.1673-3185.01921
Citation: HU J J, LI Z G, GAO Y M, et al. Investigation into abnormal heat transfer mechanism of supercritical carbon dioxide in vertically upward Tube[J]. Chinese Journal of Ship Research, 2020, 15(X): 1–11 doi: 10.19693/j.issn.1673-3185.01921
    • 超临界CO2布雷顿循环因其优异的性能,在船舶余热利用、舰艇动力推进、能源发电、分布式能源等领域均已获得应用[1- 2],并展现出广阔的应用前景。与水蒸汽相比,超临界CO2有诸多优点:临界条件易实现、稳定性好;临界点附近具有高密度,可降低压缩功;功率密度高,有利于减小循环系统的尺寸[3-4]

      循环工质的换热性能与系统稳健运行密切相关。超临界压力下,CO2的物理性质十分特殊:密度与液体相近,动力黏度与气体相近,扩散系数介于气、液之间[5]。此外,在拟临界点附近,超临界CO2的定压比热容的变化非常敏感,如压力为7.5 MPa时,比热容峰值附近的温度变化0.2°C时对应的比热容变化可高达460%。物性的特殊性使超临界CO2的流动换热规律十分特殊,而在工程实际中,超临界机组工况随运行参数变化及其复杂。此外,随着科技的进步,热交换器与制冷器正在不断向微小化、紧凑化发展,在核反应堆、空调与热泵以及航天卫星等领域均有较大的应用前景。因此,很有必要对小通道内超临界CO2在加热管内的流动和换热规律进行深入的研究,并对传热强化、传热恶化机制进行探究,才可从根本上减轻传热恶化带来的危害,保障超临界机组的稳健运行。

      早在20世纪50年代,就有学者研究了超临界CO2。Kim等[6-8]对超临界CO2d=4.5 mm的上升管与下降管内的对流换热特性进行了研究,研究表明热流密度和质量流速是影响加热管的管壁温度分布的主要因素,并基于试验数据提出了换热关联式;Zahlan等[9]在亚临界、近临界和超临界压力下,分别对内径为8和22 mm的竖直上升管中超临界CO2的流动与传热特性进行了实验研究,建立了超临界CO2流动换热的数据库,并验证了其有效性;张丽娜等[10]d=4 mm的水平圆管内超临界CO2的换热特性进行了数值模拟,研究了热流密度、质量流速、压力等参数的影响。由此可见,现有研究多针对较大通道内超临界CO2的传热规律,对小通道内超临界CO2异常传热机理的研究仍不足。同时有研究表明,管径越小,轴向速度偏心程度越大,会对换热产生影响[11],常规通道内的研究结果并不一定适用于小通道。因此,对小通道内超临界CO2的异常传热机理展开研究具有重要意义。

      近年来,数值模拟方法越来越广泛地应用于超临界流体流动换热研究。直接数值模拟(DNS)可准确地对传热进行预测,但其计算量巨大,计算时间长,目前仅用于对雷诺数较低的情况进行预测[12],不适用于工程应用。雷诺平均法(RANS)在对超临界流体流动换热的数值模拟中得到了广泛应用[13],其中SST k-ω模型对超临界流体管内流动的预测效果较好[14-15]。同时,该模型结合了k-εk-ω模型的优点,即k-ε模型适用于远离壁面区域,而k-ω模型适用于近壁面区域[16-18]

      鉴于此,本文拟对d=2 mm的竖直上升管内超临界CO2流动传热特性进行数值模拟,选取多组典型数据,对模型进行验证;采用SST kω模型,研究质量流速、热流密度、压力等因素对传热的影响规律;选取典型工况,分析不同轴向截面处各物性的径向分布情况,揭示传热强化与传热恶化的发生过程以及内在机理。

    • 为保证加热段流动充分发展,需在加热段前设置绝热段。如图1所示,超临界CO2先流经内径为d=2 mm、长度Liso=280 mm的绝热流动段,在绝热段出口达到流动充分发展的湍流后进入与绝热段直径相等、长度Lh=440 mm的加热段。绝热段入口设置为质量流入口边界,加热段出口设置为压力出口边界,加热段壁面设置为无厚度、无滑移边界,加热方式采用均匀热流法。本文工况下,入口的流体焓值iin取235 kJ·kg−1,特征雷诺数11 878≤Re≤63 694;湍流度4.01%≤I ≤4.95%,均为湍流流动。

      图  1  物理模型与边界条件设置

      Figure 1.  Physical model and the setting of boundary condition

      采用三维结构化网格,利用ICEM软件对网格进行划分,圆形截面采用O型剖分,近壁面网格加密处理,径向网格尺寸以1.15的比例增长;沿管长方向采用均匀网格,网格划分结果如图2所示。为了精确地捕捉层流底层中超临界流体的流动,近壁面第1个网格的无量纲壁面距离y+必须小于1。本文工况中,网格的最大无量纲壁面距离y+=0.06,满足计算要求。

      利用Fluent软件对控制方程进行离散求解。采用控制体积法求解连续方程、动量方程与能量方程,当上述三大方程与湍流参数方程残差小于10−4时,即认为达到收敛要求。设置SIMPLEC算法进行耦合求解;离散格式均采用二阶迎风格式。考虑到超临界CO2物性对温度变化较为敏感,且对传热的影响较为显著,为准确计算各离散点处的物性,调用了制冷剂物性查询软件NIST REFPROP对物性进行即时求解。

      图  2  网格划分

      Figure 2.  Mesh generation

    • 网格密度会直接影响计算准确度,因此在数值模拟前期需要进行网格无关性验证。本文在运行压力P=8 MPa、质量流速G=100 kg·m−2·s−1、热流密度q=30 kW·m−2的工况下,采用SST k-ω模型进行网格无关性验证,逐步增加网格数量,直到计算得到的壁温−焓值变化曲线随网格数目变化不明显,即可认为此时计算结果与网格数目无关。

      图3给出了网格数量为105×104,280×104,540×104时的计算结果。由图3可以看出,网格数量较少(105×104)时,在截面流体焓值ib<290 kJ·kg−1时,壁温单调上升,但当ib>290 kJ·kg−1后,模拟得到的壁温迅速上升,并出现壁温峰值;当网格数量增加到280×104后,在整个焓值区间内,不再出现壁温峰值;进一步增加网格数量到540×104后,与网格数为280×104时相比,计算结果无明显变化。为提高计算效率,本文选用280×104网格数量进行计算。

      图  3  网格无关性验证

      Figure 3.  Mesh sensitivity validation

    • 湍流模型的选择对超临界流体计算的影响很大。为选取合适的湍流模型,本文选择了超临界流体数值模拟常用的SST k-ω模型、RNG k-ε模型和Launder-Sharma低Re-k-ε模型,分别对Kline[19]与Song[20]的实验结果进行了计算,其中包括1组传热强化工况与2组传热恶化工况。3个模型的计算结果如图4所示。

      图  4  模型验证

      Figure 4.  Model validation

      图4(a)可以看出,强化工况下,SST k-ω模型结果与实验结果吻合较好,最大误差为2.8%;而RNG k-ε模型和Launder-Sharma模型的计算值均偏离实验值,预测精度较低,最大误差分别为12%与19.6%。由图4(b)图4(c)可以看出,恶化工况下,实验数据有明显的壁温峰值出现,此时SST k-ω模型计算得到的壁温曲线中,峰值位置与实验结果基本一致,峰值高低有所出入,最大误差为18.9%;峰值外区域壁温与实验结果吻合较好,最大误差为8.7%;而RNG k-ε模型和Launder-Sharma计算得到的壁温曲线均无峰值出现,与实验得到的传热规律不符,误差高达66.7%。综上,为保证计算结果的准确性,本文选择SST k-ω模型进行数值计算。

    • 图5~图7分别给出了质量流速G=500,750,1 000 kg·m−2·s−1P=8 MPa时壁温Tw与换热系数h随热流密度的变化规律,图中ipc为流体拟临界温度处对应的焓值,Tb为流体温度。

      图  5  G=500 kg·m−2·s−1时超临界CO2的换热特性随热流密度的变化

      Figure 5.  The variation of heat transfer characteristics of S-CO2 with heat flux under the condition of G=500 kg·m−2·s−1

      图  6  G=750 kg·m−2·s−1时超临界CO2的换热特性随热流密度的变化

      Figure 6.  The variation of heat transfer characteristics of S-CO2 with heat flux under the condition of G=750 kg·m−2·s−1

      图  7  G=1 000 kg·m−2·s−1时超临界CO2的换热特性随热流密度的变化

      Figure 7.  The variation of heat transfer characteristics of S-CO2 with heat flux under the condition of G=1000 kg·m−2·s−1

      图5可以看出,质量流速G =500 kg·m−2·s−1工况下,热流密度较低(q=50 kW·m−2q/G=0.1)时,壁温最低,变化较平缓,随流体焓值的提高单调增加,无峰值出现,相应地,其换热系数最高,且在拟临界温度点附近达到峰值,为传热强化工况;随着热流密度的升高(q=75 kW·m−2q/G=0.15),拟临界点前开始出现较为平缓的壁温峰值,换热系数随之降低,达到谷值,然后随着流体焓值的增加,换热迅速恢复,在ib≈376 kJ·kg−1处达到峰值;随着热流密度的进一步升高(q=100 kW·m−2q/G=0.2),壁温在拟临界点前出现较为尖锐的峰值,换热系数也随之迅速下降,随后随着流体焓值的增加,换热系数逐渐恢复,在ib≈433 kJ·kg−1处达到最大;热流密度较高时(q=150,200 kg·m−2q/G=0.3,0.4),壁温整体水平与壁温峰值显著升高,换热系数急剧降低,且由图5(b)可以看出,在较高热流密度下,换热系数随流体焓值的变化曲线几乎重合,都处于较低水平,且在恢复区也未见换热系数峰值。

      图6可以看出,当质量流速G=750 kg·m−2·s−1时,q=50 kW·m−2q/G=0.07)时,换热性能最好,换热系数出现较为明显的尖峰;q=100 kW·m−2q/G=0.13)时,换热在低焓值区发生较为微弱的恶化,随后换热及时恢复,并出现换热系数峰值;q=200,300 kW·m−2q/G=0.27,0.4)时,壁温出现显著的峰值,换热系数出现谷值,且恢复区换热系数仍处于较低水平。质量流速G=1 000 kg·m−2·s−1的壁温与换热系数随焓值变化趋势如图7所示,其换热规律与G=500,750 kg·m−2·s−1时的较为一致,此处不再赘述。

    • 为研究超临界CO2换热特性随运行压力的变化规律,选取1个传热强化工况(G=1 000 kg·m−2·s−1q=100 kW·m−2)与2个传热恶化工况(G=1 000 kg·m−2·s−1q=200 kW·m−2G=500 kg·m−2·s−1q=100 kW·m−2),在运行压力P分别为8,9,10 MPa时对壁温与换热系数随流体焓值的变化规律进行了对比,对比结果如图8~图10所示。

      图  8  G=1 000 kg·m−2·s−1q=100 kW·m−2时超临界CO2的换热特性随压力的变化

      Figure 8.  The variation of heat transfer characteristics of S-CO2 with pressure under the condition of G=1000 kg·m−2·s−1 and q=100 kW·m−2

      图  9  G=1 000 kg·m−2·s−1q=200 kW·m−2时超临界CO2的换热特性随压力的变化

      Figure 9.  The variation of heat transfer characteristics of S-CO2 with pressure under the condition of G=1 000 kg·m−2·s−1 and q=200 kW·m−2

      图  10  G=500 kg·m−2·s−1q=100 kW·m−2时超临界CO2的换热特性随压力的变化

      Figure 10.  The variation of heat transfer characteristics of S-CO2 with pressure under the condition of G=500 kg·m−2·s−1 and q=100 kW·m−2

      图8可以看出,在传热强化工况下,不同压力下壁温变化趋势基本一致,均随着流体焓值的升高,呈现平稳上升的趋势,且随着压力的升高,壁温略有升高;不同压力下的换热系数在拟临界点附近出现峰值,随着压力的升高,换热系数峰值有所降低,位置基本不变,传热强化强度减弱。

      图9给出了传热恶化工况G=1 000 kg·m−2·s−1q=200 kW·m−2下的换热特性随压力的变化规律。可以看出,在该工况下,不同压力下的壁温曲线趋势呈现出不同的规律:压力为8 MPa时,壁温在低焓值区出现较为平缓的凸起;压力为9,10 MPa时,壁温随流体焓值变化单调上升,无峰值出现。换热系数的趋势较为一致,均在较低焓值处出现换热系数低谷,不同压力下低谷处对应的流体焓值基本一致,然后换热逐渐恢复,出现换热系数峰值。随着压力的升高,换热系数显著升高,传热恶化程度减弱。

      图10给出了传热恶化工况G=500 kg·m−2·s−1q=100 kW·m−2下的换热特性随压力的变化规律。由图10可以看出,不同压力下的壁温曲线趋势呈现出不同的规律:压力为8 MPa时,壁温在低焓值区出现非常尖锐的峰值;压力为9,10 MPa时,壁温随流体焓值变化单调上升,无峰值出现。压力为8 MPa时,换热系数在加热段入口处急剧降低,并达到最低值,随后换热系数缓慢上升;压力为9,10 MPa时,换热系数显著升高,传热恶化程度减弱。

      综上可以得出如下结论:在传热强化工况下,随着压力的升高,传热强化效果有所削弱;传热恶化工况下,随着压力的升高,传热恶化情况也有所改善。这是由于随着压力的升高,压力离临界压力越来越远,定压比热容、密度、动力粘度等物性对温度的敏感程度降低,变化由剧烈变得平缓,因此由物性变化引起的传热强化与恶化现象均有所减弱。

    • 选取传热强化工况P=8 MPa,q =50 kW·m−2G=500 kg·m−2·s−1,进行机理分析。图11给出了上述工况下的壁温与对流换热系数随流体焓值的变化曲线。

      图  11  传热强化工况下壁温与对流换热系数随流体焓值的变化曲线

      Figure 11.  The variation of wall temperature and convective heat transfer coefficient with fluid enthalpy under the condition of heat transfer enhancement

      图11可以看出,该工况下壁温随流体焓值的升高单调增加,无明显的峰值或低谷出现。随着流体焓值的升高,管壁与流体的温度差逐渐减小,在拟临界焓值附近达到最小,然后又逐渐增大。对流换热系数的变化与之呼应,也呈现先增大后减小的趋势,在拟临界焓值附近达到最佳换热效果,即发生传热强化。

      为研究传热强化机理,分别截取绝热段ib=232 kJ·kg−1,加热段ib=275,300,325,350,375 kJ·kg−1截面,得到各截面处参数的径向分布情况如图12所示,其中横坐标为无量纲壁面距离y+,其定义式如下:

      图  12  传热强化工况下的径向物性分布

      Figure 12.  Radial physical properties distribution under the condition of heat transfer enhancement

      $${y^ + } = \frac{{y\sqrt {{{{\tau _{\rm{w}}}} / \rho }} }}{{{\mu / \rho }}}$$ (1)

      式中,τw为壁面剪切应力,Pa,可通过数值模拟软件输出。一般认为,y+<5时工质处于粘性底层区;5≤y+<30时处于过渡区;30≤y+<60时处于对数律层;y+≥60时处于湍流核心区[21]

      图12给出了管径 d =2 mm时,管内流体在传热强化工况下的物性径向分布。图中u为轴向速度;( ρbρ)g· ρ−1为浮升力,其中 ρb为截面平均密度, g为重力加速度;k为湍动能;λ为导热系数。由图12(a)可以看出,随着y+的增大,加热段流体温度逐渐降低,并在轴线处达到最低,随着温度的变化,流体物性也发生变化;绝热段温度不变,相应地,物性也不变。由图12(b)可知,在近壁面处,流体密度较低,为保持质量守恒,流体迅速向上流动,近壁面处流速提高,如图12(c)所示;同时,由于近壁面处流体密度较低也导致此处浮升力较高,如图12(d)所示,较高的浮升力也使得近壁面处流体速度有所提高。由于湍动能与雷诺切应力正相关,雷诺切应力又随速度梯度的增大而增大,因此随速度梯度先增大后减小,湍动能也呈现先增大后减小的趋势,如图12(e)所示。同时由图12(e)可知,随流体焓值的增加,湍动能不断增强,湍动能峰值也越来越高。

      ib=275~350 kJ·kg−1时,截面流体会达到或跨越拟临界温度,此时大比热容区形成。如图12(f)所示,随着ib的增加,大比热容区逐渐向远离壁面的方向扩散,大比热容区流体占据的份额越来越大,ib=350 kJ·kg−1时,大比热容区开始扩散到过渡区,大比热容区流体份额达到最大,此时换热系数达到峰值。ib>350 kJ·kg−1后,截面流体温度均已跨越拟临界温度,此时比热容cp单调变化,峰值消失,大比热容区流体所占份额减少,换热系数有所降低。

      为定量计算各截面流体平均比热容水平,引入截面平均比热容${c_{{\rm{p,ave}}}}$,其定义如下:

      $${c_{{\rm{p,ave}}}} = {{\int_{y = 0}^{y = r} {{c_{\rm{p}}}\rho {u_{\rm{c}}}{\rm{d}}A} } \Big/ {\int_{y = 0}^{y = r} {\rho u{\rm{d}}A} }}$$ (2)

      式中:r为截面半径;uc为工质流速;A为截面面积。

      由式(2)求得6个截面处的截面平均比热容如表1所示。可以看出,随着流体焓值的升高,截面平均比热容先迅速升高,在ib=350 kJ·kg−1达到最高,然后随着截面流体焓值的增加,平均比热容有所降低。这与上文的定性分析结果一致。

      表 1  各截面的平均比热容

      Table 1.  Average specific heat in each section

      截面编号ib/(kJ·kg−1)cp,ave/(kJ·kg−1·K−1)
      12322.66
      22756.98
      330013.42
      432524.99
      535026.05
      637512.64
    • 选取数值模拟工况P=8 MPa,G=500 kg·m−2·s−1q=100 kW·m−2作为传热恶化典型工况进行分析。图13给出了该工况下的壁温与对流换热系数随流体焓值的变化曲线。由图13可以看出,进入加热段后,壁温急剧升高,传热恶化现象逐渐发生,并在ib=300 kJ·kg−1附近达到峰值,然后壁温逐渐降低,变化趋于平缓,传热逐渐恢复;与此对应,进入加热段后换热系数急剧降低,在ib=300 kJ·kg−1附近达到谷值,然后换热系数缓慢升高,换热有所恢复。可以看出,ib=300 kJ·kg−1为该工况下传热恶化最为剧烈的地方,因此在接下来的机理分析中重点分析此处的物性变化特征。

      图  13  传热恶化工况下壁温与对流换热系数随流体焓值的变化曲线

      Figure 13.  The variation of wall temperature and convective heat transfer coefficient with fluid enthalpy under the condition of heat transfer deterioration

      为研究传热恶化机理,分别截取绝热段ib=232 kJ·kg−1,加热段ib=250,275,300,350,400 kJ·kg−1这6个截面,得到各截面处参数的径向分布情况如图14所示。由图14(a)可见,各截面处壁面温度均高于拟临界温度,与传热强化工况相比,传热恶化工况下,径向温度变化十分剧烈,在ib由235 kJ·kg−1增大到300 kJ·kg−1的过程中,壁温越来越高,径向温度梯度也越来越大。受到温度分布的影响,近壁面低密度区形成,并在ib=300 kJ·kg−1处达到最低值,如图14(b)所示。在低密度区域的影响下,近壁面浮升力增大使得流体加速,如图14(c)图14(d)所示。随着速度梯度的增大,粘性底层区的湍动能有所增加。随着流动发展到过渡区,流体密度的急剧上升使浮升力迅速下降,速度沿径向的变化逐渐变得平缓,此时,雷诺切应力减小,因此湍流强度逐渐减弱。由图14(c)可以看出,ib=300 kJ·kg−1时,出现速度拐点,拐点处速度梯度降为0,由于管内流场为轴对称分布,此时管内呈现“M”型速度曲线;与之对应,湍动能也达到最低,此时传热传质剧烈程度降低,造成热量积聚,发生传热恶化。ib > 300 kJ·kg−1后,“M”型速度曲线消失,速度梯度增大使湍动能升高,换热性能得到改善,壁温有所降低。

      图  14  传热恶化工况下的径向物性分布

      Figure 14.  Radial physical properties distribution under the condition of heat transfer deterioration

    • 本文利用经验证的数值模拟方法,对超临界CO2换热特性随不同实验参数的变化规律进行了研究。选取典型工况,深入分析了异常传热机理。在本文研究的参数范围内,得到如下结论:

      1) 热流密度较低(q/G<0.13)时,壁温沿管长方向单调增长,换热系数在拟临界点附近达到峰值;热流密度较高(q/G≥0.13)时,出现壁温峰值,且随着热流密度的增加,传热恢复区换热系数峰值逐渐降低,并逐步向流体焓值较大区域推移,直至消失。

      2) 随着压力远离临界压力,超临界流体物性变化剧烈程度减弱,由此引起的传热强化与传热恶化程度也随之削弱。

      3) 传热强化工况下,径向湍动能随着流体焓值的增大而增大,截面流体平均比热容随流体焓值的增大而先增大后减小,在ib=350 kJ·kg−1达到最大,此时换热达到最强,因此大比热容区流体份额显著增加是传热强化工况下换热性能变好的主要原因。

      4) 传热恶化工况下,受流体径向密度的影响,近壁面浮升力有所增大,使得流体加速流动,随着流动发展到过渡区,流体密度急剧上升使浮升力迅速下降,速度也随之下降。ib=300 kJ·kg−1时,出现“M”型速度曲线,在拐点处速度梯度降为0,湍动能也达到最低,发生传热恶化,此时速度梯度降低是发生传热恶化的主要原因。

参考文献 (21)

目录

    /

    返回文章
    返回