用多群蒙特卡罗方法计算快中子核临界keff和通量密度

用多群蒙特卡罗方法对快中子核裂变系统进行了临界计算。有效增殖因子keff的计算值与实验结果符合。计算所得中子通量密度的空间分布在球形裂变系统中随半径增大单调下降。中子通量密度的能量分布在由高浓缩铀组成的活性区内呈单一能量极大值,其对应能量对于裸球核裂变系统和具有反射层裂变系统分别为0.35MeV和0.25MeV,而在由天然铀组成的反射层中在0.1MeV附近出现能量双峰。由通量密度所得中子能谱在无反射层球形裂变系统中随半径增加变硬,在有反射层球形裂变系统中随半径增加变软。     

用多群蒙特卡罗方法计算快中子核临界K_(eff)和通量密度

用多群蒙特卡罗方法对快中子核裂变系统进行了临界计算。有效增殖因子κｅｆｆ的计算值与实验结果符合。计算所得中子通量密度的空间分布在球形裂变系统中随半径增大单调下降。中子通量密度的能量分布在由高浓缩铀组成的活性区内里单一能量极大值,其对应能量对于裸球核裂变系统和具有反射层裂变系统分别为０．３５ＭｅＶ和０．２５ＭｅＶ,而在由天然铀组成的反射层中在０．１ＭｅＶ附近出现能量双峰。由通量密度所得中子能谱在无反射层球形裂变系统中随半径增加变硬,在有反射层球形裂变系统中随半径增加变软。     

多群输运-燃耗耦合蒙特卡罗计算

介绍了三维多群中子输运-燃耗耦合P3近似蒙特卡罗程序MCMG-BURN，该程序是在连续截面蒙特卡罗程序MCNP和反应堆栅元均匀化程序WIMS基础上发展而来的，使用多群截面为模拟临界实验堆和高通量工程实验堆，取得与MCNP和实验一致的结果，MCMG-BURN具有与MCNP相同的精度，但计算时间较MCNP要少得多。     

基于蒙特卡罗方法的几何因子计算程序

应用蒙特卡罗方法求解几何因子,基于蒙特卡罗方法的几何因子计算程序使用C++语言编写,可用于任意位姿的各种尺寸的圆面探测器对圆面源几何因子的计算。该程序使用了方差减小技巧。通过与国际通用蒙特卡罗计算程序（MCNP5）的计算结果对比,该方法具有结果准确（误差较小）、计算速度快、使用方便等优点。最终使用该程序计算几何因子,与实验数据进行对比,成功验证了中子深度分布分析（NDP）能谱测量系统探测器位姿的准确性（误差5%以内）,并对其移动位置进行修正,发现电机移动20 mm大约会产生1 mm的误差。     

基于蒙特卡罗方法的几何因子计算程序

应用蒙特卡罗方法求解几何因子,基于蒙特卡罗方法的几何因子计算程序使用C++语言编写,可用于任意位姿的各种尺寸的圆面探测器对圆面源几何因子的计算。该程序使用了方差减小技巧。通过与国际通用蒙特卡罗计算程序(MCNP5)的计算结果对比,该方法具有结果准确(误差较小)、计算速度快、使用方便等优点。最终使用该程序计算几何因子,与实验数据进行对比,成功验证了中子深度分布分析(NDP)能谱测量系统探测器位姿的准确性(误差5%以内),并对其移动位置进行修正,发现电机移动20mm大约会产生1mm的误差。     

多物理耦合计算中动态输运问题高效蒙特卡罗模拟方法

多物理耦合计算在众多领域都有重要应用.如果其包含粒子输运过程,用蒙特卡罗方法模拟粒子输运常占据大部分的计算时间,因此多物理耦合计算中动态输运问题的高效蒙特卡罗模拟方法意义重大,其不可避免地依赖于大规模并行.基于动态输运问题的特点,本文提出了两种新方法:一是针对输运燃耗耦合计算的新型计数规约算法;二是动态输运计算样本数自适应算法.两种算法都能在保持计算结果基本不变的前提下使计算时间大幅减少,从而提高了效率.     

提高超高能中子计算效率的蒙特卡罗抽样技巧

