Publications

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

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.

Abstract not provided.

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.

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.

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.

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.

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.

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.

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.

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.

Hydrogen Plus Other Alternative Fuels Risk Assessment Models (HyRAM+) is a software toolkit that provides a basis for quantitative risk assessment and consequence modeling for alternative fuels infrastructure and transportation systems. HyRAM+ integrates validated, analytical models of alternative fuel behavior, statistics, and a standardized quantitative risk assessment approach to generate useful, repeatable results for the safety analysis of various alternative fuel systems. This document demonstrates how to use HyRAM+ to analyze an example system, providing tutorials of HyRAM+ features with respect to system safety analysis and risk assessment.

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.

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.

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.

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.

Final report: PO 2222794

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.

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.

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.

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.

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%.

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.

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 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.

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.

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.

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.

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 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 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.

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.

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.

A design study was conducted at the National Solar Thermal Test Facility (NSTTF) at Sandia National Laboratories in Albuquerque, NM with the objective of identifying the technical readiness level, performance limits, capital and O&M costs, and expected thermal losses of particle handling and conveyance components in a particle-based CSP plant. Key findings indicated that skips and high temperature particle conveyance technology are available for moving particles up to 615° ± 25° C. This limits the use of mechanical conveyance above the heat exchanger and suggests vertical integration of the hot storage bin and heat exchanger to facilitate direct gravity fed handling of particles. Skip rails and support structures add significant cost and must be factored into cost analysis. Chutes can be a low-cost option for particle handling but uncertainties in tower costs make it difficult to know whether they can be cost effective in areas above the receiver.

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.

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.

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.

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.

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.

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.

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.

Abstract not provided.

The “Pioneer WEC” project is targeted at developing a wave energy generator for the Coastal Surface Mooring (CSM) system within the Ocean Observatories Initiative (OOI) Pioneer Array. The CSM utilizes solar photovoltaic and wind generation systems, along with rechargeable batteries, to power multiple sensors on the buoy and along the mooring line. This approach provides continuous power for essential controller functions and a subset of instruments, and meets the full power demand roughly 70% of the time. Sandia has been tasked with designing a wave energy system to provide additional electrical power and bring the CSM up-time for satisfying the full-power demand to 100%. This project is a collaboration between Sandia and Woods Hole Oceanographic Institution (WHOI), along with Evergreen Innovations, Monterey Bay Aquarium Research Institute (MBARI), Eastern Carolina University (ECU), Johns Hopkins University (JHU), and the National Renewable Energy Laboratory (NREL). This report captures Phase I of an expected two phase project and presents project scoping and concept design results. phase project and presents project scoping and concept design results.

Abstract not provided.

Large rocket motors may violently explode when exposed to accidental fires. Even hot metal fragments from a nearby accident may penetrate the propellant and ultimately cause thermal ignition. A mechanistic understanding of heated propellants leading to thermal runaway is a major unsolved problem. Here we show that thermal ignition in propellants can be predicted using a universal cookoff model coupled to a micromechanics pressurization model. Our model predicts the time to thermal ignition in cookoff experiments with variable headspace volumes. We found that experiments with headspace volumes are more prone to deformation which distorts pores and causes increased permeability when the propellant expands into this headspace. Delayed ignition with larger headspace volume correlates with lower headspace pressures during decomposition. We found that our predictions matched experimental measurements best when the initial propellant was impermeable to gas flow rather than being permeable. Similar behavior is expected with other energetic materials with rubbery binders. Our model is validated using data from a separate laboratory. We also present an uncertainty analysis using Latin Hypercube Sampling (LHS) of thermal ignition caused by a steel fragment embedded in the propellant.

Report for LDRD project 225992 with title Cryogenic Control Circuitry for Superconducting Qubits

Abstract not provided.

Abstract not provided.

The Sound Fixing and Ranging (SOFAR) channel in the ocean allows for low frequency sound to travel thousands of kilometers, making it particularly useful for detecting underwater nuclear explosions. Suggestions that an elevated SOFAR-like channel should exist in the stratosphere date back over half a century and imply that sources within this region can be reliably sensed at vast distances. However, this theory has not been supported with evidence of direct observations from sound within this channel. Here we show that an infrasound sensor on a solar hot air balloon recorded the first infrasound detection of a ground truth airborne source while within this acoustic channel, which we refer to as the AtmoSOFAR channel. Our results support the existence of the AtmoSOFAR channel, demonstrate that acoustic signals can be recorded within it, and provide insight into the characteristics of recorded signals. Results also show a lack of detections on ground-based stations, highlighting the advantages of using balloon-borne infrasound sensors to detect impulsive sources at altitude.

This is a SAND Report on cross-correlation data collected at the Redmond Salt Mine. It discusses methods, as well as temporal variability and energy characteristics of the cross-correlation data.

As large systems of Li-ion batteries are being increasingly deployed, the safety of such systems must be assessed. Due to the high cost of testing large systems, it is important to extract key safety information from any available experiments. Developing validated predictive models that can be exercised at larger scales offers an opportunity to augment experimental data In this work, experiments were conducted on packs of three Li-ion pouch cells with different heating rates and states of charge (SOC) to assess the propagation behavior of a module undergoing thermal runaway. The variable heating rates represent slow or fast heating that a module may experience in a system. As the SOC decreases, propagation slows down and eventually becomes mitigated. It was found that the SOC boundary between propagation and mitigation was higher at a heating rate of 50 °C/min than at 10 °C/min for these cells. However, due to increased pre-heating at the lower heating rate, the propagation speed increased. Simulations were conducted with a new intra-particle diffusion-limited reaction model for a range of anode particle sizes. Propagation speeds and onset times were generally well predicted, and the variability in the propagation/mitigation boundary highlighted the need for greater uncertainty quantification of the predictions.

In this paper, an approach for 3D plasma structure diagnostics using tomographic optical emission spectroscopy (Tomo-OES) of a nanosecond pulsed atmospheric pressure plasma jet (APPJ) is presented. In contrast to the well-known Abel inversion, Tomo-OES does not require cylindrical symmetry to recover 3D distributions of plasma light emission. Instead, many 2D angular projections are measured with intensified cameras and the multiplicative algebraic reconstruction technique is used to recover the 3D distribution of light emission. This approach solves the line-of-sight integration problem inherent to optical diagnostics, allowing recovery of localized OES information within the plasma that can be used to better infer plasma parameters within complex plasma structures. Here, Tomo-OES was applied to investigate an APPJ operated with helium in ambient air and impinging on planar and structured dielectric surfaces. Surface charging caused the guided streamer from the APPJ to transition to a surface ionization wave (SIW) that propagated along the surface. The SIW experienced variable geometrical and electrical material properties as it propagated, leading to 3D configurations that were non-symmetric and spatially complex. Light emission from He, N 2 + , and N2 were imaged at ten angular projections and the respective time-resolved 3D emission distributions in the plasma were then reconstructed. The spatial resolution of each tomographic reconstruction was 7.4 µm and the temporal resolution was 5 ns, sufficient to observe the guided streamer and the effects of the structured surface on the SIW. Emission from He showed the core of the jet and emission from N 2 + and N2 indicated effects of entrainment of ambient air. Penning ionization of N2 created a ring or outer layer of N 2 + that spatially converged to form the ‘plasma bullet’ or spatially diverged across a surface as part of a SIW. The SIW entered trenches of size 150 µm, leading to decreases in plasma light emission in regions above the trenches. The plasma light emission was higher in some regions with trenches, possibly due to effects of field enhancement.

Abstract not provided.

This is a milestone report due in 2 weeks

We describe a direct magneto-optical approach to measuring the magnetic field driven by a narrow pulse width (<10 ns), 20 kA electrical current flow in the transmission line of a high energy pulsed power accelerator. The magnetic field and electrical current are among the most important operating parameters in a pulsed power accelerator and are critical to understanding the properties of the radiation output. However, accurately measuring these fields and electrical currents using conventional pulsed power diagnostics is difficult due to the strength of ionizing radiation and electromagnetic interference. Our approach uses a fiber coupled laser beam with a rare earth element sensing crystal sensor that is highly resistant to electromagnetic interference and does not require external calibration. Here, we focus on device theory, operating parameters, results from an experiment on a high energy pulsed power accelerator, and comparison to a conventional electrical current shunt sensor.

Machine learning models are promising as surrogates in optimization when replacing difficult to solve equations or black-box type models. This work demonstrates the viability of linear model decision trees as piecewise-linear surrogates in decision-making problems. Linear model decision trees can be represented exactly in mixed-integer linear programming (MILP) and mixed-integer quadratic constrained programming (MIQCP) formulations. Furthermore, they can represent discontinuous functions, bringing advantages over neural networks in some cases. We present several formulations using transformations from Generalized Disjunctive Programming (GDP) formulations and modifications of MILP formulations for gradient boosted decision trees (GBDT). We then compare the computational performance of these different MILP and MIQCP representations in an optimization problem and illustrate their use on engineering applications. We observe faster solution times for optimization problems with linear model decision tree surrogates when compared with GBDT surrogates using the Optimization and Machine Learning Toolkit (OMLT).

Abstract not provided.

Accurate targeting of radioisotope classifiers and estimators requires an understanding of the target problem space. In order to facilitate clear communication on expected model behavior and performance between practitioners and stakeholders on their problems, this questionnaire was created. Stakeholder responses form the basis of a trained model as well as the start of usage requirements for the model as it is integrated with analysis processes or detection systems. This questionnaire may also be useful to machine learning practitioners and gamma spectroscopists developing new algorithms as a starting point for characterizing their problem space, especially if they are using PyRIID.

