搜索

文章查询

x

留言板

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

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

浮升力和流动加速对超临界CO2管内流动传热影响

张海松 朱鑫杰 朱兵国 徐进良 刘欢

浮升力和流动加速对超临界CO2管内流动传热影响

张海松, 朱鑫杰, 朱兵国, 徐进良, 刘欢
PDF
HTML
导出引用
导出核心图
  • 基于单相流体的概念, 超临界流体的异常传热行为已经被研究很多年了, 但是关于其流动传热机理仍没有统一的认识. 本文通过理论分析和实验研究了超临界二氧化碳在竖直管内向上流动过程中, 浮升力和流动加速效应对其流动结构和传热过程的影响. 结果表明, 没有确凿的实验证据表明超临界流体的异常传热行为是浮升力和流动加速直接导致的, 存在的估计浮升力和流动加速效应准则均是在常物性流体的基础上, 做了大量假设得出的, 不同的研究者采用浮升力和流动加速准则分析超临界流体的传热恶化得出的结论不一致. 最后, 基于拟沸腾理论分析超临界流体的传热恶化过程, 提出超临界沸腾数区分了超临界流体正常传热与恶化传热的转换边界, 为超临界流体流动传热研究提供新思路, 超临界沸腾数对建立用于不同技术的超临界流体动力循环的最佳运行条件具有重要意义.
      通信作者: 徐进良, xjl@ncepu.edu.cn
    • 基金项目: 国家级-国家重点研发计划项目(2017YFB0601801)
    [1]

    Crespi F, Gavagnin G, Sánchez, David, Martinez, Gonzalo S 2017 Appl. Energy 195 152

    [2]

    Xu J L, Sun E H, Li M J, Liu H, Zhu B G 2018 Energy 157 227

    [3]

    Ehsan M M, Guan Z, Klimenko A Y 2018 Renewable Sustainale Energy Rev. 92 658

    [4]

    Shiralkar, B S, Griffith P 1969 J. Heat Transfer 91 27

    [5]

    Bourke P J, Pulling D J, Gill L E, Denton, W H 1970 Int. J. Heat Mass Transfer 13 1339

    [6]

    Bae Y Y 2011 Nucl. Eng. Des. 241 3164

    [7]

    Brassington D J, Cairns D N H 1977 Int. J. Heat Mass Transfer 20 207

    [8]

    Hall W B, Jackson J D 1978 Advances in Heat Transfer 7 1

    [9]

    Mceligot D M, Coon C W, Perkins H C 1970 Int. J. Heat Mass Transfer 13 431

    [10]

    Liu S H, Huang Y P, Liu G X, Wang J F, Leung L K 2017 Int. J. Heat Mass Transfer 106 1144

    [11]

    Huang D, Wu Z, Sunden B, Li W 2016 Appl. Energy 162 494

    [12]

    Dang G X, Zhong F Q, Chen L H, Chang X Y 2013 Sci. China Technol. Sci. 56 416

    [13]

    Bruch A, Bontemps A, Colasson S 2009 Int. J. Heat Mass Transfer 52 2589

    [14]

    Liao S M, Zhao T S 2002 Int. J. Heat Mass Transfer 45 5025

    [15]

    Kim D E 2011 Int. J. Heat Fluid Flow 32 176

    [16]

    徐肖肖, 吴杨杨, 刘朝, 王开正, 叶建 2015 物理学报 64 054401

    Xu X X, Wu Y Y, Liu C, Wang K Z, Ye J 2015 Acta Phys. Sin. 64 054401

    [17]

    Kurganov V A, Kaptilnyi A G 1993 Int. J. Heat Mass Transfer 36 3383

    [18]

    Simeoni G G, Bryk T, Gorelli F A Krisch M, Ruocco G, Santoro M, Scopigno T 2010 Nat. Phys. 6 503

    [19]

    Gorelli F A, Bryk T, Krisch M, Ruocco G, Santoro M, Scopigno T 2013 Sci. Rep. 3 120

    [20]

    Banuti D T 2015 J. Supercrit. Fluids 98 12

    [21]

    Zhang Q, Li H X, Lei X L, Zhang J, Kong X F 2018 Int. J. Heat Mass Transfer 127 674

    [22]

    Kandlikar S G 2004 J. Heat Transfer 126 8

  • 图 1  实验系统图

    Fig. 1.  Experiment setup.

    图 2  实验段

    Fig. 2.  Vertically positioned test tube.

    图 3  圆管中简单的二维流动几何模型

    Fig. 3.  Vertically positioned test tube.

    图 4  局部壁温Tw,in, Bu, Ac随主流焓值ib的分布关系 (a) P = 8.220 MPa, G = 200 kg/(m2·s), qw = 60 kW/m2, (b) P = 8.220 MPa, G = 520 kg/(m2·s), qw = 42 kW/m2

    Fig. 4.  Local inner wall (Tw, in), Bu, Ac distributions with bulk fluid enthalpy (ib): (a) P = 8.220 MPa, G = 200 kg/(m2·s), qw = 60 kW/m2, (b) P = 8.220 MPa, G = 520 kg/(m2·s), qw = 42 kW/m2.

    图 6  不同质量流速下(a)局部壁温Tw,in, (b) Bu, (c) Ac随主流焓值ib的分布(NHT, 正常传热; HTD, 恶化传热)

    Fig. 6.  (a) Local inner wall Tw, in, (b) Bu, (c) Ac distributions with bulk fluid enthalpy ib (NHT, normal heat transfer; HTD, heat transfer deterioration).

    图 5  局部壁温Tw,in, Bu, Ac随主流焓值ib的分布关系 (a) P = 8.220 MPa, G = 700 kg/(m2·s), qw = 245 kW/m2, (b) P = 8.220 MPa, G = 1000 kg/(m2·s), qw = 245 kW/m2

    Fig. 5.  Local inner wall (Tw, in), Bu, Ac distributions with bulk fluid enthalpy (ib): (a) P = 8.220 MPa, G = 700 kg/(m2·s), qw = 245 kW/m2, (b) P = 8.220 MPa, G = 1000 kg/(m2·s), qw = 245 kW/m2.

    图 7  不同质量流速下的GrRe2.7分布关系

    Fig. 7.  Gr and Re2.7 distribution at different mass flow rates

    图 8  基于拟沸腾理论的超临界流体径向膨胀模型

    Fig. 8.  Radial expansion model of supercritical fluids based on pseudo-boiling.

    图 9  超临界沸腾数区分两个传热机制

    Fig. 9.  Supercritical boiling number distinguishes the two regimes of heat transfer.

    表 1  测量仪器的精度和范围

    Table 1.  Accuracies and ranges of measuring instruments.

    参数范围不确定度
    压力p/MPa7.510—25.231± 1.42%
    进口温度 Tin/℃5—70± 0.75%
    出口温度 Tout/℃25—500± 0.75%
    外壁面温度 Tw,o/℃30—450± 0.75%
    质量流速 G/kg·m–2·s–1488—2000± 2.05%
    热流密度qw/kW·m–230—400.36± 8.06%
    下载: 导出CSV
  • [1]

    Crespi F, Gavagnin G, Sánchez, David, Martinez, Gonzalo S 2017 Appl. Energy 195 152

    [2]

    Xu J L, Sun E H, Li M J, Liu H, Zhu B G 2018 Energy 157 227

    [3]

    Ehsan M M, Guan Z, Klimenko A Y 2018 Renewable Sustainale Energy Rev. 92 658

    [4]

    Shiralkar, B S, Griffith P 1969 J. Heat Transfer 91 27

    [5]

    Bourke P J, Pulling D J, Gill L E, Denton, W H 1970 Int. J. Heat Mass Transfer 13 1339

    [6]

    Bae Y Y 2011 Nucl. Eng. Des. 241 3164

    [7]

    Brassington D J, Cairns D N H 1977 Int. J. Heat Mass Transfer 20 207

    [8]

    Hall W B, Jackson J D 1978 Advances in Heat Transfer 7 1

    [9]

    Mceligot D M, Coon C W, Perkins H C 1970 Int. J. Heat Mass Transfer 13 431

    [10]

    Liu S H, Huang Y P, Liu G X, Wang J F, Leung L K 2017 Int. J. Heat Mass Transfer 106 1144

    [11]

    Huang D, Wu Z, Sunden B, Li W 2016 Appl. Energy 162 494

    [12]

    Dang G X, Zhong F Q, Chen L H, Chang X Y 2013 Sci. China Technol. Sci. 56 416

    [13]

    Bruch A, Bontemps A, Colasson S 2009 Int. J. Heat Mass Transfer 52 2589

    [14]

    Liao S M, Zhao T S 2002 Int. J. Heat Mass Transfer 45 5025

    [15]

    Kim D E 2011 Int. J. Heat Fluid Flow 32 176

    [16]

    徐肖肖, 吴杨杨, 刘朝, 王开正, 叶建 2015 物理学报 64 054401

    Xu X X, Wu Y Y, Liu C, Wang K Z, Ye J 2015 Acta Phys. Sin. 64 054401

    [17]

    Kurganov V A, Kaptilnyi A G 1993 Int. J. Heat Mass Transfer 36 3383

    [18]

    Simeoni G G, Bryk T, Gorelli F A Krisch M, Ruocco G, Santoro M, Scopigno T 2010 Nat. Phys. 6 503

    [19]

    Gorelli F A, Bryk T, Krisch M, Ruocco G, Santoro M, Scopigno T 2013 Sci. Rep. 3 120

    [20]

    Banuti D T 2015 J. Supercrit. Fluids 98 12

    [21]

    Zhang Q, Li H X, Lei X L, Zhang J, Kong X F 2018 Int. J. Heat Mass Transfer 127 674

    [22]

    Kandlikar S G 2004 J. Heat Transfer 126 8

  • [1] 罗奔毅, 卢义刚. 超临界点附近二氧化碳流体的声速. 物理学报, 2008, 57(7): 4397-4401. doi: 10.7498/aps.57.4397
    [2] 卢义刚, 彭健新. 运用液体声学理论研究超临界二氧化碳的声特性. 物理学报, 2008, 57(2): 1030-1036. doi: 10.7498/aps.57.1030
    [3] 付东, 王学敏, 刘建岷. 超临界二氧化碳和模型共聚物的相平衡和成核性质研究. 物理学报, 2009, 58(5): 3022-3027. doi: 10.7498/aps.58.3022
    [4] 徐肖肖, 吴杨杨, 刘朝, 王开正, 叶建. 水平螺旋管内超临界CO2冷却换热的数值模拟. 物理学报, 2015, 64(5): 054401. doi: 10.7498/aps.64.054401
    [5] 屈年瑞, 高发明. 固态二氧化碳电子结构及性能的理论研究. 物理学报, 2011, 60(6): 067102. doi: 10.7498/aps.60.067102
    [6] 蒋国平, 肖波齐, 陈玲霞, 饶连周, 王宗篪, 魏茂金. 池沸腾传热的数学分析. 物理学报, 2009, 58(4): 2523-2527. doi: 10.7498/aps.58.2523
    [7] 付鹏涛, 韩纪锋, 牟艳红, 韩丹, 杨朝文. 瑞利散射法研究超声喷流二氧化碳团簇尺度轴向分布. 物理学报, 2011, 60(5): 053602. doi: 10.7498/aps.60.053602
    [8] 程巳阳, 徐亮, 高闽光, 金岭, 李胜, 冯书香, 刘建国, 刘文清. 直射太阳光红外吸收光谱技术遥测大气中二氧化碳柱浓度. 物理学报, 2013, 62(12): 124206. doi: 10.7498/aps.62.124206
    [9] 梁捷宁, 张镭, 张武, 史晋森. 黄土高原半干旱区地表能量不闭合及其对二氧化碳通量的影响. 物理学报, 2013, 62(9): 099203. doi: 10.7498/aps.62.099203
    [10] 吕斯骅, 赵汝光, 吴凯, 王德峥, 贾金锋, 吴思诚. 氧和一氧化碳在有序合金表面Pd{001}c(2×2)-Mn上的共吸附:生成二氧化碳的微观机制. 物理学报, 1995, 44(2): 251-258. doi: 10.7498/aps.44.251
    [11] 闫晨帅, 徐进良. 超临界压力CO2在水平圆管内流动传热数值分析. 物理学报, 2020, 69(4): 044401. doi: 10.7498/aps.69.20191513
    [12] 徐 征, 徐叙, 娄志东, 徐春祥, 于 磊, 滕 枫. 电致发光加速层二氧化硅的电子高场迁移率. 物理学报, 1998, 47(1): 139-145. doi: 10.7498/aps.47.139
    [13] 雷鹏飞, 张家忠, 王琢璞, 陈嘉辉. 非定常瞬态流动过程中的Lagrangian拟序结构与物质输运作用. 物理学报, 2014, 63(8): 084702. doi: 10.7498/aps.63.084702
    [14] 郭亚丽, 魏兰, 沈胜强, 陈桂影. 双液滴撞击平面液膜的流动与传热特性. 物理学报, 2014, 63(9): 094702. doi: 10.7498/aps.63.094702
    [15] 张程宾, 许兆林, 陈永平. 粗糙纳通道内流体流动与传热的分子动力学模拟研究. 物理学报, 2014, 63(21): 214706. doi: 10.7498/aps.63.214706
    [16] 王胜, 徐进良, 张龙艳. 非对称纳米通道内流体流动与传热的分子动力学. 物理学报, 2017, 66(20): 204704. doi: 10.7498/aps.66.204704
    [17] 管仁国, 赵占勇, 黄红乾, 连超, 钞润泽, 刘春明. 冷却倾斜板熔体处理过程边界层分布及流动传热的理论研究. 物理学报, 2012, 61(20): 206602. doi: 10.7498/aps.61.206602
    [18] 王小虎, 易仕和, 付佳, 陆小革, 何霖. 二维高超声速后台阶表面传热特性实验研究. 物理学报, 2015, 64(5): 054706. doi: 10.7498/aps.64.054706
    [19] 孙良奎, 于哲峰, 黄洁. 基于超材料的平板二维定向传热结构设计. 物理学报, 2015, 64(22): 224401. doi: 10.7498/aps.64.224401
    [20] 包芸, 何建超, 高振源. 二维湍流热对流羽流运动路径对传热特性的影响. 物理学报, 2019, 68(16): 164701. doi: 10.7498/aps.68.20190323
  • 引用本文:
    Citation:
