Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Benchmark testing”

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 1,171 records · Page 65
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Fiber-optic sensor packaging assessment and benchmarking for embedded sensing applications

Fiber optic-based sensing for non-destructive evaluation (NDE) and structural health monitoring of various infrastructure and energy assets is an increasingly important sensing scheme. Because of their immunity to electromagnetic interferences, versatility in sensing mechanisms and designs as well as capability for distributed sensing, the fiber optic sensing is of great interest for embedding within infrastructures using advanced manufacturing methods. However, the protection of fiber optic sensors during advanced manufacturing-based embedding is an important aspect for infrastructure monitoring in practice. While successful installation of fiber optic sensors has been accomplished through application of epoxy or adhesives on infrastructure surfaces, there has been only limited investigation of fiber protection packaging schemes which are compatible with in-situ repair techniques such as the metallic cold-spray. With increasing interest and importance of cold spray-based repair processes, the integration of fiber optic sensors simultaneously during the repair is equally important. In this work, we focus our efforts on investigating various types of packaging schemes for their compatibility and integrability on the surfaces of metallic structures such as oil and gas pipelines during an in-situ repair coating. First, we consider sandblasting using silica sand particles of various size for optical fibers on the surface of a metallic coupon as a useful proxy test of sensor packaging stability. Structural damage of packaging is characterized by means of optical microscopy and the optical integrity of embedded fiber optic cables is examined using optical backscatter reflectometry (OBR). In addition, initial investigation of metallic cold-spray embedded fiber optic cables under various packaging on metallic substrates is also considered with early conclusions and future directions.

Fiber-optic sensors, Sensors packaging, Metallic C↗
DOE OSTI2020

An Electrochemically Mediated Amine Regeneration Process with a Mixed Absorbent for Postcombustion CO 2 Capture

Electrochemically mediated amine regeneration (EMAR) was recently developed to avoid the use of thermal means to release CO2 captured from postcombustion flue gas in the benchmark amine process. To address concerns related to the high vapor pressure of ethylenediamine (EDA) as the primary amine used in EMAR, a mixture of EDA and aminoethylethanolamine (AEEA) was investigated. The properties of the mixed amine systems, including the absorption rates, electrolyte pH and conductivity, and CO2 capacity, were evaluated in comparison with those of solely EDA. The mixed amine system had similar properties to that of EDA, indicating no significant changes would be necessary for the future implementation of the EMAR process with mixed amines as opposed to that with just EDA. The electrochemical performance of the mixed amines in terms of the cell voltage, gas desorption rate, electron utilization, and energetics was also investigated. A 50/50 mixture of EDA and AEEA displayed the lowest energetics: ~10% lower than that of 100% EDA. With this mixture, a continuous EMAR process, in which the absorption column was connected to the electrochemical cell as the desorption stage, was tested over 100 h. The cell voltage was very stable and there was a steady gas output close to theoretical values. The desorbed gas was further analyzed and found to be 100% CO2, confirming no evaporation of the amine. The mixed absorbent composition was also characterized using titration and nuclear magnetic resonance (NMR) spectroscopy, and the results showed no amine degradation. Furthermore, these findings that demonstrate a stable, low vapor pressure absorbent with improved energetics are promising and could be a guideline for the future development of EMAR for CO2 capture from flue gas and other sources.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Techno-Economic Analysis of Green Hydrogen Energy Storage in a Cryogenic Flux Capacitor

