SEARCH · Engineering Papers
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.
This is a generic presentation on the basics of radiation shielding as well as potential issues and possible solutions facing shielding fission batteries. It is to be presented at the Lightweight Materials Workshop on Nov 8 & 9 in EROB-159.
The US Department of Energy (DOE) Office of Nuclear Energy (NE) is currently investigating the feasibility of directly disposing dual-purpose (storage and transportation) canisters (DPCs) in a spent nuclear fuel (SNF) repository. Criticality during the repository performance period (10,000 years or more) is one of the major concerns related to direct disposal of DPCs, specifically as the system undergoes degradation in the repository environment and timeframe. Oak Ridge National Laboratory (ORNL) is developing an as-loaded criticality analysis methodology using full (actinides + fission products) burnup credit that exploits the inherent criticality margin associated with actual canister-specific loading configuration. Burnup credit criticality analysis requires validation of the depletion/decay codes used to generate the burned isotopic inventory of an assembly by comparing the code-predicted inventory with experimentally determined isotopic data. Currently, isotopic measurement data for boiling water reactors (BWR) SNF are limited, and additional measurements will be highly beneficial for BWR burnup credit analysis, which is essential to demonstrate disposability of BWR DPCs. Moreover, new pressurized water reactor (PWR) samples of isotopic measurements will expand the PWR sample population and consequently will greatly improve sample statistics. In turn, this will reduce uncertainty in the computational determination of the isotopic composition of commercial SNF and eliminate additional penalties currently used for lack of data. As part of the isotopic analysis task, eight diverse, high burnup (HBU) samples from several PWR rods are being dissolved and isotopically analyzed to provide high-quality measurement data for PWR SNF to reducing the uncertainties associated with PWR isotopic depletion validation. The isotopics of interest focus on nuclides important to burnup credit, shielding, and decay heat in PWR SNF. Additionally, eight high burnup BWR fuel samples from the Limerick nuclear power plant have been identified for measurement to support reducing uncertainties associated with BWR isotopic depletion validation.
Tritium confinement is performed using different barriers to minimize releases to the environment. Primary confinement is usually performed by process piping and components, secondary confinement typically by inerted gloveboxes connected to a continuously operating tritium stripper system. Some facility designs employ tertiary tritium confinement by using process room confinement with a tritium stripper system activated after an accident to further mitigate tritium releases to the environment. Secondary confinement of some tritium systems such as radiation-hardened or shielded-cell enclosures are not easily achieved using designs for inerted glovebox. These enclosures can be challenging to seal for secondary tritium confinement due to heating, ventilation, and air conditioning (HVAC) needs for temperature control, as well as feedthroughs for equipment and instrumentation, are potential pathways for tritium leaks or loses. Additional tritium leak paths can be created in secondary confinement enclosures due to a Design Basis Accident (DBA) seismic event. These new, larger tritium leak rates from the DBA event are usually not part of the stripper system design basis so all tritium released to the secondary confinement enclosure is assumed released to the environment. Detritiation strategies for shielded enclosures and/or gloveboxes with significant leak rates are necessary to minimize off-site radiological dose consequences. To address this, the relationship between the confinement system leak rate and stripper system (recirculating) flow rate was evaluated to meet prescribed maximum allowable environmental tritium emissions from the event. The simple analysis described assumes a constant leak rate from the confinement system (i.e. shielded enclosure or glovebox) while a recirculating stripper system strips tritium from the system– a competition between tritium release and tritium recovery. At a high level, the resulting expression summarizes the relationship between system leak rate (F L ) and stripper flow rate (F S ) relative to allowable tritium release (Q A ) and initial tritium release (Q 0 ): $\frac{F_L}{F_s}$ ≈ $\frac{Q_A}{Q_0}$. This report derives the relationship between these parameters, examines the impact of finite stripping times followed by purging of the stripped volume. The analysis provides example results for up to 30 gram tritium releases: the maximum tritium inventory of a Hazard Category III (Department of Energy) Nuclear Facility. The detritiation model presented represents a high-level analysis of tritium recovery versus losses for “leaky” confinement enclosures such as shielded cells after large tritium releases. Key parameters for doing the analyses are the system leak rate, stripper flow rate, initial tritium release, and allowable tritium release values. The analyses show a proportional decrease in releases with reduction in initial tritium release but a non-linear increase in tritium recovery with increased stripper flow rate or reduced leak-rate. The analyses also apply to glovebox confinement systems which could have significant increases in leak rates after a DBA. The analysis recognizes but does not include tritium absorption followed by re-emission from the walls of the confinement volume – an analysis which could be pursued in future studies.
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.
Computing in the energy dimension is one of the greatest challenges confronting present-day deterministic neutron transport solvers. Accurately resolving the neutron flux as neutrons downscatter across resonances in the nuclear cross sections currently requires considerable computing power and suffers from approximation errors. Flux uncertainty resulting from the uncertainty of the resonance structure is the single-largest cause of reactivity uncertainty. Any additional reference solution for the critical neutron downscattering problem with resonance phenomena would be a boon to verification and validation of neutronics codes. This paper establishes a benchmark to verify the accuracy of neutron transport criticality solvers along the energy dimension. For the first time, the analytic solution of the flux amplitude is derived in the particular case of an infinite homogeneous medium with isotropic scattering in the center of mass and an arbitrary number of no-threshold, neutral particle reaction resonances (e.g., radiative capture, fission, and resonance scattering). Furthermore, original analytic expressions are established to quantify the discrepancy between the ψ k (E) and ψ α (E) flux amplitudes, respective solutions of the multiplication factor k, or the exponential time-evolution frequency α eigenproblems. The physical study of these relations led to analysis of their first-order relative difference near the criticality condition α=0. Finally, numerical solutions are provided to a benchmark problem constituted of the first resonance of 239 Pu, the 6.67-eV resonance of 238 U, and a scattering isotope with a flat cross section, allowing for the computational verification of the energy resolution of current neutron transport criticality codes. Through these novel results, this analytic benchmark can serve as a reference to verify the energy resolution and sensitivity analysis of neutron transport criticality calculations.
Nuclear microreactors (MRs) represent a new class of reactors characterized by their compactness, portability, and low power output. These features enable MRs to supply electricity and process heat to remote areas like military bases; inaccessible locations; small grids, such as on islands; or disaster impacted areas. Compared to traditional light water reactors, MRs have a unique set of attributes that need to be considered for the implementation of safeguard strategies. Current safeguard methodologies are reactor technology specific and are employed on large, stationary reactors where there is easy access by safeguards inspectors and where safeguard equipment can be easily installed and retrofitted. While there are numerous benefits to MRs, their compact size, portability, scalability, and operational lifetime create challenges to the traditional safeguard approaches, thus needing novel safeguard strategies. Here, this paper addresses the unique challenges posed by MRs to the international nuclear safeguards regime, including limited human resources, and explores how technology advancements can help mitigate these challenges. Specifically, it examines novel technologies that could contribute to establishing a comprehensive safeguards framework for MRs. These safeguards-enabling technologies encompass safeguards by design, remote sensing and monitoring technologies, applications of artificial intelligence and machine learning algorithms, utilization of digital twins, and system of systems assessments. While each of these safeguards-enabling technologies offers partial solutions to the challenges posed by MRs for the international safeguards regime, none of them alone can entirely address these challenges. Consequently, a combination of the safeguards-enabling technologies outlined in this paper is recommended to establish a robust safeguards regime for MRs.
This report presents the final design (CED-2) of IER-501, a testbed for pulsed-neutron die-away (PNDA) experiments at Lawrence Livermore National Laboratory (LLNL). The PNDA experiments will be high-quality benchmarks that serve to optimize and validate thermal neutron scattering laws (TSLs), which are high-priority nuclear data for the Department of Energy’s Nuclear Criticality Safety Program (NCSP). This CED-2 was guided by experimental results gathered from a prototype experimental set-up rather than with detailed simulations as is typically done for criticality experiments. The report presents the PNDA design and its equipment, and it documents the methodology used to postprocess the measured data. It discusses how experimental configurations influenced the decay-constant eigenvalue, or the integral parameter, of the experiment. The studied variations included the distance from the target to the generator, the orientation of the generator relative to the box, the orientation of the target to the generator, the effect of HDPE lining in the shielding box, the settings of the generator, and finally the benefits of using a low-scatter facility. From these results, a final configuration for the target, detectors, shielding box, and neutron generator was chosen. Two factors particularly influenced the chosen configuration: the reproducibility of the experiment’s geometry for benchmark modeling, and the minimization of measured neutrons scattered off the shielding box.
Each year there are millions of dollars spent on the research and production of high-resolution detectors. This research indicates that the pursuit of higher resolution detectors is not always necessary. The terminal resolution of a NaI detector, or highest detector resolution, at which identification algorithms fail to identify highly enriched uranium (HEU) was evaluated using GADRAS, Genie, and GammaVision. GADRAS employs a template matching algorithm, while Genie and GammaVision utilize a mathematical approach for peak search and identification. The NaI spectra utilized for evaluation were generated using the GADRAS Inject tab and source modeling functions. Each spectrum included terrestrial and cosmic background from Dallas, TX. The resolutions for each spectrum were increased from a default 8.92% to a point where each algorithm would fail to identify 235 U from a HEU source. Six different source configurations were used in this research: bare HEU, 50% shielded HEU, 90% shielded HEU, bare HEU with an interference source of 99 mTc, bare HEU with 99 mTc both shielded 50%, and bare HEU with 99 mTc both shielded 90%. T
Current work includes reconstructing 187-group ENDF B/VII.1 covariance matrix for comparison study to 44-group ENDF B/VII.1 matrix to determine how covariance matrix structures affect the USL calculations and investigating the exact domination mechanism of region-wise sensitivities in a loosely-coupled system. Future work will involve seeking to understand how bias distributions change with reactivity and how to calculate propagate distribution error into USL calculations, as well as how cross section perturbation studies can be extended to other types of calculations, such as shielding calculations.
Here, the target segments of the Oak Ridge National Laboratory Second Target Station (STS) neutron production facility become highly activated due to spallation reactions or nuclei transmutation by primary protons and emitted neutrons. Once the target segments are removed from their location within the core vessel, decay dose rates must be accurately quantified to determine the shielding configurations of remote-handling tools and transport casks and to aid in planning maintenance activities. For this analysis, we utilized a hybrid unstructured mesh (UM)/constructive solid geometry approach for calculating spallation products and neutron fluxes, activation calculations using the AARE package that includes the CINDER2008 activation code to calculate the decay photon source at different cooling times, and the ADVANTG code to accelerate the final decay photon transport calculation. Both Type 316 stainless steel (SS-316) and lead were investigated as candidates for shielding materials. The decay photon transport calculation through the thick SS-316 or lead shields exhibited between 25 and 30 orders-of-magnitude attenuations in the radial direction, depending on the shield. Such a difficult shielding calculation required advanced variance reduction. ADVANTG has some missing features, which limits its usability in spallation neutron source applications. It does not support volumetric sources created for MCNP6.2 UM capability. An approximate source was created for this problem. Not only was this approximate source needed for running the ADVANTG calculation to generate the weight windows, but also it was essential to develop source biasing (SB) parameters that were crucial for dramatically accelerating the decay photon transport in this problem. With this approximate source, the analysis was completed in a very reasonable computational time, and the design of the STS remote-handling equipment was finalized. This paper compares the efficiency of Monte Carlo simulations with different weight window and SB parameters calculated using different approximate ADVANTG calculations.
In the inverse radiation transport problem, a radiation signal emanating from an unknown object is used to reconstruct geometric and material parameters of the object , e.g., the mass of special nuclear material (SNM), the isotopic composition of the SNM (e.g., the enrichment of 235U), the thickness of shielding, etc. The inverse radiation transport problem is one of the most fundamental and important in the emergency response mission, and has been addressed in previous NA-22 and NA-42/84 investments by using derivative-free optimization methods. Our project was aimed at (i) consolidating previous work by applying polynomial chaos expansion (PCE) for reducing significantly the number of forward model simulations required by derivative-free optimization methods; (ii) developing a general-purpose framework and stand-alone software for solving the inverse radiation transport problem by applying a new “predictive modeling” methodology (predictive modeling of coupled multi-physics systems, PM_CMPS) which enables the fusing of experimental and computational information; and (iii) developing several benchmarks (for unscattered gamma rays) for demonstrating the quantification of non-Gaussian features (e.g., skewness and kurtosis, which characterize asymmetries and long tails) of response distributions by deterministically and efficiently propagating parameter uncertainties (covariances and higher-order moments, where available) using high-order sensitivities computed by the recently developed comprehensive adjoint sensitivity analysis methodology (C-ASAM). With regard to (i), application of PCE, the most significant accomplishment of this project was an implementation of PCE into INVERSE that significantly reduced computation time for the DiffeRential Evolution Adaptive Metropolis (DREAM) and Differential Evolution (DE) methods. For inverse transport problems with measurements of neutron leakage, neutron multiplication, and neutron-induced gamma-rays, the solution times for the DREAM and DE methods were reduced by factors of 30-300, reducing run times from hours to minutes while maintaining the same accuracy as the standard (non-PCE) DREAM and DE. Speedups were further improved by exploring sparse grid quadrature sets using the Toolkit for Adaptive Stochastic Modeling and Non-Intrusive ApproximatioN (Tasmanian) toolset. An initial exploration of sparse grids indicated that further speedups factors of 5-10 are obtainable. The implementation of PCE in INVERSE is not in a form that users can easily access. With regard to (ii), solving the inverse problem using a predictive modeling paradigm, the most significant accomplishment of this project were the derivation and implementation of adjoint-based second derivatives in neutron and photon transport theory. There are two codes that compute these derivatives. SENSPG computes them for uncollided gamma rays and SENSMG computes them for neutron leakage. SENSMG is open-source code available on github. Other significant accomplishments included a stand-alone code, MULTI-PRED, that implements the predictive modeling paradigm, along with a detailed manual. The code was implemented in the existing INVERSE code for inverse problems, but a head-to-head comparison of inverse methods was not achieved. We also derived and implemented the capability to compute sensitivities for multiplicity measurements into the SENSMG code, and this capability was used in a data assimilation application, but it was not used for solving inverse problems. With regard to (iii), developing benchmarks to demonstrate the importance of high-order sensitivity and uncertainty analysis, the most significant accomplishment of this project were two test problems that showed how important second-order sensitivities are.
The official ENDF/B nuclear data evaluation releases in the past two decades have incrementally incorporated more detailed information, new nuclide evaluations, and very often have improved the accuracy of radiation transport codes when modeling shielding, fission reactors, criticality benchmarks, and fusion systems. However, like any large collaborative data compilation, these releases have all included small errors. This report documents the small corrections to ENDF/B releases ENDF/B-VII.1, ENDF/B-VIII.0, and ENDF/B-VIII.1 that have been applied during nuclear data processing to produce data libraries for the SCALE code system.
Additive manufacturing (AM) represents a promising technique to fabricate metallic alloys with greater control of the resulting material features as compared to traditional manufacturing routes. Recently, there is greater interest in AM research on 9 wt% Cr ferritic/martensitic (F/M) steels, which are commonly studied for use in the nuclear energy industry. This work aims to prove that wire arc AM can manufacture F/M steels with adequate mechanical properties in multiple processing atmospheres and aims to study how shielding gas composition can be leveraged during fabrication to induce specific precipitation pathways. The effect of shielding gas composition on MX (M=Nb and/or V, Xdouble bondC and/or N) carbonitride precipitation in a 9 wt% Cr ferritic/martensitic (F/M) steel alloy known as Grade 91 was studied using N 2 and CO 2 gas additions to an inert Ar shielding gas atmosphere during wire arc AM. The N and C atoms present in the processing atmospheres were absorbed into the melt pools during fabrication. Due to their differing affinities for precipitate-forming reactions, the varying levels of C and N between the samples contributed to differences in final carbonitride composition and morphologies. Such precipitate behavior is of interest as carbonitrides have been shown to contribute to increased mechanical performance. This increased performance was studied via electron microscopy and tested for strength, ductility, and fracture properties.
The radiation shielding analysis of the external concrete enclosure of a feature long beamline of the Advanced Photon Source Upgrade project is described in this work. The beamline has 2 undulators that can operate in two different configurations and the source terms are estimated for the worst configuration using the SRW code in the OASYS package as well as the STAC8 code. The shielding analysis is carried out using the STAC8 code, and for the case of pink beam, with the FLUKA code also. Effective dose rates on contact with the outside surface of the concrete walls for pink, reflected full spectrum and monochromatic beams are calculated. For the monochromatic beams, 2 sets of discrete bandwidths (BW) from the Double Crystal Monochromator (DCM) and the full reflectance of the Double Multilayer Monochromator (DMM) are used. Further, due to the size of the direct beam and the energy domain of interest in the direct and scattered fields, the importance of the equivalent dose to the lens of the eye is discussed with respect to the shielding calculations of such facilities. The dose rates outside the concrete hutch in the lateral directions calculated using STAC8 is about 50% higher than the FLUKA results for the case of pink beam but may be considered as a reasonable agreement due to the differences in the calculational methodologies. Those dose rates in the lateral and forward directions are higher for the pink beam to allow full occupancy from the radiation protection point of view. For the full spectrum after 3 reflections, the effective dose rates are less than 0.5 Sv/h in the lateral direction but is higher in the forward region below about 22° with respect to the beam direction. For the DCM beams, except for the scenario with one mirror and 0.1%BW, all dose rates are less than 0.5 Sv/h outside the concrete hutch in the lateral directions. Here, the higher harmonics are what contributes to the dose outside the shield and the bandwidths assumed to derive the photon flux thus becomes important. The dose rates with the DMM are generally lower compared to the 0.1% BW DCM beam but are higher than the results obtained with XOP calculated bandwidths.
Uranium trioxide displays a complex chemical phase space, with at least six structurally distinct polymorphs accessible via different synthetic routes. Remarkably, despite its technological importance, full structural and electronic characterization of these polymorphs remains an open area of study. δ-UO 3 in particular has attracted significant theoretical attention due to its high point group and space group symmetries, having U (VI) in octahedral coordination with polyhedra interconnected through corner-sharing to build a 3-D cubic lattice with space group symmetry Pm-3m and Z = 1. Critical experimental information, such as its optical vibrational spectra, are not known. Here, we study the Raman and infrared (IR) spectra of δ-UO 3 together with the support of density functional theory (DFT) calculations for spectral interpretation. A symmetry analysis of the DFT-predicted phonon eigenmodes indicates that δ-UO 3 should have two IR active modes and no Raman active modes. Experimental results, however, indicate significant Raman scattering from δ-UO 3 . We therefore propose four potential explanations for this apparent contradiction: a possible tetragonal distortion to the cubic cell, the existence of a surface impurity layer, vacancy scattering, and structural activation of Raman signal. We use powder X-ray diffraction and confocal Raman spectroscopy with depth profiling to investigate these possibilities and suggest future experiments to explore this phenomenon in more detail. Understanding the lattice dynamics of δ-UO 3 is important for identification of technogenic U phases via Raman and infrared spectroscopy and our results indicate that the simple understanding of δ-UO 3 as a high-symmetry cubic structure should be reconsidered.
The Idaho Chemical Processing Plant (ICPP) has a number of criticality alarm systems (CAS) in the various facilities. The complexity of the plant creates some unique CAS problems. A detailed review of the CAS has been made to determine the level of compliance with ANSI 8.3-1986. A major CAS upgrade project is currently in the design stage which will expand coverage and bring all CAS monitoring unshielded areas into full compliance with ANSI 8.3-1986. Because DOE Orders do not contain an exclusion that exempts the need for CAS coverage in shielded areas as provided in ANSI 8.3-1986, a future CAS upgrade project is in the planning stage which will also upgrade the CAS coverage in shielded portions of the plant.
To validate and verify the python-based code (PySL), designed to replicate the programs used by STAYSL for Beam Correction Factor (BCF) and Self-Shielding Factor (SHIELD), a series of tests were performed. To test BCF a python script was written to generate a random flux history file and both versions of the code processed the data. The test verified matching values up to at least one decimal place, approximately 10,000 tests where run and each one passed. Isotopes began to fail the tests once neutron saturation was reached. To verify this the total time of exposure was varied the isotopes that failed were compared to a list of their half-lives. The test process for SHIELD was very similar but, in this case, the code began by producing an input file with varying thickness and device type/environment for the SHIELD input. The failure condition for this test was if any of the data points for an isotope had a difference above 3%. Approximately 40 of these tests were run and there were only 3 isotopes that had reoccurring failures but only 2% of their points were above the 3% difference. A visual comparison was conducted by plotting the results from both programs. Although the test failed, the differences between their values were minuscule, and the self-shielding factor’s shape was preserved when plotted. Next steps for this project will be validating and verifying the python-based SigPhi code and then reproducing and testing the least squares unfolding performed by STAYSL.
Screening and inspection of cargo containers are two essential methods to nondestructively examine the contents of shipment. These methods enable the detection of illicit transportation of unauthorized materials such as nuclear and radioactive materials, explosives, drugs, and so on, typically at borders or secure facilities. Although high-energy X-ray transmission is a standard system and is widely used for cargo inspection, the inherent challenges of high false-positive rates and high attenuation factors necessitate the development of complementary techniques that can increase the detection efficiency and accuracy in large and dense materials. Gamma-rays, which possess higher penetration characteristics because of their high energy, offer an alternative nonintrusive modality for cargo scanning. They represent a promising inspection method when compared to X-rays for three reasons: (1) improved ability to detect nuclear and radioactive materials, (2) higher inspection throughput rates, and (3) lower false-positive rates. Currently, there are two main gamma-ray inspection techniques, active and passive interrogation. Active interrogation can be further grouped into (1) gamma-ray transmission imaging and (2) neutron-induced gamma-ray emission detection. Gamma-ray transmission imaging utilizes differences in material densities for mapping the shipment contents and detecting anomalies. It is analogous to the X-ray transmission method; however, the high-energy photons make it more difficult to shield against, which enables more efficient performance in large and dense material inspection. Neutron-induced gamma-ray emission inspection is designed for the detection of nuclear and radioactive material because those materials emit characteristic gamma-rays when they are activated by neutron absorption. On the other hand, passive interrogation techniques rely on high-efficiency detectors to detect radiation emitted from hidden special nuclear or other radioactive materials. Similar to passive interrogation, cosmic ray muon monitoring and imaging are relatively new techniques that do not require external radioactive sources. These techniques have received attention as a potential next-generation radiographic probe to identify illicit transportation of nuclear and radioactive materials in cargo containers. Cosmic ray muons have unique features, (1) much higher energies than X-rays or gamma-rays (on the order of 10−1—104 GeV), (2) enhanced penetration capability, and (3) natural occurrence, thereby eliminating the need for induced radiation sources. These features enable cosmic ray muons to be utilized for detection of special nuclear materials in high-background-noise environments. By analyzing incoming and outgoing muon trajectories, scattering angles, and energies, it has been shown that it would be possible to locate hidden and well-shielded materials in cargo containers via three-dimensional muon tomography images or signal analysis. Gamma-rays, cosmic ray muons, and other nonintrusive cargo inspection modalities are complementary to each other, allowing them to address various cargo inspection conditions (i.e., scanning time, cost, radiation exposure level, and types of target materials). This chapter presents a detailed review of the theoretical fundamentals and technical principles behind the current gamma-ray and cosmic ray muon modalities for cargo inspection. Additionally, critical assessments and suggestions for the future directions to advance the use of gamma and muon modalities are discussed.