用金活化法测定加速器中子源中子注量及MCNP4C修正

在四川大学720所2.5MeV静电质子加速器上,由核反应7Li(p,n)7Be,T(p,n)3He产生中子,对中国工程物理研究院研制的新型中子探测器进行效率刻度实验中,需要知道探测器位置处的中子绝对注量,为此我们测量了0.165、0.352、0.576、1.400MeV四个能点的中子注量。测量方法采用的是金活化法,在实验测量中,由靶头材料、冷却水层和样品的包层材料等引起的多次散射效应及中子在样品中的自屏蔽效应等均对实验结果产生影响。这些因素在实验中不可避免,也难以通过实验方法扣除,因此用Monte Carlo程序MCNP4C对上述效应进行了修正计算。     

用活化法无损测定金属铀样品中~(238)U的含量

238U作为一种重要的裂变材料,其含量的准确测定在裂变产额数据测量中具有重要意义。在四川大学2.5 Me V质子静电加速器上,利用T(p,n)3He反应产生的483 ke V单能中子照射金属铀样品,对照射后生成放射性核素239Np的特征γ射线进行测量,利用已知的238U(n,γ)俘获截面数据实现了对238U含量的准确测量。对影响测量结果准确性的因素做了细致分析,采用蒙特卡罗方法应用软件MCNPX(Monte Carlo N-Particle e Xtended)对中子的多次散射效应和中子注量衰减效应进行了修正,对γ射线在样品中的自吸收也进行了修正,修正后的实验结果是2.884 2 g金属铀含5.712 8×1021个238U原子,实验结果的不确定度是4.1%。     

中子活化法测定~(232)Th裂变产物产额

采用中子活化法测量了~(232)Th的裂变产物及其累积产额。利用加速器T(d,n)~4He反应产生的14.9 MeV高注量中子长时间照射ThO_2样品,用高纯锗γ谱仪测量其特征γ谱,求得较长半衰期核素~(99)Mo、~(141)Ce、~(143)Ce、~(131)I、~(140)Ba等的裂变产额,实验结果的典型误差为4%。其中,利用MCNP程序对中子的多次散射效应和自屏蔽效应进行修正,同时考虑了中子注量波动及γ射线在样品中的自吸收影响。     

不同描述的中子源模型对CFETR中子学计算的影响

在托卡马克实验装置中,D-T聚变反应释放出的14 Me V高能中子,与周围部件接触会引起包层材料活化、热负载过高等一系列问题,因此在包层设计和优化过程中,相关的中子学计算显得尤为重要。为了研究不同描述的中子源模型对中国聚变工程实验堆(China Fusion Engineering Test Reactor,CFETR)中子学计算的影响,使用基于蒙特卡罗方法的MCNP(Monte Carlo N Particle Transport Code)程序来模拟中子输运过程,分别计算点源、均匀体源、基于L、H、A模约束的中子源模型对不同中子学计算的影响。结果表明,不同描述的中子源模型对氚增殖比影响较小,对中子壁负载和功率密度分布影响比较明显。     

21．6MeV中子与Be作用的次级中子双微分截面测量初步结果

采用飞行时间法,测量了21．6MeV中子与Be作用的次级中子发射谱。测量结果利用Monte—Carlo方法进行了详细的模拟分析,以进行注量衰减、多次散射和有限几何的修正。通过测量谱与模拟谱的比较来确定测量的截面。实验结果以n-p散射截面作为标准进行归一。     

HFETR辐照石墨材料中子注量验证试验

