Publications

This report summarizes important nuances in local water concerns and potential climate impacts that could influence the roll-out of technologies associated with energy transitions. To understand how water and climate dynamics could be influencing these activities for three countries.

Peridigm is a meshfree peridynamics code written in C++ for use on large-scale parallel computers. It was originally developed at Sandia National Laboratories and is currently managed as an open-source, community driven software project. Its primary features include bond-based, state-based, and non-ordinary state-based constitutive models, bond failure laws, contact, and support for explicit and implicit time integration. To date, Peridigm has been used primarily by methods developers focused on solid mechanics and material failure. Peridigm utilizes foundational software components from Sandia’s Trilinos project and was designed for extensibility. This paper provides an overview of the solution methods implemented in Peridigm, a discussion of its software infrastructure, and demonstrates the use of Peridigm for the solution of several example problems.

Carbon fiber epoxy composites are increasingly used in systems requiring a material that is both strong and light weight, as in airplanes, cars, and pressure vessels. In fire environments, carbon fiber epoxy composites are a fuel source subject to oxidation. This literature review seeks to provide material properties as well as uncertainty bounds for those properties for computational models of decomposing carbon fiber epoxy composites. The goal is to guide analysts when measurements are lacking

Measurement

Ground heat flux (G0) is a key component of the land-surface energy balance of high-latitude regions. Despite its crucial role in controlling permafrost degradation due to global warming, G0 is sparsely measured and not well represented in the outputs of global scale model simulation. In this study, an analytical heat transfer model is tested to reconstruct G0 across seasons using soil temperature series from field measurements, Global Climate Model, and climate reanalysis outputs. The probability density functions of ground heat flux and of model parameters are inferred using available G0 data (measured or modeled) for snow-free period as a reference. When observed G0 is not available, a numerical model is applied using estimates of surface heat flux (dependent on parameters) as the top boundary condition. These estimates (and thus the corresponding parameters) are verified by comparing the distributions of simulated and measured soil temperature at several depths. Aided by state-of-the-art uncertainty quantification methods, the developed G0 reconstruction approach provides novel means for assessing the probabilistic structure of the ground heat flux for regional permafrost change studies.

Current nuclear facility emergency planning zones (EPZs) are based on outdated distance-based criteria, predating comprehensive dose and risk-informed frameworks. Recent advancements in simulation tools have permitted the development of site-specific, dose, and risk-based consequence-driven assessment frameworks. This study investigated the computation of advanced reactor (AR) EPZs using two atmospheric dispersion models: a straight-line Gaussian plume model (GPM) and a semi-Lagrangian Particle in Cell (PIC). Two case studies were conducted: (1) benchmarking the NRC SOARCA study for the Peach Bottom Nuclear Generating Station and (2) analyzing an advanced INL Heat Pipe Design A microreactor's end-of-cycle inventory. The dose criteria for both cases were 10 mSv at mean weather conditions and 50 mSv at 95th percentile weather conditions at 96 h post-release. Results demonstrated that GPM and PIC estimated similar mean peak dose levels for large boiling water reactors in the farfield case, placing EPZ limits beyond current regulations. For ARs with source terms remaining in the nearfield, PIC modeling without specific nearfield considerations could result in excessively high doses and inaccurate EPZ designations. PIC dispersion demonstrated an order of magnitude higher estimate of nearfield inhalation dose contribution when compared to GPM results. Both models significantly reduced EPZ sizing within the nearfield. Thus, reductions in the AR source term may eliminate the need for a separate EPZ.

The impact of the high-power impulse magnetron sputtering (HiPIMS) pulse width on the crystallization, microstructure, and ferroelectric properties of undoped HfO2 films is investigated. HfO2 films were sputtered from a hafnium metal target in an Ar/O2 atmosphere, varying the instantaneous power density by changing the HiPIMS pulse width with fixed time-averaged power and pulse frequency. The pulse width is shown to affect the ion-to-neutral ratio in the depositing species with the shortest pulse durations leading to the highest ion fraction. In situ x-ray diffraction measurements during crystallization demonstrate that the HiPIMS pulse width impacts nucleation and phase formation, with an intermediate pulse width of 110 μs stabilizing the ferroelectric phase over the widest temperature range. Although the pulse width impacts the grain size with the lowest pulse width resulting in the largest grain size, the grain size does not strongly correlate with the phase content or ferroelectric behavior in these films. These results suggest that precise control over the energetics of the depositing species may be beneficial for forming the ferroelectric phase in this material.

The complex nature of manufacturing processes stipulates electrodes to possess high variability with increased heterogeneity during production. X-ray computed tomography imaging has proved to be critical in visualizing the complicated stochastic particle distribution of as-manufactured electrodes in lithium-ion batteries. However, accurate prediction of their electrochemical performance necessitates precise evaluation of kinetic and transport properties from real electrodes. Image segmentation that characterizes voxels to particle/pore phase is often meticulous and fraught with subjectivity owing to a myriad of unconstrained choices and filter algorithms. We utilize a Bayesian convolutional neural network to tackle segmentation subjectivity and quantify its pertinent uncertainties. Otsu inter-variance and Blind/Referenceless Imaging Spatial Quality Evaluator are used to assess the relative image quality of grayscale tomograms, thus evaluating the uncertainty in the derived microstructural attributes. We analyze how image uncertainty is correlated with the uncertainties and magnitude of kinetic and transport properties of an electrode, further identifying pathways of uncertainty propagation within microstructural attributes. The coupled effect of spatial heterogeneity and microstructural anisotropy on the uncertainty quantification of transport parameters is also understood. This work demonstrates a novel methodology to extract microstructural descriptors from real electrode images through quantification of associated uncertainties and discerning the relative strength of their propagation, thus facilitating feedback to manufacturing processes from accurate image based electrochemical simulations.

