Publications

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.

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.

Contact mechanics, or the modeling of the impenetrability of solid objects, is fundamental to computational solid mechanics (CSM) applications yet is oftentimes the most challenging in terms of computational efficiency and performance. These challenges arise from the irregularity and highly dynamic nature of contact simulation, particularly with algorithms designed for distributed memory architectures. First among these challenges is the inherent load imbalance when distributing contact load across compute nodes. This imbalance is highly problem dependent, and relates to the surface area of contact manifolds and the volume around them, rather than the distribution of the mesh over compute nodes, meaning the application load can vary drastically over different phases. The dynamic nature of contact problems motivates the use of distributed asynchronous many-tasking (AMT) frameworks to efficiently handle irregular workloads. In this paper, we present our work on distBVH, a distributed contact solution using the DARMA/vt library for asynchronous tasking that is also capable of running on-node Kokkos-based kernels. We explore how distBVH addresses the various challenges of CSM contact problems. We evaluate the use of many of DARMA/vt’s dynamic load balancers and demonstrate how our load balancing approach can provide significant performance improvements on various computational solid mechanics benchmarks. Additionally, we show how our approach can take advantage of DARMA/vt for tasking and efficient on-node kernels using Kokkos to scale over hundreds of processing elements.

Closed-loop geothermal systems (CLGSs) rely on circulation of a heat transfer fluid in a closed-loop design without penetrating the reservoir to extract subsurface heat and bring it to the surface. We developed and applied numerical models to study u-shaped and coaxial CLGSs in hot-dry-rock over a more comprehensive parameter space than has been studied before, including water and supercritical CO2 (sCO2) as working fluids. An economic analysis of each realization was performed to evaluate the levelized cost of heat (LCOH) for direct heating application and levelized cost of electricity (LCOE) for electrical power generation. The results of the parameter study, composed of 2.5 million simulations, combined with a plant and economic model comprise the backbone of a publicly accessible web application that can be used to query, analyze, and plot outlet states, thermal and mechanical power output, and LCOH/LCOE, thereby facilitating feasibility studies led by potential developers, geothermal scientists, or the general public (https://gdr.openei.org/submissions/1473). Our results indicate competitive LCOH can be achieved; however, competitive LCOE cannot be achieved without significant reductions in drilling costs. We also present a site-based case study for multi-lateral systems and discuss how our comprehensive single-lateral analyses can be applied to approximate multi-lateral CLGSs. Looking beyond hot-dry-rock, we detail CLGS studies in permeable wet rock, albeit for a more limited parameter space, indicating that reservoir permeability of greater than 250 mD is necessary to significantly improve CLGS power production, and that reservoir temperatures greater than 200 °C, achieved by going to greater depths (∼3–4 km), may significantly enhance power production.

Network interface controllers (NICs) with general-purpose compute capabilities ('SmartNICs') present an opportunity for reducing host application overheads by offloading non-critical tasks to the NIC. In addition to moving computation, offloading requires that associated data is also transferred to the NIC. To meet this need, we introduce a high-performance, general-purpose data movement service that facilitates the of-floading of tasks to SmartNICs: The SmartNIC Data Movement Service (SDMS). SDMS provides near-line-rate transfer band-widths between the host and NIC. Moreover, SDMS's In-transit Data Placement (IDP) feature can reduce (or even eliminate) the cost of serializing data on the NIC by performing the necessary data formatting during the transfer. To illustrate these capabilities, we provide an in-depth case study using SDMS to offload data management operations related to Apache Arrow, a popular data format standard. For single-column tables, SDMS can achieve more than 87% of baseline throughput for data buffers that are 128 KiB or larger (and more than 95% of baseline throughput for buffers that are 1 MiB or larger) while also nearly eliminating the host and SmartNIC overhead associated with Arrow operations.

Accuracy-optimized convolutional neural networks (CNNs) have emerged as highly effective models at predicting neural responses in brain areas along the primate ventral stream, but it is largely unknown whether they effectively model neurons in the complementary primate dorsal stream. We explored how well CNNs model the optic flow tuning properties of neurons in dorsal area MSTd and we compared our results with the Non-Negative Matrix Factorization (NNMF) model, which successfully models many tuning properties of MSTd neurons. To better understand the role of computational properties in the NNMF model that give rise to optic flow tuning that resembles that of MSTd neurons, we created additional CNN model variants that implement key NNMF constraints – non-negative weights and sparse coding of optic flow. While the CNNs and NNMF models both accurately estimate the observer's self-motion from purely translational or rotational optic flow, NNMF and the CNNs with nonnegative weights yield substantially less accurate estimates than the other CNNs when tested on more complex optic flow that combines observer translation and rotation. Despite its poor accuracy, NNMF gives rise to tuning properties that align more closely with those observed in primate MSTd than any of the accuracy-optimized CNNs. This work offers a step toward a deeper understanding of the computational properties and constraints that describe the optic flow tuning of primate area MSTd.

We consider the problem of decentralized control of reactive power provided by distributed energy resources for voltage support in the distribution grid. We assume that the reactance matrix of the grid is unknown and potentially time-varying. We present a decentralized adaptive controller in which the reactive power at each inverter is set using a potentially heterogeneous droop curve and analyze the stability and the steady-state error of the resulting system. The effectiveness of the controller is validated in simulations using a modified version of the IEEE 13-bus and a 8500-node test system.

The radar generalized image quality equation (RGIQE) is a metric used to measure both monostatic and bistatic synthetic aperture radar (BSAR) image quality, it is a function of signal-to-noise ratio (SNR) and 2-D bandwidth. The 2-D bandwidth is equal to the area of the transfer function's (TF) passband region. With the exception of side-looking monostatic geometries, almost all monostatic and bistatic geometries have skewed passband shapes when waveform frequency parameters remain unchanged from pulse to pulse. Most synthetic aperture radar (SAR) applications require a rectangular-shaped passband region, this is achieved by inscribing a rectangular region within the skewed intrinsic passband region. Increasing skewness results in less inscription area reducing 2-D bandwidth, image SNR, and thus RGIQE capacity. In this article, a waveform with frequency agility is used to rectify the skewness that degrades RGIQE capacity. By changing the waveform's center frequency and instantaneous bandwidth from pulse to pulse in a particular manner, the intrinsic passband region can be de-skewed. The de-skewed shape maximizes the inscription area thus maximizing 2-D bandwidth, image SNR, and RGIQE capacity. Three examples are given in this article, one monostatic geometry, and two bistatic geometries. RGIQE capacity is increased by 52.02%, 44.42%, and 79.09% for the three examples.

We demonstrate an InAs-based terahertz (THz) metasurface emitter that can generate and focus THz pulses using a binary-phase Fresnel zone plate concept. The metalens emitter successfully generates a focused THz beam without additional THz optics.

Multifidelity emulators have found wide-ranging applications in both forward and inverse problems within the computational sciences. Thanks to recent advancements in neural architectures, they provide significant flexibility for integrating information from multiple models, all while retaining substantial efficiency advantages over single-fidelity methods. In this context, existing neural multifidelity emulators operate by separately resolving the linear and nonlinear correlation between equally parameterized high-and low-fidelity approximants. However, many complex models ensembles in science and engineering applications only exhibit a limited degree of linear correlation between models. In such a case, the effectiveness of these approaches is impeded, i.e., larger datasets are needed to obtain satisfactory predictions. In this work, we present a general strategy that seeks to maximize the linear correlation between two models through input encoding. We showcase the effectiveness of our approach through six numerical test problems, and we show the ability of the proposed multifidelity emulator to accurately recover the high-fidelity model response under an increasing number of quasi-random samples. In our experiments, we show that input encoding produces in many cases emulators with significantly simpler nonlinear correlations. Finally, we demonstrate how the input encoding can be leveraged to facilitate the fusion of information between low-and high-fidelity models with dissimilar parametrization, i.e., situations in which the number of inputs is different between low-and high-fidelity models.

Accurate understanding of the behavior of commercial-off-the-shelf electrical devices is important in many applications. This paper discusses methods for the principled statistical analysis of electrical device data. We present several recent successful efforts and describe two current areas of research that we anticipate will produce widely applicable methods. Because much electrical device data is naturally treated as functional, and because such data introduces some complications in analysis, we focus on methods for functional data analysis.

National Security Presidential Memorandum-20 defines three tier levels for launch approval of space nuclear systems. The two main factors determining the tier level are the total quantity and type of radioactive sources and the probability of any member of the public receiving doses above certain thresholds. The total quantity of radioactive sources is compared with International Atomic Energy Agency transportation regulations. The dose probability is determined by the product of three terms: 1) the probability of a launch accident occurring; 2) the probability of a release of radioactive material given an accident; and 3) the probability of exceeding the dose threshold to any member of the public given a release. This paper provides a methodology for evaluating these values and applies this methodology to an example mission as a demonstration. For the example mission, a preliminary tier determination of Tier III was concluded.