Advancements in space exploration and sample return technology present a unique opportunity to leverage sample return capsules (SRCs) towards studying atmospheric entry of meteoroids and asteroids. Specifically engineered for the secure transport of valuable extraterrestrial samples from interplanetary space to Earth, SRCs offer unexpected benefits that reach beyond their intended purpose. As SRCs enter the Earth’s atmosphere at hypervelocity, they are analogous to naturally occurring meteoroids and thus, for all intents and purposes, can be considered artificial meteors. Furthermore, SRCs are capable of generating shockwaves upon reaching the lower transitional flow regime, and thus can be detected by strategically positioned geophysical instrumentation. NASA’s OSIRIS-REx (Origins, Spectral Interpretation, Resource Identification, and Security-Regolith Explorer) SRC is one of only a handful of artificial objects to re-enter the Earth’s atmosphere from interplanetary space since the end of the Apollo era and it will provide an unprecedented observational opportunity. This review summarizes past infrasound and seismic observational studies of SRC re-entries since the end of the Apollo era and presents their utility towards the better characterization of meteoroid flight through the atmosphere.

Nonlinear modeling and optimization is a valuable tool for aiding decisions by engineering practitioners, but programming an optimization problem based on a complex electrical, mechanical, or chemical process is a time-consuming and error-prone activity. Therefore, there is a need for model analysis and debugging tools that can detect and diagnose modeling errors. One such tool is the Dulmage–Mendelsohn decomposition, which identifies structurally under- and over-determined subsets in systems of equations and variables by partitioning the bipartite graph of the system. This work provides the necessary background to understand the Dulmage–Mendelsohn decomposition and its application to the analysis of nonlinear optimization problems, demonstrates its use in diagnosing a variety of modeling errors, and introduces software implementations for analyzing nonlinear optimization problems in the Pyomo and JuMP algebraic modeling languages.

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.

Abstract not provided.

Abstract not provided.

Improved power take-off (PTO) controller design for wave energy converters is considered a critical component for reducing the cost of energy production. However, the device and control design process often remains sequential, with the space of possible final designs largely reduced before the controller has been considered. Control co-design, whereby the device and control design are considered concurrently, has resulted in improved designs in many industries, but remains rare in the wave energy community. In this paper we demonstrate the use of a new open-source code, WecOptTool, for control co-design of wave energy converters, with the aim to make the co-design approach more accessible and accelerate its adoption. Additionally, we highlight the importance of designing a wave energy converter to maximize electrical power, rather than mechanical power, and demonstrate the co-design process while modeling the PTO's components (i.e., drive-train and generator, and their dynamics). We also consider the design and optimization of causal fixed-structure controllers. The demonstration presented here considers the PTO design problem and finds the optimal PTO drive-train that maximizes annual electrical power production. The results show a 22% improvement in the optimal controller and drive-train co-design over the optimal controller for the nominal, as built, device design.

Multifidelity uncertainty quantification (MF UQ) sampling approaches have been shown to significantly reduce the variance of statistical estimators while preserving the bias of the highest-fidelity model, provided that the low-fidelity models are well correlated. However, maintaining a high level of correlation can be challenging, especially when models depend on different input uncertain parameters, which drastically reduces the correlation. Existing MF UQ approaches do not adequately address this issue. In this work, we propose a new sampling strategy that exploits a shared space to improve the correlation among models with dissimilar parameterization. We achieve this by transforming the original coordinates onto an auxiliary manifold using the adaptive basis (AB) method (Tipireddy and Ghanem, 2014). The AB method has two main benefits: (1) it provides an effective tool to identify the low-dimensional manifold on which each model can be represented, and (2) it enables easy transformation of polynomial chaos representations from high- to low-dimensional spaces. This latter feature is used to identify a shared manifold among models without requiring additional evaluations. We present two algorithmic flavors of the new estimator to cover different analysis scenarios, including those with legacy and non-legacy high-fidelity (HF) data. We provide numerical results for analytical examples, a direct field acoustic test, and a finite element model of a nuclear fuel assembly. For all examples, we compare the proposed strategy against both single-fidelity and MF estimators based on the original model parameterization.

Abstract not provided.

Abstract not provided.

This report provides an overview of some code written to streamline aspects of the finite element simulation process using Sierra/SD. Some specific examples of common simulation tasks are 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.

Approximating differential operators defined on two-dimensional surfaces is an important problem that arises in many areas of science and engineering. Over the past ten years, localized meshfree methods based on generalized moving least squares (GMLS) and radial basis function finite differences (RBF-FD) have been shown to be effective for this task as they can give high orders of accuracy at low computational cost, and they can be applied to surfaces defined only by point clouds. However, there have yet to be any studies that perform a direct comparison of these methods for approximating surface differential operators (SDOs). The first purpose of this work is to fill that gap. For this comparison, we focus on an RBF-FD method based on polyharmonic spline kernels and polynomials (PHS+Poly) since they are most closely related to the GMLS method. Additionally, we use a relatively new technique for approximating SDOs with RBF-FD called the tangent plane method since it is simpler than previous techniques and natural to use with PHS+Poly RBF-FD. The second purpose of this work is to relate the tangent plane formulation of SDOs to the local coordinate formulation used in GMLS and to show that they are equivalent when the tangent space to the surface is known exactly. The final purpose is to use ideas from the GMLS SDO formulation to derive a new RBF-FD method for approximating the tangent space for a point cloud surface when it is unknown. For the numerical comparisons of the methods, we examine their convergence rates for approximating the surface gradient, divergence, and Laplacian as the point clouds are refined for various parameter choices. We also compare their efficiency in terms of accuracy per computational cost, both when including and excluding setup costs.

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.

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.

A summary of causal loop diagram conceptual modeling performed under the Water Intersections with Climate Systems Security Strategic Initiative.

Abstract not provided.

Abstract not provided.

As Machine Learning (ML) continues to advance, it is being integrated into more systems. Often, the ML component represents a significant portion of the system that reduces the burden on the end user or significantly improves task performance. However, the ML component represents an unknown complex phenomenon that is learned from collected data without the need to be explicitly programmed. Despite the improvement in task performance, the models are often black boxes. Evaluating the credibility and the vulnerabilities of ML models poses a gap in current test and evaluation practice. For high consequence applications, the lack of testing and evaluation procedures represents a significant source of uncertainty and risk. To help reduce that risk, here we present considerations to evaluate systems embedded with an ML component within a red-teaming inspired methodology. We focus on (1) cyber vulnerabilities to an ML model, (2) evaluating performance gaps, and (3) adversarial ML vulnerabilities.

The benefits of high-performance unidirectional carbon fiber composites are limited in many cost-driven industries due to the high cost relative to alternative reinforcement fibers. Low-cost carbon fibers have been previously proposed, but the longitudinal compressive strength continues to be a limiting factor or studies are based on simplifications that warrant further analysis. A micromechanical model is used to (1) determine if the longitudinal compressive strength of composites can be improved with noncircular carbon fiber shapes and (2) characterize why some shapes are stronger than others in compression. In comparison to circular fibers, the results suggest that the strength can be increased by 10%–13% by using a specific six-lobe fiber shape and by 6%–9% for a three-lobe fiber shape. A slight increase is predicted in the compressive strength of the study two-lobe fiber but has the highest uncertainty and sensitivity to fiber orientation and misalignment direction. The underlying mechanism governing the compressive failure of the composites was linked to the unique stress fields created by the lobes, particularly the pressure stress in the matrix. This work provides mechanics-based evidence of strength improvements from noncircular fiber shapes and insight on how matrix yielding is altered with alternative fiber shapes.

The development of additively-manufactured (AM) 316L stainless steel (SS) using laser powder bed fusion (LPBF) has enabled near net shape components from a corrosion-resistant structural material. In this article, we present a multiscale study on the effects of processing parameters on the corrosion behavior of as-printed surfaces of AM 316L SS formed via LPBF. Laser power and scan speed of the LPBF process were varied across the instrument range known to produce parts with >99 % density, and the macroscale corrosion trends were interpreted via microscale and nanoscale measurements of porosity, roughness, microstructure, and chemistry. Porosity and roughness data showed that porosity φ decreased as volumetric energy density Ev increased due to a shift in the pore formation mechanism and that roughness Sq was due to melt track morphology and partially fused powder features. Cross-sectional and plan-view maps of chemistry and work function ϕs revealed an amorphous Mn-silicate phase enriched with Cr and Al that varied in both thickness and density depending on Ev. Finally, the macroscale potentiodynamic polarization experiments under full immersion in quiescent 0.6 M NaCl showed significant differences in breakdown potential Eb and metastable pitting. In general, samples with smaller φ and Sq values and larger ϕs values and homogeneity in the Mn-silicate exhibited larger Eb. The porosity and roughness effects stemmed from an increase to the overall number of initiation sites for pitting, and the oxide phase contributed to passive film breakdown by acting as a crevice former or creating a galvanic couple with the SS.

Characterizing interface trap states in commercial wide bandgap devices using frequency-based measurements requires unconventionally high probing frequencies to account for both fast and slow traps associated with wide bandgap materials. The C − ψ S technique has been suggested as a viable quasi-static method for determining the interface trap state densities in wide bandgap systems, but the results are shown to be susceptible to errors in the analysis procedure. This work explores the primary sources of errors present in the C − ψ S technique using an analytical model that describes the apparent response for wide bandgap MOS capacitor devices. Measurement noise is shown to greatly impact the linear fitting routine of the 1 / C S ∗ 2 vs ψ S plot to calibrate the additive constant in the surface potential/gate voltage relationship, and an inexact knowledge of the oxide capacitance is also shown to impede interface trap state analysis near the band edge. In addition, a slight nonlinearity that is typically present throughout the 1 / C S ∗ 2 vs ψ S plot hinders the accurate estimation of interface trap densities, which is demonstrated for a fabricated n-SiC MOS capacitor device. Methods are suggested to improve quasi-static analysis, including a novel method to determine an approximate integration constant without relying on a linear fitting routine.

