AP1000非能動換熱器數(shù)值模擬

薛若軍，王明遠(yuǎn)，李朝君，祝 賀

（1.哈爾濱工程大學(xué)核安全與仿真技術(shù)國防重點(diǎn)學(xué)科實(shí)驗(yàn)室，150001 哈爾濱;2.海軍裝備研究院，100161 北京）

AP1000非能動換熱器數(shù)值模擬

薛若軍1，2，王明遠(yuǎn)1，李朝君1，祝 賀1

（1.哈爾濱工程大學(xué)核安全與仿真技術(shù)國防重點(diǎn)學(xué)科實(shí)驗(yàn)室，150001 哈爾濱;2.海軍裝備研究院，100161 北京）

為準(zhǔn)確認(rèn)識AP1000非能動余熱排出熱交換器的具體換熱和流動過程，針對AP1000非能動換熱器，應(yīng)用FLUENT軟件對其非穩(wěn)態(tài)傳熱特性和流動進(jìn)行研究，比較分析同一時(shí)刻不同位置溫度場和流場的分布，以及不同時(shí)刻同一位置溫度場和流場的變化，分析換熱器傳熱和流動特性隨時(shí)間的變化，同時(shí)對比分析直管形管束和彎形管束對換熱器自然對流的影響.分析結(jié)果表明:溫差、漩渦和回流使溫度場和流場相互影響，彎形管束使流體自然對流更加復(fù)雜.研究成果有助于分析換熱器的自然循環(huán)能力，為我國新建的AP1000先進(jìn)壓水堆核電機(jī)組的安全運(yùn)行提供參考.

非能動換熱器;自然循環(huán);彎管;傳熱;流動特性

非能動余熱排出換熱器(PRHR HX)是AP1000非能動余熱排出系統(tǒng)中的主要設(shè)備，在事故或停堆時(shí)利用自然循環(huán)將反應(yīng)堆內(nèi)的余熱排出［1］.本文針對AP1000彎管結(jié)構(gòu)的換熱器進(jìn)行數(shù)值模擬，分析溫度場和流場變化，并同文獻(xiàn)［2］的模擬結(jié)果進(jìn)行比較.文獻(xiàn)［2］主要是將彎形管簡化為直管后進(jìn)行模擬分析.雖然這種簡化在一定程度上可以看出換熱器的溫度場和流場變化，但是仍無法反應(yīng)真實(shí)的變化過程.

用FLUENT軟件對管殼式換熱器的模擬已經(jīng)有過許多研究，如O.Ender等［3］研究了折流板的寬度對管殼式換熱器殼側(cè)傳熱特性的影響，并且從模擬結(jié)果的對比中得出了適合管殼式換熱器數(shù)值模擬的最佳湍流模型;A.Karno和 S.Ajib等［4］研究了不同傳熱管的排列方式對換熱效果的影響;C.Uday等［5］開發(fā)了一種適用于研究管殼式換熱器殼側(cè)壓降的理論模型，并且研究了管殼式換熱器殼側(cè)壓力的變化.這些研究都充分說明了數(shù)值模擬研究管殼式換熱器的可行性和有效性.

1 換熱器模型

1.1 換熱器物理模型

圖1是AP1000非能動余熱排出換熱器(PRHR HX)的局部物理模型［6－8］，換熱器安裝在安全殼內(nèi)換料水貯存箱(IRWST)里，IRWST的高度為16 m，半徑為7 m，內(nèi)側(cè)傳熱管距IRWST壁面約為2.49 m，傳熱管束粗約為2.13 m.由于沒有研究管內(nèi)的流動情況，這樣只需對傳熱管外部進(jìn)行網(wǎng)格畫分，便可以保證網(wǎng)格數(shù)量在計(jì)算機(jī)的承受能力范圍內(nèi).在盡量不改變原模型具體結(jié)構(gòu)的基礎(chǔ)上，使模擬的結(jié)果更接近于原模型，并可同時(shí)研究彎管區(qū)的流動情況.換熱器的物理模型和管束排列方式同文獻(xiàn)［2］相同.換熱器簡化模型及坐標(biāo)方向圖如圖2所示.

圖1 PRHR HX物理模型圖

