Publications

This paper explores the concept of predictive maturity for non-linear concrete constitutive models employed in the computational prediction of the structural response of reinforced concrete structures to impact from free-flying missiles. Such concrete constitutive models are widely varied in complexity. Three constitutive models were utilized within the same finite element structural model to simulate the response of the IRIS III experiment. Each of the models were individually calibrated with available material testing data and also re-calibrated assuming limited availability of test data. When full calibration is possible, more sophisticated constitutive models appear to provide more predictive maturity; however, when this data is not available (e.g. for an existing structure where representative test specimens may not be available), the expected maturity is reduced. Indeed, this hypothesis is supported by the simulations that indicate good agreement with measured experimental response quantities from the IRIS III tests with complex constitutive models and full calibration, and accordingly poor predictions when less complex models are used or when the more sophisticated models are poorly calibrated. Thus, predictions of structural response where complete material testing data is not obtainable should be understood as less predictive.

On top of machine learning (ML) models, uncertainty quantification (UQ) functions as an essential layer of safety assurance that could lead to more principled decision making by enabling sound risk assessment and management. The safety and reliability improvement of ML models empowered by UQ has the potential to significantly facilitate the broad adoption of ML solutions in high-stakes decision settings, such as healthcare, manufacturing, and aviation, to name a few. In this tutorial, we aim to provide a holistic lens on emerging UQ methods for ML models with a particular focus on neural networks and the applications of these UQ methods in tackling engineering design as well as prognostics and health management problems. Towards this goal, we start with a comprehensive classification of uncertainty types, sources, and causes pertaining to UQ of ML models. Next, we provide a tutorial-style description of several state-of-the-art UQ methods: Gaussian process regression, Bayesian neural network, neural network ensemble, and deterministic UQ methods focusing on spectral-normalized neural Gaussian process. Established upon the mathematical formulations, we subsequently examine the soundness of these UQ methods quantitatively and qualitatively (by a toy regression example) to examine their strengths and shortcomings from different dimensions. Then, we review quantitative metrics commonly used to assess the quality of predictive uncertainty in classification and regression problems. Afterward, we discuss the increasingly important role of UQ of ML models in solving challenging problems in engineering design and health prognostics. Two case studies with source codes available on GitHub are used to demonstrate these UQ methods and compare their performance in the life prediction of lithium-ion batteries at the early stage (case study 1) and the remaining useful life prediction of turbofan engines (case study 2).

A key challenge in inverse problems is the selection of sensors to gather the most effective data. In this paper, we consider the problem of inferring the initial condition to a linear dynamical system and develop an efficient control-theoretical approach for greedily selecting sensors. Our method employs a Galerkin projection to reduce the size of the inverse problem, resulting in a computationally efficient algorithm for sensor selection. As a byproduct of our algorithm, we obtain a preconditioner for the inverse problem that enables the rapid recovery of the initial condition. We analyze the theoretical performance of our greedy sensor selection algorithm as well as the performance of the associated preconditioner. Finally, we verify our theoretical results on various inverse problems involving partial differential equations.

Drinking water infrastructure in urban settings is increasingly affected by population growth and disruptions like extreme weather events. This study explores how the integration of direct wastewater reuse can help to maintain drinking water service when the system is compromised.

A major hurdle in utilizing carbon dioxide (CO2) lies in separating it from industrial flue gas mixtures and finding suitable storage methods that enable its application in various industries. To address this issue, we utilized a combination of molecular dynamics simulations and experiments to investigate the behavior of CO2 in common room-temperature ionic liquids (RTIL) when in contact with aqueous interfaces. Our investigation of RTILs, [EMIM][TFSI] and [OMIM][TFSI], and their interaction with a pure water layer mimics the environment of a previously developed ultrathin enzymatic liquid membrane for CO2 separation. We analyzed diffusion constants and viscosity, which reveals that CO2 molecules exhibit faster mobility within the selected ILs compared to what would be predicted solely based on the viscosity of the liquids using the standard Einstein-Stokes relation. Moreover, we calculated the free energy of translocation for various species across the aqueous-IL interface, including CO2 and HCO3-. Free energy profiles demonstrate that CO2 exhibits a more favorable partitioning behavior in the RTILs compared to that in pure water, while a significant barrier hinders the movement of HCO3- from the aqueous layer. Experimental measurement of the CO2 transport in the RTILs corroborates the model. These findings strongly suggest that hydrophobic RTILs could serve as a promising option for selectively transporting CO2 from aqueous media and concentrating it as a preliminary step toward storage.

This article characterises the effects of cathode photoemission leading to electrical discharges in an argon gas. We perform breakdown experiments under pulsed laser illumination of a flat cathode and observe Townsend to glow discharge transitions. The breakdown process is recorded by high-speed imaging, and time-dependent voltage and current across the electrode gap are measured for different reduced electric fields and laser intensities. We employ a 0D transient discharge model to interpret the experimental measurements. The fitted values of transferred photoelectron charge are compared with calculations from a quantum model of photoemission. The breakdown voltage is found to be lower with photoemission than without. When the applied voltage is insufficient for ion-induced secondary electron emission to sustain the plasma, laser driven photoemission can still create a breakdown where a sheath (i.e. a region near the electrode surfaces consisting of positive ions and neutrals) is formed. This photoemission induced plasma persists and decays on a much longer time scale ( ∼ 10 s μ s) than the laser pulse length ( 30 ps). The effects of different applied voltages and laser energies on the breakdown voltage and current waveforms are investigated. The discharge model can accurately predict the measured breakdown voltage curves, despite the existence of discrepancy in quantitatively describing the transient discharge current and voltage waveforms.

Here this study investigates the nonlinear frequency response of a shaft-bearing assembly with vibro-impacts occurring at the bearing clearances. The formation of nonlinear behavior as system parameters change is examined, along with the effects of asymmetries in the nominal, inherently symmetric system. The primary effect of increasing the forcing magnitude or decreasing the contact gap sizes is the formation of grazing-induced chaotic solution branches occurring over a wide frequency range near each system resonance. The system's nominal setup has very hard contact stiffness and shows no evidence of isolas or superharmonic resonances over the frequency ranges of interest. Moderate contact stiffnesses cause symmetry breaking and introduce superharmonic resonance branches of primary resonances. Even if some primary resonances are not present due to the system's inherent symmetry, their superharmonic resonances still manifest. Branches of quasiperiodic isolas (isolated resonance branches) are also discovered, along with a cloud of isolas near a high-frequency resonance. Parameter asymmetries are found to produce a few significant changes in behavior: asymmetric linear stiffness, contact stiffness, and gap size could affect the behavior of primary resonant frequencies and isolas.

Aqueous electrolytes composed of 0.1 M zinc bis-(trifluoromethyl-sulfonyl)-imide (Zn-(TFSI)₂) and acetonitrile (ACN) were studied using combined experimental and simulation techniques. The electrolyte was found to be electrochemically stable when the ACN V% is higher than 74.4. In addition, it was found that the ionic conductivity of the mixed solvent electrolytes changes as a function of ACN composition, and a maximum was observed at 91.7 V% of ACN although the salt concentration is the same. This behavior was qualitatively reproduced by molecular dynamics (MD) simulations. Detailed analyses based on experiments and MD simulations show that at high ACN composition the water network existing in the high water composition solutions breaks. As a result, the screening effect of the solvent weakens and the correlation among ions increases, which causes a decrease in ionic conductivity at high ACN V%. Furthermore, this study provides a fundamental understanding of this complex mixed solvent electrolyte system.

A two-step solar thermochemical looping cycle based on Co₃Mo₃N/Co₆Mo₆N reduction/nitridation reactions offers a pathway for green NH₃ production that utilizes concentrated solar irradiation, H₂O, and air as feedstocks. The NH₃ production cycle steps both derive process heat from concentrated solar irradiation and encompass 1) the reduction of Co₃Mo₃N in H₂ to Co₆Mo₆N and NH₃; and 2) nitridation of Co₆Mo₆N to Co₃Mo₃N with N₂. Co₃Mo₃N reduction/nitridation reactions are examined at different H₂ and/or N₂ partial pressures and temperatures. NH₃ production is quantified in situ using liquid conductivity measurements coupled with mass spectrometry (MS). Solid-state characterization is performed to identify a surface oxygen layer that necessitates the addition of H₂ during cycling to prevent surface oxidation by trace amounts of O₂. H₂ concentrations of > 5% H₂/Ar and temperatures >500 °C are required to reduce Co₃Mo₃N to Co₆Mo₆N and form NH₃ at 1 bar. Complete regeneration of Co₃Mo₃N from Co₆Mo₆N is achieved at conditions of 700 °C under 25–75% H₂/N₂. H₂ pressure-swings are observed to increase NH₃ production during Co₃Mo₃N reduction. In conclusion, the results represent the first comprehensive characterization of and definitive non-catalytic production of NH₃ via chemical looping with metal nitrides and provide insights for technology development.

Before residential photovoltaic (PV) systems are interconnected with the grid, various planning and impact studies are conducted on detailed models of the system to ensure safety and reliability are maintained. However, these model-based analyses can be time-consuming and error-prone, representing a potential bottleneck as the pace of PV installations accelerates. Data-driven tools and analyses provide an alternate pathway to supplement or replace their model-based counterparts. In this article, a data-driven algorithm is presented for assessing the thermal limitations of PV interconnections. Using input data from residential smart meters, and without any grid models or topology information, the algorithm can determine the nameplate capacity of the service transformer supplying those customers. The algorithm was tested on multiple datasets and predicted service transformer capacity with >98% accuracy, regardless of existing PV installations. This algorithm has various applications from model-free thermal impact analysis for hosting capacity studies to error detection and calibration of existing grid models.

Access to accurate solar resource data is critical for numerous applications, including estimating the yield of solar energy systems, developing radiation models, and validating irradiance datasets. However, lack of standardization in data formats and access interfaces across providers constitutes a major barrier to entry for new users. pvlib python's iotools subpackage aims to solve this issue by providing standardized Python functions for reading local files and retrieving data from external providers. All functions follow a uniform pattern and return convenient data outputs, allowing users to seamlessly switch between data providers and explore alternative datasets. The pvlib package is community-developed on GitHub: https://github.com/pvlib/pvlib-python. As of pvlib python version 0.9.5, the iotools subpackage supports 12 different datasets, including ground measurement, reanalysis, and satellite-derived irradiance data. The supported ground measurement networks include the Baseline Surface Radiation Network (BSRN), NREL MIDC, SRML, SOLRAD, SURFRAD, and the US Climate Reference Network (CRN). Additionally, satellite-derived and reanalysis irradiance data from the following sources are supported: PVGIS (SARAH & ERA5), NSRDB PSM3, and CAMS Radiation Service (including McClear clear-sky irradiance).

The properties of electrons in matter are of fundamental importance. They give rise to virtually all material properties and determine the physics at play in objects ranging from semiconductor devices to the interior of giant gas planets. Modeling and simulation of such diverse applications rely primarily on density functional theory (DFT), which has become the principal method for predicting the electronic structure of matter. While DFT calculations have proven to be very useful, their computational scaling limits them to small systems. We have developed a machine learning framework for predicting the electronic structure on any length scale. It shows up to three orders of magnitude speedup on systems where DFT is tractable and, more importantly, enables predictions on scales where DFT calculations are infeasible. Our work demonstrates how machine learning circumvents a long-standing computational bottleneck and advances materials science to frontiers intractable with any current solutions.

Recent developments integrating micromechanics and neural networks offer promising paths for rapid predictions of the response of heterogeneous materials with similar accuracy as direct numerical simulations. The deep material network is one such approaches, featuring a multi-layer network and micromechanics building blocks trained on anisotropic linear elastic properties. Once trained, the network acts as a reduced-order model, which can extrapolate the material’s behavior to more general constitutive laws, including nonlinear behaviors, without the need to be retrained. However, current training methods initialize network parameters randomly, incurring inevitable training and calibration errors. Here, we introduce a way to visualize the network parameters as an analogous unit cell and use this visualization to “quilt” patches of shallower networks to initialize deeper networks for a recursive training strategy. The result is an improvement in the accuracy and calibration performance of the network and an intuitive visual representation of the network for better explainability.

Programmable photonic integrated circuits (PICs) are emerging as powerful tools for control of light, with applications in quantum information processing, optical range finding, and artificial intelligence. Low-power implementations of these PICs involve micromechanical structures driven capacitively or piezoelectrically but are often limited in modulation bandwidth by mechanical resonances and high operating voltages. Here we introduce a synchronous, micromechanically resonant design architecture for programmable PICs and a proof-of-principle 1×8 photonic switch using piezoelectric optical phase shifters. Our design purposefully exploits high-frequency mechanical resonances and optically broadband components for larger modulation responses on the order of the mechanical quality factor Q m while maintaining fast switching speeds. We experimentally show switching cycles of all 8 channels spaced by approximately 11 ns and operating at 4.6 dB average modulation enhancement. Future advances in micromechanical devices with high Qm, which can exceed 10000, should enable an improved series of low-voltage and high-speed programmable PICs.

Quantum cascade lasers (QCLs) have emerged as promising candidates for generating chip-scale frequency combs in mid-infrared and terahertz wavelengths. In this work, we demonstrate frequency comb formation in ring terahertz QCLs using the injection of light from a distributed feedback (DFB) laser. The DFB design frequency is chosen to match the modes of the ring cavity (near 3.3 THz), and light from the DFB is injected into the ring QCL via a bus waveguide. By controlling the power and frequency of the optical injection, we show that combs can be selectively formed and controlled in the ring cavity. Numerical modeling suggests that this comb is primarily frequency-modulated in character, with the injection serving to trigger comb formation. We also show that the ring can be used as a filter to control the output of the DFB QCL, potentially being of interest in terahertz photonic integrated circuits. Our work demonstrates that waveguide couplers are a compelling approach for injecting and extracting radiation from ring terahertz combs and offer exciting possibilities for the generation of new comb states in terahertz, such as frequency-modulated waves, solitons, and more.

Future machine learning strategies for materials process optimization will likely replace human capital-intensive artisan research with autonomous and/or accelerated approaches. Such automation enables accelerated multimodal characterization that simultaneously minimizes human errors, lowers costs, enhances statistical sampling, and allows scientists to allocate their time to critical thinking instead of repetitive manual tasks. Previous acceleration efforts to synthesize and evaluate materials have often employed elaborate robotic self-driving laboratories or used specialized strategies that are difficult to generalize. Herein we describe an implemented workflow for accelerating the multimodal characterization of a combinatorial set of 915 electroplated Ni and Ni–Fe thin films resulting in a data cube with over 160,000 individual data files. Our acceleration strategies do not require manufacturing-scale resources and are thus amenable to typical materials research facilities in academic, government, or commercial laboratories. The workflow demonstrated the acceleration of six characterization modalities: optical microscopy, laser profilometry, X-ray diffraction, X-ray fluorescence, nanoindentation, and tribological (friction and wear) testing, each with speedup factors ranging from 13–46x. In addition, automated data upload to a repository using FAIR data principles was accelerated by 64x.

Efficient conversion of pentose sugars remains a significant barrier to the replacement of petroleum-derived chemicals with plant biomass-derived bioproducts. While the oleaginous yeast Rhodosporidium toruloides (also known as Rhodotorula toruloides) has a relatively robust native metabolism of pentose sugars compared to other wild yeasts, faster assimilation of those sugars will be required for industrial utilization of pentoses. To increase the rate of pentose assimilation in R. toruloides, we leveraged previously reported high-throughput fitness data to identify potential regulators of pentose catabolism. Two genes were selected for further investigation, a putative transcription factor (RTO4_12978, Pnt1) and a homolog of a glucose transceptor involved in carbon catabolite repression (RTO4_11990). Overexpression of Pnt1 increased the specific growth rate approximately twofold early in cultures on xylose and increased the maximum specific growth by 18% while decreasing accumulation of arabitol and xylitol in fast-growing cultures. Improved growth dynamics on xylose translated to a 120% increase in the overall rate of xylose conversion to fatty alcohols in batch culture. Proteomic analysis confirmed that Pnt1 is a major regulator of pentose catabolism in R. toruloides. Deletion of RTO4_11990 increased the growth rate on xylose, but did not relieve carbon catabolite repression in the presence of glucose. Carbon catabolite repression signaling networks remain poorly characterized in R. toruloides and likely comprise a different set of proteins than those mainly characterized in ascomycete fungi.

Single-molecule stretching experiments are widely utilized within the fields of physics and chemistry to characterize the mechanics of individual bonds or molecules, as well as chemical reactions. Analytic relations describing these experiments are valuable, and these relations can be obtained through the statistical thermodynamics of idealized model systems representing the experiments. Since the specific thermodynamic ensembles manifested by the experiments affect the outcome, primarily for small molecules, the stretching device must be included in the idealized model system. Though the model for the stretched molecule might be exactly solvable, including the device in the model often prevents analytic solutions. In the limit of large or small device stiffness, the isometric or isotensional ensembles can provide effective approximations, but the device effects are missing. Here a dual set of asymptotically correct statistical thermodynamic theories are applied to develop accurate approximations for the full model system that includes both the molecule and the device. The asymptotic theories are first demonstrated to be accurate using the freely jointed chain model and then using molecular dynamics calculations of a single polyethylene chain.

System provenance forensic analysis has been studied by a large body of research work. This area needs fine granularity data such as system calls along with event fields to track the dependencies of events. While prior work on security datasets has been proposed, we found a useful dataset of realistic attacks and details that are needed for high-quality provenance tracking is lacking. We created a new dataset of eleven vulnerable cases for system forensic analysis. It includes the full details of system calls including syscall parameters. Realistic attack scenarios with real software vulnerabilities and exploits are used. For each case, we created two sets of benign and adversary scenarios which are manually labeled for supervised machine-learning analysis. In addition, we present an algorithm to improve the data quality in the system provenance forensic analysis. We demonstrate the details of the dataset events and dependency analysis of our dataset cases.

