Publications

Abstract not provided.

This document provides an overview of re-start efforts at Sandia National Laboratories (SNL), National Solar Thermal Test Facility (NSTTF), for the SNL Molten Salt Test Loop (MSTL). MSTL is one of the world’s only industrial-scale molten salt test systems, capable of testing commercial-scale thermal hydraulic equipment at an industrial scale, required for confident bankability of products. This work was facilitated as part of the Phase 1 re-start effort for facilitating forensics assessments to inform final refurbishments and upgrades costs and schedule information for revitalizing MSTL for use during a separate Phase 2 effort. This work was facilitated to further de-risk a Phase 2 effort pertaining to the tank and the pump, which were considered during Phase 1 to be high-risk items for a successful re-start. Impact from a successful re-start will allow SNL and other industrial molten salt systems companies, such as those in the Market Utilization Report provided in the Phase 1 submission package. An objective of this document is to also ensure longevity of utilization of MSTL after re start.

Size effects are an unavoidable nuisance in inertial navigation using sensors which are not co-located at the navigational point of interest. When estimating transforms between the navigation point and sensor locations, some trajectories preclude observation of all model parameters. Regularization is proposed to avoid over-fitting size-effects models. The result yields robust size effects compensation in other regions of flight.

The future of High Performance Computing (HPC) is carved out for all high-performance facilities. As computer processing increases exponentially, especially with the demand for AI training, more speed will equate to more computing density and ultimately more power draw and resource demand. Sandia National Laboratories has not seen or known of anything that would lead us to think this might change in the next 5 to 10 years.

Abstract not provided.

