SP3 (simplified P3) theory is widely used in LWR (light water reactor) analyses to partly capture the transport effect, especially for pin-by-pin core analysis with pin size homogenization. In this paper, a SP3 co...SP3 (simplified P3) theory is widely used in LWR (light water reactor) analyses to partly capture the transport effect, especially for pin-by-pin core analysis with pin size homogenization. In this paper, a SP3 code named STELLA is developed and verified at SNERDI (Shanghai Nuclear Engineering Research and Design Institute). For SP3 method, neutron transport equation can be transformed into two coupled equations in the same mathematical form as diffusion equation. In this work, SANM (semi-analytic nodal method) is used to solve diffusion-like equation, due to its easy to handle multi-group problem. Whole core nodal boundary net current coupling is used to improve convergence stability in SANM, instead of solving two-node problem. CMFD (coarse-mesh finite difference) acceleration method is employed for 0-th SP3 equation, which represents the neutron balance relationship. Three benchmarks are used to verify the SP3 code, STELLA. The first one is a self-defined one dimensional problem, which demonstrates SP3 method is extremely accurate, due to no academic approximation in one dimensional for SP3. The second one is a two dimensional one-group problem cited from Larsen's paper, which is usually used to verify and prove the SP3 code correct and accurate. And the third one is modified from 2D C5G7-MOX benchmark, whose numerical results indicate that STELLA is accurate and efficient in pin size level, compared to diffusion model.展开更多
使用JMCT(J Monte Carlo Transport Code)对来自MIT的全堆芯pin-by-pin精细建模的国际基准模型BEAVRS的热零功率(HZP)状态进行了模拟计算,并与测试数据进行了对比和分析.比较的物理量包括临界本征值、控制棒价值、反应性温度系数、轴向...使用JMCT(J Monte Carlo Transport Code)对来自MIT的全堆芯pin-by-pin精细建模的国际基准模型BEAVRS的热零功率(HZP)状态进行了模拟计算,并与测试数据进行了对比和分析.比较的物理量包括临界本征值、控制棒价值、反应性温度系数、轴向积分的全堆探测器测量值和不同位置四个组件轴向相对功率密度分布.HZP状态下不同控制棒位置插入和硼浓度的临界本征值计算,JMCT结果与理论值1.000的误差小于0.2%,控制棒价值计算结果与测量值符合.JMCT对轴向积分的探测器径向相对功率分布和四个组件的轴向归一化的探测器的计算结果与测量值进行了比较和分析,计算结果与测量值一致,同时清晰地展示了模型增加格架后,轴向功率曲线在相应位置出现下凹的现象.此外,JMCT给出了轴向积分的组件径向相对功率密度分布和轴向相对功率最大处(Z轴位置)的pin径向相对功率密度分布,并与国际知名程序MC21结果进行了对比,两个图像都符合得非常好.随着计算机与并行计算的高速发展,蒙特卡罗程序开始从传统的反应堆校验工具向反应堆设计工具转变.展开更多
文摘SP3 (simplified P3) theory is widely used in LWR (light water reactor) analyses to partly capture the transport effect, especially for pin-by-pin core analysis with pin size homogenization. In this paper, a SP3 code named STELLA is developed and verified at SNERDI (Shanghai Nuclear Engineering Research and Design Institute). For SP3 method, neutron transport equation can be transformed into two coupled equations in the same mathematical form as diffusion equation. In this work, SANM (semi-analytic nodal method) is used to solve diffusion-like equation, due to its easy to handle multi-group problem. Whole core nodal boundary net current coupling is used to improve convergence stability in SANM, instead of solving two-node problem. CMFD (coarse-mesh finite difference) acceleration method is employed for 0-th SP3 equation, which represents the neutron balance relationship. Three benchmarks are used to verify the SP3 code, STELLA. The first one is a self-defined one dimensional problem, which demonstrates SP3 method is extremely accurate, due to no academic approximation in one dimensional for SP3. The second one is a two dimensional one-group problem cited from Larsen's paper, which is usually used to verify and prove the SP3 code correct and accurate. And the third one is modified from 2D C5G7-MOX benchmark, whose numerical results indicate that STELLA is accurate and efficient in pin size level, compared to diffusion model.
文摘使用JMCT(J Monte Carlo Transport Code)对来自MIT的全堆芯pin-by-pin精细建模的国际基准模型BEAVRS的热零功率(HZP)状态进行了模拟计算,并与测试数据进行了对比和分析.比较的物理量包括临界本征值、控制棒价值、反应性温度系数、轴向积分的全堆探测器测量值和不同位置四个组件轴向相对功率密度分布.HZP状态下不同控制棒位置插入和硼浓度的临界本征值计算,JMCT结果与理论值1.000的误差小于0.2%,控制棒价值计算结果与测量值符合.JMCT对轴向积分的探测器径向相对功率分布和四个组件的轴向归一化的探测器的计算结果与测量值进行了比较和分析,计算结果与测量值一致,同时清晰地展示了模型增加格架后,轴向功率曲线在相应位置出现下凹的现象.此外,JMCT给出了轴向积分的组件径向相对功率密度分布和轴向相对功率最大处(Z轴位置)的pin径向相对功率密度分布,并与国际知名程序MC21结果进行了对比,两个图像都符合得非常好.随着计算机与并行计算的高速发展,蒙特卡罗程序开始从传统的反应堆校验工具向反应堆设计工具转变.