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

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

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

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

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

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

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

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

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

1 換熱器模型

1.1 換熱器物理模型

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

圖1 PRHR HX物理模型圖

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

1.2 控制方程

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

守恒方程:

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

能量方程:

迪貝斯－貝爾特公式

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

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

2 邊界條件設(shè)置

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

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

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

3.1 溫度場(chǎng)分析

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

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

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

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

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

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

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

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

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

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

3.2 流場(chǎng)分析

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

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

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

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

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

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

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

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

4 結(jié)論

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

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

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

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

［2］薛若軍，鄧程程，彭敏俊，等.非能動(dòng)余熱排出換熱器數(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］林誠(chéng)格，郁祖盛.非能動(dòng)安全先進(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.

(編輯 楊 波)