Computer vision models have great potential as tools for international nuclear safeguards verification activities, but off-the-shelf models require fine-tuning through transfer learning to detect relevant objects. Because open-source examples of safeguards-relevant objects are rare, and to evaluate the potential of synthetic training data for computer vision, we present the Limbo dataset. Limbo includes both real and computer-generated images of uranium hexafluoride containers for training computer vision models. We generated these images iteratively based on results from data validation experiments that are detailed here. The findings from these experiments are applicable both for the safeguards community and the broader community of computer vision research using synthetic data.

Emerging and re-emerging viral pathogens present a unique challenge for anti-viral therapeutic development. Anti-viral approaches with high flexibility and rapid production times are essential for combating these high-pandemic risk viruses. CRISPR-Cas technologies have been extensively repurposed to treat a variety of diseases, with recent work expanding into potential applications against viral infections. However, delivery still presents a major challenge for these technologies. Lipid-coated mesoporous silica nanoparticles (LCMSNs) offer an attractive delivery vehicle for a variety of cargos due to their high biocompatibility, tractable synthesis, and amenability to chemical functionalization. Here, we report the use of LCMSNs to deliver CRISPR-Cas9 ribonucleoproteins (RNPs) that target the Niemann–Pick disease type C1 gene, an essential host factor required for entry of the high-pandemic risk pathogen Ebola virus, demonstrating an efficient reduction in viral infection. We further highlight successful in vivo delivery of the RNP-LCMSN platform to the mouse liver via systemic administration.

PyApprox is a Python-based one-stop-shop for probabilistic analysis of numerical models such as those used in the earth, environmental and engineering sciences. Easy to use and extendable tools are provided for constructing surrogates, sensitivity analysis, Bayesian inference, experimental design, and forward uncertainty quantification. The algorithms implemented represent a wide range of methods for model analysis developed over the past two decades, including recent advances in multi-fidelity approaches that use multiple model discretizations and/or simplified physics to significantly reduce the computational cost of various types of analyses. An extensive set of Benchmarks from the literature is also provided to facilitate the easy comparison of new or existing algorithms for a wide range of model analyses. This paper introduces PyApprox and its various features, and presents results demonstrating the utility of PyApprox on a benchmark problem modeling the advection of a tracer in groundwater.

Lifetime-encoded materials are particularly attractive as optical tags, however examples are rare and hindered in practical application by complex interrogation methods. Here, we demonstrate a design strategy towards multiplexed, lifetime-encoded tags via engineering intermetallic energy transfer in a family of heterometallic rare-earth metal-organic frameworks (MOFs). The MOFs are derived from a combination of a high-energy donor (Eu), a low-energy acceptor (Yb) and an optically inactive ion (Gd) with the 1,2,4,5 tetrakis(4-carboxyphenyl) benzene (TCPB) organic linker. Precise manipulation of the luminescence decay dynamics over a wide microsecond regime is achieved via control over metal distribution in these systems. Demonstration of this platform’s relevance as a tag is attained via a dynamic double encoding method that uses the braille alphabet, and by incorporation into photocurable inks patterned on glass and interrogated via digital high-speed imaging. This study reveals true orthogonality in encoding using independently variable lifetime and composition, and highlights the utility of this design strategy, combining facile synthesis and interrogation with complex optical properties.

An interpretable machine learning method, physics-informed genetic programming-based symbolic regression (P-GPSR), is integrated into a continuum thermodynamic approach to developing constitutive models. The proposed strategy for combining a thermodynamic analysis with P-GPSR is demonstrated by generating a yield function for an idealized material with voids, i.e., the Gurson yield function. First, a thermodynamic-based analysis is used to derive model requirements that are exploited in a custom P-GPSR implementation as fitness criteria or are strongly enforced in the solution. The P-GPSR implementation improved accuracy, generalizability, and training time compared to the same GPSR code without physics-informed fitness criteria. The yield function generated through the P-GPSR framework is in the form of a composite function that describes a class of materials and is characteristically more interpretable than GPSR-derived equations. The physical significance of the input functions learned by P-GPSR within the composite function is acquired from the thermodynamic analysis. Fundamental explanations of why the implemented P-GPSR capabilities improve results over a conventional GPSR algorithm are provided.

R. toruloides is an oleaginous yeast, with diverse metabolic capacities and high tolerance for inhibitory compounds abundant in plant biomass hydrolysates. While R. toruloides grows on several pentose sugars and alcohols, further engineering of the native pathway is required for efficient conversion of biomass-derived sugars to higher value bioproducts. A previous high-throughput study inferred that R. toruloides possesses a non-canonical l-arabinose and d-xylose metabolism proceeding through d-arabitol and d-ribulose. In this study, we present a combination of genetic and metabolite data that refine and extend that model. Chiral separations definitively illustrate that d-arabitol is the enantiomer that accumulates under pentose metabolism. Deletion of putative d-arabitol-2-dehydrogenase (RTO4_9990) results in > 75% conversion of d-xylose to d-arabitol, and is growth-complemented on pentoses by heterologous xylulose kinase expression. Deletion of putative d-ribulose kinase (RTO4_14368) arrests all growth on any pentose tested. Analysis of several pentose dehydrogenase mutants elucidates a complex pathway with multiple enzymes mediating multiple different reactions in differing combinations, from which we also inferred a putative l-ribulose utilization pathway. Our results suggest that we have identified enzymes responsible for the majority of pathway flux, with additional unknown enzymes providing accessory activity at multiple steps. Further biochemical characterization of the enzymes described here will enable a more complete and quantitative understanding of R. toruloides pentose metabolism. These findings add to a growing understanding of the diversity and complexity of microbial pentose metabolism.

Pseudomonads are ubiquitous bacteria with importance in medicine, soil, agriculture, and biomanufacturing. We report a novel Pseudomonas putida phage, MiCath, which is the first known phage infecting P. putida S12, a strain increasingly used as a synthetic biology chassis. MiCath was isolated from garden soil under a tomato plant using P. putida S12 as a host and was also found to infect four other P. putida strains. MiCath has a ~ 61 kbp double-stranded DNA genome which encodes 97 predicted open reading frames (ORFs); functions could only be predicted for 48 ORFs using comparative genomics. Functions include structural phage proteins, other common phage proteins (e.g., terminase), a queuosine gene cassette, a cas4 exonuclease, and an endosialidase. Restriction digestion analysis suggests the queuosine gene cassette encodes a pathway capable of modification of guanine residues. When compared to other phage genomes, MiCath shares at most 74% nucleotide identity over 2% of the genome with any sequenced phage. Overall, MiCath is a novel phage with no close relatives, encoding many unique gene products.

A reduced order, nonlocal model is proposed for the contact force between initially spherical particles under compression. The model in effect provides the normal component of the interaction force between elements in the discrete element method (DEM). It is applicable to high relative density and large stress in powder compaction. It takes into account the mutual interaction between multiple points of contact, in contrast to the usual assumption in DEM of pair interactions. The mathematical form of the model is derived from a variational formulation that leads to the momentum balance for the forces on each grain. The model is calibrated mainly using detailed three dimensional peridynamic simulations of single grains under compressive loading by rigid plates that move radially with prescribed velocity. This calibration takes into account the large deformation and fracture of the grains. The interaction model also includes terms for the unloading behavior and adhesion. As validation, the model is applied to test data on the compaction of microcrystalline cellulose bulk powder.

Interface resistance has become a significant bottleneck for solid-state batteries (SSBs). Most studies of interface resistance have focused on extrinsic mechanisms such as interface reactions and imperfect contact between electrodes and solid electrolytes. Interface potentials are an important intrinsic mechanism that is often ignored. Here, we highlight Kelvin probe force microscopy (KPFM) as a tool to image the local potential at interfaces inside SSBs, examining the existing literature and discussing challenges in interpretation. Drawing analogies with electron transport in metal/semiconductor interfaces, we showcase a formalism that predicts intrinsic ionic resistance based on the properties of the contacting phases, and we emphasize that future battery designs should start from material pairs with low intrinsic resistance. We conclude by outlining future directions in the study of interface potentials through both theory and experiment. Graphic abstract: [Figure not available: see fulltext.]

Real-time time-dependent density functional theory (TDDFT) is presently the most accurate available method for computing electronic stopping powers from first principles. However, obtaining application-relevant results often involves either costly averages over multiple calculations or ad hoc selection of a representative ion trajectory. We consider a broadly applicable, quantitative metric for evaluating and optimizing trajectories in this context. This methodology enables rigorous analysis of the failure modes of various common trajectory choices in crystalline materials. Although randomly selecting trajectories is common practice in stopping power calculations in solids, we show that nearly 30% of random trajectories in an FCC aluminum crystal will not representatively sample the material over the time and length scales feasibly simulated with TDDFT, and unrepresentative choices incur errors of up to 60%. We also show that finite-size effects depend on ion trajectory via “ouroboros” effects beyond the prevailing plasmon-based interpretation, and we propose a cost-reducing scheme to obtain converged results even when expensive core-electron contributions preclude large supercells. This work helps to mitigate poorly controlled approximations in first-principles stopping power calculations, allowing 1–2 order of magnitude cost reductions for obtaining representatively averaged and converged results.

Across many industries and engineering disciplines, systems of components are designed and deployed into their operational environments. It is the desire of the engineer to be able to predict if the component or system will survive its operational environment or if the component will fail due to mechanical stresses. One method to determine if the component will survive the operational environment is to expose the component to a simulation of the environment in a laboratory. One difficulty in executing such a test is that the component may not have the same boundary condition in both the laboratory and operational configurations. This paper presents a novel method of quantifying the error in the modal domain that occurs from the impedance difference between the laboratory test fixture and the operational configuration. The error is calculated from the projection from one mode shape space to the other, and the error is in terms of each mode of the operational configuration. The error provides insight into the effectiveness of the test fixture with respect to the ability to recreate the individual mode shapes of the operational configuration. A case study is presented to show the error in the modal projection between two configurations is a lower limit for the error that can be achieved by a laboratory test.

We report on the substantial advancement of long wavelength InAs-based interband cascade lasers (ICLs) utilizing advanced waveguides formed from hybrid cladding layers and targeting the 10-12μm wavelength region. Modifications in the hole injector have improved carrier transport in these ICLs, resulting in significantly reduced threshold voltages (Vth) as low as 3.62 V at 80 K. Consequently, much higher voltage efficiencies were observed, peaking at about 73% at 10.3μm and allowing for large output powers of more than 100 mW/facet. Also, low threshold current densities (Jth) of 8.8 A/cm2 in cw mode and 7.6 A/cm2 in pulsed mode near 10μm were observed; a result of adjustments in the GaInSb hole well composition intended to reduce the overall strain accumulation in the ICL. Furthermore, an ICL from the second wafer operating at a longer wavelength achieved a peak voltage efficiency of 57% at 11.7μm, with a peak output power of more than 27 mW/facet. This ICL went on to lase beyond 12μm in both cw and pulsed modes, representing a new milestone in long wavelength coverage for ICLs with the standard W-QW active region.

Behind-the-meter (BTM) battery energy storage systems (BESS) are undergoing rapid deployment. Simple equations to estimate the installed cost of BTM BESS are often necessary when a rigorous, bottom-up cost estimate is not available or not appropriate, in applications such as energy system modeling, informing a BESS sizing decision, and cost benchmarking. Drawing on project-level data from California, I estimate several predictive regression models of the installed cost of a BTM BESS as a function of energy capacity and power capacity. The models are evaluated for in-sample goodness-of-fit and out-of-sample predictive accuracy. The results of these analyses indicate stronger empirical support for models with natural log transformations of installed cost, energy, and power as compared against widely-used models that posit a linear relationship among the untransformed versions of these variables. Building on these results, I present a logarithmic model that can predict installed cost conditional on energy capacity, power capacity, AC or DC coupling with distributed generation, customer sector, and local wages for electricians. I document how the model can be easily extrapolated to future years, either with forecasts from other sources or by re-estimating the parameters with the latest data.

Perspectives for understanding the brain vary across disciplines and this has challenged our ability to describe the brain’s functions. In this comment, we discuss how emerging theoretical computing frameworks that bridge top-down algorithm and bottom-up physics approaches may be ideally suited for guiding the development of neural computing technologies such as neuromorphic hardware and artificial intelligence. Furthermore, we discuss how this balanced perspective may be necessary to incorporate the neurobiological details that are critical for describing the neural computational disruptions within mental health and neurological disorders.

The frequency, severity, and extent of climate extremes in future will have an impact on human well-being, ecosystems, and the effectiveness of emissions mitigation and carbon sequestration strategies. The specific objectives of this study were to downscale climate data for US weather stations and analyze future trends in meteorological drought and temperature extremes over continental United States (CONUS). We used data from 4161 weather stations across the CONUS to downscale future precipitation projections from three Earth System Models (ESMs) participating in the Coupled Model Intercomparison Project Phase Six (CMIP6), specifically for the high emission scenario SSP5 8.5. Comparing historic observations with climate model projections revealed a significant bias in total annual precipitation days and total precipitation amounts. The average number of annual precipitation days across CONUS was projected to be 205 ± 26, 184 ± 33, and 181 ± 25 days in the BCC, CanESM, and UKESM models, respectively, compared to 91 ± 24 days in the observed data. Analyzing the duration of drought periods in different ecoregions of CONUS showed an increase in the number of drought months in the future (2023–2052) compared to the historical period (1989–2018). The analysis of precipitation and temperature changes in various ecoregions of CONUS revealed an increased frequency of droughts in the future, along with longer durations of warm spells. Eastern temperate forests and the Great Plains, which encompass the majority of CONUS agricultural lands, are projected to experience higher drought counts in the future. Drought projections show an increasing trend in future drought occurrences due to rising temperatures and changes in precipitation patterns. Our high-resolution climate projections can inform policy makers about the hotspots and their anticipated future trajectories.

To understand the role of the grain boundary (GB) in plasticity at small scale, a concurrently coupled mesoscale plasticity model was developed to simulate micro-bending of bicrystalline micron-sized beams. By coupling dislocation dynamics (DD) with a finite element model (FEM), a novel defect dynamics model provides the means to investigate intricate interactions between dislocations and GBs under various loading conditions. Our simulations of micro-bending agree well with corresponding micro-bending experiments, and they show that mechanical response of bicrystals could have not only hardening but also softening depending on the characters of the GB. In addition, changing the location of the GB in the microbeams results in different mechanical responses; GBs located at the neutral plane show softening compared to single crystals, while inclined GBs located halfway along the length of the beam show little effect. Simulation results could provide a clear picture on detailed dislocation-GB interactions, and quantitative resolved shear stress analysis supplemented by dislocation density distribution is used to analyze the mechanical response of bicrystalline samples.

Hierarchical optimization modeling in an algebraic modeling environment facilitates construction of large models with many interchangeable sub-models. However, for dynamic simulation and optimization applications, a flattened structure that preserves time indexing is preferred. To convert from a structure that facilitates model construction to a structure that facilitates dynamic optimization, the concept of reshaping an optimization model is introduced along with the recently developed utilities in the Pyomo algebraic modeling environment that make this possible. The application of these utilities to model predictive control simulations and partial differential equation (PDE) discretization stability analysis is discussed, and two challenging nonlinear model predictive control case studies are presented to demonstrate the advantages of this approach.

The fraction of tritium converted to the water form in a fire scenario is one of the metrics of greatest interest for radiological safety assessments. The conversion fraction is one of the prime variables contributing to the hazard assessment. This paper presents measurements of oxidation rates for the non-radioactive hydrogen isotopes (protium and deuterium) at sub-flammable concentrations that are typical of many of the most likely tritium release scenarios. These measurements are fit to a simplified 1-step kinetic rate expression, and the isotopic trends for protium and deuterium are extrapolated to produce a model appropriate for tritium. The effects of the new kinetic models are evaluated via CFD simulations of an ISO-9705 standard room fire that includes a trace release of hydrogen isotope (tritium), illustrating the high importance of the correct (measurement-based) kinetics to the outcome of the simulated conversion.

Abstract not provided.

The 2-year Puerto Rico Grid Resilience and Transition to 100% Renewable Energy Study analyzed stakeholder-driven pathways to Puerto Rico’s clean energy future. Outputs relating to electricity demand modeling were partially informed by estimates of electric vehicle adoption across all classes of medium- and heavy-duty vehicles (MHDVs), and the ensuing charging loads. To create these estimates, the team developed a transportation model for MHDVs in Puerto Rico to estimate the amount and geospatial distribution of energy used. Charging schedules for the different end uses of MHDVs were then used to construct electric load shapes assuming a portion of those vehicles would be replaced by battery electric counterparts. Study results showed that, by 2050, electric vehicles may constitute roughly 50% of the MHDV population in Puerto Rico. The resulting electrical demand curve attributable to MHDV charging showed that, for solar energy-based electrical systems with limited energy storage, this demand may create challenges unless appropriately managed either on the demand or supply side.

Abstract not provided.

A series of extensively instrumented tests was performed on the Structural Evaluation Test Unit in the early 1990s. The purpose of these tests was to determine the response of a minimally designed cask to impacts that were more severe than the design basis impact. This test series provides an excellent opportunity for benchmarking explicit dynamic finite element analysis programs for behaviors that may be experienced by casks during regulatory and extra-regulatory impact events. This report provides the parameters of the test unit, the locations of instrumentation, the locations of inspection points, and the parameters of the four tests that were conducted. A companion report provides the results of the tests.

The Strategic Petroleum Reserve (SPR) is the world’s largest supply of emergency crude oil. The reserve consists of four sites in Louisiana and Texas. Each site stores crude in deep, underground salt caverns. It is the mission of the SPR’s Enhanced Monitoring Program to examine available sensing data to inform our understanding of each site. This report discusses the monitoring data, processes, and results for each of the four sites for fiscal year 2023.

Abstract not provided.

Sandia National Laboratories is a premier United States national security laboratory which develops science-based technologies in areas such as nuclear deterrence, energy production, and climate change. Computing plays a key role in its diverse missions, and within that environment, Research Software Engineers (RSEs) and other scientific software developers utilize testing automation to ensure quality and maintainability of their work. We conducted a Participatory Action Research study to explore the challenges and strategies for testing automation through the lens of academic literature. Through the experiences collected and comparison with open literature, we identify these challenges in testing automation and then present strategies for mitigation grounded in evidence-based practice and experience reports that other, similar institutions can assess for their automation needs.

