Publications

With the amount of neuromorphic tools and frame-works growing in number, we recognize a need to increase interoperability within our field. As an illustration of this, we explore linking two independently constructed tools. Specifically, we detail the construction of an a execution backend based on STACS: Simulation Tool for Asynchronous Cortical Streams for the Fugu spiking neural algorithms framework. STACS extends the computational scope of Fugu, enabling fast simulation of large-scale neural networks. Combining these two tools is shown to be mutually beneficial, ultimately enabling more functionality than either tool on its own. We discuss design considerations, in-cluding recognizing the advantages of straightforward standards. Further, we provide some benchmark results showing drastic improvements in execution time.

The Rotor Aerodynamics, Aeroelastics, and Wake (RAAW) project's main objective was collecting data for validation of aerodynamic and aeroelastic codes for large, flexible rotors. These data come from scanning lidars of the inflow and wake, met tower, profiling lidar, blade deflection from photogrammetry, turbine SCADA data (including root bending loads), and hub-mounted SpinnerLidar inflow measurements. The goal of the present work is to analyze various methods to align the SpinnerLidar inflow data in time and space with individual blade loading. These methods would prove a way of analyzing turbine response while estimating the flowfield at each blade and provide a way of improving turbine response understanding using field data in real time, not just from simulations. The hub-mounted SpinnerLidar measures the inflow in the rotor frame meaning the locations of the blades relative to the measurement pattern do not change. The present work outlines some methods for correlating the SpinnerLidar inflow measurements with root bending loads in the rotor frame of reference accounting for both changes in wind speed and rotor speed from the measurement location one diameter upstream to each blade.

In order to characterize and benchmark computational hardware, software, and algorithms, it is essential to have many problem instances on-hand. This is no less true for quantum computation, where a large collection of real-world problem instances would allow for benchmarking studies that in turn help to improve both algorithms and hardware designs. To this end, here we present a large dataset of qubit-based quantum Hamiltonians. The dataset, called HamLib (for Hamiltonian Library), is freely available online and contains problem sizes ranging from 2 to 1000 qubits. HamLib includes problem instances of the Heisenberg model, Fermi-Hubbard model, Bose-Hubbard model, molecular electronic structure, molecular vibrational structure, MaxCut, Max-k-SAT, Max-k-Cut, QMaxCut, and the traveling salesperson problem. The goals of this effort are (a) to save researchers time by eliminating the need to prepare problem instances and map them to qubit representations, (b) to allow for more thorough tests of new algorithms and hardware, and (c) to allow for reproducibility and standardization across research studies.

The Rydberg dipole blockade has emerged as the standard mechanism to induce entanglement between neutral-Atom qubits. In these protocols, laser fields that couple qubit states to Rydberg states are modulated to implement entangling gates. Here we present an alternative protocol to implement entangling gates via Rydberg dressing and a microwave-field-driven spin-flip blockade [Y.-Y. Jau, Nat. Phys. 12, 71 (2016)1745-247310.1038/nphys3487]. We consider the specific example of qubits encoded in the clock states of cesium. An auxiliary hyperfine state is optically dressed so that it acquires partial Rydberg character. It thus acts as a proxy Rydberg state, with a nonlinear light shift that plays the role of blockade strength. A microwave-frequency field coupling a qubit state to this dressed auxiliary state can be modulated to implement entangling gates. Logic gate protocols designed for the optical regime can be imported to this microwave regime, for which experimental control methods are more robust. We show that unlike the strong dipole-blockade regime usually employed in Rydberg experiments, going to a moderate-spin-flip-blockade regime results in faster gates and smaller Rydberg decay. We study various regimes of operations that can yield high-fidelity two-qubit entangling gates and characterize their analytical behavior. In addition to the inherent robustness of microwave control, we can design these gates to be more robust to laser amplitude and frequency noises at the cost of a small increase in Rydberg decay.

Vapor-deposited PETN films undergo significant microstructure evolution when exposed to elevated temperatures, even for short periods of time. This accelerated aging impacts initiation behavior and can lead to chemical changes as well. In this study, as-deposited and aged PETN films are characterized using scanning electron microscopy and ultra-high performance liquid chromatography and compared with changes in initiation behavior measured via a high-throughput experimental platform that uses laser-driven flyers to sequentially impact an array of small explosive samples. Accelerated aging leads to rapid coarsening of the grain structure. At longer times, little additional coarsening is evident, but the distribution of porosity continues to evolve. These changes in microstructure correspond to shifts in the initiation threshold and onset of reactions to higher flyer impact velocities.

