Publications

Internal LDRD virtual poster session

Hydrogen powered locomotives are being explored to reduce emissions in rail applications. The risks of operations like refueling should be understood to ensure safe environments for workers and members of the public. Sensitivity analyses were conducted using HyRAM+ to identify major drivers of risk and compare effects of system parameters on individual risk. The consequences of jet fires from full-bore leaks dominated the risk, compared to explosions or smaller leaks. Pipe size, leak detection capability, and leak frequencies of system components greatly affected risk while overpressure modeling parameters and ambient conditions had little effect. The effects of personal protective equipment (PPE) materials on individual risk were quantified by reducing the individual’s exposure time or absorbed thermal dose. PPE only showed a risk reduction in low-risk cases. This study highlighted target areas for risk mitigation, including leak detection equipment and component maintenance, and indicated that the minimal effects of other parameters on risk may not justify prescriptive requirements for refueling operators.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Poster for CCN

The formation of magnesium chloride-hydroxide salts (magnesium hydroxychlorides) has implications for many geochemical processes and technical applications. For this reason, a thermodynamic database for evaluating the Mg(OH)2–MgCl2–H2O ternary system from 0 °C–120 °C has been developed based on extensive experimental solubility data. Internally consistent sets of standard thermodynamic parameters (ΔGf°, ΔHf°, S°, and CP) were derived for several solid phases: 3 Mg(OH)2:MgCl2:8H2O, 9 Mg(OH)2:MgCl2:4H2O, 2 Mg(OH)2:MgCl2:4H2O, 2 Mg(OH)2:MgCl2: 2H2O(s), brucite (Mg(OH)2), bischofite (MgCl2:6H2O), and MgCl2:4H2O. First, estimated values for the thermodynamic parameters were derived using a component addition method. These parameters were combined with standard thermodynamic data for Mg2+(aq) consistent with CODATA (Cox et al., 1989) to generate temperature-dependent Gibbs energies for the dissolution reactions of the solid phases. These data, in combination with values for MgOH+(aq) updated to be consistent with Mg2+-CODATA, were used to compute equilibrium constants and incorporated into a Pitzer thermodynamic database for concentrated electrolyte solutions. Phase solubility diagrams were constructed as a function of temperature and magnesium chloride concentration for comparisons with available experimental data. To improve the fits to the experimental data, reaction equilibrium constants for the Mg-bearing mineral phases, the binary Pitzer parameters for the MgOH+ — Cl− interaction, and the temperature-dependent coefficients for those Pitzer parameters were constrained by experimental phase boundaries and to match phase solubilities. These parameter adjustments resulted in an updated set of standard thermodynamic data and associated temperature-dependent functions. The resulting database has direct applications to investigations of magnesia cement formation and leaching, chemical barrier interactions related to disposition of heat-generating nuclear waste, and evaluation of magnesium-rich salt and brine stabilities at elevated temperatures.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Depleted uranium hexafluoride (UF₆), a stockpiled byproduct of the nuclear fuel cycle, reacts readily with atmospheric humidity, but the gas-phase reaction mechanism and associated chemical kinetics are poorly understood. During the performance period we undertook development of a state-of-the-art ab initio gas-phase chemical kinetics simulation workflow to model the hydrolysis of uranium hexafluroride (UF₆). In doing so, we addressed several outstanding issues in the theoretical treatment of uranium-containing systems. At the outset it was unclear how to generate accurate estimates of kinetic and thermodynamic data for U-containing chemical reactions. Generation of such data has been made routine. Prior to our work, the literature associated with UF₆ hydrolysis were disparate and inaccurate. This body of work provides a modern and comprehensive theoretical assessment of the reaction mechanism, molecular clustering towards deposition, and chemical kinetics. New methodological implementations and software integrations resulting from this work are also highlighted. As much as possible, our predictions were validated against experimental data including particle morphologies, vibrational spectroscopy, atomization enthalpies, and kinetic rate constants. Nevertheless, we were unable to reconcile kinetic measurements with high-accuracy simulations.

This document provides a technical description of algorithms used in the OpenCSP code repository that are common to deflectometry. This document is meant to be a technical reference that describes the internal algorithms of OpenCSP. An example program that uses these algorithms is SOFAST 2.0 [1].

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Poster presentation at DOE OE Peer Review 2024 in Bellevue Washington

Abstract not provided.

Abstract not provided.

The slides will be presented at a DOE CESER Peer Review for the Risk Management Tools and Technology (RMT) that provide an overview of the ESSec project. The slides discuss the containerization and orchestration tools for energy storage systems and our results along with an overview of the demonstration performed for this project. The peer review is scheduled for August 27-29.

Abstract not provided.