Over the past few years, advancements in closed-loop geothermal systems (CLGS), also called advanced geothermal systems (AGS), have sparked a renewed interest in these types of designs. CLGS have certain advantages over traditional and enhanced geothermal systems (EGS), including not requiring in-situ reservoir permeability, conservation of the circulating fluid, and allowing for different fluids, including working fluids directly driving a turbine at the surface. CLGS may be attractive in environments where water resources are limited, rock contaminants must be avoided, and stimulation treatments are not available (e.g., due to regulatory or technical reasons). Despite these advantages, CLGS have some challenges, including limited surface area for heat transfer and requiring long wellbores and laterals to obtain multi-MW output in conduction-only reservoirs. CLGS have been investigated in conduction-only systems. In this paper, we explore the impact of both forced and natural convection on the levels of heat extraction with a CLGS deployed in a hot wet rock reservoir. We bound potential benefits of convection by investigating liquid reservoirs over a range of natural and forced convective coefficients. Additionally, we investigate the effects of permeability, porosity, and geothermal temperature gradient in the reservoir on CLGS outputs. Reservoir simulations indicate that reservoir permeabilities of at least ~100 mD are required for natural convection to increase the heat output with respect to a conduction-only scenario. The impact increases with increasing reservoir temperature. When subject to a forced convection flow field, Darcy velocities of at least 10-7 m/s are required to obtain an increase in heat output.

