Search NASA⌕ Search

SEARCH · Search NASA

Results for “Code Assessment”

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 397 records · Page 22
All sourcesNASA NTRSDOE OSTI
DOE OSTI2016

Determination of the Spectral Index in the Fission Spectrum Energy Regime

Neutron reaction cross sections play a vital role in tracking the production and destruction of isotopes exposed to neutron fluence. They are central to the process of reconciling the initial and final atom inventories. Measurements of irradiated samples by radiochemical methods in tangent with an algorithm are used to evaluate the fluence a sample is exposed to over the course of the irradiation. This algorithm is the Isotope Production Code (IPC) created and used by the radiochemistry data assessment team at Los Alamos National Laboratory (LANL). An integral result is calculated by varying the total neutron fluence seen by a sample. A sample, irradiated in a critical assembly, will be exposed to a unique neutron flux defined by the neutron source and distance of the sample from the source. Neutron cross sections utilized are a function of the hardness of the neutron spectrum at the location of irradiation. A spectral index is used an indicator of the hardness of the neutron spectrum. Cross sections fit forms applied in IPC are collapsed from a LANL 30-group energy structure. Several decades of research and development have been performed to formalize the current IPC cross section library. Basis of the current fission spectrum neutron reaction cross section library is rooted in critical assembly experiments performed from the 1950’s through the early 1970’s at LANL. The focus of this report is development of the spectral index used an indicator of the hardness of the neutron spectrum in the fission spectrum energy regime.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2019

Nuclear Data and Benchmarking Program: Nuclear Data Performance Assessment for Advanced Reactors

Over the last several Evaluated Nuclear Data File (ENDF)/B releases, many updates have been introduced to the nuclear data libraries that can significantly change computational results. These changes can be particularly important for advanced reactor concepts, which are not as thoroughly investigated as those used in light water reactor (LWR) systems. Therefore, a performance assessment of the ENDF/B libraries using the SCALE 6.2.3 code package was conducted under the auspices of the U.S., Department of Energy Office of Nuclear Energy (DOE-NE). The following systems relevant to the advanced reactor community were chosen as subject for investigation in this report: sodium cooled fast reactor systems, graphite moderated high temperature gas-cooled reactors, and several molten salt reactor models. At first, the similarity of the integral performance of the ENDF/B-VII.1 and ENDF/B-VIII.0 nuclear data libraries caused by compensating errors in important differential data was demonstrated. In calculations for a shipping container of high-assay low enriched uranium, the data of individual isotopes were systematically swapped between the two ENDF/B libraries. An eigenvalue difference of up to 450 pcm was found due to the use of 235 U and 238 U from one library and 1 H and 16 O from the other library. This clearly demonstrates a cross-correlation between reaction data sets of different isotopes within a library that should be reported in the evaluations. For sodium-cooled fast reactor (SFR) fuel assemblies, significant eigenvalue (k inf ) differences (200– 450 pcm) were found between calculations using the 2011 ENDF/B-VII.1 data and 2018 ENDF/B-VIII.0 data. These differences were mainly caused by updates of 238 U and 239 Pu neutron cross sections. Results from the multigroup (MG) calculations further revealed that the group structure of the applied MG library strongly influences the MG bias due to the importance of an appropriate energy resolution of the resonances in higher energy ranges. The application of previously used MG libraries requires new verification for advanced reactor simulation in comparisons with reference continuous-energy calculations. Nuclear data uncertainty analyses of these SFR systems resulted in eigenvalue uncertainties between 1,400 and 1,800 pcm, which are three to four times higher than corresponding uncertainties in light water reactor (LWR) systems. The main contributor to this uncertainty was found to be inelastic scattering on 238 U, which shows an uncertainty of up to 50% in the fast energy range. Other important contributors are the scattering reactions of 56 Fe and 23 Na which have so far not appeared in LWR analysis. Since the largest contribution to the eigenvalue uncertainty of SFR systems is coming from scattering reactions, it is expected that uncertainties in the angular scattering distributions also have a significant impact. Those uncertainties are, however, not available for the majority of nuclides, and the capabilities to determine sensitivities to this data are currently underdeveloped. For two graphite moderated high temperature gas-cooled reactor (HTGR) benchmarks, only the ENDF/BVII. 1 calculations resulted in consistent eigenvalues (k eff ) with the corresponding experimental measurements. Eigenvalue differences of about 1,000 pcm were found between the 2006 ENDF/B-VII.0 data and the 2011 ENDF/B-VII.1 data, primarily due to an updated carbon capture cross section in the thermal energy range in the later ENDF/B release. The ENDF/B-VIII.0 eigenvalue was larger than the ENDF/B-VII.1 eigenvalue by about 300 pcm mainly due to updates in multiple 235 U neutron cross sections, with offsetting updates in the 238 U cross sections. With ENDF/B-VIII.0, graphite can be modeled as perfect crystal or with two different porosities. The choice of the graphite evaluation can have a significant influence of the eigenvalue. For a HTGR benchmark, maximum eigenvalue differences as high as 650 pcm due to different porosities were observed. More detailed studies of the impact of the porosity are necessary while considering that the graphite porosity can vary between the used materials and changes as a function of neutron fluence. xii The uncertainty of the HTGR eigenvalues due to nuclear data uncertainties was found between 500 and 600 pcm, with the top contributor being the neutron multiplicity of 235 U. A gap in the form of missing uncertainties in graphite thermal scattering data, that might have significant impact on HTGR reactors, was identified. This gap is being investigated by the DOE-NE Nuclear Data and Benchmarking Program, as documented in a separate report The fast spectrum molten salt reactor calculations revealed similar comparisons as for the SFR assemblies. The eigenvalue uncertainties were significantly influenced by inelastic scattering on 238 U, while the impact of angular scattering distributions is unknown. For the graphite and zirconium hydride moderated system, eigenvalue uncertainties of up to 700 pcm were observed, requiring in depth analyses of the cross sections and corresponding uncertainties of the included materials. In general, it was observed that the flux spectra of the various molten salt reactor systems show significant differences, even between a fresh state and a depleted state. The choice of the energy group structure of MG calculations is therefore highly relevant; the applicability of previously used MG libraries needs to be verified. Due the unavailability of pin power measurements of advanced reactor systems, pin power calculations with both ENDF/B-VII.1 and VIII.0 data were performed for a light water reactor and compared to corresponding measurements. Both calculations show good agreement with the measurements. A comparison between the ENDF/B-VII.1 and VIII.0 results revealed small differences, mostly in the range of about 0.2%.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2019