Abstract not provided.

Concentrating solar mirrors need to have correct geometry to accurately focus sunlight onto a receiver. SOFAST is a tool created at Sandia National Laboratories which measures a mirror’s surface normal at a dense sampling of points, enabling evaluation of a given mirror’s expected performance. SOFAST supports general mirrors, including concentrating solar heliostats, troughs, and dishes. The examples presented in this report predominantly focus on heliostats.

Abstract not provided.

Abstract not provided.

Abstract not provided.

A series of tension and compression tests were performed on typical transfer tubes made of 304L corrosion-resistant steel. The results were compiled and converted to true stress and true plastic strain for comparison to numerical models. The parameters of a Johnson-Cook constitutive model and fracture criterion were calibrated to mimic the test data. Comparisons were made between the converted data and the empirical formulas in the Johnson-Cook model. Simulations were performed to evaluate the implementation of the formula in finite element analysis, and the parameters were adjusted to improve the correlation.

Abstract not provided.

Abstract not provided.

As hydrogen storage facilities increase in size and capacity, it may be necessary to evaluate codes and standards that regulate hydrogen, especially in relation to allowable aggregate quantities. Existing regulations for other substances with comparable hazards such as oxygen, natural gas, and petroleum were reviewed; most fuels do not have aggregate quantity limits despite sometimes having individual tank capacity restrictions or limits for certain non-industrial, indoor, or small-scale applications. This precedent suggests that specifying a hydrogen quantity limit may not be necessary. Several possible methods for identifying an appropriate limit are presented to illustrate different approaches for a limit basis. These methods are based on overall risk, or the specific consequences inflicted by a jet fire or an explosion on people and infrastructure. However, these example metrics tend to require assumptions about potential leak sizes, suggesting that a general aggregate quantity limit would be difficult to justify without more system-specific requirements.

Abstract not provided.

The rise of grid modernization has been prompted by the escalating demand for power, the deteriorating state of infrastructure, and the growing concern regarding the reliability of electric utilities. The smart grid encompasses recent advancements in electronics, technology, telecommunications, and computer capabilities. Smart grid telecommunication frameworks provide bidirectional communication to facilitate grid operations. Software-defined networking (SDN) is a proposed approach for monitoring and regulating telecommunication networks, which allows for enhanced visibility, control, and security in smart grid systems. Nevertheless, the integration of telecommunications infrastructure exposes smart grid networks to potential cyberattacks. Unauthorized individuals may exploit unauthorized access to intercept communications, introduce fabricated data into system measurements, overwhelm communication channels with false data packets, or attack centralized controllers to disable network control. An ongoing, thorough examination of cyber attacks and protection strategies for smart grid networks is essential due to the ever-changing nature of these threats. Previous surveys on smart grid security lack modern methodologies and, to the best of our knowledge, most, if not all, focus on only one sort of attack or protection. This survey examines the most recent security techniques, simultaneous multi-pronged cyber attacks, and defense utilities in order to address the challenges of future SDN smart grid research. The objective is to identify future research requirements, describe the existing security challenges, and highlight emerging threats and their potential impact on the deployment of software-defined smart grid (SD-SG).

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.

Poster presentation at DOE OE Peer Review 2024 in Bellevue Washington

Abstract not provided.

Abstract not provided.

Naturally occurring uranium complicates monitoring for occupational exposures. There are several retroactive methods that can be used to monitor for occupational exposures, with benefits and drawbacks to each. Analysis of uranium in urine by mass spectrometry and alpha spectrometry is compared, and methods of determining an occupational exposure are presented. Furthermore, the minimum detectable concentrations from each analysis and a method for intake determination based on the analytical results are compared for various solubility types and mixtures. Mass spectrometry with radiochemical separation was found to be the most sensitive analysis for detecting occupational exposures to anthropogenic mixtures based on minimum detectable doses calculated from the proposed method for intake determination.

Nuclear power plant (NPP) risk assessment is broadly separated into disciplines of nuclear safety, security, and safeguards. Different analysis methods and computer models have been constructed to analyze each of these as separate disciplines. However, due to the complexity of NPP systems, there are risks that can span all these disciplines and require consideration of safety-security (2S) interactions which allows a more complete understanding of the relationship among these risks. A novel leading simulator/trailing simulator (LS/TS) method is introduced to integrate multiple generic safety and security computer models into a single, holistic 2S analysis. A case study is performed using this novel method to determine its effectiveness. The case study shows that the LS/TS method avoided introducing errors in simulation, compared to the same scenario performed without the LS/TS method. A second case study is then used to illustrate an integrated 2S analysis which shows that different levels of damage to vital equipment from sabotage at a NPP can affect accident evolution by several hours.