Achieving robust and efficient drilling is a critical part of reducing the cost of geothermal energy exploration and extraction. Drilling performance is often evaluated using one or more of three key metrics: depth of cut (DOC), rate of penetration (ROP), and mechanical specific energy (MSE). All three of these quantities are related to each other. DOC refers to the depth a bit penetrates into rock during drilling. This is an important quantity for estimating bit behavior. ROP is the simply the DOC multiplied by the rotational rate, and represents how quickly the drill bit is advancing through the ground. ROP is often the parameter used for drilling control and optimization. Finally, MSE provides insight into drilling efficiency and rock type. MSE calculations rely on ROP, drilling force, and drilling torque. Surface-based sensors at the top of the drill are often used to measure all these quantities. However, top-hole measurements can deviate substantially from the behavior at the bit due to lag, vibrations, and friction. Therefore, relying only on top-hole information can lead to suboptimal drilling control. In this work, we describe recent progress towards estimating ROP, DOC, and MSE using down-hole sensing. We assume down-hole measurements of torque, weight-on-bit (WOB). Our hypothesis is that these measurements can provide more rapid and accurate measures of drilling performance. We show how a multi-layer perceptron (MLP) machine learning algorithm can provide rapid and accurate performance when evaluated on experimental data taken from Sandia’s Hard Rock Drilling Facility. In addition, we implement our algorithms on an embedded system intended to emulate a bottom-hole-assembly for sensing and estimation. Our experimental results show that DOC can be estimated accurately and in real-time. These estimates when combined with measurements for rotary speed, torque, and force can provide improved estimates for ROP and MSE. These results have the potential to enable better drilling assessment, improved control, and extended component lifetimes.

This technical report serves to summarize a literature search conducted that covered confidence calibration. This report is meant to serve as a solid starting reference for individuals interested in learning more about the confidence calibration domain as well as for individuals more familiar with this work – as a summarizing document for calibration metrics is notably lacking in the literature. This report is not meant to serve as a comprehensive review of everything that has been done in this field – in fact, the reader is encouraged to look further into this domain. We describe confidence and calibration and discuss properties of good calibration metrics. We detail various calibration and calibration-tangential metrics, presenting equations, algorithms, parameters, and an analysis of strengths and weaknesses. We apply a subset of these metrics to eight proxy confidence assessment datasets. We examine the various metrics in the context of model confidence. Finally, we discuss promising future directions and outstanding questions.

Validation and verification of engineering models is important to understand potential weaknesses and issues in the model. This is accomplished through the application of constraint logic to the model. These models and the constraints put upon them can be represented through a graph structure. Here we give a visualization system to aid users understanding, locating, and fixing constraint violations in their systems. We give users several ways to narrow down on the specific errors and parts of the graph they’re interested in. Users have the opportunity to choose the types of errors that will be shown in the graph. Clustering is applied to the graph to help users narrow down their searches. Several other graph interactions are given to support discovery of constraint violations.

Abstract not provided.

A series of extensively instrumented tests was performed on the Structural Evaluation Test Unit in the early 1990s. The purpose of these tests was to determine the response of a minimally designed cask to impacts that were more severe than the design basis impact. This test series provides an excellent opportunity for benchmarking explicit dynamic finite element analysis programs for behaviors that may be experienced by casks during regulatory and extra-regulatory impact events. This report provides the results of the four tests that were conducted. It is meant to go along with a companion report that defines the benchmark problem and gives the locations for the instrumentation and inspection points.

The subject of Task F of DECOVALEX-2023 concerns performance assessment modelling of radioactive waste disposal in deep mined repositories. The primary objectives of Task F are to build confidence in the models, methods, and software used for performance assessment (PA) of deep geologic nuclear waste repositories, and/or to bring to the fore additional research and development needed to improve PA methodologies. In Task F2-(salt), these objectives have been accomplished through staged development and comparison of the models and methods used by participating teams in their PA frameworks. Coupled-process submodels and deterministic simulations of the entire PA model for a reference scenario for waste disposal in domal salt have been conducted. The task specification has been updated continuously since the initiation of the project to reflect the staged development of the conceptual repository model and performance metrics.

Abstract not provided.

Abstract not provided.

Abstract not provided.

This report is a comprehensive guide to the nonlinear viscoelastic Spectacular model, which is an isotropic, thermo-rheologically simple constitutive model for glass-forming materials, such as amorphous polymers. Spectacular is intermediate in complexity to the previous PEC and SPEC models (Potential Energy Clock and Simplified Potential Energy Clock models, respectively). The model form consists of two parts: a Helmholtz free energy functional and a nonlinear material clock that controls the rate of viscoelastic relaxation. The Helmholtz free energy is derived from a series expansion about a reference state. Expressions for the stress and entropy functionals are derived from the Helmholtz free energy following the Rational Mechanics approach. The material clock depends on a simplified expression for the potential energy, which itself is a functional of the temperature and strain histories. This report describes the thermo-mechanical theory of Spectacular, the numerical methods for time-integrating the model, model verification for its implementation in LAMÉ, a user guide for its implementation in LAMÉ, and ideas for future work. A number of appendices provide supplementary mathematical details and a description of the procedure used to derive the simplified potential energy from the full expression for the potential energy. The goal of this report is create a convenient point-of-entry for engineers who wish to learn more about Spectacular, but also to serve as a reference manual for advanced users of the model.

This report covers an inquiry into seismoacoustic array processing using infrasound arrivals combined with resulting Ground Coupled Airwaves (GCA) that are present on collocated seismic sensors. In preparation, data calibration and denoising is completed for a seismoacoustic sensor array that was deployed at the Facility for Acceptance, Calibration, and Testing on Kirtland Airforce Base from August through September of 2021. The events of interest for this study are small, local explosive sources that lead to short duration, impulsive signals on the instruments. The goal is to determine if combining infrasound signals with the corresponding GCAs on collocated seismic sensors can be used to improve the results returned by automated signal detection and characterization (e.g., back azimuth estimates). Preparation for seismic and infrasound data involves removing the instrument response so that sensors have flat power spectra over the frequency range 0.1-10 Hz, where signal from events of interest may be detected. After instrument response removal, deployment conditions specific to this array require a retrospective noise analysis to determine station emplacement characteristics. Once all data is calibrated, a manual search is performed for possible GCA arrivals across the seismoacoustic network. These arrivals are then processed through beamforming and subsequent event identification, resulting in a catalogue of seismoacoustic GCA arrivals with corresponding back azimuth and trace velocity estimations.

Redox flow batteries (RFBs) that incorporate solid energy-storing materials are attractive for high-capacity grid-scale energy storage due to their markedly higher theoretical energy densities compared to their fully liquid counterparts. However, this promise of higher energy density comes at the expense of rate capability. In this work we exploit a ZnO nanorod-decorated Ni foam scaffold to create a high surface area Li metal anode capable of rates up to 10 mA cm−2, a 10× improvement over traditional planar designs. The ZnO nanorods enhance Li metal wettability and promote uniform Li nucleation, allowing the RFB to be initially operated with a prelithiated (charged) anode, or with a safety-conscious, Li-less, fully discharged anode. 5 mgS cm−1 were cycled using a mediated S cathode, whereby redox mediators help oxidize and reduce solid S particles. At 2.4 mgS cm−2 and 10 mA cm−2, the RFB becomes limited by the mediation of solid S. Nevertheless, a respectable energy density of 20.3 Wh L−1 is demonstrated, allowing considerable increase if the S mediation rate can be further improved. Lessons learned here may be broadly applied to RFBs with alkali metal anodes, offering an avenue for safe, dense, grid-scale energy storage.

Abstract not provided.

There is currently very limited research into how experts analyze and assess potentially fraudulent content in their expertise areas, and most research within the disinformation space involves very limited text samples (e.g., news headlines). The overarching goal of the present study was to explore how an individual’s psychological profile and the linguistic features in text might influence an expert’s ability to discern disinformation/fraudulent content in academic journal articles. At a high level, the current design tasked experts with reading journal articles from their area of expertise and indicating if they thought an article was deceptive or not. Half the articles they read were journal papers that had been retracted due to academic fraud. Demographic and psychological inventory data collected on the participants was combined with performance data to generate insights about individual expert susceptibility to deception. Our data show that our population of experts were unable to reliably detect deception in formal technical writing. Several psychological dimensions such as comfort with uncertainty and intellectual humility may provide some protection against deception. This work informs our understanding of expert susceptibility to potentially fraudulent content within official, technical information and can be used to inform future mitigative efforts and provide a building block for future disinformation work.

Abstract not provided.

This report summarizes the water inputs associated with four technologies playing diverse roles in energy transitions: hydrogen, solar photovoltaics (PV), wind, and batteries. Information in this report is drawn from multiple sources, including peer-reviewed literature, industry and international agency reports, EcoInvent life cycle inventory database, and subject matter expert (SME) consultations. Where possible, insights that characterized water requirements for specific stages of the technology development (e.g., operations, manufacturing, and mining) were prioritized over broader cradle-to-gate assessment values. Furthermore, both direct and indirect water requirements (i.e., associated with associated energy inputs) were considered in this literature review.

Abstract not provided.

This report represents completion of milestone deliverable M2SF-24SN010309082 Annual Status Update for OWL due on November 30, 2023. It contains the status of fiscal year 2023 (FY2023) updates for the Online Waste Library (OWL).

Here, a review of current trends in scientific computing reveals a broad shift to open-source and higher-level programming languages such as Python and growing career opportunities over the next decade. Open-source modeling tools accelerate innovation in equation-based and data-driven applications. Significant resources have been deployed to develop data-driven tools (PyTorch, TensorFlow, Scikit-learn) from tech companies that rely on machine learning services to meet business needs while keeping the foundational tools open. Open-source equation-based tools such as Pyomo, CasADi, Gekko, and JuMP are also gaining momentum according to user community and development pace metrics. Integration of data-driven and principles-based tools is emerging. New compute hardware, productivity software, and training resources have the potential to radically accelerate progress. However, long-term support mechanisms are still necessary to sustain the momentum and maintenance of critical foundational packages.

During February 2023, a total of 32 individual distributed acoustic sensing (DAS) systems acted jointly as a global seismic monitoring network. The aim of this Global DAS Month campaign was to coordinate a diverse network of organizations, instruments, and file formats to gain knowledge and move toward the next generation of earthquake monitoring networks. During this campaign, 156 earthquakes of magnitude 5 or larger were reported by the U.S. Geological Survey and contributors shared data for 60 min after each event’s origin time. Participating systems represent a variety of manufacturers, a range of recording parameters, and varying cable emplacement settings (e.g., shallow burial, borehole, subaqueous, and dark fiber). Monitored cable lengths vary between 152 and 120,129 m, with channel spacing between 1 and 49 m. The data has a total size of 6.8 TB, and are available for free download. Finally, organizing and executing the Global DAS Month has produced a unique dataset for further exploration and highlighted areas of further development for the seismological community to address.

Computational simulation is increasingly relied upon for high/consequence engineering decisions, which necessitates a high confidence in the calibration of and predictions from complex material models. However, the calibration and validation of material models is often a discrete, multi-stage process that is decoupled from material characterization activities, which means the data collected does not always align with the data that is needed. To address this issue, an integrated workflow for delivering an enhanced characterization and calibration procedure—Interlaced Characterization and Calibration (ICC)—is introduced and demonstrated. Further, this framework leverages Bayesian optimal experimental design (BOED), which creates a line of communication between model calibration needs and data collection capabilities in order to optimize the information content gathered from the experiments for model calibration. Eventually, the ICC framework will be used in quasi real-time to actively control experiments of complex specimens for the calibration of a high-fidelity material model. This work presents the critical first piece of algorithm development and a demonstration in determining the optimal load path of a cruciform specimen with simulated data. Calibration results, obtained via Bayesian inference, from the integrated ICC approach are compared to calibrations performed by choosing the load path a priori based on human intuition, as is traditionally done. The calibration results are communicated through parameter uncertainties which are propagated to the model output space (i.e. stress–strain). In these exemplar problems, data generated within the ICC framework resulted in calibrated model parameters with reduced measures of uncertainty compared to the traditional approaches.

Bimetallic, reactive multilayers are uniformly structured materials composed of alternating sputter-deposited layers that may be ignited to produce self-propagating mixing and formation reactions. These nanolaminates are most commonly used as rapid-release heat sources. The specific chemical composition at each metal/metal interface determines the rate of mass transport in a mixing and formation reaction. The inclusion of engineered diffusion barriers at each interface will not only inhibit solid-state mixing but also may impede the self-propagating reactions by introducing instabilities to wavefront morphology. This work examines the effect of adding diffusion barriers on the propagation of reaction waves in Co/Al multilayers. The Co/Al system has been shown to exhibit a reaction propagation instability that is dependent on the bilayer thickness, which allows for the occurrence of unstable modes in otherwise stable designs from the inclusion of diffusion barriers. Based on the known stability criteria in the Co/Al multilayer system, the way in which the inclusion of diffusion barriers changes a multilayer's heat of reaction, thermal conductivity, and material mixing mechanisms can be determined. These factors, in aggregate, lead to changes in the wavefront velocity and stability.

Photonic topological insulators exhibit bulk-boundary correspondence, which requires that boundary-localized states appear at the interface formed between topologically distinct insulating materials. However, many topological photonic devices share a boundary with free space, which raises a subtle but critical problem as free space is gapless for photons above the light line. Here, we use a local theory of topological materials to resolve bulk-boundary correspondence in heterostructures containing gapless materials and in radiative environments. In particular, we construct the heterostructure’s spectral localizer, a composite operator based on the system’s real-space description that provides a local marker for the system’s topology and a corresponding local measure of its topological protection; both quantities are independent of the material’s bulk band gap (or lack thereof). Moreover, we show that approximating radiative outcoupling as material absorption overestimates a heterostructure’s topological protection. Importantly, as the spectral localizer is applicable to systems in any physical dimension and in any discrete symmetry class (i.e., any Altland-Zirnbauer class), our results show how to calculate topological invariants, quantify topological protection, and locate topological boundary-localized resonances in topological materials that interface with gapless media in general.

Here, using atomistic molecular dynamics simulations, we investigate the morphology and transport properties of a new family of fluorine-free terpolymers designed as proton-exchange membranes. Simulated random terpolymers consist of three monomers with a 5-carbon backbone with a phenylsulfonate, phenyl, or no pendant group and have ion exchange capacities (IECs) ranging from 1.06–4.14 mmol/g. At a hydration level of 9, cluster analysis reveals macrophase separation between water and terpolymers with IEC < 2.1 mmol/g and continuous, percolated hydrophilic and hydrophobic nanoscale domains at higher IECs. Channel width distribution analysis of the percolated morphologies revealed that more hydrophobic units produce less uniform channels. Decreasing the surface area per sulfonate group and increasing the fractal dimension of the hydrophilic domains correlate with increased water diffusivity, due to a more acidic interface and more isotropic water channels. Relative to the previously studied phenylsulfonate homopolymer, these terpolymers with lower IECs have only modestly lower water diffusion, and we anticipate other advantages related to processability.

Nonlinear topological insulators have garnered substantial recent attention as they have both enabled the discovery of new physics due to interparticle interactions, and may have applications in photonic devices such as topological lasers and frequency combs. However, due to the local nature of nonlinearities, previous attempts to classify the topology of nonlinear systems have required significant approximations that must be tailored to individual systems. Here, we develop a general framework for classifying the topology of nonlinear materials in any discrete symmetry class and any physical dimension. Our approach is rooted in a numerical $K$ -theoretic method called the spectral localizer, which leverages a real-space perspective of a system to define local topological markers and a local measure of topological protection. Here, our nonlinear spectral localizer framework yields a quantitative definition of topologically nontrivial nonlinear modes that are distinguished by the appearance of a topological interface surrounding the mode. Moreover, we show how the nonlinear spectral localizer can be used to understand a system's topological dynamics, i.e., the time evolution of nonlinearly induced topological domains within a system. We anticipate that this framework will enable the discovery and development of novel topological systems across a broad range of nonlinear materials.

Derived from renewable feedstocks, such as biomass, polylactic acid (PLA) is considered a more environmentally friendly plastic than conventional petroleum-based polyethylene terephthalate (PET). However, PLA must still be recycled, and its growing popularity and mixture with PET plastics at the disposal stage poses a cross-contamination threat in existing recycling facilities and results in low-value and low-quality recycled products. Hybrid upcycling has been proposed as a promising sustainable solution for mixed plastic waste, but its techno-economic and life cycle environmental performance remain understudied. Here we propose a hybrid upcycling approach using a biocompatible ionic liquid (IL) to first chemically depolymerize plastics and then convert the depolymerized stream via biological upgrading with no extra separation. We show that over 95% of mixed PET/PLA was depolymerized into the respective monomers, which then served as the sole carbon source for the growth of Pseudomonas putida, enabling the conversion of the depolymerized plastics into biodegradable polyhydroxyalkanoates (PHAs). In comparison to conventional commercial PHAs, the estimated optimal production cost and carbon footprint are reduced by 62% and 29%, respectively.

The United States Department of Energy’s (DOE) Office of Nuclear Energy’s Spent Fuel and Waste Science and Technology Campaign seeks to better understand the technical basis, risks, and uncertainty associated with the safe and secure disposition of spent nuclear fuel (SNF) and high-level radioactive waste. Commercial nuclear power generation in the United States has resulted in thousands of metric tons of SNF, the disposal of which is the responsibility of DOE (Nuclear Waste Policy Act of 1982, as amended). Any repository licensed to dispose of SNF must meet requirements regarding the long-term performance of that repository. The evaluation of long-term performance of the repository may need to consider the SNF achieving a critical configuration during the postclosure period. Of particular interest is the potential for this situation to occur in dual-purpose canisters (DPCs), which are currently licensed and being used to store and transport SNF but were not designed for permanent geologic disposal. DOE has been considering disposing of SNF in DPCs to avoid the costs and worker dose associated with repackaging the SNF currently stored in DPCs into repository-specific canisters. This report examines the consequences of postclosure criticality to provide technical support to DOE in developing a disposal plan.

