Search NASA⌕ Search

SEARCH · Search NASA

Results for “C Code”

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 415 records · Page 23

Multiphysics FISPACT-II and TENDL-2019 simulation: neutron-induced damage metrics

Nuclear interactions can be the source of atomic displacement, irradiation-induced defects and transmutation in structural materials. Such quantities are derived from, or can be correlated to, nuclear kinematic simulations of primary atomic energy distributions spectra and the quantification of the numbers of secondary defects produced per primary as a function of the available recoils, residual and emitted, energies. Recoil kinematics of neutral, residual, gas, proton, alpha particle emissions are now more rigorously treated based on recent, complete and enhanced nuclear data parsed in state-of-the-art processing tools. Defect production metrics are the starting point in the complex problem of correlating and simulating the behaviour of materials under irradiation, as experimental information is rare or scattered. Detailed, segregated primary knock-on-atom metrics are now becoming available as the starting point of further simulation processes of isolated and clustered defects in material lattices. This allows more materials, lattices, neutron incident energy ranges, and irradiation conditions to be explored with sufficient data to adequately cover both standard and novel applications and materials: the broader reactor applications landscape. The damage metrics of of materials are systematically explored under typical but different reactor's type environment. The inventory code FISPACT-II combined with the enhanced nuclear data forms of the TENDL-2019 libraries allow one to not only calculate dpa from mostly scattering events but to also properly predict gas production, nuclear heating and transmutation under the same conditions.

22 GENERAL STUDIES OF NUCLEAR REACTORS↗

Bridging the gap: Deploying AI-based Models in Real-Time Fusion Plasma Control Systems

Achieving reliable real-time control in fusion plasma experiments requires strict timing guarantees across entire control algorithms. In earlier work by Abbate et al. (2023), we demonstrated the feasibility of neural-network-based control algorithms on the DIII-D tokamak using the internally developed open-source Keras2C library for model conversion into C (Conlin et al. (2021)). However, the initial implementations relied on data buffering and branching logic outside the neural network code, causing variability in execution times. Subsequent deployments on DIII-D and KSTAR—including the RTCAKENN algorithm for kinetic profile reconstruction—proved that minimizing branching and buffering throughout the pipeline yields consistent millisecond-level cycle times under real experimental conditions (Shousha et al. (2023)). However, keeping pace with rapidly evolving AI frameworks (e.g. PyTorch) is challenging. Finally, we, therefore, propose a community-driven open-source effort to expand the tool, enabling real-time deployment across diverse systems that require strictly bounded execution times.

AI-based models↗

Fast ion relaxation in ITER mediated by Alfvén instabilities

Abstract We address the critical issue for future burning plasmas of whether high-energy fusion products or auxiliary heating-beam ions will be confined for a sufficiently long time to compensate for thermal plasma energy losses. This issue can be mitigated by one of the most deleterious collective phenomena—the instability of low, sub-cyclotron frequency Alfvén eigenmodes (AEs), such as toroidicity-induced AEs and reversed-shear AEs in the ITER steady-state scenario. Using a revised quasi-linear (QL) theory applied to energetic particle (EP) relaxation in the presence of AEs, we find that the AE instabilities can affect both neutral beam ions and alpha particles, although the resulting fast ion transport is expected to be modest if classical particle slowing down is assumed. On the other hand, the QL theory predicts that the AE amplitudes will be enhanced by the background microturbulence, although this topic remains outside our scope due to the significant numerical effort required to evaluate these effects. We report our results for EP relaxation dynamics obtained utilizing several tools: (i) a comprehensive linear stability study of the sub-cyclotron Alfvénic spectrum as computed by ideal magnetohydrodynamic NOVA simulations for the AE eigenproblem, (ii) drift kinetic NOVA-C calculations for wave–particle interaction and AE growth/damping rates, and (iii) predictive QL modeling coupled with the global transport code TRANSP to assess the EP relaxation on the equilibrium timescale.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

The Completion of Surrogate Testing for Low-Temperature TEX and a Look Towards the Future

To address the mounting need for below room temperature nuclear data validation, the Low-Temperature Thermal Epithermal eXperiments (LT-TEX) have been designed. Validation of low-temperature neutron cross sections is necessary to verify any operation at temperatures below room temperature which is typically observed in environments far from the equator. For example, a fissile material transportation truck may routinely observe ambient temperatures down to -40°C, which is the lower temperature bound of the normal conditions of transportation defined in the United States Title 10 Code of Federal Regulations §71.71c2. Additionally, sub-room temperature benchmarks can validate newly produced cross sections, that include novel thermal scattering laws, from North Carolina State University.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Prediction of tissue optical properties using the Monte Carlo modeling of photon transport in turbid media and integrating spheres

