SEARCH · Engineering Papers
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
This report is a major update of the report ”Cubical moderators for the Neutron Target Demonstator (NTD), LA-UR-23-33899”. It serves as a documentation of milestone M1 of the LDRD project 20240004DR.
This is the poster presentation prepared for 2024 Geothermal Rising Conference, October 27-30, 2024, Waikoloa, HI. The poster is part of the LDRD project CS-GE^2 to demonstrate the multiple-tier fracture imaging analysis for geothermal reservoir applications.
This report summarizes the research and capabilities developed as part of the project “Analysis and Validation of Critical Infrastructure Models using Model Order Reduction” (CI-MOR) LDRD project. CI-MOR research enables the solution of large, complex optimization models that naturally arise in national security challenges involving critical infrastructures. Specifically, CI-MOR researchers developed methods to (1) rigorously approximate complex, nonlinear optimization formulations, (2) identify alternative near-optimal solutions, (3) accelerate optimization workflows used for complex applications, and (4) rigorously integrate domain knowledge in stochastic-process models. This report provides an overview of the research done in CI-MOR, and we describe application exemplars used to illustrate CI-MOR capabilities. Furthermore, we describe the software developed by CI-MOR that researchers can leverage to analyze new applications.
Development of a radioimaging diagnostic for high-voltage component reliability testing and electrical breakdown computational model validation is described. Radioimaging has its roots in radio astronomy, where aperture synthesis (also known as synthesis imaging) has been utilized for decades to image radio sources far from Earth. Radioimaging as described herein, in contrast, seeks to image radio sources in close proximity to its receivers (i.e., in a laboratory environment). Here it is shown that corona discharge, a non-destructive precursor to catastrophic (thermal) arc discharge, electromagnetically radiates strongly within a 250 kHz – 2.5 GHz bandwidth, and is readily detected and located by postprocessing the received radio signals. The ability of radioimaging to detect both corona and arc discharge (grouped together herein as high voltage breakdown or HVB) makes it a valuable tool for 100% HVB detection in materials, components, and devices, and has the ability to indicate electrical weakness (via corona detection) prior to a destructive arc discharge event. Radioimaging enables HVB to be located both internal and external to dielectric components under test in near-real-time, with multiple and/or extended HVB events located simultaneously. In contrast, existing non-destructive diagnostics (at the time of this writing) either indicate electrical breakdown without resolving failure locations (e.g., current, voltage, and chemical measurements), locate external HVB (e.g., high-speed optical and ultraviolet (UV) measurements or photography), or locate both external and internal HVB but with low fidelity (e.g., a single HVB source can be located by existing time-of-arrival (TOA) UHF or acoustic emissions). Radioimaging instead creates a sequence of high-fidelity images similar to an optical high-speed camera but at radiofrequencies (RF), and is not limited to two-dimensions. Moreover, radioimaging has already served one internal and two external industry customers, the results of which are detailed in this report. The radioimaging results described herein were part of a three-year effort funded by the Sandia Lab Directed Research and Development (LDRD) program within the Radiation, Electromagnetic, High Energy Density Science (REHEDS) investment area.
We discuss here the possibility of using the Isotope Production Facility (IPF) at LANSCE to test the feasibility and performance of a graphite cube to create a neutron target. The idea of a neutron target enables the measurement of neutron induced reactions in inverse kinematics. This idea is part of the LANSCE strategy to stay a worldwide leader for neutron-induced research. The proof of-principle is the core of the LDRD project 20240004DR. The first steps are a series of activation experiments with different neutron energy distributions. Figs. 1 and 2 show a graphite cube built and simulated from pieces available at LANSCE. First activation experiments have already been performed at lower neutron energies utilizing (p,n) reactions on Li and Be.
The Eastern Pacific Cloud Aerosol Precipitation Experiment (EPCAPE) was aimed to enhance the understanding of cloud and aerosol properties in the region surrounding La Jolla, California. To address challenges in data collection and processing from various instruments, an edge computing device known as Waggle Sage Node (WSN) was deployed at the Ellen Browning Scripps Memorial Pier. WSN is a distributed-sensing platform designed to collect and analyze environmental data at the edge. Sage is a multi-agency-supported project that designs and builds a new kind of national-scale reusable cyberinfrastructure to enable artificial intelligence (AI) at the edge based on the Waggle platform. Sponsors include the U.S. Department of Energy (DOE) Advanced Scientific Computing Research (ASCR), DOE National Nuclear Security Administration (NNSA), DOE Biological and Environmental Research (BER) through DOE Artificial Intelligence for Earth System Predictability (AI4ESP), Argonne Laboratory-Directed Research and Development (LDRD). Sage (https://sagecontinuum.org/) is funded as a National Science Foundation Mid-Scale Research Infrastructure (MSRI) project (https://www.nsf.gov/awardsearch/showAward?AWD_ID=1935984). This robust, multi-architecture edge computing platform facilitated environmental monitoring during the campaign. This report details the scientific objectives, deployment process, and key results of integrating Waggle into the EPCAPE field campaign.
This study investigates the simulation of marine stratocumulus (stratoCu) clouds – common in many offshore and coastal wind energy regions - using Lagrangian Superdroplet Method (SDM) implemented in Energy Research and Forecasting (ERF) model. Developed by the Wind Energy Technologies Office (WETO), ERF is a high-resolution atmospheric modeling tool that relies on traditional bulk microphysics schemes. Under an LLNL-led Laboratory Directed Research and Development (LDRD) project, SDM was integrated into ERF to improve cloud process fidelity by explicitly tracking individual “superdroplets”. This enables a realistic representation of the droplet size spectrum, cloud-atmosphere interactions, and potential feedbacks from wind plants on atmospheric flow and cloud structure.
Flowable carbon slurry electrodes promise continuous, scalable electrochemical desalination and energy storage when their rheology (flow and stability) and electronic connectivity are jointly optimized. This report integrates six complementary thrusts that we executed under the Laboratory Directed Research and Development (LDRD) project: (i) a fast, ex situ, centrifuge-compatible method to quantify pellet (packed slurry) electronic conductivity, with a validated cell constant (“shape factor”) and quantitative structure–property insights for carbon mixtures; (ii) a controlled circulation study showing that peristaltic pumping can irreversibly alter particle size and morphology over laboratory-relevant durations; (iii) development of a miniature jet pump that preserves particle integrity while achieving entrainment ratios of 2–4 at practical flow rates; (iv) a comprehensive mapping of how a representative non-ionic surfactant (Tween 20) collapses the carbon gel network, separating network-percolation losses from intrinsic particle contact resistance; (v) in-situ visualization revealing how gravitational settling and flow-driven mixing compete to control charge transport, with dynamic bed layers preventing saturation while achieving 30-fold local concentration enhancement; and (vi) a simple, widely adoptable fabrication route to structured, gas-nitrided titanium nitride (TiN) titanium current collectors that dramatically lower interfacial losses and enable geometry-driven performance gains. Together, these results provide a coherent toolchain and design rules to accelerate formulation screening, de-risk device prototyping, and standardize characterization for carbon slurry electrodes in flow-electrode capacitive deionization (FCDI), electrochemical flow capacitor (EFC), and related systems.
Recently, randomized algorithms in numerical linear algebra, specifically those centered around random sketching, have gained traction in primarily theoretical research due to their potential to significantly reduce problem dimensionality at the cost of an O(1) multiplicative distortion factor. It has been assumed that this sketching can be done efficiently, but thorough investigation into how precisely to do it has been neglected. Moreover, the theory-based community has argued for sketching’s ability to reduce computational cost via complexity analysis, but has not researched how it affects the stability of the algorithms. At Sandia, efficient linear solvers that scale well on modern HPC architectures while maintaining stability are imperative for practical applications. In this LDRD, we developed a random sketching strategy that is substantially faster than existing ones, and demonstrate its superior performance in practice on a NVIDIA H100 GPU. Moreover, we show how this can be used to significantly outperform existing linear least squares solvers while improving the solver’s stability as well. Additionally, we demonstrate how this sketching strategy can be used to make a fast, stable QR factorization that can subsequently be used in s-step and block Krylov solvers. Finally, we incorporate a sketching-based block orthogonalization scheme into s-step GMRES, which is stable and faster than existing approaches on the Perlmutter supercomputer.
Atomistic defects play a pivotal role in functional and structural materials’ performance across a myriad of technology applications. Quantitative prediction of the thermodynamics and kinetics of defect formation and migration, respectively, typically requires accurate but expensive first-principles approaches, such as density functional theory (DFT). Their computational expense limits the throughput needed to perform high-throughput materials discovery/screening exercises or to perform materials modeling tasks relying on extensive sampling techniques. Therefore, in this Sandia National Laboratories Laboratory Directed Research and Development (LDRD) project (Project #229366), we developed a variety of machine learning techniques, trained on density functional theory calculations, to accelerate the discovery and modeling of materials in which vacancy and interstitial defects primarily dictate material performance. These include applications such as metal oxides for water-splitting or mixed ionic-electronic conduction, metal hydrides for hydrogen storage, and transition metal dichalcogenides for electronics, and the approaches developed herein can further be applied to many other domains that similarly depend on materials’ thermodynamic and kinetic defect properties for their desired functionality.
This SAND report summarizes work supported by an Engineering Sciences Research Foundation (ESRF) Lab Directed Research and Development (LDRD) project entitled “Beyond Magic Barrels: Digital manufacturing for crystallization, process development and optimization of explosive materials.” This SAND report is written in two parts with Part 1 discusses recrystallization of our explosive exemplar and Part 2 summarizing our work with recrystallization of resveratrol. We have studied resveratrol recrystallization with a multiscale approach combining experiments, modeling and simulation. At the single crystal scale, microscopy experiments illuminate crystal time-dependent growth rates using advanced image analysis. Bench scale experiments were carried out to look at growth of multiple particles in a small reactor creating thousands of particles and analyzing the results with microscopy and μCT. For the modeling we combine kinetic Monte Carlo (kMC) models with subscale information from density functional theory (DFT) or molecular dynamics. This work is discussed in Part 1 and can also be found in a paper from the project discussing a coarse-grained kMC model specifically developed for resveratrol. For well-mixed systems, we have population balance equations (PBE) linked with species mass conservation forming a set of ordinary differential equations that can be solved quickly. For more complicated geometries, such as the vat crystallization used throughout the complex, a coupled computational fluid dynamic (CFD)/PBE method was developed to account for gradients in temperature and concentration and differences in crystallization rates throughout the domain. These simulations are more complex and require high performance computing. We present results for two cases: 5% seed fast cool with parameters fit to the well-mixed case and 5% seed slow cool using the same parameters. We show reasonable agreement with experiments though are particles are significantly larger than the experiments.
The goal of this one-year LDRD is to study novel experimental approaches to provide benchmark spectral opacity measurements at the most extreme conditions beyond existing measurements (i.e. 2x temperature, ~5x density). These conditions correspond to those at half the solar radius and would provide a unique and novel test for stellar physics. These conditions are also relevant to advance stockpile stewardship research. The basic idea is to embed a sample inside the z-pinch dynamic hohlraum on the Z facility and obtain its opacity from either 1) its self-emission or 2) its absorption from using a capsule backlighter. This study evaluates the feasibility, limitations to these two approaches from a design approach. We found that emission-based perspective gives, at best, an opacity accuracy of about 20% for the current studied pinch conditions due to the limit of LTE approximation and radiation transport effects. The capsule-backlighter-based approach overcomes this limitation although some specific design questions remain to be answered. Promising designs are proposed for follow-on Z experiments. Future designs would help maintain Z relevancy in relation to ignition-enabled platforms at NIF.
LDRD Project 229427 aimed to determine how electrothermal instability (ETI) driven heating on a z-pinch rod pulsed with intense current evolves under mixed magnetic field (azimuthal + axial) conditions, which is pertinent to pulsed-power-driven magnetically-insulated transmission lines and physics targets. Experiments focused on diagnosing ETI-driven heating from deliberately-machined and well-characterized micron-scale surface defects (referred to as engineered defects or ED). Prior to the start of this project, understanding of how unmagnetized (B z =0) ED evolve had been obtained—simulations largely reproduce the experimentally observed high temperature spots which develop at the poles of bare/uncoated ED. Project 229427 extended the Mykonos Facility ED experimental platform to include axial field. In the first class of experiments, axial field was provided “dynamically” via a helical return can (HRC). In this case, B z and B θ rise at the same rate. Generally, the HRC generated magnetic field at a fixed polarization angle Φ B =arctan(B z /B θ )=15° on the rod's surface. In the second class of experiments, axial field was provided “statically” via a slow-rising (millisecond) external Helmholtz coil pair. In this case, B z was effectively constant/static throughout the 100 ns rise of the Mykonos current. For either case, a primary goal was to determine whether ETI provides a helical seed perturbation for the subsequent growth of the helical magneto Rayleigh-Taylor modes observed in MagLIF (static B z ) and dynamic screw pinch (DSP, dynamic B z ) experiments. When dynamic field was applied using an HRC, emissions from individual ED aligned toward Φ B , while emissions from ED within pairs elongated and preferentially merged along Φ B . These data strongly support that for a randomized defect distribution, heating from nearby current-density perturbations will favorably merge about Φ B to generate an extended seed perturbation that aligns toward the surface-field polarization, and this may impact the orientation of subsequent MRT growth on imploding liners. The results from the static field experiments were largely inconclusive, as any ETI heating rotation, if present, was obscured/overwhelmed by local/random heating from ED rim imperfections.
Predictions from physics-based computational models provide critical information to inform high consequence decisions, e.g., engineering design decisions. The ability to assess the reliability of such predictions is therefore critical. However, to date, reliability assessment rely heavily on expert judgment and qualitative arguments. This report details the efforts of LDRD 233072 to develop quantitative methods to assess reliability of model predictions, especially in the context of simplifying assumptions that can impact their reliability.
This presentation is about the recent progress of my LDRD (233082) about modeling-based assessment of multiple salt performance for underground hydrogen storage, and my talk will be at the CouFrac 2024 conference, Kyoto in Japan.
A conference presentation (oral talk for ablation workshop) of LDRD results
This report summarizes key accomplishments from a 2 year Laboratory Directed Research and Development (LDRD) project focused on advancing nanoscale characterization of metallic particulates and synthesis of laboratory-scale analogs. We successfully developed and refined nanoscale crystallography techniques with significantly improved spatial resolution and throughput compared to traditional transmission electron microscopy. We also established nanotomography capabilities for 3D particle reconstruction and began optimizing throughput for broader application. Additionally, we demonstrated the ability to relocate and characterize individual particles using atomic force microscopy (AFM), providing insights into hardness, conductivity, and magnetic behavior. On the synthesis front, we leveraged controlled and combustion-based methods to investigate annealing effects and rapid oxidation behavior of U–Fe–Ni compounds. These efforts collectively enhance our ability to characterize and interpret particulate formation and transformation under relevant conditions.
The standoff problem in pulsed-power driven IFE (Inertial Fusion Energy) refers to the problem of electrically connecting a fusion target with the driver in a way that preserves the driver/target interface in the presence of a large fusion energy release (100’s of megajoules to 1000 megajoules). High yield fusion drivers for stockpile stewardship applications have similar concerns, but without the complication of the high repetition rate needed for energy production. All proposed IFE or high yield systems have the challenge of protecting the reactor vessel and driver from the energy release, but pulsed power drivers have the additional challenge of establishing an electrical connection in a way that does not create excessive debris or require costly/massive structures that must be expendable. Fortunately, there are conceptual solutions in the form of very low mass transmission lines, usually in the form of wires or foils, since very little mass is needed to overcome the magnetic forces driving the transmission lines apart or provide a low resistance pathway given the short pulse duration of the driver. Several important aspects of the standoff problem were analyzed in the course of this LDRD project. Conceptual means of introducing the low mass transmission lines into the reactor chamber were proposed, including analysis of methods to retain the power flow gap between the electrodes in the presence of the chamber environment. The energy losses due to ohmic dissipation and magnetic acceleration of the low mass transmission lines were evaluated for candidate driver pulses and different transmission line masses, and the fluence of x-ray and neutron energy on the persistent power flow structures that must survive the pulse of energy were evaluated. Finally, means of rapidly pumping down the chamber were explored as a way to return the system to low pressure following the multi-atmosphere post-pulse pressure rise. The overall conclusion from the analysis is that, while the reactor environment creates stressful conditions for a low mass transmission line standoff system, there are concepts expending 10 kg or less of electrode mass each pulse with the potential to mitigate the stresses and successfully drive fusion targets. Whether such a system could lead to a practical and economical fusion reactor would require extensive research and development beyond the scope of the feasibility study.