Publications

Abstract not provided.

Abstract not provided.

Abstract not provided.

Role-based access control (RBAC) is adopted in the information and communication technology domain for authentication purposes. However, due to a very large number of entities within organizational access control (AC) systems, static RBAC management can be inefficient, costly, and can lead to cybersecurity threats. In this article, a novel hybrid RBAC model is proposed, based on the principles of offline deep reinforcement learning (RL) and Bayesian belief networks. The considered framework utilizes a fully offline RL agent, which models the behavioral history of users as a Bayesian belief-based trust indicator. Thus, the initial static RBAC policy is improved in a dynamic manner through off-policy learning while guaranteeing compliance of the internal users with the security rules of the system. By deploying our implementation within the smart grid domain and specifically within a Distributed Energy Resources (DER) ecosystem, we provide an end-To-end proof of concept of our model. Finally, detailed analysis and evaluation regarding the offline training phase of the RL agent are provided, while the online deployment of the hybrid RL-based RBAC model into the DER ecosystem highlights its key operation features and salient benefits over traditional RBAC models.

Abstract not provided.

Abstract not provided.

Abstract not provided.

The COVID-19 pandemic has forced many organizations—from national laboratories to private companies—to change their workforce model to incorporate remote work. This study and the summarized results sought to understand the experiences of remote workers and the ways that remote work can impact recruitment and retention, employee engagement, and career development. Sandia, like many companies, has committed to establishing a hybrid work model that will persist postpandemic, and more Sandia employees than ever before have initiated remote work agreements. This parallels the nationwide increase in remote employment and motivates this study on remote work as an enduring part of workforce models.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Metal-organic frameworks (MOFs) have recently been shown to exhibit unique mechanisms of luminescence based on charge transfer between structural units in the framework. These MOFs have the potential to be structural tuned for targeted emission with little or no metal participation. A computationally led, material design and synthesis methodology is presented here that elucidates the mechanisms of light emission in interpenetrated structures comprised of metal centers (M = In, Ga, InGa, InEu) and BTB (1,3,5-Tris(4-carboxyphenyl)benzene) linkers, forming unique luminescent M-BTB MOF frameworks. Gas phase and periodic electronic structure calculations indicate that the intensity of the emission and the wavelength are overwhelmingly controlled by a combination of the number of interacting stacked linkers and their interatomic spacings, respectively. In the MOF, the ionic radii of the metal centers primarily control the expansion or shrinkage of the linker stacking distances. Experimentally, multiple M-BTB-based MOFs are synthesized and their photoluminescence was tested. Experiments validated the modeling by confirming that shifts in the crystal structure result in variations in light emission. Through this material design method, the mechanisms of tuning luminescence properties in interpenetrated M-BTB MOFs have been identified and applied to the design of MOFs with specific wavelength emission based on their structure.

Power systems are subjected to a wide range of disturbances during daily operations. Severe disturbances, such as a loss of a large generator, a three-phase bolted fault on a generator bus, or a loss of a transmission line, can lead to the loss of synchronism of a generator or group of generators. The ability of a power system to maintain synchronism during the few seconds after being subjected to a severe disturbance is known as transient stability. Most of the modern methods of controlling transient stability involve special protection schemes or remedial action schemes. These special protection schemes sense predetermined system conditions and take corrective actions, such as generator tripping or generation re-dispatch, in real time to maintain transient stability. Another method is the use of a real-time feedback control system to modulate the output of an actuator in response to a signal. This paper provides a fundamental evaluation of the use of feedback control strategies to improve transient stability in a power system. An optimal feedback control strategy that modulates the real power injected and absorbed by distributed energy-storage devices is proposed. Its performance is evaluated on a four-machine power system and on a 34-machine reduced-order model of the Western North American Power System. The result shows that the feedback control strategy can increase the critical fault clearing time by 60%, thereby improving the transient stability of the power system.

Abstract not provided.

2021-2022 Storm Water Annual Report And Supporting Documents.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Advances in sensors and robotics have dramatically improved the diversity of experimental approaches available to the materials community. Autonomous data collection platforms, either custom-made or commercially available, provide researchers with novel tools with which to probe materials behavior and perform advanced materials characterization. The application of novel control algorithms and active learning approaches can create much more robust experimental data, or can be used to improve the performance of existing characterization tools. Five papers within this special topic focus on experimental and computational methodologies for use in automatic data collection routines for materials characterization. From novel platforms for materials discovery to new statistical frameworks for assessing the autonomous experimentation process, these five papers highlight the diverse range of applications of automation for advancing materials science.

Abstract not provided.

The Advanced Reactor Safeguards (ARS) program was established in 2020 as part of appropriations for the Advanced Reactor Demonstration Program (ARDP) through the Office of Nuclear Energy in the Department of Energy. The goal of this program is to help address near term challenges that advanced nuclear reactor vendors face in meeting domestic Material Control and Accountancy (MC&A) and Physical Protection System (PPS) requirements for U.S. construction. The technical work in the program is meant to (1) support nuclear reactor vendors with advanced MC&A and PPS designs for next generation reactors, (2) provide technical bases for the regulator, and (3) promote the integration of Safeguards and Security by Design early in the design process. Existing domestic regulations for safeguards and security, as outlined in the Code of Federal Regulations, were written for large light water reactors, and rule-making efforts are underway to develop regulations more suited to different reactor designs. The ARS program seeks to remove roadblocks in the deployment of new and advanced reactors by solving regulatory challenges, reducing safeguards and security costs, and utilizing the latest technologies and approaches for robust plant monitoring and protection. This roadmap discusses the goals of the ARS program, current research, and program plan for the next five years.

