Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Advanced Codes”

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 253 records · Page 14

Particle-in-cell Simulations of Relativistic Magnetic Reconnection with Advanced Maxwell Solver Algorithms

Abstract Relativistic magnetic reconnection is a nonideal plasma process that is a source of nonthermal particle acceleration in many high-energy astrophysical systems. Particle-in-cell (PIC) methods are commonly used for simulating reconnection from first principles. While much progress has been made in understanding the physics of reconnection, especially in 2D, the adoption of advanced algorithms and numerical techniques for efficiently modeling such systems has been limited. With the GPU-accelerated PIC code WarpX, we explore the accuracy and potential performance benefits of two advanced Maxwell solver algorithms: a nonstandard finite-difference scheme (CKC) and an ultrahigh-order pseudo-spectral method (PSATD). We find that, for the relativistic reconnection problem, CKC and PSATD qualitatively and quantitatively match the standard Yee-grid finite-difference method. CKC and PSATD both admit a time step that is 40% longer than that of Yee, resulting in a ∼40% faster time to solution for CKC, but no performance benefit for PSATD when using a current deposition scheme that satisfies Gauss’s law. Relaxing this constraint maintains accuracy and yields a 30% speedup. Unlike Yee and CKC, PSATD is numerically stable at any time step, allowing for a larger time step than with the finite-difference methods. We found that increasing the time step 2.4–3 times over the standard Yee step still yields accurate results, but it only translates to modest performance improvements over CKC, due to the current deposition scheme used with PSATD. Further optimization of this scheme will likely improve the effective performance of PSATD.

79 ASTRONOMY AND ASTROPHYSICS↗

Development and Implementation of a CTF Code Verification Suite.

CTF is a thermal hydraulic subchannel code developed to predict light water reactor (LWR) core behavior. It is a version of Coolant Boiling in Rod Arrays (COBRA) developed by Oak Ridge National Laboratory (ORNL) and North Carolina State University (NCSU) and used in the Consortium for the Advanced Simulation of LWRs (CASL). Here, the existing CTF code verification matrix is expanded, which ensures that the code is a faithful representation of the underlying mathematical model. The suite of code verification tests are mapped to the underlying conservation equations of CTF and significant gaps are addressed. As such, five new problems are incorporated: isokinetic advection, conduction, pressure drop, convection, and pipe boiling. Convergence behavior and numerical errors are quantified for each of the tests and all tests converge at the correct rate to their corresponding analytic solution. A new verification utility that generalizes the code verification process is used to incorporate these problems into the CTF automated test suite.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Coding the Computing Continuum: Fluid Function Execution in Heterogeneous Computing Environments

Advances in network technologies have greatly decreased barriers to accessing physically distributed computers. This newfound accessibility coincides with increasing hardware specialization, creating exciting new opportunities to dispatch workloads to the best resource for a specific purpose, rather than those that are closest or most easily accessible. We present Delta, a service designed to intelligently schedule function-based workloads across a distributed set of heterogeneous computing resources. Delta implements an extensible architecture in which different predictors and scheduling algorithms can be integrated to provide dynamically evolving estimates of function execution times on different resources-estimates that can be used to determine the most appropriate location for execution. We describe predictors for function runtime, data transfer time, and cold-start resource provisioning and configuration delay; dynamic learning methods that update predictor models over time; and scheduling strategies that take into account both function and endpoint information. We show that these methods can halve workload makespan when compared with a strategy that selects the fastest resource, and decrease makespan by a factor of five when compared to a round robin strategy, when deployed on a heterogeneous testbed with resources ranging from a Raspberry Pi to a GPU node in an academic cloud.

Computing continuum↗

Advanced Simulation and Experiments of Strongly Coupled Geomechanics and Flow for Gas Hydrate Deposits: Validation and Field Application. Final Report