Achieving commercially acceptable Zn-MnO2 rechargeable batteries depends on the reversibility of active zinc and manganese materials, and avoiding side reactions during the second electron reaction of MnO2. Typically, liquid electrolytes such as potassium hydroxide (KOH) are used for Zn-MnO2 rechargeable batteries. However, it is known that using liquid electrolytes causes the formation of electrochemically inactive materials, such as precipitation Mn3O4 or ZnMn2O4 resulting from the uncontrollable reaction of Mn3+ dissolved species with zincate ions. In this paper, hydrogel electrolytes are tested for MnO2 electrodes undergoing two-electron cycling. Improved cell safety is achieved because the hydrogel electrolyte is non-spillable, according to standards from the US Department of Transportation (DOT). The cycling of “half cells” with advanced-formulation MnO2 cathodes paired with commercial NiOOH electrodes is tested with hydrogel and a normal electrolyte, to detect changes to the zincate crossover and reaction from anode to cathode. These half cells achieved ≥700 cycles with 99% coulombic efficiency and 63% energy efficiency at C/3 rates based on the second electron capacity of MnO2. Other cycling tests with “full cells” of Zn anodes with the same MnO2 cathodes achieved ~300 cycles until reaching 50% capacity fade, a comparable performance to cells using liquid electrolyte. Electrodes dissected after cycling showed that the liquid electrolyte allowed Cu ions to migrate more than the hydrogel electrolyte. However, measurements of the Cu diffusion coefficient showed no difference between liquid and gel electrolytes; thus, it was hypothesized that the gel electrolytes reduced the occurrence of Cu short circuits by either (a) reducing electrode physical contact to the separator or (b) reducing electro-convective electrolyte transport that may be as important as diffusive transport.

Magnetized Liner Inertial Fusion experiments have been performed at the Z facility at Sandia National Laboratories. These experiments use deuterium fuel, which produces 2.45 MeV neutrons on reaching thermonuclear conditions. To study the spatial structure of neutron production, the one-dimensional imager of neutrons diagnostic was fielded to record axial resolved neutron images. In this diagnostic, neutrons passing through a rolled edge aperture form an image on a CR-39-based solid state nuclear track detector. Here, we present a modified generalized expectation-maximization algorithm to reconstruct an axial neutron emission profile of the stagnated fusion plasma. We validate the approach by comparing the reconstructed neutron emission profile to an x-ray emission profile provided by a time-integrated pinhole camera.

Heating of the surficial layer of the atmosphere often generates convective vortices, known as “dust devils” when they entrain visible debris. Convective vortices are common on both Earth and Mars, where they affect the climate via dust loading, contribute to wind erosion, impact the efficiency of photovoltaic systems, and potentially result in injury and property damage. However, long-duration terrestrial convective vortex activity records are rare. We have developed a high-precision and high-recall method to extract convective vortex signatures from infrasound microbarometer data streams. The techniques utilizes a wavelet-based detector to capture potential events and then a template matching system to extract the duration of the vortex. Since permanent and temporary infrasound sensors networks are present throughout the globe (many with open data), our method unlocks a vast new convective vortex dataset without requiring the deployment of specialized instrumentation. SIGNIFICANCE STATEMENT: Convective vortices, or “dust devils,” contribute to regional dust loading in Earth’s atmosphere. However, long-duration convective vortex activity records are rare. We came up with a way to autonomously detect the pressure signatures left by convective vortices striking low-frequency sound, or “infrasound,” sensors. Since permanent infrasound stations have been active for decades, our method has the potential to add ordersof-magnitude more events than previously catalogued.

A study was conducted of an intermittent binary control strategy for trailing edge flaps and leading edge spoilers installed on wind turbine blades for the purpose of load alleviation. Cost estimation models were developed for the systems to predict overall impact on levelized cost of energy over the lifecycle of the turbine system. Aeroelastic simulations of turbines with the control strategy implemented showed improved levelized cost for some, but not all cases.

The multi-harmonic balance method combined with numerical continuation provides an efficient framework to compute a family of time-periodic solutions, or response curves, for large-scale, nonlinear mechanical systems. The predictor and corrector steps repeatedly solve a sequence of linear systems that scale by the model size and number of harmonics in the assumed Fourier series approximation. In this paper, a novel Newton–Krylov iterative method is embedded within the multi-harmonic balance and continuation algorithm to efficiently compute the approximate solutions from the sequence of linear systems that arise during the prediction and correction steps. The method recycles, or reuses, both the preconditioner and the Krylov subspace generated by previous linear systems in the solution sequence. A delayed frequency preconditioner refactorizes the preconditioner only when the performance of the iterative solver deteriorates. The GCRO-DR iterative solver recycles a subset of harmonic Ritz vectors to initialize the solution subspace for the next linear system in the sequence. The performance of the iterative solver is demonstrated on two exemplars with contact-type nonlinearities and benchmarked against a direct solver with traditional Newton–Raphson iterations.

