Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “mcnp”

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 307 records · Page 17

Utilization of ACE nuclear data file toolkit ACEtk to calculate relative sensitivity coefficients of point-kinetics parameters

Sensitivity and uncertainty methods are quintessential for nuclear criticality safety and experiment design. This type of analysis relies on calculations of sensitivity coefficients; sensitivity coefficients of the effective neutron multiplication factor with respect to some nuclear data are predominantly calculated and used. As a part of the Laboratory Directed Research & Development project EUCLID (Experiments Underpinned by Computational Learning for Improvements in nuclear Data) at Los Alamos National Laboratory, sensitivity coefficients of many radiation detector measurement responses with respect to nuclear data were investigated. Specifically, this paper outlines a method to calculate point-kinetics parameters relative sensitivity coefficients with respect to nuclear data. Point-kinetics parameters such as the prompt neutron decay constant, effective delayed neutron fraction, and neutron generation time are especially important to experimenters and reactor operators designing systems with dynamic neutron populations. This method couples capabilities of the ACE (A Compact ENDF) nuclear data file toolkit, ACEtk, with the ability to load cross sections into the radiation transport code Monte Carlo N-Particle (MCNP). In conclusion, key aspects of optimizing this method for a particular application and sensitivity profiles of the Jezebel criticality experiment are examined and discussed.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

Investigation of irradiation damage and heat deposition: a comparative analysis for HEU-to-LEU conversion in HFIR

The planned conversion of the High Flux Isotope Reactor (HFIR) at Oak Ridge National Laboratory from highly enriched uranium (HEU) to low-enriched uranium (LEU) fuel requires detailed evaluation of experiment-relevant parameters to ensure continued performance for materials testing and isotope production. Here, this study presents the first comprehensive assessment of displacements per atom (dpa) and heat deposition rates in target materials within the HFIR flux trap with both HEU and candidate LEU core configurations. Seven analyses were conducted to evaluate key performance metrics, including fast neutron flux distribution, cross section response functions, cross section data, and local dpa and heat deposition rates using mesh- and cell-based tallies. Simulations employed Shift, Monte Carlo N-Particle (MCNP), and the HIFR Controller (HFIRCON) tool suite for high-fidelity transport and depletion modeling. The LEU designs—using U 3 Si 2 -Al dispersion fuel and operating at 95 MW—were compared to the current 85 MW HEU configuration. Results show that while the candidate LEU cores exhibit higher dpa rates due to a harder spectrum and extended cycle lengths, they also demonstrate reduced heat deposition rates in irradiation experiments, primarily due to increased gamma self-shielding from higher 238 U content in the core. These findings confirm that LEU conversion can maintain HFIR’s materials irradiation capabilities but may require redesigning existing experimental hardware.

HEU↗

Porosity in nuclear graphite and its impact on nuclear reactor science and criticality safety applications

Porosity in nuclear-grade graphite significantly influences its low-energy neutron scattering, yet its effect on underlying phonon properties remains debated. Here, this work integrates inelastic and small-angle neutron scattering (INS/SANS) experiments, advanced atomistic simulations with a novel machine-learned potential (DeepMD), total cross-section measurements, and neutronics calculations (SCALE, MCNP, OpenMC) to investigate porosity’s impact on neutron thermalization. INS measurements on diverse graphite grades reveal no discernible porosity effect on phonon spectra, which align with crystalline graphite. Conversely, total cross-section data below ≈10 meV show increased scattering attributable to SANS. Our DeepMD simulations demonstrate that realistic micropores do not distort phonon spectra, challenging the assumptions in current ENDF/B-VIII.1 porosity thermal scattering laws (TSLs). These TSLs, based on random atom removal, produce unphysical phonon spectra and inflate inelastic cross-sections. Augmenting a crystalline TSL with an SANS component accurately captures experimental total cross-sections. Neutronics benchmarks (ICSBEP/IRPhE) show ENDF porosity TSLs unphysically increase neutron multiplication factor, keff. Crucially, incorporating SANS physics (NCrystal/OpenMC) indicates accurately modeled porosity negligibly affects keff, reactor physics, or criticality safety.