Multiple scattering is a common phenomenon in acoustic media that arises from the interaction of the acoustic field with a network of scatterers. This mechanism is dominant in problems such as the design and simulation of acoustic metamaterial structures often used to achieve acoustic control for sound isolation, and remote sensing. In this study, we present a physics-informed neural network (PINN) capable of simulating the propagation of acoustic waves in an infinite domain in the presence of multiple rigid scatterers. This approach integrates a deep neural network architecture with the mathematical description of the physical problem in order to obtain predictions of the acoustic field that are consistent with both governing equations and boundary conditions. The predictions from the PINN are compared with those from a commercial finite element software model in order to assess the performance of the method.

Visualization of flow structures within post-detonation fireballs has been performed for benchmark validation of numerical simulations. Custom pressed PETN explosives with a 12-mm diameter hemispherical form factor were used to produce a spherically symmetric post-detonation flow with low soot yield. Hydroxyl-radical planar laser induce fluorescence (OH-PLIF) was employed to visualize the structure ranging from approximately 10μs to 35μs after shock breakout from the explosive pellet. Fireball simulations were performed using the HyBurn Computational Fluid Dynamics (CFD) package. Experimental OH-PLIF results were compared to synthetic OH-PLIF from post-processing of CFD simulations. From the comparison of experimental and synthetic OH-PLIF images, CFD is shown to replicate much of the flow structure observed in the experiments, revealing potential differences in turbulent length scales and OH kinetics. Results provide significant advancement in experimental resolution of these harsh turbulent combustion environments and validate physical models thereof.

A single Synthetic Aperture Radar (SAR) image is a 2-Dimensional projection of a 3-Dimensional scene, with very limited ability to estimate surface topography. However, with multiple SAR images collected from suitably different geometries, they may be compared with multilateration calculations to estimate characteristics of the missing dimension. The ability to employ effective multilateration algorithms is highly dependent on the geometry of the data collections, and can be cast as a least-squares exercise. A measure of Dilution of Precision (DOP) can be used to compare the relative merits of various collection geometries.

The impact of more extreme climate conditions under global warming on soil organic carbon (SOC) dynamics remains unquantified. Here we estimate the response of SOC to climate extreme shifts under 1.5 °C warming by combining a space-for-time substitution approach and global SOC measurements (0–30 cm soil). Most extremes (22 out of 33 assessed extreme types) exacerbate SOC loss under warming globally, but their effects vary among ecosystems. Only decreasing duration of cold spells exerts consistent positive effects, and increasing extreme wet days exerts negative effects in all ecosystems. Temperate grasslands and croplands negatively respond to most extremes, while positive responses are dominant in temperate and boreal forests and deserts. In tundra, 21 extremes show neutral effects, but 11 extremes show negative effects with stronger magnitude than in other ecosystems. Our results reveal distinct, biome-specific effects of climate extremes on SOC dynamics, promoting more reliable SOC projection under climate change.

Traditional electronics assemblies are typically packaged using physically or chemically blown potted foams to reduce the effects of shock and vibration. These potting materials have several drawbacks including manufacturing reliability, lack of internal preload control, and poor serviceability. A modular foam encapsulation approach combined with additively manufactured (AM) silicone lattice compression structures can address these issues for packaged electronics. These preloaded silicone lattice structures, known as foam replacement structures (FRSs), are an integral part of the encapsulation approach and must be properly characterized to model the assembly stresses and dynamics. In this study, dynamic test data is used to validate finite element models of an electronics assembly with modular encapsulation and a direct ink write (DIW) AM silicone FRS. A variety of DIW compression architectures are characterized, and their nominal stress-strain behavior is represented with hyperfoam constitutive model parameterizations. Modeling is conducted with Sierra finite element software, specifically with a handoff from assembly preloading and uniaxial compression in Sierra/Solid Mechanics to linear modal and vibration analysis in Sierra/Structural Dynamics. This work demonstrates the application of this advanced modeling workflow, and results show good agreement with test data for both static and dynamic quantities of interest, including preload, modal, and vibration response.

