Search NASA⌕ Search

SEARCH · Search NASA

Results for “code modernization”

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 163 records · Page 9

Light Output and Thermal Blast Analysis of Nuclear Fireballs

During the 1940s, 50s, and 60s, there were 210 atmospheric nuclear tests. Every test was extensively recorded with high quality motion picture films and photographs, providing a unique historic record, containing ~10,000 films, with a total of ~22 million frames of information. Much of what is known about the effects of nuclear weapons has been obtained from these tests. There are approximately 6500 scientific films that are in the process of being scanned by Lawrence Livermore National Laboratory and re-analyzed using modern computer techniques. The four major effects of a nuclear detonation are prompt radiation, shockwave, thermal blast, and fallout. The objective of this interim report is to describe recent efforts to extract benchmark data from the existing streak films using modern Figure-processing techniques. These new and, hopefully, greatly-improved data will provide our code developers with the benchmark data necessary to validate our modern computer codes that simulate nuclear weapon effects.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

The Variable Infiltration Capacity model version 5 (VIC-5): infrastructure improvements for new applications and reproducibility

The Variable Infiltration Capacity (VIC) model is a macroscale semi-distributed hydrologic model. VIC development began in the early 1990s and the model has since been used extensively for basin- to global-scale applications that include hydrologic dataset construction, trend analysis of hydrologic fluxes and states, data evaluation and assimilation, forecasting, coupled climate modeling, and climate change impact assessment. Ongoing operational applications of the VIC model include the University of Washington's drought monitoring and forecasting systems and NASA's Land Data Assimilation System. This paper documents the development of VIC version 5 (VIC-5), which includes a major reconfiguration of the legacy VIC source code to support a wider range of modern hydrologic modeling applications. The VIC source code has been moved to a public GitHub repository to encourage participation by the broader user and developer communities. The reconfiguration has separated the core physics of the model from the driver source code, whereby the latter is responsible for memory allocation, preprocessing and post-processing, and input–output (I–O). VIC-5 includes four drivers that use the same core physics modules, but which allow for different methods for accessing this core to enable different model applications. Finally, VIC-5 is distributed with robust test infrastructure, components of which routinely run during development using cloud-hosted continuous integration. The work described here provides an example to the model development community for extending the life of a legacy model that is being used extensively. The development and release of VIC-5 represents a significant step forward for the VIC user community in terms of support for existing and new model applications, reproducibility, and scientific robustness.

42 ENGINEERING↗

Calibration and Validation of NRC Fuel Performance Models

Pacific Northwest National Laboratory maintains steady-state and transient fuel performance codes for the U.S Nuclear Regulatory Commission (NRC). The NRC uses these codes in performing confirmatory calculations of licensee fuel performance codes and safety analysis methods during the review and approval of these codes and methods. The NRC also uses these codes to perform calculations to direct regulatory research and to assist in the development of new regulations and regulatory guides. Because of the safety-related uses for these codes, it is critically important that the codes are capable of accurately modeling modern fuel and operations and that the models and code predictions are well validated using data collected from conditions relevant to in-reactor conditions. FRAPCON and FRAPTRAN have an extensive database of separate effects data used to validate individual material property and fuel performance models and an integral effects database used to assess the overall code prediction within the full range of fuel operations. The key code predictions that are of regulatory interest as they relate to various specified acceptable design limits are; fuel and cladding temperature, rod internal pressure, cladding hoop strain, cladding corrosion and hydriding, and fission gas release. To assess the ability of the code to accurately predict these values, the integral assessment database consists of measurements of fuel centerline temperature, fission gas release, void volume, cladding corrosion, and cladding hoop strain. This paper will discuss the activities that are performed to ensure that FRAPCON and FRAPTRAN represent the state of the industry and will describe the separate effect and integral effect databases that are used to validate the models and are continually being expanded to ensure validation over the full range of fuel operations.

Geelhood, Kenneth J.↗

Validation of CTF Droplet Entrainment and Annular/Mist Closure Models using Riso Steam/Water Experiments

