Publications

A new approach to analytically derive constitutive stress-strain relationships from modeling the work hardening behavior of alloys was developed for assessing the strength and ductility of the Ti-6Al-4V alloy. This new approach is now successfully applied for assessing the quasi-static stress-strain behavior of an additively manufactured 304L sample. A predictive capability of this modelling approach may then be extended to model material stress-strain behavior at higher strain rates of loading.

The Sound Fixing and Ranging (SOFAR) channel in the ocean allows for low frequency sound to travel thousands of kilometers, making it particularly useful for detecting underwater nuclear explosions. Suggestions that an elevated SOFAR-like channel should exist in the stratosphere date back over half a century and imply that sources within this region can be reliably sensed at vast distances. However, this theory has not been supported with evidence of direct observations from sound within this channel. Here we show that an infrasound sensor on a solar hot air balloon recorded the first infrasound detection of a ground truth airborne source while within this acoustic channel, which we refer to as the AtmoSOFAR channel. Our results support the existence of the AtmoSOFAR channel, demonstrate that acoustic signals can be recorded within it, and provide insight into the characteristics of recorded signals. Results also show a lack of detections on ground-based stations, highlighting the advantages of using balloon-borne infrasound sensors to detect impulsive sources at altitude.

A major difficulty in the analysis of molecular-level simulations is that macroscopic flow quantities are inherently noisy due to molecular fluctuations. An important example for turbulent flows is the kinetic energy dissipation rate. Traditionally, this quantity is calculated from gradients of the macroscopic velocity field, which exacerbates the noise problem. The inability to accurately compute the dissipation rate makes meaningful comparison of molecular-level and continuum simulation results a serious challenge. Herein, we extend previously developed coarse-graining theories to derive an exact molecular-level expression for the dissipation rate, which would circumvent the need to compute gradients of noisy fields. Although the exact expression cannot feasibly be implemented in Sandia’s direct simulation Monte Carlo (DSMC) code SPARTA, we utilize an approximate “hybrid” approach and compare it to the conventional gradient-based approach for planar Couette flow and the two-dimensional Taylor-Green vortex, demonstrating that the hybrid approach is significantly more accurate. Finally, we explore the possibility of adopting a Lagrangian approach to calculate the energy dissipation rate.

Abstract not provided.

This project has identified opportunities to bring further reductions in the mass and cost of modern wind turbine blades through the use of alternative material systems and manufacturing processes. The fiber reinforced polymer material systems currently used by the wind industry have stagnated as the technology continues to mature and as a means to reduce risk while introducing new products with continually increasing blade lengths. However, as blade lengths continue to increase, the challenge of controlling blade mass becomes even more critical to enabling the associated levelized cost of energy reductions. Stiffer and stronger reinforcement fibers can help to resolve the challenges of meeting the loading demands while limiting the increase in weight, but these materials are substantially more expensive than the traditional E-glass fiber systems. One goal of this project and associated work is to identify pathways that improve the cost-effectiveness of carbon fiber such that it is the reinforcement of choice in the primary structural elements of wind blades. The use of heavy-tow textile carbon fiber material systems has been shown to reduce the blade mass by 30-31% when used in the spar cap and by up to 7% when used in edgewise reinforcement. A pultrusion cost model was developed to enable a material cost comparison that includes an accurate estimate of the intermediate manufacturing step of pultrusion for the carbon fiber composite. Material cost reductions were revealed in most cases for the heavy-tow textile carbon fiber compared to infused fiberglass. The use of carbon fiber in the edgewise reinforcement produced the most notable material cost reduction of 33% for the heavy-tow textile carbon fiber. The mass and cost savings observed when using carbon fiber in edgewise reinforcement demonstrate a clear opportunity of this design approach. A carbon fiber conversion cost model was expanded to include a characterization of manufacturing costs when using advanced conversion processes with atmospheric plasma oxidation. This manufacturing approach was estimated to reduce the cost of carbon fiber material systems by greater than 10% and can be used with textile carbon systems or traditional carbon fiber precursors. The pultrusion cost model was also used to assess the opportunity for using pultruded fiberglass in wind blades, studying conventional E-glass fiber reinforcement. When using pultruded fiberglass as the spar cap material for two design classifications, the blade weight was reduced by 6% and 9% compared to infused fiberglass. However, due to the relatively large share of the pultrusion manufacturing cost compared to fiber cost, the spar cap material cost increased by 12% and 7%. When considering the system benefits of reduced blade mass and potentially lower blade manufacturing costs for pultruded composites, there may be opportunity for pultruded E-glass in wind blade spar caps, but further studies are needed. There is a clearer outcome for using pultruded fiberglass in the edgewise reinforcement where it resulted in a blade mass reduction of 2% and associated reinforcement material cost reduction of 1% compared to infused E-glass. The use of higher performing glass fibers, such as S-glass and H-glass systems, will produce greater mass savings but a study is needed to assess the cost implications for these more expensive systems. The most likely opportunity for these high-performance glass fibers is in the edgewise reinforcement, where the increased strength will reduce the damage accumulation of this fatigue-driven component. The blade design assessments in this project characterize the controlling material properties for the primary structural components in the flapwise and edgewise directions for modern wind blades. The observed trends with low and high wind speed turbine classifications for carbon and glass fiber reinforced polymer systems help to identify where cost reductions are needed, and where improvements in mechanical properties would help to reduce the material demands.

Plasma distribution in 3D space is heavily influenced by complex surfaces and the coupling interactions between plasma properties and interfacing material properties. For example, guided streamers that transition to surface ionization waves (SIWs) and propagate over structured dielectrics experience field enhancements that can lead to localized increases in ionization rates and complex 3D configurations that are difficult to analyze. Investigating these configurations requires techniques than can provide a more complete 3D picture. To help address this capability gap, a tomographic optical emission spectroscopy (tomo-OES) diagnostic system has been developed at Sandia National Laboratories that can resolve SIWs. The system includes four intensified cameras that measure the angular projections of the plasma light emission through bandpass filters. A dot calibration target co-registers each angular projection to the same voxel grid and an algebraic reconstruction technique (ART) recovers the light intensity at each voxel. An atmospheric pressure plasma jet (APPJ), provided by Peter Bruggeman, has been investigated and representative results are shown in Figure 1. Here, a bandpass filter was used to isolate emission from the N₂ second positive system (SPS) at 337.1 nm to capture the transition of the streamer to SIW on a planar dielectric surface (relative permittivity 3.3) located 3 mm below the APPJ [3]. The surface wave velocity was 3.5x10⁴ (m/s), consistent with measurements made by Steven Shannon. Characterization of this APPJ will support the group effort of standing up a reproducible APPJ across institutions for applications such as liquid treatment, catalysis, and plasma aided combustion. Future work will investigate non-planar surfaces and eventually develop tomographic laser-induced fluorescence (tomo-LIF) approaches.

This report documents the preliminary design phase of the Critical Experiment Design (CED-1) conducted as part of integral experiment request (IER) 523. The purpose of IER-523 is to determine critical configurations of 35 weight percent (wt%) enriched uranium dioxideberyllium oxide (UO₂-BeO) material with Seven Percent Critical Experiment (7uPCX) fuels at Sandia National Laboratories (Sandia). Preliminary experiment design concepts, neutronic analysis results, and proposed paths for continuing the CED process are presented. This report builds on the feasibility and justification of experimental need report (CED-0) completed in December 2021.

Abstract not provided.

Over the past decade, cybersecurity researchers have released multiple studies highlighting the insecure nature of I&C system communication protocols. In response, standards bodies have addressed the issue by adding the ability to encrypt communications to some protocols in some cases, while control system engineers have argued that encryption within these kinds of high consequence systems is in fact dangerous. Certainly, control system information between systems should be protected. But encrypting the information may not be the best way to do so. In fact, while in IT systems vendors are concerned with confidentiality, integrity, and availability, frequently in that order, in OT systems engineers are much more concerned with availability and integrity that confidentiality. In this paper, we will counter specific arguments against encrypting control system traffic, and present potential alternatives to encryption that support nuclear OT system needs more strongly that commodity IT system needs while still providing robust integrity and availability guarantees.

Improved power take-off (PTO) controller design for wave energy converters is considered a critical component for reducing the cost of energy production. However, the device and control design process often remains sequential, with the space of possible final designs largely reduced before the controller has been considered. Control co-design, whereby the device and control design are considered concurrently, has resulted in improved designs in many industries, but remains rare in the wave energy community. In this paper we demonstrate the use of a new open-source code, WecOptTool, for control co-design of wave energy converters, with the aim to make the co-design approach more accessible and accelerate its adoption. Additionally, we highlight the importance of designing a wave energy converter to maximize electrical power, rather than mechanical power, and demonstrate the co-design process while modeling the PTO's components (i.e., drive-train and generator, and their dynamics). We also consider the design and optimization of causal fixed-structure controllers. The demonstration presented here considers the PTO design problem and finds the optimal PTO drive-train that maximizes annual electrical power production. The results show a 22% improvement in the optimal controller and drive-train co-design over the optimal controller for the nominal, as built, device design.

Abstract not provided.

We describe a direct magneto-optical approach to measuring the magnetic field driven by a narrow pulse width (<10 ns), 20 kA electrical current flow in the transmission line of a high energy pulsed power accelerator. The magnetic field and electrical current are among the most important operating parameters in a pulsed power accelerator and are critical to understanding the properties of the radiation output. However, accurately measuring these fields and electrical currents using conventional pulsed power diagnostics is difficult due to the strength of ionizing radiation and electromagnetic interference. Our approach uses a fiber coupled laser beam with a rare earth element sensing crystal sensor that is highly resistant to electromagnetic interference and does not require external calibration. Here, we focus on device theory, operating parameters, results from an experiment on a high energy pulsed power accelerator, and comparison to a conventional electrical current shunt sensor.

Abstract not provided.

As a follow-up to our more comprehensive report on Adversarial Machine Learning (AML), here we provide demonstrations of AML attacks against the Limbo image database of UF6 cylinders in a variety of orientations and amongst a variety of distractor images. We demonstrate the Carlini & Wagner AML attack against a subset of Limbo images, with 100% attack success rate; meaning all attacked images were misclassified by a highly accurate trained model, yet the image changes were imperceptible to the human eye. We also demonstrate successful attacks against segmented images (images with more than one targeted object). Finally, we demonstrated the Fast Fourier Transform countermeasure that can be used to detect AML attacks on images. The intent of this and our previous report is to inform the IAEA and stakeholders of both the promise of machine learning, which could greatly improve the efficiency of surveillance monitoring, but also of the real threat of AML and potential defenses.

Recent examples provide a significant concern for the resilience of the U.S. electric grid and represent a need for enhanced decision-making to address an increasingly wide range of complex system interactions and potential consequences. In response, this LDRD project produced a proof-of-concept evaluation called the Resilience and Hazard Assessment to Prioritize Security Operations for Decisions and Impacts (RHAPSODI) methodology as an agile and flexible analytic framework capable of addressing multiple, diverse threats to desired electric grid performance. After empirically grounding needs for the future of U.S. electric grid resilience, this project employed the systems-theoretic process analysis (STPA) to develop a systems engineering risk model. The results of a completed feasibility study of a notional high voltage transmission system demonstrate an improved ability to incorporate both spatial (e.g., geographically distributed) and temporal (e.g., dynamic and time-dependent) elements of security risk to the gird. The success of this LDRD project provides the foundation for further evolution of the systems engineering risk model for the grid; derivation of quantitative approaches to evaluate risk and resilience performance; facilitation of agile experimenting and grid sensitivity to a range of vulnerabilities; and development of tools to assist decision-makers in enhancing U.S. electrical grid resilience.

Multifidelity uncertainty quantification (MF UQ) sampling approaches have been shown to significantly reduce the variance of statistical estimators while preserving the bias of the highest-fidelity model, provided that the low-fidelity models are well correlated. However, maintaining a high level of correlation can be challenging, especially when models depend on different input uncertain parameters, which drastically reduces the correlation. Existing MF UQ approaches do not adequately address this issue. In this work, we propose a new sampling strategy that exploits a shared space to improve the correlation among models with dissimilar parameterization. We achieve this by transforming the original coordinates onto an auxiliary manifold using the adaptive basis (AB) method (Tipireddy and Ghanem, 2014). The AB method has two main benefits: (1) it provides an effective tool to identify the low-dimensional manifold on which each model can be represented, and (2) it enables easy transformation of polynomial chaos representations from high- to low-dimensional spaces. This latter feature is used to identify a shared manifold among models without requiring additional evaluations. We present two algorithmic flavors of the new estimator to cover different analysis scenarios, including those with legacy and non-legacy high-fidelity (HF) data. We provide numerical results for analytical examples, a direct field acoustic test, and a finite element model of a nuclear fuel assembly. For all examples, we compare the proposed strategy against both single-fidelity and MF estimators based on the original model parameterization.

Climate and its impacts on the natural environment, and on the ability of the natural environment to support population and the built environment, stands as a threat multiplier that impacts national and global security. The Water Intersections with Climate Systems Security (WICSS) Strategic Initiative is designed to improve understanding of water’s role in, among other topics, the connection of critical infrastructure to climate in light of competing national and global security interests (including transboundary issues and stability), and identifying research gaps aligned with Sandia, and Federal agency priorities. With this impetus in mind, the WICSS Strategic Initiative team conceptualized a causal loop diagram (CLD) of the relationship between and among climate, the natural environment, population, and the built environment, with an understanding that any such regionally focused system must have externalities that influence the system from beyond its’ control, and metrics for better understanding the consequences of the set of interactions. These are discussed in light of a series of worldviews that focus on portions of the overall systems relationship. The relationships are described and documented in detail. A set of reinforcing and balancing loops are then highlighted within the context of the model. Finally, forward-looking actions are highlighted to describe how this conceptual model can be turned into modeling to address multiple problems described under the purview of the Strategic Initiative.

Abstract not provided.

Additive manufacturing of metal components enables rapid fabrication of complex geometries. However, metal additive manufacturing also introduces new morphological and microstructural characteristics which might be detrimental to component performance. Here we report the pitting corrosion properties of wrought and additively manufactured 316L stainless steel after atmospheric exposure to coastal environments and laboratory-created environments. Qualitative visualization in combination with quantitative analysis of resulting pits provided an in-depth understanding of pitting differences between wrought and additively manufactured 316L stainless steel and between coastal and laboratory-based exposure. Optical and scanning electron microscopy were utilized for visualization, while white light interferometry measured pits across approximately 5mm x 5mm areas on each sample. Post-processing of the interferometry data enables quantification of pitting attack for each sample in terms of both pit depth and pit volume. The pitting analysis introduced herein offers a new technique to compare pitting attack between different manufacturing processes and materials.

Abstract not provided.

Abstract not provided.

The ability to track the concentrations of specific molecules in the body in real time would significantly improve our ability to study, monitor, and respond to diseases. To achieve this, we require sensors that can withstand the complex environment inside the body. Electrochemical aptamer-based sensors are particularly promising for in vivo sensing, as they are among the only generalizable sensing technologies that can achieve real-time molecular monitoring directly in blood and the living body. In this project, we first focused on extending the application space of aptamer sensors to support minimally-invasive wearable measurements. To achieve this, we developed individually-addressable sensors with commercial off-the-shelf microneedles. We demonstrated sensor function in buffer, blood, and porcine skin (a common proxy for human skin). In addition to the applied sensing project, we also worked to improve fundamental understanding of the aptamer sensing platform and how it responds to biomolecular interferents. Specifically, we explored the interfacial dynamics of biofouling – a process impacting sensors placed in complex fluids, such as blood.

Abstract not provided.

Abstract not provided.

Large rocket motors may violently explode when exposed to accidental fires. Even hot metal fragments from a nearby accident may penetrate the propellant and ultimately cause thermal ignition. A mechanistic understanding of heated propellants leading to thermal runaway is a major unsolved problem. Here we show that thermal ignition in propellants can be predicted using a universal cookoff model coupled to a micromechanics pressurization model. Our model predicts the time to thermal ignition in cookoff experiments with variable headspace volumes. We found that experiments with headspace volumes are more prone to deformation which distorts pores and causes increased permeability when the propellant expands into this headspace. Delayed ignition with larger headspace volume correlates with lower headspace pressures during decomposition. We found that our predictions matched experimental measurements best when the initial propellant was impermeable to gas flow rather than being permeable. Similar behavior is expected with other energetic materials with rubbery binders. Our model is validated using data from a separate laboratory. We also present an uncertainty analysis using Latin Hypercube Sampling (LHS) of thermal ignition caused by a steel fragment embedded in the propellant.

Current approaches to securing high consequence facilities (HCF) and critical assets are linear and static and therefore struggle to adapt to emerging threats (e.g., unmanned aerial systems) and changing environmental conditions (e.g., decreasing operational control). The pace of change in technological, organizational, societal, and political dynamics necessitates a move toward codifying underlying scientific principles to better characterize the rich interactions observed between HCF security technology, infrastructure, digital assets, and human or organizational components. The promising results of Laboratory Directed Research and Development (LDRD) 20-0373—“Developing a Resilient, Adaptive, and Systematic Paradigm for Security Analysis”—suggest that when compared to traditional security analysis, invoking multilayer network (MLN) modeling for HCF security system components captures unexpected failure cases and unanticipated interactions.

Abstract not provided.