The objectives of the proposed research are (1) to investigate geomechanical responses induced by depressurization experimentally and numerically; (2) to enhance the current numerical simulation technology in order to simulate complex physically coupled processes by depressurization and (3) to perform in-depth numerical analyses of two selected potential production test sites: one based on the deposits observed at the Ulleung basin UBGH2-6 site; and the other based on well-characterized accumulations from the westend Prudhoe Bay. To these ends, the recipient will have the following specific objectives: 1). Information obtained from multi-scale experiments previously conducted at the recipient’s research partner (the Korean Institute of Geoscience and Mineral Resources (KIGAM)) that were designed to represent the most promising known Ulleung Basin gas hydrate deposit as drilled at site UBGH2-6 will be evaluated (Task 2). These findings will be further tested by new experimental studies at Lawrence Berkeley National Laboratory (LBNL) and Texas A&M (TAMU) (Task 3) that are designed capture complex coupled physical processes between flow and geomechanics, such as sand production, capillarity, and formation of secondary hydrates. The findings of Tasks 2 and 3 will be used to further improve numerical codes. 2) Develop (in Tasks 4 through 6) an advanced coupled geomechanics and non-isothermal flow simulator (T+MAM) to account for large deformation and strong capillarity. This new code will be validated using data from the literature, from previous work by the project team, and with the results of the proposed experimental studies. The developed simulator will be applied to both Ulleung Basin and Prudhoe Bay sites, effectively addressing complex geomechanical and petrophysical changes induced by depressurization (e.g., frost-heave, strong capillarity, cryo-suction, induced fracturing, and dynamic permeability).

02 PETROLEUM↗

VERA Integration in the NEAMS Workbench

The Virtual Environment for Reactor Applications (VERA) is a cutting edge simulation code suite used to solve complex multiphysics problems for nuclear reactors. VERA makes use of high performance computing (HPC) systems to solve large problems that were not historically tractable. While extensive work has been put into the usability of VERA, it still takes significant effort to build reactor models, ensure their correctness, properly execute the calculations, and conduct the necessary analysis of the results. This difficulty is common for many complex codes such as VERA. As part of the Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, the NEAMS Workbench has been developed and maintained. The goal of the NEAMS Workbench is to provide a single interface with which to execute a wide variety of physics codes, from building the input models to analyzing the results. This document serves as a manual detailing the current integration of VERA with the NEAMS Workbench and a guide for configuring Workbench to allow remote execution of VERA.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Multi-phenomenology Yield Characterization

This report serves as the first delivery of a four-year applied science effort to transform and advance the error bounds for the yield estimate of an explosion. Each year’s delivery will be in this form, culminating in the submission of this work for peer review to a scientific journal. Importantly, the yearly progress reports can then also be viewed as expanding drafts working towards a formal journal article submission. For the first tranche of funding, we collaborated with Air Force Technical Applications Center (AFTAC) scientists to identify unclassified real-world data that demonstrate and validate our advanced error propagation methods. Collaboration includes visits to AFTAC and telecons. For this development, we illustrate the fusion of seismic, acoustic, optical, and surface effect signatures from an explosion. The mathematics and code being adapted to this specific application (Williams et al., 2021) involves physics models of multiple sensor signatures. We have also identified related physics models and have integrated them into code. Current methods of underground explosion yield estimation for the Threshold Test Ban Treaty (TTBT) have served the US treaty monitoring mission well for decades. A research objective of the Defense Nuclear Nonproliferation Research and Development (DNN R&D) office of the National Nuclear Security Administration (NNSA) has always been to provide new technical capabilities for monitoring lower thresholds. The general model and error propagation code to be developed in this project is based on significant advances in error modeling and propagation needed to analyze data at lower detection thresholds. The second tranche of funding for this project began on May 1, 2022, and planned work for the second tranche includes: i) completing the integration of physical model code into the general error model framework; this code accommodates a wide range of linear/nonlinear source models, fixed/ random effects, and frequentist/Bayesian analyses (the purpose of which is not to dictate to users how to analyze data, but instead to allow users the maximum flexibility in their work); ii) illustrative application of code to identified data, and; iii) initial planning with AFTAC researchers on delivery of code to the Common Development Environment at AFTAC, and continued writing of the planned final journal article submission (year two of this progress report), with particular emphasis on descriptions of data identified for this effort.

45 MILITARY TECHNOLOGY, WEAPONRY, AND NATIONAL DEF↗

Report on Preliminary Detailed Experimental Plan for Neutron Irradiation of A709 at ATR and HFIR

