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

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

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

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

非定常辐射输运问题的蒙特卡罗自适应偏倚抽样

用蒙特卡罗(MC)方法模拟高温、高压、多介质、大变形辐射输运问题时,由于网格体积悬殊,导致各网格通量的统计误差涨落很大,随着时间步的增加,误差积累甚至会导致计算结果失真.为此,发展了针对全局网格计算的源偏倚抽样技巧.用于源偏倚抽样的价值函数基于上个时间步各网格通量及误差,通过加权构造产生,它比传统MC通过解伴随方程获取价值的性价比要高得多.数值试验表明,全局源偏倚抽样通过自适应分配当前时间步各网格的粒子数,有效地降低了当前步重要网格通量误差. 关键词： 非定常 辐射输运 蒙特卡罗 源偏倚抽样     

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

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

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

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

JMCT在CPR1000工程设计中的验证与确认

CPR1000系列反应堆是目前国内广泛应用的第二代压水堆型号之一，蒙特卡罗程序在CPR1000系列反应堆的验证与确认是该程序实现反应堆工程设计应用的关键环节。基于某CPR1000机组实际参数，使用由国内单位研发的蒙特卡罗程序JMCT在该机组开展了粒子输运建模计算，分别进行了临界计算和固定源计算，并进行了验证与确认。对于临界计算，采用JMCT建立了全堆芯pin-by-pin模型，计算了堆芯有效增殖因子和功率分布。对于固定源计算，建立适用于屏蔽分析的反应堆模型和辐照监督管精细结构模型，计算了两个核电机组多个循环的辐照监督管探测器位置累积快中子注量。通过将JMCT的计算结果与参考程序的计算结果、反应堆实际测量值进行了对比，验证了JMCT程序在CPR1000反应堆工程设计中的实际使用效果，证明了JMCT程序具备工程级的计算精度。     

基于共轭离散纵标的减方差方法

对于深穿透类型的屏蔽计算,为了得到较为可信的统计结果,蒙特卡罗方法(MC方法)需要模拟大量的粒子,巨大的计算时间是其存在的主要问题。源偏倚和权窗技巧能够有效降低深穿透问题的计数误差。开展了基于共轭离散纵标(SN)的MC减方差方法研究,根据SN方法的共轭注量率计算并生成了源偏倚和权窗参数,编写了JMCT程序的源抽样子程序,并且在秦山一期测量值基础上进行了验证,成功应用到CAP1400压力容器快中子注量率和堆腔中子和光子剂量率计算中。数值结果表明,对于深穿透屏蔽计算问题,和无偏的MC方法相比,基于共轭SN的MC减方差方法能够在保证结果精度的前提下,提高计算效率1~2个量级。     

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

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

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

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

γ射线探测器响应函数的JMCT模拟计算

采用自主研发的JMCT软件,模拟计算了射线在NaI闪烁探测器中的响应函数。JMCT模拟计算了探测过程中所有的、完整的光子-电子耦合输运过程。其模拟计算的射线能量沉积谱、探测器响应函数,与MCNP的计算结果符合得非常好,与Berger解析法计算的结果也基本符合,表明JMCT软件可以在探测器响应函数方面得到满意的结果,JMCT将在实验核物理、核分析技术应用等方面发挥越来越重要的作用。     

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

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

JMCT优化查找策略的实现

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

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

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

PCM软件包燃料组件弯曲模型的开发及验证北大核心CSCD

燃料组件在反应堆内受压紧力等作用会发生弯曲,该弯曲会显著改变反应堆局部位置的中子慢化。基于中广核核设计软件包PCM中的组件中子截面计算软件PINE和堆芯核设计软件COCO,开发了专门的燃料组件弯曲模型,以分析燃料组件弯曲对堆芯局部功率分布的影响,并和蒙特卡罗软件JMCT做了对比验证计算。计算结果表明,PCM软件包燃料组件弯曲模型的计算结果与JMCT吻合良好,该软件包可以用于燃料组件弯曲的分析计算。燃料组件的弯曲对于堆芯的局部功率分布有显著的影响,需要在设计中予以特别关注。     

中国散裂中子源多功能反射谱仪屏蔽设计

本文介绍了利用蒙特卡罗粒子输运程序MCNPX2.5.0进行中国散裂中子源多功能反射谱仪屏蔽设计的屏蔽需求、辐射源项、计算方法和设计结果等内容.在计算中考虑慢化器泄漏源项、中子导管损失源项等不同辐射源项,使用分步计算和源项角度偏移、源项能量偏移、几何分裂等多种减方差方法,在保证计算结果精度的同时提高计算速度.在谱仪束线传输段、第二中子开关、散射室等的屏蔽计算中,通过比较了不同条件下的所需屏蔽确定最终屏蔽设计,确保谱仪屏蔽外人员可到达区域的剂量低于安全限值2.5μSv/h.     

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

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