Publications

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.

Abstract not provided.

CO₂-neutral ammonia production with concentrated solar technology is theoretically possible based on advanced solar thermochemical looping technology. The parametric analysis points to the re-oxidation temperature and the H₃ yield as the most influential parameters in the energy balance. The cycle time and the nitride cost are the most influential parameters on the CAPEX. The techno-economics analysis shows the potential of the plant to achieve a target price <125 $\$$/tonne.

Solar Thermal Ammonia Production has the potential to synthesize ammonia in a green, renewable process that can greatly reduce the carbon footprint left by the conventional Haber-Bosch reaction. Co₃Mo₃N has been identified as a potential candidate for ammonia production. It is synthesized via oxide precursor synthesis followed by nitridation under 10% H₂/N₂. The synthesis method can be extended to other candidate nitrides. The Co₃Mo₃N → Co₆Mo₆N reduction is demonstrated on TGA with rapid kinetics. The formation of NH₃ is qualitatively observed, but not quantitatively determined. The material retains crystal structure, but no secondary phases are observed in XRD. Partial re-nitridation back to CMN331 of ~35% of max nitridation is observed. Reaction parameters in TGA differ from experimental conditions in the literature. Experiments at Georgia Tech better mimic re-nitridation conditions with more sensitive, quantitative analytical techniques (GC-MS). The ASU NH₃ synthesis/re-nitridation reactor is under development and will permit experiments (reduction/re-nitridation) under precisely controlled T, pH₂.

The atomic precision advanced manufacturing (APAM) enabled vertical tunneling field effect transistor (TFET) presents a new opportunity in microelectronics thanks to the use of ultra-high doping and atomically abrupt doping profiles. We present modeling and assessment of the APAM TFET using TCAD Charon simulation. First, we show, through a combination of simulation and experiment, that we can achieve good control of the gated channel on top of a phosphorus layer made using APAM, an essential part of the APAM TFET. Then, we present simulation results of a preliminary APAM TFET that predict transistor-like current-voltage response despite low device performance caused by using large geometry dimensions. Future device simulations will be needed to optimize geometry and doping to guide device design for achieving superior device performance.

Abstract not provided.

We present an efficient self-consistent implementation of the Non-Equilibrium Green Function formalism, based on the Contact Block Reduction method for fast numerical efficiency, and the predictor-corrector approach, together with the Anderson mixing scheme, for the self-consistent solution of the Poisson and Schrödinger equations. Then, we apply this quantum transport framework to investigate 2D horizontal Si:P δ-layer Tunnel Junctions. We find that the potential barrier height varies with the tunnel gap width and the applied bias and that the sign of a single charge impurity in the tunnel gap plays an important role in the electrical current.

We propose to demonstrate the feasibility of a solar thermochemical looping technology to produce and store nitrogen (N₂) from air for the subsequent production of ammonia (NH₃) via an advanced two-stage process.

This report includes a compilation of several slide presentations: 1) Interatomic Potentials for Materials Science and Beyond–Advances in Machine Learned Spectral Neighborhood Analysis Potentials (Wood); 2) Agile Materials Science and Advanced Manufacturing through AI/ML (de Oca Zapiain); 3) Machine Learning for DFT Calculations (Rajamanickam); 4) Structure-preserving ML discovery of a quantum-to-continuum codesign stack (Trask); and 5) IBM Overview of Accelerated Discovery Technology (Pitera)

Laser beam directed energy deposition has become an increasingly popular advanced manufacturing technique for materials discovery as a result of the in situ alloying capability. In this study, we leverage an additive manufacturing enabled high throughput materials discovery approach to explore the composition space of a graded Wx(CoCrFeMnNi)100−x sample spanning 0 ≤ x ≤ 21 at%. In addition to microstructural and mechanical characterization, synchrotron high speed x-ray computer aided tomography was conducted on a W20(CoCrFeMnNi)80 composition to visualize melting dynamics, powder-laser interactions, and remelting effects of previously consolidated material. Results reveal the formation of the Fe7W6 intermetallic phase at W concentrations> 6 at%, despite the high configurational entropy. Unincorporated W particles also occurred at W concentrations> 10 at% accompanied by a dissolution band of Fe7W6 at the W/matrix interface and hardness values greater than 400 HV. The primary strengthening mechanism is attributed to the reinforcement of the Fe7W6 and W phases as a metal matrix composite. The in situ high speed x-ray imaging during remelting showed that an additional laser pass did not promote further mixing of the Fe7W6 or W phases suggesting that, despite the dissolution of the W into the Fe7W6 phase being thermodynamically favored, it is kinetically limited by the thickness/diffusivity of the intermetallic phase, and the rapid solidification of the laser-based process.

This progress report describes work performed during FY21 at Sandia National Laboratories (SNL) to assess the localized corrosion performance of canister materials used in the interim storage of spent nuclear fuel (SNF). Of particular concern is stress corrosion cracking (SCC), by which a through-wall crack could potentially form in a canister outer wall over time intervals that are shorter than possible dry storage times. In FY21, modeling and experimental work was performed that further defined our understanding of the potential chemical and physical environment present on canister surfaces at both marine and inland sites. Research also evaluated the relationship between the environment and the rate, extent, and morphology of corrosion, as well as the corrosion processes that occur. Finally, crack growth rate testing under relevant environmental conditions was initiated.

Simulations of several of the end-irradiated cylindrical photoelectron driven cavity experiments (also known as B-Dot cavities) that were fielded during the July 1 through 2, 2020 shot series at the National Ignition Facility are presented in this report with comparisons to experimental measurements. All cavity B-Dots fielded on the second, third, fourth, fifth and seventh shots were simulated using coupled Integrated Tiger Series (ITS) Monte Carlo transport codes and the Electromagnetic Plasmas in Realistic Environments (EMPIRE) electromagnetic particle-in-cell code.

A new method for computing anharmonic thermophysical properties for adsorbates on metal surfaces is presented. Classical Monte Carlo phase space integration is performed to calculate the partition function for the motion of a hydrogen atom on Cu(111). A minima-preserving neural network potential energy surface is used within the integration routine. Two different sampling schema for generating the training data are presented, and two different density functionals are used. The results are benchmarked against direct state counting results by using discrete variable representation. The phase space integration results are in excellent quantitative agreement with the benchmark results. Additionally, both the discrete variable representation and the phase space integration results confirm that the motion of H on Cu(111) is highly anharmonic. The results were applied to calculate the free energy of dissociative adsorption of H2 and the resulting Langmuir isotherms at 400, 800, and 1200 K in a partial pressure range of 0-1 bar. It shows that the anharmonic effects lead to significantly higher predicted surface site fractions of hydrogen.

Since the classical molecular dynamics simulator LAMMPS was released as an open source code in 2004, it has become a widely-used tool for particle-based modeling of materials at length scales ranging from atomic to mesoscale to continuum. Reasons for its popularity are that it provides a wide variety of particle interaction models for different materials, that it runs on any platform from a single CPU core to the largest supercomputers with accelerators, and that it gives users control over simulation details, either via the input script or by adding code for new interatomic potentials, constraints, diagnostics, or other features needed for their models. As a result, hundreds of people have contributed new capabilities to LAMMPS and it has grown from fifty thousand lines of code in 2004 to a million lines today. In this paper several of the fundamental algorithms used in LAMMPS are described along with the design strategies which have made it flexible for both users and developers. We also highlight some capabilities recently added to the code which were enabled by this flexibility, including dynamic load balancing, on-the-fly visualization, magnetic spin dynamics models, and quantum-accuracy machine learning interatomic potentials.

This project explored chemical vapor deposition as a technique to synthesis cubic boron nitride (c-BN) for electronics and coating applications. Current c-BN synthesis techniques are greatly limiting due to requiring high pressure or non-equilibrium energy sources, such as ion bombardment.

Understanding liquid hydrogen tank fluid dynamics is key for modeling liquid hydrogen systems. The tank is the source for nearly all liquid hydrogen systems. Accurate flow modeling out of the tank is needed to predict flows through downstream components. Tank contains liquid and gas that may not be at equilibrium. Questions to be addressed are: Does heat and mass transfer between liquid and vapor affect the flow rate? Is boiling an important consideration? For what conditions is a pressure relief valve (PRV) sufficient to relieve pressure and when is the burst disc needed?

Two-step solar thermochemical cycles based on reversible reactions of SrFeO3−δand (Ba,La)0.15Sr0.85FeO3−δperovskites were considered for air separation. The cycle steps encompass (1) the thermal reduction of SrFeO3−δor (Ba,La)0.15Sr0.85FeO3−δperovskites driven by concentrated solar irradiation and (2) oxidation in air to remove O2and produce N2. Rate limiting mechanisms were examined for both reactions using a combination of isothermal and non-isothermal thermogravimetry for temperature-swings between 673 and 1373 K, heating rates of 10, 20, and 50 K min−1, and O2pressure-swings between 20% O2/Ar and 100% Ar at atmospheric pressure. Evolved O2and associated lag due to transport behavior were measured with gas chromatography and used with measured sample temperatures to predict equilibrium compositions from a compound energy formalism thermodynamic model. Measured and predicted chemical equilibrium changes in deviation from stoichiometry were compared. Rapid chemical kinetics were observed as the samples equilibrated rapidly for all conditions, indicative that heat and mass transfer were the rate limiting mechanisms. The effects of bulk diffusion (or gas diffusion through the bed or pellet) were examined using pelletized and loose powdered samples and determined to have no discernable impact.

Brittle materials, such as cement, compose major portions of built infrastructure and are vulnerable to degradation and fracture from chemo-mechanical effects. Currently, methods of modeling infrastructure do not account for the presence of a reactive environment, such as water, on the acceleration of failure. Here, we have developed methodologies and models of concrete and cement fracture that account for varying material properties, such as strength, shrinkage, and fracture toughness due to degradation or hydration. The models have been incorporated into peridynamics, non-local continuum mechanics methodology, that can model intersecting and branching brittle fracture that occurs in multicomponent brittle materials, such as concrete. Through development of new peridynamic capabilities, decalcification of cement and differential shrinkage in clay-cement composites have been evaluated, along with exemplar problems in nuclear waste cannisters and wellbores. We have developed methods to simulate multiphase phenomena in cement and cement-composite materials for energy and infrastructure applications.

