遠(yuǎn)臨界點垂直管內(nèi)超臨界CO2換熱特性的數(shù)值模擬研究

摘要：超臨界CO2在核能發(fā)電、太陽能發(fā)電、低溫制冷、航空航天等領(lǐng)域有著重要應(yīng)用。目前對超臨界CO2管內(nèi)對流換熱的研究大多在臨界點溫區(qū)附近，而在遠(yuǎn)離臨界點高溫高壓條件下的超臨界CO2換熱規(guī)律尚不明晰。在高溫高壓下進(jìn)行了數(shù)值模擬研究，探究了質(zhì)量流量、入口溫度、系統(tǒng)壓力、熱流密度和管徑對對流換熱系數(shù)的影響，并分析了由這些工況變化引起的浮升力和流動加速效應(yīng)對換熱特性的影響。結(jié)果表明：隨著質(zhì)量流量、入口溫度、系統(tǒng)壓力和熱流密度的增加，對流換熱系數(shù)增大；在不同熱流密度條件下，流體的對流換熱系數(shù)差值沿流動方向逐漸擴(kuò)大；對流換熱系數(shù)隨著管徑增大而減小。相較于臨界點附近的換熱規(guī)律，熱流密度和管徑對對流換熱系數(shù)的影響存在差異?？傮w而言，壓力對對流換熱系數(shù)的影響相對較小。該研究對理解和完善超臨界流體換熱規(guī)律、指導(dǎo)高效安全換熱器設(shè)計具有重要意義和工程價值。

關(guān)鍵詞：超臨界CO2；遠(yuǎn)臨界點；垂直管；換熱特性；數(shù)值模擬

中圖分類號：TK-9""" 文獻(xiàn)標(biāo)志碼：A""" 文章編號：1002-4026（2025）01-0083-13

開放科學(xué)（資源服務(wù)）標(biāo)志碼（OSID）：

DOI：10.3976/j.issn.1002-4026.20240036【能源與動力】

收稿日期：2024-03-11

基金項目：青島博士后應(yīng)用研究項目（QDBSH20220201001）

作者簡介：趙崇鑫（1997—），男，碩士，研究方向為超臨界換熱。E-mail：2567210408@qq.com

*通信作者，何燕（1973—），女，教授，博士，泰山學(xué)者，博士生導(dǎo)師，研究方向為納米材料。E-mail：heyanqustid@163.com

魏振文（1965—），男，高級工程師，研究方向為高溫?fù)Q熱。E-mail：z3692581471214@163.com

Numerical study on heat transfer characteristics of supercritical CO2

in vertical tubes at far-critical points

ZHAO Chongxin1，CUI Jianbo2，JIN Yanchao2，HAN Yazhou2，

WU Gongpeng1，2，HE Yan1*，WEI Zhenwen2*

（1. College of Electromechanical Engineering， Qingdao University of Science and Technology， Qingdao 266000， China;

2. Qingdao Doright Energy Saving Equipment Company Limited，Qingdao 266000， China）

Abstract∶Supercritical CO2 plays an important role in many applications such as nuclear power generation， solar power generation， cryogenic refrigeration， and aerospace. Currently， the majority of studies on supercritical CO2 convective heat transfer in tubes focus on the temperature range near the critical point， while the heat transfer patterns at high temperature and pressure far from the critical point remain unclear and need to be further studied. In this study， numerical simulations were performed to analyze the effects of mass flow， inlet temperature， system pressure， heat flux density， and tube diameter on the convective heat transfer coefficient at high temperature and pressure， as well as the effects of buoyancy and flow acceleration caused by operating conditions on the heat transfer characteristics. The results show that the convective heat transfer coefficient increases with increasing mass flow， inlet temperature， system pressure， and heat flux density. The difference in convective heat transfer coefficient gradually grows along the flow direction under different heat flux densities. Convective heat transfer coefficient decreases with increasing tube diameter. Compared with the heat transfer patterns near the critical point， heat flux density and tube diameter exert different effects on the convective heat transfer coefficient. In general， the effects of pressure on the convective heat transfer coefficient are small. This study provides significant values to understand the law of supercritical fluid heat transfer and guide the design of efficient and safe heat exchanger.

Key words∶supercritical CO2; far critical point; vertical tube; heat transfer characteristics; numerical simulation

雙碳目標(biāo)對國內(nèi)的能源結(jié)構(gòu)提出新的要求，相比于水、氟利昂等介質(zhì)，超臨界CO2具有化學(xué)性質(zhì)穩(wěn)定、安全無毒、容易獲取［1］等優(yōu)點，在超臨界狀態(tài)下具有優(yōu)越的換熱特性［2］，在核能發(fā)電、太陽能發(fā)電、低溫制冷、航空航天等領(lǐng)域有著重要應(yīng)用［3］。