Magnetic reconnection is a fundamental plasma physics process ubiquitous in astrophysics, and important in both magnetic confinement fusion and space weather. The MARZ fundamental science program was recently established on Z to enable the first laboratory astrophysics platform able to access and study the strongly radiatively cooled magnetic reconnection regime. Simulations of this system have successfully used a resistive-MHD approach, but in some regions of parameter space Hall physics has the potential to be important. We describe implementation of a Hall method on a staggered grid resistive-MHD method (compatible with the approach used to model MARZ experiments. We then present a different Hall method based on cell-centered field quantities. Both approaches have been implemented in the Sandia KRAKEN code, to enable us to contrast different numerical Hall-MHD methods within the same HED code.

This document analyzes the role of Generative Artificial Intelligence (GenAI) tools in cybersecurity, particularly for red teaming. While GenAI accelerates initial security assessments, its effectiveness wanes with complexity, necessitating experienced assessors. The review critiques marketing claims, highlights ethical concerns regarding uncensored models for cybercrime, and advocates for a robust defense strategy supported by skilled professionals.

Herein, we report on the ultrafast photodissociation of nickel tetracarbonyl─a prototypical metal-ligand model system─at 197 nm. Using mid-infrared transient absorption spectroscopy to probe the bound C≡O stretching modes, we find evidence for the picosecond time scale production of highly vibronically excited nickel dicarbonyl and nickel monocarbonyl, in marked contrast with a prior investigation at 193 nm. Further spectral evolution with a 50 ps time constant suggests an additional dissociation step; the absence of any corresponding growth in signal strongly indicates the production of bare Ni, a heretofore unreported product from single-photon excitation of nickel tetracarbonyl. Thus, by probing the deep UV-induced photodynamics of a prototypical metal carbonyl, this Letter adds time-resolved spectroscopic signatures of these dynamics to the sparse literature at high excitation energies.

Photoinitiated polymerization enables spatiotemporal control of reaction conditions and can thereby generate materials with high complexity while consuming minimal energy. Where ring opening metathesis polymerization (ROMP) is concerned, photo-activated processes are typically enabled by chemical inhibition of ruthenium carbenes via the careful design of complexed ligands such that photoactivation can proceed through an isomerization or ligand dissociation event. In this contribution, we have explored a new approach to photoinitiation of ROMP based on physical inhibition through microencapsulation and controlled release of metathesis catalysts. Micron-sized particles of poly(phthalaldehyde) (PPA), catalyst, and photoacid generator were fabricated by spray drying. The particles were dispersed in dicyclopentadiene monomer, after which polymerization was initiated through temperature or UV exposure, both inducing depolymerization of the PPA particles and in situ catalyst release. The monomer/particle dispersions were found to be stable and reproducibly polymerizable with 3 weeks of storage at room temperature. Furthermore, the dispersions can be used for both photo- and thermal-initiated frontal ROMP, yielding a polymerized thermoset of equivalent properties to conventional bulk- and frontally-polymerized analogues. In conclusion, this work will ultimately enable new manufacturing techniques for ROMP-based materials, due to the modular, easily tunable nature of the underlying initiating system and its unparalleled stability.

High-entropy ceramics have garnered interest due to their remarkable hardness, compressive strength, thermal stability, and fracture toughness; yet the discovery of new high-entropy ceramics (out of a tremendous number of possible elemental permutations) still largely requires costly, inefficient, trial-and-error experimental and computational approaches. The entropy forming ability (EFA) factor was recently proposed as a computational descriptor that positively correlates with the likelihood that a 5-metal high-entropy carbide (HECs) will form the desired single phase, homogeneous solid solution; however, discovery of new compositions is computationally expensive. If you consider 8 candidate metals, the HEC EFA approach uses 49 optimizations for each of the 56 unique 5-metal carbides, requiring a total of 2744 costly density functional theory calculations. Here, we describe an orders-of-magnitude more efficient active learning (AL) approach for identifying novel HECs. To begin, we compared numerous methods for generating composition-based feature vectors (e.g., magpie and mat2vec), deployed an ensemble of machine learning (ML) models to generate an average and distribution of predictions, and then utilized the distribution as an uncertainty. We then deployed an AL approach to extract new training data points where the ensemble of ML models predicted a high EFA value or was uncertain of the prediction. Our approach has the combined benefit of decreasing the amount of training data required to reach acceptable prediction qualities and biases the predictions toward identifying HECs with the desired high EFA values, which are tentatively correlated with the formation of single phase HECs. Using this approach, we increased the number of 5-metal carbides screened from 56 to 15,504, revealing 4 compositions with record-high EFA values that were previously unreported in the literature. Our AL framework is also generalizable and could be modified to rationally predict optimized candidate materials/combinations with a wide range of desired properties (e.g., mechanical stability, thermal conductivity).

Recently, the spectral localizer framework has emerged as an efficient approach for classifying topology in photonic systems featuring local nonlinearities and radiative environments. In nonlinear systems, this framework provides rigorous definitions for concepts such as topological solitons and topological dynamics, where a system’s occupation induces a local change in its topology due to nonlinearity. For systems embedded in radiative environments that do not possess a shared bulk spectral gap, this framework enables the identification of local topology and shows that local topological protection is preserved despite the lack of a common gap. However, as the spectral localizer framework is rooted in the mathematics of C*-algebras, and not vector bundles, understanding and using this framework requires developing intuition for a somewhat different set of underlying concepts than those that appear in traditional approaches for classifying material topology. In this tutorial, we introduce the spectral localizer framework from a ground-up perspective and provide physically motivated arguments for understanding its local topological markers and associated local measure of topological protection. In doing so, we provide numerous examples of the framework’s application to a variety of topological classes, including crystalline and higher-order topology. We then show how Maxwell’s equations can be reformulated to be compatible with the spectral localizer framework, including the possibility of radiative boundary conditions. To aid in this introduction, we also provide a physics-oriented introduction to multi-operator pseudospectral methods and numerical K-theory, two mathematical concepts that form the foundation for the spectral localizer framework. Finally, we provide some mathematically oriented comments on the C*-algebraic origins of this framework, including a discussion of real C*-algebras and graded C*-algebras that are necessary for incorporating physical symmetries. Looking forward, we hope that this tutorial will serve as an approachable starting point for learning the foundations of the spectral localizer framework.

Carbon dots have attracted widespread interest for sensing applications based on their low cost, ease of synthesis, and robust optical properties. We investigate structure-function evolution on multiemitter fluorescence patterns for model carbon-nitride dots (CNDs) and their implications on trace-level sensing. Hydrothermally synthesized CNDs with different reaction times were used to determine how specific functionalities and their corresponding fluorescence signatures respond upon the addition of trace-level analytes. Archetype explosives molecules were chosen as a testbed due to similarities in substituent groups or inductive properties (i.e., electron withdrawing), and solution-based assays were performed using ratiometric fluorescence excitation-emission mapping (EEM). Analyte-specific quenching and enhancement responses were observed in EEM landscapes that varied with the CND reaction time. We then used self-organizing map models to examine EEM feature clustering with specific analytes. The results reveal that interactions between carbon-nitride frameworks and molecular-like species dictate response characteristics that may be harnessed to tailor sensor development for specific applications.

There is growing interest in material candidates with properties that can be engineered beyond traditional design limits. Compositionally complex oxides (CCO), often called high entropy oxides, are excellent candidates, wherein a lattice site shares more than four cations, forming single-phase solid solutions with unique properties. However, the nature of compositional complexity in dictating properties remains unclear, with characteristics that are difficult to calculate from first principles. Here, compositional complexity is demonstrated as a tunable parameter in a spin-transition oxide semiconductor La1− x(Nd, Sm, Gd, Y)x/4CoO3, by varying the population x of rare earth cations over 0.00≤ x≤ 0.80. Across the series, increasing complexity is revealed to systematically improve crystallinity, increase the amount of electron versus hole carriers, and tune the spin transition temperature and on-off ratio. At high a population (x = 0.8), Seebeck measurements indicate a crossover from hole-majority to electron-majority conduction without the introduction of conventional electron donors, and tunable complexity is proposed as new method to dope semiconductors. First principles calculations combined with angle resolved photoemission reveal an unconventional doping mechanism of lattice distortions leading to asymmetric hole localization over electrons. Thus, tunable complexity is demonstrated as a facile knob to improve crystallinity, tune electronic transitions, and to dope semiconductors beyond traditional means.

Harmonic and subharmonic RF injection locking is demonstrated in a terahertz (THz) quantum-cascade vertical-external-cavity surface-emitting laser (QC-VECSEL). By tuning the RF injection frequency around integer multiples and submultiples of the cavity round-trip frequency, different harmonic and subharmonic orders can be excited in the same device. Modulation-dependent behavior of the device has been studied with recorded lasing spectral broadening and locking bandwidths in each case. In particular, harmonic injection locking results in the observation of harmonic spectra with bandwidths over 200 GHz. A semiclassical Maxwell-density matrix formalism has been applied to interpret QC-VECSEL dynamics, which aligns well with experimental observations.

Here, this study investigates neutron-induced displacement damage in Bipolar Junction Transistors (BJTs) using TCAD models informed by Deep-Level-Transient-Spectroscopy (DLTS) data. These models are calibrated and validated against experimental measurements performed at various neutron fluences. Both npn and pnp transistor configurations are studied to analyze the effects of individual traps on carrier recombination and base leakage currents. In npn transistors, deep traps (0.42 eV from the conduction band) dominate at low voltages, while shallow traps (0.17 eV from the conduction band) become prominent at higher voltages. Conversely, pnp transistors have base leakage current predominantly due to deep-level traps. The study observes a notable trend in trap density versus fluence, characterized by a linear relationship on a log-log scale. These insights into defect evolution under radiation conditions are crucial for optimizing semiconductor device reliability and performance in radiation-prone environments.

The tunability of metal-organic frameworks (MOFs) makes them exceptional materials for the development of highly selective, low-power sensors for toxic gas detection. Herein, we demonstrate enhanced detection of NO2 gas by a MOF-based electrical impedance sensor made using a unique mixed metal MOF-on-MOF synthesis. A combined experimental and computational study was performed using the exemplar NixMg1-x-MOF-74 to understand the fundamental structure-property relationships behind metal mixing and MOF film synthesis methods on sensor performance. Density functional theory results indicated that the presence of Ni in Mg-MOF-74 increased framework stability and increased the electron density of states at lower energies near the HOMO, as well as enhanced the NO2-Mg adsorption interaction. Impedance data of the NixMg1-x-MOF-74 films with larger Ni contents showed greater impedance change after exposure to 1 ppm of NO2 gas. Furthermore, when synthesized through either a drop-cast or direct solvothermal film growth approach, the monometallic Ni-based sensors had the best performance. However, the mixed metal NixMg1-x-MOF-74 sensors synthesized through a MOF-on-MOF approach resulted in the highest impedance change, outperforming all monometallic Ni-based sensors. In particular, the mixed metal Ni-on-Mg-MOF-74 film was the best-performing sensor with an impedance change of 309 upon trace NO2 exposure. Change in impedance response after NO2 exposure was improved by 52% compared to the best monometallic Ni-on-Ni-MOF-74 sensor. Structural analysis of the Ni-on-Mg film showed that the first Mg-MOF-74 layer acts as a structural template controlling the structural features of the final film after metal exchange with Ni. This led to improved film quality, evidenced by the greater crystallinity and larger MOF grain sizes, and resulted in enhanced sensor performance which was not achievable through other metal mixing methods. Altogether, this study identifies structure-property relationships and synthetic templating methods that inform MOF-based sensor design, allowing for improved detection of toxic compounds.

This report details operations and maintenance (O&M) activities performed across Fiscal Years 2022 through 2024 in support of the continued capabilities of the National Solar Thermal Testing Facility (NSTTF) at Sandia National Laboratories. The NSTTF O&M project is funded by the U.S. Department of Energy Solar Energy Technologies Office (SETO) to support research activities and testing on behalf of external customers at the facility under award number CPS 38491. During the project period, the NSTTF made progress in the areas of site metrics, site maintenance and utilization tracking, and customer engagement. The O&M project also supported special initiatives including procurement of a heat exchanger for particle concentrating solar thermal processes and a scoping and cost study for refurbishment and repair of component in the NSTTF heliostat field.

Conceptual models of smectite hydration include planar (flat) clay layers that undergo stepwise expansion as successive monolayers of water molecules fill the interlayer regions. However, X-ray diffraction (XRD) studies indicate the presence of interstratified hydration states, suggesting non-uniform interlayer hydration in smectites. Additionally, recent theoretical studies have shown that clay layers can adopt bent configurations over nanometer-scale lateral dimensions with minimal effect on mechanical properties. Therefore, in this study we used molecular simulations to evaluate structural properties and water adsorption isotherms for montmorillonite models composed of bent clay layers in mixed hydration states. Results are compared with models consisting of planar clay layers with interstratified hydration states (e.g. 1W–2W). The small degree of bending in these models (up to 1.5 Å of vertical displacement over a 1.3 nm lateral dimension) had little or no effect on bond lengths and angle distributions within the clay layers. Except for models that included dry states, porosities and simulated water adsorption isotherms were nearly identical for bent or flat clay layers with the same averaged layer spacing. Similar agreement was seen with Na- and Ca-exchanged clays. In conclusion, while the small bent models did not retain their configurations during unconstrained molecular dynamics simulation with flexible clay layers, we show that bent structures are stable at much larger length scales by simulating a 41.6×7.1 nm² system that included dehydrated and hydrated regions in the same interlayer.

Theoretically and computationally describing the operation of nanodiodes requires characterizing the transitions between multiple electron emission mechanisms for nanodiodes with complicated geometries. This motivates our development of techniques to determine when simplified theories for individual mechanisms suffice compared to more complete, but more computationally expensive, models. Leveraging recent theories that define a canonical gap distance to translate planar theory to nonplanar diodes, we derive the conditions for the transitions among thermal emission, field emission, and space-charge-limited current density (SCLCD) in vacuum and with collisions for non-Cartesian coordinate systems, including spherical, cylindrical, and prolate spheroidal coordinate systems. Particle-in-cell (PIC) simulations of the current density as a function of applied voltage for a tip-to-plate geometry in vacuum agreed qualitatively with the asymptotes for thermal emission at low voltage and SCLCD at higher voltage using the canonical gap distance. As a result, this demonstrates the utility of this approach for guiding system design and suggests future extensions to save simulation time for more realistic geometries that are more computationally expensive.

Researchers are exploring adding wave energy converters to existing oceanographic buoys to provide a predictable source of renewable power. A ”pitch resonator” power take-off system has been developed that generates power using a geared flywheel system designed to match resonance with the pitching motion of the buoy. However, the novelty of the concept leaves researchers uncertain about various design aspects of the system. This work presents a novel design study of a pitch resonator to inform design decisions for an upcoming deployment of the system. The assessment uses control co-design via WecOptTool to optimize control trajectories for maximal electrical power production while varying five design parameters of the pitch resonator. Given the large search space of the problem, the control trajectories are optimized within a Monte Carlo analysis to identify optimal designs, followed by parameter sweeps around the optimum to identify trends between the design parameters. The gear ratio between the pitch resonator spring and flywheel are found to be the most sensitive design variables to power performance. The assessment also finds similar power generation for various sizes of resonator components, suggesting that correctly designing for optimal control trajectories at resonance is more critical to the design than component sizing.

The stochastic weighted particle method (SWPM) is a generalization of the Direct Simulation Monte Carlo (DSMC) method where particle weights are variable and dynamic. SWPM is backed by a strong theoretical foundation but has not been critically evaluated for problems of practical interest. A thorough assessment of SWPM for boundary-driven flows reveals significant numerical artifacts near the boundary, notably a diverging heat flux. To correct the boundary heat flux, two modifications to SWPM are proposed: separated grouping and a spatially-dependent weight transfer function. To gauge the relative efficiency of SWPM in comparison to DSMC, a high-Mach-number wheel flow which forms a strong density gradient is also simulated.

The shock response of fully-dense and porous crystalline tellurium dioxide (TeO 2 ) to the high-pressure and high-temperature fluid regime was investigated within the framework of density functional theory with Mermin’s generalization to finite temperatures. The principal and porous shock Hugoniot curves were predicted from canonical ab initio molecular dynamics (AIMD) simulations, with the phase space sampled along isotherms up to 80 000 K, for densities ranging from ρ = 3 to 17 g/cm 3 . The polymorphs investigated are α - Te O 2 paratellurite ( P 4 1 2 1 2 ), Te O 2 cotunnite ( P n m a ), and Te O 2 post-cotunnite ( P 2 1 / m ). Based on the discontinuity found in the calculated U s − u p slope of TeO 2 post-cotunnite at a shock velocity of U s ≃ 8.35 km/s and a particle velocity of u p ≃ 3.64 km/s, the shock melting temperature and pressure are predicted to be ≃ 6500 K and ≃ 170 GPa. Results from the AIMD simulations are in line with the static compression data of Te O 2 paratellurite and cotunnite, and with the recent shock Hugoniot data for single-crystal α - Te O 2 for pressures up to 85 GPa, obtained using the inclined-mirror method and the velocity interferometer system for any reflector combined with powder gun and two-stage light-gas gun.

Accurate material characterization and model calibration are pivotal for simulations used for high-consequence engineering decisions. Current characterization and calibration methods (1) use simplified test specimen geometries and global data, (2) cannot guarantee that sufficient characterization data is collected for a specific model of interest, (3) provide only mean parameter values with no uncertainty quantification, and (4) are sequential, inflexible, and time-consuming. This work developed a new paradigm—coined Interlaced Characterization and Calibration (ICC)—which drives forward the state-of-the-art in model calibration by bringing together recent advancements into one improved workflow. The ICC paradigm (1) employs tools to efficiently use full-field data to calibrate high-fidelity material models, (2) aligns the data needed with the data collected by adopting an optimal experimental design protocol, (3) provides uncertainty metrics on the calibrated model parameters, and (4) incorporates these advances into a quasi real-time feedback loop. The ICC framework was validated synthetically with both low-fidelity and high-fidelity simulations paired with several different elastoplastic material models, and was also demonstrated experimentally with an aluminum 6061 cruciform exemplar specimen. Results showed that the ICC framework—in which Bayesian optimal experimental design actively guided the experiment— resulted in calibrations with similar or better accuracy than predetermined experiments based on subject matter expertise. Moreover, the ICC framework produced a complete model calibration— with quantified uncertainties on model parameters—in 1 week, a 5 - 10× increase in efficiency over traditional approaches. Thus, the ICC paradigm improves both the calibration process and quality, by (1) improving efficiency, which increases agility of solid mechanics modeling and enables utilization of computational simulation (CompSim) at earlier stages of the design cycle and (2) providing quantified, and in some cases reduced, parameter uncertainties, which increases confidence in model predictions and supports credible decision making.

A neutron fluence map and a total ionizing dose map of the Los Alamos National Laboratory Godiva IV fast burst critical assembly was generated using passive reactor dosimetry, comprised of sulfur pellets and thermoluminescent dosimeters. Godiva IV is an unmoderated, fast burst, critical assembly constructed of approximately 65 kg of highly enriched uranium fuel alloyed with 1.5 % molybdenum for strength. [1] The mapping was performed during a single 75.6 ºC temperature rise burst operation, with the top and sides of the cylindrical Godiva-IV Top Hat covered in passive dosimetry. Dosimetry was placed in a symmetric pattern around the Top Hat, with higher concentrations near the control rods and burst rod. A specific portion of the lower quadrant of the burst rod was mapped to confirm a testing region where the neutron fluence varied by no more than ± 5%. The results will be used to assess the neutron, gamma, and total ionizing dose environment in three-dimensional space around the assembly for higher fidelity experiment placement, active dosimetry positioning, and radiation field characterization.

The characterization of the neutron, prompt gamma-ray, and delayed gamma-ray radiation fields for the White Sands Missile Range (WSMR) Fast Burst Reactor, also known as molybdenum-alloy Godiva (Molly-G) has been assessed at the 6-inch irradiation location. The neutron energy spectra, uncertainties, and common radiation metrics are presented. Code-dependent recommended constants are given to facilitate the conversion of various dosimetry readings into radiation metrics desired by experimenters. The Molly-G core was designed and configured similarly to Godiva II, as an unreflected, unmoderated, cylindrical annulus of uranium-molybdenum-alloy fuel with molybdenum loading of 10%. At the 6-inch position, the axial fluence maximum is about 2.4×1013 n/cm2 per MJ of reactor energy; about 0.1% of the neutron fluence is below 1 keV and 96% is above 100 keV. The 1-MeV Damage-Equivalent Silicon (DES) fluence is estimated at 2.2×1013 n/cm2 per MJ of reactor energy. The prompt gamma-ray dose is roughly 2.5E+03 rad(Si) per MJ and the delayed gamma-ray dose is about 1.3E+03 rad(Si) per MJ.

Fully-polarimetric synthetic aperture radar (PolSAR) data contain a rich body of elementary scattering physics information that is critically valuable for a broad range of applications and scientific purposes. However, there is a lack of available high-resolution (< 0.3048-m) data available for PolSAR phenomenology research. This article introduces a high-resolution PolSAR data set collected and provided by Sandia National Laboratories (SNL). The data sets were collected to support studying high-resolution scattering physics from different types of clutter and applications such as polarimetric-based terrain classification.

The characterization of the neutron, prompt gamma-ray, and delayed gamma-ray radiation fields in the University of Texas at Austin Nuclear Engineering Teaching Laboratory (NETL) TRIGA reactor for the beam port (BP) 1/5 free-field environment at the 128-inch location adjacent to the core centerline has been accomplished. NETL is being explored as an auxiliary neutron test facility for the Sandia National Laboratories radiation effects sciences research and development campaigns. The NETL reactor is a TRIGA Mark-II pulse and steady-state, above-ground pool-type reactor. NETL is intended as a university research reactor typically used to perform irradiation experiments for students and customers, radioisotope production, as well as a training reactor. Initial criticality of the NETL TRIGA reactor was achieved on March 12, 1992, making it one of the newest test reactor facilities in the US. The neutron energy spectra, uncertainties, and covariance matrices are presented as well as a neutron fluence map of the experiment area of the cavity. For an unmoderated condition, the neutron fluence at the center of BP 1/5, at the adjacent core axial centerline, is about 8.2×1012 n/cm2 per MJ of reactor energy. About 67% of the neutron fluence is below 1 keV and 22% above 100 keV. The 1-MeV Damage-Equivalent Silicon (DES) fluence is roughly 1.6×1012 n/cm2 per MJ of reactor energy.

The finite element method can be used to compute accurate stress intensity factors (SIFs) for cracks with complex geometries and boundary conditions. In contrast, handbook solutions act as surrogate SIF models that provide significantly faster evaluation times. However, the development of conventional surrogate SIF models relies on manual development based on low-order parameterizations. This limits surrogate model accuracy and generalizability. In this paper, we develop a framework for the automated development of mechanics-guided handbook SIF solutions by using interpretable machine learning via genetic programming for symbolic regression (GPSR). Formalizing the mechanics-based approach of Raju and Newman, SIF training data is decomposed into multiple subsets. This decomposition enables parallel GPSR model development of subfunctions, each of which accounts for specific geometrical corrections with respect to a known analytical model. Using this mechanics-based approach with GPSR allows for equations to be learned with improved accuracy and reduced complexity relative to the Raju Newman equations while maintaining the inherent interpretability of mathematical expressions. In this paper, we present equations that match the complexity of the Raju Newman equations while having reduced error, as well as equations with similar errors and reduced complexity.

High-entropy materials (HEMs) emerged as promising candidates for a diverse array of chemical transformations, including CO2 utilization. However, traditional HEMs catalysts are nonporous, limiting their activity to surface sites. Designing HEMs with intrinsic porosity can open the door toward enhanced reactivity while maintaining the many benefits of high configurational entropy. Here, a synergistic experimental, analytical, and theoretical approach to design the first high-entropy metal-organic frameworks (HEMOFs) derived from polynuclear metal clusters is implemented, a novel class of porous HEMs that is highly active for CO2 fixation under mild conditions and short reaction times, outperforming existing heterogeneous catalysts. HEMOFs with up to 15 distinct metals are synthesized (the highest number of metals ever incorporated into a single MOF) and, for the first time, homogenous metal mixing within individual clusters is directly observed via high-resolution scanning transmission electron microscopy. Importantly, density functional theory studies provide unprecedented insight into the electronic structures of HEMOFs, demonstrating that the density of states in heterometallic clusters is highly sensitive to metal composition. This work dramatically advances HEMOF materials design, paving the way for further exploration of HEMs and offers new avenues for the development of multifunctional materials with tailored properties for a wide range of applications.

The cells in battery energy storage systems are monitored, protected, and controlled by battery management systems whose sensors are susceptible to cyberattacks. False data injection attacks (FDIAs) targeting batteries’ voltage sensors affect cell protection functions and the estimation of critical battery states like the state of charge (SoC). Inaccurate SoC estimation could result in battery overcharging and over discharging, which can have disastrous consequences on grid operations. This paper proposes a three-pronged online and offline method to detect, identify, and classify FDIAs corrupting the voltage sensors of a battery stack. To accurately model the dynamics of the series-connected cells a single particle model is used and to estimate the SoC, the unscented Kalman filter is employed. FDIA detection, identification, and classification was accomplished using a tuned cumulative sum (CUSUM) algorithm, which was compared with a baseline method, the chi-squared error detector. Online simulations and offline batch simulations were performed to determine the effectiveness of the proposed approach. Throughout the batch simulations, the CUSUM algorithm detected attacks, with no false positives, in 99.83% of cases, identified the corrupted sensor in 97% of cases, and determined if the attack was positively or negatively biased in 97% of cases.

Vertically aligned nanocomposites (VANs) are unique thin films with vertical nanostructures embedded in a matrix material, allowing for the integration of two distinct materials. These nanocomposites offer novel combined physical properties, such as nanocomposite-based multiferroics and strongly coupled physical properties, such as magneto-optic coupling. Much work has been conducted in exploring different two-phase combinations and various processing conditions to achieve novel tunable properties that cannot be obtained by any single-phase material alone. In this work, the target configuration effects are explored for the growth of LaFeO3-CoNi2O4 VANs. Both mixed and pie-shaped targets are utilized to compare the target configuration effects on the phase separation, morphology tuning, and their resulting physical properties, including optical and magnetic properties. This work suggests that the target configuration is another important parameter for achieving the desired VAN morphology and can be used to design different two-phase VANs with tailorable properties.

This is the manual for the TeMPI Shim library, whose goal is to facilitate scientific apps being loosely coupled through MPI. TeMPI Shim can be used to create and manage Message Passing Interface (MPI, see [MPI]) Communicators for Single Program, Multiple Data (SPMD, see [SPMD]) and Multiple Programs, Multiple Data (MPMD, see [MPMD]) workflows. MPI provides the MPI_APPNUM command number to each of the processes within its context. This number, starting at 0, is the application number being executed. In a case where 3 applications are being run in MPMD mode, there would be MPI_APPNUM values of 0, 1, and 2. TeMPI Shim creates intra- and inter-communicators between each pair of MPI_APPNUM values. In the aforementioned case, application 0 would have intra- and inter-communicators to speak to itself, application 1, and application 2. This is replicated for each of the applications. Additionally, TeMPI Shim creates intra- and inter-communicators for the first MPI rank of each application to directly communicate only with each other. Finally, TeMPI Shim creates its own copy of the default world communicator, i.e., MPI_COMM_WORLD. In the case where there is only a single application, it will have the communicators to only speak with itself. TeMPI Shim is useful in this case since it is considered good practice for MPI applications to copy the default world communicator and reference this copy (see [Duplicate World]_) anyways. Ultimately, it provides value independent of the number of applications present.

Thermal spray processes can benefit from cooling to maintain substrate temper, reduce processing times, and manage thermally induced residual stresses. “Plume quenching” is a plume-targeted cooling technique which has been shown to reduce substrate temperatures by redirection of hot plume gases using a lateral argon curtain injected into the plume, while limiting interaction with the substrate or affecting coating properties. Here, this study explores the use of this technique for residual stress management by reducing the thermally driven component in nickel and tantalum coatings on titanium and aluminum substrates. The in-situ residual stress profiles were measured for all substrate and coating pairings during spraying and cooling, and the deposition and thermal stresses recorded. For substrate and coating pairings where the predominant component of residual stress was thermal (driven by a large difference in coefficient of thermal expansion, Δα, between coating and substrate), plume quenching reduced both the thermal stress and the final stress state of the coating. This was seen primarily in tantalum on aluminum coatings where the Δα was -17 × 10^-6 /°C, and thermal stress was reduced by 7.5% and 22.4% for the plume quenching rates of 50 and 100 slpm, respectively.

Thermal spray processes can benefit from cooling to maintain substrate temper, reduce processing times, and manage thermally induced residual stresses. “Plume quenching” is a plume-targeted cooling technique which has been shown to reduce substrate temperatures by redirection of hot plume gases using a lateral argon curtain injected into the plume, while limiting interaction with the substrate or affecting coating properties. Here, this study explores the use of this technique for residual stress management by reducing the thermally driven component in nickel and tantalum coatings on titanium and aluminum substrates. The in-situ residual stress profiles were measured for all substrate and coating pairings during spraying and cooling, and the deposition and thermal stresses recorded. For substrate and coating pairings where the predominant component of residual stress was thermal (driven by a large difference in coefficient of thermal expansion, Δα, between coating and substrate), plume quenching reduced both the thermal stress and the final stress state of the coating. This was seen primarily in tantalum on aluminum coatings where the Δα was -17 × 10^-6 /°C, and thermal stress was reduced by 7.5% and 22.4% for the plume quenching rates of 50 and 100 slpm, respectively.

Physical Review Letters

The H3D was characterized using GADRAS 19.3.5 using data sent by the IAEA. Everything related to the measurement setup was provided by the IAEA. The sources measured were detailed and the distance from the detector was also provided. The data sheets of the sources were used to determine the source activity at the time of calibration, as needed for characterization.

Analyzing X-ray radiographs is crucial for understanding target behavior in Inertial Confinement Fusion (ICF) and High Energy Density (HED) platforms. However, the density of Magneto Raleigh Taylor (MRT) bands and limitations of target materials often obscure relevant spike growth and density information. To address this issue, machine learning postprocessing techniques can be applied to remove darkened regions in radiography images. In this study, a novel method is presented for removing MRT darkened regions from z-target radiographs using a convolutional neural network (CNN). The CNN, consisting of six layers, treats the darkened regions as noise and employs a mixed loss function and end-to-end frameworks to suppress them while preserving sharpness. The six-layer architecture is designed to effectively learn features when provided with a larger volume of learning space. Each layer is optimized using a mixed loss function that combines a standard loss pixel approach with a multi-scaled structural similarity index loss, which considers luminance, contrast, and structure in local neighborhoods. This approach is particularly beneficial for capturing the stochastic structure of MRT limbs. Due to the limited availability of experimental data, training is conducted using synthetic target radiography from 3D Alegra simulations.

Salt formations are one of the potential host rocks for the final disposal of high-level radioactive waste (HLW) in deep geological repositories, both in Germany and the United States. The safe isolation of radioactive waste in these repositories relies on a multi-barrier system, combining engineered and natural barriers. The natural barrier is provided by the salt rock itself, known for its self-sealing properties and long-term stability. The engineered barrier, on the other hand, comprises sealing components strategically placed within the repository to enhance its containment capabilities. In both Germany and the United States, long-term safety assessments require demonstrating the integrity of the natural barrier for a period of up to 1 million years. Concurrently, the engineered barrier system (EBS) must maintain its structural and functional integrity until the long-term sealing, such as the granular salt backfill material, has re-consolidated to its final low porosity and permeability. Based on extensive expertise and experience with engineered barriers in salt formations, BGE TECHNOLOGY GmbH and Sandia National Laboratories have partnered to develop a robust methodology for the integrity and performance assessment of EBS in HLW repositories through the RANGERS project. This collaborative effort aims to establish a unified approach to geotechnical engineering, repository design, integrity and performance evaluation of EBS in salt repositories.

Abstract not provided.

This report investigates the integration of infrasound and seismic data to improve event localization accuracy, specifically focusing on a surface explosion at the Utah Training and Testing Range (UTTR). Utilizing the Seismoacoustic Bayesian Event Locator (SABEL) framework, we incorporated atmospheric specifications derived from Ground to Space (G2S) profiles to enhance celerity-range priors. Our analysis revealed that while the combination of infrasound and seismic observations significantly reduced localization uncertainty, challenges remained, particularly with returns at distances less than 200 km from the source and the influence of specific observations on location estimates. The results indicate that broader celerity distributions, such as those from Blom et al. (2020), facilitate better alignment with ground truth locations compared to narrower models. Overall, this work demonstrates the promise of seismoacoustic approaches in refining event localization and highlights the need for further exploration of celerity-range models to ensure reliable outcomes.

Abstract not provided.

Abstract not provided.

The use of differential scanning calorimetry (DSC) to measure the thermal behavior of individual components and electrolyte/electrode combinations is common. However, here we focus on DSC tests on an anode, cathode, and electrolyte (ACE) component combination over a temperature range that includes many of the phase transitions and key reactions (i.e., to 500 °C) that contribute to thermal runaway. This method can help quantify the complex reaction network in a full cell, thereby informing potential safety issues. Here, we used DSC heat flow data from a solid-state Li0.43CoO2+C+PVDF | LLZO | Li metal ACE sample and its components to quantify key factors affecting results. We focused on three areas: (1) ACE sample preparation and assembly in DSC pans, (2) DSC measurement parameters, and (3) heat flow analysis. Key points include the choice of component ratios (e.g., commercially relevant N:P capacity ratio), the importance of conductive carbon and binder, type of pan used, DSC ramp rate, and integration method used when dealing with broad and overlapping exothermic peaks. This work deepens the scientific basis and best practices for obtaining heat flow data from ACE samples for early-stage evaluation of solid-state and beyond battery safety.

In this report, we survey the current intersection between the fields of radar technology and artificial intelligence. Three main areas are highlighted - synthetic aperture radar automatic target detection, waveform optimization, and antenna design. Literature relevant to these applications and beyond are discussed and compiled in an annotated bibliography.

A new tool has been developed to perform variance-based global sensitivity analysis (VBGSA) on data from a set-based concurrent engineering software called Success Assured (SA). The tool is part of a digital component design app, which is currently in production as an Accelerated Digital Engineering Pathfinder at Sandia National Laboratories. When working with complex digital models, it is important to understand relationships between inputs and outputs, i.e., how “sensitive” model outputs are to changes in model inputs. After extensive research and trials of various sensitivity analysis methods, it was determined that estimation of Sobol’ indices for VBGSA with Monte Carlo simulation, paired with simple surrogate models, produces the best results for SA data. This tool increases understanding of SA models and streamlines the creation of SA datasets. This report details the methodology and implementation of this sensitivity analysis tool so others can understand it and implement it.

The Blue Canyon Dome (BCD) experiments are a series of intermediate scale subsurface chemical explosions designed to explore the geophysical relationship between scaled yields and formation damage. Phase 2 of this experiment involved explosions designed to evaluate the effect of damage, demonstrating the effect of geologic features and previous damage on the near field seismic signatures, not just in amplitude, but also in spectral content. Based on the emplacement effect on the seismic signature, the moment-based yield estimates showed error between 10s-1000s of percent. This report details an alternative approach to emplacement-based filtering techniques of the seismic signatures to improve yield estimates for explosions in damaged rock. The results will be used to develop an improved inverse yield model that leverages the information in the seismic signature for each source based on emplacement type.

The Human Readiness Level (HRL) scale is a simple nine-level scale that brings structure and consistency to the real-world application of user-centered design. It enables multidisciplinary consideration of human-focused elements during the system development process. Use of the standardized set of questions comprising the HRL scale results in a single human readiness number that communicates system readiness for human use. The Human Views (HVs) are part of an architecture framework that provides a repository for human-focused system information that can be used during system development to support the evaluation of HRL levels. This paper illustrates how HRLs and HVs can be used in combination to support user-centered design processes. A real-world example for a U.S. Army software modernization program is described to demonstrate application of HRLs and HVs in the context of user-centered design.

This paper presents a probabilistic model for various machine learning (ML) applications. While deep learning (DL) has produced state-of-the-art results in many domains, DL models are complex and over-parameterized, which leads to high uncertainty about what the model has learned, as well as its decision process. Further, DL models are not probabilistic, making reasoning about their output challenging. In contrast, the proposed model, referred to as Yet Another Discriminate Analysis(YADA), is less complex than other methods, is based on a mathematically rigorous foundation, and can be utilized for a wide variety of ML tasks including classification, explainability, and uncertainty quantification. YADA is thus competitive in most cases with many state-of-the-art DL models. Ideally, a probabilistic model would represent the full joint probability distribution of its features, but doing so is often computationally expensive and intractable. Hence, many probabilistic models assume that the features are either normally distributed, mutually independent, or both, which can severely limit their performance. YADA is an intermediate model that (1) captures the marginal distributions of each variable and the pairwise correlations between variables and (2) explicitly maps features to the space of multivariate Gaussian variables. Numerous mathematical properties of the YADA model can be derived, thereby improving the theoretic underpinnings of ML. Validation of the model can be statistically verified on new or held-out data using native properties of YADA. However, there are some engineering and practical challenges that we enumerate to make YADA more useful.

In this paper, we present a Riemannian geometric derivation of the governing equations of motion of nonholonomic dynamic systems. A geometric form of the work-energy principle is first derived. The geometric form can be realized in appropriate generalized quantities, and the independent equations of motion can be obtained if the subspace of generalized speeds allowable by nonholonomic constraints can be determined. We provide a geometric perspective of the governing equations of motion and demonstrate its effectiveness in studying dynamic systems subjected to nonholonomic constraints.

In recent years, Native tribes across the country have embraced solar energy to power their communities. However, in approximately 25-30 years after installation, the solar panels will enter end-of-life (EOL) status, meaning that they will begin to operate less efficiently than upon initial implementation. At this point, tribes can determine the outcome: they can continue to operate the array, or they can dispose of the panels. Recycling is an option that promotes the use of high-value recyclable materials in solar panels in the circular economy. There are several companies across the US that offer recycling of solar panels, but at various costs. It is important that a hazardous waste expert makes an initial assessment of the solar array to identify any solar panels that cannot be recycled. Currently, recycling solar panels is not cost-effective or environmentally sustainable, but significant research into improving the cost and carbon footprint is ongoing.

Blending hydrogen into existing natural gas pipelines is a current initiative for large-scale hydrogen deployment in the United States. Understanding changes in hazards and risk for fuel blends can help inform regulations. A differential risk assessment was conducted to understand how the safety profile of fires and unignited flammable plumes of leaking fuel would change based on blend composition. The area affected by harmful heat fluxes from jet fires was found to be larger for blends with more natural gas in them, while the distance affected by unignited flammable fuel concentrations was found to be larger for blends with more hydrogen in them. Component leak frequencies and ignition probabilities greatly impact both heat flux risk and the flammable plume risk. The collection of additional blend-specific industry data could help understand which assumptions are most valid to make when calculating risk.

EPJ Web of Conferences

Next-generation space remote sensing systems may be equipped with imaging arrays that sense data at a rate that outstrips the processing capability of any computing hardware that can operate within a satellite’s power budget. This project developed novel convolutional and recurrent neural networks to detect and estimate point-like events amid clutter, and investigated their efficient and accurate implementation on analog in-memory computing systems that are 10-1000× more energy-efficient than digital processors. This project leveraged two memory devices at different levels of technological maturity: a large-scale analog computing prototype using commercial SONOS charge-trap memory, and electrochemical memory (ECRAM) with intrinsic radiation hardness. We experimentally demonstrated end-to-end analog processing of our neural networks on SONOS and characterized the radiation response of both SONOS and ECRAM. We advanced the state-of-the-art in ECRAM precision and reliability, and developed co-design methods to enable accurate long-term operation of SONOS analog accelerators in space radiation environments.

Abstract not provided.

Hydrogen-powered locomotives may present an alternative to diesel for achieving transportation-related decarbonization, energy security, and resilience goals. As an emerging technology, hydrogen locomotives may not have adequate representation in current safety codes and standards. NFPA 2 and other appropriate codes and standards were reviewed to identify technical gaps in the code for hydrogen rail defueling, refueling, maintenance, and storage. Several technical gaps pertaining to setback distances, allowable quantities, ventilation rates, and grounding have been identified. Fueling requirements can be informed by road vehicle standards but may need to be revised for larger quantity hydrogen locomotive systems. Since there are no hydrogen-specific locomotive design and safety standards, some aspects of diesel-based locomotive standards may apply. CGA G- 5.5 currently provides the most guidance on height, placement, orientation, and design of vent systems, while other standards emphasize the importance of shielding vent stacks and provide requirements for orientation of their discharge.

The valorization and dissolution of lignin using ionic liquids (ILs) is critical for developing sustainable biorefineries and a circular bioeconomy. This review aims to critically assess the current state of computational and machine learning methods for understanding and optimizing IL-based lignin dissolution and valorization processes reported since 2022. The paper examines various computational approaches, from quantum chemistry to machine learning, highlighting their strengths, limitations, and recent advances in predicting and optimizing lignin-IL interactions. Key themes include the challenges in accurately modeling lignin’s complex structure, the development of efficient screening methodologies for ionic liquids to enhance lignin dissolution and valorization processes, and the integration of machine learning with quantum calculations. These computational advances will drive progress in IL-based lignin valorization by providing deeper molecular-level insights and facilitating the rapid screening of novel IL-lignin systems.

The deliberate introduction of stacking faults has been found to induce properties markedly distinct from those of perfectly stacked particles or sheets, with even minute deviations leading to significant alterations in material characteristics. In this review article, we discuss the effects of stacking faults – both linear and rotational – on surface structures influencing ion adsorption, particle–particle affinity governing crystal growth, friction, and electronic properties. Our comprehensive overview sheds light on the multifaceted impact of stacking faults on these phenomena, offering insights that bridge multiple disciplines, and provides numerous interdisciplinary research questions, paving the way for future innovations.

The rapid expansion of photovoltaic (PV) systems, particularly inverters, has introduced new cybersecurity challenges that threaten both local operations as well as the broader electrical grid’s stability. PV inverters, integrated into critical energy infrastructure are potential targets for cyber attacks due to vulnerabilities in firmware, remote access systems, and communication protocols. The Coordinated Vulnerability Disclosure (CVD) process, as defined by the Cybersecurity and Infrastructure Security Agency (CISA), provides a framework for identifying, reporting, and addressing these vulnerabilities in a transparent and collaborative manner. This report outlines the CVD process as it applies to PV systems, detailing the roles of key stakeholders, such as manufacturers, grid operators, and security researchers. The report also highlights specific challenges in managing vulnerabilities for new and legacy PV systems, which includes those introduced by insecure communications and third-party supply chain components. By adhering to the CVD process, the PV industry can mitigate cybersecurity risks, ensure regulatory compliance, and maintain consumer trust, while safeguarding the operational resilience of the energy grid. Ultimately, the effective coordination of vulnerability management is crucial for securing the future of PV systems within the critical electric grid infrastructure landscape.

Working together, Sandia National Laboratories, Southern California Edison (SCE) - an Investor-Owned Utility (IOU) - and the California Public Utilities Commission (CPUC) are studying how electric utilities can use equity and resilience metrics to help inform the prioritization and sequencing of resilience-driven infrastructure investments. To this end, this project evaluated “Social Burden,” an equitable resilience metric which measures the potential impact of disruptions in access to non-electric critical services on people and estimates community resilience to these disruptions. The Social Burden was expanded to incorporate SCE’s existing equity metric and applied to evaluate the potential impacts from a range of climate-informed hypothetical outage scenarios developed under SCE’s 2022 Climate Adaptation Vulnerability Assessment. One baseline (“blue-sky”) state and eight different outage scenarios were evaluated to measure the potential impacts of the outages on non-electric infrastructure, critical services, and people. Key findings include: 1) the Social Burden framework is flexible enough to adapt to and build upon existing utility equity and/or resilience metrics, 2) Social Burden results highlight the high degree of non-electric service redundancy within the SCE service area with most (6/8) hypothetical outage scenarios predicted to increase people’s Social Burden by less than 10%; however, 3) access to critical services and people’s ability to obtain them is unequal and spatially clustered, meaning that there are some hypothetical outage scenarios (2/8) that will exert a higher toll on communities directly experiencing the outage as well as some nearby communities with pre-existing vulnerabilities. The report concludes with recommendations for potential use cases of the expanded Social Burden metric and identifies priority follow-on work. Potential use cases may include incorporating equity into IOU’s prioritization of climate resilience investments. Additionally, Social Burden analysis may provide additional data and insights to augment grid planning, potentially by identifying additional needs and/or prioritizing previously identified needs.

