Next Article in Journal
Relativistic Configuration-Interaction and Perturbation Theory Calculations of the Sn XV Emission Spectrum
Previous Article in Journal
Relativistic B-Spline R-Matrix Calculations for Electron Scattering from Thallium Atoms
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analyzing the Neutron Parameters in the Accelerator Driven Subcritical Reactor Using the Mixture of Molten Pb-Bi as Both Target and Coolant

Institute of Applied Technology, Thu Dau Mot University, Thu Dau Mot City 820000, Vietnam
Submission received: 7 October 2021 / Revised: 18 November 2021 / Accepted: 19 November 2021 / Published: 21 November 2021

Abstract

:
In this paper, the Accelerator Driven Subcritical Reactor (ADSR) was simulated based on the structure of the TRIGA-Mark II reactor by the MCNPX program. The proton beam interacts on the Pb-Bi molten target with various energy levels from 0.5 GeV to 2.0 GeV. The important neutron parameters to evaluate the operability of ADSR were calculated as: the neutron yields according to various thicknesses of the target and according to the energy of the incident proton beam; the effective neutron multiplication factor for various fuel mixtures, along with its stability for some fuel mixtures; the axial and radial distributions of the neutron flux along with the height and radius of the core. The obtained results had shown a good agreement in using Pb-Bi molten as the interaction target and coolant for ADSR.

1. Introduction

The Accelerator-Driven Subcritical Reactor (ADSR) is a great research interest in the world since it was suggested back by many authors, such as C. Rubbia; K. Furukawa, C. D. Bowman et al. [1,2,3]. It can generate energy and transmute radioactive wastes cleanly and safely. The recent researches about ADSR concentrate on neutron source, designing target, and fuel.
Regarding interactive targets, the most popular studies are using solid targets to generate neutrons. In 1999, the group of authors X. Ledoux, F. Borne, A. Boudard, et al. calculated the energy spectrum of neutrons emitted at different angles when a proton beam interacts on the lead target with energies of 0.8 MeV, 1.2 MeV, and 1.6 MeV respectively [4]. Also in 1999, the group of authors S. Meigo et al. calculated the distribution of the neutron flux from a thick lead target with incident proton energies of 0.5 GeV and 1.5 GeV by the MCNP4A program [5].
Commonly used coolants are light water, molten salt, or molten Pb-Bi mixture. The group of authors Hasanzadeh et al. used the MCNP program to perform simulation of Moroccan TRIGA-Mark II Reactor and calculate Neutrons Kinetic and Parameters [6,7]. In these studies, the coolant was light water. Also from this model, Abdelaziz Darif et al. [8] calculated the dependence of the effective neutron multiplication factor (keff) on the number of reflectors in the core of the reactor. The results showed that keff can be increased without changing the size of the core. These authors also analyze neutron parameters using DRAGON.v5 and TRIVAC.v5 codes [9]. Ned Xoubi researched the neutron source for Jordan’s subcritical reactor [10]. In this studying, a proton is coupled to a subcritical assembly with a lead target; the spallation neutron source is investigated; 3-D Monte Carlo simulating of radiation transport is performed for the system. These results had shown an increase of seven orders of magnitude compared to the current JSA Pu-Be source-driven core. Xue-Chao Zhao et al. conducted an analysis of Th-U breeding capability for an accelerator-driven subcritical molten salt reactor [11]. These results had shown the subcriticality of the core has a considerable effect on the Th-U breeding; a high subcriticality is favorable to improving the conversion ratio, increasing the net U-235 production, and reducing the doubling time.
There have been many proposals and calculations using liquid lead as an interactive target and coolant [12,13,14,15,16,17].
In general, there are not many studies on neutron parameters for ADSR using molten Pb-Bi as a target and coolant. In this paper, ADSR was simulated based on TRIGA-Mark II reactor structure [6,7,8] using the MCNPX program. From this structure, the neutron parameters were calculated, including:
  • The neutron yields on molten Pb-Bi target;
  • The effective neutron multiplication factors keff;
  • The distributions of the neutron fluxes along the height of ADSR with molten Pb-Bi target, U-Th mixed fuels;
  • The distributions of the neutron fluxes along a radial of ADSR with molten Pb-Bi target, U-Th mixture fuels.

