Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Distributed estimation”

Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 865 records · Page 48

A Two-Step Time-Series Data Clustering Method for Building-Level Load Profile: Preprint

Residential and commercial buildings have huge potential to contribute value to improve grid resilience by participating grid services. To reveal the significant value, it is critical to estimate the grid service capability from these buildings. Unlike the large-scale distributed energy resources such as wind and solar farms, those buildings need to participate grid services in aggregation, not by individual. Therefore, it is important to appropriately group buildings for aggregation. In this paper, we develop a load profile clustering method to classify the building-level load profiles for grid service capability estimation. In our two-step clustering approach, we first calculate the total load consumption for each building, clustering the load profiles based on energy consumption level. Then, we further cluster the load profiles in each energy cluster based on the load shape. The parameter selection for each clustering step is discussed. The proposed method is applied on actual building-level load profiles, and the results have proved the effectiveness of this method.

advanced metering infrastructure (AMI)↗

Clutter Distributions for Tomographic Image Standardization in Ground-Penetrating Radar

Multistatic ground-penetrating radar (GPR) signals can be imaged tomographically to produce 3-D distributions of image intensities. In the absence of objects of interest, these intensities can be considered to be estimates of clutter. These clutter intensities spatially vary over several orders of magnitude and vary across different arrays, which makes a direct comparison of these raw intensities difficult. However, by gathering statistics on these intensities and their spatial variation, a variety of metrics can be determined. In this study, the clutter distribution is found to fit better to a two-parameter Weibull distribution than Gaussian or log-normal distributions. Based on the spatial variation of the two Weibull parameters, scale and shape, more information may be gleaned from these data. How well the GPR array is illuminating various parts of the ground, in depth and cross track, may be determined from the spatial variation of the Weibull scale parameter, which may in turn be used to estimate an effective attenuation coefficient in the soil. The transition in depth from clutter- to noise-limited conditions (which is one possible definition of GPR penetration depth) can be estimated from the spatial variation of the Weibull shape parameter. Lastly, the underlying clutter distributions also provide an opportunity to standardize image intensities to determine when a statistically significant deviation from background (clutter) has occurred, which is convenient for buried threat detection algorithm development that needs to be robust across multiple different arrays.

42 ENGINEERING↗

Numerical characterization of support recovery in sparse regression with correlated design

Sparse regression is employed in diverse scientific settings as a feature selection method. A pervasive aspect of scientific data is the presence of correlations between predictive features. These correlations hamper both feature selection and estimation and jeopardize conclusions drawn from estimated models. On the other hand, theoretical results on sparsity-inducing regularized regression have largely addressed conditions for selection consistency via asymptotics, and disregard the problem of model selection, whereby regularization parameters are chosen. In this numerical study, we address these issues through exhaustive characterization of the performance of several regression estimators, coupled with a range of model selection strategies. These estimators and selection criteria were examined across correlated regression problems with varying degrees of signal to noise, distributions of non-zero model coefficients, and model sparsity. Our results reveal a fundamental tradeoff between false positive and false negative control in all regression estimators and model selection criteria examined. Additionally, we numerically explore a transition point modulated by the signal-to-noise ratio and spectral properties of the design covariance matrix at which the selection accuracy of all considered algorithms degrades. Overall, we find that SCAD coupled with BIC or empirical Bayes model selection performs the best feature selection across the regression problems considered.

97 MATHEMATICS AND COMPUTING↗

Estimating the color lifetime of energetic quarks

Using a simple geometric framework with a realistic nuclear density distribution, we fit published HERMES data to determine fundamental properties of hadronization using the nuclear medium as a spatial analyzer. Our approach uses a fit to the transverse momentum broadening observable and the hadronic multiplicity ratio; the simultaneous fit to two different observables strongly constrains the outcome. Using the known sizes of the target nuclei, we extract the color lifetime, finding a z h -dependent range of values from 2 to 8 fm/c for these data. We also extract estimates for the qˆ transport coefficient characterizing the strength of the interaction between the quark and the cold nuclear medium, finding an average value of 0.035±0.011 GeV 2 /fm. With a three-parameter model we obtain satisfactory fits to the data with a goodness-of-fit parameter χ 2 /dof of 1.1 or less. In a secondary fit of the results from that model we independently find a value for the Lund String Model string tension of 1.00±0.05 GeV/fm. We evaluated the sensitivity for extracting quark energy loss and effective in-medium hadronic cross sections using four-parameter variants of the model, finding large uncertainties in both cases. Our results suggest that hadronic interaction of forming hadrons in the nuclear medium is the primary dynamical cause of meson attenuation in the HERMES data, with quark energy loss playing a more minor role.

