EXPERIMENTAL ANALYSES OF SPALLATION NEUTRONS GENERATED BY 100 MEV PROTONS AT THE KYOTO UNIVERSITY CRITICAL ASSEMBLY
- DOI : 10.5516/NET.08.2012.003
- Author: PYEON CHEOL HO, AZUMA TETSUSHI, TAKEMOTO YUKI, YAGI TAKAHIRO, MISAWA TSUYOSHI
- Organization: PYEON CHEOL HO; AZUMA TETSUSHI; TAKEMOTO YUKI; YAGI TAKAHIRO; MISAWA TSUYOSHI
- Publish: Nuclear Engineering and Technology Volume 45, Issue1, p81~88, 25 Feb 2013
Neutron spectrum analyses of spallation neutrons are conducted in the accelerator-driven system (ADS) facility at the Kyoto University Critical Assembly (KUCA). High-energy protons (100 MeV) obtained from the fixed field alternating gradient accelerator are injected onto a tungsten target, whereby the spallation neutrons are generated. For neutronic characteristics of spallation neutrons, the reaction rates and the continuous energy distribution of spallation neutrons are measured by the foil activation method and by an organic liquid scintillator, respectively. Numerical calculations are executed by MCNPX with JENDL/HE-2007 and ENDF/B-VI libraries to evaluate the reaction rates of activation foils (bismuth and indium) set at the target and the continuous energy distribution of spallation neutrons set in front of the target. For the reaction rates by the foil activation method, the C/E values between the experiments and the calculations are found around a relative difference of 10%, except for some reactions. For continuous energy distribution by the organic liquid scintillator, the spallation neutrons are observed up to 45 MeV. From these results, the neutron spectrum information on the spallation neutrons generated at the target are attained successfully in injecting 100 MeV protons onto the tungsten target.
ADS , KUCA , Spallation Neutrons , High-energy Protons , Reaction rates , Continuous Energy Distribution
The Kyoto University Research Reactor Institute (KURRI) has engaged in innovative project of research and development of an accelerator-driven system (ADS) at the Kyoto University Critical Assembly (KUCA)- and an external neutron source. The main objective of the project is to develop a next-generation neutron source with the use of the ADS with high-energy protons obtained from the fixed field alternating gradient (FFAG) accelerator,- instead of the conventional research reactor (the Kyoto University Research Reactor: KUR; 5 MW) at KURRI. The basic concept of ADS was originally proposed as an energy amplifier system that couples with a high-power accelerator and a thorium sustainable system. Another possible function of ADS was resolving the issue of nuclear transmutation of radioactive materials, including minor actinides and long-lived fission products. The current research activities on ADS involve mainly experimental feasibility studies with the use of critical assemblies and test facilities: MASURCA,- YALINA-booster and - thermal,- VENUS-1, GUINEVERE and KUCA.
In the ADS facility, a study on target design was requisite experimentally and numerically for examining the characteristics of spallation neutrons generated by injecting the high-energy protons before reactor physics experiments on feasibility studies of ADS. Especially dimensions of the target (size and thickness) needed to be investigated, before the actual injection of high-energy protons, because they are expected to vary in accordance with the spot and the energy of protons. The main characteristics of spallation neutrons analyzed experimentally and numerically at another irradiation facility have contributed notably to the present study. On the basis of those valuable results, the present experimental study of spallation neutrons generated at the tungsten target was conducted, focusing on the neutron spectrum (reaction rates and continuous energy distribution of neutrons) at the location of the target. Special attention was paid to the reaction rates (neutron yield and profile) with the use of bismuth (Bi) and indium (In) activation foils, and to the continuous energy distribution of neutrons using organic liquid scintillator. These experiments on continuous energy distribution of neutrons were emphasized to be carried out near the location of the target, whereas neutron spectrum measurements are usually conducted by the time of flight method at a location far from the neutron source. Further, the neutron spectrum measurement was considered difficult near the location of the target, because the accuracy of neutron detection is dependent on two factors: the distance and the time of flight of neutrons.
The neutron spectrum was considered an important parameter for recognizing information on neutron energy at the target; as a result, the experimental analyses at the target could significantly affect the neutronic characteristics of the core by numerical analyses in the feasibility study on ADS. Also, the neutron spectrum evaluated by reliable methodologies could contribute to accurate prediction of reactor physics parameters in the core through numerical simulations of desired precision. The objective of this study was to analyze experimentally the neutron spectrum of spallation neutrons generated by 100 MeV protons at the location of the target in the ADS facility, before the actual operation of ADS and the improvement of numerical simulations. The neutron spectrum experiments with 100 MeV protons generated by the FFAG accelerator are shown in Sec. 2 and include descriptions of the experimental settings. The results of experimental and numerical analyses by the Monte Carlo calculation code combined with nuclear data libraries are presented in Sec. 3 and the conclusions are summarized in Sec. 4.
High-energy protons were generated by the FFAG accelerator under the following parameters: 100 MeV energy; 30 pA intensity; 30 Hz repetition rate; 200 ns beam width. On the downstream of the FFAG beam line, the tungsten was set at the location of the original target (80 mm diameter and 10 mm thick); the thickness was determined on the basis of previous experimental and numerical analyses for the injection of high-energy proton beams onto the tungsten target. For the proton beam configuration modeled by numerical simulations, the size of the proton beam spot was requisite experimentally and precisely, when 100 MeV protons were injected onto the tungsten target where the spallation neutrons are generated. The Gafchromic film, which is very sensitive to the charged particles, was then used to evaluate the size of the proton beam spot injected onto the tungsten target, since a graphic image on the film is acquired quickly after the irradiation of charged particles for a short time.
The reaction rates for threshold energy of high-energy neutrons and the continuous energy distribution of the spallation neutrons at the target were acquired by the foil activation method and the organic liquid scintillator, respectively. The high-energy neutrons (spallation neutrons) of threshold reactions 209Bi(
n, xn)210- xBi ( n= 3, to 12) over 15 MeV have been generated by the injection of high-energy protons over 100 MeV. Here, to obtain the reaction rates by high-energy neutrons at the target, 209Bi was selected as an activation foil (Table 1) to cover threshold energies over 15 MeV, and 115In was selected as a normalization factor (Table 2) for monitoring the spallation neutrons at the target to cover threshold reactions 115In( n, n’)115mIn over 0.32 MeV. Foil dimensions at the target were as follows: 209Bi, 50 mm in diameter and 3 mm thick; 115In, 10 × 10 × 1 mm3, and two foils were set around the target region as shown in Fig. 1(a). Additionally, nine other 115In foils (10 × 10 × 3 mm3) were placed in a circle (100 mm radius) around the target at 30 degree intervals on an acryl plate, to investigate the angular distribution of spallation neutrons as shown in Fig. 1(b). The irradiation time of 209Bi and 115In foils was four hours for measuring the neutron yield of spallation neutrons, and their reaction rates were measured by the high-purity germanium detector (ORTEC, GEM60P). As well as in the previous study,,, the detection efficiency of the germanium detector was determined by the fitting line obtained from the energy calibration with the use of several γ-ray standard sources.
The continuous energy distribution of spallation neutrons was determined by the organic liquid scintillator (Nuclear Enterprises Ltd., NE213 Scintillator; 5” in diameter and 5” long) set directly facing the tungsten target without any reactor components as shown in Fig. 2(a). The measurement circulation of the organic liquid scintillator was as indicated in Fig. 2(b). The main advantage of the measurement system is that the two signals (rise time and light output of
γ-ray and neutron) acquired coincidently with the use of specific equipment (Laboratory Equipment Corp., Dual MCA) readily provide a two-dimensional graphic image (rise time and light output). This system allows the rise time information on two signals to accomplish easy discrimination between γ-ray and neutron, and the signals can be converted into fluorescence signals of charged particles (recoil protons generated in reaction with neutrons).
The 10 mm thick tungsten target was determined on the basis of experimental and numerical analyses from the viewpoint of the full-stop of proton beams within the tungsten target. On the other hand, the proton beam irradiation experiments were carried out to monitor the size of
the proton beam spot with the use of Gafchromic film, which is highly-sensitive to charged particles. The Gafcromic film was irradiated for two min and set on the surface (downstream beam) outside the stainless steel flange (Fig. 1)
without the tungsten target, although the influence of the stainless steel flange was slightly affected to the proton beam profile. It demonstrated that the proton strength was distributed by the downstream beams (Fig. 3(a)), with a relative distribution within an approximately 40 mm diameter spot (Fig. 3(b)). The experimental result of the Gafchro-
mic film was very useful for modeling the size of the proton beam spot in numerical simulations, therefore, the irradiation experiments were important for evaluating the size of the proton beam spot for which the tungsten target (80 mm in diameter) was considered sufficient to cover.
The experimental results of the reaction rates were obtained by measuring total counts of the peak energy of
γ-ray emissions. The saturation activity D∞ [1/ s], which is proportional to the reaction rate RR[1/ s/cm3], taking into consideration of the weight M[ g] of each foil and the density ρ[ g/cm3], was calculated using the following equations:
λindicates the decay constant, Tcthe measurement counting time, Cthe counting rate, αthe internal conversion coefficient, εthe detection efficiency, Nγthe γ-ray emission rate, Tithe irradiation time, and Twthe waiting time until the start of the measurement after the irradiation. Finally, the reaction rate RRin Eq. (2) was obtained from the saturation activity D∞ in Eq. (1).
No value for the 100 MeV proton irradiation was observed in 3
n, 9 nthrough 12 nreactions of 209Bi( n, xn)210-xBi ( x= 3 to 12), since little activation was caused by insufficient irradiation with the low proton beam intensity of 30 pA and the long half-life (38.3 y) of 3 nreactions. The MCNPX calculations with JENDL/HE-2007- for nuclear data and ENDF/B-VI for cross sections of Bi were executed by a total number of 2 × 108 histories within a statistical error of 1% to obtain the reaction rates of the irradiated 209Bi foil. The spot size 40 mm diameter of proton beams was modeled in the MCNPX calculations on the basis of the experimental results in the Gafchromic film presented in Sec. 3-1. For the irradiation experiments of the 209Bi foil, a comparison (Table 3) between the experimental and the numerical values showed agreement around a relative difference of 10% in the C/E (calculation/experiment) values, excluding the 209Bi( n, 8n)202Bi reaction. Here, from the results of the 209Bi foil irradiation, the highenergy neutrons up to 50 MeV generated by 209Bi( n, xn)210-x Bi reactions were confirmed to have been bombarded by the injection of 100 MeV protons.
The measured results of the 115In reaction rates for the angular are listed in Table 4. The highest reaction rate was observed at the center of the tungsten target shown in Table 4, and the angular distribution (Fig. 4(a) and Table 4) of reaction rates appeared slightly polarized in the upper direction, when all reaction rates were normalized by that at position 210 degrees, which was the largest along the angular. The reaction rate at 210 degrees was larger than at other positions, whereas the effect of the acryl plate was considered insufficient in the measured reaction rates along the angular. Next, to investigate the effect of acryl plate on measured reaction rates, numerical simulations were executed with the use of MCNPX and JENDL/HE-2007 with (w/ ) and without (w/o) the acryl plates. A comparison of the results (Fig. 4(b)) in the presence or absence of acryl
plates showed an apparent effect on the reaction rates: the high-energy neutron flux was attenuated ahead at the target and influenced by the acryl plate. Thus, the spallation neutrons were considered significantly spherical in the angular distribution through the results in numerical simulations, although their angular distribution was observed actually reversed at the target. Subsequently the neutron yield at the target was evaluated at (9.73±0.12)×104 1/
sover 20 MeV and (1.03±0.04)×107 1/s over 0.32 MeV from the measured reaction rates of Bi and In foils, respectively.
The main characteristics of the measurements by the organic liquid scintillator are to acquire two signals of prompt (electrons) and delayed (protons, deuterons, or α-ray, etc.) fluorescence components and to discriminate
γ-ray and neutron events caused by the two signals, respectively. Thus, the discrimination between the γ-ray and the neutron was possibly caused by the difference between their fluorescence intensities in a time-dependent manner.
A comparison between the combined (
γ-ray and neutron) and the γ-ray events showed an apparent discrimination between the γ-ray and the neutron in the experimental results (Fig. 5) of fluorescence distributions. The γ-ray events were found to be considerably large in low-fluorescence distribution and difficult to discriminate the two events of γ-ray and neutron. The amount of fluorescence by highenergy neutrons was found to be small in high-fluorescence distribution because of very small counting rates in the high-channel region. Moreover, the spallation neutrons generated from the FFAG accelerator were considered to be a group of continuous energy neutrons with an ambiguity in maximum energy, since an edge of the recoil proton corresponding to the neutron energy was not found in the measurements by the organic liquid scintillator. Thus, the neutron energy calibration- of fluorescence to the light unit was conducted with the use of the results of 22Na standard source ( γ-ray energy; 1.274 MeV).
The neutron spectrum (Fig. 6) was obtained experimentally with the use of the SCINFUL-QMD code for the matrix of response functions and with the UMG code for the unfolding of experimental results (Fig. 5) together with the matrix by SCINFUL-QMD. As a reference of the neutron spectrum, the MCNPX calculations were executed with the use of ENDF/B-VI.6 for cross section data and of LA150 for the high-energy neutron and proton libraries. A comparison between the experiment (UMG) and the calculation (MCNPX) revealed an approximate reconstruction of the neutron spectrum in the experiment, ranging between 5 and 45 MeV neutron, although the discrepancy was observed in some energy regions. In the measurement system, the amount of fluorescence was insufficiently over 50 MeV neutron, because of detector sensitivity in relating the intensity and the energy of protons. Finally, the spallation neutrons up to 45 MeV were considered successfully detected by the organic liquid scintillator, since the discrimination
γ-ray and the neutron was satisfactorily conducted.
Neutron spectrum experiments on spallation neutrons were conducted in the ADS facility at KUCA to investigate the neutronic characteristics of spallation neutrons generated at the target. The reaction rates and the continuous energy distribution of spallation neutrons were measured by the foil activation method and by the organic liquid scintillator, respectively. For the reaction rate experiments of 209Bi foil, the C/E values between the experiments and the calculations (MCNPX and ENDF/B-VI) were found well within the relative difference of 10% (209Bi(
n, xn)210- xBi; x= 4 to 7), except for some reactions ( x= 8). For continuous energy distribution experiments, the spallation neutrons were observed up to 45 MeV with the use of the organic liquid scintillator and the numerical simulations (MCNPX with JENDL/HE-2007 and ENDF/B-VI.6). From these results, the neutron spectrum information (the reaction rates and the continuous energy distribution) on the spallation neutrons generated at the tungsten target was attained successfully with the combined use of the foil activation method and the organic liquid scintillator, respectively, in the injection of 100 MeV protons obtained from the FFAG accelerator onto the tungsten target.
[Table 1.] Threshold Energy, Half-life and γ-ray Energy of 209Bi(n, xn) 210-xBi Reactions (x = 3 to 12)
[Table 2.] Threshold Energy, Half-life and γ-ray Energy of 115In (n, n’)115mIn Reactions
[Fig. 1(a).] Experimental Setting of Activation Foils (Bi and In) for Measuring the Reaction Rates at the Target Position
[Fig. 1(b).] Experimental Setting of Tungsten Target for Angular Distribution of Spallation Neutrons
[Fig. 2(a).] Experimental Setting of the Organic Liquid Scintillator (NE213)
[Fig. 2(b).] Measurement Circulation of the Organic Liquid Scintillator (NE213)
[Fig. 3(a).] Measured Result (Gafchromic Film) of Proton Strength Distribution of 100 MeV Protons at Upstream Beam of the Tungsten Target
[Fig. 3(b).] Scanning Relative Result of Proton Strength Distribution in Fig. 3(a) of 100 MeV Protons at the Downstream Beam Without the Tungsten Target
[Table 3.] C/E Values between Measured and Calculated Reaction Rates of 209Bi(n, xn)210-xBi Reactions (x = 4 to 8) for 100 MeV Proton Beams
[Table 4.] Measured Results of Indium Reaction Rates for the Angular
[Fig. 4(a).] Measured Results of Angular Distribution of Indium Reaction Rates
[Fig. 4(b).] Comparison between Calculated Results of Angular Distribution of Indium Reaction Rates with and Without the Acryl Plate
[Fig. 5.] Comparison between Measurement Results of Light Output from the Organic Liquid Scintillator Before and After the Discrimination of γ-ray and Neutron
[Fig. 6.] Comparison between Measured (Unfolding) and Calculated (MCNPX) Results of Neutron Spectrum