黎竣松，孫　皖，黃　崗,，萬靈鋒，斯俊平，孫　勝，潘良明

基于MATLAB構建反應堆輻照裝置氣隙控溫研究

黎竣松1,2，孫皖1,*，黃崗1,2,3，萬靈鋒1,2，斯俊平3，孫勝3，潘良明1,2

（1. 重慶大學低品位能源利用技術及系統(tǒng)教育部重點實驗室，重慶 400044；2. 重慶大學核工程與核技術系，重慶 400044；3. 中國核動力研究設計院，四川 成都 610213）

在核反應堆正式試運行之前，需要對各種材料進行輻照考驗，以驗證材料的耐輻照性能，其中溫度控制在試驗中尤為重要，目前輻照考驗裝置的溫度控制還存在一定的問題，為解決核反應堆中的材料輻照考驗裝置的溫度控制問題，利用商業(yè)軟件MATLAB對輻照裝置的考驗段進行了氣隙溫度控制的數值模擬，分別進行了氣隙尺寸計算、氣體組分計算、陰陽輻照計算和階梯氣隙計算，得到了指定輻照溫度需求時，通過調節(jié)氣體混合組分和氣隙大小來實現對輻照裝置的溫度控制改善方案。研究結果表明，當輻照溫度需求在200～400 ℃時，各個計算的溫度控制都能達到指定需求。因此，根據一般反應堆所設計的輻照考驗裝置進行了數值模擬，對提高輻照裝置的溫度控制問題具有一定參考性。

輻照裝置；氣隙溫度控制；陰陽輻照；MATLAB

反應堆用材料在正式使用之前，必須經過堆內輻照考驗，以此來驗證其耐輻照性是否滿足使用要求，材料的堆內輻照考驗需要在一定的技術條件下進行，其中一個重要的指標則是樣品的輻照溫度。如果輻照裝置不能有效的控制輻照溫度，會導致同一輻照罐內的樣品溫度相差過大，使同一材料樣品性能變化不同。

中國的大型壓力容器式輕水試驗堆（HFETR），采用的是非均勻電加熱和階梯氣隙控制技術，試驗溫度可控制在208.5～413 ℃[1]；美國的壓水反應堆（ATR）通過使用精確尺寸和氣體組成的氣隙來進行溫度控制，試驗溫度可控制在300 ℃左右[2]，以及瑞典高通量同位素反應堆（HFIR）通過改變支架與外殼間的氣隙大小來控制溫度[3]；日本材料試驗輕水反應堆（JMTR）和荷蘭的高通量反應堆（HFR Petten）以及法國的（JHR）輻照裝置同樣也采用了氣隙溫度控制，通過在輻照罐內部和外部的間隙中使用合適的氣體混合物來調節(jié)燃料溫度[4-6]；而比利時工程測試反應堆（BR2）則是對輻照罐部分的氣壓和電加熱的調節(jié)，輻照溫度可控制在40～250 ℃[7]，以及韓國的高通量先進中子應用反應堆（HANARO）也是由微型加熱器的輸出和輻照罐中間隙的氦氣壓力進行溫度控制[8]，而（Studsvik R2）輻照裝置是通過改變裝有樣品的輻照NaK以及二級安全殼之間的氦/氖氣體混合物來實現的[9]。

在這些輻照溫度控制方法中，氣隙控溫為主要方式，需在輻照罐內設置氣體層，通過改變氣體層的熱阻達到控制材料輻照溫度的目的。本文基于氣隙溫度控制方法，MATLAB GUI的圖像處理演示系統(tǒng)，利用GUI平臺設計人機交互的圖形用戶界面來進行模擬計算[10]。以上各種輻照試驗溫度大都只能控制在某一范圍，并且溫度控制也不夠精確，本研究依據輻照試驗需求進行模擬，將溫度控制得更為精確些，減小樣品溫差，以滿足更好的輻照試驗溫度條件，對輻照裝置的溫度精確性有一定參考意義。

1　數值計算

1.1　數學模型