计量
  • 文章访问数:  266
  • PDF下载量:  7
  • 被引次数: 0
出版历程
  • 收稿日期:  2019-10-08
  • 修回日期:  2019-12-05
  • 刊出日期:  2020-03-01

浮升力和流动加速对超临界CO2管内流动传热影响

  • 华北电力大学能源动力与机械工程学院, 低品位能源多相流与传热北京市重点实验室, 北京 102206
  • 通信作者: 徐进良, xjl@ncepu.edu.cn
    基金项目: 国家级-国家重点研发计划项目(2017YFB0601801)

摘要: 基于单相流体的概念, 超临界流体的异常传热行为已经被研究很多年了, 但是关于其流动传热机理仍没有统一的认识. 本文通过理论分析和实验研究了超临界二氧化碳在竖直管内向上流动过程中, 浮升力和流动加速效应对其流动结构和传热过程的影响. 结果表明, 没有确凿的实验证据表明超临界流体的异常传热行为是浮升力和流动加速直接导致的, 存在的估计浮升力和流动加速效应准则均是在常物性流体的基础上, 做了大量假设得出的, 不同的研究者采用浮升力和流动加速准则分析超临界流体的传热恶化得出的结论不一致. 最后, 基于拟沸腾理论分析超临界流体的传热恶化过程, 提出超临界沸腾数区分了超临界流体正常传热与恶化传热的转换边界, 为超临界流体流动传热研究提供新思路, 超临界沸腾数对建立用于不同技术的超临界流体动力循环的最佳运行条件具有重要意义.