2. Model and Calculations

The TRIGA-Mark II subcritical reactor model was simulated as the Figure 1. Some characteristics of the core are shown in Table 1. The structure of the core consists of 108 fuel rods. The core is divided into 7 rings. All of these are assumed to be placed in the molten lead-bismuth medium.
From this structure, the neutron parameters were calculated. The results were presented below.

2.1. The Neutron Yields (Yn/p) from the Molten Pb-Bi Target

The proton beam from the accelerator interacts with the molten Pb-Bi target. From there, the neutron yields were calculated. The MCNPX program was used in this case. The Bertini model in MCNPX was chosen for the calculation because it is the most suitable model according to previous studies [18,19]. The interacting target was molten Pb-Bi consisting of isotopes Pb-204, Pb-206, Pb-207, Pb-208 with ratios of 1.4%, 24.1%, 22.1%, 52.4%, respectively. The molten Pb-Bi mixture consists of 44.5% lead and 55.5% bismuth.
  • The neutron yields according to target thicknesses
In this section, the incident proton with 1 GeV energy interacts with the molten Pb-Bi target. From here, the neutron yields were calculated with different target thicknesses from 10 cm to 50 cm. These results were compared with previous works [20,21]. The authors Giao N. M. et al. [20] calculated the neutron yield by using data from the Japanese library JENDL-HE. Mariono et al. [21] performed calculations based on the Bertini model. These results were presented in Table 2 and Figure 2.
These results showed that the neutron yields on the molten Pb-Bi target by thickness were similar to that for molten lead. When the target thickness was 10 cm, the medium neutron yield was about 10 n/p. As the thickness increases, Yn/p also increases; at the 50 cm thickness, Yn/p from 25 n/p to 27 n/p. Yn/p for solid Pb-Bi was lower than that calculated by the authors N. M. Giao [20] et al. with thicknesses from 10 cm to 27 cm but higher from 27 cm to 50 cm. The difference was highest at 10 cm thickness with 10.8% and lowest at 30 cm thick with 2.7%. Yn/p for liquid Pb-Bi compared with the study of Maiorino et al. [19] on the solid lead target is always higher at all thicknesses. The highest difference was at the thickness of 50 cm with 9.5% and the lowest at the thickness of 10 cm with 8.3%. With these results, the use of Pb-Bi-liquid as a neutron-generating target in ADSR is completely appropriate. Based on these data, the designing of the ADSR core can be conducted to the requirements of each system.
  • The neutron yields according to the levels of energy of the incident proton beam
In these calculations, the proton beam had energies from 0.5 GeV to 2.0 GeV. The length of the target was given to be equivalent to the height of the reactor core, 35 cm. These results were presented in Table 3 and Figure 3.
The results have shown that the neutron yield from molten Pb-Bi target increases with the incident proton beam intensity, similar to the case of solid lead and liquid lead targets. However, Yn/p for the molten Pb-Bi was always higher than all with energies of the incident proton beam, from 9.2% to 14.2%. The highest difference is when compared with the case of proton beam intensity of 0.5 GeV, on the liquid lead target, with 14.2%. The lowest difference was when compared with the case of proton beam intensity of 2.0 GeV, on the liquid lead target, with 9.2%. These results further confirmed the possibility of using molten Pb-Bi for ADSR.

2.2. The Effective Neutron Multiplication Factors keff

For ADSR can operate and generate positive energy, the effective neutron multiplication factor keff must be greater than 0.68 [21]. In this section, the neutron was investigated for ADSR using the molten Pb-Bi and compared with the molten Pb case. The fuels used were a mixture of Th-UO with different component ratios. The stability of keff in each case was also investigated.
  • Calculating keff factor
The investigated Th-UO fuel mixture consists of 8 cases, numbered from F1 to F8, according to the increasing proportion of uranium and decreasing of thorium. Details of fuel composition ratio and the results of medium keff calculation were presented in Table 4.
The results have shown that, as the uranium ratio increased, the thorium ratio decreased, the keff increased. With F1 fuel, the keff for both Pb and Pb-Bi was low, less than 0.5. From F2 fuel, keff reaches a value higher than 0.68. The keff for molten Pb and Pb-Bi differs on average by 0.01 with the fuel mixtures. However, with the F8 fuel mixture, the average keff was greater than 1.0 for Pb, but for Pb-Bi it is still less than 1.0. These should be considered carefully because, with Pb, the reactor has reached the critical state, while with Pb-Bi, it is still in the subcritical state.
  • Stability of keff in some fuel mixtures
