反应堆蒙特卡罗临界模拟中均匀裂变源算法的改进

在反应堆pin-by-pin精细建模及蒙特卡罗模拟计算研究中, 由于不同栅元的功率密度差异较大, 导致蒙特卡罗方法临界计算的样本在不同栅元之间的分配不均衡, 由此引起栅元内的各种计数的统计误差差异较大. 为使大部分栅元内计数的统计误差降至一个合理的水平, 单纯增加总样本已不是一个高效的解决方法. 通过在特定临界计算迭代算法的基础上改进并实现均匀裂变源算法的思想, 对大亚湾压水堆pin-by-pin模型取得了具有较高效率的数值结果. 本工作为具有自主知识产权的蒙特卡罗粒子输运模拟软件JMCT最终达到反应堆pin-by-pin模型(包括一系列国际基准模型)的模拟性能要求提供了一个有效的工具.     

蒙特卡罗临界计算全局计数效率新算法研究

为提高蒙特卡罗临界计算时全局计数的整体效率, 对比分析了新提出的均匀计数密度算法、均匀径迹数密度算法和原有的均匀裂变点算法. 以大亚湾核反应堆pin-by-pin模型的全局体平均通量计数和中子沉积能计数为例, 前两种算法较均匀裂变点算法都获得了整体效率的提高. 上述算法已经在自主开发的并行蒙特卡罗输运程序JMCT上予以实现.     

关于压水堆产武器级钚的模拟计算

防止核扩散是国际社会共同努力的目标,其中武器级核材料的防扩散是重中之重.钚是反应堆的副产品,如果不计较经济效益,利用铀为核燃料的反应堆都可以生产武器级钚.本文基于日本Takahama-3压水堆建立了五个模型,并进行中子和燃耗计算,得到两种燃料棒产武器级钚的条件、燃料棒轴向的燃耗分布、组件内燃料棒燃耗的变化区间和全堆芯燃料棒径向燃耗分布.基于上述模型和计算数据给出压水堆堆芯内含有武器级钚的准确位置和UO_2燃料棒中武器级钚的产量.这种低燃耗的乏燃料给国际核不扩散带来了巨大风险,国际社会应该加强对此类乏燃料的监管.     

蒙特卡罗区域分解并行计算中确保串并行结果一致的伪随机数应用

物理建模的精细化和三维模拟给蒙特卡罗粒子输运计算的规模成千上万倍的增加,甚至超过单核内存的最大规模,仅仅依靠传统的粒子并行蒙特卡罗计算无法实现对模型模拟,区域分解并行是可能的解决方法之一.然而区域分解带来了粒子在各区域间进行迁移,导致现有的伪随机数应用方式无法确保串行计算和并行计算的结果一致.针对这种现象,本文提出赋予粒子随机数属性和动态派生次级粒子随机数的技巧来确保区域分解并行计算的串并行结果一致.     

面向领域的大规模可视建模工具JLAMT的设计与实现

为开展大规模（如反应堆全堆芯pin-by-pin）输运模拟,我们研制了中子-光子耦合输运程序JMCT,快速、方便地建立输运计算模型是缩短问题模拟周期的关键.本文介绍基于领域的可视建模工具JLAMT的开发,为解决大型复杂几何模型的快速建模问题,JLAMT开发了多种快捷建模工具和行业组件工具,实现了层次化几何树的数据结构、显隐式相结合的建模方式、基于特征的转换算法;并自动输出转换为输运计算文件格式GDML.使用JLAMT建立了包括大亚湾全堆芯模型的多个复杂装置模型,JMCT进行输运模拟,验证了JLMAT建模的正确性.     

神光Ⅲ主机装置内爆中子和伽马辐射特性的数值模拟北大核心CSCD

采用三维中子-光子耦合输运蒙特卡罗程序JMCT,对神光Ⅲ主机靶场进行精细几何建模,使用该模型通过计算获得了DD内爆聚变中子及伽马射线在靶场环境的空间分布以及能谱和飞行时间谱。根据计算结果对散射中子和伽马射线的特性进行了讨论,并定量分析了散射中子噪声和伽马噪声分别对下散中子产额测量和高能X射线照相的影响。研究结果表明:必须进行有效的辐射屏蔽设计才能使下散中子产额测量和高能X射线照相的信噪比优于10。     