Poster covering a section of Thermogravimetric Analysis using a Quartz Crystal Microbalance, highlighting drop casting as a method.

This report summarizes the international collaborations conducted by Sandia funded by the US Department of Energy Office (DOE) of Nuclear Energy (DOE-NE) Spent Fuel and Waste Science & Technology (SFWST) as part of the Sandia National Laboratories Salt R&D and Salt International work packages. This report satisfies the level-three milestone M3SF-24SN010303063. Several stand-alone sections make up this summary report, each completed by the participants. The sections discuss granular salt reconsolidation (KOMPASS/MEASURES), engineered barriers (RANGERS), numerical model comparison (DECOVALEX), an NEA Salt Club working group on the development of scenarios as part of the performance assessment development process, and progress on seal percolation experiments with German colleagues at Gesellschaft für Anlagen- und Reaktorsicherheit (GRS). Finally, we summarize events related to the US/German Workshop on Repository Research, Design and Operations.

The objective of the project is to demonstrate that critical control functions can be implemented using little resources in modern microelectronics. A finite state machine (Figure 1) is implemented onto a field programmable gate array (FPGA). The functionality of the system is demonstrated by sending binary instructions to the controller. The controller transmits patterns through an LED, controls an electromechanical device, and uses pulse-width modulation (PWM) for radar functions.

This study conducts a comparative analysis, using non-equilibrium Green’s functions (NEGF), of two state-of-the-art two-well (TW) Terahertz Quantum Cascade Lasers (THz QCLs) supporting clean 3-level systems. The devices have nearly identical parameters and the NEGF calculations with an abrupt-interface roughness height of 0.12 nm predict a maximum operating temperature (T_max) of ~ 250 K for both devices. However, experimentally, one device reaches a T_max of ~ 250 K and the other a T_max of only ~ 134 K. Both devices were fabricated and measured under identical conditions in the same laboratory, with high quality processes as verified by reference devices. The main difference between the two devices is that they were grown in different MBE reactors. Our NEGF-based analysis considered all parameters related to MBE growth, including the maximum estimated variation in aluminum content, growth rate, doping density, background doping, and abrupt-interface roughness height. From our NEGF calculations it is evident that the sole parameter to which a drastic drop in T_max could be attributed is the abrupt-interface roughness height. We can also learn from the simulations that both devices exhibit high-quality interfaces, with one having an abrupt-interface roughness height of approximately an atomic layer and the other approximately a monolayer. However, these small differences in interface sharpness are the cause of the large performance discrepancy. This underscores the sensitivity of device performance to interface roughness and emphasizes its strategic role in achieving higher operating temperatures for THz QCLs. We suggest Atom Probe Tomography (APT) as a path to analyze and measure the (graded)-interfaces roughness (IFR) parameters for THz QCLs, and subsequently as a design tool for higher performance THz QCLs, as was done for mid-IR QCLs. Our study not only addresses challenges faced by other groups in reproducing the record T_max of ~ 250 K and ~ 261 K but also proposes a systematic pathway for further improving the temperature performance of THz QCLs beyond the state-of-the-art.

Finding alloys with specific design properties is challenging due to the large number of possible compositions and the complex interactions between elements. This study introduces a multi-objective Bayesian optimization approach guiding molecular dynamics simulations for discovering high-performance refractory alloys with both targeted intrinsic static thermomechanical properties and also deformation mechanisms occurring during dynamic loading. The objective functions are aiming for excellent thermomechanical stability via a high bulk modulus, a low thermal expansion, a high heat capacity, and for a resilient deformation mechanism maximizing the retention of the BCC phase after shock loading. Contrasting two optimization procedures, we show that the Pareto-optimal solutions are confined to a small performance space when the property objectives display a cooperative relationship. Conversely, the Pareto front is much broader in the performance space when these properties have antagonistic relationships. Density functional theory simulations validate these findings and unveil underlying atomic-bond changes driving property improvements.