Multi-axis testing has become a popular test method because it provides a more realistic simulation of a field environment when compared to traditional vibration testing. However, field data may not be available to derive the multi-axis environment. This means that methods are needed to generate “virtual field data” that can be used in place of measured field data. Transfer path analysis (TPA) has been suggested as a method to do this since it can be used to estimate the excitation forces on a legacy system and then apply these forces to a new system to generate virtual field data. This chapter will provide a review of using TPA methods to do this. It will include a brief background on TPA, discuss the benefits of using TPA to compute virtual field data, and delve into the areas for future work that could make TPA more useful in this application.

This paper introduces a method for detecting stealthy false data injection attacks on the sensors of state of charge estimation algorithms used in battery management systems (BMSs). This method is based on sensor encoding, which is the active modification of sensor data streams. This method implements low-cost verification of the integrity of measurement data, allowing for the detection of stealthy additive attack vectors. It is considered that these attacks are crafted by malicious actors with knowledge of system models and who are capable of tampering with any number of measurements. The solution involves encoding all vulnerable measurements. The effectiveness of the method is demonstrated by simulations, where a stealthy attack on an encoded measurement vector captured by a BMS generates large residuals that trigger a chi-squared anomaly detector. Within the context of a defense-in-depth strategy, this method can be combined with other cybersecurity controls, such as encryption of data-in-transit, to equip cyberphysical systems with an additional line of defense against cyberattacks.

Significant research and development efforts are underway to ensure the cybersecurity of critical infrastructure and other operational technology (OT) systems. Given the high demand for safety and availability of OT systems, novel OT tools and systems must be designed and tested in a consequence-aware environment. This requires the seamless integration of high-fidelity physics simulations with emulations of OT devices and networks. This paper introduces a modular simulation environment for an oil and gas compressor station and its local networks. This environment will be used to support the development of novel tools to predict, detect, and mitigate cyber threats on critical infrastructure. This paper also describes future plans to expand the scale of the environment and its use cases.

Abstract not provided.

Measurements of the oxidation rates of various forms of carbon (soot, graphite, coal char) have often shown an unexplained attenuation with increasing temperatures in the vicinity of 2000 K, even when accounting for diffusional transport limitations and gas-phase chemical effects (e.g. CO2 dissociation). With the development of oxy-fuel combustion approaches for pulverized coal utilization with carbon capture, high particle temperatures are readily achieved in sufficiently oxygen-enriched environments. In this work, a new semi-global intrinsic kinetics model for high temperature carbon oxidation is created by starting with a previously developed 5-step mechanism that was shown to reproduce all major known trends in carbon oxidation, except for its high temperature kinetic falloff, and incorporating a recently discovered surface oxide decomposition step. The predictions of this new model are benchmarked by deploying the kinetic model in a steady-state reacting particle code (SKIPPY) and comparing the simulated results against a carefully measured set of pulverized coal char combustion temperature measurements over a wide range of oxygen concentrations in N2 and CO2 environments. The results show that the inclusion of the spontaneous surface oxide decomposition reaction step significantly improves predictions at high particle temperatures. Furthermore, the simulations reveal that O atoms released from the oxide decomposition step enhance the radical pool in the near-surface region and within the particle interior itself. Incorporation of literature rates for O and OH reactions with the carbon surface results in a reduction in the predicted radical pool concentrations and a very minor enhancement of the overall carbon oxidation rate.

Intermolecular Coulombic decay (ICD) in liquid water is a relatively novel type of nonlocal electronic decay mechanism, competing with the traditional mechanism of proton transfer between neighboring water molecules. Key features of ICD are its ultrafast non-radiative decay process and ultralong-range for excess energy transfer from the excited atom/molecule to its neighbors. Since detecting unambiguous ICD signatures in bulk liquid water is technically challenging, small water clusters have often been utilized to gain insights into ICD and other ionization processes in aqueous environment. Here, we present results from quantum mechanical calculations of the electronic structures of neutral to multiply-ionized water monomer, dimer, trimer, and tetramer. Core-level electrons of water are also considered here since recent studies demonstrated that emission site and energy of the electrons released during resonant-Auger-ICD cascade can be controlled by coupling ICD to resonant core excitation. Previous studies of ICD and electronic structures of neutral and ionized small water clusters and liquid water are briefly discussed.

