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

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

JMCT反应性微扰计算功能的实现

蒙特卡罗中子输运程序用于反应堆设计与安全分析计算及校验时,微小的反应性扰动容易被统计涨落所掩盖,因此一般通过微扰计算处理此类问题。基于微分算子方法在JMCT程序中实现了反应性微扰计算功能。为保证计算的精度,JMCT考虑了非源扰动的高阶效应以及一阶裂变源扰动效应。选取了具有全局扰动、局部扰动的快中子与热中子裂变系统对JMCT反应性微扰计算功能进行了测试。测试系统覆盖了正负反应性变化。     

JMCT优化查找策略的实现

分析连续能量点截面中子输运程序JMCT中区间查找的时间复杂度.针对主能量网格,裂变释放中子数和不可分辨共振网格进行区间映射加速,并以等差和等比序列网格为例,对比分析不同散列网格的质量.最后选择寿期初与寿期中燃料棒栅元的k_∞问题进行对比测试.结果显示,对于核素数目较多的输运问题,优化查找策略具有很好的加速效果.     

非定常输运基于香农熵的自动调整样本数策略

对于非定常输运问题提出了一种基于香农熵的自动调整样本数策略。将每一计算步的总样本数划分为若干批并逐步模拟每批中的粒子, 可以在每批粒子模拟结束后通过计算得到该时间步幸存粒子属性分布对应的香农熵值。采用在线收敛性诊断方法, 一旦通过香农熵值序列判断对应的幸存粒子属性分布已经收敛, 则可以提前结束本时间步的计算。对一个空间一维非定常输运模型的计算结果表明, 该策略可以显著减少每一计算步的实际样本数且保持最终的结果基本不变, 从而减少了计算时间, 提高了计算效率。     

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

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

组合几何蒙特卡罗粒子输运支撑软件框架JCOGIN介绍

介绍了JCOGIN支撑软件框架层次式、模块化的体系结构及其核心数据结构。这些使得JCOGIN软件框架能够支撑大规模复杂几何的蒙特卡罗（MC）粒子输运并行计算。同时介绍了在JCOGIN支撑软件框架下研发的JMCT程序,并进行了性能测试,在2万处理器核上模拟20亿粒子,并行效率达到70%。     

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

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

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

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

蒙特卡罗粒子输运软件JMCT抽样工具库设计

为通用型蒙特卡罗粒子输运程序JMCT设计了抽样工具库,通过两种技术途径提供各分布的抽样。一是针对各种常见分布提供特定抽样子程序;二是提供一个通用型的抽样子程序,可以实现任意离散分布和任意一维有限区间上连续分布的自动抽样。在设计任意一维有限区间上连续分布的自动抽样工具时利用了部分开源代码,利用其功能提供给用户最大的方便性。对抽样工具库的检验表明,其可以正确、方便地实现各种输运模拟中常见分布的抽样。     

非定常输运问题适应于消息传递并行编程环境的香农熵计算方法

在多计算步的非定常输运问题的蒙特卡罗模拟中,为自动调整每一步的样本数以获得较高的计算效率,可以有多种准则.一种可选的方法是在每一步每隔若干样本监测一次系统中未死亡粒子属性分布对应的香农熵的收敛情况以决定何时停止追加样本,此种方法需要在每一步频繁计算香农熵值.由于在MPI消息传递并行编程环境下香农熵的经典计算方法必须广播大量的数据,导致每一步的计算时间随香农熵计算频率的提高而快速增大,这显然是不能满足实际需求的.本文提出了一种适应于消息传递并行编程环境的香农熵计算新方法,该方法计算得到的香农熵值并不等价于经典方法,但二者之间的差别会随着样本数的增加而趋于零.新方法的最大优势是高频计算香农熵值的时间代价大为降低,为最终实现基于香农熵收敛判断的每步样本数的自动调整奠定了必要的基础.