Publications

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.

Anelastic strain recovery, the process of measuring the time dependent recovered strain after a core is cut at depth was utilized to make a measure of the in-situ properties stresses at depth at the FORGE (Frontier Observatory for Research in Geothermal Energy) site in Milford Utah. Core was collected from a region of well 16B at approximately 4860-4870 ft. Core was instrumented with strain gages within 10 hours of the core being cut. The relaxation of the cores was measured for approximately one month, and the results analyzed, which showed that the principal stresses were slightly off vertical, and magnitudes are close to equal.

A Marx generator module from the decommissioned RITS pulsed power machine from Sandia National Labs was modified to operate in an existing setup at Texas Tech University. This will ultimately be used as a testbed for laser triggered gas switching. The existing experimental setup at Texas Tech University consists of a large Marx tank, an oil-filled coaxial pulse forming line, an adjustable peaking gap, and load section along with various diagnostics. The setup was previously operated at a lower voltage than the new experiment, so electrostatic modeling was done to ensure viability and drive needed modifications. The oil tank will house the modified RITS Marx. This Marx contains half as many stages as the original RITS module and has an expected output of 1 MV. A trigger Marx generator consisting of 8 stages has been fabricated to trigger the RITS Marx. Charging and triggering of both Marx generators will be controlled through a fiber optic network. The output from the modified RITS Marx will be used to charge the oil-filled coaxial line acting as a low impedance pulse forming line (PFL). Once charged, the self-breaking peaking gap will close, allowing the compressed pulse to be released into the load section. For testing of the Marx module and PFL, a match 10 Ω water load was fabricated. The output pulsewidth is 55 nsec. Diagnostics include two capacitive voltage probes on either side of the peaking gap, a quarter-turn Rogowski coil for load current measurement, and a Pearson coil for calibrations purposes.

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.

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.

Deep neural networks (DNNs) achieve state-of-the-art performance in video anomaly detection. However, the usage of DNNs is limited in practice due to their computational overhead, generally requiring significant resources and specialized hardware. Further, despite recent progress, current evaluation criteria of video anomaly detection algorithms are flawed, preventing meaningful comparisons among algorithms. In response to these challenges, we propose (1) a compression-based technique referred to as Spatio-Temporal N-Gram Prediction by Partial Matching (STNG PPM) and (2) simple modifications to current evaluation criteria for improved interpretation and broader applicability across algorithms. STNG PMM does not require specialized hardware, has few parameters to tune, and is competitive with DNNs on multiple benchmark data sets in video anomaly detection.

In this work, we evaluate the usefulness of nonsmooth basis functions for representing the periodic response of a nonlinear system subject to contact/impact behavior. As with sine and cosine basis functions for classical Fourier series, which have C∞ smoothness, nonsmooth counterparts with C0 smoothness are defined to develop a nonsmooth functional representation of the solution. Some properties of these basis functions are outlined, such as periodicity, derivatives, and orthogonality, which are useful for functional series applied via the Galerkin method. Least-squares fits of the classical Fourier series and nonsmooth basis functions are presented and compared using goodness-of-fit metrics for time histories from vibro-impact systems with varying contact stiffnesses. This formulation has the potential to significantly reduce the computational cost of harmonic balance solvers for nonsmooth dynamical systems. Rather than requiring many harmonics to capture a system response using classical, smooth Fourier terms, the frequency domain discretization could be captured by a combination of a finite Fourier series supplemented with nonsmooth basis functions to improve convergence of the solution for contact-impact problems.

Gallium nitride (GaN)-based nanoscale vacuum electron devices, which offer advantages of both traditional vacuum tube operation and modern solid-state technology, are attractive for radiation-hard applications due to the inherent radiation hardness of vacuum electron devices and the high radiation tolerance of GaN. Here, we investigate the radiation hardness of top-down fabricated n-GaN nanoscale vacuum electron diodes (NVEDs) irradiated with 2.5-MeV protons (p) at various doses. We observe a slight decrease in forward current and a slight increase in reverse leakage current as a function of cumulative protons fluence due to a dopant compensation effect. The NVEDs overall show excellent radiation hardness with no major change in electrical characteristics up to a cumulative fluence of 5E14 p/cm2, which is significantly higher than the existing state-of-the-art radiation-hardened devices to our knowledge. The results show promise for a new class of GaN-based nanoscale vacuum electron devices for use in harsh radiation environments and space applications.

