Publications

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Testing of a compact Bremsstrahlung diode at the High Energy Radiation Megavolt Electron Source III (HERMES-III) was performed at Sandia National Laboratories in November, 2023. The compact diode described here is the first prototype diode in a campaign to optimize a Bremsstrahlung diode in terms of size and dose production. The goal was to test the diode at 13MV, and the experiment realized between 10-12MV at the diode. Modeling and simulation of this geometry was performed after the test, shedding insight into several phenomena seen by experimental diagnostics. Modifications to the diode and experiment are proposed for future experiments to help explain the phenomena and approach a better final design.

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.

CFD (Computational Fluid Dynamic) simulation of aerosol-laden natural convective flow and particle deposition in a spent fuel storage canister with 37 assemblies is currently computationally prohibitive. PWR (Pressurized Water Reactor) assemblies have up to 289 pins or tubes with several spacer grids to align the pins. Spacer grids with mixing vanes induce swirling during operation to increase heat transfer. Each spacer grid contains hundreds of small structures such as retaining clips, channel walls, and openings. The largest canisters store 37 PWR assemblies thus, there are numerous pins, tubes, and spacer grids for which the flow region between and around these structures needs to be determined along with the movement and deposition of aerosol particles. Because of the complicated geometry, modeling the intricate flow even for just one assembly is currently impractical. Nonetheless, we are developing techniques for a practical model to assess the natural aerosol particle deposition process in a canister in the event that a release occurs from one or more fuel pins. In the previous work it was demonstrated that CFD can model the flow through a PWR spacer grid with mixing vanes, including particle deposition, in a reasonable amount of time on a personal computer. In this work, the analysis is extended to include the bypass region between an assembly and the canister basket walls. It is shown that the flow velocity in the bypass region is about three times that of the interstitial region between the pins. The lengths before and after the spacer grid are also extended to determine when the flow becomes fully developed. In addition, the approach of computationally “stitching together” segments of an assembly is demonstrated with the plan to ultimately model a full assembly. The fraction of particles that are deposited in a segment with a spacer grid is determined as a function of particle size and flow velocity.

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.

Bonding diamond to the back side of gallium nitride (GaN) electronics has been shown to improve thermal management in lateral devices; however, engineering challenges remain with the bonding process and characterizing the bond quality for vertical device architectures. Here, integration of these two materials is achieved by room-temperature compression bonding centimeter-scale GaN and a diamond die via an intermetallic bonding layer of Ti/Au. Recent attempts at GaN/diamond bonding have utilized a modified surface activation bonding (SAB) method, which requires Ar fast atom bombardment immediately followed by bonding within the same tool under ultrahigh vacuum (UHV) conditions. The method presented here does not require a dedicated SAB tool yet still achieves bonding via a room-temperature metal-metal compression process. Imaging of the buried interface and the total bonding area is achieved via transmission electron microscopy (TEM) and confocal acoustic scanning microscopy (C-SAM), respectively. The thermal transport quality of the bond is extracted from spatially resolved frequency-domain thermoreflectance (FDTR) with the bonded areas boasting a thermal boundary conductance of >100 MW/m2·K. Additionally, Raman maps of GaN near the GaN-diamond interface reveal a low level of compressive stress, <80 MPa, in well-bonded regions. FDTR and Raman were coutilized to map these buried interfaces and revealed some poor thermally bonded areas bordered by high-stress regions, highlighting the importance of spatial sampling for a complete picture of bond quality. Overall, this work demonstrates a novel method for thermal management in vertical GaN devices that maintains low intrinsic stresses while boasting high thermal boundary conductances.

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.

Coarse-grained molecular dynamics simulations are used to study the diffusion of thin nanorods in entangled polymer melts for varying nanorod length and roughness. While prior studies observed a nanorod parallel diffusion constant scaling inversely with rod length D_∥ ~ l^–1, here, we show that this scaling is not universal and depends sensitively on the nanorod surface roughness. We observe D_∥ ~ l^–k, where k < 1 and decreases with decreasing surface roughness. The weaker scaling is driven by the non-Gaussian diffusion of nanorods due to the emergence of an intermittent hopping process that becomes more pronounced with decreasing roughness at the monomer scale. Analysis shows that the mean hop size grows for smoother rods but shows little to no variation with rod length. The mean hopping frequency shows no dependence on either rod length or roughness, suggesting it originates from the polymer melt environment. Further, our results show that the small-scale features of the nanorod surface strongly influence the large-scale and long-time transport of nanorods in polymer matrices, creating new material design opportunities for precisely engineered nanocomposites.

The First Tri-Lab Workshop on Formal Verification was held in Santa Fe, New Mexico, on December 5th, 2023. This workshop gathered staff from Sandia, Los Alamos, and Lawrence Livermore National Laboratories and NASA’s Jet Propulsion Laboratory. This report summarizes and expands on the presentations given and discussion had at this workshop. In this report, we describe the current capabilities and research needs related to formal methods at the NNSA labs. In particular, we identify medium-term and long-term research gaps in programming languages, formalization efforts of complex systems, embedded systems verification, hardware verification, cybersecurity, formal methods usability, workflows, numerical methods, the use of formal methods for artificial intelligence (and its converse, artificial intelligence for formal methods), and collaboration opportunities and considerations on these topics. We conclude with a small number of exemplar research problems related to these topics.

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.

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.

