Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “estimating”

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 1,423 records · Page 79

A direct filter method for parameter estimation

Parameter estimation is an important research topic in data assimilation. In this paper, a novel parameter estimation method is introduced, where the parameter is considered as the state process in a nonlinear filtering problem and the state model that contains the parameter is used to construct a pseudo-observation. This approach is named the direct filter method since nonlinear filtering algorithms are used to estimate the parameter directly without estimating the state model as part of the solution in the nonlinear filtering problem. Numerical experiments are carried out to examine the effectiveness and accuracy of the direct filter method.

97 MATHEMATICS AND COMPUTING↗

Ensemble models for circuit topology estimation, fault detection and classification in distribution systems

This paper presents a methodology for simultaneous fault detection, classification, and topology estimation for adaptive protection of distribution systems. The methodology estimates the probability of the occurrence of each one of these events by using a hybrid structure that combines three sub-systems, a convolutional neural network for topology estimation, a fault detection based on predictive residual analysis, and a standard support vector machine with probabilistic output for fault classification. The input to all these sub-systems is the local voltage and current measurements. A convolutional neural network uses these local measurements in the form of sequential data to extract features and estimate the topology conditions. The fault detector is constructed with a Bayesian stage (a multitask Gaussian process) that computes a predictive distribution (assumed to be Gaussian) of the residuals using the input. Since the distribution is known, these residuals can be transformed into a Standard distribution, whose values are then introduced into a one-class support vector machine. The structure allows using a one-class support vector machine without parameter cross-validation, so the fault detector is fully unsupervised. Finally, a support vector machine uses the input to perform the classification of the fault types. All three sub-systems can work in a parallel setup for both performance and computation efficiency. In conclusion, we test all three sub-systems included in the structure on a modified IEEE123 bus system, and we compare and evaluate the results with standard approaches.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Estimation of iodine nutrition and thyroid function status in late-gestation pregnant women in the United States: Development and application of a population-based pregnancy model

Previously, a deterministic biologically-based dose-response (BBDR) pregnancy model was developed to evaluate moderate thyroid axis disturbances with and without thyroid-active chemical exposure in a near-term pregnant woman and fetus. In the current study, the existing BBDR model was adapted to include a wider functional range of iodine nutrition, including more severe iodine deficiency conditions, and to incorporate empirically the effects of homeostatic mechanisms. The extended model was further developed into a population-based model and was constructed using a Monte Carlo-based probabilistic framework. In order to characterize total (T4) and free (fT4) thyroxine levels for a given iodine status at the population-level, the distribution of iodine intake for late-gestation pregnant women in the U.S was reconstructed using various reverse dosimetry methods and available biomonitoring data. The range of median (mean) iodine intake values resulting from three different methods of reverse dosimetry tested was 196.5–219.9 μg of iodine/day (228.2–392.9 μg of iodine/day). There was minimal variation in model-predicted maternal serum T4 and ft4 thyroxine levels from use of the three reconstructed distributions of iodine intake; the range of geometric mean for T4 and fT4, was 138–151.7 nmol/L and 7.9–8.7 pmol/L, respectively. The average value of the ratio of the 97.5th percentile to the 2.5th percentile equaled 3.1 and agreed well with similar estimates from recent observations in third-trimester pregnant women in the U.S. In addition, the reconstructed distributions of iodine intake allowed us to estimate nutrient inadequacy for late-gestation pregnant women in the U.S. via the probability approach. The prevalence of iodine inadequacy for third-trimester pregnant women in the U.S. was estimated to be between 21% and 44%. Taken together, the current work provides an improved tool for evaluating iodine nutritional status and the corresponding thyroid function status in pregnant women in the U.S. This model enables future assessments of the relevant risk of thyroid hormone level perturbations due to exposure to thyroid-active chemicals at the population-level. - Highlights: • A population-based thyroid function model for pregnant women was developed. • The model was used specifically study the late-gestation pregnant women in the U.S. • The prevalence of iodine inadequacy was estimated in the sub-population studied. • Developed model well predicts trimester-specific thyroid hormone reference ranges. • The model can be further used to study thyroid perturbations at a population level.