Systems engineering today faces a wide array of challenges, ranging from new operational environments to disruptive technological — necessitating approaches to improve research and development (R&D) efforts. Yet, emphasizing the Aristotelian argument that the “whole is greater than the sum of its parts” seems to offer a conceptual foundation creating new R&D solutions. Invoking systems theoretic concepts of emergence and hierarchy and analytic characteristics of traceability, rigor, and comprehensiveness is potentially beneficial for guiding R&D strategy and development to bridge the gap between theoretical problem spaces and engineering-based solutions. In response, this article describes systems–theoretic process analysis (STPA) as an example of one such approach to aid in early-systems R&D discussions. STPA—a ‘top-down’ process that abstracts real complex system operations into hierarchical control structures, functional control loops, and control actions—uses control loop logic to analyze how control actions (designed for desired system behaviors) may become violated and drive the complex system toward states of higher risk. By analyzing how needed controls are not provided (or out of sequence or stopped too soon) and unneeded controls are provided (or engaged too long), STPA can help early-system R&D discussions by exploring how requirements and desired actions interact to either mitigate or potentially increase states of risk that can lead to unacceptable losses. This article will demonstrate STPA's benefit for early-system R&D strategy and development discussion by describing such diverse use cases as cyber security, nuclear fuel transportation, and US electric grid performance. Together, the traceability, rigor, and comprehensiveness of STPA serve as useful tools for improving R&D strategy and development discussions. In conclusion, leveraging STPA as well as related systems engineering techniques can be helpful in early R&D planning and strategy development to better triangulate deeper theoretical meaning or evaluate empirical results to better inform systems engineering solutions.

Magnetic properties of more than 20 Cantor alloy samples of varying composition were investigated over a temperature range of 5 K to 300 K and in fields of up to 70 kOe using magnetometry and muon spin relaxation. Two transitions are identified: a spin-glass-like transition that appears between 55K and 190K, depending on composition, and a ferrimagnetic transition that occurs at approximately 43K in multiple samples with widely varying compositions. The magnetic signatures at 43K are remarkably insensitive to chemical composition. A modified Curie-Weiss model was used to fit the susceptibility data and to extract the net effective magnetic moment for each sample. The resulting values for the net effective moment were either diminished with increasing Cr or Mn concentrations or enhanced with decreasing Fe, Co, or Ni concentrations. Beyond a sufficiently large effective moment, the magnetic ground state transitions from ferrimagnetism to ferromagnetism. The effective magnetic moments, together with the corresponding compositions, are used in a global linear regression analysis to extract element-specific effective magnetic moments, which are compared to the values obtained by ab initio based density functional theory calculations. Finally, these moments provide the information necessary to controllably tune the magnetic properties of Cantor alloy variants.

A new capability for modeling graded density reactive flow materials in the shock physics hydrocode, CTH, is demonstrated here. Previously, materials could be inserted in CTH with graded material properties, but the sensitivity of the material was not adjusted based on these properties. Of particular interest are materials that are graded in density, sometimes due to pressing or other assembly operations. The sensitivity of explosives to both density and temperature has been well demonstrated in the literature, but to-date the material parameters for use in a simulation were fit to a single condition and applied to the entire material, or the material had to be inserted in sections and each section assigned a condition. The reactive flow model xHVRB has been extended to shift explosive sensitivity with initial density, so that sensitivity is also graded in the material. This capability is demonstrated for use in three examples. The first models detonation transfer in a graded density pellet of HNS, the second is a shaped charge with density gradients in the explosive, and the third is an explosively formed projectile.

Maximizing the production of heterologous biomolecules is a complex problem that can be addressed with a systems-level understanding of cellular metabolism and regulation. Specifically, growth-coupling approaches can increase product titers and yields and also enhance production rates. However, implementing these methods for non-canonical carbon streams is challenging due to gaps in metabolic models. Over four design-build-test-learn cycles, we rewire Pseudomonas putida KT2440 for growth-coupled production of indigoidine from para-coumarate. We explore 4,114 potential growth-coupling solutions and refine one design through laboratory evolution and ensemble data-driven methods. The final growth-coupled strain produces 7.3 g/L indigoidine at 77% maximum theoretical yield in para-coumarate minimal medium. The iterative use of growth-coupling designs and functional genomics with experimental validation was highly effective and agnostic to specific hosts, carbon streams, and final products and thus generalizable across many systems.

For reactive burn models in hydrocodes, an equilibrium closure assumption is typically made between the unreacted and product equations of state. In the CTH [1] (not an acronym) hydrocode the assumption of density and temperature equilibrium is made by default, while other codes make a pressure and temperature equilibrium assumption. The main reason for this difference is the computational efficiency in making the density and temperature assumption over the pressure and temperature one. With fitting to data, both assumptions can accurately predict reactive flow response using the various models, but the model parameters from one code cannot necessarily be used directly in a different code with a different closure assumption. A new framework is intro-duced in CTH to allow this assumption to be changed independently for each reactive material. Comparisons of the response and computational cost of the History Variable Reactive Burn (HVRB) reactive flow model with the different equilibrium assumptions are presented.

Process variations within Field Programmable Gate Arrays (FPGAs) provide a rich source of entropy and are therefore well-suited for the implementation of Physical Unclonable Functions (PUFs). However, careful considerations must be given to the design of the PUF architecture as a means of avoiding undesirable localized bias effects that adversely impact randomness, an important statistical quality characteristic of a PUF. Here in this paper, we investigate a ring-oscillator (RO) PUF that leverages localized entropy from individual look-up table (LUT) primitives. A novel RO construction is presented that enables the individual paths through the LUT primitive to be measured and isolated at high precision, and an analysis is presented that demonstrates significant levels of localized design bias. The analysis demonstrates that delay-based PUFs that utilize LUTs as a source of entropy should avoid using FPGA primitives that are localized to specific regions of the FPGA, and instead, a more robust PUF architecture can be constructed by distributing path delay components over a wider region of the FPGA fabric. Compact RO PUF architectures that utilize multiple configurations within a small group of LUTs are particularly susceptible to these types of design-level bias effects. The analysis is carried out on data collected from a set of identically designed, hard macro instantiations of the RO implemented on 30 copies of a Zynq 7010 SoC.

Organic co-crystals have emerged as a promising class of semiconductors for next-generation optoelectronic devices due to their unique photophysical properties. This paper presents a joint experimental-theoretical study comparing the crystal structure, spectroscopy, and electronic structure of two charge transfer co-crystals. Reported herein is a novel co-crystal Npe:TCNQ, formed from 4-(1-naphthylvinyl)pyridine (Npe) and 7,7,8,8-tetracyanoquinodimethane (TCNQ) via molecular self-assembly. This work also presents a revised study of the co-crystal composed of Npe and 1,2,4,5-tetracyanobenzene (TCNB) molecules, Npe:TCNB, herein reported with a higher-symmetry (monoclinic) crystal structure than previously published. Npe:TCNB and Npe:TCNQ dimer clusters are used as theoretical model systems for the co-crystals; the geometries of the dimers are compared to geometries of the extended solids, which are computed with periodic boundary conditions density functional theory. UV-Vis absorption spectra of the dimers are computed with time-dependent density functional theory and compared to experimental UV-Vis diffuse reflectance spectra. Both Npe:TCNB and Npe:TCNQ are found to exhibit neutral character in the S0 state and ionic character in the S1 state. The high degree of charge transfer in the S1 state of both Npe:TCNB and Npe:TCNQ is rationalized by analyzing the changes in orbital localization associated with the S1 transitions.

Abstract not provided.

For the first time the optimal local truncation error method (OLTEM) with 125-point stencils and unfitted Cartesian meshes has been developed in the general 3-D case for the Poisson equation for heterogeneous materials with smooth irregular interfaces. The 125-point stencils equations that are similar to those for quadratic finite elements are used for OLTEM. The interface conditions for OLTEM are imposed as constraints at a small number of interface points and do not require the introduction of additional unknowns, i.e., the sparse structure of global discrete equations of OLTEM is the same for homogeneous and heterogeneous materials. The stencils coefficients of OLTEM are calculated by the minimization of the local truncation error of the stencil equations. These derivations include the use of the Poisson equation for the relationship between the different spatial derivatives. Such a procedure provides the maximum possible accuracy of the discrete equations of OLTEM. In contrast to known numerical techniques with quadratic elements and third order of accuracy on conforming and unfitted meshes, OLTEM with the 125-point stencils provides 11-th order of accuracy, i.e., an extremely large increase in accuracy by 8 orders for similar stencils. The numerical results show that OLTEM yields much more accurate results than high-order finite elements with much wider stencils. The increased numerical accuracy of OLTEM leads to an extremely large increase in computational efficiency. Additionally, a new post-processing procedure with the 125-point stencil has been developed for the calculation of the spatial derivatives of the primary function. The post-processing procedure includes the minimization of the local truncation error and the use of the Poisson equation. It is demonstrated that the use of the partial differential equation (PDE) for the 125-point stencils improves the accuracy of the spatial derivatives by 6 orders compared to post-processing without the use of PDE as in existing numerical techniques. At an accuracy of 0.1% for the spatial derivatives, OLTEM reduces the number of degrees of freedom by 900 - 4∙10⁶ times compared to quadratic finite elements. The developed post-processing procedure can be easily extended to unstructured meshes and can be independently used with existing post-processing techniques (e.g., with finite elements).

