Publications

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

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

Abstract not provided.

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

Abstract not provided.

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

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

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

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

This SAND report collects the results from the LDRD project “SHAZAM”, which aimed to push the limits of performance for self-healing, self-assembling power systems whose sectionalizing and load-control agents rely on local measurements only (i.e., only what they can measure at their own terminals, with no data sharing between agents). This work includes self-networking microgrids. The key objectives of this work were a) to demonstrate how high the performance of local-measurement-only self-assembling power systems can be; and b) to solve certain technical problems associated with such systems, such as their inability to prevent the accidental formation of closed loops and their tendency to thermally overload some conductors. “SHAZAM” investigators a) demonstrated that the performance of such systems can be surprisingly high, b) demonstrated that such systems are quite robust to all kinds of variations, and c) developed and demonstrated solutions to several key challenges associated with this type of system.

Abstract not provided.

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

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

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

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

Both human subject experiments and computational, modeling and simulations have been used to study detection of deception. This work aims to combine these two methods by integrating empirically-derived information (from human subject experiments) into agent-based models to generate novel insights into the complex problems of detection of disinformation content. Computational experiments are used to simulate across multiple scenarios for evaluation and decision-making regarding the validity of potentially deceptive scientific documents. Factors influencing the human agent behaviors in the model were identified through a human subject experiment that was conducted to evaluate and characterize decision making related to disinformation discernment. Correlation and regression analyses were used to translate insights from the human subjects experiment to inform the parameterization of agent features and scenario development. Three scenarios were evaluated with the agent-based models to help evaluate the replicability of the simulations (validation analysis) and assess the influence of human agent and document features (sensitivity analyses). A replication of the human participant experiment demonstrated that the agent-based simulations compare favorably to empirical findings. The agent-based modeling was then used to conduct sensitivity analysis on the accuracy of deception detection as a function of document proportions and human agent features. Results indicate that precision values are adversely impacted when the proportion of deceptive documents is lower in the overall sample, whereas recall values are more sensitive to changes in human agent features. These findings indicate important nuances in accuracy evaluations that should be further considered (including consideration of potential alternate metrics) in future agent-based models of disinformation. Additional areas for future exploration include extension of simulations to consider other ways to align the agent-based model design with psychological theory and inclusion of agent-agent interactions, especially as it pertains to sharing of scientific information within an organizational context.

Sandia’s Computational Engine for Particle Transport for Radiation Effects (SCEPTRE) is a computer code that solves the linear Boltzmann transport equation, particularly targeting coupled photon-electron problems. It uses unstructured finite element meshes in space, multigroup in energy, and discrete ordinates (S_n) or other methods in angle. SCEPTRE uses an xml-based input file to specify the problem. This report documents the options and syntax of that input file.

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

This report represents the milestone deliverable M2SF-23SN010301072 “Evaluation of Nuclear Spent Fuel Disposal in Clay-Bearing Rock - Process Model Development and Experimental Studies” The report provides a status update of FY23 activities for the work package Argillite Disposal work packages for the DOE-NE Spent Fuel Waste Form Science and Technology (SFWST) Program. Clay-rich geological media (often referred as shale or argillite) are among the most abundant type of sedimentary rock near the Earth’s surface. Argillaceous rock formations have the following advantageous attributes for deep geological nuclear waste disposal: widespread geologic occurrence, found in stable geologic settings, low permeability, self-sealing properties, low effective diffusion coefficient, high sorption capacity, and have the appropriate depth and thickness to host nuclear waste repository concepts. The DOE R&D program under the Spent Fuel Waste Science Technology (SFWST) campaign has made key progress (through experiment, modeling, and testing) in the study of chemical and physical phenomena that could impact the long-term safety assessment of heat-generating nuclear waste disposition in clay/shale/argillaceous rock. International collaboration activities comprising field-scale heater tests, field data monitoring, and laboratory-scale experiments provide key information on changes to the engineered barrier system (EBS) material exposed high thermal loads. Moreover, consideration of direct disposal of large capacity dual-purpose canisters (DPCs) as part of the back-end SNF waste disposition strategy has generated interest in improving our understanding of the effects of elevated temperatures on the engineered barrier system (EBS) design concepts. Chemical and structural analyses of sampled bentonite material from laboratory tests at elevated temperatures are key to the characterization of thermal effects affecting bentonite clay barrier performance. The knowledge provided by these experiments is crucial to constrain the extent of sacrificial zones in the EBS design during the thermal period. Thermal, hydrologic, mechanical, and chemical (THMC) data collected from heater tests and laboratory experiments have been used in the development, validation, and calibration of THMC simulators to model near-field coupled processes. This information leads to the development of simulation approaches to assess issues on coupled processes involving porous media flow, transport, geomechanical phenomena, chemical interactions with barrier/geologic materials, and the development of EBS concepts. These lines of knowledge are central to the design of deep geological backfilled repository concepts where temperature plays a key role in the EBS behavior, potential interactions with host rock, and long-term performance in the safety assessment.

Pressure-shear plate impact experiments were performed to quantify flow strength of wrought, as-built additively manufactured (AM), and heat-treated and recrystallized AM 304 L stainless steel (SS304L) under combined loading. Impact velocities spanned between 0.03 and 0.24 mm/μs, resulting in corresponding pressures of 0.62–5.93 GPa. Flow strength measurements are comparable for the sample variants across the studied loading conditions; however, shear wave structures significantly differ between sample type. Microstructurally aware simulations indicate local strain differences attributed to anisotropic elastic constants of large grains (~1 mm) in the as-built and heat-treated AM may impede the ability to uniformly transmit a shear wave.

We combine crossed-beam velocity map imaging with high-level ab initio/transition state theory modeling of the reaction of S(3P) with 1,3-butadiene and isoprene under single collision conditions. For the butadiene reaction, we detect both H and H2 loss from the initial adduct, and from reaction with isoprene, we see both H loss and methyl loss. Theoretical calculations confirm these arise following intersystem crossing to the singlet surface forming long-lived intermediates. For the butadiene reaction, these lose H2 to form thiophene as the dominant channel, H to form the detected 2H-thiophenyl radical, or ethene, giving thioketene. For isoprene, additional reaction products are suggested by theory, including the observed H and methyl loss radicals, but also methyl thiophene, thioformaldehyde, and thioketene. The results for S(3P) + 1,3-butadiene, showing direct cyclization to the aromatic product and yielding few bimolecular product channels, are in striking contrast to those for the analogous O(3P) reaction.

This paper explores the utility of organizational system modeling frameworks to provide valuable insight into information flows within organizations and subsequently the opportunities for increasing resilience against disinformation campaigns targeting the system's ability to utilize information within its decision making. Disinformation is a growing challenge for many organizations and in recent years has created delay in decision making. Here the paper has utilized the viable systems model (VSM) to characterize organizational systems and used this approach to outline potential subsystem requirements to promote resilience of the system. The results of this paper can support the development of simulations and models considering the human elements within the system as well as support the development of quantitative measures of resilience.

The US Department of Energy (DOE) is investigating the use of different materials that could be used to fill the void space inside a dual-purpose canister (DPC) loaded with spent nuclear fuel (SNF) just before it is emplaced in a deep geologic repository. The purpose of adding filler material is to maintain subcritical conditions in the repository during the postclosure period, which can span up to 1,000,000 years. Several types of materials have been proposed, including metals, cements, particulates, and glass. Part of this investigation addresses how the presence of filler material inside a DPC will affect the performance of the repository with respect to the repository features; the consequences of events that may occur; and the multiple thermal, hydrologic, chemical, and mechanical processes that may occur in a deep geologic repository over long timescales. This report describes some of the filler materials that have been proposed and studied; identifies 11 features, 6 events, and 25 processes that may be affected by the presence of filler materials; and discusses the effects that may require consideration for each feature, event, or process. The results of this study can be used to direct appropriate research and to develop suitable models if the DOE decides to use fillers to maintain subcritical conditions in DPCs used to dispose of SNF.

Here, a method for the nonintrusive and structure-preserving model reduction of canonical and noncanonical Hamiltonian systems is presented. Based on the idea of operator inference, this technique is provably convergent and reduces to a straightforward linear solve given snapshot data and gray-box knowledge of the system Hamiltonian. Examples involving several hyperbolic partial differential equations show that the proposed method yields reduced models which, in addition to being accurate and stable with respect to the addition of basis modes, preserve conserved quantities well outside the range of their training data.

