楊 旭，周 濤，汝小龍，林達(dá)平，方曉璐

（華北電力大學(xué)核熱工安全與標(biāo)準(zhǔn)化研究所，北京102206）

在反應(yīng)堆運(yùn)行時，由于燃料棒、堆內(nèi)構(gòu)件等部件會受高壓過冷態(tài)的冷卻劑的腐蝕沖刷、磨蝕、FAC（管道流體加速腐蝕）效應(yīng)［1］等的影響，會產(chǎn)生許多不溶性腐蝕產(chǎn)物，腐蝕產(chǎn)物在堆芯的遷移會產(chǎn)生很多問題。腐蝕氧化物在堆芯燃料棒表面的沉積，當(dāng)沉積物達(dá)到一定厚度則會降低熱交換效率，進(jìn)而導(dǎo)致堆芯燃料過熱，嚴(yán)重時會導(dǎo)致堆芯局部熔化。同時，當(dāng)冷卻劑循環(huán)經(jīng)過堆芯時，所攜帶的腐蝕產(chǎn)物會被活化而具有放射性，具有放射性的腐蝕產(chǎn)物在一回路中的一些設(shè)備上沉積會在設(shè)備周圍形成輻射場，這是反應(yīng)堆運(yùn)行人員及維修人員吸收的放射性劑量的主要來源之一，并且會嚴(yán)重危害核反應(yīng)堆的平穩(wěn)運(yùn)行并增加運(yùn)行費(fèi)用。目前國內(nèi)外學(xué)者［2－6］的研究還是集中在壓水堆結(jié)構(gòu)材料（奧氏體不銹鋼、鎳基合金、鋯合金等）的耐蝕性和結(jié)構(gòu)材料的腐蝕開裂等方面。一些學(xué)者對細(xì)顆粒在溫度場內(nèi)的沉積規(guī)律進(jìn)行了研究，包括層流 和 湍 流 對 顆 粒 沉 積 的 影 響［7，8］。J.C.Barrett等［9］研究了放射性顆粒物在電場中的運(yùn)動規(guī)律。楊林民，周濤等［10］通過編制程序，采用熱泳沉積模型，計(jì)算嚴(yán)重事故工況下不同溫度、壓力、粒徑時安全殼內(nèi)氣溶膠的熱泳沉積效率。宣益民等［11］利用 Lattiee－Boltzmann方法建立納米流體的傳熱模型，分析含有細(xì)顆粒流體的流動與傳熱，并對這類流體和水的流動換熱參數(shù)進(jìn)行比較分析。姚軍和Michael Fairweather等［12］利用大渦模擬和拉格朗日方法研究了腐蝕產(chǎn)物（顆粒）在核反應(yīng)堆冷卻劑回路管道中的沉積現(xiàn)象。對于不溶性腐蝕產(chǎn)物在堆內(nèi)的沉積現(xiàn)象，尤其是燃料棒表面的沉積現(xiàn)象，國內(nèi)外學(xué)者研究較少，腐蝕產(chǎn)物在堆內(nèi)的沉積規(guī)律尚不清楚。通過數(shù)值模擬的方法，對不溶性粒狀腐蝕產(chǎn)物在堆芯燃料棒流域的沉積過程進(jìn)行數(shù)值模擬研究，對探究核電站反應(yīng)堆安全具有重要意義。

1 幾何模型

幾何模型如圖1所示，是依據(jù)大亞灣壓水堆核電站堆芯燃料棒區(qū)域?qū)嶋H尺寸建模，燃料包殼外徑為9.5mm，棒中心距為12.6mm，模型長度尺寸為1 000mm。圖2是模型截面網(wǎng)格示意圖。

圖1 幾何模型Fig.1 Geometric model

圖2 截面網(wǎng)格Fig.2 Section grid

采用正六面體網(wǎng)格對幾何模型進(jìn)行網(wǎng)格劃分，對燃料棒外表面邊界層區(qū)域網(wǎng)格進(jìn)行加密；進(jìn)行網(wǎng)格敏感性分析，從精度和效率上綜合考慮，確定模型網(wǎng)格節(jié)點(diǎn)數(shù)為5.0×105。

