SEARCH · Search NASA
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 AMPX nuclear data processing code reads evaluated nuclear data files and prepares libraries in continuous energy (CE) and multigroup (MG) formats. These libraries are used by the transport and depletion codes in the SCALE code system.
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.
We explore a deep connection between fracton order and product codes. In particular, we propose and analyze conditions on classical seed codes which lead to fracton order in the resulting quantum product codes. Depending on the properties of the input codes, product codes can realize either Type-I or Type-II fracton models, in both nonlocal and local constructions. For the nonlocal case, we show that a recently proposed model of lineons on nonlocal graphs can be obtained as a hypergraph product code. Interestingly, constrained mobility in this model arises only from energy barriers associated with the graph. For the local case, we introduce a novel type of classical LDPC code defined on a planar aperiodic tiling. By considering the specific example of the pinwheel tiling, we demonstrate the systematic construction of local Type-I and Type-II fracton models as product codes. Our work establishes product codes as a natural setting for exploring fracton order.
A relativistic electron probe has been developed in the DIII-D tokamak, capable of simultaneously resolving pitch angles and energies of runaway electrons (REs) for the first time. Due to the relativistic speeds of REs, their gyro-orbit size becomes comparable to those of fast deuterium with energies in the tens of keV range. This allows for the measurement of RE strike images on a phosphor plane, with their orbits being deflected by the Lorentz force as they pass through a pinhole aperture. The strike positions correspond to the energies and pitch of the incident REs. Monte Carlo N-Particle Transport Code has shown that an ultra-thin phosphor coating significantly reduces the energy deposition from γ-rays, while allowing a much greater deposition from REs, minimizing the background noise. Finally, the novel system, developed for the DIII-D tokamak, is expected to provide unprecedented insights into the phase-space dynamics of REs.
PyARC was initially developed as an open source tool to support fast reactor analyses using the Argonne Reactor Computation (ARC) code suite as a part of the Nuclear Energy Advanced Modeling and Simulation (NEAMS) Workbench initiative in FY17. The goal of this initiative is to provide a common user interface for model generation, real-time validation, execution, output processing, and visualization for all integrated codes. This is accomplished through the reliance on tools available in the Workbench framework and runtime environment. While initially developed to support the ARC codes, PyARC was extended in FY22 to wrap other NEAMS and non-ARC codes, including Griffin and OpenMC, in the supported other neutronics workflows, and support users in the adoption of NEAMS-supported high fidelity analysis codes. Most recently, NUBOW-3D, a recently adopted ARC code, was integrated to support reactor bowing calculations as well. Integration of these codes into the NEAMS Workbench directly benefits the advanced reactor modeling community by: • Providing a set of controlled, maintained, documented and validated scripts to generate inputs, which promotes best practices, reduces the learning curve, and facilitates project collaboration. • Improving the user experience: the Workbench interface provides assistance for building an input through auto-completion, real-time validation, document navigation, and geometry and results visualization. • Automating complex calculations and workflows for reactor analysis. • Helping users transition to using high-fidelity NEAMS codes along-side the ARC codes. In FY22, a progress report was published that described the state of each of the tools integrated into PyARC. Since then, there have been many enhancements and upgrades to the existing integrations as well as entirely new code integrations as well. This report details all new integrations and major developments in PyARC since the version 2.0.0 release highlighted in the FY22 report.
The ability to accurately model burn propagation in inertial confinement fusion plasmas is crucial for advancing fusion energy research. This work presents enhancements to the triforce hybrid fluid-kinetic multiphysics code, focusing on its kinetic half, which employs the particle-in-cell (PIC) method with Monte Carlo collisions (MCC). We use a moment-preserving collision model that mitigates numerical noise, particularly in spherical geometries where particle weights vary significantly. Additionally, we refine the treatment of inverse bremsstrahlung to account for electron–ion collision frequency reductions in degenerate plasmas and incorporate a blackbody radiation source to enable realistic photon injection. These improvements enable the simulation of 1-dimensional (1D) spherical fusion burn propagation in deuterium–tritium plasmas. Benchmark comparisons with the hydra radiation-hydrodynamics code confirm that triforce accurately captures the dynamics of hot-spot expansion and burn propagation, demonstrating sensitivity to ignition thresholds consistent with theoretical models. Findings show the ignition cliff to be less steep in our work compared to radiation-hydrodynamic modeling. These results highlight the role of kinetic effects in fusion ignition physics and underscore the necessity of hybrid fluid-kinetic models for advancing predictive capabilities in high-energy-density plasma systems.
This paper presents a proposed revision of the International Atomic Energy Agency transport regulations, related to the A 1 and A 2 limit values used to determine the radioactive transport classification. Based on the 'Q system', a novel methodology was introduced to derive Q A and Q B values related to scenarios involving external exposure from a distant source. These values are key parameters that respectively represent the total effective dose and total equivalent dose to the skin, from all primary and secondary particles contributing to radiation exposure. The International Working Group (WG A 1 /A 2 ) is established and associated with the TRANSSC Technical Expert Group on Radiation Protection. A review of the A 1 and A 2 values is performed in response to identified limitations within the existing Q system. The followed approach is based on Monte Carlo simulations that enabled the development of transfer functions aimed at reducing computational time and increasing the flexibility of dose evaluations for any radionuclide with known particle emission spectra. This method allows updating the Q A and Q B values to account for future data evolutions (decay data, fluence-to-dose conversion coefficients) and standardizing the calculation of regulation limits across all referenced radionuclides and scenarios related to external exposure. The transfer functions are established using three Monte Carlo simulation codes—FLUKA, Geant4, and MCNP—and address the previous limitations of the 'Q system', reflecting the latest International Commission for Radiation Protection recommendations and improvements in calculation techniques. The results of the WG show consistent agreement across the codes, with minor discrepancies observed at low primary energies due to statistical uncertainties and different handling of stopping power for electrons/positrons in the codes. This revised approach aligns with current standards and recommendations, ensuring that the radiological consequences of transport accidents are acceptable for the new A 1 and A 2 limits from a radiological protection perspective.
We present the first nonplanar hot cathode discharge simulations that capture the role of the trapped-ions plasma, elucidating new phenomena unobservable in planar geometric discharges. A discharge struck between a single emitting wire filament cathode and a bounding anode is simulated in cylindrical geometry using an axisymmetric (radial) particle-in-cell Monte-Carlo collisions code. Operating the discharge near its ionization energy threshold can lead to the formation of a two plasma mode (TPM). One plasma forms in the conventional upstream region through electron impact ionization of background neutrals. A second plasma, whose global effect on the discharge was not previously well understood, forms downstream through the trapping of cold ions in the potential well of the filament’s virtual cathode, a process enabled by ion-neutral charge exchange collisions. Three space charge regions intersperse the electrode gap—an emissive sheath between the cathode filament and trapped-ions plasma, a double layer between the two plasmas, and a classical sheath between the upstream plasma and the outer anode. Simulations exhibit mode transitions and quenching instabilities that transform the discharge between the TPM and other single-plasma sheath modes that include classical (temperature-limited), space charge limited, and inverse (anode glow) modes. The transitions are explained via “aid-and-compete” dynamics wherein the growth of one plasma enhances growth in the other while concurrently exhibiting expansion dynamics antagonistic to each other. The system exhibits strong hysteresis memory during the mode transitions. Improved understanding and control of these sheath mode transitions are expected to benefit plasma applications with hot cathodes.
The Thermo Scientific P385 Neutron Generator is a compact neutron source, producing 14 MeV neutrons through the deuterium-tritium (DT) fusion reaction. It is important to measure and understand the dependence of the neutron production rate on the accelerator current and voltage. In this study we evaluated neutron production with an absolutely calibrated liquid scintillator neutron spectrometer (BTI N-Probe), an absolutely calibrated He-3 detector surrounded by HDPE shells (Detec Nested Neutron Spectrometer, NNS), and two uncalibrated ZnS fast neutron scintillators (EJ-410), for both A3082 and A3083 sealed tubes. Here we also modeled the neutron yield using the TRIM code, which calculates the trajectory and the energy loss of deuterons and tritons within the target. Experimental results showed an essentially linear dependence on beam current, as expected. A 3.59 ±0.08 power law dependence on the operating voltage was measured, in effective agreement with the modeled value of 3.5. A series of absolute NNS and N-Probe measurements, matched against MCNP calculations, showed that the A3083 and A3082 tubes provide a maximum neutron yield of 8.2 × 10 8 n/s and 4.7 × 10 8 n/s respectively, with estimated uncertainty of ±10%.We showed, through modeling, that tritium decay is not a significant consideration for tubes, such as these, with lifetimes of less than 10 years.
This technical report documents the modeling enhancements, cross-section generation updates, and bench marking with the Shift Monte Carlo code performed under the US Department of Energy Nuclear Energy Advanced Modeling and Simulation Program in FY 2024. The work performed included several modeling enhancements, such as integration of cross-section generation in Titan and the ability to produce microscopic multigroup cross sections with Shift. Benchmarking of the cross sections produced by Shift and the two-step workflow with Griffin was performed for three problems: the Advanced Breeder Test Reactor, a generic pebble bed reactor, and a TRISO heat pipe microreactor. Comparisons of results from these benchmark problems were done with Serpent, OpenMC, and Griffin. These enhancements provide a robust foundation for applying Shift for both reference and two-step neutronics analysis for advanced reactor simulation.
The authors explore performance scalability of the open-source thermal-fluids code, NekRS, on the U.S. Department of Energy's leadership computers, Crusher, Frontier, Summit, Perlmutter, and Polaris. Particular attention is given to analyzing performance and time-to-solution at the strong-scale limit for a target efficiency of 80%, which is typical for production runs on the DOE's high-performance computing systems. Several examples of anomalous behavior are also discussed and analyzed.
Lunchtime presentation requested by Advanced Fuels Campaign colleagues for their visit with visitors from the Japan Atomic Energy Agency. The presentation provides an overview of the BISON fuel performance code and covers a few already published examples related to accident tolerant fuel research using the software.
Argonne National Laboratory and Idaho National Laboratory, through a Department of Energy Gateway for Accelerated Innovation in Nuclear Voucher, supported key analysis needs of Radiant related to (i) air jacket thermal fluid performance, (ii) evaluation of decay heat source terms defining air jacket technical requirements, and (iii) assessment of modeling methodologies employed for core analysis. All of these activities center on numerical simulation of various aspects of Kaleidos using the Multiphysics Object-Oriented Simulation Environment (MOOSE) framework, the Cardinal multiphysics application, the OpenMC Monte Carlo code, and the Nek5000 computational fluid dynamics (CFD) code. This project builds upon an earlier Nuclear Energy Advanced Modeling and Simulation (NEAMS) Thermal-Hydraulic (T/H) Center of Excellence (CoE) project focused on initial demonstration of Cardinal multiphysics simulation of High Temperature Gas Reactors (HTGRs) and now focuses on Radiant’s Kaleidos concept.
Experimental analysis and simulations with the BOUT++ code show that small edge-localised modes (ELMs) in reactor-relevant high-density regimes originate in a region close to the separatrix and only marginally perturb the pedestal structure. The measured divertor peak parallel energy fluence (ε ∥,peak ) for a database of small ELM scenarios in DIII-D and ASDEX Upgrade can be reproduced, within 40 % accuracy on average, if an ad hoc modification of the Eich peak parallel ELM energy fluence model is applied to account for the small ELM pedestal birth location. This allows for first-order extrapolation of small-ELM divertor ε ∥,peak to ITER and SPARC, resulting in values that satisfy the nominal melting threshold of tungsten monoblocks of 12 MJ m −2 . The findings reported in this study, both via modelling and direct measurements, constitute a step forward in assessing small ELMs in high edge-collisionality scenarios as a viable plasma regime for the operation of next-generation fusion machines.
The ASME code qualification effort for Alloy 709 (A709) is currently underway to qualify it for Class A construction in the American Society of Mechanical Engineers (ASME) Boiler and Pressure Vessel Code, Section III, Division 5. The United States (U.S.) Department of Energy (DOE) national laboratories are collaborating in the advanced materials development initiative to investigate the mechanical performance of A709 in support of its code qualification. As part of the A709 code qualification effort, this report summarizes ORNL’s initial experimental findings that support the integration of A709 into the elastic-perfectly plastic (EPP) strain limits code case. It also covers thermomechanical fatigue testing conducted to develop viscoelastic material models, along with the preliminary results of creep-fatigue experiments at 816°C using the Simplified Model Test (SMT) method.
Quantum error-correction codes protect information from realistic noisy channels and lie at the heart of quantum computation and communication tasks. Understanding the optimal performance and other information-theoretic properties, such as the achievable rates, of a given code is crucial, as these factors determine the fundamental limits imposed by the encoding in conjunction with the noise channel. Here, we use the transpose channel to analytically obtain the near-optimal performance of any Gottesman-Kitaev-Preskill (GKP) code under pure loss and pure amplification. We present rigorous connections between GKP code’s near-optimal performance and its dual lattice geometry and average input energy. With no energy constraint, we show that when |𝜏/(1−𝜏)| is an integer, specific families of GKP codes simultaneously achieve the loss and amplification capacity. 𝜏 is the transmissivity (gain) for loss (amplification). Our results establish GKP code as the first structured bosonic code family that achieves the capacity of loss and amplification.
This report summarizes the EMPOWERED Distributed Energy Resources (DER) Permit Accelerator Pilot project. Its primary objective was to streamline the design, installation, permitting, and inspections of targeted Distributed Energy Resource (DER) solutions. Streamlining these processes is critical to strengthening and accelerating the adoption of DER solutions nationwide. Moreover, streamlined processes help expedite the clean energy transition and expand energy resilience by improving understanding and capability among key stakeholders like property owners, code officials, and the general workforce. The project team accomplished this by developing and deploying a set of design and permitting guides for simple DER solutions, which they rolled out via a permitting pilot program. The project team developed nine permitting and inspection guides for simple DER solutions, as follows: 1. Single-Family and Duplex: Electric Vehicle Service Equipment (EVSE), and Storage, Solar + Storage (each in two code cycles) 2. Multi-family and Office: EVSE, Storage, Solar + Storage (each in two code cycles) 3. A permitting process guide to support the implementation of the permitting guides. The guides cover key code requirements via plan review and field inspection checklists for code officials and building owners to follow to streamline the permitting and inspections process for DER solutions. Following the development of the initial set of guides, the project team formed two cohorts, each made up of four jurisdictions, to participate in the Permit Accelerator Pilot. Each cohort received training and technical assistance to support them in incorporating the guides into use. The project successfully engaged jurisdictions in two disparate US regions: Chelsea, Somerville, Natick, and Norwood in metro Boston; and Pima County, Town of Gilbert, Flagstaff, and Sedona in Arizona. Feedback gathered from the project was valuable in refining the guides to better address jurisdictional needs. After the pilot period, the project team developed final versions of the guides based on extensive feedback from the cohorts and from external technical peer review. Feedback included needs and barriers related to implementation, as well as detailed technical improvements from external peer review. The final guides were presented in a series of webinars with a national reach, with 192 attendees from 35 states in the live session of the final webinar, and the final guides are posted online and available for use nationwide. Key findings: • The most receptive jurisdictions were those adjacent to other jurisdictions that had already adopted similar initiatives. The least receptive jurisdictions were those in the midst of adopting new code cycles. • The most common barrier to jurisdictional adoption of new guidelines were internal administrative and process delays, more so than technical barriers. • The main technical barrier to adoption proved to be the concern of fire hazard risks posed by energy storage systems. There is some apprehension in the code enforcement community about the safety of the batteries. Education coordinated with fire and safety services will help alleviate concerns and increase confidence in the acceptance of innovative technologies. • Opportunities to further leverage and advance the guides and related resources include outreach, education, and jurisdictional support; updates to the resources to align with code cycle changes and local requirements; and advancements supporting emerging technologies.
This report describes several bench-scale fluidization experiments that can be used to validate the CFD-DEM method as encapsulated in the MFIX-Exa code. The five cases considered are the cold-flow fluidized beds of Müller et al., Link et al. (spout-fluid), and Goldschmidt et al. (bi-disperse), the hot fluidized bed of Patil et al., and the adsorbing fluidized bed of Li et al. and Janssen. In most cases, MFIX-Exa with “standard” or “typical” CFD-DEM settings, the Gidaspow drag model, and the Gunn heat transfer provide a relatively good prediction of the quantities considered: mean void fraction profiles, mean velocity profiles, fluctuating velocity profiles, mean particle temperature and segregation index. These results, with other verification and validation tests reported elsewhere, contribute to a body of work providing confidence and credibility in CFD predictions from the MFIX-Exa code.