An optically recording velocity interferometer system has been used to measure acceleration histories and maximum velocities for laser-driven aluminum foil targets launched from the output face of optical fibers. Peak flyer velocities have been determined as a function of various parameters, including driving laser fluence, laser pulse duration and target thickness. The results at high fluences are consistent with a nearly constant efficiency of coupling optical energy into flyer kinetic energy and a small ablated mass fraction; however, the coupling efficiency falls off rapidly at fluences < 15 J-cm^-2. Measurement of the time delay between laser pulse arrival at the target and the onset of flyer motion have also been performed. Significant delays are observed at low fluences, arising from the increased time required for plasma formation at the fiber/foil interface under these conditions.

Aerosol particles are known to exist in highly viscous amorphous states at a low relative humidity and temperature. The slow diffusion of molecules in viscous particles impacts the uptake and loss of volatile and semivolatile species and the rate of heterogeneous chemistry. Recent work has demonstrated that in particles containing organic molecules and salts, the formation of two-phase gel states is possible, leading to observations of rigid particles that resist coalescence. The way that molecules diffuse and transport in gel systems is not well-characterized. In this work, we use an electrodynamic balance to levitate sample particles containing a range of organic compounds in mixtures with calcium chloride and measure the rate of water diffusion. Particles of the pure organics have been shown to form viscous amorphous states, while in mixtures with divalent salts, coalescence measurements have revealed the apparent solidification of particles, consistent with the formation of a gel state facilitated by ion-molecule interactions. We report in several cases that water transport can actually be increased in the rigid gel state relative to the pure compound that forms a viscous state under similar conditions. These measurements reveal the limitations of using viscosity as a metric for predicting molecular diffusion and that the gel structure that forms is a much stronger controlling factor in the rate of diffusion. This underscores the need for diffusion measurements as well as a deeper understanding of noncovalent molecular assembly that leads to supramolecular structures in aerosol particles.

The purpose of pvOps is to support empirical evaluations of data collected in the field related to the operations and maintenance (O&M) of photovoltaic (PV) power plants. pvOps presently contains modules that address the diversity of field data, including text-based maintenance logs, current-voltage (IV) curves, and timeseries of production information. The package functions leverage machine learning, visualization, and other techniques to enable cleaning, processing, and fusion of these datasets. These capabilities are intended to facilitate easier evaluation of field patterns and extraction of relevant insights to support reliability-related decision-making for PV sites. The open-source code, examples, and instructions for installing the package through PyPI can be accessed through the GitHub repository.

Operator learning for complex nonlinear systems is increasingly common in modeling multi-physics and multi-scale systems. However, training such high-dimensional operators requires a large amount of expensive, high-fidelity data, either from experiments or simulations. In this work, we present a composite Deep Operator Network (DeepONet) for learning using two datasets with different levels of fidelity to accurately learn complex operators when sufficient high-fidelity data is not available. Additionally, we demonstrate that the presence of low-fidelity data can improve the predictions of physics-informed learning with DeepONets. We demonstrate the new multi-fidelity training in diverse examples, including modeling of the ice-sheet dynamics of the Humboldt glacier, Greenland, using two different fidelity models and also using the same physical model at two different resolutions.

Electric fields are commonplace in plasmas and affect transport by driving currents and, in some cases, instabilities. The necessary condition for instability in collisionless plasmas is commonly understood to be described by the Penrose criterion, which quantifies a sufficient relative drift between different populations of particles that must be present for wave amplification via inverse Landau damping. For example, electric fields generate drifts between electrons and ions that can excite the ion-acoustic instability. Here, we use particle-in-cell simulations and linear stability analysis to show that the electric field can drive a fundamentally different type of kinetic instability, named the electron-field instability. This instability excites electron plasma waves with wavelengths ≳30λ_De, has a growth rate that is proportional to the electric field strength, and does not require a relative drift between electrons and ions. The Penrose criterion does not apply when accounting for the electric field. Furthermore, the large value of the observed frequency, near the electron plasma frequency, further distinguishes it from the standard ion-acoustic instability, which oscillates near the ion plasma frequency. The ubiquity of macroscopic electric fields in quasineutral plasmas suggests that this instability is possible in a host of systems, including low-temperature and space plasmas. In fact, damping from neutral collisions in such systems is often not enough to completely damp the instability, adding to the robustness of the instability across plasma conditions.

Sandia National Laboratories (SNL) and the Institut de Radioprotection et de Sûreté Nucléaire (IRSN) have collaborated on the design and execution of a set of critical experiments that explore the effects of molybdenum in water moderated fuel-rod arrays. The molybdenum is included as sleeves (tubes) on some of the fuel rods in the arrays. The fuel used in the experiments is known at Sandia as the Seven Percent Critical Experiment (7uPCX) fuel. This fuel has been used is several published benchmark evaluations in including LEU-COMP-THERM-78 and LEU-COMP THERM-080.

Characterizing and quantifying microstructure evolution is critical to forming quantitative relationships between material processing conditions, resulting microstructure, and observed properties. Machine-learning methods are increasingly accelerating the development of these relationships by treating microstructure evolution as a pattern recognition problem, discovering relationships explicitly or implicitly. These methods often rely on identifying low-dimensional microstructural fingerprints as latent variables. However, using inappropriate latent variables can lead to challenges in learning meaningful relationships. In this work, we survey and discuss the ability of various linear and nonlinear dimensionality reduction methods including principal component analysis, autoencoders, and diffusion maps to quantify and characterize the learned latent space microstructural representations and their time evolution. We characterize latent spaces by their ability to represent high-dimensional microstructural data in terms of compression achieved as a function of the number of latent dimensions required to represent the data accurately, their accuracy based on their reconstruction performance, and the smoothness of the microstructural trajectories in latent dimension. We quantify these metrics for common microstructure evolution problems in material science including spinodal decomposition of a binary metallic alloy, thin film deposition of a binary metallic alloy, dendritic growth, and grain growth in a polycrystal. This study provides considerations and guidelines for choosing dimensionality reduction methods when considering materials problems that involve high dimensional data and a variety of features over a range of lengths and time scales.

This project has progressed several physics models in the EMPIRE plasma simulation code to achieve higher fidelity simulations of high-current diode operation in pulsed-power accelerators. In this report, we present details for the following major work products: (1) a set of verification problems covering all key processes involved in gap closure physics has been designed; this suite has facilitated feature vetting and overall model maturation, (2) a new EMPIRE exemplar has been developed: the Radiographic Integrated Test Stand 6 (RITS-6) diode, and (3) An exemplar for the Saturn accelerator exemplar was enabled by the models matured under this work to self-consistently simulate further into the diode pulse than previously possible (bipolar flow regime). These developments have lead to the highest confidence EMPIRE power flow predictions of the Saturn accelerator to date. Additionally, three modeling approaches for simulating electrode plasmas have been investigated. We report on these results and provide recommendations.

Methyl formate (MF; CH₃OCHO) is the smallest representative of esters, which are common components of biodiesel. The present study characterizes the thermal dissociation kinetics of the radicals formed by H atom abstraction from MF—CH₃OCO and CH₂OCHO—through a combination of modeling, experiment, and theory. For the experimental effort, excimer laser photolysis of Cl₂ was used as a source of Cl atoms to initiate reactions with MF in the gas phase. Time-resolved species profiles of MF, Cl₂, HCl, CO₂, CH₃, CH₃Cl, CH₂O, and CH₂ClOCHO were measured and quantified using photoionization mass spectrometry at temperatures of 400–750 K and 10 Torr. The experimental data were simulated using a kinetic model, which was informed by ab initio-based theoretical kinetics calculations and included chlorine chemistry and secondary reactions of radical decomposition products. Here, we calculated the rate coefficients for the H-abstraction reactions Cl + MF → HCl + CH₃OCO (R1a) and Cl + MF → HCl + CH₂OCHO (R1b): k_1a,theory = 6.71 × 10^–15·T^1.14·exp(—606/T) cm³/molecule·s; k_1b,theory = 4.67 × 10^–18·T^2.21·exp(—245/T) cm³/molecule·s over T = 200–2000 K. Electronic structure calculations indicate that the barriers to CH₃OCO and CH₂OCHO dissociation are 13.7 and 31.6 kcal/mol and lead to CH₃ + CO₂ (R3) and CH₂O + HCO (R5), respectively. The master equation-based theoretical rate coefficients are k_3,theory (P = ∞) = 2.94 × 10⁹·T^1.21·exp(—6209/T) s^–1 and k_5,theory (P = ∞) = 8.45 × 10⁸·T^1.39·exp(—15132/T) s^–1 over T = 300–1500 K. The calculated branching fractions into R1a and R1b and the rate coefficient for R5 were validated by modeling of the experimental species time profiles and found to be in excellent agreement with theory. Additionally, we found that the bimolecular reactions CH₂OCHO + Cl, CH₂OCHO + Cl₂, and CH₃ + Cl₂ were critical to accurately model the experimental data and constrain the kinetics of MF-radicals. Inclusion of the kinetic parameters determined in this study showed a significant impact on combustion simulations of larger methyl esters, which are considered as biodiesel surrogates.

An interferometric radioimager provides real-time, high-fidelity radioimaging of high voltage breakdown (HVB) both internal and external to electrical components at sub-nanosecond and sub-millimeter resolution and has an ability to resolve multiple/spatially-extensive HVB simultaneously. Therefore, radioimaging can be used to screen for early life weakness/failure and enable non-destructive screening of defective electrical components. In particular, radioimaging can detect precursors to catastrophic HVB, allowing for early detection of weakness in critical electrical components. Radioimaging can also be used to track HVB and pinpoint defects in electrical components real time, including transformers, capacitors, cables, switches, and microelectronics.

Conceptually, single-ion polymer electrolytes (SIPE) with the anion bound to the polymer could solve major issues in Li-ion batteries, but their conductivity is too low. Experimentally, weakly interacting anionic groups have the best conductivity. To provide a theoretical basis for this result, density functional theory calculations of the optimized geometries and energies are performed for charged ligands used in SIPE. Comparison is made to neutral ligands found in dual-ion conductors, which demonstrate higher conductivity. Further, the free energy differences between adding and subtracting a ligand are small enough for the neutral ligands to have the conductivity seen experimentally. However, charged ligands have large barriers, implying that lithium transport will coincide with the slow polymer diffusion, as observed in experiments. Overall, SIPE will require additional solvent to achieve a sufficiently high conductivity. Additionally, the binding of mono- and bidentate geometries varies, providing a simple and clear reason that polarizable force fields are required for detailed interactions.

Nearly all metals, alloys, ceramics, and their associated composites are polycrystalline in nature, with grain boundaries that separate well-defined crystalline regions that influence materials properties. In all but the most pure elemental systems, intentional solutes or impurities are present and can segregate to, or less commonly away from, the grain boundaries, in turn influencing boundary behavior, their stability, and associated materials properties. In some cases, grain-boundary segregation can also trigger “phase-like” structural transitions that dramatically alter the essential nature of the boundary. With the development of advanced electron microscopy techniques, researchers can directly observe grain-boundary structures and segregation with atomic precision. Despite such spatial resolution, the underlying mechanisms governing grain-boundary segregation remain difficult to characterize. As a result, computational modeling techniques such as density functional theory, molecular dynamics, mesoscale phase-field, continuum defect theory, and others are important complementary tools to experimental observations for studying grain-boundary segregation behavior. In conclusion, these computational methods offer the ability to explore the underlying formation mechanisms of grain-boundary segregation, elucidate complex segregation behavior, and provide insights into solutions to effectively controlling microstructure.

Despite their many advantages, covalent organic frameworks (COFs) built from three-dimensional monomers are synthetically difficult to functionalize. Herein, we provide a new synthetic approach to the functionalization of a three-dimensional covalent organic framework (COF-300) by using a series of solid-state linkage transformations. By reducing the imine linkages of the framework to amine linkages, we produced a more hydrolytically stable material and liberated a nucleophilic amino group, poised for further functionalization. We then treated the amine-linked COF with diverse electrophiles to generate a library of functionalized materials, which we tested for their ability to adsorb perfluoroalkyl substances (PFAS) from water. The framework functionalized with dimethylammonium groups, COF-300-dimethyl, adsorbed more than 250 mg of perfluorooctanoic acid (PFOA) per 1 g of COF, which represents an approximately 14,500-fold improvement over that of COF-300 and underscores the importance of electrostatic interactions to PFAS adsorption performance. In conclusion, this work provides a conceptually new approach to the design and synthesis of functional three-dimensional COFs.

Sandia’s contribution to the FY23/Q4 GNG milestone was met: Demonstrate one composite hydride material with experimental data showing reversibility (50 cycles) and volumetric capacity ≥ 350 bar compressed gas.

Methyl-ethyl-substituted Criegee intermediate (MECI) is a four-carbon carbonyl oxide that is formed in the ozonolysis of some asymmetric alkenes. MECI is structurally similar to the isoprene-derived methyl vinyl ketone oxide (MVK-oxide) but lacks resonance stabilization, making it a promising candidate to help us unravel the effects of size, structure, and resonance stabilization that influence the reactivity of atmospherically important, highly functionalized Criegee intermediates. We present experimental and theoretical results from the first bimolecular study of MECI in its reaction with SO2, a reaction that shows significant sensitivity to the Criegee intermediate structure. Using multiplexed photoionization mass spectrometry, we obtain a rate coefficient of (1.3 ± 0.3) × 10-10 cm3 s-1 (95% confidence limits, 298 K, 10 Torr) and demonstrate the formation of SO3 under our experimental conditions. Through high-level theory, we explore the effect of Criegee intermediate structure on the minimum energy pathways for their reactions with SO2 and obtain modified Arrhenius fits to our predictions for the reaction of both syn and anti conformers of MECI with SO2 (ksyn = 4.42 × 1011 T-7.80exp(−1401/T) cm3 s-1 and kanti = 1.26 × 1011 T-7.55exp(−1397/T) cm3 s-1). Our experimental and theoretical rate coefficients (which are in reasonable agreement at 298 K) show that the reaction of MECI with SO2 is significantly faster than MVK-oxide + SO2, demonstrating the substantial effect of resonance stabilization on Criegee intermediate reactivity.

Abstract not provided.

Abstract not provided.

Ion diffusion in atmospheric pressure plasmas is examined and particular attention is paid to the fact that ion-ion interactions can be influenced by strong Coulomb coupling. Three regimes are identified. At low ionization fractions ( x i ≲ 10 − 6 ), standard weakly correlated ion-neutral interactions set the diffusion rate. At moderate ionization fractions ( 10 − 6 ≲ x i ≲ 10 − 2 ) there is a transition from ion-neutral to ion-ion collisions setting the diffusion rate. In this regime, the effect of strong Coulomb coupling in ion-ion collisions is accounted for by applying the mean force kinetic theory. Since both ion-neutral and ion-ion interactions contribute a comparable amount to the total diffusion rate, models (such as particle-in-cell or fluid) must account for both contributions. At high ionization fractions ( x i ≳ 10 − 2 ), strongly correlated ion-ion collisions dominate and the plasma is heated substantially by a disorder-induced heating (DIH) process associated with strong correlations. The temperature increase due to DIH strongly influences the ion diffusion rate. This effect becomes even more important, and occurs at lower ionization fractions, as the pressure increases above atmospheric pressure. In addition to ion diffusion, DIH affects the neutral gas temperature, therefore influencing the neutral diffusion rate. Model predictions are tested using molecular dynamics simulations, which included a Monte Carlo collision routine to simulate the effect of ion-neutral collisions at the lowest ionization fractions. The model and simulations show good agreement over a broad range of ionization fractions. The results provide a model for ion diffusion, on a wide range of ionization fractions and pressures, solely considering the elastic contribution to the diffusion coefficient—as an illustration of how strong Coulomb coupling influences diffusion processes in general.

Recently developed techniques for assembling moiré heterostructures have enabled unprecedented control of the interlayer interactions within the heterostructure, allowing for the realization of strongly correlated electronic states. The experimentally realized states, to date, include a variety of topological states, superconductivity, and strange metal behavior. However, the phenomenology varies significantly from sample to sample. Here we will present microscopy techniques for further characterization of both the progenitor exfoliated materials and the composite moiré heterostructure. The application of these techniques has allowed for an unprecidented degree of characterization in ambient conditions. Additionally we will explore the use of resistively detected electron spin resonance (ESR) as a novel probe for spin order in moirés. While our results do not show clearly dispersive features that would unambiguously evince spin excitations within the sample, the technique is still promising.

Abstract not provided.

Mechanical metamaterials are regularly implemented in engineering applications due to their unique properties derived from their structural geometry and material composition. This study incorporates deep reinforcement learning, a subset of machine learning that teaches an agent to complete a task through interactive experiences, into mechanical metamaterial design. The approach creates a design environment for the reinforcement learning agent to iteratively construct metamaterials with tailorable deformation and hysteretic characteristics. Validation involved producing metamaterials with a thermoplastic polyurethane (TPU) base material that exhibited the deformation response of expanded thermoplastic polyurethane (E-TPU) while maximizing or minimizing hysteresis in cyclic compression. This alignment confirmed the feasibility of tailoring deformation and energy manipulation using mechanical metamaterials. The agent's generalizability was tested by tasking it to create various metamaterials with distinct loading deformation responses and specific hysteresis goals in a simulated setting. The agent consistently delivered metamaterials that met loading curve criteria and demonstrated favorable energy return. This work demonstrates the potential of deep reinforcement learning as a rapid and effective tool for designing mechanical metamaterials with customizable traits. It ushers in the possibility of on-demand metamaterial design solutions, opening avenues across industries like footwear, wearables, and medical equipment.

American Vacuum Society National Conference presentation on ferroelectric hafnium oxide

