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

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

JMCT软件的自动源偏倚抽样功能及其应用

介绍JMCT（J Monte Carlo Transport）软件的多群中子伴随输运功能以及基于伴随通量的自动源偏倚抽样功能.对某商用压水堆屏蔽模型的模拟计算表明,采用自动源偏倚后,JMCT的模拟结果与实验值符合较好,比MCNP程序采用几何重要性方法的计算效率大幅提高.     

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

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

Euler及N-S方程的二维四边形/叉树混合网格自适应算法

描述了一种基于直角叉树网格的Euler和N-S方程自适应算法.由于考虑了粘性的作用,提出并使用了四边形叉树混合网格的方法,在几何表面附近生成贴体的四边形网格,外流场使用直角叉树网格.采用中心有限体积法,对Euler及N-S方程进行数值求解,对N-S方程的计算中加入了B-L代数湍流模型.在流场中,运用了网格自适应算法,提高了数值计算对激波、流动分离等特性的捕捉和分辨能力.采用上述方法,数值分析了单段和多段翼型的绕流问题.     

基于VENUS-III国际基准模型的JMCT程序验证

基于自主研制的三维中子-光子耦合输运蒙特卡罗通用程序JMCT（J Monte Carlo Transport Code）,采用连续点截面,对国际基准屏蔽VENUS-III模型开展精细建模和中子输运临界及屏蔽计算.临界计算得到系统k_eff、重要区域的通量及能谱.结果表明,JMCT和MCNP程序的重要区域体通量计数吻合较好,偏差均在1%以内.深穿透屏蔽计算采用外源模式,点探测器计数,JMCT计算值与实验测量值偏差在15%以内,满足屏蔽设计对误差的要求.初步验证了JMCT程序临界及屏蔽计算的可用性.     

基于JMCT和FLUENT的物理-热工耦合计算

基于蒙特卡罗程序JMCT2.2和商用CFD程序FLUENT,通过C++语言,采用外耦合的方式开发了一套耦合接口程序。利用JMCT网格和FLUENT计算域之间一一映射的方式完成物理模型和CFD模型之间的网格匹配,实现了物理模型的简单划分和CFD模型网格的精细划分。利用该耦合程序计算了压水堆单根燃料棒模型和3×3带水洞的燃料子组件模型,并与MCNP与FLUENT耦合计算结果进行对比。计算结果表明,使用本文的方法,耦合JMCT程序与FLUENT程序能够用于物理-热工耦合计算并准确提供其反馈参数。     

基于VENUS-Ⅲ国际基准模型的JMCT程序验证

基于自主研制的三维中子-光子耦合输运蒙特卡罗通用程序JMCT(J Monte Carlo Transport Code),采用连续点截面,对国际基准屏蔽VENUS-III模型开展精细建模和中子输运临界及屏蔽计算.临界计算得到系统keff、重要区域的通量及能谱.结果表明,JMCT和MCNP程序的重要区域体通量计数吻合较好,偏差均在1%以内.深穿透屏蔽计算采用外源模式,点探测器计数,JMCT计算值与实验测量值偏差在15%以内,满足屏蔽设计对误差的要求.初步验证了JMCT程序临界及屏蔽计算的可用性.     

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

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

用时序探测事件模拟提升中子多重性计算效率

提出基于时序处理探测事件的中子多重性计数统计方法。在JMCT粒子输运数值模拟程序的基础上研发了用于统计中子多重性计数的专用数值模拟程序JMCT_NMC，实现在线的中子多重性计数模拟功能。展示了利用中子多重性计数算例检验程序，比较了JMCT_NMC与传统算法计算时间消耗的结果。时序探测事件处理方法不需大量存储粒子信息，在解决中子多重性计数模拟受内存限制问题的同时，提升了计算效率。在JMCT_NMC程序中，时序探测事件模拟手段在探测器关联事件模拟、本底分析等领域还有着更广泛的应用前景。     

基于JASMIN的地下水流大规模并行数值模拟

针对具有精细网格剖分、长时间跨度特征的地下水流模拟中计算时间长、存储开销大等瓶颈问题,基于MODFLOW三维非稳定流计算方法,提出基于网格片的核心算法以及基于影像区的通信机制,并在JASMIN框架上研制了大规模地下水流并行数值模拟程序JOGFLOW.通过河南郑州市中牟县雁鸣湖水源地地下水流的模拟,对程序正确性和性能进行了验证;通过建立一个具有精细网格剖分的假想地下水概念模型对可扩展性进行测试.相对于32核的并行程序,在512以及1 024个处理机上的并行效率分别可达77.2%和67.5%.数值模拟结果表明,JOGFLOW具有较好的计算性能与可扩展性,能够有效使用数百上千计算核心,支持千万量级以上网格剖分的地下水流模型的大规模并行计算.     

基于JLAMT可视化建模的JMCT模拟计算

JMCT是中国工程物理研究院高性能数值模拟软件中心粒子输运团队自主研发的三维蒙特卡罗模拟软件,JLAMT为其前处理可视化建模工具。使用JLAMT和JMCT程序对BW,KRITZ,BEAVRS等系列基准题进行了模拟,并对有效增殖系数及计数功能进行了对比分析。其中有效增殖系数计算结果与MCNP5的最大偏差为KRITZ2装载方案19基准题的89.1 pcm,除BEAVRS基准题外的计数结果与MCNP5的偏差基本小于2%,平均偏差在1%左右;BEAVRS基准题功率分布模拟结果与MCNP5及实测值最大偏差分别为7.06%,16.6%,控制棒价值计算与MCNP5及实测值均吻合较好。,     

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结果进行了对比, 两个图像都符合得非常好. 随着计算机与并行计算的高速发展, 蒙特卡罗程序开始从传统的反应堆校验工具向反应堆设计工具转变.     

Some Weighted Averaging Methods for Gradient Recovery

We propose some new weighted averaging methods for gradient recovery, and present analytical and numerical investigation on the performance of these weighted averaging methods. It is shown analytically that the harmonic averaging yields a superconvergent gradient for any mesh in one-dimension and the rectangular mesh in two-dimension. Numerical results indicate that these new weighted averaging methods are better recovered gradient approaches than the simple averaging and geometry averaging methods under triangular mesh.     

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

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

JMCT程序临界安全基准校验计算与分析

JMCT程序是北京计算数学与应用物理研究所自主开发的蒙特卡罗方法粒子输运计算程序。从国际临界安全实验数据库中选取棒栅类、U系、Pu系、233U系实验方案对JMCT程序进行验证,验证涵盖了常见裂变元素（如U,Pu等）从低富集度到高富集度的快、热能区,同时也包含了常见的中子毒物以及反射层材料。将JMCT程序的计算结果与基准实验值进行对比,并且与MCNP,MONK程序的计算结果进行了比较。结果表明,在检验的范围内,JMCT程序具有与国际通用蒙特卡罗方法粒子输运程序一样良好的计算精度。     

JMCT光子-电子耦合输运模拟计算研究

为满足光子-电子耦合输运应用模拟需求,自主研发了JMCT光子-电子耦合输运功能,并实现了多种光电耦合输运处理方式,以便满足不同的工程需求。利用光子-电子耦合输运典型算例验证了JMCT光子-电子耦合输运功能的正确性,并应用到了探测器响应函数模拟、闪光照相模拟计算中,表明JMCT在实验核物理、核分析技术应用等方面发挥着越来越重要的作用。