A mesoscale model for the shock initiation of pentaerythritol tetranitrate (PETN) films has been utilized to elucidate changes in initiation thresholds due to aging conditions and surface roughness, as has been observed from a series of high-throughput initiation (HTI) experiments. The HTI experiment has generated a wealth of thin-pulse, sub-millimeter shock initiation data for vapor deposited PETN films with thicknesses of 67-125 μm and varying accelerated aging conditions. This is because the HTI experiment provides access to growth-to-detonation information for explosives that exhibit a shock-to-detonation transition (SDT) with length and time scales that are too short to be resolved by conventional experiments. Mesoscale modeling results using experimentally characterized PETN microstructures are able to capture the general trend observed in experiments, in that increasing flyer impact velocity increases reactions until full detonation is reached. Moreover, the varying degrees of surface roughness that were considered were found to provide only minor variances in the peak particle velocity at the explosive output. The model did not predict a shift in the initiation threshold due to aged microstructures alone, indicating that additional mesoscale model improvements are necessary.

Glass wedges are used increase the dimensionality of various optical measurements. Light refracted through the wedges can be focused to closely spaced points, lines or planes as shown in the applications herein.

Structural materials used in combustion or power generation systems need to have both environmental and temperature resistance to ensure long-term performance. As the energy sector transitions to hydrogen, there is a need to ensure compatibility of highly-alloyed austenitic steels and nickel-based alloys with hydrogen over a range of temperatures. Hydrogen embrittlement of these alloy systems is often considered most detrimental near ambient temperatures and low temperatures, although there is some evidence in the literature that hydrogen can affect creep behavior at elevated temperature. In the intermediate temperature range (e.g., 100-400C), it is uncertain whether hydrogen degradation of mechanical properties will be of concern. In this study, three alloys (304L, IN625, Hastelloy X) commonly used in power generation systems were thermally precharged with hydrogen and subsequently tensile tested to failure in air at temperatures ranging from 20°C to 200°C. At 20°C, the hydrogen-precharged condition for all materials exhibited loss in ductility with relative reduction of area ranging between 32% and 57%. The three alloys exhibited different trends with temperature but, in general, the relative reduction of area improved with increasing temperature tending towards noncharged behavior. Tests were performed at a nominal strain rate of 2 x 10-3 s-1 in order to minimize loss of hydrogen during elevated temperature testing. Hydrogen contents from the grip sections were measured both before and after testing and remained within 10% of starting content for 100°C tests and within 8-23% for 200°C tests.

Engineers are interested in the ability to compare dynamic environments for many reasons. Current methods of comparing environments compare the measured acceleration at the same physical point via a direct measurement during the two environments. Comparing the acceleration at a defined point only provides a comparison of response at that location. However, the stress and strain of the structure are defined by the global response of all the points in a structure. This chapter uses modal filtering to transform a set of measurements at physical degrees of freedom into modal degrees of freedom that quantify the global response of the structure. Once the global response of the structure is quantified, two environments can be more reliably and accurately compared. This chapter compares the response of an aerospace component in a service environment to the response of the same component in a laboratory test environment. The comparison first compares the mode shapes between the two environments. Once it is determined that the same mode shapes are present in both configurations, the modal accelerations are compared in order to determine the similarity of the global response of the component.

Photovoltaic modules undergoing laboratory hail tests were observed using high speed video to analyze the key characteristics of impact-induced glass fracture, including crack onset time, initiation location relative to the impact site, and propagation trends. Fifteen commercially representative glass-glass thin-film modules were recorded at 300,000 frames per second during hail impacts which happened to cause glass fracture. Images were processed to identify the time between impact and first plausible glass crack appearance (average 126 μs, standard deviation 59 μs) along with the time to a confirmed crack (average 158 μs, standard deviation 77μs), during the ice ball impacts which had a median kinetic energy of 47 J delivered by 55 mm diameter balls. Limiting factors for identifying glass crack timings were ice ball fragmentation obscuring the impact site and indistinct initial crack appearance, which were inherent to the images and not improved with processing. Computational simulations corresponding to each impact event showed that glass stresses were still localized to the impact site during times with definitively identifiable fracture, and even impacts which did not induce failure created local stress magnitudes exceeding stress levels associated with static glass fracture. These observations confirm that impact-induced glass failure is a time-and rate-dependent phenomena. Results from this study provide baseline metrics for developing a glass fracture criterion to predict module damage during hail impact events, which in turn allows for analysis of design features that may affect damage susceptibility.

