通用二维几何的蒙特卡罗热谱程序MCTGP—2D

ＭＣＴＧＰ－２Ｄ是多能蒙特卡罗热谱程序ＭＣＴＧＰ的新版本，用ＦＯＲＴＲＡＮ语言编制，可以在ＰＣ３８６以上的微机上运行，除了保留ＭＣＴＧＰ程序的全部功能外，每个栅元可以由圆、椭圆、任意三角形、任意四边形、矩形相组合而成。可以处理通用二维几何，栅元内几何区数任意。具备微机上栅元几何图形显示功能。对束棒栅和平板栅两个样本方案的计算结果与ＭＣＴＧＰ计算结果相符。     

NPDMC中子和光子耦合输运孔道计算三维蒙特卡罗程序

ＮＰＤＭＣ是一个我国研制的多群中子和光子耦合输运的三维孔道计算蒙特卡罗程序。该程序配备了国内制作的８７群中子和２５群光子截面数据库,在程序设计中,采用了一些独特的方法和技巧,成功地解决了蒙特卡罗方法应有和到辐射屏蔽计算中所遇到的“深穿透”和“小概率”等难题。ＮＰＤＭＣ程序可以计算各种类型的研究性反应堆和动力堆中的管道,束孔和缝隙的中子能量,能谱,剂量率和γ光子通量及剂量率等。     

利用蒙特卡罗方法计算6LiD中子源的产额与能谱分布

利用^6LiD中子源转换靶室将反应堆热中子转换成聚变谱中子，可用来进行聚变中子辐照环境下的材料性能研究。应用蒙特卡罗方法模拟聚变谱中子的产生过程，从理论上验证了这种中子源的可行性。初步计算表明：1个热中子作用在^6LiD源室外表面将在源室内腔中产生0．1314个快中子；所产生的快中子具有很好的聚变谱特点，能量集中在13．5～15．5MeV之间。     

蒙特卡罗方法用于HFETR堆芯γ释热的可行性研究

采用蒙特卡罗方法对高通量工程试验堆(HFETR)堆芯内的γ释热进行了计算,并将计算结果与实测值进行了比较.结果表明:用蒙特卡罗方法计算HFETR堆的γ释热率是可行的,具有满意的计算精度.因此,在实际工程中可采用蒙特卡罗方法来计算HFETR及堆芯内任意位置的γ释热.     

在西门子机上的多能蒙特卡罗程序系统简介

一、概况西门子机(Siemens 7760)上的多能蒙特卡罗程序系统是NOVA-840机上的”直角三角形栅元的热中子通量和热利用因子的多群蒙特卡罗计算程序”的西门子机文本。除了840机上的程序功能外,西门子机上的多能蒙特卡罗程序系统增加以下三个功能:     

蒙特卡罗计算程序高温点截面数据的增补

MCNP程序是目前在反应堆计算分析中广泛应用的蒙特卡罗计算程序。MCNP自带的截面库中的绝大多数核素只包含21℃（294K）下的截面数据。虽然也提供了一些常用核素的高温点截面数据，但温度点划分很粗，间隔跨度很大，如300、400、500、587和600K等温度下的核截面数据。     

栅元热谱程序的加速

对于 THERMOS 型的热中子空间-能量分布计算程序,当进行燃耗计算时,对每一时刻均需重复计算各区之间的首次飞行碰撞几率,花费大量计算时间。本文利用文献[1]提出的一种扰动处理方法,根据初始时刻的碰撞几率结果,能十分迅速而精确地得到其他时刻的近似值。改进后的程序起名为 FTH。     

计算正比管反冲质子谱的相关蒙特卡罗方法

中子进入充有气体的正比管,产生反冲质子。由于正比管的几何尺寸有限,在反冲质子谱的测量中必须考虑壁端效应的影响。蒙特卡罗方法可以比较精确的模拟反冲质子在几何边界的飞越状况,从而得到比较真实的反冲质子谱。本文在给出反冲质子谱解析表达式的基础上提出了一种计算反冲质子谱的相关蒙特卡罗方法,比通常采用的禁区方法有较大的优越性。计算结果与实验结果符合得比较好。相关方法的方差比禁区方法的方差有所减小,平均约减小1/3。     

基于微扰计算的蒙特卡罗临界搜索方法研究

基于微扰计算,提出一种新的蒙特卡罗临界搜索方法。该方法仅需要对初始系统进行一次临界计算,利用临界计算结果及微扰计算结果,对有效增殖因子变化曲线进行估计,并以此为基础得到临界搜索的结果。计算结果表明,该方法可有效、准确地对临界搜索问题进行计算,并有效地克服了传统方法的不足。     

COBRA程序系列中的一个新版本：COBRAIIIC／MIT—2

一、前言二十多年来,在反应堆设计和安全分析领域内,随着计算机的普及和实验研究工作的开展,各种计算机程序应运而生。分析堆芯行为的程序更是层出不穷。其中美国太平洋西北实验室发展的COBRA程序系列是较好的一个。第一个版本COBRA-Ⅰ自1967年发表至今已二十多年,中间不断公布了改进的版本,如COBRA-Ⅱ,COBRA-Ⅲ等。1973年发表的COBRA-ⅢC采用了先进的横向动量方程,较逼真地描述了物理过程,还加强     

空间质子屏蔽材料优化选择的蒙特卡罗模拟

