謝榕順,張廣旭,趙弟宏,盧功豪,洪 鋼,2,張堯立,2,*
(1.廈門(mén)大學(xué) 能源學(xué)院,福建 廈門(mén) 361105;2.福建省核能工程技術(shù)研究中心,福建 廈門(mén) 361105)
S-CO2布雷頓循環(huán)技術(shù)在高效太陽(yáng)能、地?zé)?、核能和清潔化石等能源系統(tǒng)中有很大的應(yīng)用潛能[1]。研究表明,使用S-CO2布雷頓循環(huán)發(fā)電技術(shù)可以將現(xiàn)有火電廠的熱電轉(zhuǎn)化效率從42%提高到近50%[2],顯著減少工廠規(guī)模和電力成本。因此,許多研究人員對(duì)S-CO2開(kāi)展了廣泛的研究。
對(duì)S-CO2傳熱的研究始于50年代,Pioro和Duffey[3]對(duì)通道內(nèi)流動(dòng)的二氧化碳的實(shí)驗(yàn)傳熱研究進(jìn)行了綜合評(píng)述,發(fā)現(xiàn)了3種不同的傳熱模式:一般傳熱模式(NHT)、傳熱惡化模式(HTD)和傳熱增強(qiáng)模式(THE)。而這3種傳熱模式的改變主要集中在擬臨界區(qū)域(壓力p=7.4 MPa,主流流體局部溫度Tb=30~40 ℃)附近,在高溫情況下并不顯著。為了調(diào)查這一現(xiàn)象,人們做了大量的研究,但仍然很難給出一個(gè)確定的解釋。近年來(lái),Rao等[4]、Cabeza等[5]、Ehsan等[6]、Xie等[7]對(duì)S-CO2在水平、豎直、環(huán)形管等不同幾何以及流動(dòng)狀態(tài)下的傳熱特性進(jìn)行研究,評(píng)估了各種邊界條件的影響,發(fā)現(xiàn)不同實(shí)驗(yàn)中的傳熱行為是復(fù)雜的,并且由于他們實(shí)驗(yàn)研究的溫度范圍大多集中在30~200 ℃之間,熱流密度的影響尤為顯著。McEligot等[8]研究了在較低熱流密度條件下湍流向?qū)恿鬓D(zhuǎn)變過(guò)程中的熱傳遞現(xiàn)象,并提出了再層流化理論,該理論認(rèn)為,強(qiáng)烈的流動(dòng)加速效應(yīng)會(huì)將對(duì)數(shù)律湍流剖面轉(zhuǎn)化為準(zhǔn)層流化,從而導(dǎo)致動(dòng)量傳輸減少和熱傳遞受損。Kurganov等[9]的實(shí)驗(yàn)結(jié)果表明,熱加速效應(yīng)對(duì)于中低溫下(Tb=30~80 ℃)小管徑中對(duì)傳熱的影響比較顯著。Jackson等[10]采用“等效無(wú)浮力流”概念,該概念將豎直流中的傳熱惡化歸因于由強(qiáng)烈的物理性質(zhì)變化引起的浮力改變,從而導(dǎo)致流型結(jié)構(gòu)和湍流強(qiáng)度的變化。浮升力對(duì)湍流產(chǎn)生的剪切應(yīng)力的影響也在后人的實(shí)驗(yàn)和數(shù)值研究中得到驗(yàn)證,但他們的試驗(yàn)工況大多都集中在擬臨界區(qū)域附近,在高溫區(qū)域的現(xiàn)象有待進(jìn)一步觀察。而Liu等[11]提出浮升力效應(yīng)很大程度上依賴于近壁區(qū)的流動(dòng)特性,由此發(fā)展出一種基于局部壁溫普朗特?cái)?shù)的新浮升力參數(shù),但尚未得到大量實(shí)驗(yàn)數(shù)據(jù)驗(yàn)證。
本文通過(guò)實(shí)驗(yàn)研究在高溫工況下S-CO2的傳熱特性,對(duì)比其與氣態(tài)CO2傳熱的差異,同時(shí)分析浮升力以及流動(dòng)加速效應(yīng)對(duì)于高溫工況下S-CO2傳熱的影響,最后將實(shí)驗(yàn)數(shù)據(jù)與不同類(lèi)型的現(xiàn)有經(jīng)驗(yàn)關(guān)聯(lián)式進(jìn)行比較,開(kāi)發(fā)一種新的適用于高溫工況下的S-CO2傳熱經(jīng)驗(yàn)關(guān)聯(lián)式,本文的工作能夠?yàn)镾-CO2的工程應(yīng)用提供數(shù)據(jù)支撐。
圖1為實(shí)驗(yàn)系統(tǒng)的示意圖。實(shí)驗(yàn)回路為一個(gè)閉式強(qiáng)迫循環(huán)回路,主要包括循環(huán)系統(tǒng)、冷卻系統(tǒng)、電加熱系統(tǒng)、測(cè)控系統(tǒng),主要設(shè)施有液下柱塞泵、科里奧利力質(zhì)量流量計(jì)、穩(wěn)壓器、調(diào)壓閥、雙管換熱器、工業(yè)冷水機(jī)、大功率程控直流電源和冷卻塔。為防止試驗(yàn)前不凝性氣體對(duì)試驗(yàn)結(jié)果的影響,需要經(jīng)抽真空吹掃除去不凝性氣體后,將純度99.95%的CO2注入蓄能器,經(jīng)冷凝器冷卻至高壓飽和狀態(tài)。然后液體被柱塞泵壓縮到工作壓力,在此壓力下它被注入回路。主回路中的CO2由液體柱塞泵驅(qū)動(dòng),進(jìn)入實(shí)驗(yàn)段進(jìn)行加熱。為了穩(wěn)定質(zhì)量流量,減少系統(tǒng)中的壓力脈動(dòng),在柱塞泵出口處安裝隔膜脈動(dòng)阻尼器。加熱后的CO2依次通過(guò)一次、二次換熱器冷卻后再次進(jìn)入液體柱塞泵,形成閉環(huán)。帶有冷卻塔部分的換熱器(一次換熱器)負(fù)責(zé)導(dǎo)出實(shí)驗(yàn)過(guò)程中的產(chǎn)生的大量熱量,帶有工業(yè)冷水機(jī)的換熱器(二次換熱器)負(fù)責(zé)調(diào)節(jié)實(shí)驗(yàn)段入口溫度。實(shí)驗(yàn)段本體通過(guò)可編程直流電源加熱。該電源具有自動(dòng)量程輸出特性,可在整個(gè)功率范圍內(nèi)提供寬泛的電壓和電流組合,提供大功率、穩(wěn)定的直流供電。實(shí)驗(yàn)系統(tǒng)的工況參數(shù)范圍列于表1。
表1 實(shí)驗(yàn)系統(tǒng)工況參數(shù)范圍
圖1 實(shí)驗(yàn)回路示意圖
實(shí)驗(yàn)段為內(nèi)徑10 mm、外徑12 mm的Incoloy-800H不銹鋼無(wú)縫管,通過(guò)耐高壓金屬管與系統(tǒng)主回路相連。如圖2所示,加熱段長(zhǎng)1 500 mm(L/d=150),有兩個(gè)絕熱段,分別為250 mm和300 mm,作為充分發(fā)展段。為準(zhǔn)確測(cè)量進(jìn)出口溫度和系統(tǒng)壓力,保證流體進(jìn)入加熱段的穩(wěn)定性,設(shè)置了兩個(gè)分別為500 mm以及600 mm(L/d≥50)的絕熱段作為充分發(fā)展段在受熱段的前端和后端。20個(gè)K型熱電偶(φ2 mm)排列在加熱段上,前10個(gè)每個(gè)間距50 mm,以便更加精準(zhǔn)測(cè)得擬臨界區(qū)域附近的溫度,后10個(gè)每個(gè)間距100 mm。測(cè)試段采用電加熱,兩端分別設(shè)置電絕緣法蘭。實(shí)驗(yàn)段包裹60 mm厚硅酸鋁陶瓷纖維保溫毯和100 mm厚巖棉保溫管進(jìn)行保溫。
圖2 實(shí)驗(yàn)段結(jié)構(gòu)示意圖
測(cè)量實(shí)驗(yàn)段進(jìn)、出口流體溫度的熱電偶與壁溫所用相同。使用智能壓力變送器測(cè)量系統(tǒng)壓力。實(shí)驗(yàn)段的質(zhì)量流量測(cè)量選用科氏質(zhì)量流量計(jì)。所有測(cè)量數(shù)據(jù)由計(jì)算機(jī)實(shí)時(shí)處理、顯示和保存。
為了通過(guò)實(shí)驗(yàn)獲得傳熱系數(shù)和努塞爾數(shù)(Nu),需要根據(jù)外壁溫度、進(jìn)出口溫度、壓力和熱流密度對(duì)數(shù)據(jù)進(jìn)行處理。
局部對(duì)流換熱系數(shù)h的計(jì)算公式為:
(1)
式中:Tw,i為內(nèi)管壁的局部溫度,℃;qw為內(nèi)壁面熱流密度,kW/m2;Tb根據(jù)局部流體焓Hb,x和壓力p從NIST REFPROP 9.0數(shù)據(jù)庫(kù)中獲得。直流電源提供的總熱量通過(guò)測(cè)量的電流乘以電壓計(jì)算得出。但由于受環(huán)境溫度的影響,熱量損失不可避免,因此在計(jì)算熱通量之前需要先計(jì)算熱效率η,公式如下:
(2)
式中:U和I分別為施加到測(cè)試部分的電壓(V)和電流(A);Hin和Hout分別為由相應(yīng)的壓力和溫度確定的入口和出口焓,kJ;m為質(zhì)量流量,kg/s。
實(shí)驗(yàn)段電阻產(chǎn)生的等效內(nèi)熱源qv為:
(3)
式中:Ri、Ro、Di、Do分別為管子的內(nèi)外半徑和直徑,m;Lh為加熱段長(zhǎng)度,m。
加熱段內(nèi)壁面提供給管道內(nèi)流體的真實(shí)熱流密度qw為:
(4)
由于加熱段為均勻熱流邊界條件。根據(jù)熱平衡可以計(jì)算距離實(shí)驗(yàn)段x處的主流平均焓Hb,x:
(5)
式中,dx為加熱段距入口起始點(diǎn)的距離。
鎧裝熱電偶測(cè)量的是實(shí)驗(yàn)段的外壁溫Tw,out,通過(guò)一維內(nèi)熱源導(dǎo)熱公式可以推導(dǎo)出外壁溫測(cè)點(diǎn)處對(duì)應(yīng)的內(nèi)壁溫Tw,in,具體公式如下:
(6)
式中,λ為Incoloy800H的導(dǎo)熱系數(shù),W/(m·℃)。λ會(huì)隨溫度的升高而增加。制造廠商提供的導(dǎo)熱系數(shù)計(jì)算公式如下:
λ=11.138+0.017 1T
(7)
最后,努塞爾數(shù)Nuexp的計(jì)算式為:
Nuexp=hDi/λb
(8)
式中,λb為主流溫度下的S-CO2的熱導(dǎo)率,W/(m·K)。
本實(shí)驗(yàn)依據(jù)《測(cè)量不確定度評(píng)定和表示》[14]對(duì)實(shí)驗(yàn)數(shù)據(jù)開(kāi)展B類(lèi)不確定度評(píng)定。另外,測(cè)量量分為直接測(cè)得量和間接測(cè)得量?jī)煞N,對(duì)于間接測(cè)得量,一般由多個(gè)直接測(cè)得量計(jì)算得到,假設(shè)一個(gè)間接測(cè)得量為R,它是直接測(cè)得量x1,x2,x3,…,xn的函數(shù),即R=f(x1,x2,x3,…,xn),可通過(guò)Moffat[15]給出的公式計(jì)算間接測(cè)得量的不確定度:
(9)
熱量密度q、傳熱系數(shù)h、努塞爾數(shù)Nu的不確定度計(jì)算如下所示:
(10)
(11)
(12)
根據(jù)上述不確定度傳遞公式,得到本次試驗(yàn)主要參數(shù)的不確定度,結(jié)果列于表2。
表2 實(shí)驗(yàn)測(cè)量和計(jì)算的不確定度
圖3示出了壓力為9 MPa、入口溫度為290 K時(shí)兩次實(shí)驗(yàn)外壁面溫度沿軸向管道的變化情況,從趨勢(shì)上來(lái)看,兩次給出的測(cè)試結(jié)果基本相同,相對(duì)偏差不超過(guò)1%,說(shuō)明整個(gè)系統(tǒng)的可重復(fù)性高,穩(wěn)定性好。
圖3 重復(fù)性測(cè)試
圖4為不同狀態(tài)下CO2的傳熱特性變化情況,其中黑色線代表氣態(tài)CO2,其壓力為4.3 MPa;紅色線代表S-CO2,其壓力為7.8 MPa??梢钥闯?兩者換熱特性的變化趨勢(shì)基本相同,在整個(gè)實(shí)驗(yàn)段前1/4的區(qū)域,S-CO2主流溫度變化相比于氣態(tài)更為平緩。兩者的換熱系數(shù)都沿著流動(dòng)方向降低,S-CO2換熱系數(shù)變化相比于氣態(tài)更為劇烈。在整個(gè)前1/4的實(shí)驗(yàn)段流體溫度約在300~400 K之間,S-CO2處于臨界點(diǎn)附近,因此傳熱系數(shù)的變化很劇烈,較小的溫度變化可能導(dǎo)致傳熱系數(shù)較大的差別。當(dāng)S-CO2流體平均溫度高于擬臨界溫度10 K左右時(shí),傳熱系數(shù)隨流體溫度的變化開(kāi)始變得平緩。而氣態(tài)CO2則在325 K左右時(shí)換熱系數(shù)變得較為平緩。從圖中可以發(fā)現(xiàn)S-CO2換熱系數(shù)遠(yuǎn)大于氣態(tài)CO2,特別在高溫區(qū)域(主流溫度大于473.15 K),這種現(xiàn)象隨著主流溫度的升高變得更加明顯,因此,雖然高溫工況下S-CO2換熱系數(shù)變化趨勢(shì)接近于氣態(tài)CO2,但在數(shù)值上有很大差異,在高溫流域,S-CO2的換熱明顯好于氣態(tài)CO2。
圖4 氣態(tài)CO2與S-CO2換熱特性對(duì)比
以往的研究[8-12]表明,在溫度較低時(shí),流動(dòng)加速效應(yīng)和浮升力效應(yīng)會(huì)影響超臨界流體的傳熱特性。為了確定高溫工況下熱加速和浮升力的影響,選擇目前廣泛使用的熱加速無(wú)量綱數(shù)Kv來(lái)評(píng)估熱加速效應(yīng)對(duì)傳熱的影響;選取浮升力無(wú)量綱數(shù)Bu來(lái)評(píng)價(jià)浮升力對(duì)傳熱的影響。
熱加速無(wú)量綱數(shù)Kv采用McEligot等[8]的定義,如下:
(13)
式中:βb為主流流體的局部體積膨脹系數(shù),1/K;μb為主流流體的局部動(dòng)力黏度,Pa·s;cp,b為主流溫度的比定壓熱容,J/(kg·K)。McEligot等[8]認(rèn)為,當(dāng)Kv>3×10-6時(shí),湍流減少,流動(dòng)可能過(guò)渡到層流,導(dǎo)致傳熱惡化。
圖5為在高熱流密度條件下,不同壓力下S-CO2沿軸向的換熱特性變化情況。從圖5可看出,在不同壓力下?lián)Q熱特性以及主流溫度變化差異不大,并且出口溫度高達(dá)700 K,與圖4類(lèi)似,整體S-CO2換熱系數(shù)沿軸線的變化趨勢(shì)為先減少后增加,且在高溫流體區(qū)域(即主流溫度大于473 K的區(qū)域,此時(shí)軸向距離為0.9 m,本文將在該點(diǎn)之后的區(qū)域稱為高溫流體區(qū)域)S-CO2換熱系數(shù)增大的幅度逐漸變強(qiáng)。略有差異的是在入口處有一段換熱強(qiáng)化現(xiàn)象發(fā)生,鑒于本文重點(diǎn)研究在溫度在470 K以上S-CO2換熱特性的變化情況,故對(duì)相關(guān)現(xiàn)象產(chǎn)生的原因不過(guò)多贅述。
圖5 高溫工況下S-CO2沿軸向換熱特性變化
圖6為在上述工況條件下,S-CO2沿軸向的熱加速無(wú)量綱數(shù)的變化情況,可以明顯發(fā)現(xiàn)在主流溫度較低的區(qū)域,其熱加速無(wú)量綱數(shù)Kv約為1.6×10-7,而在高溫流域該值僅為3.0×10-8,較上述研究人員所得的閾值小兩個(gè)數(shù)量級(jí),說(shuō)明在高溫流域內(nèi),熱加速效應(yīng)對(duì)于傳熱的影響可以忽略。
圖6 高溫工況下S-CO2沿軸向熱加速無(wú)量綱數(shù)變化
表征浮升力無(wú)量綱數(shù)選用Bu,用于預(yù)測(cè)S-CO2豎直方向流動(dòng)中由浮升力效應(yīng)引起的傳熱惡化或增強(qiáng)的起始點(diǎn),其定義為:
(14)
圖7為在上述工況條件下,S-CO2沿軸向的浮升力無(wú)量綱數(shù)的變化情況,可以明顯發(fā)現(xiàn)在主流溫度較低的區(qū)域,Bu在相應(yīng)的理論閾值上,同樣可以明顯發(fā)現(xiàn)在相關(guān)區(qū)域發(fā)生了傳熱的惡化以及增強(qiáng)現(xiàn)象,因此可以說(shuō)明在S-CO2較低溫度區(qū)域的傳熱與浮升力密切相關(guān)。然而在高溫流域(大于0.9 m的區(qū)域)相關(guān)的無(wú)量綱數(shù)值都低于相關(guān)參數(shù)的閾值。
圖7 高溫工況下S-CO2沿軸向浮升力無(wú)量綱數(shù)變化
上述結(jié)果說(shuō)明在高溫區(qū)域熱加速以及浮升力效應(yīng)對(duì)于S-CO2傳熱的影響幾乎可以忽略。從圖7可看出,該區(qū)域內(nèi)換熱系數(shù)沿軸向流動(dòng)方向單調(diào)上升,無(wú)明顯換熱惡化現(xiàn)象出現(xiàn)。因?yàn)榇藭r(shí)主流區(qū)為低密度的類(lèi)氣態(tài)流體,流體物性變化較為平緩,傳熱規(guī)律遵循單相強(qiáng)制對(duì)流換熱。并且浮升力影響主要針對(duì)徑向上熱物性差異,熱加速影響主要針對(duì)軸向熱物性差異上,由于高溫區(qū)域內(nèi)S-CO2的比熱和導(dǎo)熱系數(shù)等物性隨溫度的變化較為緩慢,因此徑向和軸向熱物性差別不大,浮升力及熱加速效應(yīng)對(duì)高溫區(qū)S-CO2傳熱的影響可以忽略。
為了得到適用于高溫工況下的S-CO2對(duì)流換熱經(jīng)驗(yàn)關(guān)聯(lián)式,首先參考現(xiàn)有S-CO2對(duì)流換熱經(jīng)驗(yàn)關(guān)聯(lián)式的形式,然后使用實(shí)驗(yàn)數(shù)據(jù)進(jìn)行擬合得到最終的經(jīng)驗(yàn)關(guān)聯(lián)式。
雖然S-CO2在高主流溫度區(qū)域的換熱變化趨勢(shì)類(lèi)似于氣態(tài)CO2,但在傳熱強(qiáng)度上明顯高于CO2;并且相關(guān)區(qū)域內(nèi)浮升力以及熱加速效應(yīng)對(duì)他們的影響可以忽略,因此新提出的經(jīng)驗(yàn)關(guān)聯(lián)式形式基于經(jīng)典傳熱關(guān)聯(lián)式,將直接影響流體傳熱的熱物理性質(zhì),如流體密度、動(dòng)力黏度、比熱等引入。最常用是在經(jīng)典傳熱關(guān)聯(lián)式中引入主流流動(dòng)區(qū)域和近壁流體區(qū)域流體熱物理性質(zhì)的比值。因此,選擇5個(gè)與之相關(guān)的經(jīng)驗(yàn)關(guān)聯(lián)式,與當(dāng)前測(cè)量的傳熱數(shù)據(jù)進(jìn)行比較。表3列出選定的經(jīng)驗(yàn)關(guān)聯(lián)式。
表3 預(yù)測(cè)豎直管中S-CO2傳熱的經(jīng)驗(yàn)關(guān)聯(lián)式
圖8 不同經(jīng)驗(yàn)關(guān)聯(lián)式下努塞爾數(shù)預(yù)測(cè)值與實(shí)驗(yàn)值的比較
從8e可看出,雖然K-P公式的預(yù)測(cè)能力與Fewster以及Jackson類(lèi)似,但從公式可看出其形式較為復(fù)雜,不便于在工程上直接進(jìn)行應(yīng)用,因此考慮精確性以及便捷型的角度,結(jié)合前人對(duì)傳熱關(guān)聯(lián)式的研究,確定如下形式的經(jīng)驗(yàn)關(guān)聯(lián)式:
(15)
由該關(guān)聯(lián)式形式可以發(fā)現(xiàn),當(dāng)S-CO2為常物性時(shí),物性修正項(xiàng)變?yōu)?,關(guān)聯(lián)式變?yōu)槌N镄怨?具有相容性的優(yōu)勢(shì)。為了方便實(shí)踐中使用,忽略了對(duì)換熱影響較小的修正因子,同時(shí)在參照了一些關(guān)聯(lián)式系數(shù)的基礎(chǔ)上,固定了一些物性項(xiàng)的系數(shù),而后利用豎直圓管換熱數(shù)據(jù)庫(kù),采用最小二乘法進(jìn)行擬合,可以得到高溫工況下S-CO2在豎直圓管的新傳熱關(guān)聯(lián)式:
(16)
該關(guān)聯(lián)式適用于S-CO2在豎直加熱管的換熱計(jì)算。具體使用范圍:系統(tǒng)壓力為7.40~10.22 MPa;流體溫度為450~800 K;熱流密度為200~1 000 kW/m2;質(zhì)量流密度為189.45~1 514.46 kg/(m2·s);雷諾數(shù)為2.59×104~3.28×105;普朗特?cái)?shù)為0.72~14.29。
本文提出經(jīng)驗(yàn)關(guān)聯(lián)式結(jié)果與實(shí)驗(yàn)值的對(duì)比情況如圖9所示,高溫工況下新傳熱經(jīng)驗(yàn)關(guān)聯(lián)式預(yù)測(cè)的努塞爾數(shù)除少數(shù)點(diǎn)外基本與實(shí)驗(yàn)數(shù)據(jù)吻合較好,96.95%的數(shù)據(jù)點(diǎn)在±10%以內(nèi),同時(shí)平均絕對(duì)百分比誤差為7.83%、計(jì)算誤差的標(biāo)準(zhǔn)差為3.15%、相關(guān)系數(shù)高達(dá)0.997。因此新建立的經(jīng)驗(yàn)關(guān)聯(lián)式比其他五種文中所列舉的關(guān)聯(lián)式具有更好的精度,能夠?yàn)楦邷毓r下計(jì)算S-CO2傳熱系數(shù)提供參考。
圖9 新傳熱經(jīng)驗(yàn)關(guān)聯(lián)式下努塞爾數(shù)預(yù)測(cè)值與實(shí)驗(yàn)值的比較
本文通過(guò)實(shí)驗(yàn)研究了在高溫工況下S-CO2的傳熱特性,對(duì)比了其與氣態(tài)CO2傳熱的差異,同時(shí)分析了浮升力以及流動(dòng)加速效應(yīng)對(duì)于高溫工況下S-CO2傳熱的影響,并且開(kāi)發(fā)了一種新的適用于高溫工況下的S-CO2傳熱經(jīng)驗(yàn)關(guān)聯(lián)式,主要結(jié)論如下:
1) 高溫工況(450~800 K)下S-CO2換熱系數(shù)變化趨勢(shì)接近于氣態(tài)CO2,但在數(shù)值上有很大差異,最大可以達(dá)到1 200 W(m2·K),在近似溫度情況下,S-CO2的換熱明顯好于氣態(tài)CO2;
2) 對(duì)于高溫工況下S-CO2的傳熱,熱加速效應(yīng)以及浮升力效應(yīng)其影響可以忽略不計(jì);
3) 通過(guò)考慮密度以及比熱物性的變化,提出了一種適用于高溫工況S-CO2豎直上升流動(dòng)換熱關(guān)聯(lián)式,通過(guò)與已有實(shí)驗(yàn)數(shù)據(jù)的比較,結(jié)果表明新關(guān)聯(lián)式能較好地預(yù)測(cè)換熱情況。