During fracture amorphous oxides exhibit irreversible processes, including inelastic and nonrecoverable relaxation effects in the process zone surrounding the crack tip. Here, classical molecular dynamics simulations were used with a reactive forcefield to evaluate inelastic relaxation processes in five amorphous sodium silicate compositions. Overall, the 20% Na2O-SiO2(NS20) composition exhibited the most inelastic relaxation, followed by the 15% Na2O-SiO2(NS15) composition, the 25% Na2O-SiO2(NS25) composition, and finally the 10% (NS10) and 30% (NS30) Na2O-SiO2 compositions. Coordination analysis of the Na+ ions identified that during inelastic relaxation the Na+ ions were increasingly coordinated by nonbridging oxygens (NBOs) for the NS10 and NS15 compositions, which was supported by radial analysis of the O-Na-O bond angles surrounding the crack tip. Across the sodium silicate compositional range, two different inelastic relaxation mechanism were identified based on the amount of bridging oxygens (BOs) and NBOs in the Na+ ion coordination shell. At lower (NS10) and higher (NS30) sodium compositions, the entire structured relaxed toward the crack tip. In contrast at intermediate sodium concentrations (NS20) the Na+ ion migrates toward the crack tip separately from the network structure. By developing a fundamental understanding of how modified silica systems respond to static stress fields, we will be able to predict how varying amorphous silicate systems exhibit slow crack growth.

A generalized closed-form equation for the shaded collector fraction in solar arrays on rolling or undulating terrain is provided for single-axis tracking and fixed-tilt systems. The equation accounts for different rotation angles between the shaded and shading trackers, cross-axis slope between the two trackers, and offset between the collector plane and axis of rotation. The validity of the equation is demonstrated through comparison with numerical ray-tracing simulations and remaining minor sources of error are quantified. Additionally, a simple procedure to determine backtracking rotations for each row in an array installed on the rolling terrain (varying in the direction perpendicular to the tracker axes) is provided. The backtracking equation accounts for a desired shaded fraction (including complete shade avoidance) as well as an axis-collector offset. Test cases are provided to facilitate implementation of these equations.

Laccases from white-rot fungi catalyze lignin depolymerization, a critical first step to upgrading lignin to valuable biodiesel fuels and chemicals. In this study, a wildtype laccase from the basidiomycete Fomitiporia mediterranea (Fom_lac) and a variant engineered to have a carbohydrate-binding module (Fom_CBM) were studied for their ability to catalyze cleavage of β-O-4′ ether and C–C bonds in phenolic and non-phenolic lignin dimers using a nanostructure-initiator mass spectrometry-based assay. Fom_lac and Fom_CBM catalyze β-O-4′ ether and C–C bond breaking, with higher activity under acidic conditions (pH < 6). The potential of Fom_lac and Fom_CBM to enhance saccharification yields from untreated and ionic liquid pretreated pine was also investigated. Adding Fom_CBM to mixtures of cellulases and hemicellulases improved sugar yields by 140% on untreated pine and 32% on cholinium lysinate pretreated pine when compared to the inclusion of Fom_lac to the same mixtures. Adding either Fom_lac or Fom_CBM to mixtures of cellulases and hemicellulases effectively accelerates enzymatic hydrolysis, demonstrating its potential applications for lignocellulose valorization. We postulate that additional increases in sugar yields for the Fom_CBM enzyme mixtures were due to Fom_CBM being brought more proximal to lignin through binding to either cellulose or lignin itself.

During hypersonic flight, compressional and viscous heating of the air can form a plasma layer which encases the aircraft. If the boundary layer becomes turbulent, then the electron density fluctuations can effect a parasitic modulation in microwave signals transmitted through the plasma. We developed an approach for studying the interaction of microwave signals with a turbulent, hypersonic plasma layer. The approach affords a great deal of flexibility in both the plasma layer model and the antenna configuration. We then analyzed a situation in which microwaves, transmitted from a rectangular aperture antenna, propagate through a turbulent plasma layer to a distant receiver. We characterized the first- and second-order statistics of the computed parasitic modulation and quantified the depolarization of the signal. The amplitude fluctuations are lognormally distributed at low frequencies and Rice-distributed at high frequencies. Fluctuations in the copolarized phase and amplitude of the far-field signal are strongly anticorrelated. We used a multioutput Gaussian process (MOGP) to model these quantities. The efficacy of the MOGP model is demonstrated by recovering the time evolution of the copolarized phase given the copolarized amplitude and occasional measurements of the phase.

Physical networks formed by ionizable polymers with ionic clusters as crosslinks are controlled by coupled dynamics that transcend from ionic clusters through chain motion to macroscopic response. Here, the coupled dynamics, across length scales, from the ionic clusters to the networks in toluene swollen polystyrene sulfonate networks, were directly correlated, as the electrostatic environment of the physical crosslinks was altered. The multiscale insight is attained by coupling neutron spin echo measurements with molecular dynamics simulations, carried out to times typical of relaxation of polymers in solutions. The experimental dynamic structure factor is in outstanding agreement with the one calculated from computer simulations, as the networks are perturbed by elevating the temperature and changing the electrostatic environment. In toluene, the long-lived clusters remain stable over hundreds of ns across a broad temperature range, while the polymer network remains dynamic. Though the size of the clusters changes as the dielectric constant of the solvent is modified through the addition of ethanol, they remain stable but morph, enhancing the polymer chain dynamics.

The SNL Sierra Mechanics code suite is designed to enable simulation of complex multiphysicsscenarios. The code suite is composed of several specialized applications which can operate either instandalone mode or coupled with each other. Arpeggio is a supported utility that enables loose couplingof the various Sierra Mechanics applications by providing access to Framework services that facilitatethe coupling.