Compliance with future ultra-low nitrogen oxide regulations with diesel engines requires the fastest possible heating of the exhaust aftertreatment system to its proper operating temperature upon cold starting. Late post injections are commonly integrated into catalyst-heating operating strategies. This experimental study provides insight into the complex interactions between the injection-strategy calibration and the tradeoffs between exhaust heat and pollutant emissions. Experiments are performed with certification diesel fuel and blends of diesel fuel with butylal and hexyl hexanoate. Further analyses of experimental data provide insight into fuel reactivity and oxygen content as potential enablers for improved catalyst-heating operation. A statistical design-of-experiments approach is developed to investigate a wide range of injection strategy calibrations at three different intake dilution levels. Thermodynamic and exhaust emissions measurements are taken using a new medium-duty, single-cylinder research engine. Analysis of the results provides insight into the effects of exhaust gas recirculation, oxygenated fuel blends, and fuel reactivity on exhaust heat and pollutant emissions. Late-cycle heat release is an important factor in determining exhaust temperatures. Intake dilution and fuel properties certainly affect late-cycle heat release, but the methods applied in this work are not sufficient to reproduce or explain the mechanisms by which improved fuel cetane rating promotes operation with hotter exhaust and lower pollutant emissions.

For many decades, lithium metal batteries (LMBs) have been studied for their high energy density, though they suffer from high reactivity, which safety concerns such as explosions and fires. These events can occur as a result of degradation due to redox reactions and lithium morphology that evolve during cycling. Mostly, the failure is caused by the electrolyte degradation leading to the creation of a solid electrolyte interphase (SEI), which tends to induce electron leakage, lithium dendrite growth, and capacity decay. SEI formation is practically inevitable. This work investigated the minimization of degradation events by modifying different battery configurations. Modifications were made by exchanging electrolytes, electrode surfaces, and even annexing innovative components. Our goal is to achieve a configuration where Li-ion transfer is facilitated, and electron transfer is limited which would enable a longer lithium metal battery life cycle. When performing these modifications and testing their electrochemical performance, a current reduction was observed when the electrodes were spontaneously coated with preformed SEI. Herein several classes of artificial SEI layers are explored for their ability to form passivating layers through use of a redox active molecule.

Stainless steel TPBAR components undergo neutron radiation-induced segregation and dislocation loop formation. Comparison experiments with ion beams accelerate the damage, and visualize the damage process with in-situ microscopy. In-situ Au irradiation causes defect formation, but no elemental segregation.

Technologies based on water adsorption such as water harvesting from air have tremendous potential in mitigating important global crises such as water scarcity. An important challenge to the deployment of such technologies is finding optimal adsorbent materials. Given the large materials space of available adsorbents, large-scale computational screening can be extremely helpful for this task. This work explores the methods and details associated with such screening procedures and recommends best practices. We also shed light on the limitations of traditionally used and inexpensive to compute prescreening approaches involving geometric and energetic features to predict water adsorption behavior of porous materials. Such approaches can provide general trends to predict adsorption behavior but may lead to the overlook of potentially important structures due to the complex nature of water adsorption. Finally, this study offers insights for future water adsorption simulations to facilitate the development of optimal water adsorbents.

This Analysis Report (AR) documents the determination of pH correction factors for the observed pH readings. The correction factor converts the observed pH reading recorded from the brines used in geochemical studies in support of the Waste Isolation Pilot Plant (WIPP) to a corrected pH value. The data analysis in this AR falls under AP-157 Rev. 1 Analysis Plan for Determination of pH Correction Factors in Brines (Kirkes et al, 2021). Measurement of pH in some solutions can be challenging due to numerous factors such as high ionic strength, elevated or lowered temperature, complex matrix composition, etc. (Knauss et al., 1990; and Rai et al., 1995). The measured pH can be corrected by applying the correction factor, empirically obtained from a specific test solution.

Ultradoping introduces unprecedented dopant levels into Si, which transforms its electronic behavior and enables its use as a next-generation electronic material. Commercialization of ultradoping is currently limited by gas-phase ultra-high vacuum requirements. Solvothermal chemistry is amenable to scale-up. However, an integral part of ultradoping is a direct chemical bond between dopants and Si, and solvothermal dopant-Si surface reactions are not well-developed. This work provides the first quantified demonstration of achieving ultradoping concentrations of boron (∼1e14 cm2) by using a solvothermal process. Surface characterizations indicate the catalyst cross-reacted, which led to multiple surface products and caused ambiguity in experimental confirmation of direct surface attachment. Density functional theory computations elucidate that the reaction results in direct B−Si surface bonds. This proof-of-principle work lays groundwork for emerging solvothermal ultradoping processes.

We provide further details for using Lim, et al.'s marked multidimensional Hawkes processes with dissimilar decays. We first describe what makes these different from other Hawkes processes, then describe each model hyperparameter and how to initialize it informed by the input data and any prior biases. We derived tighter bounds than Lim, et al. for faster convergence of rejection sampling. The resulting hyperparameters and bounds have been helpful against both synthetic and real-world datasets.

Complex organic-inorganic interfaces are important for device and sensing applications. Charge transfer doping is prevalent in such applications and can affect the interfacial energy level alignments (ELA), which are determined by many-body interactions. We develop an approximateab initiomany-body GW approach that can capture many-body interactions due to interfacial charge transfer. The approach uses significantly less resources than a regular GW calculation but gives excellent agreement with benchmark GW calculations on an F4TCNQ/graphene interface. We find that many-body interactions due to charge transfer screening result in gate-tunable F4TCNQ HOMO-LUMO gaps. We further predict the ELA of a large system of experimental interest—4,4′-bis(dimethylamino)bipyridine (DMAP-OED) on monolayer MoS2, where charge transfer screening results in an ∼1 eV reduction of the molecular HOMO-LUMO gap. Comparison with a two-dimensional electron gas model reveals the importance of explicitly considering the intraband transitions in determining the charge transfer screening in organic-inorganic interface systems.

This memo’s objective is to report a calibration of the J2 plasticity model with the Wilkins ductile failure criterion for 17-4 PH H1150 stainless steel under slow loading at room temperature. The calibration of the hardening function was based on uniaxial tension tests, while that of the failure model included data from tension tests on notched specimens, a butterfly specimen shear test, and a set of interrupted compression tests on shear hat specimens. The procedure was that described in, minus the rate and temperature dependence.

Abstract not provided.

The coupling of inter- and intramolecular vibrations plays a critical role in initiating chemistry during the shock-to-detonation transition in energetic materials. Herein, we report on the subpicosecond to subnanosecond vibrational energy transfer (VET) dynamics of the solid energetic material 1,3,5-trinitroperhydro-1,3,5-triazine (RDX) by using broadband, ultrafast infrared transient absorption spectroscopy. Experiments reveal VET occurring on three distinct time scales: subpicosecond, 5 ps, and 200 ps. The ultrafast appearance of signal at all probed modes in the mid-infrared suggests strong anharmonic coupling of all vibrations in the solid, whereas the long-lived evolution demonstrates that VET is incomplete, and thus thermal equilibrium is not attained, even on the 100 ps time scale. Density functional theory and classical molecular dynamics simulations provide valuable insights into the experimental observations, revealing compression-insensitive time scales for the initial VET dynamics of high-frequency vibrations and drastically extended relaxation times for low-frequency phonon modes under lattice compression. Mode selectivity of the longest dynamics suggests coupling of the N-N and axial NO2stretching modes with the long-lived, excited phonon bath.

Denoising contaminated seismic signals for later processing is a fundamental problem in seismic signals analysis. Neural network approaches have shown success denoising local signals when trained on short-time Fourier transform spectrograms. One challenge of this approach is the onerous process of hand-labeling event signals for training. By leveraging the SCALODEEP seismic event detector, we develop an automated set of techniques for labeling event data. Despite region specific challenges, training the neural network denoiser on machine curated events shows comparable performance to the neural network trained on hand curated events. We showcase our technique with two experiments, one using Utah regional data and one using regional data from the Korean peninsula.

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

Thermoelectric power plants often depend on multipurpose reservoirs to supply cooling water. Although reservoirs buffer natural hydrologic variability, severe droughts can deplete storage below critical thresholds, or to levels at which the effluent water temperature exceeds the environmental compliance requirement for cooling. This study explores the effects of projected climate change and drought on water storage at 30 major reservoirs in Texas. These reservoirs collectively provide cooling water for about two thirds of thermoelectric power capacity in the Electric Reliability Council of Texas (ERCOT) power grid. Multi-ensemble runoff projections generated from eleven downscaled hydroclimate simulations are mapped to key watersheds to create spatially correlated multi-reservoir inflow sequences. These data are used to drive reservoir storage simulations, which are linked to a metric of “capacity-at-risk” using critical reservoir thresholds. We find that projected impacts of climate change are mixed, with results indicating an increase in the occurrence of thermal disruption under only half of climate models. A critical threshold of 30% storage volume—applied to all reservoirs—results in disruption to about one fifth of ERCOT thermal generation during the most severe projected droughts. The study highlights an important role for detailed reservoir behavior simulations for capturing the effects of drought and climate change on thermoelectric plant performance.

We propose a domain decomposition method for the efficient simulation of nonlocal problems. Our approach is based on a multi-domain formulation of a nonlocal diffusion problem where the subdomains share “nonlocal” interfaces of the size of the nonlocal horizon. This system of nonlocal equations is first rewritten in terms of minimization of a nonlocal energy, then discretized with a meshfree approximation and finally solved via a Lagrange multiplier approach in a way that resembles the finite element tearing and interconnect method. Specifically, we propose a distributed projected gradient algorithm for the solution of the Lagrange multiplier system, whose unknowns determine the nonlocal interface conditions between subdomains. Several two-dimensional numerical tests on problems as large as 191 million unknowns illustrate the strong and the weak scalability of our algorithm, which outperforms the standard approach to the distributed numerical solution of the problem. Finally, this work is the first rigorous numerical study in a two-dimensional multi-domain setting for nonlocal operators with finite horizon and, as such, it is a fundamental step towards increasing the use of nonlocal models in large scale simulations.