圖2 簡化模型及坐標(biāo)方向

1.2 控制方程

控制方程如下［2，9－10］.

守恒方程:

湍流模型采用k－ε湍流模型，對連續(xù)方程，φ =1;動量方程，φ =u，v，w;湍流脈動方程，φ =k;湍流耗散率方程，φ=ε.每個(gè)控制體共有6個(gè)方程.

能量方程:

迪貝斯－貝爾特公式

對于自然對流問題，密度的處理較為關(guān)鍵.由于水箱內(nèi)流體密度變換不太大，為得到更好的收斂速度，選用Boussinesq近似，除動量方程的浮力項(xiàng)外，其他方程中的密度均看成常數(shù).動量方程中的浮力項(xiàng)可表示為

其中:ρ0為流體的(常)密度;T0為操作溫度;β為熱膨脹系數(shù)，取β≈0.000 52，以上各個(gè)方程中各項(xiàng)物理意義及經(jīng)驗(yàn)系數(shù)詳情參見文獻(xiàn)［7］.

2 邊界條件設(shè)置

因?yàn)闆]有畫管內(nèi)網(wǎng)格，所以不能設(shè)置管入口和出口邊界條件，只能設(shè)置管壁邊界條件.本文采用定壁溫的邊界條件，并將管壁溫度線性化處理來反映沿管長方向溫度的變化.圖3是將直管型換熱器壁溫分不同段數(shù)處理時(shí)溫度云圖.從云圖分析可知，分段越少，溫度分層越明顯，這是因?yàn)榉謱由俑邷囟屋^長使傳熱更加明顯;但是分段越多更能反映實(shí)際的傳熱過程.本文結(jié)合彎管結(jié)構(gòu)，將管壁溫線性化分為5段.這里要注意管內(nèi)流體的溫度并不是管外壁的溫度，必須根據(jù)相關(guān)的傳熱知識計(jì)算出管外壁溫.首先由方程(1)計(jì)算出管內(nèi)側(cè)對流換熱系數(shù)，再由對流換熱公式計(jì)算出壁溫.分別計(jì)算出管進(jìn)出口外壁溫后，再將壁溫線性化分成5份處理.水箱外壁設(shè)為絕熱壁面.水箱初始水溫48.9℃，傳熱管的入口溫度為297.2℃，出口溫度為92.8℃.對于非穩(wěn)態(tài)問題的求解，要選取合適的時(shí)間步長，還要根據(jù)殘差監(jiān)視曲線的變化趨勢調(diào)整時(shí)間步長和最大迭代步數(shù).計(jì)算時(shí)先取較小的時(shí)間步長(如0.1 s)，計(jì)算到一定的時(shí)間后再逐漸調(diào)大時(shí)間步長［11－12］.

圖3 壁溫分不同段線性化處理溫度分布云圖

3 計(jì)算結(jié)果分析

3.1 溫度場分析

3.1.1 300 s時(shí)刻的溫度場分析

在300 s時(shí)刻，換熱器內(nèi)的溫度場和速度場開始出現(xiàn)明顯的變化，所以本文取300 s為研究時(shí)間進(jìn)行分析.云圖中溫度的單位為K，速度單位m/s.

從圖4(a)可以看出，越靠近管壁的直線(x=－2.5)上各點(diǎn)的溫度越高，但在水箱上部(z＞12.8)卻出現(xiàn)x=2.5 m曲線高于x=0 m曲線，x=0 m曲線高于x= －2.5 m曲線的情況，這說明水箱中上部區(qū)域出現(xiàn)了局部的回流現(xiàn)象.從每條曲線的總體趨勢來看，都是沿z方向溫度逐漸升高，這主要是受管壁線性化溫度高低不同的影響，其次也由于熱流體在浮力驅(qū)動下向上流動導(dǎo)致的.對比3條曲線，可以發(fā)現(xiàn)在同一高度z不同水平位置x的溫度變化不同.在較低位置，水箱不同水平位置的溫度差異較小;在較高位置，水箱不同水平位置的溫度差異較大，這是由于在較高位置換熱和流動更強(qiáng)烈導(dǎo)致的.這一模擬結(jié)果和參考文獻(xiàn)［2］中相似.在圖4(b)中從3條曲線的整體分布狀態(tài)來看，直管段高度越高的水平面上溫度越高，而遠(yuǎn)離彎管區(qū)域的z=15 m曲線的溫度則低于靠近彎管區(qū)的z=12.8 m曲線，這是因?yàn)樗膶?dǎo)熱系數(shù)較小，傳遞到z=15 m處的熱量相對較小.在圖5(a)中可以看出水箱中心明顯的回流現(xiàn)象，而遠(yuǎn)離中心管束的區(qū)域這一現(xiàn)象還不明顯，這說明越靠近中心管束集中區(qū)域的流體的傳熱和對流效果越強(qiáng)，而在圖5(c)中幾乎看不出管束的影響，等溫線變化也不明顯.

