沈佳飛　王波　周秋平　王雅亮　成金東

摘要：針對超臨界二氧化碳氣冷堆核電系統(tǒng)中 Z 形通道印刷電路板式換熱器的優(yōu)化設(shè)計，通過數(shù)值模擬研究轉(zhuǎn)折角對印刷電路板式換熱器中二氧化碳流動和換熱特性的影響規(guī)律，分別擬合出摩擦因子和努塞爾數(shù)的計算關(guān)聯(lián)式，分析傳熱過程的熱阻，并討論換熱器的綜合性能。結(jié)果表明，摩擦因子和對流換熱系數(shù)均隨轉(zhuǎn)折角的增大呈拋物線規(guī)律增長，在轉(zhuǎn)折角小于 20 °時增長較慢。導(dǎo)熱熱阻占總傳熱熱阻的 4.16%～16.02%，并隨二氧化碳質(zhì)量流率和轉(zhuǎn)折角的增大而升高，在印刷電路板式換熱器的傳熱計算中，不應(yīng)被忽略。隨著轉(zhuǎn)折角的增加，通道中努塞爾數(shù)的增長幅度小于摩擦因子的增長幅度，傳遞相同熱量的泵送功率增大，但所需換熱面積減小，換熱器的制造成本下降，實際應(yīng)用中需進一步通過技術(shù)經(jīng)濟分析以選取最佳轉(zhuǎn)折角。為了控制通道中流動阻力占進口壓力的比例，轉(zhuǎn)折角以不超過 20 °為宜。

關(guān)鍵詞：印刷電路板式換熱器；轉(zhuǎn)折角；二氧化碳；流動阻力；傳熱

中圖分類號：? TK 172???????????? 文獻標(biāo)志碼：?? A

Effect of turning angle on flow and heat transfer characteristics of carbon dioxide in zigzag printed circuit heat exchangers

SHENJiafei1，WANG Bo1，ZHOU Qiuping2，WANGYaliang1，CHENG Jindong1

（1. School of Energy and Power Engineering， University of Shanghai for Science and Technology， Shanghai 200093， China;2. Zhejiang Dun'an Machinery and Electronics Technology Co.， Ltd.， Shaoxing 311800， China）

Abstract： Toachievetheoptimaldesignofzigzagprintedcircuitheatexchanger （PCHE） inthe supercritical carbon dioxide gas-cooled reactor nuclear power system， the influence of the turning angle on the flow and heat transfer characteristics of carbon dioxide in PCHEs was studied by numerical simulation，andthecorrelationsforFanningfrictionfactorandNusseltnumbercalculationwere obtained. The thermal resistance of the heat transfer process was analyzed， and the comprehensive performance of the heat exchanger was discussed. The results show that the Fanning friction factor and convective heat transfer coefficient both increase parabolically with the increase of the turning angle， and the growth is relatively slow when the turning angle is less than 20°. The thermal conductionresistance accounts for 4.16% ～? 16.02%? of the overall thermal resistance， and it increases with the increase of the mass flow rate and the turning angle. It should not be ignored in the heat transfer calculation of printed circuit heat exchanger. With the increase of the turning angle， the increase of the Nusselt number is smaller than that of the Fanning friction factor. The pumping power for transferring the same heat increases， but the required heat exchange area decreases. The manufacturing cost are reduced， and in practical applications， further techno-economic analysis is required to select the optimal turning angle. To control the proportion of the flow resistance in the channel to the inlet pressure， the turning angle should not exceed 20° .

Keywords： printed circuit heat exchanger ; turning? angle ; carbon? dioxide; flow? resistance; heat transfer