Annual Energy Savings and Peak Power Reduction Simulated for CCHP Design and Baseline HP in Cold Climate Cities

This report assesses energy savings and peak power reduction of three heat pumps in multiple US cold climate cities and building types. The three heat pumps comprise a baseline heat pump and two high efficiency cold climate heat pumps (CCHP). The DOE/ORNL Heat Pump Design Model was used to simulate the heat pumps and generate performance curves for EnergyPlus building energy simulations. The annual energy simulations were performed over seven cold climates: 1) Minneapolis, MN; 2) Duluth, MN; 3) Fairbanks, AK; 4) Helena, MT; 5) Indianapolis, IN; 6) Peoria, IL and 7) Eagle County, CO. In order to assess impact of residential building types matching different energy codes, we simulated six single-family houses in Indianapolis, IN, including : 1) single-family house having a crawlspace; 2) having a slab foundation; 3) having a heated basement; 4) having an unheated basement; 5) one single-family house having a heated basement matching IECC code of 2009, 6) and one single-family house having a heated basement matching IECC code of 2012. The former 4 houses match the IECC 2006 code.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

Initial development of a high temperature life prediction method directly accounting for variability in material properties

This report presents the initial development of a life prediction method that accounts for the variability of the material properties in Grade 91 steel. We account for material variability by both fitting a variable 3-parameter Weibull distribution to the rupture data of the material and by quantifying the effect that the creep deformation properties have on stresses via a Monte Carlo approach. To ensure a reasonable time frame when applying Monte Carlo to the creep stress, we propose a methodology that models the material as an extremely viscous Stokes flow with a non-Newtonian viscosity, therefore solving the problem in a single step. We use the final probabilistic model for structural failure developed in this document to evaluate a Flat Head Vessel designed with the Section III, Division 5 ASME Boiler Pressure and Vessel Code rules with a 100,000h design life and assess its probability of premature failure.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2020