Time-resolved spectroscopies using high-energy photons in the vacuum ultraviolet (VUV) to the X-ray region of the electromagnetic spectrum, have proven to be powerful probes of chemical dynamics. These high-energy photons can access valence and core orbitals of molecules and materials, providing key information on molecular and electronic structure and their time evolution. This report details the development of table-top sources of extreme ultraviolet (XUV) and VUV pulses at Sandia National Laboratories for use in studies of gas phase chemical dynamics. Femtosecond duration XUV pulses are produced using laser-driven high harmonic generation and their detected range span ~40-140 eV photon energies. These pulses are used in conjunction with ultraviolet pulses in a pump-probe scheme to study excited state dynamics of gas phase molecules. VUV pulses at 7.75 eV are generated using a four-wave-mixing scheme driven by 800 nm and 266 nm pulses in an argon-filled hollow-core fiber.

This is a simple model designed to run fast but still maintain the key physics and feedback mechanisms of a heat pipe. First, the capillary pressure is a function of the liquid working fluid volume fraction. Second, the boiling and condensation are based on the saturation temperature that is based on the heat pipe pressure. When the pressure goes up, the saturation temperature goes up and the vapor rains on the wick. When the pressure goes down, the saturation temperature goes down and the liquid in the entire wick boils. This is how the heat pipe adjusts to stay robust under different temperatures and heat fluxes.

Polymerization-induced phase separation enables fine control over thermoset network morphologies, yielding heterogeneous structures with domain sizes tunable over 1-100 nm. However, the controlled chain-growth polymerization techniques exclusively employed to regulate the morphology at these length scales are unsuitable for a majority of thermoset materials typically formed through step-growth mechanisms. By varying the composition of a binary curing agent mixture in a classic rubber-toughened epoxy thermoset, where the two curing agents are selected based on disparate compatibility with the rubber, we demonstrate facile tunability over morphology through a single compositional parameter. Indeed, this method yields morphologies spanning the nano-scale to the macro-scale, controlled by the relative reactivities and thermodynamic compatibility of the network components. We further demonstrate a profound connection between chain dynamics and microstructure in these materials, with the tunable morphology enabling exquisite variations in glass transition. In addition, previously unattainable control over tensile mechanical properties is realized, including atypical increase of elongation at failure while maintaining the modulus and ultimate strength.

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.

First-of-their kind datasets from a high-speed X-ray tomography system were collected, and a novel numerical effort utilizing temporal information to reduce measurement uncertainty was shown. The experimental campaign used three high-speed X-ray imaging systems to collect data at 100 kHz of a scene containing high-velocity objects. The scene was a group of known objects propelled by a 12-gauge shotgun shell reaching speeds of hundreds of meters per second. These data represent a known volume where the individual components are known, with experimental uncertainties that can be used for reconstruction algorithm validation. The numerical effort used synthetic volumes in MATLAB to produce projections along known lines of sight to perform tomographic reconstructions. These projections and reconstructions were performed on a single object at two orientations, representing two timesteps, to increase the reconstruction accuracy.

The How To Manual supplements the User’s Manual and the Theory Manual. The goal of the How To Manual is to reduce learning time for complex end to end analyses. These documents are intended to be used together. See the User’s Manual for a complete list of the options for a solution case. All the examples are part of the Sierra/SD test suite. Each runs as is. The organization is similar to the other documents: How to run, Commands, Solution cases, Materials, Elements, Boundary conditions, and then Contact. The table of contents and index are indispensable. The Geometric Rigid Body Modes section is shared with the Users Manual.

Photon sources able to emit single or entangled photon pairs are key components in quantum information systems. Semiconductor quantum dots (QDs) are promising candidates due to their high efficiencies and ease of integration with other photonic or electronic components. State-of-the-art QDs, however, are limited to certain emission wavelengths and specific applications due to material choice constraints and their randomness in shape/size. This project is focused on developing a novel in-situ patterning technique to realize QDs with a broad emission range, shape/size control and the ability to emit single/entangled photons. Our approach has two key elements: (1) In-situ patterning via arsenic-induced nanovoid etching on antimonide surfaces and (2) In-filling of nanovoids to form QDs. By closely controlling the experimental conditions, it is shown that this technique can be used to realize III-V QDs in As₂- etched nanovoids on a GaSb surface. The exposure to As₂ in terms of substrate temperature, time and flux is found to have a significant impact on the process variables such as nanovoid depth, QD dimensions etc. An in-depth analysis of the etch mechanism reveals that by controlling the As₂ exposure, uniform 3-dimensional nanostructures with varying areal densities can be obtained without an in-filling step. Preliminary optical characterization of these nanostructures shows that these QDs may be relevant for realizing emitters in the telecom wavelength range.

This document presents tests from the Sierra Structural Mechanics verification test suite. Each of these tests is run nightly with the Sierra/SD code suite and the results of the test checked versus the correct analytic result. For each of the tests presented in this document the test setup, derivation of the analytic solution, and comparison of the Sierra/SD code results to the analytic solution is provided. This document can be used to confirm that a given code capability is verified or referenced as a compilation of example problems.

Computational simulation is increasingly relied upon for high-consequence engineering decisions, and a foundational element to solid mechanics simulations is a credible material model. Our ultimate vision is to interlace material characterization and model calibration in a real-time feedback loop, where the current model calibration results will drive the experiment to load regimes that add the most useful information to reduce parameter uncertainty. The current work investigated one key step to this Interlaced Characterization and Calibration (ICC) paradigm, using a finite load-path tree to incorporate history/path dependency of nonlinear material models into a network of surrogate models that replace computationally-expensive finite-element analyses. Our reference simulation was an elastoplastic material point subject to biaxial deformation with a Hill anisotropic yield criterion. Training data was generated using either a space-filling or adaptive sampling method, and surrogates were built using either Gaussian process or polynomial chaos expansion methods. Surrogate error was evaluated to be on the order of 10^⁻5 and 10^⁻3 percent for the space-filling and adaptive sampling training data, respectively. Direct Bayesian inference was performed with the surrogate network and with the reference material point simulator, and results agreed to within 3 significant figures for the mean parameter values, with a reduction in computational cost over 5 orders of magnitude. These results bought down risk regarding the surrogate network and facilitated a successful FY22-24 full LDRD proposal to research and develop the complete ICC paradigm.

The Spent Fuel and Waste Science and Technology (SFWST) Campaign of the U.S. Department of Energy (DOE) Office of Nuclear Energy (NE), Office of Spent Fuel & Waste Disposition (SFWD) is conducting research and development (R&D) on geologic disposal of spent nuclear fuel (SNF) and highlevel nuclear waste (HLW). A high priority for SFWST disposal R&D is disposal system modeling (DOE 2012, Table 6; Sevougian et al. 2019). The SFWST Geologic Disposal Safety Assessment (GDSA) work package is charged with developing a disposal system modeling and analysis capability for evaluating generic disposal system performance for nuclear waste in geologic media.

Designing next generation thin films, tailor-made for specific applications, relies on the availability of robust processing-structure-property relationships. Traditional structure zone diagrams are limited to low-dimensional mappings, with machine-learning methods only recently attempting to relate multiple processing parameters to the final microstructure. Despite this progress, structure-processing relationships are unknown for processing conditions that vary during thin-film deposition, limiting the range of microstructures and properties achievable. In this project, we employed a phase-field computational model combined with a genetic algorithm (GA) to identify and design time-dependent processing protocols that achieve tailor-made microstructures. We simulate the physical vapor deposition of a binary-alloy thin film by employing a phase-field model, where deposition rates and diffusivities are controlled via the genetic algorithm. Our GA-guided protocols achieve targeted microstructures with lateral and vertical concentration modulations, as well as more complex, hierarchical microstructures previously not described in simple structure zone diagrams. Our algorithm provides insight to experimentalists looking for additional avenues to design novel thin-film microstructures.

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

To date, terahertz quantum-cascade vertical-external-cavity surface-emitting lasers (QC-VECSELs) have tended to oscillate in only one or two lasing modes at a time. This is due to the fact that the interaction between all of the longitudinal external cavity modes and the QC gain material is mediated through a single metasurface resonance, whose spatial overlap changes little with frequency; this suppresses spatial-hole-burning induced multi-mode operation. In this Letter, a VECSEL external cavity is demonstrated using an output coupler based upon a high-resistivity silicon etalon, which presents a periodic reflectance spectrum that is nearly matched with the external cavity mode spectrum. As the cavity length is varied, a systematic transition between a single/double-mode lasing regime and a multi-mode lasing regime is realized due to the Vernier effect. Up to nine modes lasing simultaneously with a free-spectral-range of approximately 21 GHz is demonstrated. This result provides a path toward the multi-mode operation necessary for eventual frequency comb operation.

Thin, high-density layers of dopants in semiconductors, known as δ-layer systems, have recently attracted attention as a platform for exploration of the future quantum and classical computing when patterned in plane with atomic precision. However, there are many aspects of the conductive properties of these systems that are still unknown. Here we present an open-system quantum transport treatment to investigate the local density of electron states and the conductive properties of the δ-layer systems. A successful application of this treatment to phosphorous δ-layer in silicon both explains the origin of recently-observed shallow sub-bands and reproduces the sheet resistance values measured by different experimental groups. Further analysis reveals two main quantum-mechanical effects: 1) the existence of spatially distinct layers of free electrons with different average energies; 2) significant dependence of sheet resistance on the δ-layer thickness for a fixed sheet charge density.

Sierra/SolidMechanics (Sierra/SM) is a Lagrangian, three-dimensional code for finite element analysis of solids and structures. It provides capabilities for explicit dynamic, implicit quasistatic and dynamic analyses. The explicit dynamics capabilities allow for the efficient and robust solution of models with extensive contact subjected to large, suddenly applied loads. For implicit problems, Sierra/SM uses a multi-level iterative solver, which enables it to effectively solve problems with large deformations, nonlinear material behavior, and contact. Sierra/SM has a versatile library of continuum and structural elements, and a large library of material models. The code is written for parallel computing environments enabling scalable solutions of extremely large problems for both implicit and explicit analyses. It is built on the SIERRA Framework, which facilitates coupling with other SIERRA mechanics codes. This document describes the functionality and input syntax for Sierra/SM.

Relativistic weakly collisional plasmas describe a variety of astrophysical and terrestrial plasmas, ranging from relativistic outflows from active galactic nuclei to high power microwave and magnetically insulated transmission lines. In many such systems, high fidelity kinetic models are computationally infeasible due large dynamical scales and long dynamical times. Conversely, most fluid based models such as magnetohydrodynamics (MHD) miss many relevant aspects of plasma behavior. Between these two models, two fluid methods - where the electrons and ions are evolved as separate, coupled fluids - capture many of the plasma physics of a kinetic code while remaining computational tractable for large systems.

