板式擴(kuò)散焊矩形微通道換熱器中S-CO2 流動與傳熱特性的數(shù)值模擬

李 洪， 任 燕， 章立新， 高 明， 劉婧楠， 楊其國

(1.上海理工大學(xué) 能源與動力工程學(xué)院，上海 200093； 2.上海市動力工程多相流動與傳熱重點實驗室，上海 200093)

發(fā)展新型高效發(fā)電技術(shù)是提高能源利用效率、實現(xiàn)節(jié)能環(huán)保最直接有效的方式[1]。相對于傳統(tǒng)工質(zhì)的動力循環(huán)系統(tǒng)，超臨界二氧化碳(S-CO2)布雷頓循環(huán)系統(tǒng)具有效率高、潔凈度高、靈敏度高、成本低以及結(jié)構(gòu)緊湊等優(yōu)點[2]，被認(rèn)為是新興能源利用方式中最具應(yīng)用前景的能源轉(zhuǎn)換系統(tǒng)之一，得到了國際學(xué)術(shù)界和工業(yè)界的廣泛關(guān)注[3]。

目前，S-CO2布雷頓循環(huán)已成為第四代先進(jìn)核能系統(tǒng)的備選熱力方案之一[4]。循環(huán)系統(tǒng)中的熱交換器(如預(yù)冷器、回?zé)崞鞯?是系統(tǒng)中的重要部件，直接影響著整個系統(tǒng)的高效性和緊湊性。印刷電路板式換熱器(PCHE)自面世以來，由于性能優(yōu)勢，且契合我國能源技術(shù)轉(zhuǎn)型升級和未來碳中和的需求[5]，在太陽能工程、液化天然氣工程和核反應(yīng)堆工程等領(lǐng)域得到了廣泛的關(guān)注[6]。研究表明，使用條件、流體流態(tài)及通道結(jié)構(gòu)等因素均會影響PCHE的流動與傳熱特性，且在確定的應(yīng)用背景下，通道結(jié)構(gòu)的影響最為顯著。PCHE通道結(jié)構(gòu)分為連續(xù)型和非連續(xù)型，涵蓋直線型、折線型、S形和翼型等多種形式；通道截面包括圓形、半圓形、矩形、三角形和梯形等多種形狀。受制于蝕刻加工技術(shù)，目前半圓形截面通道的研究應(yīng)用較為廣泛[7]。Geyer等[8]在恒定熱流密度的邊界條件下，對比分析了半圓形截面直線型通道與S形通道的傳熱效果，結(jié)果表明對于充分發(fā)展的層流流動而言，S形通道能夠達(dá)到強(qiáng)化傳熱的目的，同等工況下的強(qiáng)化效果可達(dá)2.5～3倍。Gupta等[9]探索了恒定熱流密度時通道截面形狀對傳熱特性的影響，結(jié)果表明在雷諾數(shù)(Re)低于200時，不同截面形狀的傳熱效果為三角形>矩形>半圓形>圓形。Lee等[10]對比了折線型通道不同截面形狀對傳熱特性和壓降特性的影響，證實矩形通道同時具有最佳的傳熱特性和最大的阻力損失。范凌灝等[11]的研究發(fā)現(xiàn)，在保持矩形微通道當(dāng)量直徑不變的情況下，隨著截面高寬比的增大，通道表面的平均努塞爾數(shù)逐漸升高，傳熱性能得到提升，但壓降也隨之增大，流動惡化，微通道換熱器的綜合性能降低。Cui等[12]對翼型翅片通道的分析結(jié)果表明，將形狀適宜的翅片交錯排布可以使得近壁邊界層周期性發(fā)展和消失，削弱邊界層對流動傳熱的影響，實現(xiàn)強(qiáng)化傳熱。綜上可知，優(yōu)化通道的截面形狀及擾流部件有助于實現(xiàn)強(qiáng)化傳熱、改善流動特性，但往往會造成生產(chǎn)工藝復(fù)雜、生產(chǎn)成本高昂等問題。因此，加工工藝便捷高效的矩形通道型PCHE在實現(xiàn)其流動與傳熱綜合性能優(yōu)化之后，將會擁有廣闊的應(yīng)用前景。