Summary Report for the NEPA Impact Analysis. Revision 1

The Waste Isolation Pilot Plant (WIPP), located in southeastern New Mexico, has been developed by the U.S. Department of Energy (DOE) for the geologic (deep underground) disposal of defense-related transuranic (TRU) waste. Containment of TRU waste at the WIPP facility is derived from standards set forth in Title 40 of the Code of Federal Regulations (CFR), Part 191. The DOE assesses compliance with the containment standards according to the Certification Criteria in Title 40 CFR Part 194 by means of Performance Assessment (PA) calculations performed by Sandia National Laboratories (SNL). WIPP PA calculations estimate the probability of radionuclide releases from the repository to the accessible environment for a regulatory period of 10,000 years after facility closure. The DOE Carlsbad Field Office (CBFO) has initiated a National Environmental Policy Act (NEPA) action for a proposal to excavate and use additional transuranic (TRU) waste disposal panels at the WIPP facility. This report documents an analysis undertaken as part of an effort to evaluate the potential environmental consequences of the proposed action. Although not explicitly required for a NEPA analysis, evaluations of a dose indicator to hypothetical members of the public after final facility closure are presented in this report. The analysis is carried out in two stages: first, Performance Assessment (PA) calculations quantify the potential releases to the accessible environment over a 10,000-year post-closure period. Second, dose was evaluated for three hypothetical exposure pathways using the conservative radionuclide concentrations assumed to be released to the accessible environment.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2020

Managing Solar Photovoltaic Integration in the Western United States: Power System Flexibility Requirements and Supply

As penetrations of variable renewable energy generation technologies such as wind and solar photovoltaics (PV) continue to increase across the United States, greater uncertainty and variability in the net load often lead to a concern about how power systems may adapt. Managing the system net load (i.e., load minus contribution from variable generation technologies) may become more challenging with increasing variable generation, as the magnitude and frequency of net ramps increase. However, there is inherent flexibility in power systems through the conventional generator fleet (under least-cost unit commitment and economic dispatch), less-conventional generation sources (e.g., storage, demand response, concentrating solar power with thermal energy storage), and imports and exports with neighbors. In this analysis, we create an open-source tool to analyze the flexibility of the results of a specific commercial unit commitment and economic dispatch tool (PLEXOS), but the code can be applied generically as well. The tool assesses the flexibility requirements (or demand) of a system through a net load analysis. The constraints and limitations of each generator are then considered to determine the availability (or supply) of flexibility. Then, the supply and demand of flexibility are compared to gain a more complete picture of potential flexibility concerns. We apply this open-source tool to high-penetration PV scenarios constructed for three focus regions in the western United States defined using the Resource Planning Model (RPM) capacity expansion modeling tool: RPM-OR, RPM-CO, and RPM-AZ. Generally, we find few flexibility concerns, as the western United States represents a large and interconnected power system with significant inherent flexibility. In addition, the PV scenarios we analyzed are overbuilt on capacity, leaving plenty of ramping ability on the system. We do find that for each focus region, the impact of imports on meeting ramping needs is essential. This means the PV integration in each focus region impacts the entire rest of the system. Each system has different dominant sources of flexibility. The conventional generator fleet (especially coal and gas combined-cycle technologies) as well as less-conventional sources such as storage are all shown to be important sources of flexibility. The scenarios evaluated here were designed to study the planning and operations impact of high solar penetration in each of three focus regions. However, none of the three focus regions likely will deploy PV in isolation, meaning the ability of imports and exports to provide flexibility may be considerably different in scenarios with strong PV deployment in every region. Overall, we intend that the framework we present here will be useful in future analysis of other system evolutions to identify whether and how flexibility may constrain the successful deployment of variable generation technologies.

14 SOLAR ENERGY↗
DOE OSTI2016