圖4 沿x方向和z方向的溫度分布

圖5 3個(gè)不同縱截面的溫度分布云圖

3.1.2 不同時(shí)刻的溫度場分析

圖6是不同時(shí)刻截面y=0.109 m(第一排管)的溫度云圖.從云圖的變化來看，開始時(shí)導(dǎo)熱占主導(dǎo)作用，管束通過熱傳導(dǎo)作用把熱量傳遞給周圍流體，使周圍流體的溫度迅速升高;但由于管壁溫度沿z正方向線性升高，所以換熱效果也沿z正方向逐漸升高;管束附近的熱流體在浮力驅(qū)動下向上流動，水箱內(nèi)水溫從下至上呈逐漸升高的趨勢，且等值線向右側(cè)發(fā)生傾斜.由于水的導(dǎo)熱系數(shù)較小，因而開始時(shí)溫度明顯升高的區(qū)域只是管束周圍的小范圍區(qū)域;但隨著時(shí)間的增長，流體之間溫差的增大使得流體產(chǎn)生自然對流現(xiàn)象，流動作用增強(qiáng)，溫度等值線發(fā)生彎曲，使管束附近集中的熱量不斷擴(kuò)散出去.從t=150 s的圖中可以看出，在水箱上部溫度等值線彎曲最大，出現(xiàn)兩邊低中間高的趨勢，這說明在水箱上部發(fā)生了強(qiáng)烈的回流現(xiàn)象，促進(jìn)了流體之間以及流體和壁面之間的換熱.在t=300 s時(shí)，隨著自然對流的作用，水箱上部等溫線彎曲減弱，而水箱中下部等溫線也發(fā)生了不同程度的彎曲.同文獻(xiàn)［2］中的直管模擬結(jié)果比較發(fā)現(xiàn)，相同時(shí)刻彎管結(jié)構(gòu)換熱器水箱中云圖變化更為明顯，說明換熱效果更好，而且相同時(shí)刻彎管換熱器的繞流現(xiàn)象也更加明顯，這對于及時(shí)排出余熱是有利的.

圖7是不同時(shí)刻z=12.8 m截面的溫度云圖.

圖6 不同時(shí)刻y=0.109 m(第1排管)截面的溫度云圖

圖7 不同時(shí)刻z=12.8 m截面的溫度云圖

從溫度云圖的變化來看，開始時(shí)導(dǎo)熱作用產(chǎn)生的管束附近的溫度梯度很明顯.隨著時(shí)間的增長，流體在熱浮力驅(qū)動下產(chǎn)生自然對流，熱流體有貼水箱側(cè)壁流動的趨勢，使得水箱壁面附近溫度較高，中間部分區(qū)域溫度反而較低，出現(xiàn)繞流現(xiàn)象，且隨時(shí)間的變化繞流區(qū)的位置也不斷變化.從圖中可以看出，剛開始只在彎管區(qū)有小范圍的繞流，在t=300 s，彎管區(qū)的左上方和右上方分別存在繞流，在t=600 s時(shí)，左上方繞流消失，而右上方繞流開始像左移動，有逐漸減小的趨勢，在t=1 000 s時(shí)，繞流繼續(xù)減小，并轉(zhuǎn)化為幾個(gè)局部小的回流，在t=2 000 s時(shí)，水箱中繞流現(xiàn)象基本消失，水溫分布趨于均勻.

3.2 流場分析

3.2.1 300 s時(shí)刻的速度場分析