Objective. Standard signal processing approaches for scintillation detectors in positron emission tomography (PET) derive accurate estimates for 511 keV photon time of interaction and energy imparted to the detection media from aggregate characteristics of electronic pulse shapes. The ultimate realization of a scintillation detector for PET is one that provides a unique timestamp and position for each detected scintillation photon. Detectors with these capabilities enable advanced concepts for three-dimensional (3D) position and time of interaction estimation with methods that exploit the spatiotemporal arrival time kinetics of individual scintillation photons. Approach. In this work, we show that taking into consideration the temporal photon emission density of a scintillator, the channel density of an analog silicon photomultiplier (SiPM) array, and employing fast electronic readout with digital signal processing, a detector that counts and timestamps scintillation photons can be realized. To demonstrate this approach, a prototype detector was constructed, comprising multichannel electronic readout for a bismuth germanate (BGO) scintillator coupled to an SiPM array. Main Results. In proof-of-concept measurements with this detector, we were able to count and provide unique timestamps for 66% of all optical photons, where the remaining 34% (two-or-more-photon pulses) are also independently counted, but each photon bunch shares a common timestamp. We show this detector concept can implement 3D positioning of 511 keV photon interactions and thereby enable corrections for time of interaction estimators. The detector achieved 17.6% energy resolution at 511 keV and 237 ± 10 ps full-width-at-half-maximum coincidence time resolution (CTR) (fast spectral component) versus a reference detector. We outline the methodology, readout, and approach for achieving this detector capability in first-ever, proof-of-concept measurements for scintillation photon counting detector with analog silicon photomultipliers. Significance. The presented detector concept is a promising design for large area, high sensitivity TOF-PET detector modules that can implement advanced event positioning and time of interaction estimators, which could push state-of-the-art performance.

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.

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.

This report summarizes the proceedings of the 13^th US/German Workshop on Salt Repository Research, Design, and Operation hosted by Sandia National Laboratories on June 20-23, 2023, in Santa Fe, New Mexico, USA. Over 60 participants attended, representing Germany, United States, the Netherlands, Australia, and the United Kingdom, along with the IAEA. The purpose of the US/German Workshop is to foster in-person collaboration and dialogue amongst salt repository researchers and nuclear waste disposal implementers across international organizations. The workshop included five sessions of topical presentations and two breakout sessions to promote additional discussion on compelling topics.

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.

The Example Problems Manual supplements the User's Manual and the Theory Manual. The goal of the Example Problems Manual is to reduce learning time for complex end to end analyses. These documents are intended to be used together. See the User's Manual for a complete list of the options for a solution case. All the examples are part of the Sierra/SD test suite. Each runs as is.

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.

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.

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.

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.

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.

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.

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.

Accurate and efficient constitutive modeling remains a cornerstone issue for solid mechanics analysis. Over the years, the LAMÉ advanced material model library has grown to address this challenge by implementing models capable of describing material systems spanning soft polymers to stiff ceramics including both isotropic and anisotropic responses. Inelastic behaviors including (visco)plasticity, damage, and fracture have all incorporated for use in various analyses. This multitude of options and flexibility, however, comes at the cost of many capabilities, features, and responses and the ensuing complexity in the resulting implementation. Therefore, to enhance confidence and enable the utilization of the LAMÉ library in application, this effort seeks to document and verify the various models in the LAMÉ library. Specifically, the broader strategy, organization, and interface of the library itself is first presented. The physical theory, numerical implementation, and user guide for a large set of models is then discussed. Importantly, a number of verification tests are performed with each model to not only have confidence in the model itself but also highlight some important response characteristics and features that may be of interest to end-users. Finally, in looking ahead to the future, approaches to add material models to this library and further expand the capabilities are presented.

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.

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.

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.

Presented in this document is a small portion of the tests that exist in the Sierra/SolidMechanics (Sierra/SM) verification test suite. Most of these tests are run nightly with the Sierra/SM code suite, and the results of the test are checked versus the correct analytical result. For each of the tests presented in this document, the test setup, a description of the analytic solution, and comparison of the Sierra/SM code results to the analytic solution is provided. Mesh convergence is also checked on a nightly basis for several of these tests. This document can be used to confirm that a given code capability is verified or referenced as a compilation of example problems. Additional example problems are provided in the Sierra/SM Example Problems Manual. Note, many other verification tests exist in the Sierra/SM test suite, but have not yet been included in this manual.

Sierra/SD provides a massively parallel implementation of structural dynamics finite element analysis, required for high-fidelity, validated models used in modal, vibration, static and shock analysis of weapons systems. This document provides a user’s guide to the input for Sierra/SD. Details of input specifications for the different solution types, output options, element types and parameters are included. The appendices contain detailed examples, and instructions for running the software on parallel platforms.

Presented in this document are tests that exist in the Sierra/SolidMechanics example problem suite, which is a subset of the Sierra/SM regression and performance test suite. These examples showcase common and advanced code capabilities. A wide variety of other regression and verification tests exist in the Sierra/SM test suite that are not included in this manual.

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.