60 APPLIED LIFE SCIENCES↗

Comparing EGF Methods for Estimating Corner Frequency and Stress Drop from P-wave Spectra

Empirical Green's functions (EGF) are widely applied to correct earthquake spectra for attenuation and other path effects in order to estimate corner frequencies and stress drops, but these source parameter estimates often exhibit poor agreement between different studies. We investigate this issue by analyzing a compact cluster of over 3000 aftershocks of the 1992 Landers earthquake. We apply and compare two different analysis and modeling methods: (1) the spectral decomposition and global EGF fitting approach, and (2) a more traditional EGF method of modeling spectral ratios. We find that spectral decomposition yields event terms that are consistent with stacks of spectral ratios for individual events, but source-parameter estimates nonetheless vary between the methods. The main source of differences comes from the modeling approach used to estimate the EGF. The global EGF-fitting approach suffers from parameter tradeoffs among the absolute stress drop, the stress-drop scaling with moment, and the high-frequency falloff rate, but has the advantage that the relative spectral shapes and stress drops among the different events in the cluster are well-resolved even if their absolute levels are not. The spectral-ratio approach solves for a different EGF for each target event without imposing any constraint on the corner frequency, f c , of the smaller events, and so can produce biased results for target-event f c . Placing constraints on the small-event f c improves the performance of the spectral-ratio method and allows the two methods to yield very similar results.

58 GEOSCIENCES↗

Kernel Density Estimation of Reaction Rates in Neutron Transport Simulations of Nuclear Reactors

Kernel density estimators (KDEs) are applied to estimate neutron scalar flux and reaction rate densities in Monte Carlo neutron transport simulations of heterogeneous nuclear reactors in continuous energy. The mean free path (MFP) KDE is introduced in order to handle the issues that arise from estimating the discontinuous reaction rate densities at material interfaces. Results show the MFP KDE is more accurate at estimating reaction rates compared with previous KDE formulations. An approximate MFP (aMFP) KDE is introduced to circumvent several practical issues presented by the MFP KDE. A volume-averaged KDE is derived and used to determine the bias introduced by the aMFP KDE. A KDE is formulated for cylindrical coordinates to better represent the geometry and capture the physics in two-dimensional reactor physics problems. The results here indicate that the cylindrical MFP KDE and cylindrical aMFP KDE are accurate tools for capturing reaction rates in heterogeneous reactor physics problems in continuous energy, with local biases of less than 1%.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Variance Estimation in Monte Carlo Eigenvalue Simulations Using Spectral Analysis Method

Monte Carlo (MC) simulation is widely used to solve the eigenvalue form of the Boltzmann transport equation that mathematically represents the neutron transport process through complex multiplying (fissionable) systems. Monte Carlo eigenvalue simulation starts with an assumed fission source distribution and uses the fission sites from the previous iteration (cycle) as the starting source in the current iteration. Important system parameters (MC tallies) such as fuel pin-power distribution are estimated over several cycles after the convergence of the fission source distribution to a stationary distribution. However, the MC fission source iteration algorithm that uses fission source sites from the previous cycle introduces a cycle-to-cycle correlation. Monte Carlo simulations that do not account for the cycle-to-cycle correlation systematically underestimate the variance of the estimated system parameters (sample mean). This paper presents the relationship between the spectral density in the frequency domain at frequency zero and the variance of the sample mean. This paper introduces a novel method in the frequency domain for the MC variance estimation. For the three test problems used in this paper, researchers have observed that the new method results in an improvement of more than one order of magnitude to the standard deviation of the sample mean. The new method also compares favorably with the previously introduced batch, bootstrap, and covariance-adjusted methods when applied to the three test problems investigated in this paper. This new method does not require modification of the MC eigenvalue algorithm (power iteration), is code agnostic, and is therefore easy to use when implementing in any existing MC code. In conclusion, the new estimate can be calculated without saving tally results of all active/stationary cycles.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Estimating Code Biases for Criticality Safety Applications with Few Relevant Benchmarks