Understanding and accurately characterizing energy dissipation mechanisms in civil structures during earthquakes is an important element of seismic assessment and design. The most commonly used model is attributed to Rayleigh. This paper proposes a systematic approach to quantify the uncertainty associated with Rayleigh's damping model. Bayesian calibration with embedded model error is employed to treat the coefficients of the Rayleigh model as random variables using modal damping ratios. Through a numerical example, we illustrate how this approach works and how the calibrated model can address modeling uncertainty associated with the Rayleigh damping model.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Topological insulator–magnetic insulator (TI–MI) heterostructures hold significant promise in the field of spintronics, offering the potential for manipulating magnetization through topological surface state–enabled spin-orbit torque. However, many TI–MI interfaces are plagued by issues such as contamination within the magnetic insulator layer and the presence of a low-density transitional region of the topological insulator. These interfacial challenges often obscure the intrinsic behavior of the TI–MI system. In this study, we addressed these challenges by depositing sputtered Bi₂⁢Te₃ on liquid phase epitaxy grown Y₃⁢Fe₅⁢O₁₂/Gd₃⁢Ga₅⁢O₁₂. The liquid phase epitaxy grown Y₃⁢Fe₅⁢O₁₂ has been previously shown to have exceptional interface quality, without an extended transient layer derived from interdiffusion processes of the substrate or impurity ions, thereby eliminating rare-earth impurity-related losses in the MI at low temperatures. At the TI–MI interface, high-resolution depth-sensitive polarized neutron reflectometry confirmed the absence of a low-density transitional growth region of the TI. By overcoming these undesirable interfacial effects, we isolate and probe the intrinsic low-temperature magnetization dynamics and transport properties of the TI–MI interface. In conclusion, our findings revealed strong spin pumping at low temperatures, accompanied by an additional in-plane anisotropy. The enhanced spin pumping at low temperatures is correlated with the observed suppression of bulk conduction and the weak antilocalization in the TI film, highlighting the interplay between the transport and spin pumping behavior in the TI–MI system.

Solvent expulsion away from an intervening region between two approaching particles plays important roles in particle aggregation yet remains poorly understood. In this work, we use metadynamics molecular simulations to study the free energy landscape of removing water molecules from gibbsite and pyrophyllite slit pores representing the confined spaces between two approaching particles. For gibbsite, removing water from the intervening region is both entropically and enthalpically unfavorable. The closer the particles approach each other, the harder it is to expel water molecules. For pyrophyllite, water expulsion is spontaneous, which is different from the gibbsite system. A smaller pore makes the water removal more favorable. When water is being drained from the intervening region, single chains of water molecules are observed in gibbsite pore, while in pyrophyllite pore water cluster is usually observed. Water-gibbsite hydrogen bonds help stabilize water chains, while water forms clusters in pyrophyllite pore to maximize the number of hydrogen bonds among themselves. This work provides the first assessment into the energetics and structure of water being drained from the intervening region between two approaching particles during oriented attachment and aggregation.

We present a comprehensive benchmarking framework for evaluating machine-learning approaches applied to phase-field problems. This framework focuses on four key analysis areas crucial for assessing the performance of such approaches in a systematic and structured way. Firstly, interpolation tasks are examined to identify trends in prediction accuracy and accumulation of error over simulation time. Secondly, extrapolation tasks are also evaluated according to the same metrics. Thirdly, the relationship between model performance and data requirements is investigated to understand the impact on predictions and robustness of these approaches. Finally, systematic errors are analyzed to identify specific events or inadvertent rare events triggering high errors. Quantitative metrics evaluating the local and global description of the microstructure evolution, along with other scalar metrics representative of phase-field problems, are used across these four analysis areas. This benchmarking framework provides a path to evaluate the effectiveness and limitations of machine-learning strategies applied to phase-field problems, ultimately facilitating their practical application.

This study investigates high performance electrochromic windows used on a passive house and residential dwelling to IECC 2021 (i.e., IECC dwelling). In the lab, the electrochromic film switches transmitted solar heat gain coefficient (SHGC) from 0.09 to 0.7 and visible transmittance from 0.15 to 0.82 with power consumption of 1.23 W/m2 during switching times less than 3 minutes. We extrapolate these results to a window assembly. Building energy models of the houses were evaluated in Santa Fe, New Mexico. A Monte Carlo analysis for 2020, 2040, 2060, and 2080 was conducted for Shared Socioeconomic Pathways 2-4.5, 3-7.0, and 5-8.5. Cases with and without the electrochromic windows and with and without electricity were used to determine energy use intensity and hours beyond thermal safety thresholds. The passive house showed 1.3-3.1% mean energy savings and the IECC dwelling 4.4-5.1% with electrochromic efficiency benefits growing into the future for both cases. Even so, overall savings decrease into the future for the passive house, due to growth in cooling load being dominant, conversely overall energy savings increase into the future for the IECC dwelling due to heating loads being dominant. For thermal resilience, the passive house exhibited a mean percent decrease of 0.02-0.31% hours in the extreme caution (i.e., > 32.2 ∘C, ≤ 39.4 ∘C) range while the IECC dwelling exhibited 0.38-4.38%. The study therefore shows that electrochromic windows will have smaller benefits for the passive house in comparison to the IECC dwelling. The relationship between electrochromic windows is shown to have a complex relationship between house efficiency and climate change by these results.