The accurate modeling of electromagnetic penetration is an important topic in computational electromagnetics. Electromagnetic penetration occurs through intentional or inadvertent openings in an otherwise closed electromagnetic scatterer, which prevent the contents from being fully shielded from external fields. To efficiently model electromagnetic penetration, aperture or slot models can be used with surface integral equations to solve Maxwell's equations. A necessary step towards establishing the credibility of these models is to assess the correctness of the implementation of the underlying numerical methods through code verification. Surface integral equations and slot models yield multiple interacting sources of numerical error and other challenges, which render traditional code-verification approaches ineffective. In this paper, we provide approaches to separately measure the numerical errors arising from these different error sources for the method-of-moments implementation of the electric-field integral equation with a slot model. We demonstrate the effectiveness of these approaches for a variety of cases.

The defect density present at the dielectric-semiconductor interface in an MOS structure directly influences the channel carrier characteristics in semiconductor devices, especially in wide bandgap material systems used in power devices. While these trap defects are typically quantified through electrical characterization of MOS-capacitor test structures, this treatment offers very little insight into the physical nature of interface defects. Such shortcomings demand a physical characterization strategy to guide fabrication optimization. X-ray photoelectron spectroscopy (XPS) is suggested as a viable technique to determine chemical data for dielectric interfaces formed using atomic layer deposition (ALD) on GaN substrates. Previously, 1-D XPS characterization has confirmed the presence of a Ga_xO_y interlayer between ALD dielectrics and the GaN substrate. In this work, XPS data is serially collected to form 2-D images of an ALD-Al₂O₃/GaN interface as a proof-of-concept experiment for in-situ XPS quality monitoring during ALD processing. The information provided by this work reveals some of the challenges for incorporating XPS characterization as an in-situ strategy during fabrication of GaN-based devices. Separately, electrical mapping of a 2-D array of ALD-Al₂O₃/GaN MOS-capacitor devices provide a means to quantify the spatial variations in interface quality across a single wafer. Physical characterization techniques, such as time-of-flight secondary ion mass spectroscopy, provide additional chemical information about the Al₂O₃/Ga_xO_y/GaN structure that complement the electrical mapping results. This analysis shows that a higher Ga_xO_y content correlates with higher interface state defects for trap energies deep in the band gap.