Criticality safety analyses rely on the availability of relevant benchmark experiments to determine justifiable margins of subcriticality. When a target application lacks neutronically similar benchmark experiments, validation studies must provide justification to the regulator that the impact of modeling and simulation limitations is well understood for the application and often must provide additional subcritical margin to ensure safe operating conditions. This study estimated the computational bias in the critical eigenvalue for several criticality safety applications supported by only a few relevant benchmark experiments. The accuracy of the following three methods for predicting computational biases was evaluated: the Upper Subcritical Limit STATisticS (USLSTATS) trending analysis method; the Whisper nonparametric method; and TSURFER, which is based on the generalized linear least-squares technique. Here, these methods were also applied to estimate computational biases and recommended upper subcriticality limits for several critical experiments with known biases and for several cases from a blind benchmark study. The methods are evaluated based on both the accuracy of their predicted computation bias and upper subcriticality limit estimates, as well as on the consistency of the methods’ estimates, as the model parameters, covariance data libraries, and set of available benchmark data were varied. Data assimilation methods typically have not been used for criticality safety licensing activities, and this study explores a methodology to address concerns regarding the reliability of such methods in criticality safety bias prediction applications.

98 NUCLEAR DISARMAMENT, SAFEGUARDS, AND PHYSICAL P↗

On MCNP Stochastic Volume Estimation Normalization

This paper resolves a perennial point of confusion regarding the source-weighting normalization factor recommended in the MCNP manual (𝜋⁢𝑟 2 ) to stochastically estimate the volume of a region within an enclosing inward-directed spherical surface source with radius 𝑟. The normalization factor arises from the relationship between a sphere’ s mean chord length, its volume, and the values estimated by MCNP track-length tallies. A brief derivation is given that relates these quantities and results in the stated normalization. The correctness of this factor is demonstrated by estimating the volume of a variety of convex and nonconvex volumes. A heuristic demonstration of how biasing the inward-directed source reduces the statistical uncertainty of the stochastic volume estimate is also given, but a rigorous analysis of this improvement is left as future work.

42 ENGINEERING↗

NONPARAMETRIC BAYESIAN ESTIMATION OF PERIODIC LIGHT CURVES

Many astronomical phenomena exhibit patterns that have periodic behavior. An important step when analyzing data from such processes is the problem of identifying the period: estimating the period of a periodic function based on noisy observations made at irregularly spaced time points. This problem is still a difficult challenge despite extensive study in different disciplines. This paper makes several contributions toward solving this problem. First, we present a nonparametric Bayesian model for period finding, based on Gaussian Processes (GPs), that does not make assumptions on the shape of the periodic function. As our experiments demonstrate, the new model leads to significantly better results in period estimation especially when the light curve does not exhibit sinusoidal shape. Second, we develop a new algorithm for parameter optimization for GP which is useful when the likelihood function is very sensitive to the parameters with numerous local minima, as in the case of period estimation. The algorithm combines gradient optimization with grid search and incorporates several mechanisms to overcome the high computational complexity of GP. Third, we develop a novel approach for using domain knowledge, in the form of a probabilistic generative model, and incorporate it into the period estimation algorithm. Experimental results validate our approach showing significant improvement over existing methods.

79 ASTRONOMY AND ASTROPHYSICS↗

Estimation of plasma properties using an extended Kalman filter with plasma global models

A physically-constrained extended Kalman filter (EKF) is applied to various zero-dimensional global models for the estimation of plasma properties using time-dependent experimental data such as the plasma density or ion flux. The capability of the EKF is demonstrated to estimate unknown system states simultaneously, such as reaction rate coefficients and the absorbed electron input power, which can be difficult, if not impossible, to measure experimentally. Global models accounting for pure argon reactions and argon-oxygen reactions are used in this work to demonstrate the ability of the filter to estimate dynamic and complex systems. Furthermore, the results obtained from the EKF plasma global model illustrate that model-data fusion techniques can be used to estimate plasma properties and processes for time-varying systems, such as pulsed discharges.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