This document summarizes the work to validate the droplet entrainment and de-entrainment models, as well as two-phase closure models in the CTF (Coolant Boiling in Rod Arrays - Two Fluids) code by comparison with experimental data obtained at Risoe National Laboratory. CTF is the modernized version of the subchannel code, COBRA-TF, being jointly developed by Oak Ridge National Laboratory and North Carolina State University for applications in the Consortium for Advanced Simulation of Light Water Reactors (CASL). The Risoe data included a series of more than 250 steam/water experiments that were performed in both tube and annulus geometries over a range of various pressures and outlet qualities. Experimental conditions were set so that the majority of cases was in the annular/mist flow regime. Measurements included liquid film flow rate, droplet flow rate, film thickness, and two-phase pressure drop. CTF was used to model 180 of the tubular geometry cases, matching experimental geometry, outlet pressure, and outlet flow quality to experimental values. CTF results were compared to the experimental data at the outlet of the test section in terms of vapor flow fractions, and entrained liquid flow fractions, and pressure drop per unit length.

42 ENGINEERING↗

Update on the Code Intercomparison and Benchmark for Muon Fluence and Absorbed Dose Induced by an 18 GeV Electron Beam After Massive Iron Shielding

In 1974, Nelson, Kase and Svensson published an experimental investigation on muon shielding around SLAC high-energy electron accelerators [1]. They measured muon fluence and absorbed dose induced by 14 and 18 GeV electron beams hitting a copper/water beamdump and attenuated in a thick steel shielding. In their paper, they compared the results with the theoretical models available at that time. In order to compare their experimental results with present model calculations, we use the modern transport Monte Carlo codes MARS15, FLUKA2011 and GEANT4 to model the experimental setup and run simulations. The results are then compared between the codes, and with the SLAC data.

Fasso, A.↗

SAS4A/SASSYS-1 Validation with the FFTF LOFWOS Tests #10-12

The FFTF liquid sodium cooled test reactor was built to assist the development and testing of advanced fuels and material for fast reactors. The goals of the FFTF Passive Safety Testing program included confirming the safety margins of FFTF as a liquid metal reactor, providing data for computer code validation, and demonstrating the inherent and passive safety benefits of its specific design features. The program included thirteen unprotected loss of flow without scram tests. The results of modeling and simulation of LOFWOS Tests #10, #11, and #12 are presented in this report to support validation of Argonne’s SAS4A/SASSYS-1 fast reactor safety analysis code. These three tests did not use the primary loop pony motors and thus transitioned to natural circulation flow rates from varying initial power levels. This validation activity leverages modeling and simulation efforts that originated under an International Atomic Energy Agency Coordinated Research Project (CRP) for benchmark analysis of the FFTF LOFWOS Test #13. The SAS4A/SASSYS-1 model originally developed for the CRP has been extended to include LOFWOS Tests #10-12 and modernized to utilize new code features and modeling practices. Predicted results from SAS4A/SASSYS-1 were compared to measured test data and differences were investigated. Overall, it was concluded that, for LOFWOS Tests #10-12, there was good agreement between the flow, power, and temperature predictions and the measured data.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Validation of CTF Droplet Entrainment and Annular/Mist Closure Models using Riso Steam/Water Experiments

This document summarizes the work to validate the droplet entrainment and de-entrainment models, as well as two-phase closure models in the CTF code [1] by comparison with experimental data obtained at Risø National Laboratory [2]. CTF is the modernized version of the subchannel code, COBRA-TF [3], being jointly developed by Oak Ridge National Laboratory and North Carolina State University for applications in the Consortium for Advanced Simulation of Light Water Reactors (CASL). The Risø data included a series of more than 250 steam/water experiments that were performed in both tube and annulus geometries over a range of various pressures and outlet qualities. Experimental conditions wereset so that the majority of cases were in the annular/mist flow regime. Measurements included liquid film flow rate, droplet flow rate, film thickness, and two-phase pressure drop. CTF was used to model 180 of the tubular geometry cases, matching experimental geometry, outlet pressure, and outlet flow quality to experimental values. CTF results were compared to the experimental data at the outlet of the test section in terms of vapor flow fractions, and entrained liquid flow fractions, and pressure drop per unit length.

Wysocki, Aaron↗

Nuclear excitation functions for medical isotope production: Targeted radionuclide therapy via nat IR$(d, x)$ 193m Pt

