79urn:lsid:arphahub.com:pub:D015427E-53F6-553C-9279-2E51CD684756Nuclear Energy and TechnologyNUCET2452-3038National Research Nuclear University MEPhI (Moscow Engineering Physics Institute)10.3897/nucet.4.2983729837Research ArticleReactor physicsSolution of neutron-transport multigroup equations system in subcritical systemsShamanin1Igor V.1Bedenko1Sergey V.1Nesterov1Vladimir N.1Lutsik1Igor O.1Prets1Anatoly A.1National Research Tomsk Polytechnic University, Russian Federation, Lenin Ave., Tomsk, 634050 Russian FederationNational Research Tomsk Polytechnic UniversityTomskRussia
Corresponding author: Sergey V. Bedenko (bedenko@tpu.ru)
Academic editor: Yuri Korovin
2018181020184179851F9A038C-4977-5341-A065-5507225B903C25351561809201819092018Igor V. Shamanin, Sergey V. Bedenko, Vladimir N. Nesterov, Igor O. Lutsik, Anatoly A. PretsThis is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.
An iteration method has been implemented to solve a neutron transport equation in a multigroup diffusion approximation. A thermoelectric generator containing plutonium dioxide, used as a source of thermal and electric power in spacecraft, was studied.
Neutron yield and multigroup diffusion approximation data was used to obtain a continuous and group distribution of neutron flux density spectra in a subcritical multiplying system.
Numerical multigroup approaches were employed using BNAB-78, a system of group constants, and other available evaluated nuclear data libraries (ROSFOND, BROND, BNAB, EXFOR and ENDSF).
The functions of neutron distribution in the zero iteration for the system of multigroup equations were obtained by approximating an extensive list of calculated and experimental data offered by the EXFOR and ENDSF nuclear data libraries. The required neutronic functionals were obtained by solving a neutron transport equation in a 28-group diffusion approximation. The calculated data was verified. The approach used is more efficient in terms of computational efforts (the values of the neutron flux density fractions converge in the third iteration). The implemented technique can be used in nuclear and radiation safety problems.
Subcritical systemneutron distribution functionneutron transportmultigroup diffusion approximation.Research status
Neutron transport calculation methods are used in reactor physics, for process monitoring during nuclear fuel fabrication and reprocessing, to determine the irradiated fuel burnup, in biological shielding design, and to improve the radiation measurement procedures in the nuclear material accounting and control system.
There are a lot of dedicated programs and codes that implement the Monte Carlo method based on selection of neutron interaction probabilities. The internationally recognized program packages include MMKKENO and MMKC (IPPE, Russia), MCU (Kurchatov Institute, Russia), MCNP (USA), KENO-3D (USA), MONK (Great Britain) and others. These programs make it possible to obtain highly accurate results thanks to a 3D geometry that models the object and simulation of real neutron transport in material.
High accuracy of verified Monte Carlo codes is achieved through a great deal of computational resources consumed and is defined primarily by the error with which the cross-sections of the neutron interaction with the material nuclei were determined.
Deterministic methods based on differential and finite-difference representation of the Boltzmann equation (Bell and Glesston 1974) do not require a great deal of computational resources. Such methods simplify the main simulation phase that consists in solution of the transport equation. The major deterministic methods include the discrete ordinate method, the characteristics method, the first collision probabilities method, the PN-method, and multigroup diffusion approximation as its simplest and most widespread form.
The accuracy of the solution (Ф(r,t), Σ(r,E), k_{eff}) obtained using multigroup diffusion approximation is not high, but the method in question offers a reasonably good approximation and is still used in Russian and foreign programs, e.g. in the USA’s MCNP and KENO and in many others.
An iteration process is implemented in the study to solve a system of multigroup diffusion equations which is used to determine the spectral and integral neutronic characteristics of a subcritical multiplying system.
The considered system is a thermoelectric generator containing 5.65 kg of PuO_{2}, used as a heat and electricity source in spacecraft.
Computational model
We shall write the Boltzmann equation (Bell and Glesston 1974) in a 28-group diffusion approximation:
(1)
where i is the number of the neutron energy group for which the equation has been written; k is the energy group number; D^{(}^{i}^{)} is the neutron diffusion constant for the i^{th} group; u^{(}^{i}^{)} is the neutron velocity corresponding to the average energy of the i^{th} group; F^{(}^{i}^{)}, F^{(}^{k}^{)} are the neutron flux densities in respective groups; Σ_{a}^{(}^{i}^{)} is the macroscopic neutron absorption cross-section for the i^{th} group; Σ_{R}^{i}^{→}^{k}, Σ_{R}^{k}^{→}^{i} are the macroscopic cross-sections of the neutron transfer from the i^{th} group to the underlying k^{th} group (from the overlying k^{th} group to the considered i^{th} group) respectively; ε_{fp}^{(}^{i}^{)} is the fraction of prompt neutrons in the spectrum χ_{fp} (E) falling within the i^{th} group (for all spectra ε^{(}^{i}^{)} = ʃχ_{fp} (E)dE); ν_{f}^{(}^{k}^{)} is the average number of prompt neutrons per the induced fission event for neutrons in the k^{th} group; Σ_{f}^{(}^{k}^{)} is the macroscopic fission cross-section for neutrons in the k^{th} group; ε_{fd}^{(}^{i}^{),}^{n} is the fraction of delayed neutrons in the i^{th} group; λ_{j}^{n} is the decay constant for the j^{th} group of delayed neutrons of the nuclide n; C_{j}^{n} (r,t) is the concentration of the delayed neutron precursor nuclei; S_{sf} is the integral yield of spontaneous fission neutrons; ε_{sf}^{(}^{i}^{)} is the fraction of neutrons in the spectrum χ_{sf} (E); ε_{αn}^{(}^{i}^{)} is the fraction of neutrons in the spectrum χ_{αn}(E); S_{sf} is the integral yield of spontaneous fission neutrons; S_{αn} is the integral yield of neutrons for the (α,n) reaction; and Q^{(}^{i}^{)}(r,t) is the intensity of the external neutron source.
Equation (1) needs to be supplemented with a boundary condition and the initial condition:
To organize the iteration process, it is required to define the spectra of induced fission neutrons, spontaneous fission neutrons and the (α,n)-reaction neutrons: χ_{f} (E) (ε^{(}^{i}^{)} = ∫χ(E)dE), χ_{sf} (E) and χ_{αn}(E) respectively. And the selection of the neutron distribution function and the spectrum S (E,t,r) = Q (r,t)∙χ(E) is an important simulation phase that makes it possible to improve considerably the accuracy of the solutions sought.
The normalized spectrum of induced fission neutrons χ_{f} (E) has been calculated in an assumption that the neutron yield can be approximated by Watt’s spectral function (Abagyan et al. 1981)
The parameters a and b depend on the type of the fissionable nuclide and the energy of the neutrons that have caused the fission and are clearly associated with the energy and fissile isotope averaged quantity ν^{*}_{f}. In the BNAB-78 and BNAB-93 systems of constants, a = 0.965(0.8 + 0.083∙ν^{*}_{f}) and b = 2.245/(0.8+0.083∙ ν^{*}_{f})^{2}. The parameters are used to determine the average fission neutron spectrum energy E^{*} = 1.5a + 0.25a^{2}b.
where ν^{*}_{sf} is the average number of spontaneous fission neutrons; R_{sf} is the spontaneous fission branching ratio; λ is the decay constant; N = N_{0}×exp(–λt) is the current inventory of the spontaneously fissile nuclide; and T_{1/2} is the half-life period.
The spectral distribution of spontaneous fission neutrons is
Spontaneous and induced fission neutrons are specific in being time-correlated unlike the (α,n)-neutrons and, besides, can be approximated by one and the same function
where a = Т relates to ν^{*}_{sf} through a = 0.48 + 0.2∙(1 + ν^{*}_{sf})^{1/2}. The parameter b is defined by the temperature Т and the average kinetic energy E_{f} per spontaneous fission product nucleon and is equal to b = 4E_{f}/a^{2} = 3.04/a^{2} (Vlaskin et al. 2015).
The average energy of the spontaneous fission neutron spectrum (Vlaskin et al. 2015, Wilson et al. 2002) is
The latter relation makes it possible to obtain the average energy of the spectrum χ_{sf} (E) equal to the experimental value for the most of fissile nuclides.
The yield of neutrons in the (α,n)-reaction (Bulanenko 1980, Vukolov and Chukreev 1987) is
where ω_{i} = (1 – ω_{j}) is the fraction of the α-emitter (material i); ω_{j} is the fraction of the material on which the reaction runs (material j); q^{α}_{i} is the yield of α-particles from the material of the type i; n_{j} is the concentration of the j-type material nuclei; σ^{αn}_{j} is the integral macroscopic cross-section of the (α,n)-reaction; ε_{ij} (E) = (–dE/dx) is the stopping power of the α-particle; E_{α} is the average spectrum energy for α-particles; and B_{αn} is the threshold of the (α,n)-reaction.
The data on the mass stopping power ε_{ij} (E) were calculated using the SRIM (The Stopping and Range of Ions in Matter) code (Ziegler et al. 2010).
The dependence σ^{αn}_{j} was obtained based on data in (Hansen et al. 1967, Bair and Gomez del Campo 1979, West and Sherwood 1982, Heaton et al. 1989, Murata and Shibata 2002) recommended for solving such type of problems, and the compilation and the analysis were performed in (Vlaskin et al. 2015, Wilson et al. 2002, Bulanenko 1980, Shamanin et al. 2010).
The spectral distribution of (α,n)-neutrons is
The normalized distribution of the neutron spectrum χ_{αn}(E) was obtained through approximation of experimental (EXFOR − Experimental Nuclear Reaction Data) and evaluated (ENDSF − Evaluated Nuclear Structure Data File) data (Gorshkov 1962, Herold 1968, Anderson and Neff 1969, Taherzadeh 1971, Taherzadeh and Gingo Peter 1972, Arkhipov et al. 1972, Bair and Butler 1973, Anderson 1980) by the function of the type
where the constants а and b depend on the type of the nuclei i and j.
The investigations conducted in (Gorshkov 1962, Herold 1968, Anderson and Neff 1969, Taherzadeh 1971, Taherzadeh and Gingo Peter 1972, Arkhipov et al. 1972, Bair and Butler 1973, Anderson 1980) had the purpose of obtaining highly accurate values of the neutron yields and spectra for ^{238}PuO_{2}. The measurement error does not exceed 5%, and the compilation and the evaluation were performed in (Vlaskin et al. 2015, Bulanenko 1980, Shamanin et al. 2010). Full texts and the results of the experiments are available in the EXFOR and ENDSF nuclear data libraries.
The reference data used for the spectra acquisition were the results obtained in (Taherzadeh and Gingo Peter 1972) which present accurate characteristics of the neutron spectra for the fuel composition PuO_{2} with a mass of 5.65 kg. The PuO_{2}-containing thermoelectric generator studied in (Taherzadeh and Gingo Peter 1972) is used as a heat and electricity source in space satellites. The measurements were performed to design the radiation protection for the spacecraft’s main components.
To be able to compare correctly the calculated data we have obtained with (Taherzadeh and Gingo Peter 1972) and to have the desired distribution as the result, the isotopic composition of Pu (Pu, %: 238 – 80.27; 239 – 15.87; 240 – 3.022; 241 – 0.643, 242 – 0.132) and the presence of light impurity elements (Li, Be, B, C, F, Na, Al, Si and others) were taken into account.
Therefore, the yield of neutrons due to the spontaneous fission of Pu isotopes is equal to approximately 2.48∙10^{3} 1/s/g of PuO_{2}; 1.18∙10^{4} 1/s/g of PuO_{2} thanks to the (α,n)-reaction on oxygen, and 9.21∙10^{3} 1/s/g of PuO_{2} thanks to reactions on impurities. The total yield of neutrons for the considered composition is ~2.35∙10^{4} 1/s/g of PuO_{2}.
Figs 1 and 2 show the sought-after distributions S (E) and spectra χ(E) used to solve the transport equation in a 28-group diffusion approximation.
Normalized distribution of the continuous (a) and the 28-group spectrum (b) of neutrons: 1 – spontaneous fission; 2 – (α,n)-reactions; 3 – total vector.
https://binary.pensoft.net/fig/236418
where с = 0.36; а = 0.84; and b = 3.63. The average spectrum energy is E^{*} = 2.02 MeV, and the average number of neutrons is ν^{*}_{sf} = 2.21.
The normalized distribution of the neutron spectrum for the ^{238}PuO_{2} (α,n)-reactions (see Fig. 2) is
where а = 0.78; b = 2.5. The average spectrum energy is E^{*} = 2.5 MeV.
Solution of a system of multigroup neutron diffusion equations by iteration method
In subcritical multiplying systems, the equation written in a diffusion approximation for neutrons in all energy groups (Bell and Glesston 1974) has the form (Golovatskiy et al. 2010)
DΔF(r,t) – Σ_{a}F(r,t) – ν_{f} Σ_{f} F(r,t) < 0,
and, hence, the value of the neutron flux density Ф(r,t) in such system is expected to become zero. The value F(r,t) in subcritical systems is known not to become zero but has a particular value characterized by neutron sources S (r,t) and by the neutronic parameters of the system.
With regard for the sources of neutrons, the above equation will become stationary
Therefore, equation (1) in a stationary form, with no regard for delayed neutrons and in the absence of an external source, can be written as follows (Golovatskiy et al. 2010):
(2)
The first term in (2), describing the neutron leakage, is determined from the relation
where B^{2(}^{i}^{)} is the geometrical parameter for the i^{th} group.
Organizing the iteration process requires the system of equations to be composed as follows
F^{(}^{i}^{)}_{j} = f (F^{(–1)}_{j}_{–1}, F^{(0)}_{j}_{–1}, ... , F^{(}^{k}^{)}_{j}_{–1}, ... , F^{(26)}_{j}_{–1}), k ≠ i, (4)
where j is the number of the iteration, beginning with the first one, by transforming system (2), with regard for (3), as follows
(5)
By expressing the value of the flux density in the i^{th} group from (5), we get
(6)
The system of equations is transformed as follows:
(6a)
All parameters in system (6a) are known, except the flux densities for the previous iteration ^{(}^{k}^{)}_{j}_{–1} and, therefore, the sum
This sum defines the number of the neutrons produced in the current generation during the fission of nuclei by all neutrons in the previous generation, except the neutrons in the i^{th} group.
To start the iteration process with the zero iteration, the flux density for neutrons in the i^{th} group is determined from (2), with regard for (3), by the relation
where the number of the neutrons produced in the current generation during the fission of nuclei by all neutrons in the previous generation is defined as equal to unity, that is
After the neutron flux spectrum in the zero iteration is determined using system (6), the iteration process is implemented. The calculation was performed for a subcritical spherical system of ^{238}PuO_{2} with a mass of 5.65 kg. The system of the BNAB-78 group constants (Abagyan et al. 1981) and recommendations in (Voropayev et al. 1979) were used, which made it possible to compare correctly the results obtained with the data presented in (Taherzadeh and Gingo Peter 1972).
Calculation results and discussions
The result of the calculation for the 28-group spectrum Ф(E), neutr./(cm^{2}∙s) is shown in Figs 3 and 4. We shall note that the values of the neutron flux density fractions converge in the third iteration, with the data obtained being fairly accurate (see Fig. 4 and Table 1).
Spectrum of the neutron flux density for the sphere of ^{238}PuO_{2} with a mass of 5.65 kg: 1 – spontaneous fission, continuous and 28-group; 2 – (α,n)-reactions; 3, 4 – total spectrum.
Spectrum of the neutron flux density at a distance of 50 cm from the PuO_{2} sphere: 1 – spontaneous fission (authors); 2 − (α,n)-reactions (authors); 3 – total spectrum (authors); 4 – total spectrum of PuO_{2} (Monte Carlo method (Taherzadeh and Gingo Peter 1972)).
https://binary.pensoft.net/fig/236420
Neutronic characteristics of the PuO_{2} fuel composition.
Parameter
Code
Diffusion approximation (Vlaskin et al. 2015, Voropayev et al. 1979)
ANISN (Taherzadeh and Gingo Peter 1972)
Diffusion approximation (Vlaskin et al. 2015, Voropayev et al. 1979)
SCALE-4.3, KENO-V.a (ENDF/B-V)
MCNP4b (JENDL-3.2)
Mass, kg
5.65
5.65
24.30
25.42
24.97
k_{eff}
0.39
0.35
1.0
1.0
1.0
1/(1−k_{eff})
1.64
1.55
0
0
0
It is shown in (Taherzadeh and Gingo Peter 1972) that, at a distance of 50 cm and more, the PuO_{2} fuel composition with a mass of 5.65 kg can be viewed as a point source, and the neutron flux decay takes place by the law ~1/R^{2}:
F(r,E) = F(E)×r^{2}/(R + r)^{2},
where F(E) is the neutron flux density spectrum; r is the radius of the PuO_{2} sphere; and R is the distance from the sphere surface to the detection point.
Fig. 3 presents the flux density spectrum for the neutrons produced by the PuO_{2} point source. It can be seen from the figure that curves 3 (authors) and 4 (Taherzadeh and Gingo Peter 1972) have a shift. The explanation is that curve 3 is a group distribution, and curve 4 is a continuous distribution (the spectra were obtained by Monte Carlo method (Taherzadeh and Gingo Peter 1972)). The same shift can be also observed in Fig. 2. Besides, the fuel composition in (Taherzadeh and Gingo Peter 1972) has a more complex geometry than that being modeled and includes components of the PuO_{2}-containing capsule structures. In the region of energies above 4 MeV, the spectrum χ_{sf} (E) is higher than used in (Taherzadeh and Gingo Peter 1972), and the contributor to curve 3 is the displaced group spectrum χ_{αn}(E).
A range of group constants was obtained in the iteration process organization which can be used to estimate k_{eff}. The expression defining k_{eff} has the form (Abagyan et al. 1981)
The constants were averaged based on relations (Abagyan et al. 1981)
Table 1 presents a comparison of the found values k_{eff} and the mass with those obtained using ANISN (Taherzadeh and Gingo Peter 1972), SCALE-4.3(KENO-V.a) and MCNP4b.
Conclusion
The continuous and group distributions of spectra and the system subcriticality value have been obtained using neutron yield data and a multigroup diffusion approximation. The results fairly accurately agree with the data in (Gorshkov 1962, Herold 1968, Anderson and Neff 1969, Taherzadeh 1971, Taherzadeh and Gingo Peter 1972, Arkhipov et al. 1972, Bair and Butler 1973, Anderson 1980) and with the results obtained in the codes ANISN (Taherzadeh and Gingo Peter 1972), MCNP (JENDL-3.2.) and SCALE-4.3 (KENO-V.a, ENDF/B-V).
The spectra S (E) and χ(E) used in the calculations can be calculated in dedicated codes such as SOURCE-4C (Vlaskin et al. 2015) and NEDIS-2m (Wilson et al. 2002) which allow finding the intensity and the spectrum of neutrons produced in the (α,n) reactions and during spontaneous fission and preparing a file with output data in such form as is convenient for solving the neutron transport equation using such codes as MCU, MCNP and SCALE.
The spectra S (E) and χ(E) were obtained by approximating an extensive list of calculated and experimental data offered by the EXFOR and ENDSF nuclear data libraries, and the required functionals (Ф(r,E), k_{eff}) were found as part of solving the neutron transport equation in a 28-group diffusion approximation.
The approach used is more efficient in terms of computational efforts and the nuclear data bank storage cost and can be used to solve applied nuclear and radiation safety problems.
ReferencesAbagyanLPBazazyantsNONikolayevMNTsibulyaAM (1981) AndersonME (1980) Neutron energy spectra of ^{238}Pu-^{18}O (α,n) source-unmoderated and polyethylene moderated.AndersonМENeffRA (1969) Neutron emission rates and energy spectra of two ^{238}Pu Power Sources. Nuclear Application 7: 62. https://doi.org/10.13182/NT69-A28386ArkhipovVAGorshkovGVGrebenskiiBS (1972) Neutron radiation of ^{238}PuO_{2} containing different amounts of ^{18}O.BairJButlerH (1973) Neutron Yield form a Small High Purity ^{238}PuO_{2} Source. Nuclear Technology 19: 202. https://doi.org/10.13182/NT73-A15883BairJKGomez del CampoJ (1979) Neutron Yields from Alpha-Particle Bombardment.BellDGlesstonS (1974) BulanenkoVI (1980) Neutron yield of (α,n) reaction on oxygen.GolovatskiyAVNesterovVNShamaninIV (2010) Organization of the iterative process in the numerical reconstruction of the neutron spectrum in a multiplying system with a graphite retarder.GorshkovVA (1962) Neutron Yield form the (α,n) Reaction in Be, B, C, O, F, Mg, Al, Si, and Granite Irradiation with Polonium a-particles.HansenLFAndersonJDMcClureJW (1967) The (α,n) cross section on ^{17}O and ^{18}O between 5 and 12,5 MeV.HeatonRLeeHSkensvedPRobertsonBC (1989) Neutron Production from Thick-Target (α,n) Reactions.HeroldTR (1968) Neutron spectrum of ^{238}PuO_{2}.MurataTShibataK (2002) Evaluation of The (α,n) Reaction Nuclear Data for Light Nuclei.ShamaninIVBulanenkoVIBedenkoSV (2010) Neutron Radiation Field of the Irradiated Ceramic Nuclear Fuel of Different Types. Izvestiya vuzov.TaherzadehM (1971) Neutron Yield form the (α,n) Reaction in the Isotope ^{18}O.TaherzadehMGingo PeterJ (1972) Neutron radiation characteristics of plutonium dioxide fuel.VlaskinGNKhomyakovYSBulanenkoVI (2015) Neutron Yield of the Reaction (α,n) on Thick Targets Comprised of Light Elements.VoropayevAIVan’kovAAVozyakovVV (1979) Group neutron cross sections for fission and radiation capture of transactinides. VANT. Ser.VukolovVAChukreevFE (1987) Neutron yield for chemical compounds of actinides. Soviet Atomic Energy.62(4): 271–276.WestDSherwoodAC (1982) Measurements of Thick-Target (α,n) Yields from Light Elements. Annals of Nuclear Energy. 1982, v. 9, 551–577.WilsonWBPerryRTCharltonWS (2002) SOURCES 4С: A Сode for Calculating (α,n), Spontaneous Fission, and Delayed Neutron Sources and Spectra. Los Alamos National Laboratory report LA UR-02-1839.ZieglerJFZieglerMDBiersackJP (2010) SRIM – The Stopping and Ranges of Ions in Matter.
* Russian text published: Izvestiya vuzov. Yadernaya Energetika (ISSN 0204-3327), 2017, n.4, pp. 38-49.