目前對超臨界CO2管內(nèi)對流換熱的研究大多集中在臨界點（T = 304 K，p = 7.5 MPa）附近［4-6］。董文志等［7］探究了入口溫度（289～299 K）對傾斜圓管內(nèi)超臨界CO2換熱特性的影響。結(jié)果表明入口溫度的升高，頂母線和底母線壁面溫差減小，可以降低傳熱惡化現(xiàn)象。朱兵國等［8］在入口溫度范圍293～296 K、壓力范圍7.5～21 MPa、熱流密度范圍50～413 kW/m2、質(zhì)量流量范圍519～1 500 kg/（m2·s）工況下，通過實驗探究了熱流密度、系統(tǒng)壓力和浮升力對垂直上升管中超臨界CO2換熱特性的影響。結(jié)果表明增加熱流密度或減小壓力會導(dǎo)致傳熱惡化，而浮升力對傳熱惡化影響較大。莊曉如等［9］在壓力范圍7.5～9 MPa、溫度范圍800～1 050 K、質(zhì)量流量范圍200～500 kg/（m2·s）、熱流密度范圍100～800 kW/m2工況下進(jìn)行了超臨界CO2對流換熱數(shù)值模擬計算，研究了上述不同工況對對流換熱的影響。結(jié)果表明壓力對對流換熱系數(shù)的影響很小，而隨著質(zhì)量流量增大及熱流密度減小，對流換熱顯著增強(qiáng)；隨流體溫度增加，對流換熱系數(shù)增大但Nu數(shù)減小。

在超臨界二氧化碳布雷頓循環(huán)系統(tǒng)中，二氧化碳在熱源中受熱形成高溫高壓的遠(yuǎn)臨界點超臨界態(tài)，此時換熱器換熱性能及結(jié)構(gòu)安全性面臨更高要求和挑戰(zhàn)；且考慮到超臨界流體自身物性及換熱影響因素的復(fù)雜性，探究遠(yuǎn)臨界點條件下超臨界CO2的換熱特征，對理解和完善超臨界流體換熱規(guī)律、指導(dǎo)高效安全換熱器設(shè)計具有重要理論和工程價值。而分析上述結(jié)果發(fā)現(xiàn)，當(dāng)前針對超臨界流體換熱規(guī)律的研究主要集中在低溫或低壓及低溫低壓的近臨界點區(qū)域，在高溫高壓條件下處于遠(yuǎn)臨界點狀態(tài)的超臨界CO2換熱特性尚不清晰［10-14］。因此，為深入揭示高溫高壓下超臨界CO2換熱規(guī)律，本文通過數(shù)值模擬方法，系統(tǒng)研究管徑和入口溫度、系統(tǒng)壓力、質(zhì)量流量、熱流密度等因素對換熱特性的影響，探究高溫高壓工況下超臨界CO2的復(fù)雜換熱特性，從而助推新型超臨界二氧化碳循環(huán)換熱系統(tǒng)的應(yīng)用。

1" 物理模型

1.1" 物理模型及網(wǎng)格劃分

本文所研究的管道模型為垂直管，超臨界CO2在管道內(nèi)向上流動并在管道內(nèi)吸熱。如圖1所示，幾何模型總長為1 100 mm的管道。為了降低進(jìn)口段和出口段的影響，并確保熱流體在管道中保持充分發(fā)展?fàn)顟B(tài)，在進(jìn)口段設(shè)置長度為160 mm、出口段設(shè)置長度為100 mm的絕熱段，中間吸熱段長度為840 mm，采用均勻恒定熱流對實驗段進(jìn)行加熱，同時忽略管壁厚度。

根據(jù)朱兵國［15］的研究結(jié)果，在均勻周向加熱條件下，超臨界CO2在垂直管內(nèi)上升流動傳熱過程，可以簡化為二維軸對稱模型。為了提高模擬準(zhǔn)確性，本文建立了三維模型進(jìn)行研究。徑向及軸向網(wǎng)格劃分結(jié)果如圖2所示。由于靠近壁面附近的流體溫度和流速變化梯度較大，所以對模型壁面進(jìn)行網(wǎng)格加密，確?？拷诿嫣幍牡谝粚泳W(wǎng)格的無量綱壁面距離（y+）小于1。經(jīng)網(wǎng)格質(zhì)量檢測，網(wǎng)格質(zhì)量均大于0.7，全部角度都大于45°，說明網(wǎng)格質(zhì)量良好。