此外，S-CO2布雷頓循環(huán)系統(tǒng)中PCHE內(nèi)的工質(zhì)為超臨界流體，超臨界流體是當(dāng)CO2的溫度和壓力均高于臨界點(31 ℃和7.4 MPa)之后的一種狀態(tài)，其基本物理性質(zhì)(包括密度、黏度及擴(kuò)散性等)均處于氣體與液體之間，具有氣液的雙重特性：其黏度和擴(kuò)散系數(shù)接近氣體，因此具有良好的流動性能、滲透性能和傳遞性能[13]。

隨著對PCHE研究的不斷深入，以S-CO2為循環(huán)工質(zhì)的實驗和理論研究在流動范圍、應(yīng)用領(lǐng)域等方面仍然存在以下不足之處：通道內(nèi)S-CO2流動與傳熱局限于低雷諾數(shù)范圍，且冷側(cè)和熱側(cè)通道內(nèi)工質(zhì)均為S-CO2的研究尚不充分等[14]。為此，筆者以S-CO2為循環(huán)工質(zhì)，對板式擴(kuò)散焊矩形微通道換熱器內(nèi)的流動與傳熱特性展開數(shù)值模擬，建立了適用的流動與傳熱數(shù)學(xué)物理模型，探索不同Re條件下擾流格柵對微通道內(nèi)流動與傳熱特性的影響機(jī)理，從微通道的流動特性、傳熱特性及綜合性能等角度進(jìn)行分析，以提高微通道換熱器的流動換熱綜合性能，并為換熱器的結(jié)構(gòu)優(yōu)化提供理論依據(jù)。

1 板式擴(kuò)散焊矩形微通道換熱器的數(shù)學(xué)物理模型

1.1 物理模型及邊界條件

傳統(tǒng)PCHE的加工是通過在金屬板片上進(jìn)行光電化學(xué)刻蝕形成通道，其工藝復(fù)雜、成本高昂[15]。筆者選用的新型板式擴(kuò)散焊矩形微通道換熱器結(jié)構(gòu)如圖1所示。該換熱器利用激光切割316L不銹鋼金屬板片，隨后以擴(kuò)散焊方式將9片金屬板按照一定次序堆疊形成一組冷熱通道，其中第1、第5、第9片為冷熱流體換熱界面，第2、第3、第4和第6、第7、第8片分別形成熱流體通道和冷流體通道(以下分別簡稱為熱通道和冷通道)，每3片金屬板的厚度為矩形微通道的高度，而冷熱通道中格柵分別位于第3和第7片，堆疊而成的逆流換熱單元如圖1(b)所示，經(jīng)換熱單元周期性重復(fù)排列，最終組成一個完整的換熱器。

換熱器包含18個換熱單元，換熱器外觀尺寸長×寬×高為411.0 mm×112.0 mm×43.6 mm，其中矩形通道斷面尺寸為3.00 mm×2.85 mm，格柵高度為0.95 mm。格柵高度是本文換熱器結(jié)構(gòu)中的最小幾何尺寸(即板片的厚度)，本文數(shù)值模擬所使用的網(wǎng)格在壁面附近需設(shè)置細(xì)密的邊界層，因此不宜直接選用完整的換熱器或者換熱單元作為數(shù)值模擬對象，以免影響數(shù)值模擬的效率和精度。為此，對換熱單元的結(jié)構(gòu)進(jìn)行了簡化，建立了如圖1(c)所示的冷熱通道單元的幾何模型，其中矩形通道與格柵等的幾何尺寸均為換熱器的實際尺寸。

