脉冲源法测量次临界系统裂变衰减特性的理论分析和数值模拟

针对点堆动力学理论解释脉冲源法测试原理时存在的问题, 基于无源中子输运方程分析次临界系统总中子数、泄漏γ射线计数率随时间的变化关系。理论分析表明: 脉冲中子源作用结束后(无源条件下), 在一定时间范围内, 泄漏γ射线计数率和总中子数近似成正比, 两者随时间变化服从近似指数衰减规律, 反映系统本身的裂变衰减特性, 可以由总中子数和γ射线计数率求解瞬发中子衰减时间常数。基于蒙特卡罗程序构造类Godiva裸铀球次临界系统, 模拟脉冲中子源作用下中子和γ射线输运过程, 计算总中子数、泄漏γ射线计数率及两者比值随时间的变化关系, 结果与理论分析一致; 利用脉冲源法由总中子数、泄漏γ射线计数率计算瞬发中子衰减时间常数α₀, 得到与α-k迭代一致的α₀。说明总中子数、泄漏γ射线计数率可以准确反映系统本身的裂变衰减特性。此外, 根据理论分析和模拟计算给出脉冲源法可用数据的时间范围, 分析泄漏γ射线计数率和总中子数比值的影响因素。     

瞬发中子衰减常数的脉冲中子源法测量系统建立