Economic Model For GDP And Tangible-Assets Loss Estimation In The MACCS Offsite Consequence Analysis Code

The MELCOR Accident Consequence Code System (MACCS) code is the NRC code used to perform probabilistic health and economic consequence assessments for atmospheric releases of radionuclides. MACCS is used by U.S. nuclear power plant license renewal applicants to support the plant specific evaluation of severe accident mitigation alternatives (SAMA) analyses as part of an applicant's environmental report for license renewal. MACCS is also used in severe accident mitigation design alternatives (SAMDA) and severe accident consequence analyses for environmental impact statements (EISs) for new reactor applications. The NRC uses MACCS in its cost-benefit assessments supporting regulatory analyses that evaluate potential new regulatory requirements for nuclear power plants. NRC regulatory analysis guidelines recommend the use of MACCS to estimate the averted "offsite property damage" cost (benefit) and the averted offsite dose cost elements.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2022

Validation of energetic particle transport models for time-dependent integrated simulations of burning plasmas (Final ST Report)

The major goal of this project is to deliver a suite of validated energetic particle (EP) modeling tools for time-dependent integrated tokamak simulations (with the code TRANSP) that can be used for modeling plasma discharges featuring a substantial EP content, such as burning plasmas. During FY20, issues were resolved with an energetic particle transport model called TGLF-EP, with improvements made to the TGLF-EP workflow, generation of input files, and interpretation of runs on the NERSC supercomputing cluster. The recipient also gave input on simulations with the code TRANSP to reproduce the reference JET discharge, and assessed the accuracy and uncertainties of the simulated B field, current evolution, thermal plasma profiles, stored energy, etc. Input was also given on interpretation of fast-ion instability calculations and recommendations for modifications to JET afterglow discharge for the upcoming DT campaign.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2024

Developing And Scaling an OpenFOAM Model to Study Turbulent Flow in a HFIR Coolant Channel

Improving the understanding of how computational fluid dynamics (CFD) direct numerical simulations (DNS) of flows in the High Flux Isotope Reactor (HFIR) perform when run in parallel using the high performance computing (HPC) platform Summit at the Oak Ridge Leadership Computing Facility (OLCF) is of particular importance to boost the computational tools used to support HFIR conversion to low enriched fuel (LEU). Evaluation of scaling performance was driven by the increasing importance of graphics processing unit (GPU) usage in HPC, which is becoming the standard for modern supercomputers such as Summit. The desired results are to obtain a strong positive correlation between the computational resources dedicated to a problem and the relative speed-up of the simulation in comparison to a benchmark. This capability will allow substantially improvement in HFIR flow analytical capabilities, specifically when predicting turbulence properties at high Reynolds numbers. The study leverages previous simulation results performed with code PHASTA (finite element) on HPC platforms Cori (NERSC) and Theta (ALCF) [1] with computing options provided in the computing platform OpenFOAM (finite volume) at OLCF. Transitioning from PHASTA to OpenFOAM will (1) eliminate dependence on third-party software for mesh generation and manipulation, (2) reduce resource needs by employing modern architectures, and (3) build expertise for future modeling of HFIR-specific problems like heat transfer in involute geometry, entrance effects, flow structure in channel corners, and so on—all important issues when defining the available thermal margins in the transition to LEU. CPUs and GPUs differ significantly in their architecture and utilization, as discussed in the literature [2]. The most important differences are in the approach to computations and their memory. A single GPU contains a large quantity of cores, enabling it to perform with a much higher throughput than a CPU, but execution requires a different approach. GPU codes execute instructions using the Single-Instruction Multiple-Thread (SIMT) approach in which a single instruction is used for groups of threads called warps. A warp typically consists of 32 threads which must execute the same set of instructions, although on separate threads. Alternately, a CPU has far fewer cores that are much more flexible in their operation, excelling at quickly performing more complex serial computations. This is why GPUs have greater throughput when properly utilized. The second important difference is seen when comparing their memory spaces. Limited memory allocations and CPU–GPU communications cause a significant bottleneck in GPU-accelerated programs. Further study was required to properly take advantage of GPU resources. A comprehensive analysis of code performance and the model-specific features of turbulence constitutes the core of this work. In this study, a DNS simulation of HFIR channel turbulence was performed with the finite volume CFD code OpenFOAM v2112 and CUDA v11.0 on Red Hat Enterprise Linux v8.2. The OpenFOAM installation had AMGx integrated to enable GPU acceleration and utilizes the PETSc4FOAM library. The computational resources and the problem size were scaled on CPU and CPU + GPU architectures to gain a better understanding of the performance of a DNS problem on modern computing hardware. The study aimed to analyze the scaling of the code exclusively on CPUs and then to examine the scaling of the codes with GPU acceleration enabled. Scaling studies included CPU and GPU acceleration on a mesh of varying resolution to analyze the impact of problem size relative to computational resources. In the course of preparing the GPU configuration on Summit, mainly using the AMGX solvers, difficulties were encountered stemming from constant changes resulting from extensive ongoing development activities and the changing environment. This resulted in the inability to complete the GPU portion of the work. The code was compiled and tested, but production runs to assess acceleration were not performed because the used discretional compute time allocation expired as year-end approached. The Summit HPC platform is scheduled for decommissioning in 2024, making it unattractive for future use with Nvidia-based GPUs. Therefore, the work will be moved onto NERSC machines in FY24. An application was prepared and submitted, and sufficient node-hours were awarded to continue the research in the next calendar year. This report summarizes work performed thus far, which mostly focused on CPU OpenFOAM computing.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2025