There is an increasing aspiration to utilize machine learning (ML) for various tasks of relevance to national security. ML models have thus far been mostly applied to tasks and domains that, while impactful, have sufficient volume of data. For predictive tasks of national security relevance, ML models of great capacity (ability to approximate nonlinear trends in input-output maps) are often needed to capture the complex underlying physics. However, scientific problems of relevance to national security are often accompanied by various sources of sparse and/or incomplete data, including experiments and simulations, across different regimes of operation, of varying degrees of fidelity, and include noise with different characteristics and/or intensity. State-of-the-art ML models, despite exhibiting superior performance on the task and domain they were trained on, may suffer detrimental loss in performance in such sparse data environments. This report summarizes the results of the Laboratory Directed Research and Development project entitled Trust-Enhancing Probabilistic Transfer Learning for Sparse and Noisy Data Environments. The objective of the project was to develop a new transfer learning (TL) framework that aims to adaptively blend the data across different sources in tackling one task of interest, resulting in enhanced trustworthiness of ML models for mission- and safety-critical systems. The proposed framework determines when it is worth applying TL and how much knowledge is to be transferred, despite uncontrollable uncertainties. The framework accomplishes this by leveraging concepts and techniques from the fields of Bayesian inverse modeling and uncertainty quantification, relying on strong mathematical foundations of probability and measure theories to devise new uncertainty-aware TL workflows.

Abstract not provided.

Accurate event locations are important for many endeavors in seismology, and understanding the factors that contribute to uncertainties in those locations is complex. In this article, we present a case study that takes an in-depth look at the accuracy and precision possible for locating nine shallow earthquakes in the Rock Valley fault zone in southern Nevada. These events are targeted by the Rock Valley Direct Comparison phase of the Source Physics Experiment, as candidates for the colocation of a chemical explosion with an earthquake hypocenter to directly compare earthquake and explosion sources. For this comparison, it is necessary to determine earthquake hypocenters as accurately as possible so that different source types have nearly identical locations. Our investigations include uncertainty analysis from different sets of phase arrivals, stations, velocity models, and location algorithms. For a common set of phase arrivals and stations, we find that epicentral locations from different combinations of velocity models and algorithms are within 600 m of one another in most cases. Event depths exhibit greater uncertainties, but focusing on the S-P times at the nearest station allows for estimates within approximately 500 m.

Abstract not provided.

Abstract not provided.

This report presents our work to model the workloads of a linear electromagnetic application based on the method of moments in the frequency domain to effectively load balance the matrix assembly. This application is particularly challenging to load balance due to its lack of persistent iterative behavior, its operation under tight memory constraint (where the matrix may fill 80% of memory on each node), and the algorithmic complexity of the computational method. This report describes the first step in our work to apply an inspector-executor approach for load balancing workloads where key parameters are exposed during the inspector phase and a pre-trained model is applied to predict relative task weights for the load balancer.

Abstract not provided.

The International Database of Reference Gamma-Ray Spectra of Various Nuclear Matter is designed to hold curated gamma spectral data and will be hosted by the International Atomic Energy Agency on its public facing web site. Currently, the database to be hosted is given to the International Atomic Energy Agency by Sandia. This document describes the application used by Sandia to load spectral data into a database.

The preliminary use of the in-situ curvature measurement technique for analyzing the planar stress evolution of controlled atmosphere plasma spray (CAPS) refractory metal deposits was performed with SNL-NM org. 1834’s CAPS system. A porous refractory metal exemplar of Ta-Nb was sprayed onto Ni-200, Ti-6Al-4V, and Al 7075-T6 substrates using a constant plasma torch parameter setting and deposition toolpath. Residual stresses of the deposits were found to be largely influenced by the substrate coefficient of thermal expansion and were calculated to be 49, 90, and -136 MPa for Ni 200, Ti-6Al-4V, and Al 7075-T6, respectively. The “Evolving stress” of the Ta-Nb deposits, which more accurately describes the mean intrinsic splat quenching stress of the spray material during deposition, was calculated to be 67, 92, and 129 MPa for Ni-200, Ti-6Al-4V, and Al 7075-T6, respectively. Notable difference in curvature measurement for the 1st coating pass for the Al 7075- T6 substrate was observed, with interface micrograph evidence suggesting potential softening and/or melting of the Al 7075-T6 substrate surface during deposition. Substrate temperature measurements prior to Ta-Nb deposition were used to calculate thermal energy absorbed from the hot gas plume by the different substrates and were found to correlate to the substrate’s thermal effusively. These calculated thermal energies were also found to be ~10 to 15% of the calculated energy output from the plasma torch’s nozzle exit for these experimental conditions.

Spectrally resolved signals in the short- to mid-wave infrared (SWIR/MWIR) bands at high-temporal resolution are critical for many national security remote sensing missions. Currently available off the shelf technology can achieve either high temporal resolution or high spectral resolution, but rugged instruments that can achieve both simultaneously remain mostly in the realm of one-off R&D projects. This report documents efforts to demonstrate a new technique for designing and building high resolution, high framerate multichannel FTIR (MC-FTIR) spectrometers that operate in the SWIR/MWIR bands. The core optical element in a MC-FTIR spectrometer is an array of statically-tuned lamellar grating interferometers (LGI). In the original MC-FTIR work these arrays were fabricated using a synchrotron x-ray lithography method. We proposed to instead fabricate these LGI arrays using multiphoton lithography (MPL), a 3D printing technique that can fabricate meso-scale structures with sub-micron precision. Although we were able to fabricate LGI arrays of sufficient size using MPL, the realized optical surfaces had unsuitably high optical form errors, precluding their use in a fieldable instrument. Further advancement in MPL technology may eventually enable fabrication of interferometer-grade LGI arrays.

Abstract not provided.

Tritium exhibits unique environmental behavior because of its potential interactions with water and organic substances. Modeling the environmental consequences of tritium releases can be relatively complex and thus an evaluation of MACCS is needed to understand what updates, if any, are needed in MACCS to account for the behavior of tritium. We examine documented tritium releases and previous benchmarking assessments to perform a model intercomparison between MACCS and state-of-practice tritium-specific codes UFOTRI and ETMOD to quantify the difference between MACCS and state of practice models for assessing tritium consequences. Additionally, information to assist an analyst in judging whether a postulated tritium release is likely to lead to significant doses is provided.

Freeplay is a common type of piecewise-smooth nonlinearity in dynamical systems, and it can cause discontinuity-induced bifurcations and other behaviors that may bring about undesirable and potentially damaging responses. Prior research has focused on piecewise-smooth systems with two or three distinct regions, but less attention is devoted to systems with more regions (i.e., multi-segmented systems). In this work, numerical analysis is performed on a dynamical system with multi-segmented freeplay, in which there are four stiffness transitions and five distinct regions in the phase space. The effects of the multi-segmented parameters are studied through bifurcation diagram evolution along with induced multi-stable behavior and different bifurcations. These phenomena are interrogated through various tools, such as harmonic balance, basins of attraction, phase planes, and Poincaré section analysis. Results show that among the three multi-segmented parameters, the asymmetry has the strongest effect on the response of the system.

Nonlinear modeling and optimization is a valuable tool for aiding decisions by engineering practitioners, but programming an optimization problem based on a complex electrical, mechanical, or chemical process is a time-consuming and error-prone activity. Therefore, there is a need for model analysis and debugging tools that can detect and diagnose modeling errors. One such tool is the Dulmage–Mendelsohn decomposition, which identifies structurally under- and over-determined subsets in systems of equations and variables by partitioning the bipartite graph of the system. This work provides the necessary background to understand the Dulmage–Mendelsohn decomposition and its application to the analysis of nonlinear optimization problems, demonstrates its use in diagnosing a variety of modeling errors, and introduces software implementations for analyzing nonlinear optimization problems in the Pyomo and JuMP algebraic modeling languages.

The “Pioneer WEC” project is targeted at developing a wave energy generator for the Coastal Surface Mooring (CSM) system within the Ocean Observatories Initiative (OOI) Pioneer Array. The CSM utilizes solar photovoltaic and wind generation systems, along with rechargeable batteries, to power multiple sensors on the buoy and along the mooring line. This approach provides continuous power for essential controller functions and a subset of instruments, and meets the full power demand roughly 70% of the time. Sandia has been tasked with designing a wave energy system to provide additional electrical power and bring the CSM up-time for satisfying the full-power demand to 100%. This project is a collaboration between Sandia and Woods Hole Oceanographic Institution (WHOI), along with Evergreen Innovations, Monterey Bay Aquarium Research Institute (MBARI), Eastern Carolina University (ECU), Johns Hopkins University (JHU), and the National Renewable Energy Laboratory (NREL). This report captures Phase I of an expected two phase project and presents project scoping and concept design results. phase project and presents project scoping and concept design results.

This project developed a novel statistical understanding of compression analytics (CA), which has challenged and clarified some core assumptions about CA, and enabled the development of novel techniques that address vital challenges of national security. Specifically, this project has yielded the development of novel capabilities including 1. Principled metrics for model selection in CA, 2. Techniques for deriving/applying optimal classification rules and decision theory to supervised CA, including how to properly handle class imbalance and differing costs of misclassification, 3. Two techniques for handling nonlocal information in CA, 4. A novel technique for unsupervised CA that is agnostic with regard to the underlying compression algorithm, 5. A framework for semisupervised CA when a small number of labels are known in an otherwise large unlabeled dataset. 6. The academic alliance component of this project has focused on the development of a novel exemplar-based Bayesian technique for estimating variable length Markov models (closely related to PPM [prediction by partial matching] compression techniques). We have developed examples illustrating the application of our work to text, video, genetic sequences, and unstructured cybersecurity log files.

Abstract not provided.

Machine learning models are promising as surrogates in optimization when replacing difficult to solve equations or black-box type models. This work demonstrates the viability of linear model decision trees as piecewise-linear surrogates in decision-making problems. Linear model decision trees can be represented exactly in mixed-integer linear programming (MILP) and mixed-integer quadratic constrained programming (MIQCP) formulations. Furthermore, they can represent discontinuous functions, bringing advantages over neural networks in some cases. We present several formulations using transformations from Generalized Disjunctive Programming (GDP) formulations and modifications of MILP formulations for gradient boosted decision trees (GBDT). We then compare the computational performance of these different MILP and MIQCP representations in an optimization problem and illustrate their use on engineering applications. We observe faster solution times for optimization problems with linear model decision tree surrogates when compared with GBDT surrogates using the Optimization and Machine Learning Toolkit (OMLT).

This is a SAND Report on cross-correlation data collected at the Redmond Salt Mine. It discusses methods, as well as temporal variability and energy characteristics of the cross-correlation data.

Abstract not provided.

In this report, we investigate the effects of conductor losses in a multipole-based cable braid magnetic penetration model. Our multipole model uses a mesh of the actual cable geometry, which enables us to model more complicated structures. After summarizing the first principles model formulation, we consider a one-dimensional array of wires, for which an analytical solution is known in the lossless case. We extend this solution to the lossy case by using a complex-valued radius. We also model this structure analytically using a conformal-mapping solution. We then compare both the self-impedance and the transfer impedance results from our first principles cable braid electromagnetic penetration model to those obtained using the analytical solutions. An analysis for various frequencies (and skin depths) usually encountered in cable modeling is reported. These results are found in good agreement up to a radius to half spacing ratio of about 0.7, demonstrating a robustness needed for many commercial and non-commercial cables.

Abstract not provided.

Accurate targeting of radioisotope classifiers and estimators requires an understanding of the target problem space. In order to facilitate clear communication on expected model behavior and performance between practitioners and stakeholders on their problems, this questionnaire was created. Stakeholder responses form the basis of a trained model as well as the start of usage requirements for the model as it is integrated with analysis processes or detection systems. This questionnaire may also be useful to machine learning practitioners and gamma spectroscopists developing new algorithms as a starting point for characterizing their problem space, especially if they are using PyRIID.

Abstract not provided.

Abstract not provided.

In this paper, an approach for 3D plasma structure diagnostics using tomographic optical emission spectroscopy (Tomo-OES) of a nanosecond pulsed atmospheric pressure plasma jet (APPJ) is presented. In contrast to the well-known Abel inversion, Tomo-OES does not require cylindrical symmetry to recover 3D distributions of plasma light emission. Instead, many 2D angular projections are measured with intensified cameras and the multiplicative algebraic reconstruction technique is used to recover the 3D distribution of light emission. This approach solves the line-of-sight integration problem inherent to optical diagnostics, allowing recovery of localized OES information within the plasma that can be used to better infer plasma parameters within complex plasma structures. Here, Tomo-OES was applied to investigate an APPJ operated with helium in ambient air and impinging on planar and structured dielectric surfaces. Surface charging caused the guided streamer from the APPJ to transition to a surface ionization wave (SIW) that propagated along the surface. The SIW experienced variable geometrical and electrical material properties as it propagated, leading to 3D configurations that were non-symmetric and spatially complex. Light emission from He, N 2 + , and N2 were imaged at ten angular projections and the respective time-resolved 3D emission distributions in the plasma were then reconstructed. The spatial resolution of each tomographic reconstruction was 7.4 µm and the temporal resolution was 5 ns, sufficient to observe the guided streamer and the effects of the structured surface on the SIW. Emission from He showed the core of the jet and emission from N 2 + and N2 indicated effects of entrainment of ambient air. Penning ionization of N2 created a ring or outer layer of N 2 + that spatially converged to form the ‘plasma bullet’ or spatially diverged across a surface as part of a SIW. The SIW entered trenches of size 150 µm, leading to decreases in plasma light emission in regions above the trenches. The plasma light emission was higher in some regions with trenches, possibly due to effects of field enhancement.

Numerical simulations are used to study the dynamics of a developing suspension Poiseuille flow with monodispersed and bidispersed neutrally buoyant particles in a planar channel, and machine learning is applied to learn the evolving stresses of the developing suspension. The particle stresses and pressure develop on a slower time scale than the volume fraction, indicating that once the particles reach a steady volume fraction profile, they rearrange to minimize the contact pressure on each particle. We consider the timescale for stress development and how the stress development connects to particle migration. For developing monodisperse suspensions, we present a new physics-informed Galerkin neural network that allows for learning the particle stresses when direct measurements are not possible. We show that when a training set of stress measurements is available, the MOR-physics operator learning method can also capture the particle stresses accurately.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Static structure factors are computed for large-scale, mechanically stable, jammed packings of frictionless spheres (three dimensions) and disks (two dimensions) with broad, power-law size dispersity characterized by the exponent -β. The static structure factor exhibits diverging power-law behavior for small wave numbers, allowing us to identify a structural fractal dimension df. In three dimensions, df≈2.0 for 2.5≤β≤3.8, such that each of the structure factors can be collapsed onto a universal curve. In two dimensions, we instead find 1.0df1.34 for 2.1≤β≤2.9. Furthermore, we show that the fractal behavior persists when rattler particles are removed, indicating that the long-wavelength structural properties of the packings are controlled by the large particle backbone conferring mechanical rigidity to the system. A numerical scheme for computing structure factors for triclinic unit cells is presented and employed to analyze the jammed packings.

Strongly charged polyelectrolytes (PEs) demonstrate complex solution behavior as a function of chain length, concentrations, and ionic strength. The viscosity behavior is important to understand and is a core quantity for many applications, but aspects remain a challenge. Molecular dynamics simulations using implicit solvent coarse-grained (CG) models successfully reproduce structure, but are often inappropriate for calculating viscosities. To address the need for CG models which reproduce viscoelastic properties of one of the most studied PEs, sodium polystyrene sulfonate (NaPSS), we report our recent efforts in using Bayesian optimization to develop CG models of NaPSS which capture both polymer structure and dynamics in aqueous solutions with explicit solvent. We demonstrate that our explicit solvent CG NaPSS model with the ML-BOP water model [Chan et al. Nat Commun 10, 379 (2019)] quantitatively reproduces NaPSS chain statistics and solution structure. The new explicit solvent CG model is benchmarked against diffusivities from atomistic simulations and experimental specific viscosities for short chains. We also show that our Bayesian-optimized CG model is transferable to larger chain lengths across a range of concentrations. Overall, this work provides a machine-learned model to probe the structural, dynamic, and rheological properties of polyelectrolytes such as NaPSS and aids in the design of novel, strongly charged polymers with tunable structural and viscoelastic properties

Abstract not provided.

A collection of MATLAB functions and class definitions called System Workflow Tools (SWFT) are available to semi-automate steps in the simulation process. Some of these steps are often simple and routine for smaller finite element models, but if done directly by an analyst can quickly become labor intensive, cumbersome, and error prone for larger, system level models. Some of SWFT’s capabilities demonstrated in this report includes writing Sierra input decks and processing Quantities of Interest (QOI) from results files. SWFT also writes scripts in order to utilize other software programs such as Cubit (separating system level CAD into subassemblies and components, creating nodesets and sidesets), DAKOTA (ensemble management), and ParaView (contour plots and animations). Detailed commands and workflows from mesh generation to report generation are provided as examples for analysts to utilize SWFT capabilities.