印刷電路板式換熱器（ printedcircuitheatexchanger ，PCHE）是一種高效緊湊式換熱器，工作溫度范圍廣，承壓能力可達60 MPa，在相同換熱條件下比傳統(tǒng)管殼式換熱器的體積減小約85%，有利于降低傳熱過程的溫差和不可逆損失，可用于太陽能熱發(fā)電、核電和船舶液化天然氣供應(yīng)等系統(tǒng)中[1-3]。PCHE 芯體由多層蝕刻有流體通道的金屬板片經(jīng)擴散焊接制成，通道型式是影響 PCHE 性能的關(guān)鍵因素，有直通道、Z 形通道、 S 形翅片通道和翼型翅片通道等結(jié)構(gòu)。 Chai 等[4]通過數(shù)值模擬研究了半圓形截面直通道 PCHE 內(nèi)超臨界二氧化碳（supercritical carbon dioxide ， S-CO2）的流動和換熱特性，發(fā)現(xiàn)局部對流放熱系數(shù)從入口處開始沿流動方向迅速下降，然后保持穩(wěn)定，存在明顯的入口效應(yīng)。李磊等[5]比較了不同通道 PCHE 中氦氣的流動特性，發(fā)現(xiàn)半圓形截面 Z 形通道 PCHE 的阻力系數(shù)大于直通道。 Chu 等[6]對 Z 形通道 PCHE 中 S-CO2和水的換熱進行了實驗研究，獲得了努塞爾數(shù)和摩擦系數(shù)的計算關(guān)聯(lián)式。謝瑤等[7]研究了 PCHE 中液化天然氣的傳熱系數(shù)與流動阻力隨質(zhì)量流率、熱通量和溫度的變化規(guī)律。 Cheng 等[8]研究了 Z 形通道 PCHE 預(yù)冷器中 S- CO2和水的換熱，分析了工質(zhì)流量和進口溫度對流動阻力、換熱系數(shù)、溫度效率的影響。 Wen 等[9] 研究了波形通道 PCHE 中，波形通道的波幅和波長對 S-CO2流動和傳熱性能的影響，結(jié)果表明，相對較小的波幅有利于提高綜合性能。宋丹等[10] 對超臨界液化天然氣在錯列 S 形通道 PCHE 中的對流換熱過程進行了數(shù)值模擬，分析了擬臨界工況附近的對流換熱特征。 Ma 等[11]研究了翅片橫向節(jié)距、縱向節(jié)距和有無端壁圓角對翼型翅片通道 PCHE 熱工水力性能的影響。針對典型應(yīng)用場景，氦氣[12-14]、超臨界甲烷[15-16]、超臨界二氧化碳[17-19]和熔鹽[20-21]等工質(zhì)在 PCHE 中的流動和換熱特性研究受到廣泛關(guān)注。

在采用 S-CO2冷卻反應(yīng)堆并直接通過 S-CO2布雷頓循環(huán)發(fā)電的新一代氣冷堆核電系統(tǒng)方案中，普遍采用 PCHE作為回?zé)崞?，其?yōu)化設(shè)計有利于提高循環(huán)效率。在不同通道結(jié)構(gòu)的 PCHE 中， Z 形通道 PCHE制造難度和成本較低，傳熱性能和應(yīng)用前景較好[22]。通道轉(zhuǎn)折角是 Z 形通道 PCHE 的重要結(jié)構(gòu)參數(shù)，對 PCHE 性能的影響尚需進一步研究。因此，本文針對 S-CO2氣冷堆核電系統(tǒng)中低溫回?zé)崞鞯膫鳠釂栴}，通過數(shù)值模擬分析不同轉(zhuǎn)折角 Z 形通道 PCHE 中二氧化碳的流動和換熱特性，為回?zé)崞鞯膬?yōu)化設(shè)計和性能實驗提供參考。

1 物理模型和數(shù)值方法

1.1 物理模型

圖1（a）為逆流換熱 PCHE 整體結(jié)構(gòu)示意圖，逆流換熱段芯體的直通道結(jié)構(gòu)如圖1（b）所示，Z 形通道結(jié)構(gòu)如圖1（c）所示。冷流體和熱流體通道均為半圓形截面，上下分層交替布置，如圖1（d）所示。通道截面半徑1.0 mm，水平方向節(jié)距3.2 mm，高度方向節(jié)距2.2 mm；每個 Z 形通道沿流動方向節(jié)距為18 mm ，共15個節(jié)距，總節(jié)距270 mm。 Z 形通道轉(zhuǎn)折角分別為0（直通道），5，10，15，20，30，45°。換熱器材料為316L不銹鋼，導(dǎo)熱系數(shù)為16.2 W/（m?K）。

1.2 網(wǎng)格模型

鑒于 Z 形通道 PCHE 通道布置的周期性特征，并假定熱流體和冷流體在各自流動通道中的流量分配均勻，選取圖2所示雙層通道換熱單元為數(shù)值模擬對象，熱通道在上，冷通道在下，不同參數(shù)的二氧化碳逆向流過上下通道進行換熱，固體域截面寬3.2 mm 、高4.4 mm、總節(jié)距270 mm。圖2（b）為其中1個節(jié)距范圍內(nèi)的三維結(jié)構(gòu)。采用 ICEM 軟件對固體域和流體域進行結(jié)構(gòu)化網(wǎng)格的劃分，半圓形截面通道采用 O 形網(wǎng)格，流體通道壁面附近進行網(wǎng)格加密。

在數(shù)值模擬過程中，對網(wǎng)格進行無關(guān)性驗證。圖3為通道轉(zhuǎn)折角=30° ， 冷通道和熱通道中二氧化碳質(zhì)量流率 G=200 kg/（m2·s），采用不同網(wǎng)格尺寸時，冷通道中流體壓降和對流換熱系數(shù)的計算結(jié)果。當(dāng)網(wǎng)格數(shù)達到6191874以后，壓降和對流換熱系數(shù)的計算結(jié)果變化十分緩慢，網(wǎng)格數(shù)進一步增加到10751736，計算結(jié)果相對偏差僅為0.16%和0.51%。因此，該工況下的模擬應(yīng)選取網(wǎng)格數(shù)量為6191874的網(wǎng)格模型。

