409 Search Results

Abstract not provided.

Neutron production can be realized with a high energy electron linear accelerator by using Bremsstrahlung and photoneutron converters. In this study, Monte Carlo N-Particle Code (MCNP) was used to evaluate potential photonuclear target designs for a high energy electron linear accelerator for applications such as neutron radiography and neutron resonance spectroscopy. A computational model was developed to inform a target design that would yield a high number of neutrons. It consists of a 20 MeV electron beam incident on a Bremsstrahlung target and a photonuclear target to generate neutrons. This computational model showed that a thickness of 0.75 inches for both tungsten and depleted uranium yields the most neutrons from photoneutron reactions. Saturation in the total number of generated neutrons was observed at over 0.75-inch thickness for both evaluated materials. Depleted uranium yielded approximately twice the number of neutrons overall compared to tungsten. The highest neutron surface flux for Depleted Uranium was 1.06 × 10-4 neutrons/cm2/source electron, and for Tungsten it was 5.12 × 10-5 neutrons/cm2/source electron. The optimal target design for this study’s application would consist of a 0.75 inch-thick block of depleted uranium with the length, width, and/or diameter varying dependent on application.

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.

In this study, we have extended the detector response function toolkit (DRiFT) to provide modeling capabilities of semiconductor sensors. DRiFT provides realistic nuclear instrumentation response by post-processing Monte-Carlo N-particle (MCNP®) radiation transport outputs. MCNP® is capable of modeling radiation transport in complex environments, but has limited detector physics and readout electronics modeling capabilities. Semiconductor detector response can be calculated with a high-fidelity for a flexible range of environments by utilizing MCNP® to simulate radiation interactions inside of detector volumes, and then using DRiFT to model charge transport and signal formation in the semiconductor, as well as the readout electronics. DRiFT models charge transport in the semiconductor, the preamplifier, shaping amplifier, pulse pile-up, and electronic noise to generate detector response. The semiconductor application in DRiFT can model a range of semiconductor materials, shapes, and sizes; and is demonstrated here for a large volume coaxial high-purity germanium (HPGe) detector. Here, we compare detector response functions of a coaxial HPGe detector with measurement of ⁶⁰Co, ¹³³Ba, and ¹³⁷Cs at varying count rates, and we conduct a parameter study to demonstrate the effect of changing parameters in the DRiFT simulation. The HPGe detector response function shows excellent agreement with measurements of difference sources with varying dead times and count rates.

The double-shell inertial confinement fusion campaign, which consists of an aluminum ablator, a foam cushion, a high-Z pusher (tungsten or molybdenum), and liquid deuterium–tritium (DT) fuel, aims for its first DT filled implosions on the National Ignition Facility (NIF) in 2024. The high-Z, high density pusher does not allow x-rays to escape the double-shell capsule. Therefore, nuclear diagnostics such as the Gamma Reaction History (GRH) diagnostic on the NIF are crucial for understanding high-Z implosion performance. To optimize the GRH measurement of fusion reaction history and the pusher’s areal density, the MCNP6.3-based forward model of the detector was built. When calculating the neutron-induced inelastic gamma ray production, the interaction of neutrons with the compressed fuel was additionally included. By folding the calculated gamma ray spectrum output and the previously calibrated GRH detector responses, the optimum set of GRH energy thresholds for measuring the pusher areal density is determined to be 2.9 and 6.3 MeV for DT double-shell experiments. In addition, the effect of the down-scattering of neutrons on the gamma ray spectrum, the minimum required yield for measurements, and the attenuation of the gamma rays through the pusher are analyzed.