Abstract not provided.

Full-field, multi-measurand diagnostics provide rich validation data necessary to improve the product life cycle time of nuclear safety components. Thermophosphor digital image correlation (TP+DIC) is a method of simultaneously measuring strain and temperature fields using patterned phosphor coatings deposited with aerosol deposition (AD). While TP+DIC produces a functional diagnostic, the coating’s reproducibility and the effect of the patterned features on the inferred temperature remains uncharacterized. This NSR&D project provided the opportunity to study two areas: 1) the tunability and repeatability of aerosol deposition and 2) the robustness of aerosol deposition phosphor on deforming substrates. The first area explores the process-property relationship of parameters elucidating the significance of each on the coating. The second area explores the relationship between the features’ characteristics (namely thickness) and the phosphor emission and inferred temperature. Together, the results will lead to the improved accuracy and functionality of TP+DIC for qualification testing of nuclear safety components.

Five alternative design configurations for a heavy-duty hydrogen refueling truck stop are detailed in this work. Each of the station concepts provides fast, 5-minute, 50 kg fills of up to 4 vehicles simultaneously, with a station capacity of 4200 kg/day. Two on-site production stations using PEM electrolysis are considered: one with off-peak production of the daily capacity; and one with on-demand production of hydrogen during vehicle refueling. Three delivered liquid hydrogen station concepts are considered: one with the same, high-pressure cascade storage system for dispensing as the electrolysis supplied stations, with low-pressure vaporization of the liquid hydrogen and pressurization via a compressor; and two with on-demand pressurization: one by low-pressure vaporization and compressors; and one with a cryogenic pump and high-pressure vaporization. Design, economic, and operational considerations for each of the components needed in these station concepts is provided. Of all the station concepts, the delivered liquid station with a low-pressure vaporizer and a cascade dispensing system has the lowest capital costs and equipment footprint, but the second highest operating costs primarily due to high costs for liquid hydrogen delivery. The lowest operating cost station is that with on-site production via PEM electrolysis at off-peak hours with a cascade delivery system. The low-pressure buffer storage system and electrolyzers have a large footprint and considerable capital costs, but could result in a low total cost of ownership, depending on the design timeline. The liquid hydrogen station with a cryo-pump has moderate capital costs, the lowest operating costs of the three delivered hydrogen stations, and the same small equipment footprint as the delivered liquid, cascade dispensing system. As cryogenic pumping technology improves and the capital costs for these pumps decreases, this station concept will become even more favorable. Three-dimensional renderings of the five station concepts provide station designers with a starting point for the development of heavy-duty refueling stations.

Abstract not provided.

The research described here was performed as part of the DOE SciDAC project Coupling Approaches for Next Generation Architectures (CANGA). A framework was developed for the derivation of novel algorithms for the multirate time integration of two-component systems coupled across an interface between spatial domains. The multirate aspect means that different time steps are allowed by each component integrator. The framework provides a way to construct multirate integrators with desirable properties related to stability, accuracy and preservation of system invariants. This report describes the framework and summarizes the major results, examples and research products.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Neural networks are becoming the cornerstone for national security prediction tasks. However, designing them requires significant research and trial/error, as they have many hyperparameters, including their computation graph (“architecture”). Neural architecture search (NAS) employs secondary optimizers to search for architectures maximizing objectives like accuracy. Evolutionary algorithms (EAs) are the most used class of optimizer for NAS. However, existing Python libraries for writing EAs limit the complexity of experiments a user can design. In this project, we built ARENA, a Python framework that encodes complex, hyper-realistic EAs. ARENA collects detailed information as it runs and is flexible enough to encode non-EA search algorithms. We tested ARENA on 4 toy optimization problems by encoding 3 search algorithms for each—random search, an EA, and simulated annealing. We also designed an EA that performs NAS on the MNIST dataset. Our experiments suggest the potential for immediate mission impact through solving lab-wide optimization problems.

This manual describes the use of the Xyce Parallel Electronic Simulator. Xyce has been designed as a SPICE-compatible, high-performance analog circuit simulator, and has been written to support the simulation needs of the Sandia National Laboratories electrical designers. This development has focused on improving capability over the current state-of-the-art in the following areas: (1) Capability to solve extremely large circuit problems by supporting large-scale parallel computing platforms (up to thousands of processors). This includes support for most popular parallel and serial computers. (2) A differential-algebraic-equation (DAE) formulation, which better isolates the device model package from solver algorithms. This allows one to develop new types of analysis without requiring the implementation of analysis-specific device models. (3) Device models that are specifically tailored to meet Sandia’s needs, including some radiation-aware devices (for Sandia users only). (4) Object-oriented code design and implementation using modern coding practices. Xyce is a parallel code in the most general sense of the phrase — a message passing parallel implementation — which allows it to run efficiently a wide range of computing platforms. These include serial, shared-memory and distributed-memory parallel platforms. Attention has been paid to the specific nature of circuit-simulation problems to ensure that optimal parallel efficiency is achieved as the number of processors grows.

The objective of this project was to evaluate material- and chemical-based solutions for hydrogen storage in rail applications as an alternative to high-pressure hydrogen gas and liquid hydrogen. Three use cases were assessed: yard switchers, long-haul locomotives, and tenders. Four storage options were considered: metal hydrides, nanoporous sorbents, liquid organic hydrogen carriers, and ammonia, using 700 bar compressed hydrogen as a benchmark. The results suggest that metal hydrides, currently the most mature of these options, have the highest potential. Storage in tenders is the most likely use case to be successful, with long-haul locomotives the least likely due to the required storage capacities and weight and volume constraints. Overall, the results are relevant for high-impact regions, such as the South Coast Air Quality Management District, for which an economical vehicular hydrogen storage system with minimal impact on cargo capacity could accelerate adoption of fuel cell electric locomotives. The results obtained here will contribute to the development of technical storage targets for rail applications that can guide future research. Moreover, the knowledge generated by this project will assist in development of material-based storage for stationary applications such as microgrids and backup power for data centers.

Abstract not provided.

Large, pulsed power and high voltage systems often employ a stack of insulators to separate a vacuum section away from water or oil sections. The size of this insulator stack often drives overall costs and feasibility of these systems. An electric breakdown along the insulator surface is a primary failure mechanism and is especially impactful if it occurs while power is still being delivered downstream. This report describes a set of experimental and modeling investigations into the cause of these breakdowns, especially focusing on the much less well-understood anode-initiated breakdowns that occur during early parts of power delivery. Additionally, new diagnostics for assessing relevant material properties and behavior of insulators are described. These results describe breakdown behavior and evolution at new temporal and spatial fidelities and provide hypotheses and some answers as to how these breakdowns can occur. This new understanding of the roles of different physics phenomena guide modifications and trade-offs in generating newer insulator stack designs that are smaller and/or have higher electrical stress thresholds.

A near net shape coating is desired to be applied to the outer surface of a capped cylinder (“cake pan”) type substrate using thermal spray technology. A capped cylinder geometry is more complex than simple coupon-level substrate substrates (e.g., flat panels, cylinders) and thus requires a more complex toolpath to deposit a uniform coating. This report documents a practical theoretical approach to calculating relative torch-to-substrate speeds for coating the cylindrical, corner, and cap region of a rotating capped cylinder based on fundamental thermal spray toolpath principles. A preliminary experimental test deposited a thermal spray coating onto a mock substrate using toolpath speeds calculated by the theoretical approach proposed. The mock substrate was metallographically inspected to assess coating uniformity across the cylindrical, corner, and cap region. Inspection of the mock substrate revealed qualitatively uniform coating microstructure and thickness where theoretically predicted, demonstrating the viability of the proposed toolpath method and associated calculations. Pathways forward to optimizing coating uniformity at the cap center are proposed as near term suggested future work.

As a follow-up to our previous report on quantum sensing for safeguards, here we delve deeper into quantum-enhanced imaging & spectroscopy and address their relevance to international safeguards. Much of the approaches rely on entangled photons, a quantum phenomenon not possible with classical physics, although just correlated photons will work for some applications, such as ghost imaging. We provide a comprehensive survey of quantum approaches, including multiple entangled photon ghost imaging and spectroscopy techniques. Entangled photons for noise reduction are also described, as well as Non-Line-Of-Sight imaging, compressive techniques, and squeezed light. Of particular interest is the generation of entangled photons with large wavelength separation, such as infrared/visible entangled photon pairs. Such entangled pairs would allow interaction with objects in the IR, such as in the molecular “fingerprint” wavelength region, while the recording device captures the visible photons, thus leveraging the high efficiency and lower cost of visible detectors. Unfortunately, entangled x-ray photons are not practical, which would have been useful for safeguards to interrogate shielded materials. Entangled gamma rays are even further beyond reason. We provide our assessment for application of quantum-enhanced imaging & spectroscopy for international safeguards, including suggested improvements to existing IAEA instruments and destructive assay measurements that are done at IAEA lab facilities.

Directed energy deposition (DED) is an attractive additive manufacturing (AM) process for large structural components. The rapid solidification and layer-by-layer process associated with DED results in non-ideal microstructures, such as large grains with strong crystallographic textures. These non-ideal microstructures can lead to severe anisotropy in the mechanical properties. Despite these challenges, DED has been identified as a potential solution for the manufacturing of near net shape Ti-6Al-4V preforms, replacing lost casting and forging capabilities. Two popular wire-based directed energy deposition (W-DED) processes were considered for the manufacturing of Ti-6Al-4V with assessments on their respective metallurgical and mechanical properties, as compared to a conventionally processed material. The two W-DED processes explored were wire arc additive manufacturing (WAAM) and electron beam additive manufacturing (EBAM). High throughput inspection and tensile testing procedures were utilized to generate statistically relevant data sets related to each process and sample orientation. The 2 AM technologies produced material with remarkably different microstructures and mechanical properties. Results revealed key differences in strength and ductility for the two disparate processes which were found to be related to differences in the metallurgical properties.

Abstract not provided.

Sandia National Laboratories (SNL) was contracted by the United States Department of Energy Environmental Management (DOE-EM), Los Alamos Field Office to perform mechanical and thermal scoping calculations as part of a study seeking to understand the ignitability risk of the Remediated Nitrate Salts (RNS) waste drums during transportation from the Waste Control Specialists (WCS) facility to Waste Isolation Pilot Plant (WIPP) and permanent disposal of the waste at WIPP. The scoping thermal simulations described in this report pertain to thermal calculations performed with a packaging system consisting of one Standard Waste Box (SWB) loaded with drums placed inside a Standard Large Box 2 (SLB2). During transportation, the SLB2 is inside Transuranic Package Transporter Model III (TRUPACT-III), which provides the third layer of the packaging. Once at the WIPP, it is assumed the SLB2 is extracted from the TRUPACT-III and maintained above ground, and then subsequently placed underground for permanent disposal. In these proposed configurations, the space between the SLB2 and the SWB is always filled by a layer of insulation consisting of air-filled glass microbubbles except for the bottom which rests directly on the SLB2. The thermal scoping calculations described in this report specifically address whether the introduction of external heat inputs, combined with the contributions from the internally generated radiolytic decay heat and chemical reactions, lead to an unstable thermal state during the time of its movement and placement in the permanent disposal location. The external heat inputs are of two forms: 1) ambient thermal irradiation (e.g., solar and ambient storage/disposal temperatures) and 2) accident-induced fire. Three scoping calculation scenarios were derived as representative, conservative scenarios: 1A) TRUPACT-III transient transportation, 1B) SLB2 48-hour outdoor storage with solar radiation, and 2) fully-engulfing fire during SLB2 handling or emplacement following a steady-state analysis in a 38 °C environment. All the simulated scenarios are conservative relative to the operational conditions expected for handling the waste package during transportation and placement in the WIPP underground disposal unit. The predictions obtained from simulating the three exposure scenarios revealed that adding the SLB2 and the air-filled glass microbubbles to the transport and storage/disposal configurations provides additional thermal protection of the drums beyond what the SWB provides alone, both during long-term above ground insolation and underground during a fire accident. Under the current transportation/storage/disposal concepts, the degree of protection provided by the packaging concept is sufficient to prevent the waste from being ignitable. The simulation results demonstrate that there is adequate margin to safely transport and place the RNS waste from WCS to the WIPP under the current operational concept.

The Photovoltaic (PV) Performance Modeling Collaborative (PVPMC) organized a blind PV performance modeling intercomparison to allow PV modelers to blindly test their models and modeling ability against real system data. Measured weather and irradiance data were provided along with detailed descriptions of PV systems from two locations (Albuquerque, New Mexico, USA, and Roskilde, Denmark). Participants were asked to simulate the plane-of-array irradiance, module temperature, and DC power output from six systems and submit their results to Sandia for processing. The results showed overall median mean bias (i.e., the average error per participant) of 0.6% in annual irradiation and −3.3% in annual energy yield. While most PV performance modeling results seem to exhibit higher precision and accuracy as compared to an earlier blind PV modeling study in 2010, human errors, modeling skills, and derates were found to still cause significant errors in the estimates.

Abstract not provided.

Infrasound observations have grown increasingly important for the monitoring of earthquakes. While large earthquakes generate infrasound that can be detected thousands of kilometers away, there are few near-field observations of infrasound generated by low-magnitude events. We describe preliminary results of the West Texas Acoustic Experiment, during which infrasound sensors collected continuous data in the Permian Basin for a six-month period spanning January—June 2023. During this time, more than 1000 earthquakes with magnitudes between 1.2 and 4.2 occurred within 50 km of the network. We used spectral analysis, array processing, and manual inspection of waveforms to evaluate arrivals of infrasound signals following 84 events with magnitudes between 2.5 and 4.2. Here, we describe eight such events and the infrasound signals associated with each. We find detections of seismic-to-acoustic infrasound signals associated with seven events. We also find strong evidence of a laterally-propagating, purely acoustic wave generated by an M2.9 earthquake.

Abstract not provided.

Abstract not provided.

Physical Chemistry Chemical Physics

This report captures the results of development and testing of a integral downhole motor and percussive hammer used for drilling in near-surface hard rock formations. The work was funded through the DOE Office of Technology Transitions Technology Commercialization Fund. It was a collaboration between Sandia National Labs and The Charles Machine Works (aka Ditch Witch). In the collaboration, Sandia developed a pneumatic motor derived from an indexing tool used in other drilling applications, and Ditch Witch developed the bearing pack tied to the output shaft of the motor as well as the angled beacon housing used for directional control.

Fast and reliable estimation of engineered fracture geometries is a key factor in controlling undesirable fractures and enhancing stimulation design. Measuring the surface deformation gradient (tilt) for engineered fractures in shallow depths (<1000 m) has been proven a reliable source of data to infer fracture geometry, thanks to the impressive resolution of tiltmeter units (in the order of nano-radians). However, solving the inverse problem requires reliable and fast forward models. In this study, we present a fast and reliable machine-learned surrogate model to estimate the ground surface tilt induced by pressurised fractures. The proposed surrogate model, based on Conditional Generative Adversarial Networks (cGAN), receives a fracture aperture map in XY and XZ planes as input and predicts the corresponding surface tilts (in X and Y directions). The surrogate model with Wasserstein loss and gradient penalty has been trained using 11,000 samples and tested for a range of input parameters such as depth, dip angles, elastic properties, fluid pressures and fracture shapes. The testing results show excellent performance of the surrogate model compared with the forward finite element model for both single and multiple pressurised fractures, while running hundreds to potentially thousands of times faster.

Ultimately, our experiment measures two quantities on an aluminum bar: motion (which modeling must predict) and temperature (which sets thermal boundary conditions). For motion, stereo DIC is a technique to use imaging data to provide displacements relative to a reference image down to 1/100th of a pixel. We use a calibrated infrared imaging method for accurate temperature measurements. We will be capturing simultaneous data and then registering temperature data in space to the same coordinate system as the displacement data. While we will later show that our experiments are repeatable, indicating that separate experiments for motion and temperature would provide similar data, the simultaneous and registered data removes test to test variability as a source of uncertainty for model calibration and reduces the number of time-consuming tests that must be performed.

Abstract not provided.

Abstract not provided.

The U.S. Strategic Petroleum Reserve (SPR) is a crude oil storage system administered by the U.S. Department of Energy. SPR injected a total of over 230 MMB of raw water into 48 caverns as part of oil sales in CY22. Leaching effects were monitored in these caverns to understand how the sales operations may impact the long-term integrity of the caverns. The leaching effects were modeled here using the Sandia Solution Mining Code, SANSMIC. The modeling results indicate that leaching-induced features do not raise concern for the majority of the caverns. In addition to 12 caverns identified in previous leaching reports, seven caverns have been identified for further monitoring based on the results of this report. Twenty-two caverns had pre- and post-leach sonars that were compared with SANSMIC results. Overall, SANSMIC was able to capture the leaching well.

Mechanical metamaterials are artificial materials with unique global properties due to the structural geometry and material composition of their unit cell. Typically, mechanical metamaterial unit cells are designed such that, when tessellated, they exhibit unique mechanical properties such as zero or negative Poisson's ratio and negative stiffness. Beyond these applications, mechanical metamaterials can be used to achieve tailorable nonlinear deformation responses. Computational methods such as gradient-based topology optimization (TO) and size/shape optimization (SSO) can be implemented to design these metamaterials. However, both methods can lead to suboptimal solutions or a lack of generalizability. Therefore, this research used deep reinforcement learning (DRL), a subset of deep machine learning that teaches an agent to complete tasks through interactive experiences, to design mechanical metamaterials with specific nonlinear deformation responses in compression or tension. The agent learned to design the unit cells by sequentially adding material to a discrete design domain and being rewarded for achieving the desired deformation response. After training, the agent successfully designed unit cells to exhibit desired deformation responses not experienced during training. This work shows the potential of DRL as a high-level design tool for a wide array of engineering applications.

