搜索

x

留言板

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

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

光纤探针测量及多尺度熵鉴别超临界类沸腾传热模式

何孝天 徐进良 程怡玮

引用本文:
Citation:

光纤探针测量及多尺度熵鉴别超临界类沸腾传热模式

何孝天, 徐进良, 程怡玮

Measurements and identification of supercritical pseudo-boiling heat transfer modes based on fiber optic probes and multiscale entropy

He Xiao-Tian, Xu Jin-Liang, Cheng Yi-Wei
PDF
HTML
导出引用
  • 亚临界沸腾包括界面蒸发和气泡动力学诱导的传热, 但超临界传热是否存在类界面蒸发和类气泡传热以及两者间的转换缺少直接的实验证据. 本文进行了超临界CO2液池传热的实验研究, 压力和液池温度分别为8—10 MPa和15 ℃. 作为加热元件和感温元件, 22 mm长和70 μm直径的镍铬丝水平放置在液池中, 光纤探针垂直放置, 其顶端高于镍铬丝200 μm. 发现随热流密度或壁面过热度的持续增大, 依次发生自然对流、类界面蒸发、类蒸发-沸腾转换、类沸腾4种传热模式. 本文重点关注类界面蒸发和类沸腾传热以及两者间的转换. 在类界面蒸发模式下, 传热系数随壁面过热度增大略有下降, 光纤输出小幅/高频信号, 不存在主频, 多尺度熵大, 表征随机信号波动. 在类蒸发-沸腾转换模式下, 光纤输出大幅/低频周期信号, 存在明显主频, 多尺度熵小, 代表有序的周期性脉动传热. 在以类气泡为特征的类沸腾模式下, 光纤信号波动幅度介于类蒸发和转换模式之间, 主频不明显, 多尺度熵也介于类蒸发和转换模式之间. 研究获得了超临界类沸腾直接的实验证据, 加深了对超临界传热机理的理解, 为后续理论研究和工程应用提供了基础.
    Subcritical boiling includes interfacial evaporation and heat transfer induced by bubble dynamics. However, for supercritical heat transfer, direct experimental evidence of the existence of pseudo-evaporation and pseudo-boiling heat transfer, as well as the conversion between them is lacking. In this work, the experimental study of supercritical carbon dioxide pool heat transfer is conducted. The pressure and temperature of the cell are 8–10 MPa and 15 ℃, respectively. As heating element and temperature sensing element, a nickel-chromium alloy wire with a length of 22 mm and diameter of 70 μm is placed horizontally in the high-pressure cell. The fiber optic probe is placed vertically, with its tip 200 μm above the wire. Four heat transfer modes, i.e. natural convection, pseudo-evaporation, transition of evaporation and boiling, and pseudo-boiling, are found to occur sequentially with the increase of heat flux density or wall superheat. Natural convection occurs when the wall temperature is below the pseudo-critical temperature. This work focuses on pseudo-evaporation and pseudo-boiling heat transfer and the transition between them. In the pseudo-evaporation mode, the heat transfer coefficient decreases slightly with wall superheat increasing. The fiber outputs a high frequency signal with small amplitude, and there is no dominant frequency. The multiscale entropy is large, characterizing random signal fluctuations. In the transition of evaporation and boiling mode, the fiber outputs a large-amplitude/low-frequency periodic signal with a significant dominant frequency and small multiscale entropy, representing an ordered periodic pulsating heat transfer. In the pseudo-boiling mode characterized by bubble-like structure, the fiber signal fluctuation amplitude and multiscale entropy are between the counterparts of the first two modes i.e. natural convection mode and pseudo-evaporation mode. The dominant frequency is not obvious. The multiscale entropies in the specific case are calculated under different key parameters, such as dimensionality, time scale factor, and length of origin data. Optimal parameters are selected based on the best separation of heat transfer modes. Finally, pseudo-boiling can be distinguished from pseudo-evaporation by multiscale entropy of 0.9, and from transition of evaporation and boiling by multiscale entropy of 0.5. In this work, direct experimental evidence of supercritical-like boiling is obtained, which deepens the understanding of the supercritical heat transfer mechanism and provides a basis for theoretical studies and engineering applications in 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]

    Saeed M, Kim M H 2022 Energy 239 121868Google Scholar

    [2]

    Zhang Z, Zhou R, Ge X, Zhang J, Wu X 2020 Energy 190 116411Google Scholar

    [3]

    Bian X Y, Xuan W, Rui W, Cai J W, Hua T, Shu G Q 2022 Energy 256 124691Google Scholar

    [4]

    Cengel Y, Boles M 2005 Thermodynamics: An Engineering Approach (5th Ed.) (SI Units) (New York: McGraw-Hill) p119

    [5]

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

    [6]

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

    [7]

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

    [8]

    Ha M Y, Yoon T J, Tlusty T, Jho Y, Lee W B 2018 J. Phys. Chem. Lett. 9 1734Google Scholar

    [9]

    Kandlikar S G 1999 Handbook of Phase Change: Boiling and Condensation (1st Ed. ) (Routledge) p63

    [10]

    Nasuti F, Pizzarelli M 2021 J. Supercrit. Fluid 168 105066Google Scholar

    [11]

    Knapp K K, Sabersky R H 1966 Int. J. Heat Mass Transf. 9 41Google Scholar

    [12]

    Knapp K K 1965 Ph. D. Dissertation (California: California Institute of Technology)

    [13]

    Hahne E, Neumann R 1981 Wärme-und Stoffübertragung 15 171

    [14]

    Rousselet Y, Warrier G R, Dhir V K 2013 J. Heat Trans. 135 022501Google Scholar

    [15]

    Holman J, Rea S, Howard C 1965 Int. J. Heat Mass Transf. 8 1095Google Scholar

    [16]

    Sakurai K, Ko H S, Okamoto K, Madarame H 2001 International Conference on Nuclear Engineering, France, April 8–12, 2001 INIS-FR-759

    [17]

    Du X, Zhu X J, Yu X, Wu F Y, Qiu Q G 2022 Nucl. Eng. Des. 386 111574Google Scholar

    [18]

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

    [19]

    张海松, 朱鑫杰, 朱兵国, 徐进良, 刘欢 2020 物理学报 69 064401Google Scholar

    Zhang H S, Zhu X J, Zhu B G, Xu J L, Liu H 2020 Acta Phys. Sin. 69 064401Google Scholar

    [20]

    Maxim F, Contescu C, Boillat P, Niceno B, Karalis K, Testino A, Ludwig C 2019 Nat. Commun. 10 1Google Scholar

    [21]

    程洁, 郭亚军, 王腾, 桂淼, 刘朝辉, 随志强 2019 化工学报 70 1375

    Cheng J, Guo Y J, Wang T, Gui M, Liu C H, Sui Z Q 2019 Ciesc Journal 70 1375

    [22]

    潘艳芝, 郑永建, 王栋, 安维峥, 牛棚满 2020 工程热物理学报 41 2756

    Pan Y Z, Zheng Y J, Wang D, An W Z, Niu P M 2020 J. Engineering Thermophysics 41 2756

    [23]

    冷洁, 胡雪音, 田瑞峰 2020 核动力工程 41 74

    Leng J, Hu X Y, Tian R F 2020 Nuclear Power Engineering 41 74

    [24]

    翟路生, 金宁德 2016 物理学报 65 010501Google Scholar

    Zhai L S, Jin N D 2016 Acta Phys. Sin. 65 010501Google Scholar

    [25]

    马敏, 郭鑫, 于洁 2022 仪器仪表学报 43 110

    Ma M, Guo X, Yu J 2022 Chin. J. Sci. Instrum. 43 110

    [26]

    刘若琳, 孔明, 曹鸣, 王道档, 单良 2022 光电子·激光 33 840

    Liu R L, Kong M, Cao M, Wang D D, Shan L 2022 J. Optoelectro. Laser 33 840

    [27]

    谭万尧, 刘晓晶, 吴德操, 罗彬彬, 赵明富, 郭裕丰, 龙俊宇, 张颖 2022 光学学报 42 52

    Tan W Y, Liu X J, Wu D C, Luo B B, Zhao M F, Guo Y F, Long J Y, Zhang Y 2022 Acta Opt. Sin. 42 52

    [28]

    Avdeev M V, Konovalov A N, Bagratashvili V N, Popov V K, Tsypina S I, Sokolova M, Ke J, Poliakoff M 2004 Phys. Chem. Chem. Phys. 6 1258Google Scholar

    [29]

    Holman J P 2012 Experimental Methods for Engineers (New York: McGraw Hill) pp64–65

    [30]

    Barney R, Nourgaliev R, McCallen R, Delplanque J P 2021 Int. J. Heat Mass Transf. 178 121572Google Scholar

    [31]

    Liu H D, Liu W, Yan P G, Chen D Q, Dong K J, Qin J, Liu H Z 2022 Int. J. Heat Mass Transf. 196 123248Google Scholar

    [32]

    Xu J L, Shen S, Gan Y H, Li Y X, Zhang W, Su Q C 2005 J. Micromech Microeng 15 1344Google Scholar

    [33]

    Taylor G I 1950 Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences 201 192

    [34]

    Shannon C E 1948 Bell Sys. Tec. J. 27 379Google Scholar

    [35]

    Pincus S M 1991 Pro. Nat. Academy Sci. 88 2297Google Scholar

    [36]

    Richman J S, Moorman J R 2000 Am. J. Physiol-Heart C 278 H2039Google Scholar

    [37]

    杨大为, 赵永东, 冯辅周, 江鹏程, 丁闯 2018 兵工学报 39 1683Google Scholar

    Yang D W, Zhao Y D, Feng F Z, Jiang P C, Ding C 2018 Acta Armamentarii 39 1683Google Scholar

    [38]

    Haryadi B, Liou J J, Wei H C, Xiao M X, Wu H T, Sun C K 2018 Comput. Meth. Prog. Bio. 166 115Google Scholar

    [39]

    Silva L E V, Duque J J, Felipe J C, Murta L O, Humeau-Heurtier A 2018 Signal Process. 147 224Google Scholar

    [40]

    Li Y B, Wang X Z, Liu Z B, Liang X H, Si S B 2018 IEEE Access 6 66723Google Scholar

  • 图 1  SCF池式传热实验系统 (a) 实验台照片; (b) 实验台原理图 (1-真空泵, 2-阀门, 3-CO2气瓶, 4-增压泵, 5-减压阀, 6-压力变送器, 7-光环行器, 8-分光器, 9-激光光源, 10-光电转换电路, 11-光纤探针, 12-高压釜, 13-LED光源, 14-冷却铜管, 15-加热丝, 16-可视窗, 17-高速相机, 18-热电偶, 19-恒温槽, 20-PID控制器, 21-直流电源, 22-标准电阻, 23-同步器, 24-高速数采仪); (c) 光纤探针封装结构

    Fig. 1.  Experiment system of SCF pool heat transfer: (a) The photo of experiment system; (b) schematic diagram of the experimental system (1-vacuum pump, 2-valve, 3-CO2 cylinder, 4-piston pump, 5-reducing valve, 6-pressure transmitter, 7-optical circulator, 8-splitter, 9-laser light source, 10-optical conversion circuit, 11-fiber optic probe, 12-high pressure vessel, 13-LED light source, 14-circulating pipe, 15-heating wire, 16-viewable window, 17-high-speed camera, 18-thermocouple, 19-thermostat, 20-PID controller, 21-DC power source, 22-standard resistor, 23-synchronizer, 24-high-speed data acquisition); (c) fiber optic probe package structure.

    图 2  探针定位和标定 (a) 光纤探针和加热丝安装; (b) 探针对齐; (c) 光纤探针标定

    Fig. 2.  Positioning and calibration of fiber optic probe: (a) Installation of fiber optic probe and heating wire; (b) probe alignment; (c) fiber optic probe calibration.

    图 3  不同压力下的传热曲线和传热模式分区 (a), (b) P分别为7.955—8.048 MPa和9.928—10.035 MPa时热流密度与过热度关系; (c), (d) P分别为7.955—8.048 MPa和9.928—10.035 MPa时传热系数与过热度关系

    Fig. 3.  QTw and hTw curves and heat transfer modes under different pressure: (a), (b) Heat flux density versus superheat for P = 7.955–8.048 MPa and 9.928–10.035 MPa, respectively; (c), (d) heat transfer coefficient versus superheat for P = 7.955– 8.048 MPa and 9.928–10.035 MPa, respectively.

    图 4  不同传热模式的图像 (a) NC模式; (b) PE模式; (c) TEB模式; (d) PB模式

    Fig. 4.  Images of different heat transfer modes: (a) NC mode; (b) PE mode; (c) TEB mode; (d) PB mode.

    图 5  探针温度信号和频谱 (a) PE模式, q = 174.1 kW/m2; (b) TEB模式, q = 359.8 kW/m2; (c) PB模式, q = 409.2 kW/m2

    Fig. 5.  Temperature signal and spectrum of probe: (a) PE mode, q = 174.1 kW/m2; (b) TEB mode, q = 359.8 kW/m2; (c) PB mode, q = 409.2 kW/m2.

    图 6  基于主频比α 的传热模式识别 (a) P = 7.955—8.048 MPa; (b) P = 9.928—10.035 MPa

    Fig. 6.  Heat transfer mode identification based on the main frequency ratio α: (a) P = 7.955–8.048 MPa; (b) P = 9.928–10.035 MPa

    图 7  多尺度熵关键参数最优选取 (a) p = 0.2; (b) N = 2000; (c) τ = 10

    Fig. 7.  Optimal selection of key parameters of multiscale entropy: (a) p = 0.2; (b) N = 2000; (c) τ = 10.

    图 8  基于多尺度熵的传热模式识别 (a) P = 7.955—8.048 MPa; (b) P = 9.928—10.035 MPa

    Fig. 8.  Heat transfer modes identification based on multiscale entropy: (a) P = 7.955–8.048 MPa; (b) P = 9.928–10.035 MPa.

  • [1]

    Saeed M, Kim M H 2022 Energy 239 121868Google Scholar

    [2]

    Zhang Z, Zhou R, Ge X, Zhang J, Wu X 2020 Energy 190 116411Google Scholar

    [3]

    Bian X Y, Xuan W, Rui W, Cai J W, Hua T, Shu G Q 2022 Energy 256 124691Google Scholar

    [4]

    Cengel Y, Boles M 2005 Thermodynamics: An Engineering Approach (5th Ed.) (SI Units) (New York: McGraw-Hill) p119

    [5]

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

    [6]

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

    [7]

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

    [8]

    Ha M Y, Yoon T J, Tlusty T, Jho Y, Lee W B 2018 J. Phys. Chem. Lett. 9 1734Google Scholar

    [9]

    Kandlikar S G 1999 Handbook of Phase Change: Boiling and Condensation (1st Ed. ) (Routledge) p63

    [10]

    Nasuti F, Pizzarelli M 2021 J. Supercrit. Fluid 168 105066Google Scholar

    [11]

    Knapp K K, Sabersky R H 1966 Int. J. Heat Mass Transf. 9 41Google Scholar

    [12]

    Knapp K K 1965 Ph. D. Dissertation (California: California Institute of Technology)

    [13]

    Hahne E, Neumann R 1981 Wärme-und Stoffübertragung 15 171

    [14]

    Rousselet Y, Warrier G R, Dhir V K 2013 J. Heat Trans. 135 022501Google Scholar

    [15]

    Holman J, Rea S, Howard C 1965 Int. J. Heat Mass Transf. 8 1095Google Scholar

    [16]

    Sakurai K, Ko H S, Okamoto K, Madarame H 2001 International Conference on Nuclear Engineering, France, April 8–12, 2001 INIS-FR-759

    [17]

    Du X, Zhu X J, Yu X, Wu F Y, Qiu Q G 2022 Nucl. Eng. Des. 386 111574Google Scholar

    [18]

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

    [19]

    张海松, 朱鑫杰, 朱兵国, 徐进良, 刘欢 2020 物理学报 69 064401Google Scholar

    Zhang H S, Zhu X J, Zhu B G, Xu J L, Liu H 2020 Acta Phys. Sin. 69 064401Google Scholar

    [20]

    Maxim F, Contescu C, Boillat P, Niceno B, Karalis K, Testino A, Ludwig C 2019 Nat. Commun. 10 1Google Scholar

    [21]

    程洁, 郭亚军, 王腾, 桂淼, 刘朝辉, 随志强 2019 化工学报 70 1375

    Cheng J, Guo Y J, Wang T, Gui M, Liu C H, Sui Z Q 2019 Ciesc Journal 70 1375

    [22]

    潘艳芝, 郑永建, 王栋, 安维峥, 牛棚满 2020 工程热物理学报 41 2756

    Pan Y Z, Zheng Y J, Wang D, An W Z, Niu P M 2020 J. Engineering Thermophysics 41 2756

    [23]

    冷洁, 胡雪音, 田瑞峰 2020 核动力工程 41 74

    Leng J, Hu X Y, Tian R F 2020 Nuclear Power Engineering 41 74

    [24]

    翟路生, 金宁德 2016 物理学报 65 010501Google Scholar

    Zhai L S, Jin N D 2016 Acta Phys. Sin. 65 010501Google Scholar

    [25]

    马敏, 郭鑫, 于洁 2022 仪器仪表学报 43 110

    Ma M, Guo X, Yu J 2022 Chin. J. Sci. Instrum. 43 110

    [26]

    刘若琳, 孔明, 曹鸣, 王道档, 单良 2022 光电子·激光 33 840

    Liu R L, Kong M, Cao M, Wang D D, Shan L 2022 J. Optoelectro. Laser 33 840

    [27]

    谭万尧, 刘晓晶, 吴德操, 罗彬彬, 赵明富, 郭裕丰, 龙俊宇, 张颖 2022 光学学报 42 52

    Tan W Y, Liu X J, Wu D C, Luo B B, Zhao M F, Guo Y F, Long J Y, Zhang Y 2022 Acta Opt. Sin. 42 52

    [28]

    Avdeev M V, Konovalov A N, Bagratashvili V N, Popov V K, Tsypina S I, Sokolova M, Ke J, Poliakoff M 2004 Phys. Chem. Chem. Phys. 6 1258Google Scholar

    [29]

    Holman J P 2012 Experimental Methods for Engineers (New York: McGraw Hill) pp64–65

    [30]

    Barney R, Nourgaliev R, McCallen R, Delplanque J P 2021 Int. J. Heat Mass Transf. 178 121572Google Scholar

    [31]

    Liu H D, Liu W, Yan P G, Chen D Q, Dong K J, Qin J, Liu H Z 2022 Int. J. Heat Mass Transf. 196 123248Google Scholar

    [32]

    Xu J L, Shen S, Gan Y H, Li Y X, Zhang W, Su Q C 2005 J. Micromech Microeng 15 1344Google Scholar

    [33]

    Taylor G I 1950 Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences 201 192

    [34]

    Shannon C E 1948 Bell Sys. Tec. J. 27 379Google Scholar

    [35]

    Pincus S M 1991 Pro. Nat. Academy Sci. 88 2297Google Scholar

    [36]

    Richman J S, Moorman J R 2000 Am. J. Physiol-Heart C 278 H2039Google Scholar

    [37]

    杨大为, 赵永东, 冯辅周, 江鹏程, 丁闯 2018 兵工学报 39 1683Google Scholar

    Yang D W, Zhao Y D, Feng F Z, Jiang P C, Ding C 2018 Acta Armamentarii 39 1683Google Scholar

    [38]

    Haryadi B, Liou J J, Wei H C, Xiao M X, Wu H T, Sun C K 2018 Comput. Meth. Prog. Bio. 166 115Google Scholar

    [39]

    Silva L E V, Duque J J, Felipe J C, Murta L O, Humeau-Heurtier A 2018 Signal Process. 147 224Google Scholar

    [40]

    Li Y B, Wang X Z, Liu Z B, Liang X H, Si S B 2018 IEEE Access 6 66723Google Scholar

  • [1] 孙辉, 刘婧楠, 章立新, 杨其国, 高明. 超临界二氧化碳类液-类气区边界线数值分析. 物理学报, 2022, 71(4): 040201. doi: 10.7498/aps.71.20211464
    [2] 孙辉, 刘婧楠, 章立新, 杨其国, 高明. 超临界CO2类液-类气区边界线数值分析. 物理学报, 2021, (): . doi: 10.7498/aps.70.20211464
    [3] 张海松, 徐进良, 朱鑫杰. 基于拟沸腾理论的超临界CO2管内传热恶化量纲分析. 物理学报, 2021, 70(4): 044401. doi: 10.7498/aps.70.20201546
    [4] 张夫一, 葛曼玲, 郭志彤, 谢冲, 杨泽坤, 宋子博. 静息态功能磁共振成像评估健康老年人认知行为的多尺度熵模型研究. 物理学报, 2020, 69(10): 108703. doi: 10.7498/aps.69.20200050
    [5] 张夫一, 葛曼玲, 郭志彤, 谢冲, 杨泽坤, 宋子博. 静息态功能磁共振成像评估健康老年人认知行为 的多尺度熵模型研究. 物理学报, 2020, (): 008700. doi: 10.7498/aps.69.20200051
    [6] 黄泽徽, 李亚安, 陈哲, 刘恋. 基于多尺度熵的Duffing混沌系统阈值确定方法. 物理学报, 2020, 69(16): 160501. doi: 10.7498/aps.69.20191642
    [7] 王艳, 徐进良, 李文, 刘欢. 超临界Lennard-Jones流体结构特性分子动力学研究. 物理学报, 2020, 69(7): 070201. doi: 10.7498/aps.69.20191591
    [8] 盛子城, 王腾, 周桂耀, 夏长明, 刘建涛, 李波瑶, 樊海霞, 陈云, 侯峙云. 基于空芯微结构光纤拉曼探针的实验研究. 物理学报, 2018, 67(18): 184211. doi: 10.7498/aps.67.20180684
    [9] 郭旭东, 唐军, 刘文耀, 郭浩, 房国成, 赵苗苗, 王磊, 夏美晶, 刘俊. 锥柱型光纤探针在表面增强拉曼散射方面的应用. 物理学报, 2017, 66(4): 044208. doi: 10.7498/aps.66.044208
    [10] 严雪过, 沈毅, 潘聪, 李鹏, 丁志华. 基于拉锥结构的全光纤型内窥OCT探针研究. 物理学报, 2016, 65(2): 024201. doi: 10.7498/aps.65.024201
    [11] 刘桂媛, 宋洪胜, 张宁玉, 程传福. 飞秒激光在锥形镀膜探针传输中相位奇异的研究. 物理学报, 2015, 64(2): 024203. doi: 10.7498/aps.64.024203
    [12] 苟竞, 刘俊勇, 魏震波, Gareth Taylor, 刘友波. 基于多尺度熵的电力能量流复杂性分析. 物理学报, 2014, 63(20): 208402. doi: 10.7498/aps.63.208402
    [13] 向郑涛, 陈宇峰, 李昱瑾, 熊励. 基于多尺度熵的交通流复杂性分析. 物理学报, 2014, 63(3): 038903. doi: 10.7498/aps.63.038903
    [14] 李勇, 刘锦超, 芦鹏飞, 杨向东. 从常温常压到超临界乙醇的分子动力学模拟. 物理学报, 2010, 59(7): 4880-4887. doi: 10.7498/aps.59.4880
    [15] 刘桂媛, 滕树云, 程传福, 宋洪胜, 刘曼. 锥形镀膜光纤探针中飞秒激光脉冲的传输. 物理学报, 2009, 58(11): 7613-7620. doi: 10.7498/aps.58.7613
    [16] 郑桂波, 金宁德. 两相流流型多尺度熵及动力学特性分析. 物理学报, 2009, 58(7): 4485-4492. doi: 10.7498/aps.58.4485
    [17] 罗奔毅, 卢义刚. 超临界点附近二氧化碳流体的声速. 物理学报, 2008, 57(7): 4397-4401. doi: 10.7498/aps.57.4397
    [18] 何 亮, 杜 磊, 庄奕琪, 李伟华, 陈建平. 金属互连电迁移噪声的多尺度熵复杂度分析. 物理学报, 2008, 57(10): 6545-6550. doi: 10.7498/aps.57.6545
    [19] 徐 凯, 潘 石, 吴世法, 孙 伟, 李银丽. 热拉伸和化学腐蚀相结合制备弯曲光纤探针. 物理学报, 2003, 52(5): 1190-1195. doi: 10.7498/aps.52.1190
    [20] 张立源, 李永贵, 王潜. 扫描近场红外显微镜光纤探针的腐蚀制法. 物理学报, 2001, 50(12): 2322-2326. doi: 10.7498/aps.50.2322