As Machine Learning (ML) continues to advance, it is being integrated into more systems. Often, the ML component represents a significant portion of the system that reduces the burden on the end user or significantly improves task performance. However, the ML component represents an unknown complex phenomenon that is learned from collected data without the need to be explicitly programmed. Despite the improvement in task performance, the models are often black boxes. Evaluating the credibility and the vulnerabilities of ML models poses a gap in current test and evaluation practice. For high consequence applications, the lack of testing and evaluation procedures represents a significant source of uncertainty and risk. To help reduce that risk, here we present considerations to evaluate systems embedded with an ML component within a red-teaming inspired methodology. We focus on (1) cyber vulnerabilities to an ML model, (2) evaluating performance gaps, and (3) adversarial ML vulnerabilities.

The benefits of high-performance unidirectional carbon fiber composites are limited in many cost-driven industries due to the high cost relative to alternative reinforcement fibers. Low-cost carbon fibers have been previously proposed, but the longitudinal compressive strength continues to be a limiting factor or studies are based on simplifications that warrant further analysis. A micromechanical model is used to (1) determine if the longitudinal compressive strength of composites can be improved with noncircular carbon fiber shapes and (2) characterize why some shapes are stronger than others in compression. In comparison to circular fibers, the results suggest that the strength can be increased by 10%–13% by using a specific six-lobe fiber shape and by 6%–9% for a three-lobe fiber shape. A slight increase is predicted in the compressive strength of the study two-lobe fiber but has the highest uncertainty and sensitivity to fiber orientation and misalignment direction. The underlying mechanism governing the compressive failure of the composites was linked to the unique stress fields created by the lobes, particularly the pressure stress in the matrix. This work provides mechanics-based evidence of strength improvements from noncircular fiber shapes and insight on how matrix yielding is altered with alternative fiber shapes.

The development of additively-manufactured (AM) 316L stainless steel (SS) using laser powder bed fusion (LPBF) has enabled near net shape components from a corrosion-resistant structural material. In this article, we present a multiscale study on the effects of processing parameters on the corrosion behavior of as-printed surfaces of AM 316L SS formed via LPBF. Laser power and scan speed of the LPBF process were varied across the instrument range known to produce parts with >99 % density, and the macroscale corrosion trends were interpreted via microscale and nanoscale measurements of porosity, roughness, microstructure, and chemistry. Porosity and roughness data showed that porosity φ decreased as volumetric energy density Ev increased due to a shift in the pore formation mechanism and that roughness Sq was due to melt track morphology and partially fused powder features. Cross-sectional and plan-view maps of chemistry and work function ϕs revealed an amorphous Mn-silicate phase enriched with Cr and Al that varied in both thickness and density depending on Ev. Finally, the macroscale potentiodynamic polarization experiments under full immersion in quiescent 0.6 M NaCl showed significant differences in breakdown potential Eb and metastable pitting. In general, samples with smaller φ and Sq values and larger ϕs values and homogeneity in the Mn-silicate exhibited larger Eb. The porosity and roughness effects stemmed from an increase to the overall number of initiation sites for pitting, and the oxide phase contributed to passive film breakdown by acting as a crevice former or creating a galvanic couple with the SS.

The United States Department of Energy’s (DOE) Office of Nuclear Energy’s Spent Fuel and Waste Science and Technology Campaign seeks to better understand the technical basis, risks, and uncertainty associated with the safe and secure disposition of spent nuclear fuel (SNF) and high-level radioactive waste. Commercial nuclear power generation in the United States has resulted in thousands of metric tons of SNF, the disposal of which is the responsibility of the DOE (Nuclear Waste Policy Act of 1982, as amended). Any repository licensed to dispose of SNF must meet requirements regarding the long-term performance of that repository. For an evaluation of the long-term performance of the repository, one of the events that may need to be considered is the SNF achieving a critical configuration during the postclosure period. Of particular interest is the potential behavior of SNF in dual-purpose canisters (DPCs), which are currently licensed and being used to store and transport SNF but were not designed for permanent geologic disposal.

Characterizing interface trap states in commercial wide bandgap devices using frequency-based measurements requires unconventionally high probing frequencies to account for both fast and slow traps associated with wide bandgap materials. The C − ψ S technique has been suggested as a viable quasi-static method for determining the interface trap state densities in wide bandgap systems, but the results are shown to be susceptible to errors in the analysis procedure. This work explores the primary sources of errors present in the C − ψ S technique using an analytical model that describes the apparent response for wide bandgap MOS capacitor devices. Measurement noise is shown to greatly impact the linear fitting routine of the 1 / C S ∗ 2 vs ψ S plot to calibrate the additive constant in the surface potential/gate voltage relationship, and an inexact knowledge of the oxide capacitance is also shown to impede interface trap state analysis near the band edge. In addition, a slight nonlinearity that is typically present throughout the 1 / C S ∗ 2 vs ψ S plot hinders the accurate estimation of interface trap densities, which is demonstrated for a fabricated n-SiC MOS capacitor device. Methods are suggested to improve quasi-static analysis, including a novel method to determine an approximate integration constant without relying on a linear fitting routine.

Magnetic properties of more than 20 Cantor alloy samples of varying composition were investigated over a temperature range of 5 K to 300 K and in fields of up to 70 kOe using magnetometry and muon spin relaxation. Two transitions are identified: a spin-glass-like transition that appears between 55K and 190K, depending on composition, and a ferrimagnetic transition that occurs at approximately 43K in multiple samples with widely varying compositions. The magnetic signatures at 43K are remarkably insensitive to chemical composition. A modified Curie-Weiss model was used to fit the susceptibility data and to extract the net effective magnetic moment for each sample. The resulting values for the net effective moment were either diminished with increasing Cr or Mn concentrations or enhanced with decreasing Fe, Co, or Ni concentrations. Beyond a sufficiently large effective moment, the magnetic ground state transitions from ferrimagnetism to ferromagnetism. The effective magnetic moments, together with the corresponding compositions, are used in a global linear regression analysis to extract element-specific effective magnetic moments, which are compared to the values obtained by ab initio based density functional theory calculations. Finally, these moments provide the information necessary to controllably tune the magnetic properties of Cantor alloy variants.

Systems engineering today faces a wide array of challenges, ranging from new operational environments to disruptive technological — necessitating approaches to improve research and development (R&D) efforts. Yet, emphasizing the Aristotelian argument that the “whole is greater than the sum of its parts” seems to offer a conceptual foundation creating new R&D solutions. Invoking systems theoretic concepts of emergence and hierarchy and analytic characteristics of traceability, rigor, and comprehensiveness is potentially beneficial for guiding R&D strategy and development to bridge the gap between theoretical problem spaces and engineering-based solutions. In response, this article describes systems–theoretic process analysis (STPA) as an example of one such approach to aid in early-systems R&D discussions. STPA—a ‘top-down’ process that abstracts real complex system operations into hierarchical control structures, functional control loops, and control actions—uses control loop logic to analyze how control actions (designed for desired system behaviors) may become violated and drive the complex system toward states of higher risk. By analyzing how needed controls are not provided (or out of sequence or stopped too soon) and unneeded controls are provided (or engaged too long), STPA can help early-system R&D discussions by exploring how requirements and desired actions interact to either mitigate or potentially increase states of risk that can lead to unacceptable losses. This article will demonstrate STPA's benefit for early-system R&D strategy and development discussion by describing such diverse use cases as cyber security, nuclear fuel transportation, and US electric grid performance. Together, the traceability, rigor, and comprehensiveness of STPA serve as useful tools for improving R&D strategy and development discussions. In conclusion, leveraging STPA as well as related systems engineering techniques can be helpful in early R&D planning and strategy development to better triangulate deeper theoretical meaning or evaluate empirical results to better inform systems engineering solutions.

Abstract not provided.

Process variations within Field Programmable Gate Arrays (FPGAs) provide a rich source of entropy and are therefore well-suited for the implementation of Physical Unclonable Functions (PUFs). However, careful considerations must be given to the design of the PUF architecture as a means of avoiding undesirable localized bias effects that adversely impact randomness, an important statistical quality characteristic of a PUF. Here in this paper, we investigate a ring-oscillator (RO) PUF that leverages localized entropy from individual look-up table (LUT) primitives. A novel RO construction is presented that enables the individual paths through the LUT primitive to be measured and isolated at high precision, and an analysis is presented that demonstrates significant levels of localized design bias. The analysis demonstrates that delay-based PUFs that utilize LUTs as a source of entropy should avoid using FPGA primitives that are localized to specific regions of the FPGA, and instead, a more robust PUF architecture can be constructed by distributing path delay components over a wider region of the FPGA fabric. Compact RO PUF architectures that utilize multiple configurations within a small group of LUTs are particularly susceptible to these types of design-level bias effects. The analysis is carried out on data collected from a set of identically designed, hard macro instantiations of the RO implemented on 30 copies of a Zynq 7010 SoC.

For reactive burn models in hydrocodes, an equilibrium closure assumption is typically made between the unreacted and product equations of state. In the CTH [1] (not an acronym) hydrocode the assumption of density and temperature equilibrium is made by default, while other codes make a pressure and temperature equilibrium assumption. The main reason for this difference is the computational efficiency in making the density and temperature assumption over the pressure and temperature one. With fitting to data, both assumptions can accurately predict reactive flow response using the various models, but the model parameters from one code cannot necessarily be used directly in a different code with a different closure assumption. A new framework is intro-duced in CTH to allow this assumption to be changed independently for each reactive material. Comparisons of the response and computational cost of the History Variable Reactive Burn (HVRB) reactive flow model with the different equilibrium assumptions are presented.

Maximizing the production of heterologous biomolecules is a complex problem that can be addressed with a systems-level understanding of cellular metabolism and regulation. Specifically, growth-coupling approaches can increase product titers and yields and also enhance production rates. However, implementing these methods for non-canonical carbon streams is challenging due to gaps in metabolic models. Over four design-build-test-learn cycles, we rewire Pseudomonas putida KT2440 for growth-coupled production of indigoidine from para-coumarate. We explore 4,114 potential growth-coupling solutions and refine one design through laboratory evolution and ensemble data-driven methods. The final growth-coupled strain produces 7.3 g/L indigoidine at 77% maximum theoretical yield in para-coumarate minimal medium. The iterative use of growth-coupling designs and functional genomics with experimental validation was highly effective and agnostic to specific hosts, carbon streams, and final products and thus generalizable across many systems.

A new capability for modeling graded density reactive flow materials in the shock physics hydrocode, CTH, is demonstrated here. Previously, materials could be inserted in CTH with graded material properties, but the sensitivity of the material was not adjusted based on these properties. Of particular interest are materials that are graded in density, sometimes due to pressing or other assembly operations. The sensitivity of explosives to both density and temperature has been well demonstrated in the literature, but to-date the material parameters for use in a simulation were fit to a single condition and applied to the entire material, or the material had to be inserted in sections and each section assigned a condition. The reactive flow model xHVRB has been extended to shift explosive sensitivity with initial density, so that sensitivity is also graded in the material. This capability is demonstrated for use in three examples. The first models detonation transfer in a graded density pellet of HNS, the second is a shaped charge with density gradients in the explosive, and the third is an explosively formed projectile.

Organic co-crystals have emerged as a promising class of semiconductors for next-generation optoelectronic devices due to their unique photophysical properties. This paper presents a joint experimental-theoretical study comparing the crystal structure, spectroscopy, and electronic structure of two charge transfer co-crystals. Reported herein is a novel co-crystal Npe:TCNQ, formed from 4-(1-naphthylvinyl)pyridine (Npe) and 7,7,8,8-tetracyanoquinodimethane (TCNQ) via molecular self-assembly. This work also presents a revised study of the co-crystal composed of Npe and 1,2,4,5-tetracyanobenzene (TCNB) molecules, Npe:TCNB, herein reported with a higher-symmetry (monoclinic) crystal structure than previously published. Npe:TCNB and Npe:TCNQ dimer clusters are used as theoretical model systems for the co-crystals; the geometries of the dimers are compared to geometries of the extended solids, which are computed with periodic boundary conditions density functional theory. UV-Vis absorption spectra of the dimers are computed with time-dependent density functional theory and compared to experimental UV-Vis diffuse reflectance spectra. Both Npe:TCNB and Npe:TCNQ are found to exhibit neutral character in the S0 state and ionic character in the S1 state. The high degree of charge transfer in the S1 state of both Npe:TCNB and Npe:TCNQ is rationalized by analyzing the changes in orbital localization associated with the S1 transitions.

Accelerators that drive z-pinch experiments transport current densities in excess of 1 MA/cm² in order to melt or ionize the target and implode it on axis. These high current densities stress the transmission lines upstream from the target, where rapid electrode heating causes plasma formation, melt, and possibly vaporization. These plasmas negatively impact accelerator efficiency by diverting some portion of the current away from the target, referred to as “current loss”. Simulations that are able to reproduce this behavior may be applied to improving the efficiency of existing accelerators and to designing systems operating at ever higher current densities. The relativistic particle-in-cell code CHICAGO^® is the primary code for modeling power flow on Sandia National Laboratories’ Z accelerator. We report here on new algorithms that incorporate vaporization and melt into the standard power-flow simulation framework. Taking a hybrid approach, the CHICAGO® kinetic/multi-fluid treatment has been expanded to include vaporization while the quasi-neutral equation-of-motion has been updated for melt at high current-densities. For vaporization, a new one-dimensional substrate model provides a more accurate calculation of electrode thermal, mass, and magnetic field diffusion as well as a means of emitting absorbed contaminants and vaporized metal ions. A quasi-fluid model has been implemented expressly to mimic the motion of imploding liners for accurate inductance histories. For melt, a multi-ion Hall-MHD option has been implemented and benchmarked against Alegra MHD. This new model is described with sufficient detail to reproduce these algorithms in any hybrid kinetic code. Physics results from the new code are also presented. A CHICAGO^® Hall-MHD simulation of a radial transmission line demonstrates that Hall physics, not included in Alegra, has no significant impact on the diffusion of electrode material. When surface contaminant desorption is mocked in as a hydrogen surface plasma, both the surface and bulk-material plasmas largely compress under the influence of the j × B force. Similar results are seen in Alegra, which also shows magnetic and material diffusion scaling with peak current. Test vaporization simulations using MagLIF and a power-flow experimental geometry show Fe⁺ ions diffuse only a few hundred µm from the electrodes, so present models of Z power flow remain valid.

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 weapons systems. This document provides a user’s guide to the input for Sierra/SD. Details of input specifications for the different solution types, output options, element types and parameters are included. The appendices contain detailed examples, and instructions for running the software on parallel platforms.

Predictive modeling typically relies on Bayesian model calibration to provide uncertainty quantification. Variational inference utilizing fully independent (“mean-field”) Gaussian distributions are often used as approximate probability density functions. This simplification is attractive since the number of variational parameters grows only linearly with the number of unknown model parameters. However, the resulting diagonal covariance structure and unimodal behavior can be too restrictive to provide useful approximations of intractable Bayesian posteriors that exhibit highly non-Gaussian behavior, including multimodality. High-fidelity surrogate posteriors for these problems can be obtained by considering the family of Gaussian mixtures. Gaussian mixtures are capable of capturing multiple modes and approximating any distribution to an arbitrary degree of accuracy, while maintaining some analytical tractability. Unfortunately, variational inference using Gaussian mixtures with full-covariance structures suffers from a quadratic growth in variational parameters with the number of model parameters. The existence of multiple local minima due to strong nonconvex trends in the loss functions often associated with variational inference present additional complications, These challenges motivate the need for robust initialization procedures to improve the performance and computational scalability of variational inference with mixture models. In this work, we propose a method for constructing an initial Gaussian mixture model approximation that can be used to warm-start the iterative solvers for variational inference. The procedure begins with a global optimization stage in model parameter space. In this step, local gradient-based optimization, globalized through multistart, is used to determine a set of local maxima, which we take to approximate the mixture component centers. Around each mode, a local Gaussian approximation is constructed via the Laplace approximation. Finally, the mixture weights are determined through constrained least squares regression. The robustness and scalability of the proposed methodology is demonstrated through application to an ensemble of synthetic tests using high-dimensional, multimodal probability density functions. Here, the practical aspects of the approach are demonstrated with inversion problems in structural dynamics.

For the first time the optimal local truncation error method (OLTEM) with 125-point stencils and unfitted Cartesian meshes has been developed in the general 3-D case for the Poisson equation for heterogeneous materials with smooth irregular interfaces. The 125-point stencils equations that are similar to those for quadratic finite elements are used for OLTEM. The interface conditions for OLTEM are imposed as constraints at a small number of interface points and do not require the introduction of additional unknowns, i.e., the sparse structure of global discrete equations of OLTEM is the same for homogeneous and heterogeneous materials. The stencils coefficients of OLTEM are calculated by the minimization of the local truncation error of the stencil equations. These derivations include the use of the Poisson equation for the relationship between the different spatial derivatives. Such a procedure provides the maximum possible accuracy of the discrete equations of OLTEM. In contrast to known numerical techniques with quadratic elements and third order of accuracy on conforming and unfitted meshes, OLTEM with the 125-point stencils provides 11-th order of accuracy, i.e., an extremely large increase in accuracy by 8 orders for similar stencils. The numerical results show that OLTEM yields much more accurate results than high-order finite elements with much wider stencils. The increased numerical accuracy of OLTEM leads to an extremely large increase in computational efficiency. Additionally, a new post-processing procedure with the 125-point stencil has been developed for the calculation of the spatial derivatives of the primary function. The post-processing procedure includes the minimization of the local truncation error and the use of the Poisson equation. It is demonstrated that the use of the partial differential equation (PDE) for the 125-point stencils improves the accuracy of the spatial derivatives by 6 orders compared to post-processing without the use of PDE as in existing numerical techniques. At an accuracy of 0.1% for the spatial derivatives, OLTEM reduces the number of degrees of freedom by 900 - 4∙10⁶ times compared to quadratic finite elements. The developed post-processing procedure can be easily extended to unstructured meshes and can be independently used with existing post-processing techniques (e.g., with finite elements).

