Publications

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

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

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

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

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

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.

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

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

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

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.

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.

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

Abstract not provided.

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.

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

Abstract not provided.

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

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

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

Abstract not provided.

Abstract not provided.

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

Abstract not provided.

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

Abstract not provided.