Publications

Abstract not provided.

Abstract not provided.

Abstract not provided.

In this article, we present a general methodology to combine the Discontinuous PetrovGalerkin (DPG) method in space and time in the context of methods of lines for transient advection-reaction problems. We first introduce a semidiscretization in space with a DPG method redefining the ideas of optimal testing and practicality of the method in this context. Then, we apply the recently developed DPG-based time-marching scheme, which is of exponential-type, to the resulting system of Ordinary Differential Equations (ODEs). We also discuss how to efficiently compute the action of the exponential of the matrix coming from the space semidiscretization without assembling the full matrix. Finally, we verify the proposed method for 1D+time advection-reaction problems showing optimal convergence rates for smooth solutions and more stable results for linear conservation laws comparing to the classical exponential integrators.

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.

Hydrogen is an important resource for many different industries throughout the world, including refining, manufacturing, and as a direct energy source. Hydrogen production, through methods such as steam methane reforming, has been developed over several decades. There is a large global demand for hydrogen from these industries and safe production and distribution are paramount for hydrogen systems. Codes and standards have been developed to reduce the risk associated with hydrogen accidents to the public. These codes and standards are similar to those in other industries in which there is inherent risk to the public, such as gasoline and natural gas production and distribution. Although there will always be a risk to the public in these types of fuels, the codes and standards are developed to reduce the likelihood of an accident occurring and reduce the severity of impact, should one occur. This report reviews the current state of hydrogen in the United States and outlines the codes and standards that ensure safe operation of hydrogen systems. The total hydrogen demand and use in different industries is identified. Additionally, the current landscape of hydrogen production and fueling stations in the United States is outlined. The safety of hydrogen systems is discussed through an overview of the purpose, methods, and content included in codes and standards. As outlined in this safety overview, the risk to the public in operation of hydrogen generation facilities and fueling stations is reduced through implementation of appropriate measures. Codes, such as NFPA 2, ensure that the risk associated with a hydrogen system is no greater than the risk presented by gasoline refueling stations.

Abstract not provided.

Sandia National Laboratories (SNL) is a multimission laboratory located in Albuquerque, New Mexico, and is one of three National Nuclear Security Administration research and development laboratories located in the United States. Recently, SNL’s Emergency Response Team (ERT) responded to an incident involving a sulfur dioxide (SO2)-fixed monitor, setting off the alarm inside a laboratory and in the adjacent hallway. The potential sources for the alarm were various experiments involving batteries and an uninterrupted power supply (UPS) in the immediate area.

Abstract not provided.

Abstract not provided.

The increasing availability of advanced metering infrastructure (AMI) data has led to significant improvements in load modeling accuracy. However, since many AMI devices were installed to facilitate billing practices, few utilities record or store reactive power demand measurements from their AMI. When reactive power measurements are unavailable, simplifying assumptions are often applied for load modeling purposes, such as applying constant power factors to the loads. The objective of this work is to quantify the impact that reactive power load modeling practices can have on distribution system analysis, with a particular focus on evaluating the behaviors of distributed photovoltaic (PV) systems with advanced inverter capabilities. Quasi-static time-series simulations were conducted after applying a variety of reactive power load modeling approaches, and the results were compared to a baseline scenario in which real and reactive power measurements were available at all customer locations on the circuit. Overall, it was observed that applying constant power factors to loads can lead to significant errors when evaluating customer voltage profiles, but that performing per-phase time-series reactive power allocation can be utilized to reduce these errors by about 6x, on average, resulting in more accurate evaluations of advanced inverter functions.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Optical fiber diagnostics are extensively used in pulsed power experiments, such as the Sandia Z machine. However, radiation produced in a pulsed power environment can significantly affect these measurements. Catastrophic fiber darkening may be mitigated with shielding, but no flexible material can stop all radiation produced by the machine and/or target. Radiation-induced refractive index modulations are particularly challenging for optical interferometry. Several approaches for radiation-tolerant photonic Doppler velocimetry are discussed here.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

A proof-of-concept tool, the Produced Water-Economic, Socio, Environmental Simulation model (PW-ESESim), was developed to support ease of analysis. The tool was designed to facilitate head-to-head comparison of alternative produced water source, treatment, and reuse water management strategies. A graphical user interface (GUI) guides the user through the selection and design of alternative produced water treatment and reuse strategies and the associated health and safety risk and economic benefits. At the highest conceptual level, alternative water strategies include the selection of a source water (locally or regionally available produced water), treatment strategy (pre-treatment, physical, chemical, biological, desalination, and post-treatment processes) and product water purpose (e.g., irrigation, industrial processing, environmental). After selection of these details, the PW-ESESim output a number of key economic, societal, environmental, public/ecological health and safety metrics to support user decision-making; specific examples include, cost of treatment, improvements in freshwater availability, human and ecologic health impacts and growth in local jobs and the economy. Through the simulation of different produced water treatment and management strategies, tradeoffs are identified and used to inform fit-for-purpose produced water treatment and reuse management decisions. While the tool was initially designed using Southeastern New Mexico (Permian Basin) as a case study, the general design of the PW-ESESim model can be extended to support other oil and gas regions of the U.S.

Soil carbon can be divided into two categories: organic and inorganic. Soil inorganic carbon (SIC) is present in carbonate minerals in the soil and is often found in dry, arid regions. Examples of SIC include calcium carbonate (CaCO₃) and magnesium carbonate (MgCO₃), both of which play important roles in soil health. Soil organic carbon (SOC) is found in fresh plant matter (available SOC) and as humus or charcoal (inert SOC). Both types of carbon act as storage in the global carbon cycle. As a carbon sink, soil carbon has the potential to store carbon that would otherwise remain in the atmosphere as CO₂, one of the primary greenhouse emissions. As such, soil is under increasing attention and research to be used as a sequestration (i.e., isolation) method to reduce the amount of carbon in the atmosphere. This type of carbon sequestration is called biological sequestration. SOC typically stores carbon for several decades (depending on decomposition rates) while SIC can store carbon for more than 70,000 years. Common sequestration techniques for SOC usually fall under the category of land management: planting perennials, keeping plant residue and composting, reducing tilling, and other agricultural practices that vary by region. SIC sequestration through carbonates naturally takes thousands of years but there have been studies to increase SIC sequestration through the addition of silicates.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

This document is a reference guide to the Xyce Parallel Electronic Simulator, and is a companion document to the Xyce Users' Guide. The focus of this document is (to the extent possible) exhaustively list device parameters, solver options, parser options, and other usage details of Xyce. This document is not intended to be a tutorial. Users who are new to circuit simulation are better served by the Xyce Users' Guide.

Abstract not provided.

Niobium (Nb)-doped lead-tin-zirconate-titanate (PSZT) ceramics near the lead-zirconate-titanate 95/5 orthorhombic AFE-rhombohedral FE morphotropic phase boundary (PSZT 13.5/81/5.5 -1.6Nb) were prepared with up to 10 mol.% of hafnium (Hf) substituted for zirconium. The ceramics were prepared by a traditional solid-state synthesis route and sintered to near full density at 1150°C for 6 h in sealed alumina crucibles with self-same material as the lead vapor source. All compositions were ∼98% dense with no detectable secondary phases by X-ray diffraction. The grain size was ∼3 μm for all compositions, consisting of equiaxed grains with intergranular porosity. The compositions exhibited remnant polarization values of ∼32 μC/cm2. Depolarization by the hydrostatic pressure-induced FE-AFE phase transition occurred at 310 MPa for all compositions, resulting in a total depolarization output of 32.4 μC/cm2 for the PSZT ceramics. Evaluation of the R3c-R3m and R3m-Pm (Formula presented.) m phase transition temperatures by impedance spectroscopy showed temperatures on heating ranging from 86 to 92°C and 186 to 182°C, respectively, for increasing nominal Hf content. Thermal hysteresis of the phase transitions was also observed in the ceramics, with the transition temperature on cooling being 1–4°C lower. The study demonstrated that the PSZT ceramics are relatively insensitive to variations in Hf content in the range of 0 to 10 mol.%.

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.

This report begins with an evaluation of cladding degradation mechanisms deemed important to assessing barrier capability. Unlike similar efforts done in the past, this evaluation accounts for the hypothetical conditions associated with direct dual-purpose canister (DPC) disposal, including conditions resulting from a postulated, in-package, steady-state criticality event. A total of 16 cladding degradation mechanisms are examined assuming direct disposal of DPCs in two different hypothetical repositories: a saturated repository in shale and an unsaturated repository in alluvium.