This paper is concerned with goal-oriented a posteriori error estimation for nonlinear functionals in the context of nonlinear variational problems solved with continuous Galerkin finite element discretizations. A two-level, or discrete, adjoint-based approach for error estimation is considered. The traditional method to derive an error estimate in this context requires linearizing both the nonlinear variational form and the nonlinear functional of interest which introduces linearization errors into the error estimate. In this paper, we investigate these linearization errors. In particular, we develop a novel discrete goal-oriented error estimate that accounts for traditionally neglected nonlinear terms at the expense of greater computational cost. We demonstrate how this error estimate can be used to drive mesh adaptivity. Here, we show that accounting for linearization errors in the error estimate can improve its effectivity for several nonlinear model problems and quantities of interest. We also demonstrate that an adaptive strategy based on the newly proposed estimate can lead to more accurate approximations of the nonlinear functional with fewer degrees of freedom when compared to uniform refinement and traditional adjoint-based approaches.

The formation of a stress corrosion crack (SCC) in the canister wall of a dry cask storage system (DCSS) has been identified as a potential issue for the long-term storage of spent nuclear fuel. The presence of an SCC in a storage system could represent a through-wall flow path from the canister interior to the environment. Modern, vertical DCSSs are of particular interest due to the commercial practice of using higher backfill pressures in the canister, up to approximately 800 kPa, compared to their horizontal counterparts. This pressure differential offers a relatively high driving potential for blowdown of any particulates that might be present in the canister. In this study, the rates of gas flow and aerosol transmission of a spent fuel surrogate through an engineered microchannel with dimensions representative of an SCC were evaluated experimentally using coupled mass flow and aerosol analyzers. The microchannel was formed by mating two gage blocks with a linearly tapering slot orifice nominally 13 μm (0.005 in.) tall on the upstream side and 25 μm (0.0010 in.) tall on the downstream side. The orifice is 12.7 mm (0.500 in.) wide by 8.86 mm (0.349 in.) long (flow length). Surrogate aerosols of cerium oxide, CeO₂, were seeded and mixed with either helium or air inside a pressurized tank. The aerosol characteristics were measured immediately upstream and downstream of the simulated SCC at elevated and ambient pressures, respectively. These data sets are intended to add to previous testing that characterized SCCs under well-controlled boundary conditions through the inclusion of testing improvements that establish initial conditions in a more consistent way. While the engineered microchannel has dimensions similar to actual SCCs, it does not reproduce the tortuous path the aerosol laden flow would have to traverse for eventual transmission. SCCs can be rapidly grown in a laboratory setting given the right conditions, and initial characterization and clean-flow testing has begun on lab grown crack samples provided to Sandia National Laboratories (SNL). Many such samples are required to produce statistically relevant transmission results, and SNL is developing a procedure to produce samples in welded steel plates. These ongoing testing efforts are focused on understanding the evolution in both size and quantity of a hypothetical release of aerosolized spent fuel particles from failed fuel to the canister interior and ultimately through an SCC.

Grid-scale battery energy storage systems (BESSs) are vulnerable to false data injection attacks (FDIAs), which could be used to disrupt state of charge (SoC) estimation. Inaccurate SoC estimation has negative impacts on system availability, reliability, safety, and the cost of operation. In this article a combination of a Cumulative Sum (CUSUM) algorithm and an improved input noise-aware extended Kalman filter (INAEKF) is proposed for the detection and identification of FDIAs in the voltage and current sensors of a battery stack. The series-connected stack is represented by equivalent circuit models, the SoC is modeled with a charge reservoir model and the states are estimated using the INAEKF. Further, the root mean squared error of the states’ estimation by the modified INAEKF was found to be superior to the traditional EKF. By employing the INAEKF, this article addresses the research gap that many state estimators make asymmetrical assumptions about the noise corrupting the system. Additionally, the INAEKF estimates the input allowing for the identification of FDIA, which many alternative methods are unable to achieve. The proposed algorithm was able to detect attacks in the voltage and current sensors in 99.16% of test cases, with no false positives. Utilizing the INAEKF compared to the standard EKF allowed for the identification of FDIA in the input of the system in 98.43% of test cases.

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.

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.

This user’s guide documents capabilities in Sierra/SolidMechanics which remain “in-development” and thus are not tested and hardened to the standards of capabilities listed in Sierra/SM 5.16 User’s Guide. Capabilities documented herein are available in Sierra/SM for experimental use only until their official release. These capabilities include, but are not limited to, novel discretization approaches such as the conforming reproducing kernel (CRK) method, numerical fracture and failure modeling aids such as the extended finite element method (XFEM) and J-integral, explicit time step control techniques, dynamic mesh rebalancing, as well as a variety of new material models and finite element formulations.

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.

Using compressive mechanical forces, such as pressure, to induce crystallographic phase transitions and mesostructural changes while modulating material properties in nanoparticles (NPs) is a unique way to discover new phase behaviors, create novel nanostructures, and study emerging properties that are difficult to achieve under conventional conditions. In recent decades, NPs of a plethora of chemical compositions, sizes, shapes, surface ligands, and self-assembled mesostructures have been studied under pressure by in-situ scattering and/or spectroscopy techniques. As a result, the fundamental knowledge of pressure-structure-property relationships has been significantly improved, leading to a better understanding of the design guidelines for nanomaterial synthesis. In the present review, we discuss experimental progress in NP high-pressure research conducted primarily over roughly the past four years on semiconductor NPs, metal and metal oxide NPs, and perovskite NPs. We focus on the pressure-induced behaviors of NPs at both the atomic- and mesoscales, inorganic NP property changes upon compression, and the structural and property transitions of perovskite NPs under pressure. We further discuss in depth progress on molecular modeling, including simulations of ligand behavior, phase-change chalcogenides, layered transition metal dichalcogenides, boron nitride, and inorganic and hybrid organic-inorganic perovskites NPs. These models now provide both mechanistic explanations of experimental observations and predictive guidelines for future experimental design. We conclude with a summary and our insights on future directions for exploration of nanomaterial phase transition, coupling, growth, and nanoelectronic and photonic properties.

Intimately intertwined atomic and electronic structures of point defects govern diffusion-limited corrosion and underpin the operation of optoelectronic devices. For some materials, complex energy landscapes containing metastable defect configurations challenge first-principles modeling efforts. Here, we thoroughly reevaluate native point defect geometries for the illustrative case of α-Al2O3 by comparing three methods for sampling candidate geometries in density functional theory calculations: displacing atoms near a naively placed defect, initializing interstitials at high-symmetry points of a Voronoi decomposition, and Bayesian optimization. We find symmetry-breaking distortions for oxygen vacancies in some charge states, and we identify several distinct oxygen split-interstitial geometries that help explain literature discrepancies involving this defect. We also report a surprising and, to our knowledge, previously unknown trigonal geometry favored by aluminum interstitials in some charge states. These new configurations may have transformative impacts on our understanding of defect migration pathways in aluminum-oxide scales protecting metal alloys from corrosion. Overall, the Voronoi scheme appears most effective for sampling candidate interstitial sites because it always succeeded in finding the lowest-energy geometry identified in this study, although no approach found every metastable configuration. Finally, we show that the position of defect levels within the band gap can depend strongly on the defect geometry, underscoring the need to conduct careful searches for ground-state geometries in defect calculations.

Pultrusion manufacturing of fiber reinforced polymers has been shown to yield some of the highest mechanical properties for unidirectional composites, having a high degree of fiber alignment with consistent performance. Pultrusions offer a low-cost manufacturing approach for producing unidirectional composites with a constant cross-section and are used in many applications, including spar caps of wind turbine blades. However, as an intermediate processing step for wind blades, the additional cost of manufacturing pultrusions must be accompanied by sufficient increases in mechanical performance and system benefits. Wind turbine blades are manufactured using vacuum-assisted resin transfer molding with infused unidirectional fiberglass or carbon pultrusions for the spar cap. Infused fiberglass composites are among the most cost-effective structural materials available and replacing this material in the cost-driven wind industry has proven challenging, where infused fiberglass spar caps are still the predominant material system in use. To evaluate alternative material systems in a pultruded composite form, it is necessary to understand the costs for this additional manufacturing step which are shown to add 33%–55% on top of the material costs. A pultrusion cost model has been developed and used to quantify cost sensitivities to various processing parameters. The mechanical performance for pultruded composites is improved versus resin-infusion manufacturing with a 17% increase in design strength at a constant fiber volume fraction, but also enables higher achievable fiber volume fractions. The cost-specific mechanical performance is compared as a function of processing parameters for pultruded composites to identify the opportunities for alternative material and manufacturing approaches for wind turbine spar caps. Finally, four materials are compared in a representative wind turbine blade model to assess the performance of pultruded carbon fiber systems and pultruded fiberglass relative to infused fiberglass, where the pultruded systems produce lower weight blades with various cost distinctions.

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.

The objective of the crystalline disposal work packages is to advance our understanding of long-term disposal of used fuel in crystalline rocks and to develop necessary experimental and computational capabilities to evaluate various disposal concepts in such media.