This paper presents a low-power staggered tunning LNA with a variable threshold limiter designed for impulse-radio ultra-wideband (IR-UWB). This amplifier has a high gain of 36.5 dB with a 3dB S21 frequency of 6.8 - 9.4 GHz while consuming 4.28 mW from a 0.8 V supply. The input return loss is better than 8 dB across the band and 10 dB from 6.9 GHz. The output return loss is better than 18 dB across the operating bandwidth. At 9 GHz, the minimum noise figure measured is 4.85 dB. The OP1dB compression point is measured at 7.8 GHz as -5.5 dBm. We also present an adjustable threshold limiter, providing additional protection over a typical diode limiter circuit. This amplifier is fabricated in a 45nm PD-SOI process. To the author's knowledge, this LNA demonstrates the highest linear figure of merit (FoM) in C and X band work.

Ionic liquid (IL) pretreatment methods show incredible promise for the efficient conversion of lignocellulosic feedstocks to fuels and chemicals. Given their low vapor pressures, distillation-based methods of extracting ionic liquids out of biomass post-pretreatment have historically been ignored in favor of alternative methods. We demonstrate a process to distill four acetate-based ionic liquids ([EthA][OAc], [PropA][OAc], [MAEthA][OAc], and [DMAEthA][OAc]) at low pressure and high purity that overcome some disadvantages of “water washing” and “one pot” recovery methods. Out of four tested ILs, ethanolamine acetate ([EthA][OAc]) is shown to have the most agreeable conversion metrics for commercial bioconversion processes achieving 73.6 % and 51.4 % of theoretical glucose and xylose yields respectively and >85 % recovery rates. Our process metrics are factored into a techno-economic analysis where [EthA][OAc] distillation is compared to other recovery methods as well as ethanolamine pretreatment at both milliliter and liter scales. Although our TEA shows [EthA][OAc] distillation underperforming against other processes, we show a step-by-step avenue to reduce sugar production cost below the wholesale dextrose price at scale.

Point-source transient events (PSTEs) - optical events that are both extremely fast and extremely small - pose several challenges to an imaging system. Due to their speed, accurately characterizing such events often requires detectors with very high frame rates. Due to their size, accurately detecting such events requires maintaining coverage over an extended field-of-view, often through the use of imaging focal plane arrays (FPA) with a global shutter readout. Traditional imaging systems that meet these requirements are costly in terms of price, size, weight, power consumption, and data bandwidth, and there is a need for cheaper solutions with adequate temporal and spatial coverage. To address these issues, we develop a novel compressed sensing algorithm adapted to the rolling shutter readout of an imaging system. This approach enables reconstruction of a PSTE signature at the sampling rate of the rolling shutter, offering a 1-2 order of magnitude temporal speedup and a proportional reduction in data bandwidth. We present empirical results demonstrating accurate recovery of PSTEs using measurements that are spatially undersampled by a factor of 25, and our simulations show that, relative to other compressed sensing algorithms, our algorithm is both faster and yields higher quality reconstructions. We also present theoretical results characterizing our algorithm and corroborating simulations. The potential impact of our work includes the development of much faster, cheaper sensor solutions for PSTE detection and characterization.

For an Energy System to be truly equitable, it should provide affordable and reliable energy services to disadvantaged and underserved populations. Disadvantaged communities often face a combination of economic, social, health, and environmental burdens and may be geographically isolated (e.g., rural communities), which systematically limits their opportunity to fully participate in aspects of economic, social, and civic life.

Explosives exposed to conditions above the Chapman-Jouget (CJ) state exhibit an overdriven response that is transient. Reactive flow models are often fit to the CJ conditions, and they transition to detonation based on inputs lower than or near CJ, but these models may also be used to predict explosive behavior in the overdriven regime. One scenario that can create a strongly overdriven state is a Mach stem shock interaction. These interactions can drive an already detonating or transitioning explosive to an overdriven state, and they can also cause detonation at the interaction location where the separate shocks may be insufficient to detonate the material. In this study, the reactive flow model XHVRB utilizing a Mie-Grüneisen equation of state (EOS) for the unreacted explosive, and a Sesame table for the reacted products, will be used to examine Mach stem interactions from multi-point detonation schemes in CTH. The effect of the overdriven response driven by PETN-based explosive pellets will be tracked to determine the transient detonation behavior, and the predicted states from the burn model will be compared to previously published data.

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.

Fault location, isolation, and service restoration of a self-healing, self-Assembling microgrid operating off-grid from distributed inverter-based resources (IBRs) can be a unique challenge because of the fault current limitations and uncertainties regarding which sources are operational at any given time. The situation can become even more challenging if data sharing between the various microgrid controllers, relays, and sources is not available. This paper presents an innovative robust partitioning approach, which is used as part of a larger self-Assembling microgrid concept utilizing local measurements only. This robust partitioning approach splits a microgrid into sub-microgrids to isolate the fault to just one of the sub-microgrids, allowing the others to continue normal operation. A case study is implemented in the IEEE 123-bus distribution test system in Simulink to show the effectiveness of this approach. The results indicate that including the robust partitions leads to less loss of load and shorter overall restoration times.

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.