Electron-ion collider↗

A practical and efficient approach for Bayesian quantum state estimation

Bayesian inference is a powerful paradigm for quantum state tomography, treating uncertainty in meaningful and informative ways. Yet the numerical challenges associated with sampling from complex probability distributions hampers Bayesian tomography in practical settings. In this article, we introduce an improved, self-contained approach for Bayesian quantum state estimation. Leveraging advances in machine learning and statistics, our formulation relies on highly efficient preconditioned Crank–Nicolson sampling and a pseudo-likelihood. We theoretically analyze the computational cost, and provide explicit examples of inference for both actual and simulated datasets, illustrating improved performance with respect to existing approaches.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

Study of Quasielastic Barrier Distributions as a Step towards the Synthesis of Superheavy Elements with Hot Fusion Reactions

The excitation functions for quasielastic scattering of Ne 22 + Cm 248 , Mg 26 + Cm 248 , and Ca 48 + U 238 are measured using a gas-filled recoil ion separator. The quasielastic barrier distributions are extracted for these systems and are compared with coupled-channel calculations. The results indicate that the barrier distribution is affected dominantly by deformation of the actinide target nuclei, but also by vibrational or rotational excitations of the projectile nuclei, as well as neutron transfer processes before capture. From a comparison between the experimental barrier distributions and the evaporation residue cross sections for Sg ( Z = 106 ), Hs (108), Cn (112), and Lv (116), it is suggested that the hot fusion reactions take advantage of a compact collision, where the projectile approaches along the short axis of a prolately deformed nucleus. Overall, a new method is proposed to estimate the optimum incident energy to synthesize unknown superheavy nuclei using the barrier distribution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Non-parametric projections of national income distribution consistent with the Shared Socioeconomic Pathways

Understanding and projecting income distributions within countries and regions is important to understanding consumption trends and the distributional consequences of climate impacts and responses. Several global, country-level projections of income distribution are available but most project only the Gini coefficient (a summary statistic of the distribution) or utilize the Gini along with the assumption of a lognormal distribution. We test the lognormal assumption and find that it typically underestimates income in the highest deciles and over-estimates it in others. We find that a new model based on two principal components of national time series data for income distribution provides a better fit to the data for all deciles, especially for the highest and lowest. We also construct a projection model in which the first principal component is driven by the Gini coefficient and the second captures deviations from this relationship. We use the model to project income distribution by decile for all countries for the five shared socioeconomic pathways. We find that inequality is consistently higher than projections based on the Gini and the lognormal functional form, with some countries reaching ratios of the highest to lowest income deciles that are almost three times their value using the lognormal assumption.

29 ENERGY PLANNING, POLICY, AND ECONOMY↗

Inverter Reliability Estimation for Advanced Inverter Functionality

In the near future, grid operators are expected to regularly use advanced distributed energy resource (DER) functions, defined in IEEE 1547-2018, to perform a range of grid-support operations. Many of these functions adjust the active and reactive power of the device through commanded or autonomous modes, which will produce new stresses on the grid-interfacing power electronics components, such as DC/AC inverters. In previous work, multiple DER devices were instrumented to evaluate additional component stress under multiple reactive power setpoints. We utilize quasi-static time-series simulations to determine voltage-reactive power mode (volt-var) mission profile of inverters in an active power system. Mission profiles and loss estimates are then combined to estimate the reduction of the useful life of inverters from different reactive power profiles. It was found that the average lifetime reduction was approximately 0.15% for an inverter between standard unity power factor operation and the IEEE 1547 default volt-var curve based on thermal damage due to switching in the power transistors. For an inverter with an expected 20-year lifetime, the 1547 volt-var curve would reduce the expected life of the device by 12 days. This framework for determining an inverter's useful life from experimental and modeling data can be applied to any failure mechanism and advanced inverter operation.

