Simulation of the response functions of an extended rangeneutron multisphere spectrometer using FLUKA

In this paper, the distribution of radiation field in the CSNS spectrometer hall at Dongguan, China, was simulated by the FLUKA program. The results show that the radiation field of the high energy proton accelerator is dominated by neutron radiation, with a broad range of neutron energies, spanning about eleven orders of magnitude. Simulation and calculation of the response functions of four Bonner spheres with a simplified model is done with FLUKA and MCNPX codes respectively, proving the feasibility of the FLUKA program for this application and the correctness of the calculation method. Using the actual model, we simulate and calculate the energy response functions of Bonner sphere detectors with polyethylene layers of different diameters, including detectors with lead layers, using the FLUKA code. Based on the simulation results, we select eleven detectors as the basic structure for an Extended Range Neutron Multisphere Spectrometer (ERNMS).     

LS301快中子探测器的效率刻度

主要介绍了利用252Cf标准裂变中子能谱刻度快脉冲LS301型液闪中子探测器的探测效率的实验方法和结果, 简要介绍了实验数据的处理过程, 得到了阈值分别为0.5, 0.7, 1.0和1.6 MeV, 以及中子能量在10 MeV以下的探测器效率, 并对测量结果进行了误差分析。同时为了检验实验结果的准确性, 将实验结果与理论计算结果进行比较, 两者在不确定度范围内符合得很好。Neutron detection efficiency of LS301 fast neutron detector was calibrated by measuring the neutron energy spectrum of 252Cf source, which has a standard fission neutron spectrum. A low mass, fast ionization chamber is used as the fission fragments detector in the time of flight(TOF) spectrometer and afforded the start signal of neutron flight. The stop signal was offered by the anode of LS301. A measured TOF spectrum was turned to the neutron energy spectrum which will be compared with the standard one. Consequently, the fast neutron detection efficiency of LS301 was obtained. Calibration procedures of experimental and data processing was given. Relative detection efficiencies were obtained separately for threshold settings of 0.5, 0.7, 1.0 and 1.6 MeV for neutron energies under 10 MeV. Experimental results accorded with theoretical efficiency curves which were calculated with the Monte Carlo code NEFF50.     

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.     

First neutron spectrometry measurement at the HL-2A Tokamak

A compact neutron spectrometer based on the liquid scintillator is presented for neutron energy spectrum measurements at the HL-2A Tokamak. The spectrometer was well characterized and a fast digital pulse shape discrimination software was developed using the charge comparison method. A digitizer data acquisition system with a maximum frequency of 1 MHz can work under an environment with a high count rate at HL-2A Tokamak. Specific radiation and magnetic shielding for the spectrometer were designed for the neutron spectrum measurement at the HL-2A Tokamak. For pulse height spectrum analysis, dedicated numerical simulation utilizing NUBEAM combined with GENESIS was performed to obtain the neutron energy spectrum. Subsequently, the transportation process from the plasma to the detector was evaluated with Monte Carlo calculations. The distorted neutron energy spectrum was folded with the response matrix of the liquid scintillation spectrometer, and good consistency was found between the simulated and measured pulse height spectra. This neutron spectrometer based on a digital acquisition system could be well adopted for the investigation of the auxiliary heating behavior and the fast-ion related phenomenon on different tokamak devices.     

PGNAA系统中D-T中子管的慢化装置优化设计

为了提高PGNAA系统中D-T中子管的中子慢化效率,获得更高的热中子产额,借助蒙特卡罗模拟,确定了以铅为中子反射层、5个聚乙烯层和铅层相互交替作为中子慢化层、碳化硼含量为3%的含硼聚乙烯作为中子吸收层以及铅作为γ屏蔽层的中子慢化装置模型。针对中子产额为3×107 n/s 的D-T中子管,该慢化装置输出面低于5 eV中子通量可达5.28×106 n/s,占总中子通量的30.8%,有效提高了中子慢化效率。经过验证模拟结果能够满足实验要求。To improve the moderating efficiency of D-T Neutron Generator in PGNAA system, and get higher thermal neutron yield, the Monte Carlo code MCNP was used to optimize the moderation setup. The lead was selected as neutron reflector and gamma absorber, 5 polyethylene layers and 4 lead layers constituted the neutron moderator and 3% boron-doping polyethylene was selected as neutron absorber. For the yield of 3107 n/s D-T Neutron Generator, this moderation setup can provide a yield of lower than 5 eV of 5.28106 n/s, accounting for 30.8% of total neutron yield, dramatically improves the moderating efficiency. It is proved that the simulation results can satisfy the requirement of PGNAA system by preliminary experimental verification.     