In this paper, we develop a nested chi-squared likelihood ratio test for selecting among shrinkage-regularized covariance estimators for background modeling in hyperspectral imagery. Critical to many target and anomaly detection algorithms is the modeling and estimation of the underlying background signal present in the data. This is especially important in hyperspectral imagery, wherein the signals of interest often represent only a small fraction of the observed variance, for example when targets of interest are subpixel. This background is often modeled by a local or global multivariate Gaussian distribution, which necessitates estimating a covariance matrix. Maximum likelihood estimation of this matrix often overfits the available data, particularly in high dimensional settings such as hyperspectral imagery, yielding subpar detection results. Instead, shrinkage estimators are often used to regularize the estimate. Shrinkage estimators linearly combine the overfit covariance with an underfit shrinkage target, thereby producing a well-fit estimator. These estimators introduce a shrinkage parameter, which controls the relative weighting between the covariance and shrinkage target. There have been many proposed methods for setting this parameter, but comparing these methods and shrinkage values is often performed with a cross-validation procedure, which can be computationally expensive and highly sample inefficient. Drawing from Bayesian regression methods, we compute the degrees of freedom of a covariance estimate using eigenvalue thresholding and employ a nested chi-squared likelihood ratio test for comparing estimators. This likelihood ratio test requires no cross-validation procedure and enables direct comparison of different shrinkage estimates, which is computationally efficient.

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.

We demonstrate evanescently coupled waveguide integrated silicon photonic avalanche photodiodes designed for single photon detection for quantum applications. Simulation, high responsivity, and record low dark currents for evanescently coupled devices are presented.

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.

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.

We investigate the kinetics and report the time-resolved concentrations of key chemical species in the oxidation of tetrahydrofuran (THF) at 7500 torr and 450-675 K. Experiments are carried out using high-pressure multiplexed photoionization mass spectrometry (MPIMS) combined with tunable vacuum ultraviolet radiation from the Berkely Lab Advanced Light Source. Intermediates and products are quantified using reference photoionization (PI) cross sections, when available, and constrained by a global carbon balance tracking approach at all experimental temperatures simultaneously for the species without reference cross sections. From carbon balancing, we determine time-resolved concentrations for the ROO˙ and ˙OOQOOH radical intermediates, butanedial, and the combined concentration of ketohydroperoxide (KHP) and unsaturated hydroperoxide (UHP) products stemming from the ˙QOOH + O2 reaction. Furthermore, we quantify a product that we tentatively assign as fumaraldehyde, which arises from UHP decomposition via H2O or ˙OH + H loss. The experimentally derived species concentrations are compared with model predictions using the most recent literature THF oxidation mechanism of Fenard et al., (Combust. Flame, 2018, 191, 252-269). Our results indicate that the literature mechanism significantly overestimates THF consumption and the UHP + KHP concentration at our conditions. The model predictions are sensitive to the rate coefficient for the ROO˙ isomerization to ˙QOOH, which is the gateway for radical chain propagating and branching pathways. Comparisons with our recent results for cyclopentane (Demireva et al., Combust. Flame, 2023, 257, 112506) provide insights into the effect of the ether group on reactivity and highlight the need to determine accurate rate coefficients of ROO˙ isomerization and subsequent reactions.

The MACCS code was created by Sandia National Laboratories for the U.S. Nuclear Regulatory Commission and has been used for emergency planning, level 3 probabilistic risk assessments, consequence analyses and other scientific and regulatory research for over half a century. Specializing in modeling the transport of nuclear material into the environment, MACCS accounts for atmospheric transport and dispersion, wet and dry deposition, probabilistic treatment of meteorology, exposure pathways, varying protective actions for the emergency, intermediate and long-term phases, dosimetry, health effects (including but not limited to population dose, acute radiation injury and increased cancer risk), and economic impacts. Routine updates and recent enhancements to the MACCS code, such as the inclusion of a higher fidelity atmospheric transport and dispersion model, the addition of a new economic impact model, and the application of nearfield modeling, have continuously increased the codes capabilities in consequence analysis. Additionally, investigations of MACCS capabilities for advanced reactor applications have shown that MACCS can provide realistic and informative risk assessments for the new generation of reactor designs. Even so, areas of improvement as well as gaps have been identified that if resolved can increase the usefulness of MACCS in any application regarding a release of nuclear material into the environment.

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.

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.