本文的輻照罐考驗段結構有兩層氣隙，外氣隙和內氣隙，氣隙里充入惰性氣體（氦氣、氬氣、氖氣或其他氣體），通過調節(jié)氬氣與氦氣的比例來控制溫度。通過在氣體進入氣隙之前進行混合。

假設模擬考驗段內部只存在熱傳導，不存在熱對流和熱輻射（氣隙外殼和內殼為定壁溫條件，氣體混流擬合為函數關系式，對流傳熱可忽略，輻射傳熱基本不存在，也可忽略）。輻照罐考驗段大多為圓柱狀，本文所構建模型也為圓柱狀，并采用MATLAB PDE Toolbox的傳熱方程。

式（1）是一個拋物線偏微分方程描述了給定時間內特定區(qū)域的溫度分布。

式中：——溫度；

——導熱系數；

——比熱容；

——內熱源。

1.2　幾何及邊界條件

輻照裝置主要用于裝載試驗樣品在試驗堆內進行輻照考驗試驗，一般采用分段結構，主要由輻照試驗段、輔助密封段、上下端頭和氣管等組件構成。將樣品固定在輻照罐內，罐內介質最簡單的為空氣，但因空氣在中子照射下會形成酸性氣體，因此一般往罐內充入惰性氣體[11]。本文所建立模型的為輻照裝置試驗段部分，模型外徑31.5 mm，內徑5 mm，模型一共由八部分組成，分別是氣隙外壁、外氣隙、輻照樣品為12.2 mm×8.0 mm的矩形（外面包裹著一層鋯-4合金）、氣管、氣管殼、內氣隙、氣隙內壁、夾塊；由內而外分別是內壁半徑1、內氣隙半徑2、夾塊內半徑3、夾塊外半徑4、外氣隙半徑5、外壁半徑6。幾何條件可根據輻照溫度所需氣隙、所需輻照樣品大小更改，通過改變3與4設置不同氣隙尺寸，具體形狀和幾何條件參數[12]如圖1和表1所示。

圖1　考驗段橫截面

表1　建模材料物性參數表

氣隙中充滿氦氣和氬氣的混流，內壁和外壁材料采用不銹鋼，夾塊部分材料為鋁，內壁和外壁附近為冷卻流道，內外流道通有溫度為50 ℃的冷卻水，模擬過程中內壁面與外壁面受冷卻水影響，溫度變化不大，采用定壁溫的邊界條件，壁面溫度近似等于冷卻水溫度。由于不銹鋼、鋁、鋯-4合金等材料的物性隨溫度變化很小，設置為常物性，模型頂部與底部設置為絕熱條件，內部設置為耦合傳熱。當單一氣體組分時氣隙條件可采用常物性，當混合氣體組分時氣體混流的熱導率和密度通過NIST REFPROP軟件查詢后使用MATLAB擬合為溫度的函數，表達式如下。

式中：——擬合后混合氣體的導熱系數；

——擬合后的系數；

——溫度。

2　計算結果分析

本文分別進行了不同氣隙組分、不同氣隙尺寸、不同釋熱率、陰陽輻照以及階梯氣隙的模擬計算，具體工況如表2所示。

氦氣導熱率0.244 W/（m·K）比氬氣導熱率0.030 W/（m·K）大，增加氣體中的氦氣含量會降低溫度，而增加氬氣含量會提高溫度。所以氬氣在氣隙混流中所占比例不能過高，并且模擬計算后發(fā)現，氦氣在氣隙混流中所占比例越高，溫度可調控范圍越廣，所以將氬氣比例控制在60%以下，為比較單一氣體組分的情況，同時計算了氣隙中全為氦氣的工況。

2.1　網格無關性驗證

基于MATLAB GUI進行網格劃分，網格劃分采用網格尺寸因子乘以，定義見式（3），為確定網格質量及網格數量對計算結果的影響，選定氣隙尺寸為0.3 mm，氣隙組分占比氦氣濃度為100%，氬氣濃度為0%，釋熱率為3 W/g的工況進行網格無關性驗證，圖2左為網格因子乘5.0，圖2右為網格因子乘2.5，網格尺寸因子越小，網格越密。驗證結果如圖3所示，結果表明網格對計算結果無影響。