Determining the thermal response of energetic materials at high densities can be difficult when pressure dependent reactions occur within the interior of the material. At high temperatures, reactive components such as hexahydro-l,3,5-tri-nitro-l,3,5-triazine (RDX), ammonium perchlorate (AP), and hydroxyl-terminated polybutadiene (HTPB) decompose and interact. The decomposition products accumulate near defects where internal pressure ultimately causes mechanical damage with closed pores transitioning into open pores. Gases are no longer confined locally; instead, they freely migrate between open pores and ultimately escape into the surrounding headspace or vent. Recently we have developed a universal cookoff model (UCM) coupled to a micromechanics pressurization (MMP) model to address pressure-dependent reactions that occur within the interior of explosives. Parameters for the UCM/MMP model are presented for an explosive and two propellants that contain similar portions of both aluminum (Al) and a binder. The explosive contains RDX and the propellants contain AP with no RDX. One of the propellants contains small amounts of curing catalysts and a burn modifier whereas the other propellant does not. We found that the cookoff behavior of the two propellants behave similarly leading and conclude that small amounts of catalysts or burn modifiers do not influence cookoff behavior appreciably. Kinetic parameters for the UCM/MMP models were obtained from the Sandia Instrumented Thermal Ignition (SITI) experiment. Validation is done with data from other laboratories.

In general, multiple-input/multiple-output (MIMO) vibration testing utilizes a response-controlled test methodology where specifications are in the form of response quantities at various locations distributed on the device under test (DUT). There are some advantages to this approach, namely that DUT response could be measured in some field environment and directly used as MIMO specifications for subsequent MIMO vibration tests on similar DUTs. However, in some cases it may be advantageous to control the MIMO vibration test at the inputs rather than the responses. One such case is free-flight environments, where the DUT is unconstrained, and all loads come from aerodynamic pressures. In this case, the force-controlled test method is much more robust to system changes such as unit-to-unit variability as compared to a response-controlled test method. This could make force-controlled MIMO test specifications more generalizable and easier to derive. This is exactly akin to transfer path analysis, where pseudo-forces are applicable in special circumstances. This paper will explore the force-controlled test concept and demonstrate it with a numerical example, comparing performance under various conditions vs. the traditional response-controlled test method.

Resonant plate shock testing techniques have been used for mechanical shock testing at Sandia for several decades. A mechanical shock qualification test is often done by performing three separate uniaxial tests on a resonant plate to simulate one shock event. Multi-axis mechanical shock activities, in which shock specifications are simultaneously met in different directions during a single shock test event performed in the lab, are not always repeatable and greatly depend on the fixture used during testing. This chapter provides insights into various designs of a concept fixture that includes both resonant plate and angle bracket used for multi-axis shock testing from a modeling and simulation point of view based on the results of finite element modal analysis. Initial model validation and testing performed show substantial excitation of the system under test as the fundamental modes drive the response in all three directions. The response also shows that higher order modes are influencing the system, the axial and transverse response are highly coupled, and tunability is difficult to achieve. By varying the material properties, changing thicknesses, adding masses, and moving the location of the fixture on the resonant plate, the response can be changed significantly. The goal of this work is to identify the parameters that have the greatest influence on the response of the system when using the angle bracket fixture for a mechanical shock test for the intent of tunability of the system.

Blind geothermal systems are believed to be common in the Basin and Range province and represent an underutilized source of renewable green energy. Their discovery has historically been by chance but more methodological strategies for exploration of these resources are being developed. One characteristic of blind systems is that they are often overlain by near-surface zones of low-resistivity caused by alteration of the overlying sediments to swelling clays. These zones can be imaged by resistivity-based geophysical techniques to facilitate their discovery and characterization. Here we present a side-by-side comparison of resistivity models produced from helicopter transient electromagnetic (HTEM) and ground-based broadband magnetotelluric (MT) surveys over a previously discovered blind geothermal system with measured shallow temperatures of ~100°C in East Hawthorne, NV. The HTEM and MT data were collected as part of the BRIDGE project, an initiative for improving methodologies for discovering blind geothermal systems. HTEM data were collected and modelled along profiles, and the results suggest the method can resolve the resistivity structure 300 - 500 m deep. A 61-station MT survey was collected on an irregular grid with ~800 m station spacing and modelled in 3D on a rotated mesh aligned with HTEM flight directions. Resistivity models are compared with results from potential fields datasets, shallow temperature surveys, and available temperature gradient data in the area of interest. We find that the superior resolution of the HTEM can reveal near-surface details often missed by MT. However, MT is sensitive to several km deep, can resolve 3D structures, and is thus better suited for single-prospect characterization. We conclude that HTEM is a more practical subregional prospecting tool than is MT, because it is highly scalable and can rapidly discover shallow zones of low resistivity that may indicate the presence of a blind geothermal system. Other factors such as land access and ground disturbance considerations may also be decisive in choosing the best method for a particular prospect. Resistivity methods in general cannot fully characterize the structural setting of a geothermal system, and so we used potential fields and other datasets to guide the creation of a diagrammatic structural model at East Hawthorne.

