Search NASA⌕ Search

SEARCH · Search NASA

Results for “System 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 559 records · Page 31
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Uncertainty Quantification in High-Low Dynamic System Coupling using RAVEN and TRANSFORM

This work demonstrates new functionality and applications stemming from the development of high-fidelity to low-fidelity (high-low) coupling for system simulations and to further explore the capabilities of the Risk Analysis Virtual Environment (RAVEN) in the performance of uncertainty quantification in this kind of high-low coupled system models. The work builds from previous work on high-low coupling that utilized COBRA-TF (CTF), the high-fidelity subchannel analysis code, with a low fidelity model built in ORNL’s TRANSFORM, the system analysis code, utilizing the Functional Mock-Up Interface (FMI). Steady-state and transient analysis examples using the high/low coupled models generated from CTF and TRANSFORM/FMI are investigated. The workflows for both steady-state and transient coupled simulations are described. A steady-state parameter sweep and uncertainty analysis of the primary flow rates and reactor power are demonstrated. Likewise, a transient pump trip and power ramp sensitivity studies are also demonstrated. This work elucidates some of the potential benefits and future needs of using RAVEN for high/low system coupling analysis of energy systems. It also shows some of the difficulties that can be encountered in coupling system simulations.

Williams, Wesley↗
DOE OSTI2025

Analysis of US Industrial Assessment Centers (IACs) implementation

Industrial energy assessments are a fundamental action toward developing a decarbonization strategy at any level. They provide an understanding of where energy efficiency opportunities exist and help to make informed business decisions about the costs and benefits of implementing sustainable policies and practices. This paper examines the effectiveness of the Department of Energy's Industrial Assessment Centers Program at providing useful energy-efficiency audits as well as the barriers faced by the program and plans for future growth and improvement. This paper presents an analysis of the program between 1981 and 2022, covering 20,290 industrial assessments and 151,198 recommendations, with $2.6 billion of recommended savings. The analysis includes a breakdown of the IAC recommendations and implemented projects based on both the industrial subsectors according to the Standard Industrial Classification code, energy type, and systems evaluated. The results include a 47 % implementation rate, a gap analysis to understand the missed opportunities, and a discussion about reasons for recommendation rejection. A comparison of the IAC Program with other country-level energy auditing or assessment programs was conducted, and suggestions for improving implementation rates were mentioned.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Probing the antiknock effect of anisole through an ignition, speciation and modeling study of its blends with isooctane

In order to unravel the reaction pathways relevant to anisole co-oxidation within a fuel blend, a detailed study of isooctane/anisole blends was performed with the ULille RCM. Ignition delays as well as mole fraction profiles were measured during a two-stage ignition delay using sampling and GC techniques. These results are used to validate a kinetic model developed from ab initio calculations for the most relevant rate constants which included H-atom abstraction reactions from anisole, and reactions on the potential energy surfaces of methoxyphenyl + O 2 and anisyl + O 2 . Pressure dependent rate constants were computed for the methoxyphenyl + O 2 and anisyl + O 2 reactive systems using master equation code analysis. The new kinetic model shows good agreement with the experimental data. Dual brute-force sensitivity analysis was performed, on both first- and second-stages of ignition, allowing the identification of the most important reactions in the prediction of both ignition delays. It was observed that while pure anisole does not show NTC behavior, a 60/40 isooctane/anisole blend displays such behavior, as well as two-stage ignition. This suggests anisole addition may not be as beneficial to knock resistance as expected from its high octane number. The kinetic modeling results demonstrate the importance of H-abstraction reactions both from the methoxy group and from the aryl ring in ortho-position and the addition of the resultant radicals to O2, mostly leading to the formation of polar or non-aromatic products.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2025

Magnetohydrodynamic simulations preliminarily predict the habitability and radio emission of TRAPPIST-1e