1.2" 湍流模型及控制方程

根據(jù)Zhang等［16］對幾種湍流模型的研究結(jié)果表明，不管是在正常傳熱工況下，還是在傳熱惡化工況下，SST k-ω模型依舊可以準(zhǔn)確預(yù)測壁面溫度等數(shù)據(jù)，而Standard k-ω、Standard k-ε和Realizable k-ε等模型會受到一些條件的制約而不能準(zhǔn)確預(yù)測所需要的數(shù)據(jù)。因此，選用SST k-ω模型進(jìn)行模擬。為了保證數(shù)值模擬計算結(jié)果的精確性，本文的超臨界CO2熱物性參數(shù)來源于美國國家標(biāo)準(zhǔn)與技術(shù)研究NIST REFPROP數(shù)據(jù)庫。收斂標(biāo)準(zhǔn)為連續(xù)性方程、動量守恒方程、能量守恒方程、湍動能方程及耗散率方程的殘差低于10-6，控制方程如下：

連續(xù)性方程

（ρu）xi=0，（1）

動量守恒方程

ρg+xjμ（ρui）xj-ρuiuj-pxi=（ρuiuj）xi，

（2）

能量守恒方程

xiμ（1Pr+μt/μPrt）ixi=（ρuii）xi，（3）

湍動能方程

Gk-Yk+xj（μ+μtσk）kxj=（ρuik）xi，（4）

耗散率方程

Gω-Yω+Dω+xj（μ+μtσω）kxj=（ρuiω）xi，（5）

式中，u為速度；ρ為密度；p為系統(tǒng)壓力；g為重力加速度；μ為黏性系數(shù)；μt為湍流黏性系數(shù)；Pr為普朗特數(shù)；Prt為湍流普朗特數(shù)；

k為湍動能，ω為湍動能耗散率;Gk為速度梯度產(chǎn)生的湍動能；Yk、Yω分別為由于湍流產(chǎn)生的k、ω的損耗；Gω為ω的產(chǎn)生項；Dω為交叉擴(kuò)散項；σk、σω分別為湍動能與耗散率對應(yīng)的普朗特數(shù)。

1.3" 邊界條件及計算式

進(jìn)口用質(zhì)量流量入口條件，出口用壓力出口條件，進(jìn)口段和出口段壁面為絕熱邊界條件，中間吸熱段壁面為恒定熱流密度邊界條件，具體工況如表1所示。

利用ANSYS FLUENT軟件對超臨界CO2在垂直吸熱管內(nèi)的換熱特性進(jìn)行數(shù)值模擬，計算得到超臨界CO2在垂直吸熱管內(nèi)的對流換熱系數(shù)、浮升力特征數(shù)和流動加速特征數(shù)等參數(shù)，進(jìn)而分析超臨界CO2的換熱特性。

對流換熱系數(shù)h（W·m-2·K-1）的計算式如下：

h=qTw-Tb， （6）

式中，Tw為壁面溫度；Tb為流體主流溫度。

1.4" 數(shù)據(jù)驗證

為了保證本文模擬數(shù)據(jù)的準(zhǔn)確性及計算效率，首先進(jìn)行網(wǎng)格無關(guān)性驗證，將模擬數(shù)據(jù)與文獻(xiàn)［17］的數(shù)據(jù)進(jìn)行了對比，選擇382萬個網(wǎng)格進(jìn)行數(shù)值模擬。從圖3中可以看出，數(shù)據(jù)取點位置到入

口的距離x和管道長

度L的比值與文獻(xiàn)數(shù)據(jù)變化趨勢保持一致，即在管入口處先迅速增加到一個峰值，然后迅速降低并趨于穩(wěn)定；同時也存在偏差，最高點出現(xiàn)的位置稍微向右偏移，這可能是由于模型中設(shè)置的絕熱段長度和截面初始流體速度分布與文獻(xiàn)中實際情況存在差異造成的。而由此計算得到的沿程壁面溫度的平均相對誤差小于10%，從而驗證了模擬結(jié)果的可靠性。

2" 結(jié)果與分析

2.1" 管徑的影響

圖4所示為不同管徑條件下，在壓力15 MPa、熱流密度60 kW/m2、質(zhì)量流量200 kg/（m2·s）、入口溫度773 K的垂直上升管中，壁溫和超臨界CO2對流換熱系數(shù)沿程的分布情況。