1.3 邊界條件

PCHE 固體域上下壁面和左右壁面設(shè)置為周期性邊界條件，進口端和出口端壁面設(shè)置為絕熱邊界條件，流體通道的壁面為流固耦合界面。流體域入口采用質(zhì)量流量入口，出口采用壓力出口。熱側(cè)二氧化碳入口溫度和壓力分別為553.15 K 和3.2 MPa，冷側(cè)二氧化碳入口溫度和壓力分別為381.15 K 和10.5 MPa。冷側(cè)和熱側(cè)的入口質(zhì)量流率相等，質(zhì)量流率的范圍為100～400 kg/（m2?s）。

1.4 物性參數(shù)

二氧化碳的熱物理性質(zhì)隨壓力和溫度變化。由于 PCHE 進口二氧化碳的壓力高，進出口壓差和進口壓力的比值較小，故近似認為熱通道和冷通道中二氧化碳的物性參數(shù)僅隨溫度變化。通過 REFPROP 軟件分別查取10.5 MPa 和3.2 MPa 壓力下，350～600 K 范圍內(nèi)二氧化碳的密度ρ、動力黏度μ、定壓比熱容 cp 和導(dǎo)熱系數(shù)λ， 并擬合為關(guān)于溫度 T 的函數(shù)，結(jié)果如表1所示，擬合優(yōu)度 R2均高于0.97。

1.5 計算模型與求解方法

在本文研究的工況范圍內(nèi)，二氧化碳在 Z 形通道中流動的雷諾數(shù)均大于2000，湍流模型采用 k-ε模型，近壁區(qū)采用標(biāo)準(zhǔn)壁面函數(shù)。采用 SIMPLE 算法實現(xiàn)速度與壓力的耦合求解，變量梯度采用 Least Squares Cell Based 方法求解，壓力方程采用標(biāo)準(zhǔn)離散格式，動量方程和能量方程的離散采用二階迎風(fēng)格式，湍流方程的離散采用一階迎風(fēng)格式，所有變量收斂精度為10?6。該數(shù)值模擬方法已廣泛用于 PCHE 中二氧化碳傳熱特性研究，數(shù)值模擬與實驗結(jié)果吻合較好[14，23]。Wang等[24]也對本文的模型進行了驗證，分析表明，計算模型比較合理。

1.6 數(shù)據(jù)處理

流體通道的水力直徑

式中： V 表示流體通道的體積；A 表示流體通道的壁面積。

范寧摩擦因子

式中：ΔPf表示流體的壓降；表示流體平均密度； L 表示流體通道的長度； u表示通道有效截面上的平均速度。

流體和通道壁面間的平均表面對流換熱系數(shù)

式中： q 表示平均熱流密度；ΔTm表示壁面和流體的平均溫差。

流體和通道壁面間對流換熱的努塞爾數(shù)

式中：λ表示流體的平均導(dǎo)熱系數(shù)。

本文研究的 PCHE 中，熱通道和冷通道結(jié)構(gòu)完全相同，大端溫差ΔTd和小端溫差ΔTx 的比值小于2，ΔTm可近似采用算術(shù)平均溫差，即

式中： Tw 表示平均壁溫；Tf表示流體平均溫度。

以任一側(cè)流體通道壁面積為基準(zhǔn)，熱流體和冷流體通過金屬壁傳熱的總傳熱系數(shù)

式中：Tf，h為熱流體平均溫度；Tf，c為冷流體平均溫度。

在數(shù)值模擬中不考慮污垢熱阻，故導(dǎo)熱熱阻

式中：hh， hc分別為熱通道與冷通道的平均表面對流換熱系數(shù)。

Z 形通道 PCHE 采用不同通道轉(zhuǎn)折角，會影響流體和通道壁面間的平均表面換熱系數(shù)，以及流體的流動阻力，所以，采用綜合評價因子ζ評價Z 形通道強化換熱的綜合性能

ζ=ΔP（ΔPfr） f = Nf（u）/（/）fref（Nur）ef（8）

式中：下標(biāo) ref 表示參比通道，本文選取轉(zhuǎn)折角為0°的直通道為參比通道。

在相同的流體進口質(zhì)量流率和溫度下，ζ可以反映傳遞相同熱量時不同轉(zhuǎn)折角 PCHE 中流體流動阻力的高低[25]。

2 結(jié)果與分析

2.1PCHE 中二氧化碳的流動特性

2.1.1 PCHE 中二氧化碳的壓力分布

為了探討轉(zhuǎn)折角對二氧化碳流動和換熱特性的影響，對直通道 PCHE 和轉(zhuǎn)折角為5，10，15，20，30，45°的每一種 Z 形通道 PCHE，分別模擬二氧化碳質(zhì)量流率為100，150，200，250，300，350，400 kg/（m2?s）的7個工況。模擬結(jié)果表明，在所有49個工況中，基于每個通道內(nèi)二氧化碳平均熱物性計算所得流動雷諾數(shù)最低為5031，最高為20370，均處于湍流狀態(tài)。