此外，建立了總長度和當(dāng)量直徑相同的半圓形直通道與矩形通道幾何模型，如圖2(a)所示，其中g(shù)為重力加速度。研究表明，在通道內(nèi)增設(shè)擾流格柵，可以起到強(qiáng)化傳熱的效果[16]。因此，在圖2(c)所示的矩形通道基礎(chǔ)上，設(shè)置了如圖2(d)～圖2(g)所示的不同的擾流格柵，探索其對換熱單元流動與傳熱特性的影響。圖2(e)中的冷熱通道幾何模型即為圖1(c)所示的冷熱通道中部設(shè)置水平格柵的換熱器幾何模型。各個模型的相關(guān)幾何尺寸見表1。

表1 冷熱通道結(jié)構(gòu)參數(shù)

利用Ansys Fluent 18.0軟件展開數(shù)值模擬研究，其中S-CO2熱物性參數(shù)來自于REFPROP V9.1，并通過引入實際氣體模型(NIST Real Gas)實現(xiàn)動態(tài)求解。數(shù)值模擬邊界條件見表2，其中計算域固體區(qū)域材料為316L不銹鋼，其密度、比定壓熱容和導(dǎo)熱系數(shù)分別為8 360 kg/m3、0.417 kJ/(kg·K)和0.021 kW/(m·K)。

表2 數(shù)值模擬邊界條件

1.2 流動與傳熱模型

板式擴(kuò)散焊矩形微通道中S-CO2流動與傳熱特性研究的數(shù)值模擬控制方程[17]如下：

連續(xù)性方程

div(ρu)=0

(1)

動量方程

(2)

ρg+Sv

(3)

(4)

能量方程

(5)

式中：ρ和μ分別為流體的密度和動力黏度；p和T分別為壓力和溫度；K和c分別為流體傳熱系數(shù)和比熱容；u、v和w分別為流體沿x、y和z方向的流速分量；S為廣義源項；μt為湍動黏度，本文選用SSTk-ω湍流模型進(jìn)行湍流計算，因此μt根據(jù)式(6)[18]計算。

μt=ρk/ω

(6)

式中：k和ω分別為湍動能和湍動頻率，計算方法詳見式(7)和式(8)。

Yk+Sk

(7)

Yω+Sω

(8)

式中：Gk為由層流流速梯度而產(chǎn)生的湍動能；Gω為由ω方程產(chǎn)生的湍動能；Γk和Γω分別為k和ω的有效擴(kuò)散率；Yk和Yω分別為k和ω的發(fā)散項；Sk和Sω為自定義源項。

數(shù)值模擬選用穩(wěn)態(tài)計算，求解方法采用基于壓力-流速耦合的SIMPLE算法，壓力項采用二階迎風(fēng)離散格式，梯度項采用格林-高斯離散格式，動量和能量均采用二階迎風(fēng)離散格式。

1.3 流動與傳熱特性評價指標(biāo)

選取努塞爾數(shù)(Nu)和流動摩擦阻力系數(shù)(f)表征S-CO2在不同通道內(nèi)的傳熱特性和流動阻力特性。Nu越大表明流體與壁面之間的傳熱系數(shù)越大；f越大表明流體在通道內(nèi)的流動損失越大[19]。Nu和f的計算方法如下：

(9)

(10)

式中：qw為監(jiān)測得到的冷熱通道壁面熱通量；De為通道的當(dāng)量直徑；λ為通道內(nèi)流體的導(dǎo)熱系數(shù)；Tf和Tw分別為流體和壁面的溫度；Δp為通道進(jìn)出口處流體之間的壓差；uc為流體流速。

為了兼顧傳熱特性和流動特性，以綜合傳熱增強(qiáng)因子(PEC)θ作為表征流動與傳熱特性的綜合指標(biāo)，計算方法見式(11)。當(dāng)PEC大于1時，表示綜合性能優(yōu)于比較基準(zhǔn)(以圖2(c)所示的無格柵矩形通道為比較基準(zhǔn)，其努塞爾數(shù)和流動摩擦阻力系數(shù)分別為Nu0和f0)。

(11)

2 模型驗證

2.1 網(wǎng)格無關(guān)性驗證