Combinatorial research, the incorporation of multiple domains in a unified research agenda, is a strong contributor to the growing corpus of scientific knowledge and technological advancements worldwide. In 2019, a study team at Sandia National Laboratories (Sandia, the Labs) used a systems approach to understand if and how combinatorial research agendas were playing out at Sandia, one of America’s premiere national security research venues. The study team used the data collection effort described in this report to ground the discussion of the broad social environment and particular organizational environments within which combinatorial research agendas are developed, as described in the full study. The team interviewed twenty-five staff members engaged in combinatorial research at Sandia in New Mexico and California during the months of June – September 2019. Analysis of this corpus of ethnographic data, combined with knowledge drawn from relevant literature, concluded that there is an individual type who would be most likely to engage in combinatoric research, described by both demographic and psychographic components. This type demonstrates both intellectual depth and the curiosity which leads to breadth. The analysis also showed that Sandia as an organization and as perceived by the respondents, set up tension for the combinatorial researcher. While Sandia was generally agnostic towards combinatorial research, that agnostic posture depended on whether the researcher was able to fulfill all her customer obligations – obligations that are structured primarily in transactional relationships with customers with relatively short time horizons. This report concludes with suggestions for additional research in the ethnographic domain.

The hydrogen sorption properties of single-phase bcc (TiVNb)100-xCrx alloys (x = 0-35) are reported. All alloys absorb hydrogen quickly at 25 °C, forming fcc hydrides with storage capacity depending on the Cr content. A thermodynamic destabilization of the fcc hydride is observed with increasing Cr concentration, which agrees well with previous compositional machine learning models for metal hydride thermodynamics. The steric effect or repulsive interactions between Cr-H might be responsible for this behavior. The cycling performances of the TiVNbCr alloy show an initial decrease in capacity, which cannot be explained by a structural change. Pair distribution function analysis of the total X-ray scattering on the first and last cycled hydrides demonstrated an average random fcc structure without lattice distortion at short-range order. If the as-cast alloy contains a very low density of defects, the first hydrogen absorption introduces dislocations and vacancies that cumulate into small vacancy clusters, as revealed by positron annihilation spectroscopy. Finally, the main reason for the capacity drop seems to be due to dislocations formed during cycling, while the presence of vacancy clusters might be related to the lattice relaxation. Having identified the major contribution to the capacity loss, compositional modifications to the TiVNbCr system can now be explored that minimize defect formation and maximize material cycling performance.

Barium titanate (BTO) is a widely researched ferroelectric useful for energy storage. While BTO’s surface chemistry is commonly studied using density functional theory, little has been published on the TiO₂ surface. Here, we determined that BTO’s surface response can be decoupled from the ferroelectric response by using a pre-optimized ferroelectric slab and allowing only the top three atomic z-layers to respond to ligand binding. Multiple favorable binding modes were identified for hydrogen, hydroxyl, water, and tert-butyl phosphonic acid on BTO’s TiO₂ surface. Of these ligands, tBuPA dominates surface binding with binding energies as low as -2.61 eV for its nine configurations.

Immune checkpoint immunotherapy (ICI) can re-activate immune reactions against neoantigens, leading to remarkable remission in cancer patients. Nevertheless, only a minority of patients are responsive to ICI, and approaches for prediction of responsiveness are needed to improve the success of cancer treatments. While the tumor mutational burden (TMB) correlates positively with responsiveness and survival of patients undergoing ICI, the influence of the subcellular localizations of the neoantigens remains unclear. Here, we demonstrate in both a mouse melanoma model and human clinical datasets of 1,722 ICI-treated patients that a high proportion of membrane-localized neoantigens, particularly at the plasma membrane, correlate with responsiveness to ICI therapy and improved overall survival across multiple cancer types. We further show that combining membrane localization and TMB analyses can enhance the predictability of cancer patient response to ICI. Our results may have important implications for establishing future clinical guidelines to direct the choice of treatment toward ICI.

Recent advances in the growth of aluminum scandium nitride films on silicon suggest that this material platform could be applied for quantum electromechanical applications. Here, we model, fabricate, and characterize microwave frequency silicon phononic delay lines with transducers formed in an adjacent aluminum scandium nitride layer to evaluate aluminum scandium nitride films, at 32% scandium, on silicon interdigital transducers for piezoelectric transduction into suspended silicon membranes. We achieve an electromechanical coupling coefficient of 2.7% for the extensional symmetric-like Lamb mode supported in the suspended material stack and show how this coupling coefficient could be increased to at least 8.5%, which would further boost transduction efficiency and reduce the device footprint. The one-sided transduction efficiency, which quantifies the efficiency at which the source of microwave photons is converted to microwave phonons in the silicon membrane, is 10% at 5 GHz at room temperature and, as we discuss, there is a path to increase this toward near-unity efficiency based on a combination of modified device design and operation at cryogenic temperatures.

This FY2023 report is the second update to the Disposal Research (DR) Research and Development (R&D) 5-year plan for the Spent Fuel and Waste Science and Technology (SFWST) Campaign DR R&D activities. In the planning for FY2020 in the U.S. Department of Energy (DOE) NE-81 SFWST Campaign, the DOE requested development of a high-level summary plan for activities in the DR R&D program for the next five (5)-year period, with periodic updates to this summary plan. The DR R&D 5-year plan was provided to the DOE based initially on the FY2020 priorities and program structure (initial 2020 version of this 5-year plan) and provides a strategic summary guide to the work within the DR R&D technical areas (Control Accounts, CA), focusing on the highest priority technical thrusts. This 5-year plan is a living document (planned to be updated periodically) that provides review of SFWST R&D accomplishments (as seen on the 2021 revision of this 5-year plan), describes changes to technical R&D prioritization based on (a) progress in each technical area (including external technical understanding) with specific accomplishments and (b) any changes in SFWST Campaign objectives and/or funding levels (i.e., Program Direction). Updates to this 5-year plan include the DR R&D adjustments to high-priority knowledge gaps to be investigated in the near-term, as well as the updated longer-term DR R&D directions for the program activities. This plan fulfills the Milestone M2SF23SN010304083 in DR Work Package (WP) SF-23SN01030408 (GDSA - Framework Development – SNL).

This report summarizes the activities performed by Sandia National Laboratories in FY23 to identify and test coating materials for the prevention, mitigation, and/or repair of potential chloride-induced stress corrosion cracking in spent nuclear fuel dry storage canisters. This work continues efforts by Sandia National Laboratories that are summarized in previous reports from FY20 through FY22 on the same topic. In FY23, Sandia National Laboratories, in collaboration with five industry partners through a memorandum of understanding, evaluated the physical, mechanical, and corrosion-resistance properties of eight different coating systems. The evaluation included thermal and radiation environments relevant to various time periods of storage for spent nuclear fuel canisters. The coating systems include polymeric (polyetherketoneketone, modified polyimide/polyurea, modified phenolic resin, epoxy), organic/inorganic ceramic hybrids (silane-based polyurethane hybrid and a quasi-ceramic sol-gel polyurethane hybrid), and coatings utilizing a Zn-rich primer applied to stainless steel coupons. The results and implications of these tests are summarized in this report. These analyses will be used to identify the most effective coatings for potential use on spent nuclear fuel dry storage canisters and to identify specific needs for further optimization of coating technologies for application on spent nuclear fuel canisters.

Here, we present a classical interatomic force field, silica-DDEC, to describe the interactions of amorphous and crystalline silica surfaces, parametrized using density functional theory-based charges. Charge schemes for silica surfaces were developed using the density-derived electrostatic and chemical (DDEC) method, which reproduces atomic charges of the periodic models as well as the electrostatic potential away from the atom sites. Lennard–Jones parameters were determined by requiring the correct description of (i) the amorphous silica density, coordination defects, and local coordination geometry, relative to experimental measurements, and (ii) water-silica interatomic distances compared with ab initio results. Deprotonated surface silanol sites are also described within the model based on DDEC charges. The result is a general electronic structure-derived model for describing fully flexible amorphous and crystalline silica surfaces and interactions of liquids with silica surfaces of varying structure and protonation state.