The purpose of this report is to document updates on the apparatus to simulate commercial vacuum drying procedures at the Nuclear Energy Work Complex at Sandia National Laboratories. Validation of the extent of water removal in a dry spent nuclear fuel storage system based on drying procedures used at nuclear power plants is needed to close existing technical gaps. Operational conditions leading to incomplete drying may have potential impacts on the fuel, cladding, and other components in the system during subsequent storage and disposal. A general lack of data suitable for model validation of commercial nuclear canister drying processes necessitates well-designed investigations of drying process efficacy and water retention. Scaled tests that incorporate relevant physics and well-controlled boundary conditions are essential to provide insight and guidance to the simulation of prototypic systems undergoing drying processes. This report documents fiscal year 2023 (FY23) updates on the Advanced Drying Cycle Simulator (ADCS). This apparatus was built to simulate commercial drying procedures and quantify the amount of residual water remaining in a pressurized water reactor (PWR) fuel assembly after drying. The ADCS was constructed with a prototypic 17×17 PWR fuel skeleton and waterproof heater rods to simulate decay heat. These waterproof heaters are the next generation design to heater rods developed and tested at Sandia National Laboratories in FY20. In FY23, a series of four simulated commercial drying tests was completed. This report presents the temperature and pressure histories of the drying tests as well as axial temperature profiles that can be compared to data from the Electric Power Research Institute (EPRI) High Burnup Demonstration TN-32B cask. Water content measurements and dew point calculations from a Hiden Analytical HPR-30 mass spectrometer are also presented in this report. Due to familiarization with this first-of-a-kind system, refinements to equipment calibration and test procedures have been identified to better match commercial drying cycles for future simulated tests. However, the presented data demonstrate the successful construction and operation of a viable research platform for quantifying residual water content closely approaching that expected in dry storage canisters during commercial drying procedures.

Work evaluating spent nuclear fuel (SNF) dry storage canister surface environments and canister corrosion progressed significantly in FY23, with the goal of developing a scientific understanding of the processes controlling initiation and growth of stress corrosion cracking (SCC) cracks in stainless steel canisters in relevant storage environments. The results of the work performed at Sandia National Laboratories (SNL) will guide future work and will contribute to the development of better tools for predicting potential canister penetration by SCC.

In this study, we demonstrate the ability of polarity inversion of sputtered aluminum scandium nitride thin films through post-fabrication processes with domain widths as small as 220 nm at a periodicity of 440 nm. An approach using photo- and electron-beam lithography to generate sub-quarter micrometer feature size with adjustable duty cycle through a lift-off process is presented. The film with a coercive field Ec+ of 5.35 MV/cm was exercised first with a 1 kHz triangular double bipolar wave and ultimately poled with a 0.5 kHz double monopolar wave using a Radiant Precision Premier II tester. The metal polar (M-polar) and nitrogen polar (N-polar) domains were identified and characterized through potassium hydroxide wet etching as well as piezoresponse force microscopy (PFM). Well-distinguished boundaries between the oppositely polarized domain regions were confirmed through the phase diagram of the PFM results. The relationship between the electrode width, poling voltage, and domain growth was experimentally studied and statistically analyzed, where 7.96 nm/V domain width broadening vs escalating poling voltage was observed. This method produces extremely high domain spatial resolution in III-nitride materials via poling and is transferable to a CMOS-compatible photolithography process. The spatial resolution of the periodically poled Al0.68Sc0.32N is suitable for second-harmonic generation of deep ultraviolet through quasi-phase-matching and RF MEMS operating in the X-Band spectrum.

Accurate fuel oxidation mechanisms can enable predictive capabilities that aid in advancing combustion technologies. High-level computational kinetics can yield reasonable rate coefficients with uncertainties, in some cases, below a factor of 2. Computed rate coefficients can be constrained further by optimizing against experimental data. Here, we explore the application of genetic algorithm (GA) optimization to constrain computed rate coefficients in complex fuel oxidation mechanisms in conjunction with temperature-dependent species mole fractions from jet-stirred reactor (JSR) measurements. Cyclohexane is a model candidate for understanding the reactivity of cyclic fuels. In this work, we optimize the rate coefficients of the most recent literature cyclohexane mechanism, which incorporates theoretically computed rate coefficients for the reaction networks stemming from the first and second O2 addition pathways, against the experimental results of two separate literature JSR studies. Optimization consistency is evaluated by carrying out three GA optimizations: fitting to the temperature-dependent species mole fractions in each JSR experiment separately and simultaneously fitting the species mole fractions in both experiments. Local sensitivity analyses are used to identify five influential low-temperature oxidation reactions for optimization. Although the three optimizations do not yield identical rate coefficients, the direction of change in all five rate coefficients is consistent among the three optimizations. Performance of the models from the three optimizations is assessed against literature ignition delay times with differences in the level of agreement observed among the different optimizations. Comparisons are made with our recent optimization work of a cyclopentane oxidation master-equation model against time-resolved species concentrations, and insights and improvements of the strategy for constraining rate coefficients using GA optimization are discussed.

Tritium permeability in zirconium-based tritium getter critically impacts tritium storage and environmental safety during operation of tritium-producing burnable absorber rods (TPBARs). Previous experiments indicated that during irradiation operation, the hydrogen equilibrium pressured is increased. Further experimental and modeling studies suggested that the enhanced tritium release observed for reactor scale assemblies might be related to a thermal diffusion known as the Soret effect. A direct measurement of the Soret factor, however, has not been performed. To improve TPBAR and other nuclear applications, here we have applied two non-equilibrium molecular dynamics methods to study thermal diffusion of hydrogen isotopes in low-concentration zirconium hydrides. One of the methods produces sufficiently converged results to distinguish crystal orientation, isotope type, and concentration effects. In conclusion, with this method, crystal orientation, isotope type, and concentration effects are discussed.

Local crystallographic features negatively affect quantum spin defects by changing the local electrostatic environment, often resulting in degraded or varied qubit optical and coherence properties. Few tools exist that enable the deterministic synthesis and study of such intricate systems on the nano-scale, making defect-to-defect strain environment quantification difficult. In this paper, we highlight state-of-the-art capabilities from the U.S. Department of Energy’s Nanoscale Science Research Centers that directly address these shortcomings. Specifically, we demonstrate how complementary capabilities of nano-implantation and nano-diffraction can be used to demonstrate the quantum relevant, spatially deterministic creation of neutral divacancy centers in 4H silicon carbide, while investigating and characterizing these systems on the ≤ 25 nm scale with strain sensitivities on the order of 1 × 10 − 6 , relevant to defect formation dynamics. This work lays the foundation for ongoing studies into the dynamics and deterministic formation of low strain homogeneous quantum relevant spin defects in the solid state.

Microneedle sensors could enable minimally-invasive, continuous molecular monitoring – informing on disease status and treatment in real-time. Wearable sensors for pharmaceuticals, for example, would create opportunities for treatments personalized to individual pharmacokinetics. Here, we demonstrate a commercial-off-the-shelf (COTS) approach for microneedle sensing using an electrochemical aptamer-based sensor that detects the high-toxicity antibiotic, vancomycin. Wearable monitoring of vancomycin could improve patient care by allowing targeted drug dosing within its narrow clinical window of safety and efficacy. To produce sensors, we miniaturize the electrochemical aptamer-based sensors to a microelectrode format, and embed them within stainless steel microneedles (sourced from commercial insulin pen needles). The microneedle sensors achieve quantitative measurements in body-temperature undiluted blood. Further, the sensors effectively maintain electrochemical signal within porcine skin. This COTS approach requires no cleanroom fabrication or specialized equipment, and produces individually-addressable, sterilizable microneedle sensors capable of easily penetrating the skin. In the future, this approach could be adapted for multiplexed detection, enabling real-time monitoring of a range of biomarkers.

The Sibling Pin test campaign is a Department of Energy (DOE) research activity within the Spent Fuel and Waste Science and Technology (SFWST) program that is tasked with characterization of high burnup (HBU) fuel in support of the High Burnup Spent Fuel Data Project. Of the 25 fuel rods in the Sibling Pin inventory, approximately 9 rod lengths have been consumed during the first phase (Phase I) of the test campaign leaving approximately 16 rod lengths for the second phase (Phase II) of testing. This plan outlines the Phase II testing and the motivations for performing these tests. Priorities for Phase II testing are based on previously identified knowledge gaps, lessons-learned from Phase I work, the original objectives of the High Burnup Spent Fuel Data Project and the Sibling Pin test campaign, and input from external stakeholders. The priorities for Phase II testing are to obtain data to characterize the effects of annealing on cladding mechanical properties and fuel rod performance, to quantify the creep behavior of cladding materials and fuel rods and the effects of creep deformations on the performance of cladding and fuel rods, and to gather data to support the final closure of the hydride reorientation and radial hydride induced embrittlement gap for HBU fuel rods.

Translating the surging interest in neuromorphic electronic components, such as those based on nonlinearities near Mott transitions, into large-scale commercial deployment faces steep challenges in the current lack of means to identify and design key material parameters. These issues are exemplified by the difficulties in connecting measurable material properties to device behavior via circuit element models. Here, the principle of local activity is used to build a model of VO2/SiN Mott threshold switches by sequentially accounting for constraints from a minimal set of quasistatic and dynamic electrical and high-spatial-resolution thermal data obtained via in situ thermoreflectance mapping. By combining independent data sets for devices with varying dimensions, the model is distilled to measurable material properties, and device scaling laws are established. The model can accurately predict electrical and thermal conductivities and capacitances and locally active dynamics (especially persistent spiking self-oscillations). The systematic procedure by which this model is developed has been a missing link in predictively connecting neuromorphic device behavior with their underlying material properties, and should enable rapid screening of material candidates before employing expensive manufacturing processes and testing procedures.