Junctions are discontinuities in flat grain boundaries that arise in all polycrystalline materials and are thought to play important roles in the response of a grain boundary network to thermal and mechanical loads. A key open question concerns the mechanisms by which solute segregation to junctions impacts properties of the grain boundary. Here, in this work, we investigate the influence of grain boundary facet junctions on solute embrittlement, and we present an analytical model that uses the hydrostatic stress field contributed by dislocations at multiple junctions to describe these effects. Specifically, we study junctions between {112} facets of various lengths in Au $\langle111\rangle$ Σ3 tilt grain boundaries. Copper and silver solutes are employed to determine if the effect of junctions on solute segregation and embrittlement is dependent on size relative to the host. Combined, atomistic simulation data and the analytical model show that Cu and Ag have opposite segregation responses to junctions due to the sign of the hydrostatic stress field induced by junctions. However, a positive shift in the embrittling potency is computed near junctions regardless of solute type or the stress state of the segregation site. Hence, for the conditions studied, junctions consistently shift the energetic landscape towards embrittlement.

Glycoboehmite (GB) materials are synthesized by a solvothermal reaction to form layered aluminum oxyhydroxide (boehmite) modified by intercalated butanediol molecules. These hybrid materials offer a platform to design materials with potentially novel sorption, wetting, and catalytic properties. Several synthetic methods have been used, resulting in different structural and spectroscopic properties, but atomistic detail is needed to determine the interlayer structure to explore the synthetic control of GB materials. Here, we use classical molecular dynamics (MD) simulations to compare the structural properties of GB interlayers containing chemisorbed butanediol molecules as a function of diol loading. Accompanying quantum (density functional theory, DFT) static calculations and MD simulations are used to validate the classical model and compute the infrared spectra of various models. Classical MD results reveal the existence of two unique interlayer environments at higher butanediol loading, corresponding to smaller (cross-linked) and expanded interlayers. DFT-computed infrared spectra reveal the sensitivity of the aluminol O-H stretch frequencies to the interlayer environment, consistent with the spectrum of the synthesized material. Insight from these simulations will aid in the characterization of the newly synthesized GB materials.

The paper deals with a new effective numerical technique on unfitted Cartesian meshes for simulations of heterogeneous elastic materials. Here, we develop the optimal local truncation error method (OLTEM) with 27- point stencils (similar to those for linear finite elements) for the 3-D time-independent elasticity equations with irregular interfaces. Only displacement unknowns at each internal Cartesian grid point are used. The interface conditions are added to the expression for the local truncation error and do not change the width of the stencils. The unknown stencil coefficients are calculated by the minimization of the local truncation error of the stencil equations and yield the optimal second order of accuracy for OLTEM with the 27-point stencils on unfitted Cartesian meshes. A new post-processing procedure for accurate stress calculations has been developed. Similar to basic computations it uses OLTEM with the 27-point stencils and the elasticity equations. The post-processing procedure can be easily extended to unstructured meshes and can be independently used with existing numerical techniques (e.g., with finite elements). Numerical experiments show that at an accuracy of 0.1% for stresses, OLTEM with the new post-processing procedure significantly (by 10⁵-10⁹ times) reduces the number of degrees of freedom compared to linear finite elements. OLTEM with the 27-point stencils yields even more accurate results than high-order finite elements with wider stencils.

Ionic assemblies, or clusters, determine the structure and dynamics of ionizable polymers and enable their many applications. Fundamental to attaining well-defined materials is controlling the balance between the van der Waals interactions that govern the backbone behavior and the forces that drive the formation of ionic clusters. Here, using small-angle neutron scattering and fully atomistic molecular dynamics simulations, the structure of a model ionomer, sulfonated polystyrene in toluene solutions, was investigated as the cluster cohesion was tweaked by the addition of ethanol. The static structure factor was measured by both techniques and correlated with the size of the ionic clusters as the polymer concentration was varied. The conjunction of SANS results and molecular insight from MD simulations enabled the determination of the structure of these inhomogeneous networks on multiple length scales. We find that across the entire concentration range studied, a network driven by the formation of ionic clusters was formed, where the size of the clusters drives the inhomogeneity of these systems. Tweaking the ionic clusters through the addition of ethanol impacts the packing of the sulfonated groups, their shape, and their size distribution, which, in turn, affects the structure of these networks.

Characterizing and identifying cells in multicellular in vitro models remain a substantial challenge. Here, we utilize hyperspectral confocal Raman microscopy and principal component analysis coupled with linear discriminant analysis to form a label-free, noninvasive approach for classifying bone cells and osteosarcoma cells. Through the development of a library of hyperspectral Raman images of the K7M2-wt osteosarcoma cell lines, 7F2 osteoblast cell lines, RAW 264.7 macrophage cell line, and osteoclasts induced from RAW 264.7 macrophages, we built a linear discriminant model capable of correctly identifying each of these cell types. The model was cross-validated using a k-fold cross validation scheme. The results show a minimum of 72% accuracy in predicting cell type. We also utilize the model to reconstruct the spectra of K7M2 and 7F2 to determine whether osteosarcoma cancer cells and normal osteoblasts have any prominent differences that can be captured by Raman. We find that the main differences between these two cell types are the prominence of the β-sheet protein secondary structure in K7M2 versus the α-helix protein secondary structure in 7F2. Additionally, differences in the CH2 deformation Raman feature highlight that the membrane lipid structure is different between these cells, which may affect the overall signaling and functional contrasts. Overall, we show that hyperspectral confocal Raman microscopy can serve as an effective tool for label-free, nondestructive cellular classification and that the spectral reconstructions can be used to gain deeper insight into the differences that drive different functional outcomes of different cells.