Pacific Northwest National Laboratory Facility Radionuclide Emission Points and Sampling Systems

Battelle–Pacific Northwest Division operates numerous research and development laboratories in Washington State. The U.S. Department of Energy (DOE) contracts to Battelle at Richland facilities on both the DOE Hanford Site and the Pacific Northwest National Laboratory (PNNL) Richland campus. These facilities have the potential for radionuclide air emissions. The PNNL contract with DOE also includes operations at the PNNL-Sequim campus in Sequim, where there is also the potential for radionuclide air emissions. This document is a periodic update that describes current PNNL facility emission units and sampling systems. The National Emission Standard for Hazardous Air Pollutants (NESHAP [40 Code of Federal Regulations 61, Subpart H]) requires an assessment of all emission units that have the potential for radionuclide air emissions. Emission units are registered with the State of Washington. Potential emissions from emission units are assessed annually by PNNL staff. Sampling, monitoring, and other regulatory compliance requirements are designated based on the potential to-emit dose criteria, a graded approach to facility-identified potential impact categories, and regulatory requirements. The purpose of this document is to describe the facility radionuclide air emission sampling program and provide current and historical facility emission unit system performance, operation, and design information. For sampled emission units, the building, exhaust unit, control technologies, and sample extraction details are provided. Additionally, applicable configuration drawings, figures, and photographs are included. For non-sampled emission units, emission estimation and radionuclide source details are provided. Site-wide permits for the lowest potential impact category are described. Deregistered/transitioned emission unit details are also provided as necessary for at least 5 years post-closure/transition. Currently, five emission units are sampled continuously for particulate radionuclides at PNNL managed facilities on the PNNL-Richland campus (3 of the 5) and on the Hanford Site (2 of the 5). Four of these units have sampling systems that comply with the American National Standards Institute/Health Physics Society (ANSI/HPS) N13.1–2011 standard for sampling from stacks and ducts of nuclear facilities, and the fifth is grandfathered and compliant with the older ANSI N13.1–1969 standard. In addition, the PNNL-managed Hanford Site 325 Building EP 325-01-S stack is sampled continuously for emissions of tritium. No emissions sampling is required for the single licensed emission unit on the PNNL-Sequim campus.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Magnetic Flux Rope Identification and Characterization from Observationally Driven Solar Coronal Models