We present a novel technique for numerically modeling relativistic magnetrons. The electrons are represented with a 5-moment relativistic fluid. Typically, the particle in cell method is used for simulated relativistic high-power microwave sources. This study considers the A6 magnetron presented by Palevsky and Bekefi [1].

C22-Gd films were deposited onto 316 stainless steel up to 2.4mm thick using three different cold spray conditions. No differences in visible coating adhesion was observed between grit-blasted and as-received 316 stainless steel substrates. All of the coatings produced had a density greater than 95% of bulk. EDS and backscattered imaging suggest that no secondary phases precipitated during coating fabrication, but potential micro-segregation of phases were likely present in the as-received material. The phase segregation was still observable in the as-deposited coatings. EDS could not resolve the differences in composition of the two phases.

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

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

We present a new geodesic-based method for geometry optimization in a basis set of redundant internal coordinates. Our method updates the molecular geometry by following the geodesic generated by a displacement vector on the internal coordinate manifold, which dramatically reduces the number of steps required to converge to a minimum. Our method can be implemented in any existing optimization code, requiring only implementation of derivatives of the Wilson B-matrix and the ability to numerically solve an ordinary differential equation.

The passivation of polycrystalline nickel surfaces against hydrogen uptake by oxygen is investigated experimentally with low energy ion scattering (LEIS), direct recoil spectroscopy (DRS), and thermal desorption spectroscopy (TDS). These techniques are highly sensitive to surface hydrogen, allowing the change in hydrogen adsorption in response to varying amounts of oxygen exposure to be measured. The chemical composition of a nickel surface during a mixed oxygen and hydrogen exposure was characterized with LEIS and DRS, while the uptake and activation energies of hydrogen on a nickel surface with preadsorbed oxygen were quantified with TDS. By and large, these measurements of how the oxygen and hydrogen surface coverage varied in response to oxygen exposure were found to be consistent with predictions of a simple site-blocking model. This finding suggests that, despite the complexities that arise due to polycrystallinity, the oxygen-induced passivation of a polycrystalline nickel surface against hydrogen uptake can be approximated by a simple site-blocking model.

Presented in this document are the theoretical aspects of capabilities contained in the Sierra/SM code. This manuscript serves as an ideal starting point for understanding the theoretical foundations of the code. For a comprehensive study of these capabilities, the reader is encouraged to explore the many references to scientific articles and textbooks contained in this manual. It is important to point out that some capabilities are still in development and may not be presented in this document. Further updates to this manuscript will be made as these capabilities come closer to production level.

A coded aperture is used to demonstrate emission spectroscopy from multiple one-dimensional measurement locations simultaneously with a single camera. The coded aperture mask has several columns of periodic apertures, each with a unique spatial frequency. Light transmitted through all mask columns is detected through an imaging spectrometer. Dispersed light from the various mask columns overlaps on the spectrometer camera but is separated using Fourier-domain filtering using the known spatial frequencies of the mask. As the coded aperture is placed at an image plane, each Fourier-filtered spectrogram comes from a unique one-dimensional measurement location. This technique represents a significant increase in the amount of spatially and spectrally resolved emission data available using a single emission spectrometer and camera at the expense of some spatial resolution due to the Fourier filtering. This instrument is particularly useful for studying transient, non-repeating events. Megahertz-rate emission spectroscopy from five one-dimensional measurement locations is demonstrated with explosive fireballs using a single camera. Optical design parameters and instrument performance characteristics are discussed.

Electrochemical models at different scales and varying levels of complexity have been used in the literature to study the evolution of the anode surface in lithium metal batteries. This includes continuum, mesoscale (phase-field approaches), and multiscale models. In this paper, using a motivating example of a moving boundary model in one dimension, we show how battery models need proper formulation for mass conservation, especially when simulated over multiple charge and discharge cycles. The article concludes with some thoughts on mass conservation and proper formulation for multiscale models.

This report summarizes the current actives in FY21 related to the effort by Sandia National Laboratories to identify and test coating materials for the prevention, mitigation, and repair of spent nuclear fuel dry storage canisters against potential chloride-induced stress corrosion cracking. This work follows up on the details provided in Sandia National Laboratories FY20 report on the same topic, which provided a detailed description of the specific coating properties desired for application and implementation on spent nuclear fuel canisters, as well as provided detail into several different coatings and their applicability to coat spent nuclear fuel canisters. In FY21, Sandia National Laboratories has engaged with private industry to create a Memorandum of Understanding and established a collaborative R&D program building off the analytical and laboratory capabilities at Sandia National Laboratories and the material design and synthesis capabilities of private industry. The resulting Memorandum of Understanding included four companies to date (Oxford Performance Materials, White Horse R&D, Luna Innovations, and Flora Coating) proposing six different coating technologies (polyetherketoneketone, modified polyimide/polyurea, modified phenolic resin, silane-based polyurethane hybrid with and without a Znrich primer, and a quasi-ceramic sol-gel polyurethane hybrid) to be tested, evaluated, and optimized for their potential use for this application. This report provides a detailed description of each of the coating systems proposed by the participating industry partners. It also provides a description of the planned experimental actives to be performed by Sandia National Laboratories including physical tests, electrochemical tests, and characterization methods. These analyses will be used to identify specific ways to further improve coating technologies toward their application and implementation on spent nuclear fuel canisters. In FY21, Sandia National Laboratories began baseline testing of the base metal material in according with activities of the Memorandum of Understanding. In FY22, Sandia National Laboratories will receive coated coupons from each of the participating industry partners and begin characterization, physical, and electrochemical testing following the test plan described herein.

The finite element method is a scheme to discretize the infinite number of degrees of freedom in continuum-level problems down to a finite number of degrees of freedom. This discretization is done in conjunction with methods that also reduce the field differential equations to sets of algebraic ones that can be solved by arithmetical operations. Therefore, solutions attained by finite element models are approximations to the exact solutions of the field equations.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The propagation of a wave pulse due to low-speed impact on a one-dimensional, heterogeneous bar is studied. Due to the dispersive character of the medium, the pulse attenuates as it propagates. This attenuation is studied over propagation distances that are much longer than the size of the microstructure. A homogenized peridynamic material model can be calibrated to reproduce the attenuation and spreading of the wave. The calibration consists of matching the dispersion curve for the heterogeneous material near the limit of long wavelengths. It is demonstrated that the peridynamic method reproduces the attenuation of wave pulses predicted by an exact microstructural model over large propagation distances.

Abstract not provided.

Achieving efficient learning for AI systems was identified as a major challenge in the DOE's recently released, AI for Science, report. The human brain is capable of efficient and low-powered learning. It is likely that implementing brain-like principles will lead to more efficient AI systems. In this LDRD, I aim to contribute to this goal by creating a foundation for implementing and studying a brain phenomenon termed short term plasticity (STP) in spiking artificial neural networks within Sandia. First, data collected by the Allen Institute for Brain Science (AIBS) was analyzed to see if STP could be classified into types using the data collected. Although the data was inadequate at the time, AIBS has updated their database and created models that could be utilized in the future. Second, I began creating a software package to assess the ability of a Boltzmann machine utilizing STP to sample from national security data.

This report provides a summary of notes for building and running the Sandia Computational Engine for Particle Transport for Radiation Effects (SCEPTRE) code. SCEPTRE is a general- purpose C++ code for solving the linear Boltzmann transport equation in serial or parallel using unstructured spatial finite elements, multigroup energy treatment, and a variety of angular treatments including discrete ordinates (Sn) and spherical harmonics (Pn). Either the first-order form of the Boltzmann equation or one of the second-order forms may be solved. SCEPTRE requires a small number of open-source Third Party Libraries (TPL) to be available, and example scripts for building these TPL are provided. The TPL needed by SCEPTRE are Trilinos, Boost, and Netcdf. SCEPTRE uses an autotools build system, and a sample configure script is provided. Running the SCEPTRE code requires that the user provide a spatial finite-elements mesh in Exodus format and a cross section library in a format that will be described. SCEPTRE uses an xml-based input, and several examples will be provided.

Abstract not provided.

Ongoing progress in synthetic biology, metabolic engineering, and catalysis continues to produce a diverse array of advanced biofuels with complex molecular structure and functional groups. In order to integrate biofuels into existing combustion systems, and to optimize the design of next-generation combustion systems, understanding connections between molecular structure and ignition at low-temperature conditions (< 1000 K) remains a priority that is addressed in part using chemical kinetics modeling. The development of predictive models relies on detailed information, derived from experimental and theoretical studies, on molecular structure and chemical reactivity, both of which influence the balance of chain reactions that occur during combustion – propagation, termination, and branching. In broad context, three main categories of reactions affect ignition behavior: (i) initiation reactions that generate a distribution of organic radicals, R˙; (ii) competing unimolecular decomposition of R˙ and bimolecular reaction of R˙ with O2; (iii) decomposition mechanisms of peroxy radical adducts (ROO˙), including isomerization via ROO˙ ⇌ Q˙OOH. All three categories are influenced by functional groups in different ways, which causes a shift in the balance of chain reactions that unfold over complex temperature- and pressure-dependent mechanisms. The objective of the present review is three-fold: (1) to provide a historical account of research on low-temperature oxidation of biofuels, including initiation reactions, peroxy radical reactions, Q˙OOH-mediated reaction mechanisms, and chain-branching chemistry; (2) to summarize the influence of functional groups on chemical kinetics relevant to chain-branching reactions, which are responsible for the accelerated production of radicals that leads to ignition; (3) to identify areas of research that are needed – experimentally and computationally – to address fundamental questions that remain. Results from experimental, quantum chemical, and chemical kinetics modeling studies are reviewed for several classes of biofuels – alcohols, esters, ketones, acyclic ethers and cyclic ethers – and are compared against analogous results in alkane oxidation. The review is organized into separate sections for each biofuel class, which include studies on thermochemistry and bond dissociation energies, rate coefficients for initiation reactions via H-abstraction and related branching fractions, reaction mechanisms and product formation from reactive intermediates, ignition delay times, and chemical kinetics modeling. Each section is then summarized in order to identify areas for which additional functional group-specific work is required. The review concludes with an outline for research directions for improving the fundamental understanding of biofuel ignition chemistry and related chemical kinetics modeling.