193m Pt is an Auger emitting radionuclide which may have therapeutic potential, particularly when labeled to the chemotherapeutic drug cisplatin. One challenge to broader explorations of its clinical potential is the need for production routes with high specific activity. As part of a larger campaign to address gaps in reaction data for emerging medical radionuclides, this work seeks to characterize the nat Ir(d,x) reactions as a potential production pathway for 193m Pt. A stacked target irradiation, consisting of natural iridium, iron, nickel, and copper foils, was performed using a 33 MeV deuteron beam at the Lawrence Berkeley National Laboratory 88-Inch Cyclotron. This measurement, along with previous experimental data, suggests an energy window between 11 to 18 MeV to maximize the production and radiopurity of 193m Pt. This experiment has yielded cross sections for 43 channels of deuteron-induced reactions from threshold to 30 MeV, including the first experimental results of nat Ir(d,x) 188m1+g,190m1+g Ir (cumulative), nat Ni(d,x) 56,57,58 m,58g Co (independent), nat Cu(d,x) 61 Co (cumulative) and nat Fe(d,x) 53 Fe, 48 V (cumulative). The results were compared with literature data, the TENDL-2023 database, and default theoretical calculations from the TALYS-2.04, CoH-3.6.0, EMPIRE 3.2.3, and ALICE-2020 reaction modeling codes. Here, this work presents another example of the lack of predictive capabilities for this set of modern nuclear-reaction modeling codes, and highlights the unsatisfactory modeling of experimental cross sections. Experimental data are important to improve the codes in general, and new experimental results can be used to improve the models. Finally, this measurement has revealed the need for an updated evaluation of the nat Cu(d,x) 63 Zn deuteron monitor reaction.

193mPt↗

The Development of the Generalized Spallation Model [Thesis]

The generalized spallation model (GSM) is a general spallation physics event generator with a physics-rich history dating back to the early 1970s during the development of the CEM. GSM estimates the outcome of high-energy hadron-nucleus and nucleus-nucleus collisions, where collision energies often exceed hundreds of MeV. Such a collision may include a neutron incident on carbon, for example. The thesis presented here outlines the continued development and modernization of GSM, whereby it was migrated from a legacy Fortran code to utilize an advanced and modern Fortran architecture. Verification and validation of GSM also took place throughout the continued development of GSM. The results of this work have transformed GSM into a highly modern and powerful event generator, where numerous benefits have been fostered and precipitated, including but not limited to portability, modularization, flexibility, robustness, and scalability. GSM is a robust and well developed software library that clients and end-users may confidently use, within the scope outlined here, as a result of this work.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

SAS4A/SASSYS-1 Validation with EBR-II Tests Performed During the SHRT Testing Program

The EBR-II sodium-cooled fast reactor supported a wide range of fast reactor research activities. The final fifteen years of operations at EBR-II were used for experiments and tests to demonstrate the importance of passive safety in liquid metal reactors and the capability of such systems to provide strong passive responses during off-normal conditions. The Shutdown Heat Removal Test program at EBR-II provided test data supporting the validation of computer codes for design, licensing, and operation of LMRs, among other objectives. The program included nearly sixty tests, most of which were protected and unprotected loss of flow and loss of heat sink tests run at various initial powers and flow rates. This report documents the results of modeling and simulation of tests from the Shutdown Heat Removal Test program to support validation of Argonne’s SAS4A/SASSYS-1 fast reactor safety analysis code. SHRT-17 was a protected loss of flow test where a loss of electrical power to all sodium coolant pumps was simulated to demonstrate the effectiveness of natural circulation cooling characteristics. SHRT-45R was an unprotected loss of flow test where the control rod scram function of the plant protection system was disabled to demonstrate the effectiveness of EBR-II’s passive reactivity feedbacks. BOP-301 and BOP-302R were unprotected loss of heat sink tests that demonstrated how reactivity feedback effects driven by changes to the core inlet temperature can shut down the fission process. This validation activity leverages modeling and simulation efforts that originated under an International Atomic Energy Agency Coordinated Research Project for benchmark analysis of the EBR-II Shutdown Heat Removal Tests. The SAS4A/SASSYS-1 models originally developed for the CRP have been modernized to utilize new code features and modeling practices. Differences between the predicted and measured test data were explored and the causes of these differences were identified as being due to modeling approximations, uncertainty in transient component behavior, and instrumentation uncertainties. Overall, it was concluded that for SHRT-17, SHRT-45R, BOP-301, and BOP-302R there was good agreement between the flow, power, and temperature predictions and the measured data.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

SAS4A/SASSYS-1 Validation with EBR-II Tests Performed During the SHRT Testing Program (Rev. 2)