基于D-T紧凑型中子源的快中子照相准直屏蔽体设计及中子束特性模拟研究

中子照相是一种重要的无损检测技术,它能用于火工产品、毒品和核燃料元件等的检测。基于紧凑型D-T中子发生器,完成了一个用于快中子照相的准直屏蔽体系统（BSA）的物理设计。根据D-T中子源的能谱和角分布建立了中子源模型,采用MCNP4C蒙特卡罗程序,模拟了准直屏蔽体系统中中子和γ射线的输运,准直中子束相对于单位源中子的中子注量可以达到9.30×10-6 cm-2,准直中子束中主要是能量大于10 MeV的快中子;在设置的样品平面直径14 cm的照射视野范围,准直束中子注量的不均匀度为4.30%,准直束中中子注量与γ注量的比值为17.20,中子通量和中子注量比值J/Φ为0.992,说明准直中子束有好的平行性;准直屏蔽体外的泄露中子注量率与准直束中子注量率相比降低了2个量级。所设计的准直屏蔽体能满足快中子照相的要求。Neutron radiography is an important nondestructive testing technique. It can be used to detect the explosive devices, drug and the nuclear fuel element, etc. A beam-shaping-assembly (BSA) based on a compact D-T neutron generator is designed for fast neutron radiography in this paper. D-T neutron source model is constructed based on the neutron energy spectrum and angular distribution data. The transportation of neutron and γ-ray in the BSA is simulated using MCNP4C code. The neutron fluence of the collimated neutron beam with respect to the neutron source of the unit source is 9.30×10-6 cm-2. The collimated neutron beams is mainly fast neutrons with energies greater than 10 MeV. In the irradiation field range with a diameter of 14 cm, the neutron fluence uniformity of the collimated beam is 4.3%, the ratio of the neutron fluence to the gamma fluence in the collimated beam is 17.20, and the neutron flux and the neutron fluence ratio (J/Φ) is 0.992 which indicates that the collimated neutron beam has good parallelism. The leakage neutron fluence in outside of BSA is two orders of magnitude lower than that of the collimated neutron beam. The designed BSA can meet the need of fast neutron radiography.     

D-D 中子发生器n-γ 密度测井的蒙特卡罗模拟研究

根据D-D 反应中子的能谱和角分布数据,建立了中子源模型;根据石灰岩地层标准刻度井群数据,建立了井模型。采用MCNP 程序模拟了井中中子和γ 射线的输运,得到了不同地层密度和不同源距处NaI 探测器中的混合γ 射线能谱和非弹γ 射线能谱。在混合γ 射线能谱2.5~5.0 MeV 能区开窗,研究了开窗区混合γ 射线相对计数随源距的变化关系,确定源距应选择在30~80 cm 范围,给出了密度与混合γ 射线计数之间的非线性关系。研究表明,可以利用D-D中子源的混合γ 射线能谱来实现n-γ 密度测井。A D-D neutron source model was developed according to the neutron spectrum and the neutron angular distribution of D-D reaction. A standard calibration well model was built. The transports of the D-D neutrons and γ-rays in the well were simulated using MCNP code. The mixed γ-spectra and inelastic γ-spectra in the NaI detector were obtained for both different distances from a neutron source and different densities. Mixed γ count in the energy range of 2.5 to 5.0 MeV as a function of distance shows that the NaI detector should be located at a distance of 3080 cm from the neutron source. The nonlinear relationships between the density and the mixed γ count were presented in this paper.It was demonstrated that the energy spectrum of mixed γ rays can realize the n-γ density logging.     