Methods of improving brain dose estimates for internally deposited radionuclides *

The US National Council on Radiation Protection and Measurements (NCRP) convened Scientific Committee 6–12 (SC 6–12) to examine methods for improving dose estimates for brain tissue for internally deposited radionuclides, with emphasis on alpha emitters. This Memorandum summarises the main findings of SC 6–12 described in the recently published NCRP Commentary No. 31, ‘Development of Kinetic and Anatomical Models for Brain Dosimetry for Internally Deposited Radionuclides’. The Commentary examines the extent to which dose estimates for the brain could be improved through increased realism in the biokinetic and dosimetric models currently used in radiation protection and epidemiology. A limitation of most of the current element-specific systemic biokinetic models is the absence of brain as an explicitly identified source region with its unique rate(s) of exchange of the element with blood. The brain is usually included in a large source region called Other that contains all tissues not considered major repositories for the element. In effect, all tissues in Other are assigned a common set of exchange rates with blood. A limitation of current dosimetric models for internal emitters is that activity in the brain is treated as a well-mixed pool, although more sophisticated models allowing consideration of different activity concentrations in different regions of the brain have been proposed. Here case studies for 18 internal emitters indicate that brain dose estimates using current dosimetric models may change substantially (by a factor of 5 or more), or may change only modestly, by addition of a sub-model of the brain in the biokinetic model, with transfer rates based on results of published biokinetic studies and autopsy data for the element of interest. As a starting place for improving brain dose estimates, development of biokinetic models with explicit sub-models of the brain (when sufficient biokinetic data are available) is underway for radionuclides frequently encountered in radiation epidemiology. A longer-term goal is development of coordinated biokinetic and dosimetric models that address the distribution of major radioelements among radiosensitive brain tissues.

61 RADIATION PROTECTION AND DOSIMETRY↗

Adaptive Linear State Estimation for Unbalanced Distribution System

The inclusion of PMU functionality in distribution relays enables the implementation of a linear state estimator (LSE) in Distribution Systems (DS). However, the unbalanced topology and phase coupling in distribution lines necessitate modifications to the LSE formulation. Additionally, the higher fault frequency in distribution systems requires a state estimation approach that is resilient to contingencies. This work proposes an adaptive linear state estimation algorithm tailored for unbalanced distribution systems with single-phase and two-phase laterals. Furthermore, a modified Optimal PMU Placement (OPP) strategy is introduced to ensure full observability in distribution systems with single-phase and two-phase buses. To maintain adaptability to topology changes, the state estimator incorporates circuit breaker status data provided by PMUs, ensuring robust performance during topology changes triggered by faults. The performance of the algorithm is verified on the IEEE 13-bus, 34-bus, and 123-bus systems.

PMUs↗

Efficient Network Partitioning: Application for Decentralized State Estimation in Power Distribution Grids

Increase in the proliferation of distributed energy resources require real-time situational awareness for efficient grid operations. State estimation plays an important role for the real-time control and management of the power grid. As the sensing infrastructure grows, aggregating and handling high volumes of data at a centralized location is extremely difficult. To address this challenge, this paper first proposes a novel and efficient hier-archical spectral clustering-based network partitioning algorithm followed by a decentralized compressive sensing (DCS)-based state estimation. The applicability of the proposed network partitioning algorithm is tested on an IEEE 123-bus network, an IEEE 8,500-node system, and a 6,000+ node distribution network. The results shows that the proposed approach efficiently divides the network into multiple sub-networks with the minimum number of edge connections among the neighbors. Then, we perform DCS-based state estimation on the 6,000+ node distribution network after dividing the network into 18 optimal partitions. Simulation results show that the DCS-based state estimation recovers the system states with high accuracy and low complexity.

alternating direction method of multipliers↗

Missing-Data Nonparametric Coherency Estimation