Abstract The Cryogenic Flux Capacitor (CFC) is a cold, dense energy storage core that is being studied in the cryo-compressed, about 300 bar and 80K, region of gaseous hydrogen (GH2) storage and liquid hydrogen (LH2) region near the normal boiling point. Hydrogen storage is improved by physically bonding the molecules within the nanoscale pores of the aerogel composite blanket material. The process of bonding or debonding is governed by principles of physical adsorption (physisorption) and thermodynamics. The large surface area afforded by the nanoporous aerogel (∼1,000 m2/g) allows its storage performance to easily exceed capacities of high-pressure GH2 storage for an equivalent volume. With the integrated aerogel, subscale tests have shown that storage is increased by about 36% over a simple tank filled with GH2 at the same operating temperature and pressure. For LH2 conditions, the CFC is shown to operate at improved densities, but testing is ongoing. For the techno-economic analysis (TEA), the source of hydrogen is compared between onsite steam methane reforming (SMR) and onsite solar photovoltaic (PV) panels providing power to electrolyzers to produce green GH2. The TEA compares pure hydrogen produced at a small scale for a 25 MW power system and at a large scale in a 500 MW power system. The system allowed for hydrogen imports and exports at a set price with a tank sized for 10 hours of power production. The two power producing technologies are a combined cycle gas turbine (CCGT) and hydrogen fuel cells. The SMR system uses natural gas as an input and includes a carbon capture and storage (CCS) system. The levelized cost of electricity (LCOE), levelized cost of hydrogen (LCOH), and levelized cost of storage (LCOS) are developed based on the capital cost and operating cost of the systems. The results are shown for current costs using a 2021 benchmark and DOE projections for cost improvements by 2030. The TEA showed that onsite hydrogen generation from SMR has an LCOH of about 1.4 to 2 USD per kg over the life of the plant and the PV hydrogen production LCOH is about 5.2 to 5.5 USD per kg. The LCOS of conventional GH2 systems is estimated to be $210/MWh and cost of storage for LH2 systems is $205/MWh for fuel cell systems and $249/MWh for CCGT systems. CFC improved the LCOS of all these systems to $198/MWh, $191/MWh and $233/MWh respectively. The LCOE also improved with conventional systems between $171/MWh and $228/MWh improved by CFC to between $167/MWh and $212/MWh. Using projections for improvement in costs following DOE’s goals by 2030, green hydrogen improved to as low as $78/MWh LCOS and LCOE for conventional cases. CFC improved over conventional storage with the lowest LCOS being $62/MWh and the lowest LCOE being $73/MWh. These results correspond to an LCOH of $2/kg. Finally, the TEA shows how LCOE is improved for hydrogen conditioning and storage over conventional systems and caverns in the 10 to 50 hour range.

08 HYDROGEN↗
DOE OSTI2021

Discrete modeling of a transformer with ALEGRA

We report progress on a task to model transformers in ALEGRA using the “Transient Magnetics” option. We specifically evaluate limits of the approach resolving individual coil wires. There are practical limits to the number of turns in a coil that can be numerically modeled, but calculated inductance can be scaled to the correct number of turns in a simple way. Our testing essentially confirmed this “turns scaling” hypothesis. We developed a conceptual transformer design, representative of practical designs of interest, and that focused our analysis. That design includes three coils wrapped around a rectangular ferromagnetic core. The secondary and tertiary coils have multiple layers. The tertiary has three layers of 13 turns each; the secondary has five layers of 44 turns; the primary has one layer of 20 turns. We validated the turns scaling of inductance for simple (one-layer) coils in air (no core) by comparison to available independent calculations for simple rectangular coils. These comparisons quantified the errors versus reduced number of turns modeled. For more than 3 turns, the errors are <5%. The magnetic field solver failed to converge (within 5000 iterations) for >10 turns. Including the core introduced some complications. It was necessary to capture the core surfaces in thin grid sheaths to minimize errors in computed magnetic energy. We do not yet have quantitative benchmarks with which to compare, but calculated results are qualitatively reasonable.

42 ENGINEERING↗
DOE OSTI2026

Monitoring of Liquid Metal Reactor Heater Zones with Recurrent Neural Network Learning of Temperature Time Series

Advanced high-temperature fluid reactors (ARs), such as sodium fast reactors (SFRs) and molten salt cooled reactors (MSCRs) utilize high-temperature fluids at ambient pressure. To melt the fluid during reactor startup and prevent fluid freezing during cooldown, the thermal–hydraulic systems of such ARs include heater zones consisting of specific heaters with controllers, temperature sensors, and thermal insulation. The failure of heater zones due to insulation material degradation or improper installation, resulting in parasitic heat losses, can lead to fluid freezing. The detection of faults using a heat-transfer model is difficult because of a lack of knowledge of the experimental details. Data-driven machine learning of heater zone temperature time series offers a viable alternative. In this study, we benchmarked the performance of recurrent neural networks (RNNs) in an analysis of heat-up transient temperature time series of heater zones installed on a liquid sodium vessel. The RNN models include long short-term memory (LSTM) and gated recurrent unit (GRU) networks, as well as their bi-directional variants, BiLSTM and BiGRU. Anomalous temperature points were designated using a percentile-based threshold applied to residual fluctuations in the detrended temperature time series. Additionally, the impact of the exponentially weighted moving average (EWMA) method on detection accuracy was examined. The RNN models’ performance was assessed using precision, recall, and F 1 score metrics. Results demonstrated that RNN models effectively detect anomalies in temperature time series with the best models for each heater zone achieving F 1 scores of over 93%. To explain the variations in RNN model performance across different heater zones, we used Kullback–Leibler (KL) divergence to quantify the relative entropy between training and testing data, and the Detrended Fluctuation Analysis (DFA) to assess long-range temporal correlations. For datasets with strong long-range correlations and minimal relative entropy between training and testing data, GRU is the best-performing model. When the data exhibits weaker long-term correlations and a significant relative entropy between training and testing distributions, BiGRU shows the best performance. For the data sets with intermediate values of both KL divergence and DFA, the best performance is obtained with LSTM and BiLSTM, respectively.

