空間核熱推進(jìn)粒子球床堆堆芯設(shè)計(jì)參數(shù)敏感性分析

姜奪玉，江新標(biāo)，王立鵬，潘孝兵，張信一

（西北核技術(shù)研究所 強(qiáng)脈沖輻射環(huán)境模擬與效應(yīng)國(guó)家重點(diǎn)實(shí)驗(yàn)室，陜西 西安 710024）

核能推進(jìn)是繼化學(xué)能推進(jìn)后人類目前所掌握的較先進(jìn)的推進(jìn)方式，太空核能推進(jìn)的開(kāi)發(fā)使人類進(jìn)行深空探索的愿望得以實(shí)現(xiàn)。空間核熱推進(jìn)（SNTP）是核能推進(jìn)的重要分支。20世紀(jì)80至90年代初，美國(guó)相繼啟動(dòng)了星球大戰(zhàn)計(jì)劃和空間探索計(jì)劃，這期間隨著材料技術(shù)的發(fā)展，固相核熱火箭推進(jìn)技術(shù)重新受到重視，并提出了一系列新的改進(jìn)型堆芯方案［1-2］，粒子球床堆（PBR）便是其代表之一。

控制鼓外徑與反射層厚度、10B 富集度及厚度、柵徑比與高徑比、燃料球體積填充率等均關(guān)系到PBR 堆芯結(jié)構(gòu)尺寸的選擇，研究它們的變化對(duì)堆芯反應(yīng)性的影響規(guī)律具有重要意義。研究結(jié)果對(duì)PBR堆芯設(shè)計(jì)參數(shù)的優(yōu)化具有參考價(jià)值。MCNP程序［3］在處理復(fù)雜幾何結(jié)構(gòu)方面有較強(qiáng)優(yōu)勢(shì)，適合PBR 小型空間堆的模擬計(jì)算。本文基于空間核熱推進(jìn)粒子球床堆物理模型，開(kāi)展堆芯設(shè)計(jì)參數(shù)的敏感性分析研究。

1 PBR堆芯模型

PBR 燃料元件模型如圖1 所示。堆芯橫向和縱向截面如圖2 所示，堆芯外圍為含有12個(gè)控制鼓的氧化鈹反射層。PBR 堆芯物理參數(shù)列于表1。燃料球基本柵元采用點(diǎn)對(duì)點(diǎn)六面柱體柵元（CHPOP）模型，實(shí)際體積填充率為56.97%。

圖1 PBR 燃料元件模型Fig.1 Model of PBR fuel element

圖2 PBR 堆芯的橫截面（a）和縱截面（b）示意圖Fig.2 Schemes of transverse section（a）and vertical section（b）for PBR

表1 PBR 堆芯主要參數(shù)［4］Table 1 Main parameters of PBR core［4］

2 敏感性分析

2.1 控制鼓外徑與反射層厚度對(duì)PBR 堆芯物理特性的影響

在PBR 堆芯中，減小反射層厚度可增加徑向中子泄漏，進(jìn)而提高控制鼓控制價(jià)值，但反射層需能支撐控制鼓，故反射層的厚度取決于控制鼓的尺寸。圖3為keff及中子泄漏率隨控制鼓外徑和反射層厚度的變化。

由圖3a可知，控制鼓全開(kāi)至全閉過(guò)程中，隨控制鼓外徑的增大，keff減小、Δkeff增大，這主要是因?yàn)榭刂乒耐鈴皆龃髸r(shí)，全開(kāi)至全閉過(guò)程中吸收體之間的相對(duì)距離增加，控制鼓價(jià)值增大，且Δkeff的增加與中子泄漏率呈正相關(guān)。由圖3b可知，keff隨反射層厚度的增大而增大，且開(kāi)始時(shí)增幅較快，這主要是因?yàn)榉瓷鋵雍穸仍龃笫怪凶有孤p少，且反射層節(jié)省是反射層厚度的正割函數(shù)，使得keff的增幅先快后趨于平緩。

圖3 keff及中子泄漏率隨控制鼓外徑和反射層厚度的變化Fig.3 keffand neutron leakage rate vs drum outer diameter and reflector thickness

2.2 10B富集度及厚度對(duì)PBR 堆芯物理特性的影響

控制鼓吸收體材料10B的變化對(duì)控制鼓的可控性能有一定影響，本文研究10B 富集度和10B厚度對(duì)控制鼓控制價(jià)值的影響，在控制鼓全開(kāi)、全閉情況下對(duì)Δkeff進(jìn)行計(jì)算，結(jié)果示于圖4。

