EBR-Ⅱ无保护失流试验SHRT-45R分析

Benchmark Analysis of EBR-Ⅱ Unprotected Loss of Flow Test SHRT-45R

  • 摘要: 为验证中国原子能科学研究院自主开发的快堆系统分析程序FASYS,对美国钠冷快堆EBR-Ⅱ的SHRT-45R无保护失流试验进行了计算分析。利用FASYS程序对试验的堆芯和一回路进行建模,以两台一回路主泵的转速、中间热交换器二次侧入口流量和温度作为计算边界条件。通过对比分析计算值与试验值发现,以堆芯功率为输入数据时,泵流量和XX09测量组件冷却剂温度计算值与试验值吻合良好,由于采用点模型模拟堆芯上腔室温度,Z形管道进口温度计算值变化较试验值快。在堆芯功率和温度耦合计算情况下,堆芯功率的计算值与实测功率总体上吻合良好,堆芯相对功率低于10%后计算值略有偏大。FASYS程序对SHRT45R试验的分析,验证了该程序的堆芯热工水力模型、一回路热工水力模型、点堆模型,特别是反应性反馈模型。

     

    Abstract: Experimental Breeder Reactor Ⅱ (EBR-Ⅱ) of Argonne National Laboratory (ANL) is a sodium cooled fast reactor. The SHRT45R unprotected loss of flow test which demonstrated the effectiveness of passive feedback in the EBRⅡ reactor was carried out on April 3, 1986. The SHRT45R test was analyzed using the FASYS code which was developed at China Institute of Atomic Energy for thermal, hydraulic, and neutronic analysis of system transients in sodiumcooled fast reactor. The FASYS code contains point reactor model, decay heat model, reactivity feedback model, singlepin model for fuelpin temperature calculation, core and circuit hydraulic model, intermediate heat exchanger model, pump model, decay heat exchanger model and sodiumtoair heat exchanger model, etc. The core and primary circuit of the SHRT45R test were modeled by the FASYS code. Besides, the primary pump speed, the intermediate heat exchanger inlet flow and inlet temperature were taken as the calculation boundary conditions. The calculation was divided into two stages. In the initial simulation stage, the reactivity feedback was not calculated, the core power was taken as the input data, and only the thermal hydraulic behaviors of the core and the primary circuit of SHRT45R were calculated. In the final simulation stage, the reactivity feedback, the core power, and the thermal and hydraulic behavior of the core and the primary circuit were calculated simultaneously. The measured data were compared with calculated results of the FASYS code. It shows that the pump flowrate and XX09 instrumented subassembly coolant temperature predicted by the FASYS code are in good agreement with the measured data when the core power is taken as the input data. The Zpipe inlet temperature is predicted faster due to modeling the upper plenum as a 0D volume. Since the influence of adjacent hightemperature subassemblies around the lowtemperature XX10 instrumented subassembly is not considered, the coolant temperature of XX10 instrumented subassembly predicted is significantly lower than the measured data. In the final simulation of core power and temperature coupled calculation, the calculated value of core power is in good agreement with measured power overall. Specifically, the total core power by the FASYS code is slightly underpredicted between 1 and 4 minutes, and slightly overpredicted after 6 minutes. The increase in total core power after 10 minutes is due to the auxiliary EM pump head increasing, which leads to lower core temperatures and a small reactivity increase. The core thermal hydraulic model, primary circuit thermal hydraulic model, point reactor model, especially the reactivity feedback model of the FASYS code were verified by analysis of SHRT45R test. When the core outlet temperature changes rapidly, the temperature change trend of core outlet plenum simulated by the 0D volume model will be faster than the actual situation. In the future, the thermal stratification model should be considered for core outlet plenum temperature simulation and the subchannel code should be used for the center lowpower subassembly and the surrounding highpower subassemblies transient thermal simulation.

     

/

返回文章
返回