This document serves as a comprehensive final report for the Single-Volume Scatter Camera (SVSC) project. Most of the work over the course of this project has been documented in journal articles, conference papers, and other reports. We therefore reference those materials for detailed presentation of our technical results. The most recent efforts on the project have not been published and are presented in detail here. We present characterizations of two neutron scatter camera prototypes; one using a monolithic geometry, and one using an optically segmented geometry. Both detectors employ plastic scintillator with SiPM-based readout. For the monolithic prototype, we present calibrations of several detector parameters made using a set of dark counts. In particular, we employ a coincidence analysis to characterize the distribution external optical crosstalk, along with the relative timing of the different SiPMs. For the optically segmented prototype we present the results of calibrations made using a tagged ²²Na source, as well as the results of an imaging measurement made using a tagged AmBe source.

Abstract not provided.

The accurate estimation of seismic wave amplitude is vital to precisely determine the yield, magnitude, and event discrimination possible for a given network – a critical element in nuclear explosion monitoring. This task is complicated by several factors, including but not limited to radiation pattern, scattering effects, and crustal variations, which can lead to the attenuation or amplification of amplitude along a given raypath. In this report, we explore the novel application of machine learning to the task of seismic amplitude estimation by training a simple Artificial Neural Network (ANN) on an S-wave amplitude dataset from Lai et al. (2019). Attributes from this dataset used as input to the ANN included event-station distances, station locations (latitude, longitude), event locations (latitude, longitude), event depths, event magnitudes, radiation patterns, signal-to noise ratio (SNR) measurements (average-amplitude, peak-to-trough, maximum peak), and signal periods. We find that the trained ANN predicts S-wave amplitudes with a modest tendency toward underestimating the actual values, as indicated by a linear regression between predicted and actual data (slope: 0.892, intercept: -0.651). These results suggest that an ANN can perform this task, with potential for significant improvements through improved datasets, architectures, and parameter tuning.

Classical short-circuit programs that linearize the power network are no longer applicable for inverter based resources (IBRs), necessitating an iterative approach. Phasor domain programs can model the IBRs using an iterative approach considering nonlinear fault responses. In phasor domain models, the IBR can be modeled as a voltage controlled current source (VCCS) in tabular form with positive, negative, and zero sequence information for balanced and unbalance short-circuit faults. In the VCCS modeling of the IBR, positive and negative incremental reactive current, also known as the k-factor, plays an important role in short circuit program convergence. In this work, a few approaches: conventional VCCS modeling with a k-factor of 2, conventional VCCS modeling with a k-factor of 2 with modified pre-fault voltages, VCCS characteristics based on the power flow solution with a k-factor of 2, and VCCS characteristics based on the power flow solution with a k-factor of 1 are investigated for short circuit program convergence under higher IBR penetration. The IEEE 39 bus New England Test System is taken as the test system, and simulations are carried out in PSS®CAPE 15.0.26 simulation software. Simulation results demonstrate that IBR penetration is higher for the VCCS model, which corresponds to the power flow solution with k-factor 1, compared to other approaches.

Pozzolans rich in silica and alumina react with lime to form cementing compounds and are incorporated into portland cement as supplementary cementitious materials (SCMs). However, pozzolanic reactions progress slower than portland cement hydration, limiting their use in modern construction due to insufficient early-age strength. Hence, alternative SCMs that enable faster pozzolanic reactions are necessary including synthetic zeolites, which have high surface areas and compositional purity that indicate the possibility of rapid pozzolanic reactivity. Synthetic zeolites with varying cation composition (Na-zeolite, H-zeolite), SiO2/Al2O3 ratio, and framework type were evaluated for pozzolanic reactivity via Ca(OH)2 consumption using ion exchange and in-situ X-ray diffraction experiments. Na-zeolites exhibited limited exchange reactions with KOH and Ca(OH)2 due to the occupancy of acid sites by Na+ and hydroxyl groups. Meanwhile, H-zeolites readily adsorbed K+ and Ca2+ from a hydroxide solution by exchanging cations with H+ at Brønsted acid sites or cation adsorption at vacant acid sites. By adsorbing cations, the H-zeolite reduced the pH and increased Ca2+ solubility to promote pozzolanic reactions in a system where Ca(OH)2 dissolution/diffusion was a rate limiting factor. High H-zeolite reactivity resulted in 0.8 g of Ca(OH)2 consumed per 1 g of zeolites after 16 h of reaction versus 0.4 g of Ca(OH)2 consumed per 1 g of Na-zeolite. The H-zeolite modulated the pore fluid alkalinity and created a low-density amorphous silicate phase via mechanisms analogous to two-step C-S-H nucleation experiments. Controlling these reaction mechanisms is key to developing next generation pozzolanic cementitious systems with comparable hydration rates to portland cement.

This report describes computational modeling in the HyRAM+ software for study of hydrogen behavior in two common scenarios. First, models of unignited plumes exiting a vent stack were considered. It was shown that entrainment and vent backpressure were major factors in plume physics. Second, HyRAM+ was extended to include wind effects on plumes by modifying the plume momentum balance and entrainment modeling. Use of the model showed that plume shape and length changed with wind speed and direction; in all cases, wind causes a shortening of the plume along the streamline. While the no-wind case in HyRAM+ has been validated and the newly developed wind model was fitted to very limited experimental data, more controlled experimental configurations would help validate the models and ensure accurate simulation of hydrogen plume behavior for vent stack releases or in wind.

Machine-learning function representations such as neural networks have proven to be excellent constructs for constitutive modeling due to their flexibility to represent highly nonlinear data and their ability to incorporate constitutive constraints, which also allows them to generalize well to unseen data. In this work, we extend a polyconvex hyperelastic neural network framework to (isotropic) thermo-hyperelasticity by specifying the thermodynamic and material theoretic requirements for an expansion of the Helmholtz free energy expressed in terms of deformation invariants and temperature. Different formulations which a priori ensure polyconvexity with respect to deformation and concavity with respect to temperature are proposed and discussed. The physics-augmented neural networks are furthermore calibrated with a recently proposed sparsification algorithm that not only aims to fit the training data but also penalizes the number of active parameters, which prevents overfitting in the low data regime and promotes generalization. The performance of the proposed framework is demonstrated on synthetic data, which illustrate the expected thermomechanical phenomena, and existing temperature-dependent uniaxial tension and tension-torsion experimental datasets.

Sandia National Labs will provide Technical Assistance and consulting to Requester on a series of Requester’s test structures and devices to characterize and aid in the determination of optimal fabrication parameters via fabrication of functional test structures and cells. Working with previously established processes and test die modules, Sandia will construct a series of test chips with varying materials properties to allow the electrical characterization and component testing of Requester’s test structures and devices. Material variations will be selected for optimization of yield and desired target operation parameters. A deliverable in the form of a report will be written. Sandia intends to deliver testable die to the Requester.

Abstract not provided.

The Engineered Barrier System (EBS) plays an important role in ensuring the long-term safety and containment of high-level waste (HLW) and spent nuclear fuel (SNF) in deep geological repositories in salt formation. As part of a multi-barrier system, the EBS works alongside the natural barrier, which is the salt formation itself and the technical barrier comprising the disposal casks. The primary function of the EBS is to maintain containment during a defined period until the backfill used in the repository made of crushed salt, develops its sealing capacity through compaction. Over the time, the backfill eventually compacts to a state of low porosity and permeability, acting as a long-term seal. However, until this process is complete, the EBS must retain its structural and functional integrity. Regulatory guidelines in Germany currently require the EBS to remain effective for up to next ice age, that is expected in 50,000 years. The significant hydro-geological and topographic changes expected during an ice age could make it impossible to accurately predict the hydro-chemical conditions within the repository system at that time. In response to these challenges, BGE TECHNOLOGY GmbH (BGE TEC) and Sandia National Laboratories (SNL) have jointly developed a comprehensive methodology for the design and safety assessment of engineered barrier systems within the scope of the RANGERS project. This methodology is tailored for repositories in salt formations. The developed methodology provides a structured approach for designing and assessing the performance of the EBS in salt-based repositories. It begins with defining a sealing concept based on the geological characteristics of the selected site and the overall repository design. The entire repository system, comprising the geological site, repository infrastructure, and EBS, is then subjected to a Features, Events, and Processes (FEP) analysis, focusing solely on those FEPs that affect the EBS. The derived FEPs help identify the loads and stresses acting on the EBS, which serve as the foundation for conducting an integrity assessment. This analysis helps predict the EBS’s evolution and performance over the regulatory time frame, feeding into integrated performance assessment simulations.

Biaxial stress is identified to play an important role in the polar orthorhombic phase stability in hafnium oxide-based ferroelectric thin films. However, the stress state during various stages of wake-up has not yet been quantified. In this work, the stress evolution with field cycling in hafnium zirconium oxide capacitors is evaluated. The remanent polarization of a 20 nm thick hafnium zirconium oxide thin film increases from 9.80 to 15.0 µC cm−2 following 106 field cycles. This increase in remanent polarization is accompanied by a decrease in relative permittivity that indicates that a phase transformation has occurred. The presence of a phase transformation is supported by nano-Fourier transform infrared spectroscopy measurements and scanning transmission electron microscopy that show an increase in ferroelectric phase content following wake-up. The stress of individual devices field cycled between pristine and 106 cycles is quantified using the sin2(ψ) technique, and the biaxial stress is observed to decrease from 4.3 ± 0.2 to 3.2 ± 0.3 GPa. The decrease in stress is attributed, in part, to a phase transformation from the antipolar Pbca phase to the ferroelectric Pca21 phase. This work provides new insight into the mechanisms controlling and/or accompanying polarization wake-up in hafnium oxide-based ferroelectrics.

Hydrogen fuel sources offer alternatives to conventional fuels in the rail transportation industry. Hydrogen powered locomotive designs utilizing either a fuel cell or an internal combustion engine can make migration to alternative fuels possible for rail transportation. Codes and standards are still in development for rail application of hydrogen and safety risks must be assessed for hydrogen locomotive applications. This report utilizes a failure mode & effects analysis framework to help qualitatively understand possible risks from a hydrogen locomotive system. Findings illustrate how a combination of three mitigations greatly reduces the risks from a hydrogen locomotive system.