Chemically robust, low-power sensors are needed for the direct electrical detection of toxic gases. Metal-organic frameworks (MOFs) offer exceptional chemical and structural tunability to meet this challenge, though further understanding is needed regarding how coadsorbed gases influence or interfere with the electrical response. To probe the influence of competitive gases on trace NO2 detection in a simulated flue gas stream, a combined structure-property study integrating synchrotron powder diffraction and pair distribution function analyses was undertaken, to elucidate how structural changes associated with gas binding inside Ni-MOF-74 pores correlate with the electrical response from Ni-MOF-74-based sensors. Data were evaluated for 16 gas combinations of N2, NO2, SO2, CO2, and H2O at 50 °C. Fourier difference maps from a rigid-body Rietveld analysis showed that additional electron density localized around the Ni-MOF-74 lattice correlated with large decreases in Ni-MOF-74 film resistance of up to a factor of 6 × 103, observed only when NO2 was present. These changes in resistance were significantly amplified by the presence of competing gases, except for CO2. Without NO2, H2O rapidly (<120 s) produced small (1-3×) decreases in resistance, though this effect could be differentiated from the slower adsorption of NO2 by the evaluation of the MOF’s capacitance. Furthermore, samples exposed to H2O displayed a significant shift in lattice parameters toward a larger lattice and more diffuse charge density in the MOF pore. Evaluating the Ni-MOF-74 impedance in real time, NO2 adsorption was associated with two electrically distinct processes, the faster of which was inhibited by competitive adsorption of CO2. Together, this work points to the unique interaction of NO2 and other specific gases (e.g., H2O, SO2) with the MOF’s surface, leading to orders of magnitude decrease in MOF resistance and enhanced NO2 detection. Understanding and leveraging these coadsorbed gases will further improve the gas detection properties of MOF materials.

Here we have observed the behavior of striations caused by ionization waves propagating in low-pressure helium DC discharges using the non-invasive laser-collision induced fluorescence (LCIF) diagnostic. To achieve this, we developed an analytic fit of collisional radiative model (CRM) predictions to interpret the LCIF data and recover quantitative two-dimensional spatial maps of the electron density, n_e, and the ratios of LCIF emission states that can be correlated with T_e with the use of accurate distribution functions at localized positions within striated helium discharges at 500 mTorr, 750 mTorr, and 1 Torr. To our knowledge, these are the first spatiotemporal, laser-based, experimental measurements of n_e in DC striations. The n_e and 447:588 ratio distributions align closely with striation theory. Constriction of the positive column appears to occur with decreased gas pressure, as shown by the radial n_e distribution. We identify a transition from a slow ionization wave to a fast ionization wave between 750 mTorr and 1 Torr. These experiments validate our analytic fit of n_e, allowing the implementation of an LCIF diagnostic in helium without the need to develop a CRM.

We present an experimental study on a terahertz quantum cascade laser (THz QCL) design that combines both two-well injector and direct-phonon scattering schemes, i.e., a so-called two-well injector direct-phonon design. As a result of the two-well injector direct-phonon scheme presented here, the lasers benefit from both a direct phonon scattering scheme for the lower laser level depopulation and a setback for the doping profile that reduces the overlap of the doped region with active laser states. Additionally, our design also has efficient isolation of the active laser levels from excited and continuum states as indicated by negative differential resistance behavior all the way up to room temperature. This scheme serves as a good platform for improving the temperature performance of THz QCLs as indicated by the encouraging temperature performance results of the device with a relatively high doping level of 7.56 × 10¹⁰ cm⁻² and T_max ∼ 167 K. With the right optimization of the molecular beam epitaxy growth and interface quality, the injection coupling strength, and the doping density and its profile, the device could potentially reach higher temperatures than the latest records reached for the maximum operating temperature (T_max) of THz QCLs.

Over the past few decades, inorganic nitride materials have grown in importance in part due to their potential as catalysts for the synthesis of NH3, a key ingredient in fertilizer and precursor to industrial chemicals. Of particular interest are the ternary (ABN) or higher-order nitrides with high metal-to-nitrogen ratios that show promise in enhancing NH3 synthesis reaction rates and yields via heterogeneous catalysis or chemical looping. Although metal nitrides are predicted to be numerous, the stability of nitrogen triple bonds found in N2, especially in comparison to the metal-nitrogen bonds, has considerably hindered synthetic efforts to produce complex nitride compounds. In this study, we present an exhaustive down-selection process to identify ternary nitrides for a promising chemical looping NH3 production mechanism. We also report on a facile and efficient two-step synthesis method that can produce well-characterized η-carbide Co3Mo3N/Fe3Mo3N or filled β-manganese Ni2Mo3N ternaries, as well as their associated quaternary, (Co,Fe)3Mo3N, (Fe,Ni)2Mo3N, and (Co,Ni)2Mo3N, solid solutions. To further explore the quaternary space, syntheses of (Co,Ni)3Mo3N (Ni ≤ 10 mol %) and Co3(Mo,W)3N (W ≤ 10 mol %) were also investigated. The structures of the nitrides were characterized via X-ray powder diffraction. The morphology and compositions were characterized with scanning electron microscopy. The multitude of chemically unique, but structurally related, nitrides suggests that properties such as nitrogen activity may be tunable, making the materials of great interest for NH3 synthesis schemes.

Cryosphere/Ocean Distributed Acoustic Sensing (CODAS) data collected from the Beaufort Sea, Alaska, using ~37.4 km of dark telecommunications fiber located at Oliktok Point, Alaska. Data were collected with a Silixa iDAS, using 10 m gauge length, 2 m spatial resolution, and 1000 Hz sample rate. Provided here are the DAS-recorded time series for the rapid refreeze event described in Baker & Abbott (2022) (see link below). This covers a date range of 2021-11-10 15:00 UTC to 2021-11-11 17:00 UTC. Data have been decimated to 100 Hz and 20 m (i.e., every 10th channel for 1831 channels, total), as used in Baker & Abbott (2022). Data have been extracted from raw format into 1-hour long .sac* files and organized into directories by channel number, spanning channels 100 to 18400. Time series units are nano-strainrate (nm/m/s). For distribution, data have been compressed into .zip files containing all time series files for 100 channels. *For information on the Seismic Analysis Code (SAC) file format: https://seiscode.iris.washington.edu/projects/sac

Short-ranged and line-gapped non-Hermitian Hamiltonians have strong topological invariants given by an index of an associated Fredholm operator. It is shown how these invariants can be accessed via the signature of a suitable spectral localizer. Here, this numerical technique is implemented in an example with relevance to the design of topological photonic systems, such as topological lasers.

Early-bird communication is a communication/computation overlap technique that combines fine-grained communication with partitioned communication to improve application run-time. Communication is divided among the compute threads such that each individual thread can initiate transmission of its portion of the data as soon as it is complete rather than waiting for all of the threads. However, the benefit of early-bird communication depends on the completion timing of the individual threads. In this paper, we measure and evaluate the potential overlap, the idle time each thread experiences between finishing their computation and the final thread finishing. These measurements help us understand whether a given application could benefit from early-bird communication. We present our technique for gathering this data and evaluate data collected from three proxy applications: MiniFE, MiniMD, and MiniQMC. To characterize the behavior of these workloads, we study the thread timings at both a macro level, i.e., across all threads across all runs of an application, and a micro level, i.e., within a single process of a single run. We observe that these applications exhibit significantly different behavior. While MiniFE and MiniQMC appear to be well-suited for early-bird communication because of their wider thread distribution and more frequent laggard threads, the behavior of MiniMD may limit its ability to leverage early-bird communication.

Traditional point-to-point communication sends data only after the entirety of the data is available. This includes situations where multiple actors (e.g., threads) contribute to the send buffer. As a result, cases where the completion times of these actors are widely distributed may be lost opportunities for optimization because data ready to be sent is waiting to be transmitted. Fine-grained communication exposes these opportunities by allowing buffers to be divided into elements that can then be sent independently (see e.g., Partitioned Communication in Message Passing Interface v4.0). While some research has been directed at exploring the utility of such 'early-bird' transmission, the overall search space for finding the best performing actor completion timings and element counts is large. In this work, we present an abstract model of fine-grained communication based on the LogGP model and a complementary benchmark. We use the model to explore actor completion timing scenarios and identify trends in communication behavior based on factors such as overall message size and delay between actor completions. We evaluate the benchmarks on three systems utilizing distinct network technologies and show that: (i) smaller numbers of elements are able to exploit most of the benefit of early-bird communication, (ii) performance benefit will depend non-trivially on application behavior, and (iii) benefits are highly network-dependent.

The plastic deformation of metals is a dissipative process. Some fraction of the plastic work is converted to heat which, given the temperature dependent response of metals, produces a thermal-mechanical coupling. In various cases, for instance when the loading is dynamic, this interaction can impact the resulting response of a material and/or system. Thus, appropriately capturing the heat generation from plastic work is necessary for various solid mechanics analysis. Determination of the fraction of work converted to heat has been long studied. Recent developments have demonstrated that the fraction is not constant but depends on various state variables. Resolving these features requires combined modeling and experimental studies. To this end, 304L stainless steel – a poor thermal conductor – was recently subjected to such an investigation. Advanced modeling capabilities were deployed to assess novel thermomechanically coupled experiments. As a complement to that study, in the current work a similar investigation is performed on copper – a good thermal conductor – to assess performance on the opposite end of the spectrum. The current document discusses these modeling efforts.