The ability to accurately predict the structure and dynamics of pool fires using computational simulations is of great interest in a wide variety of applications, including accidental and wildland fires. However, the presence of physical processes spanning a broad range of spatial and temporal scales poses a significant challenge for simulations of such fires, particularly at conditions near the transition between laminar and turbulent flow. In this study, we examine the transition to turbulence in methane pool fires using high-resolution simulations with multi-step finite rate chemistry, where adaptive mesh refinement (AMR) is used to directly resolve small-scale flow phenomena. We perform three simulations of methane pool fires, each with increasing diameter, corresponding to increasing inlet Reynolds and Richardson numbers. As the diameter increases, the flow transitions from organized vortex roll-up via the puffing instability to much more chaotic mixing associated with finger formation along the shear layer and core collapse near the inlet. These effects combine to create additional mixing close to the inlet, thereby enhancing fuel consumption and causing more rapid acceleration of the fluid above the pool. We also make comparisons between the transition to turbulence and core collapse in the present pool fires and in inert helium plumes, which are often used as surrogates for the study of buoyant reacting flows.

A new Adaptive Mesh Refinement (AMR) keyword was added to the CTH¹ hydrocode developed at Sandia National Laboratories (SNL). The new indicator keyword, "ratec*ycle", allows the user to specify the minimum number of computational cycles before an AMR block is allowed to be un-refined. This option is designed to allow the analyst to control how quickly a block is un-refined to avoid introducing anomalous waves in their solution due to information propagating across mesh resolution changes. For example, in reactive flow simulations it is often desirable to accurately capture the expansion region behind the reaction front. The effect of this new option was examined using the XHVRB^{2, 3} model for XTX8003 to model the propagation of the detonation wave in explosives in small channels, and also for a simpler explosive model driving a steel case. The effect on computational cost as a function of this new option was also examined.

The primary goal of any laboratory test is to expose the unit-under-test to conservative realistic representations of a field environment. Satisfying this objective is not always straightforward due to laboratory equipment constraints. For vibration and shock tests performed on shakers over-testing and unrealistic failures can result because the control is a base acceleration and mechanical shakers have nearly infinite impedance. Force limiting and response limiting are relatively standard practices to reduce over-test risks in random-vibration testing. Shaker controller software generally has response limiting as a built-in capability and it is done without much user intervention since vibration control is a closed loop process. Limiting in shaker shocks is done for the same reasons, but because the duration of a shock is only a few milliseconds, limiting is a pre-planned user in the loop process. Shaker shock response limiting has been used for at least 30 years at Sandia National Laboratories, but it seems to be little known or used in industry. This objective of this paper is to re-introduce response limiting for shaker shocks to the aerospace community. The process is demonstrated on the BARBECUE testbed.

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.

Quantum sensing utilizes the inherent sensitivity of a quantum system to external stimuli. Our goal is to leverage this sensitivity to develop a quantum sensor designed for the detection of ionizing radiation. Here we report on the design, fabrication, and measurement of a new quantum device for hard x-ray and gamma-ray detection. Our quantum device is based on a superconducting quantum bit (qubit) with superconducting tunnel junctions as the core device elements. We describe our experimental investigation directed toward the detection metrics of energy resolution, dynamic range, and active area. In contrast to existing superconducting detectors, the active area per qubit may be much larger than the physical area of the tunnel junctions or the physical area of the qubit device, due to the sensitivity of quantum coherence to ionizing radiation deposition within a radius on the millimeter or centimeter scale. Our experimental design enables an ionizing radiation source at room temperature to be detected by our quantum sensor at low temperature.

Operation and control of a galvanically isolated three-phase AC-AC converter for solid state transformer applications is described. The converter regulates bidirectional power transfer by phase shifting voltages applied on either side of a high-frequency transformer. The circuit structure and control system are symmetrical around the transformer. Each side operates independently, enabling conversion between AC systems with differing voltage magnitude, phase angle, and frequency. This is achieved in a single conversion stage with low component count and high efficiency. The modulation strategy is discussed in detail and expressions describing the relationship between phase shift and power transfer are presented. Converter operation is demonstrated in a 3 kW hardware prototype.

Additive manufacturing has ushered in a new paradigm of bottom-up materials-by-design of spatially non-uniform materials. Functionally graded materials have locally tailored compositions to provide optimized global properties and performance. In this letter, we propose an opportunity for the application of graded magnetic materials as lens elements for charged particle optics. A Hiperco50/Hymu80 (FeCo-2 V/Fe-80Ni-5Mo) graded magnetic alloy was successfully additively manufactured via Laser Directed Energy Deposition with spatially varying magnetic properties. The compositional gradient is then applied using computational simulations to demonstrate how a tailored material can enhance the magnetic performance of a critical, image-forming component of a transmission electron microscope.

