Salt formations have long been recognized as a highly favorable host rock for the final disposal of high-level radioactive waste (HLW) in deep geological repositories. Their unique properties, including exceptional impermeability, self-healing capabilities, and thermal conductivity, make them a reliable natural barrier for the deep disposal of radioactive waste. This report focuses on the development and application of a methodology for assessing the integrity and per formance of the Engineered Barrier System (EBS) within salt-based repositories, a critical component of the multi-barrier system ensuring safe radioactive waste disposal.
This report provides a comprehensive assessment of physical security modeling and simulation tools available for use in the vulnerability assessment (VA) process for nuclear facilities. It outlines the historical evolution of VA methodologies, emphasizing the transition from traditional layer-based approaches to a more holistic framework that integrates detection probabilities directly into combat simulations. The document details the critical components of the VA process, including the characterization of targets, threats, and protective measures, as well as the development of adversary scenarios that reflect both insider and outsider threats. It highlights the importance of performance assurance programs, emphasizing the need for continuous evaluation and testing of security systems to ensure their effectiveness against evolving threats. Additionally, the report discusses the significance of utilizing accredited modeling and simulation tools in accredited areas to accurately represent adversary actions and the corresponding responses of protective forces. By establishing a systematic approach to VA, this document aims to enhance the overall security posture of nuclear facilities, ensuring compliance with regulatory standards while effectively mitigating risks associated with potential adversarial actions.
Under direction from the DOE Office of Electricity, Sandia National Laboratories performed testing of the Viziv system to evaluate the quality of the Zenneck surface wave and potential application to long range power transfer. This report documents the test methodology as well as the test results. This includes an analysis of prior test data collected by Viziv.
This report documents a cross-site effort to identify and evaluate materials that are suitable replacements for the polyurethanes known as EN-7 and EN-8. EN-7 and EN-8 are commercial formulations that contain free toluene diisocyanate (TDI) and which are used as adhesives and as encapsulants in many weapon applications. TDI is an OSHA-regulated volatile diisocyanate and has been targeted for elimination from future weapons use, prompting the need for a TDI-free replacement material.
The purpose of this protocol is to define procedures and practices to be used by the PACT center for field testing of metal halide perovskite (MHP) photovoltaic (PV) modules. The protocol defines the physical, electrical, and analytical configuration of the tests and applies equally to mounting systems at a fixed orientation or sun tracking systems. While standards exist for outdoor testing of conventional PV modules, these do not anticipate the unique electrical behavior of perovskite cells. Further, the existing standards are oriented toward mature, relatively stable products with lifetimes that can be measured on the scale of years to decades. The state of the art for MHP modules is still immature with considerable sample to sample variation among nominally identical modules. Version 0.0 of this protocol does not define a minimum test duration, although the intent is for modules to be fielded for periods ranging for weeks to months. This protocol draws from relevant parts of existing standards, and where necessary includes modifications specific to the behavior of perovskites.
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.
This report examines narrow slot maximum versus average transmission cross sections and the resulting directivity gain. Other aperture shapes at high frequencies are also briefly discussed.
Non-von Neumann computing using neuromorphic systems based on analogue synaptic and neuronal elements has emerged as a potential solution to tackle the growing need for more efficient data processing, but progress toward practical systems has been stymied due to a lack of materials and devices with the appropriate attributes. Recently, solid state electrochemical ion-insertion, also known as electrochemical random access memory (ECRAM) has emerged as a promising approach to realize the needed device characteristics. ECRAM is a three terminal device that operates by tuning electronic conductance in functional materials through solid-state electrochemical redox reactions. This mechanism can be considered as a gate-controlled bulk modulation of dopants and/or phases in the channel. Early work demonstrating that ECRAM can achieve nearly ideal analogue synaptic characteristics has sparked tremendous interest in this approach. More recently, the realization that electrochemical ion insertion can be used to tune the electronic properties of many types of materials including transition metal oxides, layered two-dimensional materials, organic and coordination polymers, and that the changes in conductance can span orders of magnitude has further attracted interest in ECRAM as the basis for analogue synaptic elements for inference accelerators as well as for dynamical devices that can emulate a wide range of neuronal characteristics for implementation in analogue spiking neural networks. At its core, ECRAM shares many fundamental aspects with rechargeable batteries, where ion insertion materials are used extensively for their ability to reversibly store charge and energy. Computing applications, however, present drastically different requirements: systems will require many millions of devices, scaled down to tens of nanometers, all while achieving reliable electronic-state tuning at scaled-up rates and endurances, and with minimal energy dissipation and noise. In this review, we discuss the history, basic concepts, recent progress, as well as the challenges and opportunities for different types of ECRAM, broadly grouped by their primary mobile ionic charge carrier, including Li, protons, and oxygen vacancies.
A 2D phononic pseudocrystal isolator exhibiting cyclic symmetry and radial self-similarity is measured and demonstrated to block a wide range of ultrasonic vibration. Measurements of longitudinal and shear wave blocking effects are made and compared with computational results. The use of the bandgap edge ratio is recommended for quantifying suppression in very-wide-bandgap materials. The upper-to-lower suppression edge frequency ratios of 3-4 are remarkably large for shear waves and even larger for longitudinal waves upper-to-lower suppression ratio (13 at 5 dB), such that 92.5% of frequencies in that range experience ≥ 5 dB of suppression.
Permafrost thaw increases the bioavailability of ancient organic matter, facilitating microbial metabolism of volatile organic compounds (VOCs), carbon dioxide, and methane (CH4). The formation of thermokarst (thaw) lakes in icy, organic-rich Yedoma permafrost leads to high CH4 emissions, and subsurface microbes that have the potential to be biogeochemical drivers of organic carbon turnover in these systems. However, to better characterize and quantify rates of permafrost changes, methods that further clarify the relationship between subsurface biogeochemical processes and microbial dynamics are needed. In this study, we investigated four sites (two well-drained thermokarst mounds, a drained thermokarst lake, and the terrestrial margin of a recently formed thermokarst lake) to determine whether biogenic VOCs (1) can be effectively collected during winter, and (2) whether winter sampling provides more biologically significant VOCs correlated with subsurface microbial metabolic potential. During the cold season (March 2023), we drilled boreholes at the four sites and collected cores to simultaneously characterize microbial populations and captured VOCs. VOC analysis of these sites revealed “fingerprints” that were distinct and unique to each site. Total VOCs from the boreholes included > 400 unique VOC features, including > 40 potentially biogenic VOCs related to microbial metabolism. Subsurface microbial community composition was distinct across sites; for example, methanogenic archaea were far more abundant at the thermokarst site characterized by high annual CH4 emissions. The results obtained from this method strongly suggest that ∼10% of VOCs are potentially biogenic, and that biogenic VOCs can be mapped to subsurface microbial metabolisms. By better revealing the relationship between subsurface biogeochemical processes and microbial dynamics, this work advances our ability to monitor and predict subsurface carbon turnover in Arctic soils.
BiVO4 photoanodes are promising for solar water splitting, with photogenerated electrons and holes preferentially reacting at top {010} and lateral {110} facets, respectively. However, the mechanisms driving this facet-dependent reactivity remain unclear. Here, we investigate facet-dependent photocurrent and material heterogeneity using correlative scanning photoelectrochemical microscopy (SPCM), electron beam induced current (EBIC) mapping, and mid-IR scattering scanning near-field optical microscopy (s-SNOM). SPCM measurements of 62 BiVO4 particles confirmed higher photocurrents at lateral {110} facets compared to top {010} facets, but unexpectedly revealed variations in photocurrent among lateral facets within the same particle. Variations in lateral facet surface termination could explain the intraparticle-level reactivity heterogeneity, consistent with theoretical predictions. Nano-FTIR spectroscopy and Raman microspectroscopy indicated significant materials chemistry heterogeneity within individual particles and facets that could be attributed to variations in lattice vibration distortions that enhance the overlap between Bi 6s and O 2p orbitals. The increased orbital overlap is significant as it potentially increases hole mobility in the valence band and potentially explains the lateral facet-dependent charge separation efficiency observed in photocurrent maps. Facet-dependent electrical and EBIC measurements showed no space charge regions at interfacet junctions or metal-BiVO4 contacts under vacuum, suggesting that photogenerated holes beneath top {010} facets are unlikely to transport to lateral {110} facets to drive water/sulfite oxidation. These findings indicate the potential influence of distinct bulk properties and surface termination chemistries across different particles and facets, highlighting the importance of carefully controlling defects and surface chemistry during sample growth to optimize photocatalytic performance.
Continuum shock mixture models are reviewed and applied to determine the equations of state for five different compositions of Ni xAl y, as well as bulk Ni+Al reactive multilayers, by combining the fundamental property data for elemental nickel and aluminum. From the literature, we down-select and evaluate two analytical models for the mixture Hugoniot, i.e., the well-known method of kinetic energy averaging (KEA) and a recent model proposed by Jordan and Baer [J. Appl. Phys. 111, 083516 (2012)]. Fundamentally, the former method assumes pressure equilibrium, whereas the latter assumes a common particle velocity and mixture sound speed from compressible two-phase cavitating flows. Additionally, we construct thermodynamically complete equations of state by fitting Einstein oscillator series models for the specific heat at constant volume. Finally, the solid solution approximation is invoked for intermetallic compositions, which are not strictly physical mixtures. Overall, the KEA model provides a better fit to the available Ni xAl y and Ni+Al multilayer shock compression data; however, there are combinations of material properties where the performance of these two models is thought to be reversed. Moreover, the results of this work include the first analytical solution of Jordan–Baer that does not require numerical root finding, as well as proposed modifications to the Einstein oscillator series to incorporate some effects of local pressure–temperature equilibrium and reaction–diffusion. Future work is planned that will use these equations of state in mesoscale simulations to study shock-induced reaction in Ni+Al multilayers, and the intended application is illustrated with a brief 2D hydrocode example.
Production of low carbon gasoline-like fuels such as methanol-to-gasoline (MTG) is a promising approach to achieve rapid greenhouse gas emission reduction of the transportation sector. Despite the fact that gasoline that meets the ASTM D4814 standard for automotive spark-ignition engine fuel can be readily produced from these processes, it is unclear how the composition of MTG may affect engine performance and emissions. In this paper, a surrogate for an MTG is used to numerically study the effects of gasoline composition on knock propensity and on the sensitivity of knock to thermal and fuel stratification, to oxygen dilution and to nitric oxide from exhaust gas recirculation of residual gases. Simulations were performed in ANSYS CHEMKIN-PRO using a comprehensive chemical kinetic mechanism for gasoline surrogates, and results of the MTG surrogate were compared against those of a petroleum-based regular E10 gasoline, termed PACE-20. A premium-grade MTG fuel was also formulated by adding ethanol to the MTG surrogate, and results were compared against those of four premium-grade, gasoline-like fuels representative of future alternative gasoline formulations. Surrogates and mechanism were evaluated by comparison against experimental engine data, and the model showed high accuracy at stoichiometric conditions (mean absolute error of ignition timing equal to 1.46 crank angle degrees) but larger deviations at lean conditions (mean absolute error of ignition timing equal to 5.52 crank angle degrees). Despite the fact that the MTG surrogate has a RON 1.1 units higher than that of PACE-20, it may show higher knock propensity at medium temperature conditions due to a less intense NTC behavior. MTG autoignition was more temperature- and equivalence ratio-sensitive than that of PACE20, suggesting that MTG can benefit more from naturally-occurring thermal stratification or from induced fuel stratification of the end gas to mitigate knock intensity. The sensitivity of autoignition reactivity to oxygen dilution and to NO concentration was higher for MTG than for regular gasoline at medium loads, but the opposite trend was observed at high loads due to the effect of pressure on the low-temperature chemistry of regular gasoline. Approximately 14 %vol ethanol content was required to upgrade the octane rating of MTG from regular grade to premium grade. Adding 13.6 %vol ethanol made the fuel autoignition less sensitive to both oxygen dilution and NO content (ignition time varies approx. 17 % and 50 % less with oxygen dilution and NO addition, respectively, when adding ethanol at high engine loads).
An optimization-based design framework is proposed to tune the response of soft metamaterials involving both geometric instabilities and nonlinear viscoelastic material behavior. Designing the response of soft metmaterials to harness instabilities and undergo large, tailored configuration changes will enable advancements in soft robotics, shock and vibration mitigation, and flexible electronics. In line with the metamaterial concept, the response of these materials is governed to a large extent by the geometric and topological makeup of their small-scale features. However, the link between structure and response is less intuitive for soft metamaterials due to their reliance upon highly nonlinear responses triggered by geometric instabilities. This is further complicated by the effects of viscoelastic relaxation, which recent studies have shown to alter the emergence of instabilities in non-intuitive ways. hese effects are accounted for in our framework to achieve various design objectives, including tailored force–displacement response and maximized energy absorption from both geometric and material effects. To fully automate this process, it is essential to have a completely robust equation solver for forward problems involving instabilities and viscoelastic relaxation. We achieve this by casting the search for stable mechanical equilibrium — i.e. the forward problem — as a minimization problem and utilize a trust region algorithm to robustly handle instabilities and follow energetically-favorable equilibrium paths through critical points.
This study investigates the influence of tunneling contact resistances between carbon nanotubes (CNTs) on electron transport and electrical conductivity of macroscopic carbon nanofibers (CNFs), which profoundly impacts the performance of CNT thin film electronics, CNF electron emitters and cathodes, and energy conversion and storage devices. Utilizing a self-consistent electrical contact model coupling a transmission line model with tunneling current, we calculate the contact resistances of a plethora of CNT-CNT contacts within a CNF fiber, which consists of aligned, densely packed CNTs. A statistical analysis is conducted, using Gaussian distributions to account for variations in contact lengths, tunneling gap distances, and single CNT aspect ratios, to calculate the CNT-CNT contact resistance and the overall resistance of CNT fiber. By scaling our model to a macroscopic level, our results are in good agreement with experimental measurements. Our calculation suggests that while increasing the contact overlap length diminishes individual CNT-CNT contact resistance, it could paradoxically increase macroscopic CNT fiber resistance for a given constant CNF mass density, which is due to that fact that a larger overlap length allows more CNTs to pack along an electrical conduction path per unit length, leading to more tunneling contact junctions connected in series and thus less number of parallel conduction paths within the fiber cross section. Increasing tunneling gap distance increases both individual contact and overall fiber resistance. This research provides a simple design tool for tailoring CNT fiber electrical properties to promote real-world applications using CNTs or similar low-dimensional materials.
Underground caverns in a salt dome are promising geologic features to store hydrogen because of salt's extremely low permeability and self-healing behavior. The salt cavern storage community, however, has not fully understood the geomechanical behaviors of salt rock driven by quick operation cycles of injection–production, which may significantly impact the cost-effective storage-recovery performance of multiple caverns. Our field-scale generic model captures the impact of cyclic loading–unloading on the salt creep behavior and deformation under different cycle frequencies, operating pressure, and spatial order of operating cavern(s). This systematic simulation study indicates that the initial operation cycle and arrangement of multiple caverns play a significant role in the creep-driven loss of cavern volumes and cavern deformation. Our future study will develop a new salt constitutive model based on geomechanical tests of site-specific salt rock to probe the cyclic behaviors of salt precisely both beneath and above the dilatancy boundary, including reverse (inverse transient) creep, the Bauschinger effect, and damage-healing mechanism.
At micro- and nanoscales, momentum transfer between surfaces is influenced by various physical mechanisms, including quantum fluctuations, electromagnetic interactions, electric charges, and the dynamics of (rarefied) gases. Under non-isothermal conditions, rarefied gases give rise to thermal Knudsen forces whose magnitudes strongly depend on the gas species and surface characteristics. Knudsen forces are particularly relevant in nanotechnology, optical manipulation, and aerospace systems, where gas rarefaction occurs due to highly confined geometries, sub-micrometer length scales, and reduced particle densities. Despite their significance, predictive modeling of Knudsen forces is limited by a lack of comprehensive experimental data across diverse materials and surface morphologies. In this work, we present a highly sensitive and adaptable measurement platform capable of directly quantifying Knudsen forces using a suspended, interchangeable micro-cantilever within controlled rarefied helium and nitrogen environments. The system integrates optical fiber interferometry to precisely capture out-of-plane displacements at sub-micrometer resolution, driven by Knudsen forces. From the empirical data, we derive a robust correlation linking the magnitudes of Knudsen forces to energy accommodation coefficients, offering deeper insights into the underlying gas-surface interaction mechanisms.
Underground caverns in a salt dome are promising geologic features to store hydrogen because of salt's extremely low permeability and self-healing behavior. The salt cavern storage community, however, has not fully understood the geomechanical behaviors of salt rock driven by quick operation cycles of injection–production, which may significantly impact the cost-effective storage-recovery performance of multiple caverns. Our field-scale generic model captures the impact of cyclic loading–unloading on the salt creep behavior and deformation under different cycle frequencies, operating pressure, and spatial order of operating cavern(s). This systematic simulation study indicates that the initial operation cycle and arrangement of multiple caverns play a significant role in the creep-driven loss of cavern volumes and cavern deformation. Our future study will develop a new salt constitutive model based on geomechanical tests of site-specific salt rock to probe the cyclic behaviors of salt precisely both beneath and above the dilatancy boundary, including reverse (inverse transient) creep, the Bauschinger effect, and damage-healing mechanism.
Conventional methods for extracting rare earth metals (REMs) from mined mineral ores are inefficient, expensive, and environmentally damaging. Recent discovery of lanmodulin (LanM), a protein that coordinates REMs with high-affinity and selectivity over competing ions, provides inspiration for new REM refinement methods. Here, we used quantum mechanical (QM) methods to investigate trivalent lanthanide cation (Ln3+) interactions with coordination systems representing bulk solvent water and protein binding sites. Energy decomposition analysis (EDA) showed differences in the energetic components of Ln3+ interaction with representatives of solvent (water, H2O) and protein binding sites (acetate, CH3COO-), highlighting the importance of accurate description of electrostatics and polarization in computational modeling of REM interactions with biological and bioinspired molecules. Relative binding free energies were obtained for Ln3+ with coordination complexes originating from binding sites in PDB structures of a lanthanum binding peptide (PDB entry 7CCO) and LanM, with explicit consideration of the first hydration shell waters, according to quasi-chemical theory (QCT). Beyond the first shell, the bulk solvent environment was represented with an implicit continuum model. Ln3+ interactions with (H2O)9 and both binding site models became more favorable, moving down the periodic series. This trend was more pronounced with the protein binding site models than with water, resulting in affinity increasing with periodic number, except for the last REM, Lu3+, which bound less favorably than the preceding element, Yb3+. Using the truncated 7CCO binding site model, the magnitude and trend of the experimental Ln3+ relative binding free energies for the whole 7CCO peptide were reproduced. Conversely, the previously reported experimental data for LanM show a preference for the earlier lanthanides; this is likely due to longer-range interactions and cooperative effects, which are not represented by the reduced models. Using the truncated 7CCO binding site model, the magnitude and trend of the experimental Ln3+ relative binding free energies for the whole 7CCO peptide were reproduced. In contrast to the previously reported experimental data for LanM, the peptide preferentially binds the earlier lanthanides. This difference likely arises due to longer-range interactions and cooperative effects not represented by the peptide. Further investigation of Ln3+ interactions with whole proteins using polarizable molecular mechanics models with explicit solvent is warranted to understand the influence of longer-ranged interactions, cooperativity, and bulk solvent. Nevertheless, the present work provides new insights into Ln3+ interactions with biomolecules and presents an effective computational platform for designing specific single-site REM binding peptides more efficiently.
This work presents the thermal runaway propagation model LIM1TR (Lithium-ion Modeling with 1-D Thermal Runaway) as an efficient tool to predict different cell-to-cell thermal runaway propagation scenarios. Here, we explored the vent gas volume production and reaction duration highlighting the relationship between these parameters and thermal runaway propagation due to convection by the vented gases. Two metrics based on gas production rate and heating rate are utilized as good indicators of the start and end of thermal runaway. LIM1TR results are compared with and validated by experiments from the literature for single-cell and multicell array experiments of 5 Ah and 10 Ah cells. By accounting for intraparticle diffusion of reacting species in the electrodes, we were able to capture the general dynamics of thermal runaway propagation and estimate acceptable reaction durations compared with the experimental values. Simulation results further demonstrated that varying heating modes lead to distinct reaction durations, consistent with experimental observations. Vent gas volume predictions indicate the need to consider both full and partial oxidation of the electrolyte. The outcomes of this work are building blocks for further investigations of module-to-module propagation by vented gases through convective heat transfer.
The space-charge-limited current (SCLC) in a vacuum diode is given by the Child-Langmuir law (CLL), whose electric potential ϕ(x) ∝ (x/D)4/3, where x is the spatial coordinate across the gap and D is the gap separation distance. For a collisional diode, SCLC is given by the Mott-Gurney law (MGL) and ϕ(x) ∝ (x/D)3/2. Here, we apply a capacitance argument for SCLC and use the transit time from a recent exact solution for collisional SCLC to show that ϕ(x) ∝ (x/D)ξ for a general collisional gap, where 4/3 ≤ ξ ≤ 3/2 . Furthermore, ξ is strictly a function of νT, where ν is the collision frequency and T is the electron transit time. Using this definition of ξ, we estimate the spatial dependence of the electron velocity and use the gap capacitance to derive an analytic equation for collisional SCLC that agrees within ~4.5% of the exact solution that requires solving parametrically through T. This analytic equation for general ξ asymptotically recovers the CLL as ν → 0 and the MGL as ν → ∞. As a result, matching these limits shows that ξ ≈ 1.40 and V ∝ D2ν2 at the transition from a vacuum to a collisional diode for any device condition.
Low-velocity impact of 2D woven glass fiber reinforced polymer (GFRP) and carbon fiber reinforced polymer (CFRP) composite laminates was studied experimentally and numerically. Hybrid laminates containing blocked layers of GFRP/CFRP/GFRP with all plies oriented at 0° were investigated. Relatively high impact energies were used to obtain full perforation of the laminate in a low-velocity impact setup. Numerical simulations were carried out using the in-house transient dynamics finite element code, Sierra/SM, developed at Sandia National Laboratories. A three-dimensional continuum damage model was used to describe the response of a woven composite ply. Two methods for handling delamination were considered and compared: (1) cohesive zone modeling and (2) continuum damage mechanics. The reduced model size achieved by omission of the cohesive zone elements produced acceptable results at reduced computational cost. The comparison between different modeling techniques can be used to inform modeling decisions relevant to low velocity impact scenarios. The modeling was validated by comparing with the experimental results and showed good agreement in terms of predicted damage mechanisms and impactor velocity and force histories.
Coastal communities face unique challenges in maintaining continuous service from critical infrastructure. This research advances capabilities for evaluating the impact of using wave energy to desalinate water on the resilience of coastal communities. The study focuses on the feasibility of using wave energy conversion to provide drinking water to communities in need and applying resilience metrics to quantify its impact on the community. To assess the feasibility of wave-powered desalination, this research couples the open-source software Wave Energy Converter SIMulator (WEC-Sim) and Water Network Tool for Resilience (WNTR). This research explores variations in both the wave resource (location, seasonality, and duration) and the ability to maintain drinking water service during a disruption scenario by applying the simulation framework to three case studies, which are based on communities in Puerto Rico. The simulation framework provides a contextualized assessment of the ability of wave-powered desalination to improve the resilience of coastal communities, which can serve as a methodology for future studies seeking the integration of wave-powered desalination with water distribution systems.
The (a)-type screw dislocations are known to be significant mediators of plasticity in hexagonal-close-packed (HCP) metals. These dislocations have polymorphic core structures, and subtle changes in the relative energies of these core structures are known to have a large impact on the dynamics of the dislocations. This work identifies a previously neglected long-range elastic interstitial-solute/dislocation interaction that influences the core structures. Essentially, interstitial solutes induce a change in the dislocation core structure to minimize the energy of interaction between the solutes and the dislocation. Molecular dynamics simulations, continuum linear elasticity, and statistical analysis show that this long-range interaction can locally alter the dislocation cores so that many different polymorphs appear along a single dislocation not only because of direct contact between interstitials and the dislocation core but also because of this long-range elastic interaction.