The siting of nuclear waste is a process that requires consideration of concerns of the public. This report demonstrates the significant potential for natural language processing techniques to gain insights into public narratives around “nuclear waste.” Specifically, the report highlights that the general discourse regarding “nuclear waste” within the news media has fluctuated in prevalence compared to “nuclear” topics broadly over recent years, with commonly mentioned entities reflecting a limited variety of geographies and stakeholders. General sentiments within the “nuclear waste” articles appear to use neutral language, suggesting that a scientific or “facts-only” framing of “waste”-related issues dominates coverage; however, the exact nuances should be further evaluated. The implications of a number of these insights about how nuclear waste is framed in traditional media (e.g., regarding emerging technologies, historical events, and specific organizations) are discussed. This report lays the groundwork for larger, more systematic research using, for example, transformer-based techniques and covariance analysis to better understand relationships among “nuclear waste” and other nuclear topics, sentiments of specific entities, and patterns across space and time (including in a particular region). By identifying priorities and knowledge needs, these data-driven methods can complement and inform engagement strategies that promote dialogue and mutual learning regarding nuclear waste.

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.

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.

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.

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.

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.

Rotational testbeds are ubiquitous in the test and evaluation of inertial sensors and systems. However, the use of rotational testbeds is typically restricted to static states employed over long integration windows to allow for data aggregation. These methods ignore the transitions between states, and data aggregation masks potentially useful signals. In this paper, we discuss the development of modular equations for the description of the inertial inputs to a test sensor using any rotational testbed. Implementing our equations in software, specific force and angular rates can be computed from idealized table motion or measured encoder data. Results are presented using simulated data and measured data. The measured data was acquired using a three axis rate table and a MEMS IMU. The experimental results validate our model equations and demonstrate the benefits of modeling sensor inputs at the sensor rates to compensate for testbed errors.

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.

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.

This paper develops a new domain-decomposition method for solving the KKT system with heat-equation constraints.

The development of multi-axis force sensing ca-pabilities in elastomeric materials has enabled new types of human motion measurement with many potential applications. In this work, we present a new soft insole that enables mobile measurement of ground reaction forces (GRFs) outside of a lab-oratory setting. This insole is based on hybrid shear and normal force detecting (SAND) tactile elements (taxels) consisting of optical sensors optimized for shear sensing and piezoresistive pressure sensors dedicated to normal force measurement. We develop polynomial regression and deep neural network (DNN) GRF prediction models and compare their performance to ground-truth force plate data during two walking experiments. Utilizing a 4-layer DNN, we demonstrate accurate prediction of the anterior-posterior (AP), medial-lateral (ML) and vertical components of the GRF with normalized mean absolute errors (NMAE) of <5.1 %, 4.1 %, and 4.5%, respectively. We also demonstrate the durability of the hybrid SAND insole construction through more than 20,000 cycles of use.

We demonstrate high-efficiency emission at wavelengths longer than 540 nm from InGaN quantum wells regrown on periodic arrays of GaN nanostructures and explore their incorporation into nanophotonic resonators for semiconductor laser development.

We use complete polarization tomography of photon pairs generated in semiconductor metasurfaces via spontaneous parametric down-conversion to show how bound states in the continuum resonances affect the polarization state of the emitted photons.

We present a materials study of AlGaInP grown on GaAs leveraging deep-level optical spectroscopy and time resolved photoluminescence. Our materials may serve as the basis for wide-bandgap analogs of silicon photomultipliers optimized for short wavelength sensing.

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.

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.

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.

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.