Building a stronger bioeconomy requires production capabilities that are largely generated through microbial genetic engineering. Plant feedstocks can additionally be genetically engineered to generate desirable feedstock traits and provide precursors for direct microbial conversion into desired products. The oleaginous yeast Rhodosporidium toruloides is a promising organism for this type of conversion as it can grow on a wide range of deconstructed biomass and consume a variety of carbon sources. Here, we leveraged R. toruloides native p-coumaric acid consumption pathway to accumulate protocatechuate (PCA) from 4-hydroxybenzoate (4HBA) released from a sorghum feedstock line genetically engineered to overproduce 4HBA. We did so by generating and evaluating an R. toruloides strain that accumulates PCA, RSΔ12623. We then show that at two scales a cholinium lysinate pretreatment with enzymatic saccharification successfully extracts 95% of the 4HBA from the engineered sorghum biomass while producing deconstructed lignin that can be more efficiently depolymerized in a subsequent thermochemical reaction. We also demonstrate that strain RSΔ12623 can convert more than 95% of 4HBA to PCA while consuming >95% of the glucose and >80% of the xylose present in sorghum hydrolysates. Finally, to evaluate the scalability of such fermentations, we conducted the conversion of 4HBA to PCA in a 2 L bioreactor under controlled conditions. This work demonstrates the potential of purposefully producing aromatic precursors in planta that can be liberated during biomass deconstruction for direct microbial conversion to desirable bioproducts.

Metamaterial resonators have become an efficient and versatile platform in the terahertz frequency range, finding applications in integrated optical devices, such as active modulators and detectors, and in fundamental research, e.g., ultrastrong light–matter investigations. Despite their growing use, characterization of modes supported by these subwavelength elements has proven to be challenging and it still relies on indirect observation of the collective far-field transmission/reflection properties of resonator arrays. Here, we present a broadband time-domain spectroscopic investigation of individual metamaterial resonators via a THz aperture scanning near-field microscope (a-SNOM). The time-domain a-SNOM allows the mapping and quantitative analysis of strongly confined modes supported by the resonators. In particular, a cross-polarized configuration presented here allows an investigation of weakly radiative modes. These results hold great potential to advance future metamaterial-based optoelectronic platforms for fundamental research in THz photonics.

UV photofragment spectroscopy and IR-UV double resonance methods are used to determine the structure and spectroscopic responses of a three-dimensional [2.2.2]-benzocryptand cage to the incorporation of a single K+ or Ba2+ imbedded inside it (labeled as K+-BzCrypt, Ba2+-BzCrypt). We studied the isolated ion-cryptand complex under cryo-cooled conditions, brought into the gas phase by nano-electrospray ionization. Incorporation of a phenyl ring in place of the central ethyl group in one of the three N-CH2-CH2-O-CH2-CH2-O-CH2-CH2-N chains provides a UV chromophore whose S0-S1 transition we probe. K+-BzCrypt and Ba2+-BzCrypt have their S0-S1 origin transitions at 35,925 and 36,446 cm-1, respectively, blue-shifted by 174 and 695 cm-1 from that of 1,2-dimethoxybenzene. These origins are used to excite a single conformation of each complex selectively and record their IR spectra using IR-UV dip spectroscopy. The alkyl CH stretch region (2800-3000 cm-1) is surprisingly sensitive to the presence and nature of the encapsulated ion. We carried out an exhaustive conformational search of cage conformations for K+-BzCrypt and Ba2+-BzCrypt, identifying two conformations (A and B) that lie below all others in energy. We extend our local mode anharmonic model of the CH stretch region to these strongly bound ion-cage complexes to predict conformation-specific alkyl CH stretch spectra, obtaining quantitative agreement with experiment for conformer A, the gas-phase global minimum. The large electrostatic effect of the charge on the O- and N-lone pairs affects the local mode frequencies of the CH2 groups adjacent to these atoms. The localized CH2 scissors modes are pushed up in frequency by the adjacent O/N-atoms so that their overtones have little effect on the alkyl CH stretch region. However, the localized CH2 wags are nearly degenerate and strongly coupled to one another, producing an array of delocalized wag normal modes, whose highest frequency members reach up above 1400 cm-1. As such, their overtones mix significantly with the CH stretch modes, most notably involving the CH2 symmetric stretch fundamentals of the central ethyl groups in the all-alkyl chains and the CH stretches adjacent to the N-atoms and antiperiplanar to the nitrogen lone pair.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

When high-energy-density materials are subjected to thermal or mechanical insults at extreme conditions (shock loading), a coupled response between the thermo-mechanical and chemical behaviors is systematically induced. Herein we develop a reaction model for the fast chemistry of 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) at the mesoscopic scale, where the chemical behavior is determined by underlying microscopic reactive simulations. The slow carbon cluster formation is not discussed in the present work. All-atom reactive molecular dynamics (MD) simulations are performed with the ReaxFF potential, and a reduced-order chemical kinetics model for TATB is fitted to isothermal and adiabatic simulations of single crystal chemical decomposition. Unsupervised machine learning techniques based on non-negative matrix factorization are applied to MD trajectories to model the decomposition kinetics of TATB in terms of a four-component model. The associated heats of reaction are fit to the temperature evolution from adiabatic decomposition trajectories. Using a chemical species analysis, we show that non-negative matrix factorization captures the main chemical decomposition steps of TATB and provides an accurate estimation of their evolution with temperature. The final analytical formulation, coupled to a diffusion term, is incorporated into a continuum formalism, and simulation results are compared one-to-one against MD simulations of 1D reaction propagation along different crystallographic directions and with different initial temperatures. A good agreement is found for both the temporal and spatial evolution of the temperature field.

Abstract not provided.

Abstract not provided.

Thermal-Hydrologic (TH) modeling of DECOVALEX 2023, Task C has continued in FY23. This report summarizes progress in TH modeling of Step 1c, with calibration modeling and the addition of shotcrete. The work involves 3-D modeling of the full-scale emplacement experiment at the Mont Terri Underground Rock Laboratory (Nagra, 2019). While Step 1 is focused on modeling the heating phase of the FE experiment with changes in pore pressure in the Opalinus clay resulting from heating, Step 1c is focused on calibration of models using available data.

Abstract not provided.

An x-ray imaging scheme using spherically bent crystals was implemented on the Z-machine to image x rays emitted by the hot, dense plasma generated by a Magnetized Liner Inertial Fusion (MagLIF) target. This diagnostic relies on a spherically bent crystal to capture x-ray emission over a narrow spectral range (<15 eV), which is established by a limiting aperture placed on the Rowland circle. The spherical crystal optic provides the necessary high-throughput and large field-of-view required to produce a bright image over the entire, one-cm length of the emitting column of a plasma. The average spatial resolution was measured and determined to be 18 µm for the highest resolution configuration. With this resolution, the radial size of the stagnation column can be accurately determined and radial structures, such as bifurcations in the column, are clearly resolved. The success of the spherical-crystal imager has motivated the implementation of a new, two-crystal configuration for identifying sources of spectral line emission using a differential imaging technique.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract: Advantages of the 2.5D HI (Heterogeneous Integration) electronics packaging of the power electronics compared to PCB packaging will be presented. Current 2.5D packaging effort using TSV (Through Silicon Via) will be presented in terms of fabrication, microstructural analysis, reliability, and thermal simulation.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Additive manufacturing (AM) is a relatively new technological advancement that allows for rapid prototyping, development of intricate shapes, and reduction in manufacturing time. The materials of interest for this project are Ultem 1010, ABS M30, FDM Nylon 12, PC, and PPSF. However, little is known regarding the aging behavior of these AM materials. The limited aging study outlined herein was designed to compare the chemical, physical, and mechanical properties of AM parts as they experience accelerated aging at 70 °C for a total of 24 weeks. In general, ABS M30 stood out as it appeared to undergo chemical and physical changes leading to increase in density and an overall more brittle material, making this commonly used material not attractive for long-term use.

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.

A robust data-management infrastructure is a key enabler for National Security Enterprise (NSE) capabilities in artificial intelligence and machine learning. This document describes efforts from a team of researchers at Sandia National Laboratories, Los Alamos National Laboratory, and Livermore National Laboratory to complete ASC Level II milestone #8854 “Assessment of Data-Management Infrastructure Needs for Production use of Advanced Machine learning and Artificial Intelligence.”

Abstract not provided.

This work shows that the static random access memory (SRAM) error rate for a commercial 65-nm device in a dose rate environment can be highly dependent upon the integrated dose (dose rate × pulse duration). While the typical metric for such testing is dose rate upset (DRU) level in rad(Si)/s, a series of dose rate experiments at Little Mountain Test Facility (LMTF) shows dependence on the integrated dose. The error rate is also found to be dependent on the core voltage, and the preradiation value of the bits. We believe that these effects are explained by a well charge depletion caused by gamma ray photocurrent.

Abstract not provided.

Abstract not provided.

Abstract not provided.