熱通道和冷通道中二氧化碳的壓力沿流動方向都近似線性變化。以轉(zhuǎn)折角30°、質(zhì)量流率250 kg/（m2·s）的工況為例，冷通道和熱通道中二氧化碳的壓力分布分別如圖4（a）和（b）所示。從圖中可以看出，冷通道中二氧化碳的壓力從10.5 MPa 線性下降到10.4874 MPa，流動阻力12.6 kPa，流動阻力和進口壓力的比值為1.2‰；熱通道中二氧化碳的壓力從3.2 MPa線性下降到3.1546 MPa，流動阻力45.4 kPa，流動阻力和進口壓力的比值為14.2‰。

2.1.2轉(zhuǎn)折角對二氧化碳流動壓降的影響

通道中二氧化碳的壓降隨質(zhì)量流率和轉(zhuǎn)折角的增大而增大。如圖5所示，在400 kg/（m2·s）的質(zhì)量流率下，對轉(zhuǎn)折角為0，5，10，15，20、30，45°的 PCHE，冷通道二氧化碳的壓降分別為6.881，7.035，7.811，9.947，13.796，29.592，86.581 kPa，與冷通道二氧化碳進口壓力的比值分別為0.66‰，0.67‰，0.74‰，0.95‰，1.31‰，2.82‰，8.25‰；熱通道二氧化碳的壓降分別為23.889，24.483，27.234，34.952，50.552，107.179，316.088 kPa ，與熱通道二氧化碳進口壓力的比值分別為7.46‰，7.65‰，8.51‰，10.92‰，15.80‰，33.49‰，98.78‰；熱通道二氧化碳的壓降是冷通道的3.47～3.65倍。由此可見，冷通道和熱通道結(jié)構(gòu)相同、質(zhì)量流率相等的條件下，冷通道二氧化碳進口壓力高，流動阻力小，而熱通道二氧化碳進口壓力低，流動阻力大，在 PCHE 設(shè)計中可增大熱通道流通面積以控制流動阻力。

2.1.3轉(zhuǎn)折角對流動通道摩擦因子的影響

二氧化碳在 PCHE 通道中流動的摩擦因子隨轉(zhuǎn)折角和質(zhì)量流率變化的規(guī)律如圖6所示。在每個模擬工況中，熱通道和冷通道的轉(zhuǎn)折角和質(zhì)量流率相同，摩擦因子也接近，而熱通道中二氧化碳的流動阻力高，主要是由熱通道中二氧化碳的壓力低、比容大引起的。對一定轉(zhuǎn)折角的 Z 形通道，由于流動進入紊流光滑管區(qū)，摩擦因子與雷諾數(shù)負相關(guān)，質(zhì)量流率越大，雷諾數(shù)越大，摩擦因子越小。在相同的質(zhì)量流率下，Z 形通道轉(zhuǎn)折角越大，摩擦因子越高，這是因為隨著轉(zhuǎn)折角的增大，流體通過轉(zhuǎn)折角處時方向變化更急，流動截面上速度不均勻性更大，流動分離加劇，流體對前方壁面的碰撞增強，再加上分離區(qū)存在二次流，導(dǎo)致能量耗散，增加局部阻力[26]。從圖 6 中還可以看出，摩擦因子隨轉(zhuǎn)折角呈拋物線規(guī)律變化。在一定的質(zhì)量流率下，直通道的摩擦因子最??；當(dāng) Z 形通道轉(zhuǎn)折角小于 20?°時，摩擦因子隨轉(zhuǎn)折角的增大而緩慢增大；當(dāng)轉(zhuǎn)折角大于 20?°時，隨著轉(zhuǎn)折角的增大，摩擦因子大幅增長。以質(zhì)量流率400 kg/（m2·s）的工況為例，直通道 PCHE熱通道中二氧化碳平均密度34.79 kg/m3，摩擦因子為0.01176，流動阻力為23.889 kPa；轉(zhuǎn)折角20°的 Z 形通道 PCHE 熱通道中二氧化碳平均密度34.99 kg/m3，摩擦因子增大到0.02351，流動阻力增大到 50.552 kPa ；當(dāng)轉(zhuǎn)折角為45°時，熱通道中二氧化碳平均密度35.32 kg/m3，摩擦因子達到0.11165，流動阻力高達316.088 kPa。流動阻力增加的幅度顯著超過摩擦因子增加的幅度，這是因為在相同的總節(jié)距下，流動通道總長度和轉(zhuǎn)折角的余割值成正比，轉(zhuǎn)折角越大，不僅摩擦因子增加，流動通道也越長；轉(zhuǎn)折角為45°時，在相同的總節(jié)距下，其流動通道長度是直通道的1.414倍。在不同轉(zhuǎn)折角下，熱通道內(nèi)二氧化碳的平均密度略有變化，是影響熱通道流動阻力的次要因素。因此，將轉(zhuǎn)折角控制在20°以下，有利于避免過大的流動阻力。