This report contains a summary of our efforts to correlate head injury simulations predicting intracranial fluid cavitation with clinical assessments of brain injury from blunt impact to the head. Magnetic resonance imaging (MRI) data, collected on traumatic brain injury (TBI) subjects by researchers at the MIND Institute of New Mexico, was acquired for the current work. Specific blunt impact TBI case histories were selected from the TBI data for further study and possible correlation with simulation. Both group and single-subject case histories were examined. We found one single-subject case that was particularly suited for correlation with simulation. Diffusion tensor image (DTI) analysis of the TBI subject identified white matter regions within the brain displaying reductions in fractional anisotropy (FA), an indicator of local damage to the white matter axonal structures. Analysis of functional magnetic resonance image (fMRI) data collected on this individual identified localized regions of the brain displaying hypoactivity, another indicator of brain injury. We conducted high fidelity simulations of head impact experienced by the TBI subject using the Sandia head-neck-torso model and the shock physics computer code CTH. Intracranial fluid cavitation predictions were compared with maps of DTI fractional anisotropy and fMRI hypoactivity to assess whether a possible correlation exists. The ultimate goal of this work is to assess whether one can correlate simulation predictions of intracranial fluid cavitation with the brain injured sites identified by the fMRI and DTI analyses. The outcome of this effort is described in this report.

Consider a scalar field f : M → R, where M is a triangulated simplicial mesh in R^d. A level set, or contour, at value v is a connected component of f^–1 (v). As v is changed, these contours change topology, merge into each other, or split. Contour trees are concise representations of f that track this contour behavior. The vertices of these trees are the critical points of f, where the gradient is zero. The edges represent changes in the topology of contours. It is a fundamental data structure in data analysis and visualization, and there is significant previous work (both theoretical and practical) on algorithms for constructing contour trees. Suppose M has n vertices, N facets, and t critical points. A classic result of Carr, Snoeyink, and Axen (2000) gives an algorithm that takes O(n log n+Nα(N)) time (where α(·) is the inverse Ackermann function). A further improvement to O(t log t + N) time was given by Chiang et al. All these algorithms involve a global sort of the critical points, a significant computational bottleneck. Unfortunately, lower bounds of Ω(t log t) also exist. We present the first algorithm that can avoid the global sort and has a refined time complexity that depends on the contour tree structure. Intuitively, if the tree is short and fat, we get significant improvements in running time. For a partition of the contour tree into a set of descending paths, P, our algorithm runs in O($\Sigma$_pϵP |p| log |p| + tα(t) + N). This is at most O(t log D + N), where D is the diameter of the contour tree. Moreover, it is O(tα(t) + N) for balanced trees, a significant improvement over the previous complexity. Our algorithm requires numerous ideas: partitioning the contour tree into join and split trees, a local growing procedure to iteratively build contour trees, and the use of heavy path decompositions for the time complexity analysis. There is a crucial use of a family of binomial heaps to maintain priorities, ensuring that any comparison made is between comparable nodes of the contour tree. We also prove lower bounds showing that the $\Sigma$_pϵP |p| log |p| complexity is inherent to computing contour trees.

An existing shared risk framework designed for assessing and comparing threat-based risks to water utilities is being extended to incorporate electric power. An important differentiating characteristic of this framework is the use of a system-centric rather than an asset-centric approach. This approach allows anonymous sharing of results and enables comparison of assessments across different utilities within an infrastructure sector. By allowing utility owners to compare their assessments with others, they can improve their self-assessments and identification of "unknown unknowns". This document provides an approach for extension of the framework to electric power, including treatment of dependencies and interdependencies. The systems, threats, and mathematical description of associated risks used in a prototype framework are provided. The method is extensible so that additional infrastructure sectors can be incorporated. Preliminary results for a proof of concept calculation are provided.

Biological monitoring logs for Swainson's and Red-tailed Hawk's nesting behaviors during construction at SNL. Logs also include information on other observed wildlife. Three different monitoring logs are present and include images and buffer map.

Aeroseismometery is a novel, cutting edge capability that involves balloon based systems for detecting and geolocating sources of infrasound. The incident infrasound from a range of sources such as volcanos, earthquakes, explosions, supersonic aircraft impinges upon the balloon system causing it to respond dynamically. The dynamic response is post-processed to locate the infrasound source. This report documents the derivation of an analytical model that predicts the balloon dynamics. Governing equations for the system are derived as well as a transfer function relating the infrasound signal to the net force on the balloon components. Experimental measurements of the infrasound signals are convolved with the transfer function and the governing equations numerically time integrated to obtain predictions of the displacement, velocity and acceleration of the balloon system. The predictions are compared to the experimental measurements with good agreement observed. The derivation focuses only on the vertical dynamics of the balloon system. Future work will develop governing equations for the swinging response of the balloon to the incident infrasound.

QOOH radicals are key intermediates in the chain of reactions leading to the autoignition of hydrocarbons and oxygenated organic compounds. They are thought to undergo two main reactions: OH elimination to form a cyclic ether and HO2 elimination to form an alkene. However, theoretical analysis of various substituted hydroperoxyalkyl radicals has found two new pathways: OH transfer and internal H abstraction assisted OH elimination. To determine the importance of these new pathways, their barrier heights for several substituted alkanes were calculated using various quantum chemical theories and compared to those of the well-known pathways. Several cases revealed possible competition with the well-known pathways. Rate coefficients were calculated for propyl systems but further studies will need to complete rate coefficients and branching fractions for all systems analyzed to understand these new reactions’ role in autoignition.

In this paper, we propose a hybrid method that uses stochastic and deterministic search to compute the maximum likelihood estimator of a low-rank count tensor with Poisson loss via state-of-theart local methods. Our approach is inspired by Simulated Annealing for global optimization and allows for fine-grain parameter tuning as well as adaptive updates to algorithm parameters. We present numerical results that indicate our hybrid approach can compute better approximations to the maximum likelihood estimator with less computation than the state-of-the-art methods by themselves.

Drinking water has and will continue to be at the foundation of our nation’s well-being and there is a growing interest in United States (US) drinking water quality. Nearly 30% of the United States population obtained their water from community water systems that did not meet federal regulations in 2019. Given the heavy interactions between society and drinking water quality, this study integrates social constructionism, environmental injustice, and sociohydrological systems to evaluate local awareness of drinking water quality issues. By employing text analytics, we explore potential drivers of regional water quality narratives within 25 local news sources across the United States. Specifically, we assess the relationship between printed local newspapers and water quality violations in communities as well as the influence of social, political, and economic factors on the coverage of drinking water quality issues. Results suggest that the volume and/or frequency of local drinking water violations is not directly reflected in local news coverage. Additionally, news coverage varied across sociodemographic features, with a negative relationship between Hispanic populations and news coverage of Lead and Copper Rule, and a positive relationship among non-Hispanic white populations. These findings extend current understanding of variations in local narratives to consider nuances of water quality issues and indicate opportunities for increasing equity in environmental risk communication.

Electrode-scale heterogeneity can combine with complex electrochemical interactions to impede lithium-ion battery performance, particularly during fast charging. This study investigates the influence of electrode heterogeneity at different scales on the lithium-ion battery electrochemical performance under operational extremes. We employ image-based mesoscale simulation in conjunction with a three-dimensional electrochemical model to predict performance variability in 14 graphite electrode X-ray computed tomography data sets. Our analysis reveals that the tortuous anisotropy stemming from the variable particle morphology has a dominating influence on the overall cell performance. Cells with platelet morphology achieve lower capacity, higher heat generation rates, and severe plating under extreme fast charge conditions. On the contrary, the heterogeneity due to the active material clustering alone has minimal impact. Our work suggests that manufacturing electrodes with more homogeneous and isotropic particle morphology will improve electrochemical performance and improve safety, enabling electromobility.

An in situ counted ion implantation experiment improving the error on the number of ions required to form a single optically active silicon vacancy (SiV) defect in diamond 7-fold compared to timed implantation is presented. Traditional timed implantation relies on a beam current measurement followed by implantation with a preset pulse duration. It is dominated by Poisson statistics, resulting in large errors for low ion numbers. Instead, our in situ detection, measuring the ion number arriving at the substrate, results in a 2-fold improvement of the error on the ion number required to generate a single SiV compared to timed implantation. Through postimplantation analysis, the error is improved 7-fold compared to timed implantation. SiVs are detected by photoluminescence spectroscopy, and the yield of 2.98% is calculated through the photoluminescence count rate. Hanbury-Brown-Twiss interferometry is performed on locations potentially hosting single-photon emitters, confirming that 82% of the locations exhibit single photon emission statistics.