IT Service Management (ITSM) is an inimitable, ever-changing practice that is critical to all far-reaching organizations, including Sandia National Laboratories. With recent developments in technology and society – spurred by major events like the advent of ChatGPT and the pandemic – it is as important as ever to consider the implications and opportunities for ITSM. This report strategically synthesizes existing information about ITSM structures and examines the current state of Sandia’s service management entity: CCHD. Then, it looks at the emerging state of the world, culminating in CCHD-tailored recommendations for continual service improvement (CSI). Ultimately, the biggest matters to address are staff tenure, ticket documentation, and self-service facilitation. The next step might be to introduce AI to service channels in a non-system-invasive manner, namely a chatbot on the main CCHD page. All would serve to enhance end-user experiences, and by proxy Sandia’s output.

This special section of the Bulletin of the Seismological Society of America provides a broad overview on recent advances to the understanding of the seismoacoustic wavefield through 19 articles. Leveraging multiphenomenology datasets is instrumental for the continued success of future planetary missions, nuclear test ban treaty verification, and natural hazard monitoring. Progress in our theoretical understanding of mechanical coupling, advancements in coupled-media wave modeling, and developments of efficient multitechnology inversion procedures are key to fully exploiting geophysical datasets on Earth and beyond. We begin by highlighting papers describing experimental setups and instrumentation, followed by characterization of natural and anthropogenic sources of interest, and ending in new open-access datasets. Finally, we conclude with an overview of challenges that remain as well as some potential directions for future investigation within the growing multidisciplinary field of seismoacoustics.

A Machine and Deep Learning (MLDL) methodology is developed and applied to give a high fidelity, fast surrogate for 2D resistive MagnetoHydroDynamic (MHD) simulations of Magnetic Liner Inertial Fusion (MagLIF) implosions. The resistive MHD code GORGON is used to generate an ensemble of implosions with different liner aspect ratios, initial gas preheat temperatures (that is, different adiabats), and different liner perturbations. The liner density and magnetic field as functions of x, y, and z were generated. The Mallat Scattering Transformation (MST) is taken of the logarithm of both fields and a Principal Components Analysis (PCA) is done on the logarithm of the MST of both fields. The fields are projected onto the PCA vectors and a small number of these PCA vector components are kept. Singular Value Decompositions of the cross correlation of the input parameters to the output logarithm of the MST of the fields, and of the cross correlation of the SVD vector components to the PCA vector components are done. This allows the identification of the PCA vectors vis-a-vis the input parameters. Finally, a Multi Layer Perceptron (MLP) neural network with ReLU activation and a simple three layer encoder/decoder architecture is trained on this dataset to predict the PCA vector components of the fields as a function of time. Details of the implosion, stagnation, and the disassembly are well captured. Examination of the PCA vectors and a permutation importance analysis of the MLP show definitive evidence of an inverse turbulent cascade into a dipole emergent behavior. The orientation of the dipole is set by the initial liner perturbation. The analysis is repeated with a version of the MST which includes phase, called Wavelet Phase Harmonics (WPH). While WPH do not give the physical insight of the MST, they can and are inverted to give field configurations as a function of time, including field-to-field correlations.

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.

Photocathodes based on GaAs and other III-V semiconductors are capable of producing highly spin-polarized electron beams. GaAs/GaAsP superlattice photocathodes exhibit high spin polarization; however, the quantum efficiency (QE) is limited to 1% or less. To increase the QE, we fabricated a GaAs/GaAsP superlattice photocathode with a Distributed Bragg Reflector (DBR) underneath. This configuration creates a Fabry-Pérot cavity between the DBR and GaAs surface, which enhances the absorption of incident light and, consequently, the QE. These photocathode structures were grown using molecular beam epitaxy and achieved record quantum efficiencies exceeding 15% and electron spin polarization of about 75% when illuminated with near-bandgap photon energies.

Abstract not provided.

Abstract not provided.

The Xyce XDM Netlist Translator is a software solution developed by Sandia National Laboratories for translating SPICE netlists into a format readable by Xyce, Sandia’s internally developed SPICE-compatible simulation tool.

Several sources of interest often generate both low-frequency acoustic and seismic signals due to energy propagation through the atmosphere and the solid Earth. Seismic and acoustic observations are associated with a wide range of sources, including earthquakes, volcanoes, bolides, chemical and nuclear explosions, ocean noise, and others. The fusion of seismic and acoustic observations contributes to a better understanding of the source, both in terms of constraining source location and physics, as well as the seismic to acoustic coupling of energy. In this review, we summarize progress in seismoacoustic data processing, including recent developments in open-source data availability, low-cost seismic and acoustic sensors, and large-scale deployments of collocated sensors from 2010 to 2022. Similarly, we outline the recent advancements in modeling efforts for both source characteristics and propagation dynamics. Finally, we highlight the advantages of fusing multiphenomenological signals, focusing on current and future techniques to improve source detection, localization, and characterization efforts. This review aims to serve as a reference for seismologists, acousticians, and others within the growing field of seismoacoustics and multiphenomenology research.

Abstract not provided.

Abstract not provided.

The table presented below suggests the basic information that should be covered in a facility NMAC Plan for an NMAC program that is designed for nuclear security. The topics are appropriate for and should be addressed by all facilities in their NMAC Plans. They are appropriate for NMAC Plans for nuclear power plants, research reactors, fuel manufacturing facilities, facilities that produce medical isotopes, and other facilities. The difference is in the intensity with which the various measures are applied and the thoroughness of the description of the application (i.e., the program requirements). The robustness of a facility NMAC program and the content of its NMAC Plan should be graded in accordance with the type of facility and the category of its nuclear material.

Abstract not provided.

Abstract not provided.

Liquefied petroleum gas (LPG) is used in heating, cooking, and as a vehicle fuel (called autogas). A safety risk assessment may be needed to assess potential hazard scenarios and inform the regulations, codes, and standards that apply to LPG facilities, such as autogas refueling facilities. The frequency of unintended releases in an LPG system is an important aspect of a system quantitative risk assessment. This report documents estimation of leakage frequencies for individual components of LPG systems. These frequencies are described using uncertainty distributions obtained with Bayesian statistical methods, generic data, and LPG data which were publicly available. There was a lack of LPG data in the literature, so frequencies for most components were developed with generic data and should be used cautiously; without additional information about component leak frequencies in LPG systems, it is not known whether these generic frequencies may be conservative or non-conservative.

Abstract not provided.

Dominic D’Onofrio is currently a Junior studying Information and Technology at New Mexico Institute of Mining and Technology. He recently secured an Internship with NMCCoE, where he is involved with the TracerFIRE 12 project. Additionally, he is contributing to the load and security testing team by researching ways to implement pipelining and DevSecOps; this is his main project while he is at part time capacity for TracerFIRE 12. He is doing these projects to enhance his knowledge as a system administrator and gain a deeper understating of cybersecurity practices within national labs.

Abstract not provided.

The growth of helium bubbles impacts structural integrity of materials in nuclear applications. Understanding helium bubble nucleation and growth mechanisms is critical for improved material applications and aging predictions. Systematic molecular dynamics simulations have been performed to study helium bubble nucleation and growth mechanisms in Fe70Ni11Cr19 stainless steels. First, helium cluster diffusivities are calculated at a variety of helium cluster sizes and temperatures for systems with and without dislocations. Second, the process of diffusion of helium atoms to join existing helium bubbles is not deterministic and is hence studied using ensemble simulations for systems with and without vacancies, interstitials, and dislocations. We find that bubble nucleation depends on diffusion of not only single helium atoms, but also small helium clusters. Defects such as vacancies and dislocations can significantly impact the diffusion kinetics due to the trapping effects. Vacancies always increase the time for helium atoms to join existing bubbles due to the short-range trapping effect. This promotes bubble nucleation as opposed to bubble growth. Interestingly, dislocations can create a long-range trapping effect that reduces the time for helium atoms to join existing bubbles. This can promote bubble growth within a certain region near dislocations.

In this paper, we highlight how computational properties of biological dendrites can be leveraged for neuromorphic applications. Specifically, we demonstrate analog silicon dendrites that support multiplication mediated by conductance-based input in an interception model inspired by the biological dragonfly. We also demonstrate spatiotemporal pattern recognition and direction selectivity using dendrites on the Loihi neuromorphic platform. These dendritic circuits can be assembled hierarchically as building blocks for classifying complex spatiotemporal patterns.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Reversible logic schemes using flux solitons (fluxons) on long Josephson junctions (LJJs) have recently been proposed. The attraction of the fluxon is that it propagates ballistically along an LJJ until it encounters a change in the character of the LJJ, often a designed circuit element. Logic gates involve fluxons interacting with circuit elements and with other fluxons. However, testing of ballistic fluxon circuits requires other circuits outside the logic family to direct and control fluxon motion. We discuss two such non-reversible fluxon control circuits. First, the polarity filter gate is a simple non-reversible gate that allows one polarity of fluxon to pass, while reflecting the other polarity. In the off state both polarities reflect. Second, the polarity separator generalizes on the polarity filter concept and allows separation of the two fluxon polarities into different LJJs. We discuss simulations of these structures and possible applications.