利用蒙特卡罗方法模拟中子与中子碰撞及超高能中子的产生与输运,原有抽样方法需要模拟大量样本方能使与超高能中子有关统计量的计算误差达到要求.针对ICF聚变靶中子输运问题的特点,发展一种"加权赌分裂抽样方法",以增加重要区域的抽样数、减少非重要区域的抽样数,同时通过权修正保证计算结果无偏.典型模型数值模拟计算结果显示,该方法能够有效增加聚变区的中子抽样数、增加中子相互碰撞数,使超高能中子通量计算误差显著降低,达到提高超高能中子计算效率的目的.     

Monte Carlo方法计算快中子在样品中多次散射和通量衰减修正因子

本文叙述计算快中子在样品中自吸收通量衰减因子以及(n,γ)和(n,n')反应多次散射修正因子的Monte Carlo方法数学处理过程。     

HL-2A装置上DD热核聚变中子产额和通量的估算

用蒙特卡罗方法,利用MCNP3B编码计算了HL 2A装置芯部离子温度为5keV、离子密度为5×1013cm-3的DD热核聚变中子产额的空间分布,以及高场边A、低场边D、偏滤器的B和C点所在位置各种能量的中子辐射通量和平均通量。计算了由于D+D→P+T反应产生氚,在10%的氚与氘发生DT反应产生中子的条件下,热核聚变中子的产额以及在A、B、C、D各位置通量发生的变化。计算结果为发展新的聚变中子诊断系统、辐射防护和环境评估提供了依据。     

钚体源样品γ能谱计算的蒙特卡罗方法

提出了钚体源样品γ能谱分步计算的蒙特卡罗模型和方法.采用分步和直接计算两种方法对球壳形钚体源样品的γ能谱进行了模拟计算.用实验测量的γ能谱验证了计算结果的可靠性，计算结果与实验测量值在10%的误差范围符合.通过对Steve Fetter模型的模拟计算表明：对于外面包有若干屏蔽层的钚体源样品大系统，采用分步方法进行能谱计算，可以使结果快速收敛.     

蒙特-卡罗法计算黑体空腔有效发射率

介绍了蒙特-卡罗法计算黑体空腔发射率的基本思想和主要算法。首先阐述了计算黑体空腔有效发射率的理论基础,然后分别从正向光线追迹和逆向光线追迹的角度对蒙特-卡罗法的具体算法予以说明,提出逆向追迹思想能方便地计算空腔的定向有效发射率,经过运算可以得到腔积分半球有效发射率和平均垂直有效发射率。最后对比了采用蒙特-卡罗法和其它方法的计算结果,分析表明：蒙特-卡罗法在计算黑体空腔有效发射率上是准确的,劣势在于计算速度慢。     

2.8MeV中子诱发238U裂变的瞬发中子谱高能段实验测量

在中国原子能科学研究院的高压倍加器上利用中子飞行时间方法测量了2.8 MeV中子引起238U裂变的瞬发中子能谱,通过增大样品量和设计合适的屏蔽体提高了效应/本底比,使得测量数据的不确定度满足预期目标（在5.5～14 MeV能区内,能量间隔为0.5 MeV条件下能谱的不确定度小于10%）。将实验几何、中子源能量分布及角分布、探测效率、束流的时间结构等实验条件输入到MCNP程序里,模拟了出射的中子谱,模拟结果与测量结果在不确定度范围内一致,验证了在入射中子能量较低时238>U的裂变瞬发中子谱评价数据的可靠性。The Prompt Fission Neutron Spectrum (PFNS) of 238U induced by 2.8 MeV neutron was measured using the Cockcroft-Walton accelerator in China Institute of Atomic Energy (CIAE). The signal/background ratio was improved by increasing the amount of sample mass and using an appropriate shielding system. The final uncertainty of neutron energy spectrum in 5.5～14 MeV region is less than 10% with a bin size of 0.5 MeV which has reached this project's anticipation. The experimental geometry, the angular distribution and energy distribution of neutron source, the detection efficiency and time structure of deuteron beam were inputted into the MCNP code to simulate the outgoing neutron spectrum. The simulated results agree with the experimental ones within the uncertainty. The result indicates that the evaluated PFNS for 238U at low neutron energy is reliable.     

蒙卡中子输运程序JMCT和子通道热工水力程序COBRA-EN耦合计算

介绍了中子输运蒙特卡罗方法与热工水力耦合计算的流程。开发了一套蒙卡中子输运程序JMCT和子通道分析程序COBRA-EN耦合接口。通过33棒束模型的计算展示了考虑耦合计算和不考虑耦合计算的差异,论证了耦合计算在反应堆分析中的重要性。通过对反应堆组件的模拟计算,测试了耦合计算的正确性。最后分析了蒙卡计算的统计涨落和迭代计算过程中收敛标准的关系,讨论了蒙卡中子输运和热工水力耦合过程中收敛标准设置的方案和可行性。     

蒙特卡洛模拟中伪随机数的STTDMR检验

伪随机数发生器(PNG)是蒙特卡洛法概率模拟中的基本工具,但它也成为影响蒙特卡洛法精确度的最大误差源.提出了一种检验方法——二维介质辐射对称性检验法STTDMR,它可在不同条件下对各种伪随机数程序进行无限制的检验,在更精细的层次上区分伪随机数发生器的优劣.STTDMR不仅可检验伪随机数的特性,还可检验概率模型是否正确.利用STTDMR分析了多种典型的伪随机数产生法,获得一个性能上较优越的组合伪随机数ZH产生程序.     

75As中子辐射俘获截面的实验测量及蒙特卡罗修正

在四川大学2.5 MV静电加速器上， 用活化法测量了29—1 100 keV 能区中子引起的75As(n， γ)76As反应截面， 并用MCNP 4C程序对实验中由中子的多次散射、 注量率衰减效应等引起的偏差进行了修正， 最后将所得结果与已有数据进行了比较。 We have measured the neutron radiative capture cross section of 75As in the energy range of En=29—1 100 keV using the activation technique at 2.5 MV electrostatic accelerator of Sichuan University. The experimental data were then corrected for multiple scattering and self shielding effects with MCNP 4C code. We also compared our results， with the existing measurement， calculation and evaluated data in ENDF/B Ⅵ.     

X光导管传输特性的蒙特卡罗模拟

 X光在导管中的传输性能是设计X光透镜的基础。在研究X光导管传输理论基础上，建立了X光在导管中传输的射线跟踪方法，编制了蒙特卡罗模拟计算程序。对具体的模型，计算给出了导管的传输效率随X光源与导管入口面的距离、导管直径和曲率半径的变化关系。模拟表明：直圆柱导管的传输效率随导管直径增大而减小，随Ｘ光源所在平面与导管入口面之间的距离增大而增大；弯曲圆柱导管的传输效率随其曲率半径的增大而增大，当曲率半径趋向无穷大时，X光在弯圆柱导管的传输效率趋向于直圆柱导管的传输效率。利用编制程序计算得到的传输效率与实验结果符合较好。     

Monte-Carlo simulation of light scattering in fractal porous media

Monte-Carlo simulation of light propagation in a porous medium with a mass fractal morphology was carried out. It was shown that the simulation results can be used to analyze experimental data on light scattering and to study the porous medium structure. __________ Translated from Zhurnal Prikladnoi Spektroskopii, Vol. 75, No. 2, pp. 193–196, March–April, 2008.     

微波器件微放电阈值计算的蒙特卡罗方法研究

为了精准快速地计算微波器件中微放电效应的阈值,在传统蒙特卡罗方法的基础上,提出了三种不同的蒙特卡罗方法,分别对二次电子的初始能量、出射角度和初始相位等参数进行随机,结合四阶龙格-库塔法和Furman模型计算电子的运动轨迹和单次碰撞产生的二次电子发射系数,然后应用不同的方法计算有效二次电子发射系数作为微放电效应的判据.以平板传输线TEM模式为研究对象,采用四种不同的蒙特卡罗方法计算微放电阈值,并与统计模型结果进行对比.结果表明单电子-多碰撞蒙特卡罗方法误差最小,而且稳定性最好.