Nuclear magnetic resonance (NMR) and nuclear quadrupole resonance (NQR) spectroscopy of bulk quantum materials have provided insight into phenomena, such as quantum phase criticality, magnetism, and superconductivity. With the emergence of nanoscale 2D materials with magnetic phenomena, inductively detected NMR and NQR spectroscopy are not sensitive enough to detect the smaller number of spins in nanomaterials. The nitrogen-vacancy (NV) center in diamond has shown promise in bringing the analytic power of NMR and NQR spectroscopy to the nanoscale. However, due to depth-dependent formation efficiency of the defect centers, noise from surface spins, band bending effects, and the depth dependence of the nuclear magnetic field, there is ambiguity regarding the ideal NV depth for surface NMR of statistically polarized spins. In this work, we prepared a range of shallow NV ensemble layer depths and determined the ideal NV depth by performing NMR spectroscopy on statistically polarized 19F in Fomblin oil on the diamond surface. We found that the measurement time needed to achieve a signal-to-noise ratio of 3 using XY8-N noise spectroscopy has a minimum at an NV ensemble depth of 5.5 ± 1.5 nm for ensembles activated from 100 ppm nitrogen concentration. To demonstrate the sensing capabilities of NV ensembles, we perform NQR spectroscopy on the 11B of hexagonal boron nitride flakes. We compare our best diamond to previous work with a single NV and find that this ensemble provides a shorter measurement time with excitation diameters as small as 4 μm. This analysis provides ideal conditions for further experiments involving NMR/NQR spectroscopy of 2D materials with magnetic properties.

Efforts to expand the technological capability of batteries have generated increased interest in divalent cationic systems. Electrolytes used for these electrochemical applications often incorporate cyclic ethers as electrolyte solvents; however, the detailed solvation environments within such systems are not well-understood. To foster insights into the solvation structures of such electrolytes, Ca(TFSI)2and Zn(TFSI)2dissolved in tetrahydrofuran (THF) and 2-methyl-tetrahydrofuran were investigated through multi-nuclear magnetic resonance spectroscopy (17O, 43Ca, and 67Zn NMR) combined with quantum chemistry modeling of NMR chemical shifts. NMR provides spectroscopic fingerprints that readily couple with quantum chemistry to identify a set of most probable solvation structures based on the best agreement between the theoretically predicted and experimentally measured values of chemical shifts. The multi-nuclear approach significantly enhances confidence that the correct solvation structures are identified due to the required simultaneous agreement between theory and experiment for multiple nuclear spins. Furthermore, quantum chemistry modeling provides a comparison of the solvation cluster formation energetics, allowing further refinement of the preferred solvation structures. It is shown that a range of solvation structures coexist in most of these electrolytes, with significant molecular motion and dynamic exchange among the structures. This level of solvation diversity correlates with the solubility of the electrolyte, with Zn(TFSI)2/THF exhibiting the lowest degree of each. Comparisons of analogous Ca2+and Zn2+solvation structures reveal a significant cation size effect that is manifested in significantly reduced cation-solvent bond lengths and thus stronger solvent bonding for Zn2+relative to Ca2+. The strength of this bonding is further reduced by methylation of the cyclic ether ring. Solvation shells containing anions are energetically preferred in all the studied electrolytes, leading to significant quantities of contact ion pairs and consequently neutrally charged clusters. It is likely that the transport and interfacial de-solvation/re-solvation properties of these electrolytes are directed by these anion interactions. These insights into the detailed solvation structures, cation size, and solvent effects, including the molecular dynamics, are fundamentally important for the rational design of electrolytes in multivalent battery electrolyte systems.

Lithium-sulfur is a "beyond-Li-ion" battery chemistry attractive for its high energy density coupled with low-cost sulfur. Expanding to the MWh required for grid scale energy storage, however, requires a different approach for reasons of safety, scalability, and cost. Here we demonstrate the marriage of the redox-targeting scheme to the engineered Li solid electrolyte interphase (SEI), enabling a scalable, high efficiency, membrane-less Li-S redox flow battery. In this hybrid flow battery architecture, the Li anode is housed in the electrochemical cell, while the solid sulfur is safely kept in a separate catholyte reservoir and electrolyte is pumped over the sulfur and into the electrochemical cell. Electrochemically facile decamethylferrocene and cobaltocene are chosen as redox mediators to kick-start the initial reduction of solid S into soluble polysulfides and final reduction of polysulfides into solid Li2S, precluding the need for conductive carbons. On the anode side, a LiI and LiNO3pretreatment strategy encourages a stable SEI and lessens capacity fade, avoiding use of ion-selective separators. Complementary materials characterization confirms the uniform distribution of LiI in the SEI, while SEM confirms the presence of lower surface area globular Li deposition and UV-vis corroborates evolution of the polysulfide species. Equivalent areal loadings of up to 50 mgScm-2(84 mAh cm-2) are demonstrated, with high capacity and voltage efficiency at 1-2 mgScm-2(973 mAh gS-1and 81.3% VE in static cells and 1142 mAh gS-1and 86.9% VE in flow cells). These results imply that the fundamental Li-S chemistry and SEI engineering strategies can be adapted to the hybrid redox flow battery architecture, obviating the need for ion-selective membranes or flowing carbon additives, and offering a potential pathway for inexpensive, scalable, and safe MWh scale Li-S energy storage.

Aqueous aluminum (Al) batteries are posited to be a cheap and energy dense alternative to conventional Li-ion chemistries, but an aqueous electrolyte mediating trivalent aluminum cations (Al3+) warrants greater scrutiny. This study provides a rigorous examination of aqueous Al electrolytes, with the first compelling evidence for a dynamic octahedral solvation structure around the Al3+, without Al-OTf contact ion pairs, even at high concentrations. This solvation behavior and the concomitant, transient electrostatic hydrolysis of Al-OH2 ligands contrasts strongly with previously reported water-in-salt electrolytes, and occurs due to the high charge density of the Lewis acidic Al3+. Nuclear magnetic resonance spectroscopy and other physicochemical measurements quantitatively reveal how species activity evolves with concentration and temperature. This new understanding exposes practical concerns related to the corrosiveness of the acidic aqueous solutions, the degree of hydration of aluminum trifluoromethanesulfonate (Al(OTf)3) salt, and the grossly insufficient reductive stability of the proposed electrolytes (>1 V between HER onset and Al3+/Al). Collectively, these factors constitute multiple fundamental barriers to the feasibility of rechargeable aqueous Al batteries.

We demonstrate that stimulated microwave optical sideband generation using parametric frequency conversion can be utilized as a powerful technique for coherent state detection in atomic physics experiments. The technique has advantages over traditional absorption or polarization rotation-based measurements and enables the isolation of signal photons from probe photons. We outline a theoretical framework that accurately models sideband generation using a density matrix formalism. Using this technique, we demonstrate a novel intrinsic magnetic gradiometer that detects magnetic gradient fields between two spatially separated vapor cells by measuring the frequency of the beat note between sidebands generated within each cell. The sidebands are produced with high efficiency using parametric frequency conversion of a probe beam interacting with Rb87 atoms in a coherent superposition of magnetically sensitive hyperfine ground states. Interference between the sidebands generates a low-frequency beat note whose frequency is determined by the magnetic field gradient between the two vapor cells. In contrast to traditional gradiometers the intermediate step of measuring the magnetic field experienced by the two vapor cells is unnecessary. We show that this technique can be readily implemented in a practical device by demonstrating a compact magnetic gradiometer sensor head with a sensitivity of 25 fT/cm/Hz with a 4.4 cm baseline, while operating in a noisy laboratory environment unshielded from Earth's field.