The purpose of this report is to document updates on the apparatus to simulate commercial vacuum drying procedures at the Nuclear Energy Work Complex at Sandia National Laboratories. Validation of the extent of water removal in a dry spent nuclear fuel storage system based on drying procedures used at nuclear power plants is needed to close existing technical gaps. Operational conditions leading to incomplete drying may have potential impacts on the fuel, cladding, and other components in the system during subsequent storage and disposal. A general lack of data suitable for model validation of commercial nuclear canister drying processes necessitates well-designed investigations of drying process efficacy and water retention. Scaled tests that incorporate relevant physics and well-controlled boundary conditions are essential to provide insight and guidance to the simulation of prototypic systems undergoing drying processes. This report documents fiscal year 2023 (FY23) updates on the Advanced Drying Cycle Simulator (ADCS). This apparatus was built to simulate commercial drying procedures and quantify the amount of residual water remaining in a pressurized water reactor (PWR) fuel assembly after drying. The ADCS was constructed with a prototypic 17×17 PWR fuel skeleton and waterproof heater rods to simulate decay heat. These waterproof heaters are the next generation design to heater rods developed and tested at Sandia National Laboratories in FY20. In FY23, a series of four simulated commercial drying tests was completed. This report presents the temperature and pressure histories of the drying tests as well as axial temperature profiles that can be compared to data from the Electric Power Research Institute (EPRI) High Burnup Demonstration TN-32B cask. Water content measurements and dew point calculations from a Hiden Analytical HPR-30 mass spectrometer are also presented in this report. Due to familiarization with this first-of-a-kind system, refinements to equipment calibration and test procedures have been identified to better match commercial drying cycles for future simulated tests. However, the presented data demonstrate the successful construction and operation of a viable research platform for quantifying residual water content closely approaching that expected in dry storage canisters during commercial drying procedures.

Work evaluating spent nuclear fuel (SNF) dry storage canister surface environments and canister corrosion progressed significantly in FY23, with the goal of developing a scientific understanding of the processes controlling initiation and growth of stress corrosion cracking (SCC) cracks in stainless steel canisters in relevant storage environments. The results of the work performed at Sandia National Laboratories (SNL) will guide future work and will contribute to the development of better tools for predicting potential canister penetration by SCC.

The US Department of Energy (DOE) Nuclear Energy Research Initiative funded the design and construction of the Seven Percent Critical Experiment (7uPCX) at Sandia National Laboratories. The start-up of the experiment facility and the execution of the experiments described here were funded by the DOE Nuclear Criticality Safety Program. The 7uPCX is designed to investigate critical systems with fuel for light water reactors in the enrichment range above 5 % ²³⁵U. The 7uPCX assembly is a water-moderated and -reflected array of aluminum-clad U(6.90 %)O₂ fuel rods. Other critical experiments performed in the 7uPCX assembly are documented in LEU-COMP-THERM-078, LEU-COMP-THERM-080, LEU-COMP-THERM-096, LEUCOMP-THERM-097, LEU-COMP-THERM-101, and LEU-COMP-THERM-102. The purpose of these experiments was to measure the effects of molybdenum in nearly-critical systems. The molybdenum was introduced into the fuel arrays as tubular sleeves that surrounded some of the fuel rods in the fuel arrays measured. Four hundred molybdenum tubes nominally 12.7 mm outside diameter, 498 mm long, with 0.762 mm wall thickness were provided for the experiments by the Institut de Radioprotection et de Sûreté Nucléaire (IRSN). Small polyethylene adapters at each end of the tubes were used to center each tube on a fuel rod in the assembly. The critical experiments were done using a set of triangular-pitched grid plates fabricated for these experiments. The grid plate set accommodated a fuel array of a total of 1261 fuel rod positions on a pitch of 0.610 in (1.5494 cm) in a series of 20 hexagonal rings surrounding the central fuel rod. The fuel used in these experiments was fabricated using unirradiated 6.90 % enriched UO₂ fuel pellets from fuel elements designed to be used in the internal nuclear superheater section of the Pathfinder boiling water reactor operated in South Dakota by the Northern States Power Company in the 1960s. The fuel elements were obtained from The Pennsylvania State University where they had been stored for many years. The fuel pellets in those fuel elements were removed from the original Incoloy cladding and reclad in 3003 aluminum tubes and end caps for use in the experiments reported here. The five critical experiments in this series were performed in August through December 2022, in the Sandia Critical Experiments (SCX) at the Sandia Pulsed Reactor Facility. Case 1 had no molybdenum sleeves, Case 2 had 208 molybdenum sleeves clustered at the center of the array, Case 3 had 397 molybdenum sleeves clustered at the center of the array, Case 4 had 175 molybdenum sleeves in the central position and in five alternating hexagonal rings, and Case 5 had 331 molybdenum sleeves in the central position and in seven alternating hexagonal rings. All five critical experiments are judged to be acceptable as benchmark experiments.

Accurate fuel oxidation mechanisms can enable predictive capabilities that aid in advancing combustion technologies. High-level computational kinetics can yield reasonable rate coefficients with uncertainties, in some cases, below a factor of 2. Computed rate coefficients can be constrained further by optimizing against experimental data. Here, we explore the application of genetic algorithm (GA) optimization to constrain computed rate coefficients in complex fuel oxidation mechanisms in conjunction with temperature-dependent species mole fractions from jet-stirred reactor (JSR) measurements. Cyclohexane is a model candidate for understanding the reactivity of cyclic fuels. In this work, we optimize the rate coefficients of the most recent literature cyclohexane mechanism, which incorporates theoretically computed rate coefficients for the reaction networks stemming from the first and second O2 addition pathways, against the experimental results of two separate literature JSR studies. Optimization consistency is evaluated by carrying out three GA optimizations: fitting to the temperature-dependent species mole fractions in each JSR experiment separately and simultaneously fitting the species mole fractions in both experiments. Local sensitivity analyses are used to identify five influential low-temperature oxidation reactions for optimization. Although the three optimizations do not yield identical rate coefficients, the direction of change in all five rate coefficients is consistent among the three optimizations. Performance of the models from the three optimizations is assessed against literature ignition delay times with differences in the level of agreement observed among the different optimizations. Comparisons are made with our recent optimization work of a cyclopentane oxidation master-equation model against time-resolved species concentrations, and insights and improvements of the strategy for constraining rate coefficients using GA optimization are discussed.

This memo summarizes the simulation of ductile failure propagation work conducted under the ASC project “V&V of Ductile Failure” conducted during FY 23. Physically, the failure propagation consists of crack propagation in the material. In the numerical setting—specifically in a finite element model—propagation can be accomplished through element death when critical conditions occur locally at an element that is then deleted from the simulation. The validation of the finite element models is evaluated by direct comparison between the experimental and simulation results regarding the rate of crack growth and its influence on the load-deflection response of the specimens tested. This work considers two geometries that display stable crack propagation under displacement-controlled conditions. The first geometry consists of hat specimens loaded in compression with nominally identical geometries but made with three different materials: Steel A286, Al 7075-T651 and 304L stainless steel. The three materials represent a range of ductility values that affect the response and crack propagation within the specimen. The crack induced propagates under an essentially mode-II type of deformation. The second geometry consists of a pre-cracked 304L stainless steel compact tension test specimen loaded so as to induce a mode-I deformation at the crack.

In this study, we demonstrate the ability of polarity inversion of sputtered aluminum scandium nitride thin films through post-fabrication processes with domain widths as small as 220 nm at a periodicity of 440 nm. An approach using photo- and electron-beam lithography to generate sub-quarter micrometer feature size with adjustable duty cycle through a lift-off process is presented. The film with a coercive field Ec+ of 5.35 MV/cm was exercised first with a 1 kHz triangular double bipolar wave and ultimately poled with a 0.5 kHz double monopolar wave using a Radiant Precision Premier II tester. The metal polar (M-polar) and nitrogen polar (N-polar) domains were identified and characterized through potassium hydroxide wet etching as well as piezoresponse force microscopy (PFM). Well-distinguished boundaries between the oppositely polarized domain regions were confirmed through the phase diagram of the PFM results. The relationship between the electrode width, poling voltage, and domain growth was experimentally studied and statistically analyzed, where 7.96 nm/V domain width broadening vs escalating poling voltage was observed. This method produces extremely high domain spatial resolution in III-nitride materials via poling and is transferable to a CMOS-compatible photolithography process. The spatial resolution of the periodically poled Al0.68Sc0.32N is suitable for second-harmonic generation of deep ultraviolet through quasi-phase-matching and RF MEMS operating in the X-Band spectrum.

Accurate prediction of ductile failure is critical to Sandia’s NW mission, but the models are computationally heavy. The costs of including high-fidelity physics and mechanics that are germane to the failure mechanisms are often too burdensome for analysts either because of the person-hours it requires to input them or because of the additional computational time, or both. In an effort to deliver analysts a tool for representing these phenomena with minimal impact to their existing workflow, our project sought to develop modern data-driven methods that would add microstructural information to business-as-usual calculations and expedite failure predictions. The goal is a tool that receives as input a structural model with stress and strain fields, as well as a machine-learned model, and output predictions of structural response in time, including failure. As such, our project spent substantial time performing high-fidelity, three-dimensional experiments to elucidate materials mechanisms of void nucleation and evolution. We developed crystal-plasticity finite-element models from the experimental observations to enrich the findings with fields not readily measured. We developed engineering length-scale simulations of replicated test specimens to understand how the engineering fields evolve in the presence of fine-scale defects. Finally, we developed deep learning convolutional neural networks, and graph-based neural networks to encode the findings of the experiments and simulations and make forward predictions in time for structural performance. This project demonstrated the power of data-driven methods for model development, which have the potential to vastly increase both the accuracy and speed of failure predictions. These benefits and the methods necessary to develop them are highlighted in this report. However, many challenges remain to implementing these in real applications, and these are discussed along with potential methods for overcoming them.

Tritium permeability in zirconium-based tritium getter critically impacts tritium storage and environmental safety during operation of tritium-producing burnable absorber rods (TPBARs). Previous experiments indicated that during irradiation operation, the hydrogen equilibrium pressured is increased. Further experimental and modeling studies suggested that the enhanced tritium release observed for reactor scale assemblies might be related to a thermal diffusion known as the Soret effect. A direct measurement of the Soret factor, however, has not been performed. To improve TPBAR and other nuclear applications, here we have applied two non-equilibrium molecular dynamics methods to study thermal diffusion of hydrogen isotopes in low-concentration zirconium hydrides. One of the methods produces sufficiently converged results to distinguish crystal orientation, isotope type, and concentration effects. In conclusion, with this method, crystal orientation, isotope type, and concentration effects are discussed.

Local crystallographic features negatively affect quantum spin defects by changing the local electrostatic environment, often resulting in degraded or varied qubit optical and coherence properties. Few tools exist that enable the deterministic synthesis and study of such intricate systems on the nano-scale, making defect-to-defect strain environment quantification difficult. In this paper, we highlight state-of-the-art capabilities from the U.S. Department of Energy’s Nanoscale Science Research Centers that directly address these shortcomings. Specifically, we demonstrate how complementary capabilities of nano-implantation and nano-diffraction can be used to demonstrate the quantum relevant, spatially deterministic creation of neutral divacancy centers in 4H silicon carbide, while investigating and characterizing these systems on the ≤ 25 nm scale with strain sensitivities on the order of 1 × 10 − 6 , relevant to defect formation dynamics. This work lays the foundation for ongoing studies into the dynamics and deterministic formation of low strain homogeneous quantum relevant spin defects in the solid state.

This report describes specific activities in the Fiscal Year (FY) 2023 associated with the Geologic Disposal Safety Assessment (GDSA) Repository Systems Analysis (RSA) work package funded by the Spent Fuel and Waste Science and Technology (SFWST) Campaign of the U.S. Department of Energy Office of Nuclear Energy (DOE-NE), Office of Spent Fuel and Waste Disposition (SFWD). The overall objective of the GDSA RSA work package is to develop generic deep geologic repository concepts and system performance assessment (PA) models in several potential host-rock environments, and to simulate and analyze these generic repository concepts and models using the GDSA Framework toolkit and other tools as needed. The specific GDSA goal addressed in this report is reference case development, simulation, and analysis for the unsaturated alluvium (UZ), one of the four potential host-rocks considered by the GDSA. Further, we aim to exercise the simulation tools and methodologies under development by GDSA for PA modelling.

Aging of silicon photovoltaic (PV) module packaging is one of the greatest limiters of PV module service lifetimes. Module characterization typically focuses on power degradation metrics, which do not convey the complexities of often simultaneous degradation mechanisms. In this work, PV modules with pristine references and known fielding histories were investigated by non-destructive and destructive methods. Modules from Canadian Solar, Mission Solar, and Hanwha Q-Cells were fielded for up to three years; select modules were removed from fielding each year for coring to allow for characterization of the encapsulant. Modules are commonly encapsulated with two protective layers of partially-crystalline ethylene vinyl acetate (EVA) polymer that must undergo a crosslinking reaction to achieve desired properties. The extent of crystallinity of the encapsulants as studied by differential scanning calorimetry showed differences between manufacturers and over time. Some encapsulants showed different magnitudes of crystal sizes which changed after fielding; encapsulants with the monodisperse crystal sizes did not change with fielding. This is due to differences in thermal history. These results have implications for stress development during module aging, since EVA crystal melting and crosslinking reactions can result in encapsulant density changes.

Microneedle sensors could enable minimally-invasive, continuous molecular monitoring – informing on disease status and treatment in real-time. Wearable sensors for pharmaceuticals, for example, would create opportunities for treatments personalized to individual pharmacokinetics. Here, we demonstrate a commercial-off-the-shelf (COTS) approach for microneedle sensing using an electrochemical aptamer-based sensor that detects the high-toxicity antibiotic, vancomycin. Wearable monitoring of vancomycin could improve patient care by allowing targeted drug dosing within its narrow clinical window of safety and efficacy. To produce sensors, we miniaturize the electrochemical aptamer-based sensors to a microelectrode format, and embed them within stainless steel microneedles (sourced from commercial insulin pen needles). The microneedle sensors achieve quantitative measurements in body-temperature undiluted blood. Further, the sensors effectively maintain electrochemical signal within porcine skin. This COTS approach requires no cleanroom fabrication or specialized equipment, and produces individually-addressable, sterilizable microneedle sensors capable of easily penetrating the skin. In the future, this approach could be adapted for multiplexed detection, enabling real-time monitoring of a range of biomarkers.

The Sibling Pin test campaign is a Department of Energy (DOE) research activity within the Spent Fuel and Waste Science and Technology (SFWST) program that is tasked with characterization of high burnup (HBU) fuel in support of the High Burnup Spent Fuel Data Project. Of the 25 fuel rods in the Sibling Pin inventory, approximately 9 rod lengths have been consumed during the first phase (Phase I) of the test campaign leaving approximately 16 rod lengths for the second phase (Phase II) of testing. This plan outlines the Phase II testing and the motivations for performing these tests. Priorities for Phase II testing are based on previously identified knowledge gaps, lessons-learned from Phase I work, the original objectives of the High Burnup Spent Fuel Data Project and the Sibling Pin test campaign, and input from external stakeholders. The priorities for Phase II testing are to obtain data to characterize the effects of annealing on cladding mechanical properties and fuel rod performance, to quantify the creep behavior of cladding materials and fuel rods and the effects of creep deformations on the performance of cladding and fuel rods, and to gather data to support the final closure of the hydride reorientation and radial hydride induced embrittlement gap for HBU fuel rods.

Translating the surging interest in neuromorphic electronic components, such as those based on nonlinearities near Mott transitions, into large-scale commercial deployment faces steep challenges in the current lack of means to identify and design key material parameters. These issues are exemplified by the difficulties in connecting measurable material properties to device behavior via circuit element models. Here, the principle of local activity is used to build a model of VO2/SiN Mott threshold switches by sequentially accounting for constraints from a minimal set of quasistatic and dynamic electrical and high-spatial-resolution thermal data obtained via in situ thermoreflectance mapping. By combining independent data sets for devices with varying dimensions, the model is distilled to measurable material properties, and device scaling laws are established. The model can accurately predict electrical and thermal conductivities and capacitances and locally active dynamics (especially persistent spiking self-oscillations). The systematic procedure by which this model is developed has been a missing link in predictively connecting neuromorphic device behavior with their underlying material properties, and should enable rapid screening of material candidates before employing expensive manufacturing processes and testing procedures.

Motivated by recent experimental results we calculate from first-principles the lifetime of low-energy quasiparticles in bilayer graphene (BLG). Here, we take into account the scattering rate arising from electron-electron interactions within the GW approximation for the electron self-energy and consider several p-type doping levels ranging from 0 to ρ ≈ 2.4 × 10¹² holes/cm². In the undoped case we find that the average inverse lifetime scales linearly with energy away from the charge neutrality point, with values in good agreement with experiments. The decay rate is approximately three times larger than in monolayer graphene, a consequence of the enhanced screening in BLG. In the doped case, the dependence of the inverse lifetime on quasiparticle energy acquires a non-linear component due to the opening of an additional decay channel mediated by acoustic plasmons.

