Measurement of k_eff by delayed neutron multiplication in subcritical systems

Measurement of k eff by delayed neutron multiplication in subcritical systems 1 2 first-author 1 2 corresp Corresponding author. E-mail Address:zhengchun@caep.cn Corresponding author. E-mail Address:zhengchun@caep.cn zhengchun@caep.cn College of Physical Science and Technology, Sichuan University, Chengdu 610064, China College of Physical Science and Technology, Sichuan University, Chengdu 610064, China Institute of Nuclear Physics and Chemistry, China Academy of Engineering Physics, Mianyang 621900, China Institute of Nuclear Physics and Chemistry, China Academy of Engineering Physics, Mianyang 621900, China In this paper, we build on the concept of equivalent fundamental-mode source to propose using delayed neutrons as a neutron source in multiplication experiments to acquire the effective multiplication factor k eff of subcritical systems, which is difficult to acquire directly from conventional neutron source multiplication method. We analyzed the difference between a fundamental-mode fission source and delayed neutron source, then adopted a factor to convert delayed neutron distribution to an equivalent fundamental-mode source distribution and employed Monte Carlo code to acquire this factor. The delayed neutron multiplication measurement method was established for the first time, and corresponding experiments were conducted in subcritical systems. The multiplication of delayed neutrons was measured based on Chinese Fast Burst Reactor-II (CFBR-II) at subcritical states, and k eff was acquired from delayed neutron multiplication successfully (0.9921 and 0.9969, respectively). The relative difference between k eff obtained by the new method and previous values acquired by the positive period method is less than 1% for these two studied cases. k eff Delayed neutrons Neutron source multiplication Subcritical system 1 Introduction 1 Subcriticality is an important concept for monitoring the safety of reactor operation in critical or near-critical states. The neutron source multiplication (NSM) method is widely used to acquire nuclear system reactivity in nuclear safety analysis [ 1 1 - 7 7 1-7 ]. Experimental NSM investigations were carried out in several facilities worldwide in recent years, including the Kyoto University Critical Assembly "KUCA" facility [ 8 8 - 10 10 8-10 ], the YALINA thermal subcritical assembly [ 11 11 , 12 12 ], and the MASURCA reactor through the series of MUSE experiments of the fast subcritical mock-up [ 13 13 - 15 15 13-15 ]. These studies investigated the nuclear characteristics in NSM experiments and examined the neutron properties of subcritical systems. Neutrons involved in those experiments were 14 MeV pulsed neutrons, 252 Cf spontaneous fission neutrons, deuterium–deuterium (D–D) and deuterium–tritium (D–T) neutrons. Different neutron source distributions (distributed in space, energy, and angle) have different effects on the multiplication of the nuclear systems. The NSM method is a simple measurement technique which may feasibly be carried out in real-time. This method does not require any special equipment other than a stationary external neutron source, and an ordinary neutron detector. Additionally, the NSM method is based on steady-state analysis, so that this technique is very suitable for real-time measurement. Despite these advantages, the absolute value of subcriticality cannot be measured directly by the NSM method. One of the important reasons is that extra neutron sources used in the NSM experiments are very different from eigen-distribution neutron source (fundamental-mode fission sources) in space, energy, and angle distribution. It is impossible for extra neutron sources to distribute uniformly throughout the system volume. When the neutron sources are eigen-distribution neutron sources, the effective multiplication factor k eff can be directly related to the neutron multiplication factor M by M =1/(1- k eff ). However, the fundamental-mode source is a fictitious source distribution that does not exist in reality. In order to establish the relationship between the ordinary neutron source and fundamental-mode source, G.D.Spriggs et al. [ 16 16 ] have developed the theory of equivalent fundamental-mode sources, and demonstrated a factor that could convert any arbitrary source distribution to an equivalent fundamental-mode source distribution. They also demonstrated a method for calculating this factor in subcritical systems. In China, Jin-Feng Du and Yan-Peng Yin [ 17 17 - 18 18 17-18 ] developed the Monte Carlo analysis method to calculate the equivalent fundamental-mode source strength of 14 MeV neutrons and spontaneous fission neutron sources in CFBR-II reactors, respectively. Two neutron source distributions (distributed in space, energy and angle) producing the same number of neutrons will not necessarily contribute equally toward the multiplication of the given system. For this reason, equally sized units added to the system will have different effects on the multiplication of the system. Fundamental-mode neutrons and delayed neutrons that are all produced from fission events in critical system should be identically distributed in space and angle, varying only in energy spectrum. Therefore, we could introduce a factor to convert delayed neutron distribution to an equivalent fundamental-mode source distribution, i.e., introduce a factor to eliminate the energy distribution difference between eigen-distribution neutron sources and delayed neutron sources in the NSM experiments to acquire k eff of the nuclear system. In this paper, we build on the concept of equivalent fundamental-mode source to propose using delayed neutrons as an extra neutron source in NSM experiments to acquire k eff of the subcritical system. The paper is structured to begin with a statement of basic theory; thereafter, our use of the Monte Carlo analysis method to calculate the equivalent fundamental-mode source strength of delayed neutrons is demonstrated; then, the multiplication experiments we conducted in subcritical systems using delayed neutrons as a driving source to measure k eff are described; and finally, conclusions are drawn, and future work is outlined. 2 Delayed neutron multiplication 1 In this section, we describe a novel method to determine k eff of subcritical system using delayed neutron multiplication. Delayed neutrons are born from fission events and spread all over the system. They can induce fission events and proliferate neutrons like other ordinary neutrons. Delayed neutron multiplication can be measured by an appropriate method. The determination of k eff from delayed neutron multiplication is predicated based on an understanding of the concept of equivalent fundamental-mode source. In a subcritical system, if the source S was distributed as a fundamental-mode fission source, the fundamental-mode multiplication, M 0 , would be related to the effective multiplication factor of the system, k eff , as follows: <math display="block" id="M1"><mrow><msub><mtext>M</mtext><mn>0</mn></msub><mo>=</mo><mfrac><mn>1</mn><mrow><mn>1</mn><mo>−</mo><msub><mi>k</mi><mrow><mi>e</mi><mi>f</mi><mi>f</mi></mrow></msub></mrow></mfrac><mo>.</mo></mrow></math> However, the equivalent fundamental-mode source is a fictitious theoretical source. In experiments, arbitrary sources that have been placed in or near the assembly (such as an external startup neutron source) will produce system multiplication, M , that can differ significantly from the fundamental-mode multiplication, M 0 . Because it is customary in reactor physics to express most quantities in terms of k eff , it is necessary to modify Eq.(1) by including a factor g * that allows us to express the actual multiplication produced by an arbitrary source distribution in terms of the fundamental-mode multiplication. That is, <math display="block" id="M2"><mrow><mi>M</mi><mo>=</mo><msup><mi>g</mi><mo>∗</mo></msup><msub><mi>M</mi><mn>0</mn></msub><mo>=</mo><mfrac><mrow><msup><mi>g</mi><mo>∗</mo></msup></mrow><mrow><mn>1</mn><mo>−</mo><msub><mi>k</mi><mrow><mtext>eff</mtext></mrow></msub></mrow></mfrac><mo>,</mo></mrow></math> <math display="block" id="M3"><mrow><msup><mi>g</mi><mo>∗</mo></msup><mo>=</mo><mfrac><mrow><mover accent="true"><mrow><msubsup><mi>ψ</mi><mi>s</mi><mo>∗</mo></msubsup></mrow><mo stretchy="true">¯</mo></mover></mrow><mrow><mover accent="true"><mrow><msubsup><mi>ψ</mi><mi>f</mi><mo>∗</mo></msubsup></mrow><mo stretchy="true">¯</mo></mover></mrow></mfrac><mo>,</mo></mrow></math> where <math id="M4"><mrow><mover accent="true"><mrow><msubsup><mi>ψ</mi><mi>s</mi><mo>∗</mo></msubsup></mrow><mo stretchy="true">¯</mo></mover></mrow></math> is the average importance of a source neutron, and <math id="M5"><mrow><mover accent="true"><mrow><msubsup><mi>ψ</mi><mi>f</mi><mo>∗</mo></msubsup></mrow><mo stretchy="true">¯</mo></mover></mrow></math> is the Average importance of a fission neutron. This factor <math id="M6"><mrow><msup><mi>g</mi><mo>∗</mo></msup></mrow></math> can be defined as the ratio of the fixed-source multiplication to the fundamental-mode multiplication, and is used to relate a given source strength to its equivalent fundamental-mode source strength [ 16 16 ]. Therefore, <math id="M7"><mrow><msup><mi>g</mi><mo>∗</mo></msup></mrow></math> is multiplied by the source strength of the neutron source to yield a fundamental-mode source strength, Q , defined as <math id="M8"><mrow><mi>Q</mi><mo>=</mo><msup><mi>g</mi><mo>∗</mo></msup><mi>S</mi></mrow></math> . Therefore, the product, <math id="M9"><mrow><msup><mi>g</mi><mo>∗</mo></msup><mi>S</mi></mrow></math> , represents the equivalent fundamental-mode source. A method for determining k eff by delayed neutron multiplication will now be described. Delayed neutrons need to be acquired first. To produce these, a nuclear device is operated in critical mode for a moment and delayed neutron precursors are produced. Delayed neutrons precursors distribute throughout the system, and then, the device is shut down to stop the fission events. Subsequently the existent precursors in the system decay. At this point, only delayed neutrons exist in the system (and no more prompt neutrons). A few minutes later, the reactor enters a subcritical state from shut-down mode; thus, delayed neutron multiplication is realized. Then, k eff of the system can be acquired according to the relationship between k eff and delayed neutron multiplication ( <math id="M10"><mrow><msub><mi>M</mi><mrow><mi>d</mi><mi>n</mi></mrow></msub></mrow></math> ). The multiplication of delayed neutrons in every step above is now described. While the reactor is running, delayed neutron precursors are all continuously being produced and decaying. Precursors release neutrons called delayed neutrons. When the reactor operates in critical mode for more than ten minutes, precursors′ generation and decay will reach dynamic equilibrium according to the characteristics of delayed neutron precursors. The neutron detectors count rate, <math id="M11"><mrow><msub><mi>n</mi><mn>1</mn></msub></mrow></math> , in the system should be: <math display="block" id="M12"><mrow><msub><mi>n</mi><mn>1</mn></msub><mo>=</mo><mi>ε</mi><mo stretchy="false">(</mo><mi>ν</mi><mo>−</mo><mi>α</mi><mo>−</mo><mn>1</mn><mo stretchy="false">)</mo><mover><mi>F</mi><mo>.</mo></mover><mo>,</mo></mrow></math> where <math id="M13"><mi>ε</mi></math> is the neutron detector efficiency, <math id="M14"><mi>ν</mi></math> is the Average yield of neutrons per prompt fission, <math id="M15"><mi>α</mi></math> is the ratio of prompt neutrons capture to fission cross section, and <math id="M16"><mrow><mover><mi>F</mi><mo>.</mo></mover></mrow></math> is the fission rate of reactor in critical mode (s -1 ), which can be calculated from the running power: <math id="M17"><mrow><mover><mi>F</mi><mo>.</mo></mover><mo>=</mo><mi>P</mi><mo>×</mo><mn>3.1</mn><mo>×</mo><msup><mrow><mn>10</mn></mrow><mrow><mn>10</mn></mrow></msup></mrow></math> . When the density of delayed neutron precursors reaches saturation, the reactor is shut-down. Fission events cease and delayed neutron precursor production also ceases. The precursors that existed in the whole system subsequently decay. The decay properties of precursors determine that the shorter a precursor’s half-life, the faster it disappears. According to six groups of delayed neutron model, the 1-th precursor with the longest half-life ( 87 Br: T 1/2 =55.6 s in ENDF/B-VI) has the slowest decay behavior. Thereafter, 300 s after the fission events stop, only 1-th delayed neutrons exist in the system [ 19 19 ]. The emission rate (which can be defined as neutron source intensity) of the 1-th delayed neutrons at time t beginning from the shut-down moment, <math id="M18"><mrow><msub><mi>S</mi><mrow><mi>d</mi><mi>n</mi></mrow></msub></mrow></math> , can be given by <math display="block" id="M19"><mrow><msub><mi>S</mi><mrow><mi>d</mi><mi>n</mi></mrow></msub><mo>=</mo><mover><mi>F</mi><mo>.</mo></mover><mo>×</mo><mo stretchy="false">(</mo><mn>1</mn><mo>−</mo><msup><mi>e</mi><mrow><mo>−</mo><mi>λ</mi><mi>T</mi></mrow></msup><mo stretchy="false">)</mo><mo>×</mo><msub><mi>ν</mi><mi>d</mi></msub><mo>×</mo><msub><mi>a</mi><mn>1</mn></msub><mo>×</mo><msup><mi>e</mi><mrow><mo>−</mo><mi>λ</mi><mi>t</mi></mrow></msup><mo>,</mo></mrow></math> where t is the time, start from the reactor shutdown(s); T is the time duration while reactor run in critical mode (s); <math id="M20"><mrow><msub><mi>ν</mi><mtext>d</mtext></msub></mrow></math> the total yield of delayed neutron precursors per fission; a 1 is the relative abundance of 1-th delayed neutron precursor; λ is the radioactive decay constant of 1-th delayed neutron precursor(s -1 ); and <math id="M21"><mrow><msub><mi>S</mi><mrow><mi>d</mi><mi>n</mi></mrow></msub></mrow></math> is the delayed neutrons source intensity at time t (s -1 ). If the time duration, T , that the system is kept at in a critical state is greater than ten half-lives of 87 Br, all the precursors could reach saturation, giving: <math id="M22"><mrow><mn>1</mn><mo>−</mo><msup><mi>e</mi><mrow><mo>−</mo><mi>λ</mi><mi>T</mi></mrow></msup></mrow></math> ≈1. Then, Eq. (5) should be <math display="block" id="M23"><mrow><msub><mi>S</mi><mrow><mi>d</mi><mi>n</mi></mrow></msub><mo>=</mo><mover><mi>F</mi><mo>.</mo></mover><mo>×</mo><msub><mi>ν</mi><mtext>d</mtext></msub><mo>×</mo><msub><mi>a</mi><mn>1</mn></msub><mo>×</mo><msup><mi>e</mi><mrow><mo>−</mo><mi>λ</mi><mi>t</mi></mrow></msup><mo>.</mo></mrow></math> Alternatively, according to the theory of equivalent fundamental-mode source, 1-th delayed neutrons fundamental-mode source strength ( <math id="M24"><mrow><msub><mi>Q</mi><mrow><mi>d</mi><mi>n</mi></mrow></msub></mrow></math> ) could be written as <math display="block" id="M25"><mrow><msub><mi>Q</mi><mrow><mi>d</mi><mi>n</mi></mrow></msub><mo>=</mo><msubsup><mi>g</mi><mrow><mi>d</mi><mi>n</mi></mrow><mo>*</mo></msubsup><mo>×</mo><msub><mi>S</mi><mrow><mi>d</mi><mi>n</mi></mrow></msub><mo>,</mo></mrow></math> where <math id="M26"><mrow><msub><mi>g</mi><mrow><mtext>dn</mtext></mrow></msub><msup><mrow/><mo>∗</mo></msup></mrow></math> is defined as the ratio of the delayed neutron multiplication to the fundamental-mode multiplication, as previously mentioned. More than 300 s after the fission events stop, the reactor steps into a subcritical state, the relationship between multiplication ( M ), k eff , the multiplication of the system, leakage of the neutrons and efficiency of the detector, and the neutron detector counting rate <math id="M27"><mrow><msub><mi>n</mi><mn>2</mn></msub></mrow></math> should be <math display="block" id="M28"><mrow><msub><mi>n</mi><mn>2</mn></msub><mo>=</mo><mi>ε</mi><mo>×</mo><mfrac><mrow><mi>ν</mi><mo>−</mo><mi>α</mi><mo>−</mo><mn>1</mn></mrow><mi>ν</mi></mfrac><mo>×</mo><msub><mi>Q</mi><mrow><mi>d</mi><mi>n</mi></mrow></msub><mo>×</mo><mi>M</mi><mo>,</mo></mrow></math> where all the parameters in the Eq. (8) are the same as above mentioned. Substituting Eq. (6) and Eq. (7) into Eq. (8), <math id="M29"><mrow><msub><mi>n</mi><mn>2</mn></msub></mrow></math> can be expressed as <math display="block" id="M30"><mrow><msub><mi>n</mi><mn>2</mn></msub><mo>=</mo><mi>ε</mi><mo>×</mo><mfrac><mrow><mi>ν</mi><mo>−</mo><mi>α</mi><mo>−</mo><mn>1</mn></mrow><mi>ν</mi></mfrac><mo>×</mo><msub><mi>g</mi><mrow><mi>d</mi><mi>n</mi></mrow></msub><msup><mrow/><mo>∗</mo></msup><mo>×</mo><mover><mi>F</mi><mo>.</mo></mover><mo>×</mo><msub><mi>ν</mi><mtext>d</mtext></msub><mo>×</mo><msub><mi>a</mi><mn>1</mn></msub><mo>×</mo><msup><mi>e</mi><mrow><mo>−</mo><mi>λ</mi><mi>t</mi></mrow></msup><mo>×</mo><mfrac><mn>1</mn><mrow><mn>1</mn><mo>−</mo><msub><mi>k</mi><mrow><mtext>eff</mtext></mrow></msub></mrow></mfrac><mo>.</mo></mrow></math> Combined with the Eq. (4) and Eq. (9), Eq.(10) is obtained by <math display="block" id="M31"><mrow><mfrac><mrow><msub><mi>n</mi><mn>2</mn></msub></mrow><mrow><msub><mi>n</mi><mn>1</mn></msub></mrow></mfrac><mo>=</mo><mfrac><mrow><msub><mi>a</mi><mn>1</mn></msub><mo>×</mo><msub><mi>ν</mi><mi>d</mi></msub><mo>×</mo><msup><mi>e</mi><mrow><mo>−</mo><mi>λ</mi><mi>t</mi></mrow></msup><mo>×</mo><msub><mi>g</mi><mrow><mi>d</mi><mi>n</mi></mrow></msub><msup><mrow/><mo>∗</mo></msup></mrow><mrow><mo stretchy="false">(</mo><mn>1</mn><mo>−</mo><msub><mi>k</mi><mrow><mtext>eff</mtext></mrow></msub><mo stretchy="false">)</mo><mo>×</mo><mi>ν</mi></mrow></mfrac><mo>.</mo></mrow></math> Then, <math id="M32"><mrow><msub><mi>k</mi><mrow><mtext>eff</mtext></mrow></msub></mrow></math> can be obtained as follows <math display="block" id="M33"><mrow><msub><mi>k</mi><mrow><mtext>eff</mtext></mrow></msub><mo>=</mo><mn>1</mn><mo>−</mo><mo stretchy="false">(</mo><mfrac><mrow><msub><mi>a</mi><mn>1</mn></msub><mo>×</mo><msub><mi>ν</mi><mi>d</mi></msub><mo>×</mo><msup><mi>e</mi><mrow><mo>−</mo><mi>λ</mi><mi>t</mi></mrow></msup><mo>×</mo><msub><mi>g</mi><mrow><mi>d</mi><mi>n</mi></mrow></msub><msup><mrow/><mo>∗</mo></msup></mrow><mi>ν</mi></mfrac><mo stretchy="false">)</mo><mo>×</mo><mfrac><mrow><msub><mi>n</mi><mn>1</mn></msub></mrow><mrow><msub><mi>n</mi><mn>2</mn></msub></mrow></mfrac><mo>.</mo></mrow></math> Eq. (11) is the basic equation of this work. With the parameters <math id="M34"><mrow><msub><mi>a</mi><mn>1</mn></msub></mrow></math> 、 <math id="M35"><mrow><msub><mi>ν</mi><mi>d</mi></msub></mrow></math> 、 <math id="M36"><mi>λ</mi></math> 、 <math id="M37"><mi>ν</mi></math> , when <math id="M38"><mrow><msub><mi>n</mi><mn>1</mn></msub></mrow></math> and <math id="M39"><mrow><msub><mi>n</mi><mn>2</mn></msub></mrow></math> are measured and <math id="M40"><mrow><msubsup><mi>g</mi><mrow><mi>d</mi><mi>n</mi></mrow><mo>*</mo></msubsup></mrow></math> is calculated, the value of <math id="M41"><mrow><msub><mi>k</mi><mrow><mtext>eff</mtext></mrow></msub></mrow></math> can be obtained. 3 Monte Carlo analysis of the effective coefficient g dn * 1 In order to obtain the factor g dn * of 1-th delayed neutrons, the Monte Carlo numerical method was adopted to simulate the reaction and transport process of delayed neutrons and eigen-distribution neutrons, respectively. As we mentioned above in Eq. (3), the effective coefficient g * is defined as the ratio of the fixed-source multiplication to the fundamental-mode multiplication. In this work, a neutron source multiplication method was used to acquire g dn * of 1-th delayed neutrons. Here, the enriched uranium bare sphere assembly model with 8.805 cm radius and 93.71% enrichment 235 U is built. A neutron detector was put outside of the system. The detector was spherical and closely surrounded the system, improving the detector efficiency. The resultant efficiency of leaked neutron detection was 100%. Additionally, the F1 tally card (current integrated over a surface) of MCNP5 [ 20 20 ] was adopted to acquire the average importance of every neutron emitted from the neutron source. The calculated formula is: <math id="M42"><mrow><msub><mi>F</mi><mn>1</mn></msub><mo>=</mo><mstyle displaystyle="true"><mrow><msub><mo>∫</mo><mi>A</mi></msub><mrow><mstyle displaystyle="true"><mrow><msub><mo>∫</mo><mi>μ</mi></msub><mrow><mstyle displaystyle="true"><mrow><msub><mo>∫</mo><mi>t</mi></msub><mrow><mstyle displaystyle="true"><mrow><msub><mo>∫</mo><mi>E</mi></msub><mrow><mi>J</mi><mo stretchy="false">(</mo><mi>r</mi><mo>,</mo><mi>E</mi><mo>,</mo><mi>t</mi><mo>,</mo><mi>μ</mi><mo stretchy="false">)</mo><mtext>d</mtext><mi>E</mi><mtext>d</mtext><mi>t</mi><mtext>d</mtext><mi>μ</mi><mtext>d</mtext><mi>A</mi></mrow></mrow></mstyle></mrow></mrow></mstyle></mrow></mrow></mstyle></mrow></mrow></mstyle><mi>，</mi></mrow></math> where <math id="M43"><mrow><mi>J</mi><mo stretchy="false">(</mo><mi>r</mi><mo>,</mo><mi>E</mi><mo>,</mo><mi>t</mi><mo>,</mo><mi>μ</mi><mo stretchy="false">)</mo><mo>=</mo><mrow><mo>|</mo><mi>μ</mi><mo>|</mo></mrow><mtext>Φ</mtext><mo stretchy="false">(</mo><mi>r</mi><mo>,</mo><mi>E</mi><mo>,</mo><mi>t</mi><mo stretchy="false">)</mo><mi>A</mi></mrow></math> is the neutron surface current. The neutron source distribution produces primary neutrons, and simulates the neutron transport process in the subcritical system. The leaked neutrons are statistically integrated the over the detector surface, and normalized to delayed neutrons or eigen neutrons. The normalized value represents the average importance of a source neutron. The average importance of delayed neutrons and eigen neutrons is acquired, allowing the effective coefficient, g * , of delayed neutrons to be obtained from their ratio. In this work, we used 1-th delayed neutrons to analyze the average importance of 1-th delayed neutrons in a Monte Carlo numerical simulation. The distribution of delayed neutron precursors was in accordance with the relative fission rate in the critical nuclear system volume, as what that of the 1-th delayed neutron precursors. Relative delayed neutron emission rates in the volume are required to determine the relative fission rate along the radius. The FMn tally multiplier card in MCNP5 was used to calculate fission rate with the form: <math id="M44"><mrow><mi>C</mi><mstyle displaystyle="true"><mrow><mo>∫</mo><mrow><mi>φ</mi><mo stretchy="false">(</mo><mi>E</mi><mo stretchy="false">)</mo></mrow></mrow></mstyle><mi>R</mi><mo stretchy="false">(</mo><mi>E</mi><mo stretchy="false">)</mo><mtext>d</mtext><mi>E</mi><mo>,</mo></mrow></math> where φ ( E ) <math id="M45"><mrow/></math> is the energy-dependent fluency (particles/cm 2 ) and can be acquired from F5 card; R ( E ) is an operator of additive and/or multiplicative response functions from the MCNP cross-section libraries or specially designated quantities; the constant C is any arbitrary scalar quantity that can be used for normalization. Here, R ( E ), the reaction cross-sections are microscopic (with units in barns) corresponding to reaction numbers "−6" (total fission cross section) on an FM card. 19 observation points are set at every 0.5 cm interval along the radius. Relative fission rate distribution along the radius is tallied statistically, and the results, which can be used to sample the relative intensity of delayed neutron distribution in the next step, are shown in Fig.1 Fig.1 . 1-th delayed neutron distributions are sampled as follows: (1) A DNB card is used to turned off delayed neutron production in MCNP5; (2) enriched uranium bare sphere assembly is first divided into 18 concentric spherical shells along the radius; (3) relative intensity distribution of 1-th delayed neutrons at the initial time in each fissile spherical shell is according to relative fission rate as shown in Fig.1 Fig.1 ; (4) 1-th delayed neutrons, with the average energy of 250 keV, emitted homogeneously in the whole volume; (5) the detector counts of delayed neutrons, n dn , is given via the F1 tally card. As to eigen neutrons, they are sampled as following: (1) critical calculation is performed to produce eigen neutrons via a KCODE card; (2) an SSW card is used and eigen neutron information is written in a WSSA file; (3) an SSR card is used to read the location, direction, and energy distribution of eigen neutrons from the RSSA file, which called the WSSA file in the former step; (4) the calculation is executed and detector counts for eigen neutrons, n en , are recorded with the use of an F1 tally card. After performing calculations on delayed neutrons and eigen neutrons, the detector counts for these two neutron types are obtained, and their ratio, n dn / n en , is calculated. This is the effective coefficient of delayed neutrons ( g dn * ), which is 1.22. Statistical uncertainty has been controlled within 1%. This result is consistent with expectations. From Eq. (3), the 1-th delayed neutrons’ effective coefficient ( g dn * ) and neutron importance are connected. Neutrons’ importance, determined by the neutrons properties, include direction, energy, and distribution. In other words, when neutrons are injected into subcritical and/or critical system, the neutron population will increase in terms of neutrons′ properties. For eigen-distribution neutrons, a kind of fission neutron produced in a critical system, its "relative neutron importance" should be 1. For 235 U nucleus, a low-energy neutron has a greater probability to induce fission. Delayed neutrons have a softer spectrum than prompt neutrons, so they induce fission events more easily than prompt neutrons. Thus, it is reasonable to have an effective coefficient greater than 1 for delayed neutrons. We can also obtain this conclusion from effective delayed neutron fraction β eff , which is larger than the delayed neutron fraction, β. Furthermore, comparative studies were conducted to demonstrate that the 1-th delayed neutrons’ effective coefficient, g dn * , is larger than the effective coefficient of the total delayed neutrons. In [ 18 18 ] and [ 21 21 ], Yanpeng Yin calculated the effective delayed neutron fraction β eff of CFBR-II, which is 704 pcm. Gregory experimentally studied the effective delayed neutron fraction β eff of fast critical assembly XIX-1, determining it to be 737 pcm. Thus, drawn from experimental results, we could obtain g * of these two according to Eq.(12). The g dn * values obtained in this work and comparative results of the two reference works are listed in Table 1 Table 1 . 1 1 Results in this work 1 1 Reference value [ 18 18 ] 1 1 Reference value [ 21 21 ] 1 1 1.22 1 1 1.083 1 1 1.13 <math display="block" id="M46"><mrow><msub><mi>g</mi><mrow><mi>d</mi><mi>n</mi></mrow></msub><msup><mrow/><mo>∗</mo></msup><mo>=</mo><mfrac><mrow><msub><mi>β</mi><mrow><mtext>eff</mtext></mrow></msub></mrow><mrow><mstyle displaystyle="true"><munderover><mo>∑</mo><mrow><mi>i</mi><mo>=</mo><mn>1</mn></mrow><mn>6</mn></munderover><mrow><msub><mi>β</mi><mi>i</mi></msub></mrow></mstyle></mrow></mfrac><mo>=</mo><mfrac><mrow><msub><mi>β</mi><mrow><mtext>eff</mtext></mrow></msub></mrow><mrow><mn>0.0065</mn></mrow></mfrac><mo>.</mo></mrow></math> The average energy of 1-th delayed neutrons (250 keV) is lower than the average energy of total delayed neutrons (about 450 keV), so it is easier for 1-th delayed neutrons to induce fission. It can be concluded that it is reasonable for the effective coefficient g * of 1-th delayed neutrons to be larger than that of the total delayed neutrons. Owing to this, we propose that using the Monte Carlo method based on the MCNP code to calculate the effective coefficient of 1-th delayed neutrons is practical. G.D. Spriggs et al .calculated the parameter g * of 252 Cf spontaneous fission neutrons [ 21 21 ], and 235 U and 238 U spontaneous fission neutrons [ 16 16 ], but did not establish g * of delayed neutrons. As such, our achievement in this work is an effective complement to their theory of equivalent fundamental-mode sources. 4 Experiment apparatus review 1 Our experimental procedure was performed using the Chinese Fast Burst Reactor-II (CFBR-II). The CFBR-II is a spherical assembly consisting of cast-enriched uranium. Two hemispherical enriched-uranium blocks, referred to as upper and lower safety blocks, are driven mechanically. Between the two safety blocks, there is a 31-cm diameter and 5.2-cm-thick stainless-steel plate. There are three horizontal circular rod slots in the plate, which hold three control rods: an auto-adjustment rod, composition rod, and pulse rod to pass through. Control rods are all made from cast-enriched uranium. The reactivity of the assembly can be controlled by adjusting the depth to which control rods are inserted. The assembly can work at subcritical, delayed critical, and super prompt critical states. Two BF 3 proportional counters were used to record the leaked neutrons (proportional to the neutron density in the system) during measurement. A BF 3 proportional counter, SZJ-1 type, the production of Beijing Nuclear Instrument Factory, was put in the paraffin barrel to construct the Hanson long counter, along with a preamplifier, a high-voltage power supply, a main amplifier, and a multi-scalar to constitute the measurement system. The front of the detectors was 180 cm away from the center of the system, and the dead time of the detectors was 3 μs. Fig.2 Fig.2 demonstrates the system structure and the location of detectors. 45 0 45 0 Stainless steel plate180cmBF 3 180cmBF 3 CFBR-IIUpper blockLower block(a) 5 Experiments and results analysis 1 We performed experiments as follows: (1) the reactor was kept in critical state for more than 10 min and the delayed neutron precursors reached saturation, at which point the neutron detectors count rate ( <math id="M47"><mrow><msub><mi>n</mi><mn>1</mn></msub></mrow></math> ) in experiments was 2.07e7 s -1 (case 1: 50 W running power) and 1.25e7 s -1 (case 2: 30 W running power); (2) we drew out a control rod (made of enriched uranium) and separated the two safety blocks promptly, indicating the system would be in a shut-down state and the fission events stopped, precursors could not produce, and the ones that already existed in the system decayed subsequently; 300 s after the reactor shut down, only 1-th delayed neutrons existed in the system and other groups of delayed neutrons disappeared, as we outlined in section 3; (3) upper and lower safety blocks of the assembly connected together rapidly, taking reactor from a shut-down state to a subcritical state. To this point, the multiplication of 1-th delayed neutrons had been realized. Two neutron detectors counted the leaked neutrons of the experimental system. Taking into consideration the fluctuation and noise effect on the count rate, the delayed neutrons’ behaviors after the system shut-down are shown in Fig.3 Fig.3 . Curves were fitted with an exponential function, and the half-life of 87 Br was obtained as 55.45 s and 54.79 s, respectively, in these two experiments. The results are in good agreement with the half-life of 87 Br, 55.6 s(ENDF/B-VI). The supplementary specifications for the neutron count rate n 1 and n 2 are as follows. These rates are acquired from a common, current method, which is often used when the neutron count rate is too large for detectors to handle. A brief description of the method is as follows. First, the relationship between electric current and neutron count rate is determined from the ionization chamber and the neutron detector, working simultaneously at low power. The detector count rate when the nuclear assembly is working at high power can be determined by extrapolation from readings obtained via the ionization chamber. With the constants <math id="M48"><mrow><msub><mi>a</mi><mn>1</mn></msub></mrow></math> (0.0395 for 235U), <math id="M49"><mrow><msub><mi>ν</mi><mi>d</mi></msub></mrow></math> (0.0163 for 235U), <math id="M50"><mi>λ</mi></math> (0.0127 s -1 for 87Br), <math id="M51"><mi>ν</mi></math> (2.59 for 235 U prompt fission) [ 22 22 ], detector count rates n 1 and <math id="M52"><mrow><msub><mi>n</mi><mn>2</mn></msub></mrow></math> , and the calculated <math id="M53"><mrow><msubsup><mi>g</mi><mrow><mi>d</mi><mi>n</mi></mrow><mo>*</mo></msubsup></mrow></math> (1.22), the average value of k eff from all points in Fig.3 Fig.3 (a) and (b) can be obtained based on Eq.(11). The average value of k eff was 0.9921 and 0.9969, respectively, in these two experiments. The comparison between the results obtained through the reference method (positive period method: identified by <math id="M54"><mrow><msubsup><mi>k</mi><mrow><mtext>eff</mtext></mrow><mrow><mtext>Ref</mtext></mrow></msubsup></mrow></math> ) and our proposed method (identified by k eff ) is calculated by this formula: <math id="M55"><mrow><mi>Δ</mi><mi>ε</mi><mo>=</mo><mo>|</mo><mtext>(</mtext><msub><mi>k</mi><mrow><mtext>eff</mtext></mrow></msub><mo>-</mo><msubsup><mi>k</mi><mrow><mtext>eff</mtext></mrow><mrow><mtext>Ref</mtext></mrow></msubsup><mo stretchy="false">)</mo><mo>/</mo><msubsup><mi>k</mi><mrow><mtext>eff</mtext></mrow><mrow><mtext>Ref</mtext></mrow></msubsup><mtext>|</mtext><mo>×</mo><mtext>100%</mtext></mrow></math> . The comparison results are listed in Table 2 Table 2 . The accuracy of the results can be attributed to the power of the nucleus device when kept in steady mode. The reactor ran at a power level of 50 W in case 1 and 30 W in case 2. The greater the power level, the more fission neutrons are produced in the system, the more neutrons leak from reactor, and more easily the neutron detectors catch the neutrons. This means that the count rate of detectors is influenced by the power of the nuclear system. Thus, the statistical uncertainty of detectors is lower in case 1 and its results are closer to the positive period method results. 1 1 case 1 1 Results in this work( k eff ) 1 1 positive period method( <math id="M56"><mrow><msubsup><mi>k</mi><mrow><mtext>eff</mtext></mrow><mrow><mtext>Ref</mtext></mrow></msubsup></mrow></math> ) 1 1 ∆ ε % 1 1 1 1 1 0.9921 1 1 0.9922 1 1 0.01 1 1 2 1 1 0.9969 1 1 0.9939 1 1 0.30 Here is the complementary to the reference method (the positive period method). The control rods of CFBR-II were calibrated using the positive period method. The curves describing the relationship between reactivity and the position of control rods are well established. When the reactivity of the system is changed by adjusting the position of the control rods, the reactivity of the system can be read off the curves. 6 Conclusion 1 A methodology to use delayed neutrons as the neutron source in neutron multiplication experiments to acquire k eff was first proposed and investigated via preliminary NSM experiments on CFBR-II. The evaluated results were approximately equal to the reference method (the positive period method). An experimental delayed neutron multiplication method was established, and the Monte-Carlo method was adopted to analyze the effective coefficient g * of delayed neutrons. The effective neutron multiplication factor k eff was obtained experimentally in two subcritical states. The results were approximately equal to the reference method (the positive period method) and were not heavily dependent on detection efficiency. We posit that running the reactor at an appropriate power level would reduce the relative differences in experimental results. We conclude that measurement of k eff by delayed neutron multiplication in subcritical systems, as we outlined in this paper, is feasible for simple subcritical apparatus, such as CFBR-II. References Y.Q. Shi , Q.F. Zhu , P. Xia et al ., Neutron source multiplication method research in reactor physics experimen . Nucl. Sci. Eng . 25 , 14 - 19 ( 2005 ). doi: 10.3321/j.issn:0258-0918.2005.01.002 http://doi.org/10.3321/j.issn:0258-0918.2005.01.002 Neutron source multiplication method research in reactor physics experimen G. Truchet , W.F.G. van Rooijen , Y. Shimazu et al ., Application of the modified neutron source multiplication method to the prototype FBR Monju . Ann. Nucl. Energy . 51 , 94 - 106 ( 2013 ). doi: 10.1016/j.anucene.2012.07.040 http://doi.org/10.1016/j.anucene.2012.07.040 Application of the modified neutron source multiplication method to the prototype FBR Monju G.M. Sun , M.S. Cheng , Z.M. Dai , Preliminary analysis of fuel management for a small modular molten salt fast reactor . Nucl. Sci. Tech . 39 , 070603 ( 2016 ). doi: 10.11889/j.0253-3219.2016.hjs.39.070603 http://doi.org/10.11889/j.0253-3219.2016.hjs.39.070603 Preliminary analysis of fuel management for a small modular molten salt fast reactor K. Yang , J.G. Chen , X.Z. Cai , Using the neutron balance method to access the feed fuel requirements for CANDLE . Nucl. Sci. Tech . 39 , 060601 ( 2016 ). doi: 10.11889/j.0253-3219.2016.hjs.39.060601 http://doi.org/10.11889/j.0253-3219.2016.hjs.39.060601 Using the neutron balance method to access the feed fuel requirements for CANDLE T. Endo , A. Yamamoto , Y. Yamane , Detected-neutron multiplication factor measured by neutron source multiplication method . Ann. Nucl. Energy . 38 , 2417 - 2427 ( 2011 ). doi: 10.1016/j.anucene.2011.01.007 http://doi.org/10.1016/j.anucene.2011.01.007 Detected-neutron multiplication factor measured by neutron source multiplication method B. Ye , L.N. Zeng , Y.P. Yin et al ., Measurement of keff with an improved neutron source multiplication method based on numerical analysis . Nucl. Sci. Tech . 25 , 020602 ( 2014 ). doi: 10.13538/j.1001-8042/nst.25.020602 http://doi.org/10.13538/j.1001-8042/nst.25.020602 Measurement of keff with an improved neutron source multiplication method based on numerical analysis H.W. Xiao , G.M. Liu , H. Yao et al ., Influence of extended low power operation on neutronics parameters of CNP600 . Nucl. Sci. Tech . 25 , 110602 ( 2016 ). doi: 10.11889/j.0253-3219.2016.hjs.39.110602 http://doi.org/10.11889/j.0253-3219.2016.hjs.39.110602 Influence of extended low power operation on neutronics parameters of CNP600 H. Shahbunder , C.H. Pyeon , T. Misawa et al ., Experimental analysis for neutron multiplication by using reaction rate distribution in accelerator-driven system . Ann. Nucl. Energy . 37 , 592 - 597 ( 2010 a). doi: 10.1016/j.anucene.2009.12.022 http://doi.org/10.1016/j.anucene.2009.12.022 Experimental analysis for neutron multiplication by using reaction rate distribution in accelerator-driven system H. Shahbunder , C.H. Pyeon , T. Misawa et al ., Subcritical multiplication factor and source efficiency in accelerator-driven system . Ann. Nucl. Energy . 37 , 1214 - 1222 ( 2010 b). doi: 10.1016/j.anucene.2010.04.010 http://doi.org/10.1016/j.anucene.2010.04.010 Subcritical multiplication factor and source efficiency in accelerator-driven system H. Shahbunder , C.H. Pyeon , T. Misawa et al ., Effects of neutron spectrum and external neutron source on neutron multiplication parameters in accelerator-driven system . Ann. Nucl. Energy . 37 , 1785 - 1791 ( 2010c ). doi: 10.1016/j.anucene.2010.07.003 http://doi.org/10.1016/j.anucene.2010.07.003 Effects of neutron spectrum and external neutron source on neutron multiplication parameters in accelerator-driven system C.M. Persson , P. Seltborg , A. Ahlander et al ., Analysis of reactivity determination methods in the subcritical experiment Yalina . Nucl. Instrum. Meth. A . 554 , 374 - 383 ( 2005 ). doi: 10.1016/j.nima.2005.07.058 http://doi.org/10.1016/j.nima.2005.07.058 Analysis of reactivity determination methods in the subcritical experiment Yalina C.M. Persson , A. Fokau , I. Serafimovich et al ., Pulsed neutron source measurements in the subcritical ADS experiment YALINA-Booster . Ann. Nucl. Energy . 35 , 2357 - 2364 ( 2008 ). doi: 10.1016/j.anucene.2008.07.011 http://doi.org/10.1016/j.anucene.2008.07.011 Pulsed neutron source measurements in the subcritical ADS experiment YALINA-Booster D.L. Chichester , M.T. Kinlaw , The MARVEL assembly for neutron multiplication . Appl. Radiat. Isotopes . 80 , 42 - 48 ( 2013 ). doi: 10.1016/j.apradiso.2013.05.012 http://doi.org/10.1016/j.apradiso.2013.05.012 The MARVEL assembly for neutron multiplication R. Soule , W. Assal , P. Chaussonnet et al ., Neutronic Studies in Support of Accelerator-Driven Systems: The MUSE Experiments in the MASURCA Facility . Nucl. Sci. Eng . 148 , 124 - 152 ( 2004 ). doi: 10.13182/NSE01-13C http://doi.org/10.13182/NSE01-13C Neutronic Studies in Support of Accelerator-Driven Systems: The MUSE Experiments in the MASURCA Facility A. Billebaud , R. Brissot , C. LeBrun et al ., Prompt multiplication factor measurements in subcritical systems: From MUSE experiment to a demonstration ADS . Prog. Nucl. Energ . 49 , 142 - 160 ( 2007 ). doi: 10.1016/j.pnucene.2006.11.001 http://doi.org/10.1016/j.pnucene.2006.11.001 Prompt multiplication factor measurements in subcritical systems: From MUSE experiment to a demonstration ADS G.D. Spriggs , R.D. Busch , T. Sakurai et al ., The Equivalent Fundamental-Mode Source . Ann. Nucl. Energy . 26 , 237 - 264 ( 1999 ). doi: 10.1016/S0306-4549(98)00048-614 http://doi.org/10.1016/S0306-4549(98)00048-614 The Equivalent Fundamental-Mode Source J.F. Du , X.Q. Fan , Calculation of effective coefficient of spontaneity fission neutron source in CFBR-Ⅱreactor . Nucl. Power. Eng . 31 , 131 - 133 ( 2010 ). doi: 10.7538/yzk.2013.47.04.0619 http://doi.org/10.7538/yzk.2013.47.04.0619 Calculation of effective coefficient of spontaneity fission neutron source in CFBR-Ⅱreactor Y.P. Yin , C. Zheng , P. Huang , Calculation of CFBR-Ⅱ's multiplication-reactivity conversion factor . Chin. J. Comput. Phys . 27 , 799 - 804 ( 2010 ). doi: 10.3969/j.issn.1001-246X.2010.06.002 http://doi.org/10.3969/j.issn.1001-246X.2010.06.002 Calculation of CFBR-Ⅱ's multiplication-reactivity conversion factor C. Zheng , Q. Wang , M. Li et al ., Measurement of total fission number by delayed neutron method for burst reactor . Atom. Energy. Sci. Technol . 41 , 699 - 701 ( 2007 ). doi: 10.3969/j.issn.1000-6931.2007.06.014 http://doi.org/10.3969/j.issn.1000-6931.2007.06.014 Measurement of total fission number by delayed neutron method for burst reactor X-5 Monte Carlo Team Los Alamos National Laboratory , MCNP-A General Monte Carlo N-Particle Transport Code , Version 5 . LA-UR-03-1987 ( 2008 ). MCNP-A General Monte Carlo N-Particle Transport Code G.D. Spriggs , T. Sakurai , S. Okajima , Rossi-α and βeff measurements in a fast critical assembly . Prog. Nucl. Energy . 35 , 169 - 181 ( 1999 ). doi: 10.1016/S0149-1970(99)00010-4 http://doi.org/10.1016/S0149-1970(99)00010-4 Rossi-α and βeff measurements in a fast critical assembly D.E. Gremyachkin , V.M. Piksaikin , K.V. Mitrofanov et al ., Verification of the evaluated fission product yields data from the neutron induced fission of 235U, 238U and 239Pu based on the delayed neutron characteristics . Prog. Nucl. Energy . 83 , 13 - 25 ( 2015 ). doi: 10.1016/j.pnucene.2015.02.016 http://doi.org/10.1016/j.pnucene.2015.02.016 Verification of the evaluated fission product yields data from the neutron induced fission of 235U, 238U and 239Pu based on the delayed neutron characteristics 10.1007/s41365-018-0355-7 This work was supported by the National Natural Science Foundation of China (Nos.11175164, 91326109) 2017-03-03 2017-03-24 2017-03-31 2018-02-07 2018-02-01 2021-04-11 © China Science Publishing & Media Ltd. (Science Press), Shanghai Institute of Applied Physics, the Chinese Academy of Sciences, ChineseNuclear Society and Springer Nature Singapore Pte Ltd. 2018 This is a post-peer-review, pre-copyedit version of an article published in Nuclear Science and Techniques. The final authenticated version is available online at: http://dx.doi.org/10.1007/s41365-018-0355-7 http://dx.doi.org/10.1007/s41365-018-0355-7 . 2018 Nuclear Science and Techniques 02 29