Plenoptic background-oriented schlieren is a diagnostictechnique that enables the measure-ment of three-dimensional refractive gradients by a combination of background-oriented schlieren and a plenoptic light field camera. This plenoptic camera is a modification of a traditional camera via the insertion of an array of microlenses between the imaging lens and digital sensor. This allows the collection of both spatial and angular information on the incoming light rays and therefore provides three-dimensional information about the imaged scene. Background-oriented schlieren requires a relatively simple experimental configurationincludingonlyacameraviewing a patterned background through the density field of interest. By using a plenoptic camera to capture background-oriented schlieren images the optical distortion created by density gradients in three dimensions can be measured. This chapter is intended to review critical developments in plenoptic background-oriented schlieren imaging and provide an outlook for future applications of this measurement technique.

We present a novel design of a III-V-on-silicon heterogeneously integrated tunable ring laser, achieving >80 nanometers of tuning bandwidth, the widest conceived using only two rings, fostering many applications such as spectroscopy and beam steering.

Spatial navigation involves the formation of coherent representations of a map-like space, while simultaneously tracking current location in a primarily unsupervised manner. Despite a plethora of neurophysiological experiments revealing spatially-tuned neurons across the mammalian neocortex and subcortical structures, it remains unclear how such representations are acquired in the absence of explicit allocentric targets. Drawing upon the concept of predictive learning, we utilize a biologically plausible learning rule which utilizes sensory-driven observations with internally-driven expectations and learns through a contrastive manner to better predict sensory information. The local and online nature of this approach is ideal for deployment to neuromorphic hardware for edge-applications. We implement this learning rule in a network with the feedforward and feedback pathways known to be necessary for spatial navigation. After training, we find that the receptive fields of the modeled units resemble experimental findings, with allocentric and egocentric representations in the expected order along processing streams. These findings illustrate how a local and self-supervised learning method for predicting sensory information can extract latent structure from the environment.

I experimentally investigated and modeled a proposed frequency-domain method for detecting and tracking cislunar spacecraft and near-earth asteroids using heliostat fields at night. Unlike imaging, which detects spacecraft and asteroids by their streak in star-fixed long-exposure photographs, the proposed detection method oscillates the orientation of heliostats concentrating light from the stellar field and measures the light's photocurrent power spectrum at sub-milliHertz resolution. If heliostat oscillation repetitively traces out a closed loop fixed to the stars, light from spacecraft or asteroids moving along that loop produce photocurrent at a frequency slightly shifted from starlight. The frequency shift is proportional to the spacecraft or asteroid's apparent angular rate relative to sidereal. Relative phase corresponds to relative angular position, enabling tracking. Since heliostats are inexpensive compared to an astronomical observatory and otherwise unused at night, the proposed method may cost-effectively augment observatory systems such as NASA's Asteroid Terrestrial-impact Last Alert System (ATLAS).

The Explosive Destruction System (EDS) V31 containment vessel was procured by the US Army Recovered Chemical Materiel Directorate (RCMD) as a third-generation system used to destroy chemical munitions. It is the fifth individual EDS vessel to be fabricated under Code Case 2564 of the 2019 ASME Boiler and Pressure Vessel Code, which provides rules for the design of impulsively loaded vessels. The explosive rating for this vessel, based on the code case, is twenty-four (24) pounds TNT-equivalent for up to 1092 detonations. This report documents the results of explosive tests that were performed on the vessel at Sandia National Laboratories in Albuquerque, New Mexico to qualify the vessel for field operations use. There were three design basis configurations for qualification testing. Qualification test (1) consisted of a simulated M55 rocket motor and warhead assembly of 24lbs of Composition C-4 (30 lb TNT equivalent). This test was considered the maximum load case, based on modeling and simulation methods performed by Sandia prior to the vessel design phase. Qualification test (2) consisted of a regular, right circular cylinder, unitary charge, located central to the vessel interior of 19.2 lb of Composition C-4 (24 lb TNT equivalent). Qualification test (3) consisted of a 12-pack of regular, right circular cylinders, distributed evenly inside the vessel (totaling 19.2 lb of C-4, or 24 lb TNT equivalent). The ASME certification was based exclusively on the analytical predictions because the data required for certification cannot be obtained through testing. Strains through the thickness of the wall and on the inside surface of the cylinder were required and could only be obtained through analysis. Strain gages were placed on the outside of the vessel in three locations and the displacement of the door was measured using a Photonic Doppler Velocimetry (PDV) system. These measured values are compared to analytical predictions to help ensure the accuracy of the predicted strains and displacements throughout the rest of the model.

