Search NASA⌕ Search

SEARCH · Search NASA

Results for “MCNP Model”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 199 records · Page 11

Thermal Hydraulics Safety Analyses of TRIGA Reactor Upon Fuel Reconfiguration

Safety analyses of the Training, Research, Isotope, General Atomics (TRIGA) reactor at Kansas State University(KSU) is performed using thermal hydraulic simulations based on changes to reactivity or the volumetric heat generation rate due to changes in the positioning of the enriched fuel elements in the reactor core. If the reactivity or neutronics parameters are known, then they can be used to simulate the thermal energy state of the material which is crucial to assess the safety of the nuclear fuel inside a nuclear reactor. Several thermal hydraulics models and codes have been developed in the past to assess reactor safety. RELAP (Reactor Excursion and Leak Analysis Program) is one such tool widely used throughout the world for safety analyses. Although RELAP is considered a reliable tool by regulators, vendors, and operators, it lacks the numerical fidelity provided by 3-D computational fluid dynamics (CFD) models. This work involves the critical task of assessing safety of a research reactor upon fuel reconfiguration based on models developed in RELAP and comparing them with a commercially available 3-D CFD program. The TRIGA Mark-II reactor, designed by General Atomics and located at Kansas State University, is a class of small pool type inherently safe research reactor licensed to operate at steady state power of up to 1.25 MW. It has a cylindrical core measuring 0.45 meters diametrically and 0.38 meters axially positioned in a tank measuring 1.98 meters diametrically and 6.25 meters axially. It has uranium zirconium hydride (UZrH) fuel rods housed in stainless steel cladding material arranged in five concentric rings (B,C,D,E and F rings, in that order, moving outwards from the center). TRIGA fuel is nominally 20% enriched, but the reactivity of the fuel can be changed by using burnable poisons (not used at KSU), or by varying the total weight percent of uranium in the UZrH fuel meat. The TRIGA reactor is designed to be an inherently safe reactor and is designed to be cooled by natural convection. Thus it is of paramount importance to avoid the film boiling regime; failing to do so will subject the fuel to severe conditions, possibly resulting in its failure due to fuel/clad interactions. A high DNBR ratio is a primary consideration in the reactor core design and thus any modification in the core should ensure that this criteria is always met. In order to fully utilize the capabilities of KSU TRIGA reactor facility, it should be operated at its full capacity, i.e. 1.25 MW, which requires increasing the fuel loading in some of the E- and F-ring elements from 8.5% U (nominal) to 12.0% U (nominal). Only four elements are to be switched to 12% loading, causing local power peaking in these elements. The degree of power peaking to be expected in these 12%-loaded fuel elements was calculated based on a detailed MCNP (Monte Carlo N-Particle code) model of the reactor core. Some preliminary results of the calculations are presented in this paper based on changes in peaking factors and power levels. To predict reactor thermal behavior, axial changes in the peaking factors obtained from MCNP models were used to calculate volumetric energy density functions for a single channel one dimensional RELAP model and a three dimensional COMSOL model of the TRIGA MARK-II reactor. (authors)

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Processing MCNP Elemental Edit Outputs

The Monte Carlo N-Particle (MCNP) transport code version 6 (also known as MCNP6) has the capability for tracking particles on unstructured mesh (UM) geometry models embedded into constructive solid geometry (CSG) cells. A UM geometry is a collection of elements representing a solid geometry. The first step of MCNP UM modeling is using other software packages to create a finite element mesh representation of a solid 3D geometry. Computer-aided design (CAD) or computer-aided manufacturing (CAM) software is typically used to create a solid geometry model, which is later imported into mesh generation software to create a UM model. The MCNP UM feature was originally designed for models generated by the Abaqus/CAE software. The MCNP code version 6.0 and later can process UM models formatted as Abaqus input files. MCNP can process a UM model consisting of several different element types including linear tetrahedral or hexahedral elements and calculate quantities of interest such as flux and energy deposition at elements. An MCNP UM simulation provides high-fidelity elemental edit (i.e., tally) outputs, which can be further used in multiphysics calculations. The MCNP UM feature was used for multiphysics simulations where quantities of interest calculated by MCNP are used as inputs for heat transfer calculations in Abaqus. MCNP6.3 can produce two types of elemental edit output (EEOUT) file formats: ASCII and HDF5. An EEOUT file type must be requested on an EMBED card while output type (flux or energy deposition) must be requested on an EMBEE card. We wrote Python3 scripts to extract energy deposition values in an ASCII or HDF5 EEOUT file and compute a heat flux profile for an Abaqus heat transfer calculation.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