從圖4（a）中可以看出，在入口附近，不同管徑下壁溫的差別較小，而在xgt;0.28 m范圍內(nèi)，壁溫隨管徑的增大而減小，并且壁溫差值沿流動方向不斷增大，D = 4 mm的壁溫增幅最大。由圖5所示的截面流體溫度和密度分布可知，隨著管徑增加，流體平均溫度降低而密度增加，導(dǎo)致從管壁吸收的熱量增加，從而使得壁面溫度上升緩慢。與Bae等［18］的研究相反，由圖4（b）得出在遠(yuǎn)離臨界點條件下，對流換熱系數(shù)隨管徑增加而減小，當(dāng)管徑從4 mm增加到8 mm時，換熱系數(shù)的降幅明顯大于管徑從8 mm增加到12 mm時換熱系數(shù)的變化幅度，這說明超臨界CO2的換熱性能隨管徑增加而減弱，且換熱系數(shù)與管徑存在非線性關(guān)系。一方面，恒定流量條件下，管徑增加時，流速降低而密度增加，而從圖10可以看出，流體動力黏度在溫度變化范圍內(nèi)雖然降低，但量級較小，最終導(dǎo)致表征流體湍流程度的雷諾數(shù)隨管徑增加降低，進(jìn)而對流體換熱系數(shù)造成影響。另一方面，管徑增大導(dǎo)致流體的浮升力作用加強(qiáng)，傳熱惡化的現(xiàn)象開始凸顯［13］，而且從圖10知，流體的熱導(dǎo)率隨溫度降低而減小，流體傳熱性能被削弱。在以上綜合作用下，增大管徑削弱了流體的換熱性能。

2.2" 質(zhì)量流量的影響

圖6所示為壓力15 MPa、熱流密度60 kW/m2、管徑4 mm、入口溫度773 K的垂直上升管中，不同質(zhì)量流量條件下超臨界CO2流動換熱過程壁溫和對流換熱系數(shù)沿程的分布情況。

由圖6可知，隨著質(zhì)量流量的增加，壁溫逐漸降低，而且質(zhì)量流量在G = 100 kg/（m2·s）工況時，壁溫沿軸向增加幅度明顯高于其他兩個工況，根據(jù)文獻(xiàn)［19-21］可知，這是因為流體流速的增加會導(dǎo)致湍流擴(kuò)散速率的增大，流體與壁面間的熱量傳遞加快，導(dǎo)致壁面溫度降低，而在低質(zhì)量流量工況下的流體流速較小，湍流擴(kuò)散速率小，熱量在壁面易積聚，所以壁溫沿軸向變化幅度大。與壁面溫度變化相反，對流換熱系數(shù)則隨質(zhì)量流量增加而增大；與壁溫軸向變化相對應(yīng)，低質(zhì)量流量時，對流換熱系數(shù)的軸向變化更加明顯。質(zhì)量流量的變化導(dǎo)致流速的變化，進(jìn)而影響換熱過程。當(dāng)質(zhì)量流量增加時，流體流速增加，使得雷諾數(shù)增加的同時壁面處流體邊界層厚度減小，可以抑制類液膜厚度的增加［20］，導(dǎo)致導(dǎo)熱熱阻減小，增強(qiáng)了流體的換熱。

2.3" 入口溫度的影響

圖7所示為質(zhì)量流量200 kg/（m2·s）、壓力15 MPa、熱流密度60 kW/m2、管徑4 mm的垂直上升管中，不同入口溫度條件下超臨界CO2流動換熱過程壁溫和對流換熱系數(shù)沿程的分布情況。

由圖7可知，壁面溫度和對流換熱系數(shù)隨入口溫度的升高而增加，而且壁面溫度和對流換熱系數(shù)隨入口流體溫度均呈現(xiàn)相對均勻的變化。在恒定熱流密度工況下，對于相同幅度的流體溫升，當(dāng)入口溫度升高時，管內(nèi)主流流體從壁面吸收的熱量減少，從而使得壁溫升高。隨著溫度的增加，超臨界CO2密度逐漸減小，恒定質(zhì)量流量下流體流速增加，使得流體湍流強(qiáng)度增強(qiáng)，導(dǎo)致對流換熱系數(shù)增加。此外，由文獻(xiàn)［10］可知，入口溫度對換熱效果有密切聯(lián)系，當(dāng)入口溫度較低時，在壁面附近，低密度的蒸汽狀流體會持續(xù)擴(kuò)散，導(dǎo)致壁面上的低密度層流體的厚度增加，從而在內(nèi)壁上形成了一層較厚的氣膜，進(jìn)而增加了導(dǎo)熱熱阻。當(dāng)入口溫度升高時，低密度的蒸汽狀流體會被抑制膨脹，使得更多的高密度液體狀流體靠近壁面，傳熱過程得到強(qiáng)化，因此對流換熱系數(shù)隨入口溫度的升高而增大。