钍基熔盐堆石墨材料辐照考验目标为:中子注量为5×10~(20)cm~(-2)(±15%)(E>0.1 Me V),堆内辐照试验温度650℃(允许偏差±50℃)。为了满足辐照考验要求,在高通量工程试验反应堆(HFETR)第92-I炉的K07孔道进行辐照验证试验。该验证试验辐照装置采用分段构成的型式,主要由辅助密封段、辐照试验段、气管组件3部分构成,辐照罐外围为去离子水,辐照罐内为惰性气体用于控制辐照试验温度。使用MCNP程序对各样品中子注量进行预示计算,同时在辐照装置阳面和阴面都布置了探测器进行中子注量测量。试验表明:在辐照试验过程中,在辐照装置调气系统最佳导热模式下辐照温度略高于上限700℃;利用MCNP程序预示计算中子注量结果为5.7×10~(20)cm~(-2)(E>0.1 Me V),而中子注量测量结果为4.83×10~(20)cm~(-2)(E>0.1 Me V),基本满足石墨材料辐照考验中子注量要求。     

用于(n,γ)辐射俘获研究的大液体闪烁计数器的γ射线探测效率计算

本文叙述了用蒙特卡罗(Monte Carlo)方法计算大液体闪烁计数器探测γ射线的效率和中子在样品中的多次散射修正项的计算步骤,最后给出一些计算结果和讨论。     

1.4MeV-5MeV中子诱发~(238)U裂变产额测量

~(238)U裂变产额测量工作在核数据测量中有着重要意义,本工作利用2.5MeV质子静电加速器产生的1.4MeV-5MeV单能中子诱发238U裂变,通过对裂变产物放射性的测量对裂变产物核素~(135)I、~(133)I、~(105)Ru和~(91)Sr的产额进行了测定。照射过程中中子通量用活化法确定。分析了影响实验测量的多个因素,包括用MCNPX程序对中子在靶头及样品中的多次散射和自屏蔽效应进行了修正,对γ射线在样品中的自吸收进行修正等。得到产额数据典型误差为3.5%,最后把测量结果与已有的裂变产额数据进行比对。     

DT中子下聚乙烯球中~(232)Th(n,2n)反应率的测量与分析

~(232)Th(n,2n)数据对于钍基反应堆研究十分重要,有必要针对其开展中子积分实验,因此,建立了聚乙烯球基准宏观装置,利用活化法在DT中子源下开展了测量~(232)Th(n,2n)反应率的中子积分实验。实验中分别采用了钍粉末以及钍片两种样品形态,以消除实验样品状态对实验结果的影响,将样品置于与D离子入射方向成0°的位置进行辐照,利用金硅面α探测器进行中子产额监测以及中子注量波动监测。辐照结束后,利用高纯锗谱仪离线测量反应产物231Th发射的能量为84.2 ke V的γ射线,得到~(232)Th(n,2n)反应率值。同时使用MCNP(Monte Carlo N Particle Transport Code)在数据库ENDF/B-Ⅶ.1、ENDF/B-Ⅶ.0、JENDL4.0下对实验进行了精确模拟,数据库JENDL4.0下的反应率计算结果与实验符合较好。     

反冲质子探测系统中子灵敏度的Monte Carlo算法

在强流脉冲中子测量工作中,常利用PIN探测器测量中子与聚乙烯中的H核发生弹性散射产生的反冲质子.本文利用Monte Carlo技术进行快速计算的反冲质子探测系统中子灵敏度表达式,给出了算法和计算流程.     

基于D-T中子源的板状铌样品的评价数据基准检验实验

利用中国原子能科学研究院核数据国家重点实验室的脉冲化氘氚聚变中子源产生的145 MeV单能中子,通过飞行时间法,测量了5、10、15 cm厚度板状铌（Nb）样品在与60°和120°两个方向上的泄漏中子飞行时间谱。利用蒙特卡罗模拟软件MCNP 4C进行了泄漏中子飞行时间谱的模拟计算,分别获得了CENDL 31、ENDF/B Ⅷ0和JENDL 40 3个数据库中Nb评价数据的模拟结果。通过各数据库不同能区的模拟结果与实验结果的比值（C/E）,对3个数据库中93Nb与145 MeV中子作用的角分布和双微分截面等相关评价数据进行了检验,重点分析了CENDL 31库的数据。结果表明,CENDL 31数据库的模拟结果在弹性散射能区、非弹性散射能区以及(n,2n)反应能区与实验结果均存在一定的偏差。而JENDL 40数据库除在120°弹性散射能区有高估现象,其他能区的模拟结果与实验结果均符合较好。ENDF/B Ⅷ0数据库的模拟结果除在60°方向弹性散射峰偏低外,其他能量范围的模拟结果均高于实验。     