As-Run ATR Physics Evaluation for the UCF-1 HSIS Experiment

This Engineering Calculations and Analysis Report (ECAR) documents the Advanced Test Reactor (ATR) detailed MCNP full core model physics analysis performed to support the irradiation of the University of Central Florida (UCF-1) Hydraulic Shuttle Irradiation System Experiment in the ATR B-7 position. The purpose of this analysis is to calculate: • The as-run neutron and gamma heating rates for the experiment components (W/g); • The as-run DPA for the experiment specimens after 20 (test train 1) and 2 (test train 2) days of irradiation in the ATR B-7 position (neutrons/cm2-sec); • The radionuclide source term after 2 days and 20 days of irradiation in the ATR B-7 position; • The decay heat associated with the UCF-1 radionuclide source term. These calculations were performed using the computer codes MCNP and ORIGEN2. The neutronics model, analysis details, and results are presented in the attached report.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Characterization of Uranium Foil Irradiations at the WSU TRIGA Reactor Using a New Reactor Model in SCALE [Slides]

A depleted uranium foil and a highly enriched uranium foil that were irradiated at the Washington State University TRIGA reactor were used to measure the relative production of 237 U and 239 U, two key isotopes of interest in nuclear forensics, as well as their decay daughters 237 Np and 239 Pu. Total fissions in the two foils were also estimated using 99 Mo measurement. A reactor model of the TRIGA was built in the SCALE neutron transport code using dimensions obtained from an existing MCNP geometry. The model was successfully used to simulate foil fissions and the reaction rate ratios in 238 U, using the laboratory-measured isotope ratios as a benchmark. The flux distribution in the foils was determined. Finally, using the continuous energy Tsunami-3d module of SCALE, a sensitivity/uncertainty analysis was performed to determine the effects on the reaction ratio caused by other reactions in the reactor and by uncertainties in the known cross-sections.

07 ISOTOPE AND RADIATION SOURCES↗

Methodology to Generate the Decay Gamma Source for a Second Target Station Target Wedge

This report details the specifics of applying the Position-Averaged Method methodology for the calculation of the rotating lasagna-style target wedge activation and the resulting decay gamma source terms. The methodology uses MCNP® Code Version 6.2.0 with the RNUCS patch coupled with CINDER2008 from the AARE V1.0 package. The process begins with model preparations that include ensuring impurities are included in the materials in the regions of interest along with segmentation of the geometry in order to capture the decay gamma source gradients. This report also details the physics models and other MCNP®, CINDER2008, and ADVANTG 3.2.0 related options used to calculate the decay gamma source terms. ADVANTG 3.2.0 is used to generate a weight window set to aid in convergence of the neutron fluxes in the initial MCNP® protons-on-target particle transport calculation to calculate the neutron fluxes and spallation products in the regions of interest. The neutron fluxes and spallation products are then used in the CINDER2008 transmutation calculation, and the decay gamma sources are generated by the AARE_GAMMA_SOURCE_SCRIPT as a part of the AARE V1.0 package. These decay gamma sources include the energy distributions as well as the necessary distributions to sample the location of the source using the cell rejection method. The gamma sources that are produced using the methodology may be used for a wide variety of analysis involving the target wedge. This report includes the material specifications, physics options, proton source term description, geometrical configuration, and post-processing details.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Monte Carlo N-Particle forward modeling for density reconstruction of double shell capsule radiographs