English Abstract

    • 以超临界流体为工作媒介的能量循环, 在热力循环效率、能量容量、安全性和经济性以及使用性能等方面具有很大的提高潜力, 因此, 被广泛地应用于工业领域中[1]. 最近, Xu等[2]提出了一个1000 MWe S-CO2的发电设计方案, 主蒸汽参数为620 ℃/30 MPa, 功率效率可达到48.37%, 这优于以超临界水为工作媒介的水蒸气朗肯循环. 毫无疑问, 这些新系统的设计和优化都需要深入地了解S-CO2独特的热工水力特性, 尤其是一定条件下发生的传热恶化现象[3].

      较早发现超临界流体传热危机现象的是Shiralkar和Griffith[4], 他们在竖直光管中发现热流密度高于一定值后, 加热壁面壁温会突然升高至863 K, 并且热流密度越高, 壁温峰值越大. 此后, 各国学者对超临界流体的传热恶化现象展开了研究, 主要包括H2O, CO2及有机工质等流体, 并提出了诸多传热恶化的发生机理和预测传热恶化发生模型[5-16]. Bourke等[5]通过对超临界二氧化碳传热实验研究发现, 高热流密度下, 传热恶化仅发生在向上流动过程中, 认为传热恶化是由浮升力引起. Bae[6]通过实验发现向下和向上流动加热过程中均存在传热恶化现象, 向下流动恶化程度要比向上恶化程度低. Brassington和Cairns[7]对超临界氦进行了研究, 发现对于向上流动, 可以用浮升力诱导传热恶化很好地解释, 但是一些壁面峰值没有出现在临界点附近, 认为浮升力并不是唯一的影响因素. Hall和Jackson[8]根据常物性的假设, 结合自然对流与强制对流解释了超临界流体局部传热恶化现象, 认为浮升力是导致恶化的原因, 并提出了一个可以预测传热恶化发生的理论模型. Mceligot等[9]提出了一个无量纲参数, 描述了由于加热而产生的流动加速度的影响. Liu等[10]通过理论分析和实验研究了浮力和加速度对混合对流的影响, 并发展了估计浮升力和流动加速效应的新准则. Huang等[11]回顾了公开文献中的浮升力和流动加速准则, 发现存在的区分强制对流和混合对流的准则估算与实验数据有较大的差异, 不同研究者采用相同的准则得到的临界值不一致. 除此之外, 有研究者发现浮升力准则临界值还与流体流动方向、工质种类有关[12]. 从公开的文献来看, 缺少确凿的实验证据, 证明超临界流体传热恶化发生时, 浮升力或流动加速是直接原因, 研究者仅仅通过流动方向和常物性假设去解释超临界流体传热异常行为, 认为浮升力或流动加速是导致超临界流体传热恶化的原因, 得出的结论仅在各自的实验范围内适用性较好, 这些无量纲准则不具有通用性. 除此之外, 研究者进行实验和理论研究时, 可能忽略了一些因素, 导致不同研究者得出了不同的结论[13-17]. 近年来, 超临界流体的单相概念不断受到物理学家的质疑, 2010年Simeoni等[18]通过非弹性X射线散射和超临界氩分子动力学证明了Widom线(WL)可以将超临界区域划分为两个不同的动力学机制: 类汽和类液区域. 2013年, Gorelli等[19]通过分子动力学模拟表明, 超临界状态下的水可以通过WL划分为类汽和类液两个区域. 2015年, Banuti[20]采用拟沸腾概念解释了超临界流体跨越WL, 从类汽到类液转变这一过程, 并给出了临界点以上汽液共存的拟沸腾线, 这个转变过程的WL和临界点的汽液共存线类似.

      本文研究了超临界CO2在垂直向上流动过程中的流动传热特性, 检查了浮升力和流动加速效应对局部传热恶化的影响, 工况参数范围为: 管子内径din = 8 mm, 质量流速G = 100—1500 kg/(m2·s), 热流密度qw = 20—400 kW/m2, 进口压力Pin = 8.220—20.822 MPa. 最后, 根据实验结果, 依据拟沸腾理论对超临界流体的恶化进行解释和研究, 为超临界布雷顿循环设计及优化提供理论依据和参考.

    • 实验系统装置如图1所示, 实验回路主要包括主循环回路系统、冷却水循环系统以及二氧化碳充液回收系统, 实验系统设计压力和温度分别为30 MPa和500 ℃. 实验所使用的工质为纯度高于99.0%的工业CO2, 工质从高压泵出来分为两路: 一路流过缓冲罐, 用以平衡系统的流量和压力脉动, 然后进入预热器并加热到实验所需要的入口温度, 最后工质流入实验段; 另一路在实验过程中作为旁路调节, 可以结合高压泵对实验流量进行精确的调节, 压力通过背压阀进行调节, 其中电加热系统采用低压大电流通过不锈钢管本身的电阻直接发热, 通过连续调节电压, 实现电加热功率连续改变, 整个系统的最大加热功率为160 kW. 实验中涉及到的主要仪器及相关参数的不确定度如表1所列.

      图  1  实验系统图

      Figure 1.  Experiment setup.

      参数范围不确定度
      压力p/MPa7.510—25.231± 1.42%
      进口温度 Tin/℃5—70± 0.75%
      出口温度 Tout/℃25—500± 0.75%
      外壁面温度 Tw,o/℃30—450± 0.75%
      质量流速 G/kg·m–2·s–1488—2000± 2.05%
      热流密度qw/kW·m–230—400.36± 8.06%

      表 1  测量仪器的精度和范围

      Table 1.  Accuracies and ranges of measuring instruments.

      实验段采用规格为$ \phi12\times 2~\rm mm $的Cr18Ni9Ti不锈钢圆管, 总长3600 mm, 其中加热段长度为2000 mm, 实验段几何结构和温压测点布置如图2所示. 实验段的两端均布置了800 mm长的流动稳定段, 作为不加热处理和出、入口截面. 两铜电极板之间的距离为有效加热长度, 通过连接铜编线接入直流电源, 在管壁上形成均匀的热流. 实验段主流温度分别由布置在进、出口处的两个铠装热电偶测量. 在实验段上取了39个等距截面, 每个截面布置一个K型(NiCr-NiSi)热电偶, 用以获得实验段外壁面温度, 整个实验段的管道包裹50 mm厚度硅酸铝保温棉, 减少实验段与外界环境的热交换, 实验段的实际加热功率Pc由下式计算:

      图  2  实验段

      Figure 2.  Vertically positioned test tube.

      ${P_{\rm{c}} } = \eta \cdot {P_{\rm{e}} },$

      式中Pe为电加热功率; η为热平衡效率, 热平衡效率的定义为实验段达到稳态热平衡时, 流过整个加热段流体的焓升与电加热功率的比值,

      $\eta = \frac{{m\left( {{i_{{\rm{b}},{\rm{out}}}} - {i_{{\rm{b}},{\rm{in}}}}} \right)}}{{UI}},$

      其中m为质量流量; ib,outib,in分别为流体的出口和进口温度对应的焓值, 这两个值分别通过NIST物性软件查得; UI分别为电压和电流. 给定入口温度和质量流速, 假设主流焓升沿加热段轴向线性增加, 根据能量守恒, 可计算任意轴向位置的流体主流焓值:

      ${i_{\rm{b}} }\left( x \right) = {i_{{\rm{b}},{\rm{in}}}} + \frac{{4{q_{\rm{w}} }x}}{{G{d_{{\rm{in}} }}}}.$

      垂直光滑圆管内壁温可简化为具有内热源的一维稳态导热问题, 圆柱坐标下的导热微分方程描述为

      $\frac{{{{\rm{d}} ^2}t}}{{{\rm{d}} {r^2}}} + \frac{1}{r}\frac{{{\rm{d}} t}}{{{\rm{d}} r}} + \frac{\varphi }{{{\lambda _{\rm{w}} }}} = 0,$

      式中φ为内热源, λw为壁面导热系数. 给定边界条件: r = R, t = tw,o, dt/dr = 0, 求解得

      ${t_{{\rm{w}},{\rm{i}}}} = {t_{{\rm{w}},{\rm{o}}}} - \frac{{{q_{\rm{w}} }{d_{{\rm{in}} }}}}{{2{\lambda _{\rm{w}} }}}\left( {\frac{1}{2} - \frac{{d_{{\rm{in}}}^2}}{{d_{{\rm{out}} }^2 - d_{{\rm{in}}}^2}}{\rm{In}} \frac{{{d_{{\rm{out}} }}}}{{{d_{{\rm{in}} }}}}} \right),$

      式中dindout分别为管子的内、外径; qw为表面热流密度; tw,itw,o分别为管子的内、外壁温. 实验段内表面热流密度可由下式计算:

      ${q_{\rm{w}} } = \frac{{{P_{\rm{c}} }}}{{\text{π} {d_{{\rm{in}} }}L}}.$

    • 考虑一半径为R的充分发展圆管管道内流动, 二维坐标系下的简单流动几何模型如图3所示. 对动量方程从壁面到通道中心积分有

      图  3  圆管中简单的二维流动几何模型

      Figure 3.  Vertically positioned test tube.

      $\frac{{\partial \bar u}}{{\partial x}}\int_0^R {\overline {\rho u} } {\rm{d}} y = - \frac{{{\rm{d}} \overline p }}{{{\rm{d}} x}}R - {\tau _{\rm{w}} } - g{\overline \rho _{\rm{b}} }R,$

      这里$ {\overline \rho _{\rm{b}} }$为流体的平均密度, τw为壁面剪切力, g为重力加速度, 为了估计由于流动加速而产生的轴向动量变化, 假设$\partial \bar u/\partial x$在通道内是均匀的, 且$\partial {\bar u_{\rm{b}} }/\partial x = 0$, 可以得到:

      $0 = \overline G \frac{{{\rm{d}} {{\bar u}_{\rm{b}} }}}{{{\rm{d}} x}} + \frac{{{\tau _{\rm{w}} }}}{R} + \frac{{\partial \tau }}{{\partial y}} + \left( {{{\overline \rho }_{\rm{b}} } - \overline \rho } \right)g,$

      式中右边第一项是流体受热膨胀引起的加速项, 第二项是壁面产生的剪切力项, 第三项是流体剪切力项, 最后一项是由于密度差引起的浮升力项. 因此, 对于竖直管道中向上流动的单相流体, 任意位置y处的剪切力可以表达为

      $\tau = {\tau _{\rm{w}} }\frac{{\left( {R - \delta } \right)}}{R} - \left( {{{\overline \rho }_{\rm{b}} } - \overline \rho } \right)g{\delta _{{\rm{Bu}} }} - \overline G \frac{{{\rm{d}} {{\bar u}_{\rm{b}} }}}{{{\rm{d}} x}}{\delta _{{\rm{Ac}} }},$

      其中δBuδAc分别为与浮升力和流动加速有关的厚度, 从(9)式可以看出, 浮升力和流动加速总是使剪切力趋于减小. 为了得到方程(9)的解, 将任意位置的剪应力积分:

      $\begin{split}\tau =\; & {\tau _{\rm{w}} }\frac{R}{{R - y}} - \frac{1}{{y - R}}\\ & \times\int_0^y {\left[ {\left( {{{\overline \rho }_{\rm{b}} } - \overline \rho } \right)g + \overline G \frac{{{\rm{d}} {{\bar u}_{\rm{b}} }}}{{{\rm{d}} x}}} \right]} \left( {y - R} \right){\rm{d}} y.\end{split}$

      边界条件为$y = 0$$\tau = {\tau _{\rm{w}} }$, ${\overline \rho _{\rm{b}} } = \overline \rho $, 则可得圆管内的圆柱形流体元的力平衡:

      $0 = \frac{{2{\tau _{\rm{w}} }}}{R} + \overline G \frac{{{\rm{d}} {{\bar u}_{\rm{b}} }}}{{{\rm{d}} x}},$

      代入(10)式, 可得

      $\frac{\tau }{{{\tau _{\rm{w}} }}} = \frac{{R - y}}{R} - \frac{1}{{{\tau _{\rm{w}} }\left( {y - R} \right)}}\int_0^y {\left( {{{\overline \rho }_{\rm{b}} } - \overline \rho } \right)g} \left( {y - R} \right){\rm{d}} y.$

      这里定义无量纲${R^ + } = R/\delta $, ${Y^ + } = y/\delta $, 上述方程的无量纲形式为

      $\begin{split}\frac{\tau }{{{\tau _{\rm{w}}}}} =\;& \frac{{{R^ + } - {Y^ + }}}{{{R^ + }}} - Bu\frac{1}{{({Y^ + } - {R^ + })}}\\ & \times\int_0^{{Y^ + }} {\frac{{({{\overline \rho }_{\rm{b}}} - \overline \rho )}}{{({{\overline \rho }_{\rm{b}}} - {{\overline \rho }_{\rm{w}}})}}} ({Y^ + } - {R^ + }){\rm{d}}{Y^ + },\end{split}$

      其中

      $Bu = \frac{{({{\overline \rho }_{\rm{b}}} - {{\overline \rho }_{\rm{w}}})g{\delta _{{\rm{Bu}} }}}}{{{\tau _{\rm{w}} }}},$

      Bu为浮升力参数, 表征了浮升力与壁面剪切应力的相对大小, 决定了右侧第二项的重要性. 如果Bu $\ll 1 $, 浮力的影响可以忽略不计; 当Bu足够大时, 剪切力分布会发生改变, 在极端情况下, 剪切力会变为负值.

    • 当流体接收到壁面的热量后, 由于热膨胀会使流体产生加速效应, 考虑物性剧烈变化的单相流体沿x轴竖直向上流动, 应用能量和质量守恒及与温度相关的流体性质, 假设流体的密度变化主要是由轴向温度变化引起的, 密度沿流动方向的变化可以表示为

      $\frac{{{\rm{d}} \rho }}{{{\rm{d}} x}} = - \frac{{\rho \beta }}{{{c_{\rm{P}} }}}\frac{{4{q_{\rm{w}} }}}{{G{d_{{\rm{in}} }}}}.$

      定义无量纲量ρ* = ρ/ρpc, X = x/din, ρpc为拟临界温度对应的密度, 则方程(15)变为

      ${\rho ^*}\frac{{{\rm{d}} {\rho ^*}}}{{{\rm{d}} X}} = - \frac{{4\beta }}{{{c_{\rm{P}} }}}\frac{{{q_{\rm{w}} }}}{G},$

      式中的β为膨胀系数, cP为定压比热. 管内流动过程中的任意横截面均满足质量守恒, 这里ρu存在如下关系:

      $\rho \frac{{{\rm{d}} u}}{{{\rm{d}} x}} = - u\frac{{{\rm{d}} u}}{{{\rm{d}} x}},$

      定义无量纲速度为

      $U = {u}/{{{u_0}}},$

      其中du/dx的值可以通过流体的膨胀率估计, 应用理想气体假设和能量守恒、质量守恒定律, 可以近似得到

      $\frac{{{\rm{d}} U}}{{{\rm{dX}} }} = - \frac{{{\rm{d}} {\rho ^*}}}{{{\rm{d}} X}}.$

      因此, 在边界层区域, 由于低密度层的存在, 壁面附近流体的速度比核心区流体的速度低, 因此, 主流流体加速度产生一个压力梯度, 该压力梯度均匀作用于边界层, 由于密度梯度发生变化, 导致速度梯度发生变化, 最终改变了近壁面剪切力分布. 为了估计流动加速对剪切力的改变, 忽略浮升力的影响, 因此, 可以得到类似上述浮升力的无量纲表达式

      $\frac{\tau }{{{\tau _{\rm{w}}}}} = \frac{{{R^ + } - {Y^ + }}}{{{R^ + }}} - Ac\frac{1}{{{Y^ + } - {R^ + }}}\int_0^{{Y^ + }} {\frac{1}{2}} ({Y^ + } - {R^ + }){\rm{d}}{Y^ + },$

      这里

      $Ac = \frac{{G{\rm{d}}\overline {{u_{\rm{b}}}} /{\rm{d}} x{\delta _{{{Ac}}}}}}{{{\tau _{\rm{w}}}}}.$

      (20)式右侧第二项的重要性将由Ac决定, Ac为流动加速参数, 表示流动加速度和壁面切应力的相对大小, 当Ac $\ll 1 $时, 流动加速效应可以被忽略. 实际上, 由于热膨胀使流体加速导致湍流边界层的“分层”现象, 对超临界流体的传热过程有着重要意义, 但是没有引起足够的重视.

      公开文献中用来估计浮升力和流动加速对剪切力改变的准则均做了大量假设, 大都是分别考虑了浮升力和加速效应, 实际上, 浮升力和加速效应总是同时存在的. 对(10)式变形, 并无量纲化, 可以得到一个同时包含浮生力项和流动加速项的无量纲方程:

      $\begin{split}\frac{\tau }{{{\tau _{\rm{w}} }}} =\; & \frac{{{R^ + } - {Y^ + }}}{{{R^ + }}} - Bu\frac{1}{{{Y^ + } - {R^ + }}}\\ & \times\int_0^{{Y^ + }} {\frac{{({{\overline \rho }_{\rm{b}} } - \overline \rho )}}{{({{\overline \rho }_{\rm{b}}} - {{\overline \rho }_{\rm{w}}})}}} ({Y^ + } - {R^ + }){\rm{d}}{Y^ + } \\ &- Ac\frac{1}{{{Y^ + } - {R^ + }}}\int_0^{{Y^ + }} {\frac{1}{2}} ({Y^ + } - {R^ + }){\rm{d}}{Y^ + }.\end{split}$

      从以上推导来看, 得到的准则均是无量纲准则, 这些准则仅在某些情况下近似满足假设条件, 并且都假定了壁面剪应力与切应力之间的一般关系. 事实上, 当传热恶化发生时, 并没有确凿的实验证据表明浮升力和流动加速效应是直接原因, 这些假设与实际情况偏离较大. 从公开的文献来看, 不同的研究者得出的结论也不一致. 除此之外, 这些模型均在“两区”模型的假设基础上进行理论分析, 这需要一个合理的特征温度, 但是“两区”模型并不是总是成立的, 模型的有效性很大程度上取决于壁面热流密度的水平. 本文选取两个典型准则即浮升力准则[8]和流动加速准则[9], 通过实验数据检验存在的经典准则的适用性, 即:

      $Bu = \frac{{G{r_{}}}}{{R{e^{2.7}}}},$

      $Ac = \frac{{4{q_{\rm{w}} }{d_{{\rm{in}}}}\beta }}{{R{e^2}\mu {c_{\rm{p}}}}},$

      式中的Gr为格拉晓夫数, Re为雷诺数.

    • 不同质量流速下的Tw,in, Bu以及Ac随主流焓值分布, 如图4图6所示, 所有传热恶化均发生在主流温度Tb < Tpc < Tw, Tpc为拟临界温度. 当质量流速为200 kg/(m2·s)的某一热流密度时(如图4(a)所示), 壁面温度急剧飞升, 然后恢复. 此时, Bu最高可达1 × 10–3的量级, 整个传热过程均远大于Bu的临界值. 对于质量流速为520 kg/(m2·s)的工况(如图4(b)所示), 恶化起始点Bu大于1 × 10–5, 但是当Bu小于其临界值时, 传热异常早已完全恢复. 两个质量流速下, 虽然壁温发生异常时, Bu临界值均大于1 × 10–5, 但是并不和浮升力准则存在单值性关系. 两个工况随着主流温度增加, Ac均先增加后减小, 壁面温度与Ac最大值没有对应关系.

      图  4  局部壁温Tw,in, Bu, Ac随主流焓值ib的分布关系 (a) P = 8.220 MPa, G = 200 kg/(m2·s), qw = 60 kW/m2, (b) P = 8.220 MPa, G = 520 kg/(m2·s), qw = 42 kW/m2

      Figure 4.  Local inner wall (Tw, in), Bu, Ac distributions with bulk fluid enthalpy (ib): (a) P = 8.220 MPa, G = 200 kg/(m2·s), qw = 60 kW/m2, (b) P = 8.220 MPa, G = 520 kg/(m2·s), qw = 42 kW/m2.

      图  6  不同质量流速下(a)局部壁温Tw,in, (b) Bu, (c) Ac随主流焓值ib的分布(NHT, 正常传热; HTD, 恶化传热)

      Figure 6.  (a) Local inner wall Tw, in, (b) Bu, (c) Ac distributions with bulk fluid enthalpy ib (NHT, normal heat transfer; HTD, heat transfer deterioration).

      当质量流速为700 kg/(m2·s)的某一热流密度下时(如图5(a)所示), 恶化与BuAc准则不存在一定的对应关系. 随着质量流速进一步增大到1000 kg/(m2·s), 发现Ac数的最大值与壁面温度最值发生了对应关系, 浮升力完全不能预测恶化起始点的发生, 如图5(b)所示. 因此, 有研究者认为流动加速是导致传热恶化的主要原因[15], 实际上, 这个准则大小仅与流体的主流温度有关, 对于大质量流速, 当流体主流温度在拟临界点附近时, Ac数通常就会达到最大值.

      图  5  局部壁温Tw,in, Bu, Ac随主流焓值ib的分布关系 (a) P = 8.220 MPa, G = 700 kg/(m2·s), qw = 245 kW/m2, (b) P = 8.220 MPa, G = 1000 kg/(m2·s), qw = 245 kW/m2

      Figure 5.  Local inner wall (Tw, in), Bu, Ac distributions with bulk fluid enthalpy (ib): (a) P = 8.220 MPa, G = 700 kg/(m2·s), qw = 245 kW/m2, (b) P = 8.220 MPa, G = 1000 kg/(m2·s), qw = 245 kW/m2.

      图6给出了不同质量流速下的壁温、BuAc随焓值分布, 其中▲和□分别对应质量流速为1500和1250 kg/(m2·s). 图6(a)给出了两者在相同压力不同热流密度下的壁温分布, 质量流速为1500 kg/(m2·s)在靠近拟临界区发生了局部恶化, 壁面温度出现峰值, 然后传热恢复, 其中a为壁温恶化起始点, b为壁温恶化完全恢复点; 质量流速为1250 kg/(m2·s)为非恶化工况, 壁面温度光滑、连续地上升, a'点、b'点分别为与1500 kg/(m2·s)质量流速相同主流温度下对应的特征壁温. 两个质量流速下的浮升力准则值沿主流温度分布如图6(b)所示, 相同的主流温度下正常传热的Bu值始终大于恶化传热, 在恶化起始点处, 大质量流速的Bu值是小质量流速的二倍, 根据Hall和Jackson[8]的理论, 此时浮升力对大质量流速的剪切力改变要大于小质量流速剪切力的改变, 传热恶化应该更严重, 但是结果恰恰相反. 流动加速的影响如图6(c)所示, Ac在靠近拟临界焓值附近同时达到了最大值, 整个过程二者在数值上并没有大的差距, 尤其是质量流速为1500 kg/(m2·s)恶化起始点对应的主流温度下, 二者的Ac分别为2.26 × 10–11和2.15 × 10–11, 质量流速为1500 kg/(m2·s)稍大一点, 但是不足以构成不同寻常传热行为的决定因素.

      Kurganov和Kaptilnyi[17]在相对较大的质量流速的实验研究过程中也发现了Hall和Jackson[8]的浮升力模型不适用, 但是他们认为该浮升力模型假定在y+ = 30时, τ/τw < 0.9, 浮升力和加速度的影响是明显的, 这过高地估计了壁面热流密度, 作者认为τ/τw < 0.97是合理的. 实际上, 在大质量流速下, 管道横截面依然会出现较大的密度梯度. 图7给出了所有数据的GrRe2.7的分布关系, □代表恶化传热, ▲代表正常传热. Hall和Jackson[8]认为当Gr > Re2.7时, 浮升力会减小剪切力可能会导致传热发生恶化, 并从理论推导出一个临界值1 × 10–5. 从图7可以看出, 这个临界值并不适用于所有数据, 仅凭GrRe2.7的相对大小不能区分出正常和恶化. 除此之外, 正常传热和恶化传热也难以区分, 这在其他的文献也被报道过[21].

      图  7  不同质量流速下的GrRe2.7分布关系

      Figure 7.  Gr and Re2.7 distribution at different mass flow rates

    • 将超临界传热类比为亚临界压力下的过冷沸腾, 对于亚临界沸腾过程, 气泡在壁面上成核并生长, 这些气泡脱离壁面与核心区的过冷液体融合并凝结, 因此, 蒸汽不会在管道的核心区中积聚. 相反, 如果气泡不离开壁面, 而是合并形成一个蒸汽层, 结果会导致相应的壁面温度急剧升高. Kandlikar[22]认为蒸汽层的形成与一个重要的无量纲参数K1有关, 对于气泡受到的蒸发动量与惯性力的大小由K1表征:

      ${K_1} = \frac{{{F_{{\rm{M}} '}}}}{{{F_{{\rm{I}}'}}}} = {\left( {\frac{{{q_{\rm{w}} }}}{{G{i_{{\rm{fg}} }}}}} \right)^2}\frac{{{\rho _{\rm{f}} }}}{{{\rho _{\rm{g}} }}} = \frac{{{\rho _{\rm{f}} }}}{{{\rho _{\rm{g}} }}}B{o^2},$

      其中Bo为沸腾数, Bo = qw/(Gifg), ifg为潜热, ρfρg分布为液相密度和汽相密度. 对于给定工质的密度比, 沸腾数对确定传热机理具有重要意义. 较大的K1值表示蒸汽充分膨胀, 使气泡附着在壁上, 当气泡的结合形成了一个蒸汽层, 导致壁面温度的飞升, 称为传热恶化. 现在我们类比亚临界沸腾过程, 来处理超临界流体的异常行为, 如图8(a)图8(b)所示, 当壁面的温度达到拟临界温度, 类汽流体首先在壁面附近产生, 并且类液态流体逐渐向类汽态流体转变. 此时, K1数变成

      图  8  基于拟沸腾理论的超临界流体径向膨胀模型

      Figure 8.  Radial expansion model of supercritical fluids based on pseudo-boiling.

      ${K_1} = \frac{{{F_{{\rm{M}} '}}}}{{{F_{{\rm{I}}'}}}} = {\left( {\frac{{{q_{\rm{w}} }}}{{G\Delta i}}} \right)^2}\frac{{{\rho _{\rm{f}} }}}{{{\rho _{\rm{g}} }}}.$

      对于亚临界压力下的沸腾, qw/ifg表示汽液界面上的蒸发质量通量, 这个蒸发过程满足汽、液平衡条件, 对于特定的工质, 潜热ifg仅与压力有关. 但是, 对于超临界流体的类液与类汽转化过程, 亚临界下的潜热的概念不再适用, 从类液到类汽的相变焓由∆i量化, ∆i被认为与拟临界焓值ipc有关, 因此, 在WL上的由类液态向类汽态转变的蒸发通量可以由qw/kipc来衡量, 为了方便计算, K1数最后变成

      ${K_1} = \frac{{{F_{{\rm{M}} '}}}}{{{F_{{\rm{I}}'}}}} = {\left( {\frac{{{q_{\rm{w}} }}}{{G{i_{{\rm{pc}} }}}}} \right)^2}\frac{{{\rho _{\rm{f}} }}}{{{\rho _{\rm{g}} }}}.$

      K1值的大小代表壁面蒸汽层的膨胀力与惯性力对比, 其中qw/Gipc为超临界沸腾数(SBO)数或拟沸腾数, ipc为给定压力下的拟临界焓值. 根据Banuti[20]的拟沸腾理论, 超临界流体跨过“相变点”吸收的能量, 不仅用于克服分子间的引力, 而且还要提供流体本身温度的升高, 拟沸腾温度随着压力升高而增大, 而ipc是流体本身的物理属性, 随着压力增大而增大, 也就是拟沸腾相变焓随着压力增大而增大, 这与Banuti[20]的拟沸腾理论一致.

      为了检验SBO数的有效性, 通过实验获得了8 mm管径下大量S-CO2流动传热实验数据, 共89个工况, 参数范围涵盖$ G\!= \!\! $100—1500 kg/(m2·s), qw = 20—400 kW/m2, pin = 8.220—20.822 MPa. 这些数据对应的壁温恶化和SBO的分布关系如图9所示, 其中红色代表恶化工况, 黑色代表正常工况. 几乎所有正常和恶化工况均可以通过5.204 × 10–4这一临界值区分开. 图中A点工况为8 mm管径下正常传热外壁温随ib的分布, 没有明显的壁面温度飞升, 壁面温度光滑分布, 其中SBO = 3.685 × 0–4, 这个值小于临界值, 而B点对应的工况明显壁面温度飞升, 为恶化传热, 其中∆T = 103 ℃为壁面温度飞升值, 其中SBO = 1.012 × 10–3, 这个值大于临界值.

      图  9  超临界沸腾数区分两个传热机制

      Figure 9.  Supercritical boiling number distinguishes the two regimes of heat transfer.

    • 采用传统的单相流体概念对超临界流体的流动传热过程进行研究已经持续了近80年了, 这些研究大都是在常物性的基础上, 做了大量假设, 考虑浮升力和流动加速效应对其流动传热的影响, 得出了大量的无量纲准则来判断其发生传热恶化的起始点, 但是这些准则不具有通用性, 单相的理想流体概念很难处理超临界流体的异常流动传热过程. 关于超临界流体有关密度的相变过程一直都是争议的焦点, 本文采用拟沸腾理论处理超临界流体的异常流动传热过程, 通过类比超临界压力与亚临界压力之间的传热, 提出超临界沸腾数SBO, 对于8 mm管径超临界二氧化碳从正常传热向恶化传热转变的临界值为5.204 × 10–4, 在一定条件下, 这个无量纲数表示低密度层膨胀引起的动量力与惯性力对比, 成功地采用拟沸腾数区分了超临界二氧化碳正常传热和恶化传热的边界, 为理解超临界流体的流动传热机理开辟了新的途径. 除此之外, 我们后续的工作会继续探讨管径对超临界二氧化碳传热的影响, 关注超临界状态下类汽和类液相互转变的微观和宏观过程的关联, 尤其是这种转变过程中的相变焓如何确定, 我们将继续采用分子动力学模拟对其进行研究.

参考文献 (22)

目录

    /

    返回文章
    返回