2,4,dimethyloxetane is an important cyclic ether intermediate that is produced from hydroperoxyalkyl (QOOH) radicals in the low-temperature combustion of n-pentane. However, the reaction mechanisms and rates of consumption pathways remain unclear. In the present work, the pressure- and temperature-dependent kinetics of seven cyclic ether peroxy radicals, which stem from 2,4,dimethyloxetane via H-abstraction and O2 addition, were determined. The automated kinetic workflow code, KinBot, was used to model the complexity of the chemistry in a stereochemically resolved manner and solve the resulting master equations from 300-1000 K and from 0.01-100 atm. The main conclusions from the calculations include (i) diastereomeric cyclic ether peroxy radicals show significantly different reactivities, (ii) the stereochemistry of the peroxy radical determines which QOOH isomerization steps are possible, (iii) conventional QOOH decomposition pathways, such as cyclic ether formation and HO2 elimination, compete with ring-opening reactions, which primarily produce OH radicals, the outcome of which is sensitive to stereochemistry. Ring-opening reactions lead to unique products, such as unsaturated, acyclic peroxy radicals, that form direct connections with species present in other chemical kinetics mechanisms through "cross-over" reactions that may complicate the interpretation of experimental results from combustion of n-pentane and, by extension, other alkanes. For example, one cross-over reaction involving 1-hydroperoxy-4-pentanone-2-yl produces 2-(hydroperoxymethyl)-3-butanone-1-yl, which is an iso-pentane-derived ketohydroperoxide (KHP). At atmospheric pressure, the rate of chemical reactions of all seven peroxy radicals compete with that of collisional stabilization, resulting in well-skipping reactions. However, at 100 atm, only one out of seven peroxy radicals undergoes significant well-skipping reactions. The rates produced from the master equation calculations provide the first foundation for the development of detailed sub-mechanisms for cyclic ether intermediates. In addition, analysis of the complex reaction mechanisms of 2,4-dimethyloxetane-derived peroxy radicals provides insights into the effects of stereoisomers on reaction pathways and product yields.

All-solid-state batteries are often assumed to be safer than conventional Li-ion ones. In this work, we present the first thermodynamic models to quantitatively evaluate solid-state and Li-ion battery heat release under several failure scenarios. The solid-state battery analysis is carried out with an Li7La3Zr2O12 solid electrolyte but can be extended to other configurations using the accompanying spreadsheet. We consider solid-state batteries that include a relatively small amount of liquid electrolyte, which is often added at the cathode to reduce interfacial resistance. While the addition of small amounts of liquid electrolyte increases heat release under specific failure scenarios, it may be small enough that other considerations, such as manufacturability and performance, are more important commercially. We show that short-circuited all-solid-state batteries can reach temperatures significantly higher than conventional Li-ion, which could lead to fire through flammable packaging and/or nearby materials. Our work highlights the need for quantitative safety analyses of solid-state batteries.

Finite-temperature Kohn-Sham density functional theory (KS-DFT) is a widely-used method in warm dense matter (WDM) simulations and diagnostics. Unfortunately, full KS-DFT-molecular dynamics models scale unfavourably with temperature and there remains uncertainty regarding the performance of existing approximate exchange-correlation (XC) functionals under WDM conditions. Of particular concern is the expected explicit dependence of the XC functional on temperature, which is absent from most approximations. Average-atom (AA) models, which significantly reduce the computational cost of KS-DFT calculations, have therefore become an integral part of WDM modeling. In this paper, we present a derivation of a first-principles AA model from the fully-interacting many-body Hamiltonian, carefully analyzing the assumptions made and terms neglected in this reduction. We explore the impact of different choices within this model—such as boundary conditions and XC functionals—on common properties in WDM, for example equation-of-state data, ionization degree and the behavior of the frontier energy levels. Furthermore, drawing upon insights from ground-state KS-DFT, we discuss the likely sources of error in KS-AA models and possible strategies for mitigating such errors.

The goal of this project is to reconsider core market and reliability processes that can potentially yield to transformative advances in power grid security, reliability, and efficiency. Current electric power market designs are strongly a function of computing capabilities and limitations that were available in the mid-to-late 1990s, circa deregulation. This includes constructs such as: (1) a 2-tiered day-ahead/real-time market construct; and (2) linearized (“DC”) real power flow approximations in dispatch and pricing. At that time, state-of-the-art computational capabilities could at the limit address deterministic mixed-integer programming formulations of unit commitment (UC) and linear programming formulations of economic dispatch (ED) at limited fidelity and scale. Such constraints forced limited look-ahead time-horizons, crude approximations of AC power flow physics and operations, and artificial partitioning between day-ahead markets, hour(s)-ahead reliability processes, and real-time markets. Consequently, these limitations have resulted in limited security and reliability with increasing out-of-market payments, particularly as uncertainty associated with renewables and distributed energy resources grows.

This guide is meant to assist communities – from residents to energy experts to decision makers – in developing a conceptual microgrid design that meets site-specific energy resilience goals. Using the framework described in this guidebook, stakeholders can come together and start to quantify site-specific vulnerabilities, identify the most significant risks to delivery of electricity, and establish electric outage tolerances across the community. In addition to establishing minimum service needs, this framework encourages communities to consider broader sustainability goals and policy constraints and begin to estimate up-front costs associated with the installation of alternative microgrid solutions. The framework guides a community through data collection and a high-level assessment of its needs, constraints, and priorities, prior to engaging engineers, vendors, and contractors. The first sections of this guidebook provide a high-level primer on electric systems. The latter sections include guidance for step-by-step data gathering and analysis of site conditions. The ultimate product resulting from the stepwise approach is a conceptual microgrid design. A conceptual design is defined as an initial design (10%-20% complete) that considers the specific threats, needs, limitations, and investment options for a given location. Going through this exercise and developing the conceptual microgrid design as a community ensures the same community members who will ultimately live with the solution are the developers of its foundational design. Often, these are also the very same people who understand system tolerances and needs the best and are therefore the ideal candidates for establishing these criteria. Especially when it comes to evaluating critical infrastructure, it is the community that best understands the most critical services. The framework is intended to facilitate a systematic approach to planning for resilience and provide a deeper understanding of how to use a framework to make decisions around microgrid solutions. Like many processes where tradeoffs need to be considered, this is often an iterative process. If this guide serves to help educate and empower communities who are beginning the process of deploying a microgrid, it has met the goal of its authors.

High-pressure storage and cyclic (de)pressurization of hydrogen gas is known to result in degradation and failure of gas canisters, hoses, linings, and O-rings as the relatively small hydrogen molecule can readily permeate most materials. Hence, identifying material compositions that are less susceptible to hydrogen-induced damage is of significant importance to the hydrogen energy infrastructure. Here, we use classical atomistic molecular dynamics simulations to study hydrogen exposed ethylene-propylene-diene monomer (EPDM) rubber, an elastomer typically used in O-rings. We make chemical modifications to the model by adjusting the crosslink density and report on gas solubility, diffusivity, and molecular restructuring in response to rapid decompression. Our simulations indicate that increases in crosslink density can reduce volumetric expansion during decompression and result in smaller free volume pore sizes. However, these favorable properties for sealing materials come with a tradeoff. At pressure, crosslinks introduce extra free volume, providing potential sites for gas localization, the precursor to cavitation-induced failure.

Recent high profile cyber attacks on critical infrastructures have raised awareness about the severe and widespread impacts that these attacks can have on everyday life. This awareness has spurred research into making industrial control systems and other cyber-physical systems more resilient. A plethora of cyber resilience metrics and frameworks have been proposed for cyber resilience assessments, but these approaches typically assume that data required to populate the metrics is readily available, an assumption that is frequently not valid. This paper describes a new cyber experimentation platform that can be used to generate relevant data and to calculate resilience metrics that quantify how resilient specified industrial control systems are to specified threats. Demonstration of the platform and analysis process are illustrated through a use case involving the control system for a pressurized water reactor.

Social systems are uniquely complex and difficult to study, but understanding them is vital to solving the world’s problems. The Ground Truth program developed a new way of testing the research methods that attempt to understand and leverage the Human Domain and its associated complexities. The program developed simulations of social systems as virtual world test beds. Not only were these simulations able to produce data on future states of the system under various circumstances and scenarios, but their causal ground truth was also explicitly known. Research teams studied these virtual worlds, facilitating deep validation of causal inference, prediction, and prescription methods. The Ground Truth program model provides a way to test and validate research methods to an extent previously impossible, and to study the intricacies and interactions of different components of research.

A computationally efficient radiative transport model is presented that predicts a camera measurement and accounts for the light reflected and blocked by an object in a scattering medium. The model is in good agreement with experimental data acquired at the Sandia National Laboratory Fog Chamber Facility (SNLFC). The model is applicable in computational imaging to detect, localize, and image objects hidden in scattering media. Here, a statistical approach was implemented to study object detection limits in fog.

A major challenge in shape optimization is the coupling of finite element method (FEM) codes in a way that facilitates efficient computation of shape derivatives. This is particularly difficult with multiphysics problems involving legacy codes, where the costs of implementing and maintaining shape derivative capabilities are prohibitive. The volume and boundary methods are two approaches to computing shape derivatives. Each has a major drawback: the boundary method is less accurate, while the volume method is more invasive to the FEM code. We introduce the strip method, which computes shape derivatives on a strip adjacent to the boundary. The strip method makes code coupling simple. Like the boundary method, it queries the state and adjoint solutions at quadrature nodes, but requires no knowledge of the FEM code implementations. At the same time, it exhibits the higher accuracy of the volume method. As an added benefit, its computational complexity is comparable to that of the boundary method, that is, it is faster than the volume method. We illustrate the benefits of the strip method with numerical examples.