Here, it is demonstrated that single-crystalline and highly doped GaAs membranes are excellent candidates for realizing infrared-transparent shields of electromagnetic interference at millimeter frequencies. Measured optical transmittance spectra for the semiconductor membranes show resonant features between 750 and 2500 nm, with a 100% maximum transmittance. The shielding effectiveness of the membranes is extracted from measured scattering parameters between 65 and 85 GHz. Selected GaAs membranes and membranes/polyamide films exhibit shielding effectiveness ranging from 22 to 40 dB, which are suitable values to ensure the safe operation of infrared devices for commercial applications. Theoretical calculations based on a plane wave model show that the interplay of primary reflection and multiple internal reflections of the radio-frequency waves results in broadband shielding capabilities of the membrane between 10 and 300 GHz.

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

The laser system is the most complex component of a light-pulse atom interferometer (LPAI), controlling frequencies and intensities of multiple laser beams to configure quantum gravity and inertial sensors. Its main functions include cold-atom generation, state preparation, state-selective detection, and generating a coherent two-photon process for the light-pulse sequence. To achieve substantial miniaturization and ruggedization, we integrate key laser system functions onto a photonic integrated circuit. Our study focuses on a high-performance silicon photonic suppressed-carrier single-sideband (SC-SSB) modulator at 1560 nanometers, capable of dynamic frequency shifting within the LPAI. By independently controlling radio frequency (RF) channels, we achieve 30-decibel carrier suppression and unprecedented 47.8-decibel sideband suppression at peak conversion efficiency of -6.846 decibels (20.7%). We investigate imbalances in both amplitudes and phases between the RF signals. Using this modulator, we demonstrate cold-atom generation, state-selective detection, and atom interferometer fringes to estimate gravitational acceleration, g ≈ 9.77 ± 0.01 meters per second squared, in a rubidium (87Rb) atom system.

Helium-4-based scintillation detector technology is emerging as a strong alternative to pulse-shape discrimination-capable organic scintillators for fast neutron detection and spectroscopy, particularly in extreme gamma-ray environments. The ⁴He detector is intrinsically insensitive to gamma radiation, as it has a relatively low cross-section for gamma-ray interactions, and the stopping power of electrons in the ⁴He medium is low compared to that of ⁴He recoil nuclei. Consequently, gamma rays can be discriminated by simple energy deposition thresholding instead of the more complex pulse shape analysis. The energy resolution of ⁴He scintillation detectors has not yet been well-characterized over a broad range of energy depositions, which limits the ability to deconvolve the source spectra. In this work, an experiment was performed to characterize the response of an Arktis S670 ⁴He detector to nuclear recoils up to 9 MeV. The ⁴He detector was positioned in the center of a semicircular array of organic scintillation detectors operated in coincidence. Deuterium–deuterium and deuterium–tritium neutron generators provided monoenergetic neutrons, yielding geometrically constrained nuclear recoils ranging from 0.0925 to 8.87 MeV. The detector response provides evidence for scintillation linearity beyond the previously reported energy range. Finally, the measured response was used to develop an energy resolution function applicable to this energy range for use in high-fidelity detector simulations needed by future applications.

In global sensitivity analysis (GSA) of a model, a proper convergence analysis of metrics is essential for ensuring a level of confidence or trustworthiness in sensitivity results obtained, yet is somewhat deficient in practice. The level of confidence in sensitivity measures, particularly in relation to their influence and support for decisions from scientific, social and policy perspectives, is heavily reliant on the convergence of GSA. We review the literature and summarize the available methods for monitoring and assessing convergence of sensitivity measures based on application purposes. The aim is to expose the various choices for convergence assessment and encourage further testing of available methods to clarify their level of robustness. Furthermore, the review identifies a pressing need for comparative studies on convergence assessment methods to establish a clear hierarchy of effectiveness and encourages the adoption of systematic approaches for enhanced robustness in sensitivity analysis.

Individual lanthanide elements have physical/electronic/magnetic properties that make each useful for specific applications. Several of the lanthanides cations (Ln3+) naturally occur together in the same ores. They are notoriously difficult to separate from each other due to their chemical similarity. Predicting the Ln3+ differential binding energies (ΔΔE) or free energies (ΔΔG) at different binding sites, which are key figures of merit for separation applications, will help design of materials with lanthanide selectivity. We apply ab initio molecular dynamics (AIMD) simulations and density functional theory (DFT) to calculate ΔΔG for Ln3+ coordinated to ligands in water and embedded in metal-organic frameworks (MOFs), and ΔΔE for Ln3+ bonded to functionalized silica surfaces, thus circumventing the need for the computational costly absolute binding (free) energies ΔG and ΔE. Perturbative AIMD simulations of water-inundated simulation cells are applied to examine the selectivity of ligands towards adjacent Ln3+ in the periodic table. Static DFT calculations with a full Ln3+ first coordination shell, while less rigorous, show that all ligands examined with net negative charges are more selective towards the heavier lanthanides than a charge-neutral coordination shell made up of water molecules. Amine groups are predicted to be poor ligands for lanthanide-binding. We also address cooperative ion binding, i.e., using different ligands in concert to enhance lanthanide selectivity.

