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.
The presentation discusses the scope of EMPIRE, EMPIRE's convenience, the EMPIRE-3.2 (Malta) Nuclear Reaction Model Code, reaction models, and what is needed to improve EMPIRE's predictive power.
Many evaluation, assessment, and modeling tasks at nuclear power plants require spatial in-formation; this often requires physical visits to locations within the facility because the 2D or 3D schematics and current models do not contain enough detail or do not capture as-built and real-world conditions. These visits require extensive manual labor for not only the requesting party, but also support groups, such as security. Light Detection and Ranging (LiDAR) mapping is trying to solve that problem by providing very detailed 3D models for low costs. However, the use of these models can be very limited because either component reference information is missing and too costly to add or there is no way to extract specific spatial data needed for other tools. This report outlines two main efforts. First, to reduce the effort of “Tagging” data in large 3D models, a general Application Programming Interface (API) was developed to import a variety of existing plant database information into a 3D-visualization engine. Filters allow the user to have only zone-specific items listed; then, they can simply click and assign the information to a specific spot or component in the 3D model. The second part of the work is the development of an interface for importing pieces from the3D-LiDAR model into other systems needed for modeling and simulation, outlined around fire modeling. This interface allows for the retrieval of item location and boundaries, enabling the auto generation of models for varying tools.
Improving and adapting industrial systems to timely meet changing programmatic and market demands is an important goal to achieve, including when operating and maintaining complex nuclear processes and facilities. However, changes to these complex systems are costly, particularly when they are already in place and bounded to stringent requirements and constraints such as when handling radioactive material and contaminated equipment. These conditions often exist when treating spent nuclear fuel remotely within shielded nuclear radiation chambers, commonly referred as hot cells, to condition nuclear material and/or fabricate products for utilization in other nuclear enterprises such as in the manufacture of advanced nuclear fuel. The illustrative case considered here is the production of high assay low enriched uranium (HALEU) products supporting the deployment of advanced nuclear reactors. For the HALEU program, resources invested were and are being systematically analyzed so that these investments are maximized in a facility that is nearly 60 years old. A methodology that has effectively enabled optimized and improvements in the Spent Fuel Treatment (SFT) program, and consequently the HALEU program, involves discrete event simulation as addressed in this article. Here, the quantification of multiple productivity metrics, including material processing rates, cycle times, bottlenecks, number of material transfers as well as equipment, workstation, and material handling utilization, has resulted in a myriad of diverse discoveries and data-informed decisions regarding process layout and constituent, labor levels and schedules, selection of new process units, storage needs, and other critical process configurations. This article describes such a computational capability being applied for decision-making, illustrates its application to an actual process and program, provides illustrative results, and argues how computational methods for the modeling, analysis, and optimization of complex processes and facilities does lead to informed decisions derived from data and not only from intuition.
Whether or not femto-scale droplets of quark-gluon plasma (QGP) are formed in so-called small systems at high-energy colliders is a pressing question in the phenomenology of the strong interaction. For proton-proton or proton-nucleus collisions the answer is inconclusive due to the large theoretical uncertainties plaguing the description of these processes. While upcoming data on collisions of 16 O nuclei may mitigate these uncertainties in the near future, here we demonstrate the unique possibilities offered by complementing 16 O + 16 O data with collisions of 20 Ne ions. We couple both nuclear lattice effective field theory (NLEFT) and projected generator coordinate method (PGCM) ab initio descriptions of the structure of 20 Ne and 16 O to hydrodynamic simulations of 16 O + 16 O and 20 Ne + 20 Ne collisions at high energy. We isolate the imprints of the bowling-pin shape of 20 Ne on the collective flow of hadrons, which can be used to perform quantitative tests of the hydrodynamic QGP paradigm. In particular, we predict that the elliptic flow of 20 Ne + 20 Ne collisions is enhanced by as much as 1.174(8) stat (31) syst for NLEFT and 1.139(6) stat (39) syst for PGCM relative to 16 O + 16 O collisions for the 1% most central events. At the same time, theoretical uncertainties largely cancel when studying relative variations of observables between two systems. This demonstrates a method based on experiments with two light-ion species for precision characterizations of the collective dynamics and its emergence in a small system.
Abstract not provided.
Abstract not provided.
The Full Law Analysis Scattering System Hub (FLASSH) is a modern, advanced code which evaluates the thermal scattering law (TSL) along with accompanying cross sections. FLASSH features generalized methods which accommodate any material structure. Historical approximations including the incoherent and cubic approximations have been removed. Instead, the latest release of FLASSH features advanced physics options including distinct corrections (1-phonon contributions) and non-cubic formulations. The non-cubic elastic and inelastic contributions are necessary to accurately evaluate 1-phonon contributions. Both non-cubic and 1-phonon calculations require high-density sampling of the various scattering directions. Optimization and parallelization of these routines were therefore necessary to produce results in a reasonable timeframe. With these notable improvements to the generalized TSL, FLASSH 1.0 meets benchmark requirements, demonstrating noticeable agreement with experiment for both TSLs and the resulting integrated cross sections. Additional features including a graphical user interface (GUI), plotting diagnostics, and formatted output options including ACE files allow users to complete a TSL evaluation with minimal input and maximum flexibility. The user GUI creates input files for FLASSH, reducing user error and also providing built-in error checks. Autofill options and suggested input values help make TSL evaluation accessible to novice users. The FLASSH code is compiled to run on both Windows and Linux platforms with automatic parallelization. (authors)
Near the critical point in the QCD phase diagram, hydrodynamics breaks down at a momentum where the frequency of the fastest hydrodynamic mode becomes comparable with the decay rate of the slowest nonhydrodynamic mode. Hydro+ was developed as a framework which extends the range of validity of hydrodynamics beyond that momentum value. This was achieved through coupling the hydrodynamic modes to the slowest nonhydrodynamic mode. In this work, analyzing the spectrum of linear perturbations in Hydro+, we find that a slow mode falls out of equilibrium if its momentum is greater than a characteristic momentum value. That characteristic momentum turns out to be set by the branch points of the dispersion relations. These branch points occur at the critical momenta of so-called spectral curves and are related to the radius of convergence of the derivative expansion. The existence of such a characteristic momentum scale suggests that a particular class of slow modes has no remarkable effect on the flow of the plasma. Based on these results, we show that there is an enhancement of the speed of sound due to the Hydro+ corrections compared to the known thermodynamic corrections near the critical point. This enhancement can not exceed a temperature-dependent upper bound which we derive near the critical point in the QCD phase diagram.
The speed of sound squared (SSS) s 2 in massive neutron stars (NSs) characterizes not only the stiffness of supradense neutron-rich matter within but also equivalently properties of the curved geometry due to the strong-field gravity and matter-geometry coupling. A peaked density or radius profile of s 2 has been predicted for massive NSs using various NS equation of state (EOS) models. However, the nature, cause, location and size of the peak in s 2 profiles are still very EOS model dependent. In this work, we investigate systematically s 2 profiles in massive NSs in a new approach that is independent of the nuclear EOS model and without any presumption about the NS structure and/or composition. In terms of the small quantities (reduced radius, the energy density and pressure scaled by their central values), we perform double-element perturbative expansions in solving perturbatively the scaled Tolman-Oppenheimer-Volkoff (TOV) equations and analyzing s 2 profiles from the Newtonian limit to the general relativistic (GR) case. The GR term in the TOV equations plays a twofold role: it compresses NS matter and modifies the pressure/energy density ratio from small values in Newtonian stars showing no s 2 peak to large ones for massive NSs possessing a peak in their s 2 profiles, and eventually takes away the peak in extremely compact/massive NSs approaching the causality limit. In particular, the peaked behavior in s 2 is expected to emerge near the center of massive NSs like PSR J0740+6620, while a sharp phase transition is unlikely to occur there. Furthermore, these features revealed from our analyses are universal as they are intrinsic properties of the GR stellar structure equations independent of the still very uncertain EOS of supradense neutron-rich matter in NSs.
This report from Idaho National Laboratory (INL) summarizes the key modeling and simulation activities conducted under the Department of Energy (DOE) Molten Salt Reactor (MSR) Campaign during the Fiscal Year 2023 (FY23). The focus of the work was to leverage state-of-the-art modeling capabilities from the DOE Nuclear Energy Advanced Modeling and Simulation (NEAMS) codes to enable novel multiphysics and multiscale modeling and simulation of MSRs. Through collaboration with NEAMS code developers, advanced multiphysics analysis capabilities for MSR systems were demonstrated by coupling depletion, thermal-hydraulics, and thermochemistry into an innovative framework for chemical species transport in MSRs. As a result, the framework can track nuclides throughout their lifetimes in the core, from production (depletion) to advection throughout the salt volume (thermal-hydraulics) and off-gassing or precipitation outside of the salt (thermochemistry). This work supports the near-term deployment of MSRs by integrating the synergistic efforts between the DOE’s MSR Campaign and NEAMS program. The resulting framework will help better connect system design modelers with experimentalists to better understand and predict complex physical behaviors in MSRs. Researchers and MSR developers alike can now leverage these new modeling and simulation capabilities to perform novel analyses with applications including: • MSR dynamics during normal operational transients and accident scenarios • Off-gas system design and performance for fuel cycle and depletion analysis • Corrosion and active chemistry control for reactor component health and lifetime determination • Source term, decay heat and activity determination in accident scenarios • Special nuclear material accountancy and chemical forensic analysis for safeguards • Digital twin development of experiments and experimental reactor demonstrations • Measurement requirements for instrumentation and control design • Uncertainty and sensitivity analysis of missing data to inform future experimental data collection.
A detailed study of high-energy processes in solar flares is reported, including the production of neutrons and pions, and incorporating isobaric and scaling models and a recent compilation of pion production data (Dermer, 1986). The broad-band gamma-ray spectrum resulting from the decay of neutral pions, the bremsstrahlung of positrons and electrons from the decay of charged pions, and the annihilation in flight of positrons is evaluated. Also evaluated is the 0.511 MeV gamma-ray line resulting from the annihilation of the positrons which survive annihilation in flight. Calculations were based on an isotropic, thick-target model using the best available nuclear data and models. Results are compared with extensive observation of the June 3, 1982 flare (10-120 MeV gamma rays), 0.511 and 2.2 MeV line emission, nuclear line emission, high-energy neutrons, and interplanetary charged particles.
Large-scale models of nuclear structure are currently the only way to provide consistent datasets for the many properties of thousands of exotic nuclei that are required by nucleosynthesis simulations. In [W. Ryssens et al., Eur. Phys. J. A 58, 246 (2022)], we recently presented the new BSkG2 model based on an energy density functional of the Skyrme type. Relying on a flexible three-dimensional coordinate representation of the nucleus, the model takes into account both triaxial deformation and time-reversal symmetry breaking. BSkG2 achieves a state-of-the-art global description of nuclear ground state (g.s.) properties and reproduces in particular the known masses with a root-mean-square (rms) deviation of 678 keV. Moving beyond g.s. properties, the model also reproduces all empirical values for the primary and secondary barriers as well as isomer excitation energies of actinide nuclei with rms deviations below 500 keV, i.e. with unprecedented accuracy. Here we discuss in detail the extension of our framework to the calculation of the fission barriers of 45 actinide nuclei, including odd-mass and odd-odd systems. We focus in particular on the impact of symmetry breaking which is key to the accuracy of the model: we allow systematically for axial, reflection and time-reversal symmetry breaking. The effect of the latter on the fission properties of odd-mass and odd-odd nuclei is small, but we find that allowing for shapes with triaxial or octupole deformation, as well as shapes with both, is crucial to achieving this accuracy. Here, the numerical accuracy of our coordinate space approach, the variety of nuclear configurations explored and the simultaneous successful description of fission properties and known masses makes BSkG2 the tool of choice for the large-scale study of nuclear structure.
Various aspects of nuclear electric propulsion (NEP) systems analysis and modeling are discussed. The following specific topics are covered: (1) systems analysis challenges; (2) goals for NEP systems analysis; (3) the Nuclear Propulsion Office approach; and (4) NEP subsystem model development. The discussion is presented in vugraph form.
This work presents the latest improvements to, and investigations performed with, the Fluoride-Cooled High-Temperature Pebble-Bed Reactor reference plant models for the United States Nuclear Regulatory Commission. These models, developed with the Comprehensive Reactor Analysis Bundle, or BlueCRAB, serve as the foundation for the future development of detailed design evaluation models based on license applications. BlueCRAB is the code suite proposed for non-light-water reactor systems safety analyses, and it incorporates various simulation tools developed by the Nuclear Energy Advance Modeling and Simulation program, including the Griffin code for reactor physics, the Pronghorn and SAM codes for core thermal fluids, the BISON code for solid conduction and fuel performance, and the SAM code for system analysis. The primary objective of this work is to assess the level of readiness of BlueCRAB for modeling fluoride-cooled high-temperature pebble-bed reactors. We first developed numerical models in BlueCRAB that include the key physics for this technology, ensuring an adequate level of fidelity for modeling the core performance during accident scenarios. This was followed by simulation of transient scenarios, two loss-of-forced-cooling events (one protected and one unprotected), and two control rod withdrawal events (one delayed and one prompt supercritical reactivity insertion). The analysis includes comparisons between the 2-D thermal fluid porous media models in Pronghorn and SAM, comparisons between coupled Pronghorn-Griffin and coupled SAM-Griffin models for two loss-of-forced cooling events and one control rod withdrawal event, and comparisons between SAM single-solve and domain-overlapping approaches for multi-scale thermal fluid coupling. In addition, we performed comparisons between 3-D, 2-D, and 0-D neutronic models for the two control rod withdrawal scenarios with Pronghorn-Griffin. The results show that the BlueCRAB models led to physically intuitive solutions for the scenarios examined. The changes in the various scalar and vector fields such as the neutron flux, power, temperatures, densities, pressures, and velocities are all within the expected ranges, and their distributions can be explained from the system response of the transients and the geometric and material variations. Several comparisons suggest that the porous media models in Pronghorn and SAM can lead to similar solutions, even though they are based on different methodologies. The simulations demonstrate that there are differences between the various levels of fidelity, and it is advisable to have flexible tools that can cover the breadth and depth of needs that may arise in future technical evaluations. We believe that BlueCRAB’s capabilities represent a significant asset for confirmatory analyses aimed at resolving important safety questions.
Background: The nucleus 32 Mg (N=20 and Z=12) plays a central role in the so-called “island of inversion,” where in the ground states sd-shell neutrons are promoted to the fp-shell orbitals across the shell gap, resulting in the disappearance of the canonical neutron magic number N=20. Purpose: The primary goals of this work are to extend the level scheme of 32Mg, provide spin-parity assignments to excited states, and discuss the microscopic structure of each state through comparisons with theoretical calculations. Method: In-beam γ-ray spectroscopy of 32 Mg was performed using two direct-reaction probes: one-neutron (two-proton) knockout reactions on 33 Mg ( 34 Si). Final-state exclusive cross sections and parallel momentum distributions were extracted from the experimental data and compared with eikonal-based reaction model calculations combined with shell-model overlap functions. Results: Owing to the remarkable selectivity of the one-neutron and two-proton knockout reactions, a significantly updated level scheme for 32 Mg, which exhibits negative-parity intruder and positive-parity normal states, was constructed. The experimental results were confronted with four different nuclear structure models. Conclusions: In some of these models, different aspects of 32 Mg and the transition into the island of inversion are well described. However, unexplained discrepancies remain, and, even with the help of these state-of-the-art theoretical approaches, the structure of this key nucleus is not yet fully captured.
In this work, we experimentally investigated the batch-to-glass conversion kinetics of two low-activity and two high-level nuclear waste melter feeds that represent widely different behavior during melting; tests comprised thermogravimetric analysis, the feed expansion test, evolved gas analysis, and x-ray diffraction. Our study indicates that the conversion processes, which include gas-evolving reactions, dissolution of crystalline phases, and primary foam expansion and collapse, are interconnected and influenced by the feed’s thermal history. The fraction of dissolved silica was adopted as a measure of the extent of conversion inside the cold cap—the floating layer of feed material on the top surface of the melt in an electric melter. The interface between the cold cap and the free-flowing melt below was associated with a constant fraction of dissolved silica. The estimated viscosity of glass-forming melt at this interface lies within a narrow range for the feeds studied.
This work presents the latest improvements to, and investigations performed with, the pebble-bed high-temperature gas-cooled reactor (PB-HTGR) reference plant models for the United States Nuclear Regulatory Commission. These models serve as the foundation for the future development of detailed design evaluation models based on license applications. The reference plant models have been developed with the Comprehensive Reactor Analysis Bundle, or BlueCRAB, which is the code suite proposed for non-light-water reactor systems safety analysis. It incorporates various simulation tools developed by the Nuclear Energy Advanced Modeling and Simulation program, including the Griffin code for reactor physics, the Pronghorn and SAM codes for core thermal fluids, the BISON code for solid conduction and fuel performance, and the SAM code for system analysis. The primary objective of the work that was performed was to assess BlueCRAB’s level of readiness for modeling a PB-HTGR. To do so, we first developed numerical models in BlueCRAB that include the key physics for this technology to ensure an adequate level of fidelity for modeling PB-HTGR core performance and for performing multiphysics simulations for equilibrium core conditions and different accident scenarios. Then we simulated transient scenarios, including depressurized and pressurized loss of forced cooling accidents, over-cooling, and control rod withdrawal events with delayed and prompt supercritical reactivity insertions. The analysis in this report includes comparisons of the 2D thermal fluid porous media models in Pronghorn and SAM, and comparisons of coupled SAM/Griffin/SAM and coupled Pronghorn/Griffin for depressurized and pressurized loss of forced cooling, over-cooling, and control rod withdrawal events. In addition, we compare 3D, 2D, and 0D/PKE neutronic models for the two control rod withdrawal scenarios with coupled Pronghorn/Griffin. The comparisons show that the BlueCRAB models lead to physically intuitive solutions for the scenarios examined. The changes in the various scalar and vector fields, such as neutron flux, power, temperature, density, pressure, and velocity, are within the expected ranges, and their distributions can be explained by the system response of the transients and the geometric and material variations. Several comparisons suggest that the porous media models in Pronghorn and SAM can lead to similar solutions, even though they are based on different methodologies. This work further highlights the need for flexible tools with various levels of fidelity to cover the breadth and depth of needs that may arise in future technical evaluations of the PB-HTGR. We believe that the BlueCRAB capabilities will be a significant asset for confirmatory analyses in order to resolve important safety questions.
The Galactic Cosmic Ray spectra in Mars orbit were generated with the recently expanded HZETRN (High Z and Energy Transport) and QMSFRG (Quantum Multiple-Scattering theory of nuclear Fragmentation) model calculations. These model calculations are compared with the first eighteen months of measured data from the MARIE (Martian Radiation Environment Experiment) instrument onboard the 2001 Mars Odyssey spacecraft that is currently in Martian orbit. The dose rates observed by the MARIE instrument are within 10% of the model calculated predictions. Model calculations are compared with the MARIE measurements of dose, dose-equivalent values, along with the available particle flux distribution. Model calculated particle flux includes GCR elemental composition of atomic number, Z = 1-28 and mass number, A = 1-58. Particle flux calculations specific for the current MARIE mapping period are reviewed and presented.