From the results above, some fuel mixtures were selected to evaluate the stability of the keff in some cycles. The selected mixtures were F2, F3, F5, and F7. With these mixtures, the ADSR generates positive energy and operates in a subcritical state. Critical calculations skipped the first 100 cycles, selects 16 keff values of cycles from 100 to 250 to evaluate. The results were shown in Figure 4.
These results have shown that the values of keff were relatively stable throughout all cycles. For the case of the F2 fuel, the medium keff is 0.708338 then the range is from 0.69541 to 0.72226. With the F3 fuel, the medium keff is 0.902307, the values vary from 0.88528 to 0.91899. In the case of the F5 fuel, the value of keff is 0.902307, ranging from 0.94279 to 0.97565. With the F7 fuel, the mean keff is 0.979434 and ranges from 0.9621 to 1.00029. In particular, for the case of the F7 fuel, the maximum of keff is sometimes greater than 1.0. Therefore, careful studies and calculations are needed if using this fuel at this ratio.

2.3. The Radial and Axial Distributions of the Neutron Fluxes

In this section, the neutron fluxes were calculated using the F2, F3, F5, F7 fuel mixtures. The MCNP5 program and F4mesh tally were used. The neutron fluxes were calculated from the F4mesh tally in MCNPX based on Equation (1):
F 4 = C Φ ( E ) R ( E ) d E
with normalization for the current of proton beam I (mA) and neutron yield Yn/p. Where C is a multiplication, (ncm−2s−1) is the neutron flux and R(E) is any combination of sums and products of energy-dependent quantities known to MCNPX [22]. The neutron fluxes have been determined according to Equation (2) [7]:
Φ = 2 × 1 0 3 C / s m A × 1 p 1.6 × 1 0 19 C × F 4 × Y n / p
with Yn/p (n/p) is the neutron yield.
  • The axial distributions of the neutron fluxes
The neutron fluxes were calculated at 20 locations along with the height of the fuel rod inside the core. The results were shown in Figure 5.
These results have shown that the distributions of the neutron fluxes with different fuel cases have the same form. The neutron flux was maximum value in the center of the core and decreases towards the sides. As the proportion of uranium in the mixture increased, as the proportion of thorium decreased, the maximum neutron flux increased. However, the exact maximum locations of the neutron flux were also varying. With the F2 fuel, the maximum flux was 5.15903 × 1014 (ncm−2s−1) at position h = 23.438 cm; with the F3 fuel, the maximum flux was 1.25981 × 1015 (ncm−2s−1) at position h = 21.562 cm. Particularly for the two cases F5 and F7, the maximum flux position was the same, at the same position 17.812 cm, with the maximum flux being 2.62486 × 1015 (ncm−2s−1) and 5.6457 × 1015 (ncm−2s−1), respectively. These results had shown that the greater the percentage of uranium, the closer the position of the maximum neutron flux was to the center of the core.
  • The radial distributions of the neutron fluxes
The neutron fluxes were calculated at 20 locations along the radius of the core. These results were shown in Figure 6.
These results also showed that the maximum neutron flux at the center of the core and gradually decreases to the outside. As the proportion of uranium in the mixture increased, the maximum neutron flux increased. With the F7 fuel mixture, the maximum neutron flux reached more than 8 × 1015 (ncm−2s−1). The neutron fluxes decreased rapidly from the center of the core to a radius of 8 cm, then it began to decrease slowly, and from 8 cm to 12 cm, the neutron fluxes were almost unchanged. From 12 cm, the neutron fluxes continued to decrease but slower than at the beginning.
These results showed that the neutron fluxes for the molten Pb-Bi target using the mixed fuel cases F2, F3, F5, F7 were very high and suitable for the operation of ADSR. Depending on the actual conditions, requirements, and operational purposes of each ADSR, it is necessary to choose the appropriate fuel component ratio.