Measurement and Analysis of Neutron Spectra in Lithium Fluoride and Polytetrafluoroethylene Piles

For the assessment of neutron cross section data for fluorine, angular neutron spectra in the lithium fluoride (LiF) and polytetrafluoroethylene ((CF₂)n) piles were measured in the energy range from a few keV to a few MeV by the time-of-flight method with an electron linac, and the results were compared with those calculated by using nuclear data from JENDL-2 and ENDF/B-IV. Spatial distributions of neutron and X-ray fluxes were also measured in the test piles by the activation method, and the influence of photoneutrons generated in the sample material on the neutron spectrum in each pile was estimated. As a result, it was found that their influence on the neutron spectrum shape below 1 MeV was not so large as was necessary to be taken into account for the present assessment.

The calculated spectra using the JENDL-2 data and the ENDF/B-IV data show generally good agreement with those measured in both piles. However, both calculations underestimate the neutron fluxes around several 100 keV, and overestimate those below 100 keV, when they are normalized in the energy range of 10 keV～1 MeV. Large discrepancies are found between the shapes of the measured and calculated spectra around the resonances of fluorine cross section below 100 keV. The present measurements and analyses suggest that the reevaluations of the inelastic and elastic scattering cross sections below 1MeV and the resonance cross sections below 100 keV are necessary to reduce the observed discrepancies.     

CENDL-3.2与ENDF/B-Ⅷ.0的56Fe评价截面对屏蔽计算影响研究

CENDL-3.2评价库对56Fe非弹性散射截面进行了更新,为了验证其与ENDF/B-Ⅷ.0评价库中截面以及屏蔽计算能力的差异,通过NJOY2016程序对56Fe共振重造后的非弹性散射、总截面等微观截面进行了比较;并制作了多群截面,在56Fe非弹性散射能量范围对以56Fe为主要核素的3个系列屏蔽基准题ILL-Fe、OKTAVIAN-Fe、IPPE-Fe进行了屏蔽计算性能的比较。结果表明,CENDL-3.2评价库的非弹性散射截面在4~12 MeV能量范围内低于ENDF/B-Ⅷ.0评价库的结果;多群截面基准题验证表明,CENDL-3.2评价库计算结果与实验值总体符合较好;对于OKTAVIAN-Fe基准题,在0.1~1 MeV能量范围内两评价库计算结果吻合较好。此外,所有基准题验证结果都有共同的现象,即在56Fe非弹性散射截面占主要贡献的1~10 MeV能量范围内,CENDL-3.2的计算结果比ENDF/B-Ⅷ.0的计算结果偏高。      

中子与58,60Ni反应的理论计算

根据中子与天然Ni及其同位素反应的总截面、去弹截面和弹性散射角分布的实验数据，得到中子的光学模型势参量。应用得到的光学模型势参量，根据光学模型、统一的Hauser-Feshbach和激子模型理论以及扭曲波玻恩近似理论，系统计算和分析了中子与^58,60Ni反应的非弹散射角分布和双微分截面，理论结果与实验很好地一致。     

利用瞬发γ射线法测量209Bi的中子非弹性散射截面