gated recurrent unit↗
DOE OSTI2024

Development of Low-Cost, High-Performance, Easy-To-Apply, Non-Flammable, Inorganic Phase Change Material (PCM) Technology (Project Final Report)

This report describes a 45-months long research program focused on the development of novel, easy-to-apply, non-flammable, and high-performance inorganic phase change materials (PCMs) for building and industrial applications. The University of Massachusetts Lowell (UML) formed a world-class team consisting of researchers form InsolCorp (only N. American manufacturer of inorganic PCM systems for building applications), and a group of industrial advisors, to develop a universal/multipurpose, simple-to-manufacture and cost-effective PCM technology. The project team expects that the results of this work will spur in the future the adoption of thermal storage materials – a key building energy saving technology as identified by DOE BTO – for a variety of building envelope applications. The main goal of this project was to demonstrate a suite of low-cost, multipurpose, and durable inorganic PCM formulations with phase transition temperatures encompassing typical building applications (between +5 o C and +55 o C). The first objective was to design, fabricate, and experimentally validate a performance of inexpensive, durable, highly efficient, non-flammable, and easy to manufacture PCMs. To allow a variety of building applications, the project team focused on formulations that exhibit repeatable phase transitions between +5 o C and +55 o C. To follow the DOE BTO cost efficiency target without compromising thermal performance, our work was based on inorganic compounds (mostly salt hydrates) and their blends, which represent a fraction of the cost of most of organic PCMs with about twice as high density as well as significantly higher thermal conductivity and phase change enthalpy. The second objective was to develop easy-to-manufacture and -install packaging/encapsulation designs that are 1) a superior barrier to current state-of-the-art macro-packaging, which significantly reduces the risk of loss of hydration water and PCM leak, and 2) optimal in enhancing the heat exchange rates with the surroundings and within the PCM core to ensure complete charging/discharging of the entire PCM within the product. Finally, the project’s intend was to scale-up the fabrication process to demonstrate installation on system-scale applications, and to validate the performance under field conditions. This work aimed at developing low-cost, high-energy storage, and reliable latent heat storage technology for building applications. This development was realized by formulating and integrating the following two technology components: 1) inorganic salt hydrate based PCMs that have high latent enthalpies and are low-cost and durable, and 2) PCM encapsulation (packaging) technology that maximizes PCM concentration and enhances heat transport characteristics in the product and with the external environment/materials. High thermal storage capacity, low cost and fire resistance are key to the building market entry for PCM technology. Therefore, the project’s focus was on salt-hydrate-based formulations which satisfy all these criteria. Packaging and/or encapsulation of PCM is a key processing step. The project team recognized that a low-cost and simple-to-manufacture salt hydrate-based PCM technology holds the best chance to be successful in the building construction market, a market which is traditionally extremely sensitive to cost and where commodity thermal insulations are the benchmark for envelope-related energy saving measures. That is why, in this project, the main intention was to minimize the production cost and maximize the product energy storage density without sacrificing the PCM performance. It was achieved through: 1. Minimizing the non-PCM components (plastics, additives, packaging/encapsulation materials, etc.) because they are significantly more expensive than salt hydrates, 2. Using highly thermally conductive and lightweight PCM carrier (packaging material) to facilitate more complete phase cycling, and 3. Optimizing the thickness and minimizing air spaces in product design (such as in pouched PCM). For this purpose, our approach was to enable an easy system design, including selection of the PCM operating temperatures, optimizing the necessary heat storage capacity (by stacking together several layers of PCM products), and if needed, a synchronized usage of PCM products of different temperatures. A specially designed, robust, highly thermally conducting and highly impermeable packaging (to retain salt hydrate water during phase transition cycles) was designed and tested to increase the overall system thermal performance and durability. All PCM products developed during this project were tested in both lab scale and in full scale field conditions. It is expected that, after further developments and commercialization, the developed PCM technologies may be also applied in space conditioning, energy storage technologies, and heat transfer applications.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2021

Initial Assessment of TeF 6 Adsorption

Research and development that supports the management of off-gases from nuclear fuel reprocessing has historically been focused on the off-gas streams that arise from aqueous reprocessing technology. However, as Gen-IV reactor development pathways move toward deployment, alternative spent nuclear fuel (SNF) processing and disposition pathways have been considered more actively. This work is focused upon aspects of fluoride volatility (FV) processing. The versatility of this method for deployment against multiple types of spent fuel encourages the continued advancement of both the primary separations processes and the secondary processes, including waste treatment, material control and accountability, and engineering designs. Recent work noted that TeF 6 , the most highly volatile fluorinated compound produced during FV processing, did not have a clear abatement technology recommended in the literature. Thus, this work performed scoping tests on activated alumina and copper shot to assess whether they could be used to remove TeF 6 from gas streams that bear F 2 . Previous work on this topic was not well described in the literature and was not performed with excess F 2 in the stream as would be typical of spent fuel processing via FV. To support an understanding of the concentration of TeF 6 contacting the adsorbent beds, a series of preliminary testing identified the TeF 6 production rate and the equilibrium concentration of TeF 6 in the gas stream contacting the adsorbent. Excess F 2 was determined to not affect the ability of activated alumina to remove TeF 6 from the gas stream quickly and completely. A determination of whether excess F 2 affected the distribution depth of Te in the sorbent bed is still pending analysis of the used sorbent. Literature suggests that when activated alumina is near saturation, TeF 6 could migrate from the sorbent bed. Future testing should investigate this possibility. Copper metal did not adsorb TeF 6 in the presence of F 2 across the sorbent temperature range of 50 to 335°C. F 2 was fully removed by the copper bed. Although preliminary thermodynamics would indicate adsorption to be energetically favorable, other factors that impact adsorption (e.g., slow kinetics, excess fluorine on the copper surface, an unfavorable transition state) are likely preventing the adsorption of TeF 6 at an easily measurable rate. The work also investigated TeF 6 production from multiple forms of Te in the temperature range of 100 to 250°C. No previous study had assessed the initial reaction rates for this process. A carefully designed study allowed determination of the activation energy for the production of TeF 6 from Te metal. The substantial amount of data collected during this study merits analysis beyond what is described here. A more in-depth kinetic analysis will be pursued. Additional analytical results will provide the ability to benchmark adsorption coefficients for TeF 6 , understand the distribution of TeF 6 within the alumina bed, and better understand the effect of F 2 partial pressure on Te fluorination. The data from this report, as supplemented by these additional analyses, will be submitted to a peer-reviewed journal.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2020

The good, the bad, and the ugly: Data-driven load profile discord identification in a large building portfolio

Reducing the overall energy consumption and associated greenhouse gas emissions in the building sector is essential for meeting our future sustainability goals. Recently, smart energy metering facilities have been deployed to enable monitoring of energy consumption data with hourly or subhourly temporal resolution. This unprecedented data collection has created various opportunities for advanced data analytics involving load profiles (e.g., building energy benchmarking programs, building-to-grid integration, and calibration of urban-scale energy models). These applications often need preprocessing steps to detect daily load profile discords, such as: 1) outliers due to system malfunctions (the bad) and 2) irregular energy consumption patterns, such as those resulting from holidays (the ugly) compared to normal consumption patterns (the good). However, current preprocessing methods predominantly focus on filtering using statistical threshold values, which fail to capture the contextual discords of daily profiles. In addition, discord detection algorithms in building research are often aimed at finding individual building-level discords, which are not suitable at a large scale. Thus, here, we develop a method for automated load profile discord identification (ALDI) in a large portfolio of buildings (more than 100 buildings). Specifically, ALDI 1) uses the matrix profile (MP) method to quantify the similarities of daily subsequences in time series meter data, 2) compares daily MP values with typical-day MP distributions using the Kolmogorov-Smirnov test, and 3) identifies daily load profile discords in a large building portfolio. We evaluate ALDI using the metering data of both an academic campus and a residential neighborhood. Our results demonstrate that ALDI efficiently discovers measurement errors by system malfunctions and low energy consumption days in the academic campus portfolio, and it detects unique load shape patterns likely driven by occupant behavior and extreme weather conditions in the residential neighborhood.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2020

3D CFD Model Validation Using Benchmark Data of 1/16th Scaled VHTR Upper Plenum and Development of Wall Heat-Transfer Correlation For Laminar Flow