BEAVRS基准模型热零功率状态的JMCT分析

使用JMCT (J Monte Carlo Transport Code)对来自MIT的全堆芯pin-by-pin精细建模的国际基准模型BEAVRS的热零功率(HZP)状态进行了模拟计算, 并与测试数据进行了对比和分析. 比较的物理量包括临界本征值、控制棒价值、反应性温度系数、轴向积分的全堆探测器测量值和不同位置四个组件轴向相对功率密度分布. HZP状态下不同控制棒位置插入和硼浓度的临界本征值计算, JMCT结果与理论值1.000 的误差小于0.2%, 控制棒价值计算结果与测量值符合. JMCT对轴向积分的探测器径向相对功率分布和四个组件的轴向归一化的探测器的计算结果与测量值进行了比较和分析, 计算结果与测量值一致, 同时清晰地展示了模型增加格架后, 轴向功率曲线在相应位置出现下凹的现象. 此外, JMCT给出了轴向积分的组件径向相对功率密度分布和轴向相对功率最大处(Z轴位置)的pin径向相对功率密度分布, 并与国际知名程序MC21结果进行了对比, 两个图像都符合得非常好. 随着计算机与并行计算的高速发展, 蒙特卡罗程序开始从传统的反应堆校验工具向反应堆设计工具转变.     

JMCT程序在线多普勒展宽研究

在核反应堆堆芯物理分析中, 考虑燃料核素的多普勒效应对反应性的影响是必要的. 基于最小二乘拟合法, 在自主蒙特卡罗粒子输运程序JMCT中开发了在线多普勒展宽功能. 在300-3000 K的温度范围内, 利用核数据处理程序(NJOY)制作了系列温度点下的连续能量核数据. 在统一的能量网格上, 进行最小二乘拟合, 得到拟合参数. JMCT利用这些拟合参数, 快速计算得到任意温度的核数据, 用于输运计算. 计算结果表明, 拟合得到的核数据与NJOY制作的核数据符合很好, 满足精度要求. JMCT临界模型计算结果表明, 拟合参数计算有效增殖因数k_eff与直接采用NJOY加工得到核数据的计算结果吻合.     

蒙特卡罗粒子输运软件JMCT的网格计数功能设计与实现

基于蒙特卡罗输运的网格计数能够统计指定区域内密集的分布计数量.为实现JMCT统计局部计数量分布,设计并实现了JMCT网格计数功能,该功能支持三种正交几何（xyz直角几何、rθz柱几何、rθφ球几何）的均匀及不均匀剖分;介绍xyz直角几何网格计数的算法;基于大亚湾反应堆pin-by-pin模型、Venus模型以及ITER次临界包层模型等初步验证了网格计数功能;应用铀阵列模型,采用若干种网格剖分办法,在单机上对比测试了xyz直角几何下JMCT和MCNP5网格计数功能的串行性能,结果显示JMCT的计数时间较短,具有较高的效率.     

非齐次燃耗方程数值解法

非齐次燃耗方程常用于描述具有显著核素迁移效应的核能系统中核素含量随时间的变化规律.国内外许多燃耗计算程序无法求解方程非齐次项含时的情况.本文在方程非齐次项能够被有限阶关于时间的泰勒展开逼近这一前提下,研究了非齐次项含时情况下方程的两种解法.首先通过Laplace变换推导出了方程基于线性子链方法的解析解形式,然后使用Carathéodory-Fejér方法计算出了方程矩阵级数解的近最佳Chebyshev有理逼近式.将两种方法在燃耗计算程序JBURN中实现,并进行了数值计算,绝大部分计算结果符合很好,部分结果在较长有效数字内仍能保持一致,验证了方法的正确性和精度.同时为求解具有其他非齐次项形式的燃耗方程提供了一种思路.