2.4" 壓力的影響

圖8所示為熱流密度60 kW/m2、質(zhì)量流量200 kg/（m2·s）、管徑4 mm、入口溫度773 K的垂直上升管中，不同壓力條件下超臨界CO2流動換熱過程壁溫和對流換熱系數(shù)沿程的分布情況。

由圖8可知，壓力的升高對壁溫的影響較小，而在同一主流溫度下，壁溫隨壓力的升高而降低，這與文獻(xiàn)［16］觀察到的現(xiàn)象一致，在相同溫度下提高壓力，使得流體的比體積減小，比熱容增加，流體升溫需要的熱量更多，導(dǎo)致壁溫減小。對流換熱系數(shù)隨著壓力的升高而增加，且在流動方向上保持了一致的變化?？紤]到超臨界CO2的密度隨著壓力的增大而增大，在恒定質(zhì)量流量下，使得流體流速減小，湍流強(qiáng)度減弱，對流換熱能力減小。但另一方面，隨著壓力的增大，CO2的熱導(dǎo)率增大，從而強(qiáng)化傳熱，兩者綜合作用使得對流換熱系數(shù)變大，說明超臨界CO2的熱導(dǎo)率在遠(yuǎn)臨界區(qū)域內(nèi)對換熱特性影響較大。

2.5" 熱流密度的影響

圖9所示壓力15 MPa、質(zhì)量流量200 kg/（m2·s）、管徑4 mm、入口溫度773 K的垂直上升管中，不同熱流密度條件下壁溫和對流換熱系數(shù)的沿程分布情況。

由圖9可知，壁溫隨壁面熱流密度的增加而增大，這是因為熱流密度增大，提供給管壁的熱量增多，從而使得壁溫升高。對流換熱系數(shù)表現(xiàn)出與壁溫相同的變化趨勢，超臨界CO2的熱導(dǎo)率隨溫度的升高而增大，靠近壁面的流體溫度因壁溫的升高而升高，隨著溫度增加，超臨界CO2密度逐漸減小，流體流速增加，流體湍流強(qiáng)度相應(yīng)增強(qiáng)，對流換熱系數(shù)增加。在不同熱流密度條件下，壁溫和對流換熱系數(shù)差值沿流體流動方向逐漸擴(kuò)大，并在出口處達(dá)到最大，這說明熱流密度對遠(yuǎn)臨界點區(qū)域超臨界CO2管內(nèi)換熱的影響隨流動發(fā)展愈發(fā)明顯。流體在流動過程中逐漸吸熱升溫，熱流密度增加會導(dǎo)致流體升溫速率和幅度相應(yīng)增加，而隨著流動發(fā)展，這種差異引起的熱量傳遞過程的差別也隨著擴(kuò)大，從而使得壁溫和對流換熱系數(shù)的差值增加。

2.6" 流動換熱效應(yīng)分析

根據(jù)文獻(xiàn)［22-29］的研究得出，浮升力和流動加速效應(yīng)能夠影響管道內(nèi)超臨界CO2的組成結(jié)構(gòu)，進(jìn)而影響超臨界CO2的換熱特性。不過，目前對這兩種效應(yīng)的判別依據(jù)，各種文獻(xiàn)提出了不同的看法，本文引入Jackson等［30］與McEligot等［31］提出的判別依據(jù)，浮升力（Bu）和流動加速參數(shù)（Ac）的計算公式具體如下：

Bu=GrRe2.7，（7）

Ac=4qβDμbcpRe2，（8）

其中

Gr=（ρb-ρw）ρbgD3μ2b，（9）

Re=GDμb，（10）

式中，Gr是格拉曉夫數(shù)，Re是雷諾數(shù)；β是體積膨脹系數(shù)；D是管徑；μb是主流流體黏度；Cp是定壓比熱；ρb是主流流體密度，ρw是近壁面流體密度；G是流體流速。

為了更好理解兩種效應(yīng)對傳熱的影響機(jī)理，圖10給出了不同超臨界壓力下CO2的熱物性

對比。從圖10（a）中可以看出，壓力越大，超臨界CO2密度也越大，但都隨著溫度的升高而減少；從圖10（b）圖可以看出，壓力越大，超臨界CO2的定壓比熱也越大，并且隨溫度的升高，p=8 MPa和p=15 MPa下的定壓比熱增加，p=22 MPa下的定壓比熱在700 K附近出現(xiàn)最低值；從圖10（c）圖和圖10（d）圖可以看出，壓力越大，超臨界CO2的熱導(dǎo)率和黏度也越大，而且都隨溫度的升高而增加。由此看出，隨壓力的增大，超臨界CO2的各熱物性均增大，而且壓力對密度和定壓比熱的影響更大。