component degradation↗

Dynamic Modeling, Parameter Estimation, and Data Reconciliation of a Supercritical Pulverized Coal-Fired Boiler

A high-fidelity, spatially distributed dynamic model of a supercritical pulverized coal-fired boiler was developed for calculating transient thermal profiles of the boiler tubes, which are often not measured or unmeasurable in industrial settings due to the harsh operating conditions. In this work, rigorous models of the subcritical and supercritical steam properties, which can be highly nonlinear, were used to capture the transition between two-phase and one-phase flow, respectively, as supercritical boilers frequently transition through the critical point during load-following operation. The transient boiler model was also used for data reconciliation and parameter estimation for the purpose of model validation against industrial data. An approach was developed to reduce the size of large-scale reconciliation problem by focusing on the implementation of bias terms for the reconciled variables, which can be estimated in place of the reconciled variables themselves. In addition, a functional approximation of the reconciliation problem was defined so that it can be more simply optimized. Combining this functional approximation optimization with reformulation of the reconciliation problem considering measurement biases can significantly improve the tractability of many such large-scale problems. The validated boiler model was then used to study transient responses under load-following operation with specific attention on the calculation of unmeasurable variables.

42 ENGINEERING↗

Simulation Evaluation of a Large-Scale Implementation of Virtual-Phase Link-Based Model Predictive Control

Traffic congestion is a serious problem in the US, and traffic signal control is one of the effective solutions to congestion. Previous research on model predictive control (MPC)-based traffic signal control showed substantial benefits over conventional methods. This study focused on implementing MPC over a large-scale network with complex intersections and the impact of cycle length, network size, and imperfect state estimation on performances. This study implemented a virtual phase link (VPL)-based model predictive control method which used the number of vehicles in each VPL as input state variables and was suitable for National Electrical Manufacturing Association (NEMA) ring-barrier control. To test the impact of network size, the performance of distributed MPC (36 intersections in the network are divided into five subnetworks) was compared with that of MPC over the full network for a set of cycle lengths. To test the impact of imperfect state estimation, we synthetically infused estimation error and developed two scenarios, MPC-error and MPC-error narrow, which had higher and lower estimation errors, respectively. The performance of these MPC methods was compared with that of the existing time-of-day (TOD) method and an offline method that used Webster's method for split and MULTIBAND for cycle length and offset optimization. Trajectory and linkwise signal performance measures were collected from the simulation to evaluate performance. The distributed MPC method with perfect state estimation had the lowest delay and highest energy efficiency of all the methods. The performance of MPC decreased as the prediction inaccuracy increased. MPC-error had 7% and 11% more delay than MPC-error narrow in the morning and evening peaks, respectively. Overall, simulation results suggest that even with imperfect state estimation, MPC methods will outperform offline methods significantly.

large-scale simulation↗

Improved multifidelity Monte Carlo estimators based on normalizing flows and dimensionality reduction techniques

Here, we study the problem of multifidelity uncertainty propagation for computationally expensive models. In particular, we consider the general setting where the high-fidelity and low-fidelity models have a dissimilar parameterization both in terms of number of random inputs and their probability distributions, which can be either known in closed form or provided through samples. We derive novel multifidelity Monte Carlo estimators which rely on a shared subspace between the high-fidelity and low-fidelity models where the parameters follow the same probability distribution, i.e., a standard Gaussian. We build the shared space employing normalizing flows to map different probability distributions into a common one, together with linear and nonlinear dimensionality reduction techniques, active subspaces and autoencoders, respectively, which capture the subspaces where the models vary the most. We then compose the existing low-fidelity model with these transformations and construct modified models with an increased correlation with the high-fidelity model, which therefore yield multifidelity estimators with reduced variance. A series of numerical experiments illustrate the properties and advantages of our approaches.

97 MATHEMATICS AND COMPUTING↗

Compressed Natural Gas Component Leak Frequency Estimation