计量
  • 文章访问数:  383
  • PDF下载量:  16
  • 被引次数: 0
出版历程
  • 收稿日期:  2022-08-24
  • 修回日期:  2022-09-28
  • 上网日期:  2023-01-05
  • 刊出日期:  2023-03-05

光纤探针测量及多尺度熵鉴别超临界类沸腾传热模式

  • 1. 华北电力大学, 低品位能源多相流与传热北京市重点实验室, 北京 102206
  • 2. 华北电力大学, 电站能量传递转化与系统教育部重点实验室, 北京 102206
  • 通信作者: 徐进良, xjl@ncepu.edu.cn
    基金项目: 国家自然科学基金重点项目(批准号: 52130608)和国家自然科学基金创新研究群体科学基金(批准号: 51821004)资助的课题

摘要: 亚临界沸腾包括界面蒸发和气泡动力学诱导的传热, 但超临界传热是否存在类界面蒸发和类气泡传热以及两者间的转换缺少直接的实验证据. 本文进行了超临界CO2液池传热的实验研究, 压力和液池温度分别为8—10 MPa和15 ℃. 作为加热元件和感温元件, 22 mm长和70 μm直径的镍铬丝水平放置在液池中, 光纤探针垂直放置, 其顶端高于镍铬丝200 μm. 发现随热流密度或壁面过热度的持续增大, 依次发生自然对流、类界面蒸发、类蒸发-沸腾转换、类沸腾4种传热模式. 本文重点关注类界面蒸发和类沸腾传热以及两者间的转换. 在类界面蒸发模式下, 传热系数随壁面过热度增大略有下降, 光纤输出小幅/高频信号, 不存在主频, 多尺度熵大, 表征随机信号波动. 在类蒸发-沸腾转换模式下, 光纤输出大幅/低频周期信号, 存在明显主频, 多尺度熵小, 代表有序的周期性脉动传热. 在以类气泡为特征的类沸腾模式下, 光纤信号波动幅度介于类蒸发和转换模式之间, 主频不明显, 多尺度熵也介于类蒸发和转换模式之间. 研究获得了超临界类沸腾直接的实验证据, 加深了对超临界传热机理的理解, 为后续理论研究和工程应用提供了基础.