在向没有任何缓发中子先驱核存在的次临界系统中注入一束脉冲中子后，系统中的瞬发中子强度随空间和能量分布的形状在经过很短一段时间的调整后，将A=Aoexp(λ，t）衰减，A为瞬发中子强度，s^-1；A0为初始瞬发中子强度，s^-1；t为瞬发中子衰减时间，s；λ为瞬发中子衰减常数，s^-1。     

散射中子对超瞬发临界系统动态行为的影响北大核心CSCD

在点反应堆理论的基础上,建立了研究散射中子影响超瞬发临界系统动态行为的模拟计算方法,利用CFBR-Ⅱ堆的脉冲实验验证了该方法的正确性。通过分析计算表明,在超瞬发临界状态,散射中子的影响包括增加反应性、使脉冲波形后沿展宽和增加脉冲裂变产额。建立的模拟计算方法为预测脉冲堆外各种物体产生的散射中子的动力学效应提供了途径,对提高脉冲堆的核安全有应用价值。     

基于脉冲中子源法的CiADS反应性测量蒙卡模拟分析

脉冲中子源法(PNS)是加速器驱动次临界系统反应性测量的一种重要技术.利用蒙卡软件建立CiADS次临界反应堆模型,模拟注入脉冲质子束的过程,获得的中子通量的时间演化谱.采用Python语言编程实现脉冲叠加过程,给出稳定缓发中子本底,实现脉冲中子源法的次临界反应堆的反应性测量模拟,给出脉冲周期注入下堆芯不同位置处的中子变...     

西安脉冲堆127I嬗变计算方法与实验验证

为了开展129I的热中子嬗变的研究,在西安脉冲堆上开展了127I靶件辐照实验。以探索实验条件,对127I靶件的嬗变率进行了蒙特卡罗计算,并与实验测量值进行了比对。利用NJOY程序,以ENDF/B VII.0库为基础,制作了127I在西安脉冲堆堆芯辐照温度下的MCNP格式截面库,与MCNP自带库(ENDF/B VI.2)同温度下截面库进行了比较,在不可分辨共振区做了改进,使用新制的截面库,利用MCNP程序对ORIGEN2数据库中的127I辐射俘获截面进行了修正,结合ORIGEN2程序分析了127I靶件在西安脉冲堆实际辐照后的嬗变率和核素的变化,研究了中子能谱和辐照时间对靶件嬗变计算的影响。使用MCNPX自带的燃耗模块CINDER90对127I靶件的嬗变情况进行模拟,结果与ORIGEN2基本一致,与实验数值有2%~3%的偏差,主要原因是MCNP计算中子通量密度存在误差。     

求解中子动力学方程的加权蒙特卡罗方法

为了实现基于蒙特卡罗方法的中子动力学计算,在传统的直接蒙特卡罗动力学方法的基础上,提出了一种加权蒙特卡罗动力学方法。该方法通过引入粒子权重的概念,隐式考虑中子俘获反应和裂变反应过程中中子数目的变化,避免了模拟粒子的数目随时间的变化,降低了统计偏差,消除了程序计算过程中粒子的存库操作,提高了计算精度。基于单能点堆模型,开发了中子动力学计算程序NECP-Dandi,进行了大量数值验证与分析,包括无缓发中子、单组缓发中子、六组缓发中子、正阶跃反应性引入、负阶跃反应性引入、正脉冲反应性、负脉冲反应性和正线性反应性引入等情况。数值结果表明,相比于直接蒙特卡罗动力学方法,加权蒙特卡罗动力学方法在计算结果的精度和计算效率上有较为明显的改进,程序结构更为简洁。     

源强对Rossi-α噪声分析的影响

通过数值模拟手段研究了源强对Rossi-测量的影响。为了定量研究中子源强对Rossi-测量的影响,基于MCNP软件;开发了用于计算Rossi-分布的数值模拟工具。反应堆模型和加速器驱动次临界系统的示意模型的数值模拟结果被用于展示瞬发中子衰减常数与源强之间的关联。数值模拟结果表明,入射源强源中子强度对Rossi-测量结果有显著影响。对于处于次临界状态的反应堆和加速器驱动次临界系统模型,在入射源强较小时,Rossi-方法可以正确给出反应堆中子衰减常数,但当入射源强较大时,测量不能给出正确的中子衰减常数。通过研究源强和测量结果的关系,可以找到能够给出正确测量结果的最大可用源强。通过选择一个可用源强范围内的较强中子源,可以减少测量所需时间。     

300#研究堆首炉中央孔道中子通量密度计算

基于MCNP程序对300#研究堆首炉堆芯进行精细建模,通过并行计算方式得到了实验临界棒位下堆芯的有效增殖因数为1.002 29,与临界值之间的相对误差为0.229%,验证了物理模型的正确性。探讨并解决了并行计算的中断与接续问题,提出了体通量计数与点探测器计数应用中的合理化建议,即对大体积空间计数时尽量使用体通量计数。计算值与实验值对比结果表明:两者在3 MW功率水平下热中子通量密度相差4.6％,符合得较好。     

MCMGP3三维多群P3近似蒙特卡罗中子输运程序基准检验

三维多群P3近似Monte Carlo中子输运程序MCMGP3是为反应堆临界安全计算设计的,它是从连续点截面中－－光孙耦合输运Monte Carlo程序MCNP发展而来,程序用多群截面代替了MCNP程序的连续点截面,但保留了MCNP程序的几何处理能力,计数能力和降低方差技巧及图形功能。能群数可扩展,使用宏观截面或微观截面均可,中子解分布采用P3近似和广义Gaussia求积。多个基准问题结果显示,MCMGP3程序结果与其它方法计算结果符合很好,计算还表明在同样计算精度下,MCMGP3程序的计处时间较MCNP程序少得多。此外,MCMGP3程序还实验了与WIMS程序的连算,可作反应堆全数值模拟。     

氘-氚中子源辐照旋转靶剂量演化模拟计算

首次将蒙特卡罗自动建模系统(MCAM)、蒙特卡罗粒子输运程序(MCNP)及自主研发的活化程序BURNDOT相结合,实现了中子输运、材料活化、光子输运模拟计算的耦合。对14MeV氘-氚(D-T)中子源旋转靶室剂量演化分析,计算了氘-氚中子源辐照旋转靶室的瞬发中子、γ三维辐射剂量场分布及连续辐照8小时后缓发γ剂量变化情况,并考察了材质、栅元、主要核素对缓发γ剂量贡献的影响,得到了旋转靶的剂量时空演化规律,把计算结果与欧洲活化程序FISPACT-2007进行了对比。     

Calculation of the reactor neutron time of flight spectrum by convolution technique

It is a very complex and time-consuming process to simulate the nuclear reactor neutron spectrum from the reactor core to the export channel by applying a Monte Carlo program. This paper presents a new method to calculate the neutron spectrum by using the convolution technique which considers the channel transportation as a linear system and the transportation scattering as the response function. It also applies Monte Carlo Neutron and Photon Transport Code （MCNP） to simulate the response function numerically. With the application of convolution technique to calculate the spectrum distribution from the core to the channel, the process is then much more convenient only with the simple numerical integral numeration. This saves computer time and reduces some trouble in re-writing of the MCNP program.     

基于小型加速器的编码中子源成像研究

编码中子源成像可以在对中子注量率影响不大的情况下大大提高成像的准直比, 从而提高成像质量.北京大学开展了基于小型加速器的编码中子源成像技术研究工作. 不同于已有的基于反应堆的小面积编码板的研究工作, 北京大学建立了基于小型加速器的大面积编码板的编码中子源成像实验平台, 并对加速器中子源上的实验方法和数据处理进行了探索, 对比了重建算法, 获得了初步的重建照片.研究工作表明, 编码中子源成像技术可用于加速器中子源, 但重建图像质量仍须提高. 关键词： 加速器中子源 中子成像 编码源成像 图像重建     

D-T中子与平板型聚乙烯材料作用的次级中子角度谱实验研究

基于反冲质子法建立了一种测量D-T中子与平板型宏观样品作用的次级中子角度谱的实验方法.为保证探测器的能量线性并在较低的中子有效测量下阈(0.5 MeV)情况下获得好的中子-伽马射线甄别性能,采用高、低能段分别测量的方法.采用事件记录法,同时记录了次级中子和伴随伽马射线的脉冲形状甄别和脉冲幅度二维信息,利用基于ROOT数据分析平台编写的离线数据分析程序,完成了伴随伽马射线的挑选和扣除,以及高、低两能段反冲质子谱的拼接,并成功的将神经网络技术应用于中子能谱的解谱,获得了D-T中子与9和18 cm厚平板型聚乙烯材料作用的0.5-15 MeV的次级中子角度谱实验结果.实验模型的MC模拟由MCNP5完成,数据库采用ENDF-VI,实验结果和MC计算结果在实验不确定度范围内一致. 关键词： D-T中子 积分中子学 反冲质子法 次级中子能谱     

墙壁反射中子对脉冲堆波形的影响

墙壁的反射中子会对快脉冲堆的波形产生明显的影响.堆芯中子泄漏后,经过墙壁的反射有一定的概率返回堆芯,由于能量的差异,泄漏中子的返回时间是一个连续的分布.传统的双区模型只考虑了相互作用概率,而没有时间信息,尽管可以很好地解决稳态问题,而无法解决瞬态问题.本文采用等效的方法,把泄漏中子等效为时间相关的堆芯本征源,建立了含有反射效应的时间关联双区模型.求解得到的脉冲波形与CFBR-Ⅱ的实验结果一致,从而合理解释了脉冲波形后沿衰减变慢和坪功率提高的实验现象.     

Chapman-Enskog-maximum entropy method on time-dependent neutron transport equation

The time-dependent neutron transport equation in semi and infinite medium with linear anisotropic and Rayleigh scattering is proposed. The problem is solved by means of the flux-limited, Chapman-Enskog-maximum entropy for obtaining the solution of the time-dependent neutron transport. The solution gives the neutron distribution density function which is used to compute numerically the radiant energy density E(x,t), net flux F(x,t) and reflectivity R_f. The behaviour of the approximate flux-limited maximum entropy neutron density function are compared with those found by other theories. Numerical calculations for the radiant energy, net flux and reflectivity of the proposed medium are calculated at different time and space.     

A conceptual high flux reactor design with scope for use in ADS applications

A 100 MWt reactor design has been conceived to support flux level of the order of 10¹⁵ n/cm²/s in selected flux trap zones. The physics design considers high enriched metallic alloy fuel in the form of annular plates placed in a D₂O moderator tank in a hexagonal lattice arrangement. By choosing a tight lattice pitch in the central region and double the lattice pitch in the outer region, it is possible to have both high fast flux and thermal flux trap zones. By design the flux level in the seed fuel has been kept lower than in the high flux trap zones so that the burning rate of the seed is reduced. Another important objective of the design is to maximize the time interval of refueling. As against a typical refueling interval of a few weeks in such high flux reactor cores, it is desired to maximize this period to as much as six months or even one year. This is possible to achieve by eliminating the conventional control absorbers and replacing them with a suitable amount of fertile material loading in the reactor. Requisite number of seedless thorium-aluminum alloy plates are placed at regular lattice locations vacated by seed fuel in alternate fuel layers. It is seen that these thorium plates are capable of acquiring asymptotic fissile content of 14 g/kg in about 100 days of irradiation at a flux level of 8 × 10¹⁴ n/cm²/s. In summary, the core has a relatively higher fast flux in the central region and high thermal flux in the outer region. The present physics design envisages a flat core excess reactivity for the longest possible cycle length of 6 months to one year. It is also possible to modify the design for constant subcriticality for about the same period or longer duration by considering neutron spallation source at the centre and curtailing the power density in the inner core region by shielding it with a layer of thoria fuel loading.      

Phonon density of states of tetracyanoethylene from coherent inelastic neutron scattering at Dhruva reactor

Inelastic neutron scattering experiments to determine phonon density of states of coherent scattering samples of polycrystalline complex solids are generally intensity-limited and therefore are feasible only at high flux facilities. Phonon density of states of the monoclinic phase of tetracyanoethylene at 300 K, obtained using the medium resolution triple axis spectrometer at the new Indian medium flux reactor Dhruva are reported here. The raw data is converted to the “neutron weighted” phonon density of states by applying suitable corrections. Comparison made with results from a theoretical calculation based on a semirigid molecule model of lattice dynamics is fair. Results from Dhruva are also consistent with that obtained (to be published) at the high flux pulsed neutron source (ISIS) of the Rutherford Appleton Laboratory in United Kingdom.     

Study of a nTHGEM-based thermal neutron detector

With new generation neutron sources, traditional neutron detectors cannot satisfy the demands of the applications, especially under high flux. Furthermore, facing the global crisis in ~3He gas supply, research on new types of neutron detector as an alternative to ~3He is a research hotspot in the field of particle detection. GEM(Gaseous Electron Multiplier) neutron detectors have high counting rate, good spatial and time resolution, and could be one future direction of the development of neutron detectors. In this paper, the physical process of neutron detection is simulated with Geant 4 code, studying the relations between thermal conversion efficiency, boron thickness and number of boron layers. Due to the special characteristics of neutron detection, we have developed a novel type of special ceramic n THGEM(neutron THick GEM) for neutron detection. The performance of the n THGEM working in different Ar/CO_2 mixtures is presented, including measurements of the gain and the count rate plateau using a copper target X-ray source. A detector with a single n THGEM has been tested for 2-D imaging using a ~(252)Cf neutron source. The key parameters of the performance of the n THGEM detector have been obtained, providing necessary experimental data as a reference for further research on this detector.     

Measurements of effective delayed neutron fraction in a fast neutron reactor using the perturbation method

A perturbation method is proposed to obtain the effective delayed neutron fraction β_(eff) of a cylindrical highly enriched uranium reactor.Based on reactivity measurements with and without a sample at a specified position using the positive period technique,the reactor reactivity perturbation △ρ of the sample in β_(eff) units is measured.Simulations of the perturbation experiments are performed using the MCNP program.The PERT card is used to provide the difference dκ of effective neutron multiplication factors with and without the sample inside the reactor.Based on the relationship between the effective multiplication factor and the reactivity,the equation β~(eff)=dκ/△ρ is derived.In this paper,the reactivity perturbations of 13 metal samples at the designable position of the reactor are measured and calculated.The average β_(eff) value of the reactor is given as 0.00645,and the standard uncertainty is 3.0%.Additionally,the perturbation experiments for β_(eff) can be used to evaluate the reliabilities of the delayed neutron parameters.This work shows that the delayed neutron data of ~(235)U and ~(238)U from G.R.Keepin's publication are more reliable than those from ENDF-B6.0.ENDF-B7.0,JENDL3.3 and CENDL2.2.