The frequency of unintended releases in a compressed natural gas system is an important aspect of the system quantitative risk assessment. The frequencies for possible release scenarios, along with engineering models, are utilized to quantify the risks for compressed natural gas facilities. This report documents component leakage frequencies representative of compressed natural gas components that were estimated as a function of the normalized leak size. A Bayesian statistical method was used which results in leak frequency distributions for each component which represent variation and uncertainty in the leak frequency. The analysis shows that there is high uncertainty in the estimated leak frequencies due to sparsity in compressed natural gas data. These leak frequencies may still be useful in compressed natural gas system risk assessments, as long as this high uncertainty is acknowledged and considered appropriately.

03 NATURAL GAS↗

A Two-Step Time-Series Data Clustering Method for Building-Level Load Profile

Residential and commercial buildings have huge potential to contribute value to improve grid resilience by participating grid services. To reveal the significant value, it is critical to estimate the grid service capability from these buildings. Unlike the large-scale distributed energy resources such as wind and solar farms, those buildings need to participate grid services in aggregation, not by individual. Therefore, it is important to appropriately group buildings for aggregation. The load profiles in the same group will have similar characteristics at the same time step, so grid operators can send the grid service signal to the customer group with a higher chance to respond at that time step. In this paper, we develop a load profile clustering method to classify the building-level load profiles for grid service capability estimation and the results have proved the effectiveness of this method.

AMI↗

Predicting climate change impacts on poikilotherms using physiologically guided species abundance models

Poikilothermic animals comprise most species on Earth and are especially sensitive to changes in environmental temperatures. Species conservation in a changing climate relies upon predictions of species responses to future conditions, yet predicting species responses to climate change when temperatures exceed the bounds of observed data is fraught with challenges. We present a physiologically guided abundance (PGA) model that combines observations of species abundance and environmental conditions with laboratory-derived data on the physiological response of poikilotherms to temperature to predict species geographical distributions and abundance in response to climate change. The model incorporates uncertainty in laboratory-derived thermal response curves and provides estimates of thermal habitat suitability and extinction probability based on site-specific conditions. We show that temperature-driven changes in distributions, local extinction, and abundance of cold, cool, and warm-adapted species vary substantially when physiological information is incorporated. Notably, cold-adapted species were predicted by the PGA model to be extirpated in 61% of locations that they currently inhabit, while extirpation was never predicted by a correlative niche model. Failure to account for species-specific physiological constraints could lead to unrealistic predictions under a warming climate, including underestimates of local extirpation for cold-adapted species near the edges of their climate niche space and overoptimistic predictions of warm-adapted species.

54 ENVIRONMENTAL SCIENCES↗

LTAU-FF: Loss Trajectory Analysis for Uncertainty in atomistic Force Fields

LTAU (Loss Trajectory Analysis for Uncertainty) is a technique for estimating the uncertainty in a model's predictions for any regression task by approximating the cumulative distribution function (CDF) of the model's errors over the course of training. The approximated CDF, combined with a similarity search algorithm in the model's descriptor space, can then be used to estimate the likelihood that the model's predictions on any given test point will fall below a chosen threshold. LTAU-FF (LTAU in atomistic Force Fields) is the application of LTAU specifically for use with atomistic force fields.

Vita, JoshuaA↗

ARPA-E PNDIODES "PN Junctions by Ion Implantation" Project (Final Report)