With support from the U.S. Department of Energy-Office of Nuclear Energy’s (DOE-NE’s) Nuclear Energy Advanced Modeling and Simulation (NEAMS) program, an effort has been pursued to support high-temperature gas-cooled reactor (HTGR) technology development and its modeling and simulation needs. There is a particular need for advanced modeling and simulation tools to predict thermal-fluid behavior in the nuclear reactor primary system, especially in the core and the lower and upper plena, during safety-related transients. In this report, two main such activities are presented relevant to the HTGRs: (1) three-dimensional (3D) computational fluid dynamics (CFD) validation using benchmark data from the upper plenum of Texas A&M University’s 1/16th scaled very-high-temperature gas-cooled reactor (VHTR), and (2) development of wall heat-transfer correlation for laminar flow in a wall-heated pipe. The CFD tool validation exercises can be helpful to choose the models and CFD tools to simulate and design specific components of the HTRGs such as upper plenum where jet mixing is a complex phenomenon. In a loss of forced circulation event, the laminar flow can be observed during the development of natural circulation flow. This work includes the development and validation of heat transfer correlations for laminar flow using the Nek5000 CFD code due to limited available experimental data for laminar flow conditions to guide low-order models (1D). In this report, the flow characteristics of a single isothermal jet discharging into the upper plenum was investigated using the Nek5000 Large-Eddy Simulation (LES) CFD tool. Several numerical simulations were performed for various jet-discharged Reynolds numbers ranging from 3,413 to 12,819. A grid-independent study was performed. The numerical results of mean velocity, root-mean-square fluctuating velocity, and Reynolds stress were compared against the benchmark data. Good agreement was obtained between simulated and measured data for axial mean velocities, except near the upper plenum hemisphere. The maximum predicted errors for axial mean velocities at various normalized coolant channel diameter heights of 1, 5, and 10 are 1.56%, 1.88%, and 3.82%, respectively. In addition, the predicted root-mean-square fluctuating velocity and Reynolds stress are qualitatively in agreement with the experimental data. The Nek5000 code was used to develop wall-heat transfer correlation for laminar flow in a cylindrical tube. Several simulations were performed for various Reynolds flow and wall-heat fluxes. A new heat transfer correlation was developed using data from Nek5000 simulation results and regression functions in Matlab. The developed heat transfer correlation is valid for various Reynolds flows from 200 to 2000. The predicted R² value for model fit was 0.875, which ensures that 87.5% of the model data lies on the Nek5000 data. Moreover, a machine learning (ML) tool was used to train and test the Nek5000 data. A good fit of the ML-based model was observed with the test data.

21 SPECIFIC NUCLEAR REACTORS AND ASSOCIATED PLANTS↗
DOE OSTI2020

The Strong Impact of Weak Horizontal Convergence on Continental Shallow Convection

Idealized large-eddy simulation (LES) is a basic tool for studying three-dimensional turbulence in the planetary boundary layer. LES is capable of providing benchmark solutions for parameterization development efforts. However, real small-scale atmospheric flows develop in heterogeneous and transient environments with locally varying vertical motions inherent to open multiscale interactive dynamical systems. These variations are often too subtle to detect them by state-of-the-art remote and in situ measurements, and are typically excluded from idealized simulations. The present study addresses the impact of weak [i.e., O(10-6) s-1] short-lived low-level large-scale convergence/divergence perturbations on continental shallow convection. The results show a strong response of shallow nonprecipitating convection to the applied weak large-scale dynamical forcing. Evolutions of CAPE, mean liquid water path, and cloud-top heights are significantly affected by the imposed convergence/divergence. In contrast, evolving cloud-base properties, such as the area coverage and mass flux, are only weakly affected. To contrast those impacts with microphysical sensitivity, the baseline simulations are perturbed assuming different observationally based cloud droplet number concentrations and thus different rainfall. For the tested range of microphysical perturbations, the imposed convergence/divergence provides significantly larger impact than changes in the cloud microphysics. Simulation results presented here provide a stringent test for convection parameterizations, especially important for large-scale models progressing toward resolving some nonhydrostatic effects.

Meteorology & Atmospheric Sciences↗
DOE OSTI2022

Oak Ridge National Laboratory Pilot Demonstration of an Attestation and Anomaly Detection Framework using Distributed Ledger Technology for Power Grid Infrastructure