Advanced nuclear power technologies will use higher temperatures as a means to extract energy at a higher efficiency than current plants and therefore put a larger demand on the structural materials. Improved performance of structural materials could enable greater safety margins, longer plant lifetimes, and reduce maintenance costs. Alloy 709 is championed as the next generation of austenitic alloys for advanced nuclear reactors. In parallel to the ASME code case pursued, the AMMT program is initiating a neutron irradiation campaign to provide first-of-a-kind engineering data to establish operational design parameters and how the mechanical response is modified by environmental factors. This document refines the AMMT neutron irradiation campaign to a 4 year program to support the generation of creep knockdown factors for Alloy 709 and welded Alloy 709. The campaign is divided among two national laboratories, Oak Ridge National Laboratory and Idaho National Laboratory, to use the strengths of each laboratory. Through a cooperative plan, time-independent properties and time-dependent properties will be obtained across a large temperature window, nominally 300°C to 800°C, damage levels up to 10 dpa, and with and without the impacts of transmutation produced helium.

99 GENERAL AND MISCELLANEOUS↗

Quantitative metrics for evaluating thermonuclear design codes and physics models applied to the National Ignition Campaign

In this study, physics models and design codes for hot dense plasmas undergoing thermonuclear burn are evaluated objectively using statistical metrics that compare the difference between calculations and data relative to the experimental uncertainties. The analysis is applied to the National Ignition Campaign (NIC) because it is relevant, comprehensive, and well documented. The statistics confirm that a key process afflicting NIC performance is mix driven by hydrodynamic instabilities as approximated here using the KL model [G. Dimonte and R. Tipton, Phys. Fluids 18, 085101 (2006)]. New physics models are also presented for instability-driven magnetic fields [B. Srinivasan et al., Phys. Rev. Lett. 108, 165002 (2012)] and the Coulomb logarithm for electron–ion thermal relaxation [G. Dimonte and J. Daligault, Phys. Rev. Lett. 101, 135001 (2008)]. The plasma-generated magnetic fields improve the agreement between code and data in a statistically significant manner by reducing the electron thermal conduction in the hot-spot via the Hall term. The Coulomb logarithm is presented in a numerically practical form that incorporates recent theoretical advances. However, it does not improve the agreement between the code and data because the thermal relaxation is so fast in non-ignited NIC plasmas that the typical 20% error does not change the tight coupling between electrons and ions. Even with these improved physics models, the one-dimensional simulations presented here are not able to describe the high-convergence NIC implosions in a statistically acceptable manner. Of the heroic multi-dimensional simulations of Clark et al. [Phys. Plasmas 23, 056302 (2016)], only those in three-dimensions satisfy the statistical acceptance criterion.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Development of an Advanced Multiphysics Simulation Capability for Radiant's Microreactor Design

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.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Elevated-temperature cyclic properties of advanced manufactured materials

Microreactor developers at the 2019 GAIN Microreactors workshop expressed a need for advanced manufacturing to fabricate microreactor components. Powder metallurgy (PM) hot isostatic pressing (HIP) is more mature than other advanced manufacturing techniques for the following reasons: 1. a code case exists for Section III, Division 1, Subsection NB components of light water reactors to be manufactured using PM HIP Alloy 316L, and 2. tensile and creep properties have been shown to be equivalent to wrought material for Type 316 stainless steel and Grade 91. Thus, PM HIP is more readily deployable in the near term. Elevated-temperature cyclic properties of PM-HIP materials, however, have not been established. The scope of testing for Code qualification could be significantly reduced compared to a new material for the following reasons: 1. only a limited amount of time consuming creep testing would be required, and 2. testing of fatigue and creep-fatigue behavior could be reduced if it can be demonstrated that the properties are comparable to wrought material. The objective of this work package is to develop an understanding of the relationship of advanced materials processing on the material characteristics as they relate to elevated-temperature component design and construction. This will be achieved by the following: • Elevated-temperature fatigue and creep-fatigue testing of PM HIP and wrought Alloy 316L. • Fracture and microstructural characterization of the PM HIP and wrought Alloy 316L. • Analysis of the elevated-temperature fatigue and creep-fatigue properties of PM HIP Alloy 316L. • Draft a report evaluating the elevated-temperature fatigue and creep-fatigue properties of PM HIP Alloy 316L. • Pursue procurement of Alloy 316H and Grade 91 manufactured using PM HIP. • Initiate scoping studies of the procured Alloy 316H or Grade 91 manufactured using PM HIP contingent upon successful material procurement.

36 MATERIALS SCIENCE↗

NuGraph2 with context-aware inputs: physics-inspired improvements in semantic segmentation