Jackson等［30］與McEligot等［31］指出，當(dāng)浮升力Bult;1×10-5和流動加速

參數(shù)Aclt;3×10-6時，可以忽略浮升力和流動加速效應(yīng)。本文針對不同參數(shù)下的Bu和Ac進(jìn)行了計算，如圖11所示。從圖11可以看出，沿流體流動方向，Bu和Ac在各工況下均減小，這與流體溫度沿著x/D增大的方向升高有關(guān)。如圖11（a）所示，Bu和Ac隨著質(zhì)量流量的增加而減小。這歸因于增大質(zhì)量流量導(dǎo)致管內(nèi)流體流速增加，使得Re數(shù)增大，從而使Bu和Ac值降低。在G=200和300 kg/（m2·s）工況時，Bu沿流動方向減小且小于其閾值，說明此時浮升力效應(yīng)對超臨界CO2換熱特性影響較小。而在G=100 kg/（m2·s）工況下的Bu與其他質(zhì)量流量工況下的Bu的差值最大，且沿流動方向差值不斷減小，在x/Dlt;125時大于其閾值，而在x/Dgt;125時小于其閾值，此處對流換熱系數(shù)很小，換熱能力不強(qiáng)。在3種不同質(zhì)量流量工況下的Ac小于其閾值，說明此時流動加速效應(yīng)對超臨界CO2換熱特性影響不大，可以忽略。

如圖11（b）所示，Bu和Ac隨著溫度的升高而減小，由圖10（a）可知，超臨界CO2的密度隨溫度的升高而減少，在定質(zhì)量流量（G=ρu）下，流體流速增大，使Re數(shù)增大，從而導(dǎo)致Bu和Ac均減小，這與質(zhì)量流量對Bu和Ac的影響效果是一致的。在這3種溫度工況下的Bu和Ac都小于其閾值，說明此時浮升力和流動加速效應(yīng)對超臨界CO2換熱特性影響不大，可以忽略。

圖11（續(xù)）

如圖11（c）所示，Bu隨著壓力的增大而增大，由圖10（d）可知，超臨界CO2的黏度隨壓力的升高而變大，從而使得（7）式中的Bu增大。而Ac隨壓力的增大而減小，由圖10（b）和（d）可知，超臨界CO2的定壓比熱和黏度隨壓力的升高而增大，從而使得（8）式中的Ac減小，而且p=8 MPa下的定壓比熱的斜率最大，p=15 MPa和p=22 MPa下的定壓比熱的斜率緩慢增加，從而使得Ac的差值產(chǎn)生變化。但在這3種壓力工況下的Bu和Ac遠(yuǎn)小于其閾值，說明此時浮升力和流動加速效應(yīng)對超臨界CO2換熱特性影響不大，可以忽略。

如圖11（d）所示，Bu和Ac都隨著熱流密度的增加而增大，因為熱流密度增加，使超臨界CO2的溫度升高，導(dǎo)致密度減小，從而浮升力和流動加速效應(yīng)影響變大，進(jìn)而減小了換熱過程。但在這3種熱流密度工況下的Bu和Ac都小于其閾值，說明此時浮升力和流動加速效應(yīng)對超臨界CO2換熱特性影響不大，可以忽略。

3" 結(jié)論

本文利用數(shù)值模擬的方法對超臨界CO2在垂直吸熱管內(nèi)遠(yuǎn)離臨界點的換熱特性進(jìn)行了研究，分析了管徑、質(zhì)量流量、溫度、壓力和熱流密度以及浮升力和流動加速效應(yīng)對超臨界CO2換熱特性的影響，得出如下結(jié)論：

（1）在入口段，管徑對壁面溫度的影響較小，而在xgt;0.28范圍內(nèi)，隨著管徑的增加，壁溫開始明顯減小。在遠(yuǎn)臨界點區(qū)域，相較于質(zhì)量流量和熱流密度，入口溫度對壁溫的影響在整個軸向隨管徑增大表現(xiàn)為均勻增加，而壓力對壁溫的影響在4～12 mm的管徑范圍內(nèi)則顯著減弱。