Chave recently proposed an estimator for multitaper spectral density where the time series contains missing values. In this article we generalize this technique to a multitaper estimator of coherence and phase and show that one can also obtain bootstrapped confidence intervals. Additionally, we give two examples. The first is a toy example in which the true coherence is known. In the second example we show that the multitaper missing-data coherence estimator computed on real data with a single gap comprising 11% of the data outperforms the Daniell-smoothed coherence estimator where there are no gaps. The case where the two time series have different missing indices is also discussed.

42 ENGINEERING↗

Recent Development of Frequency Estimation Methods for Future Smart Grid

The frequency estimated by the Phasor Measurement Unit (PMU) is a critical index of power system status and supports many smart grid applications. The future smart grid features high penetration of renewables and more fast-moving power electronics inverters but raises challenges to the reliable frequency estimation. This article presents three methods to address these challenges. First, an enhanced zero-crossing algorithm was developed to track the fast-changing frequency in system dynamics. Second, we propose a technology that can tolerate the system transient and suppress the outliers. Third, an algorithm was developed to export high time-resolution frequency estimations with minimum computational effort. All of the proposed methods are realized in hardware and compared with classical frequency estimation methods. The testing results indicate that the proposed methods have excellent performance. They can be used in future PMUs and provide reliable and high time resolution data for smart grid applications.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Online Output-Based Inertia Estimation of Modern Power Systems

The overall inertia of modern power systems is drastically impacted by the penetration level of converter-based resources (CBRs). The intermittent output of CBRs and their dependence on the weather conditions affirm the need for a realtime inertia estimation approach. This can help the transmission system operators(TSOs) to be continuously aware of the changes in the system inertia and to take suitable control actions. This paper is a step forward to a totally data-driven framework for power system inertia estimation using the available measurements from the phasor measurement units (PMUs). This inertia estimation is done based on the spectral analysis of the so-called Koopman linear operator. Koopman operator considers the relationship between the system's inertia and its dynamical response to slight disturbances in the load. Simulation results show that Koopman operator can approximate the system behavior in a reconstructed linear space hence inertia constant is estimated accurately.

data-driven↗

Online State Estimation for Time-Varying Systems

The paper investigates the problem of estimating the state of a time-varying system with a linear measurement model; in particular, the paper considers the case where the number of measurements available can be smaller than the number of states. In lieu of a batch linear least-squares (LS) approach well-suited for static networks, where a sufficient number of measurements could be collected to obtain a full-rank design matrix the paper proposes an online algorithm to estimate the possibly time-varying state by processing measurements as and when available. The design of the algorithm hinges on a generalized LS cost augmented with a proximal-point-type regularization. With the solution of the regularized LS problem available in closed-form, the online algorithm is written as a linear dynamical system where the state is updated based on the previous estimate and based on the new available measurements. Conditions under which the algorithmic steps are in fact a contractive mapping are shown, and bounds on the estimation error are derived for different noise models. Numerical simulations are provided to corroborate the analytical findings.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Data-Driven Day-Ahead PV Estimation Using Autoencoder-LSTM and Persistence Model

Inherent variability in photovoltaic (PV) and associated impacts on power systems is a challenging problem for both the PV owners and the grid operators. Existing statistical and machine learning algorithms typically work well for weather conditions similar to historical data. Furthermore, uncertain weather conditions pose a great challenge to the estimation accuracy of the estimation models. With the enhanced integration of intelligent electronic devices and the realization of associated automation in the power grid, renewable energy data is becoming more accessible, which can be utilized by deep learning models and improve the PV power generation estimation accuracy. In this paper, a hybrid deep learning model driven by external weather data is proposed to do day-ahead PV output forecasting at 15-minute-interval. The proposed model is motivated by the recent advancement of Long-Short-Term-Memory (LSTM) networks and AutoEncoder (AE), which estimates uncertainties in sequence while making the prediction for complex weather conditions. Meanwhile, the persistence model (PM) is used to predict continuous sunny weather conditions. The forecasting result is validated with data from multiple locations

42 ENGINEERING↗