Motivated by recent experimental results we calculate from first-principles the lifetime of low-energy quasiparticles in bilayer graphene (BLG). Here, we take into account the scattering rate arising from electron-electron interactions within the GW approximation for the electron self-energy and consider several p-type doping levels ranging from 0 to ρ ≈ 2.4 × 10¹² holes/cm². In the undoped case we find that the average inverse lifetime scales linearly with energy away from the charge neutrality point, with values in good agreement with experiments. The decay rate is approximately three times larger than in monolayer graphene, a consequence of the enhanced screening in BLG. In the doped case, the dependence of the inverse lifetime on quasiparticle energy acquires a non-linear component due to the opening of an additional decay channel mediated by acoustic plasmons.

Here, we study the impact of deep-level defects on trap-assisted Auger–Meitner recombination in c-plane InGaN/GaN LEDs using a small-signal electroluminescence (SSEL) method and deep-level optical spectroscopy (DLOS). Carrier dynamics information, including carrier lifetime, recombination rate, and carrier density, is obtained from SSEL, while DLOS is used to obtain the deep-level defect density. Through fitting the nonradiative recombination rates of wafers with different deep-level defect densities, we obtain the Shockley–Read–Hall (SRH) and trap-assisted Auger–Meitner recombination (TAAR) coefficients. We show that defect-related nonradiative recombination, including both SRH and TAAR, accounts for a relatively small fraction of the total nonradiative recombination, which is dominated by intrinsic Auger–Meitner recombination. The interplay between carrier localization and Coulomb enhancement has a different impact on radiative and intrinsic Auger–Meitner recombination. Evidence is presented that the imbalance between the change of radiative and intrinsic Auger–Meitner recombination is the primary cause of the efficiency droop at high carrier densities in the samples studied.

We study the impact of deep-level defects on trap-assisted Auger-Meitner recombination in c-plane InGaN/GaN LEDs using a small-signal electroluminescence (SSEL) method and deep-level optical spectroscopy (DLOS). Carrier dynamics information, including carrier lifetime, recombination rate, and carrier density, is obtained from SSEL, while DLOS is used to obtain the deep-level defect density. Through fitting the nonradiative recombination rates of wafers with different deep-level defect densities, we obtain the Shockley-Read-Hall (SRH) and trap-assisted Auger-Meitner recombination (TAAR) coefficients. We show that defect-related nonradiative recombination, including both SRH and TAAR, accounts for a relatively small fraction of the total nonradiative recombination, which is dominated by intrinsic Auger-Meitner recombination. The interplay between carrier localization and Coulomb enhancement has a different impact on radiative and intrinsic Auger-Meitner recombination. Evidence is presented that the imbalance between the change of radiative and intrinsic Auger-Meitner recombination is the primary cause of the efficiency droop at high carrier densities in the samples studied.

Cast Monel alloys are used in applications requiring a combination of good mechanical properties and excellent resistance to corrosion. Despite prevalent industrial use, relatively few studies have been conducted to investigate the relationships between composition, solidification behavior, and microstructure. Given that these alloys are used in the cast and welded conditions, these factors have a significant influence over the material properties. Here, in this work, microstructural characterization, electron probe microanalysis, X-ray diffraction, and differential scanning calorimetry were used to study how changes in Si and Nb concentrations affected the solidification path and microstructure of Monel alloys. It was found that increasing Nb concentration stabilized higher amounts of MC carbides and suppressed graphite formation during solidification. It was also found that the high nominal concentration and segregation of Si to the liquid led to the formation of Ni₃₁Si₁₂ and other silicides via terminal eutectic reactions at the end of solidification. A pseudo-binary solidification diagram was constructed using experimental data and was applied to predict the mass fraction of solidified eutectic as a function of composition. The modeled microstructures were found to be in good agreement with experimentally measured phase fractions.

Here, we perform all-atom molecular dynamics simulations of lithium triflate in 1,2-dimethoxyethane using six different literature force fields. This system is representative of many experimental studies of lithium salts in solvents and polymers. We show that multiple historically common force fields for lithium ions give qualitatively incorrect results when compared with those from experiments and quantum chemistry calculations. We illustrate the importance of correctly selecting force field parameters and give recommendations on the force field choice for lithium electrolyte applications.

Alkene ozonolysis generates short-lived Criegee intermediates that are a significant source of hydroxyl (OH) radicals. This study demonstrates that roaming of the separating OH radicals can yield alternate hydroxycarbonyl products, thereby reducing the OH yield. Specifically, hydroxybutanone has been detected as a stable product arising from roaming in the unimolecular decay of the methyl-ethyl-substituted Criegee intermediate (MECI) under thermal flow cell conditions. The dynamical features of this novel multistage dissociation plus a roaming unimolecular decay process have also been examined with ab initio kinetics calculations. Experimentally, hydroxybutanone isomers are distinguished from the isomeric MECI by their higher ionization threshold and distinctive photoionization spectra. Moreover, the exponential rise of the hydroxybutanone kinetic time profile matches that for the unimolecular decay of MECI. A weaker methyl vinyl ketone (MVK) photoionization signal is also attributed to OH roaming. Complementary multireference electronic structure calculations have been utilized to map the unimolecular decay pathways for MECI, starting with 1,4 H atom transfer from a methyl or methylene group to the terminal oxygen, followed by roaming of the separating OH and butanonyl radicals in the long-range region of the potential. Roaming via reorientation and the addition of OH to the vinyl group of butanonyl is shown to yield hydroxybutanone, and subsequent C-O elongation and H-transfer can lead to MVK. A comprehensive theoretical kinetic analysis has been conducted to evaluate rate constants and branching yields (ca. 10-11%) for thermal unimolecular decay of MECI to conventional and roaming products under laboratory and atmospheric conditions, consistent with the estimated experimental yield (ca. 7%).

As a part of NASA's efforts in space, options are being examined for an Artemis moon base project to be deployed. This project requires a system of interconnected, but separate, DC microgrids for habitation, mining, and fuel processing. This in-place use of power resources is called in-situ resource utilization (ISRU). These microgrids are to be separated by 9-12 km and each contains a photovoltaic (PV) source, energy storage systems (ESS), and a variety of loads, separated by level of criticality in operation. The separate microgrids need to be able to transfer power between themselves in cases where there are generation shortfall, faults, or other failures in order to keep more critical loads running and ensure safety of personnel and the success of mission goals. In this work, a 2 grid microgrid system is analyzed involving a habitation unit and a mining unit separated by a tie line. A set of optimal controls that has been developed, including power flow controls on the tie line, dispatch of PV generation, and dispatch of non-critical loads, is analyzed, and validated in hardware on the Secure Scalable Microgrid Testbed (SSMTB). This testbed includes hardware emulators for a variety of energy sources, energy storage devices, pulsed loads, and other loads.

This report describes specific activities in the Fiscal Year (FY) 2023 associated with the Geologic Disposal Safety Assessment (GDSA) Repository Systems Analysis (RSA) work package funded by the Spent Fuel and Waste Science and Technology (SFWST) Campaign of the U.S. Department of Energy Office of Nuclear Energy (DOE-NE), Office of Spent Fuel and Waste Disposition (SFWD).

Experimental studies and ab initio quantum chemistry calculations were combined to investigate the process by which a Fenton reaction breaks down polystyrene sulfonate. The experimental results show that both molecular weight reduction and loss of aromaticity occur nearly simultaneously, a finding that is supported by the calculations. The results show that more than half of the material is broken down to low molecular weight compounds (< 500 g/mol) with two molar equivalents of H2O2 per styrene monomer. The calculations provide insights into the reaction pathways and indicate that at least two hydroxyl radicals are required to cleave backbone C–C bonds or to eliminate aromaticity. The calculations also show that, of the aromatic carbons, hydroxyl radical is most likely to add to the carbon bonded to sulfur. This finding explains the loss of hydrogen sulfite anion early in the process and also the efficient reduction of Fe(III) to Fe(II) through semiquinone formation. Taken together the experimental and computational results indicate that the reaction is very efficient and that very little H2O2 is lost to unproductive reactions. This high efficiency is attributed to the close association of Fe atoms with the sulfonate group such that hydroxyl radicals are generated near the polymer chains.

Abstract not provided.

A direct numerical simulation (DNS) campaign is deployed for a series of confined downward oriented, non-isothermal turbulent impinging jet configurations. A baseline Reynolds number of 9960 is obtained through a precursor DNS pipe flow simulation (Reτ=505). Three jet temperature configurations (confinement height to nozzle diameter of three) enter a cylindrical domain that share ambient and impingement plate temperatures (298.15K). The range of jet temperatures are crafted such that the ratio of inlet to ambient density varies from unity to 0.52, showcasing the effect of density disparity on flow characteristics such as core collapse, radial mixing of momentum and energy, near-wall stagnation behavior, wall-jet profiles, and large-scale vortical structures. Surface quantities provided include mean radial heat flux and wall-shear stress profiles, and heat flux histograms at select radial stations. Results showcase increased radial normal stresses for higher temperature jets that support increased mixing, resulting in large-scale recirculation structures that are smaller, while retaining similar normalized radial wall profiles for shear stress, heat flux and pressure. Radial plots for wall shear stress and Nusselt number showcase strong radial decay as compared to previous configurations that share similar jet and ambient temperatures. For the 373.15 K case, a Gaussian-like histogram for heat fluxes at the impingement plate transitions to a log-normal profile as radial distances increase. In contrast, the 573.15 K configuration displays a bi-modal heat flux characteristic at the impingement plate, and in similar manner to the moderate temperature counterpart, transitions to a log-normal profile at larger radial distances.