In the Double Shell Inertial Confinement Fusion concept, characterizing the shape asymmetry of imploding metal shells is vital for understanding energy-efficient compression and radiative losses of the thermonuclear fuel. The Monte Carlo N-Particle MCNP® code forward models radiography of Double Shell capsule implosions using the Advanced Radiographic Capability at the National Ignition Facility. A procedure is developed for using MCNP to reconstruct density profiles from the radiograph image intensity. For a given Double Shell imploding target geometry, MCNP radiographs predict image contrast, which can help guide experimental design. Finally, in future work, the calculated MCNP synthetic radiographs will be compared with experimental radiographs to determine the radial and azimuthal density profiles of the Double Shell capsules.

47 OTHER INSTRUMENTATION↗

TH-CD-207A-08: Simulated Real-Time Image Guidance for Lung SBRT Patients Using Scatter Imaging

Purpose: To develop a comprehensive Monte Carlo-based model for the acquisition of scatter images of patient anatomy in real-time, during lung SBRT treatment. Methods: During SBRT treatment, images of patient anatomy can be acquired from scattered radiation. To rigorously examine the utility of scatter images for image guidance, a model is developed using MCNP code to simulate scatter images of phantoms and lung cancer patients. The model is validated by comparing experimental and simulated images of phantoms of different complexity. The differentiation between tissue types is investigated by imaging objects of known compositions (water, lung, and bone equivalent). A lung tumor phantom, simulating materials and geometry encountered during lung SBRT treatments, is used to investigate image noise properties for various quantities of delivered radiation (monitor units(MU)). Patient scatter images are simulated using the validated simulation model. 4DCT patient data is converted to an MCNP input geometry accounting for different tissue composition and densities. Lung tumor phantom images acquired with decreasing imaging time (decreasing MU) are used to model the expected noise amplitude in patient scatter images, producing realistic simulated patient scatter images with varying temporal resolution. Results: Image intensity in simulated and experimental scatter images of tissue equivalent objects (water, lung, bone) match within the uncertainty (∼3%). Lung tumor phantom images agree as well. Specifically, tumor-to-lung contrast matches within the uncertainty. The addition of random noise approximating quantum noise in experimental images to simulated patient images shows that scatter images of lung tumors can provide images in as fast as 0.5 seconds with CNR∼2.7. Conclusions: A scatter imaging simulation model is developed and validated using experimental phantom scatter images. Following validation, lung cancer patient scatter images are simulated. These simulated patient images demonstrate the clinical utility of scatter imaging for real-time tumor tracking during lung SBRT.

60 APPLIED LIFE SCIENCES↗

Study of Whisper Statistics for Criticality Safety Code Validation

Whisper is computational software designed to assist with validation studies with the Monte Carlo particle transport package MCNP. Historical approaches rely on expert based judgment to select benchmarks while Whisper uses sensitivity/uncertainty (S/U) methods to select relevant benchmarks to a particular application model. Whisper will also derive a calculational margin and margin of subcriticality (MOS) that may be used to generate an upper subcritical limit (USL) for a particular application model. Whisper finds benchmark experiments whose isotope sensitivity profiles align best with the MCNP model in question, and assign a correlation coefficient ('ck' in Whisper parlance) to each selected benchmark as a measure of how closely those profiles agree. When ck values are high, there is less uncertainty between model and benchmarks and a higher USL is possible, relative to cases when ck values are low. The default Whisper settings recommends a 100,000 particles per cycle run for 500 active cycles, resulting in 50 M active neutron histories per application run. This is computationally intensive and may not be necessary to produce accurate results. A Whisper statistical study was done varying the neutron histories per cycle from 10,000 to the default 100,000. A lower number of neutron histories will allow Whisper to be used on low memory systems such as laptops. The intent of this paper is to present the changes in the Whisper results as the number of neutrons histories run is varied.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