Measurement of the neutron spectrum of a Pu-C source with a liquid scintillator

The neutron response function for a BC501A liquid scintillator （LS） has been measured using a series of monoenergetic neutrons produced by the p-T reaction. The proton energies were chosen such as to produce neutrons in the energy range of 1 to 20 MeV. The principles of the technique of unfolding a neutron energy spectrum by using the measured neutron response function and the measured Pulse Height （PH） spectrum is briefly described. The PH spectrum of neutrons from the Pu-C source, which will be used for the calibration of the reactor antineutrino detectors for the Daya Bay neutrino experiment, was measured and analyzed to get the neutron energy spectrum. Simultaneously the neutron energy spectrum of an Am-Be source was measured and compared with other measurements as a check of the result for the Pu-C source. Finally, an error analysis and a discussion of the results are given.     

8.19 MeV中子与9Be作用的次级中子双微分截面测量

用中国原子能科学研究院HI13串列加速器上的多探测器快中子飞行时间谱仪, 测量了8.19 MeV中子与9Be作用时, 从20°到160°区间26个角度的次级中子双微分截面。测量截面以np散射截面作为标准进行归一。实验结果用Monte Carlo方法进行了中子注量率衰减、多次散射和有限几何修正, 并用MCNP4C程序对所用的Monte Carlo程序进行了验证。测量结果与评价数据以及其它实验室的数据进行了比较。The secondary neutron emission double differential cross section of 9Be induced by 8.19 MeV neutron was measured at 26 different angles from 20°to 150°by using the multi detector fast neutron TOF spectrometer at the HI 13 Tandem Accelerator at China Institute of Atomic Energy(CIAE). The results were normalized to np scattering measurement. A special Monte Carlo code which was validated with the MCNP 4C code was employed to analyze the measured data for the corrections of neutron flux attenuation, multiple scattering and finite geometry. The measured results were compared with the evaluated data and the other measurements.      

(n,γ)反应实验研究中的中子屏蔽设计

用射线全吸收型装置(Gamma-ray Total Absorption Facility,GTAF),可以对中子俘获反应截面进行高精度测量。为了降低实验本底,实验中需要对源中子进行准直和屏蔽,还要对被样品散射的中子进行吸收以减少它们进入探测器后所形成的干扰。采用MCNP对中子的准直器、屏蔽体和中子吸收体进行了模拟设计,中子准直屏蔽体材料选用含硼聚乙烯(BC4 的质量分数为3%) 和铅。准直孔直径为13 mm,长度为500mm,经准直后样品处中子束斑坪顶直径为21 mm。中子吸收体材料选用聚乙烯和碳化硼,吸收体球壳内腔半径30 mm,聚乙烯壳层厚度60 mm,碳化硼壳层厚度10 mm,被样品散射的中子经吸收体后衰减93.7%。Neutron capture cross section can be measured by Gamma-ray Total Absorption Facility (GTAF) with high precision. To reduce the background of experiments, the neutron source must be collimated and shielded, and the neutrons scattered from the sample must be absorbed to minimise interference after they go into the detector. The shield, collimator and absorber were simulated and designed with MCNP code. Boron-ontainingpolyethylene with 3% BC4 and lead are used as the materials for the neutron collimator and shield. The diameter of the collimating aperture is 13 mm, and the length of the collimator is 500 mm. After being collimated, the diameter of neutron beam plateau at the sample position is 21 mm. The neutron absorber is made of polyethylene and BC4, and the thickness of polyethylene shell and BC4 shell are 60 and 10 mm, respectively. The simulated result shows that neutrons scattered from the sample can decay 93.7% through the neutron absorber.     

Simulation of the response functions of an extended range neutron multisphere spectrometer using FLUKA