The Spent Fuel Waste Disposition (SFWD) program is planning to conduct a full-scale seismic shake table test on the dry storage systems of spent nuclear fuel (SNF) to close the gap related to seismic loads on fuel assemblies in dry storage systems. This test will allow for quantifying the strains and accelerations on surrogate fuel assembly hardware and cladding during earthquakes of different magnitudes and frequency content. Full-scale testing is needed because a dry storage system is a complex and highly nonlinear system making it hard to predict (model) the responses to seismic excitations. The non-linearity arises from the multiple spatial gaps in the system – between fuel rods and the basket, between the basket and dry storage canister, between the dry storage canister and the storage cask (overpack), and ventilation gaps. The non-linearities pose significant limitations on the value of tests with scaled systems.

Here, we study the impact of deep-level defects on trap-assisted Auger–Meitner recombination in c-plane InGaN/GaN LEDs using a small-signal electroluminescence (SSEL) method and deep-level optical spectroscopy (DLOS). Carrier dynamics information, including carrier lifetime, recombination rate, and carrier density, is obtained from SSEL, while DLOS is used to obtain the deep-level defect density. Through fitting the nonradiative recombination rates of wafers with different deep-level defect densities, we obtain the Shockley–Read–Hall (SRH) and trap-assisted Auger–Meitner recombination (TAAR) coefficients. We show that defect-related nonradiative recombination, including both SRH and TAAR, accounts for a relatively small fraction of the total nonradiative recombination, which is dominated by intrinsic Auger–Meitner recombination. The interplay between carrier localization and Coulomb enhancement has a different impact on radiative and intrinsic Auger–Meitner recombination. Evidence is presented that the imbalance between the change of radiative and intrinsic Auger–Meitner recombination is the primary cause of the efficiency droop at high carrier densities in the samples studied.

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.

We study the impact of deep-level defects on trap-assisted Auger-Meitner recombination in c-plane InGaN/GaN LEDs using a small-signal electroluminescence (SSEL) method and deep-level optical spectroscopy (DLOS). Carrier dynamics information, including carrier lifetime, recombination rate, and carrier density, is obtained from SSEL, while DLOS is used to obtain the deep-level defect density. Through fitting the nonradiative recombination rates of wafers with different deep-level defect densities, we obtain the Shockley-Read-Hall (SRH) and trap-assisted Auger-Meitner recombination (TAAR) coefficients. We show that defect-related nonradiative recombination, including both SRH and TAAR, accounts for a relatively small fraction of the total nonradiative recombination, which is dominated by intrinsic Auger-Meitner recombination. The interplay between carrier localization and Coulomb enhancement has a different impact on radiative and intrinsic Auger-Meitner recombination. Evidence is presented that the imbalance between the change of radiative and intrinsic Auger-Meitner recombination is the primary cause of the efficiency droop at high carrier densities in the samples studied.

Cast Monel alloys are used in applications requiring a combination of good mechanical properties and excellent resistance to corrosion. Despite prevalent industrial use, relatively few studies have been conducted to investigate the relationships between composition, solidification behavior, and microstructure. Given that these alloys are used in the cast and welded conditions, these factors have a significant influence over the material properties. Here, in this work, microstructural characterization, electron probe microanalysis, X-ray diffraction, and differential scanning calorimetry were used to study how changes in Si and Nb concentrations affected the solidification path and microstructure of Monel alloys. It was found that increasing Nb concentration stabilized higher amounts of MC carbides and suppressed graphite formation during solidification. It was also found that the high nominal concentration and segregation of Si to the liquid led to the formation of Ni₃₁Si₁₂ and other silicides via terminal eutectic reactions at the end of solidification. A pseudo-binary solidification diagram was constructed using experimental data and was applied to predict the mass fraction of solidified eutectic as a function of composition. The modeled microstructures were found to be in good agreement with experimentally measured phase fractions.

Heterogenous materials under shock compression can be expected to reach different shock states throughout the material according to local differences in microstructure and the history of wave propagation. Here, a compact, multiple-beam focusing optic assembly is used with high-speed velocimetry to interrogate the shock response of porous tantalum films prepared through thermal-spray deposition. The distribution of particle velocities across a shocked interface is compared to results obtained using a set of defocused interferometric beams that sampled the shock response over larger areas. The two methods produced velocity distributions along the shock plateau with the same mean, while a larger variance was measured with narrower beams. The finding was replicated using three-dimensional, mesoscopically resolved hydrodynamics simulations of solid tantalum with a pore structure mimicking statistical attributes of the material and accounting for radial divergence of the beams, with agreement across several impact velocities. Accounting for pore morphology in the simulations was found to be necessary for replicating the rise time of the shock plateau. The validated simulations were then used to show that while the average velocity along the shock plateau could be determined accurately with only a few interferometric beams, accurately determining the width of the velocity distribution, which here was approximately Gaussian, required a beam dimension much smaller than the spatial correlation lengthscale of the velocity field, here by a factor of ∼30×, with implications for the study of other porous materials.

Here, we perform all-atom molecular dynamics simulations of lithium triflate in 1,2-dimethoxyethane using six different literature force fields. This system is representative of many experimental studies of lithium salts in solvents and polymers. We show that multiple historically common force fields for lithium ions give qualitatively incorrect results when compared with those from experiments and quantum chemistry calculations. We illustrate the importance of correctly selecting force field parameters and give recommendations on the force field choice for lithium electrolyte applications.

Objective. Our objective is to formulate the problem of the magnetoencephalographic (MEG) sensor array design as a well-posed engineering problem of accurately measuring the neuronal magnetic fields. This is in contrast to the traditional approach that formulates the sensor array design problem in terms of neurobiological interpretability the sensor array measurements. Approach. We use the vector spherical harmonics (VSH) formalism to define a figure-of-merit for an MEG sensor array. We start with an observation that, under certain reasonable assumptions, any array of m perfectly noiseless sensors will attain exactly the same performance, regardless of the sensors' locations and orientations (with the exception of a negligible set of singularly bad sensor configurations). We proceed to the conclusion that under the aforementioned assumptions, the only difference between different array configurations is the effect of (sensor) noise on their performance. We then propose a figure-of-merit that quantifies, with a single number, how much the sensor array in question amplifies the sensor noise. Main results. We derive a formula for intuitively meaningful, yet mathematically rigorous figure-of-merit that summarizes how desirable a particular sensor array design is. We demonstrate that this figure-of-merit is well-behaved enough to be used as a cost function for a general-purpose nonlinear optimization methods such as simulated annealing. We also show that sensor array configurations obtained by such optimizations exhibit properties that are typically expected of 'high-quality' MEG sensor arrays, e.g. high channel information capacity. Significance. Our work paves the way toward designing better MEG sensor arrays by isolating the engineering problem of measuring the neuromagnetic fields out of the bigger problem of studying brain function through neuromagnetic measurements.

Alkene ozonolysis generates short-lived Criegee intermediates that are a significant source of hydroxyl (OH) radicals. This study demonstrates that roaming of the separating OH radicals can yield alternate hydroxycarbonyl products, thereby reducing the OH yield. Specifically, hydroxybutanone has been detected as a stable product arising from roaming in the unimolecular decay of the methyl-ethyl-substituted Criegee intermediate (MECI) under thermal flow cell conditions. The dynamical features of this novel multistage dissociation plus a roaming unimolecular decay process have also been examined with ab initio kinetics calculations. Experimentally, hydroxybutanone isomers are distinguished from the isomeric MECI by their higher ionization threshold and distinctive photoionization spectra. Moreover, the exponential rise of the hydroxybutanone kinetic time profile matches that for the unimolecular decay of MECI. A weaker methyl vinyl ketone (MVK) photoionization signal is also attributed to OH roaming. Complementary multireference electronic structure calculations have been utilized to map the unimolecular decay pathways for MECI, starting with 1,4 H atom transfer from a methyl or methylene group to the terminal oxygen, followed by roaming of the separating OH and butanonyl radicals in the long-range region of the potential. Roaming via reorientation and the addition of OH to the vinyl group of butanonyl is shown to yield hydroxybutanone, and subsequent C-O elongation and H-transfer can lead to MVK. A comprehensive theoretical kinetic analysis has been conducted to evaluate rate constants and branching yields (ca. 10-11%) for thermal unimolecular decay of MECI to conventional and roaming products under laboratory and atmospheric conditions, consistent with the estimated experimental yield (ca. 7%).

A wind rotor is disclosed that produces energy optimally for a given thrust overturning moment. By designing rotors with suboptimal aerodynamic efficiency, they can have optimal thrust performance, which will reduce the substructure cost and/or enable greater energy capture for a given substructure.

As a part of NASA's efforts in space, options are being examined for an Artemis moon base project to be deployed. This project requires a system of interconnected, but separate, DC microgrids for habitation, mining, and fuel processing. This in-place use of power resources is called in-situ resource utilization (ISRU). These microgrids are to be separated by 9-12 km and each contains a photovoltaic (PV) source, energy storage systems (ESS), and a variety of loads, separated by level of criticality in operation. The separate microgrids need to be able to transfer power between themselves in cases where there are generation shortfall, faults, or other failures in order to keep more critical loads running and ensure safety of personnel and the success of mission goals. In this work, a 2 grid microgrid system is analyzed involving a habitation unit and a mining unit separated by a tie line. A set of optimal controls that has been developed, including power flow controls on the tie line, dispatch of PV generation, and dispatch of non-critical loads, is analyzed, and validated in hardware on the Secure Scalable Microgrid Testbed (SSMTB). This testbed includes hardware emulators for a variety of energy sources, energy storage devices, pulsed loads, and other loads.

This report describes specific activities in the Fiscal Year (FY) 2023 associated with the Geologic Disposal Safety Assessment (GDSA) Repository Systems Analysis (RSA) work package funded by the Spent Fuel and Waste Science and Technology (SFWST) Campaign of the U.S. Department of Energy Office of Nuclear Energy (DOE-NE), Office of Spent Fuel and Waste Disposition (SFWD).

Abstract not provided.

A critical mission need exists to develop new materials that can withstand extreme environments and multiple sequential threats. High entropy materials, those containing 5 or more metals, exhibit many exciting properties which would potentially be useful in such situations. However, a particularly hard challenge in developing new high entropy materials is determining a priori which compositions will form the desired single phase material. The project outlined here combined several modeling and experimental techniques to explore several structure-property-relationships of high entropy ceramics in an effort to better understand the connection between their compositional components, their observed properties, and stability. We have developed novel machine learning algorithms which rapidly predict stable high entropy ceramic compositions, identified the stability interplay between configurational entropy and cation defects, and tested the mechanical stability of high entropy oxides using the unique capabilities at the Dynamic Compression Sector facility and the Saturn accelerator.

Developed in 2018, PathTrace is a software package built with the intention of making path analysis simple and intuitive. PathTrace is a top-down pathway analysis software where a user is able to explore vulnerable pathways into a facility. The intention of utilizing a software tool like PathTrace is to characterize an existing physical protection system (PPS) and to upgrade the system to achieve a high level of response interruption, or probability of interruption (P_I) of the adversary. There are four steps for conducting path analysis using PathTrace. The first step is to identify an image to use to build the model and scale the model within PathTrace using a section of known distance (wall or fence perimeter, for example). The scaling process will produce a grid of cells through which the user is able to build a model. The second step is to fill out the grid of cells with four categories of materials: Barriers, Detection Areas, Jumps, and Targets. These materials apply associated delay and detection values to the cells in which they are applied. The third step is to represent the adversary and response forces. The adversaries are represented by their capabilities in interacting with the materials identified in step two, and the response is represented by how quickly they will be able to respond to an adversary attack. Finally, the user is able to take all of the information from the previous three steps and perform a Most Vulnerable Path (MVP) analysis. In this stage, the user is able to visualize vulnerable adversary pathways and reason about how to upgrade these pathways to provide a high level of P_I.

Abstract not provided.

Accurate simulation of the evolution of polar ice-sheets requires a massive amount of computational power. In order to take advantage of the newest generation of supercomputing clusters, the Albany Land Ice code has been modernized for performance portability across a variety of parallel architectures, with a focus on enabling end-to-end GPU capability. Albany uses a multigrid preconditioning approach for solving linear systems via performance portable smoothers from the Trilinos package Ifpack2. Since the Albany Land Ice code is constantly evolving and both Albany and Trilinos are in constant development, it is likely that the optimal choice of solver parameters will change over time. It is therefore critical to have an automatic performance tuning framework to ensure that the best possible performance is maintained. Toward this effect, we have developed an automatic performance tuning framework to determine the best fine- and coarse-grid smoothing algorithms and parameters. We treat the underlying performance model of the linear solve as a black box and use the python-based GPTune Bayesian optimization library to determine the optimal smoother choice and parameters. Using this approach, we have found smoothers and their corresponding parameters that result in, on average, 1.2 times faster, and up to 1.5 times faster solve-times than our manually-tuned parameters. We also show that the proposed auto-tuning approach produces reliably better parameters than naive black box optimization techniques like random search for a given function evaluation budget. By implementing our tuning framework in the Python-based workflow management tool parsl, we also ensure that we efficiently use available computing resources during the tuning process and avoid unnecessary long wait times in computing cluster job queues.

Traditionally, electric grid planning seeks to maintain safe, reliable, efficient, and affordable service for current and future customers. As policies, expectations of the energy system, and the threat landscape evolve, additional objectives for power system planners are emerging, including decarbonization, resilience, and equity. Renewable and clean energy goals, especially in the context of deep decarbonization strategies, are changing the mix of resources on the electric grid and prompting new considerations for grid architecture. The increased frequency and severity of extreme weather events over the last two decades, coupled with cybersecurity concerns, have elevated resilience as a key system need. More recently, there has been greater focus on equity and energy justice in grid planning to ensure that disadvantaged communities are not adversely affected by grid modernization and have equal access to its benefits. In response, new thinking around multi-objective decision planning is exploring improvements in grid planning processes to better integrate approaches to meet decarbonization, resilience, and equity objectives. To provide a foundation for this work, a series of white papers was produced to summarize these emerging objectives.

Superconducting qubits have reached the point where system designers are worried about the heat that control wiring brings into the cryostat. To continue scaling cryogenic quantum systems, control solutions that work inside the cold space must be explored. One possibility is to use control electronics that is native to superconductivity, so called single-flux-quantum (SFQ) circuitry, to form an interface between qubits and whatever other electronics is needed to control eventual quantum systems. To begin exploring the utility of SFQ as control circuitry, we performed modeling and experiments on qubit readout using ballistic fluxons which are SFQ in the limit of ballistic fluxon transport. Our modeling results show that a flavor of qubit, the fluxonium, can be read out using ballistic fluxons. We designed test samples to prove some of the key concepts needed for such a readout but were ultimately unable to getting a working demonstration. The lack of testing success was due to challenges in fabrication and running short of time to perform testing rather than a fundamental problem with our analysis.

Abstract not provided.

We demonstrate the discrimination of ground-state hyperfine manifolds of a cesium atom in an optical tweezer using a simple probe beam with Formula Presented% detection fidelity and 0.9(2)% detection-driven loss of bright-state atoms. Our detection infidelity of Formula Presented% is an order of magnitude better than previously published low-loss readout results for alkali-metal atoms in optical tweezers. We achieve these results by identifying and mitigating an extra depumping mechanism due to stimulated Raman transitions induced by trap light in the presence of probe light. In this work, complex optical systems and stringent vacuum pressures are not required, enabling straightforward adoption of our techniques on contemporary experiments.

Traditionally, electric grid planning seeks to maintain safe, reliable, efficient, and affordable service for current and future customers. As policies, expectations of the energy system, and the threat landscape evolve, additional objectives for power system planners are emerging, including decarbonization, resilience, and equity. Renewable and clean energy goals, especially in the context of deep decarbonization strategies, are changing the mix of resources on the electric grid and prompting new considerations for grid architecture. The increased frequency and severity of extreme weather events over the last two decades, coupled with cybersecurity concerns, have elevated resilience as a key system need. More recently, there has been greater focus on equity and energy justice in grid planning to ensure that disadvantaged communities are not adversely affected by grid modernization and have equal access to its benefits. In response, new thinking around multi-objective decision planning is exploring improvements in grid planning processes to better integrate approaches to meet decarbonization, resilience, and equity objectives. To provide a foundation for this work, a series of white papers was produced to summarize these emerging objectives.

The broad goal of this project was to develop new analytical and numerical insights for how important photonic processes can be improved using recently discovered principles in topological and non-Hermitian physics. In particular, there are two recent discoveries that we aimed to harness to achieve this goal. First, it was discovered in condensed matter physics that crystalline symmetries can protect low-dimensional topologically protected states in lattices without the need for breaking time-reversal symmetry. These so-called ‘higher-order’ topological systems represent an important development for photonic systems, where it is very difficult to break time-reversal symmetry, and which had been previously thought necessary to realize topological phenomena. Second, the last decade has seen a significant amount of interest in phenomena which are unique to non-Hermitian systems, i.e., systems which do not conserve energy. For example, spatially patterned gain and loss can be used to realize exceptional points, which are degeneracies in a system’s spectrum where the system becomes defective, while the existence of radiative losses also enables a new route to confinement through bound states in the continuum. For such non-Hermitian phenomena, photonics again represents a critical platform, as photonic systems naturally lose light to their radiative environments, making them generally non-Hermitian, and it is also possible to incorporate additional gain or loss. Based on these broad principles, we pursued a range of projects to harness bound states in the continuum in a variety of different systems and architectures, develop real-space methods for classifying topological systems to yield better photonic design principles, and a novel Brillouin-based fiber laser for sensing strain.