表2　計算工況表

圖2　考驗段建模網格示意圖（左5.0G，右2.5G）

圖3　網格無關性驗證結果

2.2　不同釋熱率條件計算結果

以輻照需求為200～400 ℃為需求，選定釋熱率為1～6 W/g，當釋熱率為3 W/g時，氣隙尺寸為0.1 mm、0.3 mm、0.5 mm，氣隙組分氦氣占100%，輻照罐截面的溫度云圖和樣品中心線的徑向變化如圖4、圖5所示。

圖4　不同氣隙尺寸下的樣品溫度

從圖4可以看出，隨著氣隙的增大，輻照樣品的溫度逐漸升高，氣隙尺寸為0.1 mm、0.3 mm、0.5 mm時輻照樣品最高溫度分別為140.93 ℃、253.12 ℃、347.38 ℃，最低溫度分別為130.86 ℃、242.79 ℃、336.75 ℃。三種氣隙結構下的樣品溫差大致相同，分別為10.07 ℃、10.33 ℃、10.63 ℃。其他釋熱率工況下，樣品輻照溫度隨著釋熱率增大而上升，樣品溫差也逐漸增大，分別如圖6和圖7所示。

從圖7可以看出樣品溫差在不同釋熱率條件下以3 ℃左右的差別逐漸增長的趨勢，但每種氣隙尺寸的樣品溫差相差不大。當內外氣隙采用相同尺寸時，外流道面積大于內流道面積，大部分熱量都從外流道表面排出，使外流道側的樣品溫度低于內流道側的樣品溫度，導致輻照樣品內側與外側溫度分布不均勻，也就導致樣品溫差增大，為了減小由于內外溫度分布不均勻而造成樣品溫差，可采用不同的內外氣隙尺寸。當內氣隙為 0.05 mm，外氣隙為 0.1 mm時，對比內氣隙為0.1 mm，外氣隙為0.1 mm，低釋熱率3 W/g的樣品溫差從 10.07 ℃下降至7.76 ℃，高釋熱率 6 W/g的樣品溫差從 20.09 ℃下降至15.52 ℃，如圖8所示。

圖5　不同氣隙尺寸的溫度云圖

圖6　不同釋熱率下的樣品平均溫度

圖7　不同氣隙尺寸下的樣品溫差

圖8　不同內外氣隙尺寸的樣品溫度

2.3　不同氣隙組分條件計算結果

氣隙尺寸一般設為0.1 mm至0.5 mm，當氣隙尺寸過小時，達不到所需的輻照溫度，通常需額外采用其他輔助加熱裝置，當氣隙尺寸過大時，輻照溫度過高可能會導致燒毀的風險，因此不能只從釋熱率條件和氣隙尺寸條件去進行輻照試驗，還應從氣隙組分方面考慮，通常氣隙混流中含有兩種或兩種以上的氣體（如氦氣、氬氣、氖氣等）。本文選取了氣隙尺寸為0.3 mm，釋熱率分別為3 W/g和6 W/g的不同氣隙組分計算工況，從圖9中可以看出，隨著氬氣濃度在氣隙組分中增加，樣品溫度逐漸升高，但樣品溫差基本保持不變，低釋熱率時的樣品溫差維持在10.81 ℃，高釋熱率時的樣品溫差維持在20.87 ℃。且隨著氬氣在氣隙混流中所占比例的增高，線段斜率越小，氬氣濃度為60%時比氬氣濃度為50%時的溫度高出118 ℃，氬氣濃度為50%時比氬氣濃度為40%時的溫度高出125 ℃，氬氣濃度為40%時比氬氣濃度為30%時的溫度高出134 ℃；由此類推可說明，當氬氣在氣體混流中所占比例越高，溫度的可調節(jié)范圍越大。反之氦氣濃度所占比例越高，溫度的可調控范圍越小，在此基礎上更改氣隙尺寸和釋熱率，可進一步增加樣品溫度的控制范圍。能夠反向控制氦氣或氬氣濃度來達到所需要的輻照溫度。