TRAPPIST-1e, an Earth-sized exoplanet in the habitable zone of the nearby M dwarf TRAPPIST-1, may experience magnetospheric responses that vary with stellar space weather, which could potentially influence both its habitability and radio emissions. Our objective is to investigate how different Earth-like magnetospheric configurations of TRAPPIST-1e – specifically variations in dipolar magnetic field strength and axial tilt – respond to diverse stellar space weather conditions, including events analogous to coronal mass ejections (CMEs), and to assess their implications for potential habitability and expected radio emissions. We conducted 3D magnetohydrodynamic simulations of the TRAPPIST-1e system using the PLUTO code in spherical coordinates. The planetary magnetic field was modelled as dipolar, with equatorial strengths from Earth-like to several times stronger. The dipole axis spans a representative range of axial tilts. We investigate four stellar wind environments, from sub-Alfvénic flow to CME-like disturbances. Planetary shielding was quantified based on the magnetopause standoff distance, and radio powers were estimated via empirical scaling laws. Our simulations show that both shielding and radio power depend strongly on the magnetic configuration. Stronger fields increase protection, while larger tilts reduce it. Radio power increases with both field strength and tilt across all wind regimes. An Earth-like magnetic field can provide effective shielding even under intense CMEs, whereas high tilts require stronger fields. Predicted radio powers reach ~10 20 erg s −1 during CMEs, making bursts from close-in, magnetised planets more detectable. However, for TRAPPIST-1e, the maximum cyclotron frequency lies below the Earth’s ionospheric cutoff (~10 MHz), making ground-based detection currently infeasible.

Sun: magnetic fields↗
DOE OSTI2025

RELAP5-3D Capabilities for High-Temperature Gas-Cooled Reactor Analysis

RELAP5-3D is the premier systems thermal-hydraulic code in the world. The code has been used significantly for high-temperature gas-cooled reactor (HTGR) analyses over the last 10 years, but there has not been a comprehensive summary of RELAP5-3D capabilities and methods for HTGR modeling. In this paper, we provide a guide for how novice users can build models for prismatic HTGRs. This guide includes best practices for developing heat structures and setting up radiation and conduction heat transfer. We also provide a summary of relevant validation studies. These studies are divided into three categories: separate effects tests, integral effects tests, and nuclear reactor studies. We present a comprehensive review of RELAP5-3D HTGR applications for reactor modeling over the last 10 years. These applications include work for the Modular High-Temperature Gas-Cooled Reactor benchmark, design activities for a prismatic HTGR as part of the U.S. Department of Energy’s Advanced Demonstration and Test Reactor Options Study, and studies performed to design gas-cooled microreactors. Following the comprehensive literature review, we present a gap analysis highlighting capabilities of RELAP5-3D that are available and missing for a complete analysis of HTGR thermal hydraulics. We also discuss the limitations of existing code capabilities and potential improvements that could improve analyses. Current research efforts have yet to fully validate RELAP5-3D for a prismatic HTGR analysis; however, some of this is attributable to confounding factors in the available experimental data as opposed to code capabilities. In conclusion, this paper presents the most comprehensive review to date on RELAP5-3D capabilities, limitations, and applications for HTGR analysis.

22 - GENERAL STUDIES OF NUCLEAR REACTORS↗
DOE OSTI2021

Anomaly detection in the Zwicky Transient Facility DR3

We present results from applying the SNAD anomaly detection pipeline to the third public data release of the Zwicky Transient Facility (ZTF DR3). The pipeline is composed of three stages: feature extraction, search of outliers with machine learning algorithms, and anomaly identification with followup by human experts. Our analysis concentrates in three ZTF fields, comprising more than 2.25 million objects. A set of four automatic learning algorithms was used to identify 277 outliers, which were subsequently scrutinized by an expert. From these, 188 (68 per cent) were found to be bogus light curves – including effects from the image subtraction pipeline as well as overlapping between a star and a known asteroid, 66 (24 per cent) were previously reported sources whereas 23 (8 per cent) correspond to non-catalogued objects, with the two latter cases of potential scientific interest (e.g. one spectroscopically confirmed RS Canum Venaticorum star, four supernovae candidates, one red dwarf flare). Moreover, using results from the expert analysis, we were able to identify a simple bi-dimensional relation that can be used to aid filtering potentially bogus light curves in future studies. We provide a complete list of objects with potential scientific application so they can be further scrutinised by the community. These results confirm the importance of combining automatic machine learning algorithms with domain knowledge in the construction of recommendation systems for astronomy. Our code is publicly available.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2020