Critical benchmarks↗

Demonstration of RELAP5-3D for transient analysis of a dual coolant lead lithium fusion blanket concept

Our research aims to build a multiphysics framework for transient analysis of the Dual- Coolant Lead-Lithium (DCLL) blanket design of the proposed Fusion Nuclear Science Facility (FNSF) for all Onticipated Operational Occurrences (AOOs) using Reduced Order Models (ROMs). This paper presents an assessment of using RELAP5-3D for transient thermal-hydraulic blanket analysis of the DCLL design. Preliminary validation studies for the RELAP5-3D properties were conducted using a representative vertical flow loop model and a comparison study using experimental data from the thermal convection corrosion loops at ORNL. These studies inspected the thermal hydraulic response of Lead-Lithium Eutectic (PbLi) systems within the RELAP5-3D code; mainly pressure drop and heat transfer, as compared to analytical calculations and experimental data. The Magnetohydrodynamic (MHD) pressure drop effect for liquid metal under the effects of a uniform magnetic field was implemented into our models using an equivalent forms loss coefficient based on the phasic momentum equations. We verified our implementation of MHD pressure drop within RELAP5-3D using a systematic approach consisting of two studies. Models were developed to perform comparison studies based on experimental data from the Argonne Liquid metal EXperiment (ALEX) facility and the Magnetohydrodynamic PbLi Experiment (MaPLE) facility. Using the developed validation basis, we built a simplified DCLL channel model within RELAP5-3D, inspected MHD pressure drop and performed transient temperature analysis. The DCLL channel model used radial heating data from the International Thermonuclear Experimental Reactor (ITER) which was normalized to the FNSF design. Here, a representative startup transient was performed using a relative power curve from the initial ramp up of a long-term power pulse of the Tore Supra reactor. This study demonstrated that RELAP5-3D ROMs are capable of accurate transient analysis of the DCLL blanket. We have since utilized MCNP to develop heating profiles based on the FNSF design and are working toward implementing them into our DCLL model.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Virtual prototyping of liquid metal blanket performance in fusion pilot plant

Liquid metal blanket is a dominant design option for the next step fusion devices responsible for harvesting energy from fusion reaction, and simultaneously producing fuel for the same reaction through tritium breeding. Liquid metal blankets introduce additional complexity to the design due to fluid motion, fluid structure interaction, and magnetohydrodynamic (MHD) effects arising from the motion of the conducting fluid through the magnetic field. They are also directly affected by the plasma heat flux and neutronic fluence. PPPL is currently developing a virtual prototyping system for numerical analysis of the liquid metal blankets for future fusion devices. The system has a customized 3D computational fluid dynamics (CFD) code in its core, allowing MHD flow and conjugate heat transfer analysis in blankets fluids and solids. The code was successfully used before for dual coolant blanket analysis [A. Khodak et al., Fusion Eng. and Des. 137 (2018)]. Recently the same code was modified to allow verified simulation of MHD flows at high Hartmann numbers of several thousand typical for blanket applications. CFD code receives volumetric heat source distribution from the neutronic analysis based on MCNP code. In addition, direct tritium breeding simulation will be performed allowing optimization of the blanket performance. 2D axisymmetric version of neutronics code will be used for rapid optimization, with 3D version employed for detailed analysis. The surface heat distribution on the plasma facing wall will be defined by the software HEAT allowing 3D modeling of the heat flux based on the magnetic field distribution including gyro-orbit effects. Results of thermal analysis are imported into structural analysis code also included in the system. Finally, direct import of CAD geometry will be used for analyzing all components and as a result design option can be efficiently optimized.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

FAST irradiations and initial post irradiation examinations – Part I