AGR 1 and 2 Isotopic Depletion Validation

This presentation includes work performed developing a computational benchmark for the AGR 1 and 2 experiments. These experiments involve irradiation of TRISO particle compact in the Advanced Test Reactor (ATR). Many post-irradiation measurements were made from these experiments, including isotopic concentrations of fission products. The benchmark involves a MCNP-ORIGEN coupled model of the ATR and AGR experiment capsules to calculate a comprehensive inventory of fission product and actinide concentrations in the compacts. This presentation describes the experiments, the measurements obtained, the benchmark model, and results of the model evaluation.

97 - MATHEMATICS AND COMPUTING↗

Evolving Adjustments to External (Gamma) Slope Factors for CERCLA Risk and Dose Assessments: The MCNP Years - 16307

To model the external exposure pathway in risk and dose assessments of radioactive contamination at Superfund sites, the U.S. Environmental Protection Agency (EPA) uses slope factors (SFs), also known as risk coefficients, and dose conversion factors (DCFs). Without any adjustment these external radiation exposure pathways effectively assumes that an individual is exposed to a source geometry that is effectively an infinite slab. The concept of an 'infinite slab' means that the thickness of the contaminated zone and its aerial extent are so large that it behaves as if it were infinite in its physical dimensions. EPA has been making increasingly complex adjustments to account for the extent of the contamination and its corresponding radiation field to provide more accurate risk and dose assessment modeling when using its calculators. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

MO-FG-CAMPUS-TeP3-02: Benchmarks of a Proton Relative Biological Effectiveness (RBE) Model for DNA Double Strand Break (DSB) Induction in the FLUKA, MCNP, TOPAS, and RayStation™ Treatment Planning System

Purpose: Validate implementation of a published RBE model for DSB induction (RBEDSB) in several general purpose Monte Carlo (MC) code systems and the RayStation™ treatment planning system (TPS). For protons and other light ions, DSB induction is a critical initiating molecular event that correlates well with the RBE for cell survival. Methods: An efficient algorithm to incorporate information on proton and light ion RBEDSB from the independently tested Monte Carlo Damage Simulation (MCDS) has now been integrated into MCNP (Stewart et al. PMB 60, 8249–8274, 2015), FLUKA, TOPAS and a research build of the RayStation™ TPS. To cross-validate the RBEDSB model implementation LET distributions, depth-dose and lateral (dose and RBEDSB) profiles for monodirectional monoenergetic (100 to 200 MeV) protons incident on a water phantom are compared. The effects of recoil and secondary ion production ({sub 2}H{sub +}, {sub 3}H{sub +}, {sub 3}He{sub 2+}, {sub 4}He{sub 2+}), spot size (3 and 10 mm), and transport physics on beam profiles and RBEDSB are examined. Results: Depth-dose and RBEDSB profiles among all of the MC models are in excellent agreement using a 1 mm distance criterion (width of a voxel). For a 100 MeV proton beam (10 mm spot), RBEDSB = 1.2 ± 0.03 (− 2–3%) at the tip of the Bragg peak and increases to 1.59 ± 0.3 two mm distal to the Bragg peak. RBEDSB tends to decrease as the kinetic energy of the incident proton increases. Conclusion: The model for proton RBEDSB has been accurately implemented into FLUKA, MCNP, TOPAS and the RayStation™TPS. The transport of secondary light ions (Z > 1) has a significant impact on RBEDSB, especially distal to the Bragg peak, although light ions have a small effect on (dosexRBEDSB) profiles. The ability to incorporate spatial variations in proton RBE within a TPS creates new opportunities to individualize treatment plans and increase the therapeutic ratio. Dr. Erik Traneus is employed full-time as a Research Scientist at RaySearch Laboratories. The research build of the RayStation used in the study was made available to the University of Washington free of charge. RaySearch Laboratories did not provide any monetary support for the reported studies.

60 APPLIED LIFE SCIENCES↗