The initial target design used for Pu-238 production at Idaho National Laboratory was designed by Oak Ridge National Laboratory to optimize the production of Pu-238 in the High Flux Isotope Reactor (HFIR) and are referred to as HFIR GEN II targets. To take advantage of the Advanced Test Reactor’s (ATR) taller active core region a redesign of the HFIR GEN II targets was needed. It was proposed to stack two HFIR GEN II targets nose to nose about the core center line; however, this resulted in excessive neutron and photon heating in the pellets located in the center. This peak heating was not desirable so three alternative designs were investigated for the ATR GEN I targets. The python-based code, MCNP to ORIGEN2 in Python (MOPY), was used to calculate the heating rates after 40 days of irradiation to capture the effects of each configuration. The purpose of this paper is to document the details of these conceptual design calculations and comparisons for the ATR GEN I targets.

This paper describes an experimental method to acquire high resolution energy- and spatially- resolved neutron beam spots using the time-gated neutron imaging system with Teledyne Pi-MAX4 camera. These experimental data offer a unique opportunity for benchmarking beam spot simulations. High-quality simulations depend significantly on a high-fidelity geometry model, which can be challenging for legacy facilities. We informed our MCNPX geometry model by latest metrology survey employing Leica laser tracker ATS600. It gave us a high fidelity description of our facility geometry. Such a robust integration of novel tools and methods yields a previously unattainable level of accuracy in both predicting and capturing neutron beam spots, marking a notable advancement over traditional methods reliant on static image plates. Here, to demonstrate the practical application of these tools, we are showing a non-uniform beam spot challenge at our Flight Path 14 (FP14) at the Los Alamos Neutron Science Center (LANSCE). Our precise MCNPX prediction of beam spot shifting as function of neutron energy was confirmed by experimental beam spot measured with extremely high level of detail. Results of this research demonstrate a significant leap in neutron beam optimization at LANSCE and set a new benchmark in beam spot characterization. The advanced methods presented here have potential for adoption at similar research facilities worldwide, aiming at substantial improvement in neutron beam delivery for experiments.

– Gas-filled neutron detectors have numerous applications across the nuclear engineering and nuclear physics fields. The ability to accurately model and simulate these detectors is important for those applications but is currently limited by the lack of readily-useable detector response software. Recently, the capabilities of DRiFT, a Detector Response Function Toolkit, were expanded to model gas-filled, He-3 and BF3, neutron detectors so that, combined with the radiation transport capabilities of the MCNP code, a high-fidelity treatment of gas-filled neutron detectors can be obtained. Further, this model has been validated by an experiment carried out with the Epithermal Neutron Multiplicity Counter and its capabilities have been demonstrated in two additional experiments. This work shows that utilizing DRiFT to post-process MCNP outputs produces more accurate results than using the MCNP code alone, reducing the difference between experimental and simulated results for measurements taken near the end of a He-3 tube, where the MCNP code struggles to model inactive regions of the detector, from a maximum of 35% with the MCNP code alone to 15% with the MCNP code plus DRiFT. DRiFT's diagnostic capabilities are also demonstrated with measurements for scenarios when pulse pileup or room return effects are significant and must be considered. Altogether, these measurements underpin the ability of DRiFT to accurately model and predict the behavior of gas-filled neutron detectors, making it a valuable tool for the design and testing of systems and experiments that utilize these detectors.

Standard MCNP particle tracking takes place along straight-line trajectories from interaction point to interaction point. There is a feature within MCNP that is planned for deprecation that provides surface boundary conditions for approximating gravity for planetary cases, but this feature is not applicable to a cold neutron beam. A new extension has been developed to track particles along parabolic trajectories with a constant acceleration. MCNP contains 1st and 2nd-order surfaces as well as a special case of 4th-order surfaces for simple tori, and the intersection of parabolic trajectories with these surfaces becomes 2nd, 4th, and 8th-order equations in time, respectively. Solving these equations utilizes a fast algorithm for finding the roots of polynomials. Finally, the theory, MCNP input card, and examples of using this new feature will be discussed.

MontePy is an open-source python software library for reading, editing, and writing MCNP input files. This presentation was given to a university research group interested in it. MontePy provides an object-oriented interface for working with these MCNP input files. It uses a lexer and parser system in order to create a concrete syntax tree representing the input from the input file.