通過Ansys ICEM在幾何模型計算域生成六面體結(jié)構(gòu)化網(wǎng)格，網(wǎng)格質(zhì)量在0.9以上。為準(zhǔn)確反映通道內(nèi)流體真實的流動與傳熱狀態(tài)，在冷熱通道壁面處設(shè)置了細(xì)密的邊界層網(wǎng)格。邊界層網(wǎng)格的初始高度為0.001 6 mm，增長率為1.2，以確?？拷诿娴谝粚泳W(wǎng)格節(jié)點處的y+小于1。

為驗證網(wǎng)格無關(guān)性，選取65萬、91萬、126萬和214萬4種網(wǎng)格數(shù)展開數(shù)值模擬。其中，冷、熱通道S-CO2的入口質(zhì)量流量分別為2.022 4×10-3kg/s和1.305 2×10-3kg/s，其余參數(shù)見表2，數(shù)值模擬結(jié)果如圖3所示。

圖3 不同網(wǎng)格數(shù)下冷熱通道內(nèi)的流動與傳熱特性

由圖3可知，在相同的工況條件下，冷熱通道中的壓降和溫差均隨著網(wǎng)格數(shù)的增加而變化，并在網(wǎng)格數(shù)超過91萬時基本趨于穩(wěn)定，表明此時網(wǎng)格數(shù)增加對計算結(jié)果的影響可以忽略。因此，為提高計算效率，選用126萬網(wǎng)格方案進(jìn)行數(shù)值模擬。

2.2 模擬準(zhǔn)確性驗證

以Meshram等[20]的研究結(jié)果為參考標(biāo)準(zhǔn)，在相同的邊界條件下驗證所建物理模型的準(zhǔn)確性，幾何模型采用半圓形直通道模型(即model1)。驗證結(jié)果見表3，其中誤差為(模擬值-參考值)/參考值。

表3 模擬準(zhǔn)確性驗證結(jié)果

由表3可知，數(shù)值模擬得到的冷熱通道內(nèi)的壓降及溫差與參考值基本一致，壓降的誤差范圍為-7.52%～13.23%，溫差的誤差范圍為-4.30%～8.79%。綜上所述，所建立的冷熱通道換熱單元數(shù)學(xué)物理模型可以滿足S-CO2流動與傳熱特性數(shù)值模擬的準(zhǔn)確性要求。

3 結(jié)果與分析

3.1 無格柵時直通道內(nèi)S-CO2的流動與傳熱特性

首先研究了通道截面為半圓形(圖2(b))和矩形(圖2(c))時冷熱通道換熱單元中S-CO2的傳熱特性和流動阻力特性，Nu和f隨Re的變化情況如圖4所示。由圖4可知，當(dāng)5 000≤Re≤30 000時，隨著Re的增大，無格柵的冷熱通道內(nèi)Nu增大，而f則減小。這是由于當(dāng)Re增大時，通道內(nèi)流速升高、湍流加劇，增強(qiáng)了S-CO2與壁面之間的對流換熱，因而表現(xiàn)出更好的傳熱特性；同時，流速升高使得近壁面附近的層流邊界層變薄，通道中流動摩擦阻力相應(yīng)減小。對比相同Re工況下的數(shù)值模擬結(jié)果不難發(fā)現(xiàn)，截面為半圓形的直通道傳熱性能較弱但流動摩擦阻力較小，而截面為矩形的直通道傳熱性能較好，但流動摩擦阻力較大。

圖4 無格柵時直通道內(nèi)Nu和f隨Re的變化情況

圖5給出了無格柵半圓形直通道PEC隨Re的變化情況。由圖5可以看出，當(dāng)Re>5 000時，在無格柵情況下，半圓形冷熱直通道的PEC均略大于1，且在相同Re條件下，熱通道內(nèi)流速大于冷通道內(nèi)流速，因而熱通道的PEC均大于冷通道的PEC。此外，半圓形直通道內(nèi)的PEC隨著Re的增大先升高后降低，當(dāng)Re為25 000左右時，綜合傳熱效果最好。

圖5 無格柵半圓形直通道PEC隨Re的變化情況