In this paper, the distribution of radiation field in the CSNS spectrometer hall at Dongguan, China, was simulated by the FLUKA program. The results show that the radiation field of the high energy proton accelerator is dominated by neutron radiation, with a broad range of neutron energies, spanning about eleven orders of magnitude.Simulation and calculation of the response functions of four Bonner spheres with a simplified model is done with FLUKA and MCNPX codes respectively, proving the feasibility of the FLUKA program for this application and the correctness of the calculation method. Using the actual model, we simulate and calculate the energy response functions of Bonner sphere detectors with polyethylene layers of different diameters, including detectors with lead layers, using the FLUKA code. Based on the simulation results, we select eleven detectors as the basic structure for an Extended Range Neutron Multisphere Spectrometer(ERNMS).     

Experimental validation of a Tokamak neutron spectrometer based on Bonner spheres

In order to realize real-time fusion neutron spectrum diagnosis for the HL 2A Tokamak, a Bonner Sphere Spectrometer(BSS) array has been developed, consisting of eight polyethylene spheres(PS) with embedded3 He proportional counters.To validate its spectrometric capability, spectrum measurement of an241Am-Be neutron source was carried out and is described.The Monte Carlo code Geant4 was used to calculate the response functions,taking this interference into consideration.Finally, the neutron spectrum was unfolded in the energy range from10-9MeV to 20 MeV.The unfolded spectrum has remarkable consistency with the ISO 8529-1 standard241Am-Be neutron spectrum which is a preliminary demonstration that this BSS is reliable and practical.     

Calibration of a Bonner sphere extension (BSE) for high-energy neutron spectrometry

In a recent work, we constructed modular multisphere system which expands upon the design of an existing, commercially available Bonner sphere system by adding concentric shells of copper, tungsten, or lead. Our modular multisphere system is referred to as the Bonner Sphere Extension (BSE). The BSE was tested in a high energy neutron beam (thermal to 800 MeV) at Los Alamos Neutron Science Center and provided improvement in the measurement of the neutron spectrum in the energy regions above 20 MeV when compared to the standard BSS (Burgett, 2008 and Howell et al., 2009).However, when the initial test of the system was carried-out at LANSCE, the BSE had not yet been calibrated. Therefore the objective of the present study was to perform calibration measurements. These calibration measurements were carried out using monoenergetic neutron ISO 8529-1 reference beams at the Physikalisch-Technische Bundesanstalt (PTB), Braunschweig, Germany. The following monoenergetic reference beams were used for these experiments: 14.8 MeV, 1.2 MeV, 565 keV, and 144 keV. Response functions for the BSE were calculated using the Monte Carlo N-Particle Code, eXtended (MCNPX). The percent difference between the measured and calculated responses was calculated for each sphere and energy. The difference between measured and calculated responses for individual spheres ranged between 7.9 % and 16.7 % and the arithmetic mean for all spheres was (10.9 ± 1.8) %. These sphere specific correction factors will be applied for all future measurements carried-out with the BSE.     

Iterative unfolding for Bonner sphere spectrometers using the MSANDB code – Sensitivity analysis and dose calculation

Since many decades, Bonner sphere spectrometers (BSSs) are routinely used for assessment of neutron spectra over a wide energy range from some meV to GeV. Typically, a spectrometer consists of a neutron detector sensitive to thermal neutrons located inside moderating polyethylene (PE) spheres of various sizes. Based on the response functions of the detector/sphere systems that must be calculated by Monte Carlo codes, an unfolding procedure is applied to deduce the present neutron spectrum from the count rates of the detectors. To start the unfolding procedure, a guess spectrum is required that includes some prior information on the physics of the investigated neutron spectrum, and that is iteratively modified to match the measured count rates. For the present investigation, a BSS-system consisting of an ³He proportional counter and seventeen spheres containing PE or a combination of PE and lead is used. The system is used for example to measure secondary neutrons from cosmic radiation at mountain altitudes, or from high-energy accelerators outside the radiation shielding. A systematic study was performed to quantify the influence of the chosen guess spectrum and the number of iteration steps on the unfolded neutron spectrum, and on integral quantities deduced such as total neutron fluence or ambient dose equivalent. It turned out that none of these changes resulted in dose quantities that were more than a few percent different to those deduced when the optimised start spectrum was used. Similarly, use of the two different response matrices available for our BSS system provided similar values for the ambient dose equivalent.     