The Advanced Fuels Campaign Fission Accelerated Steady-state Test (FAST) at Idaho National Laboratory (INL) completed its first irradiation cycle within the Advanced Test Reactor (ATR). The test focused on the irradiation of alloy fuel forms for use in sodium fast reactors. Tests included the use of fuel additives (Sn, Sb, and Pd) for reducing the gettering of lanthanide fission products, zirconium liners to mitigate fuel-cladding chemical interaction, and annular fuel geometry for sodium-free fuel designs. The first cycle of FAST testing was completed and a series of rodlets were removed for the initial post irradiation examination (PIE) of FAST rodlets. The rodlet irradiation conditions was evaluated using Monte Carlo N-Particle (MCNP) for as-run power history and COMSOL for temperature analysis. These rodlets include a set of low burnups (~2.5 % heavy metal [%FIMA]) control rodlets and a helium bonded annular rodlet (4.7 %FIMA). Initial non-destructive PIE has been completed and includes, visual inspection, neutron radiography and gamma scanning of the FAST rodlets. Radiography confirmed the integrity of the experiments as well as showing the complete filling of the annulus in the annular fuel at a modest burnup (4.7 %FIMA) and potential slumping of the cooler rodlets at lower burnup. Precision gamma scanning indicated mostly usual fission product behavior with the exception of the ceasium in the He-bonded annular fuel. Future destructive PIE will be necessary in order to fully understand the effects of accelerated irradiation on U-Zr metallic fuel behavior.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

New evaluation on angular distributions and energy spectra for neutron-induced charged-particle measurements

Recent development on neutron-induced charged particle reaction measurements can provide high-quality double differential cross sections at Los Alamos Neutron Science Center. With newly measured cross section data, we found lacking nuclear outputs for neutron-induced charged particle reactions in the latest ENDF/B-VIII.0 library, when compared with Monte Carlo simulations. Here, this paper presents how we improved the missing angular distributions and energy spectra for these reactions and the validation of our new evaluations with the measurements performed with the LENZ instrument at LANSCE. The impact of this additional evaluation is demonstrated by comparing the MCNP simulations with the LENZ experimental data.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Confirming the absence of nuclear warheads via passive gamma-ray measurements

Arms-control agreements between the United States and Russia negotiated after the end of the Cold War have imposed limits on the number of deployed strategic nuclear weapons. Verification of these agreements has relied on onsite inspections, sometimes supported by radiation detection techniques to confirm the absence of a nuclear warhead when ambiguities arise. So far, these measurements have sought to detect neutron emissions associated with the presence of plutonium, but they would be inadequate for uranium devices. In an effort to offer instruments that could be used to confirm the absence of both plutonium and uranium weapons, here we propose an inspection system that uses only passive gamma radiation detection techniques. Such a system would be particularly valuable for next-generation arms-control agreements that limit total numbers of weapons and would involve containerized items in storage. We conducted extensive Monte Carlo simulations to support the development of a verification protocol and detection algorithm. Here, we demonstrate the viability of the technique using standard laboratory check sources and MCNP simulations for simplified configurations of special nuclear material.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

A holistic approach to the optimization of neutron beam transport at the LANSCE facility

An advanced neutron beam collimation system has been developed as part of an effort to optimize neutron beam transport at the Weapons Neutron Research (WNR) Facility, located within the Los Alamos Neutron Science Center (LANSCE). The goal of this work was to develop and demonstrate techniques to tailor neutron delivery in order to provide maximum available flux on sample with a specific beam profile, while simultaneously reducing unwanted background. A holistic approach was taken, upgrading the facility spallation target, facility metrology infrastructure, and flight path shutter insert to support the implementation of this advanced collimation system. Here, modern instruments and software were employed to conduct facility surveys, characterization of as-built geometry of critical components, 3D ray tracing and neutron transport calculations. Beam images and flux measurements were acquired after installation to evaluate the performance of the collimation system and to demonstrate consistent agreement with MCNP simulations which showed an 87% increase of flux on sample while providing suppression of background neutron impingement on the sample frame by 10 3 .

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Lithium indium diselenide — An advanced material for neutron detection