Abstract not provided.

We present an elementary and self-contained proof that there are exactly four symmetry classes of the elasticity tensor in two dimensions: oblique, rectangular, square, and isotropic. In two dimensions, orthogonal transformations are either reflections or rotations. The proof is based on identification of constraints imposed by reflections and rotations on the elasticity tensor, and it simply employs elementary tools from trigonometry, making the proof accessible to a broad audience. For completeness, we identify the sets of transformations (rotations and reflections) for each symmetry class and report the corresponding equations of motions in classical linear elasticity.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

We propose developing numerical methods informed by novel experimental diagnostics that transition from solid-to-fluid, while accurately predicting the stress and deformation regardless of phase.

This report documents the updated seismic shake table test plan. The report describes the shake table inputs (ground motions), test hardware, shake table facility, friction experiment, and proposed instrumentation.

This report summarizes the international collaborations conducted by Sandia funded by the US Department of Energy Office (DOE) of Nuclear Energy Spent Fuel and Waste Science & Technology (SFWST) as part of the Sandia National Laboratories Salt R&D and Salt International work packages. This report satisfies the level-three milestone M3SF-22SN010303063. Several stand-alone sections make up this summary report, each completed by the participants. The sections discuss international collaborations on geomechanical benchmarking exercises (WEIMOS), granular salt reconsolidation (KOMPASS), engineered barriers (RANGERS), numerical model comparison (DECOVALEX) and an NEA Salt Club working group on the development of scenarios as part of the performance assessment development process. Finally, we summarize events related to the US/German Workshop on Repository Research, Design and Operations. The work summarized in this annual update has occurred during the COVID-19 pandemic, and little international or domestic travel has occurred. Most of the collaborations have been conducted via email or as virtual meetings, but a slow return to travel and in-person meetings has begun.

This SAND Report Guide offers support to authors, technical writers, principal investigators, and others involved in the process of creating, formatting, or refining a SAND Report. It details what you need to know before you begin compiling a SAND Report, directs you to the SAND Report templates, outlines the order of elements in a SAND Report, and explains what to do when your report is completed and ready for Review and Approval and subsequent distribution. Supporting information is provided in the appendix, such as where to get technical assistance, trademarks, Microsoft Word, and equations.

CMOS-based computing systems that employ the von Neumann architecture are relatively limited when it comes to parallel data storage and processing. In contrast, the human brain is a living computational signal processing unit that operates with extreme parallelism and energy efficiency. Although numerous neuromorphic electronic devices have emerged in the last decade, most of them are rigid or contain materials that are toxic to biological systems. In this work, we report on biocompatible bilayer graphene-based artificial synaptic transistors (BLAST) capable of mimicking synaptic behavior. The BLAST devices leverage a dry ion-selective membrane, enabling long-term potentiation, with ~50 aJ/µm² switching energy efficiency, at least an order of magnitude lower than previous reports on two-dimensional material-based artificial synapses. The devices show unique metaplasticity, a useful feature for generalizable deep neural networks, and we demonstrate that metaplastic BLASTs outperform ideal linear synapses in classic image classification tasks. With switching energy well below the 1 fJ energy estimated per biological synapse, the proposed devices are powerful candidates for bio-interfaced online learning, bridging the gap between artificial and biological neural networks.

It has been demonstrated that grid cells in the brain are encoding physical locations using hexagonally spaced, periodic phase-space representations. We explore how such a representation may be computationally advantageous for related engineering applications. Theories of how the brain decodes from a phase-space representation have been developed based on neuroscience data. However, theories of how sensory information is encoded into this phase space are less certain. Here we show a method for how a navigation-relevant input space such as elevation trajectories may be mapped into a phase-space coordinate system that can be decoded using previously developed theories. We also consider how such an algorithm may then also be mapped onto neuromrophic systems. Just as animals can tell where they are in a local region based on where they have been, our encoding algorithm enables the localization to a position in space by integrating measurements from a trajectory over a map. In this paper, we walk through our approach with simulations using a digital elevation model.

It has been demonstrated that grid cells in the brain are encoding physical locations using hexagonally spaced, periodic phase-space representations. We explore how such a representation may be computationally advantageous for related engineering applications. Theories of how the brain decodes from a phase-space representation have been developed based on neuroscience data. However, theories of how sensory information is encoded into this phase space are less certain. Here we show a method for how a navigation-relevant input space such as elevation trajectories may be mapped into a phase-space coordinate system that can be decoded using previously developed theories. We also consider how such an algorithm may then also be mapped onto neuromrophic systems. Just as animals can tell where they are in a local region based on where they have been, our encoding algorithm enables the localization to a position in space by integrating measurements from a trajectory over a map. In this paper, we walk through our approach with simulations using a digital elevation model.