本文考虑在地球内辐射带、轨道高度为556 km、轨道与赤道倾角为90°的飞行器空间的质子环境,对不同的单质屏蔽材料(铝、钨、钽、铅等)和含钨复合涂层,利用SHIELD程序进行蒙特卡罗(M-C)模拟,分析了经过屏蔽材料层后在硅材料中产生的能量沉积、非电离能损( Non-Ionizing Energy Losses,NIEL)、透射质子能谱、二次粒子及其能量和核反应生成的核素等计算结果,并进行了屏蔽材料组分和排序的优化选择.     

ABAQUS程序在微机上的移植和开发

概述了大型力学非线性有限元结构分析通用程序ＡＢＡＱＵＳ微机版的移植和开发工作，解决了不同计算机浮点数格式差别、数据传输中的死机等特殊移植技术难题，并通过二次开发给程序增加了可视化输入输出功能，既降低了分析工作的劳动强度，又保证了计算的可靠性。微机版ＡＢＡＱＵＳ程序通过了原版所附标准例题的验算，结果完全符合。用微机版ＡＢＡＱＵＳ对中国实验快堆（ＣＥＦＲ）壳体进行了强度和变形计算，计算结果亦与ＡＤＩＮＡ程序的计算结果完全一致。     

γ射线自吸收修正因子精确计算方法研究

对具有一定几何尺寸的样品进行γ谱分析时,样品自身对γ射线的吸收影响对核素含量的精确测量。本文在对比国内外关于自吸收修正因子计算方法的基础上,分析了被广泛采用的简化计算模型存在的问题及对修正结果的影响。基于混合蒙特卡罗模拟的思想,提出了自吸收修正因子的精确计算模型,并使用FORTRAN程序进行了随机抽样和积分计算,得到精确的自吸收修正因子。通过加标样品及不同质量标准源的对比测量,将精确计算模型与简化计算模型和无源效率校准软件计算结果进行了对比分析。结果表明,简化计算模型过高评估了自吸收干扰,而精确计算模型计算结果与实验测量结果及无源效率校准软件计算结果具有较好的一致性,相对偏差不大于5%。最后针对核电厂周围环境中主要关注的γ核素,使用精确计算模型得到了不同γ核素在土壤中的自吸收修正因子。     

A power spectrum approach to tally convergence in Monte Carlo criticality calculation

In Monte Carlo criticality calculation, confidence interval estimation is based on the central limit theorem (CLT) for a series of tallies from generations in equilibrium. A fundamental assertion resulting from CLT is the convergence in distribution (CID) of the interpolated standardized time series (ISTS) of tallies. In this work, the spectral analysis of ISTS has been conducted in order to assess the convergence of tallies in terms of CID. Numerical results obtained indicate that the power spectrum of ISTS is equal to the theoretically predicted power spectrum of Brownian motion for tallies of effective neutron multiplication factor; on the other hand, the power spectrum of ISTS of a strongly correlated series of tallies from local powers fluctuates wildly while maintaining the spectral form of fractional Brownian motion. The latter result is the evidence of a case where a series of tallies are away from CID, while the spectral form supports normality assumption on the sample mean. It is also demonstrated that one can make the unbiased estimation of the standard deviation of sample mean well before CID occurs.     

蒙特卡罗方法在中子能谱研究中的应用

针对快中子能谱实验测量数据的修正 ,开发了FAMS MC蒙特卡罗计算程序 ,对 5 9、6 4和1 4 1MeV中子在Be核上产生的次级中子双微分截面的实验结果进行了修正计算 ,并与用MCNP蒙特卡罗程序修正的结果进行了比较。用FAMS MC程序进行中子能谱实验测量数据修正得到满意结果。     

一个低本底中子探测器的效率计算

一、引言自从M.Fleishman和S.Pons报告了他们实现冷核聚变后,引起了科学工作者的极大兴趣并掀起了热潮,我们于1989年4月也开始了这方面的研究工作。我们采用液体闪烁体耦合到光电倍增管来探测中子。ST-451液体闪烁体为φ10.5cm×5cm的圆柱体,用n-γ甄别技术排除γ射线讯号以降低本底。在设置探测阈窗从1.5～3.5 MeV中子能量时,它的平均本底计数为0.025中子/s,在中子源强度为5中子/s时,     

利用LabWindows/CVI实现γ能谱的蒙特卡罗模拟

基于蒙特卡罗模拟思想,利用虚拟仪器开发工具LabWindows/CVI编制软件模拟了γ光子在NaI(T1)闪烁晶体中的输运过程.直接以γ光子在NaI(T1)闪烁晶体中的沉积能量近似对应实际γ能谱仪的脉冲幅度输出,从而得到NaI(T1)探测器中γ射线的响应谱线.结果表明,蒙特卡罗模拟得到的能谱能较好地拟合实验测量得到的能谱.     

X射线荧光分析中原级谱分布的计算

X射线荧光分析中,X射线管产生的原级谱的分布对荧光分析的影响很大。尤其是在元素间吸收增强效应的校正过程中,基本参数法要求准确获得X射线原级谱的强度分布。使用MCNP程序模拟不同加速电压、不同靶材料、不同铍窗厚度等条件下电子打靶后的X射线能谱分布,为X射线荧光分析仪研制过程优化X光管靶材、管压等提供依据,实现高精度的X射线荧光分析。