Single particle levitation methods are a powerful subset of aerosol instrumentation that allow a wide range of particle properties and processes to be explored. One of the most common forms of single particle levitation uses electric fields and is generally referred to as an electrodynamic balance (EDB). There are many different kinds of EDB's that have been designed with different applications in mind, and a corresponding array of analytical tools have been developed to characterize particles held in these traps. In this tutorial, we review the design and development of the EDB and discuss a range of analytical methods, including electrostatic analysis, light scattering, spectroscopy, and imaging, that allow for measurements of hygroscopic growth, volatility, surface tension and viscosity, diffusion, and phase and morphology. We go on to review recent advanced analytical methods using mass spectrometry to probe particle composition. This review is intended to provide readers with the basic knowledge to set up an EDB platform, design measurement protocols based on the available analytical tools, and run experiments to probe the fundamental properties of aerosol particles relevant to their role in the atmosphere, impacts on clouds and climate, effects on air quality, role in health and disease, and applications in industrial processes.

Uncertainty in severe accident evolution and outcome is driven by event bifurcations that represent distinctive challenges to defensive layers and tend to promote the emergence of discrete classes of core damage and accident risk. This discrete set of "attractor" states arise from the complex networks of competing physical phenomena and conditional event cascades occurring as the overall system degrades – a process that yields increasing degrees of freedom and accident progression pathways. Characterization of these event spaces has proven elusive to more traditional data interrogation methods, but proves tractable by application of more advanced data collection and machine learning approaches. Through application of these approaches we demonstrate a conceptual framework that enables real-time/robust, risk-informed decision-making support to improve accident mitigation and encourage “graceful exits” during low probability, extreme events limiting accident consequences. In this analysis, we simulated over 8,000 short-term station blackout (STSBO) accidents with the state-of-the-art integral severe accident code, MELCOR, and demonstrate the potential for ML approaches to predict simulation outcomes. We chose to pair ML tools with interpretable and mechanistic event trees for the considered STSBO accident space to predict the likelihood of future event paths along the tree. In addition to the current state of the system, we use information from recent trajectories of temperature, pressure, and other physical features, combining both the current state and past trajectories to forecast future event paths. Finally, we simulate the random injection of variable amounts of water to quantify the efficacy of available actions at reducing risks along the many branches in the event tree. We identify scenarios and windows of opportunity to mitigate risk as well as scenarios in which such actions are unlikely to alter the accident end-state.

Aquaculture systems require careful consideration of location, which determines water conditions, pollution impacts, and hazardous conditions. Mobility may be able to address these factors while also supporting the targeting of renewable energy sources such as wind, wave, and solar power throughout the year. In this paper, a purpose-built mobile aquaculture ship is identified and modeled with a combination of renewable energy harvesting capabilities as a case study with the objective of assessing the potential benefits of targeting high renewable energy potentials to power aquaculture operations. A route optimization algorithm is created and tuned to simulate the mobility of the aquaculture platform and cost-basis comparisons are made to a stationary system. The small spatial variability in renewable energy potential when combining multiple resources significantly limits the benefits of a mobile, renewable-targeting aquaculture system. On the other hand, the consistent energy harvest from a blend of renewable energy types (13 kW installed wind capacity, 661 m2 installed solar, and 1 m characteristic width wave-energy converter) suggests that the potential benefits of a mobile platform for offshore aquaculture (mitigation of environmental and social concerns, any potential positive impact on yields, hazard avoidance, etc.) can likely be pursued without significant increases in energy harvester costs.

Abstract not provided.

The Ghareb Formation is a shallowly buried porous chalk in southern Israel that is being considered as a host rock for a geologic nuclear waste repository. Setup and operation of a repository will induce significant mechanical, hydrological and chemical perturbations in the Ghareb. Developing a secure repository requires careful characterization of the rock behavior to different loads. To characterize hydromechanical behavior of the Ghareb, several short- and long-term deformation experiments were conducted. Hydrostatic loading tests were conducted both dry and water-saturated, using different setups to measure elastic properties, time-dependent behavior, and permeability. A set of triaxial tests were conducted to measure the elastic properties and rock strength under differential loading at dry and water-saturated conditions. The hydrostatic tests showed the Ghareb began to deform inelastically around 12–15 MPa, a relatively low effective pressure. Long-term permeability measurements demonstrated that permeability declined with increasing effective pressure and was permanently reduced by ~ 1 order of magnitude after unloading pressure. Triaxial tests showed that water saturation significantly degrades the rock properties of the Ghareb, indicating water-weakening is a significant risk during repository operation. Time-dependent deformation is observed during hold periods of both the hydrostatic and triaxial tests, with deformation being primarily visco-plastic. The rate of deformation and permeability loss is strongly controlled by the effective pressure as well. Additionally, during holds of both hydrostatic and triaxial tests, it is observed that when water-saturated, radial strain surpassed axial strain when above effective pressures of 13–20 MPa. Thus, deformation anisotropy may occur in situ during operations even if the stress conditions are hydrostatic when above this pressure range.

Peat fires are a major contributor to greenhouse gas emissions. The estimates of these emissions currently contain major uncertainties, due to the difficulty of determining the mass of peat burned in a fire. To address these uncertainties, we develop a computational physics-based peat smoldering model, which will be leveraged for high-fidelity quantitative estimates of peat fire emissions relevant to climate change. We present the verification of the 2-D axisymmetric model, a first step towards developing a full 3-D model. Verification includes the solution verification against a literature model for the 0-D smoldering case and verification of the heat transfer problem in 1-D and 2-D. Also presented is the effect of reaction mechanism on the smoldering model, for which we found a relatively simple three-step reaction mechanism is able to capture key behavior. These verification results provide the foundation for moving forward with validation against experimental data of the 2-D model.

We present methods to estimate parameters for models for the incidence angle modifier for simulating irradiance on a photovoltaic array. The incidence angle modifier quantifies the fraction of direct irradiance that is reflected away at the array’s face, as a function of the direct irradiance’s angle of incidence. Parameters can be estimated from data and the fitting method can be used to convert between models. We show that the model conversion procedure results in models that produce similar annual insolation on a fixed plane.

This special memorial issue pays tribute to James (Jim) A. Miller, a giant of combustion science who died in 2021, with a celebration of his enormous influence on the field. We were touched by the responses we received after we sent out the invitations for it. Jim inspired several generations of scientists, who viewed him as a mentor, a father figure, and a friend. Together with Nils Hansen and Peter Glarborg, we have written a detailed account on his life and work. Furthermore, it appeared in this journal shortly after his death; and so here we focus on the scientific areas he had interest in and influence on, and how they relate to the 34 papers in this issue. The topics of these papers span a variety of Jim's interests including nitrogen chemistry, polycyclic aromatic hydrocarbon (PAH) chemistry, oxidation chemistry, energy transfer, prompt dissociations, and codes to facilitate combustion chemistry simulations.

Abstract not provided.

Abstract not provided.

The Daedalus ultrafast x-ray imager is the latest generation in Sandia’s hybrid CMOS detector family. With three frames along an identical line of sight, 1 ns minimum integration time, a higher full well than Icarus, and added features, Daedalus brings exciting new capabilities to diagnostic applications in inertial confinement fusion and high energy density science. In this work, we present measurements of time response, dynamic range, spatial uniformity, pixel cross-talk, and absolute x-ray sensitivity using pulsed optical and x-ray sources. We report a measured 1.5 Me− full well, pixel sensitivity at 9.58 × 10−7 V/e−, and an estimate of spatial uniformity at ∼5% across the sensor array.

Marine aerosol injections are a key component in further understanding of both the potentials of deliberate injection for marine cloud brightening (MCB), a potential climate intervention (CI) strategy, and key aerosol-cloud interaction behaviors that currently form the largest uncertainty in global climate model (GCM) predictions of our climate. Since the rate of spread of aerosols in a marine environment directly translates to the effectiveness and ability of aerosol injections in impacting cloud radiative forcing, it is crucial to understand the spatial and temporal extent of injected-aerosol effects following direct injection into marine environments. The ubiquity of ship-injected aerosol tracks from satellite imagery renders observational validation of new parameterizations possible in 2D, however, 3D compatible data is more scarce, and necessary for the development of subgrid scale parameterizations of aerosol-cloud interactions in GCMs. This report introduces two novel parameterizations of atmospheric aerosol injection behavior suitable for both 3D (GCM-compatible) and 2D (observation-related) modeling. Their applicability is highlighted using a wealth of different observational data: small and larger scale salt-aerosol injection experiments conducted at SNL, 3D large eddy simulations of ship-injected aerosol tracks and 2D satellite images of ship tracks. The power of experimental data in enhancing knowledge of aerosol-cloud interactions is in particular emphasized by studying key aerosol microphysical and optical properties as observed through their mixing in cloud-like environments.

Moving target defenses (MTDs) are widely used as an active defense strategy for thwarting cyberattacks on cyber-physical systems by increasing diversity of software and network paths. Recently, machine Learning (ML) and deep Learning (DL) models have been demonstrated to defeat some of the cyber defenses by learning attack detection patterns and defense strategies. It raises concerns about the susceptibility of MTD to ML and DL methods. In this article, we analyze the effectiveness of ML and DL models when it comes to deciphering MTD methods and ultimately evade MTD-based protections in real-time systems. Specifically, we consider a MTD algorithm that periodically randomizes address assignments within the MIL-STD-1553 protocol - a military standard serial data bus. Two ML and DL-based tasks are performed on MIL-STD-1553 protocol to measure the effectiveness of the learning models in deciphering the MTD algorithm: 1) determining whether there is an address assignments change i.e., whether the given system employs a MTD protocol and if it does 2) predicting the future address assignments. The supervised learning models (random forest and k-nearest neighbors) effectively detected the address assignment changes and classified whether the given system is equipped with a specified MTD protocol. On the other hand, the unsupervised learning model (K-means) was significantly less effective. The DL model (long short-term memory) was able to predict the future addresses with varied effectiveness based on MTD algorithm's settings.

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

Abstract not provided.

Here, a line list for the N second positive system, $B^3Π_g—C^3Π_u$, has been compiled using the PGOPHER spectral simulation software. The line list extends the number of vibrational states of the $B^3Π_g$ up to v=29 and a maximum rotational state of J=150 for simulation temperatures up to 7000 K. New electronic–vibrational transition moments were calculated using refined potential energy curves and a transition dipole moment with the DUO software. Comparisons to experimental data and the SPECAIR software have been used to validate the new line list. The results are available in ASCII ExoMol .state and .trans files and as a PGOPHER input file for use in spectral analysis.

Commercial generation of energy by nuclear power plants in the United States (U.S.) has produced thousands of metric tons of spent nuclear fuel (SNF), the disposal of which is the responsibility of the U.S. Department of Energy (DOE). Utilities typically utilize the practice of storing this SNF in dual-purpose canisters (DPCs). DPCs were designed, licensed, and loaded to meet Nuclear Regulatory Commission (NRC) requirements that preclude the possibility of a criticality event during SNF storage and transport, but were not designed or loaded to preclude the possibility of a criticality event during the regulated post-closure period following disposal, which could be up to 1,000,000 years (Price, 2019). There are several options being investigated that could facilitate the disposal of SNF stored in DPCs in a geologic repository (Hardin et al., 2015; SNL 2020b; SNL 2021b). These include: (1) repackage the SNF into canisters that are designed to prevent criticality during the regulated post-closure period following disposal, but with an increased disposal cost estimated at approximately $\$$20B in United States dollars (USD) (Freeze et al., 2019); (2) analysis of the probability and consequences of criticality from the direct disposal of DPCs during a 1,000,000-year post-closure period in several geologic disposal media (Price, 2019); and (3) filling the void space of a DPC with a material before its disposal that significantly limits the potential for criticality over the post-closure regulatory period. This report further investigates the third option, filling DPC already containing SNF with a material to limit the potential for criticality over the post-closure regulatory period.

The Clifford spectrum is a form of joint spectrum for noncommuting matrices. This theory has been applied in photonics, condensed matter and string theory. In applications, the Clifford spectrum can be efficiently approximated using numerical methods, but this only is possible in low dimensional example. In this paper we examine the higher-dimensional spheres that can arise from theoretical examples. We also describe a constructive method to generate five real symmetric almost commuting matrices that have a K-theoretical obstruction to being close to commuting matrices. For this, we look to matrix models of topological electric circuits.

Opacity-on-NIF has obtained opacity data under conditions similar to those achieved by the entirely different Opacity-on-Z platform. From low- and high-Z elements at different anchor points, rigorously compare the opacity data between the laboratories and to multiple opacity theory models. Compare and assess the data acquisition and processing methods for obtaining opacities and for measuring/inferring sample conditions. Explain, or develop hypotheses for, any discrepancies. Map progress to the National Opacity Strategy and define future directions.

We report a spontaneous and hierarchical self-assembly mechanism of carbon dots prepared from citric acid and urea into nanowire structures with large aspect ratios (>50). Scattering-type scanning near-field optical microscopy (s-SNOM) with broadly tunable mid-IR excitation was used to interrogate details of the self-assembly process by generating nanoscopic chemical maps of local wire morphology and composition. s-SNOM images capture the evolution of wire formation and the complex interplay between different chemical constituents directing assembly over the nano- to microscopic length scales. We propose that residual citrate promotes tautomerization of melamine surface functionalities to produce supramolecular shape synthons comprised of melamine-cyanurate adducts capable of forming long-range and highly directional hydrogen-bonding networks. This intrinsic, heterogeneity-driven self-assembly mechanism reflects synergistic combinations of high chemical specificity and long-range cooperativity that may be harnessed to reproducibly fabricate functional structures on arbitrary surfaces.

Nonlocal models allow for the description of phenomena which cannot be captured by classical partial differential equations. The availability of efficient solvers is one of the main concerns for the use of nonlocal models in real world engineering applications. Here, we present a domain decomposition solver that is inspired by substructuring methods for classical local equations. In numerical experiments involving finite element discretizations of scalar and vectorial nonlocal equations of integrable and fractional type, we observe improvements in solution time of up to 14.6x compared to commonly used solver strategies.

The heat generated by high-level radioactive waste can pose numerical and physical challenges to subsurface flow and transport simulators if the liquid water content in a region near the waste package approaches residual saturation due to evaporation. Here, residual saturation is the fraction of the pore space occupied by liquid water when the hydraulic connectivity through a porous medium is lost, preventing the flow of liquid water. While conventional capillary pressure models represent residual saturation using asymptotically large values of capillary pressure, here, residual saturation is effectively modeled as a tortuosity effect alone. Treating the residual fluid as primarily dead-end pores and adsorbed films, relative permeability is independent of capillary pressure below residual saturation. To test this approach, PFLOTRAN is then used to simulate thermal-hydrological conditions resulting from direct disposal of a dual-purpose canister in unsaturated alluvium using both conventional asymptotic and revised, smooth models. Importantly, while the two models have comparable results over 100 000 years, the number of flow steps required is reduced by approximately 94%.

The basic building block of a distributed-memory cluster or supercomputer is a node. Each node includes a host, which is a processor (xPU) + memory hierarchy. The host can communicate with other hosts via its NIC (network interface controller). A network connects the nodes. The nodes may be arranged in some topology, which determines the network’s carrying capacity and cost.

The finite element method (FEM) is widely used to simulate a variety of physics phenomena. Approaches that integrate FEM with neural networks (NNs) are typically leveraged as an alternative to conducting expensive FEM simulations in order to reduce the computational cost without significantly sacrificing accuracy. However, these methods can produce biased predictions that deviate from those obtained with FEM, since these hybrid FEM-NN approaches rely on approximations trained using physically relevant quantities. In this work, an uncertainty estimation framework is introduced that leverages ensembles of Bayesian neural networks to produce diverse sets of predictions using a hybrid FEM-NN approach that approximates internal forces on a deforming solid body. The uncertainty estimator developed herein reliably infers upper bounds of bias/variance in the predictions for a wide range of interpolation and extrapolation cases using a three-element FEM-NN model of a bar undergoing plastic deformation. This proposed framework offers a powerful tool for assessing the reliability of physics-based surrogate models by establishing uncertainty estimates for predictions spanning a wide range of possible load cases.

The assembly of ultra-complex structures from simple building units remains a long-term challenge in chemistry. Using small molecular building blocks (MBBs) in a mixed-ligand approach permitted the assembly of unprecedented metal-organic frameworks (MOFs), M-kum-MOF-1 (M = Y, Tb), exhibiting extra-large mesoporous cavities with small access windows. The ultra-complex cage of M-kum-MOF-1 consists of 240 vertices bridged by 432 edges, leading to a 194 faces-containing tile. This tile exhibits more faces than in any periodic structures (zeolites, MOFs, metal-organic polyhedra [MOPs], etc.) known to date. M-kum-MOF-1 not only possess zeolitic features (anionic framework), but they also contain an underlying wse zeolitic topology, which is observed for the first time.