Abstract not provided.

When designing measures to control infectious disease spread, it is crucial to understand the structure of the population for which interventions are being implemented. Recent work has highlighted the need for models that incorporate demographic heterogeneity not just in age structure but also by socioeconomic status (SES). Appropriately capturing additional sources of population heterogeneity requires considerable data and model development. To understand the potential disagreement between SES-explicit or SES-agnostic disease models, we adapted Sandia’s Adaptive Recovery Model (ARM) model to consider differences in contact structure and mortality by Social Vulnerability Index (SVI) on a theoretical network. We also incorporated an Average network that did not consider SVI. By exploring disparities in vaccine and PPE uptake by SES and comparing to Average networks, as well as analyzing the influence of global vs. local contact, we found that the two model constructions often predicted different outcomes. Whether these differences are truly reflective of incorporating SES, and which model most closely represents reality, merits further investigation.

There is growing interest to extend low-rank matrix decompositions to multi-way arrays, or tensors. One fundamental low-rank tensor decomposition is the canonical polyadic decomposition (CPD). The challenge of fitting a low-rank, nonnegative CPD model to Poisson-distributed count data is of particular interest. Several popular algorithms use local search methods to approximate the global maximum likelihood estimator from local minima. Simultaneously, a recent trend in theoretical computer science and numerical linear algebra leverages randomization to solve very large, hard problems. The typical approach is to use randomization for a fast approximation and determinism for refinement to yield effective algorithms with theoretical guarantees. Two popular algorithms for Poisson CPD reflect that emergent dichotomy: CP Alternating Poisson Regression is a deterministic algorithm and Generalized Canonical Polyadic decomposition makes use of stochastic algorithms in several variants. This work extends recent work to develop two new methods that leverage randomized and deterministic algorithms for improved accuracy and performance.

This report summarizes the fiscal year 2023 (FY23) status of the second phase of a series of borehole heater tests in salt at the Waste Isolation Pilot Plant (WIPP) funded by the Disposal Research and Development (R&D) program of the Spent Fuel & Waste Science and Technology (SFWST) office at the US Department of Energy’s Office of Nuclear Energy’s (DOE-NE) Office in the Spent Fuel and Waste Disposition (SFWD) program.

Nuclear power plants (NPPs) are considering flexible plant operations to take advantage of excess thermal and electrical energy. One option for NPPs is to pursue hydrogen production through high temperature electrolysis as an alternate revenue stream to remain economically viable. The intent of this study is to investigate the risk of a hydrogen production facility in close proximity to an NPP. A 100 MW, 500 MW, and 1,000 MW facility are evaluated herein. Previous analyses have evaluated preliminary designs of a hydrogen production facility in a conservative manner to determine if it is feasible to co-locate the facility within 1 km of an NPP. This analysis specifically evaluates the risk components of different hydrogen production facility designs, including the likelihood of a leak within the system and the associated consequence to critical NPP targets. This analysis shows that although the likelihood of a leak in an HTEF is not negligible, the consequence to critical NPP targets is not expected to lead to a failure given adequate distance from the plant.

Factor Fixing (FF) is a common method for reducing the number of model parameters to lower computational cost. FF typically starts with distinguishing the insensitive parameters from the sensitive and pursues uncertainty quantification (UQ) on the resulting reduced-order model, fixing each insensitive parameter at a fixed value. There is a need, however, to expand such a common approach to consider the effects of decision choices in the FF-UQ procedure on metrics of interest. Therefore, to guide the use of FF and increase confidence in the resulting dimension-reduced model, we propose a new adaptive framework consisting of four principles: (a) re-parameterize the model first to reduce obvious non-identifiable parameter combinations, (b) focus on decision relevance especially with respect to errors in quantities of interest (QoI), (c) conduct adaptive evaluation and robustness assessment of errors in the QoI across FF choices as sample size increases, and (d) reconsider whether fixing is warranted. The framework is demonstrated on a spatially-distributed water quality model. The error in estimates of QoI caused by FF can be estimated using a Polynomial Chaos Expansion (PCE) surrogate model. Built with 70 model runs, the surrogate is computationally inexpensive to evaluate and can provide global sensitivity indices for free. For the selected catchment, just two factors may provide an acceptably accurate estimate of model uncertainty in the average annual load of Total Suspended Solids (TSS), suggesting that reducing the uncertainty in these two parameters is a priority for future work before undertaking further formal uncertainty quantification.

Earthquakes have repeatedly been shown to produce inaudible acoustic signals (< 20 Hz), otherwise known as infrasound. These signals can propagate hundreds to thousands of kilometers and still be detected by ground-based infrasound arrays depending on the source strength, distance between source and receiver, and atmospheric conditions. Another type of signal arrival at infrasound arrays is the seismic induced motion of the sensor itself, or ground-motion-induced sensor noise. Measured acoustic and seismic waves produced by earthquakes can provide insight into properties of the earthquake such as magnitude, depth, and focal mechanism, as well as information about the local lithology and atmospheric conditions. Large earthquakes that produce strong acoustic signals detected at distances greater than 100 km are the most commonly studied; however, more recent studies have found that smaller magnitude earthquakes (Mw < 2:0) can be detected at short ranges. In that vein, this study will investigate the ability for a long-term deployment of infrasound sensors (deployed as part of the Source Physics Experiments [SPE] from 2014 to 2020) to detect both seismic and infrasonic signals from earthquakes at local ranges (< 50 km). Methods used include a combination of spectral analysis and automated array processing, supported by U.S. Geological Survey earthquake bulletins. This investigation revealed no clear acoustic detections for short range earthquakes. However, secondary infrasound from an Mw 7.1 earthquake over 200 km away was detected. Important insights were also made regarding the performance of the SPE networks including detections of other acoustic sources such as bolides and rocket launches. Finally, evaluation of the infrasound arrays is performed to provide insight into optimal deployments for targeting earthquake infrasound.

Abstract not provided.

This presentation introduces the use of a novel methodology for analyzing emerging technology development through the case study of atomic clocks.

Commercially available thermal interface materials are investigated in this study. First, eight different samples were thermal cycled for 30 to 201 days from ambient to 125°C or ambient to 175°C in an accelerated aging test. Compared to the new samples, the aged samples exhibited varying optical and chemical changes. The most deteriorated samples were eliminated from the study and the remaining five samples underwent further analysis. A transient plane source measurement was conducted on four of the five remaining samples to verify their thermal conductivities, an off-gassing study was performed to detect volatile compounds, and a salt fog corrosion test was conducted to explore corrosion effects. The silicone pads and gels with diboron trioxide or aluminum oxide filler materials exhibited the best overall performance with the least amount of visual and spectroscopic degradation after thermal cycling. While the offgassing results were inconclusive, the silicone-based materials released approximately 1000 to 25000 ppb of volatile organic compounds, with Fujipoly’s GR130A gel performing the best, and Laird’s TFLEX760 and Parker Chomerics Therm-A-Gap Gel 30 releasing the most volatile compounds. None of the downselected samples exhibited corrosion effects in salt fog.

The classical Drude model provides an accurate description of the plasma resonance of three-dimensional materials, but only partially explains two-dimensional systems where quantum mechanical effects dominate such as P:δ layers - atomically thin sheets of phosphorus dopants in silicon that induce electronic properties beyond traditional doping. Previously it was shown that P:δ layers produce a distinct Drude tail feature in ellipsometry measurements. However, the ellipsometric spectra could not be properly fit by modeling the δ layer as a discrete layer of classical Drude metal. In particular, even for large broadening corresponding to extremely short relaxation times, a plasma resonance feature was anticipated but not evident in the experimental data. In this work, we develop a physically accurate description of this system, which reveals a general approach to designing thin films with intentionally suppressed plasma resonances. Our model takes into account the strong charge-density confinement and resulting quantum mechanical description of a P:δ layer. We show that the absence of a plasma resonance feature results from a combination of two factors: (i) the sharply varying charge-density profile due to strong confinement in the direction of growth; and (ii) the effective mass and relaxation time anisotropy due to valley degeneracy. The plasma resonance reappears when the atoms composing the δ layer are allowed to diffuse out from the plane of the layer, destroying its well-confined two-dimensional character that is critical to its distinctive electronic properties.

Abstract not provided.

Abstract not provided.