Electricity accounts for nearly 40% of the energy generated in the United States, of which 65% is lost to electrical system inefficiency before reaching end-use sectors. Power electronics will play a growing role in this distribution and end-use, with up to 80% of electricity estimated to pass through power electronics from generation to consumption by 2030. This gives technical improvements in power electronics enormous potential to drive efficiency gains throughout the U.S. economy. Current Si-based power architectures including MOSFETs and IGBTs are limited by high losses, low switching frequency, and poor high-temperature performance. Given the material properties of GaN, devices can exhibit a given turn-on resistance and breakdown voltage with a minute fraction of the area of a Si device—all with higher switching speeds. While n-type doping of III-nitrides has reached a high level of maturity and integration into commercialized products, p-type doping presents distinct challenges that requires continued research. The success of the next generation of GaN-based power devices depends in large part on the ability to form selective-area p-type regions, which can be accomplished by ion implantation of the prevailing acceptor dopant, Mg. This approach is relatively inexpensive and is a fab-ready processing step for creation of current-generation power electronic devices including BJTs, IGBTs, and diode rectifiers. However, ion implantation induces lattice damage and creates point defects within the material, which negatively impacts dopant activation annealing. These forms of damage may be removed by high-temperature annealing, however short annealing time and higher pressures are required to prevent decomposition of the material. We approach this challenge through an innovative combination of approaches toward obtaining high-conductivity p-type GaN. By the novel employment of a gyrotron, a microwave source that can be shaped into a highly localized beam to dissipate very high power in a semiconductor target, a wafer of GaN can be annealed in short (<2 s) ‘pulses’ at temperatures of 1450 °C (2640 °F) to achieve p-type activation without appreciable surface degradation. The p-type activation of Mg-implanted GaN can be further aided by co-implantation with N, believed to enhance the removal of point defects that compensate the acceptors dopant. Elevation to such high temperatures for even a short time induces degradation of GaN-to overcome, we have implemented additional measures to retain high-quality crystalline GaN while removing implant damage and activating Mg dopant. A bilayer cap of AlN is deposited by a combination of chemical and physical vapor deposition to protect the GaN surface. To prevent decomposition of GaN loss of N during annealing, a custom vessel—to contain pressurized N while allowing injection of microwave energy—is employed. By detailed study of the interaction of various measures that are implemented to protect the sample surface (GaN), activate the implanted dopants, and mitigate damage of the lattice heating damage during annealing, we are able to achieve p-type GaN and observe recovery of the crystal and removal of substantial density of defects created during the implantation process.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Conservative Estimation of Tail Probabilities from Limited Sample Data

Several sparse-sample uncertainty quantification (UQ) methods are compared for conservative but not overly conservative estimation of small tail probabilities involving responses that lay beyond specified thresholds in the tails of probability distributions. Sixteen very differently shaped distributions (or probability density functions, PDFs) and tail probability magnitudes ranging from 10 -5 to 10 -1 are considered in order for the study to be relevant to a wide range of risk analysis and quantification of margins and uncertainty (QMU) problems. The emphasis of the study is on limited data regimes ranging from N = 2 to 20 samples, reflective of most experimental and some expensive computational situations. Relatively simple sparse-sample UQ methods tested for this regime involve statistical tolerance interval "Equivalent Normal and related "Ensemble of Normals" and "Superdistribution (SD) approaches. (The independently derived SD is effectively equivalent to the Bayesian posterior predictive distribution given the assumptions of the derivation.) The performance of the methods was generally improved for N ≥ 5 samples with a generalized Jackknife resampling technique, which determines a tail probability estimate by averaging estimates from smaller sub-samples. Several quantitative metrics for method conservatism and accuracy of tail probability estimation are used to assess and rank the methods' performance over many random trials for each test PDF and probability magnitude. A variant of Bootstrap resampling was also tried, but did not significantly improve tail probability estimates in most cases. Detailed results are presented from over 100-million tests over the above factors that provide useful granular information on which methods or combination of methods perform best in various areas of the factor space.

97 MATHEMATICS AND COMPUTING↗

Neural simulation-based inference of the neutron star equation of state directly from telescope spectra

Neutron stars provide a unique opportunity to study strongly interacting matter under extreme density conditions. The intricacies of matter inside neutron stars and their equation of state are not directly visible, but determine bulk properties, such as mass and radius, which affect the star's thermal X-ray emissions. However, the telescope spectra of these emissions are also affected by the stellar distance, hydrogen column, and effective surface temperature, which are not always well-constrained. Uncertainties on these nuisance parameters must be accounted for when making a robust estimation of the equation of state. In this study, we develop a novel methodology that, for the first time, can infer the full posterior distribution of both the equation of state and nuisance parameters directly from telescope observations. This method relies on the use of neural likelihood estimation, in which normalizing flows use samples of simulated telescope data to learn the likelihood of the neutron star spectra as a function of these parameters, coupled with Hamiltonian Monte Carlo methods to efficiently sample from the corresponding posterior distribution. Our approach surpasses the accuracy of previous methods, improves the interpretability of the results by providing access to the full posterior distribution, and naturally scales to a growing number of neutron star observations expected in the coming years.

79 ASTRONOMY AND ASTROPHYSICS↗