Abstract not provided.

As the seismic monitoring community advances toward detecting, identifying, and locating ever-smaller natural and anthropogenic events, the need is constantly increasing for higher resolution, higher fidelity data, models, and methods for accurately characterizing events. Local-distance seismic data provide robust constraints on event locations, but also introduce complexity due to the significant geologic heterogeneity of the Earth’s crust and upper mantle, and the relative sparsity of data that often occurs with small events recorded on regional seismic networks. Identifying the critical characteristics for improving local-scale event locations and the factors that impact location accuracy and reliability is an ongoing challenge for the seismic community. Using Utah as a test case, we examine three data sets of varying duration, finesse, and magnitude to investigate the effects of local earth structure and modeling parameters on local-distance event location precision and accuracy. We observe that the most critical elements controlling relocation precision are azimuthal coverage and local-scale velocity structure, with tradeoffs based on event depth, type, location, and range.

In this paper we analyze the noise in macro-particle methods used in plasma physics and fluid dynamics, leading to approaches for minimizing the total error, focusing on electrostatic models in one dimension. We begin by describing kernel density estimation for continuous values of the spatial variable x, expressing the kernel in a form in which its shape and width are represented separately. The covariance matrix C(x,y) of the noise in the density is computed, first for uniform true density. The bandwidth of the covariance matrix is related to the width of the kernel. A feature that stands out is the presence of constant negative terms in the elements of the covariance matrix both on and off-diagonal. These negative correlations are related to the fact that the total number of particles is fixed at each time step; they also lead to the property ∫C(x,y)dy=0. We investigate the effect of these negative correlations on the electric field computed by Gauss's law, finding that the noise in the electric field is related to a process called the Ornstein-Uhlenbeck bridge, leading to a covariance matrix of the electric field with variance significantly reduced relative to that of a Brownian process. For non-constant density, ρ(x), still with continuous x, we analyze the total error in the density estimation and discuss it in terms of bias-variance optimization (BVO). For some characteristic length l, determined by the density and its second derivative, and kernel width h, having too few particles within h leads to too much variance; for h that is large relative to l, there is too much smoothing of the density. The optimum between these two limits is found by BVO. For kernels of the same width, it is shown that this optimum (minimum) is weakly sensitive to the kernel shape. We repeat the analysis for x discretized on a grid. In this case the charge deposition rule is determined by a particle shape. An important property to be respected in the discrete system is the exact preservation of total charge on the grid; this property is necessary to ensure that the electric field is equal at both ends, consistent with periodic boundary conditions. We find that if the particle shapes satisfy a partition of unity property, the particle charge deposited on the grid is conserved exactly. Further, if the particle shape is expressed as the convolution of a kernel with another kernel that satisfies the partition of unity, then the particle shape obeys the partition of unity. This property holds for kernels of arbitrary width, including widths that are not integer multiples of the grid spacing. We show results relaxing the approximations used to do BVO optimization analytically, by doing numerical computations of the total error as a function of the kernel width, on a grid in x. The comparison between numerical and analytical results shows good agreement over a range of particle shapes. We discuss the practical implications of our results, including the criteria for design and implementation of computationally efficient particle shapes that take advantage of the developed theory.

Abstract not provided.

Abstract not provided.

With the rapid proliferation of additive manufacturing and 3D printing technologies, architected cellular solids including truss-like 3D lattice topologies offer the opportunity to program the effective material response through topological design at the mesoscale. The present report summarizes several of the key findings from a 3-year Laboratory Directed Research and Development Program. The program set out to explore novel lattice topologies that can be designed to control, redirect, or dissipate energy from one or multiple insult environments relevant to Sandia missions, including crush, shock/impact, vibration, thermal, etc. In the first 4 sections, we document four novel lattice topologies stemming from this study: coulombic lattices, multi-morphology lattices, interpenetrating lattices, and pore-modified gyroid cellular solids, each with unique properties that had not been achieved by existing cellular/lattice metamaterials. The fifth section explores how unintentional lattice imperfections stemming from the manufacturing process, primarily sur face roughness in the case of laser powder bed fusion, serve to cause stochastic response but that in some cases such as elastic response the stochastic behavior is homogenized through the adoption of lattices. In the sixth section we explore a novel neural network screening process that allows such stocastic variability to be predicted. In the last three sections, we explore considerations of computational design of lattices. Specifically, in section 7 using a novel generative optimization scheme to design novel pareto-optimal lattices for multi-objective environments. In section 8, we use computational design to optimize a metallic lattice structure to absorb impact energy for a 1000 ft/s impact. And in section 9, we develop a modified micromorphic continuum model to solve wave propagation problems in lattices efficiently.

Although many software teams across the laboratories comply with yearly software quality engineering (SQE) assessments, the practice of introducing quality into each phase of the software lifecycle, or the team processes, may vary substantially. Even with the support of a quality engineer, many teams struggle to adapt and right-size software engineering best practices in quality to fit their context, and these activities aren’t framed in a way that motivates teams to take action. In short, software quality is often a “check the box for compliance” activity instead of a cultural practice that both values software quality and knows how to achieve it. In this report, we present the results of our 6600 VISTA Innovation Tournament project, "Incentivizing and Motivating High Confidence and Research Software Teams to Adopt the Practice of Quality." We present our findings and roadmap for future work based on 1) a rapid review of relevant literature, 2) lessons learned from an internal design thinking workshop, and 3) an external Collegeville 2021 workshop. These activities provided an opportunity for team ideation and community engagement/feedback. Based on our findings, we believe a coordinated effort (e.g. strategic communication campaign) aimed at diffusing the innovation of the practice of quality across Sandia National Laboratories could over time effect meaningful organizational change. As such, our roadmap addresses strategies for motivating and incentivizing individuals ranging from early career to seasoned software developers/scientists.

Abstract not provided.

The properties of materials can change dramatically at the nanoscale new and useful properties can emerge. An example is found in the paramagnetism in iron oxide magnetic nanoparticles. Using magnetically sensitive nitrogen-vacancy centers in diamond, we developed a platform to study electron spin resonance of nanoscale materials. To implement the platform, diamond substrates were prepared with nitrogen vacancy centers near the surface. Nanoparticles were placed on the surface using a drop casting technique. Using optical and microwave pulsing techniques, we demonstrated T1 relaxometry and double electron-electron resonance techniques for measuring the local electron spin resonance. The diamond NV platform developed in this project provides a combination of good magnetic field sensitivity and high spatial resolution and will be used for future investigations of nanomaterials and quantum materials.

In this project, our goal was to develop methods that would allow us to make accurate predictions about individual differences in human cognition. Understanding such differences is important for maximizing human and human-system performance. There is a large body of research on individual differences in the academic literature. Unfortunately, it is often difficult to connect this literature to applied problems, where we must predict how specific people will perform or process information. In an effort to bridge this gap, we set out to answer the question: can we train a model to make predictions about which people understand which languages? We chose language processing as our domain of interest because of the well- characterized differences in neural processing that occur when people are presented with linguistic stimuli that they do or do not understand. Although our original plan to conduct several electroencephalography (EEG) studies was disrupted by the COVID-19 pandemic, we were able to collect data from one EEG study and a series of behavioral experiments in which data were collected online. The results of this project indicate that machine learning tools can make reasonably accurate predictions about an individual?s proficiency in different languages, using EEG data or behavioral data alone.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

There is a dearth in the literature on how to capture the uncertainty generated by material surface evolution in thermal modeling. This leads to inadequate or highly variable uncertainty representations for material properties, specifically emissivity when minimal information is available. Inaccurate understandings of prediction uncertainties may lead decision makers to incorrect conclusions, so best engineering practices should be developed for this domain. In order to mitigate the aforementioned issues, this study explores different strategies to better capture the thermal uncertainty response of engineered systems exposed to fire environments via defensible emissivity uncertainty characterizations that can be easily adapted to a variety of use cases. Two unique formulations (one physics-informed and one mathematically based) are presented. The formulations and methodologies presented herein are not exhaustive but more so are a starting point and give the reader a basis for how to customize their uncertainty definitions for differing fire scenarios and materials. Finally, the impact of using this approach versus other commonly used strategies and the usefulness of adding rigor to material surface evolution uncertainty is demonstrated.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Whether calibrating equipment or inspecting products on the factory floor, metrology requires many complicated statistical calculations to achieve a full understanding and evaluation of measurement uncertainty and quality. In order to assist its workforce in performing these calculations in a consistent and rigorous way, the Primary Standards Lab at Sandia National Laboratories (SNL) has developed a free and open-source software package for computing various metrology calculations from uncertainty propagation to risk analysis. In addition to propagating uncertainty through a measurement model using the well-known Guide to Expression of Uncertainty in Measurement or Monte Carlo approaches, evaluating the individual Type A and Type B uncertainty components that go into the measurement model often requires other statistical methods such as analysis of variance or determining uncertainty in a fitted curve. Once the uncertainty in a measurement has been calculated, it is usually evaluated from a risk perspective to ensure the measurement is suitable for making a particular conformance decision. Finally, SNL’s software can perform all these calculations in a single application via an easy-to-use graphical interface, where the different functions are integrated so the results of one calculation can be used as inputs to another calculation.

Identification and characterization of underground events from surface or remote data requires a thorough understanding of the rock material properties. However, material properties usually come from borehole data, which is expensive and not always available. A potential alternative is to use topographic characteristics to approximate the strength, but this has never been done before quantitatively. Here we present the results from the first steps towards this goal. We have found that there are strong correlations between compressive and tensile strengths and slopes, but these correlations vary depending on data analysis details. Rugosity may be better correlated to strength than slope values. More comprehensive analyses are needed to fully understand the best method of predicting strength from topography for this area. We also found that misalignment of multiple GIS datasets can have a large influence on the ability to make interpretations. Lastly, these results will require further study in a variety of climatic conditions before being applicable to other sites.

This report summarizes GADRAS methods and gamma spectrometry rodeo uranium isotopic analysis results for high energy resolution H3D M400 cadmium zinc telluride (CZT) and ORTEC Micro Detective high-purity germanium (HPGe) spectra of uranium isotopic standards collected at Oak Ridge National Laboratory (ORNL) and Lawrence Livermore National Laboratory (LLNL) over a two-year measurement campaign. During the campaign, measurements were performed with the detectors unshielded, side shielded, and collimated. In addition, measurements of the uranium isotopic standards were performed unshielded and shielded.