Dirac semimetals have attracted a great deal of current interests due to their potential applications in topological quantum computing, low-energy electronic devices, and single photon detection in the microwave frequency range. Herein are results from analyzing the low magnetic (B) field weak-antilocalization behaviors in a Dirac semimetal Cd₃As₂ thin flake device. At high temperatures, the phase coherence length lΦ first increases with decreasing temperature (T) and follows a power law dependence of l_Φ ∝ T^–0.4. Below ~3 K, l_Φ tends to saturate to a value of ~180 nm. Another fitting parameter α, which is associated with independent transport channels, displays a logarithmic temperature dependence for T > 3 K, but also tends to saturate below ~3 K. The saturation value, ~1.45, is very close to 1.5, indicating three independent electron transport channels, which we interpret as due to decoupling of both the top and bottom surfaces as well as the bulk. This result, to our knowledge, provides first evidence that the surfaces and bulk states can become decoupled in electronic transport in Dirac semimetal Cd₃As₂.

The principle of least action is the cornerstone of classical mechanics, theory of relativity, quantum mechanics, and thermodynamics. Here, we describe how a neural network (NN) learns to find the trajectory for a Lennard-Jones (LJ) system that maintains balance in minimizing the Onsager-Machlup (OM) action and maintaining the energy conservation. The phase-space trajectory thus calculated is in excellent agreement with the corresponding results from the "ground-truth" molecular dynamics (MD) simulation. Furthermore, we show that the NN can easily find structural transformation pathways for LJ clusters, for example, the basin-hopping transformation of an LJ38from an incomplete Mackay icosahedron to a truncated face-centered cubic octahedron. Unlike MD, the NN computes atomic trajectories over the entire temporal domain in one fell swoop, and the NN time step is a factor of 20 larger than the MD time step. The NN approach to OM action is quite general and can be adapted to model morphometrics in a variety of applications.

We are currently witnessing the dawn of hydrogen (H₂) economy, where H₂ will soon become a primary fuel for heating, transportation, and longdistance and long-term energy storage. Among diverse possibilities, H₂ can be stored as a pressurized gas, a cryogenic liquid, or a solid fuel via adsorption onto porous materials. Metal–organic frameworks (MOFs) have emerged as adsorbent materials with the highest theoretical H₂ storage densities on both a volumetric and gravimetric basis. However, a critical bottleneck for the use of H₂ as a transportation fuel has been the lack of densification methods capable of shaping MOFs into practical formulations while maintaining their adsorptive performance. Here, we report a high-throughput screening and deep analysis of a database of MOFs to find optimal materials, followed by the synthesis, characterization, and performance evaluation of an optimal monolithic MOF (monoMOF) for H₂ storage. After densification, this _monoMOF stores 46 g L^–1 H₂ at 50 bar and 77 K and delivers 41 and 42 g L^–1 H₂ at operating pressures of 25 and 50 bar, respectively, when deployed in a combined temperature– pressure (25–50 bar/77 K → 5 bar/160 K) swing gas delivery system. This performance represents up to an 80% reduction in the operating pressure requirements for delivering H₂ gas when compared with benchmark materials and an 83% reduction compared to compressed H₂ gas. Our findings represent a substantial step forward in the application of high-density materials for volumetric H₂ storage applications.

The structural and dynamical properties of nanoconfined solutions can differ dramatically from those of the corresponding bulk systems. Understanding the changes induced by confinement is central to controlling the behavior of synthetic nanostructured materials and predicting the characteristics of biological and geochemical systems. A key outstanding issue is how the molecular-level behavior of nanoconfined electrolyte solutions is reflected in different experimental, particularly spectroscopic, measurements. This is addressed here through molecular dynamics simulations of the OH stretching infrared (IR) spectroscopy of NaCl, NaBr, and NaI solutions in isotopically dilute HOD/D2O confined in hydroxylated amorphous silica slit pores of width 1-6 nm and pH ∼2. In addition, the water reorientation dynamics and spectral diffusion, accessible by pump-probe anisotropy and two-dimensional IR measurements, are investigated. The aim is to elucidate the effect of salt identity, confinement, and salt concentration on the vibrational spectra. It is found that the IR spectra of the electrolyte solutions are only modestly blue-shifted upon confinement in amorphous silica slit pores, with both the size of the shift and linewidth increasing with the halide size, but these effects are suppressed as the salt concentration is increased. This indicates the limitations of linear IR spectroscopy as a probe of confined water. However, the OH reorientational and spectral diffusion dynamics are significantly slowed by confinement even at the lowest concentrations. The retardation of the dynamics eases with increasing salt concentration and pore width, but it exhibits a more complex behavior as a function of halide.

