The Cramér-Rao Lower Bound (CRLB) is used as a classical benchmark to assess estimators. Online algorithms for estimating modal properties from ambient data, i.e., mode meters, can benefit from accurate estimates of forced oscillations. The CRLB provides insight into how well forced oscillation parameters, e.g., frequency and amplitude, can be estimated. Previous works have solved the lower bound under a single-channel PMU measurement; thus, this paper extends works further to study CRLB under two-channel PMU measurements. The goal is to study how correlated/uncorrelated noise can affect estimation accuracy. Interestingly, these studies shows that correlated noise can decrease the CRLB in some cases. This paper derives the CRLB for the two-channel case and discusses factors that affect the bound.
In accident scenarios involving release of tritium during handling and storage, the level of risk to human health is dominated by the extent to which radioactive tritium is oxidized to the water form (T2O or THO). At some facilities, tritium inventories consist of very small quantities stored at sub-atmospheric pressure, which means that tritium release accident scenarios will likely produce concentrations in air that are well below the lower flammability limit. It is known that isotope effects on reaction rates should result in slower oxidation rates for heavier isotopes of hydrogen, but this effect has not previously been quantified for oxidation at concentrations well below the lower flammability limit for hydrogen. This work describes hydrogen isotope oxidation measurements in an atmospheric tube furnace reactor. These measurements consist of five concentration levels between 0.01% and 1% protium or deuterium and two residence times. Oxidation is observed to occur between about 550°C and 800°C, with higher levels of conversion achieved at lower temperatures for protium with respect to deuterium at the same volumetric inlet concentration and residence time. Computational fluid dynamics simulations of the experiments were used to customize reaction orders and Arrhenius parameters in a 1-step oxidation mechanism. The trends in the rates for protium and deuterium are extrapolated based on guidance from literature to produce kinetic rate parameters appropriate for tritium oxidation at low concentrations.
Integrating recent advancements in resilient algorithms and techniques into existing codes is a singular challenge in fault tolerance - in part due to the underlying complexity of implementing resilience in the first place, but also due to the difficulty introduced when integrating the functionality of a standalone new strategy with the preexisting resilience layers of an application. We propose that the answer is not to build integrated solutions for users, but runtimes designed to integrate into a larger comprehensive resilience system and thereby enable the necessary jump to multi-layered recovery. Our work designs, implements, and verifies one such comprehensive system of runtimes. Utilizing Fenix, a process resilience tool with integration into preexisting resilience systems as a design priority, we update Kokkos Resilience and the use pattern of VeloC to support application-level integration of resilience runtimes. Our work shows that designing integrable systems rather than integrated systems allows for user-designed optimization and upgrading of resilience techniques while maintaining the simplicity and performance of all-in-one resilience solutions. More application-specific choice in resilience strategies allows for better long-term flexibility, performance, and - importantly - simplicity.
This work explores deriving transmissibility functions for a missile from a measured location at the base of the fairing to a desired location within the payload. A pressure on the outside of the fairing and the rocket motor’s excitation creates an acceleration at a measured location and a desired location. Typically, the desired location is not measured. In fact, it is typical that the payload may change, but measured acceleration at the base of the fairing is generally similar to previous test flights. Given this knowledge, it is desired to use a finite-element model to create a transmissibility function which relates acceleration from the previous test flight’s measured location at the base of the fairing to acceleration at a location in the new payload. Four methods are explored for deriving this transmissibility, with the goal of finding an appropriate transmissibility when both the pressure and rocket motor excitation are equally present. These methods are assessed using transient results from a simple example problem, and it is found that one of the methods gives good agreement with the transient results for the full range of loads considered.
Adaptive protection is defined as a real-time system that can modify the protective actions according to the changes in the system condition. An adaptive protection system (APS) is conventionally coordinated through a central management system located at the distribution system substation. An APS depends significantly on the communication infrastructure to monitor the latest status of the electric power grid and send appropriate settings to all of the protection relays existing in the grid. This makes an APS highly vulnerable to communication system failures (e.g., broken communication links due to natural disasters as well as wide-range cyber-attacks). To this end, this paper presents the addition of local adaptive modular protection (LAMP) units to the protection system to guarantee its reliable operation under extreme events when the operation of the APS is compromised. LAMP units operate in parallel with the conventional APS. As a backup, if APS fails to operate because of an issue in the communication system, LAMP units can accommodate a reliable fault detection and location on behalf of the protection relay. The performance of the proposed APS is verified using IEEE 123 node test system.
A high-speed, two-color pyrometer was developed and employed to characterize the temperature of the ejecta from pyrotechnic igniters. The pyrometer used a single objective lens, beamsplitter, and two high-speed cameras to maximize the spatial and temporal resolutions. The pyrometer used the integrated intensity of under-resolved particles to maintain a large region of interest to capture more particles. The spectral response of the pyrometer was determined based on the response of each optical component and the total system was calibrated using a black body source to ensure accurate intensity ratios over the range of interest.
A method for measuring the velocity of a passing object by recording the voltage induced by a magnet fixed to that object is explored in this article. Both theory and applications are discussed along with the strengths and limitations of the technique.
This paper presents a literature review on current practices and trends on cyberphysical security of grid-connected battery energy storage systems (BESSs). Energy storage is critical to the operation of Smart Grids powered by intermittent renewable energy resources. To achieve this goal, utility-scale and consumer-scale BESS will have to be fully integrated into power systems operations, providing ancillary services and performing functions to improve grid reliability, balance power and demand, among others. This vision of the future power grid will only become a reality if BESS are able to operate in a coordinated way with other grid entities, thus requiring significant communication capabilities. The pervasive networking infrastructure necessary to fully leverage the potential of storage increases the attack surface for cyberthreats, and the unique characteristics of battery systems pose challenges for cyberphysical security. This paper discusses a number of such threats, their associated attack vectors, detection methods, protective measures, research gaps in the literature and future research trends.
This paper studies the differences in a synthetic inertia controller of using two different feedback measurements: (i) an estimate of the rate of change of frequency from local voltage measurements, and (ii) a remote machine acceleration from a generator nearby to the actuator. The device that provides the synthetic inertia action is a converter interfaced generator (CIG). The paper carries out analysis in the frequency domain, using Bode plots, to show that synthetic inertia control using frequency estimates is more prone to instabilities than for the case where a machine speed is used. The paper then proposes a controller (or a filter) to mitigate these effects. In addition, the paper shows the effects that a delay of the machine speed signal of the nearby generator has on the synthetic inertia control of the system and how a controller is also needed in this case. Finally, the paper shows the difference in performance of a synthetic inertia controller when using these different measurement signals with simulations in time domain a electromagnetic transient program platform.
OpenMP 5.0 added support for reductions over explicit tasks. This expands the previous reduction support that was limited primarily to worksharing and parallel constructs. While the scope of a reduction operation in a worksharing construct is the scope of the construct itself, the scope of a task reduction can vary. This difference requires syntactical means to define the scope of reductions, e.g., the task_reduction clause, and to associate participating tasks, e.g., the in_reduction clause. Furthermore, the disassociation of the number of threads and the number of tasks creates space for different implementations in the OpenMP runtime. In this work, we provide insights into the behavior and performance of task reduction implementations in GCC/g++ and LLVM/Clang. Our results indicate that task reductions are well supported by both compilers, but their performance differs in some cases and is often determined by the efficiency of the underlying task management.
Incorrect modeling of control characteristics for inverter-based resources (IBRs) can affect the accuracy of electric power system studies. In many distribution system contexts, the control settings for behind-the-meter (BTM) IBRs are unknown. This paper presents an efficient method for selecting a small number of time series samples from net load meter data that can be used for reconstructing or classifying the control settings of BTM IBRs. Sparse approximation techniques are used to select the time series samples that cause the inversion of a matrix of candidate responses to be as well-conditioned as possible. We verify these methods on 451 actual advanced metering infrastructure (AMI) datasets from loads with BTM IBRs. Selecting 60 15-minute granularity time series samples, we recover BTM control characteristics with a mean error less than 0.2 kVAR.
Limited access to transmission lines after a major contingency event can inhibit restoration efforts. After Hurricane Maria, for example, flooding and landslides damaged roads and thus limited travel. Transmission lines are also often situated far from maintained roadways, further limiting the ability to access and repair them. Therefore, this paper proposes a methodology for assessing Puerto Rico's infrastructure (i.e., roads and transmission lines) to identify potentially hard to reach areas due to natural risks or distance to roads. The approach uses geographic information system (GIS) data to define vulnerable areas, that may experience excessive restoration times. The methodology also uses graph theory analysis to find transmission lines with high centrality (or importance). Comparison of these important transmission lines with the vulnerability results found that many reside near roads that are at risk for landslides or floods.
In this paper, we address the problem of convergence of sequential variational inference filter (VIF) through the application of a robust variational objective and H∞-norm based correction for a linear Gaussian system. As the dimension of state or parameter space grows, performing the full Kalman update with the dense covariance matrix for a large-scale system requires increased storage and computational complexity, making it impractical. The VIF approach, based on mean-field Gaussian variational inference, reduces this burden through the variational approximation to the covariance usually in the form of a diagonal covariance approximation. The challenge is to retain convergence and correct for biases introduced by the sequential VIF steps. We desire a frame-work that improves feasibility while still maintaining reasonable proximity to the optimal Kalman filter as data is assimilated. To accomplish this goal, a H∞-norm based optimization perturbs the VIF covariance matrix to improve robustness. This yields a novel VIF-H∞ recursion that employs consecutive variational inference and H∞ based optimization steps. We explore the development of this method and investigate a numerical example to illustrate the effectiveness of the proposed filter.
Performance assessment is an important tool to estimate the long-term safety for a nuclear waste repository. Performance assessment simulations are subject to multiple kinds of uncertainty including stochastic uncertainty, state of knowledge uncertainty, and model uncertainty. Task F1 of the DECOVALEX project involves comparison of the models and methods used in post-closure performance assessment of deep geologic repositories in fractured crystalline rock, providing an opportunity to compare the effects of different sources of uncertainty. A generic reference case for a mined repository in fractured crystalline rock was put together by participating teams, where each team was responsible for determining how best to represent and implement the model. This work presents the preliminary crystalline reference case results for the Department of Energy (DOE) team.
The penetration of renewable energy resources (RER) and energy storage systems (ESS) into the power grid has been accelerated in recent times due to the aggressive emission and RER penetration targets. The Integrated resource planning (IRP) framework can help in ensuring long-term resource adequacy while satisfying RER integration and emission reduction targets in a cost-effective and reliable manner. In this paper, we present pIRP (probabilistic Integrated Resource Planning), an open-source Python-based software tool designed for optimal portfolio planning for an RER and ESS rich future grid and for addressing the capacity expansion problem. The tool, which is planned to be released publicly, with its ESS and RER modeling capabilities along with enhanced uncertainty handling make it one of the more advanced non-commercial IRP tools available currently. Additionally, the tool is equipped with an intuitive graphical user interface and expansive plotting capabilities. Impacts of uncertainties in the system are captured using Monte Carlo simulations and lets the users analyze hundreds of scenarios with detailed scenario reports. A linear programming based architecture is adopted which ensures sufficiently fast solution time while considering hundreds of scenarios and characterizing profile risks with varying levels of RER and ESS penetration levels. Results for a test case using data from parts of the Eastern Interconnection are provided in this paper to demonstrate the capabilities offered by the tool.
Many teams struggle to adapt and right-size software engineering best practices for quality assurance to fit their context. Introducing software quality is not usually framed in a way that motivates teams to take action, thus resulting in it becoming a "check the box for compliance"activity instead of a cultural practice that values software quality and the effort to achieve it. When and how can we provide effective incentives for software teams to adopt and integrate meaningful and enduring software quality practices? We explored this question through a persona-based ideation exercise at the 2021 Collegeville Workshop on Scientific Software in which we created three unique personas that represent different scientific software developer perspectives.
The goal of this paper is to present a set of measurements from a benchmark structure containing two bolted joints to support future efforts to predict the damping due to the joints and to model nonlinear coupling between the first two elastic modes. Bolted joints introduce nonlinearities in structures, typically causing a softening in the natural frequency and an increase in damping because of frictional slip between the contact interfaces within the joint. These nonlinearities pose significant challenges when characterizing the response of the structure under a large range of load amplitudes, especially when the modal responses become coupled, causing the effective damping and natural frequency to not only depend on the excitation amplitude of the targeted mode, but also the relative amplitudes of other modes. In this work, two nominally identical benchmark structures, known in some prior works as the S4 beam, are tested to characterize their nonlinear properties for the first two elastic modes. Detailed surface measurements are presented and validated through finite element analysis and reveal distinct contact interactions between the two sets of beams. The free-free test structures are excited with an impact hammer and the transient response is analyzed to extract the damping and frequency backbone curves. A range of impact amplitudes and drive points are used to isolate a single mode or to excite both modes simultaneously. Differences in the nonlinear response correlate with the relative strength of the modes that are excited, allowing one to characterize mode coupling. Each of the beams shows different nonlinear properties for each mode, which is attributed to the different contact pressure distributions between the parts, although the mode coupling relationship is found to be consistent between the two. The test data key finding are presented in this paper and the supporting data is available on a public repository for interested researchers.
Proceedings of the Nuclear Criticality Safety Division Topical Meeting, NCSD 2022 - Embedded with the 2022 ANS Annual Meeting
Salazar, Alex
The postclosure criticality safety assessment for the direct disposal of dual-purpose canisters (DPCs) in a geologic repository includes considerations of transient criticality phenomena. The power pulse from a hypothetical transient criticality event in an unsaturated alluvial repository is evaluated for a DPC containing 37 spent pressurized water reactor (PWR) assemblies. The scenario assumes that the conditions for baseline criticality are achieved through flooding with groundwater and progressive failure of neutron absorbing media. A preliminary series of steady-state criticality calculations is conducted to characterize reactivity feedback due to absorber degradation, Doppler broadening, and thermal expansion. These feedback coefficients are used in an analysis with a reactor kinetics code to characterize the transient pulse given a positive reactivity insertion for a given length of time. The time-integrated behavior of the pulse can be used to model effects on the DPC and surrounding barriers in future studies and determine if transient criticality effects are consequential.
A 0.2-2 GHz digitally programmable RF delay element based on a time-interleaved multi-stage switched-capacitor (TIMS-SC) approach is presented. The proposed approach enables hundreds of ns of broadband RF delay by employing sample time expansion in multiple stages of switched-capacitor storage elements. The delay element was implemented in a 45 nm SOI CMOS process and achieves a 2.55-448.6 ns programmable delay range with < 0.12% delay variation across 1.8 GHz of bandwidth at maximum delay, 2.42 ns programmable delay steps, and 330 ns/mm2 area efficiency. The device achieves 24 dB gain, 7.1 dB noise figure, and consumes 80 mW from a 1 V supply with an active area of 1.36 mm2.
This paper describes how the performance of motion primitive-based planning algorithms can be improved using reinforcement learning. Specifically, we describe and evaluate a framework that autonomously improves the performance of a primitive-based motion planner. The improvement process consists of three phases: exploration, extraction, and reward updates. This process can be iterated continuously to provide successive improvement. The exploration step generates new trajectories, and the extraction step identifies new primitives from these trajectories. These primitives are then used to update rewards for continued exploration. This framework required novel shaping rewards, development of a primitive extraction algorithm, and modification of the Hybrid A* algorithm. The framework is tested on a navigation task using a nonlinear F-16 model. The framework autonomously added 91 motion primitives to the primitive library and reduced average path cost by 21.6 s, or 35.75% of the original cost. The learned primitives are applied to an obstacle field navigation task, which was not used in training, and reduced path cost by 16.3 s, or 24.1%. Additionally, two heuristics for the modified Hybrid A* algorithm are designed to improve effective branching factor.
The proper coordination of power system protective devices is essential for maintaining grid safety and reliability but requires precise knowledge of fault current contributions from generators like solar photovoltaic (PV) systems. PV inverter fault response is known to change with atmospheric conditions, grid conditions, and inverter control settings, but this time-varying behavior may not be fully captured by conventional static fault studies that are used to evaluate protection constraints in PV hosting capacity analyses. To address this knowledge gap, hosting capacity protection constraints were evaluated on a simplified test circuit using both a time-series fault analysis and a conventional static fault study approach. A PV fault contribution model was developed and utilized in the test circuit after being validated by hardware experiments under various irradiances, fault voltages, and advanced inverter control settings. While the results were comparable for certain protection constraints, the time-series fault study identified additional impacts that would not have been captured with the conventional static approach. Overall, while conducting full time-series fault studies may become prohibitively burdensome, these findings indicate that existing fault study practices may be improved by including additional test scenarios to better capture the time-varying impacts of PV on hosting capacity protection constraints.
With the increase in penetration of inverter-based resources (IBRs) in the electrical power system, the ability of these devices to provide grid support to the system has become a necessity. With standards previously developed for the interconnection requirements of grid-following inverters (GFLI) (most commonly photovoltaic inverters), it has been well documented how these inverters 'should' respond to changes in voltage and frequency. However, with other IBRs such as grid-forming inverters (GFMIs) (used for energy storage systems, standalone systems, and as uninterruptable power supplies) these requirements are either: not yet documented, or require a more in deep analysis. With the increased interest in microgrids, GFMIs that can be paralleled onto a distribution system have become desired. With the proper control schemes, a GFMI can help maintain grid stability through fast response compared to rotating machines. This paper will present an experimental comparison of commercially available GFMIand GFLI ' responses to voltage and frequency deviation, as well as the GFMIoperating as a standalone system and subjected to various changes in loads.
A forward analytic model is required to rapidly simulate the neutron time-of-flight (nToF) signals that result from magnetized liner inertial fusion (MagLIF) experiments at Sandia’s Z Pulsed Power Facility. Various experimental parameters, such as the burn-weighted fuel-ion temperature and liner areal density, determine the shape of the nToF signal and are important for characterizing any given MagLIF experiment. Extracting these parameters from measured nToF signals requires an appropriate analytic model that includes the primary deuterium-deuterium neutron peak, once-scattered neutrons in the beryllium liner of the MagLIF target, and direct beamline attenuation. Mathematical expressions for this model were derived from the general-geometry time- and energy-dependent neutron transport equation with anisotropic scattering. Assumptions consistent with the time-of-flight technique were used to simplify this linear Boltzmann transport equation into a more tractable form. Models of the uncollided and once-collided neutron scalar fluxes were developed for one of the five nToF detector locations at the Z-Machine. Numerical results from these models were produced for a representative MagLIF problem and found to be in good agreement with similar neutron transport simulations. Twenty experimental MagLIF data sets were analyzed using the forward models, which were determined to only be significantly sensitive to the ion temperature. The results of this work were also found to agree with values obtained separately using a zero scatter analytic model and a high-fidelity Monte Carlo simulation. Inherent difficulties in this and similar techniques are identified, and a new approach forward is suggested.
We propose the use of balanced iterative reducing and clustering using hierarchies (BIRCH) combined with linear regression to predict the reduced Young's modulus and hardness of highly heterogeneous materials from a set of nanoindentation experiments. We first use BIRCH to cluster the dataset according to its mineral compositions, which are derived from the spectral matching of energy-dispersive spectroscopy data through the modular automated processing system (MAPS) platform. We observe that grouping our dataset into five clusters yields the best accuracy as well as a reasonable representation of mineralogy in each cluster. Subsequently, we test four types of regression models, namely linear regression, support vector regression, Gaussian process regression, and extreme gradient boosting regression. The linear regression and Gaussian process regression provide the most accurate prediction, and the proposed framework yields R2 = 0.93 for the test set. Although the study is needed more comprehensively, our results shows that machine learning methods such as linear regression or Gaussian process regression can be used to accurately estimate mechanical properties with a proper number of grouping based on compositional data.
Quantifying gas-surface interactions for hypersonic reentry applications remains a challenging and complex problem where credible models are needed to design and analyze thermal protection systems. A flexible sensitivity analysis approach is demonstrated to analyze finite-rate ablation models to identify reaction parameters and mechanisms of influence on predicted quantities of interest. Simulations of hypersonic flow over a sphere-cone are presented using parameterized Park, Zhluktov and Abe (ZA), and MURI finite-rate models that describe the oxidation and sublimation of carbon. The results presented in this study emphasize the importance of characterizing model inputs that are shown to have a high impact on predicted quantities and build evidence to assess credibility of these models.
Primary frequency control in power systems is being challenged by the large-scale integration of inverter-based resources (IBRs) because they do not typically respond to frequency fluctuations. This paper suggests introducing new reserve products into the electricity market that provide incentive for IBRs to contribute to primary frequency control in ways that take advantage of their fast-acting capabilities. In addition to a Primary Frequency Response (PFR) reserve product, which accommodates standard droop control, we suggest introducing a Fast Frequency Response (FFR) reserve product, a reserve product for Virtual Inertia (VI), which is also known as synthetic inertia, and an inertia product. We adopt a reserve requirement that guarantees sufficient primary frequency response reserve to adequately arrest frequency decline in response to a large generator outage within a certain margin. We place this reserve requirement into a real-time co-optimization problem, derive prices for each product and analyze the incentives provided to IBRs.
Modern distribution systems can accommodate different topologies through controllable tie lines for increasing the reliability of the system. Estimating the prevailing circuit topology or configuration is of particular importance at the substation for different applications to properly operate and control the distribution system. One of the applications of circuit configuration estimation is adaptive protection. An adaptive protection system relies on the communication system infrastructure to identify the latest status of power. However, when the communication links to some of the equipment are outaged, the adaptive protection system may lose its awareness over the status of the system. Therefore, it is necessary to estimate the circuit status using the available healthy communicated data. This paper proposes the use of machine learning algorithms at the substation to estimate circuit configuration when the communication to the tie breakers is compromised. Doing so, the adaptive protection system can identify the correct protection settings corresponding to the estimated circuit topology. The effectiveness of the proposed approach is verified on IEEE 123 bus test system.
We examine coupling into azimuthal slots on an infinite cylinder with a infinite length interior cavity operating both at the fundamental cavity modal frequencies, with small slots and a resonant slot, as well as higher frequencies. The coupling model considers both radiation on an infinite cylindrical exterior as well as a half space approximation. Bounding calculations based on maximum slot power reception and interior power balance are also discussed in detail and compared with the prior calculations. For higher frequencies limitations on matching are imposed by restricting the loads ability to shift the slot operation to the nearest slot resonance; this is done in combination with maximizing the power reception as a function of angle of incidence. Finally, slot power mismatch based on limited cavity load quality factor is considered below the first slot resonance.
Geothermal energy has been underutilized in the U.S., primarily due to the high cost of drilling in the harsh environments encountered during the development of geothermal resources. Drilling depths can approach 5,000 m with temperatures reaching 170 C. In situ geothermal fluids are up to ten times more saline than seawater and highly corrosive, and hard rock formations often exceed 240 MPa compressive strength. This combination of extreme conditions pushes the limits of most conventional drilling equipment. Furthermore, enhanced geothermal systems are expected to reach depths of 10,000 m and temperatures more than 300 °C. To address these drilling challenges, Sandia developed a proof-of-concept tool called the auto indexer under an annual operating plan task funded by the Geothermal Technologies Program (GTP) of the U.S. Department of Energy Geothermal Technologies Office. The auto indexer is a relatively simple, elastomer-free motor that was shown previously to be compatible with pneumatic hammers in bench-top testing. Pneumatic hammers can improve penetration rates and potentially reduce drilling costs when deployed in appropriate conditions. The current effort, also funded by DOE GTP, increased the technology readiness level of the auto indexer, producing a scaled prototype for drilling larger diameter boreholes using pneumatic hammers. The results presented herein include design details, modeling and simulation results, and testing results, as well as background on percussive hammers and downhole rotation.
Using the power balance method we estimate the maximum electric field on a conducting wall of a cavity containing an interior structure supporting eccentric coaxial modes in the frequency regime where the resonant modes are isolated from each other.
