搜索

x

留言板

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

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

超临界CO2池式传热流固耦合传热特性数值模拟

于博文 何孝天 徐进良

引用本文:
Citation:

超临界CO2池式传热流固耦合传热特性数值模拟

于博文, 何孝天, 徐进良

Numerical simulation of fluid-structure coupled heat transfer characteristics of supercritical CO2 pool heat transfer

Yu Bo-Wen, He Xiao-Tian, Xu Jin-Liang
PDF
HTML
导出引用
  • 超临界类沸腾传热已有初步研究, 对类气-类液界面的定义仍不统一. 采用层流模型对超临界CO2池式传热特性进行流固耦合数值模拟研究, d = 70 μm的铂丝为加热件, 热流密度qw范围为0—2000 kW/m2, 压力P范围为8—10 MPa, 加热丝附近采用多尺度网格进行建模, 模拟值与实验数据吻合较好. 通过与亚临界传热类比通过计算导热占比Qcon/Qt沿r方向变化规律, 将超临界划分为3个区, 当T < TL时为类液区 (liquid-like, LL), TL < T < TM时为类液向类气转变的过度区 (two-phase-like, TPL), T > TM为类气区 (vapor-like, VL), 与用超临界类沸腾干度x分区规律相同; 根据通过计算导热占比确定的类气膜内平均导热λave和导热热阻RG随热流密度qw变化曲线可很好地解释类蒸发区换热系数hqw的变化规律. 本文提出了确定超临界池式传热类气-类液界面的新方法, 较好地解释了类蒸发区传热机理, 为后续超临界类沸腾池式传热的发展提供了理论基础.
    The heat transfer of supercritical pseudo-boiling has been preliminarily studied, but the definition of gas-liquid interface is still not unified. The fluid-structure coupling numerical simulation of heat transfer characteristics in supercritical CO2 pool is carried out by using laminar flow model. Platinum wire is the heating element, with diameter d = 70 μm. The heat flux density qw is in a range of 0–2000 kW/m2, and the pressure P is in a range of 8–10 MPa. Multi-scale mesh is used to model the heating wire, and simulation values accord well with the experimental data. The results show that due to the increase of the circumferential average Rayleigh number Raave of the heating filament with qw, the characteristic of the natural convection zone is that h increases with qw. The temperatures of the four characteristic working conditions in the evaporation-like zone show a downward trend along the r direction. Through analogy with subcritical heat transfer and by calculating the thermal conductivity ratio Qcon/Qt, the supercritical is divided into three regions, T < TL is liquid-like region (LL), TL < T < TM is two-phase-like region (TPL), T > TM is vapor-like region (VL). The rule is the same as that of x partition according to supercritical pseudo-boiling dryness. According to the curves of average thermal conductivity λave and thermal resistance RG versus heat flux qw, determined by calculating thermal conductivity ratio, the variation law of heat transfer coefficient h with qw in evaporation-like region can be well explained, as qw increases, the thermal conductivity thermal resistance RG increases, and the heat from the heating filament is difficult to transfer to the fluid outside the vapor-like membrane, leading the heat transfer coefficient h to decrease when qA < qw < qC, and a significant increase in λave when qw > qC, and the recovery of heat transfer when h rises again. In this paper, a new method of determining the gas-liquid interface of supercritical pool heat transfer is proposed. This method can effectively explain the heat transfer mechanism in the evaporation-like zone, and provide a theoretical basis for developing supercritical pool heat transfer in the future.
      通信作者: 徐进良, xjl@ncepu.edu.cn
    • 基金项目: 国家自然科学基金重点项目(批准号: 52130608)和国家自然科学基金创新研究群体科学基金(批准号: 51821004)资助的课题.
      Corresponding author: Xu Jin-Liang, xjl@ncepu.edu.cn
    • Funds: Project supported by the Key Program of the National Natural Science Foundation of China (Grant No. 52130608) and the Science Fund for Creative Research Groups of the National Natural Science Foundation of China (Grant No. 51821004).
    [1]

    Tsai W C, Wang Y D 2019 Prog. Polym. Sci. 98 101161Google Scholar

    [2]

    Knez Z, Markocic E, Leitgeb M, Primozic M, Hrncic M K, Skerget M 2014 Energy 77 235Google Scholar

    [3]

    Sarkar J 2015 Renew Sust. Energ Rev. 48 434Google Scholar

    [4]

    Liu Y P, Wang Y, Huang D G 2019 Energy 189 115900Google Scholar

    [5]

    Kafengauz N L, Fedorov M I 1966 J. Eng. Phys. 11 99Google Scholar

    [6]

    Ackerman J W 1970 J. Heat Trans. 92 490Google Scholar

    [7]

    Knapp K K, Sabersky R H 1965 Int. J. Heat Mass Transfer 9 41Google Scholar

    [8]

    Abadzic E, Goldstein R J 1970 Int. J. Heat Mass Transfer 13 1163Google Scholar

    [9]

    Tamba J, Takahashi T, Ohara T, Aihara T 1998 Exp. Therm. Fluid Sci. 17 248Google Scholar

    [10]

    Zhu B G, Xu J L, Yan C S, Xie J 2020 Int. J. Heat Mass Transfer 148 119080Google Scholar

    [11]

    Xu J L, Zhang H S, Zhu B G, Xie J 2020 Sol. Energy 195 27Google Scholar

    [12]

    Zhua B G, Xu J L, Wua X M, Xie J, Li M J 2019 Int. J. Therm. Sci. 136 254Google Scholar

    [13]

    Huang D, Wu Z, Sunden B, Li W 2016 Appl. Energ. 162 494Google Scholar

    [14]

    Rousselet Y, Warrier G 2012 J. Heat Transfer 135 022501Google Scholar

    [15]

    Warrier G, Rousselet Y 2012 J. Heat Transfer 135 022502Google Scholar

    [16]

    张海松, 徐进良, 朱鑫杰 2021 物理学报 70 044401Google Scholar

    Zhang H S, Xu J L, Zhu X J 2021 Acta Phys. Sin. 70 044401Google Scholar

    [17]

    Cheng L Y, Wang Q Y, Xu J L, 2024 Int. J. Heat Mass Transfer 220 124953Google Scholar

    [18]

    He X T, Xu J L, Yu X J, Xie J 2023 Int. J. Heat Mass Transfer 214 124417Google Scholar

    [19]

    何孝天, 徐进良, 程怡玮 2023 物理学报 72 057801Google Scholar

    He X T, Xu J L, Cheng Y W 2023 Acta Phys. Sin. 72 057801Google Scholar

    [20]

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

    [21]

    Maxi F, Contescu C, Boillat P, Niceno B, Karalis K, Testino A, Ludwig C H 2019 Nat. Commun. 19 12117Google Scholar

    [22]

    Banuti D T 2015 J. Supercrit. Fluid 98 12Google Scholar

    [23]

    Xu J L, Wang Y, Ma X J 2021 Phys. Rev. 104 014142Google Scholar

    [24]

    Wang Q Y, Ma X J, Xu J L, Li M J, Wang Y 2021 Int. J. Heat Mass Transfer 181 121875Google Scholar

    [25]

    闫晨帅, 徐进良 2020 物理学报 4 044401Google Scholar

    Yan C S, Xu J L 2020 Acta Phys. Sin. 4 044401Google Scholar

    [26]

    Tripathi P M, Basu S 2021 Phys. Fluids 33 043304Google Scholar

    [27]

    Li X L, Yu X Y, Liu P T, Fan Y H, Yang D L, Tang G H 2023 Int. J. Heat Mass Transfer 216 124563Google Scholar

    [28]

    Wang J T, Li Z L, Zhai Y L, Wang H 2023 Int. J. Heat Mass Transfer 201 123571Google Scholar

    [29]

    于博文, 何孝天, 徐进良 2023 中国科学: 技术科学 3 25Google Scholar

    Yu B W, He X T, Xu J L 2023 Sci. Sin. Tech. 3 25Google Scholar

    [30]

    杨世铭, 陶文铨 2011 传热学 (北京: 高等教育出版社)

    Yang S M, Tao W Q 2011 Heat Transfer Science (Beijing: Higher Education Press

    [31]

    He X T, Xu J L, Xie J 2024 Phys. Fluids 36 013343Google Scholar

    [32]

    Liu M, Liu S H, Xi D P, Huang S F, Huang Y P 2023 Nucl. Energy 104 2391Google Scholar

  • 图 1  几何模型和网格划分示意图

    Fig. 1.  Geometric model and grid diagram.

    图 2  模型验证 (a)热流密度随Tw-Tpc变化; (b)平均瑞利数Raave随热流密度变化; (c)热流密度随壁温变化; (d)换热系数随热流密度变化规律

    Fig. 2.  Model verification: (a) Heat flux varies with Tw-Tpc, (b) Raave varies with heat flux; (c) heat flux varies with wall temperature; (d) heat transfer coefficient varies with heat flux.

    图 3  不同工况温度和密度沿r方向变化规律 (a), (b) Case A, qw = 87.73 kW/m2; (c), (d) P = 10 MPa, Tb = 288 K; θ = 90°

    Fig. 3.  Temperature and density change in r direction under different working conditions: (a), (b) Case A, qw = 87.73 kW/m2; (c), (d) P = 10 MPa, Tb = 288 K; θ = 90°.

    图 4  Qcon计算网格二维示意图

    Fig. 4.  Two-dimensional diagram of Qcon computing grid.

    图 5  qw = 233.33 kW/m2, θ = 90°时, 亚临界与超临界导热占比Qcon/Qt沿r方向变化类比 (a) 亚临界P = 6 MPa; (b) 超临界P = 10 MPa

    Fig. 5.  Subcritical and supercritical thermal conductivity ratio Qcon/Qt changes along the r direction with qw = 233.33 kW/m2, θ = 90°: (a) Subcritical P = 6 MPa; (b) supercritical P = 10 MPa.

    图 6  不同压力下速度矢量图 (a) 亚临界P = 6 MPa; (b) 超临界P = 10 MPa

    Fig. 6.  Velocity vector diagram at different pressures: (a) Subcritical P = 6 MPa; (b) supercritical P = 10 MPa.

    图 7  Qcon/Qtr变化规律 (a) Case A工况; (b) 同一角度不同工况下Qcon/Qtr变化规律

    Fig. 7.  Qcon/Qt changes with r: (a) Case A working condition; (b) Qcon/Qt changes with r under the same angle and different working conditions.

    图 8  不同工况下Qcon/Qt随类沸腾干度x变化曲线

    Fig. 8.  Qcon/Qt curve with boiling dryness x under different working conditions.

    图 9  (a) 导热热阻RG随热流密度变化曲线; (b) 平均导热λave随热流密度变化曲线

    Fig. 9.  (a) Curve of thermal resistance RG with heat flux; (b) average thermal conductivity λave with heat flux curve.

    表 1  网格独立性验证

    Table 1.  Grid independence verification.

    Grid No. 径向节点数 周向节点数 Tw/K 网格数 误差/%
    1 100 21 310.2 407.7×104 –3.71
    2 100 41 316.0 82.6×104 –1.88
    3 100 81 321.3 169.8×104 –0.26/0.23
    4 100 121 322.1 261.4×104 0
    5 40 81 328.7 156.3×104 2.08
    6 70 81 326.1 163.1×104 1.29
    7 130 81 322.0 176.5×104 0
    下载: 导出CSV
  • [1]

    Tsai W C, Wang Y D 2019 Prog. Polym. Sci. 98 101161Google Scholar

    [2]

    Knez Z, Markocic E, Leitgeb M, Primozic M, Hrncic M K, Skerget M 2014 Energy 77 235Google Scholar

    [3]

    Sarkar J 2015 Renew Sust. Energ Rev. 48 434Google Scholar

    [4]

    Liu Y P, Wang Y, Huang D G 2019 Energy 189 115900Google Scholar

    [5]

    Kafengauz N L, Fedorov M I 1966 J. Eng. Phys. 11 99Google Scholar

    [6]

    Ackerman J W 1970 J. Heat Trans. 92 490Google Scholar

    [7]

    Knapp K K, Sabersky R H 1965 Int. J. Heat Mass Transfer 9 41Google Scholar

    [8]

    Abadzic E, Goldstein R J 1970 Int. J. Heat Mass Transfer 13 1163Google Scholar

    [9]

    Tamba J, Takahashi T, Ohara T, Aihara T 1998 Exp. Therm. Fluid Sci. 17 248Google Scholar

    [10]

    Zhu B G, Xu J L, Yan C S, Xie J 2020 Int. J. Heat Mass Transfer 148 119080Google Scholar

    [11]

    Xu J L, Zhang H S, Zhu B G, Xie J 2020 Sol. Energy 195 27Google Scholar

    [12]

    Zhua B G, Xu J L, Wua X M, Xie J, Li M J 2019 Int. J. Therm. Sci. 136 254Google Scholar

    [13]

    Huang D, Wu Z, Sunden B, Li W 2016 Appl. Energ. 162 494Google Scholar

    [14]

    Rousselet Y, Warrier G 2012 J. Heat Transfer 135 022501Google Scholar

    [15]

    Warrier G, Rousselet Y 2012 J. Heat Transfer 135 022502Google Scholar

    [16]

    张海松, 徐进良, 朱鑫杰 2021 物理学报 70 044401Google Scholar

    Zhang H S, Xu J L, Zhu X J 2021 Acta Phys. Sin. 70 044401Google Scholar

    [17]

    Cheng L Y, Wang Q Y, Xu J L, 2024 Int. J. Heat Mass Transfer 220 124953Google Scholar

    [18]

    He X T, Xu J L, Yu X J, Xie J 2023 Int. J. Heat Mass Transfer 214 124417Google Scholar

    [19]

    何孝天, 徐进良, 程怡玮 2023 物理学报 72 057801Google Scholar

    He X T, Xu J L, Cheng Y W 2023 Acta Phys. Sin. 72 057801Google Scholar

    [20]

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

    [21]

    Maxi F, Contescu C, Boillat P, Niceno B, Karalis K, Testino A, Ludwig C H 2019 Nat. Commun. 19 12117Google Scholar

    [22]

    Banuti D T 2015 J. Supercrit. Fluid 98 12Google Scholar

    [23]

    Xu J L, Wang Y, Ma X J 2021 Phys. Rev. 104 014142Google Scholar

    [24]

    Wang Q Y, Ma X J, Xu J L, Li M J, Wang Y 2021 Int. J. Heat Mass Transfer 181 121875Google Scholar

    [25]

    闫晨帅, 徐进良 2020 物理学报 4 044401Google Scholar

    Yan C S, Xu J L 2020 Acta Phys. Sin. 4 044401Google Scholar

    [26]

    Tripathi P M, Basu S 2021 Phys. Fluids 33 043304Google Scholar

    [27]

    Li X L, Yu X Y, Liu P T, Fan Y H, Yang D L, Tang G H 2023 Int. J. Heat Mass Transfer 216 124563Google Scholar

    [28]

    Wang J T, Li Z L, Zhai Y L, Wang H 2023 Int. J. Heat Mass Transfer 201 123571Google Scholar

    [29]

    于博文, 何孝天, 徐进良 2023 中国科学: 技术科学 3 25Google Scholar

    Yu B W, He X T, Xu J L 2023 Sci. Sin. Tech. 3 25Google Scholar

    [30]

    杨世铭, 陶文铨 2011 传热学 (北京: 高等教育出版社)

    Yang S M, Tao W Q 2011 Heat Transfer Science (Beijing: Higher Education Press

    [31]

    He X T, Xu J L, Xie J 2024 Phys. Fluids 36 013343Google Scholar

    [32]

    Liu M, Liu S H, Xi D P, Huang S F, Huang Y P 2023 Nucl. Energy 104 2391Google Scholar

  • [1] 赵明慧, 刘忠军, 姬帅, 刘晨, 敖庆波. 超临界氮气在单壁碳纳米管内吸附行为的GCMC模拟研究. 物理学报, 2022, 71(22): 220201. doi: 10.7498/aps.71.20220765
    [2] 牛越, 包为民, 李小平, 刘彦明, 刘东林. 大功率热平衡感应耦合等离子体数值模拟及实验研究. 物理学报, 2021, 70(9): 095204. doi: 10.7498/aps.70.20201610
    [3] 陈国华, 石科军, 储进科, 吴昊, 周池楼, 肖舒. 环形磁场金属等离子体源冷却流场的数值模拟与优化. 物理学报, 2021, 70(7): 075203. doi: 10.7498/aps.70.20201368
    [4] 邹雄, 漆小波, 张涛先, 高章帆, 黄卫星. 惯性约束聚变靶丸内杂质气体抽空流洗过程的数值模拟. 物理学报, 2021, 70(7): 075207. doi: 10.7498/aps.70.20201491
    [5] 叶欣, 单彦广. 疏水表面振动液滴模态演化与流场结构的数值模拟. 物理学报, 2021, 70(14): 144701. doi: 10.7498/aps.70.20210161
    [6] 庄晓如, 徐心海, 杨智, 赵延兴, 余鹏. 高温吸热管内超临界CO2传热特性的数值模拟. 物理学报, 2021, 70(3): 034401. doi: 10.7498/aps.70.20201005
    [7] 左娟莉, 杨泓, 魏炳乾, 侯精明, 张凯. 气力提升系统气液两相流数值模拟分析. 物理学报, 2020, 69(6): 064705. doi: 10.7498/aps.69.20191755
    [8] 胡兵, 郁殿龙, 刘江伟, 朱付磊, 张振方. 流固耦合声子晶体管路冲击振动特性研究. 物理学报, 2020, 69(19): 194301. doi: 10.7498/aps.69.20200414
    [9] 闫晨帅, 徐进良. 超临界压力CO2在水平圆管内流动传热数值分析. 物理学报, 2020, 69(4): 044401. doi: 10.7498/aps.69.20191513
    [10] 喻明浩. 非平衡感应耦合等离子体流场与电磁场作用机理的数值模拟. 物理学报, 2019, 68(18): 185202. doi: 10.7498/aps.68.20190865
    [11] 吴晓笛, 刘华坪, 陈浮. 基于浸入边界-多松弛时间格子玻尔兹曼通量求解法的流固耦合算法研究. 物理学报, 2017, 66(22): 224702. doi: 10.7498/aps.66.224702
    [12] 徐肖肖, 吴杨杨, 刘朝, 王开正, 叶建. 水平螺旋管内超临界CO2冷却换热的数值模拟. 物理学报, 2015, 64(5): 054401. doi: 10.7498/aps.64.054401
    [13] 蒋勇, 贺少勃, 袁晓东, 王海军, 廖威, 吕海兵, 刘春明, 向霞, 邱荣, 杨永佳, 郑万国, 祖小涛. CO2激光光栅式扫描修复熔石英表面缺陷的实验研究与数值模拟. 物理学报, 2014, 63(6): 068105. doi: 10.7498/aps.63.068105
    [14] 史冬岩, 王志凯, 张阿漫. 任意复杂流-固边界的格子Boltzmann处理方法. 物理学报, 2014, 63(7): 074703. doi: 10.7498/aps.63.074703
    [15] 殷鹏飞, 张蓉, 熊江涛, 李京龙. 搅拌摩擦焊准稳态热力耦合过程数值模拟研究. 物理学报, 2013, 62(1): 018102. doi: 10.7498/aps.62.018102
    [16] 王新鑫, 樊丁, 黄健康, 黄勇. 双钨极耦合电弧数值模拟. 物理学报, 2013, 62(22): 228101. doi: 10.7498/aps.62.228101
    [17] 孙东科, 项楠, 陈科, 倪中华. 格子玻尔兹曼方法模拟弯流道中粒子的惯性迁移行为. 物理学报, 2013, 62(2): 024703. doi: 10.7498/aps.62.024703
    [18] 孙健, 刘伟强. 疏导式结构在头锥热防护中的应用. 物理学报, 2012, 61(17): 174401. doi: 10.7498/aps.61.174401
    [19] 靳冬欢, 刘文广, 陈星, 陆启生, 赵伊君. 三股互击式喷注器及燃烧室流场的数值模拟. 物理学报, 2012, 61(6): 064206. doi: 10.7498/aps.61.064206
    [20] 刘东戎, 桑宝光, 康秀红, 李殿中. 考虑固相移动的大尺寸钢锭宏观偏析数值模拟. 物理学报, 2009, 58(13): 104-S111. doi: 10.7498/aps.58.104
计量
  • 文章访问数:  1677
  • PDF下载量:  90
  • 被引次数: 0
出版历程
  • 收稿日期:  2023-12-11
  • 修回日期:  2024-03-17
  • 上网日期:  2024-04-02
  • 刊出日期:  2024-05-20

/

返回文章
返回