Modelling of the anti-neutrino production and spectra from a Magnox reactor - 57134

Full text of publication follows: The anti-neutrino source properties of a fission reactor are governed by the production and beta decay of the radionuclides present and the summation of their individual anti-neutrino spectra. The fission product radionuclide production changes during reactor operation and different fissioning species give rise to different product distributions. It is thus possible to determine some details of reactor operation, such as power, from the anti-neutrino emission to confirm safeguards records. Also according to some published calculations, it may be feasible to observe different anti-neutrino spectra depending on the fissile contents of the reactor fuel and thus determine the reactor's fissile material inventory during operation which could considerable improve safeguards. In mid-2014 the University of Liverpool deployed a prototype anti-neutrino detector at the Wylfa R1 station in Anglesey, United Kingdom based upon plastic scintillator technology developed for the T2K project. The deployment was used to develop the detector electronics and software until the reactor was finally shutdown in December 2015. To support the development of this detector technology for reactor monitoring and to understand its capabilities, the National Nuclear Laboratory modelled this graphite moderated and natural uranium fuelled reactor with existing codes used to support Magnox reactor operations and waste management. The 3D multi-physics code PANTHER was used to determine the individual powers of each fuel element (8*6152) during the year and a half period of monitoring based upon reactor records. The WIMS/TRAIL/FISPIN code route was then used to determine the radionuclide inventory of each nuclide on a daily basis in each element. These nuclide inventories were then used with the BTSPEC code to determine the anti-neutrino spectra and source strength using JEFF-3.1.1 data. Finally the anti-neutrino flux at the detector was determined, and using the anti-neutrino cross-section on hydrogen a preliminary estimate of the expected reaction rate in the detector for each day during the year and a half of monitored reactor operation calculated. The results of the preliminary calculations are shown and limitations in the methods and data discussed including the use of different methods to calculate the anti-neutrino spectra such as the recent BetaShape code. (authors)