MCNP6.3 Unstructured Mesh Verification: GodivR and CANDU Models

A geometric cell of the Monte Carlo N-Particle (MCNP)1 transport code is traditionally created by using Boolean operators on defined surfaces. This constructive solid geometry (CSG) capability has been available in the MCNP code since its beginning. However, a CSG model approach is limited when it comes to constructing a representative geometry for a complex model in its ability to capture a correct model representation. Starting with the version 6.0, the MCNP code has the ability of embedding an unstructured mesh (UM) model into a CSG cell to create a hybrid geometry [1]. The MCNP UM feature provides the flexibility of defining very complex geometries because computer aided design (CAD) and mesh generation software packages can be utilized to construct UM models for MCNP simulations.

97 MATHEMATICS AND COMPUTING↗

ACRRF High-Bay Dose Calculations using MCNP (Part A)

Analytical tools and models have been developed as a starting point for directly assessing dose in the Annular Core Research Reactor Facility (ACRRF) due to reactor operation. Key results include peak dose along the Central Cavity (CC) centerline (beamline) at the cavity level, dose throughout the High-Bay (HB), and dose on the facility roof for partially-shielded reactor operation where the 4” insert is removed from the CC Shield Plug (SP). Model results in the beamline are benchmarked against measured doses from passive dosimetry evaluations. Personnel total (neutron and gamma) dose in the ACRRF HB is calculated using Monte Carlo N-Particle (MCNP). Various CC and SP configurations are analyzed, including unshielded (no SP) and partially shielded (SP installed but 4” insert removed). Novel application of Variance Reduction (VR) techniques, namely the Surface Source Write (SSW) and Surface Source Read (SSR) capabilities in MCNP, enable impressive resolution (in a Monte Carlo modeling sense) of dose throughout much the facility. The VR techniques reduce stochastic error for challenging tallies, with more advanced techniques explored in the companion to this report (Part B) [1]. Supplementary studies (including a verification analysis) and pedagogic evaluations in Part B involve neutron spectra, angular distributions, and the dose impact of facility characteristics. With the SP 4” insert removed and the Lead-Boron (44”) Bucket (LB–44) in the reactor cavity, Total Effective Dose (TED) within the CC beamline is ≈140 rem per 300 MJ of reactor yield (or 3900 rem per hour at 100% Steady-State (SS) power). With no SP (unshielded) and a Free-Field (FF) cavity, TED within the beamline is ≈610 rem per 300 MJ (or 17000 rem per hour at 100% SS power). Due to the predicted collimation of radiation by the reactor pool (and partial SP, if present), beamline dose is much greater than the scattered radiation field surrounding the cavity and reactor tank. Comparisons are made to beamline dosimetry measurements to validate the model. Model predictions agree reasonably well (⪅10%) with measured quantities of neutron fluence, gamma fluence, and spectral metrics. Away from the beamline, comparisons made to previous dose measurements in the HB agreement within an order of magnitude.

61 RADIATION PROTECTION AND DOSIMETRY↗

MCNP6 and SERPENT2 Performances with Unstructured Mesh Geometry