Nanothermite NiO-Al is a promising material system for low gas emission heat sources; yet, its reactive properties are highly dependent on material processing conditions. In the current study, sputter deposition is used to fabricate highly controlled nanolaminates comprised of alternating NiO and Al layers. Films having an overall stoichiometry of 2Al to 3NiO were produced with different bilayer thicknesses to investigate how ignition and self-sustained, high temperature reactions vary with changes to nanometer-scale periodicity and preheat conditions. Ignition studies were carried out with both hot plate and laser irradiation and compared to slow heating studies in hot-stage x-ray diffraction. Ignition behavior has bilayer thickness and heating rate dependencies. The 2Al/3NiO with λ ≤ 300 nm ignited via solid/solid diffusion mixing (activation energy, Ea = 49 ± 3 kJ/mole). Multilayers having λ≥ 500 nm required a more favorable mixing kinetics of solid/liquid dissolution into molten Al (Ea = 30 ± 4 kJ/mole). This solid/liquid dissolution Ea is a factor of 5 lower than that of the previously reported powder compacts due to the elimination of a passivating Al oxide layer present on the powder. The reactant mixing mechanism between 300 and 500 nm bilayer thicknesses was dependent on the ignition source's heating rate. The self-propagating reaction velocities of 2Al/3NiO multilayers varied from 0.4 to 2.5 m/s. Pre-heating nanolaminates to temperatures below the onset reaction temperatures associated with forming intermediate nickel aluminides at multilayer interfaces led to increased propagation velocities, whereas pre-heating samples above the onset temperatures inhibited subsequent attempts at laser ignition.

The selective amorphization of SiGe in Si/SiGe nanostructures via a 1 MeV Si+ implant was investigated, resulting in single-crystal Si nanowires (NWs) and quantum dots (QDs) encapsulated in amorphous SiGe fins and pillars, respectively. The Si NWs and QDs are formed during high-temperature dry oxidation of single-crystal Si/SiGe heterostructure fins and pillars, during which Ge diffuses along the nanostructure sidewalls and encapsulates the Si layers. The fins and pillars were then subjected to a 3 × 1015 ions/cm2 1 MeV Si+ implant, resulting in the amorphization of SiGe, while leaving the encapsulated Si crystalline for larger, 65-nm wide NWs and QDs. Interestingly, the 26-nm diameter Si QDs amorphize, while the 28-nm wide NWs remain crystalline during the same high energy ion implant. This result suggests that the Si/SiGe pillars have a lower threshold for Si-induced amorphization compared to their Si/SiGe fin counterparts. However, Monte Carlo simulations of ion implantation into the Si/SiGe nanostructures reveal similar predicted levels of displacements per cm3. Molecular dynamics simulations suggest that the total stress magnitude in Si QDs encapsulated in crystalline SiGe is higher than the total stress magnitude in Si NWs, which may lead to greater crystalline instability in the QDs during ion implant. The potential lower amorphization threshold of QDs compared to NWs is of special importance to applications that require robust QD devices in a variety of radiation environments.

A study was conducted to investigate the outgassing characteristics of four thermal filler materials. The purpose of this screening was to identify any outgassing products that might be considered reactive, specifically compounds that could result in corrosion in the systems where these materials are used. A range of compounds was observed in the sample headspaces, though most do not stand out as being known reactive species of concern. However, several halogenated compounds and sulfurous compounds- classes compounds known to facilitate corrosion reactions under certain conditions- were observed in low concentrations. The TFLEX 760 exhibited the highest total outgassing, while the GR130 had the lowest. Therm-a-Gap75 and the Si thermal grease exhibited very similar outgassing profiles. It is difficult to predict the extent to which any given compound observed in an analysis of this type might pose a risk in an actual system; factors such as temperature, system geometry, concentration, and gas conductance all play a role in the kinetics governing chemical reactions. It is recommended that the results of these analyses are shared with pertinent materials SMEs familiar with the system(s) in question to evaluate potential risks.

Next generation ion traps will likely need to support tens if not hundreds of ions in order to achieve several logical qubits. As we scale to those sizes, the same problems we face now – rf dissipation, control I/O, and optical access – will only grow and become more complicated. While many of these challenges can potentially be solved with technology integration, independently researching the feasibility of that integration and other solutions may help reduce the time and risk in scaling up to larger traps, by testing on smaller less complex devices. We should also consider other fabrication techniques that may help scale to larger devices, such as: through-substrate-vias (TSVs), different metal coatings, exotic rf routing, on chip laser sources, or even a secondary macroscopic trap to reload ions from. To have these technologies ready for full scale integration when we need them, ion traps with some of these capabilities need to be produced now. Developing the rigorous fabrication methods for producing reliable traps takes time and experimentation. We propose developing larger ion traps and reliable integrated technology in conjunction to make both available faster.

This summary report, prepared by Condor Country Consulting Inc. (CCCI), covers all construction monitoring of the Swainson’s hawk (Buteo swainsoni, SWHA) nest that occurred from April 2022 to July 2022 at the Sandia National Laboratories, California (SNL/CA). SNL/CA personnel first observed the SWHA pair near the pine tree along East Avenue that contained last year’s nest on April 19th, 2022. SWHA are a threatened species in the state of California and require a quarter-mile non-disturbance buffer. SNL/CA received approval from the California Department of Fish and Wildlife (CDFW) to continue work within the buffer area with a biological monitor on site to observe the SWHA for construction related disturbances. Monitoring began on April 20th, 2022 and concluded on July 1st, 2022.

This report deals with the development and evaluation of a numerical model to examine applied coating to a metal substrate subjected to a ballistic impact. The numerical model will be used to examine the benefit of the coating in resisting penetration due to the impact. For a detailed examination the Retch-Ipson curve is used as a metric. The numerical data is plotted and then fit to the Retch-Ipson curve and error calculations are used to compare the difference between the numerical output and the experimental data. This initial study is an examination of a few shortcomings of the standard material models used, and demonstrate the future work that is needed to understand the ballistic behavior of materials.