Superconducting topological systems formed by a strong 3D topological insulator (TI) in proximity to a conventional s-wave superconductor (SC) have been intensely studied, as they may host Majorana zero modes. However, there are limited experimental realizations of TI-SC systems in which robust superconducting pairing is induced on the surface states of the TI and a topological superconducting state is established. Here, we fabricate a TI-SC system by depositing, via a focused ion beam, tungsten (W) nanoscale clusters on the surface of TI Bi0.91Sb0.09. We find that the resulting heterostructure supports phase-slip lines (PSLs) that act as effective Josephson junctions (JJs). We probe the response of the system to microwave radiation. We find that for some ac frequencies, and powers, the resulting Shapiro steps’ structure of the voltage-current characteristic exhibits a missing first step and an unexpectedly wide second Shapiro step. The theoretical analysis of the measurements shows that the unusual Shapiro response arises from the interplay between a static JJ and a dynamic one and allows us to identify the conditions under which the missing first step can be attributed to the topological nature of the JJs formed by the PSLs. Our results suggest an approach to induce superconductivity in a TI, a route to realizing highly transparent topological JJs, and show how the response of superconducting systems to microwave radiation can be used to infer the dynamics of PSLs. Highly transparent topological junctions are promising candidates to realize vector field sensors with very high sensitivity. In addition, due to the nontrivial Berry phase of the TI’s surface states such junctions can be in a topological state which is ideal to create topologically protected qubits.

This presentation includes a look into Sandia critical experiments including the 7uPCX, BUCCX, and assembly design. This presentation touches on the completion of IER 305 with CED-3b, CED-4a, and CED-4b. Finally, there are preparations to perform IER 441 including new hardware, critical configurations, and next steps.

This presentation includes a IER 441 assembly overview and the difference between IER 305 AND IER 441 with central test region assembly and hex pitch. Next this presentation looks at IER 441 procurement issues and delays and new hardware. additionally conducted was a SPRF/CX: IER 441 hardware test fit (Success)). This presentation concludes with lessons learned and acknowledgements.

Testing is pivotal for early identification of disease and subsequent infection control. Pathogens’ nucleic acid sequence can change due to naturally-occurring genetic drift or intentional modification. Because of the reliance on molecular assays for human, animal, and plant disease diagnosis, we must understand how nucleotide mutations affect test accuracy. Primers designed against original lineages of a pathogen may be less efficient at detecting variants with genetic changes in priming regions. Here, we made single- and multi-point mutations in priming regions of a model SARS-CoV-2 template that was used as input for a loop-mediated isothermal amplification (LAMP) assay. We found that many of the modifications impacted assay sensitivity, amplification speed, or both. Further research exploring mutations at every position in each of the eight priming regions should be conducted to evaluate trends and determine generalizability.

This presentation titled "Updates on UO₂-BeO Experiment (IER 523)" covers experiment status, experiment motivation, CED-1 summary, current efforts (CED-2), and includes a concluding summary.

Wurtzite ferroelectrics are an emerging material class that expands the functionality and application space of wide bandgap semiconductors. Promising physical properties of binary wurtzite semiconductors include a large, reorientable spontaneous polarization, direct band gaps that span from the infrared to ultraviolet, large thermal conductivities and acoustic wave velocities, high mobility electron and hole channels, and low optical losses. The ability to reverse the polarization in ternary wurtzite semiconductors at room temperature enables memory and analog type functionality and quasi-phase matching in optical devices and boosts the ecosystem of wurtzite semiconductors, provided the appropriate combination of properties can be achieved for any given application. In this article, advances in the design, synthesis, and characterization of wurtzite ferroelectric materials and devices are discussed. Highlights include: the direct and quantitative observation of polarization reversal of ∼135 μC/cm2 charge in Al1−xBxN via electron microscopy, Al1−xBxN ferroelectric domain patterns poled down to 400 nm in width via scanning probe microscopy, and full polarization retention after over 1000 h of 200 °C baking and a 2× enhancement relative to ZnO in the nonlinear optical response of Zn1−xMgxO. The main tradeoffs, challenges, and opportunities in thin film deposition, heterostructure design and characterization, and device fabrication are overviewed.

In this study, x-ray imaging addresses many challenges with visible light imaging in extreme environments, where optical diagnostics such as digital image correlation (DIC) and particle image velocimetry (PIV) suffer biases from index of refraction changes and/or cannot penetrate visibly occluded objects. However, conservation of intensity—the fundamental principle of optical image correlation algorithms—may be violated if ancillary components in the X-ray path besides the surface or fluid of interest move during the test. The test series treated in this work sought to characterize the safe use of fiber-epoxy composites in aerospace and aviation industries during fire accident scenarios. Stereo X-ray DIC was employed to measure test unit deformation in an extreme thermal environment involving a visibly occluded test unit, incident surface heating to temperatures above 600°C, and flames and soot from combusting epoxy decomposition gasses. The objective of the current work is to evaluate two solutions to resolve the violation of conservation of intensity that resulted from both the test unit and the thermal chamber deforming during the test. The first solution recovered conservation of intensity by pre-processing the path-integrated X-ray images to isolate the DIC pattern of the test unit from the thermal chamber components. These images were then correlated with standard, optical DIC software. The second solution, called Path-Integrated Digital Image Correlation (PI-DIC), modified the fundamental matching criterion of DIC to account for multiple, independently-moving components contributing to the final image intensity. The PI-DIC algorithm was extended from a 2D framework to a stereo framework and implemented in a custom DIC software. Both solutions accurately measured the cylindrical shape of the undeformed test unit, recovering radii values of R = 76.20±0.12 mm compared to the theoretical radius of R_theor = 76.20 mm. During the test, the test unit bulged asymmetrically as decomposition gasses pressurized the interior and eventually burned in a localized jet. Both solutions measured the heterogeneous radius of this bulge, which reached a maximum of approximately R = 91 mm, with a discrepancy of 2–3% between the two solutions. Two solutions that resolve the violation of conservation of intensity for path-integrated X-ray images were developed. Both were successfully employed in a stereo X-ray DIC configuration to measure deformation of an aluminum-skinned, fiber-epoxy composite test unit in a fire accident scenario.