Analytical and semi–analytical models for stream depletion with transient stream stage drawdown induced by groundwater pumping are developed to address a deficiency in existing models, namely, the use of a fixed stream stage condition at the stream–aquifer interface. Here field data are presented to demonstrate that stream stage drawdown does indeed occur in response to groundwater pumping near aquifer–connected streams. A model that predicts stream depletion with transient stream drawdown is developed based on stream channel mass conservation and finite stream channel storage. The resulting models are shown to reduce to existing fixed–stage models in the limit as stream channel storage becomes infinitely large, and to the confined aquifer flow with a no–flow boundary at the streambed in the limit as stream storage becomes vanishingly small. The model is applied to field measurements of aquifer and stream drawdown, giving estimates of aquifer hydraulic parameters, streambed conductance, and a measure of stream channel storage. The results of the modeling and data analysis presented herein have implications for sustainable groundwater management.

Boron nitride nanotubes (BNNTs) are high-strength, high-modulus nanotubes with high thermal and oxidative stabilities. Two hybrid composites were prepared with satin weave carbon fiber (CF) and resole-type phenolic resin: one with surface layers of BNNTs and one with alternating interlayers of BNNTs. The samples were subjected to hot jet tests that simulate realistic high-pressure-temperature conditions to understand the behavior of BNNTs under high-pressure erosion. Adding BNNTs to CF/phenolic laminates enhanced the ablation resistance by reinforcing the char material and mitigated localized thermal damage. Hybrid laminates exhibited up to 14% lower weight loss, 55% increase in flexural modulus, higher thermal diffusivity, and improved char yield and microstructure compared to CF/phenolic samples. The surface layer hybrid had many surviving nanotubes reinforcing the char and crystalline oxide structures that could mitigate further oxygen diffusion. Further, various characterization methods were used to deduce possible mechanisms and their products, indicating that BNNTs could serve as growth templates for direct crystalline boron oxide formation. Overall, hybrid BNNT/CF/phenolic laminates displayed better ablation resistance and favorable microstructure evolution under high-pressure conditions.

Materials simulations based on direct numerical solvers are accurate but computationally expensive for predicting materials evolution across length- and time-scales, due to the complexity of the underlying evolution equations, the nature of multiscale spatiotemporal interactions, and the need to reach long-time integration. We develop a method that blends direct numerical solvers with neural operators to accelerate such simulations. This methodology is based on the integration of a community numerical solver with a U-Net neural operator, enhanced by a temporal-conditioning mechanism to enable accurate extrapolation and efficient time-to-solution predictions of the dynamics. We demonstrate the effectiveness of this hybrid framework on simulations of microstructure evolution via the phase-field method. Such simulations exhibit high spatial gradients and the co-evolution of different material phases with simultaneous slow and fast materials dynamics. We establish accurate extrapolation of the coupled solver with large speed-up compared to DNS depending on the hybrid strategy utilized. This methodology is generalizable to a broad range of materials simulations, from solid mechanics to fluid dynamics, geophysics, climate, and more.

This dataset is comprised of a library of atomistic structure files and corresponding X-ray diffraction (XRD) profiles and vibrational density of states (VDoS) profiles for bulk single crystal silicon (Si), gold (Au), magnesium (Mg), and iron (Fe) with and without disorder introduced into the atomic structure and with and without mechanical loading. Included with the atomistic structure files are descriptor files that measure the stress state, phase fractions, and dislocation content of the microstructures. All data was generated via molecular dynamics or molecular statics simulations using the Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) code. This dataset can inform the understanding of how local or global changes to a materials microstructure can alter their spectroscopic and diffraction behavior across a variety of initial structure types (cubic diamond, face-centered cubic (FCC), hexagonal close-packed (HCP), and body-centered cubic (BCC) for Si, Au, Mg, and Fe, respectively) and overlapping changes to the microstructure (i.e., both disorder insertion and mechanical loading).