通道內(nèi)二氧化碳流動的摩擦因子與雷諾數(shù)的關(guān)系如圖7所示。在本文研究的每一個工況下，冷通道和熱通道中二氧化碳的質(zhì)量流率相同，盡管冷通道中二氧化碳壓力高、溫度低，熱通道中二氧化碳壓力低、溫度高，兩側(cè)流體的動力黏度卻很接近，相對偏差都在3%以內(nèi)，因而兩側(cè)流動的雷諾數(shù)也非常接近。在相同的雷諾數(shù)下，冷通道中二氧化碳流動的摩擦因子比熱通道高約3%～14%；轉(zhuǎn)折角越大，兩側(cè)的摩擦因子相差越小。這可能是因為冷通道中二氧化碳溫度沿流動方向不是線性變化，進口端二氧化碳的溫度上升較快，出口端溫度上升較緩，按平均溫度查取密度后計算得到的平均速度和動壓略微偏小，使摩擦因子偏大一點。

2.1.4流動通道摩擦因子與雷諾數(shù)的關(guān)聯(lián)式

對摩擦因子和雷諾數(shù)進行非線性擬合

f = C1ReC2（9）

得到不同轉(zhuǎn)折角 Z 形通道 PCHE 中冷通道和熱通道摩擦因子與雷諾數(shù)之間關(guān)聯(lián)式的系數(shù) C1和 C2，如表2所示。

2.2PCHE 中二氧化碳的換熱特性

2.2.1 PCHE 中二氧化碳的溫度分布

各工況下熱通道中二氧化碳的溫度沿流動方向近似線性下降，冷通道中二氧化碳在入口端沿流動方向溫度上升較快，然后平穩(wěn)增長；冷通道入口端的傳熱溫差較大，往出口端逐漸減小，并在較大范圍內(nèi)保持基本穩(wěn)定。以轉(zhuǎn)折角30°、質(zhì)量流率250 kg/（m2·s）條件下 Z 形通道 PCHE 內(nèi)的流動為例，如圖8所示，冷通道中二氧化碳溫度從 381.15 K 上升到532.35K，熱通道中二氧化碳溫度從553.15 K 下降到430.41K，高溫端傳熱溫差20.8 K，低溫端傳熱溫差49.26 K。由于冷側(cè)進口二氧化碳比定壓熱容為1.490 kJ/（kg?K），熱側(cè)出口二氧化碳比定壓熱容為1.035 kJ/（kg?K），兩側(cè)質(zhì)量流率相等，熱側(cè)和冷側(cè)交換一定的熱量后，引起的熱側(cè)二氧化碳溫降，應(yīng)高于冷側(cè)二氧化碳的溫升，而圖中溫度分布表明低溫端冷側(cè)二氧化碳溫度反而上升快，傳熱溫差大，這意味著在 PCHE 中由于流動通道的水力直徑很小，冷熱通道之間的金屬壁從高溫端向低溫端的軸向?qū)嶙饔貌豢珊雎?。結(jié)果表明，在回?zé)崞髟O(shè)計中如要進一步減小傳熱溫差，還需增加通道長度。

2.2.2轉(zhuǎn)折角對二氧化碳和通道壁面間對流換熱的影響

二氧化碳與 PCHE通道壁面間的對流換熱系數(shù)隨通道轉(zhuǎn)折角和質(zhì)量流率的增大而增大。從圖9中可以看出，對流換熱系數(shù)隨轉(zhuǎn)折角呈拋物線規(guī)律增大，在一定的質(zhì)量流率下，直通道的對流換熱系數(shù)最??；當(dāng) Z 形通道轉(zhuǎn)折角小于20°時，對流換熱系數(shù)隨轉(zhuǎn)折角的增大而緩慢增大；當(dāng)轉(zhuǎn)折角大于20°時，隨著轉(zhuǎn)折角的增大，對流換熱系數(shù)的增幅加大。以質(zhì)量流率400 kg/（m2·s）的工況為例，直通道 PCHE熱通道中二氧化碳和通道壁面間對流換熱系數(shù)為3268.4 W/（m2·K），轉(zhuǎn)折角20°的 Z 形通道 PCHE熱通道中二氧化碳和通道壁面間對流換熱系數(shù)增大到4004.6 W/（m2·K），轉(zhuǎn)折角45°的 Z 形通道 PCHE熱通道中二氧化碳和通道壁面間對流換熱系數(shù)進一步提高到6724.8 W/（m2·K）。這是由于轉(zhuǎn)折角增大時，流體經(jīng)過轉(zhuǎn)折角處的流向改變更大，流動分離加劇，對前方壁面的沖刷增強，通過流體擾動更充分地破壞邊界層，從而增強與壁面的對流換熱。結(jié)合圖6可知，轉(zhuǎn)折角增大到45°， 對流換熱系數(shù)增長幅度遠遠低于摩擦因子增長的幅度，這會引起泵送功率的提高。