Traditionally, Monte Carlo codes model a nuclear reactor by combinatorial geometry. In this model, volume primitives (e,g, spheres, cylinders, parallelepipeds...) are combined together by Boolean operations (e.g. union, intersection, and complement) and replicated by lattice definitions to build the core geometry. Recently, high performances computing platforms allowed using unstructured mesh geometry in Monte Carlo codes. When this option is chosen, any type of core geometry can be modeled since the restriction of using volume primitives is removed. In an unstructured mesh geometry, Boolean operations and lattice definitions are not used since each volume is defined by triangular surfaces (facets). This study focuses on two Monte Carlo codes: MCNP (version 6.1.1 beta) developed at Los Alamos National Laboratory since 1957 and written in FORTRAN90, C, and C++ (267,000 total source code lines), and SERPENT (version 2.1.27) developed at VTT Technical Research Centre since 2004 and written in C (195,000 total source code lines). These two codes follow different approaches for the unstructured mesh geometry modeling: MCNP uses ABAQUS format, whereas SERPENT uses the STL format. In the STL (STereoLithographic) geometry model only the surfaces of the volume are meshed; the interior part of the volumes is not meshed. In the STL geometry model, the whole geometry is defined by the sequence of SOLID and ENDSOLID sections. Each volume of the geometry must have a SOLID and ENDSOLID section. The STL file contains no information specifying the materials of the volumes. MCNP and SERPENT can model any kind of complicated geometry using unstructured meshes. MCNP and SERPENT use the ABAQUS and STL file formats, respectively. The STL file format does not have tetrahedrons and has only triangular facets on the volumes boundary. This allows SERPENT simulations to run 10 times faster than MCNP. However, if Monte Carlo results have to be coupled to thermal-hydraulics codes, the interior of the geometry volumes must be meshed, which may reduce SERPENT performances.

97 MATHEMATICS AND COMPUTING↗

Development of MCNP Training Modules for Safeguards Practitioners [Abstract]

The Monte-Carlo N-Particle (MCNP) software developed at LANL is the most widely used neutron transport code in the world. It is an essential tool for a variety of applications including detector development and design, nuclear fuel burnup simulation, criticality safety, and nondestructive assay system optimization. For this reason, it is indispensable within the safeguards and materials control & accountability (MC&A) communities. Multiple MCNP training courses have been created and taught over the last several decades by the MCNP development team at LANL, however there are no existing courses that cover specialized topics considered fundamental to NDA and safeguards models. To fill this gap, the MCNP team and Safeguards Science and Technology group at LANL have co-created a set of training modules customized to meet the specialized needs of the safeguards and MC&A communities. The basic modules cover concepts such as NDA system optimization, He-specific and other capture tallies, and tools for improved theoretical understanding. An advanced module was also created to cover topics including variance reduction for active interrogation simulations, use of the LANL MCNPTools post-processor, PTRAC (particle tracking) and list-mode data simulations, and fuel burnup simulations. The training modules teach to the latest and most state-of-the-art MCNP features and tools released by the development team at LANL and are intended to be taught jointly by the developers and safeguards experts. Ultimately, we hope that creation of these modules will serve to capture and convey the safeguards modeling and MCNP expertise at LANL, and that we will be able to share the modules more broadly with the MC&A and safeguards communities.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Zero Power Reactor Database (ZPRD) Development Plan

Past sodium-cooled fast reactors (SFR) were built with an active experimental program in place to support the design and development work. Most of the experimental facilities in the United States that were important for SFR design were shutdown in the 1980s and 1990s. Reactor licensing and construction requires any reactor design to be verified against existing reactor facilities or experimental measurements. With the absence of those experimental facilities, modern SFR projects must rely on historical measurements to demonstrate that the engineering modeling software and data being used for the new reactor design work are reliable. There has been a considerable push in the last 6 years by both DOE and commercial companies to obtain historical experimental measurements that are relevant for SFRs, in particular those with features that are important for the new reactor designs of interest. The zero power reactor experiments carried out at Argonne National Laboratory’s critical facilities (ZPR-3, ZPR-6, ZPR-9, and ZPPR) from the 1950s to the 1980s are some of the best reactor physics experiments on SFR technology that are available today. Of particular interest today are the ZPPR-15 measurements done at the ZPPR facility for the Integral Fast Reactor project in the 1980s as they are in line with most commercial and DOE interests today. In the past 10 years, the measurements done on ZPPR-15 have been processed into both Monte Carlo (MCNP) and deterministic models (MC2-3 and DIF3D) useable for validating the engineering modeling software for key parts of the SFR design work. To achieve this, a detailed model description must be created for the experiment and the experimental measurement that the engineering modeling software is to reproduce. Then, an assessment of the uncertainty on the measured quantity which considers all of the sources of uncertainty in defining the model must be obtained and documented. The models created for ZPPR-15 provide the best validation basis available today for neutronics modeling software. Reference 2 is a good resource to understand how these models were built and how the uncertainties on the measured quantities were derived. The intention of the Zero Power Reactor Database (ZPRD), hosted at frdb.ne.anl.gov, is to make available the experimental measurements and models that have been constructed to-date. Though ZPPR-15 measurements are the primary data requested for validation needs, other measurements on ZPPR, ZPR-6, and ZPR-9 in support of the Clinch River Breeder Reactor (CRBR) and Fast Test Reactor (FFTF) should also be considered important for future software validation needs. In this manuscript, the details of available measurements on ZPR-3, ZPR-6, ZPR-9, and ZPPR facilities are summarized, and a general organization of the web interface is displayed. Many of the documents associated with the measurements are export controlled information so access to the database will also have to be controlled.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Development, Testing and Validation of a Waste Assay System for the Measurement and Characterisation of Active Spent Fuel Element Debris From UK Magnox Reactors - 12533