Auto-magnetizing (AutoMag) liners are cylindrical tubes that employ helical current flow to produce strong internal axial magnetic fields prior to radial implosion on ~100 ns timescales. AutoMag liners have demonstrated strong uncompressed axial magnetic field production (>100 T) and remarkable implosion uniformity during experiments on the 20 MA Z accelerator. However, both axial field production and implosion morphology require further optimization to support the use of AutoMag targets in magnetized liner inertial fusion (MagLIF) experiments. Data from experiments studying the initiation and evolution of dielectric flashover in AutoMag targets on the Mykonos accelerator have enabled the advancement of magnetohydrodynamic (MHD) modeling protocols used to simulate AutoMag liner implosions. Implementing these protocols using ALEGRA has improved the comparison of simulations to radiographic data. Specifically, both the liner in-flight aspect ratio and the observed width of the encapsulant-filled helical gaps during implosion in ALEGRA simulations agree more closely with radiography data compared to previous GORGON simulations. Although simulations fail to precisely reproduce the measured internal axial magnetic field production, improved agreement with radiography data inspired the evaluation of potential design improvements with newly developed modeling protocols. Three-dimensional MHD simulation studies focused on improving AutoMag target designs, specifically seeking to optimize the axial magnetic field production and enhance the cylindrical implosion uniformity for MagLIF. Importantly, by eliminating the driver current prepulse and reducing the initial inter-helix gap widths in AutoMag liners, simulations indicate that the optimal 30–50 T range of precompressed axial magnetic field for MagLIF on Z can be accomplished concurrently with improved cylindrical implosion uniformity.

Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the. U.S. Department of Energy’s National Nuclear Security Administration. The National Nuclear Security Administration’s Sandia Field Office administers the contract and oversees contractor operations at the Sandia National Laboratories Kaua'i Test Facility in Hawai'i. Activities at the site are conducted in support of U.S. Department of Energy weapons programs., and the site has operated as a rocket preparation launching and tracking facility since 1962. The U.S. Department of Energy and its management and operating contractor are committed to safeguarding the environment, assessing sustainability practices, and ensuring the validity and accuracy of the monitoring data presented in this annual site environmental report. This report summarizes the environmental protection, restoration, and monitoring programs in place at Sandia National Laboratories, Kaua'i Test Facility, during calendar year 2022. Environmental topics include cultural resource management, chemical management, air quality, meteorology, ecology, oil storage, site sustainability, terrestrial surveillance, waste management, water quality, wastewater discharge, and implementation of the National Environmental Policy Act. This report is prepared in accordance with and as required by DOE O 231.1B, Admin Change 1, Environment, Safety and Health Reporting and has been approved for public distribution.

Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the. U.S. Department of Energy’s National Nuclear Security Administration. The National Nuclear Security Administration’s Sandia Field Office administers the contract and oversees contractor operations at Sandia National Laboratories, Tonopah Test Range. Activities at the site are conducted in support of U.S. Department of Energy weapons programs and have operated at the site since 1957. The U.S. Department of Energy and its management and operating contractor are committed to safeguarding file environment, assessing sustainability practices, and ensuring the validity and accuracy of the monitoring data presented in this annual site environmental report. This report summarizes the environmental protection, restoration, and monitoring programs in place at Sandia National Laboratories, Tonopah Test Range during calendar year 2022. Environmental topics include cultural resource management, chemical management, air quality, ecology, environmental restoration, oil storage, site sustainability, terrestrial surveillance, waste management, water quality, wastewater discharge, and implementation of the National Environmental Policy Act. This report is prepared in accordance with and as required by DOE 0 231.IB, Admin Change 1, Environment, Safety and Health Reporting and has been approved for public distribution.

Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration. The National Nuclear Security Administration’s Sandia Field Office administers the contract and oversees contractor operations at Sandia National Laboratories, California. Activities at this multiprogram engineering and science laboratory support the nuclear weapons stockpile program, energy and environmental research, homeland security, micro- and nanotechnologies, and basic science and engineering research. The U.S. Department of Energy and its management and operating contractor are committed to safeguarding the environment, assessing sustainability practices, and ensuring the validity and accuracy of the monitoring data presented in this annual site environmental report. This report provides a summary of environmental monitoring information and compliance activities that occurred at Sandia National Laboratories, California during calendar year 2022 unless noted otherwise. General site and environmental program information is also included. This report was prepared in accordance with DOE O 231.1B, Environment, Safety and Health Reporting.

Sandia National Laboratories is a multimission laboratory managed and operated by National Technology & Engineering Solutions of Sandia, LLC, a wholly owned subsidiary of Honeywell International Inc., for the U.S. Department of Energy’s National Nuclear Security Administration. The National Nuclear Security Administration’s Sandia Field Office administers the contract and oversees contractor operations at Sandia National Laboratories, New Mexico. Activities at the site support research and development programs with a wide variety of national security missions, resulting in technologies for nonproliferation, homeland security, energy and infrastructure, and defense systems and assessments. The U.S. Department of Energy and its management and operating contractor are committed to safeguarding the environment, assessing sustainability practices, and ensuring the validity and accuracy of the monitoring data presented in this annual site environmental report. This report summarizes the environmental protection and monitoring programs in place at Sandia National Laboratories, New Mexico, during calendar year 2022. Environmental topics include cultural resource management, chemical management, air quality, ecology, environmental restoration, oil storage, site sustainability, terrestrial surveillance, waste management, water quality, and implementation of the National Environmental Policy Act. This report is prepared in accordance with and as required by DOE O 231.1B, Admin Change 1, Environment, Safety and Health Reporting, and has been approved for public distribution.

There is a common need in the advancement of optical diagnostic techniques to increase the dimensionality of measurements. For example, point measurements could be improved to multi-point, line, planar, volumetric, or time-resolved volumetric measurements. In this work, a unique optical element is presented to enable multidimensional measurements, namely, an array of glass wedges. A light source is passed through the wedges, and different portions of the illumination are refracted by different amounts depending on the glass wedge angle. Subsequent optics can be used to focus the light to multiple points, lines, or planes. Basic characterization of a glasswedge array is presented. Additionalwedge-array configurations are discussed, including the use of a periodic intensity mask for multi-planar measurements via structured illumination. The utility of this optical element is briefly demonstrated in (a) multi-planar flame particulate measurements, (b) multi-point femtosecond-laser electronic excitation tagging for flow velocimetry, and (c) multi-line nitric oxide molecular tagging velocimetry in a hypersonic shock-tunnel. One significant advantage of this optical component is its compatibility with highenergy laser sources, which may be a limiting factor with other beam-splitting or beam-forming elements such as some diffractive optics. Additionally, an array of glass wedges is simple and easily customizable compared to other methods for forming multiple closely spaced illumination patterns. Suggestions for further development and applications are discussed.

The interactions of carboxylate anions with water and cations are important for a wide variety of systems, both biological and synthetic. Here, in order to gain insight on properties of the local complexes, we apply density functional theory, to treat the complex electrostatic interactions, and investigate mixtures with varied numbers of carboxylate anions (acetate) and waters binding to monovalent cations, Li⁺, Na⁺ and K⁺. The optimal structure with overall lowest free energy contains two acetates and two waters such that the cation is four-fold coordinated, similar to structures found earlier for pure water or pure carboxylate ligands. More generally, the complexes with two acetates have the lowest free energy. In transitioning from the overall optimal state, exchanging an acetate for water has a lower free energy barrier than exchanging water for an acetate. In most cases, the carboxylates are monodentate and in the first solvation shell. As water is added to the system, hydrogen bonding between waters and carboxylate O atoms further stabilizes monodentate structures. These structures, which have strong electrostatic interactions that involve hydrogen bonds of varying strength, are significantly polarized, with ChelpG partial charges that vary substantially as the bonding geometry varies. Overall, these results emphasize the increasing importance of water as a component of binding sites as the number of ligands increases, thus affecting the preferential solvation of specific metal ions and clarifying Hofmeister effects. Finally, structural analysis correlated with free energy analysis supports the idea that binding to more than the preferred number of carboxylates under architectural constraints are a key to ion transport.

This milestone work baselines electromagnetic particle-in-cell capability of the EMPIRE plasma simulation code to model key processes germane to the physics of electrode plasmas arising in magnetically-insulated transmission lines operating at or near 20 MA. This evaluation is done so through the provision of benchmark verification problems designed to exercise the individual and combined physics models on a small-scale surrogate geometry for the final-feed-to-load region of the Z accelerator under representative operating conditions. In this report, we overview our test designs, and present a portfolio of simulation results along with performance assessments which altogether establish state-of-the-art. In particular, two main verification categories are covered this report: (1) Z-relevant desorption physics (Temkin isotherm), and (2) two approaches to simulate electrode plasma creation and dynamics (automatic creation versus self-consistent creation through direct simulation Monte Carlo collisions).

The natural assemblage of a symbiotic bacterial microbiome (bacteriome) with microalgae in marine ecosystems is now being investigated as a means to increase algal productivity for industry. When algae are grown in open pond settings, biological contamination causes an estimated 30% loss of the algal crop. Therefore, new crop protection strategies that do not disrupt the native algal bacteriome are needed to produce reliable, high-yield algal biomass. Bacteriophages offer an unexplored solution to treat bacterial pathogenicity in algal cultures because they can eliminate a single species without affecting the bacteriome. To address this, we identified a highly virulent pathogen of the microalga Nannochloropsis gaditana, the bacterium Bacillus safensis, and demonstrated rescue of the microalgae from the pathogen using phage. 16S rRNA amplicon sequencing showed that phage treatment did not alter the composition of the bacteriome. It is widely suspected that the algal bacteriome could play a protective role against bacterial pathogens. To test this, we compared the susceptibility of a bacteriome-attenuated N. gaditana culture challenged with B. safensis to a N. gaditana culture carrying a growth-promoting bacteriome. We showed that the loss of the bacteriome increased the susceptibility of N. gaditana to the pathogen. Transplanting the microalgal bacteriome to the bacteriome-attenuated culture reconstituted the protective effect of the bacteriome. Finally, the success of phage treatment was dependent on the presence of beneficial bacteriome. This study introduces two synergistic countermeasures against bacterial pathogenicity in algal cultures and a tractable model for studying interactions between microalgae, phages, pathogens, and the algae microbiome.

The Python-Cubit^® enhancement code base is intended to be used as an extension to already existing Cubit^® functionality. It provides the user with a number of functionalities that are either currently outside the realm of the python functions which Cubit^® supplies internally (such as vector math), or that are comprised of commonly used combinations of already existing python functionalities (such as removing a full round from a slot cut). The foreseen style of use for many of these scripts is to utilize volume names and geometric data such as surface area, surface type, etc. as a way to filter out geometries, and provide a powerful id-less method. These filters combined with a number of already existing python functionalities such as the set() operator and zip() function can be used to operate on many geometries at a single time without a need for the user to manually select them or use their ids. Please refer to the example given in the documents examples section for a demonstration of the work flow.

The dislocation core structure has a significant role in determining the dominant slip plane and the magnitude of the Peierls stress for a dislocation. An important challenge when studying dislocation cores is to determine the stable and metastable core morphologies, and then relate these structures to the dynamics of the dislocations. ere this study introduces a method for identifying core structures that are metastable at zero temperature. Application of this method to $\langle$a$\rangle$-type screw dislocations in α-Ti (as described using an empirical potential) reveals a multitude of (meta)stable nonplanar cores. Molecular dynamics studies show how the competing metastable core structures determine the properties of the dislocations at temperature and under a range of non-Schmid stresses.

There are growing needs to understand how extreme weather events impact the electrical grid. Renewable energy sources such as solar photovoltaics are expanding in use to help sustainably meet electricity demands. Wildfires and, notably, the widespread smoke resulting from them, are one such extreme event that can impair the performance of solar photovoltaics. However, isolating the impact that smoke has on photovoltaic energy production, separate from ambient conditions, can be difficult. In this work, we seek to understand and quantify the impacts of wildfire smoke on solar photovoltaic production within the Western United States. Our analysis focuses on the construction of a random forest regression model to predict overall solar photovoltaic production. The model is used to separate and quantify the impacts of wildfire smoke in particular. To do so, we fuse historical weather, solar photovoltaic energy production, and PM2.5 particulate matter (primary smoke pollutant) data to train and test our model. The additional weather data allows us to capture interactions between wildfire smoke and other ambient conditions, as well as to create a more powerful predictive model capable of better quantifying the impacts of wildfire smoke on its own. We find that solar PV energy production decreases 8.3% on average during high smoke days at PV sites as compared to similar conditions without smoke present. This work allows us to improve our understanding of the potential impact on photovoltaic-based energy production estimates due to wildfire events and can help inform grid and operational planning as solar photovoltaic penetration levels continue to grow.

The Xyce™ Parallel Electronic Simulator has been written to support the simulation needs of Sandia National Laboratories’ electrical designers. Xyce™ is a SPICE-compatible simulator with the ability to solve extremely large circuit problems on large-scale parallel computing platforms, but also includes support for most popular parallel and serial computers.

A series of MD and DFT simulations were performed to investigate hydrogen self-clustering and retention in tungsten. Using a newly develop machine learned interatomic potential, spontaneous formation of hydrogen platelets was observed after implanting low-energy hydrogen into tungsten at high fluxes and temperatures. The platelets formed along low miller index orientations and neighboring tetrahedral and octahedral sites and could grow to over 50 atoms in size. High temperatures above 600 K and high hydrogen concentrations were needed to observe significant platelet formation. A critical platelet size of six hydrogen atoms was needed for long term stability. Platelets smaller than this were found to be thermally unstable within a few nanoseconds. To verify these observations, characteristic platelets from the MD simulations were simulated using large-scale DFT. DFT corroborated the MD results in that large platelets were also found to be dynamically stable for five or more hydrogen atoms. The LDOS from the DFT simulated platelets indicated that hydrogen atoms, particularly at the periphery of the platelet, were found to be at least as stable as hydrogen atoms in bulk tungsten. In addition, electrons were found to be localized around hydrogen atoms in the platelet itself and that hydrogen atoms up to 4.2 Å away within the platelet were found to share charge suggesting that the hydrogen atoms are interacting across longer distances than previously suggested. These results reveal a self-clustering mechanisms for hydrogen within tungsten in the absence of radiation induced or microstructural defects that could be a precursor to blistering and potentially explain the experimentally observed high hydrogen retention particularly in the near surface region.

This work presents a spectral micromechanical formulation for obtaining the full-field and homogenized response of elastic micropolar composites. The algorithm relies on a coupled set of convolution integral equations for the micropolar strains, where periodic Green’s operators associated with a linear homogeneous reference medium are convolved with functions of the Cauchy and couple stress fields that encode the material’s heterogeneity, as well as any potential material nonlinearity. Such convolution integral equations take an algebraic form in the reciprocal Fourier space that can be solved iteratively. In this vein, the fast Fourier transform (FFT) algorithm is leveraged to accelerate the numerical solution, resulting in a mesh-free formulation in which the periodic unit cell representing the heterogeneous material can be discretized by a regular grid of pixels in two dimensions (or voxels in three dimensions). For verification, the numerical solutions obtained with the micropolar FFT solver are compared with analytical solutions for a matrix with a dilute circular inclusion subjected to plane strain loading. The developed computational framework is then used to study length-scale effects and effective (micropolar) moduli of composites with various topological configurations.

The Z accelerator at Sandia National Laboratories conducts z-pinch experiments at 26 MA in support of DOE missions in stockpile stewardship, dynamic materials, fusion, and other basic sciences. Increasing the current delivered to the z-pinch would extend our reach in each of these disciplines. To achieve increases in current and accelerator efficiency, a fraction of Z’s shots are set aside for research into transmission-line power flow. These shots, with supporting simulations and theory, are incorporated into this Advanced Diagnostics milestone report. The efficiency of Z is reduced as some portion of the total current is shunted across the transmission-line gaps prior to the load. This is referred to as “current loss”. Electrode plasmas have long been implicated in this process, so the bulk of dedicated power-flow experiments are designed to measure the plasma environment. The experimental analyses are enhanced by simulations conducted using realistic hardware and Z voltage pulses. In the same way that diagnostics are continually being improved for sensitivity and resolution, the modeling capability is continually being improved to provide faster and more realistic simulations. The specifics of the experimental hardware, diagnostics, simulations, and algorithm developments are provided in this report. The combined analysis of simulation and data confirms that electrode plasmas have the most detrimental impact on current delivery. Experiments over the last three years have tested the theoretical current-loss mechanisms of enhanced ion current, plasma gap closure, and Hall-related current. These mechanisms are not mutually exclusive and may be coincident in the final feed as well as in upstream transmission lines. The final-feed geometries tested here, however, observe lower-density plasmas without dominant ion currents which is consistent with a Hall-related current. The picture of plasma formation and transport formed from experiment and simulation is informing hardware designs being fielded on Z now and being proposed for the Next-Generation Pulsed Power (NGPP) facility. In this picture, the strong magnetic fields that heat the electrodes above particle emission thresholds also confine the charged particles near the surface. Some portion of the plasmas thus formed is transported into the transmission-line gap under the force of the electric field, with aid from plasma instabilities. The gap plasmas are then transported towards the load by a cross-field drift, where they accumulate and contribute to a likely Hall-related cross-gap current. The achievements in experimental execution, model validation, and physical analysis presented in this report set the stage for continued progress in power flow and load diagnostics on Z. The planned shot schedule for Z and Mykonos will provide data for extrapolation to higher current to ensure the predicted performance and efficiency of a NGPP facility.

The United Sates Department of Energy (DOE) Generation 3 Concentrated Solar Power (CSP) program is interested in higher efficiency power systems at lower costs, potentially with systems utilizing chloride molten salts. Ternary chloride molten salts are corrosive and need to be held at high temperatures to achieve higher power system efficiencies. However, materials and cost of manufacturing of such a facility can be very expensive, particularly using exotic materials that are not always readily available. Materials that can withstand the harsh corrosive and thermal-mechanical environments of high-temperature molten salt systems (>700 ℃) are needed. High temperature systems offer greater thermodynamic efficiency but must also make cost efficient use of corrosion-resistant alloys. To ensure reliable high-performance operation for molten salt power plant designs confidence in materials compatibility with CSP Gen 3 halide salts must be established. This paper will present an analysis of Inconel 625 as an alternative to the costly Haynes 230 at 760℃ for 500 hours. Both metals were tested in an unaltered state as well as a homogenous weld. Each sample was weighed pre- and post-test, with a final composition analysis using Scanning Electron Microscopy (SEM) and Energy Dispersive X-Ray Spectroscopy (EDS). Preliminary findings suggest that Haynes 230 outperformed Inconel 625, but more research at longer durations, 1,000 hours will be required for full reliable assessment.

This paper summarizes findings from a small, mixed-method research study examining industry perspectives on the potential for new forms of automation to invigorate the concentrating solar power (CSP) industry. In Fall 2021, the Solar Energy Technologies Office (SETO) of the United States Department of Energy (DOE) funded Sandia National Laboratories to elicit industry stakeholder perspectives on the potential role of automated systems in CSP operations. We interviewed eleven CSP professionals from five countries, using a combination of structured and open comment response modes. Respondents indicated a preference for automated systems that support heliostat manufacturing and installation, calibration, and responsiveness to shifting weather conditions. This pilot study demonstrates the importance of engaging industry stakeholders in discussions of technology research and development, to promote adoptable, useful innovation.