Efficient carbon capture requires engineered porous systems that selectively capture CO2 and have low energy regeneration pathways. Porous liquids (PLs), solvent-based systems containing permanent porosity through the incorporation of a porous host, increase the CO2 adsorption capacity. A proposed mechanism of PL regeneration is the application of isostatic pressure in which the dissolved nanoporous host is compressed to alter the stability of gases in the internal pore. This regeneration mechanism relies on the flexibility of the porous host, which can be evaluated through molecular simulations. Here, the flexibility of porous organic cages (POCs) as representative porous hosts was evaluated, during which pore windows decreased by 10-40% at 6 GPa. POCs with sterically smaller functional groups, such as the 1,2-ethane in the CC1 POC resulted in greater imine cage flexibility relative to those with sterically larger functional groups, such as the cyclohexane in the CC3 POC that protected the imine cage from the application of pressure. Structural changes in the POC also caused CO2 adsorption to be thermodynamically unfavorable beginning at ∼2.2 GPa in the CC1 POC, ∼1.1 GPa in the CC3 POC, and ∼1.0 GPa in the CC13 POC, indicating that the CO2 would be expelled from the POC at or above these pressures. Energy barriers for CO2 desorption from inside the POC varied based on the geometry of the pore window and all the POCs had at least one pore window with a sufficiently low energy barrier to allow for CO2 desorption under ambient temperatures. The results identified that flexibility of the CC1, CC3, or CC13 POCs under compression can result in the expulsion of captured gas molecules.

Abstract not provided.

Abstract not provided.

5 min panel presentation for the EPRI connector workshop.

Abstract not provided.

Abstract not provided.

Abstract not provided.

This is the seminar I will present at WCCM conference highlighting our latest research work on incorporating genetic programming to obtain data-driven strength models for complex materials.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

This is the seminar I will present at WCCM conference highlighting our latest research work on incorporating genetic programming to obtain data-driven strength models for complex materials.

The panel discussion titled "Innovations in Advanced Cybersecurity Technologies for Safety Systems Within Nuclear Power Plants (NPPs)" took place at the American Nuclear Society (ANS) Annual Meeting in Las Vegas on Tuesday, June 18, 2024. The session was moderated by Lee Maccarone of Sandia National Laboratories (SNL) and featured four experts: Benjamin Karch from Sandia National Laboratories, Ismael Garcia from the U.S. Nuclear Regulatory Commission (NRC), Trevor Taylor from Paragon Energy Solutions (ES), and Forrest Shriver from Sentinel Devices.

Abstract not provided.

In order to make design decisions, engineers may seek to identify regions of the design domain that are acceptable in a computationally efficient manner. A design is typically considered acceptable if its reliability with respect to parametric uncertainty exceeds the designer’s desired level of confidence. Despite major advancements in reliability estimation and in design classification via decision boundary estimation, the current literature still lacks a design classification strategy that incorporates parametric uncertainty and desired design confidence. To address this gap, this works offers a novel interpretation of the acceptance region by defining the decision boundary as the hypersurface which isolates the designs that exceed a user-defined level of confidence given parametric uncertainty. This work addresses the construction of this novel decision boundary using computationally efficient algorithms that were developed for reliability analysis and decision boundary estimation. The proposed approach is verified on two physical examples from structural and thermal analysis using Support Vector Machines and Efficient Global Optimization-based contour estimation.

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.

Presentation for the GeoMIP 2024 Meeting

Abstract not provided.

Presentation for the OPTICA Laser Applications to Chemical Security and Environmental Analysis

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.

Molecular dynamics (MD) has been an incredibly useful tool to model physical processes that were synthesized experimentally but not fully understood. MD, through the use of semi-empirical inter-atomic potentials, has allowed understanding of different physical processes in materials science. Yet as well as providing useful insights into materials science, molecular dynamics has a wider range of usability. In this report, I will be detailing how MD can be used to study graphene formation from a carbon liquid which requires high temperatures and pressures. Beyond this, I will describe the usefulness of MD for understanding the physics for phonon transport quantum sensors. To do this, MD was employed to determine the dynamical matrix by treating atoms as coupled oscillators. An accurate understanding of the dynamical matrix of a system is required to calculate the non-equilibrium Green’s function used to describe the phonon transport within phonon wave-guides. I found that, across multiple pressures and temperatures, randomly placed carbon atoms will show evidence of pent-first formation with semi-empirical models. Density functional theory (DFT), on the other hand, was too computationally expensive to use for full scale MD simulations, but we have the possibility of training a machine learned interatomic potential to approximate DFT for carbon in the environments being studied for pent-first graphene sheet formation.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

A presentation on MRPIB for the 2024 INMM meeting.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Poster to be presented at NAMBE to share results of recent research

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

MAD³ (Material Data Driven Design) is a novel and unique software solution that provides initial plastic anisotropy of polycrystalline metals using crystallographic texture information, developed at Sandia National Laboratories. In this document, we describe the structure and functionality of the current MAD³ software (v1.01).

