Publications

Abstract not provided.

Abstract not provided.

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.

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.

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.

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.

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.

Abstract not provided.

Abstract not provided.

Abstract not provided.

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.

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.

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.

Abstract not provided.

Abstract not provided.

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.

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.

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 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.

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.

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.

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 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.