Abstract not provided.

A simple approach to simulate contact between deformable objects is presented which relies on levelset descriptions of the Lagrangian geometry and an optimization-based solver. Modeling contact between objects remains a significant challenge for computational mechanics simulations. Common approaches are either plagued by lack of robustness or are exceedingly complex and require a significant number of heuristics. In contrast, the levelset contact approach presented herein is essentially heuristic free. Furthermore, the presented algorithm enables resolving and enforcing contact between objects with a significant amount of initial overlap. Examples demonstrating the feasibility of this approach are shown, including the standard Hertz contact problem, the robust removal of overlap between two overlapping blocks, and overlap-removal and pre-load for a bolted configuration.

Arithmetic Coding (AC) using Prediction by Partial Matching (PPM) is a compression algorithm that can be used as a machine learning algorithm. This paper describes a new algorithm, NGram PPM. NGram PPM has all the predictive power of AC/PPM, but at a fraction of the computational cost. Unlike compression-based analytics, it is also amenable to a vector space interpretation, which creates the ability for integration with other traditional machine learning algorithms. AC/PPM is reviewed, including its application to machine learning. Then NGram PPM is described and test results are presented, comparing them to AC/PPM.

Abstract not provided.

Physical protection of public buildings has long been a concern of police and security services where a balance of facility security and personnel safety is vital. Due to the nature of public spaces, the use of permanently installed and deploy-on-demand physical barrier systems must be safe for the legitimate occupants and visitors of that space. Such systems must seek to mitigate the personal and organizational consequences of unintentionally seriously injuring or killing an innocent bystander by slamming a heavy, rigid, and quick-deploying barrier into place. Consideration and implementation of less-than-lethal technologies is necessary to reduce risk to visitors and building personnel. One potential barrier solution is a fast-acting, high-strength, composite airbag barrier system for doorways and hallways to quickly deploy a less-than-lethal barrier at entry points as well as isolate intruders who have already gained access. This system is envisioned to be stored within an architecturally attractive selectively frangible shell that could be permanently installed at a facility or installed in remote or temporary locations as dictated by risk. The system would be designed to be activated remotely (hardwired or wireless) from a Central Alarm Station (CAS) or other secure location.

State-of-the-art engineering and science codes have grown in complexity dramatically over the last two decades. Application teams have adopted more sophisticated development strategies, leveraging third party libraries, deploying comprehensive testing, and using advanced debugging and profiling tools. In today's environment of diverse hardware platforms, these applications also desire performance portability-avoiding the need to duplicate work for various platforms. The Kokkos EcoSystem provides that portable software stack. Based on the Kokkos Core Programming Model, the EcoSystem provides math libraries, interoperability capabilities with Python and Fortran, and Tools for analyzing, debugging, and optimizing applications. In this article, we overview the components, discuss some specific use cases, and highlight how codesigning these components enables a more developer friendly experience.

The twenty-seven critical experiments in this series were performed in 2020 in the SCX at the Sandia Pulsed Reactor Facility. The experiments are grouped by fuel rod pitch. Case 1 is a base case with a pitch of 0.8001 cm and no water holes in the array. Cases 2 through 6 have the same pitch as Case 1 but contain various configurations with water holes, providing slight variations in the fuel-to-water ratio. Similarly, Case 7 is a base case with a pitch of 0.854964 cm and no water holes in the array. Cases 8 through 11 have the same pitch as Case 7 but contain various configurations with water holes. Cases 12 through 15 have a pitch of 1.131512 cm and differ according to the number of water holes in the array, with Case 12 having no water holes. Cases 16 through 19 have a pitch of 1.209102 cm and differ according to number of water holes in the array, with Case 16 having no water holes. Cases 20 through 23 have a pitch of 1.6002 cm and differ according to number of water holes in the array, with Case 20 having no water holes. Cases 24 through 27 have a pitch of 1.709928 cm and differ according to number of water holes in the array, with Case 24 having no water holes. As the experiment case number increases, the fuel-to-water volume ratio decreases.

Abstract not provided.

Typical approaches to classify scenes from light convert the light field to electrons to perform the computation in the digital electronic domain. This conversion and downstream computational analysis require significant power and time. Diffractive neural networks have recently emerged as unique systems to classify optical fields at lower energy and high speeds. Previous work has shown that a single layer of diffractive metamaterial can achieve high performance on classification tasks. In analogy with electronic neural networks, it is anticipated that multilayer diffractive systems would provide better performance, but the fundamental reasons for the potential improvement have not been established. In this work, we present extensive computational simulations of two - layer diffractive neural networks and show that they can achieve high performance with fewer diffractive features than single layer systems.

Metal oxides have been an attractive option for a range of applications, including hydrogen sensors, microelectronics, and catalysis, due to their reactivity and tunability. The properties of metal oxides can vary greatly on their precise surface structure; however, few surface science techniques can achieve atomistic-level determinations of surface structure, and fewer yet can do so for insulator surfaces. Low energy ion beam analysis offers a potential insulator-compatible solution to characterizing the surface structure of metal oxides. As a feasibility study, we apply low energy ion beam analysis to investigate the surface structure of a magnetite single crystal, Fe₃O₄(100). We obtain multi-angle maps using both forward-scattering low energy ion scattering (LEIS) and backscattering impact-collision ion scattering spectroscopy (ICISS). Both sets of experimental maps have intensity patterns that reflect the symmetries of the Fe₃O₄(100) surface structure. However, analytical interpretation of these intensity patterns to extract details of the surface structure is significantly more complex than previous LEIS and ICISS structural studies of one-component metal crystals, which had far more symmetries to exploit. To gain further insight into the surface structure, we model our experimental measurements with ion-trajectory tracing simulations using molecular dynamics. Our simulations provide a qualitative indication that our experimental measurements agree better with a subsurface cation vacancy model than a distorted bulk model.

Abstract not provided.

This project developed prototype germanium telluride switches, which can be used in RF applications to improve SWAP (size, weight, and power) and signal quality in RF systems. These switches can allow for highly reconfigurable systems, including antennas, communications, optical systems, phased arrays, and synthetic aperture radar, which all have high impact on current National Security goals for improved communication systems and communication technology supremacy. The final result of the project was the demonstration of germanium telluride RF switches, which could act as critical elements necessary for a single chip RF communication system that will demonstrate low SWAP and high reconfigurability

Abstract not provided.

Abstract not provided.

Ultra-low voltage drop tunnel junctions (TJs) were utilized to enable multi-active region blue light emitting diodes (LEDs) with up to three active regions in a single device. The multi-active region blue LEDs were grown monolithically by metal-organic chemical vapor deposition (MOCVD) without growth interruption. This is the first demonstration of a MOCVD grown triple-junction LED. Optimized TJ design enabled near-ideal voltage and EQE scaling close to the number of junctions. This work demonstrates that with proper TJ design, improvements in wall-plug efficiency at high output power operation are possible by cascading multiple III-nitride based LEDs.

Currently a set of 71 radionuclides are accounted for in off-site consequence analysis for LWRs. Radionuclides of dose consequence are expected to change for non-LWRs, with radionuclides of interest being type-specific. This document identifies an expanded set of radionuclides that may need to be accounted for in multiple non-LWR systems: high temperature gas reactors (HTGRs); fluoride-salt-cooled high-temperature reactors (FHRs); thermal-spectrum fluoride-based molten salt reactors (MSRs); fast-spectrum chloride-based MSRs; and, liquid metal fast reactors with metallic fuel (LMRs) Specific considerations are provided for each reactor type in Chapter 2 through Chapter 5, and a summary of all recommendations is provided in Chapter 6. All identified radionuclides are already incorporated within the MACCS software, yet the development of tritium-specific and carbon-specific chemistry models are recommended.

Abstract not provided.

Recent developments in the renewable energy sector have seen an unprecedented growth in residential photovoltaic (PV) installations. However, high PV penetration levels often lead to overvoltage problems in low-voltage (LV) distribution feeders. Smart inverter control such as active power curtailment (APC)-based overvoltage control can be implemented to overcome these challenges. The APC technique utilizes a constant droop-based approach which curtails power rigidly, which can lead to significant energy curtailment in the LV distribution feeders. In this paper, different variations of the APC technique with linear, quadratic, and exponential droops have been analyzed from the point-of-view of energy curtailment for a LV distribution network in North America. Further, a combinatorial approach using various droop-based APC methods in conjunction with adaptive dynamic programming (ADP) as a supplementary control scheme has also been proposed. The proposed approach minimizes energy curtailment in the LV distribution network by adjusting the droop gains. Simulation results depict that ADP in conjunction with exponential droop reduces the energy curtailment to approximately 50% compared to using the standard linear droop.

This report summarizes the activities performed as part of the Science and Engineering of Cybersecurity by Uncertainty quantification and Rigorous Experimentation (SECURE) Grand Challenge LDRD project. We provide an overview of the research done in this project, including work on cyber emulation, uncertainty quantification, and optimization. We present examples of integrated analyses performed on two case studies: a network scanning/detection study and a malware command and control study. We highlight the importance of experimental workflows and list references of papers and presentations developed under this project. We outline lessons learned and suggestions for future work.

Leakage pathways through caprock lithologies for underground storage of CO2 and/or enhanced oil recovery (EOR) include intrusion into nano-pore mudstones, flow within fractures and faults, and larger-scale sedimentary heterogeneity (e.g., stacked channel deposits). To assess multiscale sealing integrity of the caprock system that overlies the Morrow B sandstone reservoir, Farnsworth Unit (FWU), Texas, USA, we combine pore-to-core observations, laboratory testing, well logging results, and noble gas analysis. A cluster analysis combining gamma ray, compressional slowness, and other logs was combined with caliper responses and triaxial rock mechanics testing to define eleven lithologic classes across the upper Morrow shale and Thirteen Finger limestone caprock units, with estimations of dynamic elastic moduli and fracture breakdown pressures (minimum horizontal stress gradients) for each class. Mercury porosimetry determinations of CO2 column heights in sealing formations yield values exceeding reservoir height. Noble gas profiles provide a “geologic time-integrated” assessment of fluid flow across the reservoir-caprock system, with Morrow B reservoir measurements consistent with decades-long EOR water-flooding, and upper Morrow shale and lower Thirteen Finger limestone values being consistent with long-term geohydrologic isolation. Together, these data suggest an excellent sealing capacity for the FWU and provide limits for injection pressure increases accompanying carbon storage activities.