There has always been a desire to port high-fidelity reactive flow models from one code to another. For example, the AWE reactive burn model known as CREST has been or is being implemented in several of the U.S. Department of Energy hydrocodes. Those involved with reactive burn model implementation recognize the challenges immediately, e.g., Eulerian versus Lagrangian frameworks, the form of the equation of state, the closure relations, etc. In this work, we report the development of the CREST reactive burn model in CTH, a multidimensional, multi-material hydrocode developed by Sandia National Laboratories, following an earlier implementation shown at the last International Detonation Symposium. Results include code-to-code comparisons between CTH and the AWE hydrocode PERUSE, focusing on the simulated particle velocity histories during a shock-to-detonation transition, and corresponding to previous gas gun impact experiments as well as new model verification studies. Lessons learned are provided, including discussions of the numerical accuracy, in addition to the role of artificial viscosity and artificial viscous work. Finally, simulation results are shown to compare the Snowplough versus P-Alpha porosity model options.

Conceptual models of smectite hydration include planar (flat) clay layers that undergo stepwise expansion as successive monolayers of water molecules fill the interlayer regions. However, X-ray diffraction (XRD) studies indicate the presence of interstratified hydration states, suggesting non-uniform interlayer hydration in smectites. Additionally, recent theoretical studies have shown that clay layers can adopt bent configurations over nanometer-scale lateral dimensions with minimal effect on mechanical properties. Therefore, in this study we used molecular simulations to evaluate structural properties and water adsorption isotherms for montmorillonite models composed of bent clay layers in mixed hydration states. Results are compared with models consisting of planar clay layers with interstratified hydration states (e.g. 1W–2W). The small degree of bending in these models (up to 1.5 Å of vertical displacement over a 1.3 nm lateral dimension) had little or no effect on bond lengths and angle distributions within the clay layers. Except for models that included dry states, porosities and simulated water adsorption isotherms were nearly identical for bent or flat clay layers with the same averaged layer spacing. Similar agreement was seen with Na- and Ca-exchanged clays. While the small bent models did not retain their configurations during unconstrained molecular dynamics simulation with flexible clay layers, we show that bent structures are stable at much larger length scales by simulating a 41.6×7.1 nm2 system that included dehydrated and hydrated regions in the same interlayer.

This letter presents a one-step predictive run-to-run controller (R2R-MPC) for the automation of mechanical serial sectioning (MSS), a destructive material analysis process. To address the inherent uncertainty and disturbances in the MSS process, a robust closed-loop approach is presented. The robust R2R-MPC models the uncertainty of the MSS process using a linear differential inclusion. As an analytical model of the MSS process is unavailable, the differential inclusion is identified from historical data. The R2R-MPC is posed as an optimization problem that computes incremental changes to the control input which minimize the worst-case material removal errors. This optimization-based controller is combined with a run-to-run controller to provide integral action that rejects constant disturbances and tracks constant reference removal rates. To demonstrate the efficacy of our robust R2R-MPC, we present simulation results which compare the presented controller with a conventional non-robust R2R.

Neuromorphic computing uses brain-inspired concepts to accelerate and efficiently execute a wide range of applications, such as mimicking biological circuits, solving NPhard optimization problems and accelerating machine learning at the edge. In particular, neuromorphic architectures to efficiently execute Spiking Neural Networks (SNNs) have gained popularity. SNNs extend artificial neural networks (ANNs) by encoding information in time as either rates or delays between spiking events, shared between neurons via their weighted connections. SNN-based platforms are event-driven, resulting in naturally sparse, noise-tolerant and power-efficient computation. In this tutorial, we present the state-of-the-art in scalable digital and analog spiking neuromorphic system architectures, and discuss current research trends within the neuromorphic architecture field at the system level. We further introduce our SANA-FE tool for Simulation of Advanced Neuromorphic Architectures for Fast Exploration, which has been developed as part of a collaboration between the University of Texas at Austin and Sandia National Laboratories. SANA-FE allows for modeling and performance-power prediction of different spiking hardware architectures executing SNN applications to support rapid, early system-level design-space exploration, hardware-aware application development and system architecture co-design. The tutorial includes a hands-on component in which SANA-FE's capabilities are demonstrated and used to perform system design and application mapping case studies.