Averaged cluster approach to including chemical short-range order in KKR-CPA

The single-site Korringa-Kohn-Rostoker coherent potential approximation (KKR-CPA) ignores short-range ordering present in disordered metallic systems. Here in this paper, we establish a technique to fix this shortcoming by embedding an averaged cluster that displays chemical short-range order (SRO). The degree of SRO can be tuned by externally defined order parameters. This averaged cluster can be embedded in the single-site CPA medium or a self-consistently obtained effective medium that contains SRO information. The validity of this method is demonstrated by applying it to two alloy systems—the CuZn body-centered cubic (BCC) solid solution, and AlCrTiV, a four-element BCC high entropy alloy. A comparison between the non-self-consistent and self-consistent modes is also provided for the two above-mentioned systems. We make the code available on the internet. Planned extensions to this paper are discussed.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2022

Differentiable Preisach Modeling for Characterization and Optimization of Particle Accelerator Systems with Hysteresis

Future improvements in particle accelerator performance are predicated on increasingly accurate online modeling of accelerators. Hysteresis effects in magnetic, mechanical, and material components of accelerators are often neglected in online accelerator models used to inform control algorithms, even though reproducibility errors from systems exhibiting hysteresis are not negligible in high precision accelerators. Here, we combine the classical Preisach model of hysteresis with machine learning techniques to efficiently create nonparametric, high-fidelity models of arbitrary systems exhibiting hysteresis. We experimentally demonstrate how these methods can be used in situ, where a hysteresis model of an accelerator magnet is combined with a Bayesian statistical model of the beam response, allowing characterization of magnetic hysteresis solely from beam-based measurements. Finally, we explore how using these joint hysteresis-Bayesian statistical models allows us to overcome optimization performance limitations that arise when hysteresis effects are ignored.

43 PARTICLE ACCELERATORS↗
DOE OSTI2020

Quantum Dynamics at Scale: Ultrafast Control of Emergent Functional Materials

Confluence of extreme-scale quantum dynamics simulations (i.e. quantum@scale) and cutting-edge x-ray free-electron laser experiments are revolutionizing materials science. An archetypal example is the exciting concept of using picosecond light pulses to control emergent material properties on demand in atomically-thin layered materials. This paper describes efforts to scale our quantum molecular dynamics engine toward the United States' first exaflop/s computer, under an Aurora Early Science Program project named "Metascalable layered material genome". Key algorithmic and computing techniques incorporated are: (1) globally-scalable and locally-fast solvers within a linear-scaling divide-conquer-recombine algorithmic framework; (2) algebraic `BLASification' of computational kernels; and (3) data alignment and loop restructuring, along with register and cache blocking, for enhanced vectorization and efficient memory access. The resulting weak-scaling parallel efficiency was 0.93 on 131,072 Intel Xeon Phi cores for a 56.6 million atom (or 169 million valence-electron) system, whereas the various code transformations achieved 5-fold speedup. The optimized simulation engine allowed us for the first time to establish a significant effect of substrate on the dynamics of layered material upon electronic excitation.

Algebraic BLASification↗
DOE OSTI2021

PSSE Network Reduction

This code can reduce power system network models in PSSE. The user selects single-port or two-port subnetworks for reduction, and the software reduces the subnetwork as described in the paper "Power System Network Reduction for Power Hardware-in-the-Loop Simulation" (submitted).

Wang, Bin↗
DOE OSTI2019

Design and Evaluation of a Pixelated PSD-capable Scintillator Detector with SiPM Readout

We present the detector response comparison between a 10x10 pixellated array of scintillator read out with Anger logic using four 2" Hamamatsu R7724-100 super bialkali photomultiplier tubes (PMT) and a custom Silicon photomultiplier (SiPM) board consisting of 100 C-series 6x6 mm SiPMs from SensL. An array of these pixellated detectors are currently used in the Neutron Coded Aperture (NCA) imaging system. The energy, timing and pulse shape discrimination response using both readout schemes are presented, along with an analysis of multiple scatter events occurring within the block. An evaluation of the impact of photodetector readout on the overall detection efficiency and imaging accuracy is presented.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2020