（2）在遠(yuǎn)臨界點區(qū)域，對流換熱系數(shù)隨著管徑增加而表現(xiàn)出非線性減小趨勢，管徑由4 mm增加到8 mm時的變化幅度明顯大于從8 mm增加到12 mm時的變化幅度；隨著入口溫度和熱流密度的增加，超臨界CO2的密度減小，湍流流動狀態(tài)增強(qiáng)，從而提高了對流換熱系數(shù)；增大質(zhì)量流量使得流體流速增加，削弱了傳熱惡化；與壁溫的影響一致，壓力對對流換熱系數(shù)的影響最小。

（3）浮升力和流動加速效應(yīng)在入口段最為明顯，并隨著流動發(fā)展逐漸減弱。在遠(yuǎn)臨界點區(qū)域內(nèi)，相同條件下浮升力作用大于流動加速效應(yīng)，特別地，在G=100 kg/（m2·s）的低質(zhì)量流量條件下，入口段浮升力作用不可忽略。而由于遠(yuǎn)臨界點區(qū)域熱物性參數(shù)對工況條件的敏感性降低，浮升力和流動加速效應(yīng)對超臨界CO2換熱特性的影響普遍較小。

參考文獻(xiàn)：

［1］LIANG Y C， SUN Z L， DONG M R， et al. Investigation of a refrigeration system based on combined supercritical CO2 power and transcritical CO2 refrigeration cycles by waste heat recovery of engine［J］. International Journal of Refrigeration， 2020， 118： 470-482. DOI： 10.1016/j.ijrefrig.2020.04.031.

［2］SINGH A S， CHOUDHARY T， SANJAY S. Thermal analysis of aircraft auxiliary power unit： potential of super-Critical CO2 brayton cycle［C］. Aero Tech Americas.2019. DOI：10.4271/2019-01-1391.

［3］錢中. 微型換熱器瞬態(tài)傳熱分析［J］. 壓力容器， 2011， 28（9）： 26-29. DOI： 10.3969/j.issn.1001-4837.2011.09.006.

［4］CABEZA L F， DE GRACIA A， FERNNDEZ A I， et al. Supercritical CO2 as heat transfer fluid： A review［J］. Applied Thermal Engineering， 2017， 125： 799-810. DOI： 10.1016/j.applthermaleng.2017.07.049.

［5］XIE J Z， LIU D C， YAN H B， et al. A review of heat transfer deterioration of supercritical carbon dioxide flowing in vertical tubes： Heat transfer behaviors， identification methods， critical heat fluxes， and heat transfer correlations［J］. International Journal of Heat and Mass Transfer， 2020， 149： 119233. DOI： 10.1016/j.ijheatmasstransfer.2019.119233.

［6］BOVARD S， ABDI M， NIKOU M R K， et al. Numerical investigation of heat transfer in supercritical CO2 and water turbulent flow in circular tubes［J］. The Journal of Supercritical Fluids， 2017， 119： 88-103. DOI： 10.1016/j.supflu.2016.09.010.

［7］董文志， 韋武， 周亭羽， 等. 入口溫度對傾斜圓管內(nèi)超臨界CO2的流動傳熱影響研究［J］. 計算機(jī)與數(shù)字工程， 2023， 51（9）： 2165-2170. DOI： 10.3969/j.issn.1672-9722.2023.09.041.

［8］朱兵國， 吳新明， 張良， 等. 垂直上升管內(nèi)超臨界CO2流動傳熱特性研究［J］. 化工學(xué)報， 2019， 70（4）： 1291-1299. DOI： 10.11949/j.issn.0438-1157.20180695.

［9］莊曉如， 徐心海， 楊智， 等. 高溫吸熱管內(nèi)超臨界CO2傳熱特性的數(shù)值模擬［J］. 物理學(xué)報， 2021， 70（3）： 170-182. DOI： 10.7498/aps.70.20201005.

［10］朱兵國， 張海松， 孫恩慧， 等. 超高參數(shù)CO2在垂直管中的傳熱分析［J］. 化工進(jìn)展， 2019， 38（11）： 4880-4889. DOI： 10.16085/j.issn.1000-6613.2019-0582.

［11］QIU Y， LI M J， HE Y L， et al. Thermal performance analysis of a parabolic trough solar collector using supercritical CO2 as heat transfer fluid under non-uniform solar flux［J］. Applied Thermal Engineering， 2017， 115： 1255-1265. DOI： 10.1016/j.applthermaleng.2016.09.044.

［12］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.

［13］靳遵龍， 劉東來， 劉敏珊， 等. 超臨界CO2冷卻條件下水平微圓管中對流換熱特性［J］. 壓力容器， 2012， 29（7）： 9-13. DOI： 10.3969/j.issn.1001-4837.2012.07.002.