3. Conclusions

In this paper, the ADSR was simulated based on the structure of the TRIGA-Mark II reactor. From this structure, the neutron parameters have been calculated by MCNPX and MCNP5 programs. Here, the molten Pb-Bi mixtures were used as both a target and coolant. The results included the neutron yield on molten lead and molten Pb-Bi targets; the effective neutron multiplication factor (keff) with different ratios of fuel mixtures containing thorium and uranium; evaluate the stability of the keff in several fuel mixtures over several cycles; calculated the radial and axial distribution of the neutron fluxes. The results have shown that Molten Pb-Bi was very suitable to be used both as an interactive target and as a coolant in ADSR. There are only a few remarks about the proportions of uranium in the mixture that needs to be carefully calculated so that the ADSR generates positive energy and remains subcritical. These results contribute more data for the researching and designing of ADSR most effectively and economically.

Funding

This research received no external funding.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

Not applicable.

Conflicts of Interest

The author declare that there is no conflict of interests regarding the publication of this paper.

References

  1. Roche, C.; Revol, J.P.; Mandrillon, P.; Klapisch, R.; Kadi, Y.; Gelès, C.; Galvez, J.; Fiétier, N.; Carminati, F.; Buono, S.; et al. Conceptual Design of a Fast Neutron Operated High Power Energy Amplifier (No. CERN-AT-95-44-ET); Los Alamos National Laboratory: Los Alamos, NM, USA, 2005. [Google Scholar]
  2. Furukawa, K.; Kato, Y.; Ohmichi, T.; Ohno, H. Combined system of accelerator molten-salt breeder (AMSB) apd molten-salt converter reactor (MSCR). At. Tekhnika Rub. 1983, 15, 23–29. [Google Scholar]
  3. Bowman, C.; Arthur, E.; Lisowski, P.; Lawrence, G.; Jensen, R.; Anderson, J.; Blind, B.; Cappiello, M.; Davidson, J.; England, T.; et al. Nuclear energy generation and waste transmutation using an accelerator-driven intense thermal neutron source. Nucl. Instrum. Methods Phys. Res. Sect. A Accel. Spectrometers Detect. Assoc. Equip. 1992, 320, 336–367. [Google Scholar] [CrossRef]
  4. LeDoux, X.; Borne, F.; Boudard, A.; Brochard, F.; Crespin, S.; Drake, D.; Duchazeaubeneix, J.C.; Durand, D.; Durand, J.M.; Fréhaut, J.; et al. Spallation Neutron Production by 0.8, 1.2, and 1.6 GeV Protons on Pb Targets. Phys. Rev. Lett. 1999, 82, 4412. [Google Scholar] [CrossRef]
  5. Meigo, S.; Takada, H.; Chiba, S.; Nakamoto, T.; Ishibashi, K.; Matsufuji, N.; Maehata, K.; Shigyo, N.; Watanabe, Y.; Numajiri, M. Measurements of neutron spectra produced from a thick lead target bombarded with 0.5- and 1.5-GeV protons. Nucl. Instrum. Methods Phys. Res. Sect. A Accel. Spectrometers Detect. Assoc. Equip. 1999, 431, 521–530. [Google Scholar] [CrossRef]
  6. Hassanzadeh, M. Optimization and analysis of the effects of physical parameters in a TRIGA-ADSR. Kerntechnik 2017, 82, 334–338. [Google Scholar] [CrossRef]
  7. Hassanzadeh, M.; Feghhi, S.A.H. Sensitivity analysis of core neutronic parameters in accelerator driven subcritical reactors. Ann. Nucl. Energy 2014, 63, 228–232. [Google Scholar] [CrossRef]
  8. Darif, A.; Chetaine, A.; Mghar, M.; Kabach, O. Investigation of the influence of various reflector elements material types on the core nuclear parameters of the 2 MW Triga Mark-II research reactor using the Monte Carlo code. Turk. J. Phys. 2018, 42, 41–46. [Google Scholar] [CrossRef]
  9. Darif, A.; Chetaine, A.; Kabach, O.; Mghar, M.; Saidi, A. Neutronic Analysis of Moroccan TRIGA MARK-II Research Reactor using the DRAGON. v5 and TRIVAC. v5 codes. Phys. AUC 2017, 27, 41–49. [Google Scholar]
  10. Xoubi, N. Neutronic design study of accelerator driven system (ADS) for Jordan subcritical reactor as a neutron source for nuclear research. Appl. Radiat. Isot. 2018, 131, 71–76. [Google Scholar] [CrossRef] [PubMed]
  11. Zhao, X.-C.; Cui, D.-Y.; Cai, X.-Z.; Chen, J.-G. Analysis of Th-U breeding capability for an accelerator-driven subcritical molten salt reactor. Nucl. Sci. Tech. 2018, 29, 121. [Google Scholar] [CrossRef]
  12. Tien, T.M. Distributions of neutron flux from (p, n) reaction on the liquid lead target for accelerator driven subcritical reactor (ADSR). J. Phys. Conf. Ser. 2019, 1324, 012061. [Google Scholar] [CrossRef]
  13. Tien, T.M. Distribution of Neutrons from the Reaction (p, n) on the Liquid Lead Target in the Accelerator Driven System Reactor. J. Phys. Conf. Ser. 2019, 1172, 012066. [Google Scholar] [CrossRef]
  14. Tien, T.M.; Dung, T.Q. Calculation of the neutron flux distribution in the accelerator driven subcritical reactor with (Th-233U) O2 and (Th-235U) O2 mix fuel. J. Phys. Conf. Ser. 2020, 1451, 012009. [Google Scholar] [CrossRef]
  15. Tien, T.M.; Khanh, N.K.; Hien, B.T.T.; Luong, N.T.T.; Phung, N.H.T.; Thi, N.T.M. K effective factor in the ADSR using liquid lead target and (Th233U) O2, (Th235UO)2, (Th238U) O2 fuel mixture. J. Phys. Conf. Ser. 2020, 1706, 012009. [Google Scholar] [CrossRef]
  16. Tien, T.M.; Khanh, N.K.; Ngan, N.K.; Nhi, N.T.T. Radioactive decay of thorium and uranium in the liquid lead and molten salt. IOP Conf. Ser. Mater. Sci. Eng. 2021, 1070, 012024. [Google Scholar] [CrossRef]
  17. Tien, T.M.; Phung, N.H.T.; Hien, B.T.T. Effect of reflector materials to the neutron flux and k effective in the accelerator driven subcritical reactor. IOP Conf. Ser. Mater. Sci. Eng. 2021, 1070, 012025. [Google Scholar] [CrossRef]
  18. İskender, D.; Arasoğlu, A.; Tel, E. Neutron multiplicity with 1.0 and 1.2 GeV proton-induced spallation reactions on thin targets. Chin. J. Phys. 2008, 46, 124–134. [Google Scholar]
  19. Maiorino; Jose, R.; Mongelli, S.T.; Santos, A.; Anefalos, S.; Deppman, A.; Carluccio, T. A review of models and codes for neutron source (spallation) calculation for ADS application. In Proceedings of the International Nuclear Atlantic Conference-INAC 2005, Santos, Brazil, 28 August–2 September 2005; p. 85-99141-01-5. [Google Scholar]
  20. Nguyen; Giao, M.; Le Thi, T.T.; Nguyen, T.A.T. Screening Effect in (pn) Reactions on Heavy Element Targets 82 206 Pb, 92 238 U, 74 184 W, 79 197 Au. No. IC-2010/057. Abdus Salam Int. Cent. Theor. Phys. 2010, 43, 43048729. [Google Scholar]
  21. Nifenecker, H.; David, S.; Loiseaux, J.M.; Meplan, O. Basics of accelerator driven subcritical reactors. Nucl. Instrum. Methods Phys. Res. Sect. A Accel. Spectrometers Detect. Assoc. Equip. 2001, 463, 428–467. [Google Scholar] [CrossRef]
  22. Pelowitz, Desise B. “MCNPX User’s Manual, Version 2.5.0, April 2005.” LACP05-0369, 2005. Available online: https://www.osti.gov/etdeweb/servlets/purl/602571 (accessed on 21 August 2021).