This report describes the proposed surface sampling techniques and plan for the multi-year Canister Deposition Field Demonstration (CDFD). The CDFD is primarily a dust deposition test that will use three commercial 32PTH2 NUHOMS welded stainless steel storage canisters in Advanced Horizontal Storage Modules, with planned exposure testing for up to 10 years at an operating ISFSI site. One canister will be left at ambient condition, unheated; the other two will have heaters to achieve canister surface temperatures that match, to the degree possible, spent nuclear fuel (SNF) loaded canisters with heat loads of 10 kW and 40 kW. Surface sampling campaigns for dust analysis will take place on a yearly or bi-yearly basis. The goal of the planned dust sampling and analysis is to determine important environmental parameters that impact the potential occurrence of stress corrosion cracking on SNF dry storage canisters. Specifically, measured dust deposition rates and deposited particle sizes will improve parameterization of dust deposition models employed to predict the potential occurrence and timing of stress corrosion cracks on the stainless steel SNF canisters. The size, morphology, and composition of the deposited dust and salt particles will be quantified, as well as the soluble salt load per unit area and the rate of deposition, as a function of canister surface temperature, location, time, and orientation. Previously, a preliminary sampling plan was developed, identifying possible sampling locations on the canister surfaces and sampling intervals; possible sampling methods were also described. Further development of the sampling plan has commenced through three different tasks. First, canister surface roughness, a potentially important parameter for air flow and dust deposition, was characterized at several locations on one of the test canisters. Second, corrosion testing to evaluate the potential lifetime and aging of thermocouple wires, spot welds, and attachments was initiated. Third, hand sampling protocols were developed, and initial testing was carried out. The results of those efforts are presented in this report. The information obtained from the CDFD will be critical for ongoing efforts to develop a detailed understanding of the potential for stress corrosion cracking of SNF dry storage canisters.

Systems subjected to dynamic loads often require monitoring of their vibrational response, but limitations on the total number and placement of the measurement sensors can hinder the data-collection process. This paper presents an indirect approach to estimate a system's full-field dynamic response, including all uninstrumented locations, using response measurements from sensors sparsely located on the system. This approach relies on Bayesian inference that utilizes a system model to estimate the full-field response and quantify the uncertainty in these estimates. By casting the estimation problem in the frequency domain, this approach utilizes the modal frequency response functions as a natural, frequency-dependent weighting scheme for the system mode shapes to perform the expansion. This frequency-dependent weighting scheme enables an accurate expansion, even with highly correlated mode shapes that may arise from spatial aliasing due to the limited number of sensors, provided these correlated modes do not have natural frequencies that are closely spaced. Furthermore, the inherent regularization mechanism that arises in this Bayesian-based procedure enables the utilization of the full set of system mode shapes for the expansion, rather than any reduced subset. This approach can produce estimates when considering a single realization of the measured responses, and with some modification, it can also produce estimates for power spectral density matrices measured from many realizations of the responses from statistically stationary random processes. A simply supported beam provides an initial numerical validation, and a cylindrical test article excited by acoustic loads in a reverberation chamber provides experimental validation.

A major challenge in shape optimization is the coupling of finite element method (FEM) codes in a way that facilitates efficient computation of shape derivatives. This is particularly difficult with multiphysics problems involving legacy codes, where the costs of implementing and maintaining shape derivative capabilities are prohibitive. The volume and boundary methods are two approaches to computing shape derivatives. Each has a major drawback: the boundary method is less accurate, while the volume method is more invasive to the FEM code. We introduce the strip method, which computes shape derivatives on a strip adjacent to the boundary. The strip method makes code coupling simple. Like the boundary method, it queries the state and adjoint solutions at quadrature nodes, but requires no knowledge of the FEM code implementations. At the same time, it exhibits the higher accuracy of the volume method. As an added benefit, its computational complexity is comparable to that of the boundary method, that is, it is faster than the volume method. We illustrate the benefits of the strip method with numerical examples.

The crystallization of an amorphous Ag–In–Sb–Te (AIST) phase change material (PCM) is studied using multiple in situ imaging techniques to directly quantify crystal growth rates over a broad range of temperatures. The measurable growth rates span from ≈ 10–9 to ≈ 20 m/s. Recent results using dynamic transmission electron microscopy (TEM), a photoemission TEM technique, and TEM with sub-framed imaging are reported here and placed into the context of previous growth rate measurements on AIST. Dynamic TEM experiments show a maximum observed crystal growth rate for as-deposited films to be > 20 m/s. It is shown that crystal growth above the glass transition can be imaged in a TEM through use of subframing and a high-frame-rate direct electron detection camera. Challenges associated with the determination of temperature during in situ TEM experiments are described. Preliminary nanocalorimetry results demonstrate the feasibility of collecting thermodynamic data for crystallization of PCMs with simultaneous TEM imaging. Graphical abstract: [Figure not available: see fulltext.]

We employ ultrafast mid-infrared transient absorption spectroscopy to probe the rapid loss of carbonyl ligands from gas-phase nickel tetracarbonyl following ultraviolet photoexcitation at 261 nm. Here, nickel tetracarbonyl undergoes prompt dissociation to produce nickel tricarbonyl in a singlet excited state; this electronically excited tricarbonyl loses another CO group over tens of picoseconds. Our results also suggest the presence of a parallel, concerted dissociation mechanism to produce nickel dicarbonyl in a triplet excited state, which likely dissociates to nickel monocarbonyl. Mechanisms for the formation of these photoproducts in multiple electronic excited states are theoretically predicted with one-dimensional cuts through the potential energy surfaces and computation of spin-orbit coupling constants using equation of motion coupled cluster methods (EOM-CC) and coupled cluster theory with single and double excitations (CCSD). Bond dissociation energies are calculated with CCSD, and anharmonic frequencies of ground and excited state species are computed using density functional theory (DFT) and time-dependent density functional theory (TD-DFT).

Swelling clay hydration/dehydration is important to many environmental and industrial processes. Experimental studies usually probe equilibrium hydration states in an averaged manner and thus cannot capture the fast water transport and structural change in interlayers during hydration/dehydration. Using molecular simulations and thermogravimetric analyses, we observe a two-stage dehydration process. The first stage is controlled by evaporation at the edges: water molecules near hydrophobic sites and the first few water molecules of the hydration shell of cations move fast to particle edges for evaporation. The second stage is controlled by slow desorption of the last 1-2 water molecules from the cations and slow transport through the interlayers. The two-stage dehydration is strongly coupled with interlayer collapse and the coordination number changes of cations, all of which depend on layer charge distribution. This mechanistic interpretation of clay dehydration can be key to the coupled chemomechanical behavior in natural/engineered barriers.

Abstract not provided.

Chemical interactions on the surface of a functional nanoparticle are closely related to its crystal facets, which can regulate the corresponding energy storage properties like hydrogen absorption. In this study, we reported a one-step growth of magnesium (Mg) particles with both close- and nonclose-packed facets, that is, {0001} and {21¯ 1¯ 6} planes, on atomically thin reduced graphene oxide (rGO). The detailed microstructures of Mg/rGO hybrids were revealed by X-ray diffraction, selected-area electron diffraction, high-resolution transmission electron microscopy, and fast Fourier transform analysis. Hydrogen storage performance of Mg/rGO hybrids with different orientations varies: Mg with preferential high-index {21¯ 1¯ 6} crystal surface shows remarkably increased hydrogen absorption up to 6.2 wt % compared with the system exposing no preferentially oriented crystal surfaces showing inferior performance of 5.1 wt % within the first 2 h. First-principles calculations revealed improved hydrogen sorption properties on the {21¯ 1¯ 6} surface with a lower hydrogen dissociation energy barrier and higher stability of hydrogen atoms than those on the {0001} basal plane, supporting the hydrogen uptake experiment. In addition, the hydrogen penetration energy barrier is found to be much lower than that of {0001} because of low surface atom packing density, which might be the most critical process to the hydrogenation kinetics. The experimental and calculation results present a new handle for regulating the hydrogen storage of metal hydrides by controlled Mg facets.