從圖10中可以看出，努塞爾數(shù)隨通道轉(zhuǎn)折角和雷諾數(shù)的增大而增大。在每個工況中，冷通道和熱通道的轉(zhuǎn)折角相等，結(jié)構(gòu)完全相同，但是，在相同的雷諾數(shù)下，熱通道中努塞爾數(shù)比冷通道高出40%～150%，而且雷諾數(shù)越小，轉(zhuǎn)折角越大，熱通道中努塞爾數(shù)比冷通道中努塞爾數(shù)增加的幅度越大。

對努塞爾數(shù)和雷諾數(shù)、普朗特數(shù)進行非線性擬合：

Nu = C1ReC2Pr n（10）

式中：Pr為普朗特數(shù)； n 為系數(shù)，對熱通道， n 取0.3，對冷通道， n 取0.4。

不同轉(zhuǎn)折角 Z 形通道 PCHE 中冷通道和熱通道中努塞爾數(shù)關(guān)聯(lián)式的系數(shù)，如表3所示。

2.2.3轉(zhuǎn)折角對導(dǎo)熱熱阻的影響

本文研究的 PCHE 中上下2個相鄰的半圓形截面冷熱通道之間的壁厚是漸變的，可采用等效壁厚法計算冷熱通道之間的導(dǎo)熱熱阻[1，27]。等效壁厚

tw = t -（11）

式中： t 為冷熱通道之間的最大壁厚，對本文研究的 PCHE，其值為0.0022 m；r 為半圓形截面的半徑。

導(dǎo)熱熱阻

rw =（12）

式中：λw 為金屬的導(dǎo)熱系數(shù)。

采用等效壁厚法計算得到冷熱通道之間金屬壁的導(dǎo)熱熱阻為8.7321×10?5 m2·K/W。

通過數(shù)值模擬獲得的總傳熱系數(shù)隨轉(zhuǎn)折角和質(zhì)量流率的變化規(guī)律如圖11所示。

根據(jù)總傳熱系數(shù)、熱通道的平均表面對流換熱系數(shù)和冷通道的平均表面對流換熱系數(shù)，按照式（7）計算可得傳熱過程的導(dǎo)熱熱阻。同一轉(zhuǎn)折角下，導(dǎo)熱熱阻隨質(zhì)量流率的變化很小，近似為一個定值。對于轉(zhuǎn)折角0，5，10，15，20，30，45°的Z 形通道 PCHE，計算所得導(dǎo)熱熱阻分別為7.8522×10?5，7.8283×10?5，7.7639×10?5，7.7213×10?5，7.7565×10?5，7.8764×10?5，7.9214×10?5 m2·K/W，與采用等效壁厚法計算所得導(dǎo)熱熱阻的相對偏差分別為10.08%，10.35%，11.09%，11.58%，11.17%，9.80%，9.28%。等效壁厚法計算所得導(dǎo)熱熱阻偏高，其原因可能是等效熱阻法是簡化計算方法，沒有考慮轉(zhuǎn)折角對導(dǎo)熱的影響，也沒有考慮冷熱通道之間的金屬壁從高溫端向低溫端的軸向?qū)嶙饔谩?/p>

在相同質(zhì)量流率下，冷通道中二氧化碳的對流換熱熱阻在總傳熱熱阻中所占比重較大。以質(zhì)量流率為100 kg/（m2·s）的工況為例，轉(zhuǎn)折角為0，5，10，15，20，30，45°時，冷通道對流換熱熱阻分別占總傳熱熱阻的60.04%，60.07%，60.23%，60.88%，61.53%，63.17%，65.31%。隨著質(zhì)量流率的增大，冷通道熱阻在總熱阻中所占比重略有減小。如圖12所示，質(zhì)量流率和轉(zhuǎn)折角越大，冷通道和熱通道的對流換熱熱阻越小，導(dǎo)熱熱阻占總傳熱熱阻的比例就越高。在所有的工況中，導(dǎo)熱熱阻占總傳熱熱阻的比例在4.16%～16.02%之間，表明在 PCHE 的傳熱計算中，不宜忽略導(dǎo)熱熱阻。2.3 轉(zhuǎn)折角對 PCHE 綜合性能的影響

通過上述分析可知，隨著 Z 形通道轉(zhuǎn)折角的增大， PCHE 換熱性能增強，同時流動阻力也會增大。為了綜合評價轉(zhuǎn)折角對 Z 形通道 PCHE 性能的影響，以直通道 PCHE 作為參比結(jié)構(gòu)，按式（8）計算綜合評價因子ζ。