Figure 1. Structure of the TRIGA-Mark II core by MCNP code.
Figure 1. Structure of the TRIGA-Mark II core by MCNP code.
Atoms 09 00095 g001
Figure 2. The neutron yields from molten Pb-Bi target with different thicknesses from 10 cm to 50 cm.
Figure 2. The neutron yields from molten Pb-Bi target with different thicknesses from 10 cm to 50 cm.
Atoms 09 00095 g002
Figure 3. The neutron yields from the proton beam with the levels of energy from 0.5 GeV to 2.0 GeV.
Figure 3. The neutron yields from the proton beam with the levels of energy from 0.5 GeV to 2.0 GeV.
Atoms 09 00095 g003
Figure 4. The keff from 100 to 250 cycles.
Figure 4. The keff from 100 to 250 cycles.
Atoms 09 00095 g004
Figure 5. The axial distributions of the neutron fluxes with molten Pb-Bi target.
Figure 5. The axial distributions of the neutron fluxes with molten Pb-Bi target.
Atoms 09 00095 g005
Figure 6. The radial distributions of the neutron flux with molten Pb-Bi target.
Figure 6. The radial distributions of the neutron flux with molten Pb-Bi target.
Atoms 09 00095 g006
Table 1. Structure of the core of the TRIGA- Mark II reactor.
Table 1. Structure of the core of the TRIGA- Mark II reactor.
CharacteristicsDetails
Fuels
Density of the fuel
Coolant
Reflectors
Density of the reflectors
Control rods
Diameter of the fuels
Height of the fuels
Diameter of the core
Height of the core
UThO
10.32 g cm−3
Molten lead-bismuth
Graphite
2.25 g cm−3
B4C
3.73 cm
38.1 cm
56 cm
72 cm
Table 2. The neutron yields from molten Pb-Bi target with different thicknesses from 10 cm to 50 cm.
Table 2. The neutron yields from molten Pb-Bi target with different thicknesses from 10 cm to 50 cm.
Thickness (cm)Neutron Yield (n/p)
Pb-Solid [14]Pb-Solid [15]Pb-Bi-Liquid (This Work)
1010.9689.0809.899
2018.85216.08717.616
3022.21520.85922.839
4023.92623.48025.862
5024.85624.95427.561
Table 3. The neutron yields from the proton beam with the levels of energy from 0.5 GeV to 2.0 GeV.
Table 3. The neutron yields from the proton beam with the levels of energy from 0.5 GeV to 2.0 GeV.
Proton Energy (GeV)Neutron Yield (n/p)
Pb(Solid)Pb (Liquid)Pb-Bi (Liquid)
0.58.35408.53109.7400
0.816.89617.03519.193
1.021.95221.96824.526
1.532.80432.75836.569
2.042.59542.19246.909
Table 4. Details of the ratio of components in the fuel mixture and the results of keff.
Table 4. Details of the ratio of components in the fuel mixture and the results of keff.
Fuel T 90 232 h U 92 233 O 8 16 keff (Pb)keff (Pb-Bi)
F10.8000.1000.1000.453450.430634
F20.7000.2000.1000.715460.708338
F30.6100.2900.1000.907190.902307
F40.6000.3000.1000.927640.919357
F50.5800.3200.1000.961230.959649
F60.5700.3300.1000.979840.974404
F70.5650.3350.1000.987510.979434
F80.5600.3400.1001.000440.983260
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tran Minh, T. Analyzing the Neutron Parameters in the Accelerator Driven Subcritical Reactor Using the Mixture of Molten Pb-Bi as Both Target and Coolant. Atoms 2021, 9, 95. https://0-doi-org.brum.beds.ac.uk/10.3390/atoms9040095

AMA Style

Tran Minh T. Analyzing the Neutron Parameters in the Accelerator Driven Subcritical Reactor Using the Mixture of Molten Pb-Bi as Both Target and Coolant. Atoms. 2021; 9(4):95. https://0-doi-org.brum.beds.ac.uk/10.3390/atoms9040095

Chicago/Turabian Style

Tran Minh, Tien. 2021. "Analyzing the Neutron Parameters in the Accelerator Driven Subcritical Reactor Using the Mixture of Molten Pb-Bi as Both Target and Coolant" Atoms 9, no. 4: 95. https://0-doi-org.brum.beds.ac.uk/10.3390/atoms9040095

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop