Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Performance benchmark”

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,711 records · Page 95
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Measurement of the Neutron Magnetic Form Factor at Large Momentum Transfer Using the Super-Bigbite Apparatus in Jefferson Lab Hall-A

The nucleon elastic electromagnetic form factors help us study the electromagnetic structure of the nucleon, benchmark theoretical models, and improve our understanding of non-perturbative quantum chromodynamics and confinement. The Nobel Prize-winning electron-nucleon scattering experiments by Robert Hofstadter and collaborators in the 1950s at Stanford High Energy Physics Lab were the first nucleon form factor measurements performed using leptonic probes. The Super Bigbite Spectrometer (SBS) program at Hall-A of Jefferson Lab represents the latest efforts to measure nucleon form factors. This ambitious program aims to significantly extend the current data set in terms of square momentum transfer (Q2) with high precision. The advent of novel detector technologies, like Gas Electron Multipliers (GEM), which provide excellent position resolution (< 100 ?m) while withstanding high background particle rates (several hundred MHz/cm2) over a large active area, has paved the way for open-geometry, moderate solid angle spectrometers, which are central to all form factor experiments in the SBS program. The first experimental run group in SBS ran successfully between September 2021 and February 2022, collecting data for the measurement of the magnetic form factor of the neutron Gn M at five squared momentum-transfer values: 3.0, 4.5, 7.5, 9.8, and 13.5 (GeV /c)2. This extends the existing high-precision data for Gn M by about a factor of four. The ratio technique was used, which involved the simultaneous measurement of exclusive quasielastic scattering of D(e,e?n)p and D(e,e?p)n from a deuterium target. Pre-preliminary results for D(e,e?n)p and D(e,e?p)n quasi-elastic ratio, and the neutron magnetic form factor Gn M , for Q2 points 3.0, 9.8, and 13.5 (GeV /c)2 are presented.

Mudiyanselage, Anuruddha Rathnayake↗
DOE OSTI2024

Polysilaketals: High‐Performance Polyether‐Based Electrolytes with Tunable Disubstituted Silane Linkers

Abstract Polymer electrolytes exhibit higher energy density and improved safety in lithium‐ion batteries relative to traditionally used liquid electrolytes but are currently limited by their lower electrochemical performance. Aiming to access polymer electrolytes with competitive electrochemical properties, we developed the anionic ring‐opening polymerization (AROP) of cyclic silaketals to synthesize amorphous silicon‐containing polyether‐based electrolytes with varying substituent bulk of the general formula [OSi(R) 2 (CH 2 CH 2 O) 2 ] n (R=alkyl, phenyl). As opposed to previously reported uncontrolled polycondensation routes toward low molecular weight polysilaketals, AROP allows access to targeted molecular weights above the entanglement threshold of the polymers. The polysilaketal with the lowest steric bulk (P(OSi Me,Me ‐2EO)) exceeds the conductivity of poly(ethylene oxide) (PEO), a leading polymer electrolyte. To the best of our knowledge, this is the first solid polymer electrolyte to achieve this benchmark. Steric bulk in polysilaketals was found to impart stability and two bulkier polysilaketals, P(OSi Et,Et ‐2EO) and P(OSi Me,Ph ‐2EO), exhibited higher current fractions than PEO over a wide range of salt loadings. Moreover, the efficacy of P(OSi Et,Et ‐2EO) was competitive with that of PEO. Taken together, the tunable and competitive electrochemical properties of polysilaketals validate the systematic incorporation of silyl groups as a strategy to access high performance polymer electrolytes.

Rugh, Haley J.↗
DOE OSTI2023

High-Fidelity Numerical Investigation on Elucidating Sodium Heat Transfer Characteristics for 37-Pin Wire-Wrapped Fuel Bundle in the PLANDTL Facility

This study involved a Reynolds-averaged Navier-Stokes- (RANS-) based computational fluid dynamics (CFD) analysis of the 37-pin wire-wrapped fuel bundle of the PNC Plant dynamics test loop (PLANDTL) facility. Previously, mainly the hydrodynamic phenomena of the wire-wrapped fuel bundle were analyzed, but the present study additionally included heat transfer analysis through conjugate heat transfer. The main purpose of the study was to benchmark the experimental data of the PLANDTL 37-pin wire-wrapped fuel bundle to investigate the heat transfer phenomena. In addition, the aim was to verify the accuracy of the RANS-based CFD analysis method using the STAR-CCM+ simulation software in comparison with the experimental data. The grid used for verification was an innovative grid system consisting of hexahedra using Fortran-based code. The development of the RANS-based CFD methodology included grid sensitivity analysis, turbulence model sensitivity analysis, and turbulent Prandtl number sensitivity analysis. Information on the temperature, mass flow rate, and area of the CFD results for each subchannel was provided for the top of the heated section and is expected to serve as a reference for future studies aiming to perform the validation and verification of a PLANDTL facility. In addition, the dependence of the peak temperature on the azimuth angle of each pin was analyzed.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2022

FORCE Development Status Update: Vertical Integration and Benchmarking of System Dynamics

Recent efforts to establish effective models for grid energy analysis, especially given the increase in variable renewable energy (VRE) sources and the economic challenges faced by traditional nuclear energy, have generated new technological considerations. One effort to improve the economic viability of nuclear power involves investigating integrated energy systems (IES) which include secondary energy systems that introduce flexibility and secondary market possibilities to existing and perceived future nuclear energy generation technologies. To analyze the technical and economic potential of IES, the Framework for Optimization of Resources and Economics (FORCE) tool suite was developed through a collaboration among national laboratories. Within the FORCE tool suite, the Holistic Energy Resource Optimization Network (HERON) provides algorithms for analyzing the long-term viability of potential IES technologies, while HYBRID provides algorithms and models to achieve high-resolution analysis of coupling physics over a short time period. Continued maturing of the FORCE tool suite requires further interconnections between the various tools in the suite in order to ensure consistent analysis. Analyses performed by applying HYBRID to transient process modeling should be easily harvestable as inputs to HERON analyses. The first item in this status update is a demonstration of an automated data pipeline for loading data from HYBRID into HERON analyses. Application of HYBRID results to HERON, as part of using the FORCE tool suite, relies on robust modeling assumptions for the various models included in HYBRID. The second result of this status update is the benchmarking and validation of cost and operational data, with a particular focus on natural gas energy generators. These generators are benchmarked with a focus on contrasting them with proposed thermal energy storage (TES) technologies.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗
DOE OSTI2021

$\mathrm{PPT}$-Multicore: performance prediction of Open$\mathrm{MP}$ applications using reuse profiles and analytical modeling

In this report we present PPT-Multicore, an analytical model embedded in the Performance Prediction Toolkit (PPT) to predict parallel applications’ performance running on a multicore processor. PPT-Multicore builds upon our previous work towards a multicore cache model. We extract LLVM basic block labeled memory trace using an architecture-independent LLVM-based instrumentation tool only once in an application’s lifetime. The model uses the memory trace and other parameters from an instrumented sequentially executed binary. We use probabilistic and computationally efficient reuse profiles to predict the cache hit rates and runtimes of OpenMP programs’ parallel sections. We model Intel’s Broadwell, Haswell, and AMD’s Zen2 architectures and validate our framework using different applications from PolyBench and PARSEC benchmark suites. The results show that PPT-Multicore can predict cache hit rates with an overall average error rate of 1.23% while predicting the runtime with an error rate of 9.08%.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

A Sparse Distributed Gigascale Resolution Material Point Method

In this paper, we present a four-layer distributed simulation system and its adaptation to the Material Point Method (MPM). The system is built upon a performance portable C++ programming model targeting major High-Performance-Computing (HPC) platforms. A key ingredient of our system is a hierarchical block-tile-cell sparse grid data structure that is distributable to an arbitrary number of Message Passing Interface (MPI) ranks. We additionally propose strategies for efficient dynamic load balance optimization to maximize the efficiency of MPI tasks. Our simulation pipeline can easily switch among backend programming models, including OpenMP and CUDA, and can be effortlessly dispatched onto supercomputers and the cloud. Finally, we construct benchmark experiments and ablation studies on supercomputers and consumer workstations in a local network to evaluate the scalability and load balancing criteria. We demonstrate massively parallel, highly scalable, and gigascale resolution MPM simulations of up to 1.01 billion particles for less than 323.25 seconds per frame with 8 OpenSSH-connected workstations.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2025

Adaptive Interface-PINNs (AdaI-PINNs) for transient diffusion: Applications to forward and inverse problems in heterogeneous media

We model transient diffusion in heterogeneous materials using a novel physics-informed neural networks framework (PINNs) termed Adaptive interface physics-informed neural networks or AdaI-PINNs (Roy et al. arXiv preprint arXiv:2406.04626, 2024). AdaI-PINNs utilize different activation functions with trainable slopes tailored to each material region within the computational domain, allowing for a fully automated and adaptive PINNs approach to model interface problems with strongly and weakly discontinuous solutions. To enhance its performance in highly heterogeneous transient diffusion systems, we prescribe a suite of robust practices, including appropriate non-dimensionalization of equations, a biased sampling method, Glorot initialization, and the hard enforcement of boundary and initial conditions. Here we evaluate the efficacy of the proposed method on several benchmark forward and inverse problems. Comparative studies on one-dimensional and two-dimensional benchmark problems reveal that the modified AdaI-PINNs outperform its unmodified counterpart, achieving root-mean-square errors that are at least two orders of magnitude better in forward problems. For inverse problems, the maximum errors in the approximated diffusion coefficients by modified AdaI-PINNs are four orders of magnitude better than those of the unmodified version. Additionally, modified AdaI-PINNs demonstrate improved stability in problems with large material mismatches.

42 ENGINEERING↗
DOE OSTI2025

Convergence study of wakefield simulations with GdfidL and ECHO3D

The interaction of charged particle beams with vacuum chamber components gives rise to electromagnetic wakefields, whose frequency-domain representation is known as beam coupling impedance. Geometric impedance arising from discontinuities and transitions in the vacuum chamber is the focus of this study. Minimizing this impedance is essential to mitigate adverse collective effects in modern storage rings operating with high-intensity particle beams. Accurate and reliable impedance simulations is a key factor of the vacuum chamber design. This paper presents the results of a convergence study of two widely used electromagnetic solvers, GdfidL and ECHO3D, applied to key vacuum-chamber components of the National Synchrotron Light Source II (NSLS-II) storage ring. Detailed comparisons are performed for several geometries, including flange absorbers, RF bellows, button-type beam position monitors, and an in-vacuum undulator (IVU). The results show notable differences in convergence and computational efficiency between the two codes. While GdfidL provides highly resolved results and serves as a common benchmark tool, ECHO3D yields consistent results with coarser meshes, significantly reducing simulation time and memory demands. Simulations with a full-geometry IVU model demonstrate that simplified taper-transition models can miss important impedance contributions. In conclusion, these findings provide practical guidelines for efficient and accurate impedance modeling to optimize design of vacuum chamber components for accelerators.

36 MATERIALS SCIENCE↗
DOE OSTI2025

The accuracy of multi-group models for nonlocal electron transport in magnetized plasmas

In the extreme conditions of inertial confinement fusion experiments, heat flow plays a vital role, but local diffusive models frequently break down and overestimate the heat flow. The situation becomes more complicated again in the significant magnetic fields generated during laser–plasma interactions or in magnetized fusion schemes. Accurate non-local and magnetized heat flow computations can be carried out using Vlasov–Fokker–Planck (VFP) simulations, but these are computationally expensive. There is, therefore, significant interest in using faster multi-group models to accurately calculate the non-local heat flow in magnetized plasmas. We benchmark two such multi-group models for calculating the heat flow, M1 and hybrid-AWBS-BGK, against diffusive models and full VFP simulations, before applying the models to realistic example test cases, both magnetized and unmagnetized. We find that the multi-group models generally perform very well for moderate non-localities up to kλmfp∼0.01, but the computational cost increases dramatically. hybrid-AWBS-BGK performs more effectively than M1 at high non-localities, up to kλmfp∼1, due to its adaptive solver and robust P1 closure, but tends to fail in very strong magnetic fields. Both codes are much faster than VFP simulations but are still slow in steep temperature gradients.