Disposal of commercial spent nuclear fuel in a geologic repository is studied. In situ heater experiments in underground research laboratories provide a realistic representation of subsurface behavior under disposal conditions. This study describes process model development and modeling analysis for a full-scale heater experiment in opalinus clay host rock. The results of thermal-hydrology simulation, solving coupled nonisothermal multiphase flow, and comparison with experimental data are presented. The modeling results closely match the experimental data.

The explosive BTF (benzotrifuroxan) is an interesting molecule for sub-millimeter studies of initiation and detonation. It has no hydrogen, thus no water in the detonation products and a subsequently high temperature in the reaction zone. The material has impact sensitivity that is comparable or less than that of PETN (pentaerythritol tetranitrate) and slightly greater than RDX, HMX, and CL-20. Physical vapor deposition (PVD) can be used to grow high-density films of pure explosives with precise control over geometry, and we apply this technique to BTF to study detonation and initiation behavior as a function of sample thickness. The geometrical effects on detonation and corner turning behavior are studied with the critical detonation thickness experiment and the micromushroom test, respectively. Initiation behavior is studied with the high-throughput initiation experiment. Vapor-deposited films of BTF show detonation failure, corner turning, and initiation consistent with a heterogeneous explosive. Scaling of failure thickness to failure diameter shows that BTF has a very small failure diameter.

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.

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.

Deep neural networks for automatic target recognition (ATR) have been shown to be highly successful for a large variety of Synthetic Aperture Radar (SAR) benchmark datasets. However, the black box nature of neural network approaches raises concerns about how models come to their decisions, especially when in high-stake scenarios. Accordingly, a variety of techniques are being pursued seeking to offer understanding of machine learning algorithms. In this paper, we first provide an overview of explainability and interpretability techniques introducing their concepts and the insights they produce. Next we summarize several methods for computing specific approaches to explainability and interpretability as well as analyzing their outputs. Finally, we demonstrate the application of several attribution map methods and apply both attribution analysis metrics as well as localization interpretability analysis to six neural network models trained on the Synthetic and Measured Paired Labeled Experiment (SAMPLE) dataset to illustrate the insights these methods offer for analyzing SAR ATR performance.

Visualization of flow structures within post-detonation fireballs has been performed for benchmark validation of numerical simulations. Custom pressed PETN explosives with a 12-mm diameter hemispherical form factor were used to produce a spherically symmetric post-detonation flow with low soot yield. Hydroxyl-radical planar laser induce fluorescence (OH-PLIF) was employed to visualize the structure ranging from approximately 10μs to 35μs after shock breakout from the explosive pellet. Fireball simulations were performed using the HyBurn Computational Fluid Dynamics (CFD) package. Experimental OH-PLIF results were compared to synthetic OH-PLIF from post-processing of CFD simulations. From the comparison of experimental and synthetic OH-PLIF images, CFD is shown to replicate much of the flow structure observed in the experiments, revealing potential differences in turbulent length scales and OH kinetics. Results provide significant advancement in experimental resolution of these harsh turbulent combustion environments and validate physical models thereof.

A new Adaptive Mesh Refinement (AMR) keyword was added to the CTH¹ hydrocode developed at Sandia National Laboratories (SNL). The new indicator keyword, "ratec*ycle", allows the user to specify the minimum number of computational cycles before an AMR block is allowed to be un-refined. This option is designed to allow the analyst to control how quickly a block is un-refined to avoid introducing anomalous waves in their solution due to information propagating across mesh resolution changes. For example, in reactive flow simulations it is often desirable to accurately capture the expansion region behind the reaction front. The effect of this new option was examined using the XHVRB^{2, 3} model for XTX8003 to model the propagation of the detonation wave in explosives in small channels, and also for a simpler explosive model driving a steel case. The effect on computational cost as a function of this new option was also examined.

Our present electric power grid maximizes spinning inertia of fossil fuel generators (inherent energy storage) to meet stability and performance requirements. Our goal is to begin to investigate the replacement of the large spinning inertia of fossil fuel generators with energy storage systems (ESS) including information flow as a necessary part of the renewable energy sources (RES) and subject to certain criteria. General criteria metrics include: energy storage, information flow, estimation, communication links, central versus decentralized, etc. Our focus is on evaluating the Fisher Information Equivalency (FIE) metric as a multi-criteria trade-off cost function for the minimization of ESS options and information flow. This paper begins with a formal conceptual definition of an infinite bus. Then a simple example of a One Machine Infinite Bus (OMIB) system with a Unified Power Flow Controller (UPFC) to demonstrate the FIE-based approach to minimize the ESS. A second more detailed example of several spinning machines are included with representative power electronic and ESS for RES that are attached to the electric power grid. A simple trade-study begins to highlight requirements to support large penetration of RES. Keep in mind for a large scale high penetration of RES will require large investments in ESS which we want to minimize.