Requirements Description of the REBUS Software for the Versatile Test Reactor

This report presents the modeling and simulation capabilities of Argonne National Laboratory’s fuel cycle analysis code, the REactor BUrnup System (REBUS), that will be relied on for the VTR project. These capabilities will then be used to establish the set of REBUS verification tasks necessary to verify REBUS for usage on the VTR project. A similar path was followed for the DIF3D software where its’ requirements and verification tasks were established. REBUS has been maintained by Argonne since the early 1960s to support its reactor design mission. That software transitioned from the original REBUS to REBUS-2 in the mid 1970s and to REBUS-3 in the mid 1980s. Since then REBUS has gone through many revisions to the current REBUS-11. Note that this version numbering is consistent with the progression of DIF3D, the base flux solver that REBUS is built upon. The name REBUS refers to a pictorial based puzzle as the original developers were inspired by having to track thousands of unique fuel assemblies as they are inserted into the reactor, depleted, shuffled, discharged, and reprocessed.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2023

Neutronics Calculations for EBR-II Run 137E

Neutronics calculations were performed for EBR-II Run 137E using MC2 -3/TWODANT, DIF3D, and PERSENT codes. Parameters required for system safety analysis such as the reactivity feedback coefficients and core power distribution were evaluated. The detailed results are provided in the accompanying Excel spreadsheet Run_137E_Neutronics_Results.xls. This report describes the core model of EBR-II at the beginning of Run 137E, the computation methodology, and how to interpret the results listed in the Excel spreadsheet.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2021

GOTHIC Aerosol Source Depletion Studies

Pacific Northwest National Laboratory has continued work to develop an aerosol-laden flow modeling capability with the Generation of Thermal Hydraulic Information in Containment (GOTHIC TM ) computer code to perform simulations for thermal hydraulic conditions and aerosol transport and deposition in spent fuel casks. This report describes our recent work to expand our model that was originally developed in 2019 to allow for thermal characterization, carrier gas flow determination, and tracking of particulate behavior throughout the entire canister volume. This was achieved through conversion of decay heat source in all fuel tubes within the model from heaters to thermal conductors in GOTHIC as well as remeshing the internal volume of the canister. These model improvements allowed for simulations of base case scenarios for comparison of similar efforts by our collaborators with alternative tools at Sandia and Oak Ridge. Several characteristics of the GOTHIC code were elucidated by this effort and were documented for consideration in how the code is used for this effort moving forward. For instance, treatment of particle size distributions, initial spatial distributions, mesh gradients, and minimum volume concentrations were identified as important issues for consideration. Ultimately, the current version of the code is capable of tracking temperatures, flow rates, and particle behavior throughout the canister internal volume. This report presents predicted depletion times and preferred deposition patterns within the canister. Additionally, this report documents code performance and insights for further model development and use for comparisons to other codes and for predicting system behavior in experiments being performed and planned by our collaborators.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2021

Two-Phase Turbulence Statistics from High Fidelity Dispersed Droplet Flow Simulations in a Pressurized Water Reactor (PWR) Sub-Channel with Mixing Vanes

In the dispersed flow film boiling regime (DFFB), which exists under post-LOCA (loss-of-coolant accident) conditions in pressurized water reactors (PWRs), there is a complex interplay between droplet dynamics and turbulence in the surrounding steam. Experiments have accredited particular significance to droplet collision with the spacer-grids and mixing vane structures and their consequent positive feedback to the heat transfer recorded in the immediate downstream vicinity. Enabled by high-performance computing (HPC) systems and a massively parallel finite element-based flow solver—PHASTA (Parallel Hierarchic Adaptive Stabilized Transient Analysis)—this work presents high fidelity interface capturing, two-phase, adiabatic simulations in a PWR sub-channel with spacer grids and mixing vanes. Selected flow conditions for the simulations are informed by the experimental data found in the literature, including the steam Reynolds number and collision Weber number (Wec={40,80}), and are characteristic of the DFFB regime. Data were collected from the simulations at an unprecedented resolution, which provides detailed insights into the continuous phase turbulence statistics, highlighting the effects of the presence of droplets and the comparative effect of different Weber numbers on turbulence in the surrounding steam. Further, axial evolution of droplet dynamics was analyzed through cross-sectionally averaged quantities, including droplet volume, surface area and Sauter mean diameter (SMD). The downstream SMD values agree well with the existing empirical correlations for the selected range of Wec. The high-resolution data repository from the simulations herein is expected to be of significance to guide model development for system-level thermal hydraulic codes.