圖9　不同氣隙組分的樣品溫度

2.4　陰陽輻照條件計算結果

在進行輻照試驗時，由于輻照裝置與反應堆的結構特點，通常會導致輻照裝置的試驗段一面正對堆芯（陽面），另一面背對堆芯（陰面），在試驗過程中，陽面所接收的射線強度比陰面的高，陰陽兩面的發(fā)熱量不同，導致陰陽兩面輻照樣品的溫度存在一定差異，使材料的性能變化不同，如圖10所示。

圖10　不同輻照方式計算的陰陽面溫度差異

圖中外氣隙尺寸為0.3 mm，內氣隙尺寸為0.1 mm，陰面釋熱率為 3 W/g，陽面釋熱率為6 W/g，陰陽面氣隙組分全為氦氣；陰陽兩面最高溫度分別為255.45 ℃、287.48 ℃，最低溫度分別為246.61 ℃、273.62 ℃；陰面樣品最高溫差8.84 ℃，陽面樣品最高溫差13.86 ℃，陰陽面最高處溫差為32.03 ℃，最低處溫差為27.01 ℃，若采用常規(guī)輻照計算，陰陽面的溫度差異過大，對試驗影響太大，試驗結果不理想，為解決這種影響，可將陰陽兩面采用不同夾塊形式，如以半圓形分開設置氣隙尺寸，或充入不同氣體組分或加入電加熱棒等方式來彌補溫度差異，計算結果如圖10所示。

保持其他條件不變，改變陰陽兩面的氣隙組分，陽面氣隙組分全為氦氣，陰面氣隙組分為80%氦氣、20%氬氣。陰陽兩面最高溫度分別為330.35 ℃、338.75 ℃，最低溫度分別為324.24 ℃、325.63 ℃，此時對比常規(guī)輻照計算陰面樣品溫差為6.11 ℃，陽面樣品溫差為13.12 ℃，陰陽兩面最高處樣品溫差為8.4 ℃，最低處樣品溫差為1.39 ℃，極大地縮減了陰陽面的溫度差異。

2.5　階梯氣隙控溫計算結果

由于γ射線的強度在縱向上呈余弦分布，同時有限長度的裝置必然會產生端部導熱的影響[13]。其軸向釋熱相差大，造成考驗段軸向溫度梯度分布不均勻，為降低軸向溫度梯度，對軸向距離較長的裝置可通過設置軸向階梯型氣隙，改變不同位置的氣隙熱阻來匹配相應的熱流密度。將考驗段劃分為十段，每一段采用不同的氣隙尺寸計算。根據不同氣隙尺寸下樣品溫度與釋熱率的關系，以全是氦氣組分時樣品平均溫度為200 ℃作為輻照試驗需求，求解出每段的氣隙尺寸及各氣隙段的輻照樣品最大溫差見表3所示。

表3　階梯氣隙溫度控制

由表3可見，考驗段氣隙尺寸最小點在L4段，最大點在L10段，氣隙尺寸從中間部分向兩端逐漸增大，每段活性區(qū)域的溫度都穩(wěn)定在200 ℃左右，對于軸向方向較長的輻照罐的溫度差異可控制在5 ℃以內，極大地縮小了軸向方向的溫度差異。

3　結論

（1）建立了輻照罐溫度控制有限元模型，可助于不同輻照溫度要求下的氣隙結構、陰陽輻照罐結構及階梯氣隙結構設計，并可快速預測變輻照溫度條件下的所需調節(jié)氣隙組分；

（2）在200～400 ℃輻照溫度區(qū)間內，通過計算可實現溫度控制達到10 ℃以內；

（3）針對陰陽輻照條件，可通過引入陰陽輻照罐控溫技術控制陰陽面輻照溫度的均勻性，溫度控制可達5 ℃左右；

（4）針對軸向長度較長的輻照裝置，通過選配階梯氣隙，可使溫度控制達到5 ℃以內。

［1］ 劉洋，張亮，楊文華.HFETR材料高溫輻照試驗溫度控制技術研究［J］. 核動力工程，2016，37（S1）：23-25.