由 圖4可 知：1）改 變10B 富 集 度 與10B 厚度，得到的最大控制價(jià)值相近，Δkeff均近似為0.055；2）隨著10B厚度的增大，Δkeff存在1個(gè)峰值（圖4b），這是因?yàn)殡S著10B 厚度的增大，中子吸收增加，控制鼓價(jià)值增大，與此同時(shí)，控制鼓全開(kāi)到全閉過(guò)程中，吸收體相對(duì)距離減小，故10B 厚度達(dá)到某一臨界值時(shí)，控制價(jià)值 反 而 下 降；3）10B 富 集 度 和10B 厚 度 變 化 對(duì)Δkeff影響較小，相應(yīng)Δkeff最多變化約0.009 87（10B富集度變化范圍19.8%～100%，10B厚度變化范圍0.2～5cm），相反，控制鼓外徑變化對(duì)Δkeff的影響較明顯，相應(yīng)Δkeff最多變化約0.037（圖3a），約為前者的3.7倍。

2.3 組件柵徑比與堆芯高徑比對(duì)PBR 堆芯物理特性的影響

1）組件柵徑比對(duì)PBR堆芯物理特性的影響

保持燃料區(qū)體積不變，取柵徑比x＝1.2、1.3、1.4、1.5、1.75、2.0，得到keff示于圖5，圖6為徑向中子注量率不均勻因子隨徑向距離的變化（x 的變化主要影響徑向中子分布）。

由圖5可知，隨x 增大，keff先增大后減小，且有1個(gè)明顯峰值。這主要因?yàn)椋憾研救剂蠀^(qū)體積不變，當(dāng)x 較小時(shí)，相應(yīng)慢化劑量較少，導(dǎo)致中子慢化不充分，堆芯處于欠慢化區(qū)；相反，x 較大時(shí)，中子過(guò)度慢化，堆芯處于過(guò)慢化區(qū)。另一方面，增加慢化劑將導(dǎo)致堆芯中子能譜軟化。理論上，小的x 使得堆芯結(jié)構(gòu)尺寸較小，更符合小堆芯的理念，但要兼顧堆芯臨界條件。由圖6可知，隨x 增大，慢化劑增多，組件柵距增大，中子注量率不均勻因子變化范圍增大，不均勻性增大。

圖4 Δkeff隨10B富集度和10B厚度的變化Fig.4 Δkeffvs 10B enrichment and 10B thickness

2）堆芯高徑比對(duì)PBR堆芯物理特性的影響保持柵距、燃料區(qū)體積不變，取高徑比z＝1.016、1.1、1.2、1.3，得到控制鼓全開(kāi)與全閉時(shí)的有效增殖因數(shù)差值Δkeff示于圖7。由圖7可知，隨z增大，堆芯整個(gè)體積增大，慢化區(qū)增大，而燃料區(qū)體積不變，使中子泄漏減少，Δkeff減小。需要說(shuō)明的是，對(duì)于給定的堆芯，慢化劑厚度是一定的，這里未考慮慢化劑的實(shí)際厚度，實(shí)際上z 增大可增加徑向中子泄漏，對(duì)控制鼓的可控性是有益的。圖8為不同堆芯高徑比下的軸向中子注量率不均勻因子分布（堆芯高徑比變化主要影響軸向中子分布）。由圖8可知，隨z增大，軸向中子注量率的不均勻性變化范圍增大，堆芯不均勻性增大。

圖5 keff隨x 的變化Fig.5 keffvs x

圖6 徑向中子注量率不均勻因子隨徑向距離的變化Fig.6 Radial neutron fluence rate inhomogeneous factor vs radial distance

2.4 燃料球體積填充率對(duì)PBR 堆芯物理特性的影響

為簡(jiǎn)化計(jì)算過(guò)程，選取簡(jiǎn)單立方體柵元（理論體積填充率52.36%，用RPP 表示）和點(diǎn)對(duì)點(diǎn)六面柱體柵元（理論體積填充率60.46%，用HEX 表示）兩種模型［5-6］。

圖7 Δkeff隨z的變化Fig.7 Δkeffvs z

圖8 不同堆芯高徑比下軸向中子注量率不均勻因子分布Fig.8 Axial neutron fluence rate inhomogeneous factor at different z

1）體積填充率對(duì)keff的影響

研究［7］表明，采用隨機(jī)方式填充燃料球時(shí)，一般填充系數(shù)為60.5%時(shí)出現(xiàn)的概率最高，這也是選擇HEX 為主要模型的原因。使用MCNP程序?qū)BR 全堆芯進(jìn)行統(tǒng)計(jì)，結(jié)果列于表2。由表2可看出，隨體積填充率的增大，keff增大，可控性能增大。

2）堆芯燃料區(qū)能譜

選取全堆芯的燃料區(qū)作為統(tǒng)計(jì)對(duì)象，PBR堆芯中子能譜示于圖9a，為清晰表示快區(qū)的計(jì)算結(jié)果，放大0.01～1.0 MeV 能區(qū)上的能譜（圖9b）。由圖9a可知，隨體積填充率的增大，熱能區(qū)的中子份額減小，這是因?yàn)槿剂锨蛑g空隙是充當(dāng)冷卻工質(zhì)的氫，其有慢化能力，故間隙減小時(shí)，氫的量減少，導(dǎo)致中子慢化能力減弱，能譜變硬；由圖9b可看出，隨體積填充率的增大，快中子份額略有增加。