This report summarizes the design and pilot demonstration of a framework called Grid Guard that was created to provide increased data and device trustworthiness to electric grid devices by leveraging distributed ledger technology (DLT), specifically blockchain. Grid Guard contains a combination of core cryptographic methods such as the secure hash algorithm (SHA), and asymmetric cryptography, private permissioned blockchain, baselining configuration data, consensus algorithm (Raft) and the Hyperledger Fabric (HLF) framework. The system implements a low energy, fast, and robust enhancement to system trustworthiness within and across electric grid systems such as substations, control centers and metering infrastructures. Blockchain is a distributed database structured that provides a practically unalterable (immutable) timeline of stored transactions. By relying on hashing and the Raft consensus algorithm, if an entity tries to illegitimately alter a record at one instance of the database the other ledger nodes are not altered. They work to cross-reference each other and easily locate any incorrectly added data and remove it. The bulk raw data is stored in an off-chain storage (outside of the blockchain ledger) and a hash of this baseline data is stored in the Blockchain ledger via hashing windows of time-series and configuration data, after aggregation and filtering. The bulk off-chain data repository is then considered to be trust-anchored using the hashes stored in the blockchain. To secure the electric grid testbed devices and data, device configuration baselines were compared to those baselines that had been previously stored in the ledger. Statistical baselines for device configurations, network communication patterns, and high-speed sensor data are calculated and then stored off-chain and hashes stored in the ledger. Measurements such as three-phase voltage and current, frequency, breaker status, protection scheme settings, network configuration settings (and other device configuration artifacts) and network traffic features (packet interarrival times) are compared every minute or other selected time windows. During phase 1 of the Grid Guard DLT project different DLT technologies were studies, and an assessment was performed on DLT technology vulnerabilities, uses, and key characteristics. DLT consensus protocols were studies (e.g., RAFT, named after Reliable, Replicated, Redundant, And Fault-Tolerant). Also, cryptography, public, private and permissioned or permissionless systems were assessed. Grid Guard implements a permissioned private DLT. Consensus algorithm selection and choice of DLT implementation depended heavily on the use-case. For this use-case, parameters were selected to measure performance and existing tools for assessment. Benchmarking was performed theoretically and practically. During phase 2 hashed transactions/blocks were inserted into the ledger every second. During phase 2 of the Grid Guard DLT project, a prototype framework was developed and demonstrated for attestation of critical substation devices and data using precision timing systems that use PTP and IRIG-B protocols) on a testbed of operational devices that emulated a distribution substation, control center, and power metering infrastructure using real Operational Technology (OT). The testbed includes OT devices such as protective relays, human machine interfaces (HMI), and power meters. To determine when to collect and compare system and network baselines, an initial examination of an anomaly detection capability to identify malicious manipulation of data streams was conducted. The resulting anomaly detection was demonstrated in a set of experiments and leveraged to trigger device artifact attestation checks. Attestation checks occur against device configuration baselines when compared with the immutable blockchain-stored baselines, which provided a cryptographically supported means by which to store baselines. The electrical substation-grid testbed was created to test the Grid Guard framework. The testbed emulates the operations of a portion of a power grid and SCADA systems as closely as possible. The testbed integrates real protocols, mainly IEC 61850 standard protocols, such as the Sampled Value (SV) and the GOOSE protocols. The testbed also supports DNP3 and other layer 2 and layer 3 protocols such as Telnet, SSH, SFTP/FTP and other proprietary protocols needed to connect to industrial control system equipment. The testbed emulates real power conditions using the OpalRT hardware-in-the-loop (HIL) device which can create fault situations that cannot be easily tested on real systems. The electrical substation-grid testbed was created using real measurement, communication, and protection devices that electrical utilities commonly use.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2025

Investigation of tearing mode stability near ideal stability boundaries via asymptotic matching techniques

A number of improvements to the TJ toroidal tearing mode code [Fitzpatrick, Phys. Plasmas 31, 102507 (2024)] are documented. The TJ code is also successfully benchmarked against the STRIDE toroidal tearing mode code [Glasser and Koleman, Phys. Plasmas 25, 082502 (2018)]. Finally, the new capabilities of the TJ code are used to investigate the stability of tearing modes in tokamak plasmas as an ideal stability boundary, associated with either an external-kink or an internal-kink mode, is approached. All elements of the tearing stability matrix are found to tend to infinity as an ideal stability boundary is approached. Furthermore, as the stability boundary is approached, the eigenfunctions of the various tearing modes in the plasma, which are decoupled by sheared plasma rotation, are all found to morph into that of the marginally stable ideal mode. However, the growth-rates and real frequencies of the various “ideal-tearing-modes” are different from one another. Moreover, the growth-rate of the ideal-tearing-mode that reconnects magnetic flux at the rational surface that lies closest to the edge of the plasma is the one that tends to a very large value as the stability boundary is approached. A relatively simple test for ideal stability that is capable of detecting stability boundaries for external-kink and internal-kink modes, even in the presence of a very close-fitting ideal wall, is described and verified.

FOS: Physical sciences↗
DOE OSTI2023

Evaluation of OpenAI Codex for HPC Parallel Programming Models Kernel Generation