We present a simple and powerful technique for finding a good error model for a quantum processor. The technique iteratively tests a nested sequence of models against data obtained from the processor, and keeps track of the best-fit model and its wildcard error (a metric of the amount of unmodeled error) at each step. Each best-fit model, along with a quantification of its unmodeled error, constitutes a characterization of the processor. We explain how quantum processor models can be compared with experimental data and to each other. We demonstrate the technique by using it to characterize a simulated noisy two-qubit processor.

Understanding the fundamental mechanisms underpinning shock initiation is critical to predicting energetic material (EM) safety and performance. Currently, the timescales and pathways by which shock-excited lattice modes transfer energy into specific chemical bonds remains an open question. Towards understanding these mechanisms, our group has previously measured the vibrational energy transfer (VET) pathways in several energetic thin films using broadband, femtosecond transient absorption spectroscopy. However, new technologies are needed to move beyond these thin film surrogates and measure broadband VET pathways in realistic EM morphologies. Herein, we describe a new broadband, femtosecond, attenuated total reflectance spectroscopy apparatus. Performance of the system is benchmarked against published data and the first VET results from a pressed EM pellet are presented. This technology enables fundamental studies of VET dynamics across sample configurations and environments (pressure, temperature, etc .) and supports the potential use of VET studies in the non-destructive surveillance of EM components.

Abstract not provided.

A series of nanopillar compression tests were performed on tungsten as a function of temperature using in situ transmission electron microscopy with localized laser heating. Surface oxidation was observed to form on the pillars and grow in thickness with increasing temperature. Deformation between 850◦C and 1120◦C is facilitated by long-range diffusional transport from the tungsten pillar onto adjacent regions of the Y2O3-stabilized ZrO2 indenter. The constraint imposed by the surface oxidation is hypothesized to underly this mechanism for localized plasticity, which is generally the so-called whisker growth mechanism. The results are discussed in context of the tungsten fuzz growth mechanism in He plasma-facing environments. The two processes exhibit similar morphological features and the conditions under which fuzz evolves appear to satisfy the conditions necessary to induce whisker growth.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The Spent Fuel and Waste Science and Technology (SFWST) Campaign of the U.S. Department of Energy Office of Nuclear Energy, Office of Spent Fuel and Waste Disposition (SFWD), has been conducting research and development on generic deep geologic disposal systems (i.e., geologic repositories). This report describes specific activities in the Fiscal Year (FY) 2021 associated with the Geologic Disposal Safety Assessment (GDSA) Repository Systems Analysis (RSA) work package within the SFWST Campaign. 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 host-rock environments, and to simulate and analyze these generic repository concepts and models using the GDSA Framework toolkit, and other tools as needed.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Virtual prototyping in engineering design rely on modern numerical models of contacting structures with accurate resolution of interface mechanics, which strongly affect the system-level stiffness and energy dissipation due to frictional losses. High-fidelity modeling within the localized interfaces is required to resolve local quantities of interest that may drive design decisions. The high-resolution finite element meshes necessary to resolve inter-component stresses tend to be computationally expensive, particularly when the analyst is interested in response time histories. The Hurty/Craig-Bampton (HCB) transformation is a widely used method in structural dynamics for reducing the interior portion of a finite element model while having the ability to retain all nonlinear contact degrees of freedom (DOF) in physical coordinates. These models may still require many DOF to adequately resolve the kinematics of the interface, leading to inadequate reduction and computational savings. This study proposes a novel interface reduction method to overcome these challenges by means of system-level characteristic constraint (SCC) modes and properly orthogonal interface modal derivatives (POIMDs) for transient dynamic analyses. Both SCC modes and POIMDs are computed using the reduced HCB mass and stiffness matrices, which can be directly computed from many commercial finite element analysis software. Comparison of time history responses to an impulse-type load in a mechanical beam assembly indicate that the interface-reduced model correlates well with the HCB truth model. Localized features like slip and contact area are well-represented in the time domain when the beam assembly is loaded with a broadband excitation. The proposed method also yields reduced-order models with greater critical timestep lengths for explicit integration schemes.

Abstract not provided.

Abstract not provided.

Abstract not provided.

This project focused on providing a fundamental physico-chemical understanding of the coupling mechanisms of corrosion- and radiation-induced degradation at material-salt interfaces in Ni-based alloys operating in emulated Molten Salt Reactor(MSR) environments through the use of a unique suite of aging experiments, in-situ nanoscale characterization experiments on these materials, and multi-physics computational models. The technical basis and capabilities described in this report bring us a step closer to accelerate the deployment of MSRs by closing knowledge gaps related to materials degradation in harsh environments.

Abstract not provided.

Abstract not provided.

This report is a preliminary test plan of the seismic shake table test. The final report will be developed when all decisions regarding the test hardware, instrumentation, and shake table inputs are made. A new revision of this report will be issued in spring of 2022. The preliminary test plan documents the free-field ground motions that will be used as inputs to the shake table, the test hardware, and instrumentation. It also describes the facility at which the test will take place in late summer of 2022.

Abstract not provided.

Machine-learned models, specifically neural networks, are increasingly used as “closures” or “constitutive models” in engineering simulators to represent fine-scale physical phenomena that are too computationally expensive to resolve explicitly. However, these neural net models of unresolved physical phenomena tend to fail unpredictably and are therefore not used in mission-critical simulations. In this report, we describe new methods to authenticate them, i.e., to determine the (physical) information content of their training datasets, qualify the scenarios where they may be used and to verify that the neural net, as trained, adhere to physics theory. We demonstrate these methods with neural net closure of turbulent phenomena used in Reynolds Averaged Navier-Stokes equations. We show the types of turbulent physics extant in our training datasets, and, using a test flow of an impinging jet, identify the exact locations where the neural network would be extrapolating i.e., where it would be used outside the feature-space where it was trained. Using Generalized Linear Mixed Models, we also generate explanations of the neural net (à la Local Interpretable Model agnostic Explanations) at prototypes placed in the training data and compare them with approximate analytical models from turbulence theory. Finally, we verify our findings by reproducing them using two different methods.

Magnetic microscopy with high spatial resolution helps to solve a variety of technical problems in condensed-matter physics, electrical engineering, biomagnetism, and geomagnetism. In this work we used quantum diamond magnetic microscope (QDMM) setups, which use a dense uniform layer of magnetically-sensitive nitrogen-vacancy (NV) centers in diamond to image an external magnetic field using a fluorescence microscope. We used this technique for imaging few-micron ferromagnetic needles used as a physically unclonable function (PUF) and to passively interrogate electric current paths in a commercial 555 timer integrated circuit (IC). As part of the QDMM development, we also found a way to calculate ion implantation recipes to create diamond samples with dense uniform NV layers at the surface. This work opens the possibility for follow-up experiments with 2D magnetic materials, ion implantation, and electronics characterization and troubleshooting.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The purpose of this report is to document improvements in the simulation of 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. A general lack of data suitable for model validation of commercial nuclear canister drying processes necessitates additional, 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.

Abstract not provided.

Abstract not provided.

While computer systems, software applications, and operational technology (OT)/Industrial Control System (ICS) devices are regularly updated through automated and manual processes, there are several unique challenges associated with distributed energy resource (DER) patching. Millions of DER devices from dozens of vendors have been deployed in home, corporate, and utility network environments that may or may not be internet-connected. These devices make up a growing portion of the electric power critical infrastructure system and are expected to operate for decades. During that operational period, it is anticipated that critical and noncritical firmware patches will be regularly created to improve DER functional capabilities or repair security deficiencies in the equipment. The SunSpec/Sandia DER Cybersecurity Workgroup created a Patching Subgroup to investigate appropriate recommendations for the DER patching, holding fortnightly meetings for more than nine months. The group focused on DER equipment, but the observations and recommendations contained in this report also apply to DERMS tools and other OT equipment used in the end-to-end DER communication environment. The group found there were many standards and guides that discuss firmware lifecycles, patch and asset management, and code-signing implementations, but did not singularly cover the needs of the DER industry. This report collates best practices from these standards organizations and establishes a set of best practices that may be used as a basis for future national or international patching guides or standards.

A six-month research effort has advanced the hybrid kinetic-fluid modeling capability required for developing non-thermal warm x-ray sources on Z. The three particle treatments of quasi-neutral, multi-fluid, and kinetic are demonstrated in 1D simulations of an Ar gas puff. The simulations determine required resolutions for the advanced implicit solution techniques and debug hybrid particle treatments with equation-of-state and radiation transport. The kinetic treatment is used in preliminary analysis of the non-Maxwellian nature of a gas target. It is also demonstrates the sensitivity of the cyclotron and collision frequencies in determining the transition from thermal to non-thermal particle populations. Finally, a 2D Ar gas puff simulation of a Z shot demonstrates the readiness to proceed with realistic target configurations. The results put us on a very firm footing to proceed to a full LDRD which includes continued development transition criteria and x-ray yield calculation.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Current methods for stochastic media transport are either computationally expensive or, by nature, approximate. Moreover, none of the well-developed, benchmarked approximate methods can compute the variance caused by the stochastic mixing, a quantity especially important to safety calculations. Therefore, we derive and apply a new conditional probability function (CPF) for use in the recently developed stochastic media transport algorithm Conditional Point Sampling (CoPS), which 1) leverages the full intra-particle memory of CoPS to yield errorless computation of stochastic media outputs in 1D, binary, Markovian-mixed media, and 2) leverages the full inter-particle memory of CoPS and the recently developed Embedded Variance Deconvolution method to yield computation of the variance in transport outputs caused by stochastic material mixing. Numerical results demonstrate errorless stochastic media transport as compared to reference benchmark solutions with the new CPF for this class of stochastic mixing as well as the ability to compute the variance caused by the stochastic mixing via CoPS. Using previously derived, non-errorless CPFs, CoPS is further found to be more accurate than the atomic mix approximation, Chord Length Sampling (CLS), and most of memory-enhanced versions of CLS surveyed. In addition, we study the compounding behavior of CPF error as a function of cohort size (where a cohort is a group of histories that share intra-particle memory) and recommend that small cohorts be used when computing the variance in transport outputs caused by stochastic mixing.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Wellbore integrity is a significant problem in the U.S. and worldwide, which has serious adverse environmental and energy security consequences. Wells are constructed with a cement barrier designed to last about 50 years. Indirect measurements and models are commonly used to identify wellbore damage and leakage, often producing subjective and even erroneous results. The research presented herein focuses on new technologies to improve monitoring and detection of wellbore failures (leaks) by developing a multi-step machine learning approach to localize two types of thermal defects within a wellbore model, a prototype mechatronic system for automatically drilling small diameter holes of arbitrary depth to monitor the integrity of oil and gas wells in situ, and benchtop testing and analyses to support the development of an autonomous real-time diagnostic tool to enable sensor emplacement for monitoring wellbore integrity. Each technology was supported by experimental results. This research has provided tools to aid in the detection of wellbore leaks and significantly enhanced our understanding of the interaction between small-hole drilling and wellbore materials.