3）堆芯中子注量率分布

選取堆芯燃料區(qū)作為統(tǒng)計(jì)區(qū)域，圖10為不同體積填充率下軸向和徑向中子注量率不均勻因子分布。

表2 HEX和RPP兩種情況下的keff Table 2 keffin HEX and RPP

圖9 不同體積填充率下的PBR 堆芯中子能譜Fig.9 Neutron spectra as a function of pebble picking rate for PBR core

由圖10可知，堆芯底部軸向中子注量率不均勻因子較大，這主要是因?yàn)槟Ｐ偷撞糠瓷鋵又行目讖捷^大，中子泄漏較多。此外，體積填充率的變化對(duì)堆芯軸向、徑向中子注量率不均勻因子分布影響較小，這是因?yàn)槿剂锨蚓o鄰燃料區(qū)邊界，所有球與球之間的距離會(huì)整體變化。

圖10 不同體積填充率下軸向和徑向中子注量率不均勻因子分布Fig.10 Radial and axial distributions of neutron fluence rate inhomogeneous factor in different pebble picking rates

3 結(jié)論

本文基于空間核熱推進(jìn)粒子球床堆的設(shè)計(jì)模型，研究了設(shè)計(jì)參數(shù)的變化對(duì)PBR 堆芯反應(yīng)性的影響，得到以下結(jié)論：

1）反射層厚度取決于控制鼓外徑的選擇，在PBR 堆芯設(shè)計(jì)中應(yīng)保持增加徑向堆芯中子泄漏，提高控制鼓價(jià)值，增強(qiáng)安全性這一原則。

2）10B富集度與10B 厚度變化對(duì)控制鼓可控性影響較小，且10B 厚度變化時(shí)，Δkeff存在1個(gè)最優(yōu)化峰值。

3）對(duì)于PBR 堆芯，存在最優(yōu)化柵徑比，考慮小堆芯設(shè)計(jì)思想和可控性能優(yōu)化，在其他條件允許的情況下，可適當(dāng)降低柵徑比的取值；增大堆芯高徑比可增加堆芯徑向中子泄漏，提高控制鼓的控制價(jià)值。

4）隨燃料球體積填充率的增大，堆芯keff增大，中子泄漏減少，能譜硬化，燃料球體積填充率的變化對(duì)堆芯軸向、徑向中子注量率不均勻因子影響較小。

［1］ 廖宏圖.核熱推進(jìn)技術(shù)綜述［J］.火箭推進(jìn)，2011，37（4）：1-11.LIAO Hongtu.Overview of nuclear thermal propulsion technologies［J］.Journal of Rocket Propulsion，2011，37（4）：1-11（in Chinese）.

［2］ 何偉鋒，向紅軍，蔡國(guó)飆.核火箭原理、發(fā)展及應(yīng)用［J］.火箭推進(jìn)，2005，31（2）：37-43.HE Weifeng，XIANG Hongjun，CAI Guobiao.The fundamentals developments and applications of nuclear rocket propulsion［J］.Journal of Rocket Propulsion，2005，31（2）：37-43（in Chinese）.

［3］ BREISMEISTER J F.MCNP—A general Monte Carlo natural particle transport code，LA-UR-03-1987［R］.USA：Los Alamos National Laboratory，2003.

［4］ LUDEWIG H，POWELL J R，TODOSOW M，et al.Design of particle bed reactors for the space nuclear thermal propulsion program［J］.Progress in Nuclear Energy，1996，30（1）：8-19.

［5］ 經(jīng)滎清，楊永偉，許云林.蒙特卡羅方法用于HTR-10首次臨界燃料裝料預(yù)估的校算［J］.核動(dòng)力工程，2005，26（1）：28-34.JING Xingqing，YANG Yongwei，XU Yunlin.Application of Monte Carlo method for verification calculation in fuel loading prediction for first criticality of HTR-10［J］.Nuclear Power Engineering，2005，26（1）：28-34（in Chinese）.

［6］ 常鴻，楊永偉，經(jīng)滎清，等.球床式高溫氣冷堆初次臨界物理計(jì)算的蒙特卡羅方法模型分析［J］.核動(dòng)力工程，2005，26（5）：419-424.CHANG Hong，YANG Yongwei，JING Xingqing，et al.Model analysis of Monte-Carlo method for first criticality physics calculation in pebble bed high temperature gas-cooled reactor［J］.Nuclear Power Engineering，2005，26（5）：419-424（in Chinese）.

［7］ FRATONI M.Development and applications of methodologies for the neutronic design of the pebble bed advanced high temperature reactor（PB-AHTR）［R］.USA：University of California，Berkeley，2008.