Salt offers an optimal medium for the permanent isolation of heat-producing radioactive waste due to its impermeability, high thermal conductivity, and ability to close fractures through creep. A thorough understanding of the thermal-hydrological-mechanical (THM) processes, encompassing brine migration, is fundamental for secure radioactive waste disposal within salt formations. At the Waste Isolation Pilot Plant (WIPP), we conducted joint in situ geophysical monitoring experiments during active heating to investigate brine migration near excavations. This experiment incorporated electrical resistivity tomography (ERT) alongside high-resolution fiber-optic-based distributed temperature sensing within a controlled heating experiment. Additionally, discrete element model (DEM) based numerical simulations were conducted to simulate THM processes during heating, providing a more mechanistic understanding of the coupled processes leading to the observed changes in geophysical measurements. During heating, resistivity shifts near the heater were reasonably explained by temperature effects. However, in more distant, cooler regions, the resistivity decrease exceeded predictions based solely on temperature. DEM simulations highlighted brine migration, propelled by pore pressure gradients, as the likely primary factor contributing to the additional resistivity decline beyond temperature effects. The comparison between the predicted ERT responses and observations was much improved when considering the effects of brine migration based on the DEM simulations. These geophysical and simulation findings shed light on brine migration in response to salt heating, enhancing our understanding of the coupled THM processes in salt for safe radioactive waste disposal.

In this project, we experimented the focused ion beam (FIB) based fabrications of semiconductor quantum dots (QDs) by using metal nano particles (NPs) (e.g., Al) on semiconductor as a template and by means of the FIB induced direct metal-to-QD conversion. We have examined effect of the experimental conditions, including Ga⁺ ion energy and dose as well as substrate temperature. The results of experiments have shown AlGaSb QD formation on GaSb substrate can be achieved under certain conditions but there are many challenges about the techniques, including compositional nonuniformity of the QDs formed, partial conversion of the metal NP to QD, and high defect concentration in the QDs.

Abstract not provided.

Distribution system model calibration is a key enabling task for incipient failure identification within the distribution system. This report summarizes the work and publications by Sandia National Laboratories on the GMLC project titled “Incipient Failure Identification for Common Grid Asset Classes”. This project was a joint effort between Sandia National Laboratories, Lawrence Livermore National Laboratory, National Energy Technology Laboratory, and Oak Ridge National Laboratory. The included work covers distribution system topology identification, transformer groupings, phase identification, regulator and tap position estimation, and the open-source release and implementation of the developed algorithms.

We propose a novel statistical inference methodology for multiway count data that is corrupted by false zeros that are indistinguishable from true zero counts. Our approach consists of zero-truncating the Poisson distribution to neglect all zero values. This simple truncated approach dispenses with the need to distinguish between true and false zero counts and reduces the amount of data to be processed. Inference is accomplished via tensor completion that imposes low-rank tensor structure on the Poisson parameter space. Our main result shows that an N-way rank-R parametric tensor M ∈ (0, ∞)I×.....×I generating Poisson observations can be accurately estimated by zero-truncated Poisson regression from approximately IR2 log22(I) non-zero counts under the nonnegative canonical polyadic decomposition. Our result also quantifies the error made by zero-truncating the Poisson distribution when the parameter is uniformly bounded from below. Therefore, under a low-rank multiparameter model, we propose an implementable approach guaranteed to achieve accurate regression in under-determined scenarios with substantial corruption by false zeros. Several numerical experiments are presented to explore the theoretical results.

Abstract not provided.

Abstract not provided.

A direct numerical simulation (DNS) campaign is deployed for a series of confined downward oriented, non-isothermal turbulent impinging jet configurations. A baseline Reynolds number of 9960 is obtained through a precursor DNS pipe flow simulation (Reτ=505). Three jet temperature configurations (confinement height to nozzle diameter of three) enter a cylindrical domain that share ambient and impingement plate temperatures (298.15K). The range of jet temperatures are crafted such that the ratio of inlet to ambient density varies from unity to 0.52, showcasing the effect of density disparity on flow characteristics such as core collapse, radial mixing of momentum and energy, near-wall stagnation behavior, wall-jet profiles, and large-scale vortical structures. Surface quantities provided include mean radial heat flux and wall-shear stress profiles, and heat flux histograms at select radial stations. Results showcase increased radial normal stresses for higher temperature jets that support increased mixing, resulting in large-scale recirculation structures that are smaller, while retaining similar normalized radial wall profiles for shear stress, heat flux and pressure. Radial plots for wall shear stress and Nusselt number showcase strong radial decay as compared to previous configurations that share similar jet and ambient temperatures. For the 373.15 K case, a Gaussian-like histogram for heat fluxes at the impingement plate transitions to a log-normal profile as radial distances increase. In contrast, the 573.15 K configuration displays a bi-modal heat flux characteristic at the impingement plate, and in similar manner to the moderate temperature counterpart, transitions to a log-normal profile at larger radial distances.

Traditionally, electric grid planning seeks to maintain safe, reliable, efficient, and affordable service for current and future customers. As policies, expectations of the energy system, and the threat landscape evolve, additional objectives for power system planners are emerging, including decarbonization, resilience, and equity. Renewable and clean energy goals, especially in the context of deep decarbonization strategies, are changing the mix of resources on the electric grid and prompting new considerations for grid architecture. The increased frequency and severity of extreme weather events over the last two decades, coupled with cybersecurity concerns, have elevated resilience as a key system need. More recently, there has been greater focus on equity and energy justice in grid planning to ensure that disadvantaged communities are not adversely affected by grid modernization and have equal access to its benefits. In response, new thinking around multi-objective decision planning is exploring improvements in grid planning processes to better integrate approaches to meet decarbonization, resilience, and equity objectives. To provide a foundation for this work, a series of white papers was produced to summarize these emerging objectives.

In this report we present our findings and outcomes of the NNRDS (analysis of Neural Networks as Random Dynamical Systems) project. The work is largely motivated by the analogy of a large class of neural networks (NNs) with a discretized ordinary differential equation (ODE) schemes. Namely, residual NNs, or ResNets, can be viewed as a discretization of neural ODEs (NODEs) where the NN depth plays the role of the time evolution. We employ several legacy tools from ODE theory, such as stiffness, nonlocality, autonomicity, to enable regularization of ResNets thus improving their generalization capabilities. Furthermore, armed with NN analysis tools borrowed from the ODE theory, we are able to efficiently augment NN predictions with uncertainty overcoming wellknown dimensionality challenges and adding a degree of trust towards NN predictions. Finally, we have developed a Python library QUiNN (Quantification of Uncertainties in Neural Networks) that incorporates improved-architecture ResNets, besides classical feed-forward NNs, and contains wrappers to PyTorch NN models enabling several major classes of uncertainty quantification methods for NNs. Besides synthetic problems, we demonstrate the methods on datasets from climate modeling and materials science.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Calibrated measurements of lightning optical emissions are critical for both quantifying the impacts of lightning in our atmosphere and devising detection instruments with sufficient dynamic range capable of yielding close to 100% detection efficiency. However, to date, there is only a limited number of investigations that have attempted to take such calibrated measurements. In this work, we report the power radiated by lightning in both visible and infrared bands, assuming isotropic emission, and accounting for atmospheric absorption. More precisely, we report peak radiated power and total radiated energy in the combined visible plus near-infrared range (VNIR, 0.34–1.1 μm), around the Hα line (652–667 nm), and for the 2–2.5 μm infrared band. The estimated peak power and total energy radiated by negative cloud-to-ground return strokes in the VNIR range is 130 MW and 20 kJ, respectively. Additionally, we detected peak radiated powers of 12 and 0.19 MW in the Hα and infrared bands, respectively. We cross-reference the optical data set with peak current reported by a lightning detection network. The resulting trend is that optical power emitted around the Hα line scales with peak return stroke current according to a power law with exponent equal to 1.25. This trend, which should be approximately true across the entire visible spectrum, can be attributed to the plasma negative differential resistance of the lightning return stroke channel. We conclude by discussing the challenges in performing calibrated measurements of lightning optical power in different bands and comparing the results with previously-collected data with different experimental setups, observation conditions, and calibration methods.

Abstract not provided.

Abstract not provided.

Abstract not provided.

This project applies methods in Bayesian inference and modern statistical methods to quantify the value of new experimental data, in the form of new or modified diagnostic configurations and/or experiment designs. We demonstrate experiment design methods that can be used to identify the highest priority diagnostic improvements or experimental data to obtain in order to reduce uncertainties on critical inferred experimental quantities and select the best course of action to distinguish between competing physical models. Bayesian statistics and information theory provide the foundation for developing the necessary metrics, using two high impact experimental platforms on Z as exemplars to develop and illustrate the technique. We emphasize that the general methodology is extensible to new diagnostics (provided synthetic models are available), as well as additional platforms. We also discuss initial scoping of additional applications that began development in the last year of this LDRD.

Abstract not provided.

Human activities involving subsurface reservoirs—resource extraction, carbon and nuclear waste storage—alter thermal, mechanical, and chemical steady-state conditions in these systems. Because these systems exist at lithostatic pressures, even minor chemical changes can cause chemically assisted deformation. Therefore, understanding how chemical effects control geomechanical properties is critical to optimizing engineering activities. The grand challenge in predicting the effect of chemical processes on mechanical properties lays in the fact that these phenomena take place at molecular scales, while they manifest all the way to reservoir scales. To address this fundamental challenge, we investigated chemical effects on deformation in model and real systems spanning molecular- to centimeter scales. We used theory, experiment, molecular dynamics simulation, and statistical analysis to (1) identify the effect of simple reactions, such as hydrolysis, on molecular structures in interfacial regions of stressed geomaterials; (2) quantify chemical effects on the bulk mechanical properties, fracture and displacement for granular rocks and single crystals; (3) develop initial understanding of universal scaling for individual displacement events in layered geomaterials; and (4) develop analytic approximations for the single-chain mechanical response utilizing asymptotically correct statistical thermodynamic theory. Taken together, these findings advance the challenging field of chemo-mechanics.

The injection and storage of anthropogenic CO2 in the subsurface is being deployed as a climate change mitigation tool; however, diagenetic-paragenetic heterogeneity in sandstone reservoirs often contributes to interval specific chemomechanical changes that affect injection and can increase leakage risk. Here, we address reservoir heterogeneities’ impact on chemomechanical changes in a macroporous-dominated lithofacies of Morrow B sandstone, a formation containing several diagenetically-distinct hydraulic facies while undergoing enhanced oil recovery (EOR) and carbon dioxide (CO2) sequestration. We performed three flow-through experiments using a CO2-charged or uncharged formation water combined with four indirect tensile strength tests per post-test sample. We then used the microstructure and paragenetic sequence to understand chemomechanical weakening with key observations as follows: dissolution of carbonates and feldspars changed porosity; increased permeability led to reclassifying each sample in a different hydraulic flow unit; decreased ultrasonic velocity; and did not lead to a loss of tensile strength. Tensile strength maintenance occurred due to the low abundance and minor dissolution of siderite, the stability of quartz, and the relative position of diagenetic ankerite within feldspar. This macroporous-dominated lithofacies is the primary reservoir for the Morrow B Sandstone, and is analogous to other porous sandstone reservoirs. It represents an end-member of a chemomechanically low-risk siliceous CO2 sequestration and CO2-EOR reservoir.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The ion velocity distribution functions of thermonuclear plasmas generated by spherical laser direct drive implosions are studied using deuterium-tritium (DT) and deuterium-deuterium (DD) fusion neutron energy spectrum measurements. A hydrodynamic Maxwellian plasma model accurately describes measurements made from lower temperature (<10 keV), hydrodynamiclike plasmas, but is insufficient to describe measurements made from higher temperature more kineticlike plasmas. The high temperature measurements are more consistent with Vlasov-Fokker-Planck (VFP) simulation results which predict the presence of a bimodal plasma ion velocity distribution near peak neutron production. These measurements provide direct experimental evidence of non-Maxwellian ion velocity distributions in spherical shock driven implosions and provide useful data for benchmarking kinetic VFP simulations.

This report describes a Late-Start LDRD project on a compact neutron time-of-flight (nTOF) diagnostic, based on the Ultra-Fast Pixel Array Camera (UPAC) platform. By leveraging advancements in microelectronics, the UPAC nTOF is designed to address existing challenges of current nTOFs at the Z-machine that rely on large scintillators and photomultiplier tubes using bulky, discrete digitizing units hindering scalability with signal routing presenting an additional challenge. The report describes the UPAC diagnostic platform, outlines the calibration of the diagnostic and the latest results of fielding the UPAC as a ride-a-long diagnostic on the Magnetized liner inertial fusion (MagLIF) campaigns at the Z-machine supported by the Late-Start LDRD, and finally discusses conclusions and proposed future work.

Cybersecurity is a persistent concern to the safety and security of Nuclear Power Plants (NPPs), but has lacked data-driven, evidence-based research. Rigorous cybersecurity analysis is critical for the licensing of advanced reactors using a performance-based approach. One tool that enables cybersecurity analysis is modeling and simulation. The nuclear industry makes extensive use of modeling and simulation throughout the decision process but lacks a method to incorporate cybersecurity analysis with existing models. To meet this need, the Advanced Reactor Cyber Analysis and Development Environment (ARCADE) was developed. ARCADE is a suite of publicly available tools that can be used to develop emulations of industrial control system devices and networks and integrate those emulations with physics simulators. This integration of cyber emulations and physics models enables rigorous cyber-physical analysis of cyber-attacks on NPP systems. This report provides an overview of key considerations for using ARCADE with existing physics models and demonstrates ARCADE’s capabilities for cybersecurity analysis. Using a model of the Small Modular Advanced High Temperature Reactor (SmAHTR), ARCADE was able to determine the sensitivity of the primary heat exchangers (PHX) to coordinated cyber-attacks. The analysis determined that while the PHX’s failures cause disruption to the reactor, they did not cause any safety limits to be exceeded because of the plant design, including passive safety features. Further development of ARCADE will enable rigorous, repeatable, and automated cyber-physical analysis of advanced reactor control systems. These efforts will also help reduce regulatory uncertainty by presenting similar types of cybersecurity analyses in a common format, driving standard approaches and reporting.

Abstract not provided.

Abstract not provided.

Computational and mathematical models are essential to understanding complex systems and phenomena. However, when developing such models, limited knowledge and/or resources necessitates the use of simplifying assumptions. It is therefore crucial to quantify the impact of such simplifying assumptions on the reliability and accuracy of resulting model predictions. This work develops a first-of-its-kind approach to quantify the impact of physics modeling assumptions on predictions. Here, we leverage the emerging field of model-form uncertainty (MFU) representations, which are parameterized modifications to modeling assumptions, in combination with grouped Sobol’ indices to quantitatively measure an assumption’s importance. Specifically, we compute the grouped Sobol’ index for the MFU representation’s parameters as a single importance measure of the assumption for which the MFU representation characterizes uncertainty. To ensure this approach is robust to the subjective choice of how to parameterize a MFU representation, we establish bounds for the difference between sensitivity results for two different MFU representations based on differences in model prediction statistics. The capabilities associated with this approach are demonstrated on three exemplar problems: an upscaled subsurface contaminant transport problem, ablation modeling for hypersonic flight, and nuclear waste repository modeling. We found that our grouped approach is able to assess the impact of modeling assumptions on predictions and offers computational advantages over classical Sobol’ index computation while providing more interpretable results.

Abstract not provided.

This primary purpose of this project was to evaluate alternative gas mixtures to sulfur hexafluoride (SF6) developed for high voltage power delivery applications for use in high voltage spark gap switches. These SF6 alternatives lower global warming potential emissions and enable improvements to the pressure-voltage design space. A combined experimental, computational, and theoretical study was used to quantify the impact of persistent breakdown products on the breakdown distribution of SF6-replacement gas mixtures. Viable SF6 replacements suitable for use in spark gap switches were studied to enable performance and agility improvements for next-generation pulsed power research relevant to national security missions. Experimental campaign included establishing parameters of switch gases as function of concentration. Various concentrations and pressures were tested for trends in breakdown voltage, repeatability, and durability, and breakdown constituents. A zero-dimensional plasma global model was used to simulate the plasma arc decay and recombination process in spark-gap switches relevant to the Z machine. Finally, a complete and consistent set of electron-neutral collision cross-sections for the novel insulating gas C4F7N is reported.

Modeling of phenomena such as anomalous transport via fractional-order differential equations has been established as an effective alternative to partial differential equations, due to the inherent ability to describe large-scale behavior with greater efficiency than fully resolved classical models. In this review article, we first provide a broad overview of fractional-order derivatives with a clear emphasis on the stochastic processes that underlie their use. We then survey three exemplary application areas — subsurface transport, turbulence, and anomalous materials — in which fractional-order differential equations provide accurate and predictive models. For each area, we report on the evidence of anomalous behavior that justifies the use of fractional-order models, and survey both foundational models as well as more expressive state-of-the-art models. We also propose avenues for future research, including more advanced and physically sound models, as well as tools for calibration and discovery of fractional-order models.