The resurgence of interest in a hydrogen economy and the development of hydrogen-related technologies has initiated numerous research and development efforts aimed at making the generation, storage, and transportation of hydrogen more efficient and affordable. Solar thermochemical hydrogen production (STCH) is a process that potentially exhibits numerous benefits such as high reaction efficiencies, tunable thermodynamics, and continued performance over extended cycling. Although CeO₂ has been the de facto standard STCH material for many years, more recently 12R-Ba₄CeMn₃O₁₂ (BCM) has demonstrated enhanced hydrogen production at intermediate H₂/H₂O conditions compared to CeO₂, making it a contender for large-scale hydrogen production. However, the thermo-reduction stability of 12R-BCM dictates the oxygen partial pressure (pO₂) and temperature conditions optimal for cycling. In this study, we identify the formation of a 6H-BCM polytype at high temperature and reducing conditions, experimentally and computationally, as a mechanism and pathway for 12R-BCM decomposition. 12R-BCM was synthesized with high purity and then controllably reduced using thermogravimetric analysis (TGA). Synchrotron X-ray diffraction (XRD) data is used to identify the formation of a 6H-Ba₃Ce_0.75Mn_2.25O₉ (6H-BCM) polytype that is formed at 1350 degrees C under strongly reducing pO₂. Density functional theory (DFT) total energy and defect calculations show a window of thermodynamic stability for the 6H-polytype consistent with the XRD results. These data provide the first evidence of the 6H-BCM polytype and could provide a mechanistic explanation for the superior water-splitting behaviors of 12R-BCM.

The Synchronic Web is a network of information that is locked into a single global view of history. When clients notarize their data to the Synchronic Web, they gain the ability to irrefutably prove the following statement to the rest of the world: "I commit to this information—and only this information—at this moment in time." Much like encryption or digital signatures, this capability has the potential to bolster the integrity of public cyberspace at a foundational level and scale.

In April 5-7, 2022, Sandia National Laboratories hosted a second virtual workshop to further explore the potential for developing AI-enhanced co-design for microelectronics (AICoM). This second piece in an ongoing workshop series again brought together two themes. The first theme, co-design for next generation microelectronics, was drawn from the 2018 Department of Energy Office of Science (DOE SC) “Basic Research Needs for Microelectronics” (BRN) report (DOE/SC, 2018, 2021), which called for a “fundamental rethinking” of the traditional design approach to microelectronics, in which subject matter experts (SMEs) in each microelectronics discipline (materials, devices, circuits, algorithms, etc.) work near-independently. Instead, the BRN called for a non-hierarchical, egalitarian vision of co-design, wherein “each scientific discipline informs and engages the others” in “parallel but intimately networked efforts to create radically new capabilities.” The second theme, exploiting and advancing artificial intelligence (AI) to support co-design for microelectronics, acknowledges the continuing breakthroughs in AI that are currently enhancing and accelerating solutions to traditional design problems in materials synthesis and processing, circuit design, and electronic design automation (EDA).

Interactions of ceramic proton conductors with the environment under operating conditions play an essential role on material properties and device performance. It remains unclear how the chemical environment of material, as modulated by the operating condition, affects the proton conductivity. Combining near-ambient pressure X-ray photoelectron spectroscopy and impedance spectroscopy, we investigate the chemical environment changes of oxygen and the conductivity of BaZr0.9Y0.1O3-δunder operating condition. Changes in O 1s core level spectra indicate that adding water vapor pressure increases both hydroxyl groups and active proton sites at undercoordinated oxygen. Applying external potential further promotes this hydration effect, in particular, by increasing the amount of undercoordinated oxygen. The enhanced hydration is accompanied by improved proton conductivity. This work highlights the effects of undercoordinated oxygen for improving the proton conductivity in ceramics.

Ion trap quantum computing utilizes electronic states of atomic ions such as Ca+ to encode information on to a qubit. To explore the fundamental properties of Ca+ inside molecular cavities, we describe here a computational study of Ca+ bound inside neutral [n]-cycloparaphenylenes (n = 5-12), often referred to as “nanohoops”. This ab initio study characterizes optimized structures, harmonic vibrational frequencies, potential energy surfaces, and ion molecular orbital distortion as functions of increasing nanohoop size. The results of this work provide a first step in guiding experimental studies of the spectroscopy of these ion-molecular cavity complexes.

Aria is a Galerkin nite element based program for solving coupled-physics problems described by systems of PDEs and is capable of solving nonlinear, implicit, transient and direct-to-steady state problems in two and three dimensions on parallel architectures. The suite of physics currently supported by Aria includes thermal energy transport, species transport, and electrostatics as well as generalized scalar, vector and tensor transport equations. Additionally, Aria includes support for manufacturing process flows via the incompressible Navier-Stokes equations specialized to a low Reynolds number (Re < 1) regime. Enhanced modeling support of manufacturing processing is made possible through use of either arbitrary Lagrangian-Eulerian (ALE) and level set based free and moving boundary tracking in conjunction with quasi-static nonlinear elastic solid mechanics for mesh control. Coupled physics problems are solved in several ways including fully-coupled Newton’s method with analytic or numerical sensitivities, fully-coupled Newton-Krylov methods and a loosely-coupled nonlinear iteration about subsets of the system that are solved using combinations of the aforementioned methods. Error estimation, uniform and dynamic ℎ-adaptivity and dynamic load balancing are some of Aria’s more advanced capabilities.

The SIERRA Low Mach Module: Fuego, henceforth referred to as Fuego, is the key element of the ASC fire environment simulation project. The fire environment simulation project is directed at characterizing both open large-scale pool fires and building enclosure fires. Fuego represents the turbulent, buoyantly-driven incompressible ow, heat transfer, mass transfer, combustion, soot, and absorption coefficient model portion of the simulation software. Using MPMD coupling, Scefire and Nalu handle the participating-media thermal radiation mechanics. This project is an integral part of the SIERRA multi-mechanics software development project. Fuego depends heavily upon the core architecture developments provided by SIERRA for massively parallel computing, solution adaptivity, and mechanics coupling on unstructured grids.

Aria is a Galerkin finite element based program for solving coupled-physics problems described by systems of PDEs and is capable of solving nonlinear, implicit, transient and direct-to-steady state problems in two and three dimensions on parallel architectures. The suite of physics currently supported by Aria includes thermal energy transport, species transport, and electrostatics as well as generalized scalar, vector and tensor transport equations. Additionally, Aria includes support for manufacturing process flows via the incompressible Navier-Stokes equations specialized to a low Reynolds number ( $Re$ < 1 ) regime. Enhanced modeling support of manufacturing processing is made possible through use of either arbitrary Lagrangian-Eulerian (ALE) and level set based free and moving boundary tracking in conjunction with quasi-static nonlinear elastic solid mechanics for mesh control. Coupled physics problems are solved in several ways including fully-coupled Newton's method with analytic or numerical sensitivities, fully-coupled Newton-Krylov methods and a loosely-coupled nonlinear iteration about subsets of the system that are solved using combinations of the aforementioned methods. Error estimation, uniform and dynamic $h$-adaptivity and dynamic load balancing are some of Aria's more advanced capabilities.

The SIERRA Low Mach Module: Fuego, henceforth referred to as Fuego, is the key element of the ASC fire environment simulation project. The fire environment simulation project is directed at characterizing both open large-scale pool fires and building enclosure fires. Fuego represents the turbulent, buoyantly-driven incompressible flow, heat transfer, mass transfer, combustion, soot, and absorption coefficient model portion of the simulation software. Using MPMD coupling, Scefire and Nalu handle the participating-media thermal radiation mechanics. This project is an integral part of the SIERRA multi-mechanics software development project. Fuego depends heavily upon the core architecture developments provided by SIERRA for massively parallel computing, solution adaptivity, and mechanics coupling on unstructured grids.

SIERRA/Aero is a compressible fluid dynamics program intended to solve a wide variety compressible fluid flows including transonic and hypersonic problems. This document describes the commands for assembling a fluid model for analysis with this module, henceforth referred to simply as Aero for brevity. Aero is an application developed using the SIERRA Toolkit (STK). The intent of STK is to provide a set of tools for handling common tasks that programmers encounter when developing a code for numerical simulation. For example, components of STK provide field allocation and management, and parallel input/output of field and mesh data. These services also allow the development of coupled mechanics analysis software for a massively parallel computing environment.

Presented in this document is a portion of the tests that exist in the Sierra Thermal/Fluids verification test suite. Each of these tests is run nightly with the Sierra/TF code suite and the results of the test checked under mesh refinement against the correct analytic result. For each of the tests presented in this document the test setup, derivation of the analytic solution, and comparison of the code results to the analytic solution is provided. This document can be used to confirm that a given code capability is verified or referenced as a compilation of example problems.