Here, this paper describes the synthesis, crystal growth, detector fabrication, radiation hardening studies, MCNP modeling, and characterization of lithium indium diselenide or LiInSe 2 . This newly-developed room temperature thermal neutron detector has semiconducting and scintillating properties and it is suitable for neutron detection application. LiInSe 2 was synthesized starting from elemental Li, In, Se in two steps due to high reactivity of Li. A single crystal of LiInSe 2 was grown using the Vertical Bridgman method. The room temperature band gap was found to be 2.8 eV using optical absorption measurements. Bulk resistivity was measured at ∼5 × 10 11 Ωcm. Photoconductivity measurements of LiInSe 2 wafers identified a peak in the photocurrent around 445 nm. Nuclear radiation detectors were fabricated from single crystal wafer and the responses to alpha particles at various biases were measured. The mobility-lifetime product was estimated. Gamma irradiation studies were performed with calculated absorbed doses ranging from 0.2126 to 21,262 Gy. The characterization of the two wafers for their scintillator performance was conducted after each irradiation. The gamma irradiation produced a reduction of the light yield that translated to a lower channel number for the centroid of alpha detection spectra. It also showed a considerable reduction of the decay time after the first irradiation. These are the first studies on gamma radiation hardening with this material.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Experimental demonstration and modeling of a robotic neutron detector with spectral and directional sensitivity for treaty verification

International safeguards and arms control agreements often require labor-intensive, intrusive onsite inspections to perform verification tasks. The ability to localize a neutron source and/or characterize a neutron field may be imperative for identifying anomalies. We are interested in the role of autonomous mobile robots, which, if designed properly, may be more effective and efficient and less intrusive than their human counterparts. Toward developing such a capability, we, for this study, have constructed the N-SpecDir Bot, comprised of three boron-coated straw detectors azimuthally-distributed within a cylinder of high-density polyethylene, which is mounted on an omni-directional robotic platform. Our N-SpecDir Bot is specifically designed to provide spectral and directional sensitivity, in addition to gross counts, by utilizing the signals from the three detectors. The detection system has been extensively characterized by MCNP modeling, which has been benchmarked to experiments conducted at the Princeton Plasma Physics Laboratory. We demonstrate the spectral and directional sensitivity experimentally and in simulation, and provide a simple yet robust method for directional measurements.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

Spherical time-encoded radiation imaging simulations

Radiation source localization is important for nuclear nonproliferation and can be obtained using time-encoded imaging systems with unsegmented detectors. A scintillation crystal can be used with a moving coded-aperture mask to vary the detected count rate produced from radiation sources in the far field. The modulation of observed counts over time can be used to reconstruct an image with the known coded-aperture mask pattern. Current time-encoded imaging systems incorporate cylindrical coded-aperture masks and have limits to their fully coded imaging field-of-view. This work focuses on expanding the field-of-view to 4π by using a novel spherical coded-aperture mask. A regular icosahedron is used to approximate a spherical mask. This icosahedron consists of 20 equilateral triangles; the faces of which are each subdivided into four equilateral triangle-shaped voxels which are then projected onto a spherical surface, creating an 80-voxel coded-aperture mask. Furthermore, these polygonal voxels can be made from high-Z materials for gamma-ray modulation and/or low-Z materials for neutron modulation. In this work, we present Monte Carlo N-Particle (MCNP) simulations and simple models programmed in Mathematica to explore image reconstruction capabilities of this 80-voxel coded-aperture mask.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Applied radiation tags for warhead dismantlement transparency