Abstract not provided.

Update to prior 5.14 user manual. I think updates are minor and mostly in the Johnson-cook section. In there those updates are more writing and less on technical changes.

Crystal plasticity finite element model (CPFEM) is a powerful numerical simulation in the integrated computational materials engineering toolboxes that relates microstructures to homogenized materials properties and establishes the structure–property linkages in computational materials science. However, to establish the predictive capability, one needs to calibrate the underlying constitutive model, verify the solution and validate the model prediction against experimental data. Bayesian optimization (BO) has stood out as a gradient-free efficient global optimization algorithm that is capable of calibrating constitutive models for CPFEM. In this paper, we apply a recently developed asynchronous parallel constrained BO algorithm to calibrate phenomenological constitutive models for stainless steel 304 L, Tantalum, and Cantor high-entropy alloy.

Calibrated measurements of lightning optical emissions are critical for both quantifying the impacts of lightning in our atmosphere and devising detection instruments with sufficient dynamic range capable of yielding close to 100% detection efficiency. However, to date, there is only a limited number of investigations that have attempted to take such calibrated measurements. In this work, we report the power radiated by lightning in both visible and infrared bands, assuming isotropic emission, and accounting for atmospheric absorption. More precisely, we report peak radiated power and total radiated energy in the combined visible plus near-infrared range (VNIR, 0.34–1.1 μm), around the Hα line (652–667 nm), and for the 2–2.5 μm infrared band. The estimated peak power and total energy radiated by negative cloud-to-ground return strokes in the VNIR range is 130 MW and 20 kJ, respectively. Additionally, we detected peak radiated powers of 12 and 0.19 MW in the Hα and infrared bands, respectively. We cross-reference the optical data set with peak current reported by a lightning detection network. The resulting trend is that optical power emitted around the Hα line scales with peak return stroke current according to a power law with exponent equal to 1.25. This trend, which should be approximately true across the entire visible spectrum, can be attributed to the plasma negative differential resistance of the lightning return stroke channel. We conclude by discussing the challenges in performing calibrated measurements of lightning optical power in different bands and comparing the results with previously-collected data with different experimental setups, observation conditions, and calibration methods.

We propose a novel statistical inference methodology for multiway count data that is corrupted by false zeros that are indistinguishable from true zero counts. Our approach consists of zero-truncating the Poisson distribution to neglect all zero values. This simple truncated approach dispenses with the need to distinguish between true and false zero counts and reduces the amount of data to be processed. Inference is accomplished via tensor completion that imposes low-rank tensor structure on the Poisson parameter space. Our main result shows that an N-way rank-R parametric tensor M ∈ (0, ∞)I×.....×I generating Poisson observations can be accurately estimated by zero-truncated Poisson regression from approximately IR2 log22(I) non-zero counts under the nonnegative canonical polyadic decomposition. Our result also quantifies the error made by zero-truncating the Poisson distribution when the parameter is uniformly bounded from below. Therefore, under a low-rank multiparameter model, we propose an implementable approach guaranteed to achieve accurate regression in under-determined scenarios with substantial corruption by false zeros. Several numerical experiments are presented to explore the theoretical results.

Previous studies have shown that selective laser flash sintering (SLFS) can be initiated in dielectrics that exhibit ionic or electronic conduction at high temperature. These materials required high laser powers to reach the temperatures where electrical conduction is sufficient to initiate SLFS. In this study, SLFS in lanthanum chromite (LC), an intrinsic electronic conductor with high conductivity, and lanthanum strontium chromite (LSC), which is doped to further increase electronic conductivity, were investigated with a focus on understanding the initiation mechanisms. Results show that the initiation of SLFS in LC and LSC occurs when electronic charge carriers are activated and flow to the electrode where the current is measured. A combination of carriers produced at the electrode, temperature-activated intrinsic charge carriers, and extrinsic charge carriers present in LSC due to doping are responsible for the facile initiation of SLFS.

Abstract not provided.

Abstract not provided.

Achieving brain-like efficiency in computing requires a co-design between the development of neural algorithms, brain-inspired circuit design, and careful consideration of how to use emerging devices. The recognition that leveraging device-level noise as a source of controlled stochasticity represents an exciting prospect of achieving brain-like capabilities in probabilistic neural algorithms, but the reality of integrating stochastic devices with deterministic devices in an already-challenging neuromorphic circuit design process is formidable. Here, we explore how the brain combines different signaling modalities into its neural circuits as well as consider the implications of more tightly integrated stochastic, analog, and digital circuits. Further, by acknowledging that a fully CMOS implementation is the appropriate baseline, we conclude that if mixing modalities is going to be successful for neuromorphic computing, it will be critical that device choices consider strengths and limitations at the overall circuit level.

Code verification plays an important role in establishing the credibility of computational simulations by assessing the correctness of the implementation of the underlying numerical methods. In computational electromagnetics, the numerical solution to integral equations incurs multiple interacting sources of numerical error, as well as other challenges, which render traditional code-verification approaches ineffective. In this paper, we provide approaches to separately measure the numerical errors arising from these different error sources for the method-of-moments implementation of the combined-field integral equation. We demonstrate the effectiveness of these approaches for cases with and without coding errors.

Global climate change has prompted many national plans for rapid emissions reductions. For example, the United States recently committed to transitioning to 100% carbon-free electricity by 2035 and net-zero emissions economy-wide by 2050. Parallel to conversations surrounding emissions reductions is the call for energy justice, or the demand for more equitable distribution of energy-related burdens and benefits among communities. To date, energy justice has evolved as a mostly academic conversation, which may limit its utility to praxis. In response, we offer an interdisciplinary framework that aims to organize existing knowledge and lessons learned from energy development. Specifically, we developed the Meaningful Marine Renewable Energy (MRE) Development Framework and conducted a literature review using MRE as a case study. MRE was chosen because it is a nascent renewable energy technology in the US with projects mostly in demonstration stages and no commercial deployment, making it a useful case study to apply lessons learned from other energy sectors and other countries. Discussion of current resources being developed among the MRE community and their implications for furthering energy justice priorities are also explored. We conclude the review with a compiled list of questions meant to support stakeholders in translating theoretical concepts of Meaningful MRE Development to practice. Although the Meaningful MRE framework was developed using MRE as a use case, our interdisciplinary theoretical framework can be applied beyond MRE to other sustainable and renewable energy projects.

Fracture prognosis and characterization efforts require knowledge of crack tip position and the Stress Intensity Factors (SIFs) acting in the vicinity of the crack. Here, we present an efficient numerical approach to infer both of these characteristics under a consistent theoretical framework from noisy, unstructured displacement data. The novel approach utilizes the separability of the asymptotic linear elastic fracture mechanics fields to expedite the search for crack tip position and is particularly useful for noisy displacement data. The manuscript begins with an assessment of the importance of accurately locating crack tip position when quantifying the SIFs from displacement data. Next, the proposed separability approach for quickly inferring crack tip position is introduced. Comparing to the widely used displacement correlation approach, the performance of the separability approach is assessed. Cases involving both noisy data and systematic deviation from the asymptotic linear elastic fracture mechanics model are considered, e.g. inelastic material behavior and finite geometries. An open source python implementation of the proposed approach is available for use by those doing field and laboratory work involving digital image correlation and simulations, e.g. finite element, discrete element, molecular dynamics and peridynamics, where the crack tip position is not explicitly defined.

Abstract not provided.

This study presents a deep learning based methodology for both remote sensing and design of acoustic scatterers. The ability to determine the shape of a scatterer, either in the context of material design or sensing, plays a critical role in many practical engineering problems. This class of inverse problems is extremely challenging due to their high-dimensional, nonlinear, and ill-posed nature. To overcome these technical hurdles, we introduce a geometric regularization approach for deep neural networks (DNN) based on non-uniform rational B-splines (NURBS) and capable of predicting complex 2D scatterer geometries in a parsimonious dimensional representation. Then, this geometric regularization is combined with physics-embedded learning and integrated within a robust convolutional autoencoder (CAE) architecture to accurately predict the shape of 2D scatterers in the context of identification and inverse design problems. An extensive numerical study is presented in order to showcase the remarkable ability of this approach to handle complex scatterer geometries while generating physically-consistent acoustic fields. The study also assesses and contrasts the role played by the (weakly) embedded physics in the convergence of the DNN predictions to a physically consistent inverse design.

Abstract not provided.

Abstract not provided.

Abstract not provided.

In polymer-filled granular composites, damage may develop in mechanical loading prior to material failure. Damage mechanisms such as microcracking or plastic deformation in the binder phase can substantially alter the material's mesostructure. For energetic materials, such as solid propellants and plastic bonded explosives, these mesostructural changes can have far reaching effects including degraded mechanical properties, potentially increased sensitivity to further insults, and changes in expected performance. Unfortunately, predicting damage is nontrivial due to the complex nature of these composites and the entangled interactions between inelastic mechanisms. In this work, we assess the current literature of experimental knowledge, focusing on the pressure-dependent shear response, and propose a simple simulation framework of bonded particles to study four limiting-case material formulations at both meso- and macro-scales. To construct the four cases, we systematically vary the relative interfacial strength between the polymer binder and granular filler phase and also vary the polymer's glass transition temperature relative to operating temperature which determines how much the binder can plastically deform. These simulations identify key trends in global mechanical response, such as the emergence of strain hardening or softening regimes with increasing pressure which qualitatively resemble experimental results. By quantifying the activation of different inelastic mechanisms, such as bonds breaking and plastically straining, we identify when each mechanism becomes relevant and provide insight into potential origins for changes in mechanical responses. The locations of broken bonds are also used to define larger, mesoscopic cracks to test various metrics of damage. We primarily focus on triaxial compression, but also test the opposite case of triaxial extension to highlight the impact of Lode angle on mechanical behavior.

