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.7.7415674156Research ArticleReactor physicsOn the determination of neutron multiplication by the Rossi-alpha methodGrabezhnoyVladimir A.1DulinViktor A.1DulinVitaliy V.1vitdulin@yandex.ruMikhailovGennady M.1JSC SSC RF-IPPE n.a. A.I. Leypunsky, 1 Bondarenko Sq., 249033 Obninsk, Kaluga Reg., RussiaJSC SSC RF-IPPE n.a. A.I. LeypunskyObninskRussia
Corresponding author: Vitaliy V. Dulin (dulin@ippe.ru)
Academic editor: Georgy Tikhomirov
202123092021732532578186FFC3-3532-5A0B-92EC-FD7C04C29E1A0912202001072021Vladimir A. Grabezhnoy, Viktor A. Dulin, Vitaliy V. Dulin, Gennady M. MikhailovThis 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.
Introduction. This work contains the results of determining the prompt neutron multiplication factor in the subcritical state of a one-core BFS facility, obtained by the neutron coincidence method, for which the influence of the error in the β_{eff} in determining the multiplication factor turned out to be insignificant. The core of the facility consisted of rods filled with pellets of metallic depleted uranium, 37% enriched uranium dioxide and 95% enriched plutonium, sodium, stainless steel and Al_{2}O_{3}. Stainless steel served as a reflector.
Methods. In contrast to the inverse kinetics equation solving (IKES) method, which is convenient for determining reactor subcritical states, the neutron coincidence method practically does not depend on the error in the value of the effective fraction of delayed neutrons β_{eff}. If in the IKES method the reactivity value is obtained in fractions of β_{eff}, i.e., from the measurement of delayed neutrons, the neutron coincidence method is based on the direct measurement of the value (1 – k_{σ}_{p})^{2}, where is the effective multiplication factor by prompt neutrons. The total multiplication factor is defined as k_{eff} = k_{σ}_{p} + β_{eff}. If, for example, k_{eff} ≈ 0.9 (which is typical for determining the fuel burnup campaign), then it is the error in determining k_{σ}_{p} that is the main one in comparison with the error in β_{eff}. Thus, a 10% error in β_{eff} of 0.003–0.004 (typical for plutonium breeders) will make a contribution to the error 1 – k_{eff} equal to 1 – k_{σ}_{p} + β_{eff} ≈ 0.00035, i.e., approximately 0.35%, but not 10%, as in the IKES method.
Rossi-alpha measurements were carried out using two ^{3}He counters and a time analyzer. The measurement channel width Δt was 1.0 μs. From these measurements, the value of the prompt neutron multiplication factor was obtained. In this case, the space-isotope correlation factor for the medium with a source was calculated using the following values: Φ(x) – solutions of the inhomogeneous equation for the neutron flux and Φ^{+}(x) – solutions of the ajoint inhomogeneous equation.
Results. The authors also present a comparison of the results of the Rossi-alpha experiment and measurements of the BFS-73 subcritical facility by the standard IKES method in determining the multiplication factor value. The data of the IKES method differ insignificantly from the results of the Rossi-alpha method over the entire range of changes in the subcriticality with an increase in the subcriticality of the BFS-73 one-core facility.
Conclusion. It was impossible to apply the neutron coincidence method to fast reactors; however, the method turned out to be quite workable on their models created at the BFS facility, which was successfully demonstrated in this study.
The multiplication factor, estimated by the standard methods adopted at the BFS facility, i.e., the IKES method (Kazansky and Matusevich 1994) at shallow subcriticality (~ β_{eff}), and the inverse multiplication method, taking into account the change in the neutron registration efficiency at deeper subcriticality (~ 10 β_{eff}), is directly related to the accuracy of measurements of the β_{eff} value itself. Previously implemented international programs for measuring β_{eff} (see, for example, (Dulin VA and Dulin VV 1998)) in fairly simple single-zone cores with dense reflectors made of uranium dioxide made it possible to determine the error in β_{eff} ≈ 3–4%. Nevertheless, it is of interest to determine deep subcriticalities by a fundamentally different method, without using β_{eff}.
The work contains the results of determining the multiplication factor in the subcritical state, obtained by the neutron coincidence method (Dulin and Matvienko 2002), for which the influence of the error in β_{eff} in determining the multiplication factor is insignificant (tens of times less than in methods where the value of reactivity is inversely proportional to the value of β_{eff}, as in the IKES method (Kazansky and Matusevich 1994)).
BFS facility measurements
The effective neutron multiplication factor in the subcritical state of a single-core BFS facility was determined using both of the above methods. The core of the facility consisted of rods filled with pellets of metallic depleted uranium, 37% enriched uranium dioxide and 95% enriched plutonium, sodium, stainless steel and Al_{2}O_{3}. Stainless steel served as a reflector. The height of the core was 45 cm, the equivalent diameter ≈ 65 cm. The subcritical state was achieved by removing from the core 11 fuel rods, of which the controls and CPS were composed.
Measurements by the neutron coincidence method were carried out using two SNM-18 ^{3}He counters placed in the central channel (instead of a fuel rod) and a time analyzer (TA). The helium pressure in the counters was about two atmospheres. One of the counters was located 1 cm above the central plane of the critical facility core, and the other was below it. The length of the counter was approximately 25 cm. The first was used as a starting one for the TA, and the second was used as a counting one. The width of the VA channel Δt = 1 μs. The delay T in the counting channel was measured and turned out to be 40.1 μs.
As an example, Figure 1 presents the result of one of such measurements and Figure 2 shows the temporal distribution of its correlated part obtained after the constant background was subtracted (determined in the region of 125–250 channels). Along the abscissa, the TA channels (1 μs wide) are located; the ordinate shows the count in them.
In (Mikhailov and Dulin 1995), the analysis method is described in detail, which makes it possible, using the measured correlated part of the distributions (see Fig. 2), the average constant background N̄ (see Fig. 1, to the right of the 125^{th} channel), the integral of readings above zero in the right of the correlated part of the distribution (see Fig. 2) J_{right} and the source of spontaneous fissions Q_{sp}, to obtain the value of the prompt neutron multiplication factor (1 – k_{σ}_{p}_{exp})^{2}:
where Φ(x) is the solution to the inhomogeneous equation for the neutron flux; Φ^{+}(x) is the fission spectrum-averaged solution to the ajoint inhomog eneous equation, which is the probability for the detector to register a neutron that appears at the point x (Marchuk and Orlov 1961). The numbers in the denominator Ξ are the ν values for fissions of ^{235}U, ^{239}Pu, ^{238}U and the source of spontaneous fissions of ^{240}Pu , and those in the numerator Θ are the 〈ν(ν – 1)〉 values for themselves, respectively.
The processing of Rossi-α temporal distributions was carried out taking into account the asymmetry of the temporal distribution of the background relative to the center of the correlated signal, which was at 40.1 µs of the TA scale. For this, the signal of a periodic pulse generator with a duration of 0.05 µs was applied simultaneously to the starting and counting (through the delay line T) TA channels with a channel width of 0.1 µs, and the count was recorded in the 401^{st} channel, which corresponded to the delay value T = 40.1 µs.
This type of background asymmetry was found experimentally on a fast critical assembly (Tsutomu 1968), and a quantitative description of the time dependence of the background N (t) was given in (Mikhailov and Dulin 1995). Having carried out measurements with exact knowledge of the delay T, but not assuming, as was done in (Grabezhnoy et al. 2015), that it is determined by the position of the peak maximum visible on the screen, we compared the results of both approaches and experimentally verified the theoretical conclusions (Mikhailov and Dulin 1995), according to which:
Nleft(t)=N1-ACαα+C2αα-C·eC(t-T)-α+Cα-C·eα(t-T); t < T, (2)
Nright(t)=N1-ACαα+C·e-α(t-T); t ≥ T. (3)
Note that the temporal distribution of the background is continuous, i.e., N_{left}(T) = N_{right}(T).
The correlated part of the distribution at t ≥ T obtained after the constant mean background is subtracted, equal to [ACα/(α + C)]e^{−α(}^{t} – ^{T}^{)} (the right side visible in the experiment), according to (Mikhailov and Dulin 1995),has an amplitude in α/(α + C) times less than the theoretical value. Adding the ‘lost’ part of the background (3), we obtain the true distribution of the right side:
R_{cor}(t) = AC_{count}Ce^{−α(}^{t} – ^{T}^{)}; t ≥ T. (4)
Correspondingly, the area under the right correlated component (integral over t ≥ T) is
J_{right} = AC_{count}C/α. (5)
To obtain the true measured correlated areas J_{cor} from those ‘visible’ above the constant background level N, it is necessary to additionally measure the speed of incoming counts C at the start and then find the α value from the experiment. According to (7), it is on the right side of the equality R_{cor}^{exp}(t) (see Figs 2, 3), since the theoretical temporal behavior of the correlated component and the background ‘under it’ on the right have the same exponential dependence. Figure 3 shows that the decay time of the neutron density τ is approximately 8 μs.
Time distribution of the correlated part obtained after the constant background is subtracted N̄ (Axis Y).
https://binary.pensoft.net/fig/592156
Knowing C and α, we can calculate the value of the addition to the area AC/α to the right of t ≥ T:
ΔJ_{right} = (AC/α)×[C/(α + C)]. (6)
When calculating ΔJ_{left}(t), the integration, according to (Mikhailov and Dulin 1995), is taken for all t < T. Then the addition to the left side of the square will be
ΔJleft=-ACαC(α+C)2α2(α-C)N~-α+C(α-C)α (7)
Table 1 shows the results of the experiment (the first four columns), Table 2 uses the account of experimental additions in formulas (6), (7) and the values of the left side of formula (1), where J_{right} and J_{left} are taken from the fifth and sixth columns of Table 1. It can be seen that these experimental additions confirm the description of the time-background dependence N (t) obtained in (Mikhailov and Dulin 1995) and are essential in terms of determining the total area J.
Experimental results, taking into account the experimental additives
Files
Starts C, s –1
Channel background on the right
α×10^{5}, s^{–1}
J_{left}
J_{right}
RT-1
1544 ± 44
8666
1.191
36811
73436
RT-21
1816 ± 46
13836
1.002
56401
121326
RT-3
1796 ± 15
14993
1.119
62041
130442
Comparison of the results of the Rossi-α and IKES methods
Files
J_{left} + ΔJ_{left}
J_{right} + ΔJ_{right}
N/(2J_{right}ΔtQ)
N/[(J_{right}+J_{left})ΔtQ]
RT-1
73150
74376
0.0859
0.0858
RT-21
111799
123485
0.0820
0.0869
RT-3
123102
132503
0.0831
0.0859
0.0837 ± 0.0017
0.0862 ± 0.0005
k_{eff} (Rossi-α) – k_{eff} (IKES)
0.9624 ± 0.0008 –
0.9630 ± 0.0003 –
– 0.9595 ± 0.0030 =
– 0.9595 ± 0.0030 =
= 0.0029
= 0.0035
The source of spontaneous fissions Q_{sp} was the ^{240}Pu source, which was included in the plutonium fuel. Using the calculated factor D = 0.0472 and the effective source of spontaneous fissions ^{240}Pu Q_{sp} = (7.64 ± 0.21)×10^{6} n/s (1σ), we obtain the prompt neutron multiplication factor 1/(1 – k_{σ}_{p}_{exp}) and the values k_{eff}(φ^{+}_{σ}) = k_{σ}_{p}_{exp} + β_{eff}. In this method the absolute error β_{eff} but not the relative error, as in the |ρ_{$}|-subcriticality (in units of β_{eff}) of the assembly. The difference between the k_{eff} values obtained by different methods is shown in the last row of the Table 2.
Table 3 shows the results of using the Rossi-α experiment for the left-hand side of (1) and the result of measurements of the same subcritical assembly by the IKES method (Grabezhnoy et al. 2015).
The results of using the experiment of Rossi-α and measurements of the subcritical assembly by the IKES method
A similar method was used to determine the subcritical state of a single-core BFS-73 facility, consisting of metallic uranium (18% enrichment), steel and sodium, in proportions close to the composition of a fast reactor (Dulin VA and Dulin VV 1998).
The multiplication factor, estimated by the standard methods adopted at the BFS facilities (by the IKES method (Kazansky and Matusevich 1994) and by the method of inverse neutron multiplication), was k = 0.90 ± 0.01.
The processing of these experiments carried out in (Marchuk and Orlov 1961), taking into account the results of (Mikhailov and Dulin 1995), made it possible to more correctly calculate the space-isotope correlation factor for a medium with a source; however, when calculating J_{left}(t), the integration of C_{left}(t) from (2) over t was carried out in the range from 0 to Т, but not from –∞ to Т, as follows from (Mikhailov and Dulin 1995). Now this is done by means of integration similar to (2), (3) and (6), (7), analyzing the previously measured files with an accurate knowledge of the delay T, but not assuming, as was done in (Grabezhnoy et al. 2015), that it is determined by the position of the peak maximum visible on the screen.
The results of using the Rossi-α experiment and measurements of the subcritical assembly by the IKES method presented in Table 3 are similar to the process described above. Since, in this case, the counters were at the border of the core, and the Cf sources were 4–10 times smaller than the source of neutrons of spontaneous ^{240}Pu fissions, when files of type S were measured, covers made of polyethylene 1 mm thick were put on the counters and, in the case of files of type R, the covers were 4.1 mm thick.
The data of the IKES method differ insignificantly from the results of the Rossi-α method with an increase in the subcriticality of the assembly (the difference did not exceed the errors). Perhaps, this was due to the simplicity of the geometry of the subcritical single-core BFS-73 facility (Marchuk and Orlov 1961) and, accordingly, the small error in calculating the spatial corrections in both methods. Since in the IKES method, the error in β_{eff} directly enters the error in k_{eff} (IKES) (as an error in the denominator), and the error in k_{eff} (Rossi-α) — only as an error in the sum k_{eff}(φ^{+}_{σ}) = k_{σ}_{p} + β_{eff,} then (if it is true) a small value of the difference k_{eff}(Rossi-α) – k_{eff}(IKES) indicates a small error in β_{eff} in the IKES method.
ReferencesDulinVADulinVV (1998) On the measurement of the subcritical states of the breeding media by the Rossi-alpha method. Izvestia Vysshikh Uchebnykh Zavedeniy. Yadernaya Energetika [News of Higher Educational Institutions.DulinVAMatvienkoIP (2002) Definition of deep subcriticality of breeding media states by Rossi-α method. Izvestia Vysshikh Uchebnykh Zavedeniy. Yadernaya Energetika [News of Higher Educational Institutions.GrabezhnoyVADulinVADulinVV (2015) On spatial corrections in Rossi-α experiments. Izvestia Vysshikh Uchebnykh Zavedeniy. Yadernaya Energetika [News of Higher Educational Institutions.KazanskyYuAMatusevichEC (1994) Experimental Physics of Reactors. Energoatomizdat Publ., Moscow, 116–119. [in Russian]MarchukGIOrlovVV (1961) On the theory of conjugate functions. In: Neutron Physics. Gosatomizdat Publ., Moscow, 31–34. [in Russian]MikhailovGMDulinVA (1995) Measurement of β_{eff} by the Rossi-α Method.TakeshiSShigeakiO (1970) Experimental core for experiments of effective delayed fraction β_{eff}.TsutomuI (1968) On the Background Counts in the Rossi-alpha Experiment. Nukleonik, Band 11, Heft 3: 157.
* Russian text published: Izvestiya vuzov. Yadernaya Energetika (ISSN 0204-3327), 2021, n. 2, pp. 50–58.