2 計(jì)算模型與方法

2.1 標(biāo)準(zhǔn)k－ε模型

在湍動能k方程的基礎(chǔ)上，引入一個湍動耗散率ε的方程，形成k－ε雙方程模型，稱為標(biāo)準(zhǔn)k－ε模型。該模型中的ε定義為：

湍流黏度μt可表示成k與ε的函數(shù)：

標(biāo)準(zhǔn)k－ε模型的輸運(yùn)方程為：

式中，μ為流體動力黏度，Pa·s；ρ為流體密度，kg／m3；ui為時均速度，m／s；Gk為層流速度梯度產(chǎn)生的湍流動能，kg·m2／s2；Gb為浮力產(chǎn)生的湍流動能，kg·m2／s2；YM為可壓縮湍流中脈動擴(kuò)張項(xiàng)，kg·m2／s2；Cμ、C1ε、C2ε和C3ε為經(jīng)驗(yàn)常量；αk和αε分別為k方程和ε方程的湍流普朗特?cái)?shù)；Sk和Sε為k方程和ε方程的用戶自定義源項(xiàng)，kg·m2／s2；xi、xj、xk是方向坐標(biāo)，x、y、z分別記作x1、x2、x3，其中i＝1、2、3；j＝1、2、3；k＝1、2、3。

2.2 離散相軌道計(jì)算模型

FLUENT是通過積分拉氏坐標(biāo)系下的顆粒作用力微分方程來求解離散相顆粒的軌道。顆粒的平衡方程在笛卡爾坐標(biāo)系下的形式（x方向）為：

式中，F(xiàn)D（u－up）為顆粒的單位質(zhì)量曳力，N／kg；u為流體相速度，m／s；up為顆粒速度，m／s；ρp為顆粒密度（骨架密度），kg／m3；gx為顆粒在x方向的重力加速度，m／s2；Fx為顆粒平衡方程中的其他作用力，N／kg。對于粒徑為1～10μm的顆粒，Stokes′曳力公式是適用的。這種情況下，F(xiàn)D定義為：FD＝18μ／dp2ρpCc，μ為流體動力黏度，Pa·s；dp為顆粒直徑，m；Cc為Stokes′曳力公式的Cunningham修正常數(shù)（考慮稀薄氣體力學(xué)的顆粒壁面速度滑移的修正），Cc＝1＋2λ／d（1.257＋0.4E－（1.1dp／2λ）），λ為氣體分子平均自由程，m。

2.3 計(jì)算方法

模擬使用的顆粒材料為碳酸鈣，粒徑為1μm。平均熱流密度為62W／m2，流場入口溫度為279.4℃，進(jìn)口設(shè)置為速度入口，流速為3.34m／s，出口設(shè)置為自由出口?？紤]顆粒在流場中所受的重力、熱泳力、布朗運(yùn)動、及沙夫曼力等。采用離散相模型，在歐拉坐標(biāo)系下考察流體相的運(yùn)動，在拉格朗日坐標(biāo)系下考察離散顆粒相的運(yùn)動。為了計(jì)算逃逸的顆粒，離散相進(jìn)出口設(shè)置為逃逸面（Escape）；離散相壁面條件設(shè)置為吸收（Trap）類型，即假定顆粒與壁面接觸后被捕獲，顆粒運(yùn)動終止。計(jì)算得到流場入口雷諾數(shù)超過2 320，從而判定流態(tài)為湍流。因此，對氣相采用標(biāo)準(zhǔn)模型來預(yù)測通道內(nèi)流場的湍流變化。

3 數(shù)值模擬結(jié)果與分析

3.1 截面溫度分布

圖3為流場入口溫度為279.4℃時，流場充分發(fā)展段堆芯流域溫度分布。

圖3 溫度分布Fig.3 Temperature distribution

由圖3可以看出，堆芯流域流場充分發(fā)展段溫度分布，整體上呈中心對稱分布，燃料棒外表面溫度最高，中心溫度最低，這是由于燃料棒對流體的加熱作用引起的。此外，相比于圓形管道中流場溫度分布，堆芯流域結(jié)構(gòu)對流場溫度分布有擠壓作用，使溫度分布有畸變。