圖13為不同轉(zhuǎn)折角 PCHE 的綜合評價因子隨質(zhì)量流率的變化規(guī)律。從圖中可以看出，所有轉(zhuǎn)折角下 PCHE 的綜合評價因子均小于1，表明摩擦因子的增長幅度大于努塞爾數(shù)的增長幅度。在相同的冷熱流體質(zhì)量流率和進出口溫度下，與直通道 PCHE 相比，采用 Z 形通道 PCHE 時實現(xiàn)相同傳熱量，克服流動阻力所消耗的功率隨著轉(zhuǎn)折角的增大而升高，導(dǎo)致運行成本增加。但是，隨著轉(zhuǎn)折角的增大，努塞爾數(shù)提高，所需換熱面積減小，換熱器結(jié)構(gòu)更緊湊，重量更輕，這將減小換熱器的制造成本。在 PCHE 設(shè)計時，需綜合考慮制造成本和運行成本，以及對發(fā)電循環(huán)效率的影響。結(jié)合圖5所示冷通道和熱通道流動阻力占冷通道進口壓力的比例，轉(zhuǎn)折角以不超過20°為宜，否則過大的流動阻力會降低布雷頓發(fā)電循環(huán)的效率。

3 結(jié)論

針對 S-CO2氣冷堆核電系統(tǒng)中低溫回?zé)崞鞯膫鳠?，通過數(shù)值模擬分析不同轉(zhuǎn)折角 Z 形通道 PCHE 中二氧化碳的流動和換熱特性，得出如下結(jié)論：

a.冷通道和熱通道內(nèi)二氧化碳流動和換熱過程中，摩擦因子和對流換熱系數(shù)均隨轉(zhuǎn)折角的增大呈拋物線規(guī)律變化，當(dāng) Z 形通道轉(zhuǎn)折角小于20°時，摩擦因子和對流換熱系數(shù)隨轉(zhuǎn)折角的增大而緩慢增大；當(dāng)轉(zhuǎn)折角大于20°時，隨著轉(zhuǎn)折角的增大，摩擦因子和對流換熱系數(shù)增長幅度加大，且對流換熱系數(shù)增長幅度顯著低于摩擦因子增長幅度。

b.在 PCHE冷通道和熱通道結(jié)構(gòu)相同、質(zhì)量流率相等的條件下，熱通道中二氧化碳的流動阻力和對流換熱系數(shù)顯著高于冷通道，可考慮通過增加熱側(cè)二氧化碳流通截面積，降低熱通道中二氧化碳流速和阻力。

c.質(zhì)量流率和轉(zhuǎn)折角越大，導(dǎo)熱熱阻占總傳熱熱阻的比例越高。在本文研究的工況中，導(dǎo)熱熱阻占總傳熱熱阻的比例在4.16%～16.02%之間，在 PCHE 的傳熱計算中，不宜忽略導(dǎo)熱熱阻。采用等效壁厚法計算所得的導(dǎo)熱熱阻相對偏大。

d.不同轉(zhuǎn)折角下 PCHE 的綜合評價因子均小于1，意味著在相同的質(zhì)量流率和進出口溫度下，隨著轉(zhuǎn)折角的增加，傳遞相同熱量的泵送功率增大，但所需換熱面積減小，換熱器體積、重量和制造成本下降，在實際應(yīng)用中應(yīng)通過技術(shù)經(jīng)濟分析以選取最佳轉(zhuǎn)折角。結(jié)合通道中流動阻力占進口壓力的比例，轉(zhuǎn)折角以不超過20°為宜，以免流動阻力過大導(dǎo)致布雷頓發(fā)電循環(huán)效率下降。

參考文獻：

[1]徐哲， 張明輝， 段天應(yīng)， 等.超臨界二氧化碳在印刷電路板式換熱器內(nèi)的流動換熱特性研究[J].原子能科學(xué)技術(shù)， 2021， 55（5）：849–855.

[2] KLEMES? J J. Compact heat exchangers for energy transfer intensification： low grade heat and fouling mitigation[M]. Boca Raton： CRC Press， 2016：101–104.

[3] NIKITINK，KATOY，NGOL. Printedcircuitheat exchanger thermal-hydraulic performanceinsupercritical CO2experimentalloop[J].InternationalJournalofRefrigeration， 2006， 29（5）：807–814.

[4] CHAIL，TASSOUSA. Numericalstudyofthe thermohydraulicperformanceofprintedcircuitheat exchangersforsupercriticalCO2Braytoncycle applications[J]. Energy Procedia， 2019， 161：480–488.

[5]李磊， 楊劍， 馬挺，等.印刷電路板通道的高溫傳熱和阻力特性研究[J].工程熱物理學(xué)報， 2014， 35（5）：931–934.

[6] CHU W X， LI X H， CHEN Y T， et al. Experimental study on small scale printed circuit heat exchanger with zigzag channels[J]. HeatTransferEngineering， 2021， 42（9）：723–735.