The EBR-II sodium-cooled fast reactor supported a wide range of fast reactor research activities. The final fifteen years of operations at EBR-II were used for experiments and tests to demonstrate the importance of passive safety in liquid metal reactors and the capability of such systems to provide strong passive responses during off-normal conditions. The Shutdown Heat Removal Test program at EBR-II provided test data supporting the validation of computer codes for design, licensing, and operation of LMRs, among other objectives. The program included nearly sixty tests, most of which were protected and unprotected loss of flow and loss of heat sink tests run at various initial powers and flow rates. This report documents the results of modeling and simulation of tests from the Shutdown Heat Removal Test program to support validation of Argonne’s SAS4A/SASSYS-1 fast reactor safety analysis code. SHRT-17 and SHRT-6 were protected loss of flow tests where a loss of electrical power to all sodium coolant pumps was simulated to demonstrate the effectiveness of natural circulation cooling characteristics. SHRT-45R, SHRT-39, and SHRT-43R were unprotected loss of flow tests where the control rod scram function of the plant protection system was disabled to demonstrate the effectiveness of EBR-II’s passive reactivity feedbacks. BOP-301 and BOP-302R were unprotected loss of heat sink tests that demonstrated how reactivity feedback effects driven by changes to the core inlet temperature can shut down the fission process. PICT-5/6 demonstrated how core power could be increased and decreased via changes in the intermediate sodium and feedwater flow rates. This validation activity leverages modeling and simulation efforts that originated under an International Atomic Energy Agency Coordinated Research Project for benchmark analysis of the EBR-II Shutdown Heat Removal Tests. The SAS4A/SASSYS-1 models originally developed for the CRP have been modernized to utilize new code features and modeling practices. Differences between the predicted and measured test data were explored and the causes of these differences were identified as being due to modeling approximations, uncertainty in transient component behavior, and instrumentation uncertainties. Overall, it was concluded that for all tests simulated there was good agreement between the flow, power, and temperature predictions and the measured data.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Coupled reactors analysis: New needs and advances using Monte Carlo methodology

Coupled reactors and the coupling features of large or heterogeneous core reactors can be investigated with the Avery theory that allows a physics understanding of the main features of these systems. However, the complex geometries that are often encountered in association with coupled reactors, require a detailed geometry description that can be easily provided by modern Monte Carlo (MC) codes. This implies a MC calculation of the coupling parameters defined by Avery and of the sensitivity coefficients that allow further detailed physics analysis. The results presented in this paper show that the MC code SERPENT has been successfully modifed to meet the required capabilities.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Data-driven modeling of coarse mesh turbulence for reactor transient analysis using convolutional recurrent neural networks

Advanced nuclear reactors often exhibit complex thermal-fluid phenomena during transients. To accurately capture such phenomena, a coarse-mesh three-dimensional (3-D) modeling capability is desired for modern nuclear-system code. In the coarse-mesh 3-D modeling of advanced-reactor transients that involve flow and heat transfer, accurately predicting the turbulent viscosity is a challenging task that requires an accurate and computationally efficient model to capture the unresolved fine-scale turbulence. In this work, we propose a data-driven coarse-mesh turbulence model based on local flow features for the transient analysis of thermal mixing and stratification in a sodium-cooled fast reactor. The model has a coarse mesh setup to ensure computational efficiency, while it is trained by fine-mesh computational fluid dynamics (CFD) data to ensure accuracy. A novel neural network architecture, combining a densely connected convolutional network and a long-short-term-memory network, is developed that can efficiently learn from the spatial temporal CFD transient simulation results. The neural network model was trained and optimized on a loss-of flow transient and demonstrated high accuracy in predicting the turbulent viscosity field during the whole transient. The trained model's generalization capability was also investigated on two other transients with different inlet conditions. The study demonstrates the potential of applying the proposed data-driven approach to support the coarse-mesh multi-dimensional modeling of advanced reactors.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Optimizing the hit finding algorithm for liquid argon TPC neutrino detectors using parallel architectures

Neutrinos are particles that interact rarely, so identifying them requires large detectors which produce lots of data. Processing this data with the computing power available is becoming even more difficult as the detectors increase in size to reach their physics goals. Liquid argon time projection chamber (LArTPC) neutrino experiments are expected to grow in the next decade to have 100 times more wires than in currently operating experiments, and modernization of LArTPC reconstruction code, including parallelization both at data- and instruction-level, will help to mitigate this challenge. The LArTPC hit finding algorithm is used across multiple experiments through a common software framework. In this paper we discuss a parallel implementation of this algorithm. Using a standalone setup we find speedup factors of two times from vectorization and 30–100 times from multi-threading on Intel architectures. The new version has been incorporated back into the framework so that it can be used by experiments. On a serial execution, the integrated version is about 10 times faster than the previous one and, once parallelization is enabled, more speedups comparable to the standalone program are achieved.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗

Towards Resiliency Evaluation of Vector Programs

The systems resilience research community has developed methods to manually insert additional source-program level assertions to trap errors, and also devised tools to conduct fault injection studies for scalar program codes. In this work, we contribute the first vector oriented LLVM-level fault injector VULFI to help study the effects of faults in vector architectures that are of growing importance, especially for vectorizing loops. Using VULFI, we conduct a resiliency study of nine real-world vector benchmarks using Intel’s AVX and SSE extensions as the target vector instruction sets, and offer the first reported understanding of how faults affect vector instruction sets. We take this work further toward automating the insertion of resilience assertions during compilation. This is based on our observation that during intermediate (e.g., LLVM-level) code generation to handle full and partial vectorization, modern compilers exploit (and explicate in their code-documentation) critical invariants. These invariants are turned into error-checking code. We confirm the efficacy of these automatically inserted low-overhead error detectors for vectorized for-loops.

fault injection↗

A Statistical Interpolation Code for Ocean Analysis and Forecasting

Abstract We present a data assimilation package for use with ocean circulation models in analysis, forecasting, and system evaluation applications. The basic functionality of the package is centered on a multivariate linear statistical estimation for a given predicted/background ocean state, observations, and error statistics. Novel features of the package include support for multiple covariance models, and the solution of the least squares normal equations either using the covariance matrix or its inverse—the information matrix. The main focus of this paper, however, is on the solution of the analysis equations using the information matrix, which offers several advantages for solving large problems efficiently. Details of the parameterization of the inverse covariance using Markov random fields are provided and its relationship to finite-difference discretizations of diffusion equations are pointed out. The package can assimilate a variety of observation types from both remote sensing and in situ platforms. The performance of the data assimilation methodology implemented in the package is demonstrated with a yearlong global ocean hindcast with a 1/4° ocean model. The code is implemented in modern Fortran, supports distributed memory, shared memory, multicore architectures, and uses climate and forecasts compliant Network Common Data Form for input/output. The package is freely available with an open source license from www.tendral.com/tsis/ .

Srinivasan, Ashwanth↗

Evolution of a minimal parallel programming model

Here, we take a historical approach to our presentation of self-scheduled task parallelism, a programming model with its origins in early irregular and nondeterministic computations encountered in automated theorem proving and logic programming. We show how an extremely simple task model has evolved into a system, asynchronous dynamic load balancing (ADLB), and a scalable implementation capable of supporting sophisticated applications on today’s (and tomorrow’s) largest supercomputers; and we illustrate the use of ADLB with a Green’s function Monte Carlo application, a modern, mature nuclear physics code in production use. Our lesson is that by surrendering a certain amount of generality and thus applicability, a minimal programming model (in terms of its basic concepts and the size of its application programmer interface) can achieve extreme scalability without introducing complexity.

97 MATHEMATICS AND COMPUTING↗

Update of the Nuclear Criticality Slide Rule Calculations: Plutonium systems – Delayed Fission Gamma

IRSN (France), LLNL (USA) and ORNL (USA) began a long-term collaboration effort in 2015 to update the nuclear criticality Slide Rule for the emergency response to a nuclear criticality accident. The Slide Rule permits the estimation of neutron and gamma dose rates and integrated doses based upon estimated fission yields, as a function of distance from the fission source, and time after criticality accidents for different critical systems. This paper presents results from the fourth phase of the current update of the Slide Rule project, in which delayed fission-product gamma (DFG) dose rates of unreflected plutonium critical systems were compared by several modern 3D radiation transport codes (MCNP, COG, SCALE), using updated flux-to-dose conversion factors. Dose rates are calculated for fissile material at five moderation ratios (H/Pu), at 1 m above the ground as a function of distance (between 30 cm and 1,200 m) from the external surface of the source to the center of the detector, and for periods between 1 s and 1,000 min after the critical instantaneous event. Further efforts have been devoted to the determination of the delayed gamma source, by comparing the time-dependent energy spectra obtained from several methods. Overall, DFG dose rates calculated by each participant led to consistent results. Extra effort is under way to identify the cause of the remaining differences, by comparing precisely the gamma source, and particularly nuclides inventories.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