This paper describes the development, testing and validation of a waste measurement instrument for characterising active remote handled radioactive waste arising from the operation of Magnox reactors in the United Kingdom. Following operation in UK Magnox gas cooled reactors and a subsequent period of cooling, parts of the magnesium-aluminium alloy cladding were removed from spent fuel and the uranium fuel rods with the remaining cladding were removed to Sellafield for treatment. The resultant Magnox based spent fuel element debris (FED), which constitutes active intermediate level waste (ILW) has been stored in concrete vaults at the reactor sites. As part of the decommissioning of the FED vaults the FED must be removed, measured and characterised and placed in intermediate storage containers. The present system was developed for use at the Trawsfynydd nuclear power station (NPS), which is in the decommissioning phase, but the approach is potentially applicable to FED characterisation at all of the Magnox reactors. The measurement system consists of a heavily shielded and collimated high purity Germanium (HPGe) detector with electromechanical cooling and a high count-rate preamplifier and digital multichannel pulse height analyser. The HPGe based detector system is controlled by a software code, which stores the measurement result and allows a comprehensive analysis of the measured FED data. Fuel element debris is removed from the vault and placed on a tray to a uniform depth of typically 10 cm for measurement. The tray is positioned approximately 1.2 meters above the detector which views the FED through a tungsten collimator with an inverted pyramid shape. At other Magnox sites the positions may be reversed with the shielded and collimated HPGe detector located above the tray on which the FED is measured. A comprehensive Monte Carlo modelling and analysis of the measurement process has been performed in order to optimise the measurement geometry and eliminate interferences from radioactive sources and FED in the immediate vicinity of the measurement position. The detector system has been calibrated and high activity radioactive sources of Cs-137, Co-60 and Na-22 have been used to validate the measurement process. The data acquisition and analysis software code has been tested and validated in keeping with the software quality assurance requirements of both ISO:9001-2008 - TICK-IT in the UK and NQA-1. The measurement and analysis system has been comprehensively tested with high activity sources, is flexible and may be applicable to a wide range of remote handled radioactive waste measurement applications. It is due to be installed at Trawsfynydd NPS later this year. This paper describes the Waste Tray Assay System (WTAS) that has been developed for the measurement of Magnox FED waste. The WTAS has been tested with a range of radioactive sources and its operation has been simulated with benchmarked MCNP Monte Carlo calculations. The measurement software has been validated as has the operation of the system for a range of strong radioactive sources. A system based on the design is due for installation and operation in 2012. The system has application to the measurement of Magnox Fuel Element Debris (FED) waste at other Magnox reactor sites. The major design objective of the WTAS that has been achieved is the ability of the assay system to determine the content of Cs-137, and in turn to enable the fissile burden to be assessed using a radionuclide fingerprint, in the presence of higher and highly variable quantities of Co-60, typically from nimonic springs. The approach can be used in other Magnox FED waste configurations where the detector is located above the FED waste sorting tray and where the collimation is fixed below the detector and at a distance above the tray. In this case, which has also been investigated, there are different shielding problems and mechanical support issues. The extensive use of MCNP Monte Carlo modelling to simulate the geometry of the sorting cell and the distribution of radioactive sources has helped to ensure that all of the detector shielding requirements are addressed and suitable Cs-137 and Co-60 discrimination can be achieved. The WTAS in its present form or in other configurations has relevance to the measurement of other active ILW and highly active RH waste. Examples include high activity RH LLW and RH TRU (Transuranic) waste as defined in the United States arising from both commercial nuclear and Department of Energy (DOE) operations. The analysis is able to analyse a range of radionuclides beyong those expected in the Magnox FED cases. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗

A novel method for modeling the neutron time of flight (nTOF) detector response in current mode to inertial confinement fusion experiments.

There are several machines in this country that produce short bursts of neutrons for various applications. A few examples are the Z machine, operated by Sandia National Laboratories in Albuquerque, NM ; the OMEGA Laser Facility at the University of Rochester in Rochester, NY ; and the National Ignition Facility (NIF) operated by the Department of Energy at Lawrence Livermore National Laboratory in Livermore, California . They all incorporate neutron time of flight (nTOF) detectors which measure neutron yield, and the shapes of the waveforms from these detectors contain germane information about the plasma conditions that produce the neutrons. However, the signals can also be %E2%80%9Cclouded%E2%80%9D by a certain fraction of neutrons that scatter off structural components and also arrive at the detectors, thereby making analysis of the plasma conditions more difficult. These detectors operate in current mode - i.e., they have no discrimination, and all the photomultiplier anode charges are integrated rather than counted individually as they are in single event counting. Up to now, there has not been a method for modeling an nTOF detector operating in current mode. MCNP PoliMi was developed in 2002 to simulate neutron and gamma ray detection in a plastic scintillator, which produces a collision data output table about each neutron and photon interaction occurring within the scintillator; however, the post processing code which accompanies MCNP PoliMi assumes a detector operating in single event counting mode and not current mode. Therefore, the idea for this work had been born: could a new post processing code be written to simulate an nTOF detector operating in current mode? And if so, could this process be used to address such issues as the impact of neutron scattering on the primary signal? Also, could it possibly even identify sources of scattering (i.e., structural materials) that could be removed or modified to produce %E2%80%9Ccleaner%E2%80%9D neutron signals? This process was first developed and then applied to the axial neutron time of flight detectors at the Z Facility mentioned above. First, MCNP PoliMi was used to model relevant portions of the facility between the source and the detector locations. To obtain useful statistics, variance reduction was utilized. Then, the resulting collision output table produced by MCNP PoliMi was further analyzed by a MATLAB postprocessing code. This converted the energy deposited by neutron and photon interactions in the plastic scintillator (i.e., nTOF detector) into light output, in units of MeVee%D1%84 (electron equivalent) vs time. The time response of the detector was then folded into the signal via another MATLAB code. The simulated response was then compared with experimental data and shown to be in good agreement. To address the issue of neutron scattering, an %E2%80%9CIdeal Case,%E2%80%9D (i.e., a plastic scintillator was placed at the same distance from the source for each detector location) with no structural components in the problem. This was done to produce as %E2%80%9Cpure%E2%80%9D a neutron signal as possible. The simulated waveform from this %E2%80%9CIdeal Case%E2%80%9D was then compared with the simulated data from the %E2%80%9CFull Scale%E2%80%9D geometry (i.e., the detector at the same location, but with all the structural materials now included). The %E2%80%9CIdeal Case%E2%80%9D was subtracted from the %E2%80%9CFull Scale%E2%80%9D geometry case, and this was determined to be the contribution due to scattering. The time response was deconvolved out of the empirical data, and the contribution due to scattering was then subtracted out of it. A transformation was then made from dN/dt to dN/dE to obtain neutron spectra at two different detector locations.

Nelson, Alan J.↗