Determining the thermal response of energetic materials at high densities can be difficult when pressure dependent reactions occur within the interior of the material. At high temperatures, reactive components such as hexahydro-l,3,5-tri-nitro-l,3,5-triazine (RDX), ammonium perchlorate (AP), and hydroxyl-terminated polybutadiene (HTPB) decompose and interact. The decomposition products accumulate near defects where internal pressure ultimately causes mechanical damage with closed pores transitioning into open pores. Gases are no longer confined locally; instead, they freely migrate between open pores and ultimately escape into the surrounding headspace or vent. Recently we have developed a universal cookoff model (UCM) coupled to a micromechanics pressurization (MMP) model to address pressure-dependent reactions that occur within the interior of explosives. Parameters for the UCM/MMP model are presented for an explosive and two propellants that contain similar portions of both aluminum (Al) and a binder. The explosive contains RDX and the propellants contain AP with no RDX. One of the propellants contains small amounts of curing catalysts and a burn modifier whereas the other propellant does not. We found that the cookoff behavior of the two propellants behave similarly leading and conclude that small amounts of catalysts or burn modifiers do not influence cookoff behavior appreciably. Kinetic parameters for the UCM/MMP models were obtained from the Sandia Instrumented Thermal Ignition (SITI) experiment. Validation is done with data from other laboratories.

Dendrites enable neurons to perform nonlinear operations. Existing silicon dendrite circuits sufficiently model passive and active characteristics, but do not exploit shunting inhibition as an active mechanism. We present a dendrite circuit implemented on a reconfigurable analog platform that uses active inhibitory conductance signals to modulate the circuit's membrane potential. We explore the potential use of this circuit for direction selectivity by emulating recent observations demonstrating a role for shunting inhibition in a directionally-selective Drosophila (Fruit Fly) neuron.

Before residential photovoltaic (PV) systems are interconnected with the grid, various planning and impact studies are conducted on detailed models of the system to ensure safety and reliability are maintained. However, these model-based analyses can be time-consuming and error-prone, representing a potential bottleneck as the pace of PV installations accelerates. Data-driven tools and analyses provide an alternate pathway to supplement or replace their model-based counterparts. In this article, a data-driven algorithm is presented for assessing the thermal limitations of PV interconnections. Using input data from residential smart meters, and without any grid models or topology information, the algorithm can determine the nameplate capacity of the service transformer supplying those customers. The algorithm was tested on multiple datasets and predicted service transformer capacity with >98% accuracy, regardless of existing PV installations. This algorithm has various applications from model-free thermal impact analysis for hosting capacity studies to error detection and calibration of existing grid models.

Frequency stability issues are common in low-inertia microgrids due to the dominance of power electronics converter-based resources and comparatively low levels of inertia of synchronous generation. While many methods have been proposed for frequency support in these systems, it is challenging to ensure both stability and an adequate level of support without directly modeling the system. Perturbation-based extremum-seeking control (PESC) is a model-free adaptive control strategy that optimally sets the system's performance measure without requiring a mathematical system model. However, PESC offers poor transient performance due to design limitations, such as an averaging theory. Some modifications have been made to address these limitations; nevertheless, the modified design of PESC operates more as a model-based control. The Quasi-Newton method is a popular class of numerical optimizers attributed with a design procedure similar to model-free control that uses the gradient and an approximate inverse Hessian of the performance measure to run an optimization loop. This paper presents the design and compares the performance of the Quasi-Newton method and a model-based PESC for frequency support of microgrids. The simulation results illustrate the comparable performance of both control schemes and show the model-free control capability of the numerical optimization method for a class of nonlinear dynamic systems.

Tabulated chemistry models are widely used to simulate large-scale turbulent fires in applications including energy generation and fire safety. Tabulation via piecewise Cartesian interpolation suffers from the curse-of-dimensionality, leading to a prohibitive exponential growth in parameters and memory usage as more dimensions are considered. Artificial neural networks (ANNs) have attracted attention for constructing surrogates for chemistry models due to their ability to perform high-dimensional approximation. However, due to well-known pathologies regarding the realization of suboptimal local minima during training, in practice they do not converge and provide unreliable accuracy. Partition of unity networks (POUnets) are a recently introduced family of ANNs which preserve notions of convergence while performing high-dimensional approximation, discovering a mesh-free partition of space which may be used to perform optimal polynomial approximation. We assess their performance with respect to accuracy and model complexity in reconstructing unstructured flamelet data representative of nonadiabatic pool fire models. Our results show that POUnets can provide the desirable accuracy of classical spline-based interpolants with the low memory footprint of traditional ANNs while converging faster to significantly lower errors than ANNs. For example, we observe POUnets obtaining target accuracies in two dimensions with 40 to 50 times less memory and roughly double the compression in three dimensions. We also address the practical matter of efficiently training accurate POUnets by studying convergence over key hyperparameters, the impact of partition/basis formulation, and the sensitivity to initialization.