Abstract not provided.

Abstract not provided.

Overview of 9321's Computing Capabilities related to Machine Learning

There is interest in increasing the size and range of the electron beams while maintaining their high dose rate. To accomplish this a Pulsed Power Electron Gun (PPEG) has been developed based around a pulsed (30ns) 300 kV Marx Generator (Applied Physical Electronics L.C.).

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.

Presentation slides

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Poster presentation at DOE OE Peer Review 2024 in Bellevue Washington

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

A pulsed laser testing system for rapid screening of single-event latch-up (SEL) in microelectronics with real-time visualization of sensitive nodes via photo-emission is presented. We provide a detailed overview of the system, followed by SEL threshold measurements in an analog-to-digital converter (AD9240). We find that SEL thresholds are comparable to measurements taken with a low-repetition rate, fixed beam system, and demonstrate full-die scanning speeds can be an order of magnitude faster

Abstract not provided.

Abstract not provided.

High-throughput image segmentation of atomic resolution electron microscopy data poses an ongoing challenge for materials characterization. In this paper, we investigate the application of the polyhedral template matching (PTM) method, a technique widely employed for visualizing three-dimensional (3D) atomistic simulations, to the analysis of two-dimensional (2D) atomic resolution electron microscopy images. This technique is complementary with other atomic resolution data reduction techniques, such as the centrosymmetry parameter, that use the measured atomic peak positions as the starting input. Furthermore, since the template matching process also gives a measure of the local rotation, the method can be used to segment images based on local orientation. We begin by presenting a 2D implementation of the PTM method, suitable for atomic resolution images. We then demonstrate the technique's application to atomic resolution scanning transmission electron microscopy images from close-packed metals, providing examples of the analysis of twins and other grain boundaries in FCC gold and martensite phases in 304 L austenitic stainless steel. Finally, we discuss factors, such as positional errors in the image peak locations, that can affect the accuracy and sensitivity of the structural determinations.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Stereo high-speed video of photovoltaic modules undergoing laboratory hail tests was processed using digital image correlation to determine module surface deformation during and immediately following impact. The purpose of this work was to demonstrate a methodology for characterizing module impact response differences as a function of construction and incident hail parameters. Video capture and digital image analysis were able to capture out-of-plane module deformation to a resolution of ±0.1 mm at 11 kHz on an in-plane grid of 10 × 10 mm over the area of a 1 × 2 m commercial photovoltaic module. With lighting and optical adjustments, the technique was adaptable to arbitrary module designs, including size, backsheet color, and cell interconnection. Impacts were observed to produce an initially localized dimple in the glass surface, with peak deflection proportional to the square root of incident energy. Subsequent deformation propagation and dissipation were also captured, along with behavior for instances when the module glass fractured. Natural frequencies of the module were identifiable by analyzing module oscillations postimpact. Limitations of the measurement technique were that the impacting ice ball obscured the data field immediately surrounding the point of contact, and both ice and glass fracture events occurred within 100 μs, which was not resolvable at the chosen frame rate. Increasing the frame rate and visualizing the back surface of the impact could be applied to avoid these issues. Applications for these data include validating computational models for hail impacts, identifying the natural frequencies of a module, and identifying damage initiation mechanisms.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Junctions are fundamental elements that support qubit locomotion in two-dimensional ion trap arrays and enhance connectivity in emerging trapped-ion quantum computers. In surface ion traps they have typically been implemented by shaping radio frequency (RF) electrodes in a single plane to minimize the disturbance to the pseudopotential. However, this method introduces issues related to RF lead routing that can increase power dissipation and the likelihood of voltage breakdown. Here, we propose and simulate a novel two-layer junction design incorporating two perpendicularly rotoreflected (rotated, then reflected) linear ion traps. The traps are vertically separated, and create a trapping potential between their respective planes. The orthogonal orientation of the RF electrodes of each trap relative to the other provides perpendicular axes of confinement that can be used to realize transport in two dimensions. While this design introduces manufacturing and operating challenges, as now two separate structures have to be precisely positioned relative to each other in the vertical direction and optical access from the top is obscured, it obviates the need to route RF leads below the top surface of the trap and eliminates the pseudopotential bumps that occur in typical junctions. In this paper the stability of idealized ion transfer in the new configuration is demonstrated, both by solving the Mathieu equation analytically to identify the stable regions and by numerically modeling ion dynamics. Our novel junction layout has the potential to enhance the flexibility of microfabricated ion trap control to enable large-scale trapped-ion quantum computing.

Abstract not provided.

Abstract not provided.

Abstract not provided.