Arran, C. (ORCID:0000000286448118)↗
DOE OSTI2024

Development of Decay Heat Sensitivity Analysis Capability in SCALE/ORIGEN

In this work, a decay heat sensitivity analysis capability was developed and implemented in the ORIGEN code of the SCALE nuclear modeling and simulation suite. This capability introduces improved numerical integration schemes, which overcome the challenges associated with accurately modeling the behavior of adjoint nuclide amounts during coarse time steps for both nuclide amount and decay heat sensitivity calculations. This capability significantly improves the accuracy of calculations without compromising computational efficiency compared to the existing method. Extensive verification was conducted for various benchmark problems, including a 238 Pu decay and an irradiation problem involving 135 Xe, evaluated with both coarse and fine time grids. The results show excellent agreement with reference direct perturbation solutions, reaffirming the computational accuracy of the newly proposed numerical integration methods. Furthermore, sensitivity analyses were performed for fission product inventories ( 147 Sm, 150 Sm, 155 Gd) in pressurized water reactor UO 2 and MOX fuel assemblies. These analyses demonstrated that the ORIGEN sensitivity analysis capability can capture detailed sensitivity coefficients and underlying physics in real applications. Additionally, a decay heat sensitivity analysis for high-assay low-enriched uranium fuel, including various initial 235 U enrichment and burnup points, highlights the extended capabilities of SCALE/ORIGEN in comprehensively assessing the factors influencing total decay heat. These advancements in ORIGEN offer valuable insights for reactor analysis, fuel design, and safety assessments, especially in the context of advanced nuclear fuel development and design changes.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2026

FAIR Data and Interpretable AI Framework for Architectured Metamaterials (Final Report)

This research program established a transformative framework for the discovery and design of mechanical metamaterials, which are architected structures engineered to control physical phenomena like sound and vibration in ways natural materials cannot. To overcome the traditional reliance on trial-and-error, the project developed an interpretable Artificial Intelligence (AI) framework that moves beyond "black box" models to reveal the specific geometric patterns—such as "unit-cell templates"—that govern a material’s performance. A major breakthrough was the development of a hierarchical design method, which allows a single material to block vibrations across multiple frequency ranges simultaneously by layering patterns at different scales without them interfering with one another. This was further expanded to include irregular, graph-based designs that use spanning tree algorithms to ensure structural connectivity while allowing for customized, direction-dependent properties like stiffness and acoustic impedance. Beyond design, the project addressed the practicalities of real-world production by developing uncertainty quantification techniques that account for manufacturing defects and material variability, reducing the need for expensive physical testing by orders of magnitude. To speed up the discovery process, the team implemented Gaussian Process Regression and other surrogate models that provide accurate performance predictions at a fraction of the traditional computational cost. The AI-generated designs were successfully validated through fabrication of physical samples and wave propagation experiments, confirming their ability to accurately guide or reflect waves as predicted. By contributing these tools and high-quality FAIR benchmark datasets to the wider scientific community, this work provides a scalable foundation for advancing technologies in aerospace vibration control, medical imaging, and noise reduction.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Integration of Online Cross-Section Generation Capability with Depletion and Transient Solvers in Griffin