This paper summarizes findings from a small, mixed-method research study examining industry perspectives on the potential for new forms of automation to invigorate the concentrating solar power (CSP) industry. In Fall 2021, the Solar Energy Technologies Office (SETO) of the United States Department of Energy (DOE) funded Sandia National Laboratories to elicit industry stakeholder perspectives on the potential role of automated systems in CSP operations. We interviewed eleven CSP professionals from five countries, using a combination of structured and open comment response modes. Respondents indicated a preference for automated systems that support heliostat manufacturing and installation, calibration, and responsiveness to shifting weather conditions. This pilot study demonstrates the importance of engaging industry stakeholders in discussions of technology research and development, to promote adoptable, useful innovation.

In this study, deformations inflicted in diamond could help design stronger materials.

A third-generation chloride salt tank system was designed for a 1 MWth pilot-scale system to be investigated at the National Solar Thermal Test Facility (NSTTF) in Albuquerque, NM, USA. This prototype Gen 3, concentrating solar power (CSP) system was designed to facilitate a minimum of 6 hrs. of thermal energy storage (TES) with operational nominal temperatures of 500°C and 720°C for a cold and hot tank respectively. For this investigation, the researchers developed steady and transient computational fluid mechanics (CFD) circulation models to assess thermal-fluid behavior within the tanks, and their respective interactions with environmental heat transfer. The models developed for this novel CSP system design included unique chloride molten salt thermodynamic properties and correlations. The results of this investigation suggest thermal gradients for the steady flow model less 1oC with overall circulation velocities as high as approximately 2.1 m/s. Higher steady flow rates of salt passing into and out of the tanks resulted in smaller thermal gradients than the slower flow rates as the molten salt mixes better (an increase of around 120% in the heat transfer coefficient) at the higher velocities associated with the higher flow rate. The port spacing of 3.85 m was found to have a highly uniform temperature distribution. For the unsteady model, nitrogen flow was found to become appreciably steady after approximately 10 minutes, and resultant molten salt flow was found to increase slowly as the overall salt level rose.

The need for reliable, cost-effective, utility scale energy storage that is universally applicable across different regions is becoming evident with the global transition towards non-polluting renewable energy resources. The operations and management of these energy storage technologies introduces a unique challenge that is inherently different from the conventional energy storage in the form of fossil fuel. The investigation into the business model, value proposition and economic viability of a utility scale thermal energy storage was part of a program sponsored by the United States Department of Energy, called Energy I-Corps. During this program, the project team reached out to a series of industry stakeholders to conduct interviews on the topic of thermal energy storage for utility scale power generation. Specific focus was placed on the business model based on the market needs in the context of the power grid in the United States. The utilization and re-use of infrastructure at existing thermo-electric power plants yielded the most viable business model for the implementation of the form of thermal energy storage discussed here.

The Sandia National Laboratories (SNL) National Solar Thermal Test Facility (NSTTF) conducted efficacy testing on a shut-off isolation valve for use with molten ternary chloride salt. A ball valve was tested under controlled N2 ullage gas pressure and connected with flanged fittings that featured a spiral-wound gasket. The valve assembly consisted of boronized nickel coated SS316 components, with design features that greatly reduce the cost of overall valve assembly. Testing results showed that the valve did not leak, and post-test analysis demonstrated that the ball, seat, packing, and body all survived both the heat loads and the relative corrosive environment. Spiral-wound gaskets for flanged connections used in the system also functioned nominally, with no leaks or signs of failures during post-test analysis. However, testing was ultimately forced to rapidly stop after testing between 500-530°C as the actuator used on the valve failed in the heat, preventing the valve from sealing in the closed position. In addition, salt plugs and salt vapor plating also prevented the test from continuing.

Falling particle receivers are a promising receiver design to couple with particle-based concentrating solar power to help meet future levelized cost of electricity targets in next generation systems. The thermal performance of receivers is critical to the economics of the overall system, and accurate models of particle receivers are necessary to predict the performance in all conditions. A model validation study was performed using falling particle receiver data recently collected at the National Solar Thermal Test Facility at Sandia National Laboratories. The particle outlet temperature, the thermal efficiency of the receiver, and the wind speed and direction around the receiver were measured in 26 steady-state experiments and compared to a corresponding receiver model. The results of this study showed improved agreement with the experimental data over past validation efforts but did not fully meet all predefined validation metrics. Future model improvements were identified to continue to strengthen the modeling capabilities.

Here, we develop a Stockmayer fluid model that accounts for the dielectric responses of polar solvents (water, MeOH, EtOH, acetone, 1-propanol, DMSO, and DMF) and NaCl solutions. These solvent molecules are represented by Lennard-Jones (LJ) spheres with permanent dipole moments and the ions by charged LJ spheres. The simulated dielectric constants of these liquids are comparable to experimental values, including the substantial decrease in the dielectric constant of water upon the addition of NaCl. Moreover, the simulations predict an increase in the dielectric constant when considering the influence of ion translations in addition to the orientation of permanent dipoles.

Sandia National Laboratories and the Institut de Radioprotection et de Sûreté Nucléaire have collaborated on the design and execution of a set of critical experiments that explore the effects of molybdenum in water-moderated fuel-rod arrays. The molybdenum was included as sleeves on some of the fuel rods in the critical experiment fuel arrays. Approach-to-critical experiments were performed on five configurations of fuel and molybdenum sleeves using the 7uPCX fuel in core hardware that set the triangular fuel rod pitch at 15.494 mm. The experiments are evaluated as benchmark critical experiments for the 2023 edition of the International Criticality Safety Benchmark Evaluation Project (ICSBEP) Handbook as LEU-COMP-THERM-111.

The Integrated Tiger Series (ITS) generates a database containing energy deposition data. This data, when stored on an Exodus file, is not typically suitable for analysis within Sierra Mechanics for finite element analysis. The its2sierra tool maps data from the ITS database to the Sierra database. This document provides information on the usage of its2sierra.

The purpose of this proposal is to design a new integral critical experiment to investigate the effects of beryllium oxide and high assay low-enriched uranium fuels. this proposal considers using several existing resources at Sandia: (1) the Critical Experiments (SCX) facility and water tank, (2) spare UO₂-BeO fuel for the Annular Core Research Reactor (ACRR), and 7uPCX fuel rods from previous benchmark experiments.

NasGen provides a path for migration of structural models from Nastran bulk data format (BDF) into both an Exodus mesh file and an ASCII input file for Sierra Structural Dynamics (Salinas) and Solid Mechanics (Adagio). Many tools at Sandia National Labs (SNL) use the Exodus format. This document describes capabilities and limitations of the NasGen translation software.

The effective management of plastic waste streams to prevent plastic land and water pollution is a growing problem that is also one of the most important challenges in polymer science today. Polymer materials that are stable over their lifetime and can also be cheaply recycled or repurposed as desired could more easily be diverted from waste streams. However, this is difficult for most commodity plastics. It is especially difficult to conceive this with intractable, cross-linked polymers such as rubbers. In this work, we explore the utility of microencapsulated Grubbs’ catalysts for the in-situ depolymerization and reprocessing of polybutadiene (PB) rubber. Second-generation Hoveyda-Grubbs catalyst (HG2) contained within glassy thermoplastic microspheres can be dispersed in PB rubber below the microsphere’s glass transition temperature (Tg) without adverse depolymerization, evidenced by rubber with and without these microspheres obtaining similar shear storage moduli of ≈16 and ≈28 kPa, respectively. The thermoplastic’s Tg can be used to tune the depolymerization temperature, via release of HG2 into the rubber matrix. For example, using poly(lactic acid) (PLA) vs polysulfone results in an 85 and 162 °C depolymerization temperature, respectively. Liquefaction of rubber to a mixture of small molecules and oligomers is demonstrated using a 0.01 mol % catalyst loading using PLA as the encapsulant. At that same catalyst loading, depolymerization occurs to a greater extent in comparison to two ex-situ approaches, including a conventional solvent-assisted method, where it occurs at roughly twice the extent at each given catalyst loading. In addition, depolymerization of the microsphere-loaded rubbers was demonstrated for samples stored under nitrogen for 23 days. Lastly, we show that the depolymerized products can be reprocessed back into solid rubber with a shear storage modulus of ≈32 kPa. Thus, we envision that this approach could be used to recycle and reuse cross-linked rubbers at the end of their product lifetime.

This report provides documentation for the Sandia Toolkit (STK) modules. STK modules are intended to provide infrastructure that assists the development of computational engineering software such as finite-element analysis applications. STK includes modules for unstructured-mesh data structures, reading/writing mesh files, geometric proximity search, and various utilities. This document contains a chapter for each module, and each chapter contains overview descriptions and usage examples. Usage examples are primarily code listings which are generated from working test programs that are included in the STK code-base. A goal of this approach is to ensure that the usage examples will not fall out of date.

A digital twin has intelligent modules that continuously monitor the condition of the individual components and the whole of a system. Digital twins can provide nuclear power plants (NPP) operators an unprecedented level of monitoring, control, supervision, and security by contributing a greater volume of data for more comprehensive data analysis and increased accuracy of insights and predictions for decision making throughout the entire NPP lifecycle. NPP operators and managers have historically relied on limited, second hand or incomplete data. With proper implementation, digital twins can provide a central hub of all intel that allows for a multidisciplinary view of an NPP. This equips operators and managers with the ability to have more information, context, and intel that can be used for greater granularity during planning and decision making. Digital twins can be used in many activities as the technology has many different concepts surrounding it. From the various definitions of a digital twin within the industry, digital twins can be differentiated by levels of integration/automation. The three main models include digital model, digital shadow, and digital twin. Digital twins offer many potential advancements to the nuclear industry that could reduce costs, improve designs, provide safer operation, and improve their overall security.

Abstract not provided.

Freeplay is a common type of piecewise-smooth nonlinearity in dynamical systems, and it can cause discontinuity-induced bifurcations and other behaviors that may bring about undesirable and potentially damaging responses. Prior research has focused on piecewise-smooth systems with two or three distinct regions, but less attention is devoted to systems with more regions (i.e., multi-segmented systems). In this work, numerical analysis is performed on a dynamical system with multi-segmented freeplay, in which there are four stiffness transitions and five distinct regions in the phase space. The effects of the multi-segmented parameters are studied through bifurcation diagram evolution along with induced multi-stable behavior and different bifurcations. These phenomena are interrogated through various tools, such as harmonic balance, basins of attraction, phase planes, and Poincaré section analysis. Results show that among the three multi-segmented parameters, the asymmetry has the strongest effect on the response of the system.

High Energy Arcing Faults (HEAFs) are hazardous events in which an electrical arc leads to the rapid release of energy in the form of heat, vaporized metal, and mechanical force. In Nuclear Power Plants, these events are often accompanied by loss of essential power and complicated shutdowns. To confirm the probabilistic risk analysis (PRA) methodology in NUREG/CR-6850, which was formulated based on limited observational data, the NRC led an international experimental campaign from 2014 to 2016. The results of these experiments uncovered an unexpected hazard posed by aluminum components in or near electrical equipment and the potential for unanalyzed equipment failures. Sandia National Laboratories (SNL), in support of the NRC work, collaborated with NIST, BSI, KEMA, and NRC to support the full-scale HEAF test campaign in 2022. SNL provided high speed visible and infrared video/data of ten tests that collected data from HEAFs originated on copper and aluminum buses inside switchgears and bus ducts. Part of the SNL scope was to place cameras with high-speed data collection at different vantage points within the test facility to provide NRC a more complete and granular view of the test events.

Abstract not provided.

Abstract not provided.

In the dynamic landscape of Operational Technology (OT), and specifically the emerging landscape for Advanced Reactors, the establishment of trust between digital assets emerges as a challenge for cybersecurity modernization. This report reviews existing approaches to authentication in Enterprise environments, and proposed methods for authentication in OT, and analyzes each for its applicability to future Advanced Reactor digital networks. Principles of authentication ranging from underlying cryptographic mechanisms to trust authorities are evaluated through the lens of OT. These facets emphasize the importance of mutual authentication in real-time environments, enabling a paradigm shift from the current approach of strong boundaries to a more malleable network that allows for flexible operation. This work finds that there is a need for evaluation and decision making by industry stakeholders, but current technologies and approaches can be adapted to fit needs and risk tolerances.

Abstract not provided.

Hydrogen diffusion in metals and alloys plays an important role in the discovery of new materials for fuel cell and energy storage technology. While analytic models use hand-selected features that have clear physical ties to hydrogen diffusion, they often lack accuracy when making quantitative predictions. Machine learning models are capable of making accurate predictions, but their inner workings are obscured, rendering it unclear which physical features are truly important. To develop interpretable machine learning models to predict the activation energies of hydrogen diffusion in metals and random binary alloys, we create a database for physical and chemical properties of the species and use it to fit six machine learning models. Our models achieve root-mean-squared errors between 98-119 meV on the testing data and accurately predict that elemental Ru has a large activation energy, while elemental Cr and Fe have small activation energies. By analyzing the feature importances of these fitted models, we identify relevant physical properties for predicting hydrogen diffusivity. While metrics for measuring the individual feature importances for machine learning models exist, correlations between the features lead to disagreement between models and limit the conclusions that can be drawn. Instead grouped feature importance, formed by combining the features via their correlations, agree across the six models and reveal that the two groups containing the packing factor and electronic specific heat are particularly significant for predicting hydrogen diffusion in metals and random binary alloys. This framework allows us to interpret machine learning models and enables rapid screening of new materials with the desired rates of hydrogen diffusion.

Static structure factors are computed for large-scale, mechanically stable, jammed packings of frictionless spheres (three dimensions) and disks (two dimensions) with broad, power-law size dispersity characterized by the exponent -β. The static structure factor exhibits diverging power-law behavior for small wave numbers, allowing us to identify a structural fractal dimension df. In three dimensions, df≈2.0 for 2.5≤β≤3.8, such that each of the structure factors can be collapsed onto a universal curve. In two dimensions, we instead find 1.0df1.34 for 2.1≤β≤2.9. Furthermore, we show that the fractal behavior persists when rattler particles are removed, indicating that the long-wavelength structural properties of the packings are controlled by the large particle backbone conferring mechanical rigidity to the system. A numerical scheme for computing structure factors for triclinic unit cells is presented and employed to analyze the jammed packings.

Abstract not provided.

Tritium exhibits unique environmental behavior because of its potential interactions with water and organic substances. Modeling the environmental consequences of tritium releases can be relatively complex and thus an evaluation of MACCS is needed to understand what updates, if any, are needed in MACCS to account for the behavior of tritium. We examine documented tritium releases and previous benchmarking assessments to perform a model intercomparison between MACCS and state-of-practice tritium-specific codes UFOTRI and ETMOD to quantify the difference between MACCS and state of practice models for assessing tritium consequences. Additionally, information to assist an analyst in judging whether a postulated tritium release is likely to lead to significant doses is provided.

Numerical simulations are used to study the dynamics of a developing suspension Poiseuille flow with monodispersed and bidispersed neutrally buoyant particles in a planar channel, and machine learning is applied to learn the evolving stresses of the developing suspension. The particle stresses and pressure develop on a slower time scale than the volume fraction, indicating that once the particles reach a steady volume fraction profile, they rearrange to minimize the contact pressure on each particle. We consider the timescale for stress development and how the stress development connects to particle migration. For developing monodisperse suspensions, we present a new physics-informed Galerkin neural network that allows for learning the particle stresses when direct measurements are not possible. We show that when a training set of stress measurements is available, the MOR-physics operator learning method can also capture the particle stresses accurately.

Abstract not provided.

The methodology described in this article enables a type of holistic fleet optimization that simultaneously considers the composition and activity of a fleet through time as well as the design of individual systems within the fleet. Often, real-world system design optimization and fleet-level acquisition optimization are treated separately due to the prohibitive scale and complexity of each problem. This means that fleet-level schedules are typically limited to the inclusion of predefined system configurations and are blind to a rich spectrum of system design alternatives. Similarly, system design optimization often considers a system in isolation from the fleet and is blind to numerous, complex portfolio-level considerations. In reality, these two problems are highly interconnected. To properly address this system-fleet design interdependence, we present a general method for efficiently incorporating multi-objective system design trade-off information into a mixed-integer linear programming (MILP) fleet-level optimization. This work is motivated by the authors' experience with large-scale DOD acquisition portfolios. However, the methodology is general to any application where the fleet-level problem is a MILP and there exists at least one system having a design trade space in which two or more design objectives are parameters in the fleet-level MILP.

This report documents the preliminary design phase of the Critical Experiment Design (CED-1) conducted as part of integral experiment request (IER) 523. The purpose of IER-523 is to determine critical configurations of 35 weight percent (wt%) enriched uranium dioxideberyllium oxide (UO₂-BeO) material with Seven Percent Critical Experiment (7uPCX) fuels at Sandia National Laboratories (Sandia). Preliminary experiment design concepts, neutronic analysis results, and proposed paths for continuing the CED process are presented. This report builds on the feasibility and justification of experimental need report (CED-0) completed in December 2021.

A major difficulty in the analysis of molecular-level simulations is that macroscopic flow quantities are inherently noisy due to molecular fluctuations. An important example for turbulent flows is the kinetic energy dissipation rate. Traditionally, this quantity is calculated from gradients of the macroscopic velocity field, which exacerbates the noise problem. The inability to accurately compute the dissipation rate makes meaningful comparison of molecular-level and continuum simulation results a serious challenge. Herein, we extend previously developed coarse-graining theories to derive an exact molecular-level expression for the dissipation rate, which would circumvent the need to compute gradients of noisy fields. Although the exact expression cannot feasibly be implemented in Sandia’s direct simulation Monte Carlo (DSMC) code SPARTA, we utilize an approximate “hybrid” approach and compare it to the conventional gradient-based approach for planar Couette flow and the two-dimensional Taylor-Green vortex, demonstrating that the hybrid approach is significantly more accurate. Finally, we explore the possibility of adopting a Lagrangian approach to calculate the energy dissipation rate.

Abstract not provided.