English Abstract

    • 超临界流体(supercritical fluid, SCF)是现代发电动力循环的主要工作介质, 由于压力、温度高, 可有效提升系统循环效率, 提高系统紧凑性[1], 如朗肯循环中使用超临界水[2]、布雷顿循环中的超临界二氧化碳(supercritical carbon dioxide, sCO2)[3]. 经典教科书将SCF视为单相流体[4], 依据浮升力和加速效应提出关联式, 但不能准确预测传热恶化和强化等异常现象[5]. 现代物理学家通过X射线散射[6]和分子动力学模拟[7,8]等方法发现SCF有类似两相的物质结构, 在Widom线两侧可区分为类气和类液两种流体. 超临界和亚临界流体的这种相似性提供了研究SCF传热的新角度.

      亚临界流体吸热相变时有蒸发和沸腾两种模式, 二者定义清晰、互有异同. 蒸发是在已有相界面上从液体到气体的质量传递; 沸腾是在固体壁面或过热液体中核化并形成气泡, 液气相变产生新的两相界面[9]. 在超临界领域, 20世纪60年代就已提出“类沸腾”概念: 传热实验中, 类沸腾概念仅代表部分工况下的传热异常和性能相似性; 热力学中, 类液体升温成为类气体的过程被统称为类沸腾[10]. 传热模式的模糊定义阻碍了超临界类沸腾相关研究.

      在过去的研究中, SCF传热表征手段有限, 主要通过视窗观察流体现象与传热数据之间的联系. Knapp和Sabersky[11]使用镍铬合金丝在高压容器中加热sCO2, 通过纹影法发现, 产生类气泡时qwT曲线斜率变陡, 与亚临界沸腾传热起始点相似. Knapp[12]试验了多种加热方式, 只在水平金属丝上拍摄到类气泡, 金属条或垂直金属丝上没有产生. Hahne和Neumann[13]的实验表明高压釜中形成类气泡与金属丝材质有关, 铂丝上不产生类气泡, 也没有传热强化. 2013年Rousselet等[14]总结了前人实验结果, 当Tw < Tpc时传热类似单相自然对流, Tw > Tpc时类似亚临界膜态沸腾, 其中Tw为壁温, Tpc为拟临界温度. SCF管内强制对流实验研究较多, 但针对类两相分布的很少. Holman等[15]观察到超临界氟利昂12在玻璃管中加热时有气态轨迹, 沿程壁温先升后降. Sakurai等[16]用阴影法拍摄矩形管中单侧加热的sCO2流场, 随热流升高, 壁面出现密度波的黑色阴影, 传热系数下降. Du等[17]设计了类似实验段并直接拍摄流场, 传热恶化时壁面有白色阴影. 本课题组[18,19]通过管外壁温数据推测管内流体状态, 提出类气膜的分布状态是造成传热恶化的主要原因.

      总之, 过去的研究已经发现SCF传热和类两相分布之间存在关联, 如类气泡造成传热强化, 类气膜引起传热恶化, 但只能通过可视化实验和传热数据确认传热模式, 手段单一、缺乏定量判据. Maxim等[20]发展了SCF类沸腾中子射线测量技术, 可以有效分辨类两相分布, 但需在通道中填充多孔碳纤维, 成本高且难以推广. 亚临界流体相分布测量技术相当成熟, 发展出射线法[21,22]、超声波[23]、电容/电导[24,25]和激光光学[26,27]等多种测量技术, 有效识别流动传热模式, 如泡状流、段塞流、环状流等. 然而, 这些技术大多依赖气液两相物性的突变 (密度、声速、电导率、折射率等), 不能直接用于测量物性连续变化的SCF.

      本文在高压容器中用镍铬金属丝加热sCO2, 采用改进的光纤探针测量sCO2温度场连续波动, 研究不同传热模式的现象和转换机理, 实验参数覆盖P = 8—10 MPa, q = 0—1800 kW/m2, Tb = 15 ℃. 随热流密度升高, 传热模式从自然对流转换为类蒸发、随后过渡为内沸腾, 对应于传热系数先升后降, 最后再升高. 引入频谱和熵分析光纤探针信号特征, 提出基于多尺度熵的SCF传热模式定量判据, 获得了超临界类沸腾直接的实验证据.

    • SCF池式传热系统和测量原理如图1所示, 主要包括高压釜、温度压力控制管路、铂丝电路、光纤探针系统、高速摄影及数采. 实验工质选用纯度99.999%的CO2, 临界点为Pcr = 7.377 MPa, Tcr = 31.2 ℃. 工质被气动活塞泵增压注入高压釜, 由入口减压阀和排气阀控制釜内压力P. 高压釜为304不锈钢材质筒体, 内径100 mm, 长为140 mm, 两端面有直径30 mm蓝宝石玻璃视窗. 通过釜中纯铜盘管循环水温PID控制池温Tb, 本文中Tb均为15 ℃.

      图  1  SCF池式传热实验系统 (a) 实验台照片; (b) 实验台原理图 (1-真空泵, 2-阀门, 3-CO2气瓶, 4-增压泵, 5-减压阀, 6-压力变送器, 7-光环行器, 8-分光器, 9-激光光源, 10-光电转换电路, 11-光纤探针, 12-高压釜, 13-LED光源, 14-冷却铜管, 15-加热丝, 16-可视窗, 17-高速相机, 18-热电偶, 19-恒温槽, 20-PID控制器, 21-直流电源, 22-标准电阻, 23-同步器, 24-高速数采仪); (c) 光纤探针封装结构

      Figure 1.  Experiment system of SCF pool heat transfer: (a) The photo of experiment system; (b) schematic diagram of the experimental system (1-vacuum pump, 2-valve, 3-CO2 cylinder, 4-piston pump, 5-reducing valve, 6-pressure transmitter, 7-optical circulator, 8-splitter, 9-laser light source, 10-optical conversion circuit, 11-fiber optic probe, 12-high pressure vessel, 13-LED light source, 14-circulating pipe, 15-heating wire, 16-viewable window, 17-high-speed camera, 18-thermocouple, 19-thermostat, 20-PID controller, 21-DC power source, 22-standard resistor, 23-synchronizer, 24-high-speed data acquisition); (c) fiber optic probe package structure.

      池式传热加热元件为镍铬合金丝, 牌号Cr20Ni35, 直径70 μm, 长22 mm. 镍铬丝、标准电阻和直流电源(Itech 6132B)组成串联电路, 控制加热功率. 镍铬丝同时也是测温元件, 实验前标定电阻率和温度的关系, 实验中测量电阻计算加热丝壁温. 镍铬丝水平悬挂在高压釜中心, 高速相机透过视窗拍摄现象(Keyence VW-9000), 拍摄帧率4000 frames/s, 分辨率1.9 μm.

    • Avdeev等[28]首先制备适合超临界环境下的反射式光纤探针, 命名为Fibre Optic Reflectometer (FOR), 用于测量超临界纯流体和混合物的密度. 本文将FOR用于SCF传热过程, 测量流场密度和温度变化, 获取传热过程瞬态特征. 光纤探针测量原理如图1(b)所示. 光源发出的1310 nm激光被分光器分为测量光束和参考光束. 参考光束经光电二极管和跨阻放大器直接转换为参考电压信号Uref. 测量光束经光环行器进入探针, 在探针端面反射, 入射光强和反射光强分别为I0IR. 反射光沿原光路返回光环行器, 传递至光电回路转换为测量电压信号Usig. 光纤端面入射光强I0和反射光强IR的比值与光纤折射率n0和流体折射率n有关, 且与电压信号成正比[28]:

      $ \frac{{{I_{\text{R}}}}}{{{I_0}}} = {\left( {\frac{{n - {n_0}}}{{n + {n_0}}}} \right)^2} = \kappa \frac{{{U_{{\text{sig}}}}}}{{{U_{{\text{ref}}}}}} , $

      其中n0为光纤材料SiO2的折射率, n0 = 1.45, 比例系数κ与光纤性质和沿程损失有关, 与测量流体无关. 在空气中测得的电压信号Usigair/Urefair满足关系

      $ {\left( {\frac{{{n_{{\text{air}}}} - {n_0}}}{{{n_{{\text{air}}}} + {n_0}}}} \right)^2} = \kappa \frac{{{U_{{\text{sigair}}}}}}{{{U_{{\text{refair}}}}}} , $

      (1)式和(2)式可消去κ, 得到光纤端面流体折射率n[28]

      $ n = {n_0}\frac{{\sqrt {\dfrac{{{U_{{\text{sigair}}}}}}{{{U_{{\text{refair}}}}}}\dfrac{{{U_{{\text{ref}}}}}}{{{U_{{\text{sig}}}}}}} \left( {{n_0} + {n_{{\text{air}}}}} \right) - \left( {{n_0} - {n_{{\text{air}}}}} \right)}}{{\sqrt {\dfrac{{{U_{{\text{sigair}}}}}}{{{U_{{\text{refair}}}}}}\dfrac{{{U_{{\text{ref}}}}}}{{{U_{{\text{sig}}}}}}} \left( {{n_0} + {n_{{\text{air}}}}} \right) + \left( {{n_0} - {n_{{\text{air}}}}} \right)}} , $

      其中nair = 1.00027 (20 ℃, 1 atm, 1 atm = 1.01 × 105 Pa). 流体折射率n和密度ρ之间满足Lorentz-Lorenz方程[28]

      $ A\rho = \frac{{{n^2} - 1}}{{{n^2} + 2}}. $

      对于非极性分子CO2, A = 1.42×10–4 m3/kg, 是与温度无关的常数. 根据sCO2的密度ρ和压力P, 在REFPROP软件中可查得流体温度Tg.

      光纤探针密封结构如图1(c)所示. 探针本体为普通通信光纤, 端面直径125 μm, 纤芯直径9 μm, 用密封胶封装在长200 mm、直径0.6 mm不锈钢毛细管中. 毛细管和耐压套管间填充锥形聚四氟, 由顶部螺栓挤压密封. 耐压套管外加工有齿距350 μm超细螺纹, 配合角度盘纵向调节探针高度. 实验中探针位于加热丝上方200 μm处, 如图2(a)所示. 装配时将光纤光源替换为红色可见光, 调节加热丝下方手动位移台, 当加热丝上光斑最亮时光纤位置对正, 如图2(b)所示.

      图  2  探针定位和标定 (a) 光纤探针和加热丝安装; (b) 探针对齐; (c) 光纤探针标定

      Figure 2.  Positioning and calibration of fiber optic probe: (a) Installation of fiber optic probe and heating wire; (b) probe alignment; (c) fiber optic probe calibration.

      为验证光纤探针精度, 在高压釜中与铠装热电偶同时测量sCO2温度, 测量压力范围8—12 MPa, 温度5—40 ℃. 结果如图2(c)所示, 二者平均相对误差eA = 3.57%, 平均相对误差eR = 4.67%, 平均标准差eS = 7.49%, 光纤探针可准确测量sCO2流体温度. 与普通热电偶相比, 光纤没有热响应引起的时间延迟, 可实时反应温度波动.

    • 本文在8 MPa和10 MPa压力下开展sCO2传热实验, 实测压力P范围分别为7.955—8.048 MPa和9.928—10.035 MPa, 后文以8 MPa和10 MPa指代. 池温为(15±0.1) ℃, 热流密度范围0—1800 kW/m2. 恒河数采(Yokogawa DL750)测量电源电压Ut和标准电阻电压Us, 加热丝热流密度q和电阻率$ {\overline \rho _{\text{w}}} $分别为

      $ q = \frac{1}{{{\text{π }}dL}}\frac{{{U_{\text{t}}} - {U_{\text{s}}}}}{{{R_{\text{s}}}}}{U_{\text{s}}} \text{, } $

      $ {\overline \rho _{\text{w}}} = \frac{{{\text{π }}{d^2}}}{{4L}}\frac{{{U_{\text{t}}} - {U_{\text{s}}}}}{{{U_{\text{s}}}}}{R_{\text{s}}} \text{, } $

      其中dL分别为加热丝直径和长度. 提前标定的镍铬丝电阻率和温度关系为

      $ {\overline \rho _{{\text{CrNi}}}} = 1.1061 \times {10^{ - 6}} \times \left( {1 + {{4.2155}^{ - 4}}{T_{\text{w}}}} \right) \text{, } $

      根据壁温Tw计算传热系数h:

      $ h = {q}/({{{T_{\text{w}}} - {T_{\text{b}}}}}) . $

      热电偶经过精密水银温度计校准, 测量Tb不确定度为0.1 ℃, 压力和电压不确定度分别为0.075%和0.5%. 不能直接测量的参数测量不确定度按下列公式传递[29]:

      $ {\text{δ}}r = \sqrt {\sum\limits_{i = 1}^N {{{\left( {\frac{{\partial r}}{{\partial {x_i}}}{\text{δ}} {x_i}} \right)}^2}} } \text{, } $

      其中xi为直接测得量的不确定度, 则热流密度、壁温和传热系数不确定度分别为1.74%, 3.71%和4.84%.

    • SCF在指定压力下定压比热容cp存在最大值, 此时的温度为类临界温度Tpc. 在Tpc附近, 流体密度ρ、黏度μ和导热系数λ都随温度发生剧烈变化. Maxim等[20]和Barney等[30]认为这种显著的物性变化类似亚临界流体在沸点处的相变, 当T > Tpc时, 流体性质类似气体; 反之, T < Tpc时, 性质类似液体. 定义加热丝壁面过热度$ \Delta {T_{\text{w}}} = {T_{\text{w}}} - {T_{{\text{pc}}}} $, 加热丝传热规律和现象见图3图4. 在不同压力下, 随ΔTw增高传热系数h都先增后减, 最后在振荡中保持上升趋势. 在宽广的工况范围内存在4种传热模式: 自然对流(natural convection, NC), 类蒸发(pseudo evaporation, PE), 类蒸发-沸腾转换(transition of evaporation and boiling, TEB), 以及类沸腾(pseudo boiling, PB).

      图  3  不同压力下的传热曲线和传热模式分区 (a), (b) P分别为7.955—8.048 MPa和9.928—10.035 MPa时热流密度与过热度关系; (c), (d) P分别为7.955—8.048 MPa和9.928—10.035 MPa时传热系数与过热度关系

      Figure 3.  QTw and hTw curves and heat transfer modes under different pressure: (a), (b) Heat flux density versus superheat for P = 7.955–8.048 MPa and 9.928–10.035 MPa, respectively; (c), (d) heat transfer coefficient versus superheat for P = 7.955– 8.048 MPa and 9.928–10.035 MPa, respectively.

      图  4  不同传热模式的图像 (a) NC模式; (b) PE模式; (c) TEB模式; (d) PB模式

      Figure 4.  Images of different heat transfer modes: (a) NC mode; (b) PE mode; (c) TEB mode; (d) PB mode.

    • 在ΔTw < 0时, sCO2近似单相液体, 物性没有发生显著变化, 因此传热规律与亚临界单相自然对流一致, h随ΔTw上升呈线性增长. NC模式高速图像如图4(a)所示, 只能看到加热丝本体, 直接光学拍摄不能反映自然对流中微弱的物性变化.

    • 当ΔTw > 0时, sCO2从类液体转变为类气体, 密度和导热系数迅速减小, 换热能力差, 传热系数h随ΔTw先下降, 之后基本维持恒定. PE模式下图像如图4(b)所示, 加热丝周围出现均匀黑色阴影, 且热流密度越大, 阴影范围越宽. 这是因为类气体密度低, 折射率远小于外围的类液体, 光线穿透时折射角度大, 在图像中形成黑色阴影.

      这种现象与亚临界膜态沸腾时的界面蒸发十分相似. 亚临界沸腾时, 高热流密度下气泡大量合并为气膜, 加热丝附近传热以气膜导热为主, 导致沸腾危机[31]. 此时液体通过界面蒸发, 传热系数与热流密度几乎无关[32]. 超临界压力下, 类液体被加热升温, 膨胀为类气体, 传热系数下降, 且受热流密度的影响降低, 类似亚临界膜态沸腾时液体蒸发, 因此称为类蒸发传热模式, 尽管超临界时没有气液界面和表面张力.

    • 当ΔTw接近150 ℃时, hTw曲线出现拐点, 传热显著强化. 从图4(c)的时间序列图像可知, TEB模式中加热丝上下振荡, 伴随出现类气泡. T = 0 ms时, 加热丝周围仍然被均匀类气膜包裹, 上方波纹状阴影是上个循环的遗留. 70 ms时, 加热丝有向下移动趋势, 类气膜开始失稳, 厚度分布不均, 部分区域阴影减弱, 可以看到加热丝本体, 类气膜的扰动影响加热丝上方尾流, 出现波纹阴影. 96.25 ms时, 加热丝向下移动至最低点, 类气膜完全失稳形成类气泡, 在加热丝下方呈半球形排布. 128.75 ms时, 加热丝从最低点向上移动, 周围类气体恢复稳定, 类气泡逐渐消散. 在一个振荡周期中, 加热丝传热在类蒸发和类沸腾之间转换, 称为蒸发沸腾过渡传热模式. 8 MPa时TEB模式范围较窄, 10 MPa时TEB模式范围宽广且难以预测, 甚至出现随热流密度增加, 壁温反而下降的现象, 与亚临界核态沸腾起始点(ONB)十分相似.

      图4(c)所示, 从0—128.75 ms的图像序列可以看出, 类气泡的产生与加热丝振动直接相关. 加热丝向下运动时, 推动底部类气体向类液体加速运动. Taylor等[33]通过线性失稳理论证明, 低密度流体向高密度流体加速运动时, 二者界面必然失稳, 不稳定性的发展速度与$\sqrt {\dfrac{{{\rho _2} - {\rho _1}}}{{{\rho _2} + {\rho _1}}}} $成正比, 其中密度${\rho _2} > {\rho _1}$, 即密度差越大传播越快. 这种不稳定性与表面张力无关, 即使超临界压力下不存在气液界面, 加速度也会导致类气和类液过渡区失稳, 形成气泡状结构. 加热丝振动有两个原因, 一是金属丝热膨胀变长, 振动无需克服金属弹性模量, 二是sCO2自然对流拖曳力足够与加热丝自重对抗. 传热模式转换时拖曳力时大时小, 因此加热丝忽上忽下. 加热丝运动和类气膜失稳相互促进, 初始的微小扰动最终发展为加热丝大幅振荡, 热边界层受到剧烈扰动, 厚度减薄, 传热强化.

    • 随过热度升高, h持续增长, 但h−ΔTw曲线斜率变缓. 此时传热过渡为类沸腾模式, 如图4(c)所示, 加热丝下侧持续产生类气泡. 由于没有表面张力束缚, 类气泡中的流体在浮升力作用下向上流动, 形成柱状流. 随着温度下降, 向上流动的类气体与周围类液体的边界逐渐模糊.

      PB模式与TEB模式的显著区别是加热丝不再上下振动, 产生的类气泡尺寸也略大于TEB模式. 这说明PB模式中类气膜失稳机理转变, 不再依赖加热丝运动. 假设加热丝初始温度沿轴向分布均匀, 对流流场随机小扰动可能引起局部壁温略高于均值. 从(7)式可知, 镍铬金属温度越高电阻越大, 在恒压串联电路中局部电阻Ri与局部功率Qi的关系为

      $ {Q_i} = \frac{{{R_i}U_{\text{t}}^2}}{{{{\left( {{R_{\text{r}}} + {R_i}} \right)}^2}}} \text{, } $

      其中Ut为电路总电压, Rr为标准电阻和加热丝其余部分电阻之和, 均为定值. 对Ri求导:

      $ \frac{{\partial Q}}{{\partial {R_i}}} = U_{\text{t}}^2\frac{{{R_{\text{r}}} - {R_i}}}{{{{\left( {{R_{\text{r}}} + {R_i}} \right)}^3}}} . $

      由于Rr $\gg $ Ri, 导数大于0, 因此局部温度越高发热量越大, 反之温度越低发热越小, 这导致加热丝表面温度总是趋于不均匀. 另一方面, 金属导热系数较大, 热量沿轴线传递, 会使温度分布趋于均匀. 加热丝上这两个因素相互竞争, 当热流密度较小时, 如NC, PE和TEB模式, 金属导热占主导, 温度分布均匀; 当热流密度较大时, 如PB模式, 非均匀热源占主导, 使加热丝温度分布不均. 发热强的位置类液向类气转化量更大, 聚集形成类气泡. 加热丝实际的温度分布有待未来测量证实.

      以上实验结果和分析表明, 尽管SCF没有界面和张力, 传热过程依然呈现多种截然不同的模式, 与亚临界流体自然对流向核态沸腾和膜态沸腾的转变类似.

    • 图3可以发现, NC模式仅在热流密度很小、过热度为负的条件下出现, 提高热流很快进入PE模式, 因此本文主要关注PE, TEB和PB模式流型识别, 特征工况的光纤温度信号如图5所示. PE模式加热丝表面被类气膜覆盖, 加热丝上方对流流动稳定, 探针测量的温度波动很小, 如图5(a)所示. TEB模式加热丝周期振荡, 与探针距离时远时近, 因此温度也大幅度周期波动, 如图5(b)所示. PB模式加热丝不再振动, 类气泡上方柱状流掠过探针针头, 测量温度出现随机波动, 如图5(c)所示.

      图  5  探针温度信号和频谱 (a) PE模式, q = 174.1 kW/m2; (b) TEB模式, q = 359.8 kW/m2; (c) PB模式, q = 409.2 kW/m2

      Figure 5.  Temperature signal and spectrum of probe: (a) PE mode, q = 174.1 kW/m2; (b) TEB mode, q = 359.8 kW/m2; (c) PB mode, q = 409.2 kW/m2.

      三种传热模式的探针信号波动性质有较明显区别, 进行快速傅里叶变换(FFT)后幅频特性如图5右栏所示. PE模式幅值较小, 频率分布均匀, 没有明显主频. TEB模式有明显周期特征, 主频为5.67 Hz, 幅值显著突出. PB模式主频为0.67 Hz, 频率升高幅值逐渐降低, 没有明显突跳. 可以发现, 三种模式主频fm突出程度不同, 为定量表征主频fm在频谱中的占比, 定义主频比α

      $ \alpha = \frac{{{A^2}\left( {{f_{m} }} \right)}}{{\displaystyle\sum\limits_{i = 1}^n {{A^2}\left( {{f_i}} \right)} }}, $

      其中, A(fi)为频率fi对应的幅值.

      8 MPa和10 MPa各工况点探针信号频谱计算的α值, 如图6所示. TEB模式α值均大于0.3, PB模式整体大于0.1, PE模式略小于0.1. 主频比α可以有效区分TEB模式, 但PE和PB的分界线相对模糊. 这说明PE和PB模式中温度波动随机性都比较强, PB的周期性并不突出, 因此仅采用主频比不能有效区分三种传热模式.

      图  6  基于主频比α 的传热模式识别 (a) P = 7.955—8.048 MPa; (b) P = 9.928—10.035 MPa

      Figure 6.  Heat transfer mode identification based on the main frequency ratio α: (a) P = 7.955–8.048 MPa; (b) P = 9.928–10.035 MPa

    • 在热力学理论中, 熵作为热力系统混乱程度的衡量参数, 由克劳修斯在19世纪提出. 1948年, Shannon[34]首次提出了信息熵的概念, 利用时间序列中不同状态的概率分布来评估其不确定性. 如果不同状态的概率值相似, 则很难确定未来状态, 时间序列不确定性较大, 具有较大熵值; 如果某状态的概率远大于其他状态, 则说明信号可预测性强, 时间序列较为有序, 具有较小熵值. 1991年, 受信息熵的启发, Pincus[35]引入了近似熵来量化时间序列的不规则性和自相似性. 然而, 近似熵对数据长度有着较高的要求, 如果数据长度很短, 则得到的值通常小于实际值. Richman和Moorman[36]于2000年提出了样本熵, 对数据长度要求较低且一致性更强. 样本熵的算法基于单时间尺度, 而忽略了时间序列中可能存在不同的时间尺度. 由于多尺度熵很好地考虑了潜在的多尺度属性, 在设备故障识别[37]、医疗诊断[38]和图形处理[39]等方面得到广泛应用, 其基本计算步骤如下.

      对长度为N的原始信号{x(i)}做粗粒化处理, 按下式构造出长度为N/τ的新信号{y(j)}:

      $ y(j) = \frac{1}{\tau }\sum\limits_{i = 1 + (j - 1)\tau }^{j\tau } {x(i)} ,{\text{ }}j = 1,{\text{ }}2,{\text{ }} \cdots ,{\text{ }}\frac{N}{\tau }\text{, } $

      其中, τ 为粗粒化处理的尺度因子. 定义容许偏差r,

      $ r = p \times {\text{STD}}\left( {\left\{ y \right\}} \right) \text{, } $

      其中STD({y})为信号{y(j)}的标准差, p为常系数. 使用新信号重构一组m维空间向量Y(1), Y(2), ··· , Y(N/τ – m + 1),

      $\begin{split} & {\boldsymbol{Y}}(i) = \left\{ {{y}(i),y(i + 1), \cdots ,y(i + m - 1)} \right\},\\ & ~~~~~~~~~~~~~~~~~ i = 1,2, \cdots ,{N \mathord{\left/{\vphantom {N \tau }} \right. } \tau } - m + 1.\end{split}$

      定义向量Y(i)和Y(j)之间的距离为两向量所有对应元素之差绝对值的最大值, 表示为

      $\begin{split} & d\left[ {Y(i),Y(j)} \right] = \max \left( {\left| {y(i + k) - y(j + k)} \right|} \right),\\ & ~~~~~~~~~~~~~~~ k = 0,\;1,\; \cdots ,\; m - 1 .\end{split} $

      按照(16)计算向量Y(i)与其余所有向量的距离, 并统计其中距离小于容许偏差r的数目, 记为Bi, 并计算其与距离总数之比

      $ B_i^m\left( r \right) = {{{B_i}}}/({{{N/\tau } - m}}). $

      对所有i做上述计算, 并求平均值

      $ {B^m}\left( r \right) = \frac{1}{{{N \mathord{\left/ {\vphantom {N \tau }} \right. } \tau } - m + 1}}\sum\limits_{i = 1}^{{N \mathord{\left/ {\vphantom {N \tau }} \right. } \tau } - m + 1} {B_i^m\left( r \right)} . $

      将维数增加至m+1, 重复上述步骤, 得到

      $ {B^{m + 1}}\left( r \right) = \frac{1}{{{N \mathord{\left/ {\vphantom {N \tau }} \right. } \tau } - m}}\sum\limits_{i = 1}^{{N \mathord{\left/ {\vphantom {N \tau }} \right. } \tau } - m} {B_i^{m + 1}\left( r \right)} . $

      可得对应τ尺度的多尺度熵MSE(τ)

      $ {\text{MSE}}\left( \tau \right) = - \ln \frac{{{B^{m + 1}}\left( r \right)}}{{{B^m}\left( r \right)}} . $

    • 由多尺度熵的计算过程可知, 计算结果受4个关键参数影响: 重构向量组时的维数m、常系数p、原始数据的长度N、对原始信号做粗粒化处理时的尺度因子τ. 需要分析参数的最佳选值. 现有研究表明, 取m = 1或m = 2时, 所得结果对数据长度的要求较低, 且统计特征比较合理[40], 本文选用m = 2进行计算.

      针对图5三种典型光纤信号计算MSE值, 比较不同p值对计算结果的影响. 结果见图7(a), 长度N = 2000, τ = 10. 当p值较低时, PE与PB模式曲线出现交叉, 区分效果较差; 当p = 0.2时, 三信号的MSE值区分较好; 而p取值过大时, PB和TEB相互靠近, 且可能造成系统信息丢失. 因此本文选用p = 0.2.

      图  7  多尺度熵关键参数最优选取 (a) p = 0.2; (b) N = 2000; (c) τ = 10

      Figure 7.  Optimal selection of key parameters of multiscale entropy: (a) p = 0.2; (b) N = 2000; (c) τ = 10.

      当原始信号长度N值过小时, 计算结果无法准确反应信号特征; 当N过大时, 会造成计算量突增. 同样计算图5光纤信号MSE值, 比较N值的影响. 计算结果如图7(b)所示, τ = 10. 当N = 500时, PE与PB模式MSE值十分接近. 随着N的增长, 差异逐渐增大. 当N达到2000后, TEB与PB模式MSE值基本稳定, PE模式增速放缓. 此时三种信号MSE值已有足够的区分度, 继续增加N值, 区分程度提升较小, 且增大了计算量, 因此选用N = 2000.

      多尺度熵能够针对不同时间尺度对信号特征进行提取, 这是相对于样本熵的优势. 为选择最优的时间尺度因子τ, 计算了三种典型光纤信号在不同τ 下的MSE值, 如图7(c)所示. 整体上, 各尺度下PE模式MSE值最大, TEB模式最小. τ < 10时, PB和TEB模式的MSE值接近; τ > 10时, PE和PB模式相互靠近. 因此, 取τ = 10时, MSE值能较好地区分不同传热模式.

      使用上文确定的参数值(m = 2, p = 0.2, N = 2000, τ = 10), 计算8 MPa与10 MPa下各工况光纤信号MSE值, 结果如图8所示. MSE可以有效区分三种传热模式, 且两种压力下分界值相同, MSE > 0.9时为PE模式, 0.5 < MSE < 0.9时为PB模式, MSE < 0.5时为TEB模式. 仅10 MPa个别工况点PB与PE接近.

      图  8  基于多尺度熵的传热模式识别 (a) P = 7.955—8.048 MPa; (b) P = 9.928—10.035 MPa

      Figure 8.  Heat transfer modes identification based on multiscale entropy: (a) P = 7.955–8.048 MPa; (b) P = 9.928–10.035 MPa.

      PE模式信号相对稳定, 其他两种信号存在着大幅波动, 但计算结果却是PE模式MSE值最大, 似乎与直观印象不符. 这是由于MSE的计算中将信号的波动幅度归一化, 研究波动的规律性与周期性, 波动幅值对熵的影响较小. MSE实际体现系统的不确定性: PE模式信号波动主要来自电气测量噪声, 随机性较强, 规律性较差, 难以预测, 因此MSE值较大; TEB信号呈现明显的周期性波动, 规律性强, MSE值较小; 从图5图6可以看出, PB模式单个类气泡的产生和消散具有周期性, 但柱状流的摆动和消散又具有随机性, 因此信号周期性较弱, MSE值介于前两者之间.

    • 本文通过实验研究了sCO2在池内镍铬丝上的类沸腾传热特性, 实验在超临界压力高压釜中进行, 热流密度覆盖0—1800 kW/m2宽广范围. 采用电阻测温、高速相机、光纤探针等多种方法测量sCO2传热过程, 分析光纤探针信号特征, 并提出基于光纤探针和多尺度熵的SCF传热模式识别方法, 得到以下结论.

      1) 在宽广的工况范围类, 池内镍铬丝上sCO2hTw传热曲线先升后降, 随后在振荡中再次上升直至相对稳定, 其中存在4种传热模式. 当Tw低于Tpc时, 传热与亚临界自然对流相同. 当Tw超过Tpc后, 类气体包裹加热丝导致传热恶化, 类似亚临界膜态沸腾时气膜缓慢蒸发. 较高温度下镍铬丝受力失衡开始上下振动, 类气膜在Taylor不稳定性作用下失稳形成类气泡, 传热模式在类蒸发和类沸腾之间周期转换. 温度极高时镍铬丝温度分布不均, 持续产生类气泡, 完全进入类沸腾模式.

      2) 类气泡是类气膜失稳的产物, 形成机理与亚临界核态沸腾气泡完全不同, 但二者都通过搅动热边界层使传热增强, 是单相流传热中不会自发出现的现象. 这再次说明SCF在传热中不能被视作单相流, 应充分考虑类两相流型对传热模式的影响.

      3) 光纤探针热响应时间可视为0, 能准确测量镍铬丝上方流场温度波动. 本文分析了不同工况光纤信号频谱特征, TEB模式温度周期波动, 采用主频比可有效识别; 但PE和PB模式周期性都比较弱, 主频比容易混淆. 采用多尺度熵分析光纤信号可以有效识别PE, TEB和PB传热模式. PE模式信号中包含较多电气设备噪音, 随机性强难预测, MSE值高于0.9. TEB模式信号周期性强, 复杂度低, MSE值小于0.5. PB模式类气泡有弱周期性, MSE值介于二者之间.

      SCF密度连续变化, 没有相界面, 流型识别相对困难. 本文使用光纤探针可连续测量流体密度, 并利用多尺度熵成功识别类沸腾传热模式和流型, 尚未见其他文献有类似报道. 未来可将本文采用的技术和方法推广至SCF管内流动传热, 测量管内流体分布和流型特征, 为SCF类沸腾传热机理研究提供坚实依据.

参考文献 (40)

目录

    /

    返回文章
    返回