Griffin is a Multiphysics Object-Oriented Simulation Environment (MOOSE)-based reactor multiphysics analysis application jointly developed by Argonne and Idaho National Laboratories under the DOENE Nuclear Energy Advanced Modeling and Simulation (NEAMS) program. In FY25, an online crosssection generation capability based on the Self-Shielding Application Programming Interface (SSAPI) was demonstrated for TRISO-fueled reactor problems under steady-state conditions. This fiscal year, that capability was extended to support depletion and transient multiphysics calculations, enabling high-fidelity analyses that generate self-shielded cross sections on the fly from the actual evolving composition and temperature states rather than from pre-tabulated libraries. For depletion, a two-way coupling was established in which SSAPI computes compact-averaged self-shielded cross sections that the depletion solver then uses to advance the Bateman equations, with the updated compositions returned to SSAPI at each step; the depletion module was refactored to support both library-based and SSAPI-based cross sections, and additional logic was added to track daughter isotopes and to exclude minor isotopes for efficiency. For transient analysis, the SSAPI multigroup library was extended with the kinetics data required for time-dependent calculations, the Improved Quasi-Static (IQS) scheme was coupled with SSAPI, and several supporting capabilities were implemented, including a self-shielding treatment that lets control rods and drums move within a self-shielded model, which had previously been impossible and had ruled out rod- and drum-movement transients with on-the-fly cross sections altogether, a new mixing scheme for delayed-neutron precursor decay constants, a checkpoint-based restart workflow, and performance improvements such as pointwise cross-section interpolation and the bypassing of unnecessary Dancoff factor calculations. The implemented capabilities were verified against Serpent Monte Carlo solutions. For depletion, a prismatic pin-cell problem based on a Next Generation Nuclear Plant (NGNP) Very High Temperature Reactor benchmark showed excellent agreement, with eigenvalue differences within 200 pcm over the entire burnup range (up to 140 MWD/kgU) and fission-product and actinide inventories agreeing to within 0.8% and 2.5%, respectively; a heat-pipe microreactor assembly problem with a much higher fuel loading confirmed the same behavior and quantified the bias introduced when the multigroup equivalence effect is neglected. For transient analysis, a pin-cell problem with a step reactivity insertion and temperature feedback reproduced the analytically expected asymptotic power and showed close agreement between the direct and IQS solutions, and a two-dimensional microreactor core problem with control-drum rotation exercised the new moving-drum self-shielding treatment and demonstrated successful coupling of the online crosssection generation with both the direct and IQS transient methods. The capability was further exercised on a full-core pebble-bed problem, in which Griffin was coupled with the System Analysis Module (SAM) to simulate load-following operation of the gPBR with the Doppler feedback resolved at the TRISO fuel kernel temperature. These developments in Griffin provide a convenient, high-fidelity approach to cross-section generation for advanced thermal reactors with geometrically complex and highly heterogeneous configurations, including TRISO-fueled prismatic and pebble-bed systems, and support steady-state, depletion, and transient multiphysics calculations. They also enable self-shielded cross sections to be evaluated directly at the actual coupled state of the system, thereby establishing a foundation for high-fidelity, fully coupled multiphysics analysis of advanced reactors

Park, H.↗
DOE OSTI2025

ATF Cladding Mechanical Properties Report: Capability Demonstration

This report documents mechanical testing capability demonstration activities performed in fiscal year (FY)2025 at Oak Ridge National Laboratory(ORNL)on chromium-coated (Cr-coated) and uncoated advanced zirconium alloy claddings irradiated in the High Flux Isotope Reactor (HFIR) to approximately 4 displacements per atom (dpa), corresponding to ~13GWd/t burnup. Specimens were prepared in axial tension (ATT) and ring tension (RTT) geometries, and passive silicon carbide(SiC)thermometry (TM) was employed to determine irradiation temperatures, which averaged 38–43 °C below the 330 °C design target. Mechanical testing at ambient temperature demonstrated the expected irradiation-induced hardening, with yield strength(YS)and ultimate tensile strength(UTS)values increasing substantially relative to unirradiated counterparts. However, this strengthening was accompanied by a reduction in ductility, as indicated by lower uniform and total elongations(UE and TE). Both coated and uncoated claddings exhibited similar mechanical response, though Cr-coated specimens showed surface cracking perpendicular to the loading direction, attributable to the hardness mismatch between the coating and substrate. Fracture in all cases remained ductile, and no coating spallation was observed following HFIR irradiation. Complementary efforts were directed toward the fabrication of test specimens from commercially irradiated cladding (rod 47I, ~31.1 GWd/t average burnup). Axial sectioning and computer numerical control(CNC)machining successfully produced ATT geometries suitable for benchmarking against HFIR-irradiated specimens. This capability enables direct comparison of cladding behavior between test reactor and commercial reactor environments, thereby supporting the validation of HFIR as a surrogate irradiation platform for accident tolerant fuel (ATF) development. Once HFIR irradiations are completed in FY26,the relevant comparison tests will be completed. Collectively, the FY 2025 PIE campaign has provided mechanical performance data for irradiated advanced claddings The demonstrated capabilities support the framework for mechanical testing and further evaluations in subsequent years. These efforts will represent an important contribution toward the licensing and deployment of Cr-coated zirconium alloy cladding as a near-term ATF solution.