SIERRA/Aero is a compressible fluid dynamics program intended to solve a wide variety compressible fluid flows including transonic and hypersonic problems. This document describes the commands for assembling a fluid model for analysis with this module, henceforth referred to simply as Aero for brevity. Aero is an application developed using the SIERRA Toolkit (STK). The intent of STK is to provide a set of tools for handling common tasks that programmers encounter when developing a code for numerical simulation. For example, components of STK provide field allocation and management, and parallel input/output of field and mesh data. These services also allow the development of coupled mechanics analysis software for a massively parallel computing environment.

The SNL Sierra Mechanics code suite is designed to enable simulation of complex multiphysics scenarios. The code suite is composed of several specialized applications which can operate either in standalone mode or coupled with each other. Arpeggio is a supported utility that enables loose coupling of the various Sierra Mechanics applications by providing access to Framework services that facilitate the coupling. More importantly Arpeggio orchestrates the execution of applications that participate in the coupling. This document describes the various components of Arpeggio and their operability. The intent of the document is to provide a fast path for analysts interested in coupled applications via simple examples of its usage.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The nuclear accident consequence analysis code MACCS has traditionally modeled dispersion during downwind transport using a Gaussian plume segment model. MACCS is designed to estimate consequence measures such as air concentrations and ground depositions, radiological doses, and health and economic impacts on a statistical basis over the course of a year to produce annualaveraged output measures. The objective of this work is to supplement the Gaussian atmospheric transport and diffusion (ATD) model currently in MACCS with a new option using the HYSPLIT model. HYSPLIT/MACCS coupling has been implemented, with HYSPLIT as an alternative ATD option. The subsequent calculations in MACCS use the HYSPLIT-generated air concentration, and ground deposition values to calculate the same range of output quantities (dose, health effects, risks, etc.) that can be generated when using the MACCS Gaussian ATD model. Based on the results from the verification test cases, the implementation of the HYSPLIT/MACCS coupling is confirmed. This report contains technical details of the HYSPLIT/MACCS coupling and presents a benchmark analysis using the HYSPLIT/MACCS coupling system. The benchmark analysis, which involves running specific scenarios and sensitivity studies designed to examine how the results generated by the traditional MACCS Gaussian plume segment model compare to the new, higher fidelity HYSPLIT/MACCS modeling option, demonstrates the modeling results that can be obtained by using this new option. The comparisons provided herein can also help decision-makers evaluate the potential benefit of using results based on higher fidelity modeling with the additional computational burden needed to perform the calculations. Three sensitivity studies to investigate the potential impact of alternative modeling options, regarding 1) input meteorological data set, 2) method to estimate stability class, and 3) plume dispersion model for larger distances, on consequence results were also performed. The results of these analyses are provided and discussed in this report.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Wind energy can provide renewable, sustainable electricity to rural Native homes and power schools and businesses. It can even provide tribes with a source of income and economic development. The purpose of this research is to determine the potential for deploying community and utility-scale wind renewable technologies on Turtle Mountain Band of Chippewa tribal lands. Ideal areas for wind technology development were investigated, based on wind resources, terrain, land usage, and other factors. This was done using tools like the National Renewable Energy Laboratory Wind Prospector, in addition to consulting tribal members and experts in the field. The result was a preliminary assessment of wind energy potential on Turtle Mountain lands, which can be used to justify further investigation and investment into determining the feasibility of future wind technology projects.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The objective of this project is the demonstration, and validation of hydrogen fuel cells in the marine environment. The prototype generator can be used to guide commercial development of a fuel cell generator product. Work includes assessment and validation of the commercial value proposition of both the application and the hydrogen supply infrastructure through third-party hosted deployment as the next step towards widespread use of hydrogen fuel cells in the maritime environment.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Silicon spin qubits satisfy the necessary criteria for quantum information processing. However, a demonstration of high-fidelity state preparation and readout combined with high-fidelity single- and two-qubit gates, all of which must be present for quantum error correction, has been lacking. We use a two-qubit Si/SiGe quantum processor to demonstrate state preparation and readout with fidelity greater than 97%, combined with both singleand two-qubit control fidelities exceeding 99%. The operation of the quantum processor is quantitatively characterized using gate set tomography and randomized benchmarking. Our results highlight the potential of silicon spin qubits to become a dominant technology in the development of intermediate-scale quantum processors.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

We demonstrate SONOS (silicon-oxide-nitride-oxide-silicon) analog memory arrays that are optimized for neural network inference. The devices are fabricated in a 40nm process and operated in the subthreshold regime for in-memory matrix multiplication. Subthreshold operation enables low conductances to be implemented with low error, which matches the typical weight distribution of neural networks, which is heavily skewed toward near-zero values. This leads to high accuracy in the presence of programming errors and process variations. We simulate the end-To-end neural network inference accuracy, accounting for the measured programming error, read noise, and retention loss in a fabricated SONOS array. Evaluated on the ImageNet dataset using ResNet50, the accuracy using a SONOS system is within 2.16% of floating-point accuracy without any retraining. The unique error properties and high On/Off ratio of the SONOS device allow scaling to large arrays without bit slicing, and enable an inference architecture that achieves 20 TOPS/W on ResNet50, a > 10× gain in energy efficiency over state-of-The-Art digital and analog inference accelerators.

There has been ever-growing interest and engagement regarding net-zero and carbon neutrality goals, with many nations committing to steep emissions reductions by mid-century. Although water plays critical roles in various sectors, there has been a distinct gap in discussions to date about the role of water in the transition to a carbon neutral future. To address this need, a webinar was convened in April 2022 to gain insights into how water can support or influence active strategies for addressing emissions activities across energy, industrial, and carbon sectors. The webinar presentations and discussions highlighted various nuances of direct and indirect water use both within and across technology sectors (Figure ES-1). For example, hydrogen and concrete production, water for mining, and inland waterways transportation are all heavily influenced by the energy sources used (fossil fuels vs. renewable sources) as well as local resource availabilities. Algal biomass, on the other hand, can be produced across diverse geographies (terrestrial to sea) in a range of source water qualities, including wastewater and could also support pollution remediation through nutrient and metals recovery. Finally, water also influences carbon dynamics and cycling within natural systems across terrestrial, aquatic, and geologic systems. These dynamics underscore not only the critical role of water within the energy-water nexus, but also the extension into the energy-watercarbon nexus.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The Health Management Clinic (HMC) is a worksite specialty clinic designed to provide an exceptional level of health care for Sandia employees with diabetes, cholesterol and blood pressure disorders, and for those employees that need help with smoking cessation, depression, anxiety, sleep disorders, or weight management. With a unified commitment to the best care practices available, the HMC is Sandia’s interface to workplace healthcare and health plan services. The HMC provides Sandia employees access to onsite screenings, health care exams, preventative health education, disease management education, care management, periodic laboratory testing, immunizations, podiatry services, and behavioral, fitness, and nutrition counseling/education. Our multidisciplinary team of health professionals consists of physicians, nurses, medical assistants, certified diabetes educators, dietitians, health educators, and exercise specialists. Services offered by the Health Management clinic have been designed to reduce further complications from disease states and promote healthy behavior changes for Sandia employees.

Abstract not provided.

The inverse methods team provides a set of tools for solving inverse problems in structural dynamics and thermal physics, and also sensor placement optimization via Optimal Experimental Design (OED). These methods are used for designing experiments, model calibration, and verfication/validation analysis of weapons systems. This document provides a user's guide to the input for the three apps that are supported for these methods. Details of input specifications, output options, and optimization parameters are included.

Abstract not provided.

The Department of Energy (DOE) is the owner and part operator of multiple facilities in Northern California. The facilities include those located at Lawrence Livermore National Laboratory (LLNL), Lawrence Berkeley National Laboratory (LBNL), Sandia National Laboratories/California (SNL/CA) and SLAC National Accelerator Laboratory (SLAC) among other sites. Through their operations, the facilities generate hazardous waste and, thereby, are subject to the requirements of Chapter 31 of the Title 22 California Code of Regulations, Waste Minimization. The Northern California sites are primarily research and development facilities in the areas relating to national security, high-energy physics, engineering, bioscience and environmental health and safety disciplines. As mentioned above these DOE facilities are primarily research and development facilities. The hazardous wastes generated may be associated with operations that range in size from small, bench scale R&D to major maintenance and operations waste streams. Therefore, even though this document breaks down the waste streams based on California Waste Codes (CWC), the quantities of waste within one waste code category could be from many different locations and dissimilar processes. Because of the nature of the work at the sites, it is not economically feasible to try to implement source reduction measures for every process that generates a portion of the waste stream. This document identifies the processes that generate the major portion of the waste within an identified major waste stream and reports on progress made toward source reduction. In accomplishing the mission, it is DOE’s goal to eliminate waste generation and emissions giving priority to those that may present the greatest risk to human health and the environment.