Graph neural networks have recently shown strong promise for event reconstruction tasks in Liquid Argon Time Projection Chambers, yet their performance remains limited for underrepresented classes of particles, such as Michel electrons. In this work, we investigate physics-informed strategies to improve semantic segmentation within the NuGraph2 architecture. We explore three complementary approaches: (i) enriching the input representation with context-aware features derived from detector geometry and track continuity, (ii) introducing auxiliary decoders to capture class-level correlations, and (iii) incorporating energy-based regularization terms motivated by Michel electron energy distributions. Experiments on MicroBooNE public datasets show that physics-inspired feature augmentation yields the largest gains, particularly boosting Michel electron precision and recall by disentangling overlapping latent space regions. In contrast, auxiliary decoders and energy-regularization terms provided limited improvements, partly due to the hit-level nature of NuGraph2, which lacks explicit particle- or event-level representations. Our findings highlight that embedding physics context directly into node-level inputs is more effective than imposing task-specific auxiliary losses, and suggest that future hierarchical architectures such as NuGraph3, with explicit particle- and event-level reasoning, will provide a more natural setting for advanced decoders and physics-based regularization. The code for this work is publicly available on Github at https://github.com/vitorgrizzi/nugraph_phys/tree/main_phys.

Other Experiments↗

Modeling, analysis, and code/data validation of DIII-D tokamak divertor experiments on ELM and non-ELM plasma tungsten sputtering erosion

Abstract We analyzed recent DIII-D tokamak tungsten divertor probe experiments using advanced, coupled, sputter erosion/redeposition, plasma, and surface response code packages. Modeling is done for ELMing H-mode, and L-mode plasmas, impinging on various size tungsten deposits on Divertor Material Evaluation System (DiMES) carbon probes. The simulations compute 3D, full kinetic, sub-gyromotion, impurity sputtering and transport, including changes in tungsten surface composition and response due to mixed deuterium and carbon ions irradiation. Per our analysis, ELM (edge localized mode) plasma sputtering in DIII-D mostly involves free-streaming high energy (∼500–1000 eV) D + and C +6 ions, with high near-surface plasma density. L-Mode sputtering is due to impurity sputtering (C, W) only, with lower density. All cases show complete redeposition of tungsten on the divertor, with significant redeposition on the tungsten spots themselves, and low self-sputtering. Comparison of ELM plasma gross tungsten erosion simulation results with in-situ spectroscopic data is good, as are code/data comparisons of net erosion using post-exposure Rutherford backscattering (RBS) data for the L-mode probes. The analysis, extrapolated to a full tungsten divertor, implies low net erosion and negligible plasma contamination from sputtering. These results support the use of high-Z plasma facing surfaces in ITER and beyond.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Fault-tolerant optical interconnects for neutral-atom arrays

We analyze the use of photonic links to enable large-scale fault-tolerant connectivity of locally error-corrected modules based on neutral atom arrays. Our approach makes use of recent theoretical results showing the robustness of surface codes to boundary noise and combines recent experimental advances in atom-array quantum computing with logical qubits with optical quantum networking techniques. We find the conditions for fault tolerance can be achieved with local two-qubit Rydberg gate and nonlocal Bell-pair errors below 1% and 10%, respectively, without requiring distillation or space-time overheads. Realizing the interconnects with a lens, a single optical cavity, or an array of cavities enables—with sufficient multiplexing—a Bell-pair generation rate in the 1–50 MHz range. When directly interfacing logical qubits, this rate translates to error-correction cycles in the 25–2000 kHz range, satisfying all requirements for fault tolerance and in the upper range fast enough for 100 kHz logical clock cycles. Published by the American Physical Society 2025

Sinclair, Josiah (ORCID:0000000215238295)↗

Sockeye Validation Support Using the SPHERE Facility

The engineering-scale heat pipe application Sockeye, being developed through the Nuclear En- ergy Advanced Modeling and Simulation program is concurrently undergoing code development and verification and validation activities. The Single Primary Heat Extraction and Removal Emu- lator (SPHERE) facility at Idaho National Laboratory has performed numerous sodium heat pipe experiments, one of which is modeled in this work to support Sockeye validation. Heat pipe data produced by the SPHERE facility was found to be very useful not only for Sockeye validation but also as a useful resource for model development. A Sockeye simulation performed using Sockeye’s conduction model modeled an experimental run at SPHERE and was found to be able to accurately reproduce the steady temperature distribution but revealed discrepancies in transient behavior after shutoff of the power. This effort also revealed a number of challenges of modeling an experimental setup such as SPHERE, for example, the importance of details in the setup such as possible gaps in insulation. Also, continuous collaboration between Sockeye developers and SPHERE experimentalists has been established, which allows future experiments to be refined in response to questions and challenges found by simulating earlier experiments.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Preliminary Study on TRISO Fuel Cross Section Generation