一种宽能谱多球中子谱仪能量响应的MC模拟

利用蒙特卡罗程序FLUKA模拟计算了聚乙烯慢化球和辅助材料慢化球对低能中子到高能中子的响应函数曲线。结果表明,对纯聚乙烯球来说,随着聚乙烯层厚度的增加,响应曲线峰逐步右移,峰值在高能区有所下降,对20 Me V以上的中子,无论纯聚乙烯球的尺寸有多大,其响应均下降到很低的程度;对辅助材料慢化球来说,中子能量小于1 Me V时,辅助材料慢化球与聚乙烯慢化球的响应曲线相似,但当中子能量大于20 Me V时,中子与辅助材料层发生(n,xn)反应,慢化球的响应呈显著上升趋势。分析计算结果,最终能够确定宽能谱多球中子谱仪的尺寸组合。     

Neutron,photon and proton energy spectra at high altitude measured using a phoswich-type neutron detector

Neutron energy spectrum from 7 to 180 MeV, photon energy spectrum from 4 to 50 MeV and proton energy spectrum from 94 to 145 MeV were measured simultaneously using a phoswich-type neutron detector with particle discrimination methods at atmospheric depth of 249 g/cm², a vertical cut-off rigidity of 10.2 GV and at a heliocentric potential of 312 MV. We compared our results with other measured and calculated particle energy spectra. Our measured results give a large, sharp neutron peak around 70 MeV, although Bonner balls show a broad peak around 100 MeV due to low energy resolution. The measured photon and proton spectra are between the calculated energy spectra. This onboard study provides the first experimental neutron energy spectrum over 10 MeV with a high-energy resolution.     

Monte Carlo calculations and experimental calibrations of Bonner sphere systems with a new cylindrical helium-3 proportional counter

The response functions to mono-energetic neutrons from thermal energy to 20 MeV of eight polyethylene moderating spheres with a cylindrical proportional counter were calculated with the MCNP-4B code for each moderator and confirmed by measurements in wide neutron fields (ISO), as Am–Be and ²⁵²Cf, and with the CANEL and SIGMA IPSN Cadarache (France) facilities, which simulate realistic field and thermal neutron response. With the described spectrometry system, we have carried out neutron measurements at different points from the annular zone inside the containment building of the Vandellòs II Nuclear Power Plant (Tarragona, Spain). The aim of this work is to present the results obtained during the sphere spectrometer development, with reference to the experimental and theoretical aspects and to show that our calculated spectra correspond to points of measurement compared with those obtained using the MORELPA code.     

聚变中子能谱测量系统脉冲中子灵敏度的实验研究

为多种复杂环境下的稳态和脉冲DT聚变中子能谱测量建立了一种灵敏度优化反冲质子磁谱仪. 使用成像板和同位素α源测量了谱仪的反冲质子能量-位置投影关系. 利用稳态加速器中子源平台、通过单粒子计数方法结合三维带电粒子输运程序模拟,研究了谱仪脉冲中子灵敏度能量响应. 通过高探测效率参数设置使谱仪对DT中子的探测效率达到2×10^-5 cm²水平,从而在较弱中子源上获得了较高统计精度实验数据. 程序模拟结果与谱仪α粒子刻度和DT中子标定实验结果取得了良好的一致性,可由此发展精细解谱技术,以提高脉冲中子能谱测量的灵敏度和能量分辨. 关键词： 聚变中子能谱 磁反冲质子 脉冲中子灵敏度 粒子输运     

用于中子能谱测量的反冲质子磁谱仪

根据神光-Ⅲ原型内爆实验条件下初级中子能谱诊断的要求,提出了反冲质子磁谱仪原型的设计方案,并通过蒙特卡罗程序Geant4对设计方案进行了仿真模拟。该质谱仪的测量能区在6~30 MeV,测量产额范围大于1012,能够测量内爆实验的全中子能谱。模拟结果表明:设计的谱仪对于14 MeV的中子探测效率大于10-11,能量分辨达到1%,信噪比大于10,满足神光-Ⅲ原型装置内爆实验离子温度诊断的要求。