Large-scale CO2 sequestration into geological formations has been suggested to reduce CO2 emissions from industrial activities. However, much like enhanced geothermal stimulation and wastewater injection, CO2 sequestration has a potential to induce earthquake along weak faults, which can be considered a negative impact on safety and public opinion. This study shows the physical mechanisms of potential seismic hazards along basement faults driven by CO2 sequestration under variation in geological and operational constraints. Specifically we compare the poroelastic behaviors between multiphase flow and single-phase flow cases, highlighting specific needs of evaluating induced seismicity associated with CO2 sequestration. In contrast to single-phase injection scenario, slower migration of the CO2 plume than pressure pulse may delay accumulation of pressure and stress along basement faults that may not be mitigated immediately by shut-in of injection. The impact of multiphase flow system, therefore, needs to be considered for proper monitoring and mitigation strategies.

Erosion and flooding impacts to Arctic coastal environments are intensifying with nearshore oceanographic conditions acting as a key environmental driver. Robust and comprehensive assessment of the nearshore oceanographic conditions require knowledge of the following boundary conditions: incident wave energy, water level, incident wind energy, ocean temperature and salinity, bathymetry, and shoreline orientation. The number of offshore oceanographic boundary conditions can be large, requiring a significant computational investment to reproduce nearshore conditions. This present study develops location-independent typologies to reduce the number of boundary conditions needed to assess nearshore oceanographic environments in both a Historical (2007–2019) and Future (2020–2040) timespan along the Alaskan North Slope. We used WAVEWATCH III® and Delft3D Flexible Mesh model output from six oceanographic sites located along a constant ∼50 m bathymetric line spanning the Chukchi to Beaufort Seas. K-means clustering was applied to the energy-weighted joint-probability distribution of significant wave height (Hs) and peak period (Tp). Distributions of wave and wind direction, wind speed, and water level associated with location-independent centroids were assigned single values to describe a reduced order, typological rendition of offshore oceanographic conditions. Reanalysis data (e.g., ASRv2, ERA5, and GOFS) grounded the historical simulations while projected conditions were obtained from downscaled GFDL-CM3 forced under RCP8.5 conditions. Location-dependence for each site is established through the occurrence joint-probability distribution in the form of unique scaling factors representing the fraction of time that the typology would occupy over a representative year. As anticipated, these typologies show increasingly energetic ocean conditions in the future. They also enable computationally efficient simulation of the nearshore oceanographic environment along the North Slope of Alaska for better characterization of coastal processes (e.g., erosion, flooding, or sediment transport).

The mesostructure of porous electrodes used in lithium-ion batteries strongly influences cell performance. Accurate imaging of the distribution of phases in these electrodes would allow this relationship to be better understood through simulation. However, imaging the nanoscale features in these components is challenging. While scanning electron microscopy is able to achieve the required resolution, it has well established difficulties imaging porous media. This is because the flat imaging planes prepared using focused ion beam milling will intersect with the pores, which makes the images hard to interpret as the inside walls of the pores are observed. It is common to infiltrate porous media with resin prior to imaging to help resolve this issue, but both the nanoscale porosity and the chemical similarity of the resins to the battery materials undermine the utility of this approach for most electrodes. In this study, a technique is demonstrated which uses in situ infiltration of platinum to fill the pores and thus enhance their contrast during imaging. Reminiscent of the Japanese art of repairing cracked ceramics with precious metals, this technique is referred to as the kintsugi method. The images resulting from applying this technique to a conventional porous cathode are presented and then segmented using a multi-channel convolutional method. We show that while some cracks in active material particles were empty, others appear to be filled (perhaps with the carbon binder phase), which will have implications for the rate performance of the cell. Energy dispersive X-ray spectroscopy was used to validate the distribution of phases resulting from image analysis, which also suggested a graded distribution of the binder relative to the carbon additive. The equipment required to use the kintsugi method is commonly available in major research facilities and so we hope that this method will be rapidly adopted to improve the imaging of electrode materials and porous media in general.

Understanding the potential risk of stress corrosion cracking of spent nuclear fuel dry storage canisters has been identified as a knowledge gap for determining the safety of long-term interim storage of spent nuclear fuel. To address this, the DOE is funding a multi-lab DOE effort to understand the timing, occurrence, and consequences of potential canister SCC. Sandia National Laboratories has developed a probabilistic model for canister penetration by SCC. This model has been continuously updated at SNL since 2014. Model uncertainties are treated using a nested loop structure, where the outer loop accounts for uncertainties due to lack of data and the inner aleatoric loop accounts for uncertainties due to variation in nature. By separating uncertainties into these categories, it is possible to focus future work on reducing the most influential epistemic uncertainties. Several experimental studies have already been performed to improve the modeling approach through expanded process understanding and improved model parameterization. The resulting code is physics-based and intended to inform future work by identifying (1) important modeling assumptions, (2) experimental data needs, and (3) necessary model developments. In this document, several of the sub-models in the probabilistic SCC model have been exercised, and the intermediate results, as the model progresses from one sub-model to the next, are presented. Evaluating the sub-models in this manner provides a better understanding of sub-model outputs and has identified several unintended consequences of model assumptions or parameterizations, requiring updates to the modeling approach. The following updates have been made, and future updates have been identified.