3.2 流場分布

圖4和圖5分別為流場入口溫度為279.4℃時，流場充分發(fā)展段堆芯流域軸向速度分布和徑向速度分布。

圖4 軸向速度分布Fig.4 Axial average velocity distribution

由圖4可以看出，堆芯流域軸向速度整體上呈中心對稱分布，和圖3溫度分布類似，軸向速度在堆芯流域中心區(qū)域速度最大，在靠近燃料棒區(qū)域軸向速度快速下降，形成這種分布的原因是堆芯流域結(jié)構(gòu)對流場擠壓造成的。

圖5 徑向速度分布Fig.5 Radial average velocity distribution

從圖5可以看出，堆芯流域徑向速度分布也為中心對稱分布，但徑向速度分布在靠近燃料棒表面附近形成了四個對稱分布的峰值區(qū)域，這是由于堆芯流域中心區(qū)域湍流較強(qiáng)，導(dǎo)致流體橫向流動較強(qiáng)引起的；而在靠近對稱面的區(qū)域則形成了四個對稱分布谷值，這是由于此區(qū)域流場速度較小，流場分布已經(jīng)變成層流，流體橫向流動較弱引起的。

3.3 湍動能分布

圖6為流場入口溫度為279.4 ℃時，流場充分發(fā)展段堆芯流域湍動能分布。

圖6 湍動能分布Fig.6 Turbulent kinetic energy distribution

由圖6可以看出，與前面流場分布類似，流場湍動能分布呈中心對稱分布，在堆芯流域中心和對稱面中心形成谷值，這是由于這些區(qū)域流體的湍流速度漲落較小，導(dǎo)致流場中心區(qū)域出現(xiàn)湍動能谷值；而在燃料棒外表面近壁面區(qū)域形成峰值，這是由于該區(qū)域流場徑向速度出現(xiàn)峰值，湍流速度漲落較大，導(dǎo)致湍流增強(qiáng)引起的。

3.4 顆粒濃度分布

圖7～圖9分別為流場入口溫度為279.4℃時，進(jìn)出口處不溶性顆粒濃度分布和燃料棒外表面顆粒沉積分布。

圖7 進(jìn)口處顆粒濃度分布Fig.7 Mass concentration distribution at entrance

圖8 出口處顆粒濃度分布Fig.8 Mass concentration distribution at exit

由圖7和圖8可以看出，不溶性腐蝕產(chǎn)物顆粒在堆芯進(jìn)出口處靠近燃料棒表面區(qū)域形成了高濃度區(qū)域，尤其在對稱面附近顆粒濃度更高，而在堆芯中心區(qū)域顆粒濃度很低，這主要是由于在湍流較弱和溫度梯度較大區(qū)域，顆粒的熱泳效應(yīng)致使顆粒在此區(qū)域形成高濃度區(qū)域。但從圖7和圖8中也明顯可以看出，進(jìn)口處顆粒濃度明顯高于出口處，這主要是由于顆粒的沉降引起的。從圖8中還可看出，燃料棒表面高濃度區(qū)域，在靠近燃料棒表面附近，出現(xiàn)了一個低濃度區(qū)域，這主要是由于在燃料棒近壁面區(qū)域，溫度梯度較大，流場湍動能出現(xiàn)峰值，不利于顆粒在此區(qū)域形成穩(wěn)定的高濃度區(qū)域。

圖9 燃料包殼表面顆粒沉積分布Fig.9 Mass concentration distribution at fuel cladding

由圖9可以看出，顆粒在燃料包殼出口段表面形成點(diǎn)狀沉積，在包殼表面點(diǎn)狀沉積區(qū)域會出現(xiàn)包殼點(diǎn)蝕現(xiàn)象，點(diǎn)蝕區(qū)域也會引起局部傳熱惡化，破壞包殼完整性；顆粒在燃料包殼進(jìn)口段表面形成大面積附著沉積，這會改變包殼材料的熱導(dǎo)率和堆芯中子通量分布，引起堆芯軸向功率偏移。針對這種現(xiàn)象，應(yīng)該加強(qiáng)一回路水化學(xué)控制，停堆換料期間對顆粒沉積量大的局部區(qū)域進(jìn)行定點(diǎn)清除，并且應(yīng)當(dāng)建立堆內(nèi)構(gòu)件損傷追蹤計(jì)劃，及時更換損傷嚴(yán)重的堆內(nèi)構(gòu)件。

