Billebaud, A.

Brissot, R.

Chaussonnet, P.

Destouches, C.

Heuer, D.

Kerveno, M.

Laurens, Jean-Marc

Le Brun, C.

Liatard, E.

Loiseaux, J.M.

Meplan, O.

Merle, E.

Perdu, F.

Vollaire, J.

HAL-IN2P3

The MUSE-4 experiment : prompt reactivity andneutron spectrum measurementsA. Billebaud, R. Brissot, D. Heuer, M. Kerveno, C. Le Brun, E. Liatard, J.M.Loiseaux, O. Meplan, E. Merle, F. Perdu, et al.To cite this version:A. Billebaud, R. Brissot, D. Heuer, M. Kerveno, C. Le Brun, et al.. The MUSE-4 experi-ment : prompt reactivity and neutron spectrum measurements. PHYSOR 2002 InternationalConference on the New Frontiers of Nuclear Technology Reactor Physics. Safety and High-Performance Computing, Oct 2002, Seoul, South Korea. American Nuclear Society, 13C-01-(9p.), 2002. <in2p3-00011957>HAL Id: in2p3-00011957http://hal.in2p3.fr/in2p3-00011957Submitted on 15 Oct 2002HAL is a multi-disciplinary open accessarchive for the deposit and dissemination of sci-entific research documents, whether they are pub-lished or not. The documents may come fromteaching and research institutions in France orabroad, or from public or private research centers.L’archive ouverte pluridisciplinaire HAL, estdestine´e au de´poˆt et a` la diffusion de documentsscientifiques de niveau recherche, publie´s ou non,e´manant des e´tablissements d’enseignement et derecherche franc¸ais ou e´trangers, des laboratoirespublics ou prive´s.PHYSOR 2002, Seoul, Korea, October 7-10, 2002 THE MUSE-4 EXPERIMENT: PROMPT REACTIVITY AND NEUTRON SPECTRUM MEASUREMENTS   A.Billebaud, R.Brissot, D.Heuer, M.Kerveno, C.LeBrun, E.Liatard,  J.-M.Loiseaux, O.Méplan, E.Merle, F.Perdu, and J.Vollaire. Institut des Sciences Nucléaires CNRS-IN2P3/UJF 53, avenue des Martyrs, F-38026 Grenoble Cedex, FRANCE billebaud@isn.in2p3.fr   C.Destouches, P.Chaussonnet,  and J.-M.Laurens. CEA Cadarache DER/SPEx/LPE F-13108 Saint-Paul-lez-Durance Cedex, FRANCE   On behalf of the MUSE European collaboration   ABSTRACT   In the case of the use of ADS for incineration of nuclear waste and energy production a better knowledge of the ways to monitor the different reactor parameters is required. To do so, the MASURCA fast neutron reactor in different subcritical configurations has been coupled to a deuteron accelerator producing a pulse of neutrons in the middle of the core. We performed both dynamic and spectroscopic measurements. Firstly, the time response of the core after a pulse is measured with a 235U fission chamber and a 3He proportional counter. We show that for a core close to criticality the neutron population behaves as predicted by the point kinetics theory. On the other hand, for subcriticality level relevant for ADS, the evolution of the neutron population is not a pure exponential and thus the point kinetics cannot be used to deduce the prompt multiplication factor. So we propose a new approach based partly on MCNP simulations and with less restrictive  assumptions than those of the point kinetics. This method allows a determination of different prompt k values which are in good agreement with the expected ones.  Spectroscopic measurements are also performed with a 3He proportional counter. From the specific response of the detector to monoenergetic neutrons it is possible to reconstruct the neutron energy spectrum at the detector location, which is found in good agreement with the simulated one below 0.6 MeV. Above this energy some studies are in progress to improve the detector response.  These very preliminary results will be completed by further tests and experiments planned till the end of 2003.     1. INTRODUCTION   In the field of the new interest for accelerator driven systems (ADS), new experiments are taking place at the MASURCA facility (CEA Cadarache, France) which consists of a small fast neutron reactor coupled to a deuteron accelerator (called GENEPI) creating a pulsed neutron source inside the reactor core. This reactor, operating at different subcritical levels allows the investigation of neutronic properties of such a medium and more particularly of its dynamic behaviour.  In ADS  two  multiplication factors can be defined. ks is the multiplication factor of  the source neutrons and is depending both on the source and the reactor caracteristics: it can be deduced from  1 PHYSOR 2002, Seoul, Korea, October 7-10, 2002 the total power of the reactor, but not so easily regarding to the tricky definition of a source neutron. The other factor which can be considered is the effective multiplication factor keff  corresponding to the multiplication  of the stabilized  fission distribution many generations after the source has been turned off. It is an intrinsic characteristic of the reactor and it governs its safety. Operating a subcritical reactor will require to monitor this keff  value. We explore here a way to determine the subcriticality level of the core from dynamic measurements performed in the MASURCA core submitted to source pulses and from related simulations. More precisely the measurement concerns the prompt multiplication factor kp as the βeff  value can be measured separately and is less essential here.  Finally we also show preliminary results concerning an attempt of neutron spectrometry in the reactor.    2. EXPERIMENTAL SET-UP   2.1 REACTOR   The reactor consists of about 21x23 modular vertical 4"x4" tubes. The central zone (~60 cm in diameter) is filled with MOX (UO2, PuO2) and Na rodlets, and is surrounded by a 25 cm Na and stainless steel reflector and a thick shield. The deuteron pulsed beam provided by the accelerator (250 kV) is brought in the middle of the core by mean of a horizontal glove finger ended by a TiD or TiT target (figure 1). The beam peak intensity is about 50 mA for a width of less than 1µs, and the repetition rate can vary from  a few Hz up to 5 kHz, providing (with the TiD target) about 4.104 neutrons per pulse (~4.107 n/s at 1 kHz). The core configuration for the measurements corresponds to keff = 0.994 with all control rods up (value obtained from MCNP simulations made with a detailed description of the reactor [1]), and when one rod is inserted it has been measured to 0.960±0.004 (MSA measurements [2]), the detailed composition of the control rod being not enough accurately known at the moment to perform the corresponding simulation.  Figure 1. MASURCA reactor  facility coupled to the pulsed neutron source "GENEPI".    2 PHYSOR 2002, Seoul, Korea, October 7-10, 2002 2.2 DETECTORS   Time and energy spectra are measured by mean of a fission chamber and a 3He proportional counter inserted in vertical 1" channels in the core. Signals are recorded with a 12 bit flash-ADC, working with a 10 MHz clock (100 ns/channel) and synchronized on the neutron pulsed source. The fission chamber has a U5 deposit (about 1 mg). The 3He detector is filled with about 10 mb of 3He and 6 b of Ar (and a few mb of CO2). For spectrometry in the fuel zone this detector is shielded from γ-rays with a 1" blanket of  Pb over the whole fuel zone length. Detector locations are shown in figure 2 on the transversal section (median plane).  Figure 2. Detector locations in the core: transversal section in the median plane.   3. DYNAMIC RESULTS   3.1 LIMIT OF THE POINT KINETICS MODEL  When a pulsed neutron source is injected into the core those neutrons initialise fission chains and thus the neutron flux is suddenly increased above its intrinsic value. After a sufficiently large number of collisions these neutrons represent an additional flux with a spectral and spatial distribution close to the flux of the already existing neutron population. One can assume that in those new chains each neutron creates kp neutrons in the next generation after an average time l  called the generation time. The value of l  can be estimated through a MCNP simulation: for a stabilized fission source l  is found equal to 0.58 µs. Under this assumption the decay of the neutron population follows a pure exponential decrease    N(t) = N0 exp (- αt)  (1)  with α = (1-kp)/l.  3 PHYSOR 2002, Seoul, Korea, October 7-10, 2002 The decrease of the neutron population can be obtained  by measuring the reaction rate of a detector located in the core. Figure 3 shows the logarithmic time spectra recorded with a fission chamber  in the fuel zone for subcriticality levels of 0.994 and 0.960 for a repetition rate of the pulsed source of 2 kHz. Are also shown the polynomial fits of the spectra.  When the reactor is close to the criticality (kef f= 0.994), after a few microseconds the curve behaves as a pure exponential with a slope α close to 0.016 µs-1. From this value we can calculate kp = 1-αl  = 0.990 ± 0.001 to be compared to the kp value also simulated kp = 0.991.  For a subcriticality level relevant for ADS, namely 0.960, shown in figure 3, the slope of the logarithmic decrease is clearly time dependent. This behaviour shows the limit of validity of the point kinetics model. 20 40 60 80 100 120 140012345678  Figure 3. Logarithmic time spectra obtained with a U5 fission chamber in the fuel zone of MASURCA for two subcriticality levels ke ff  = 0.994 and keff  = 0.960. Plain curves are polynomial fits of the spectra.  3.2 A MORE SOPHISTICATED MODEL     To explain the failure of  such a simple model for subcritical levels, we must pay attention to the distribution of neutron generations after several microseconds. When the multiplication factor is low, neutrons from the first generations become relatively more important than the ones of the later generations since these last generations induce less and less neutrons. Subsequently the lifetime of the early neutrons, much greater than the average generation time l , has to be taken into account. We propose a more sophisticated model where we do not assume anymore that each neutron creates another one at an average time l , but gives birth to a neutron in the next fission generation with a time distribution called P(τ), τ being the time elapsed since its own birth by fission. That distribution can be easily obtained by Monte-Carlo simulation for a stabilized neutron source and is shown in figure 4.             4 PHYSOR 2002, Seoul, Korea, October 7-10, 2002 0 20 40 60 8010-610-510-410-310-210-1110102100 Figure 4. Fission generation time distribution P(τ) provided by a MCNP simulation.  From that definition we deduce that     (2) pkdP =∫ ττ )( and thus we can normalise P(τ) to kp = 1. With that normalised distribution  P'(τ), we have access to the number of neutrons in the core at any time for any kp value, summing the contribution of each generation:  N kp (t) = kp P'(τ) + kp 2 P'(τ)*P'(τ) + kp 3 P'(τ)*P'(τ)*P'(τ) + ... (3) where * denotes the convolution operator. The decrease rate  α kp (t) can then be calculated for different kp values from the logarithmic derivative:   (4) dtdNNtkp1)( =α and compared to the α (t) obtained from the experimental N(t) spectrum, the one fitting the best the experiment giving the kp value of the reactor. This method has been applied to the spectra shown in figure 3 but also to the time spectrum measured with the 3He counter at keff = 0.960 shown in figure 5. The α (t) curves obtained from the fits of these spectra are shown together with relevant calculated αkp(t) curves in figure 6. With the U5 chamber the kp values of the core are found bounded by 0.990  and 0.993 for the assembly close to the criticality, and by 0.955 and 0.960 in the other case. They are  in rather good agreement with respectively the kp simulated value 0.991 and the one calculated with  kp = keff  (1 - βeff ) = 0.960(1 - 0.003) = 0.957.         5 PHYSOR 2002, Seoul, Korea, October 7-10, 2002 10 20 30 40 50 60 70 800246810 Figure 5. Logarithmic time spectrum obtained with a 3He counter in the fuel zone of MASURCA for keff  = 0.960. Plain curve is a polynomial fit of the spectrum.  The  α (t) obtained from the 3He counter exhibits a slight different behaviour compared to the U5 data. This can be explained by the location of the counters: the U5 chamber is in the middle of the fuel zone while the 3He counter is close to the reflector and thus submitted to more low energy neutrons. To improve measurements in such locations the use of detectors with energy threshold (like U8 chamber for instance) is required.  This method is very promising as it can give the kp value from a measurement made in the first tens of microseconds after the pulse, where the counting rate is still high. It is not necessary to reach an asymptotic behaviour of α(t), which may be not accessible to experiment. Moreover it is clearly seen that the α value given by the point kinetics model (0.073 µs-1) is not in agreement with the behaviour of the experimental α(t).     Figure 6. Logarithmic derivatives of the fits of the experimental time spectra compared to several α kp (t) calculated by the proposed method.  More details about this method and its robustness can be found in [3].   6 PHYSOR 2002, Seoul, Korea, October 7-10, 2002  4. NEUTRON SPECTROMETRY IN THE REACTOR   Neutron spectrometry in fast neutron reactors is not trivial considering spectra are continuous, ranging up to a few MeV and especially because the neutron flux is very large even in a small research reactor (about 106 n/cm2/s). Classical time of flight methods cannot be applied here essentially because neutron transport outside the core is most of the time impossible, and the use  of a time-energy relation between the neutron energy En and its tof is not as simple in a reactor as in an inert medium. Spectrometry based upon the detection of recoil protons from elastic scattering processes with hydrogen in plastic or organic scintillators can be done but it requires space which is sorely lacking in measurement channels. The recoil proton energies can be directly measured in a proportional gas counter. But a neutron of a given energy will induce a recoil proton energy distribution due to the recoil angle distribution. Moreover the presence of hydrogen in a fast reactor is not desirable. Finally, in proportional counters it is also possible to measure the energy of particles induced by neutron on light nuclei reactions, the particles sharing the energy released by the reaction Q and the incident neutron energy. We decided to use a 3He proportional counter whose energy calibration is easier due to the thermal neutron peak at Q = 764 keV. But the response of such counters to monoenergetic neutrons is not really ideal. The peak corresponding to the collection of the total kinetic energy of  the induced proton and triton is associated to a plateau due to wall effect in the counter and subsequent partial energy collection. This wall effect is reduced by increasing the counter pressure, but cannot be entirely suppressed for fast neutrons. In the response, one has to consider also the structure arising from a competing reaction occurring in the counter which is the elastic scattering of neutrons by 3He. The probability of this reaction becomes greater than that of the (n,p) reaction  above 100 keV, and gives a recoil 3He distribution ending at 3/4En which cannot be neglected above 1 MeV as both responses superimpose. The typical response of such a counter is illustrated in figure 7 for 1 MeV incident neutrons.   Figure 7. Typical response spectrum of 3He proportional counter to 1 MeV monoenergetic neutrons.  To extract a neutron energy spectrum from the deposited energy spectrum obtained with this counter one has to modelize the response as a function of neutron energy and to build the associated matrix M for the energy range of interest. To achieve this, spectra with monoenergetic neutrons have been obtained at a Van de Graaff accelerator (CEN-Bordeaux-Gradignan, France) from 200 keV up to 1.5 MeV. Detector responses have been fitted (and combined to simulations for recoil 3He distributions)  7 PHYSOR 2002, Seoul, Korea, October 7-10, 2002 and parametrized, the two reaction contributions being weighted according to their respective cross sections. The matrix Mij elements (i  for the neutron energy and j for the deposited energy) with a 10 keV energy binning are then calculated (by extrapolation above 1.5 MeV). Figure 7 represents Mi for i=100. The neutron spectrum E whose size is imax can be extracted from the deposited energy spectrum R whose size is jmax by an iterative subtraction method,   (R)1 = (R)0 - Mimax . Eimax (5)  where (R)0 is the experimental spectrum, (R)1 the spectrum obtained after the first iteration, and Eimax satisfies Mimax,jmax . Eimax = (Rjmax)1. At rank k we have:  (R)k = (R)k-1 - Mi . Ei (6)  with i= imax-k, j=jmax-k, and  Ei satisfying Mi,j . Ei = (Rj)k. This method has been applied to a measurement performed in the fuel zone for keff=0.97. The result, shown in figure 8 is compared to the neutron spectrum simulated with the Monte Carlo code MCNP-4C. The overall agreement is rather good up to 600 keV, which represents about 75% of the flux. However around  450 keV an underestimation is seen which is inherent to our subtraction method which fails to fully reproduce the flux depression due to the oxygen resonance around 400 keV: it induces an overestimation for a few following lower energy bins. The extraction  method is still being improved. We aim to be less dependent upon the flux structures and to be able to measure neutron spectra up to higher energies. This technique will allow us to validate of our simulations in a subcritical medium.   Figure 8. Neutron spectrum in the fuel zone simulated with MCNP for keff  = 0.97 compared to the spectrum measured with a 3He counter.   CONCLUSIONS   After having shown the failure of the point kinetics model to reproduce the kp value typical of an ADS system, we proposed a more sophisticated model to obtain kp  from experiment by introducing the generation time distribution P(τ). From this simulated distribution and successive convolutions it is possible to calculate the time evolution of the neutron population for any value of kp and to compare their logarithmic derivative to the experimental reactor response to  a pulsed source. With this method the measured kp values are found in good agreement with the expected ones. Especially it allows to  8 PHYSOR 2002, Seoul, Korea, October 7-10, 2002 obtain kp  in a time range where the decrease rate α(t) of the neutron population has not yet reached its asymptotic value. We also highlighted the importance of the choice of the detector in such measurements, which should be sensitive to fast neutrons, when its location  in the reactor is too close to the reflector. Concerning the neutron spectrometry, promising results are found with a 3He proportional counter and the reconstruction of its response in the major part of the spectrum (i.e. below 600 keV). A work is in progress to improve the flux extraction above this energy where the detector response has been extrapolated. These are very first results but all these techniques will  be tested and improved in several different experimental conditions in the immediate future during the following of the experimental MUSE programme planned till the end of 2003.   ACKNOWLEDGEMENTS    We would like to thank the GENEPI accelerator technical team from I.S.N. for its availability and help during the running of experiments and also the CEA MASURCA reactor team. One of us (M.K.) acknowledges a post-doctoral grant from the "C.Fréjacques-L.Michel Foundation" ("Société de Secours des Amis des Sciences"). This work is partially supported by the 5th framework programme of the European Commission.   REFERENCES  1. R.Soule, W.Assal, D.Villamarin, “Geometrical and physical data associated to the MUSE-4 reference configuration,” CEA-SPEX/LPE/NT-01-015, (2001).  2. R.Soule et al, CEA, Private communication, (2002).  3. F.Perdu, J.M.Loiseaux, A.Billebaud, R.Brissot, D.Heuer, C.LeBrun, E.Liatard, O.Méplan, E.Merle, H.Nifenecker, J.Vollaire, “Prompt reactivity determination in a subcritical assembly through the response to a dirac pulse,” Progress in Nuclear Energy, submitted, (2002).    9 

The MUSE-4 experiment : prompt reactivity and neutron spectrum measurements

Hal - Université Grenoble Alpes

The MUSE-4 experiment : prompt reactivity and neutron spectrum measurements

Abstract

Similar works

Full text

Available Versions

HAL-IN2P3

Hal - Université Grenoble Alpes