Accurate diagnosis of areal density ($\textit{ρR}$) is critical for the inference of performance metrics in inertial confinement fusion implosions. One potential source of error in this diagnosis is the existence of low mode perturbations in the imploding target, which lead to asymmetries in the inference of the $\textit{ρR}$ from different lines of sight. Here, the error accrued as a result of limited coverage of the sphere due to a finite number of detectors is quantified, and the development of a forward scatter measurement from the OMEGA neutron time-of-flight detectors is motivated. In this work, a method by which the 1D-equivalent 4 π-averaged $\langleρR\rangle$ can be reconstructed, if accurate mode information can be diagnosed by other means, is validated.

Multi-phase problems have so many more unknowns, we’d like to have a tool to constrain some open questions related to microstructure and twin & dislocation behavior. We want an atomistic scale perspective on aspects of strength. Some multi-scale questions accessible to atomistic study: What lattice-specific behavior influences dislocation production/mobility and/or twinning? Do the phase transformations wipe-out, modify or preserve grain size and orientation? Does plastic strain reset at phase transition? If so under what conditions? Tin is the material chosen for the effort because it is non-hazardous and has multiple accessible solid phases at relatively low pressures.

Infrasound is generated by a variety of natural and anthropogenic sources. Infrasonic waves travel through the dynamic atmosphere, which can change on the order of minutes to hours. Infrasound propagation largely depends on the wind and temperature structure of the atmosphere. Numerical weather prediction models are available to provide atmospheric specifications, but uncertainties in these models exist and they are computationally expensive to run. Machine learning has proven useful in predicting tropospheric weather using Long Short-Term Memory (LSTM) networks. An LSTM network is utilized to make atmospheric specification predictions up to ~30 km for three different training and testing scenarios: (a) the model is trained and tested using only radiosonde data from the Albuquerque, NM, USA station, (b) the model is trained on radiosonde stations across the contiguous US, excluding the Albuquerque, NM, USA station, which was reserved for testing, and (c) the model is trained and tested on radiosonde stations across the contiguous US. Long Short-Term Memory predictions are compared to a state-of-the-art reanalysis model and show cases where the LSTM outperforms, performs equally as well, or underperforms in comparison to the state-of-the-art. Regional and temporal trends in model performance across the US are also discussed. Results suggest that the LSTM model is a viable tool for predicting atmospheric specifications for infrasound propagation modeling.

We report peridynamics is a relatively new field in continuum mechanics that has developed over the past 20 years. This paper studies the evolution of collaborations in the field of peridynamics since its inception using social network analysis. For this purpose, we construct a network for each year from 2000 to 2019 describing co-authorship between scientists in peridynamics. In these networks, each node represents a scientist and each link connects two co-authoring scientists with a link weight representing the frequency and strength of their collaboration; each network as a whole can be thought of as a graph representation of the peridynamics community for the given year. By constructing a network for each year, we are able to analyze the evolution of the network in time and discuss the implications of this evolution for the peridynamics community. Our study demonstrates that the peridynamics community has been growing exponentially in size in recent years. Centrality metrics are also used to identify the most collaborative scientists in the community. Moreover, we compute link recommendations based on both elevating a scientist’s position in the network with respect to certain centrality metrics or closing structural holes in the network identified with persistent homology. We further extend the analysis to higher-order networks whose nodes represent groups of scientists in the community and whose links connect collaborating groups. In some sense, our work studies the past, present, and future of the peridynamics community.

Some past and future supercomputer nodes incorporate High- Bandwidth Memory (HBM). Compared to standard DRAM, HBM has similar latency, higher bandwidth and lower capacity. In this paper, we evaluate algorithms for managing High- Bandwidth Memory automatically. Previous work suggests that, in the worst case, performance is extremely sensitive to the policy for managing the channel to DRAM. Prior theory shows that a priority-based scheme (where there is a static strict priority-order among p threads for channel access) is O(1)-competitive, but FIFO is not, and in the worst case is ?(p) competitive. Following this theoretical guidance would be a disruptive change for vendors, who currently use FIFO variants in their DRAMcontroller hardware. Our goal is to determine theoretically and empirically whether we can justify recommending investment in priority-based DRAM controller hardware. In order to experiment with DRAM channel protocols, we chose a theoretical model, validated it against real hardware, and implemented a basic simulator. We corroborated the previous theoretical results for the model, conducted a parameter sweep while running our simulator on address traces from memory bandwidth-bound codes (GNU sort and TACO sparse matrix-vector product), and designed better channel-access algorithms.

Symbolic regression (SR) with a multi-gene genetic program has been used to elucidate new empirical equations describing diffusion in Lennard-Jones (LJ) fluids. Examples include equations to predict self-diffusion in pure LJ fluids and equations describing the finite-size correction for self-diffusion in binary LJ fluids. The performance of the SR-obtained equations was compared to that of both the existing empirical equations in the literature and to the results from artificial neural net (ANN) models recently reported. It is found that the SR equations have improved predictive performance in comparison to the existing empirical equations, even though employing a smaller number of adjustable parameters, but show an overall reduced performance in comparison to more extensive ANNs.