Saini, Nadish↗
DOE OSTI2019

Non-Intrusive Inference Reduced Order Model for Fluids Using Deep Multistep Neural Network

In this effort we propose a data-driven learning framework for reduced order modeling of fluid dynamics. Designing accurate and efficient reduced order models for nonlinear fluid dynamic problems is challenging for many practical engineering applications. Classical projection-based model reduction methods generate reduced systems by projecting full-order differential operators into low-dimensional subspaces. However, these techniques usually lead to severe instabilities in the presence of highly nonlinear dynamics, which dramatically deteriorates the accuracy of the reduced-order models. In contrast, our new framework exploits linear multistep networks, based on implicit Adams–Moulton schemes, to construct the reduced system. The advantage is that the method optimally approximates the full order model in the low-dimensional space with a given supervised learning task. Moreover, our approach is non-intrusive, such that it can be applied to other complex nonlinear dynamical systems with sophisticated legacy codes. We demonstrate the performance of our method through the numerical simulation of a two-dimensional flow past a circular cylinder with Reynolds number Re = 100. The results reveal that the new data-driven model is significantly more accurate than standard projection-based approaches.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Abundances and Transients from Neutron Star–White Dwarf Mergers

Abstract We systematically investigate the mergers of neutron star–white dwarf binaries from beginning to end, with a focus on the properties of the inflows and outflows in accretion disks and their electromagnetic emissions. Using population synthesis models, we determine a subset of these binaries in which the white dwarf companion undergoes unstable mass transfer and complete tidal disruption, forming a large accretion disk around the neutron star. The material evolves according to a one-dimensional advection-dominated accretion-disk model with nuclear burning, neutrino emissions, and disk-surface wind ejection. The extreme dynamics of the entire process have proven difficult to analyze, and thus currently, the properties are poorly understood. The outflows from the mergers are iron- and nickel-rich, giving rise to optical and infrared emissions powered by the decay of the radioactive iron-type isotopes, calculated via the SuperNu light-curve code. We find these systems capable of powering bright, yet short-lived, optical transients with the potential to power gamma-ray bursts.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Comparing the Performance of Julia on CPUs versus GPUs and Julia-MPI versus Fortran-MPI: a case study with MPAS-Ocean (Version 7.1)

Abstract. Some programming languages are easy to develop at the cost of slow execution, while others are fast at runtime but much more difficult to write. Julia is a programming language that aims to be the best of both worlds – a development and production language at the same time. To test Julia's utility in scientific high-performance computing (HPC), we built an unstructured-mesh shallow water model in Julia and compared it against an established Fortran-MPI ocean model, the Model for Prediction Across Scales–Ocean (MPAS-Ocean), as well as a Python shallow water code. Three versions of the Julia shallow water code were created: for single-core CPU, graphics processing unit (GPU), and Message Passing Interface (MPI) CPU clusters. Comparing identical simulations revealed that our first version of the Julia model was 13 times faster than Python using NumPy, where both used an unthreaded single-core CPU. Further Julia optimizations, including static typing and removing implicit memory allocations, provided an additional 10–20× speed-up of the single-core CPU Julia model. The GPU-accelerated Julia code was almost identical in terms of performance to the MPI parallelized code on 64 processes, an unexpected result for such different architectures. Parallelized Julia-MPI performance was identical to Fortran-MPI MPAS-Ocean for low processor counts and ranges from 2× faster to 2× slower for higher processor counts. Our experience is that Julia development is fast and convenient for prototyping but that Julia requires further investment and expertise to be competitive with compiled codes. We provide advice on Julia code optimization for HPC systems.

54 ENVIRONMENTAL SCIENCES↗