Distributed Acoustic Sensing (DAS) can record acoustic wavefields at high sampling rates and with dense spatial resolution difficult to achieve with seismometers. Using optical scattering induced by cable deformation, DAS can record strain fields with spatial resolution of a few meters. However, many experiments utilizing DAS have relied on unused, dark telecommunication fibers. As a result, the geophysical community has not fully explored DAS survey parameters to characterize the ideal array design. This limits our understanding of guiding principles in array design to deploy DAS effectively and efficiently in the field. A better quantitative understanding of DAS array behavior can improve the quality of the data recorded by guiding the DAS array design. Here we use steered response functions, which account for DAS fiber’s directional sensitivity, as well as beamforming and back-projection results from forward modelling calculations to assess the performance of varying DAS array geometries to record regional and local sources. A regular heptagon DAS array demonstrated improved capabilities for recording regional sources over other polygonal arrays, with potential improvements in recording and locating local sources. These results help reveal DAS array performance as a function of geometry and can guide future DAS deployments.

Uncertainty Quantification (UQ) is vital to safety-critical model-based analyses, but the widespread adoption of sophisticated UQ methods is limited by technical complexity. In this paper, we introduce UM-Bridge (the UQ and Modeling Bridge), a high-level abstraction and software protocol that facilitates universal interoperability of UQ software with simulation codes. It breaks down the technical complexity of advanced UQ applications and enables separation of concerns between experts. UM-Bridge democratizes UQ by allowing effective interdisciplinary collaboration, accelerating the development of advanced UQ methods, and making it easy to perform UQ analyses from prototype to High Performance Computing (HPC) scale. In addition, we present a library of ready-to-run UQ benchmark problems, all easily accessible through UM-Bridge. These benchmarks support UQ methodology research, enabling reproducible performance comparisons. We demonstrate UM-Bridge with several scientific applications, harnessing HPC resources even using UQ codes not designed with HPC support.

Welding processes used in the production of pressure vessels impart residual stresses in the manufactured component. Computational modeling is critical to predicting these residual stress fields and understanding how they interact with notches and flaws to impact pressure vessel durability. Here, in this work, we present a finite element model for a resistance forge weld and validate it using laboratory measurements. Extensive microstructural changes, near-melt temperatures, and large localized deformations along the weld interface pose significant challenges to Lagrangian finite element modeling. The proposed modeling approach overcomes these roadblocks in order to provide a high-fidelity simulation that can predict the residual stress state in the manufactured pressure vessel; a rich microstructural constitutive model accounts for material recrystallization dynamics, a frictional-to-tied contact model is coordinated with the constitutive model to represent interfacial bonding, and adaptive remeshing is employed to alleviate severe mesh distortion. An interrupted-weld approach is applied to the simulation to facilitate comparison to displacement measures. Several techniques are employed for residual stress measurement in order to validate the finite element model: neutron diffraction, the contour method, and the slitting method. Model-measurement comparisons are supplemented with detailed simulations that reflect the configurations of the residual-stress measurement processes themselves. The model results show general agreement with experimental measurements, and we observe some similarities in the features around the weld region. Factors that contribute to model-measurement differences are identified. Finally, we conclude with some discussion of the model development and residual stress measurement strategies, including how to best leverage the efforts put forth here for other weld problems.

To decrease operating costs associated with fungal infections in algal crops used for biofuel production, we developed bacterial consortia that displayed antifungal properties. These bacteria were grown in culture with algae species without any additional operating costs or need for re-inoculation with bacteria. These co-cultures maintained their antifungal properties for the entirety of the project period and increased mean time to failure (MTTF) by up to 350% when challenged with high levels of fungal pests. Multiple fungal and fungus-like pests were tested and the consortia showed efficacy against three species.

CrMnFeCoNi, also called the Cantor alloy, is a well-known high-entropy alloy whose magnetic properties have recently become a focus of attention. Here, we present a detailed muon spin relaxation study of the influence of chemical composition and sample processing protocols on the magnetic phase transitions and spin dynamics of several different Cantor alloy samples. Specific samples studied include a pristine equiatomic sample, samples with deficient and excess Mn content, and equiatomic samples magnetized in a field of 9 T or plastically deformed in pressures up to 0.5 GPa. The results confirm the sensitive dependence of the transition temperature on composition and demonstrate that post-synthesis pressure treatments cause the transition to become significantly less homogeneous throughout the sample volume. In addition, we observe critical spin dynamics in the vicinity of the transition in all samples, reminiscent of canonical spin glasses and magnetic materials with ideal continuous phase transitions. Application of an external magnetic field suppresses the critical dynamics in the Mn-deficient sample, while the equiatomic and Mn-rich samples show more robust critical dynamics. The spin-flip thermal activation energy in the paramagnetic phase increases with Mn content, ranging from 3.1⁢(3) × 10^-21 J for 0% Mn to 1.2⁢(2) × 10^-20 J for 30% Mn content. These results shed light on critical magnetic behavior in environments of extreme chemical disorder and demonstrate the tunability of spin dynamics in the Cantor alloy via chemical composition and sample processing.

In this article, global climate change has the potential to alter soil organic carbon (SOC) stocks in rice paddies, because increases in temperature and atmospheric carbon dioxide concentration ([CO₂]) both influence the primary input (i.e., net primary production, NPP) and output (i.e. heterotrophic respiration) of carbon (C). We used two types of open-top chambers representing present conditions (+0°C, +0 ppm CO₂) and projected climate change conditions (+2°C, +200 ppm CO₂) to investigate the net effect of climate change on SOC stock in rice paddy. Additional chambers with elevated temperature only (+2°C, +0 ppm CO₂) allowed us to quantify the individual effects of temperature and [CO₂]. We calculated changes in SOC stock using net ecosystem C balance (NECB) analysis (i.e., the balance between C inputs and outputs). Compared to present conditions, projected climate change did not change grain yield due to a trade-off between the effects of warming and [CO₂] on grain yield components. NPP during the fallow season significantly decreased under combined warming and CO₂, as the impact of warming outweighed that of elevated [CO₂]. However, rice NPP remained unchanged during the cropping season. Warming plus elevated CO₂ increased SOC mineralization by 157–429 %, particularly through warming-induced soil CO₂ emission during the fallow season. Consequently, climate change conditions decreased (119–271 %) NECB values compared to present conditions, primarily through the response to warming. Our findings demonstrate that rice paddies represent positive feedback on climate change, because accelerated C release from warmed soils will override C gains from NPP under elevated CO₂. Reducing SOC depletion in rice paddy agriculture under a changing climate therefore requires conservative soil management practices during the fallow season.

Here, polymerization-induced phase separation is a useful method for the construction of heterogeneous epoxy networks with properties exceeding their homogeneous counterparts. In this work, we examine the static and dynamic thermomechanical properties of phase-separated epoxy networks salient to their application as encapsulants. Three heterogeneous epoxy-amine networks with nano-, meso-, and macro-phase-separated morphologies comprised of hard and soft domains are compared to a rigid, unstructured network. The glass transition profiles of the heterogeneous networks are complex, spanning many decades in the frequency domain. The nanophase-separated morphology leads to higher coefficient of thermal expansion, yet surprisingly is characterized by reduced residual stress. Under both quasi-static and dynamic compression (strain rates of order 10^–3 and 10³ s^–1, respectively), the nanophase-separated network also exhibits higher modulus and strength. In split-Hopkinson bar experiments, the energy dissipation characteristics of the epoxy networks were nearly identical. Curiously, however, the Hugoniot response of the macro-phase-separated network determined by ballistic shockwave analysis indicates a remarkable ability of this material to mitigate shockwave propagation in comparison to many homogeneous and heterogeneous polymer materials. Collectively, this work reveals several previously unreported phenomena with respect to structure–property relationships in phase-separated epoxy networks, illustrating the potential value of systematically tuned microstructures for optimization of application-specific physical properties.

The pitfalls of numerical computations using floating-point numbers are well known. Existing static techniques for floating-point error analysis provide a single upper bound across all potential values for a floating-point variable. We present a new abstract domain for floating-point error analysis which describes error as a function of each variable’s value. This domain accurately models the nature of floating-point error as dependent on the magnitude of its operands. We use this domain to effectively handle exceptional values (e.g., NaN), branch instability, and binade boundaries. The granular analysis provides users with a detailed understanding of forward error. We implement the abstract domain in a tool that supports analyzing a subset of C including conditionals, arrays, and arithmetic operators. We compare our implementation with Fluctuat and show how our analysis can improve the error bounds for subranges of possible outputs.

Data-consistent inversion is designed to solve a class of stochastic inverse problems where the solution is a pullback of a probability measure specified on the outputs of a quantities of interest (QoI) map. Here, this work presents stability and convergence results for the case where finite QoI data result in an approximation of the solution as a density. Given their popularity in the literature, separate results are proven for three different approaches to measuring discrepancies between probability measures: f-divergences, integral probability metrics, and L^p metrics. In the context of integral probability metrics, we also introduce a pullback probability metric that is well-suited for data-consistent inversion. This fills a theoretical gap in the convergence and stability results for data-consistent inversion that have mostly focused on convergence of solutions associated with approximate maps. Numerical results are included to illustrate key theoretical results with intuitive and reproducible test problems that include a demonstration of convergence in the measure-theoretic "almost" sense.

Redox flow batteries (RFBs) are an attractive choice for stationary energy storage of renewables such as solar and wind. Non-aqueous redox flow batteries (NARFBs) have garnered broad interest due to their high voltage operation compared to their aqueous counterparts. Further, the utilization of bipolar redox-active molecules (BRMs) is a practical way to alleviate crossover faced by asymmetric RFBs. In this work, ferrocene (Fc) and phthalimide (PI) are covalently linked with various tethering groups which vary in structure and length. The compiled results suggest that the length and steric shielding ability of the linker group can greatly influence the stability and overall performance of Fc-n-PI BRM-based NARFBs. Primary sources of capacity loss are found to be BRM degradation for straight chain spacers <6 carbons and membrane (Nafion) fouling. Fc-hexyl-PI provided the most stable battery cycling and coulombic efficiencies of >98 % over 100 cycles (~13 days). NARFB using Fc-hexyl-PI as an active material exhibited high working voltage (1.93 V) and maximum capacity (1.28 Ah L−1). Additionally, this work highlights rational strategies to improve cycling stability and optimize NARFB performance.

Recent points of emphasis in the Library of Advanced Materials for Engineering (LAMÉ) have been to enable flexibility in formulations via the adoption of a variety of modular frameworks. While more established phenomenologies such as plasticity and viscoelasticity have been considered, elastically orthotropic models (e.g. elastic_3D_orthotropic) have not. For the elasticity component, not much can be modularized. However, a potential feature of interest would be the evaluation of failure criteria to consider the possibility of damage. Many such forms exist in the literature providing a good basis for modularity.

Presented in this document is a small portion of the tests that exist in the Sierra/SolidMechanics (Sierra/SM) verification test suite. Most of these tests are run nightly with the Sierra/SM code suite, and the results of the test are checked versus the correct analytical result. For each of the tests presented in this document, the test setup, a description of the analytic solution, and comparison of the Sierra/SM code results to the analytic solution is provided. Mesh convergence is also checked on a nightly basis for several of these tests. This document can be used to confirm that a given code capability is verified or referenced as a compilation of example problems. Additional example problems are provided in the Sierra/SM Example Problems Manual. Note, many other verification tests exist in the Sierra/SM test suite, but have not yet been included in this manual.

The performance and reliability of many structures and components depend on the integrity of interfaces between dissimilar materials. Interfacial toughness Γ is the key material parameter that characterizes resistance to interfacial crack growth, and Γ is known to depend on many factors including temperature. For example, previous work showed that the toughness of an epoxy/aluminum interface decreased 40 % as the test temperature was increased from −60 °C to room temperature (RT). Interfacial integrity at elevated temperatures is of considerable practical importance. Recent measurements show that instead of continuing to decrease with increasing temperature, Γ increases when test temperature is above RT. Cohesive zone finite element calculations of an adhesively bonded, asymmetric double cantilever beam specimen of the type used to measure Γ suggest that this increase in toughness may be a result of R-curve behavior generated by plasticity-enhanced toughening during stable subcritical crack growth with interfacial toughness defined as the critical steady-state limit value. In these calculations, which used an elastic-perfectly plastic epoxy model with a temperature-dependent yield strength, the plasticity-enhanced increase in Γ above its intrinsic value Γo depended on the ratio of interfacial strength σ* to the yield strength σyb of the bond material. There is a nonlinear relationship between Γ/Γo and σ*/σyb with the value Γ/Γo increasing rapidly above a threshold value of σ*/σyb. The predicted increase in toughness can be significant. For example, there is nearly a factor of two predicted increase in Γ/Γo during micrometer-scale crack-growth when σ*/σyb = 2 (a reasonable choice for σ*/σyb). Furthermore, contrary to other reported results, plasticity-enhanced toughening can occur prior to crack advance as the cohesive zone forms and the peak stress at the tip of the original crack tip translates to the tip of the fully formed cohesive zone. These results suggest that plasticity-enhanced toughening should be considered when modeling interfaces at elevated temperatures.

Tailorable discontinuous fiber composite laminates provide relative formability beyond that of continuous fiber laminates, while achieving improved mechanical performance over comparable stochastic systems. Here, in this work, the notch sensitivity of engineered prepreg platelet molded composite (PPMC) laminates is investigated using the open-hole tension (OHT) test and compared to available data for stochastic PPMCs and continuous fiber laminates made with the same material. The press-formed thermoplastic composites (AS4/PEKK) were molded with a quasi-isotropic stacking sequence. The discontinuous PPMC laminate was found to be notch insensitive with OHT strengths ranging from 145.4 MPa (CV $=$ 7%) for d/w $=$ 0.5 to 229.3 MPa (CV $=$ 9%) for d/w $=$ 0.25. The highly ordered meso-structure of the engineered PPMC laminate yields comparatively excellent mechanical properties for relatively thin laminates in contrast to stochastic systems. Both net- and gross-section failures were observed for d/w $=$ 0.25, which suggests that the engineered PPMC laminates studied here maintain a degree of inherent, internal stress concentrations that compete with those caused by geometric features such as a circular hole. Computational simulations of the OHT tests with explicitly represented platelets were found to be in good agreement with experimental measurements. The progressive failure analysis was used to conduct a numerical investigation of the stacking sequence and platelet meso-morphology.

Presented in this document is a portion of the tests that exist in the Sierra Thermal/Fluids verification test suite. Each of these tests is run nightly with the Sierra/TF code suite and the results of the test checked under mesh refinement against the correct analytic result. For each of the tests presented in this document the test setup, derivation of the analytic solution, and comparison of the code results to the analytic solution is provided.