We evaluate AI-assisted generative capabilities on fundamental numerical kernels in high-performance computing (HPC), including AXPY, GEMV, GEMM, SpMV, Jacobi Stencil, and CG. We test the generated kernel codes for a variety of language-supported programming models, including (1) C++ (e.g., OpenMP [including offload], OpenACC, Kokkos, SyCL, CUDA, and HIP), (2) Fortran (e.g., OpenMP [including offload] and OpenACC), (3) Python (e.g., numpy, Numba, cuPy, and pyCUDA), and (4) Julia (e.g., Threads, CUDA.jl, AMDGPU.jl, and KernelAbstractions.jl). We use the GitHub Copilot capabilities powered by the GPT-based OpenAI Codex available in Visual Studio Code as of April 2023 to generate a vast amount of implementations given simple + + prompt variants. To quantify and compare the results, we propose a proficiency metric around the initial 10 suggestions given for each prompt. Results suggest that the OpenAI Codex outputs for C++ correlate with the adoption and maturity of programming models. For example, OpenMP and CUDA score really high, whereas HIP is still lacking. We found that prompts from either a targeted language such as Fortran or the more general purpose Python can benefit from adding code keywords, while Julia prompts perform acceptably well for its mature programming models (e.g., Threads and CUDA.jl). We expect for these benchmarks to provide a point of reference for each programming model's community. Overall, understanding the convergence of large language models, AI, and HPC is crucial due to its rapidly evolving nature and how it is redefining human-computer interactions.

Godoy, William↗
DOE OSTI2022

Large‐Eddy Simulations of a Convection Cloud Chamber: Sensitivity to Bin Microphysics and Advection

Abstract Bin microphysics schemes are useful tools for cloud simulations and are often considered to provide a benchmark for model intercomparison. However, they may experience issues with numerical diffusion, which are not well quantified, and the transport of hydrometeors depends on the choice of advection scheme, which can also change cloud simulation results. Here, an atmospheric large‐eddy simulation model is adapted to simulate a statistically steady‐state cloud in a convection cloud chamber under well‐constrained conditions. Two bin microphysics schemes, a spectral bin method and the method of moments, as well as several advection methods for the transport of the microphysical variables are employed for model intercomparison. Results show that different combinations of microphysics and advection schemes can lead to considerable differences in simulated cloud properties, such as cloud droplet number concentration. We find that simulations using the advection scheme that suffers more from numerical diffusion tends to have a smaller droplet number concentration and liquid water content, while simulation with the microphysics scheme that suffers more from numerical diffusion tends to have a broader size distribution and thus larger mean droplet sizes. Sensitivities of simulations to bin resolution, spatial resolution, and temporal resolution are also tested. We find that refining the microphysical bin resolution leads to a broader cloud droplet size distribution due to the advection of hydrometeors. Our results provide insight for using different advection and microphysics schemes in cloud chamber simulations, which might also help understand the uncertainties of the schemes used in atmospheric cloud simulations.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2021

Developing and Testing a Novel Stochastic Ice Microphysics Parameterization for Cloud and Climate Models Using ARM Field Campaign Data (Final Progress Report)

The major goals of this project were: 1) to use field campaign measurements from DOE’s Atmospheric Radiation Measurement (ARM) program to characterize variability of important parameters describing properties of ice particles in the atmosphere; 2) based on this observational analysis, to develop a parameterization scheme for weather and climate models that stochastically varies these parameters, and implement the new scheme into a weather model called the Weather Research and Forecasting model (WRF); 3) to use WRF coupled with the new stochastic scheme to simulate ARM field campaign thunderstorm cases and analyze how accounting for this parameter variability affects the model simulations. This work was performed jointly between the National Center for Atmospheric Research, University of Oklahoma, and University of Utah. To accomplish these goals, we extended an approach previously developed to characterize the variability in the size distribution of ice particles to parameters that are explicitly represented in models (i.e., relationships between ice particle mass and size, and between particle fall velocity and size). Our project was, to our knowledge, the first to apply observationally-constrained estimates of this parameter variability describing mass-size and fall velocity-size in a modeling framework. Our results showed efficacy of the approach, evaluated using ARM observations. Similarly, to our knowledge, work in this project was the first to propose and evaluate in detail a stochastic approach for unresolved turbulent mixing in high-resolution model simulations against detailed, benchmark large eddy simulations and ARM observations. Results showed some promising behavior, particularly with increased mixing and dilution of air in thunderstorm cores with surrounding environmental air, bringing the stochastic simulations closer to the benchmark large eddy simulations; however, results were somewhat degraded using stochastic mixing compared to observations from the AMIE/DYNAMO field campaign. This project also further refined and applied a modeling methodology called “piggybacking” that can robustly separate dynamical and thermodynamic impacts of model changes, and comparison studies of different models based on cases developed from ARM observations. Finally, this project directly supported three graduate students who completed their PhDs as well as a postdoctoral research fellow.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Solar neutrino detection in liquid xenon detectors via charged-current scattering to excited states