A semi-analytic fluid model has been developed for characterizing relativistic electron emission across a warm diode gap. Here we demonstrate the use of this model in (i) verifying multi-fluid codes in modeling compressible relativistic electron flows (the EMPIRE-Fluid code is used as an example; see also Ref. 1), (ii) elucidating key physics mechanisms characterizing the influence of compressibility and relativistic injection speed of the electron flow, and (iii) characterizing the regimes over which a fluid model recovers physically reasonable solutions.

Wall impingement and fuel film deposition in gasoline direct injection engines under cold start conditions are major concerns for emissions reduction. However, it is challenging to study the dynamics of film deposition under realistic conditions because of the difficulty of measuring the thicknesses of these microscale films. Low-coherence interferometry provides a quantitative optical film thickness measurement technique that can be applied to study this problem. This work presents the first high-speed spectral low-coherence interferometry measurements of impinging gasoline direct injection sprays. The feasibility and practical concerns associated with high-speed low-coherence interferometry systems are explored. Two approaches to spectral low-coherence interferometry: Michelson interferometry and Fizeau interferometry, were implemented and are compared. The results show that Fizeau interferometry is the better option for measurements of impinging sprays in closed spray vessels. The high-speed low-coherence interferometry system was applied in the Fizeau configuration to measure time-resolved film thickness of impinging sprays under engine-relevant conditions to demonstrate its capabilities.

Security engineering approaches can often focus on a particular domain—physical security, cyber security, or personnel security, for example. Yet, security systems engineering consistently faces challenges requiring socio-technical solutions to address evolving and dynamic complexity. While some drivers of this complexity stem from complex risk environments, innovative adversaries, and disruptive technologies, other drivers are endogenous and emerge from the interactions across security engineering approaches. In response, INCOSE's Systems Security Working Group identified the need to better coordinate “disparate security solutions [that] operate independently” as one of eleven key concepts in their IS21 FuSE Security Roadmap. From this perspective, this need for “security orchestration” aligns with the perspective that security is a property that emerges from interactions within complex systems. Current efforts at Sandia National Laboratories are developing a systems security engineering approach that describes high consequence facility (HCF) security as a multidomain set of interacting layers. The result is a multilayered network (MLN)-based approach that captures the interactions between infrastructure, physical components, digital components, and humans in nuclear security systems. This article will summarize the MLN-based approach to HCF security and describe two preliminary results demonstrating potential benefits from incorporating interactions across disparate security solutions. Here, leveraging the logical structure of networks, this MLN model-based approach provides an example of how security orchestration provides enhanced systems security engineering solutions.

Abstract not provided.

The purpose of this report is to document the probable origin and cause of the 2022 Lurance Fire. As weather was mostly clear skies and lightning was not a factor, the fire had to be the result of human activity.

Spark-ignition (SI) engine efficiency can be increased by operating lean and with increased compression ratio (CR), but both of these measures tend to increase the propensity for undesirable acoustic knock generation. It is well known that increased CR makes the engine more prone to knock due to increased combustion pressures and temperatures, but it may be less well understood why lean operation would exacerbate knock generation. For typical gasoline-range fuels, the laminar flame speed becomes very low (roughly only 20% compared to stoichiometric conditions) for an air-excess ratio (λ) of 2. Indirectly, this exacerbates the knock challenge in two ways; a) it may necessitate operation with a combustion phasing near Top Dead Center (TDC) to complete the combustion before expansion cooling occurs, b) it increases cycle-to-cycle variations, making it more challenging to operate near the knock limits. In addition, the high intake pressure required for lean operation (nearly a factor of two higher for λ = 2 compared to λ = 1) increases the oxygen concentration which promotes end-gas autoignition and knock generation. Towards overcoming these challenges of lean combustion, this study aims to provide a better understanding of fuel autoignition under various conditions. First, to reveal the octane appetite under lean conditions, this experimental work utilized fuels of varying Research Octane Number (RON) and octane sensitivity (S). It was found that lean operation favored fuels that have high RON and high S since those were less knock limited. However, two compositionally different fuels with similarly high RON and S exhibited notable difference in knock limits under lean operation, indicating that RON and S may fail to accurately rank order fuels' knock propensity. Second, the experiments show that under boosted conditions end-gas autoignition becomes sensitive to the level of trapped residual nitric oxide (NO), which in turn is very sensitive to variations of both actual λ and combustion phasing, among other factors. The results suggest that strong knock-suppression benefits could be realized if single-ppm NO mole fraction can be consistently maintained in the reactants. Finally, it is noted that maintaining knock-free operation is particularly important for lean operation because the lower peak combustion temperatures lower the speed of sound, which in turn shifts the frequency content of the in-cylinder knock to a lower frequency range. Lower knock frequencies can increase the transmission efficiency from the combustion chamber to the outer surfaces of the engine, potentially increasing engine noise levels if knock occurs.

