基于RMC的CABRI-like基准题计算研究

RMC-based computational study of CABRI-like benchmark

  • 摘要: CABRI研究堆是位于CEA Cadarache的一座池式压水堆,该堆配备独特的He-3瞬态棒系统,可通过快速减压引入剧烈反应性变化,产生与压水堆反应性引入事故(RIA)特征高度相似的极端功率脉冲,是国际上开展燃料瞬态行为研究的重要实验平台。本研究采用国产自主化蒙特卡罗粒子输运程序RMC,针对CABRI-like基准题堆芯结构开展了CSG几何建模与瞬态核热耦合计算,旨在提供一套基于RMC的独立数值结果,并与国际同类型程序的计算数据进行交叉对比,以验证基准题设置的合理性及结果的一致性。计算中,采用RMC时空动力学模块模拟He-3减压引入反应性所引发的功率与温度随时间演化及径向分布,通过Python脚本将实时功率转化为燃料温度变化,并结合内置的TMS温度-截面反馈算法实现核热耦合,从而捕捉反应性引入期间因功率骤升导致的燃料温升、截面变化以及多普勒效应引发的功率偏移。

     

    Abstract:
    Background Reactivity-initiated accidents (RIAs) are critical safety concerns for pressurized water reactors. The CABRI research reactor, equipped with a unique He-3 transient rod system, can generate extreme power pulses closely resembling RIA conditions, serving as a key experimental platform for fuel transient behavior studies. The CABRI-like benchmark problem, established under an international framework, is designed to test transient simulation capabilities of numerical codes, with several internationally recognized codes such as TRIPOLI-4® and Serpent participating in cross-comparisons. The domestically developed Monte Carlo code RMC possesses built-in transient kinetics and temperature-feedback modules, offering an independent computational tool for such benchmark analyses.
    Purpose This study aims to perform independent neutronics/thermal-hydraulics coupled transient calculations for the CABRI-like benchmark using RMC, and to cross-validate the results against those obtained from international peer codes, thereby assessing the rationality of the benchmark and the consistency of the numerical results.
    Methods A detailed CSG (Constructive Solid Geometry) model of the benchmark core was constructed with RMC, incorporating all 40 assemblies including fuel assemblies, BCS control assemblies, and He3 transient assemblies. The LOW transient scenario, driven by a prescribed rapid decrease in He3 density, was simulated for the first 0.15 s, discretized into 150 time steps of 0.001 s each. The RMC space-time kinetics module based on the predictor-corrector quasi-static method was employed. Neutronics/thermal-hydraulics coupling was realized via a Python script that updates the fuel temperature of each assembly at every time step according to the deposited energy, while the built-in Target Motion Sampling (TMS) algorithm provides on-the-fly temperature-dependent cross sections, capturing the Doppler feedback effect. Criticality calculations were performed at the beginning and end states of the LOW, INTER, and HIGH transients for comparison of keff and reactivity insertion against TRIPOLI-4® and Serpent results. The fuel temperature reactivity coefficient was evaluated by varying the fuel temperature from 300 K to 1 800 K.
    Results Criticality comparisons show a maximum keff deviation of less than 0.001 and a reactivity insertion deviation of less than 0.000 4. For the LOW transient, the relative power increases to approximately 6 times the initial 1 MW, and the average fuel temperature rises by about 0.26 K, both in good agreement with the reference data. The computed temperature feedback coefficient is approximately -0.000 02/K, consistent with the compared results. Assembly-level distributions of power, energy deposition, and temperature at t = 0.15 s are presented; the highest temperature (300.32 K) occurs in assembly No. 37, and the three columns of assemblies adjacent to the central experimental zone exhibit above-average temperature increases. A clear correlation between the normalized energy deposition rate and the assembly-averaged temperature is observed.
    Conclusions The RMC code successfully reproduces the fast transient and neutronics/thermal-hydraulics coupling behavior of the CABRI-like benchmark, with results in good agreement with those from international Monte Carlo codes, confirming both the capability of RMC for such analyses and the validity of the benchmark specification. Future work will introduce axial discretization of the fuel rods to obtain axial power and temperature profiles, thereby providing more complete reference data for the benchmark.

     

/

返回文章
返回