Formed through magnetic field shearing and reconnection in the solar corona, magnetic flux ropes are structures of twisted magnetic field, threaded along an axis. Their evolution and potential eruption are of great importance for space weather. Here we describe a new methodology for the automated detection of flux ropes in simulated magnetic fields, utilizing field-line helicity. Our Flux Rope Detection and Organization (FRoDO) code, which measures the magnetic flux and helicity content of pre-erupting flux ropes over time, as well as detecting eruptions, is publicly available. As a first demonstration, the code is applied to the output from a time-dependent magnetofrictional model, spanning 1996 June 15–2014 February 10. Over this period, 1561 erupting and 2099 non-erupting magnetic flux ropes are detected, tracked, and characterized. For this particular model data, erupting flux ropes have a mean net helicity magnitude of 2.66×10{sup 43} Mx{sup 2}, while non-erupting flux ropes have a significantly lower mean of 4.04×10{sup 42} Mx{sup 2}, although there is overlap between the two distributions. Similarly, the mean unsigned magnetic flux for erupting flux ropes is 4.04×10{sup 21} Mx, significantly higher than the mean value of 7.05×10{sup 20} Mx for non-erupting ropes. These values for erupting flux ropes are within the broad range expected from observational and theoretical estimates, although the eruption rate in this particular model is lower than that of observed coronal mass ejections. In the future, the FRoDO code will prove to be a valuable tool for assessing the performance of different non-potential coronal simulations and comparing them with observations.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2022

Quantitative Performance Assessment of Proxy Apps and Parents (Report for ECP Proxy App Project Milestone ADCD-504-28)

The ECP Proxy Application Project has an annual milestone to assess the state of ECP proxy applications and their role in the overall ECP ecosystem. Our FY22 March/April milestone (ADCD- 504-28) proposed to: Assess the fidelity of proxy applications compared to their respective parents in terms of kernel and I/O behavior, and predictability. Similarity techniques will be applied for quantitative comparison of proxy/parent kernel behavior. MACSio evaluation will continue and support for OpenPMD backends will be explored. The execution time predictability of proxy apps with respect to their parents will be explored through a carefully designed scaling study and code comparisons. Note that in this FY, we also have quantitative assessment milestones that are due in September and are, therefore, not included in the description above or in this report. Another report on these deliverables will be generated and submitted upon completion of these milestones. To satisfy this milestone, the following specific tasks were completed: Study the ability of MACSio to represent I/O workloads of adaptive mesh codes. Re-define the performance counter groups for contemporary Intel and IBM platforms to better match specific hardware components and to better align across platforms (make cross-platform comparison more accurate). Perform cosine similarity study based on the new performance counter groups on the Intel and IBM P9 platforms. Perform detailed analysis of performance counter data to accurately average and align the data to maintain phases across all executions and develop methods to reduce the set of collected performance counters used in cosine similarity analysis. Apply a quantitative similarity comparison between proxy and parent CPU kernels. Perform scaling studies to understand the accuracy of predictability of the parent performance using its respective proxy application. This report presents highlights of these efforts.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

High-Temperature Gas-Cooled Pebble-Bed Reactors Running In And Transient Modeling Capabilities Demonstration

This study presents a comprehensive benchmarking and verification effort of several thermal-hydraulic and multiphysics capabilities for high-temperature gas-cooled reactor (HTGR) applications. The first part of this effort focuses on the running-in verification of Griffin's multiphysics capabilities, specifically for simulating the evolution of Pebble Bed reactor cores from startup to equilibrium. In the absence of validation data, code-to-code comparisons are conducted with Kugelpy, showing good agreement for key quantities like maximum power density and fresh core k-eigenvalue predictions. However, discrepancies in equilibrium core predictions suggest potential issues with cross sections, underscoring the need for further refinement and evaluation. The HTTF system analysis code benchmark involves RELAP5-3D, SAM, and GAMMA+ to assess their predictive capabilities for HTTF behavior under both normal operation and pressurized conduction cooldown (PCC) transient conditions. While there is good agreement in predicting major parameters such as coolant temperature, solid temperature, and flow distribution, discrepancies in transient behavior highlight differences in modeling approaches, nodalizations, and heat transfer models. The HTTF lower plenum CFD benchmark employs nekRS to simulate flow mixing phenomena, successfully capturing relevant flow physics and demonstrating mesh independence in complex geometries. Preliminary results suggest a relatively uniform temperature field but significant unsteadiness in the flow, requiring time-averaging analyses. The GPBR200 system analysis code benchmark uses SAM's core channel and porous media models, incorporating an RCCS loop for decay heat removal. During steady-state and transient conditions, including protected de-pressurized and pressurized loss of forced cooling (DLOFC and PLOFC), both models show good agreement in predicting temperature profiles and key parameters. Notably, while the core channel model underpredicts convective heat transfer effects, both models maintain temperatures well below the TRISO fuel safety limit. These benchmarking efforts collectively enhance the predictive capabilities of the tools used in HTGR design and safety analysis, guiding developments to improve their accuracy and applicability.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2024

