Publications

Abstract not provided.

Abstract not provided.

Abstract not provided.

The Sea Change Project, managed and financed by SWITCH Maritime, produced the first zero-emission 100% hydrogen fuel-cell powered commercial ferry in the world. The MV Sea Change is a 70-foot catamaran equipped with a 360-kW fuel cell powertrain and a 250-bar 246 kg capacity hydrogen storage system. The 75-passenger vessel began public passenger service in July 2024 as part of the San Francisco Bay Ferry system. Summarized here are aspects of the design of the vessel, its construction, and vessel bunkering with high-pressure hydrogen. Results from a fuel cell performance analysis are also presented, along with estimates of the greenhouse gas (GHG) and criterial pollutant emissions assuming different sourcings of hydrogen. Lessons learned from the project will be briefly reviewed.

FLEXO (Flux-Limited Extended-MHD Ohm's Law) is a production-line multiphysics code developed at Sandia to enable more predictive modeling of target physics on pulsed-power devices. FLEXO uses an extended magnetohydrodynamics (XMHD) model which includes a generalized Ohm's law (GOL), an electron inertia term, and Hall physics. This report describes the code's numerical methods, its computational performance, and test problems of interest.

Drill rig parameter measurements are routinely used during deep well construction to monitor and guide drilling conditions for improved performance and reduced costs. While insightful into the drilling process, these measurements are of reduced value without a standard to aid in data evaluation and decision making. In the main body of this work (Volume 1), a method is demonstrated whereby rock reduction model constraints are used to interpret drilling response parameters; the method could be applied in real-time to improve decision-making in the field and to further discern technology performance during post-drilling evaluations. Drilling parameters are evaluated using laboratory-validated rock reduction models for predicting the phenomenological response of drag bits (Detournay and Defourny, 1992) in computational algorithms. The method presented has applicability to development of advanced analytics on future geothermal wells using real-time electronic data recording for improved performance and reduced drilling costs. A drilling cost model is also used to show the tradeoff between rate of penetration and bit life and the influence on interval drilling costs. Details of the bit specifications and performance are cataloged in an independent volume, documented under separate cover, for each of the four wells, and include Volume 2: Utah FORGE 16A(78)-32; Volume 3: Utah FORGE 56-32; Volume 4: Utah FORGE 78B-32 and Volume 5: Utah FORGE 16B(78)-32.

Increasing the fluence of z-pinch x-ray radiation sources above ∼ 10 keV has been a long-standing goal for scientists at Sandia National Laboratories’ Z Machine. Optimizing sources for non-thermal “cold Kα” emission in higher atomic-number materials appears to be a promising path to increase warm x-ray yield. However, this emission is generated by supra-thermal electrons, which are not treated in the magnetohydrodynamic (MHD) codes that are typically used in z-pinch source development. MHD codes do not allow for charge separation or space-charge-generated electric fields, and constrain particle kinematics to Maxwellian distributions. The kinetic codes which do accommodate discrete, non-thermal energy distributions are computationally prohibitive when modeling plasmas near solid density and when modeling/tracking higher ionization states. Thus, modeling non-thermal z-pinch sources requires a new simulation tool. In this report, we present a new hybrid modeling capability that uses the fast features of MHD-type particles to the greatest extent possible, then transitions to the slower but more complete kinetic particle treatment to correctly capture the particle energy spectra that generate non-thermal emission. This capability is founded on the fully-relativistic particle-in-cell code Chicago, which already includes fluid particle treatments. The governing equations and hybrid methodology presented here are applied in simulations of an argon gas-puff and a molybdenum wire-array to provide preliminary code validation. The argon simulation is compared to measured implosion times and yields from Jones et al., Phys. Plasmas 22, 020706 (2015). The simulated x-ray yield is within 25% of measurements and the implosion times agree within a few percent. The molybdenum wire array simulation captures the implosion timing reported in Hansen et al., Phys. Plasmas 21, 031202 (2014), but work is needed to verify the available EOS table. These exemplar simulations represents the type of non-thermal sources that will be developed using the hybrid code capability going forward.

Abstract not provided.

Remote observation of infrasound resonant signals emanating from underground tunnel structures could potentially allow the remote quantification of the geometry of the underground structures in which the signals were generated. However, the sensitivity of these observations to tunnel geometry and changes in that geometry are unknown. In this report we outline a numerical simulation study with the following three objectives: 1. Can we model infrasound resonance using numerical simulations? 2. What is the sensitivity of the modeled observations to changes in tunnel geometry or boundary conditions? And 3. Can we accurately model resonant observations of explosions occurring in Redmond Salt Mine in Central Utah. In this report we outline affirmative answers to the first two objectives, but we were unable to accurately model the Redmond Explosions due to numerical instability in the model of the complex structure of the mine. Recommendations for future work emphasize the need to acquire additional datasets and to explore a more data-based approach in which changes in data signatures are detected as a first step towards developing a method which inverts resonant infrasound signals for tunnel geometry.

Drill rig parameter measurements are routinely used during deep well construction to monitor and guide drilling conditions for improved performance and reduced costs. While insightful into the drilling process, these measurements are of reduced value without a standard to aid in data evaluation and decision making. In the main body of this work (Volume 1), a method is demonstrated whereby rock reduction model constraints are used to interpret drilling response parameters; the method could be applied in real-time to improve decision-making in the field and to further discern technology performance during post-drilling evaluations. Drilling parameters are evaluated using laboratory-validated rock reduction models for predicting the phenomenological response of drag bits (Detournay and Defourny, 1992) in computational algorithms. The method presented has applicability to development of advanced analytics on future geothermal wells using real-time electronic data recording for improved performance and reduced drilling costs. A drilling cost model is also used to show the tradeoff between rate of penetration and bit life and the influence on interval drilling costs. Details of the bit specifications and performance are cataloged in an independent volume, documented under separate cover, for each of the four wells, and include Volume 2: Utah FORGE 16A(78)-32; Volume 3: Utah FORGE 56-32; Volume 4: Utah FORGE 78B-32 and Volume 5: Utah FORGE 16B(78)-32.

Abstract not provided.

Nitrogen-containing organic molecules play essential roles in nearly every aspect of chemistry. The importance of these species within gas-phase chemistry has increased in recent years, due to their emissions from wildfires, crude biofuels, ammonia combustion, and CO₂ capture facilities. However, there is a lack of detailed relevant mechanistic and kinetic studies of N-compounds, resulting in poor representation of these pathways within chemical models. In this report, we detail our exploration of the gas-phase oxidation mechanisms of pyrrole, imidazole, pyrrolidine, methylamine, and dimethylamine, done by employing multiplexed VUV photoionization mass spectrometry coupled with KinBot, a computational tool that automatically explores multi-well potential energy surfaces. Overall, we vastly expanded the known gas-phase oxidation mechanisms and chemical kinetics of these species. This work will improve the accuracy and completeness of atmospheric and combustion chemistry models, ultimately leading to better-informed decision-makers in the energy and environmental sectors.

MELCOR has been used extensively to facilitate virtual investigations into severe nuclear accidents for light-water reactors (LWRs). Non-light water reactors (non-LWRs) render some LWR-centric approaches potentially unsuitable. MELCOR has been instrumental in analyzing source terms for LWRs and has recently expanded its applicability to non-LWRs. To simplify radionuclide (RN) tracking, MELCOR currently groups elements into 17 classes, each containing representative species. This grouping, optimized for LWRs, is not appropriate for non-LWRs due to the different chemistry. This necessitates a reevaluation of radionuclide transport modeling. This report introduces a new class scheme for MELCOR tailored to MSR modeling, expanding the current 17 classes to 32 and are explained in the context of a UF₄ fueled FLiBe carrier MSR. It provides a discussion and justification for the new groupings and outlines a methodology for discovering and defining additional classes in MELCOR using a sample calculated RN inventory and a Gibbs energy minimizer (GEM).

Abstract not provided.

Slides for a 2 hours presentation at a summer school on batteries. All the material in the presentation has been previously published.

Achieving maximum performance in a Z-pinch Magnetized Liner Inertial Fusion (MagLIF) implosion requires high levels of symmetry to reach high pressures, high areal densities, and high efficiency energy coupling between the compressed fuel and the imploding shell. In this work, we present the development of a nonlinear model to describe an asymmetric MagLIF liner undergoing magneto-Rayleigh-Taylor (MRT) instabilities in the z = 0 plane. We investigate and present results from our model that show the degradation of the pressure at stagnation under various kinds of asymmetries that can arise during MagLIF experiments.

Abstract not provided.

Abstract not provided.

A novel approach is presented for parametric analysis of remotely-sensed ground and cloud clutter. A spatial-frequency-domain clutter model is generated from an extensive, one-year database of weather imagery and statistics are given for each spatial frequency. This approach is useful for the analysis and design of spatial and temporal clutter-rejection filters, which can also be analyzed in this domain.

Bugs in digital logic have led to some significant security vulnerabilities. Hardware bugs are particularly troublesome since they cannot be easily patched. Additionally, if the bug is in the root of trust, all trust built upon it can be vulnerable. Traditional testing either require a deep knowledge of the system, creative attack vectors and lots of human interaction. This is not scalable as there are very few engineers that can wear the hat of a designer, a verification engineer, and a cybersecurity expert. Hardware fuzzing is a relatively new research area in dynamic hardware testing. It has proven to be an effective method for discovering bugs, unexpected behaviors, and security vulnerabilities in software. While hardware fuzzing is new to the hardware domain, it has a strong track record in software testing. Fuzzing is a testing technique that randomly mutates the input data to uncover bugs or vulnerabilities in the design. It is especially good at finding corner cases that test engineers can not envision. Another advantage over other dynamic testing techniques is that, if done well, deep knowledge of the design is not required. Additionally, fuzzing scales well. If the system is set up correctly, it can run unsupervised for weeks if necessary. In this work, we propose using hardware fuzzing to improve the input vector generation for an information flow tracking tool. To get reasonable throughput of test vectors, an emulator is targeted as the execution platform. Efficient emulator execution has some specific requirements.

Abstract not provided.

U.S. nuclear power facilities face increasing challenges in meeting dynamic security requirements caused by evolving and expanding threats while keeping costs reasonable to make nuclear energy competitive. The past approach has often included implementing security features after a facility has been designed and without attention to optimization, which can lead to cost overruns. Incorporating security into the design process can provide robust, cost-effective, and sufficient physical protection systems. The purpose of this report is to capture lessons learned by the Advanced Reactor Safeguards and Security (ARSS) program that may be beneficial for other advanced and small modular reactor (SMR) vendors to use when developing security systems and postures. This report will capture relevant information that can be used in the security-by-design (SeBD) process for SMR and microreactor vendors.

We study the problem of multifidelity uncertainty propagation for computationally expensive models. In particular, we consider the general setting where the high-fidelity and low-fidelity models have a dissimilar parameterization both in terms of number of random inputs and their probability distributions, which can be either known in closed form or provided through samples. We derive novel multifidelity Monte Carlo estimators which rely on a shared subspace between the high-fidelity and low-fidelity models where the parameters follow the same probability distribution, i.e., a standard Gaussian. We build the shared space employing normalizing flows to map different probability distributions into a common one, together with linear and nonlinear dimensionality reduction techniques, active subspaces and autoencoders, respectively, which capture the subspaces where the models vary the most. We then compose the existing low-fidelity model with these transformations and construct modified models with an increased correlation with the high-fidelity model, which therefore yield multifidelity estimators with reduced variance. A series of numerical experiments illustrate the properties and advantages of our approaches.

Abstract not provided.

Seismic imaging methods are critical for Global Security and Energy & Homeland Security missions and activities that rely on subsurface characterization, but traditional methods remain computationally expensive and require significant labor hours and expertise to execute. Within the past few years, machine learning (ML), namely deep learning (DL), has been used to develop data-driven end-to-end full waveform inversion (FWI) methods to estimate 2D P-wave velocity (Vp) models in a fraction of the time as conventional FWI. These methods, however, are trained on simplistic acoustic wave seismic data and Vp models that are not realistic nor representative of real-world observations, leaving a large gap between the state-of-the-art and deployable, feasible, and practical DL FWI methods. Here, we generate a synthetic active-source, 3D, elastic wave seismic data set and a variety of Vp models with realistic geologic structure for training DL FWI methods. We evaluate six different methods that have performed well for acoustic DL FWI or medical imaging tasks using our more realistic dataset. We find that these six trained models do not match the performance of published acoustic end-to-end DL FWI methods, indicating more training data may be needed, physics may need to be incorporated to achieve good accuracy at the sacrifice of the end-to-end advantage, and/or novel methods need to be developed to enable end-to-end DL FWI methods to perform well for real-world seismic data.

Beyond the title, this talk is about improving the computational modeling of germanium hole spins in lithographic quantum dots to work towards a better understanding of how potential qubit properties of these systems are affected by the effective shell filling of these holes in the electric potential.

For multi-scale multi-physics applications e.g., the turbulent combustion code Pele, robust and accurate dimensionality reduction is crucial to solving problems at exascale and beyond. A recently developed technique, Co-Kurtosis based Principal Component Analysis (CoK-PCA) which leverages principal vectors of co-kurtosis, is a promising alternative to traditional PCA for complex chemical systems. To improve the effectiveness of this approach, we employ Artificial Neural Networks for reconstructing thermo-chemical scalars, species production rates, and overall heat release rates corresponding to the full state space. Our focus is on bolstering confidence in this deep learning based non-linear reconstruction through Uncertainty Quantification (UQ) and Sensitivity Analysis (SA). UQ involves quantifying uncertainties in inputs and outputs, while SA identifies influential inputs. One of the noteworthy challenges is the computational expense inherent in both endeavors. To address this, we employ the Monte Carlo methods to effectively quantify and propagate uncertainties in our reduced spaces while managing computational demands. Our research carries profound implications not only for the realm of combustion modeling but also for a broader audience in UQ. By showcasing the reliability and robustness of CoK-PCA in dimensionality reduction and deep learning predictions, we empower researchers and decision-makers to navigate complex combustion systems with greater confidence.

Abstract not provided.

Computer programs are often viewed as collections of functions – each function has parameters (inputs) and computes a return value, and each has potential side effects that modify program state (outputs). In this research, a Sandia symbolic execution tool designed to support “human-in-the-loop” analysis was modified to automatically create “function summaries,” and a new tool, “diaboolical,” was created to support enhancing readability of the summary using a novel approach to bit-vector simplification that leverages logic synthesis and tautologies. For this effort, students at Auburn University created several finite state machines (FSMs) to serve as exemplars for high-consequence systems. Function summaries for each of the machines were obtained, and then portions of the summaries were simplified using both diaboolical and the simplification procedure of a popular SMT solver. A comparison of the results shows that diaboolical can often produce smaller function summaries, with expression length improvements over the unsimplified function summaries ranging from 0% to 90% for diaboolical and 0% to 65% for the SMT solver, though diaboolical had a significantly greater cost in time. Diaboolical was evaluated against a collection of “arbitrary” C-code as well as FSM exemplars, and for both datasets it achieved an approximately 10% improvement in expression length compared to simplifications that could be obtained using existing techniques. Function summaries can assist assurance efforts that evaluate existing systems and their executable code. A smaller function summary is likely easier for humans to understand and could thus increase the ability and efficacy of assurance practices centered around the analysis of executable artifacts.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

Abstract not provided.

This poster describes progress in training a variety of neural network architectures on a typical detection task. The poster will be presented at the Sandia MLDL Workshop

Goes over a simple software library (Python) for utilizing sparse autoencoders for more models than just language models.

Abstract not provided.

Abstract not provided.

The Spent Fuel & Waste Science and Technology (SFWST) Campaign of the Office of Spent Fuel & Waste Disposition of U.S. Department of Energy Office of Nuclear Energy (DOE-NE) is conducting research and development on geologic disposal of spent nuclear fuel (SNF) and high-level nuclear waste (HLW). This report describes fiscal year 2024 accomplishments in the Geologic Disposal Safety Assessment (GDSA) PFLOTRAN Development work package, which is charged with developing subsurface simulation software for postclosure performance assessment of deep geologic disposal of SNF and HLW.

Poster to be presented at the Applied Superconductivity Conference, 2024 in Salt Lake City, Utah in September. This poster details our calculations on the propagation of ballistic fluxon solitons in long Josephson junctions.

Presentation on content from this publication: https://www.sciencedirect.com/science/article/pii/S0952197623017992

Abstract not provided.

Radial basis functions (RBFs) are a powerful tool for constructing high-order accurate reduced representations of scattered data in arbitrary dimension and on manifolds. We present a method of constructing data approximations in which we utilize a functional tail to capture a global background profile and a RBF neural network (NN) to capture the smaller-scale features. In the RBF NN the RBF centers, matrix shape parameters were selected adaptively for each RBF. We also utilized a geodesic notion of distance on the manifold on which the data lies, e.g., the spherical geodesic for data on the sphere. Although each of these ideas have been been investigated separately in previous works, their combination into a single algorithm is novel. We defined a machine learning problem in which these properties are learned to minimize the data reduction error. We demonstrate the algorithm for applications of scattered data reduction in the plane and on the sphere.

Abstract not provided.

Abstract not provided.

The effects of raw water leaching at the U.S. SPR were modeled using SANSMIC. In addition to 18 caverns identified in previous leaching reports, six caverns have been identified for further monitoring based on the results of this report. SANSMIC continues to serve as a useful tool for monitoring changes in cavern shape. However, some caverns where string cuts were made in the middle of drawdowns are showing overprediction in the SANSMIC model when compared with recent sonars.

In this letter, we present interfacial fracture toughness data for a polymer-metal interface where tests were conducted at various test temperatures T and loading rates δ˙. An adhesively bonded asymmetric double cantilever beam (ADCB) specimen was utilized to measure toughness. ADCB specimens were created by bonding a thinner, upper adherend to a thicker, lower adherend (both 6061 T6 aluminum) using a thin layer of epoxy adhesive, such that the crack propagated along the interface between the thinner adherend and the epoxy layer. The specimens were tested at T from 25 to 65 °C and δ˙ from 0.002 to 0.2 mm/s. The measured interfacial toughness Γ increased as both T and δ˙ increased. For an ADCB specimen loaded at a constant δ˙, the energy release rate G increases as the crack length a increases. For this reason, we defined rate effects in terms of the rate of change in the energy release rate G˙. Although not rigorously correct, a formal application of time–temperature superposition (TTS) analysis to the Γ data provided useful insights on the observed dependencies. In the TTS-shifted data, Γ decreased and then increased for monotonically increasing G˙. Thus, the TTS analysis suggests that there is a minimum value of Γ. This minimum value could be used to define a lower bound in Γ when designing critical engineering applications that are subjected to T and δ˙ excursions.

Embedded machine-learned models (EMLMs) have the promise to improve the predictive accuracy of engineering simulators in environments of national interest. EMLMs often comprise complex input-output maps (e.g., neural networks), which make them unamenable to rigorous analysis and generally difficult to interpret. In the face of decades of theory, this lack of interpretability is a significant barrier to building confidence in these models. This work outlines an approach to interpret EMLMs using sparse polynomial regression for comparison with theoretical understanding. To do so, we build on the concept of Locally Interpretable Model-agnostic Explanations (LIME) using physics-informed clustering, prototype selection, and library construction. While general, we demonstrate our method on tensor-basis neural networks used in Reynolds-Averaged Navier-Stokes simulations of hypersonic fluid flows. Results are presented for a simulated toy model and for direct numerical simulations (DNS) of turbulent flows over a flat plate.