Salt offers an optimal medium for the permanent isolation of heat-producing radioactive waste due to its impermeability, high thermal conductivity, and ability to close fractures through creep. A thorough understanding of the thermal-hydrological-mechanical (THM) processes, encompassing brine migration, is fundamental for secure radioactive waste disposal within salt formations. At the Waste Isolation Pilot Plant (WIPP), we conducted joint in situ geophysical monitoring experiments during active heating to investigate brine migration near excavations. This experiment incorporated electrical resistivity tomography (ERT) alongside high-resolution fiber-optic-based distributed temperature sensing within a controlled heating experiment. Additionally, discrete element model (DEM) based numerical simulations were conducted to simulate THM processes during heating, providing a more mechanistic understanding of the coupled processes leading to the observed changes in geophysical measurements. During heating, resistivity shifts near the heater were reasonably explained by temperature effects. However, in more distant, cooler regions, the resistivity decrease exceeded predictions based solely on temperature. DEM simulations highlighted brine migration, propelled by pore pressure gradients, as the likely primary factor contributing to the additional resistivity decline beyond temperature effects. The comparison between the predicted ERT responses and observations was much improved when considering the effects of brine migration based on the DEM simulations. These geophysical and simulation findings shed light on brine migration in response to salt heating, enhancing our understanding of the coupled THM processes in salt for safe radioactive waste disposal.

Abstract not provided.

Deflection and stress calculated from an experimentally validated, high-fidelity finite element model (FEM) of a photovoltaic module experiencing mechanical load was compared to results from a simplified FEM treating the module laminate as a homogenized composite using a rule of mixtures approach, and further compared to analytical calculations treating the module as a Kirchoff-Love flat plate. The goal of this study was to determine the error incurred by analyzing module mechanics with varying levels of simplification, since resolving the aspect ratios of a module is computationally expensive. Homogenized FEMs were found to underpredict peak deflection under a 1.0 kPa load by between 13 and 19% for lower and upper bound application of the rule of mixtures. However, module shape was captured, implying that a useful replication of a resolved model could be achieved with a reduced, calibrated material stiffness. Homogenized stress results captured glass layer tensile stress components to within 46 to 52% at a sample location of interest, though agreement was poor through the remainder of the laminate due to the lack of material resolution. For plate theory, deflection was overpredicted by 45 to 67% for upper and lower bound homogenizations, and frame-adjacent module shapes were not adequately replicated. Stress results mirrored FEM trends but magnitudes were not well correlated to resolved model values. These results support the use of homogenized laminate models for module shape derivation, though resolved models remain necessary for stress analyses. The accuracy of plate theory was found to be inadequate for most applications.

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

Abstract not provided.

Abstract not provided.

Efficient operation of battery energy storage systems requires that battery temperature remains within a specific range. Current techno-economic models neglect the parasitic loads heating and cooling operations have on these devices, assuming they operate at constant temperature. In this work, these effects are investigated considering the optimal sizing of battery energy storage systems when deployed in cold environments. A peak shaving application is presented as a linear programming problem which is then formulated in the PYOMO optimization programming language. The building energy simulation software EnergyPlus is used to model the heating, ventilation, and air conditioning load of the battery energy storage system enclosure. Case studies are conducted for eight locations in the United States considering a nickel manganese cobalt oxide lithium ion battery type and whether the power conversion system is inside or outside the enclosure. The results show an increase of 42% to 300% in energy capacity size, 43% to 217% in power rating, and 43% to 296% increase in capital cost dependent on location. This analysis shows that the heating, ventilation, and air conditioning load can have a large impact on the optimal sizes and cost of a battery energy storage system and merit consideration in techno-economic studies.

The carbon phase diagram is rich with polymorphs which possess very different physical and optical properties ideal for different scientific and engineering applications. An understanding of the dynamically driven phase transitions in carbon is particularly important for applications in inertial confinement fusion, as well as planetary and meteorite impact histories. Experiments on the Z Pulsed Power Facility at Sandia National Laboratories generate dynamically compressed high-pressure states of matter with exceptional uniformity, duration, and size that are ideal for investigations of fundamental material properties. X-ray diffraction (XRD) is an important material physics measurement because it enables direct observation of the strain and compression of the crystal lattice, and it enables the detection and identification of phase transitions. Several unique challenges of dynamic compression experiments on Z prevent using XRD systems typically utilized at other dynamic compression facilities, so novel XRD diagnostics have been designed and implemented. We performed experiments on Z to shock compress carbon (pyrolytic graphite) samples to pressures of 150–320 GPa. The Z-Beamlet Laser generated Mn-Heα (6.2 keV) X-rays to probe the shock-compressed carbon sample, and the new XRD diagnostics measured changes in the diffraction pattern as the carbon transformed into its high-pressure phases. Quantitative analysis of the dynamic XRD patterns in combination with continuum velocimetry information constrained the stability fields and melting of high-pressure carbon polymorphs.

The nitrogen-vacancy (NV) color center in diamond has demonstrated great promise in a wide range of quantum sensing. Recently, there have been a series of proposals and experiments using NV centers to detect spin noise of quantum materials near the diamond surface. This is a rich complex area of study with novel nano-magnetism and electronic behavior, that the NV center would be ideal for sensing. However, due to the electronic properties of the NV itself and its host material, getting high quality NV centers within nanometers of such systems is challenging. Band bending caused by space charges formed at the metal-semiconductor interface force the NV center into its insensitive charge states. Here, we investigate optimizing this interface by depositing thin metal films and thin insulating layers on a series of NV ensembles at different depths to characterize the impact of metal films on different ensemble depths. We find an improvement of coherence and dephasing times we attribute to ionization of other paramagnetic defects. The insulating layer of alumina between the metal and diamond provide improved photoluminescence and higher sensitivity in all modes of sensing as compared to direct contact with the metal, providing as much as a factor of 2 increase in sensitivity, decrease of integration time by a factor of 4, for NV T 1 relaxometry measurements.

Abstract not provided.

Abstract not provided.

Quantum point contacts (QPC) are the building blocks of quantum dot qubits and semiconducting quantum electrical metrology circuits. QPCs also make highly sensitive electrical amplifiers with the potential to operate in the quantum-limited regime. Though the inherent operational bandwidth of QPCs can eclipse the THz regime, the impedance mismatch with the external circuitry limits the operational frequency to a few kHz. Lumped-element impedance-matching circuits are successful only up to a few hundreds of MHz in frequency. QPCs are characterised by a complex impedance consisting of quantized resistance, capacitance, and inductance elements. Characterising the complex admittance at higher frequencies and understanding the coupling of QPC to other circuit elements and electromagnetic environments will provide valuable insight into its sensing and backaction properties. In this work, we couple a QPC galvanically to a superconducting stub tuner impedance matching circuit realised in a coplanar waveguide architecture to enhance the operation frequency into the GHz regime and investigate the electrical amplification and complex admittance characteristics. The device, operating at ~1.96 $GHz$ exhibits a conductance sensitivity of 2.92 X 10^-5(e²/h)/$\sqrt{Hz}$ with a bandwidth of 13 $MHz$. Besides, the RF reflected power unambiguously reveals the complex admittance characteristics of the QPC, shining more light on the behaviour of quantum tunnel junctions at higher operational frequencies.

The ion velocity distribution functions of thermonuclear plasmas generated by spherical laser direct drive implosions are studied using deuterium-tritium (DT) and deuterium-deuterium (DD) fusion neutron energy spectrum measurements. A hydrodynamic Maxwellian plasma model accurately describes measurements made from lower temperature (<10 keV), hydrodynamiclike plasmas, but is insufficient to describe measurements made from higher temperature more kineticlike plasmas. The high temperature measurements are more consistent with Vlasov-Fokker-Planck (VFP) simulation results which predict the presence of a bimodal plasma ion velocity distribution near peak neutron production. These measurements provide direct experimental evidence of non-Maxwellian ion velocity distributions in spherical shock driven implosions and provide useful data for benchmarking kinetic VFP simulations.

The goal of quantum tracking control is to identify shaped fields to steer observable expectation values along designated time-dependent tracks. The fields are determined via an iteration-free procedure, which is based on inverting the underlying dynamical equations governing the controlled observables. In this paper, we generalize the ideas in [Phys. Rev. A 98, 043429 (2018)2469-992610.1103/PhysRevA.98.043429] to the task of orienting symmetric top molecules in three dimensions. To this end, we derive equations for the control fields capable of directly tracking the expected value of the three-dimensional dipole orientation vector along a desired path in time. We show this framework can be utilized for tracking the orientation of linear molecules as well, and present numerical illustrations of these principles for symmetric-top tracking control problems.

A bonded particle model is used to explore how variations in the material properties of brittle, isotropic solids affect critical behavior in fragmentation. To control material properties, a model is proposed which includes breakable two- and three-body particle interactions to calibrate elastic moduli and mode I and mode II fracture toughnesses. In the quasistatic limit, fragmentation leads to a power-law distribution of grain sizes which is truncated at a maximum grain mass that grows as a nontrivial power of system size. In the high-rate limit, truncation occurs at a mass that decreases as a power of increasing rate. A scaling description is used to characterize this behavior by collapsing the mean-square grain mass across rates and system sizes. Consistent scaling persists across all material properties studied, although there are differences in the evolution of grain size distributions with strain as the initial number of grains at fracture and their subsequent rate of production depend on Poisson's ratio. This evolving granular structure is found to induce a unique rheology where the ratio of the shear stress to pressure, an internal friction coefficient, decays approximately as the logarithm of increasing strain rate. The stress ratio also decreases at all rates with increasing strain as fragmentation progresses and depends on elastic properties of the solid.