Large-scale Hydrogen Storage Risk Assessment

This project investigated risks involved in deploying a large-scale hydrogen storage system at the Port of Seattle (hereafter, the Port) for its on-terminal and maritime applications in an urban industrial setting. Alongside, the project attempted to address some of the barriers to risk assessment such as need for an exact system design for a systematic investigation, direct access to surrounding communities to gauge their perceptions, and an integrated software required to undertake a full-fledged risk assessment. These barriers were overcome using illustrative reference station designs, engaging with community-facing agencies through Port support, and pooling national laboratory capabilities available for risk assessments. The project identified relevant public safety risk metrics, compared various hydrogen carriers, engaged with community-facing agencies, and explored potential gaps in existing safety codes and standards. The primary impacts of this project include the development of risk assessment guidance for ports and utilities, informing them of the trade-offs in the choice of hydrogen carriers, and the ability to increase public capacity for dialog and engagement. This paves the way toward decarbonization of the Port activities, bringing about awareness around jobs in the market for risk assessments, and the need to ramp up community engagement long before any hydrogen system deployment is undertaken.

08 HYDROGEN↗
DOE OSTI2024

VAN-DAMME: GPU-accelerated and symmetry-assisted quantum optimal control of multi-qubit systems

We present an open-source software package, VAN-DAMME (Versatile Approaches to Numerically Design, Accelerate, and Manipulate Magnetic Excitations), for massively-parallelized quantum optimal control (QOC) calculations of multi-qubit systems. To enable large QOC calculations, the VAN-DAMME software package utilizes symmetry-based techniques with custom GPU-enhanced algorithms. This combined approach allows for the simultaneous computation of hundreds of matrix exponential propagators that efficiently leverage the intra-GPU parallelism found in high-performance GPUs. In addition, to maximize the computational efficiency of the VAN-DAMME code, we carried out several extensive tests on data layout, computational complexity, memory requirements, and performance. These extensive analyses allowed us to develop computationally efficient approaches for evaluating complex-valued matrix exponential propagators based on Padé approximants. To assess the computational performance of our GPU-accelerated VAN-DAMME code, we carried out QOC calculations of systems containing 10 - 15 qubits, which showed that our GPU implementation is 18.4× faster than the corresponding CPU implementation. Our GPU-accelerated enhancements allow efficient calculations of multi-qubit systems, which can be used for the efficient implementation of QOC applications across multiple domains.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2019

R and DE Robotic Sensor Intern Team: Remote 3D Mapping Overview

