For turbulent reacting flow systems, identification of low-dimensional representations of the thermo-chemical state space is vitally important, primarily to significantly reduce the computational cost of device-scale simulations. Principal component analysis (PCA), and its variants, are a widely employed class of methods. Recently, an alternative technique that focuses on higher-order statistical interactions, co-kurtosis PCA (CoK-PCA), has been shown to effectively provide a low-dimensional representation by capturing the stiff chemical dynamics associated with spatiotemporally localized reaction zones. While its effectiveness has only been demonstrated based on a priori analyses with linear reconstruction, in this work, we employ nonlinear techniques to reconstruct the full thermo-chemical state and evaluate the efficacy of CoK-PCA compared to PCA. Specifically, we combine a CoK-PCA-/PCA-based dimensionality reduction (encoding) with an artificial neural network (ANN) based reconstruction (decoding) and examine, a priori, the reconstruction errors of the thermo-chemical state. In addition, we evaluate the errors in species production rates and heat release rates, which are nonlinear functions of the reconstructed state, as a measure of the overall accuracy of the dimensionality reduction technique. We employ four datasets to assess CoK-PCA/PCA coupled with ANN-based reconstruction: zero-dimensional (homogeneous) reactor for autoignition of an ethylene/air mixture that has conventional single-stage ignition kinetics, a dimethyl ether (DME)/air mixture which has two-stage (low and high temperature) ignition kinetics, a one-dimensional freely propagating premixed ethylene/air laminar flame, and a two-dimensional dataset representing turbulent autoignition of ethanol in a homogeneous charge compression ignition (HCCI) engine. Results from the analyses demonstrate the robustness of the CoK-PCA based low-dimensional manifold with ANN reconstruction in accurately capturing the data, specifically from the reaction zones.
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.
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.
In low inertia grids, significant frequency deviations can occur as a result of changes in power (load, generation, etc.), These deviations may activate various protection schemes designed to safeguard the system, potentially leading to blackouts. Therefore, assessing the frequency stability of the power system is crucial. The Frequency Security Index (FSI) serves as a metric for evaluating system stability. However, computing the FSI for a specific load change necessitates actual load changes on the system, which is often impractical. This paper introduces a method for calculating the FSI without requiring load changes for all values. A mathematical expression for the FSI is derived, which uses the values of microgrid parameters (such as inertia and damping constant) to compute the FSI for any load change. Subsequently, the parameters that most significantly affect the FSI are identified. Then, the paper introduces a Moving Horizon Estimation (MHE)-based parameter estimation approach, which leverages small perturbations from an energy storage system to estimate the most influential parameters for the FSI. The results show that the FSI calculation with the estimated parameters is more accurate (compared to COI averaged parameters), enabling a more effective state of health monitoring of the microgrid.
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.
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.
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.
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.
A mesoscale model for the shock initiation of pentaerythritol tetranitrate (PETN) films has been utilized to elucidate changes in initiation thresholds due to aging conditions and surface roughness, as has been observed from a series of high-throughput initiation (HTI) experiments. The HTI experiment has generated a wealth of thin-pulse, sub-millimeter shock initiation data for vapor deposited PETN films with thicknesses of 67-125 μm and varying accelerated aging conditions. This is because the HTI experiment provides access to growth-to-detonation information for explosives that exhibit a shock-to-detonation transition (SDT) with length and time scales that are too short to be resolved by conventional experiments. Mesoscale modeling results using experimentally characterized PETN microstructures are able to capture the general trend observed in experiments, in that increasing flyer impact velocity increases reactions until full detonation is reached. Moreover, the varying degrees of surface roughness that were considered were found to provide only minor variances in the peak particle velocity at the explosive output. The model did not predict a shift in the initiation threshold due to aged microstructures alone, indicating that additional mesoscale model improvements are necessary.
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.
Conceptual models of smectite hydration include planar (flat) clay layers that undergo stepwise expansion as successive monolayers of water molecules fill the interlayer regions. However, X-ray diffraction (XRD) studies indicate the presence of interstratified hydration states, suggesting non-uniform interlayer hydration in smectites. Additionally, recent theoretical studies have shown that clay layers can adopt bent configurations over nanometer-scale lateral dimensions with minimal effect on mechanical properties. Therefore, in this study we used molecular simulations to evaluate structural properties and water adsorption isotherms for montmorillonite models composed of bent clay layers in mixed hydration states. Results are compared with models consisting of planar clay layers with interstratified hydration states (e.g. 1W–2W). The small degree of bending in these models (up to 1.5 Å of vertical displacement over a 1.3 nm lateral dimension) had little or no effect on bond lengths and angle distributions within the clay layers. Except for models that included dry states, porosities and simulated water adsorption isotherms were nearly identical for bent or flat clay layers with the same averaged layer spacing. Similar agreement was seen with Na- and Ca-exchanged clays. While the small bent models did not retain their configurations during unconstrained molecular dynamics simulation with flexible clay layers, we show that bent structures are stable at much larger length scales by simulating a 41.6×7.1 nm2 system that included dehydrated and hydrated regions in the same interlayer.
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.
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.
We demonstrate for the first time waveguide integrated cascaded germanium photodetector arrays operated as photocells. We characterize several different array designs, and discuss their effects on voltage and photocurrent performance parameters.
This experimental work investigates the flow field associated with a tandem expansion-compression geometry in the Sandia Trisonic Wind Tunnel. PIV measurements of the boundary layer before the expansion characterized properties of the incoming boundary layer. Schlieren and oil-flow experiments were conducted at Mach 1.5 and 2 for a range of stagnation pressures. Further PIV experiments were conducted across the expansion and compression corners to observe the post-expansion changes in the boundary layer and its influence on the shock/boundary-layer interaction at the compression corner. Velocity profiles qualitatively matched well with RANS simulations and showed rapid growth of the boundary layer following the expansion, tapering off to a slower rate of growth with distance. Turbulence intensity, exemplified by the streamwise component of turbulent normal stress, diminished substantially after the expansion corner leading to the belief that relaminarization processes were occurring. Comparison with analysis from the literature suggests that the distortion of the boundary layer due to expansion leads to a separation length 30% larger than for an equilibrium turbulent boundary layer.
This work presents dose rate testing of an 8-bit microcontroller which executes code during radiation pulses. Multiple different failure modes are observed, which depend on the dose rate and integrated dose received by the device.
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.
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.
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.
This article aims at discovering the unknown variables in the system through data analysis. The main idea is to use the time of data collection as a surrogate variable and try to identify the unknown variables by modeling gradual and sudden changes in the data. We use Gaussian process modeling and a sparse representation of the sudden changes to efficiently estimate the large number of parameters in the proposed statistical model. The method is tested on a realistic dataset generated using a one-dimensional implementation of a Magnetized Liner Inertial Fusion (MagLIF) simulation model, and encouraging results are obtained.
Evaluation of the hydrogen storage capacity of porous rocks is crucial for underground hydrogen storage. Using 1H nuclear magnetic resonance (NMR) spectroscopy, we successfully characterized the hydrogen responses and identified storage mechanisms in Berea sandstone under varying water saturation. The results indicate that the injected hydrogen behaves as a free gas phase and is capable of occupying the empty pore volume regardless of the saturation state. No hysteresis was observed during injection and production cycles.
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.
Multi-axis testing has become a popular test method because it provides a more realistic simulation of a field environment when compared to traditional vibration testing. However, field data may not be available to derive the multi-axis environment. This means that methods are needed to generate “virtual field data” that can be used in place of measured field data. Transfer path analysis (TPA) has been suggested as a method to do this since it can be used to estimate the excitation forces on a legacy system and then apply these forces to a new system to generate virtual field data. This chapter will provide a review of using TPA methods to do this. It will include a brief background on TPA, discuss the benefits of using TPA to compute virtual field data, and delve into the areas for future work that could make TPA more useful in this application.
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.
Battery systems are typically equipped with state of charge (SoC) estimation algorithms. Sensor measurements used to estimate SoC are susceptible to false data injection attacks (FDIAs) that aim to disturb state estimation and, consequently, damage the system. In this paper, SoC estimation methods are re-purposed to detect FDIAs targeting the current and voltage sensors of a battery stack using a combination of an improved input noise aware unscented Kalman filter (INAUKF) and a cumulative sum detector. The root mean squared error of the states estimated by the INAUKF was at least 85% lower than the traditional unscented Kalman filter for all noise levels tested. The proposed method was able to detect FDIA in the current and voltage sensors of a series-connected battery stack in 99.55% of the simulations.
