This year, we focused on completing the light squeezing and building the imaging station. In this report, we present a detailed description of a quantum imaging experiment utilizing squeezed light. The entire experimental setup has two parts, namely, the squeezing station where we produce quantum-noise squeezed light where a light quadrature (either the amplitude of the phase) has reduced quantum error below the shot noise of coherent light, and the imaging station where the squeezed light is used to image an object. The squeezing station consists of an optical parametric oscillator operating below the laser threshold. We provide the status quo and the plans for the squeezing imaging experiment.
The kinetic codes used to model the coupled dynamics of electromagnetic fields and charged particle transport have requirements for spatial, temporal, and charge resolution. These requirements may vary by the solution technique and scope of the problem. In this report, we investigate the resolution limits in the energy-conserving implicit particle-in-cell code CHICAGO. This report has the narrow aim of determining the maximum acceptable grid spacing for the dense plasmas generated in models of z-pinch target gases and power-flow electrode plasmas. In the 2D sample problem, the plasma drifts without external forces with velocity of 10 cm/µs. Simulations are scaled by plasma density to maintain uniform strides across the plasma and from the plasma to the boundaries. Additionally, the cloud-in-cell technique is used with 400 particles per cell and Δt = 0.85× the Courant limit. For the linear cloud distribution, the criterion for conserving energy is ΔE/Etot < 0.01 for 50,000 time steps. The grid resolution is determined to crudely be Δx ≲ 3ls, where ls is the electron collisionless skin depth. For the second-order cloud distribution the criterion is ΔE/Etot < 0.005 yielding Δx ≤ 15ls. These scalings are functions of the chosen vd, Δt, particles-per-cell, and number of steps.
Cyberattacks against industrial control systems have increased over the last decade, making it more critical than ever for system owners to have the tools necessary to understand the cyber resilience of their systems. However, existing tools are often qualitative, subject matter expertise-driven, or highly generic, making thorough, data-driven cyber resilience analysis challenging. The ADROC project proposed to develop a platform to enable efficient, repeatable, data-driven cyber resilience analysis for cyber-physical systems. The approach consists of two phases of modeling: computationally efficient math modeling and high-fidelity emulations. The first phase allows for scenarios of low concern to be quickly filtered out, conserving resources available for analysis. The second phase supports more detailed scenario analysis, which is more predictive of real-world systems. Data extracted from experiments is used to calculate cyber resilience metrics. ADROC then ranks scenarios based on these metrics, enabling prioritization of system resources to improve cyber resilience.
Machine learning-based data-driven modeling can allow computationally efficient time-dependent solutions of PDEs, such as those that describe subsurface multiphysical problems. In this work, our previous approach (Kadeethum et al., 2021d) of conditional generative adversarial networks (cGAN) developed for the solution of steady-state problems involving highly heterogeneous material properties is extended to time-dependent problems by adopting the concept of continuous cGAN (CcGAN). The CcGAN that can condition continuous variables is developed to incorporate the time domain through either element-wise addition or conditional batch normalization. Moreover, this framework can handle training data that contain different timestamps and then predict timestamps that do not exist in the training data. As a numerical example, the transient response of the coupled poroelastic process is studied in two different permeability fields: Zinn & Harvey transformation and a bimodal transformation. The proposed CcGAN uses heterogeneous permeability fields as input parameters while pressure and displacement fields over time are model output. Our results show that the model provides sufficient accuracy with computational speed-up. This robust framework will enable us to perform real-time reservoir management and robust uncertainty quantification in poroelastic problems.
Additive manufactured Ti-5Al-5V-5Mo-3Cr (Ti-5553) is being considered as an AM repair material for engineering applications because of its superior strength properties compared to other titanium alloys. Here, we describe the failure mechanisms observed through computed tomography, electron backscatter diffraction (EBSD), and scanning electron microscopy (SEM) of spall damage as a result of tensile failure in as-built and annealed Ti-5553. We also investigate the phase stability in native powder, as-built and annealed Ti-5553 through diamond anvil cell (DAC) and ramp compression experiments. We then explore the effect of tensile loading on a sample containing an interface between a Ti-6Al-V4 (Ti-64) baseplate and additively manufactured Ti-5553 layer. Post-mortem materials characterization showed spallation occurred in regions of initial porosity and the interface provides a nucleation site for spall damage below the spall strength of Ti-5553. Preliminary peridynamics modeling of the dynamic experiments is described. Finally, we discuss further development of Stochastic Parallel PARticle Kinteic Simulator (SPPARKS) Monte Carlo (MC) capabilities to include the integration of alpha (α)-phase and microstructural simulations for this multiphase titanium alloy.
X-ray stereo digital image correlation (DIC) measurements were performed at 10 kHz on the internal surface of a jointed structure in a shock tube at a shock Mach number of 1.42 and compared with optical stereo DIC measurements on the outer, visible surface of the structure. The shock tube environment introduces temperature and density gradients in the gas through which the structure was imaged, resulting in spatial and temporal index of refraction variations. These variations cause bias errors in optical DIC measurements due to beam-steering but have minimal influence on x-ray DIC measurements. These results demonstrate the utility of time-resolved x-ray DIC measurements in complicated environments where optical measurements suffer severe errors and/or are precluded by lack of optical access.
The frequency of unintended releases in a compressed natural gas system is an important aspect of the system quantitative risk assessment. The frequencies for possible release scenarios, along with engineering models, are utilized to quantify the risks for compressed natural gas facilities. This report documents component leakage frequencies representative of compressed natural gas components that were estimated as a function of the normalized leak size. A Bayesian statistical method was used which results in leak frequency distributions for each component which represent variation and uncertainty in the leak frequency. The analysis shows that there is high uncertainty in the estimated leak frequencies due to sparsity in compressed natural gas data. These leak frequencies may still be useful in compressed natural gas system risk assessments, as long as this high uncertainty is acknowledged and considered appropriately.
The time dependence of phase diagrams and how to model rate dependent transitions remains one of the key unanswered questions in physics. When a material is loaded dynamically through equilibrium phase boundaries, it is the kinetics that determines the real time expression of a phase transition. Here we report the atomic and nanosecond-scale quantification of kinetics of shock-driven phase transition in multiple materials. We uniquely make use of a both a simple shock as well as shock-and-hold loading pathways compress different crystalline solids and induce structural phase transitions below melt. Coupling shock loading with time-resolved synchrotron x-ray diffraction (DXRD), we probe the structural transformations of these solids in the short-lived high pressure and temperature states generated. The novelty and power of using DXRD for the assessment of kinetics of phase transitions lies in the ability to discover and identify new phases and to examine kinetics without prior knowledge of a material's phase diagram. Our results provide a quantified expression and a physics model of kinetics of formation of high-pressure phases under shock loading: transition incubation time, evolution, completion time and crystallization rate.
This report discusses the progress on the collaboration between Sandia National Laboratories (Sandia) and Japan Atomic Energy Agency (JAEA) on the sodium fire research in fiscal year (FY) 2022 and is a continuation of the FY 2021 progress report. We only report the changes made to the current sodium pool fire model in MELCOR. We modified and corrected many control functions to enhance the fraction of oxygen consumed that reacts to form monoxide (FO2) parameter in the current model from the FY2021 report. This year's enhancements relate to better agreement of the suspended aerosol measurement from JAEA's F7 series tests. Staff from Sandia and JAEA conducted the validation studies of the sodium pool fire model in MELCOR. To validate this pool fire model with the latest enhancement, JAEA sodium pool fire experiments (F7-1 and F7-2) were used. The results of the calculation, including the code-to-code comparisons are discussed as well as suggestions for further model improvement. Finally, recommendations are made for new MELCOR simulations for FY 2023.
We present the results of an LDRD project, funded by the Nuclear Deterrence IA, to develop capabilities for quantitative assessment of pyrotechnic thermal output. The thermal battery igniter is used as an exemplar system. Experimental methodologies for thermal output evaluation are demonstrated here, which can help designers and engineers better specify pyrotechnic components , provide thermal output guidelines for new formulations, and generate new metrics for assessing component performance and margin given a known failure condition. A heat-transfer analysis confirms that the dominant mode of energy transfer from the pyrotechnic output plume to the heat pellet is conduction via deposition of hot titanium particles. A simple lumped-parameter model of titanium particle heat transfer and a detailed multi-phase model of deposition heat transfer are discussed. Pyrotechnic function, as defined by "go/no-go" standoff testing of a heat pellet, is correlated with experimentally measured igniter plume temperature, titanium metal particle temperature, and energy deposition. Three high-speed thermal diagnostics were developed for this task. A three-color imaging pyrometer, acquiring 100k images per second on three color channels, is deployed for measurement of titanium particle temperatures. Complimentary measurements of the overall igniter plume emission ("color") temperature were conducted using a transmission-grating spectrograph in line-imaging mode. Heat flux and energy deposition to a cold wall at the heat-pellet location were estimated using an eroding thermocouple probe, with a frequency response of ~5 kHz. Ultimate "go/no-go" function in the igniter/heat-pellet system was correlated with quantitative thermal metrics, in particular surface energy deposition and plume color temperature. Titanium metal-particle and plume color temperatures both experience an upper bound approximated by the 3245-K boiling point of TiO2. Average metal-particle temperatures remained nearly constant for all standoff distances at T = 2850 K, ± 300 K, while plume color temperature and heat flux decay with standoff—suggesting that heat-pellet failure results from a drop in metal-particle flux and not particle temperature. At 50% likelihood of heat-pellet failure, peak time-resolved plume color temperatures drop well below TiO2 boiling to ~2000 - 2200 K, near the TiO2 melting point. Estimates of peak heat flux decline from up to 1 GW/m2 for near-field standoffs to below 320 MW/m2 at 50% failure likelihood.
High-enthalpy hypersonic flight represents an application space of significant concern within the current national-security landscape. The hypersonic environment is characterized by high-speed compressible fluid mechanics and complex reacting flow physics, which may present both thermal and chemical nonequilibrium effects. We report on the results of a three-year LDRD effort, funded by the Engineering Sciences Research Foundation (ESRF) investment area, which has been focused on the development and deployment of new high-speed thermochemical diagnostics capabilities for measurements in the high-enthalpy hypersonic environment posed by Sandia's free-piston shock tunnel. The project has additionally sponsored model development efforts, which have added thermal nonequilibrium modeling capabilities to Sandia codes for subsequent design of many of our shock-tunnel experiments. We have cultivated high-speed, chemically specific, laser-diagnostic approaches that are uniquely co-located with Sandia's high-enthalpy hypersonic test facilities. These tools include picosecond and nanosecond coherent anti-Stokes Raman scattering at 100-kHz rates for time-resolved thermometry, including thermal nonequilibrium conditions, and 100-kHz planar laser-induced fluorescence of nitric oxide for chemically specific imaging and velocimetry. Key results from this LDRD project have been documented in a number of journal submissions and conference proceedings, which are cited here. The body of this report is, therefore, concise and summarizes the key results of the project. The reader is directed toward these reference materials and appendices for more detailed discussions of the project results and findings.
Many materials of interest to Sandia transition from fluid to solid or have regions of both phases coexisting simultaneously. Currently there are, unfortunately, no material models that can accurately predict this material response. This is relevant to applications that "birth stress" related to geoscience, nuclear safety, manufacturing, energy production and bioscience. Accurately capturing solidification and residual stress enables fully predictive simulations of the evolving front shape or final product. Accurately resolving flow of proppants or blood could reduce environmental impact or lead to better treatments for heart attacks, thrombosis, or aneurism. We will address a science question in this proposal: When does residual stress develop during the critical transition from liquid to solid and how does it affect material deformation? Our hypothesis is that these early phases of stress development are critical to predictive simulation of material performance, net shape, and aging. In this project, we use advanced constitutive models with yield stress to represent both fluid and solid behavior simultaneously. The report provides an abbreviated description of the results from our LDRD "Stress Birth and Death: Disruptive Computational Mechanics and Novel Diagnostics for Fluid-to-Solid Transitions," since we have written four papers that document the work in detail and which we reference. We give highlights of the work and describe the gravitationally driven flow visualization experiment on a model yield stress fluid, Carbopol, at various concentrations and flow rates. We were able to collapse the data on a single master curve by showing it was self-similar. We also describe the Carbopol rheology and the constitutive equations of interest including the Bingham-Carreau-Yasuda model, the Saramito model, and the HB-Saramito model including parameter estimation for the shear and oscillatory rheology. We present several computational models including the 3D moving mesh simulations of both the Saramito models and Bingham-Carreau-Yasuda (BCY) model. We also show results from the BCY model using a 3D level set method and two different ways of handling reduced order Hele-Shaw modeling for generalized Newtonian fluids. We present some first ever two-dimensional results for the modified Jeffries Kamani-Donley-Rogers constitutive equation developed during this project. We include some recent results with a successful Saramito-level set coupling that allows us to tackle problems with complex geometries like mold filling in a thin gap with an obstacle, without the need for remeshing or remapping. We report on some experiments for curing systems where fluorescent particles are used to track material flow. These experiments were carried out in an oven on Sylgard 184 as a model polymerizing system. We conclude the report with a summary of accomplishments and some thoughts on follow-on work.
In recent years, infections and damage caused by malware have increased at exponential rates. At the same time, machine learning (ML) techniques have shown tremendous promise in many domains, often out performing human efforts by learning from large amounts of data. Results in the open literature suggest that ML is able to provide similar results for malware detection, achieving greater than 99% classifcation accuracy [49]. However, the same detection rates when applied in deployed settings have not been achieved. Malware is distinct from many other domains in which ML has shown success in that (1) it purposefully tries to hide, leading to noisy labels and (2) often its behavior is similar to benign software only differing in intent, among other complicating factors. This report details the reasons for the diffcultly of detecting novel malware by ML methods and offers solutions to improve the detection of novel malware.
Today as well as tomorrows spaceborne assets impact almost all areas of national and nuclear security. Spaceborne assets can not only collect and disseminate valuable data, well beyond just the visual, but also track terrestrial-based mobile assets in real-time, and active spaceborne platforms potentially pose serious risk to vulnerable earth-based systems and infrastructures. The capability to defend national spaceborne assets from attack/interference is critical for security interests. This effort supports this mission through the cost-effective preeminent detection of approaching threats to our nation’s vital resources, in order to help secure and trust these high-value assets against the threats of tomorrow. This project develops novel fabrication techniques for conformal, low-profile and lightweight leakywave antenna (LWA) detection/imaging systems, which fuses technical embroidery (TE) and laser ablation (LA) processes with LWA design. Technical embroidery is an emerging field in additive textile manufacturing where flexible materials and functionalized fabrics are created for a wide variety of uses and purposes, while laser ablation is the process of removing material from a solid surface by irradiating it with a laser beam. Here, thin, conformal antenna designs are designed, modeled and fabricated using both TE and LA, to create lightweight, flexible and conformal object detection and imaging radars. This novel development ensures our nation’s ability to field advanced lightweight and conformal technologies to protect spaceborne assets.
Oblapenko, G.; Goldstein, D.; Varghese, P.; Moore, C.
We propose a new scheme for simulation of collisions with multiple possible outcomes in variable-weight DSMC computations. The scheme is applied to a 0-D ionization rate coefficient computation, and 1-D electrical breakdown simulation. We show that the scheme offers a significant (up to an order of magnitude) improvement in the level of stochastic noise over the usual acceptance-rejection algorithm, even when controlling for the slight additional computational costs. The benefits and performance of the scheme are analyzed in detail, and possible extensions are proposed.
As part of the development process, scaled testing of wave energy converter devices are necessary to prove a concept, study hydrodynamics, and validate control system approaches. Creating a low-cost, small, lightweight data acquisition system suitable for scaled testing is often a barrier for wave energy converter developers’ ability to test such devices. This paper outlines an open-source solution to these issues, which can be customized based on specific needs. This will help developers with limited resources along a path toward commercialization.
X-ray computed tomography is generally a primary step in characterization of defective electronic components, but is generally too slow to screen large lots of components. Super-resolution imaging approaches, in which higher-resolution data is inferred from lower-resolution images, have the potential to substantially reduce collection times for data volumes accessible via x-ray computed tomography. Here we seek to advance existing two-dimensional super-resolution approaches directly to three-dimensional computed tomography data. Multiple scan resolutions over a half order of magnitude of resolution were collected for four classes of commercial electronic components to serve as training data for a deep-learning, super-resolution network. A modular python framework for three-dimensional super-resolution of computed tomography data has been developed and trained over multiple classes of electronic components. Initial training and testing demonstrate the vast promise for these approaches, which have the potential for more than an order of magnitude reduction in collection time for electronic component screening.
Milestone accomplishments staged the ExaWind team for successful completion of KPP-2 challenge problem in FY23, which requires the simulation on Frontier of at least four MW-scale turbines in an atmospheric boundary layer with at least 20B gridpoints. The ExaWind project and software stack is many faceted, with team members working on multiple areas, including linear-system solvers (Trilinos, hypre, AMReX), overset meshes, turbulence modeling, and in situ visualization, all with an aim for high fidelity predictions and performance portability. This milestone marks significant improvements on many fronts and provides the team with a pathway to exascale wind farm simulations in FY23.
Bidadi, Shreyas B.; Brazell, Michael B.; Brunhart-Lupo, Nicholas B.; Henry de Frahan, Marc T.; Lee, Dong H.; Hu, Jonathan J.; Melvin, Jeremy M.; Mullowney, Paul M.; Vijayakumar, Ganesh V.; Moser, Robert D.; Rood, Jon R.; Sakievich, Philip S.; Sharma, Ashesh S.; Williams, Alan B.; Sprague, Michael A.
The goal of the ExaWind project is to enable predictive simulations of wind farms comprised of many megawatt-scale turbines situated in complex terrain. Predictive simulations will require computational fluid dynamics (CFD) simulations for which the mesh resolves the geometry of the turbines, capturing the thin boundary layers, and captures the rotation and large deflections of blades. Whereas such simulations for a single turbine are arguably petascale class, multi-turbine wind farm simulations will require exascale-class resources.
The Sandia National Laboratories, California (SNL/CA) site comprises approximately 410 acres and is located in the eastern portion of Livermore, Alameda County, California. The property is owned by the United States Department of Energy and is being managed and operated by National Technology & Engineering Solutions of Sandia, LLC. The facility location is shown on the Site Map(s) in Appendix A. This Stormwater Pollution Prevention Plan (SWPPP) is designed to comply with California’s General Permit for Stormwater Discharges Associated with Industrial Activities (General Permit) Order No. 2015-0122-DWQ (NPDES No. CAS000001) issued by the State Water Resources Control Board (State Water Board) (Ref. 6.1). This SWPPP has been prepared following the SWPPP Template provided on the California Stormwater Quality Association Stormwater Best Management Practice Handbook Portal: Industrial and Commercial (CASQA 2014). In accordance with the General Permit, Section X.A, this SWPPP contains the following required elements: Facility Name and Contact Information; Site Map; List of Significant Industrial Materials; Description of Potential Pollution Sources; Assessment of Potential Pollutant Sources; Minimum BMPs; Advanced BMPs, if applicable; Monitoring Implementation Plan (MIP); Annual Comprehensive Facility Compliance Evaluation (Annual Evaluation); and, Date that SWPPP was Initially Prepared and the Date of Each SWPPP Amendment, if Applicable.
Disastrous consequences can result from defects in manufactured parts—particularly the high consequence parts developed at Sandia. Identifying flaws in as-built parts can be done with nondestructive means, such as X-ray Computed Tomography (CT). However, due to artifacts and complex imagery, the task of analyzing the CT images falls to humans. Human analysis is inherently unreproducible, unscalable, and can easily miss subtle flaws. We hypothesized that deep learning methods could improve defect identification, increase the number of parts that can effectively be analyzed, and do it in a reproducible manner. We pursued two methods: 1) generating a defect-free version of a scan and looking for differences (PandaNet), and 2) using pre-trained models to develop a statistical model of normality (Feature-based Anomaly Detection System: FADS). Both PandaNet and FADS provide good results, are scalable, and can identify anomalies in imagery. In particular, FADS enables zero-shot (training-free) identification of defects for minimal computational cost and expert time. It significantly outperforms prior approaches in computational cost while achieving comparable results. FADS’ core concept has also shown utility beyond anomaly detection by providing feature extraction for downstream tasks.
Single photon detection (SPD) plays an important role in many forefront areas of fundamental science and advanced engineering applications. In recent years, rapid developments in superconducting quantum computation, quantum key distribution, and quantum sensing call for SPD in the microwave frequency range. We have explored in this LDRD project a new approach to SPD in an effort to provide deterministic photon-number-resolving capability by using topological Josephson junction structures. In this SAND report, we will present results from our experimental studies of microwave response and theoretical simulations of microwave photon number resolving detector in topological Dirac semimetal Cd3As2. These results are promising for SPD at the microwave frequencies using topological quantum materials.
Projection-based model order reduction allows for the parsimonious representation of full order models (FOMs), typically obtained through the discretization of a set of partial differential equations (PDEs) using conventional techniques (e.g., finite element, finite volume, finite difference methods) where the discretization may contain a very large number of degrees of freedom. As a result of this more compact representation, the resulting projection-based reduced order models (ROMs) can achieve considerable computational speedups, which are especially useful in real-time or multi-query analyses. One known deficiency of projection-based ROMs is that they can suffer from a lack of robustness, stability and accuracy, especially in the predictive regime, which ultimately limits their useful application. Another research gap that has prevented the widespread adoption of ROMs within the modeling and simulation community is the lack of theoretical and algorithmic foundations necessary for the “plug-and-play” integration of these models into existing multi-scale and multi-physics frameworks. This paper describes a new methodology that has the potential to address both of the aforementioned deficiencies by coupling projection-based ROMs with each other as well as with conventional FOMs by means of the Schwarz alternating method [41]. Leveraging recent work that adapted the Schwarz alternating method to enable consistent and concurrent multiscale coupling of finite element FOMs in solid mechanics [35, 36], we present a new extension of the Schwarz framework that enables FOM-ROM and ROM-ROM coupling, following a domain decomposition of the physical geometry on which a PDE is posed. In order to maintain efficiency and achieve computation speed-ups, we employ hyper-reduction via the Energy-Conserving Sampling and Weighting (ECSW) approach [13]. We evaluate the proposed coupling approach in the reproductive as well as in the predictive regime on a canonical test case that involves the dynamic propagation of a traveling wave in a nonlinear hyper-elastic material.