Monte Carlo methods are an established technique for simulating light transport in biological tissue. Integrating spheres make experimental measurements of the reflectance and transmittance of a sample straightforward and inexpensive. This work presents an extension to existing Monte Carlo photon transport methods to simulate integrating sphere experiments. Crosstalk between spheres in dual-sphere experiments is accounted for in the method. Analytical models, previous works on Monte Carlo photon transport, and experimental measurements of a synthetic tissue phantom validate this method. We present two approaches for using this method to back-calculate the optical properties of samples. Experimental and simulation uncertainties are propagated through both methods. Both back-calculation methods find the optical properties of a sample accurately and precisely. Our model is implemented in standard Python 3 and CUDA C++ [J. Nickolls, I. Buck, M. Garland, and K. Skadron, ACM Queue 6 , 40 ( 2008 ) ] and is publicly available in Code 1.

Cook, Patrick D. (ORCID:0000000279345428)↗

MSTAR2019 Code Description and User's Manual

This report describes an ideal cascade model of uranium enrichment that includes side feed and side product streams and flexible input options. It is based on the original MSTAR model developed by Ed Von Halle and implemented in a Visual Basic code. The current version allows the user to specify integer numbers of stages or the stage numbers of external flows instead of assays in those flows. The computational engine is written in FORTRAN-90 and is invoked by a user-friendly GUI written in C++. This version of MSTAR has been demonstrated to operate on Linux, Mac, and Windows platforms, and has undergone significant testing and quality analysis. A number of examples are presented to illustrate the operation of the code and guide the user. The code is extremely fast, and results are returned immediately. The input and output have been specially configured for analysis by the environmental sampling team of the International Atomic Energy Agency (IAEA). A number of additional output features are included to assist the user in visualizing computational results and downloading data to files for use in other analysis software.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗

TIGER, A thermodynamics Equilibrium Tool for Explosives Update: Improved Solver

The thermodynamic equilibrium code known as TIGER has been in use since the 1970s and is designed to calculate the performance of energetic materials during detonation at high temperatures and pressures. The original TIGER code utilized a limited database consisting of 12 gaseous and 3 condensed constituents, which were composed of the elements carbon (C), hydrogen (H), nitrogen (N), oxygen (O), and aluminum (Al). In contrast, the more modern JCZS3 database features a significantly larger dataset, including 756 gaseous species, 189 positive and negative ions, and 496 condensed constituents derived from 62 different elements. While this expanded product species database allows for the exploration of a broader range of problems relevant to our laboratory, it also introduces challenges related to convergence, particularly when both liquid and solid species coexist under a vapor dome. In this work, we present an improved TIGER solution technique aimed at accurately determining the equilibrium state for systems that can produce a diverse array of products, including both solid and liquid condensed species. This report presents the status of the TIGER solver as of the end of FY25. Ongoing efforts are focused on enhancing the solver, and the report outlines several planned improvements. We anticipate that further developments will require additional documentation in the future.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Corrosion of SS 316, Inconel 617, Inconel 625 in molten chloride salt with emphasis on the speciation of Cr, Ni, and Co in molten salt matrix by XAFS

Dr. Hartmann’ report summarizes the corrosion behavior of materials used in the construction and operation of molten salt reactors which is is of great interest to the molten salt reactor (MSR) community. In the present study we intent to decipher the very elemental phenomena on the corrosion of austenitic chromium-containing stainless steel SS316, as well as Inconel 625 and ASME-coded high temperature Alloy (Inconel) 617. We initiated and are in continuation of conducting ampoule type corrosion testing at 550 °C and 650 °C to determine chromium diffusivity and the related activation energy as well as the speciation of the corrosion products.

36 MATERIALS SCIENCE↗

A Look Towards the Execution of the Low-Temperature TEX Experimental Campaign

To address the mounting need for below roomtemperature nuclear data validation, the Low-Temperature Thermal Epithermal eXperiments (LT-TEX) have been designed. Validation of low-temperature neutron cross sections is necessary to verify any operation at temperatures below room temperature which is typically observed in environments far from the equator. For example, a fissile material transportation truck may routinely observe ambient temperatures down to -40°C, which is the lower temperature bound of the normal conditions of transportation defined in the United States Title 10 Code of Federal Regulations §71.71c2. Additionally, sub-room temperature benchmarks can validate newly produced cross sections, that include novel thermal scattering laws, from North Carolina State University.