The Integrated Tiger Series (ITS) generates a database containing energy deposition data. This data, when stored on an Exodus file, is not typically suitable for analysis within Sierra Mechanics for finite element analysis. The its2sierra tool maps data from the ITS database to the Sierra database. This document provides information on the usage of its2sierra.

Accurate and efficient constitutive modeling remains a cornerstone issue for solid mechanics analysis. Over the years, the LAMÉ advanced material model library has grown to address this challenge by implementing models capable of describing material systems spanning soft polymers to stiff ceramics including both isotropic and anisotropic responses. Inelastic behaviors including (visco)plasticity, damage, and fracture have all incorporated for use in various analyses. This multitude of options and flexibility, however, comes at the cost of many capabilities, features, and responses and the ensuing complexity in the resulting implementation. Therefore, to enhance confidence and enable the utilization of the LAMÉ library in application, this effort seeks to document and verify the various models in the LAMÉ library. Specifically, the broader strategy, organization, and interface of the library itself is first presented. The physical theory, numerical implementation, and user guide for a large set of models is then discussed. Importantly, a number of verification tests are performed with each model to not only have confidence in the model itself but also highlight some important response characteristics and features that may be of interest to end-users. Finally, in looking ahead to the future, approaches to add material models to this library and further expand the capabilities are presented.

NasGen provides a path for migration of structural models from Nastran bulk data format (BDF) into both an Exodus mesh file and an ASCII input file for Sierra Structural Dynamics (Salinas) and Solid Mechanics (Adagio). Many tools at Sandia National Labs (SNL) use the Exodus format. This document describes capabilities and limitations of the NasGen translation software.

This report provides society information that may include news, reviews or technical notes that should be of interest to practitioners and researchers.

Presented in this document are tests that exist in the Sierra/SolidMechanics example problem suite, which is a subset of the Sierra/SM regression and performance test suite. These examples showcase common and advanced code capabilities. A wide variety of other regression and verification tests exist in the Sierra/SM test suite that are not included in this manual.

Sierra/SD provides a massively parallel implementation of structural dynamics finite element analysis, required for high fidelity, validated models used in modal, vibration, static and shock analysis of structural systems. This manual describes the theory behind many of the constructs in Sierra/SD. For a more detailed description of how to use Sierra/SD, we refer the reader to User’s Manual. Many of the constructs in Sierra/SD are pulled directly from published material. Where possible, these materials are referenced herein. However, certain functions in Sierra/SD are specific to our implementation. We try to be far more complete in those areas. The theory manual was developed from several sources including general notes, a programmer_notes manual, the user’s notes and of course the material in the open literature.

This is an addendum to the Sierra/SolidMechanics 5.22 User’s Guide that documents additional capabilities available only in alternate versions of the Sierra/SolidMechanics (Sierra/SM) code. These alternate versions are enhanced to provide capabilities that are regulated under the U.S. Department of State’s International Traffic in Arms Regulations (ITAR) export control rules. The ITAR regulated codes are only distributed to entities that comply with the ITAR export control requirements. The ITAR enhancements to Sierra/SM include material models with an energy-dependent pressure response (appropriate for very large deformations and strain rates) and capabilities for blast modeling. This document is an addendum only; the standard Sierra/SolidMechanics 5.22 User’s Guide should be referenced for most general descriptions of code capability and use.

This report describes an approach to utilizing phasor measurement unit (PMU) data from multiple Intelligent Electronics Devices (IEDs) in a low-voltage network to produce a differential scheme for protecting the medium-voltage feeder and low-voltage network transformers. The proposed protection scheme is designed and prototyped on a real-time automation controller. Its performance is evaluated using real-time controller hardware-in-the-loop simulation. Lab testing results indicate that the proposed protection scheme allows significant distributed energy resources (DER) backfeed and enables selective and fast protection of medium voltage feeders.

Photonic computing has the potential to harness the full degrees of freedom (DOFs) of the light field, including the wavelength, spatial mode, spatial location, phase quadrature, and polarization, to achieve a higher level of computing parallelism and scalability than digital electronic processors. While multiplexing using the wavelength and other DOFs can be readily integrated on silicon photonics platforms with compact footprints, conventional mode-division multiplexed (MDM) photonic designs occupy areas exceeding tens to hundreds of microns for a few spatial modes, significantly limiting their scalability. Here, we utilize inverse design to demonstrate an ultracompact photonic computing core that calculates vector dot products based on MDM coherent mixing. Our dot-product core integrates the functionalities of two-mode multiplexers and one multimode coherent mixer within a nominal footprint of 5 μm x 3 μm . We have experimentally demonstrated computing examples on the fabricated dot-product core, including complex number multiplication and motion estimation using optical flow. The compact dot-product core design enables large-scale on-chip integration in a parallel photonic computing primitive cluster for high-throughput scientific computing and computer vision tasks.

Exploding bridgewire detonators (EBWs) containing pentaerythritol tetranitrate (PETN) exposed to high temperatures may not function following discharge of the design electrical firing signal from a charged capacitor. Knowing functionality of these arbitrarily facing EBWs is crucial when making safety assessments of detonators in accidental fires. Orientation effects are only significant when the PETN is partially melted. The melting temperature can be measured with a differential scanning calorimeter. Nonmelting EBWs will be fully functional provided the detonator never exceeds 406 K (133 °C) for at least 1 h. Conversely, EBWs will not be functional once the average input pellet temperature exceeds 414 K (141 °C) for a least 1 min which is long enough to cause the PETN input pellet to completely melt. Functionality of the EBWs at temperatures between 406 and 414 K will depend on orientation and can be predicted using a stratification model for downward facing detonators but is more complex for arbitrary orientations. A conservative rule of thumb would be to assume that the EBWs are fully functional unless the PETN input pellet has completely melted.

Abstract not provided.

An overview of recent development of our MEG system. Most of the included data/figures are already published in the open literature.

The Department of Energy (DOE)’s Technical Standard DOE-STD-1194-2019 (dated September 2019), Nuclear Materials Control and Accountability, provides key guidance for the determination of Special Nuclear Material (SNM) attractiveness levels. Attractiveness levels are a key component in the security categorization of SNM processed, used, and stored at DOE facilities. Upon review, the writing team identified specific components relating to the determination of attractiveness levels that could be modified to improve clarity, reduce burden on sites, and/or better align with the graded safeguards principle that is central to the DOE’s nuclear security program. The report outlines the original verbiage, issues with implementation of that verbiage, proposed new verbiage, and the expected benefits thus serving as the technical basis for the proposed changes.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The synthesis of Lanthanum Cobalt Oxide (LCO) via Sol-Gel method provides a potentially low-cost method of production for a P-type photocatalytic. LCO was prepped from Lanthanum and Nitrate Precursors in Aqueous solution. After a significant amount of water is evaporated, the solution is deposited onto SiO2 substrate via spin-deposition method. After annealing, the samples produce thin film LCO that display thickness of sub-500 nanometers. The samples material profile is confirmed by both Raman spectroscopy and X-Ray Diffraction spectroscopy (XRD). Additionally, two-point probing tested the conductivity of the samples. The thin film samples display characteristics of a P-type photocatalytic and may potentially be used in the formation of a P-N junction for user in water-splitting applications.

The adoption of machine learning approaches for gamma-ray spectroscopy has received considerable attention in the literature. Many studies have investigated the deployment of various algorithm architectures to a specific task. However, little attention has been afforded to the development of the datasets leveraged to train the models. Such training datasets typically span a set of environmental or detector parameters to encompass a problem space of interest to a user. Variations in these measurement parameters will also induce fluctuations in the detector response, including expected pile-up and ground scatter effects. Fundamental to this work is the understanding that 1) the underlying spectral shape varies as the measurement parameters change and 2) the statistical uncertainties associated with two spectra impact their level of similarity. While previous studies attribute some arbitrary discretization to the measurement parameters for the generation of their synthetic training data, this work introduces a principled methodology for efficient spectral-based discretization of a problem space. A signal-to-noise ratio (SNR) respective spectral comparison measure and a Gaussian Process Regression (GPR) model are used to predict the spectral similarity across a range of measurement parameters. This innovative approach effectively showcased its capability by dividing a problem space, ranging from 5 cm to 100 cm standoff distances and 5 μCi–100 μCi of 137Cs, into three unique combinations of measurement parameters. The findings from this work will aid in creating more robust datasets, which incorporate many possible measurement scenarios, reduce the number of required experimental test set measurements, and possibly enable experimental training data collection for gamma-ray spectroscopy.

Abstract not provided.