Here, this study explores the feasibility of applying radiation tags as a Chain of Custody tool to confirm the dismantlement of individual warheads by investigating time dependent emissions from special nuclear material. Precedents of photon and neutron sources that have been used to irradiate warhead components are presented. Measurements of irradiated Highly Enriched Uranium (HEU) are used to benchmark simulations with the CINDER90 library for MCNP. Delayed gamma emissions are simulated as a function of time since irradiation of moderated HEU and Weapons Grade Plutonium with photon and neutron sources. Models of commercial 9 MV and 15 MV linacs were found to induce ~10 9 fissions/g and ~10 11 fissions/g, respectively, creating 3 to 4 orders of magnitude more delayed gammas than neutron sources, modeled at 10 10 n/s. However, a DT neutron source was found to induce sufficient fissions, ~10 7 fissions/g, to create quantifiable signatures relative to intrinsic emissions. Applied radiation tagging is determined to be technically feasible, but its operational context and practical implementation requires significant development.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

Characterizing the Experiment for Calibration with Uranium (Excalibur) neutron source for use in warhead verification

Neutron sources can play a variety of roles in warhead verification. For transmission radiography, a source of directed high energy neutrons is required, while for applications to detect fissile isotopes, sub-MeV neutrons are preferred. The Excalibur (Experiment for Calibration with Uranium) neutron source has been built and used in a variety of verification-related experiments. Excalibur is based on a commercial deuterium-tritium neutron generator specified and measured to be capable of producing 14 MeV neutrons at rates of up to 8.2 × 10 8 neutrons/s. Here, the generator is enclosed in a carbon-steel 32" diameter, 23.62" high carbon-steel cylinder that moderates the mean neutron energy to under 500 keV. This, in turn, is encased in 5%-borated polyethylene such that the entire assembly is a 48" x 48" box that is 30" tall. For radiographic applications, a narrow, tapered channel in the steel and polyethylene allows 14 MeV neutrons to stream directly from the generator to a test object. Its collimating capability is demonstrated by measuring the neutron flux profile. In the moderated mode of operation, the generator is fully enclosed in the steel, but a large section of the polyethylene is removed, providing a flux of sub-MeV neutrons from a wide range of angles. Neutron angular and spectral measurements using both a nested neutron spectrometer and a commercial liquid scintillator coupled with a 3 He detector show the expected softer neutron spectrum in moderated mode in good agreement with MCNP6 calculations. The gamma-ray spectrum from Excalibur is also in good agreement with MCNP modeling. Based on these findings, the future application of Excalibur in its two configurations is discussed.

Fissile material detection↗

Optimizing a detection system for fissile material in nuclear disarmament verification

In arms control treaties, verification plays a crucial role in detecting non-compliance, deterring future violations, and building trust between state parties. Neutron interrogation that induces fission reactions in fissile isotopes and measures the resulting fission neutrons, could be employed for this purpose. This study aims to develop a system which can determine the presence of fissile material while intrinsically protecting information. In this paper, we focus on optimizing the system for discriminating between an enriched uranium block from a depleted uranium (DU) block. The system was built and we report on benchmark measurements with DU and 16% enriched uranium blocks. Furthermore, the Excalibur (Experiment for Calibration with Uranium) neutron source, a neutron spectrometer (redBubble Technology Industries (BTI) N-Probe), and superheated droplet detectors were used for these measurements. MCNP simulations provided insights into detector responses to fissile materials with varying isotopic compositions, confirming that the system functioned as designed.

Active neutron interrogation↗

Methodology for inferring reactor core power distribution from an optical fiber based gamma thermometer array