The Department of Energy (DOE) nuclear complex includes many hazardous and challenging environments that must be inspected remotely and routinely for structural integrity evaluation and long-term planning purposes. The H-Canyon Air Exhaust (HCAEX) Tunnel biannual inspection is an example of such an environment. Presently inspections are performed by teleoperated robotic crawlers provide low resolution video feedback to the structural engineers. Advanced sensors such as Lidar and panospheric hi-res cameras have the potential to provide additional valuable data in future inspections. UT-Austin Lidar work in support of an H-Canyon deployment began with the DOE EM sponsored H-Canyon Integrated Research Project. The R and DE Robotic Sensor intern team includes two of the UT-Austin Masters students working on Lidar development to facilitate the pulling of university developed solutions to the field. The goal of this project is to develop a suite of sensors that could be used to gather advanced data in future HCAEX tunnel inspections. Of particular interest is the ability to map the tunnel walls at an accuracy that will allow for the comparison of two scans performed at different time periods to estimate the concrete erosion rate. Additionally, the development of the 3D point cloud in the GPS denied environment is the first step in being able to collect and overlay other 3D data such as radiation mapping. Project specific goals include: - Design and develop 3D mapping sensor suite. Merge Lidar point cloud with high res camera images Incorporate Inertial Measurement Unit (IMU) to assist with positioning. - Develop sensor software and user GUI Sensor data collection and integration software User GUI for reach back system. - Test solution on crawler in tunnel test bed Use crawler base similar to H-Canyon planned crawler Design and build test bed simulating tunnel conditions. To facilitate rapid sensor integration, software coding and user GUI development, the Robot Operating System (ROS) was employed as the development platform. ROS is a set of open source software libraries and tools that assist with the building of robot applications. ROS includes drivers for many sensors, state-of-the-art algorithms, and powerful developer tools. Sensors selected for the summer intern solution included assessment as to whether they possessed ROS drivers and associated sample code. Additionally, the ROS 3D Robot Visualizer (rviz) tool has been employed.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2016

Towards a Consolidated Approach for the Assessment of Evaluation Models of Nuclear Power Reactors

The STARS project at the Paul Scherrer Institut (PSI) has adopted the TRACE thermal-hydraulic (T-H) code for best-estimate system transient simulations of the Swiss Light Water Reactors (LWRs). For analyses involving interactions between system and core, a coupling of TRACE with the SIMULATE-3K (S3K) LWR core simulator has also been developed. In this configuration, the TRACE code and associated nuclear power reactor simulation models play a central role to achieve a comprehensive safety analysis capability. Thus, efforts have now been undertaken to consolidate the validation strategy by implementing a more rigorous and structured assessment approach for TRACE applications involving either only system T-H evaluations or requiring interfaces to e.g. detailed core or fuel behavior models. The first part of this paper presents the preliminary concepts of this validation strategy. The principle is to systematically track the evolution of a given set of predicted physical Quantities of Interest (QoIs) over a multidimensional parametric space where each of the dimensions represent the evolution of specific analysis aspects, including e.g. code version, transient specific simulation methodology and model "nodalisation". If properly set up, such environment should provide code developers and code users with persistent (less affected by user effect) and quantified information (sensitivity of QoIs) on the applicability of a simulation scheme (codes, input models, methodology) for steady state and transient analysis of full LWR systems. Through this, for each given transient/accident, critical paths of the validation process can be identified that could then translate into defining reference schemes to be applied for downstream predictive simulations. In order to illustrate this approach, the second part of this paper presents a first application of this validation strategy to an inadvertent blowdown event that occurred in a Swiss BWR/6. The transient was initiated by the spurious actuation of the Automatic Depressurization System (ADS). The validation approach progresses through a number of dimensions here: First, the same BWR system simulation model is assessed for different versions of the TRACE code, up to the most recent one. The second dimension is the "nodalisation" dimension, where changes to the input model are assessed. The third dimension is the "methodology" dimension. In this case imposed power and an updated TRACE core model are investigated. For each step in each validation dimension, a common set of QoIs are investigated. For the steady-state results, these include fuel temperatures distributions. For the transient part of the present study, the evaluated QoIs include the system pressure evolution and water carry-over into the steam line.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2017

Modification of the SAS4A Safety Analysis Code for Integration with the ADAPT Discrete Dynamic Event Tree Framework.

It is difficult to assess the consequences of a transient in a sodium-cooled fast reactor (SFR) using traditional probabilistic risk assessment (PRA) methods, as numerous safety-related sys- tems have passive characteristics. Often there is significant dependence on the value of con- tinuous stochastic parameters rather than binary success/failure determinations. One form of dynamic PRA uses a system simulator to represent the progression of a transient, tracking events through time in a discrete dynamic event tree (DDET). In order to function in a DDET environment, a simulator must have characteristics that make it amenable to changing physical parameters midway through the analysis. The SAS4A SFR system analysis code did not have these characteristics as received. This report describes the code modifications made to allow dynamic operation as well as the linking to a Sandia DDET driver code. A test case is briefly described to demonstrate the utility of the changes.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