Over the last few years, crystalline topology has been used in photonic crystals to realize edge- and corner-localized states that enhance light-matter interactions for potential device applications. However, the band-theoretic approaches currently used to classify bulk topological crystalline phases cannot predict the existence, localization, or spectral isolation of any resulting boundary-localized modes. While interfaces between materials in different crystalline phases must have topological states at some energy, these states need not appear within the band gap, and thus may not be useful for applications. Here, we derive a class of local markers for identifying material topology due to crystalline symmetries, as well as a corresponding measure of topological protection. As our real-space-based approach is inherently local, it immediately reveals the existence and robustness of topological boundary-localized states, yielding a predictive framework for designing topological crystalline heterostructures. In conclusion, beyond enabling the optimization of device geometries, we anticipate that our framework will also provide a route forward to deriving local markers for other classes of topology that are reliant upon spatial symmetries.

Surrogate fuels that reproduce the characteristics of full-boiling range fuels are key tools to enable numerical simulations of fuel-related processes and ensure reproducibility of experiments by eliminating batch-to-batch variability. Within the PACE initiative, a surrogate fuel for regular-grade E10 (10%vol ethanol) gasoline representative of a U.S. market gasoline, termed PACE-20, was developed and adopted as baseline fuel for the consortium. Although extensive testing demonstrated that PACE-20 replicates the properties and combustion behavior of the full-boiling range gasoline, several concerns arose regarding the purity level required for the species that compose PACE-20. This is particularly important for cyclo-pentane, since commercial-grade cyclo-pentane typically shows 60%-85% purity. In the present work, the effects of the purity level of cyclo-pentane on the properties and combustion characteristics of PACE-20 were studied. Chemical kinetic simulations were performed to predict the effects of cyclo-pentane impurities on the properties, octane rating, and autoignition reactivity under homogeneous charge compression-ignition conditions of PACE-20. From the numerical results, cyclo-pentane with 85% purity or higher is required to reasonably match both the research octane number and motor octane number of the target gasoline. Finally, homogeneous charge compression-ignition engine simulations show that impurities have only a modest effect on reactivity at naturally aspirated conditions, but cyclo-pentane purity is critical to properly replicate the pressure dependency of the reactivity.

Motivated by reducing errors in the energy budget related to enthalpy fluxes within the Energy Exascale Earth System Model (E3SM), we study several physics-dynamics coupling approaches. Using idealized physics, a moist rising bubble test case, and the E3SM's nonhydrostatic dynamical core, we consider unapproximated and approximated thermodynamics applied at constant pressure or constant volume. With the standard dynamics and physics time-split implementation, we describe how the constant-pressure and constant-volume approaches use different mechanisms to transform physics tendencies into dynamical motion and show that only the constant-volume approach is consistent with the underlying equations. Using time step convergence studies, we show that the two approaches both converge but to slightly different solutions. We reproduce the large inconsistencies between the energy flux internal to the model and the energy flux of precipitation when using approximate thermodynamics, which can only be removed by considering variable latent heats, both when computing the latent heating from phase change and when applying this heating to update the temperature. Finally, we show that in the nonhydrostatic case, for physics applied at constant pressure, the general relation that enthalpy is locally conserved no longer holds. In this case, the conserved quantity is enthalpy plus an additional term proportional to the difference between hydrostatic pressure and full pressure.

Pyrimidine has two in-plane CH(δ+)/N̈(δ−)/CH(δ+) binding sites that are complementary to the (δ−/2δ+/δ−) quadrupole moment of CO2. We recorded broadband microwave spectra over the 7.5-17.5 GHz range for pyrimidine-(CO2)n with n = 1 and 2 formed in a supersonic expansion. Based on fits of the rotational transitions, including nuclear hyperfine splitting due to the two 14N nuclei, we have assigned 313 hyperfine components across 105 rotational transitions for the n = 1 complex and 208 hyperfine components across 105 rotational transitions for the n = 2 complex. The pyrimidine-CO2 complex is planar, with CO2 occupying one of the quadrupolar binding sites, forming a structure in which the CO2 is stabilized in the plane by interactions with the C-H hydrogens adjacent to the nitrogen atom. This structure is closely analogous to that of the pyridine-CO2 complex studied previously by ( Doran, J. L. J. Mol. Struct. 2012, 1019, 191-195 ). The fit to the n = 2 cluster gives rotational constants consistent with a planar cluster of C2v symmetry in which the second CO2 molecule binds in the second quadrupolar binding pocket on the opposite side of the ring. The calculated total binding energy in pyrimidine-CO2 is −13.7 kJ mol-1, including corrections for basis set superposition error and zero-point energy, at the CCSD(T)/ 6-311++G(3df,2p) level, while that in pyrimidine-(CO2)2 is almost exactly double that size, indicating little interaction between the two CO2 molecules in the two binding sites. The enthalpy, entropy, and free energy of binding are also calculated at 300 K within the harmonic oscillator/rigid-rotor model. This model is shown to lack quantitative accuracy when it is applied to the formation of weakly bound complexes.