It is essential to Sandia National Laboratory’s continued success in scientific and technological advances and mission delivery to embrace a hybrid workforce culture under which current and future employees can thrive. This report focuses on the findings of the Hybrid Work Team for the Center for Computing Research, which met weekly from March to June 2023 and conducted a survey across the Center at Sandia. Conclusions in this report are drawn from the 9 authors of this report, which comprises the Hybrid Work Team, and 15 responses to a center-wide survey, as well as numerous conversations with colleagues. A major finding was widespread dissatisfaction with the quantity, execution, and tooling surrounding formal meetings with remote participants. While there was consensus that remote work enables people to produce high quality individual and technical work, there was also consensus that there was widespread social disconnect, with particular concern about hires that were made after the onset of the Covid-19 pandemic. There were many concerns about tooling and policy to facilitate remote collaboration both within Sandia and with its external collaborators. This report includes recommendations for mitigating these problems. For problems for which obvious recommendations cannot be made, ideas of what a successful solution might look like are presented.

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.

An aircraft commander needs to be aware of weather phenomena that might be hazardous to his aircraft and mission. An important tool for this is airborne weather (WX) detection radar. The airborne WX radar needs to map weather for the aircraft commander that might be relevant to the safety of the aircraft, which involves both detecting a weather phenomenon, and to some extent seeing through it to detect weather phenomena behind it. Many factors influence the performance of an airborne WX radar

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

We report flow statistics and visualizations from molecular-gas-dynamics simulations using the direct simulation Monte Carlo (DSMC) method for turbulent Couette flow in a minimal domain where the lower wall is replaced by an idealized permeable fibrous substrate representative of thermal-protection-system materials for which the Knudsen number is O(10-1). Comparisons are made with smooth-wall DSMC simulations and smooth-wall direct numerical simulations (DNS) of the Navier-Stokes equations for the same conditions. Roughness, permeability, and noncontinuum effects are assessed. In the range of Reynolds numbers considered herein, the scalings of the skin friction on the permeable substrate and of the mean flow within the substrate suggest that they are dominated by viscous effects. While the regenerative cycle characteristic of smooth-wall turbulence remains intact for all cases considered, we observe that the near-wall velocity fluctuations are modulated by the permeable substrate with a wavelength equal to the pore spacing. Additionally, the flow within the substrate shows significant rarefaction effects, resulting in an apparent permeability that is 13% larger than the intrinsic permeability. In contrast, the smooth-wall DSMC and DNS simulations exhibit remarkably good agreement for the statistics examined, despite the Knudsen number based on the viscous length scale being as large as O(10-1). This latter result is at variance with classical estimates for the breakdown of the continuum assumption and calls for further investigations into the interaction of noncontinuum effects and turbulence.

Waveform cross-correlation is a sensitive phase-matched filtering technique that can detect seismic events for nuclear explosion monitoring. However, there are outstanding challenges with correlation detectors, most notably a direct dependence on the completeness of the waveform template library. To ameliorate these challenges, we investigate how dynamic time warping (DTW) may make waveform correlation more robust. DTW analyzes the differences between two time series and attempts to “warp” one time series relative to another in a recursive manner. We apply DTW to synthetic earthquake and recorded explosion templates to expand the capability of correlation detectors. We explore what conditions (e.g., source, station distance, frequency bands) and/or DTW algorithms generate stronger correlation scores. We show that DTW performs well on noisy signals and can dramatically improve the cross-correlation coefficient between a template and data-stream waveform. We conclude with recommendations on how to utilize DTW in nuclear monitoring detection.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

An additive manufacturing approach combining aerosol jet printing (AJP) and electrodeposition opens a new pathway to the production of lightweight coreless flyback transformer devices for power electronics. AJP of seed layers with resolution on the order of 30μm is combined with electrodeposition of Cu and Ni for decreased resistance. This combined approach addresses known shortcomings of AJP and electrodeposition. Nanoparticle inks used in AJP of metals have low conductivity versus bulk materials due to their high grain boundary resistance. There is a lack of readily available high-resolution patterning techniques for electrodeposition outside of expensive clean-room-based lithography techniques. Combining these two techniques enables the patterning of high-resolution, high-conductivity components. In this manuscript, we report on the construction of coreless flyback transformers consisting of two-layer primary and two-layer secondary spiral inductors separated by layers of a printed UV-curable dielectric. An input voltage of 17 V at 400 kHz was amplified to an output of 1250 V corresponding to a gain of 73.5. COMSOL modeling at the individual inductor level and at the transformer level was used to compare expected inductance, equivalent series resistance (ESR), and coupling with experimentally measured values.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The goal of this report is to provide insight to the development of vt-tv, a C++ HPC visualization tool designed for insightful analysis of load-balancing metrics in the DARMA toolkit. In particular, it delves into its modular data model and diverse usage scenarios, emphasizing adaptability and efficiency.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Thrifty Array Foram (TAF) files store numeric data in a binary format, minimizing storage requirements while preserving quick read access. Real data of any size and dimensionality can be stored in this format at varying degrees of numeric precision. Implicit array are associated with each dimension, eliminating the need to explicitly store uniformly-spaced grid vectors. Unlimited text comments may be included with the array for user documentation, and every file begins with a text synopsis of the binary structure. The format is deliberately designed for memory mapping, where portions of the array can be read without loading the entire file at once.

Abstract not provided.

We present a graph neural network approach that fully automates the prediction of defect formation enthalpies for any crystallographic site from the ideal crystal structure, without the need to create defected atomic structure models as input. Here we used density functional theory reference data for vacancy defects in oxides, to train a defect graph neural network (dGNN) model that replaces the density functional theory supercell relaxations otherwise required for each symmetrically unique crystal site. Interfaced with thermodynamic calculations of reduction entropies and associated free energies, the dGNN model is applied to the screening of oxides in the Materials Project database, connecting the zero-kelvin defect enthalpies to high-temperature process conditions relevant for solar thermochemical hydrogen production and other energy applications. The dGNN approach is applicable to arbitrary structures with an accuracy limited principally by the amount and diversity of the training data, and it is generalizable to other defect types and advanced graph convolution architectures. It will help to tackle future materials discovery problems in clean energy and beyond.

Abstract not provided.

Machine learning methodologies can provide insight into Brønsted-Guggenheim-Scatchard specific ion interaction theory (SIT) parameter values where experimental data availability may be limited. This study develops and executes machine learning frameworks to model the SIT interaction coefficient, ε. Key findings include successful estimations of ε via artificial neural networks using clustering and value prediction approaches. Applicability to other chemical parameters is also assessed briefly. Models developed here provide support for a use-case of machine learning in geologic nuclear waste disposal research applications, namely in predictions of chemical behaviors of high ionic strength solutions (i.e., subsurface brines).

Abstract not provided.

Abstract not provided.

Abstract not provided.

In many areas of constrained optimization, representing all possible constraints that give rise to an accurate feasible region can be difficult and computationally prohibitive for online use. Satisfying feasibility constraints becomes more challenging in high-dimensional, non-convex regimes which are common in engineering applications. A prominent example that is explored in the manuscript is the security-constrained optimal power flow (SCOPF) problem, which minimizes power generation costs, while enforcing system feasibility under contingency failures in the transmission network. In its full form, this problem has been modeled as a nonlinear two-stage stochastic programming problem. In this work, we propose a hybrid structure that incorporates and takes advantage of both a high-fidelity physical model and fast machine learning surrogates. Neural network (NN) models have been shown to classify highly non-linear functions and can be trained offline but require large training sets. In this work, we present how model-guided sampling can efficiently create datasets that are highly informative to a NN classifier for non-convex functions. We show how the resultant NN surrogates can be integrated into a non-linear program as smooth, continuous functions to simultaneously optimize the objective function and enforce feasibility using existing non-linear solvers. Overall, this allows us to optimize instances of the SCOPF problem with an order of magnitude CPU improvement over existing methods.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Threaded fastener behavior can be an important aspect of complex component and system behavior, but there is no one-size-fits-all finite element analysis technique. Proper modeling of threaded fastener joints requires careful consideration of many details, from test setup and data acquisition to constitutive modeling and uncertainty quantification approaches. This report details analysis of a “mini-radax” bolted-joint exemplar where a Discrete-Direct uncertainty quantification approach is employed to evaluate margin of the component. The mini-radax geometry is tested to failure on a drop table, and single-coupon tests of individual fasteners serve as foundational data for the analysis. Analysis predictions complement the test data well and provide additional context for engineering decision-making.

This paper describes a summer internship project undertaken at Sandia National Labs (SNL), both current status and future work. The project was to explore various machine learning approaches for use on turbulent flow data. Specifically, unsupervised classification of turbulent flow data was explored. First, the usage of models in this field is discussed, and several issues in the common usage of the models are identified. Solutions to these issues are then proposed, in the form of a Bayesian filtering approach which probabilistically incorporates multiple sources of data to improve confidence in a result. Several types of sensors are suggested for this method, the incorporation of which range from semi-supervised learning approaches to fully unsupervised. These approaches are then tested on several turbulent flow cases.

Abstract not provided.

Abstract not provided.