圖8(a)是中心面上的3條平行的豎直線沿z方向的速度分布曲線.從整體來看，流速隨高度增大不斷升高，越靠近管壁的區(qū)域這種上升趨勢越明顯(如x= －2.5 m曲線)，而較遠(yuǎn)位置處這種上升趨勢較平緩.這是因?yàn)榭拷鼙趨^(qū)的流體受管壁的影響越大傳熱效果越好，產(chǎn)生的熱浮力驅(qū)動作用最大;而遠(yuǎn)離管壁處溫度較低，產(chǎn)生的浮升力較小.在中上部位置，3條曲線由于在彎管上部區(qū)域，溫差較大，自然對流比較劇烈，流速相對彎管下部區(qū)域大很多，而隨后在此區(qū)域由于溫差和自然對流的作用出現(xiàn)了繞流現(xiàn)象，導(dǎo)致流速出現(xiàn)先增大后減小又增大的過程.從x= －2.5 m曲線可以看出，由于彎管區(qū)結(jié)構(gòu)發(fā)生突變使流動的復(fù)雜性增強(qiáng)，導(dǎo)致在上部彎管區(qū)出現(xiàn)了速度先降低后升高，隨后結(jié)構(gòu)變得開闊，速度又逐漸平穩(wěn)升高.圖8(b)為中心面上3條平行的水平線沿x方向的速度分布曲線.z=12.8 m曲線的形狀反映了管束附近區(qū)域的速度梯度分布.而z=1.24 m和z=12.8 m兩條曲線的彎管區(qū)域溫度的波動較同一位置直管區(qū)更劇烈.在z=15 m曲線中，在彎管區(qū)左側(cè)流速出現(xiàn)了先降低又升高的趨勢，這是因?yàn)樵趶澒艿淖笊戏匠霈F(xiàn)了局部的回流.而從z=15 m和z=12.8 m兩條曲線中可以看出，在彎管區(qū)右側(cè)，由于出現(xiàn)了繞流，流速先增大后減小，且繞流已經(jīng)開始向左移動，越靠近彎管壁面，換熱越明顯，向左移動得越快.

圖9是IRWST內(nèi)部分區(qū)域的流線圖.從圖中可以看到:有熱流體沿管壁上升而后向下縱向繞流的較規(guī)則的運(yùn)動趨勢;產(chǎn)生了漩渦，漩渦的攪動有利于對流換熱;局部的繞流和回流現(xiàn)象不僅僅是縱向和橫向作用的，而是相互影響，形成了斜向的漩渦.這說明水箱內(nèi)各流體的運(yùn)動規(guī)律是相互影響的復(fù)雜運(yùn)動.從圖中還可以清晰地看到彎管左上方、右上方、彎管內(nèi)側(cè)以及彎管外側(cè)局部的回流，同時(shí)還可以看到流體沿著水箱筒體壁面的繞流.

圖8 沿z和x方向的速度分布

圖9 IRWST內(nèi)部分區(qū)域流線圖

3.2.2 不同時(shí)刻的速度場分布

由圖10可知，在換熱器上部彎管區(qū)換熱很強(qiáng)烈，150 s時(shí)該部分流體流動的速度已經(jīng)很快，而且出現(xiàn)了漩渦的趨勢，隨著換熱的繼續(xù)進(jìn)行，出現(xiàn)了明顯的繞流現(xiàn)象，產(chǎn)生了更加明顯的漩渦.這種流動形式破壞了管束的邊界層使換熱效果更強(qiáng).同直管式換熱器相比，彎管換熱器的換熱效果更好，繞流現(xiàn)象更加明顯，在水箱的上部出現(xiàn)了大量的漩渦，加劇了換熱，而且在兩彎管之間的區(qū)域換熱也明顯增強(qiáng).

圖10 不同時(shí)刻y=0.109 m(第1排管)截面的速度云圖

4 結(jié)論

1)由于溫差、漩渦和繞流的存在，換料水箱中溫度場和速度場產(chǎn)生相互影響，在彎管處溫度場和速度場的變化較劇烈，有利于換熱.