36 MATERIALS SCIENCE↗
DOE OSTI2022

An independent analysis of bias sources and variability in wind plant pre-construction energy yield estimation methods

The wind resource assessment community has long had the goal of reducing the bias between wind plant pre-construction energy yield assessment (EYA) and the observed annual energy production (AEP). This comparison is typically made between the 50% probability of exceedance (P50) value of the EYA and the long-term corrected operational AEP (hereafter OA P50), and is known as the P50 bias. The industry has critically lacked an independent analysis of bias reduction investigated across multiple consultants to identify the greatest sources of uncertainty and variance in the EYA process and the best opportunities for uncertainty reduction. The present study addresses this gap by benchmarking consultant methodologies against each other and against operational data at a scale not seen before in industry collaborations. We consider data from 10 wind plants and evaluate discrepancies between eight consultancies in the steps taken from estimates of gross to net energy. Consultants tend to overestimate the gross energy produced at the turbines and then compensate by further overestimating downstream losses, leading to a mean P50 bias near zero, still with significant variability among the individual wind plants. Within our data sample, we find that consultant estimates of all loss categories, except environmental losses, tend to reduce the project-to-project variability of the P50 bias. The disagreement between consultants, however, remains flat throughout the addition of losses. Finally, we find that differences in consultants’ estimates of project performance can lead to differences up to $10/MWh in the levelized cost of energy for a wind plant.

Todd, Austin C.↗
DOE OSTI2023

Insights into Waterflooding in Hydrocarbon-Bearing Nanochannels of Varying Cross Sections from Mesoscopic Multiphase Flow Simulations

Waterflooding is one of the geotechnique used to recover fuel sources from nanoporous geological formations. The scientific understanding of the process that involves the multiphase flow of nanoconfined fluids, however, has lagged, mainly due to the complex nanopore geometries and chemical compositions. To enable benchmarked flow of nanoconfined fluids, the architected geomaterials, such as synthesized mesoporous silica with tunable pore shapes and surface chemical properties, are used for designing and conducting experiments and simulations. This work uses a modified many-body dissipative particle dynamics (mDPD) model with accurately calibrated parameters to perform parametric flow simulations for studying the influences of waterflooding driven power, pore shape, surface roughness, and surface wettability on the multiphase flow in heptane-saturated silica nanochannels. Remarkably, up to 80\% reduction in the effective permeability is found for water-driven heptane flow in a baseline 4.5 nm-wide slit channel, when compared with the Hagen–Poiseuille equation. In the 4.5 nm-wide channels with architected surface roughness, the flow rate is found either higher or lower than the baseline case, depending on the shape and size of cross-sections. High wettability of the solid surface to water is essential for achieving high recovery of heptane, regardless of surface roughness. When the solid surface is less wetting or non-wetting to water, the existence of an optimal waterflooding driven power is found to allow for the highest possible recovery. A detailed analysis on the evolution of the transient water-heptane interface in those nanochannels is presented to elucidate the underlying mechanisms that impact or dictate the multiphase flow behaviors.

58 GEOSCIENCES↗
DOE OSTI2025

Systematic softening in universal machine learning interatomic potentials