This SAND report fulfills the completion requirements for the ASC Physics and Engineering Modeling Level 2 Milestone 7836 during Fiscal Year 2021. The Sandia Simplified potential energy clock (SPEC) non-linear viscoelastic constitutive model was developed to predict a whole host of polymer glass physical behaviors in order to provide a tool to assess the effects of stress on these materials over their lifecycle. Polymer glasses are used extensively in applications such as electronics packaging, where encapsulants and adhesives can be critical to device performance. In this work, the focus is on assessing the performance of the model in predicting material evolution associated with long-term physical aging, an area that the model has not been fully vetted in. These predictions are key to utilizing models to help demonstrate electronics packaging component reliability over decades long service lives, a task that is very costly and time consuming to execute experimentally. The initiating hypothesis for the work was that a model calibration process can be defined that enables confidence in physical aging predictions under ND relevant environments and timescales without sacrificing other predictive capabilities. To test the hypothesis, an extensive suite of calibration and aging data was assembled from a combination of prior work and collaborating projects (Aging and Lifetimes as well as the DoD Joint Munitions Program) for two mission relevant epoxy encapsulants, 828DGEBA/DEA and 828DGEBA/T403. Multiple model calibration processes were developed and evaluated against the entire set of data for each material. A qualitative assessment of each calibration's ability to predict the wide range of aging responses was key to ranking the calibrations against each other. During this evaluation, predictions that were identified as non-physical, i.e., demonstrated something that was qualitatively different than known material behavior, were heavily weighted against the calibration performance. Thus, unphysical predictions for one aspect of aging response could generate a lower overall rating for a calibration process even if that process generated better quantitative predictions for another aspect of aging response. This insurance that all predictions are qualitatively correct is important to the overall aim of utilizing the model to predict residual stress evolution, which will depend on the interplay amongst the different material aging responses. The DSC-focused calibration procedure generated the best all-around aging predictions for both materials, demonstrating material models that can qualitatively predict the whole host of different physical aging responses that have been measured. This step forward in predictive capability comes from an unanticipated source, utilization of calorimetry measurements to specify model parameters. The DSC-focused calibration technique performed better than compression-focused techniques that more heavily weigh measurements more closely related to the structural responses to be predicted. Indeed, the DSC-focused calibration procedure was only possible due to recent incorporation of the enthalpy and heat capacity features into SPEC that was newly verified during this L2 milestone. Fundamentally similar aspects of the two material model calibrations as well as parametric studies to assess sensitives of the aging predictions are discussed within the report. A perspective on the next steps to the overall goal of residual stress evolution predictions under stockpile conditions closes the report.

Abstract not provided.

Abstract not provided.

Over the past four years, an informal working group has developed to investigate existing sensitivity analysis methods, examine new methods, and identify best practices. The focus is on the use of sensitivity analysis in case studies involving geologic disposal of spent nuclear fuel or nuclear waste. To examine ideas and have applicable test cases for comparison purposes, we have developed multiple case studies. Four of these case studies are presented in this report: the GRS clay case, the SNL shale case, the Dessel case, and the IBRAE groundwater case. We present the different sensitivity analysis methods investigated by various groups, the results obtained by different groups and different implementations, and summarize our findings.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Detonations and flames are characterized by three-dimensional (3D) temperature fields, yet state-of- the-art temperature measurement techniques yield information at a point or along a line. The goal of the research documented here was to combine ultrafast laser spectroscopy and structured illumination to deliver an unprecedented measurement capability—three-dimensional, instantaneous temperature measurements in a gas-phase volume. To achieve this objective, different parts of the proposed technique were developed and tested independently. Structured illumination was used to image particulate matter (soot) in a turbulent flame at multiple planes using a single laser pulse and a single camera. Emission spectroscopy with structured detection was demonstrated for emission- based measurements of explosives with enhance dimensionality. Finally, an instrument for multi- planar laser-based temperature measurement technique was developed. Structured illumination techniques will continue to be developed for multi-dimensional and multi-parameter measurements. These new measurement capabilities will be important for heat transfer and fluid dynamic research areas.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

This report describes the risk-informed technical elements that will contribute to a defense-in-depth assessment for cybersecurity. Risk-informed cybersecurity must leverage the technical elements of a risk-informed approach appropriately in order to evaluate cybersecurity risk insights. HAZCADS and HAZOP+ are suitable methodologies to model the connection between digital harm and process hazards. Risk assessment modeling needs to be expanded beyond HAZCADS and HAZOP+ to consider the sequence of events that lead to plant consequences. Leveraging current practices in PRA can lead to categorization of digital assets and prioritizing digital assets commensurate with the risk. Ultimately, the culmination of cyber hazard methodologies, event sequence modeling, and digital asset categorization will facilitate a defense-in-depth assessment of cybersecurity.

Abstract not provided.

Scientific applications run on high-performance computing (HPC) systems are critical for many national security missions within Sandia and the NNSA complex. However, these applications often face performance degradation and even failures that are challenging to diagnose. To provide unprecedented insight into these issues, the HPC Development, HPC Systems, Computational Science, and Plasma Theory & Simulation departments at Sandia crafted and completed their FY21 ASC Level 2 milestone entitled "Integrated System and Application Continuous Performance Monitoring and Analysis Capability." The milestone created a novel integrated HPC system and application monitoring and analysis capability by extending Sandia's Kokkos application portability framework, Lightweight Distributed Metric Service (LDMS) monitoring tool, and scalable storage, analysis, and visualization pipeline. The extensions to Kokkos and LDMS enable collection and storage of application data during run time, as it is generated, with negligible overhead. This data is combined with HPC system data within the extended analysis pipeline to present relevant visualizations of derived system and application metrics that can be viewed at run time or post run. This new capability was evaluated using several week-long, 290-node runs of Sandia's ElectroMagnetic Plasma In Realistic Environments ( EMPIRE ) modeling and design tool and resulted in 1TB of application data and 50TB of system data. EMPIRE developers remarked this capability was incredibly helpful for quickly assessing application health and performance alongside system state. In short, this milestone work built the foundation for expansive HPC system and application data collection, storage, analysis, visualization, and feedback framework that will increase total scientific output of Sandia's HPC users.

Abstract not provided.

Abstract not provided.

Abstract not provided.

We present a new experimental methodology for detailed experimental investigations of depolymerization reactions over solid catalysts. This project aims to address a critical need in fundamental research on chemical upcycling of polymers – the lack of rapid, sensitive, isomerselective probing techniques for the detection of reaction intermediates and products. Our method combines a heterogeneous catalysis reactor for the study of multiphase (gas/polymer melt/solid) systems, coupled to a vacuum UV photoionization time-of-flight mass spectrometer. This apparatus draws on our expertise in probing complex gas-phase chemistry and enables highthroughput, detailed chemical speciation measurements of the gas phase above the catalyst, providing valuable information on the heterogeneous catalytic reactions. Using this approach, we investigated the depolymerization of high-density polyethylene (HDPE) over Ir-doped zeolite catalysts. We showed that the product distribution was dominated by low-molecular weight alkenes with terminal C=C double bonds and revealed the presence of many methyl-substituted alkenes and alkanes, suggesting extensive methyl radical chemistry. In addition, we investigated the fundamental reactivity of model oligomer molecules n-butane and isobutane over ZSM-5 zeolites. We demonstrated the first direct detection of methyl radical intermediates, confirming the key role of methyl in zeolite-catalyzed activation of alkanes. Our results show the potential of this experimental method to achieve deep insight into the complex depolymerization reactions and pave the way for detailed mechanistic studies, leading to increased fundamental understanding of key processes in chemical upcycling of polymers.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The effect of extreme waves on the coastal community includes inundation, loss of habitats, increasing shoreline erosion, and increasing risks to coastal infrastructures (e.g., ports, breakwaters, oil and gas platforms), important for supporting coastal resilience. The coastal communities along the US Gulf of Mexico are very low-lying, which makes the region particularly vulnerable to impacts of extreme waves generated by storm events. We propose assessing and mapping the risks from extreme waves for the Gulf of Mexico coast to support coastal resiliency planning. The risks will be assessed by computing n-year recurring wave height (e.g., 1, 5, 50, 100-year) using 32-year wave hindcast data and various extreme value analysis techniques including Peak- Over-Threshold and Annual Maxima method. The characteristics of the extreme waves, e.g., relations between the mean and extreme wave climates, directions associated with extreme waves, will be investigated. Hazard maps associated with extreme wave heights at different return periods will be generated to help planners identify potential risks and envision places that are less susceptible to future storm damage.

Probabilistic and Bayesian neural networks have long been proposed as a method to incorporate uncertainty about the world (both in training data and operation) into artificial intelligence applications. One approach to making a neural network probabilistic is to leverage a Monte Carlo sampling approach that samples a trained network while incorporating noise. Such sampling approaches for neural networks have not been extensively studied due to the prohibitive requirement of many computationally expensive samples. While the development of future microelectronics platforms that make this sampling more efficient is an attractive option, it has not been immediately clear how to sample a neural network and what the quality of random number generation should be. This research aimed to start addressing these two fundamental questions by examining basic “off the shelf” neural networks can be sampled through a few different mechanisms (including synapse “dropout” and neuron “dropout”) and examine how these sampling approaches can be evaluated both in terms of evaluating algorithm effectiveness and the required quality of random numbers.

Abstract not provided.

Abstract not provided.

Abstract not provided.