Cross section self-shielding methodologies for TRISO fuel were assessed to provide accurate multigroup cross sections for a high-fidelity reactor physics code so that the code is able to accurately model and simulate advanced reactors with TRISO fuel. Initially, the two existing methodologies (the SCALE method and the Sanchez-Pomraning method) were studied and implemented to MC2-3 for detailed performance tests. Additionally, a new spatial self-shielding method, named the iterative local spatial self-shielding (ILSS) method, for particulate fuels was developed based on the disadvantage factor and implemented to MC2-3 as well. The new method approximately accounts for the effect of randomly distributed particles on the particle shadowing effect using a homogenized compact region surrounding a particle of interest at the center. The self-shielded cross sections of the particle at the center are determined iteratively since the cross sections of the homogenized compact region are calculated using them. For the energy range above 100 keV where the fuel-to-moderator ratio is more important than the random distribution of particles, a single particle unit-cell model is used by preserving the average amount of moderator per fuel particle in the system. The three self-shielding methods implemented in MC2-3 were tested using numerical benchmark problems made based on fuel compact problems of a prismatic-type very high temperature reactor. Test results indicated that the ILSS method produced slightly better results than the SCALE and Sanchez-Pomraning methods, compared to the Serpent-2 Monte Carlo results obtained with 25 independent random particle configurations. The SCALE and Sanchez-Pomraning methods tend to underestimate the heterogeneity effect by 150 and 100 pcm, respectively, while the new ILSS method overestimates the heterogeneity effect by 70 pcm. In future, the new self-shielding method will be extended to perform pebble calculations and compare results with those from the SCALE and Sanchez-Pomraning methods. Furthermore, the new method will be optimized for practical applications to on-the-fly resonance treatment for lattice or whole-core calculations for advanced reactors.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

Modeling of Prismatic High Temperature Reactors in Pronghorn

Pronghorn is a Multiphysics Object-Oriented Simulation Environment (MOOSE) based thermal-hydraulics code developed at Idaho National Laboratory (INL) for advanced reactor analysis. It has been previously applied to model Pebble-Bed High Temperature Reactors (High Temperature Reactor (HTR)s), Liquid-Metal Cooled Reactors, and Molten Salt Reactors, among others. This work applies the coarse-mesh thermal hydraulics capabilities in Pronghorn to model Prismatic-Core HTRs. In particular, the Oregon State University (OSU)’s High Temperature Test Facility (HTTF) is modeled with Pronghorn. The HTTF is a 1:4 height scaled-down facility of General Atomics’ Modular High Temperature Gas-cooled Reactor.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

NREL Hydrogen Sensor Testing Laboratory and Smart Hydrogen Wide Area Monitoring

The NREL Hydrogen sensor laboratory aims to ensure that hydrogen sensor technology is available to meet end-user needs and to foster the proper use of sensors by advancing next generation sensing and analysis techniques, supporting codes and standards development, and improving component reliability systems.

08 HYDROGEN↗

Report on Preliminary Detailed Experimental Plan for Neutron Irradiation of A709 at ATR and HFIR

Advanced nuclear power technologies will use higher temperatures as a means to extract energy at a higher efficiency than current plants and therefore put a larger demand on the structural materials. Improved performance of structural materials could enable greater safety margins, longer plant lifetimes, and reduce maintenance costs. Alloy 709 is championed as the next generation of austenitic alloys for advanced nuclear reactors. In parallel to the ASME code case pursued, the AMMT program is initiating a neutron irradiation campaign to provide first-of-a-kind engineering data to establish operational design parameters and how the mechanical response is modified by environmental factors. This document refines the AMMT neutron irradiation campaign to a 4 year program to support the generation of creep knockdown factors for Alloy 709 and welded Alloy 709. The campaign is divided among two national laboratories, Oak Ridge National Laboratory and Idaho National Laboratory, to use the strengths of each laboratory. Through a cooperative plan, time-independent properties and time-dependent properties will be obtained across a large temperature window, nominally 300°C to 800°C, damage levels up to 10 dpa, and with and without the impacts of transmutation produced helium.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