Abstract not provided.

Abstract not provided.

Although visualizations are a useful tool for helping people to understand information, they can also have unintended effects on human cognition. This is especially true for uncertain information, which is difficult for people to understand. Prior work has found that different methods of visualizing uncertain information can produce different patterns of decision making from users. However, uncertainty can also be represented via text or numerical information, and few studies have systematically compared these types of representations to visualizations of uncertainty. We present two experiments that compared visual representations of risk (icon arrays) to numerical representations (natural frequencies) in a wildfire evacuation task. Like prior studies, we found that different types of visual cues led to different patterns of decision making. In addition, our comparison of visual and numerical representations of risk found that people were more likely to evacuate when they saw visualizations than when they saw numerical representations. These experiments reinforce the idea that design choices are not neutral: seemingly minor differences in how information is represented can have important impacts on human risk perception and decision making.

Abstract not provided.

Abstract not provided.

The injection and storage of anthropogenic CO2 in the subsurface is being deployed as a climate change mitigation tool; however, diagenetic-paragenetic heterogeneity in sandstone reservoirs often contributes to interval specific chemomechanical changes that affect injection and can increase leakage risk. Here, we address reservoir heterogeneities’ impact on chemomechanical changes in a macroporous-dominated lithofacies of Morrow B sandstone, a formation containing several diagenetically-distinct hydraulic facies while undergoing enhanced oil recovery (EOR) and carbon dioxide (CO2) sequestration. We performed three flow-through experiments using a CO2-charged or uncharged formation water combined with four indirect tensile strength tests per post-test sample. We then used the microstructure and paragenetic sequence to understand chemomechanical weakening with key observations as follows: dissolution of carbonates and feldspars changed porosity; increased permeability led to reclassifying each sample in a different hydraulic flow unit; decreased ultrasonic velocity; and did not lead to a loss of tensile strength. Tensile strength maintenance occurred due to the low abundance and minor dissolution of siderite, the stability of quartz, and the relative position of diagenetic ankerite within feldspar. This macroporous-dominated lithofacies is the primary reservoir for the Morrow B Sandstone, and is analogous to other porous sandstone reservoirs. It represents an end-member of a chemomechanically low-risk siliceous CO2 sequestration and CO2-EOR reservoir.

Living document updating the 5.14 version, SAND2023-03603 of May 2023.

Frontal polymerization involves the propagation of a thermally driven polymerization wave through a monomer solution to rapidly generate high-performance polymeric materials with little energy input. The balance between latent catalyst activation and sufficient reactivity to sustain a front can be difficult to achieve and often results in systems with poor storage lives. This is of particular concern for frontal ring-opening metathesis polymerization (FROMP) where gelation occurs within a single day of resin preparation due to the highly reactive nature of Grubbs-type catalysts. In this report we demonstrate the use of encapsulated catalysts to provide remarkable latency to frontal polymerization systems, specifically using the highly active dicyclopentadiene monomer system. Negligible differences were observed in the frontal velocities or thermomechanical properties of the resulting polymeric materials. FROMP systems with encapsulated catalyst particles are shown with storage lives exceeding 12 months and front rates that increase over a well-characterized 2 month period. Moreover, the modularity of this encapsulation method is demonstrated by encapsulating a platinum catalyst for the frontal polymerization of silicones by using hydrosilylation chemistry.

We explore the use of reduced physics models for efficient kinetic particle simulations of space charge limited (SCL) emission in inner magnetically insulated transmission lines (inner MITLs), with application to Sandia National Laboratories' Z machine. We propose a drift kinetic (guiding center) model of electron motion in place of a fully kinetic model and electrostatic-magnetostatic fields in place of electromagnetic fields. The validity of these approximations is suggested by the operational parameters of the Z machine, namely, current pulse lengths of order 100 ns compared with Larmor periods typically smaller than 10-11 s, typical Larmor radii of a few (tens) of microns (magnetic fields of tens to hundreds of Tesla) compared with MITL dimensions of a few centimeters, and transient time of light waves along the inner MITL of order a fraction of a nanosecond. Guiding center orbits eliminate the fast electron gyromotion, which enables the use of tens to hundreds of times larger time steps in the numerical particle advance. Electrostatic-magnetostatic fields eliminate the Courant-Friedrichs-Lewy (CFL) numerical stability limit on the time step and allow the use of higher grid resolutions or, alternatively, larger time steps in the fields advance. Overall, potential computational cost savings of tens to hundreds of times exists. The applicability of the reduced physics models is examined on two problems. First, in the simulation of space charge limited emission of electrons from the cathode surface due to high electric fields in a radial inner MITL geometry with a short load. In particular, it is shown that a drift kinetic-based particle-in-cell (PIC) model with electrostatic-magnetostatic fields is able to accurately reproduce well-known physics of electron vortex formation, spatially and temporally. Second, deeper understanding is gained of the mechanism behind vortex formation in this MITL geometry by considering an exemplar problem of an electron block of charge. This simpler setup reveals that the main mechanism of vortex formation can be attributed to pure drift motion of the electrons, that is, the (fully kinetic) gyromotion of the electrons is inessential to the process. This exemplar problem also suggests a correlation of the spatial dimensions of vortices to the thickness of the electron layer, as observed in SCL simulations. It also confirms that the electromagnetic nature of the fields does not play an essential role. Finally, an improved hybrid fully kinetic and drift kinetic model for electron motion is proposed, as means of capturing finite Larmor radius (FLR) effects; the particular FLR physics that is missed by the drift kinetic model is the particle-wall interaction. By initializing SCL emitted electrons as fully kinetic and later transitioning them to drift kinetic, according to simple criteria, the accuracy of SCL simulations can be improved, while preserving the potential for computational efficiency.

Abstract not provided.

The wave energy control competition established a benchmark problem which was offered as an open challenge to the wave energy system control community. The competition had two stages: In the first stage, competitors used a standard wave energy simulation platform (WEC-Sim) to evaluate their controllers while, in the second stage, competitors were invited to test their controllers in a real-time implementation on a prototype system in a wave tank. The performance function used was based on converted energy across a range of standard sea states, but also included aspects related to economic performance, such as peak/average power, peak force, etc. This paper compares simulated and experimental results and, in particular, examines if the results obtained in a linear system simulation are borne out in reality. Overall, within the scope of the device tested, the range of sea states employed, and the performance metric used, the conclusion is that high-performance WEC controllers work well in practice, with good carry-over from simulation to experimentation. However, the availability of a good WEC mathematical model is deemed to be crucial.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

We demonstrate the discrimination of ground-state hyperfine manifolds of a cesium atom in an optical tweezer using a simple probe beam with Formula Presented% detection fidelity and 0.9(2)% detection-driven loss of bright-state atoms. Our detection infidelity of Formula Presented% is an order of magnitude better than previously published low-loss readout results for alkali-metal atoms in optical tweezers. We achieve these results by identifying and mitigating an extra depumping mechanism due to stimulated Raman transitions induced by trap light in the presence of probe light. In this work, complex optical systems and stringent vacuum pressures are not required, enabling straightforward adoption of our techniques on contemporary experiments.

Abstract not provided.

This project applies methods in Bayesian inference and modern statistical methods to quantify the value of new experimental data, in the form of new or modified diagnostic configurations and/or experiment designs. We demonstrate experiment design methods that can be used to identify the highest priority diagnostic improvements or experimental data to obtain in order to reduce uncertainties on critical inferred experimental quantities and select the best course of action to distinguish between competing physical models. Bayesian statistics and information theory provide the foundation for developing the necessary metrics, using two high impact experimental platforms on Z as exemplars to develop and illustrate the technique. We emphasize that the general methodology is extensible to new diagnostics (provided synthetic models are available), as well as additional platforms. We also discuss initial scoping of additional applications that began development in the last year of this LDRD.

This paper describes a process for forming a buried field shield in GaN by an etch-and-regrowth process, which is intended to protect the gate dielectric from high fields in the blocking state. GaN trench MOSFETs made at Sandia serve as the baseline to show the limitations in making a trench gated device without a method to protect the gate dielectric. Device data coupled with simulations show device failure at 30% of theoretical breakdown for devices made without a field shield. Implementation of a field shield reduces the simulated electric field in the dielectric to below 4 MV/cm at breakdown, which eliminates the requirement to derate the device in order to protect the dielectric. For realistic lithography tolerances, however, a shield-to-channel distance of 0.4 μm limits the field in the gate dielectric to 5 MV/cm and requires a small margin of device derating to safeguard a long-term reliability and lifetime of the dielectric.

Abstract not provided.

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

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

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

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

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

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

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

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

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

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

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

Interlocking metasurfaces (ILMs) are a new class of mechanical metasurfaces built from architected arrays of interlocking features that can serve as a nonpermanent, robust joining technology. An ILM's strength is governed by the structural material, orientation, and topology of its latching unit cells. The presented work optimized the topologies of ILM unit cells to maximize strength in tensile and shear loading using gradient-based parametric optimization and genetic algorithms. Experimental validation confirmed that the optimized designs achieved considerable strength increases compared to a human intuitive design. In several cases, the optimized designs were approximately double the effective interfacial strength compared to that achieved via expert intuition alone. The strength improvement was seen for isolated unit cells and arrays of interacting unit cells (metasurfaces). An analysis of the topologies of the optimized designs showed that tall dendritic geometric features with large contact surfaces yield robust solutions in tension, while short and broad geometric features with large contact surfaces yield better results in shear loading. This study revealed the importance of shape optimization to maximize ILM effectiveness under single- and multi-objective scenarios.

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

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

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

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

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

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