Compressed hydrogen stored at cryogenic temperatures has a much higher density than room-temperature storage, which enables large-scale hydrogen storage and transport. An understanding of the release of cryogenic hydrogen from pressurized vessels is needed to evaluate the risk and safety concerns with the use of this fuel. The present work extends the analysis of previous experimental studies that measured the gas concentrations of cryo-compressed hydrogen jets and methane jets using a laser Raman scattering diagnostic system. Since the Raman signals are very small, a denoising algorithm was applied to significantly reduce the noise to enable statistical analysis of the data. The transient features of the turbulent jets were characterized by their concentration intermittencies and probability density functions (PDFs). A two-part PDF was developed to predict the bimodal features of the jet concentration distributions. Then, the flammability factors of the cryogenic jets were calculated based on the intermittency and the PDF.

Experiments have shown that ducted fuel injection (DFI) effectively reduces soot emissions from direct-injection diesel engines. Although many computational studies have evaluated DFI's spray development and soot reduction mechanisms in constant volume chambers, only limited computational work on internal combustion engines exists. The DFI duct assembly changes the engine's in-cylinder flow, spray, and combustion development. Therefore, current production engine designs might not be optimal for achieving the best engine performance with DFI. This work conducted an extensive numerical study to evaluate how parameter changes affect DFI performance. The parameters include swirl ratio, piston geometry, compression ratio (CR), number of injector orifices, split injection strategy, and exhaust gas recirculation (EGR) in a heavy-duty diesel engine utilizing DFI. The combustion and soot emission data from the Sandia compression ignition optical research engine were used for model validation. Simulations showed that an increased swirl ratio resulted in more intense jet flame-piston interaction, slowing down the combustion heat release during the late combustion stage and leading to lower indicated thermal efficiency (ITE) due to higher exhaust losses. A piston-bowl design with a reentrant inner piston edge yielded the highest thermal efficiency, due to the reduced cylinder head heat transfer loss. Additional injector orifices led to higher efficiency owing to a more advanced combustion phasing. Nevertheless, the maximum pressure rise rate (MPRR) and oxides of nitrogen (NOx) emissions also increased with the number of injector orifices due to more rapid heat release and higher combustion temperature. Implementation of a split injection strategy combined with a higher EGR rate effectively inhibited the excessive MPRR and NOx formation. In general, the study concluded that DFI is not sensitive to most parameter changes but will benefit from future parameter optimization.

Thermal spray deposition is an inherently stochastic manufacturing process used for generating thick coatings of metals, ceramics and composites. The generated coatings exhibit hierarchically complex internal structures that affect the overall properties of the coating. The deposition process can be adequately simulated using rules-based process simulations. Nevertheless, in order for the simulation to accurately model particle spreading upon deposition, a set of predefined rules and parameters need to be calibrated to the specific material and processing conditions of interest. The calibration process is not trivial given the fact that many parameters do not correspond directly to experimentally measurable quantities. This work presents a protocol that automatically calibrates the parameters and rules of a given simulation in order to generate the synthetic microstructures with the closest statistics to an experimentally generated coating. Specifically, this work developed a protocol for tantalum coatings prepared using air plasma spray. The protocol starts by quantifying the internal structure using 2-point statistics and then representing it in a low-dimensional space using Principal Component Analysis. Subsequently, our protocol leverages Bayesian optimization to determine the parameters that yield the minimum distance between synthetic microstructure and the experimental coating in the low-dimensional space.

Femtosecond laser electronic excitation tagging (FLEET) velocimetry is an important diagnostic technique for seedless velocimetry measurements particularly in supersonic and hypersonic flows. Typical FLEET measurements feature a single laser line and camera system to achieve one-component velocimetry along a line, although some multiple-spot and multiple-component configurations have been demonstrated. In this work, tomographic imaging is used to track the three-dimensional location of many FLEET spots. A quadscope is used to combine four unique views onto a single high-speed image intensifier and camera. Tomographic reconstructions of the FLEET emission are analyzed for three-component velocimetry from multiple FLEET spots. Glass wedges are used to create many (nine) closely spaced FLEET spots with less than 10% transmission losses. These developments lead to a significant improvement in the dimensionality and spatial coverage of a FLEET instrument with some increases in experimental complexity and data processing. Multiple-point three-component FLEET velocimetry is demonstrated in an underexpanded jet.

In many applications, physical domains are geometrically complex making it challenging to perform coarse-scale approximation. A defeaturing process is often used to simplify the domain in preparation for approximation and analysis at the coarse scale. Herein, a methodology is presented for constructing a coarse-scale reproducing basis on geometrically complex domains given an initial fine-scale mesh of the fully featured domain. The initial fine-scale mesh can be of poor quality and extremely refined. The construction of the basis functions begins with a coarse-scale covering of the domain and generation of weighting functions with local support. Manifold geodesics are used to define distances within the local support for general applicability to non-convex domains. Conventional moving least squares is used to construct the coarse-scale reproducing basis. Applications in quasi-interpolation and linear elasticity are presented.