Sandia National Laboratories (SNL) has completed a comparative evaluation of three design assessment approaches for a 2-liter (2L) capacity containment vessel (CV) of a novel plutonium air transport (PAT) package designed to survive the hypothetical accident condition (HAC) test sequence defined in Title 10 of the United States (US) Code of Federal Regulations (CFR) Part 71.74(a), which includes a 129 meter per second (m/s) impact of the package into an essentially unyielding target. CVs for hazardous materials transportation packages certified in the US are typically designed per the requirements defined in the American Society of Mechanical Engineers (ASME) Boiler and Pressure Vessel Code (B&PVC) Section III Division 3 Subsection WB “Class TC Transportation Containments.” For accident conditions, the level D service limits and analysis approaches specified in paragraph WB-3224 are applicable. Data derived from finite element analyses of the 129 m/s impact of the 2L-PAT package were utilized to assess the adequacy of the CV design. Three different CV assessment approaches were investigated and compared, one based on stress intensity limits defined in subparagraph WB-3224.2 for plastic analyses (the stress-based approach), a second based on strain limits defined in subparagraph WB-3224.3, subarticle WB-3700, and Section III Nonmandatory Appendix FF for the alternate strain-based acceptance criteria approach (the strain-based approach), and a third based on failure strain limits derived from a ductile fracture model with dependencies on the stress and strain state of the material, and their histories (the Xue-Wierzbicki (X-W) failure-integral-based approach). This paper gives a brief overview of the 2L-PAT package design, describes the finite element model used to determine stresses and strains in the CV generated by the 129 m/s impact HAC, summarizes the three assessment approaches investigated, discusses the analyses that were performed and the results of those analyses, and provides a comparison between the outcomes of the three assessment approaches.

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.

Characterizing shielding effectiveness (SE) of enclosures is important in aerospace, military, and consumer applications. Direct SE measurement of an enclosure or chassis may be considered an exact characterization, but there are several sources of possible variability in such measurements, e.g., mechanical tolerances, the absence of components during test that exist in a final assembly, movement of components and cables, and perturbations due to probes and associated cabling. In [1] , internal stirrers were investigated as a way to sample the variation of SE of small enclosures when populated with random metallic objects. Here, we explore this idea as a way to quantify variability and sensitivity of an SE measurement, not only indicating the uncertainty of the SE measurement, but also delineating frequency ranges where either deterministic or statistical simulations should be applied.

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.

The diesel-piloted dual-fuel compression ignition combustion strategy is well-suited to accelerate the decarbonization of transportation by adopting hydrogen as a renewable energy carrier into the existing internal combustion engine with minimal engine modifications. Despite the simplicity of engine modification, many questions remain unanswered regarding the optimal pilot injection strategy for reliable ignition with minimum pilot fuel consumption. The present study uses a single-cylinder heavy-duty optical engine to explore the phenomenology and underlying mechanisms governing the pilot fuel ignition and the subsequent combustion of a premixed hydrogen-air charge. The engine is operated in a dual-fuel mode with hydrogen premixed into the engine intake charge with a direct pilot injection of n-heptane as a diesel pilot fuel surrogate. Optical diagnostics used to visualize in-cylinder combustion phenomena include high-speed IR imaging of the pilot fuel spray evolution as well as high-speed HCHO* and OH* chemiluminescence as indicators of low-temperature and high-temperature heat release, respectively. Three pilot injection strategies are compared to explore the effects of pilot fuel mass, injection pressure, and injection duration on the probability and repeatability of successful ignition. The thermodynamic and imaging data analysis supported by zero-dimensional chemical kinetics simulations revealed a complex interplay between the physical and chemical processes governing the pilot fuel ignition process in a hydrogen containing charge. Hydrogen strongly inhibits the ignition of pilot fuel mixtures and therefore requires longer injection duration to create zones with sufficiently high pilot fuel concentration for successful ignition. Results show that ignition typically tends to rely on stochastic pockets with high pilot fuel concentration, which results in poor repeatability of combustion and frequent misfiring. This work has improved the understanding on how the unique chemical properties of hydrogen pose a challenge for maximization of hydrogen's energy share in hydrogen dual-fuel engines and highlights a potential mitigation pathway.