4 結(jié)論

（1）堆芯流域溫度分布、流場分布和湍動能分布都成近似中心對稱分布。

（2）由于熱泳效應(yīng)和湍流效應(yīng)的綜合影響，顆粒在對稱面附近形成高濃度區(qū)域；由于顆粒的沉降作用，入口段顆粒濃度比出口段高。

（3）堆芯流域顆粒物在包殼入口段表面呈大面積附著沉積，這會改變包殼材料的熱導(dǎo)率和堆芯中子通量分布，引起堆芯軸向功率偏移；而在包殼出口段表面呈點(diǎn)狀沉積，這可能導(dǎo)致包殼出現(xiàn)點(diǎn)蝕現(xiàn)象，點(diǎn)蝕區(qū)域也可能會引起傳熱惡化，破壞包殼完整性。

（4）堆內(nèi)腐蝕緩解措施，首先應(yīng)該加強(qiáng)一回路水化學(xué)控制，停堆換料期間對顆粒沉積量大的局部區(qū)域進(jìn)行定點(diǎn)清除，其次應(yīng)當(dāng)建立堆內(nèi)構(gòu)件損傷追蹤計(jì)劃，及時更換損傷嚴(yán)重的堆內(nèi)構(gòu)件。

［1］ 束國剛，薛飛，遆文新，等 .核電廠管道的流體加速腐蝕及其老化管理［J］.腐蝕與防護(hù) .2006，27（2）：72－76.

［2］ Rokuro Nishimura，Yasuaki Maeda.SCC evaluation of type 304and 316austenitic stainless steels in acidic chloride solutions using the slow strain rate technique［J］.Corrosion Science.2004，46：769－785.

［3］ R.Nishimura，A.Sulaiman，Y.Maeda.Stress corrosion cracking of austenitic stainless steel Type 316in acid solutions and intergranular SCC mechanism：effects of anion species （Cl－and SO42－and sensitizing temperature）［J］.Environment－Induced Cracking of Materials.2008，1：419－433.

［4］ Pradeep Kumar Samantaroy，Girija Suresh，Panita Paul，etc.Corrosion behavior of Alloy 690and Alloy 693im simulated nuclear high level waste medium［C］.Journal of Nuclear Materials.2011，418：27－37.

［5］ 邢建東，高義民，張國賞 .不銹鋼與碳鋼的液固兩相流沖刷腐蝕磨損研究［J］.西安交通大學(xué)學(xué)報(bào)，2004，38（5）：469－473.

［6］ 吳成紅，甘復(fù)興 .金屬在兩相流動水體中的沖刷腐蝕［J］.材料保護(hù)，2000，33（4）：33－35.

［7］ WALKER K L，HOMSY G M，GEYLING R T.Thermophoretic deposition of small particles in laminar tube flow［J］.Journal of Colloid and Interface Science.1979，69（1）：138－147.

［8］ NISHIO G， KITANI S， TAKAHASHI K.Thermophoretic deposition of aerosol particles in a heatexchanger pipe［J］.Industrial ＆Engineering Chemistry Process Design and Development.1974，13 （4）：408－415.

［9］ J.C.Barrett，C.F.Clement，A.B.S.Virdee The removal of radioactive aerosols by electric fields［J］.Journal of Aerosol Science.2009，40（3）：185－192.

［10］ 楊林民，周濤，陸道綱 .壓水堆嚴(yán)重事故下氣溶膠熱映沉積規(guī)律［J］.原子能科學(xué)技術(shù) .2008，42（1）：63－66.

［11］ 宣益民，余凱，等 .基于Lattiee－Boltzmann方法的納米流體流動與傳熱分析［J］.工程熱物理學(xué)報(bào).2004，25：1022－1024.

［12］ 姚軍，Michael Fairweather，李寧 .粒狀腐蝕產(chǎn)物在管道中的沉積機(jī)理［J］.金屬學(xué)報(bào)，2011，47（7）：804－808.