LT-TEX↗

Using a Large Language Model as a Building Block to Generate Usable Validation and Verification Suite for OpenMP

In the HPC area, both hardware and software move quickly. Often new hardware is developed and deployed, the corresponding software stack, including compilers and other tools, are under active development while leading edge software developers are working to port and tune their applications, all at the same time. While the software ecosystem is in flux, one of the key challenges for users is obtaining insight into the state of implementation of key features in the programming languages and models their applications are using – whether they have been implemented, and whether the implementation conforms to the specification, especially for newly implemented features (less tested by widespread use). OpenMP is one of the most prominent shared memory programming models used for on-node programming in HPC. With the shift towards accelerators (such as GPUs and FPGAs) and heterogeneous programming OpenMP features are getting more complex. It is natural to ask whether generative AI approaches, and large language models (LLMs) in particular, can help in producing validation and verification test suites to allow users better and faster insights into the availability and correctness of OpenMP features of interest. In this work, we explore the use of ChatGPT-4 to generate a suite of tests for OpenMP features. We have chosen a set of directives and clauses, a total of 78 combinations, which first appeared in OpenMP 3.0 (released in May 2008) but are also relevant for accelerators. We prompted ChatGPT to generate tests in the C and Fortran languages, for both host (CPU) and device (accelerator). On the Summit super-computer using the GNU implementation, we found that, of the 78 generated tests 67 C tests and 43 Fortran tests compiled successfully and fewer than those executed to completion. On further analysis we show that not all generated tests are valid. We document the process, results, and provide detailed analysis regarding the quality of tests generated. With the aim of providing input to a production quality validation and verification suite, we manually implement the corrections required to make the tests valid according to the current OpenMP specification. We quantify this effort as small, medium, or large, and record the lines of code changed to correct the invalid tests. With the corrected tests we validate recent implementations from HPE, AMD, and GNU on the Frontier supercomputer. Our experiment and subsequent analysis show that although LLMs are capable of producing HPC specific codes, they are limited by their understanding of the deeper semantics and restrictions of programming models such as OpenMP. Unsurprisingly more commonly used features have better support, while some OpenMP 3.0 directives such as sections and tasking are not universally supported on accelerators. We demonstrate that successful compilation and execution to completion are inadequate metrics for evaluating generated code and that, at this time, commodity LLMs require expert intervention for code verification. This points to gaps in the training data that is currently available for HPC. We demonstrate that with "small" effort 37% of generated invalid C tests and 63% of generated invalid Fortran tests could be corrected. This improves productivity of test generation as we circumvent writing from scratch and the common programming errors associated with it.

Pophale, Swaroop [ORNL] (ORCID:0000000185446367)↗

Assessment of RELAP5-3D Code with Molten Salt Heat Transfer Experiments

The accurate thermal-hydraulic assessment of thermal storage systems is crucial for addressing the integrity and performance of the storage system design. This is particularly important for thermal storage systems using molten salts as a heat transport and storage medium, as the chemical corrosion and erosion by the high-temperature molten salts can cause serious damage to the structural components of the storage system. To understand these effects for TerraPower’s Natrium® Demonstration Reactor design, the system thermal-hydraulic analysis code, RELAP5-3D, is utilized to analyze the thermal storage system. This study evaluates two heat transfer correlations implemented in the RELAP5-3D code—the Dittus-Boelter and Gnielinski correlations—using selected benchmark cases from two molten salt heat transfer experiments conducted at Xi’an Jiao Tong University and the German Aerospace Center. The Nusselt numbers calculated using the Gnielinski correlation agree with the experimental data within ±5% of the relative deviations at 300°C and 400°C. With the increasing Reynolds numbers, the Dittus-Boelter correlation underestimates the Nusselt numbers by up to 22% compared to the Gnielinski correlation. Comparison of the RELAP5-3D calculation results with experimental data at a high temperature of 550°C revealed that the thermo-chemical characteristics of solar salt at high temperatures surpassing the decomposition temperatures of nitrate salts reduce the accuracy of heat transfer model of the RELAP5-3D code. Specifically, it has been observed that the temperatures of the wall surface and liquid near the wall can surpass the decomposition temperatures, even though the bulk temperature of molten salts remains significantly lower than these decomposition temperatures. In summary, this study recommends the use of the Gnielinski correlation for the heat transfer analysis of molten salt, rather than the Dittus-Boelter correlation.