［14］洪瑞， 袁寶強(qiáng)， 杜文靜. 垂直上升管內(nèi)超臨界二氧化碳傳熱惡化機(jī)理分析［J］. 化工學(xué)報， 2023， 74（8）： 3309-3319. DOI： 10.11949/0438-1157.20230472.

［15］朱兵國. 超臨界二氧化碳垂直管內(nèi)對流換熱研究［D］. 北京： 華北電力大學(xué)， 2020.

［16］ZHANG Q， LI H X， LIU J L， et al. Numerical investigation of different heat transfer behaviors of supercritical CO2 in a large vertical tube［J］. International Journal of Heat and Mass Transfer， 2020， 147： 118944. DOI： 10.1016/j.ijheatmasstransfer.2019.118944.

［17］LI D， XU X X， CAO Y， et al. The characteristics and mechanisms of self-excited oscillation pulsating flow on heat transfer deterioration of supercritical CO2 heated in vertical upward tube［J］. Applied Thermal Engineering， 2022， 202： 117839. DOI： 10.1016/j.applthermaleng.2021.117839.

［18］BAE Y Y， KIM H Y， KANG D J. Forced and mixed convection heat transfer to supercritical CO2 vertically flowing in a uniformly-heated circular tube［J］. Experimental Thermal and Fluid Science， 2010， 34（8）： 1295-1308. DOI： 10.1016/j.expthermflusci.2010.06.001.

［19］閆晨帥. 超臨界二氧化碳流動傳熱數(shù)值模擬研究［D］. 北京： 華北電力大學(xué)， 2021.

［20］尹少軍. 圓管內(nèi)超臨界二氧化碳傳熱特性數(shù)值模擬［D］. 北京： 華北電力大學(xué)， 2021.

［21］朱鑫杰. 超臨界CO2垂直上升和下降對流傳熱特性實驗研究［D］.北京： 華北電力大學(xué)， 2021.

［22］ZHANG Q， LI H X， KONG X F， et al. Special heat transfer characteristics of supercritical CO2 flowing in a vertically-upward tube with low mass flux［J］. International Journal of Heat and Mass Transfer， 2018， 122： 469-482. DOI： 10.1016/j.ijheatmasstransfer.2018.01.112.

［23］LEI Y C， CHEN Z Q. Numerical study on cooling heat transfer and pressure drop of supercritical CO2 in wavy microchannels［J］. International Journal of Refrigeration， 2018， 90： 46-57. DOI： 10.1016/j.ijrefrig.2018.03.023.

［24］XIANG M R， GUO J F， HUAI X L， et al. Thermal analysis of supercritical pressure CO2 in horizontal tubes under cooling condition［J］. The Journal of Supercritical Fluids， 2017， 130： 389-398. DOI： 10.1016/j.supflu.2017.04.009.

［25］XU R N， LUO F， JIANG P X. Buoyancy effects on turbulent heat transfer of supercritical CO2 in a vertical mini-tube based on continuous wall temperature measurements［J］. International Journal of Heat and Mass Transfer， 2017， 110： 576-586. DOI： 10.1016/j.ijheatmasstransfer.2017.03.063.

［26］張海松， 朱鑫杰， 朱兵國， 等. 浮升力和流動加速對超臨界CO2管內(nèi)流動傳熱影響［J］. 物理學(xué)報， 2020， 69（6）： 136-145. DOI： 10.7498/aps.69.20191521.

［27］劉光旭， 黃彥平， 王俊峰， 等. 浮升力和流動加速效應(yīng)對超臨界CO2傳熱影響研究［J］. 核動力工程， 2016， 37（2）： 48-51. DOI： 10.13832/j.jnpe.2016.02.0048.

［28］LIU S， HUANG Y P， LIU G X， et al. Improvement of buoyancy and acceleration parameters for forced and mixed convective heat transfer to supercritical fluids flowing in vertical tubes［J］.International Journal of Heat and Mass Transfer， 2017， 106：1144-1156.

［29］黃宇， 段倫博. 超臨界流體流動加速效應(yīng)及其判別式研究進(jìn)展［J］. 動力工程學(xué)報， 2022， 42（1）： 94-100. DOI： 10.19805/j.cnki.jcspe.2022.01.012.

［30］JACKSON J D， HALL W B. Influences of buoyancy on heat transfer to fluids flowing in vertical tubes under turbulent conditions［J］.Institution of Mechanical Engineers， Conference Publications， 1979， 2： 613-640.

［31］MCELIGOT D M， COON C W， PERKINS H C. Relaminarization in tubes［J］. International Journal of Heat and Mass Transfer， 1970， 13（2）： 431-433. DOI： 10.1016/0017-9310（70）90118-3.