Determining the depths of small crustal earthquakes is challenging in many regions of the world, because most seismic networks are too sparse to resolve trade-offs between depth and origin time with conventional arrival-time methods. Precise and accurate depth estimation is important, because it can help seismologists discriminate between earthquakes and explosions, which is relevant to monitoring nuclear test ban treaties and producing earthquake catalogs that are uncontaminated by mining blasts. Here, we examine the depth sensitivity of several physics-based waveform features for ∼8000 earthquakes in southern California that have well-resolved depths from arrival-time inversion. We focus on small earthquakes (2 < ML < 4) recorded at local distances (< 150 km), for which depth estimation is especially challenging. We find that differential magnitudes (Mw= ML–Mc) are positively correlated with focal depth, implying that coda wave excitation decreases with focal depth. We analyze a simple proxy for relative frequency content, Φ≡ log10 (M0)+3log10 (fc (,and find that source spectra are preferentially enriched in high frequencies, or “blue-shifted,” as focal depth increases. We also find that two spectral amplitude ratios Rg 0.5–2 Hz/Sg 0.5–8 Hz and Pg/Sg at 3–8 Hz decrease as focal depth increases. Using multilinear regression with these features as predictor variables, we develop models that can explain 11%–59% of the variance in depths within 10 subregions and 25% of the depth variance across southern California as a whole. We suggest that incorporating these features into a machine learning workflow could help resolve focal depths in regions that are poorly instrumented and lack large databases of well-located events. Some of the waveform features we evaluate in this study have previously been used as source discriminants, and our results imply that their effectiveness in discrimination is partially because explosions generally occur at shallower depths than earthquakes.

Abstract not provided.

Abstract not provided.

Abstract not provided.

This paper details a computational framework to produce automated, graphical workflows, and how this framework can be deployed to support complex modeling problems like those in nuclear engineering. Key benefits of the framework include: automating previously manual workflows; intuitive construction and communication of workflows through a graphical interface; and automated file transfer and handling for workflows deployed across heterogeneous computing resources. This paper demonstrates the framework's application to probabilistic post-closure performance assessment of systems for deep geologic disposal of nuclear waste. However, the framework is a general capability that can help users running a variety of computational studies.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

A series of drained and undrained water-saturated constant mean-stress tests were performed to investigate the strength, elasticity, and poroelastic response of a water-saturated high porosity nonwelded tuff. Drained strengths are found to increase with increasing effective confining pressures. Elastic moduli increase with increasing mean stress. Undrained strengths are small due to development of high pore pressures that generate low effective confining pressures. Skempton’s values are pressure dependent and appear to reflect the onset of inelastic deformation. Permeabilities decrease after deformation from ∼ 10–14 to ∼ 10–16 m2 and are a function of the applied confining pressure. Deformation is dominated by pore collapse, compaction, and intense microfracturing, with the undrained tests favoring microfracture-dominant deformation and the drained tests favoring compaction-dominant deformation. These property determinations and observations are used to develop/parameterize physics-based models for underground explosives testing.

The computational tool described in this report helps identify promising biological pathways that produce SAF platform molecules (either a drop-in SAF, or a precursor that can be easily converted to a drop-in SAF). The workflow the computational tool follows first identifies possible biological pathways from a user-defined metabolite. These pathways may, or may not lead to a SAF platform molecule, thus the second step involves insilico testing of the end product of each pathway to assess whether it is, or is not, a SAF platform molecule. The identification of biological pathways performed in the first step is facilitated by linking the metabolite to a biological reaction database. Pathways are found by identifying pathways in the reaction database that include the metabolite. The computational tool includes an alternative way to find pathways. The alternative way develops a Flux Balanced Analysis (FBA), and modifying the FBA to include reactions that transform the metabolite. These modifications serve as a basis for understanding, in a semi-quantitative way, if there is an increase in the flux to desirable products. The second step, in silico testing of the end-products, is accomplished by estimating key physical properties relevant to SAF. When good models are available, we have integrated those models into the computational tool. In a few instances, we have developed our own models. In all instances, we have validated the models against available measured data. Finally, we have evaluated the effectiveness of our computational tool by genetically engineering Rhodosporidium toruloides. Validation occurred without the use of a FBA, and further validation is required.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

This report summarizes the research and capabilities developed as part of the project “Analysis and Validation of Critical Infrastructure Models using Model Order Reduction” (CI-MOR) LDRD project. CI-MOR research enables the solution of large, complex optimization models that naturally arise in national security challenges involving critical infrastructures. Specifically, CI-MOR researchers developed methods to (1) rigorously approximate complex, nonlinear optimization formulations, (2) identify alternative near-optimal solutions, (3) accelerate optimization workflows used for complex applications, and (4) rigorously integrate domain knowledge in stochastic-process models. This report provides an overview of the research done in CI-MOR, and we describe application exemplars used to illustrate CI-MOR capabilities. Furthermore, we describe the software developed by CI-MOR that researchers can leverage to analyze new applications.

Abstract not provided.

The current instrumentation for observing the complex flow fields in and around wind plants struggles to match the fidelity of existing simulation tools. As a result, these measurement limitations create a hurdle for validating and assessing the quality of the wind plant numerical models. This roadmap for instrumentation development recommendations was created to offer guidance on narrowing the gap between measurement and simulation fidelity. A process was established to identify where gaps in instrumentation exist for wind energy test campaigns by analyzing the capabilities of instrumentation for capturing the various important phenomena at the necessary resolution for both the science goal and validation objectives. To this end, a multi-disciplinary team of experts on instrumentation, wind energy, and atmospheric science was assembled to identify these significant instrumentation needs. A recommendation for instrumentation to be developed is provided, and the framework developed through this process is expected to be useful to the design of future test campaigns. The mapping tools developed for this process will be distributed as part of a future International Energy Agency Wind Technology Collaboration Program task on instrumentation development.

We investigate hydrodynamic fluctuations in the flow past a circular cylinder near the critical Reynolds number Rec for the onset of vortex shedding. Starting from the fluctuating Navier-Stokes equations, we perform a perturbation expansion around Rec to derive analytical expressions for the statistics of the fluctuating lift force. Molecular-level simulations using the direct simulation Monte Carlo method support the theoretical predictions of the lift power spectrum and amplitude distribution. Notably, we have been able to collect sufficient statistics at distances Re/Rec-1=O(10-3) from the instability that confirm the appearance of non-Gaussian fluctuations, and we observe that they are associated with intermittent vortex shedding. These results emphasize how unavoidable thermal-noise-induced fluctuations become dramatically amplified in the vicinity of oscillatory flow instabilities and that their onset is fundamentally stochastic.

Abstract not provided.

Future bulk power systems are expected to operate with a high penetration of inverter-based resources (IBR). This is anticipated to fundamentally change the system’s short-circuit behaviors and demand a change in existing transmission line protection settings or schemes. Phasor-domain short circuit programs (e.g., ASPEN, CAPE, and CYME) remain the mainstream tools for developing protection settings. However, the system protection community has shared concerns about the accuracy of IBR models in short-circuit programs. This report compares the performance of existing generic IBR models in short-circuit programs with detailed electro-magnetic transient (EMT) models provided by several different IBR manufacturers in PSCAD, to quantify the short-circuit program’s IBR model accuracy, assess their ability to replicate individual IBR controls such as negative sequence current injections, and identify specific gaps in the existing generic phasor-domain IBR models that are currently available.

Abstract not provided.

Presentation at IEEE Brain Discovery and Neurotechnology Workshop https://brain.ieee.org/2024-ieee-brain-discovery-neurotechnology-workshop/

In this project we considered the initial stages of helium bubble nucleation via the proposed mechanism of self-interstitial atom nucleation. By calculating the energy barrier to self-interstitial atom nucleation in a range of Fe-Ni-Cr alloys we identified the most important energetic contributions to the phenomenon: the Frenkel-pair energy barrier in the absence of helium and the difference of insertion energy for a He cluster into a perfect lattice and vacancy. From this observation, we developed a simple model of helium-assisted self-interstitial atom nucleation.

Abstract not provided.

Liquid organic hydrogen carrier (LOHC) systems are an excellent alternative to pressurized gas and liquid hydrogen storage technologies due to their high volumetric storage capacities and straightforward adaptation to existing infrastructure. Here, we investigate various molecular and heterogenized Ru, Mn, and Fe catalysts for the reversible (de)hydrogenation of polyols. Mn catalysts, Mn-MACHO(Ph) (1) and Mn-MACHO(iPr)BH₄ (2-BH₄) were found to maintain catalytic activity for hydrogen production comparable to the Ru analogs, with greater than 98% conversion of 1,4 butanediol and quantitative hydrogen production. Furthermore, to assess the viability of utilizing heterogenized molecular catalysts for 1,4 butanediol (de)hydrogenation, molecular catalysts, Mn-MACHO(Ph), Ru-MACHO(Ph), and Ru-9, were polymerized to form Mn-MACHO-Poly, Ru-MACHO-Poly, and Ru-9-Poly respectively. These catalysts were then used to assess (de)hydrogenation of polyols, ethylene glycol and 1,4 butanediol. These studies reveal that Ru-MACHO-Poly is an efficient dehydrogenation catalyst with 99% conversion and a hydrogen percent yield of 96%. In addition, Ru-MACHO-Poly is a competent hydrogenation catalyst with approximately 98% conversion back to 1,4 butanediol based on quantitative NMR measurements. Overall, the data suggest these molecular and heterogenized catalysts have potential for practical use in polyalcohol-based LOHC systems.

Methane (CH₄), an abundant greenhouse gas, is the second largest contributor to global warming after carbon dioxide (CO₂). In comparison to CO₂, CH₄ has a larger warming effect over a much shorter lifetime. While technologies to radically reduce global carbon dioxide emissions are materializing, rapid reductions in methane emissions are needed to limit near-term warming. Methane is primarily emitted as a byproduct from agricultural activities and energy extraction/utilization and is currently monitored via bottom-up (i.e., activity level) or top-down (via airborne or satellite retrievals) approaches. However, significant methane leaks remain undetected, and emission rates are challenging to characterize with current monitoring frameworks. In this report, we study methane leaks from oil and gas infrastructure using a tiered monitoring approach that combines bottom-up and top-down approaches in an integrated framework. We describe the individual advantages of bottom-up and top-down sensors in both stationary and mobile settings before characterizing how a fully integrated framework can improve predictions and uncertainties of potential leak locations and their emission rates. Further, we study the impact of different atmospheric (wind) conditions on integrated methane monitoring and develop a probabilistic approach to optimal sensor placement, thereby shortening detection times and improving monitoring capabilities. Last, we discuss how biogenic flux modeling can be used to improve assessment of background methane concentrations needed to fully assess the sensitivity of a tiered monitoring system.

Commercial high-density polyethylene (HDPE) is lightweight with excellent mechanical properties but has limited thermal stability. To enhance this, lanthanide-based nanomaterials like CeO₂ were mixed with HDPE. Using a twin-screw extruder improved dispersion. This study examines the effects of neutron, X-rays and gamma radiation on the crystallinity of polymer of CeO₂-HDPE composites.

Infrasound sensing plays a critical role in the detection and analysis of bolides, offering passive, cost-effective global monitoring capabilities. Key objectives include determining the timing, location, and yield of these events. Achieving these goals requires a robust approach to detect, analyze, and interpret rapidly moving elevated sources such as bolides (also re-entry). In light of advancements in infrasonic methodologies, there is a need for a comprehensive overview of the characteristics that distinguish bolides from other infrasound sources and methodologies for bolide infrasound analysis. This paper provides a focused review of key considerations and presents a unified framework to enhance infrasound processing approaches specifically tailored for bolides. Three representative case studies are presented to demonstrate the practical application of infrasound processing methodologies and deriving source parameters while exploring challenges associated with bolide-generated infrasound. These case studies underscore the effectiveness of infrasound in determining source parameters and highlight interpretative challenges, such as variations in signal period measurements across different studies. Future research should place emphasis on improving geolocation and yield accuracy. This can be achieved through rigorous and systematic analyses of large, statistically significant samples of such events, aiming to resolve interpretative inconsistencies and explore the causes for variability in signal periods and back azimuths. The topic described here is also relevant to space exploration involving planetary bodies with atmospheres, such as Venus, Mars, and Titan.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The Glassy Amorphous Polymer (GAP) model is a viscoelastic/plastic model developed at Los Alamos National Laboratory to accurately model a variety of polymers across a wide range of conditions and loading rates, including shock loading. In the present report we introduce and assess this model, newly implemented in the ALEGRA shock and multiphysics code, using a series of verification and application-related validation problems. We describe the mathematical and theoretical formulation of the model, as well as its implementation in ALEGRA, in detail. We provide verification results that assess the model implementation against published computational results, as well as validation results which we compare to existing experimental results when possible. These comparisons instill confidence in the implementation and indicate that the addition of the GAP model to the ALEGRA code provides users with a high-fidelity polymer modeling capability that is capable of recreating complex polymer phenomena.

This an abstract for the fall Electrochemical Society Conference (ECS)

Although there are increasing numbers of distributed energy resources (DERs) and microgrids being deployed, current IEEE and utility standards generally strictly limit their interconnection inside secondary networks. Secondary networks are low-voltage meshed (non-radial) distribution systems that create redundancy in the path from the main grid source to each load. This redundancy provides a high level of immunity to disruptions in the distribution system, and thus extremely high reliability of electric power service. There are two main types of secondary networks, called grid and spot secondary networks, both of which are used worldwide. In the future, primary networks in distribution systems that might include looped or meshed distribution systems at the primary-voltage (medium-voltage) level may also become common as a means for improving distribution reliability and resilience.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Rapid, time-varying, three-dimensional physics underpin numerous engineering challenges. Often, these physics occur within opaque environments, internal to a component, severely limiting applicable diagnostics. Development of novel diagnostics is necessary to understand and predict transient three-dimensional (3D) phenomena within opaque environments. This report highlights progress in four key areas leading to advancements in high-speed X-ray radiography and tomography. The first area is enabling MHz-rate imaging of energetics at the Advanced Photon Source at Argonne National Laboratory. The second is modeling a high-flux, rotating-anode X-ray source to understand the heat loads on the anode. The third effort was to develop a novel reconstruction algorithm that is validated by ground experimental tomography data and synthetic tomography data. The fourth is the development of a novel approach to two-color X-ray imaging.

Sandia National Laboratories (SNL) is in the process of creating Inspecta (International Nuclear Safeguards Personal Examination and Containment Tracking Assistant), an Artificial Intelligence (AI)-powered smart digital assistant (SDA) with robotic capabilities, aimed at enhancing the effectiveness, efficiency, and safety of international nuclear safeguards inspections. This innovative tool is designed to assist inspectors on-site by supporting or automating tasks that are typically mundane, hazardous, or susceptible to errors. In 2021, the development team established the specifications for Inspecta by analyzing International Atomic Energy Agency (IAEA) documents and consulting with former IAEA inspectors and subject matter experts. This process involved aligning in-field inspection tasks with existing commercial or open-source technologies to outline a roadmap for the initial prototype of Inspecta, while also identifying areas needing further research and development. From 2022 – 2024, the focus has shifted to integrating a critical inspection activity, the examination of seals, into an early version of Inspecta. This has involved developing both the software and hardware capabilities necessary for this task. This report outlines the ongoing advancements in Inspecta's functionalities, specifically those supporting the seal examination process.

A novel algorithm for explicit temporal discretization of the variable-density, low-Mach Navier-Stokes equations is presented here. Recognizing there is a redundancy between the mass conservation equation, the equation of state, and the transport equation(s) for the scalar(s) which characterize the thermochemical state, and that it destabilizes explicit methods, we demonstrate how to analytically eliminate the redundancy and propose an iterative scheme to solve the resulting transformed scalar equations. The method obtains second-order accuracy in time regardless of the number of iterations, so one can terminate this subproblem once stability is achieved. Hence, flows with larger density ratios can be simulated while still retaining the efficiency, low cost, and parallelizability of an explicit scheme. The temporal discretization algorithm is used within a pseudospectral direct numerical simulation which extends the method of Kim, Moin, and Moser for incompressible flow [17] to the variable-density, low-Mach setting, where we demonstrate stability for density ratios up to ∼25.7.

I am going to present my work on multifidelity timeseries models at MS&T in Pittsburgh. We develop efficient machine learning methodologies to accelerate time-series predictions from a hierarchy of complex physics-based models.

Abstract not provided.

Abstract not provided.

Rattlesnake is a combined-environments, multiple input/multiple output control system for dynamic excitation of structures under test. It provides capabilities to control multiple responses on the part using multiple exciters using various control strategies. Rattlesnake is written in the Python programming language to facilitate multiple input/multiple output vibration research by allowing users to prescribe custom control laws to the controller. Rattlesnake can target multiple hardware devices, or even perform synthetic control to simulate a test virtually. Rattlesnake has been used to execute control problems with up to 200 response channels and 24 shaker drives. This document describes the functionality, architecture, and usage of the Rattlesnake controller to perform combined environments testing.

Abstract not provided.

Abstract not provided.

The Hydrogen Plus Other Alternative Fuels Risk Assessment Models (HyRAM+) software has seen various improvements and additional physics capabilities since validation against experimental data was last published for HyRAM v3.1. Notably, HyRAM+ now includes four models allowing for the calculation of overpressure resulting from vapor cloud explosions from unconfined jet releases. As with the previous HyRAM validation report, validation data was gathered from available published literature and tested against HyRAM+ capabilities. The validation comparisons include tank blowdown, unignited dispersion jet plume, ignited jet flame, and enclosed accumulation and overpressure. The unconfined overpressure calculations in HyRAM+ v5.1.1 generally show good agreement with many of the experimental data sets for all four unconfined overpressure models, though HyRAM+ overpredicts the experimental data for small and cryogenic hydrogen releases. The comparisons for the other HyRAM+ physics models are largely unchanged from the previously published validation report.

Hydrogen storage systems are becoming more widely deployed throughout the country, and as their presence continues to grow, it is possible that individual and interconnected systems will be exposed to cyber-attacks. These events can cause physical and financial harm to employees, people in the vicinity, and to the company that owns the facility. The two main mechanisms malicious actors may access information or control from a hydrogen storage facility are through information technology and operations technology devices, the former of which refers to data and information from networked devices and the latter of which refers to onsite controls for the physical system. Both types of entryways into the system should be considered when facility managers conduct cyber risk assessments and when regulators develop or revise relevant codes and standards. This report analyzes cybersecurity risks applicable to a wide variety of hydrogen storage systems by outlining the system's purpose and the importance of its cybersecurity. The hydrogen storage system architecture and communication protocols are provided to understand potential cyber vulnerabilities. Later, an event tree analysis is performed on hydrogen operation to identify system weaknesses by outlining potential attack scenarios. This report also identifies critical cyber assets related to different hydrogen operations followed by an examination of potential threats, and the impact of cyber assets on those operational assets.

Abstract not provided.

Modern concepts for next generation pulsed power (NGPP) are slated to deliver up to ten times the energy of Z today. An increase of this magnitude is concerning insofar that Z currently exhibits sizable amounts of inner magnetically insulated transmission line (MITL) loss current on the order of 5-10%. Loss phenomenon in these systems are complex and electrode heating and subsequent thermal desorption are a leading cause. Rapid heat-driven thermal desorption of contaminants scales as the square of the current. Therefore, even a modest doubling of drive current would yield an ~ 4X in non-linear surface electrode heating, quickening thermal desorption-based current loss. Exacerbating these physics is a current inability to measure ultra fast heating rates (>20°C/ns), which are paramount to benchmarking and code validation critical to NGPP design – as an empirical approach is not viable. Therefore, Ultrafast Photoluminescent Surface Heating Optical Thermometry (UP-SHOT) was developed as a new diagnostic for measurement of GHz-scale electrode heating. The discovery of UP-SHOT leveraged expertise in Engineering Science, Material Science, Pulsed-Power, and the Center for Integrated Nanotechnologies. This report includes information on: 1) The preparation of zinc oxide (ZnO) films, characterization, post-deposition treatments 2) Time-resolved photoluminescence at elevated temperatures and thermographic sensitivity

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

We numerically investigate the mechanisms that resulted in induced seismicity occurrence associated with CO2 injection at the Illinois Basin–Decatur Project (IBDP). We build a geologi-cally consistent model that honors key stratigraphic horizons and 3D fault surfaces inter-preted using surface seismic data and microseismicity locations. We populate our model with reservoir and geomechanical properties estimated using well-log and core data. We then performed coupled multiphase flow and geomechanics modeling to investigate the impact of CO2 injection on fault stability using the Coulomb failure criteria. We calibrate our flow model using measured reservoir pressure during the CO2 injection phase. Our model results show that pore-pressure diffusion along faults connecting the injection inter-val to the basement is essential to explain the destabilization of the regions where micro-seismicity occurred, and that poroelastic stresses alone would result in stabilization of those regions. Slip tendency analysis indicates that, due to their orientations with respect to the maximum horizontal stress direction, the faults where the microseismicity occurred were very close to failure prior to injection. These model results highlight the importance of accurate subsurface fault characterization for CO2 sequestration operations.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The SIREN (Scaling Ion Traps by Requiring iNnovative heterogenous integration) project explores the feasibility of heterogeneous integration (HI) as a transformative approach to scaling ion traps, a critical technology for advancing quantum computers and atomic clocks. Traditional ion trap architectures face significant challenges in scalability due to limitations in optical access, fabrication techniques, and material constraints. SIREN addresses these challenges by leveraging HI, which combines different materials and fabrication processes to create more complex and efficient ion trap structures. HI integrated structures can be manufactured without compromising the process to maintain compatibility to ion traps. This project focuses on integrating a separately fabricated waveguide with a fully functional ion trap. The respective alignment between the pieces needs to be accurate to less than 2 µm to ensure that the light from the waveguide can overlap with the trapping region. The fine alignment must also be maintained through an ultra-high vacuum bake, a critical step in preparing an ion trap experiment. The project's outcomes suggest that heterogeneous integration is a promising pathway for overcoming current scalability barriers, paving the way for the next generation of quantum technologies. SIREN's findings contribute significantly to the field, offering a scalable solution that could accelerate the development of practical quantum computers and highly accurate atomic clocks.