2)隨著時(shí)間的增長，熱流體沿著水箱壁面和彎管束內(nèi)側(cè)進(jìn)行繞流，并且強(qiáng)度、位置和方向不斷的變化.由于受到了彎管束結(jié)構(gòu)的影響，使得流體在換料水箱中的自然對流作用更加復(fù)雜，而流體的流速和方向的多變會更有利于熱量的傳遞，增加了該換料水箱的非能動性.

3)管束附近產(chǎn)生了較強(qiáng)的湍流，在彎管區(qū)繞流和回流現(xiàn)象更加明顯.由于彎管結(jié)構(gòu)增加了流體的擾動性，所以同直管換熱器相比彎管換熱器的換熱效果更好，更能及時(shí)的排出反應(yīng)堆內(nèi)的余熱，保證反應(yīng)堆的安全.

［1］西屋電氣公司.西屋公司的AP1000先進(jìn)非能動型核電廠［J］.現(xiàn)代電力，2006，25(3):60－61.

［2］薛若軍，鄧程程，彭敏俊，等.非能動余熱排出換熱器數(shù)值模擬［J］.原子能科學(xué)技術(shù)，2010，44(4):429 －435.

［3］ENDER O，LLKER T.Shell side CFD analysis of a small shell-and-tubes heat exchanger［J］.Energy Conversion and Management，2010，51:1004 －1014.

［4］KARNO A，AJIB S.Effect of tubes pitch on heat transfer in shell and tubes heat exchangers-new simulation software［J］.Heat Mass Transfer，2006，29:263 －70.

［5］UDAY C.Modeling for shell-side pressure drop for liquid flow in shell-and-tubes heat exchanger［J］.International Journal of Heat and Mass Transfer，2006，49:601－610.

［6］林誠格，郁祖盛.非能動安全先進(jìn)核電廠 AP1000［M］.北京:原子能出版社，2008:138－142，157－161.

［7］RRYES J，WOODS B.Testing of passive safety system performance for higher power advanced reactors［R］.［S.l.］:Oregon State University，2004.

［8］SCHULZ T L.Westinghouse AP1000 advanced passive plant［J］.Nuclear Engineering and Design，2006，236(5):1547－1557.

［9］陶文銓.數(shù)值傳熱學(xué)［M］.第2版.西安:西安交通大學(xué)出版社，2001:333－370.

［10］楊世銘，陶文銓.傳熱學(xué)［M］.北京:高等教育出版社，2006:237－243.

［11］ANSYS Fluent Version 6.3.Fluent 6.3 user's guide［M］.Lebanon:ANSYS Inc，2006.

［12］DONG Q W，WANG Y Q，LIU M S.Numerical and experimental investigation of shell side characteristics for ROD baffle heat exchanger［J］.Applied Thermal Engineering，2008，28:651－660.

Numerical simulation of AP1000 passive residual heat removal heat exchanger

XUE RUO-jun1，2，WANG Ming-yuan1，LI Chao-jun1，ZHU He1

(1.National Defense Key Subject Laboratory for Nuclear Safety and Simulation Technology，Harbin Engineering University，150001 Harbin，China;2.Navy Equipment Research Institute，100161 Beijing，China)

To accurately know heat transfer and flow processes of AP1000 passive residual heat removal heat exchanger(PRHR HX)，and provide a reference for the safe operation of our new AP1000 advanced pressurized water reactor nuclear power plant，F(xiàn)LUENT software is used.By comparative analysis of the distribution of temperature-field and flow-field in different locations at the same time，and the change of temperature-field and flow-field in the same location at different times，we can know the heat transfer and flow characteristics which are changed with time more intuitive，and analyze the effect of straight tubes and elbow tubes on the natural circulation.The results show that the temperature difference，vortex，flow around make the interaction of temperature field and flow field.Due to the influence of elbow tubes structure，the fluid natural convection in the IRWST is more complex.

passive residual heat removal heat exchanger;natural circulation;elbow tubes;heat-transferring;flow characteristic

TL33

A

0367－6234(2012)11－0133－06

2012－03－23.

高等學(xué)校創(chuàng)新引智計(jì)劃資助項(xiàng)目(B08047).

薛若軍(1974—)，男，副教授.

薛若軍，xueruojun@hrbeu.edu.cn.

(編輯 楊 波)