Machine learning interatomic potentials (MLIPs) have introduced a new paradigm for atomic simulations. Recent advancements have led to universal MLIPs (uMLIPs) that are pre-trained on diverse datasets, providing opportunities for universal force fields and foundational machine learning models. However, their performance in extrapolating to out-of-distribution complex atomic environments remains unclear. In this study, we highlight a consistent potential energy surface (PES) softening effect in three uMLIPs: M3GNet, CHGNet, and MACE-MP-0, which is characterized by energy and force underprediction in atomic-modeling benchmarks including surfaces, defects, solid-solution energetics, ion migration barriers, phonon vibration modes, and general high-energy states. The PES softening behavior originates primarily from the systematically underpredicted PES curvature, which derives from the biased sampling of near-equilibrium atomic arrangements in uMLIP pre-training datasets. Our findings suggest that a considerable fraction of uMLIP errors are highly systematic, and can therefore be efficiently corrected. We argue for the importance of a comprehensive materials dataset with improved PES sampling for next-generation foundational MLIPs.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Proton surface exchange kinetics of perovskite triple conducting thin films for protonic ceramic electrolysis cells: BaPr 0.9 Y 0.1 O 3–δ (BPY) vs. Ba 1–x Co 0.4 Fe 0.4 Zr 0.1 Y 0.1 O 3–δ (BCFZY)

Protonic ceramic electrolysis cells (PCECs) are an attractive green H 2 production technology, given their intermediate-temperature operating range and ability to produce dry H 2 . However, PCECs will benefit from development of more efficient and durable “triple conducting” anodes where steam is split, H incorporated, and oxygen evolved. In this work, we evaluated the kinetics of the steam-splitting/H incorporation reaction on BaPr 0.9 Y 0.1 O 3–δ (BPY) in comparison to the benchmark Ba 1–x Co 0.4 Fe 0.4 Z r0.1 Y 0.1 O 3–δ (BCFZY) composition, replacing most of the transition metal elements (Co, Fe, Zr) with the lanthanide Pr. We prepared geometrically well-defined perovskite BPY and BCFZY thin films by pulsed laser deposition and performed simultaneous optical transmission relaxation and electrical conductivity relaxation measurements at 400–500 °C in 0.21 atm O 2 during switching of the steam partial pressure to isolate and compare their proton surface exchange coefficients (k). The k values of BPY were comparable to those of BCFZY and more stable over time. According to angle-resolved XPS and STEM-EDS mapping of FIB cross-sections, the surface of BPY exhibited Ba enrichment, Pr deficiency, and Si contamination. In contrast, BCFZY exhibited Ba deficiency throughout, no obvious surface segregation, and less Si contamination. The Ba segregation on the BPY film appears to have promoted steam splitting/H incorporation kinetics even though the more basic surface reacted with the acidic environmental SiO x H y . Faster kinetics observed on stoichiometric BCFZY vs. Ba-deficient BCFZY confirmed the benefit of a high A-site Ba concentration. This result contrasts with most work on perovskites applied in solid oxide electrolysis cell anodes, in which A-site segregation is considered deleterious for surface reaction kinetics.

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

Accelerating the density-functional tight-binding method using graphical processing units

Acceleration of the density-functional tight-binding (DFTB) method on single and multiple graphical processing units (GPUs) was accomplished using the MAGMA linear algebra library. Herein two major computational bottlenecks of DFTB ground-state calculations were addressed in our implementation: the Hamiltonian matrix diagonalization and the density matrix construction. The code was implemented and benchmarked on two different computer systems: (1) the SUMMIT IBM Power9 supercomputer at the Oak Ridge National Laboratory Leadership Computing Facility with 1–6 NVIDIA Volta V100 GPUs per computer node and (2) an in-house Intel Xeon computer with 1–2 NVIDIA Tesla P100 GPUs. The performance and parallel scalability were measured for three molecular models of 1-, 2-, and 3-dimensional chemical systems, represented by carbon nanotubes, covalent organic frameworks, and water clusters.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