Hail poses a significant threat to photovoltaic (PV) systems due to the potential for both cell and glass cracking. This work experimentally investigates hail-related failures in Glass/Backsheet and Glass/Glass PV modules with varying ice ball diameters and velocities. Post-impact Electroluminescence (EL) imaging revealed the damage extent and location, while high-speed Digital Image Correlation (DIC) measured the out-of-plane module displacements. The findings indicate that impacts of 20 J or less result in negligible damage to the modules tested. The thinner glass in Glass/Glass modules cracked at lower impact energies (-25 J) than Glass/Backsheet modules (-40 J). Furthermore, both module types showed cell and glass cracking at lower energies when impacted at the module's edges compared to central impacts. At the time of presentation, we will use DIC to determine if out-of-plane displacements are responsible for the impact location discrepancy and provide more insights into the mechanical response of hail impacted modules. This study provides essential insights into the correlation between impact energy, impact location, displacements, and resulting damage. The findings may inform critical decisions regarding module type, site selection, and module design to contribute to more reliable PV systems.

In this work, the frequency response of a simplified shaft-bearing assembly is studied using numerical continuation. Roller-bearing clearances give rise to contact behavior in the system, and past research has focused on the nonlinear normal modes of the system and its response to shock-type loads. A harmonic balance method (HBM) solver is applied instead of a time integration solver, and numerical continuation is used to map out the system’s solution branches in response to a harmonic excitation. Stability analysis is used to understand the bifurcation behavior and possibly identify numerical or system-inherent anomalies seen in past research. Continuation is also performed with respect to the forcing magnitude, resulting in what are known as S-curves, in an effort to detect isolated solution branches in the system response.

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.

Abstract not provided.

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.

Interim dry storage of spent nuclear fuel involves storing the fuel in welded stainless-steel canisters. Under certain conditions, the canisters could be subjected to environments that may promote stress corrosion cracking leading to a risk of breach and release of aerosol-sized particulate from the interior of the canister to the external environment through the crack. Research is currently under way by several laboratories to better understand the formation and propagation of stress corrosion cracks, however little work has been done to quantitatively assess the potential aerosol release. The purpose of the present work is to introduce a reliable generic numerical model for prediction of aerosol transport, deposition, and plugging in leak paths similar to stress corrosion cracks, while accounting for potential plugging from particle deposition. The model is dynamic (changing leak path geometry due to plugging) and it relies on the numerical solution of the aerosol transport equation in one dimension using finite differences. The model’s capabilities were also incorporated into a Graphical User Interface (GUI) that was developed to enhance user accessibility. Model validation efforts presented in this paper compare the model’s predictions with recent experimental data from Sandia National Laboratories (SNL) and results available in literature. We expect this model to improve the accuracy of consequence assessments and reduce the uncertainty of radiological consequence estimations in the remote event of a through-wall breach in dry cask storage systems.

Underground caverns in salt formations are promising geologic features to store hydrogen (H2) because of salt's extremely low permeability and self-healing behavior.Successful salt-cavern H2 storage schemes must maximize the efficiency of cyclic injection-production while minimizing H2 loss through adjacent damaged salt.The salt cavern storage community, however, has not fully understood the geomechanical behaviors of salt rocks driven by quick operation cycles of H2 injection-production, which may significantly impact the cost-effective storage-recovery performance.Our field-scale generic model captures the impact of combined drag and back stressing on the salt creep behavior corresponding to cycles of compression and extension, which may lead to substantial loss of cavern volumes over time and diminish the cavern performance for H2 storage.Our preliminary findings address that it is essential to develop a new salt constitutive model based on geomechanical tests of site-specific salt rock to probe the cyclic behaviors of salt both beneath and above the dilatancy boundary, including reverse (inverse transient) creep, the Bauschinger effect and fatigue.

This paper provides a summary of planning work for experiments that will be necessary to address the long-term model validation needs required to meet offshore wind energy deployment goals. Conceptual experiments are identified and laid out in a validation hierarchy for both wind turbine and wind plant applications. Instrumentation needs that will be required for the offshore validation experiments to be impactful are then listed. The document concludes with a nominal vision for how these experiments can be accomplished.

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.