Gas molecule clustering within nanopores holds significance in the fields of nanofluidics, biology, gas adsorption/desorption, and geological gas storage. However, the intricate roles of nanoconfinement and surface chemistry that govern the formation of gas clusters remain inadequately explored. In this study, through free energy calculation in molecular simulations, we systematically compared the tendencies of H2 and CO2 molecules to aggregate within hydrated hydrophobic pyrophyllite and hydrophilic gibbsite nanopores. The results indicate that nanoconfinement enhances gas dimer formation in the nanopores, irrespective of surface chemistry. However, surface hydrophilicity prohibits the formation of gas clusters larger than dimers, while large gas clusters form easily in hydrophobic nanopores. Despite H2 and CO2 both being non-polar, the larger quadrupole moment of CO2 leads to a stronger preference for dimer/cluster formation compared to H2. Our results also indicate that gases prefer to enter the nanopores as individual molecules, but exit the nanopores as dimers/clusters. This investigation provides a mechanistic understanding of gas cluster formation within nanopores, which is relevant to various applications, including geological gas storage.

We present a high-level architecture for how artificial intelligences might advance and accumulate scientific and technological knowledge, inspired by emerging perspectives on how human intelligences advance and accumulate such knowledge. Agents advance knowledge by exercising a technoscientific method—an interacting combination of scientific and engineering methods. The technoscientific method maximizes a quantity we call “useful learning” via more-creative implausible utility (including the “aha!” moments of discovery), as well as via less-creative plausible utility. Society accumulates the knowledge advanced by agents so that other agents can incorporate and build on to make further advances. The proposed architecture is challenging but potentially complete: its execution might in principle enable artificial intelligences to advance and accumulate an equivalent of the full range of human scientific and technological knowledge.

In lithium-metal batteries, grains of lithium can become electrically isolated from the anode, lowering battery performance. We state that experiments reveal that rest periods after battery discharge might help to solve this problem.

Replacement of conventional petroleum fuels with renewable fuels reduces net emissions of carbon and greenhouse gases, and affords opportunities for increased domestic energy security. Here, we present alkyl dialkoxyalkanoates (or DAOAs) as a family of synthetic diesel and marine fuel candidates that feature ester and ether functionality. These compounds employ pyruvic acid and fusel alcohols as precursors, which are widely available as metabolic intermediates at high titer and yield. DAOA synthesis proceeds in high yield using a simple, mild chemical transformation performed under air that employs bioderived and/or easily recovered reagents and solvent. The scalability of the synthetic protocol was proven in continuous flow with in situ azeotropic water removal, yielding 375 g of isolated product. Chemical stability of DAOAs against aqueous 0.01 M H2SO4 and accelerated oxidative conditions is demonstrated. The isolated DAOAs were shown to meet or exceed widely accepted technical criteria for sustainable diesel fuels. In particular, butyl 2,2-dibutoxypropanoate (DAOA-2) has indicated cetane number 64, yield soot index 256 YSI per kg, lower heating value 30.9 MJ kg−1 and cloud point < −60 °C and compares favorably to corresponding values for renewable diesel, biodiesel and petroleum diesel.

Quantifying the radioactive sources present in gamma spectra is an ever-present and growing national security mission and a time-consuming process for human analysts. While machine learning models exist that are trained to estimate radioisotope proportions in gamma spectra, few address the eventual need to provide explanatory outputs beyond the estimation task. In this work, we develop two machine learning models for a NaI detector measurements: one to perform the estimation task, and the other to characterize the first model’s ability to provide reasonable estimates. To ensure the first model exhibits a behavior that can be characterized by the second model, the first model is trained using a custom, semi-supervised loss function which constrains proportion estimates to be explainable in terms of a spectral reconstruction. The second auxiliary model is an out-of-distribution detection function (a type of meta-model) leveraging the proportion estimates of the first model to identify when a spectrum is sufficiently unique from the training domain and thus is out-of-scope for the model. In demonstrating the efficacy of this approach, we encourage the use of meta-models to better explain ML outputs used in radiation detection and increase trust.

Living document updating the 5.16 version, SAND2023-09919 of September 2023.

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

As machine learning models for radioisotope quantification become more powerful, likewise the need for high-quality synthetic training data grows as well. For problem spaces that involve estimating the relative isotopic proportions of various sources in gamma spectra it is necessary to generate training data that accurately represents the variance of proportions encountered. In this report, we aim to provide guidance on how to target a desired variance of proportions which are randomly when using the PyRIID Seed Mixer, which samples from a Dirichlet distribution. We provide a method for properly parameterizing the Dirichlet distribution in order to maintain a constant variance across an arbitrary number of dimensions, where each dimension represents a distinct source template being mixed. We demonstrate that our method successfully parameterizes the Dirichlet distribution to target a specific variance of proportions, provided that several conditions are met. This allows us to follow a principled technique for controlling how random mixture proportions are generated which are then used downstream in the synthesis process to produce the final, noisy gamma spectra.

NasGen provides a path for migration of structural models from NASTRAN bulk data format (BDF) into both an Exodus mesh file and an ASCII input file for Sierra Structural Dynamics (Salinas) and Solid Mechanics (Presto).

Tests from the Sierra Structural Dynamics verification test suite are reviewed. Each is run nightly and the results of the test checked versus the correct analytic result. For each of the tests presented in this document the test setup, derivation of the analytic solution, and comparison of the Sierra code results to the analytic solution is provided. This document can be used to confirm that a given code capability is verified or referenced as a compilation of example problems.