Melanized species of filamentous fungi isolated from high radiation environments have been reported to exhibit radiotropism, defined as the directed growth toward a source of ionizing radiation. Inconsistencies in the experimental approaches and results have impeded our understanding of the key factors involved in radiotropism. In the present study, we assessed radiotropism in four isolates of fungi: Aspergillus niger, A. calidoustus JC-1043, Paecilomyces variotii SNL-1, and P. variotii IMV-00236. Of these fungi, only P. variotii IMV-00236 had been previously reported to exhibit radiotropic behavior. Plates of each fungus were placed in equivalent proximity to a ¹³⁷Cs source, with a primary gamma emission of 662 keV, and differences in the rate and direction of mycelia growth were measured over a seven-day period. Significant differences were not observed in the rate or direction of growth of the different fungi based on exposure to gamma radiation, which suggested a lack of measurable radiotropism in these experiments. Additional studies varying parameters such gamma emission rates and energies, as well as other types of ionizing radiation (e.g., alpha and beta particles, neutrons) are necessary to gain further insights to the factors critical to the expression of radiotropic behavior in filamentous fungi.

Rattlesnake is a combined-environments, multiple input/multiple output control system for dynamic excitation of structures under test. It provides capabilities to control multiple responses on the part using multiple exciters using various control strategies. Rattlesnake is written in the Python programming language to facilitate multiple input/multiple output vibration research by allowing users to prescribe custom control laws to the controller. Rattlesnake can target multiple hardware devices, or even perform synthetic control to simulate a test virtually. Rattlesnake has been used to execute control problems with up to 200 response channels and 24 shaker drives. This document describes the functionality, architecture, and usage of the Rattlesnake controller to perform combined environments testing.

Abstract not provided.

Abstract not provided.

A series of reactive-transport models of Enhanced Geothermal Systems (EGS) were constructed using the reactive transport code PFLOTRAN to examine the effect of matrix thermal contraction and mineral dissolution/precipitation on fracture flow in the context of grid cell size and model complexity. It was found that for thermal drawdown at production well, the impact of fracture zone grid cell size is negligible.

Abstract not provided.

Although there are increasing numbers of distributed energy resources (DERs) and microgrids being deployed, current IEEE and utility standards generally strictly limit their interconnection inside secondary networks. Secondary networks are low-voltage meshed (non-radial) distribution systems that create redundancy in the path from the main grid source to each load. This redundancy provides a high level of immunity to disruptions in the distribution system, and thus extremely high reliability of electric power service. There are two main types of secondary networks, called grid and spot secondary networks, both of which are used worldwide. In the future, primary networks in distribution systems that might include looped or meshed distribution systems at the primary-voltage (medium-voltage) level may also become common as a means for improving distribution reliability and resilience.

Abstract not provided.

Microchannel heat exchanger technology is being pursued for next generation CSP concepts for primary power cycle heat addition and power cycle heat recuperation due to the high heat transfer coefficients and pressure containment advantages of small sCO₂ channels. The economics of future CSP plants as dictated by the SETO 2020 or 2030 targets depend on a heat exchanger with a 30-year lifetime (resisting creep, fatigue, corrosion, erosion) and operational characteristics such as fast ramping and the ability to withstand thermal shock. However, the lifetime and operational limits of microchannel heat exchangers operating at high-temperatures, particularly those constructed from high-nickel alloys, are not well known. This uncertainty has resulted in heat exchanger vendors not being able to accurately forecast heat exchanger lifetime as required by customers, specify operational limits as required by process engineers to prevent premature heat exchanger failure, or overdesign heat exchanger which leads to higher cost than necessary. Our goal is to evaluate heat exchanger lifetime and operational limits for the manufacturing and prototype design for next-generation CSP heat exchanger technology through a combination of collecting experimental data and modeling studies.

Abstract not provided.

Abstract not provided.

Development of a radioimaging diagnostic for high-voltage component reliability testing and electrical breakdown computational model validation is described. Radioimaging has its roots in radio astronomy, where aperture synthesis (also known as synthesis imaging) has been utilized for decades to image radio sources far from Earth. Radioimaging as described herein, in contrast, seeks to image radio sources in close proximity to its receivers (i.e., in a laboratory environment). Here it is shown that corona discharge, a non-destructive precursor to catastrophic (thermal) arc discharge, electromagnetically radiates strongly within a 250 kHz – 2.5 GHz bandwidth, and is readily detected and located by postprocessing the received radio signals. The ability of radioimaging to detect both corona and arc discharge (grouped together herein as high voltage breakdown or HVB) makes it a valuable tool for 100% HVB detection in materials, components, and devices, and has the ability to indicate electrical weakness (via corona detection) prior to a destructive arc discharge event. Radioimaging enables HVB to be located both internal and external to dielectric components under test in near-real-time, with multiple and/or extended HVB events located simultaneously. In contrast, existing non-destructive diagnostics (at the time of this writing) either indicate electrical breakdown without resolving failure locations (e.g., current, voltage, and chemical measurements), locate external HVB (e.g., high-speed optical and ultraviolet (UV) measurements or photography), or locate both external and internal HVB but with low fidelity (e.g., a single HVB source can be located by existing time-of-arrival (TOA) UHF or acoustic emissions). Radioimaging instead creates a sequence of high-fidelity images similar to an optical high-speed camera but at radiofrequencies (RF), and is not limited to two-dimensions. Moreover, radioimaging has already served one internal and two external industry customers, the results of which are detailed in this report. The radioimaging results described herein were part of a three-year effort funded by the Sandia Lab Directed Research and Development (LDRD) program within the Radiation, Electromagnetic, High Energy Density Science (REHEDS) investment area.

A method for estimating covariance matrices which capture the uncertainties in calculated reactor spectra has been developed. This method is based on perturbing the parameters of a physics-based analytic model fitted to a calculated spectrum. The covariance of the perturbed analytic spectra imposes energy-dependent correlations due to the physics of the neutron processes in the reactor, i.e., a fission component, a 1/E down-scatting component, and a thermal Maxwellian component. An analytic model is developed which is shown to produce good fits to several reactor environments. The covariance matrices produced via this method are then used as the prior spectrum in STAYSL least squares spectrum adjustment where it is combined with integral metrics, such as activation measurements, to produce a high-fidelity neutron spectrum characterization. It was concluded that the methodology showed agreeable results for the ACRR free-field spectrum adjustment in STAYSL resulting in a 𝜒² value of 2.21 (per degree of freedom), but further work is needed to describe scattering and interface regions.

Abstract not provided.

Abstract not provided.

This manual describes the use of the Xyce Parallel Electronic Simulator. Xyce has been designed as a SPICE-compatible, high-performance analog circuit simulator, and has been written to support the simulation needs of the Sandia National Laboratories electrical designers. This development has focused on improving capability over the current state-of-the-art in the following areas: • Capability to solve extremely large circuit problems by supporting large-scale parallel computing platforms (up to thousands of processors). This includes support for most popular parallel and serial computers. • A differential-algebraic-equation (DAE) formulation, which better isolates the device model package from solver algorithms. This allows one to develop new types of analysis without requiring the implementation of analysis-specific device models. • Device models that are specifically tailored to meet Sandia’s needs, including some radiation-aware devices (for Sandia users only). • Object-oriented code design and implementation using modern coding practices. Xyce is a parallel code in the most general sense of the phrase — a message passing parallel implementation — which allows it to run efficiently a wide range of computing platforms. These include serial, shared-memory and distributed-memory parallel platforms. Attention has been paid to the specific nature of circuit-simulation problems to ensure that optimal parallel efficiency is achieved as the number of processors grows.

Abstract not provided.

Abstract not provided.

The SNL Sierra Mechanics code suite is designed to enable simulation of complex multiphysics scenarios. The code suite is composed of several specialized applications which can operate either in standalone mode or coupled with each other. Arpeggio is a supported utility that enables loose coupling of the various Sierra Mechanics applications by providing access to Framework services that facilitate the coupling.

Abstract not provided.

Abstract not provided.

In this study, we address the challenge of enhancing image quality and spatial resolution in computed tomography (CT) imaging by introducing simulation and fabrication of high aspect ratio, point-like transmission targets. Utilizing advanced electroplating techniques, traditionally employed in the fabrication of Through Substrate Via (TSV) interconnects for CMOS circuitry, we successfully embed copper targets within silicon substrates. This method allows us to create high-aspect-ratio features specifically designed for X-ray transmission targets, resulting in micro targets that exhibit a volume increase compared to conventional evaporated surface targets. Furthermore, we present simulation results of the X-ray spectrum generated by these targets, demonstrating their potential to significantly improve both image quality and spatial resolution in CT applications. Our findings suggest that leveraging advanced fabrication techniques can open new avenues for the development of enhanced imaging technologies in medical diagnostics and beyond.

Nematic liquid crystal elastomers (LCEs) are a unique class of network polymers with the potential for enhanced mechanical energy absorption and dissipation capacity over conventional network polymers because they exhibit both conventional viscoelastic behavior and soft-elastic behavior (nematic director changes under shear loading). This additional inelastic mechanism makes them appealing as candidate damping materials in a variety of applications from vibration to impact. The lattice structures made from the LCEs provide further mechanical energy absorption and dissipation capacity associated with packing out the porosity under compressive loading. Understanding the extent of mechanical energy absorption, which is the work per unit mass (or volume) absorbed during loading, versus dissipation, which is the work per unit mass (or volume) dissipated during a loading cycle, requires measurement of both loading and unloading response. In this study, a bench-top linear actuator was employed to characterize the loading-unloading compressive response of polydomain and monodomain LCE polymers and polydomain LCE lattice structures with two different porosities (nominally, 62% and 85%) at both low and intermediate strain rates at room temperature. As a reference material, a bisphenol-A (BPA) polymer with a similar glass transition temperature (9 °C) as the nematic LCE (4 °C) was also characterized at the same conditions for comparing to the LCE polymers. Based on the loading-unloading stress-strain curves, the energy absorption and dissipation for each material at different strain rates (0.001, 0.1, 1, 10 and 90 s-1) were calculated with considerations of maximum stress and material mass/density. The strain-rate effect on the mechanical response and energy absorption and dissipation behaviors was determined. The energy dissipation ratio was also calculated from the resultant loading and unloading stress-strain curves. All five materials showed significant but different strain rate effects on energy dissipation ratio. The solid LCE and BPA materials showed greater energy dissipation capabilities at both low (0.001 s−1) and high (above 1 s−1) strain rates, but not at the strain rates in between. The polydomain LCE lattice structure showed superior energy dissipation performance compared with the solid polymers especially at high strain rates.

In remote sensing systems, the capabilities of the system are constrained by the complex interactions between size, weight, and power (SWAP) of potential designs. In electro-optical (EO) systems, examples of these critical parameters include the system’s sensitivity and resolution. Those parameters can be increased by ever larger optical apertures and focal planes but at the cost of more SWAP. Multi-image super resolution (MISR) techniques allow resolution to be enhanced via computation rather than more sophisticated optical hardware. These algorithms combine multiple images together into a single, higher resolution image, trading temporal resolution and computation for spatial resolution. Fielded MISR techniques, such as Drizzle, can require several hundred images to create a single super resolved image, implying reduced temporal resolution, increased data acquisition load, and limiting mission applications. Iterative techniques, such as model-based image reconstruction and compressive sensing, have been shown to create super resolved images using fewer images than Drizzle. They do this by posing an optimization problem that balances accuracy between a highly accurate physical model and an image model. In the case of super resolution, the physical model is defined by the relation between low resolution input images and the desired high resolution output image. The image model encodes some assumptions about the super resolved image. These assumptions are meant to suppress reconstruction artifacts that arise due to deterministic physical model error, stochastic measurement noise, and potential undersampling. In practice, the performance of iterative methods are limited by imaging models compatible with optimization. Deep learning-based methods can effectively learn image models of arbitrary complexity, but lack the theoretical explainability and robustness of iterative techniques. Consensus equilibrium (CE) generalizes the iterative techniques beyond optimization, enabling blackbox algorithms such as traditional and neural image denoisers to be used as the image model. CE-based approaches retain much of the explainability and robustness of iterative techniques while allowing the expressiveness of machine learning image models to be used. Additionally, by unrolling iterations of CE with an embedded image denoiser, the image denoiser can be further trained and specialized to the specific application with potentially higher quality reconstructions. Under this project, we demonstrated the feasibility of training an unrolled neural network based upon CE. While we didn’t train one, we showed that the CE process is differentiable and its gradient can be tractably computed. We also explored the usage of a variants of CE akin to generative neural works. Most importantly, we applied the CE framework to a number of problems including non-blind deconvolution, upsampling, single-image super resolution, MISR, event-based sensing, and saturated deconvolution. Our MISR prototype creates high quality reconstructions with an order of magnitude fewer images than previous approaches and, critically, produces these reconstructions fast enough for practical usage.

We build upon recent work on the use of machine-learning models to estimate Hamiltonian parameters using continuous weak measurement of qubits as input. We consider two settings for the training of our model: (1) supervised learning, where the weak-measurement training record can be labeled with known Hamiltonian parameters, and (2) unsupervised learning, where no labels are available. The first has the advantage of not requiring an explicit representation of the quantum state, thus potentially scaling very favorably to a larger number of qubits. The second requires the implementation of a physical model to map the Hamiltonian parameters to a measurement record, which we implement using an integrator of the physical model with a recurrent neural network to provide a model-free correction at every time step to account for small effects not captured by the physical model. We test our construction on a system of two qubits and demonstrate accurate prediction of multiple physical parameters in both the supervised context and the unsupervised context. We demonstrate that the model benefits from larger training sets, establishing that it is "learning,"and we show robustness regarding errors in the assumed physical model by achieving accurate parameter estimation in the presence of unanticipated single-particle relaxation.

Abstract not provided.