綜上所述，半圓形直通道的綜合傳熱性能雖然優(yōu)于矩形通道，但優(yōu)勢并不明顯。因此，考慮到半圓形直通道加工工藝復(fù)雜、生產(chǎn)成本高昂，優(yōu)化矩形通道結(jié)構(gòu)、實現(xiàn)通道強(qiáng)化傳熱、提升綜合傳熱性能更具實際意義。

3.2 設(shè)置格柵時矩形通道內(nèi)S-CO2的流動與傳熱特性

為了分析通道內(nèi)熱流場情況，沿幾何模型Z軸方向選取了4個取樣面，取樣面與熱通道入口之間的距離分別為60 mm、175 mm、393 mm和508 mm，即第2個取樣面位于熱通道格柵下游3 mm處，第3個取樣面位于冷通道格柵下游3 mm處。圖6給出了Re為20 000時，3種擾流格柵設(shè)置方式對矩形通道內(nèi)溫度分布的影響。由圖6可知，在相同Re條件下，熱通道內(nèi)的溫度梯度大于冷通道內(nèi)的溫度梯度，且對比Z=175 mm和Z=393 mm處的溫度云圖不難發(fā)現(xiàn)，熱通道內(nèi)擾流格柵對溫度場的影響更加顯著。

圖6 矩形通道內(nèi)不同位置處的溫度分布云圖

Fig.6 Contours of temperature distribution at different positions in the rectangular channel

此外，由圖2(b)可以看出，將擾流格柵設(shè)置在通道中部，對通道壁面處溫度邊界層的影響較小；而將擾流格柵設(shè)置在通道兩側(cè)，對通道壁面處溫度邊界層的影響則較大。將圖6(a)與圖6(d)、圖6(e)對比可知，S-CO2經(jīng)過位于通道兩側(cè)的水平格柵后，其對側(cè)近壁處溫度邊界層變薄，近壁處流體與壁面之間的溫度梯度減小，實現(xiàn)了熱通道與冷通道之間的強(qiáng)化換熱，其中圖6(d)及圖6(e)同為兩側(cè)水平格柵模型，不同的是圖6(d)反映的是遠(yuǎn)離水平格柵換熱界面處的溫度分布，圖6(e)反映的是緊靠水平格柵換熱界面處的溫度分布?？紤]到湍動能是衡量流體混合能力的重要指標(biāo)，冷熱通道內(nèi)S-CO2在流經(jīng)擾流格柵后湍動能會發(fā)生急劇變化，以熱通道為例，圖7給出了通道中心截面處3種擾流格柵設(shè)置方式對熱通道內(nèi)S-CO2湍動能分布的影響。 以圖7(a)為比較基準(zhǔn)可知，圖7(b)及圖7(c)中在通道中設(shè)置擾流格柵可以顯著增強(qiáng)流體擾動、促進(jìn)流體混合，且圖7(d)中擾流格柵下游流體的湍動能峰值最高、湍動程度最強(qiáng)、湍動范圍最大。因此，將擾流格柵水平設(shè)置在通道兩側(cè)(即model5)時，擾流效果最好，可以有效改善近壁面流體與壁面之間的流動與傳熱。

圖7 矩形熱通道內(nèi)S-CO2的湍動能分布云圖

Fig.7 Contours of turbulent kinetic energy distribution of S-CO2 in the hot rectangular channel

圖8給出了Z=175 mm處矩形熱通道內(nèi)S-CO2的流速分布。由圖8可知，無格柵時矩形通道(即model2)內(nèi)S-CO2在X方向和Y方向的流速分布均呈梯形狀，近似于管內(nèi)完全發(fā)展的湍流流速分布；設(shè)置格柵后，隨著格柵設(shè)置方式的變化，矩形通道內(nèi)S-CO2在不同方向的流速分布呈不同V形駝峰狀。在相同邊界條件下，將擾流格柵水平設(shè)置在通道兩側(cè)(即model5)時，矩形通道內(nèi)S-CO2在X方向的平均流速和Y方向的峰值流速分別為2.71 m/s和3.38 m/s，均為最大值，可以有效促進(jìn)流體換熱，但同時也會影響矩形通道內(nèi)的阻力特性。