[7]謝瑤， 李劍銳， 胡海濤.印刷電路板式換熱器內(nèi)超臨界甲烷流動換熱特性模擬[J].化工學(xué)報 ， 2021， 72（S1）：203–209.

[8] CHENGKY，ZHOUJZ，ZHANGHZ，etal. Experimentalinvestigationofthermal-hydraulic characteristics of a printed circuit heat exchanger used as a pre-coolerforthesupercriticalCO2Braytoncycle[J]. Applied Thermal Engineering， 2020， 171：115116.

[9] WEN Z X， LV Y G， LI Q， et al. Numerical study on heat transfer behavior of wavy channel supercritical CO2 printed circuitheatexchangerswithdifferentamplitudeand wavelengthparameters[J]. InternationalJournalofHeat and Mass Transfer， 2020， 147：118922.

[10]宋丹， 蔣慶峰， 馮國增， 等.超臨界 LNG 在錯列 S 形翅片微通道的流動傳熱特性研究[J].低溫與超導(dǎo) ， 2021，49（1）：16–21.

[11] MA T， XIN F， LI L， et al. Effect offin-endwall fillet on thermal hydraulic performance of airfoil printed circuit heat exchanger[J]. AppliedThermalEngineering， 2015， 89：1087–1095.

[12] MCCORMACK D. The application of printed circuit heat exchanger technology in the pebble bed modular reactor demonstration plant[C]//Proceedings of ASME Turbo Expo 2001： Power for Land， Sea， and Air. New Orleans： ASME， 2001：2001-GT-0473.

[13] KIMIH，NOHC. Physicalmodeldevelopmentand optimal design of PCHE for intermediate heat exchangers in HTGRs[J]. Nuclear Engineering and Design， 2012， 243：243–250.

[14] KIMDE，KIMMH，CHAJE，etal. Numerical investigationonthermal –hydraulicperformanceofnew printedcircuitheatexchangermodel[J].Nuclear Engineering and Design， 2008， 238（12）：3269–3276.

[15]李瑋哲， 林文勝.超臨界甲烷在印刷電路板換熱器中加熱過程模擬[J].低溫工程， 2017（5）：60–64，74.

[16]賈丹丹， 趙忠超， 張永， 等.超臨界 LNG 在印刷板式汽化器微細流道內(nèi)的流動與換熱性能數(shù)值研究[J].船舶工程， 2017， 39（5）：35–40.

[17] BAEKS，KIMJH，JEONGS，etal. Developmentof highlyeffectivecryogenic printedcircuitheatexchanger （PCHE） with low axial conduction[J]. Cryogenics， 2012， 52（7/9）：366–374.

[18]張海燕， 郭江峰， 淮秀蘭， 等. PCHE 內(nèi)軸向?qū)釋植繐Q熱性能的影響研究 [J].化工學(xué)報 ， 2019， 70（12）：4590–4598.

[19]王績德， 馮巖， 韓東江.超臨界二氧化碳動力循環(huán)在鈉冷快堆中的應(yīng)用綜述[J].核科學(xué)與工程， 2019， 39（2）：289–297.

[20] SHI H Y， LI M J， WANG W Q， et al. Heat transfer and friction of molten salt and supercritical CO2 flowing in an airfoilchannelofaprintedcircuitheatexchanger[J].InternationalJournalofHeatandMassTransfer， 2020， 150：119006.

[21]何雅玲， 王文奇， 邱羽， 等.熔鹽在復(fù)雜換熱結(jié)構(gòu)內(nèi)的對流換熱特性實驗研究及進展[J].科學(xué)通報 ， 2019， 64（28/29）：3007–3019.

[22] KATO Y， NITAWAKI T， MUTO Y. Medium temperaturecarbon dioxide gas turbine reactor[J]. Nuclear Engineering and Design， 2004， 230（1/3）：195–207.

[23]高毅超 ， 夏文凱 ， 龍穎 ， 等.管徑和轉(zhuǎn)折角對 Z 型 PCHE 換熱及壓降影響的研究[J].熱能動力工程， 2019，34（2）：94–100.

[24] WANGB，SHENJF，CHENGJD，etal. Numerical improvement using flow and heat transfer calculations of thezigzaggeometryforcarbondioxidePCHEs[J]. Energies， 2022， 15（8）：2831.

[25]顧維藻， 神家銳， 馬重芳， 等.強化傳熱[M].北京：科學(xué)出版社， 1990：10.

[26]史以奇， 李凌.螺紋管管內(nèi)流動與傳熱的數(shù)值模擬[J].上海理工大學(xué)學(xué)報， 2016， 38（2）：133–136，141.

[27] KIMIH，NOHC，LEEJI，etal. Thermalhydraulic performance analysis of the printed circuit heat exchanger usingaheliumtestfacilityandCFDsimulations[J]. NuclearEngineeringandDesign，2009，239（11）：2399–2408.

（編輯：石瑛）