Brillouin based distributed fiber sensors present a unique set of characteristics amongst fiber sensing architectures. They are able to measure absolute strain and temperature over long distances, with high spatial resolution, and very large dynamic range in off-the-shelf fiber. However, Brillouin sensors traditionally provide only modest sensitivity due to the weak dependence of the Brillouin frequency on strain and the high signal to noise ratio required to identify the resonance’s peak frequency to within a small fraction of its linewidth. Recently, we introduced a technique which substantially improves the precision of Brillouin fiber sensors by exciting a series of lasing modes in a fiber loop cavity that experience Brillouin amplification at discrete locations in the fiber. The narrow-linewidth and high intensity of the lasing modes enabled ultra-low noise Brillouin sensors with large dynamic range. However, our initial demonstration was only modestly distributed: measuring strain at 40, non-contiguous positions along a 400 m fiber. In this work, we greatly extend this methodology to enable fully distributed sensing at 1000 contiguous locations along 3.5 km of fiber—an order of magnitude increase in sensor count and range. This highly-multiplexed Brillouin fiber laser sensor provides a strain noise as low as 34 nε/√Hz and we analyze the limiting factors in this approach.

Sandia National Laboratories has tested and evaluated a new version of the Chaparral 64S infrasound sensor, designed and manufactured by Chaparral Physics. The purpose of this infrasound sensor evaluation is to measure the performance characteristics in such areas as power consumption, sensitivity, full scale, self-noise, dynamic range, response, passband, sensitivity variation due to changes in barometric pressure and temperature, and sensitivity to acceleration. The Chaparral 64S infrasound sensors are being evaluated for use in the International Monitoring System (IMS) of the Preparatory Commission to the Comprehensive Nuclear-Test-Ban Treaty Organization (CTBTO).

Abstract not provided.

In the near future, grid operators are expected to regularly use advanced distributed energy resource (DER) functions, defined in IEEE 1547-2018, to perform a range of grid-support operations. Many of these functions adjust the active and reactive power of the device through commanded or autonomous operating modes which induce new stresses on the power electronics components. In this work, an experimental and theoretical framework is introduced which couples laboratory-measured component stress with advanced inverter functionality and derives a reduction in useful lifetime based on an applicable reliability model. Multiple DER devices were instrumented to calculate the additional component stress under multiple reactive power setpoints to estimate associated DER lifetime reductions. A clear increase in switch loss was demonstrated as a function of irradiance level and power factor. This is replicated in the system-level efficiency measurements, although magnitudes were different—suggesting other loss mechanisms exist. Using an approximate Arrhenius thermal model for the switches, the experimental data indicate a lifetime reduction of 1.5% when operating the inverter at 0.85 PF—compared to unity PF—assuming the DER failure mechanism thermally driven within the H-bridge. If other failure mechanisms are discovered for a set of power electronics devices, this testing and calculation framework can easily be tailored to those failure mechanisms.

Abstract not provided.

Abstract not provided.

Sandia National Laboratories has tested and evaluated an updated version of the MB3a infrasound sensor, designed by CEA and manufactured by SeismoWave. The purpose of this infrasound sensor evaluation is to measure the performance characteristics in such areas as power consumption, sensitivity, full scale, self-noise, dynamic range, response, passband, sensitivity variation due to changes in barometric pressure and temperature, and sensitivity to acceleration. The MB3a infrasound sensors are being evaluated for use in the International Monitoring System (IMS) of the Preparatory Commission to the Comprehensive Nuclear-Test-Ban Treaty Organization (CTBTO).

Abstract not provided.

Abstract not provided.

Targeted doping of grain boundaries is widely pursued as a pathway for combating thermal instabilities in nanocrystalline metals. However, certain dopants predicted to produce grain-boundary-segregated nanocrystalline configurations instead form small nanoprecipitates at elevated temperatures that act to kinetically inhibit grain growth. Here, thermodynamic modeling is implemented to select the Mo–Au system for exploring the interplay between thermodynamic and kinetic contributions to nanostructure stability. Using nanoscale multilayers and in situ transmission electron microscopy thermal aging, evolving segregation states and the corresponding phase transitions are mapped with temperature. The microstructure is shown to evolve through a transformation at lower homologous temperatures (<600 °C) where solute atoms cluster and segregate to the grain boundaries, consistent with predictions from thermodynamic models. An increase in temperature to 800 °C is accompanied by coarsening of the grain structure via grain boundary migration but with multiple pinning events uncovered between migrating segments of the grain boundary and local solute clustering. Direct comparison between the thermodynamic predictions and experimental observations of microstructure evolution thus demonstrates a transition from thermodynamically preferred to kinetically inhibited nanocrystalline stability and provides a general framework for decoupling contributions to complex stability transitions while simultaneously targeting a dominant thermal stability regime.

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.

Sandia National Laboratories performed tests to address the potential vulnerability concerns of a coupled High-Altitude Electromagnetic Pulse (HEMP) inducing secondary coupling onto critical instrumentation and control cables in a nuclear power plant, with specific focus on early-time HEMP. Three types of receiving cables in nine configurations were tested to determine transfer functions between two electrically separated cables referenced to the common mode input current on the transmitting cable. One type of transfer function related the input short circuit current and resulting open circuit voltage on the receiving cable. The other transfer function related the input short circuit current and the resulting short circuit current on the receiving cable. A 500 A standard HEMP waveform was input into the transfer functions to calculate peak coupling values on the receiving cables. The highest level of coupling using the standard waveform occurred when cables were in direct contact, with a peak short circuit current of 85 A and open circuit voltage of 9.8 kV, while configurations with separated cables predicted coupling levels of less than 5 A or 500 V.

Abstract not provided.

Abstract not provided.

Abstract not provided.

A collection of x-ray computed tomography scans of specimens from the Museum of Southwestern Biology.