We investigate the prospects for real-time detection of solar neutrinos via the charged-current neutrino-nucleus scattering process in liquid xenon time projection chambers. We use a nuclear shell model, benchmarked with experimental data, to calculate the cross sections for populating specific excited states of the cesium nuclei produced by neutrino capture on Xe 131 and Xe 136 . The shell model is further used to compute the decay schemes of the low-lying 1 + excited states of Cs 136 , for which there is sparse experimental data. We explore the possibility of tagging the characteristic deexcitation γ rays/conversion electrons using two techniques: spatial separation of their energy deposits using event topology and their time separation using delayed coincidence. The efficiencies in each case are evaluated within a range of realistic detector parameters. We find that the topological signatures are likely to be dominated by radon backgrounds, but that a delayed-coincidence signature from long-lived states predicted in Cs 136 may enable background-free detection of CNO neutrino interactions in next-generation experiments with smaller uncertainty than current measurements. We also estimate the sensitivity as a function of exposure for detecting the solar-temperature-induced line shift in Be 7 neutrino emission, which may provide a new test of solar models.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

Establishment of a Vertically Integrated Domestic Manufacturing Process for Production of Substrates Needed for Manufacture of Gas Diffusion Layers

In this project, AvCarb, LLC evaluated the baseline performance metrics of commercial carbon veils and their corresponding Gas Diffusion Layers (GDLs) with the goal of establishing an optimized, vertically integrated production system for wet-laid nonwoven substrates used in gas diffusion media for electrochemical energy storage and conversion devices. Mechanical testing and microstructural characterization were conducted and used to develop a multiscale computational model capable of simulating and predicting the performance of GDLs in fuel cells. Although the project successfully generated foundational transport and modeling data, it was terminated prior to identifying the critical GDL design parameters necessary for full optimization. The program aimed to improve carbon veil fabrication through enhanced fiber dispersion, fiber-fiber adhesion control, and improved web formation, enabling the production of high-quality, uniform substrates. Simulations were intended to guide mixing and solution delivery system design and process conditions, followed by production-scale trials to evaluate fiber dispersion, web uniformity, and mechanical robustness. At full deployment, the proposed production line would have been capable of producing approximately 650,000 m² of carbon veil annually. This capability remains strategically important, as the United States currently lacks a domestic source of wet-laid nonwoven carbon substrates that satisfy the stringent quality requirements for fuel-cell GDLs and electrolyzers representing an ongoing supply-chain vulnerability. Beyond supply-chain benefits, the project established a robust benchmarking dataset for existing commercial carbon veils while advancing next-generation material concepts targeting improved performance and manufacturing consistency.

Olson, Cynthia Lemay↗
DOE OSTI2020

First-principles calculation of the configurational energy density of states for a solid-state ion conductor with a variant of the Wang and Landau algorithm

In this work, a variant of the Wang and Landau algorithm for calculation of the configurational energy density of states is proposed. The algorithm was developed for the purpose of using first-principles simulations, such as density functional theory, to calculate the partition function of disordered sublattices in crystal materials. The expensive calculations of first-principles methods make a parallel algorithm necessary for a practical computation of the configurational energy density of states within a supercell approximation of a solid-state material. The algorithm developed in this work is tested with the two-dimensional (2d) Ising model to bench mark the algorithm and to help provide insight for implementation to a materials science application. Tests with the 2d Ising model revealed that the algorithm has good performance compared to the original Wang and Landau algorithm and the 1/$\textit{t}$ algorithm, in particular the short iteration performance. Further, a proof of convergence is presented within an adiabatic assumption, and the analysis is able to correctly predict the time dependence of the modification factor to the density of states. The algorithm was then applied to the lithium and lanthanum sublattice of the solid-state lithium ion conductor Li 0.5 La 0.5 TiO 3 . This was done to help understand the disordered nature of the lithium and lanthanum. The results find, overall, that the algorithm performs very well for the 2d Ising model and that the results for Li 0.5 La 0.5 TiO 3 are consistent with experiment while providing additional insight into the lithium and lanthanum ordering in the material. The primary result is that the lithium and lanthanum become more mixed between layers along the c axis for increasing temperature. In part, the simulation of the disordered Li 0.5 La 0.5 TiO 3 system serves as a benchmark for what size systems are currently and in the near future practical to calculate with density functional theory methods.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