The Department of Energy (DOE) is the owner of multiple facilities in Northern California. The facilities include Lawrence Livermore National Laboratory (LLNL), Lawrence Berkeley National Laboratory (LBNL), Sandia National Laboratories/California (SNL/CA) and SLAC National Accelerator Laboratory (SLAC) among other sites. Through their operations, the facilities generate hazardous waste and, thereby, are subject to the requirements of Chapter 31 of the Title 22 California Code of Regulations, Waste Minimization. The Northern California sites are primarily research and development facilities in the areas relating to national security, high-energy physics, bioscience and the environment.

Abstract not provided.

A straight fiber with nonlocal forces that are independent of bond strain is considered. These internal loads can either stabilize or destabilize the straight configuration. Transverse waves with long wavelength have unstable dispersion properties for certain combinations of nonlocal kernels and internal loads. When these unstable waves occur, deformation of the straight fiber into a circular arc can lower its potential energy in equilibrium. The equilibrium value of the radius of curvature is computed explicitly.

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.

High-speed, optical imaging diagnostics are presented for three-dimensional (3D) quantification of explosively driven metal fragmentation. At early times after detonation, Digital Image Correlation (DIC) provides non-contact measures of 3D case velocities, strains, and strain rates, while a proposed stereo imaging configuration quantifies in-flight fragment masses and velocities at later times. Experiments are performed using commercially obtained RP-80 detonators from Teledyne RISI, which are shown to create a reproducible fragment field at the benchtop scale. DIC measurements are compared with 3D simulations, which have been ‘leveled’ to match the spatial resolution of DIC. Results demonstrate improved ability to identify predicted quantities-of-interest that fall outside of measurement uncertainty and shot-to-shot variability. Similarly, video measures of fragment trajectories and masses allow rapid experimental repetition and provide correlated fragment size-velocity measurements. Measured and simulated fragment mass distributions are shown to agree within confidence bounds, while some statistically meaningful differences are observed between the measured and predicted conditionally averaged fragment velocities. Together these techniques demonstrate new opportunities to improve future model validation.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Small molecules encoded by biosynthetic pathways mediate cross-species interactions and harbor untapped potential, which has provided valuable compounds for medicine and biotechnology. Since studying biosynthetic gene clusters in their native context is often difficult, alternative efforts rely on heterologous expression, which is limited by host-specific metabolic capacity and regulation. Here, in this work, we describe a computational-experimental technology to redesign genes and their regulatory regions with hybrid elements for cross-species expression in Gram-negative and -positive bacteria and eukaryotes, decoupling biosynthetic capacity from host-range constraints to activate silenced pathways. These synthetic genetic elements enabled the discovery of a class of microbiome-derived nucleotide metabolites—tyrocitabines—from Lactobacillus iners. Tyrocitabines feature a remarkable orthoester-phosphate, inhibit translational activity, and invoke unexpected biosynthetic machinery, including a class of “Amadori synthases” and “abortive” tRNA synthetases. Our approach establishes a general strategy for the redesign, expression, mobilization, and characterization of genetic elements in diverse organisms and communities.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Operability thresholds that differentiate between functional RP-87 exploding bridge wire (EBW) detonators and nonfunctional RP-87 EBW detonators (duds) were determined by measuring the time delay between initiation and early wall movement (function time). The detonators were inserted into an externally heated hollow cylinder of aluminum and fired with current flow from a charged capacitor using an exploding bridge wire (EBW initiated). Functioning detonators responded like unheated pristine detonators when the function time was 4 μs or less. The operability thresholds of the detonators were characterized with a simple decomposition cookoff model calibrated using a modified version of the Sandia Instrumented Thermal Ignition (SITI) experiment. These thresholds are based on the calculated state of the PETN when the detonators fire. The operability threshold is proportional to the positive temperature difference (ΔT) between the maximum temperature within the PETN and the onset of decomposition (∼406 K). The temperature difference alone was not sufficient to define the operability threshold. The operability threshold was also proportional to the time that the PETN had been at elevated temperatures. That is, failure was proportional to both temperature and reaction rate. The reacted gas fraction is used in the current work for the reaction correlation. Melting of PETN also had a significant effect on the operability threshold. Detonator failure occurred when the maximum temperature exceeded the nominal melting point of PETN (414 K) for 45±5 s or more.

A collection of x-ray computed tomography scans of specimens from the Bailey-Matthews National Shell Museum.

Corynebacterium glutamicum has been successfully employed for the industrial production of amino acids and other bioproducts, partially due to its native ability to utilize a wide range of carbon substrates. We demonstrated C. glutamicum as an efficient microbial host for utilizing diverse carbon substrates present in biomass hydrolysates, such as glucose, arabinose, and xylose, in addition to its natural ability to assimilate lignin-derived aromatics. As a case study to demonstrate its bioproduction capabilities, L-lactate was chosen as the primary fermentation end product along with acetate and succinate. C. glutamicum was found to grow well in different aromatics (benzoic acid, cinnamic acid, vanillic acid, and p-coumaric acid) up to a concentration of 40 mM. Besides, ¹³C-fingerprinting confirmed that carbon from aromatics enter the primary metabolism via TCA cycle confirming the presence of β-ketoadipate pathway in C. glutamicum. ¹³C-fingerprinting in the presence of both glucose and aromatics also revealed coumarate to be the most preferred aromatic by C. glutamicum contributing 74 and 59% of its carbon for the synthesis of glutamate and aspartate respectively. ¹³C-fingerprinting also confirmed the activity of ortho-cleavage pathway, anaplerotic pathway, and cataplerotic pathways. Finally, the engineered C. glutamicum strain grew well in biomass hydrolysate containing pentose and hexose sugars and produced L-lactate at a concentration of 47.9 g/L and a yield of 0.639 g/g from sugars with simultaneous utilization of aromatics. Succinate and acetate co-products were produced at concentrations of 8.9 g/L and 3.2 g/L, respectively. Our findings open the door to valorize all the major carbon components of biomass hydrolysate by using C. glutamicum as a microbial host for biomanufacturing.

Abstract not provided.

Abstract not provided.

This project is part of a multi-lab consortium that leverages U.S. research expertise and facilities at national labs and universities to significantly advance electric drive power density and reliability, while simultaneously reducing cost. The final objective of the consortium is to develop a 100 kW traction drive system that achieves 33 kW/L, has an operational life of 300,000 miles, and a cost of less than $\$$6/kW. One element of the system is a 100 kW inverter with a power density of 100 kW/L and a cost of $\$$2.7/kW. New materials such as wide bandgap semiconductors, soft magnetic materials, and ceramic dielectrics, integrated using multi-objective co optimization design techniques, will be utilized to achieve these program goals. This project focuses on a subset of the power electronics work within the consortium, specifically the design, fabrication, and evaluation of vertical GaN power devices suitable for automotive applications.

Fuel costs and emissions in maritime ports are an opportunity for transportation energy efficiency improvement and emissions reduction efforts. Ocean-going vessels, harbor craft, and cargo handling equipment are still major contributors to air pollution in and around ports. Diesel engine costs continually increase as tighter criteria pollutant regulations come into effect and will continue to do so with expected introduction of carbon emission regulations. Diesel fuel costs will also continue to rise as requirements for cleaner fuels are imposed. Both aspects will increase the cost of diesel-based power generation on the vessel and on shore. Although fuel cells have been used in many successful applications, they have not been technically or commercially validated in the port environment. One opportunity to do so was identified in Honolulu Harbor at the Young Brothers Ltd. wharf. At this facility, barges sail regularly to and from neighboring islands and containerized diesel generators provide power for the reefers while on the dock and on the barge during transport, nearly always at part load. Due to inherent efficiency characteristics of fuel cells and diesel generators, switching to a hydrogen fuel cell power generator was found to have potential emissions and cost savings. Deployment in Hawaii showed the unit needed greater reliability in the start-up sequence, as well as an improved interface to the end-user, thereby presenting opportunities for repairing/upgrading the unit for deployment in another locale. In FY2018, the unit was repaired and upgraded based on the Hawaii experience, and another deployment site was identified for another 6-month deployment of the 100 kW MarFC.

Abstract not provided.

Abstract not provided.

Abstract not provided.