［2］ Grover S B，Furstenau R V.The Advanced Test Reactor Irradiation Facilities and Capabilities［R］. Idaho National Lab.（INL），Idaho Falls，ID（United States），2007.

［3］ Qualls A L，Muroga T.A varying temperature irradiation experiment for operation in HFIR［J］. Journal of nuclear materials，1998，258：407-412.

［4］ Muroga T，Nagasaka T，Watanabe H，et al.The effect of final heat treatment temperature on radiation response of V-4Cr-4Ti［J］. Journal of nuclear materials，2011，417（1-3）：310-313.

［5］ Ahlf J，Conrad R，Cundy M，et al.Irradiation experiments on high temperature gas-cooled reactor fuels and graphites at the high flux reactor petten［J］. Journal of nuclear materials，1990，171（1）：31-36.

［6］ Bamber，R. Colin，C. Bignan，G. et al. Conceptual design of test devices for the JHR tailored to the needs of the nuclear fusion community［J］. Fusion Engineering and Design 2017，122：113-123.

［7］ Le Pape Y，Giorla A，Sanahuja J.Combined effects of temperature and irradiation on concrete damage［J］. Journal of Advanced Concrete Technology，2016，14（3）：70-86.

［8］ Choo K N，Kang Y H，Choi M H，et al.Development of Irradiation Capsule Technology in HANARO［C］. Proceeding of the Korean Nuclear Society Spring Meeting，2006.

［9］ Grounes，M. Research reactors for power reactor fuel and materials testingStudsvik’s experience［R］. 1998.

［10］孫鵬崴，王俊，王樹軍，等. 基于MATLAB GUI的圖像處理系統(tǒng)的設計［J］. 計算機技術與發(fā)展，2022，32（04）：215-220.

［11］楊越，單潤華. 輻照技術研究［J］. 核動力工程，1980（04）：61-67+27.

［12］孫勝，楊文華，童明炎，等. HFETR高注量率區(qū)輻照裝置研制與應用［J］. 核動力工程，2016，37（03）：99-102.

［13］李耀鑫，俞爾俊. 輻照樣品的溫度控制［J］. 核科學與工程，1989（01）：57-65+6.

Study on the Gas Gap Temperature Control of Material Irradiation Test Based on MATLAB

LI Junsong1,2，SUN Wan1,*，HUANG Gang1,2,3，WAN Lingfeng1,2， SI Junping3，SUN Sheng3，PAN Liangming1,2

（1. Key Laboratory of Low-Grade Energy Utilization Technologies and Systems，Ministry of Education，Chongqing University，Chongqing 400044，China；2. Department of Nuclear Engineering and Technology，Chongqing University，Chongqing 400044，China；3. Nuclear Power Institute of China，Chengdu of Sichuan Prov. 610213，China）

Before the official operation of nuclear reactors, irradiation tests need to be conducted to verify the irradiation resistance of materials. Precise temperature control is particularly important in the test and there are some problems in the temperature control of the irradiation test device yet. In order to solve the temperature control problem, a numerical simulation of gas gap temperature control in the test section of the irradiator was performed by MATLAB, including gas gap size calculation, gas component calculation, irradiation of front-to-core and back-to-core, and stepped gas gap calculation. According to the simulation, precise temperature control of the irradiator could be achieved by adjusting the gas component and the gas gap size to reach the specific temperature. The results of the study indicated that the temperature control of each calculation can meet the requirements of the irradiation temperature demand within the range of 200～400 ℃. Therefore, numerical simulations based on irradiation test devices of general reactors are informative for the temperature control problem of the irradiators.

Irradiation device; Gas gap temperature control; Irradiation of front-to-core and back-to-core; MATLAB

TL334

A

0258-0918（2023）05-0971-08

2022-10-04

國家重點研發(fā)計劃資助項目（2018YFB1900400）

黎竣松（1999—），男，重慶巫溪人，碩士研究生，現主要從事材料輻照溫度控制方法研究

孫 皖，E-mail：sunwan@cqu.edu.cn