21 - SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLAN↗

Depletion capabilities in the OpenMC Monte Carlo particle transport code

A depletion solver has been implemented in OpenMC and is described herein. The depletion solver is implemented in Python and interfaces with OpenMC’s transport solver through a C++ application programming interface, which enables an in-memory transport-depletion coupling. Multiple integration methods for advancing in time have been implemented and exhibit tradeoffs in cost, accuracy, and memory use. For all time integration methods, evaluation of the matrix exponential is performed by using the incomplete partial fraction form of the Chebyshev rational approximation method. Simulations of a pressurized water reactor (PWR) pincell and a sodium-cooled fast reactor (SFR) assembly were carried out with OpenMC and Serpent. For both problems, the use of a high-fidelity depletion chain results in predictions of k eff that agree within 20–30 pcm between OpenMC and Serpent. Predicted actinide concentrations were found to agree to a fraction of a percent, and most fission product concentrations were found to agree within 1%. Here, the few cases where larger differences were observed can be attributed either to differences in how the energy dependence of fission product yields is handled or deficiencies in the nuclear data used. OpenMC simulations of the PWR and SFR problems using a simplified 228-nuclide depletion chain demonstrate that it achieves accuracy close to that of the full, high-fidelity depletion chain with respect to the studied responses.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

In silico evolution of oscillatory gene dynamics formerly Code sharing for a manuscript to be submitted to the journal iScience

The disclosure is about the numerical simulations and visualization implemented in general-purpose language Python and C++. The simulations were in-silico evolution of hypothetical gene circuits that potentially admit oscillatory dynamics. The aim of the study, which was invited to be submitted to a special issue "Timekeeping, noise and development across kingdoms" in the journal iScience, is to demonstrate that in contrast to common belief, noise in gene expression dynamics could significantly facilitate the evolutionary process for searching oscillatory dynamics, which are essential for developing circadian rhythm in living organisms.

Lin, Yen Ting↗

Deployment of the HFIRCON transport and depletion tool for plutonium-238 production studies

Irradiation of {sup 237}Np-bearing targets in Oak Ridge National Laboratory's (ORNL) High Flux Isotope Reactor (HFIR) results in the efficient production of {sup 238}Pu, which, in the form of heat source PuO{sub 2}, is used as a reliable power source for deep-space and planetary NASA missions. A technology demonstration subproject was initiated at ORNL in 2011 to develop and implement the technology required to establish a {sup 238}Pu supply chain. A systematic progression of NpO{sub 2}/Al cermet (20 vol.% NpO{sub 2}) activities to date has successfully demonstrated target fabrication, irradiation, and chemical recovery processes. Recent program tasks have included the development of the HFIRCON transport and depletion tool for efficient reactor physics analyses and the evaluation of increased NpO{sub 2} loadings (i.e., beyond 20 vol.%) and NpN-based targets. This paper documents the deployment of the HFIRCON code to assess various Np concentrations in NpO{sub 2}- and NpN-based targets in HFIR's inner small vertical experiment facilities. Results indicate that {sup 238}Pu production and quality can be enhanced with increased Np loadings; however, target conversion rates are reduced. The results recorded in this paper, thermal and material balance evaluations, and testing requirement planning will be used to determine whether increased NpO{sub 2} loadings or NpN-based targets will be further considered. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗

Power distribution estimation method for SMR using ex-core detectors: experimental demonstration by plural control rod patterns at KUCA

The power distribution estimation method based on the ex-core detectors, or PHOEBE, was demonstrated at Kyoto University Critical Assembly, KUCA. Generally, core monitoring systems use in-core neutron detectors. Since inside the core is a harsh environment, the maintainability and reliability of the detectors are deteriorated. On the other hand, the environment outside the core is milder: core monitoring by the ex-core detectors improves the maintainability and reliability especially for small modular reactors and micro-reactors. However, neutron information from the inner region of the core is lost at the ex-core detectors. To recover the information, the authors proposed to utilize the power correlation between the fuel regions. PHOEBE concept was demonstrated at KUCA with distorted power distributions simulated by control rod patterns. The relative power distribution estimated by PHOEBE agreed with that calculated by Monte Carlo simulation code MVP. PHOEBE approximately reproduced the trend of the distorted power distributions calculated by MVP code. By contrast, the case of without power correlation between fuel regions produced significant different power distributions. Therefore, the advantage of considering the power correlation between fuel regions was also demonstrated. (authors)

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