In polymer-filled granular composites, damage may develop in mechanical loading prior to material failure. Damage mechanisms such as microcracking or plastic deformation in the binder phase can substantially alter the material's mesostructure. For energetic materials, such as solid propellants and plastic bonded explosives, these mesostructural changes can have far reaching effects including degraded mechanical properties, potentially increased sensitivity to further insults, and changes in expected performance. Unfortunately, predicting damage is nontrivial due to the complex nature of these composites and the entangled interactions between inelastic mechanisms. In this work, we assess the current literature of experimental knowledge, focusing on the pressure-dependent shear response, and propose a simple simulation framework of bonded particles to study four limiting-case material formulations at both meso- and macro-scales. To construct the four cases, we systematically vary the relative interfacial strength between the polymer binder and granular filler phase and also vary the polymer's glass transition temperature relative to operating temperature which determines how much the binder can plastically deform. These simulations identify key trends in global mechanical response, such as the emergence of strain hardening or softening regimes with increasing pressure which qualitatively resemble experimental results. By quantifying the activation of different inelastic mechanisms, such as bonds breaking and plastically straining, we identify when each mechanism becomes relevant and provide insight into potential origins for changes in mechanical responses. The locations of broken bonds are also used to define larger, mesoscopic cracks to test various metrics of damage. We primarily focus on triaxial compression, but also test the opposite case of triaxial extension to highlight the impact of Lode angle on mechanical behavior.

Abstract not provided.

Abstract not provided.

Deep learning (DL) models have enjoyed increased attention in recent years because of their powerful predictive capabilities. While many successes have been achieved, standard deep learning methods suffer from a lack of uncertainty quantification (UQ). While the development of methods for producing UQ from DL models is an active area of current research, little attention has been given to the quality of the UQ produced by such methods. In order to deploy DL models to high-consequence applications, high-quality UQ is necessary. This report details the research and development conducted as part of a Laboratory Directed Research and Development (LDRD) project at Sandia National Laboratories. The focus of this project is to develop a framework of methods and metrics for the principled assessment of UQ quality in DL models. This report presents an overview of UQ quality assessment in traditional statistical modeling and describes why this approach is difficult to apply in DL contexts. An assessment on relatively simple simulated data is presented to demonstrate that UQ quality can differ greatly between DL models trained on the same data. A method for simulating image data that can then be used for UQ quality assessment is described. A general method for simulating realistic data for the purpose of assessing a model’s UQ quality is also presented. A Bayesian uncertainty framework for understanding uncertainty and existing metrics is described. Research that came out of collaborations with two university partners are discussed along with a software toolkit that is currently being developed to implement the UQ quality assessment framework as well as serve as a general guide to incorporating UQ into DL applications.

We explore the use of reduced physics models for efficient kinetic particle simulations of space charge limited (SCL) emission in inner magnetically insulated transmission lines (inner MITLs), with application to Sandia National Laboratories' Z machine. We propose a drift kinetic (guiding center) model of electron motion in place of a fully kinetic model and electrostatic-magnetostatic fields in place of electromagnetic fields. The validity of these approximations is suggested by the operational parameters of the Z machine, namely, current pulse lengths of order 100 ns compared with Larmor periods typically smaller than 10-11 s, typical Larmor radii of a few (tens) of microns (magnetic fields of tens to hundreds of Tesla) compared with MITL dimensions of a few centimeters, and transient time of light waves along the inner MITL of order a fraction of a nanosecond. Guiding center orbits eliminate the fast electron gyromotion, which enables the use of tens to hundreds of times larger time steps in the numerical particle advance. Electrostatic-magnetostatic fields eliminate the Courant-Friedrichs-Lewy (CFL) numerical stability limit on the time step and allow the use of higher grid resolutions or, alternatively, larger time steps in the fields advance. Overall, potential computational cost savings of tens to hundreds of times exists. The applicability of the reduced physics models is examined on two problems. First, in the simulation of space charge limited emission of electrons from the cathode surface due to high electric fields in a radial inner MITL geometry with a short load. In particular, it is shown that a drift kinetic-based particle-in-cell (PIC) model with electrostatic-magnetostatic fields is able to accurately reproduce well-known physics of electron vortex formation, spatially and temporally. Second, deeper understanding is gained of the mechanism behind vortex formation in this MITL geometry by considering an exemplar problem of an electron block of charge. This simpler setup reveals that the main mechanism of vortex formation can be attributed to pure drift motion of the electrons, that is, the (fully kinetic) gyromotion of the electrons is inessential to the process. This exemplar problem also suggests a correlation of the spatial dimensions of vortices to the thickness of the electron layer, as observed in SCL simulations. It also confirms that the electromagnetic nature of the fields does not play an essential role. Finally, an improved hybrid fully kinetic and drift kinetic model for electron motion is proposed, as means of capturing finite Larmor radius (FLR) effects; the particular FLR physics that is missed by the drift kinetic model is the particle-wall interaction. By initializing SCL emitted electrons as fully kinetic and later transitioning them to drift kinetic, according to simple criteria, the accuracy of SCL simulations can be improved, while preserving the potential for computational efficiency.

This report summarizes the measurement and analysis results of initial performance demonstrator wafers patterned at LumArray in January 2023 in accordance with the project Statement of Work.

Previous studies have shown that selective laser flash sintering (SLFS) can be initiated in dielectrics that exhibit ionic or electronic conduction at high temperature. These materials required high laser powers to reach the temperatures where electrical conduction is sufficient to initiate SLFS. In this study, SLFS in lanthanum chromite (LC), an intrinsic electronic conductor with high conductivity, and lanthanum strontium chromite (LSC), which is doped to further increase electronic conductivity, were investigated with a focus on understanding the initiation mechanisms. Results show that the initiation of SLFS in LC and LSC occurs when electronic charge carriers are activated and flow to the electrode where the current is measured. A combination of carriers produced at the electrode, temperature-activated intrinsic charge carriers, and extrinsic charge carriers present in LSC due to doping are responsible for the facile initiation of SLFS.

Achieving brain-like efficiency in computing requires a co-design between the development of neural algorithms, brain-inspired circuit design, and careful consideration of how to use emerging devices. The recognition that leveraging device-level noise as a source of controlled stochasticity represents an exciting prospect of achieving brain-like capabilities in probabilistic neural algorithms, but the reality of integrating stochastic devices with deterministic devices in an already-challenging neuromorphic circuit design process is formidable. Here, we explore how the brain combines different signaling modalities into its neural circuits as well as consider the implications of more tightly integrated stochastic, analog, and digital circuits. Further, by acknowledging that a fully CMOS implementation is the appropriate baseline, we conclude that if mixing modalities is going to be successful for neuromorphic computing, it will be critical that device choices consider strengths and limitations at the overall circuit level.

Code verification plays an important role in establishing the credibility of computational simulations by assessing the correctness of the implementation of the underlying numerical methods. In computational electromagnetics, the numerical solution to integral equations incurs multiple interacting sources of numerical error, as well as 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 combined-field integral equation. We demonstrate the effectiveness of these approaches for cases with and without coding errors.

A discussion of systems like ChatGPT, what their legal issues may be, how they are affecting society, and what the ethical considerations of their existence and use are.

We develop methods for previously introduced NMR hardware that enables sensitive, noninvasive spectroscopic study of coin cell batteries, by overcoming the magnetic field distortions intrinsic to the battery limited chemical resolution. Specifically, we employ an interferometric NMR technique to overcome field heterogeneity, reducing peak widths by two orders of magnitude and obtaining chemically resolved solution-state spectra.

Abstract not provided.

Abstract not provided.

Distributed acoustic sensing (DAS) has a demonstrated potential for wide-scale and continuous in situ monitoring of near-surface environmental and anthropogenic processes. DAS is attractive for development as a multi-geophysical observatory due to the prevalence of existing fiber infrastructure in regions with environmental, cultural, or strategic significance. To evaluate the efficacy of this technology for monitoring of polar environmental processes, we collected DAS data from a 37-km long section of seafloor telecommunications fiber located on the continental shelf of the Beaufort Sea, Alaska. This experiment spanned eight, one-week, seasonally-distributed periods across two years. This was the first ever deployment of seafloor DAS beneath sea ice, and the first deployment in any marine environment to span multiple seasons. We recorded a variety of environmental and anthropogenic signals with demonstrable utility for the study of sea ice dynamics and tracking of ocean vessels and ice-traversing vehicles.

Experimental studies and ab initio quantum chemistry calculations were combined to investigate the process by which a Fenton reaction breaks down polystyrene sulfonate. The experimental results show that both molecular weight reduction and loss of aromaticity occur nearly simultaneously, a finding that is supported by the calculations. The results show that more than half of the material is broken down to low molecular weight compounds (< 500 g/mol) with two molar equivalents of H2O2 per styrene monomer. The calculations provide insights into the reaction pathways and indicate that at least two hydroxyl radicals are required to cleave backbone C–C bonds or to eliminate aromaticity. The calculations also show that, of the aromatic carbons, hydroxyl radical is most likely to add to the carbon bonded to sulfur. This finding explains the loss of hydrogen sulfite anion early in the process and also the efficient reduction of Fe(III) to Fe(II) through semiquinone formation. Taken together the experimental and computational results indicate that the reaction is very efficient and that very little H2O2 is lost to unproductive reactions. This high efficiency is attributed to the close association of Fe atoms with the sulfonate group such that hydroxyl radicals are generated near the polymer chains.

Scientific software (SciSoft) is complex, often containing a mixture of production capabilities co-mingled with features under active research and development. Furthermore, SciSoft is often developed over decades by non-computer scientists who may not have a strong background in or prioritize software architecture design, testing, and quality (e.g., test coverage). These conditions lead to difficulty in understanding which software components or functions implement what user-facing features and therefore those features’ software quality pedigree. This lack of understanding poses challenges in assessing readiness and credibility of user features, and often relies on a SciSoft subject matter expert’s (SME) laborious investigation and assertion. This final report of a one-year Computing and Information Sciences Lab Directed Research and Development project presents a general framework for modeling SciSoft architecture as a direct relationship between user features and the software components/functions that implement them. Our approach leverages automated labeling of the SciSoft’s regression test suite and employs machine learning algorithms to construct the architecture model. We demonstrate this framework on the Solid Mechanics component of the SIERRA multi-physics engineering analysis suite developed at Sandia National Laboratories.

The SALSA3D tomographic model provides a crucial community resource for improving the quality (in terms of both accuracy and precision) of predictions of the traveltimes of seismic waves and therefore improving our ability to locate anthropogenic or natural seismic events. Constructing the requisite tomographic model requires addressing the challenges implied by a massive and growing dataset of traveltime measurements. This study explores one approach to tackle this challenge: the use of summary rays, which average traveltime measurements from sources within evenly spaced cells, thereby eliminating redundant data.

This report summarizes FY23 activities to improve mechanistic source term modeling for MSR concepts. Relevant MELCOR capability enhancements made during FY23 are summarized including development of a flexible python-based EOS generator (MELEOS), porous domain modeling capabilities for validation applications, and development of a MELCOR model for the LSTL facility in anticipation of upcoming molten salt experiments.

Abstract not provided.

This report documents the generation of a skeletal chemical reaction mechanism for use with hemispherical pentaerythritol tetranitrate charges. Skeletal mechanisms can substantially reduce computation time while maintaining accuracy. The methodology within uses faster running sample simulations to build a representative thermodynamic state space. These thermodynamic states are used with a constant-volume reactor analysis and a reaction flow analysis to remove unimportant species and reactions from a full chemical reaction mechanism. For the given test case, this results in a 6x speedup in computation time for directly comparable simulations in 2D axisymmetric simulations. We see a 30x speedup in simulations in 3D Cartesian coordainates when compared to a prior full kinetics simulation. There is strong agreement between temperature and species mass fraction profiles between the full and skeletal chemical reaction mechanisms. These methodologies can be applied to any explosive, given the availability of sample simulations.

Abstract not provided.

The Canada-US Blended Cyber-Physical Exercise was a successful, first of its kind, multiorganization and multi-laboratory exercise that culminated years of complex system development and planning. The project aimed to answer three driving research questions, (1) How do cyberattacks support malicious acts leading to theft or sabotage [at a nuclear site]? (2) What are aspects of an effective combined cyber-physical response? (3) How to evaluate effectiveness of that response? Which derived the following primary objectives, 1. The May 2023 Cyber-Physical Exercise shall present a cyber-attack scenario that supports malicious acts leading to theft or sabotage. 2. The May 2023 Cyber-Physical Exercise shall define aspects of an effective combined cyber-physical response. 3. Analysis of the May 2023 Cyber-Physical Exercise shall evaluate the effectiveness of the incident response against pre-established exercise evaluation criteria. 4. Analysis of the May 2023 Cyber-Physical Exercise shall assess the effectiveness of the evaluation criteria itself. 5. Exercises shall be performed in a real-life environment. The team believes these objectives were met, and the evidence will be presented in this report. Due to the novelty of the exercise, there were several lessons learned that will be presented in this report.

During this LDRD project, our team developed a technology which enables the fabrication of novel nanostructures replicating seashell – “nature’s toughest material”. The resulting coatings exhibit high thermal stability up to 1650°C, which exceeds the hardness of Spectra® by ~44%, as well as the compressive strength of aluminum by ~57%. Coatings made with this technology are stronger, environmentally friendly, more sustainable, and more versatile than other comparable materials. Beryllium wafers, the current, most favorable shielding material in terms of thermal and mechanical properties, are very toxic and cost hundreds of times more than the new material developed in this project. The coatings on silicon wafer and stainless steel, respectively, have been tested as ride-along on the Z machine and clearly outperform the bare substrate. Use of this technology will have a profound global impact for pulsed power and fusion energy development, debris mitigation for spacecraft and satellites, durability of drill bits used in deep well drilling and tunnel boring operations, thermal protection of aircraft and manned spacecraft, and various other thermal and mechanical protection applications.

Abstract not provided.

Experimental measurements of low density plasmas forming in Magnetically Insulated Transmission Line (MITL) regions are desired to improve our understanding of current loss and power flow. Therefore, a new optical interferometer diagnostic was commissioned via this LDRD project. To measure the expected 10¹³ - 10¹⁷ cm^-3 electron densities inside the 0.5 - 6 mm Anode-Cathode (A-K) gaps, a colinear SHOP interferometer diagnostic was constructed. The diagnostic was initially fielded on the University of New Mexico (UNM) Helicon-Cathode (HelCat) plasma device which provided a highly repeatable and well understood plasma source for which the colinear SHOP interferometer’s functionality could be verified and measured. Utilizing the highly repeatable plasma source and shot averaging, the interferometer was able to achieve an areal density sensitivity of 1×10¹⁴ cm^-2. This work at UNM lead to a Review of Scientific Instruments (RSI) publication [20], DOI:10.1063/5.0101687. After the diagnostic’s capability was proven at UNM, the colinear SHOP interferometer was commissioned for use on the Sandia National Laboratories (SNL) Mykonos accelerator. Here, it provided the first temporal areal density measurements of plasma formation in a parallel plate MITL. The diagnostic was able to achieve a single shot (no multi-shot averaging like at UNM) areal density sensitivity of 1×10¹⁵ cm^-2 along a ~ 2mm probing path length, which provided adequate capability to conduct fundamental physics research of MITL plasma formation. CHICAGO and ALEGRA simulations support the diagnostics experimental findings. More experimental and computational work will continue, likely leading to another publication(s). The smaller scale Mykonos accelerator work has also provided justification that the colinear SHOP interferometer is a capable diagnostic for measuring plasma areal densities in the inner MITL and convolute regions of larger TW-class accelerators like SNL’s Z machine.

In this brief report we document algorithmic choices and updates to our code related to the earthquake relocation portion of our tomographic imaging algorithm. We show results of these improvements by relocating over 40,000 events located within 20-30 km of the Rock Valley Direct Comparison (RV/DC) site using both absolute and differential arrival times within the context of two different 3-D Earth models. Accurate hypocentral locations and Earth models are important to the ultimate goals of the RV/DC program, which will co-locate a chemical explosion with a shallow earthquake within Rock Valley, southern Nevada, to investigate differences between the source types and improve our analysis algorithms for both types (Snelson et al., 2022). Our improvements to our relocation algorithms comprise just one step toward achieving these goals

This is the presentation, that goes with a previously-approved, Information-Released abstract of the same title,

In this work, thermogravimetric analysis (TGA) was performed on samples of a carbon fiber epoxy composite, a glass fiber epoxy composite, and a mixed carbon fiber/glass fiber epoxy composite, as well on each constituent material (polymer epoxy, carbon fibers and glass fibers). TGA was conducted for heating rates from 1-20 C/min with purified purge gases of nitrogen and dry air. For the fiberglass composite, we find that ~70% of the material remains after heating in air to 1200 C. For the carbon fiber epoxy composite, we observe greater mass loss as the carbon fibers can oxidize, leaving little material by the end of the test. The mixed composite, which has a 2:1 ratio of glass fibers to carbon fibers, experienced a total mass loss between the two other composites. By determining the relationship between the thermal decomposition of a composite material and its constituent materials, we can predict the fire behavior of novel composites during the material design phase.