以無格柵矩形通道為比較基準(zhǔn)，對比分析了不同格柵設(shè)置方式對通道內(nèi)流動與傳熱特性的影響，如圖9所示。由圖9可知，在本文研究范圍內(nèi)，設(shè)置格柵后矩形通道內(nèi)的Nu隨著Re的增大而逐漸增大，而f則隨著Re的增大而逐漸減小，且隨著Re的持續(xù)增大，f的變化將趨于平緩，與圖4中無格柵時矩形通道內(nèi)Nu和f隨Re的變化趨勢相同。此外，隨著Re的增大，矩形通道內(nèi)的Nu/Nu0幾乎保持不變，f/f0則緩慢增大。當(dāng)Re相同時，將擾流格柵水平設(shè)置在通道兩側(cè)(即model5)，Nu/Nu0和f/f0相對較大，這意味著該通道的傳熱性能較好，但流動摩擦阻力較大；將擾流格柵水平或者豎直設(shè)置在通道中部(即model3或model4)，Nu/Nu0和f/f0差別不大，此時矩形通道的傳熱性能較差，但流動摩擦阻力也較小。

為了綜合對比不同通道中的流動與傳熱特性，設(shè)置格柵后矩形通道內(nèi)的PEC如圖10所示。由圖10可知，3種擾流格柵設(shè)置方式使得冷熱通道內(nèi)部的PEC均小于1，且隨著Re的增大，PEC整體上呈緩慢下降趨勢。與無格柵的矩形通道相比，設(shè)置格柵后雖然可以增強(qiáng)傳熱，但同時也會增大流動摩擦阻力，導(dǎo)致綜合傳熱性能降低，且流速越高則綜合傳熱性能越差。因此，矩形通道內(nèi)擾流格柵的優(yōu)化設(shè)計需在保障強(qiáng)化傳熱的基礎(chǔ)上，盡可能降低對流動摩擦阻力的影響。對比圖6～圖10不難看出，在本文研究范圍內(nèi)，將擾流格柵水平設(shè)置在矩形通道兩側(cè)的擾流效果最好，在單獨考慮傳熱特性時，其傳熱性能最優(yōu)。因此，之后將在此基礎(chǔ)上，進(jìn)一步探索擾流格柵間隔、個數(shù)以及形狀尺寸等因素對流動和傳熱特性的影響機(jī)理，實現(xiàn)通道綜合傳熱性能以及換熱器整體性能的優(yōu)化。

(a) Nu和Nu/Nu0隨Re的變化情況

圖10 設(shè)置格柵后矩形通道內(nèi)PEC隨Re的變化情況

4 結(jié) 論

(1) 隨著Re的增大，板式擴(kuò)散焊矩形微通道換熱器換熱單元冷熱直通道內(nèi)的S-CO2流速增大、湍流加劇，Nu隨之增大，f反而減小。

(2) 未設(shè)置擾流格柵時，與矩形通道相比，在相同Re條件下，半圓形直通道的傳熱性能較弱，同時流動摩擦阻力較小，綜合傳熱增強(qiáng)因子大于1，即相同當(dāng)量直徑的半圓形直通道的綜合傳熱性能優(yōu)于矩形通道。

(3) 設(shè)置擾流格柵后，在相同Re條件下，擾流格柵水平設(shè)置在矩形通道兩側(cè)時的擾流效果和傳熱性能最優(yōu)，綜合傳熱性能相對最高。

(4) 與半圓形直通道和無格柵矩形通道相比，目前的擾流格柵設(shè)置方式在增強(qiáng)傳熱的同時會增大流動摩擦阻力，導(dǎo)致綜合傳熱增強(qiáng)因子降低。為兼顧換熱器流動與傳熱特性以及生產(chǎn)制造成本，需以擾流格柵水平設(shè)置在通道兩側(cè)的換熱單元為基礎(chǔ)，進(jìn)一步深入探索擾流格柵的優(yōu)化設(shè)計。