Here, a 4-step data analytic methodology has been devised for the purpose of inferring the distribution of power in a reactor core, based on the response of an array of optical fiber based gamma thermometers (OFBGTs). This data analytic methodology is crucial for the development of a system of OFBGTs for the purpose of calibrating local power range monitors in boiling water reactors. Such a system would be an improvement to the present calibration system in boiling water reactors, in terms of safety, efficiency, and permanence. The first step of this methodology is to establish an energy balance method. In this method, one uses MCNP to determine response functions, which allow one to convert from gamma thermometer response to power. The gamma thermometers and the reactor core are segmented, such that each gamma thermometer segment provides an estimate of power for each reactor core segment. The estimates of power for the reactor core segments (hereafter referred to as fuel assembly segments) are calculated based on the measured response of the gamma thermometers. The second step of the methodology is to employ a weighting scheme to combine the estimates of the power of the various fuel assembly segments, based on the response of the various OFGBT segments and the incremental dose rates of the various fuel assembly segments to the OFBGT segments. The third step is to iteratively calculate the estimates of the power of the fuel assembly segments, until a convergence criterion is met, which indicates that the calculation has converged. The fourth and final step of the methodology is to estimate the uncertainty of the power for each fuel assembly segment. The mathematical basis for this step is not the focus of this paper. We have used a 3-D homogeneous reactor model to demonstrate the data analytic methodology; and have found that the data analytic methodology operates as intended.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

An approach towards plastic scintillators from thermally activated delayed fluorescent dyes and cross-linkable bismuth compounds

Cross-linked polyvinyltoluene (PVT) based plastic scintillators loaded with dyes exhibiting thermally activated delayed fluorescence (TADF) and various amounts of (caproate)di-(methacrylate) bismuth (CMB) were developed. The co-polymerization of vinyl toluene monomer loaded with different weight ratios of CMB and 1 wt% of TADF dye yielded highly transparent cross-linked plastic scintillators with intense green photoluminescence. The thermal and mechanical properties of the scintillators were investigated using thermogravimetric analysis (TGA), differential scanning calorimetry (DSC) and micro indentation techniques. Here, their performance for radiation detection was assessed by performing gamma spectroscopy using a Cs-137 source and benchmarked to commercial plastic scintillators. The loading of CMB was guided by Monte Carlo N-Particle (MCNP) simulations. The results of this study show that the cross-linking approach and the use of TADF dyes yield scintillators with good mechanical properties with a uniaxial yield strength up to ~66 MPa for compositions loaded with 40 wt% of CMB and high optical quality. However, the light yield in gamma spectroscopy experiments is found to be reduced with the addition of CMB, suggesting a trade-off between mechanical yield strength and light yield in gamma spectroscopy.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

An Improved Energy Deposition Model in MPACT and Explicit Heat Generation Coupling with CTF

The default energy deposition model in the CASL neutronics code MPACT assumes all fission energy is deposited locally in fuel rods. Furthermore, equilibrium delayed energy release is assumed for both steady-state and transient calculations. These approximations limit the accurate representation of the heat generation distribution in space and its variations over time, which are essential for power distribution and thermal-hydraulic coupling calculations. In this paper, an improved energy deposition model is presented in both the spatial and time domains. Spatially, the energy deposition through fission, neutron capture, and slowing-down reactions are explicitly modeled to account for the heat generation from all regions of a reactor core, and a gamma smearing scheme is developed that utilizes the gamma sources from neutron fission and capture. In the time domain, the delayed energy release is modeled by solving an additional equation of delayed heat emitters, similar to the equation of delayed neutron precursors.To allow the explicit heat generation coupling, the interfaces between MPACT and CTF were updated to transfer separate heat sources for different material regions (fuel, clad, moderator and guide tube). The results show that the distributions of the energy deposition between MPACT and MCNP agree very well for various 2-D assembly and quarter-core problems without TH feedback. The MPACT/CTF coupled calculation for the hot full power quarter-core case exhibited a reduced peak pin power by 2.3% and a reduced peak fuel centerline temperature by 17 K when using the explicit energy deposition and heat transfer. The new model also shows a maximum 100 pcm keff effect on assembly depletion problems and an increased overall energy release by 7% in a PWR reactivity-initiated accident (RIA) problem.

Liu, Yuxuan↗