Although visualizations are a useful tool for helping people to understand information, they can also have unintended effects on human cognition. This is especially true for uncertain information, which is difficult for people to understand. Prior work has found that different methods of visualizing uncertain information can produce different patterns of decision making from users. However, uncertainty can also be represented via text or numerical information, and few studies have systematically compared these types of representations to visualizations of uncertainty. We present two experiments that compared visual representations of risk (icon arrays) to numerical representations (natural frequencies) in a wildfire evacuation task. Like prior studies, we found that different types of visual cues led to different patterns of decision making. In addition, our comparison of visual and numerical representations of risk found that people were more likely to evacuate when they saw visualizations than when they saw numerical representations. These experiments reinforce the idea that design choices are not neutral: seemingly minor differences in how information is represented can have important impacts on human risk perception and decision making.

The primary purpose of this document is to outline the progress made on the LDRD titled “Identifying and Explaining Anomalous Activity in Surveillance Video with Compression Algorithms” in FY22 and FY23. In this LDRD, we explored the usage of compression-based analytics to identify anomalous activity in video. We developed a novel algorithm, Spatio-Temporal N-Gram PPM (STNG PPM) that accounts for spatially and temporally aware anomalies in video. We extracted features using motions vectors from video as well as operating on the raw features. STNG PPM is comparable to many deep learning approaches but does not require specialized hardware (GPUs) to run efficiently. We also examine the evaluation metrics and propose novel measures addressing faults in the current evaluation measures.

Abstract not provided.

In this report, we document the process related to developing a regional geologic model of a 605 x 1334 km area centered around Utah and encompassing surrounding states. This model is developed to test the effect that composition of a model has on the generation of synthetic data with the intent of using this information to improve upon full waveform moment tensor inversions. We compare observed data from three seismic events and five stations to the synthetic data generated by a preliminary model derived from a geologic framework model (GFM) developed by the USGS. The synthetic data and observed data comparisons indicate that our preliminary model performs well at smaller offset distances in the northern and central sections of the model. However, the southern stations consistently display synthetic data P- and S-wave arrival times that do not match the observed data arrival times, indicating that the velocity structure of the southern part of the model especially is inaccurate.

Laser powder bed fusion (LPBF) Additive Manufacturing (AM) has the potential to enable the production of components with novel designs and material properties unachievable otherwise. However, process repeatability is a challenge, making qualification ill-defined and greatly reducing the utility of what could be an important manufacturing technology. In this work, a combination of modeling, uncertainty quantification (UQ), and experimentation are used in an effort to predict and bound the range of possible outcomes of the LPBF process. Quantities of interest predicted are melt pool dimensions, microstructure features, and mechanical distortions. A combination of high fidelity thermal-fluid models, microstructure growth models, and reduced fidelity, rapid thermal and mechanical models are used. Uncertainty propagation techniques are used to predict probability distributions of quantities of interest from estimates of process uncertainties. Repeated experiments are done to quantify observed probability distributions and compared to predicted distributions to determine if predictions are precise and accurate. Novel modeling methods are microstrucutre characterization techniques are also discussed. It is found that high fidelity models do a generally good job bounding experimentally observed melt pool morphologies for both bead-on-plate and powder bed cases. Microstructure models are able to bound a number of experimentally observed microstructure statistics, but with low precision due to challenges with calibrating the microstructure growth model parameters. A developed modified inherent strain distortion model does not accurately predict observed distortions. A lumped laser distortion model shows promise in being both accurately and precisely bounding observed outcomes from the deflection comb build, but requires further evaluation on more builds and geometries.

Quantum point contacts (QPC) are the building blocks of quantum dot qubits and semiconducting quantum electrical metrology circuits. QPCs also make highly sensitive electrical amplifiers with the potential to operate in the quantum-limited regime. Though the inherent operational bandwidth of QPCs can eclipse the THz regime, the impedance mismatch with the external circuitry limits the operational frequency to a few kHz. Lumped-element impedance-matching circuits are successful only up to a few hundreds of MHz in frequency. QPCs are characterised by a complex impedance consisting of quantized resistance, capacitance, and inductance elements. Characterising the complex admittance at higher frequencies and understanding the coupling of QPC to other circuit elements and electromagnetic environments will provide valuable insight into its sensing and backaction properties. In this work, we couple a QPC galvanically to a superconducting stub tuner impedance matching circuit realised in a coplanar waveguide architecture to enhance the operation frequency into the GHz regime and investigate the electrical amplification and complex admittance characteristics. The device, operating at ~1.96 $GHz$ exhibits a conductance sensitivity of 2.92 X 10^-5(e²/h)/$\sqrt{Hz}$ with a bandwidth of 13 $MHz$. Besides, the RF reflected power unambiguously reveals the complex admittance characteristics of the QPC, shining more light on the behaviour of quantum tunnel junctions at higher operational frequencies.

Securing satellite groundstations against cyber-attacks is vital to national security missions. However, these cyber threats are constantly evolving. As vulnerabilities are discovered and patched, new vulnerabilities are discovered and exploited. In order to automate the process of discovering existing vulnerabilities and the means to exploit them, a reinforcement learning framework is presented in this report. We demonstrate that this framework can learn to successfully navigate an unknown network and detect nodes of interest despite the presence of a moving target defense. The agent then exfiltrates a file of interest from the node as quickly as possible. This framework also incorporates a defensive software agent that learns to impede the attacking agents progress. This setup allows for the agents to work against each other and improve their abilities. We anticipate that this capability will help uncover unforeseen vulnerabilities and the means to mitigate them. The modular nature of the framework enables users to swap out learning algorithms and modify the reward functions in order to adapt the learning tasks to various use cases and environments. Several algorithms, viz., tabular Q learning, deep Q networks, proximal policy optimization, advantage actor-critic, generative adversarial imitation learning, are explored for the agents and the results highlighted. The agent learns to solve the tasks in a light-weight abstract environment. Once the agent learns to perform sufficiently well, it can be deployed in a minimega virtual machine environment (or a real network) with wrappers that map abstract actions to software commands. The agent also uses a local representation of the actions called a ‘slot-mechanism’. This allows the agent to learn in a certain network and generalize it to different networks. The defensive agent learns to predict the actions taken by an offensive agent and uses that information to anticipate the threat. This information can then either be used to raise an alarm or to take actions to thwart the attack. We believe that with the appropriate reward design, a representative environment, and action set, this framework can be generalized to tackle other cybersecurity tasks. By sufficiently training these agents, we can anticipate vulnerabilities leading to robust future designs. We can also deploy automated defensive agents that can help secure satellite groundstation and their vital national security missions.

This paper describes a process for forming a buried field shield in GaN by an etch-and-regrowth process, which is intended to protect the gate dielectric from high fields in the blocking state. GaN trench MOSFETs made at Sandia serve as the baseline to show the limitations in making a trench gated device without a method to protect the gate dielectric. Device data coupled with simulations show device failure at 30% of theoretical breakdown for devices made without a field shield. Implementation of a field shield reduces the simulated electric field in the dielectric to below 4 MV/cm at breakdown, which eliminates the requirement to derate the device in order to protect the dielectric. For realistic lithography tolerances, however, a shield-to-channel distance of 0.4 μm limits the field in the gate dielectric to 5 MV/cm and requires a small margin of device derating to safeguard a long-term reliability and lifetime of the dielectric.

Abstract not provided.

The wave energy control competition established a benchmark problem which was offered as an open challenge to the wave energy system control community. The competition had two stages: In the first stage, competitors used a standard wave energy simulation platform (WEC-Sim) to evaluate their controllers while, in the second stage, competitors were invited to test their controllers in a real-time implementation on a prototype system in a wave tank. The performance function used was based on converted energy across a range of standard sea states, but also included aspects related to economic performance, such as peak/average power, peak force, etc. This paper compares simulated and experimental results and, in particular, examines if the results obtained in a linear system simulation are borne out in reality. Overall, within the scope of the device tested, the range of sea states employed, and the performance metric used, the conclusion is that high-performance WEC controllers work well in practice, with good carry-over from simulation to experimentation. However, the availability of a good WEC mathematical model is deemed to be crucial.

Using DTW for HPC / LDMS data clustering

Fracture prognosis and characterization efforts require knowledge of crack tip position and the Stress Intensity Factors (SIFs) acting in the vicinity of the crack. Here, we present an efficient numerical approach to infer both of these characteristics under a consistent theoretical framework from noisy, unstructured displacement data. The novel approach utilizes the separability of the asymptotic linear elastic fracture mechanics fields to expedite the search for crack tip position and is particularly useful for noisy displacement data. The manuscript begins with an assessment of the importance of accurately locating crack tip position when quantifying the SIFs from displacement data. Next, the proposed separability approach for quickly inferring crack tip position is introduced. Comparing to the widely used displacement correlation approach, the performance of the separability approach is assessed. Cases involving both noisy data and systematic deviation from the asymptotic linear elastic fracture mechanics model are considered, e.g. inelastic material behavior and finite geometries. An open source python implementation of the proposed approach is available for use by those doing field and laboratory work involving digital image correlation and simulations, e.g. finite element, discrete element, molecular dynamics and peridynamics, where the crack tip position is not explicitly defined.

The modern global economy relies heavily on carbon-based products that are derived from petroleum, which presents sustainability, resource management, and greenhouse gas exacerbated climate change challenges. Due to these challenges, there is the need for a global industrial transition towards green and sustainable production. Microbial production of valuable chemicals from renewable biomass represents one promising route. However, high-volume low-value products such as commodity chemicals are still difficult to make profitable. One fundamental bottleneck is a waste of more than 1/3 of the feedstock carbon as CO₂ in the fermentation process. Here the project focuses on fundamentally reconfiguring the metabolism to reduce CO₂ loss in central metabolic pathways thereby also improving bioproduct yields. Here we present technologies to prevent CO₂ loss and balance reducing equivalents within the cell to enable complete conversion of glucose from renewable feedstocks into bioproducts.

Frontal polymerization involves the propagation of a thermally driven polymerization wave through a monomer solution to rapidly generate high-performance polymeric materials with little energy input. The balance between latent catalyst activation and sufficient reactivity to sustain a front can be difficult to achieve and often results in systems with poor storage lives. This is of particular concern for frontal ring-opening metathesis polymerization (FROMP) where gelation occurs within a single day of resin preparation due to the highly reactive nature of Grubbs-type catalysts. In this report we demonstrate the use of encapsulated catalysts to provide remarkable latency to frontal polymerization systems, specifically using the highly active dicyclopentadiene monomer system. Negligible differences were observed in the frontal velocities or thermomechanical properties of the resulting polymeric materials. FROMP systems with encapsulated catalyst particles are shown with storage lives exceeding 12 months and front rates that increase over a well-characterized 2 month period. Moreover, the modularity of this encapsulation method is demonstrated by encapsulating a platinum catalyst for the frontal polymerization of silicones by using hydrosilylation chemistry.

This study presents a deep learning based methodology for both remote sensing and design of acoustic scatterers. The ability to determine the shape of a scatterer, either in the context of material design or sensing, plays a critical role in many practical engineering problems. This class of inverse problems is extremely challenging due to their high-dimensional, nonlinear, and ill-posed nature. To overcome these technical hurdles, we introduce a geometric regularization approach for deep neural networks (DNN) based on non-uniform rational B-splines (NURBS) and capable of predicting complex 2D scatterer geometries in a parsimonious dimensional representation. Then, this geometric regularization is combined with physics-embedded learning and integrated within a robust convolutional autoencoder (CAE) architecture to accurately predict the shape of 2D scatterers in the context of identification and inverse design problems. An extensive numerical study is presented in order to showcase the remarkable ability of this approach to handle complex scatterer geometries while generating physically-consistent acoustic fields. The study also assesses and contrasts the role played by the (weakly) embedded physics in the convergence of the DNN predictions to a physically consistent inverse design.

Decarbonization efforts highlight hydrogen as an attractive alternative to fossil fuels, but its tendency to embrittle structural metals demands careful consideration when designing hydrogen infrastructure. Moreover, the mechanisms by which hydrogen degrades these materials are still being elucidated. The current work develops new computational tools to quantify the different contributions of hydrogen to the energy barrier of cross-slip, a key deformation mechanism. Novel features are implemented to a line tension model, which include the use of non-singular dislocation interactions, character-dependent dislocation energies and simulations of the constriction configurations. A new molecular dynamics technique is developed to calculate the interaction energy between the partials of a dissociated dislocation via fixing the centers of mass of the regions below and above the Shockley partials and performing time-averaged calculations. Hydrogen is found to impact the stacking fault width of dislocations in different ways depending on their characters: it decreases for dislocations with a character θ > 30°, remains unchanged for θ = 30° and increases for θ < 30°. The latter regime is a newly identified mechanism by which hydrogen inhibits cross-slip. Moreover, formation of nano-hydrides is predicted to occur around screw dislocations for high hydrogen concentrations, a phenomenon previously identified only in dislocations with an edge component. If nano-hydrides develop, their influence extending the equilibrium stacking fault width and increasing both the constriction and cross-slip energy barriers dominate over all other hydrogen contributions. The theory and tools developed will pave the way towards a comprehensive understanding of hydrogen-dislocation interactions in structural metals.

The carbon phase diagram is rich with polymorphs which possess very different physical and optical properties ideal for different scientific and engineering applications. An understanding of the dynamically driven phase transitions in carbon is particularly important for applications in inertial confinement fusion, as well as planetary and meteorite impact histories. Experiments on the Z Pulsed Power Facility at Sandia National Laboratories generate dynamically compressed high-pressure states of matter with exceptional uniformity, duration, and size that are ideal for investigations of fundamental material properties. X-ray diffraction (XRD) is an important material physics measurement because it enables direct observation of the strain and compression of the crystal lattice, and it enables the detection and identification of phase transitions. Several unique challenges of dynamic compression experiments on Z prevent using XRD systems typically utilized at other dynamic compression facilities, so novel XRD diagnostics have been designed and implemented. We performed experiments on Z to shock compress carbon (pyrolytic graphite) samples to pressures of 150–320 GPa. The Z-Beamlet Laser generated Mn-Heα (6.2 keV) X-rays to probe the shock-compressed carbon sample, and the new XRD diagnostics measured changes in the diffraction pattern as the carbon transformed into its high-pressure phases. Quantitative analysis of the dynamic XRD patterns in combination with continuum velocimetry information constrained the stability fields and melting of high-pressure carbon polymorphs.

This project pursued a novel strategy for incorporating multiple qubits per ion into ion-trap based quantum computing (ITQC) involving Ca⁺ and Ba⁺. By forming molecular complexes of these cations with molecular-scale cages, we hypothesized that molecular energy levels could be incorporated into quantum computing while retaining key properties of the atomic ions intact. We experimented with a variety of molecular cages and found that Na⁺, K⁺, Rb⁺, Ca²⁺, Sr²⁺, and Ba²⁺ could be captured and brought into the gas phase efficiently by imbedding them inside [2.2.2]-benzocryptand. IR and UV spectra of these cage complexes are sensitive to the size and charge state of the ion, reporting on the structures and binding properties of the cage complexes. UV photofragmentation of the Ba²⁺-Acetate^-1-BzCrypt complex produces Ba⁺-BzCrypt, the complex targeted for exploration in the original hypothesis. Follow-on funding is needed to pursue the spectroscopy of this complex as a target for ITQC.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The understanding of power flow plasmas is important as we look towards next generation pulsed power (NGPP) as current losses could prohibit the goals of that facility. Therefore, it is important to have accurate diagnostics of the plasma parameters on the current machines, which can be used to help inform and improve simulations. Having these plasma parameters will help validate models and simulations to provide confidence when they are expanded to conditions relevant to NGPP. One important plasma parameter that can be measured is the electron density, which can be measured by photonic Doppler velocimetry (PDV). A PDV system has several key advantages over other interferometers by measuring relatively low densities (> 1 × 10¹⁵ cm^-2) with both spatial and temporal resolution. Experiments were performed on the Mykonos pulsed power machine, which is a 1 MA sub scale machine in which recent platforms have been developed to explore current densities relevant to the inner magnetically insulated transmission line (MITL) on the Z machine. Experiments were performed on two different platforms, the thin foil platform and the Mykonos parallel plate platform (MP3). In addition, a combination of both single-point and multi-point measurements were used. The single-point measurements proved to be very promising, providing a clear increase in density at about 70 ns into the current rise on thin foil experiments up to about 5 × 10¹⁷ cm^-3 before the probe stopped providing signal. While we did also see returns from multi-point measurements on both platforms, the signals were not as easy to interpret due to strong background effects. However, they do show initial promise for this diagnostic to measure density at several points across a 1 mm gap. These measurements provide insights in how to improve the diagnostic so that it can provide useful information on power flow relevant experiments.

This report documents the generation of a mechanism to predict the inclusion of carbon soot particles in a high explosive flow. The mechanism includes gasification and oxidation reactions, formation, sublimation, radiation, and agglomeration. Each part of the mechanism is derived from properties in the literature. The influence of each part of the mechanism is explored using simple, example simulations consisting of a 12 mm diameter 2,4,6-Trinitrotoluene charge detonated in ambient air. The mechanism has not been quantitatively compared to experiments. Additional efforts will be required to tune and validate it, which will require continued advancements in experimental diagnostics and simulation techniques.