在铅铋快堆、空间堆等先进反应堆中,铋作为冷却剂和慢化剂材料被大量使用,其中子核反应截面,尤其是中子非弹性散射截面的准确性对这些核装置的安全性和经济性等具有重要的影响。基于中国原子能科学研究院HI 13串列加速器瞬发γ射线实验平台,通过瞬发γ射线法测量了209Bi在90、105和120 MeV 3个能点的中子非弹性散射截面。在相对于中子束30°、70°、110°和150°方向放置4个Clover探测器测量中子与样品相互作用产生的γ射线。实验采用相对测量,通过测量中子与48Ti发生非弹性散射发射的9835 keV γ射线的产生截面来确定209Bi的截面。209Bi金属样品的尺寸为50 mm×4 mm,参考样品为1块50 mm×1 mm的天然钛金属样品。将实验测量结果与已发表的实验数据、ENDF/B Ⅷ.0、JEFF 33、JENDL 40、ROSFOND 2010和CENDL 31等评价库数据以及Talys 195程序默认参数的计算结果进行对比,发现趋势一致,90、105 MeV能点的测量结果与Talys 195程序的计算结果符合得更好,120 MeV能点的测量结果与ROSFOND 2010评价库数据符合得更好。     

Neutron Nuclear Data Evaluation for Thorium-232

An evaluation was made on the neutron cross sections, resonance parameters and average neutron yield in fission for ²³²Th in the energy range from thermal energy to 20 MeV. The fission and capture cross sections were evaluated on the basis of the experimental data by converting the relative ratio data into cross section values by making use of recent evaluations for reference cross sections. The total cross section was determined from experimental data in the region from 24 keV to 15 MeV and then extrapolated to lower and higher energies by using the optical model whose parameters had been adjusted as so to reproduce the measured data. The elastic and inelastic scattering, (n, 2n) and (n, 3n) reaction cross sections were calculated by means of the statistical model combined with the optical model. A set of resonance parameters were recommended in the energy range below 3.5 keV and average resonance parameters were deduced in the unresolved resonance region. A value of 7.40 b was chosen for the capture cross section at 0.025 eV, and the picket-fence negative-energy levels were introduced so as to reproduce the non-l/v behavior of the capture cross section in the epithermal region.

The results were incorporated in the Japanese Evaluated Nuclear Data Library, Version 2 (JENDL-2). Comparison was made between the present and other evaluations such as ENDF/B-V and possible reasons for the discrepancy were discussed.     

Neutron Self-Shielding Factor of Scattering-Predominant Resonance Foils

The neutron self-shielding factor of ⁵⁹Co resonance foil as an example of foils whose scattering cross section predominate over their absorption cross sections was obtained by both Monte Carlo method (analog) and the collision probability method for various thicknesses of the foil. Also, the transmission and reflection probabilities of neutrons which have various energies near the resonance energy were obtained, and the effects of multiple scattering of neutrons on the neutron self-shielding factor are discussed.

The neutron self-shielding factors obtained by the Monte Carlo method and by the collision probability method agreed well with each other in the cases Σ t ～ 4.0, in which the Monte Carlo method requires considerably longer machine time. Although for the cases of large Σ t (～4.0) the agreement is not always good because of the flat flux approximation in the collision probability method, the calculation time by Monte Carlo is conveniently short. A combination of both methods is useful in obtaining the neutron self-shielding factor of resonance foils.     

Cross Section Assessment by Measurement of Time Dependent Neutron Spectra from Lithium Assemblies

Time dependent neutron spectra from lithium assemblies were measured to assess the neutron cross sections of ⁷Li in ENDF/B-IV, which is important nuclide for the D-T fusion reactor blanket material. Pulsed neutrons produced by D-D or D-T reaction were used to measure leakage neutron spectra from cubical lithium assemblies as a function of time by the use of NE213 liquid scintillator. Calculations of time dependent neutron spectra were carried out by the Monte Carlo code SIMON, which was prepared for this study. The group constants used in these calculations were processed from ENDF/B-IV data. The calculated and the measured neutron spectra were compared for the following three; a stationary spectrum, spectra at each time interval and decay curves for specified energy groups. Discrepancies between the measured and the calculated neutron spectra were found in these comparisons. In order to assure the cause of these discrepancies, some calculations were carried out with recently measured cross sections of inelastic scattering which excite 0.478 and 4.63 MeV level of ⁷Li. It was concluded that some of the neutron cross section data of ⁷Li in ENDF/B-IV should be ameliorated.