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 1,081 records · Page 60

Revealing the Statistics of Extreme Events Hidden in Short Weather Forecast Data

Extreme weather events have significant consequences, dominating the impact of climate on society. While high-resolution weather models can forecast many types of extreme events on synoptic timescales, long-term climatological risk assessment is an altogether different problem. A once-in-a-century event takes, on average, 100 years of simulation time to appear just once, far beyond the typical integration length of a weather forecast model. Therefore, this task is left to cheaper, but less accurate, low-resolution or statistical models. But there is untapped potential in weather model output: despite being short in duration, weather forecast ensembles are produced multiple times a week. Integrations are launched with independent perturbations, causing them to spread apart over time and broadly sample phase space. Collectively, these integrations add up to thousands of years of data. We establish methods to extract climatological information from these short weather simulations. Using ensemble hindcasts by the European Center for Medium-range Weather Forecasting archived in the subseasonal-to-seasonal (S2S) database, we characterize sudden stratospheric warming (SSW) events with multi-centennial return times. Consistent results are found between alternative methods, including basic counting strategies and Markov state modeling. By carefully combining trajectories together, we obtain estimates of SSW frequencies and their seasonal distributions that are consistent with reanalysis-derived estimates for moderately rare events, but with much tighter uncertainty bounds, and which can be extended to events of unprecedented severity that have not yet been observed historically. These methods hold potential for assessing extreme events throughout the climate system, beyond this example of stratospheric extremes.

58 GEOSCIENCES↗

Estimation of Snow Parameters from Dual-Wavelength Airborne Radar

Estimation of snow characteristics from airborne radar measurements would complement In-situ measurements. While In-situ data provide more detailed information than radar, they are limited in their space-time sampling. In the absence of significant cloud water contents, dual-wavelength radar data can be used to estimate 2 parameters of a drop size distribution if the snow density is assumed. To estimate, rather than assume, a snow density is difficult, however, and represents a major limitation in the radar retrieval. There are a number of ways that this problem can be investigated: direct comparisons with in-situ measurements, examination of the large scale characteristics of the retrievals and their comparison to cloud model outputs, use of LDR measurements, and comparisons to the theoretical results of Passarelli(1978) and others. In this paper we address the first approach and, in part, the second.

Liao, Liang↗

Understanding Climate Policy Data Needs. NASA Carbon Monitoring System Briefing: Characterizing Flux Uncertainty, Washington D.C., 11 January 2012

Climate policy in the United States is currently guided by public-private partnerships and actions at the local and state levels. This mitigation strategy is made up of programs that focus on energy efficiency, renewable energy, agricultural practices and implementation of technologies to reduce greenhouse gases. How will policy makers know if these strategies are working, particularly at the scales at which they are being implemented? The NASA Carbon Monitoring System (CMS) will provide information on carbon dioxide fluxes derived from observations of earth's land, ocean and atmosphere used in state of the art models describing their interactions. This new modeling system could be used to assess the impact of specific policy interventions on CO2 reductions, enabling an iterative, results-oriented policy process. In January of 2012, the CMS team held a meeting with carbon policy and decision makers in Washington DC to describe the developing modeling system to policy makers. The NASA CMS will develop pilot studies to provide information across a range of spatial scales, consider carbon storage in biomass, and improve measures of the atmospheric distribution of carbon dioxide. The pilot involves multiple institutions (four NASA centers as well as several universities) and over 20 scientists in its work. This pilot study will generate CO2 flux maps for two years using observational constraints in NASA's state-of -the-art models. Bottom-up surface flux estimates will be computed using data-constrained land and ocean models; comparison of the different techniques will provide some knowledge of uncertainty in these estimates. Ensembles of atmospheric carbon distributions will be computed using an atmospheric general circulation model (GEOS-5), with perturbations to the surface fluxes and to transport. Top-down flux estimates will be computed from observed atmospheric CO2 distributions (ACOS/GOSAT retrievals) alongside the forward-model fields, in conjunction with an inverse approach based on the CO2 model of GEOS ]Chem. The forward model ensembles will be used to build understanding of relationships among surface flux perturbations, transport uncertainty and atmospheric carbon concentration. This will help construct uncertainty estimates and information on the true spatial resolution of the top-down flux calculations. The relationship between the top-down and bottom-up flux distributions will be documented. Because the goal of NASA CMS is to be policy relevant, the scientists involved in the flux modeling pilot need to understand and be focused on the needs of the climate policy and decision making community. If policy makers are to use CMS products, they must be aware of the modeling effort and begin to design policies that can be evaluated with information. Improving estimates of carbon sequestered in forests, for example, will require information on the spatial variability of forest biomass that is far more explicit than is presently possible using only ground observations. Carbon mitigation policies being implemented by cities around the United States could be designed with the CMS data in mind, enabling sequential evaluation and subsequent improvements in incentives, structures and programs. The success of climate mitigation programs being implemented in the United States today will hang on the depth of the relationship between scientists and their policy and decision making counterparts. Ensuring that there is two-way communication between data providers and users is important for the success both of the policies and the scientific products meant to support them..

Brown, Molly E.↗

Estimating risk for earth-satellite attenuation prediction

Annual cumulative distributions of attenuation measurements and of rain-rate measurements were obtained from several locations in Europe and the United States. They were analyzed to estimate the year-to-year variability to be associated with a prediction of the expected cumulative distribution. Two models were constructed to estimate variability, an ad hoc model that summarized a number of observations when compared to model predictions and a probabilistic model that applied the ideas of order statistics to the prediction problem when the number of independent attenuation or rain-rate events in a sample year could be estimated. Based on these models, the statistical uncertainty in a model prediction may be estimated. The estimation procedure also provides an answer to the question of the number of years of observation needed to provide an estimate of the empirical distribution with a specified statistical uncertainty.

Crane, Robert K.↗

Plasma-parameter dependence of ro-vibrational temperatures for H 2 in LHD divertor

We analyzed a thousand visible spectra of Fulcher-α band measured for divertor plasmas in Large Helical Device. With a coronal model and Bayesian inference, the population distribution of hydrogen molecule in the electronic ground state was estimated. The non-thermal population distribution was recovered with a two-temperature model, which has two sets of rotational and vibrational temperatures, as well as their mixture coefficient. The lower rotational temperature significantly changes according to the plasma parameters. Its nearly linear dependence on the electron density was found, which is consistent with previous works. The lower vibrational temperature also shows a small density dependence, as reported by a previous work. On the other hand, the higher rotational and vibrational temperatures as well as the mixture coefficient only show slight changes over the broad range of plasma parameters. These population parameters show a significant correlation; with higher electron density, all the temperatures and the fraction of the higher-temperature component increase simultaneously. In conclusion, this suggests that the electron impact plays an important role to determine the population distribution.

42 ENGINEERING↗

Determination of aquatic hazard concentrations for the oil spill response product class of surface washing agents using species sensitivity distributions

Surface washing agents (SWAs) are a diverse class of oil spill response products intended to facilitate removal of stranded oil from shorelines. This class of agents has high application rates relative to other categories of spill response products, but global toxicity data is generally limited to two standard test species: inland silverside and mysid shrimp. Here, in this study, we provide a framework to maximize the utility of limited toxicity data across a product class. To characterize species sensitivity to SWAs, the toxicity of three agents spanning a range of chemical and physical properties were tested in eight species. The relative sensitivity of mysids shrimp and inland silversides as surrogate test organisms was determined. Toxicity normalized species sensitivity distributions (SSDn) were used to estimate fifth centile hazard concentration (HC5) values for SWAs with limited toxicity data. Chemical toxicity distributions (CTD) of SWA HC5 values were used to compute a fifth centile chemical hazard distribution (HD5) to provide a more comprehensive assessment of hazard across a spill response product class with limited toxicity data than traditional single species or single agent approaches can give.

63 RADIATION, THERMAL, AND OTHER ENVIRON. POLLUTAN↗

The Poisson tensor completion non-parametric differential entropy estimator

We introduce the Poisson tensor completion (PTC) estimator, a non-parametric differential entropy estimator. The PTC estimator leverages inter-sample relationships to compute a low-rank Poisson tensor decomposition of the frequency histogram. Our crucial observation is that the histogram bins are an instance of a space partitioning of counts and thus can be identified with a spatial Poisson process. The Poisson tensor decomposition leads to a completion of the intensity measure over all bins—including those containing few to no samples—and leads to our proposed PTC differential entropy estimator. A Poisson tensor decomposition models the underlying distribution of the count data and guarantees non-negative estimated values and so can be safely used directly in entropy estimation. Our estimator is the first tensor-based estimator that exploits the underlying spatial Poisson process related to the histogram explicitly when estimating the probability density with low-rank tensor decompositions for the purpose of tensor completion. Furthermore, we demonstrate that our PTC estimator is a substantial improvement over standard histogram-based estimators for sub-Gaussian probability distributions because of the concentration of norm phenomenon.

42 ENGINEERING↗

Stochastic Filtering of Rain Profiles Using Radar, Surface-Referenced Radar, or Combined Radar-Radiometer Measurements

This paper describes a computationally efficient nearly optimal Bayesian algorithm to estimate rain (and drop size distribution) profiles, given a radar reflectivity profile at a single attenuating wavelength. In addition to estimating the averages of all the mutually ambiguous combinations of rain parameters that can produce the data observed, the approach also calculates the n-ns uncertainty in its estimates (this uncertainty thus quantifies "the amount of ambiguity" in the "solution"). The paper also describes a more general approach that can make estimates based on a radar reflectivity profile together with an approximate measurement of the path-integrated attenuation, or a radar reflectivity profile and a set of passive microwave brightness temperatures. This more general "combined" algorithm is currently being adapted for the Tropical Rainfall Measuring Mission.

Haddad, Ziad S.↗

Probabilistic Study Conducted on Sensor-Based Engine Life Calculation

Turbine engine life management is a very complicated process to ensure the safe operation of an engine subjected to complex usage. The challenge of life management is to find a reasonable compromise between the safe operation and the maximum usage of critical parts to reduce maintenance costs. The commonly used "cycle count" approach does not take the engine operation conditions into account, and it oversimplifies the calculation of the life usage. Because of the shortcomings, many engine components are regularly pulled for maintenance before their usable life is over. And, if an engine has been running regularly under more severe conditions, components might not be taken out of service before they exceed their designed risk of failure. The NASA Glenn Research Center and its industrial and academic partners have been using measurable parameters to improve engine life estimation. This study was based on the Monte Carlo simulation of 5000 typical flights under various operating conditions. First a closed-loop engine model was developed to simulate the engine operation across the mission profile and a thermomechanical fatigue (TMF) damage model was used to calculate the actual damage during takeoff, where the maximum TMF accumulates. Next, a Weibull distribution was used to estimate the implied probability of failure for a given accumulated cycle count. Monte Carlo simulations were then employed to find the profiles of the TMF damage under different operating assumptions including parameter uncertainties. Finally, probabilities of failure for different operating conditions were analyzed to demonstrate the importance of a sensor-based damage calculation in order to better manage the risk of failure and on-wing life.

Guo, Ten-Huei↗

Ionizing Radiation Environment on the International Space Station: Performance vs. Expectations for Avionics and Material

The role of structural shielding mass in the design, verification, and in-flight performance of International Space Station (ISS), in both the natural and induced orbital ionizing radiation (IR) environments, is reported. Detailed consideration of the effects of both the natural and induced ionizing radiation environment during ISS design, development, and flight operations has produced a safe, efficient manned space platform that is largely immune to deleterious effects of the LEO ionizing radiation environment. The assumption of a small shielding mass for purposes of design and verification has been shown to be a valid worst-case approximation approach to design for reliability, though predicted dependences of single event effect (SEE) effects on latitude, longitude, SEP events, and spacecraft structural shielding mass are not observed. The Figure of Merit (FOM) method over predicts the rate for median shielding masses of about 10g/cm(exp 2) by only a factor of 3, while the Scott Effective Flux Approach (SEFA) method overestimated by about one order of magnitude as expected. The Integral Rectangular Parallelepiped (IRPP), SEFA, and FOM methods for estimating on-orbit (Single Event Upsets) SEU rates all utilize some version of the CREME-96 treatment of energetic particle interaction with structural shielding, which has been shown to underestimate the production of secondary particles in heavily shielded manned spacecraft. The need for more work directed to development of a practical understanding of secondary particle production in massive structural shielding for SEE design and verification is indicated. In contrast, total dose estimates using CAD based shielding mass distributions functions and the Shieldose Code provided a reasonable accurate estimate of accumulated dose in Grays internal to the ISS pressurized elements, albeit as a result of using worst-on-worst case assumptions (500 km altitude x 2) that compensate for ignoring both GCR and secondary particle production in massive structural shielding.

Koontz, Steven L.↗

Proportion estimation using prior cluster purities

The prior distribution of CLASSY component purities is studied, and this information incorporated into maximum likelihood crop proportion estimators. The method is tested on Transition Year spring small grain segments.

Terrell, G. R.↗

Tractable Learning in Underexcited Power Grids

Estimating the structure of physical flow networks, such as power grids, is critical to secure delivery of energy. This article discusses statistical structure estimation in power grids in the “underexcited” regime, where a subset of internal nodes has zero injection fluctuations. Prior estimation algorithms based on nodal voltages fail for such grids as the voltage covariance matrix is not invertible. Here, we propose a novel topology learning algorithm for learning underexcited general networks. Our algorithm uses physics-informed conservation laws to first identify the zero-injection buses and their neighbors, and then estimates the remaining edges in the grid. We prove the asymptotic correctness of our algorithm for grids with nonadjacent internal zero-injection nodes. More important, we theoretically analyze our algorithm’s efficacy under noisy measurements, and determine bounds on maximum noise under which asymptotically correct recovery is guaranteed. Our approach is validated through simulations with voltage samples generated on test distribution grids with real injection data and nonlinear power flow models.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Three-dimensional migration behavior of juvenile salmonids in reservoirs and near dams

To acquire 3-D tracking data on juvenile salmonids, Juvenile Salmon Acoustic Telemetry System (JSATS) cabled hydrophone arrays were deployed in the forebays of two dams on the Snake River and at a mid-reach reservoir between the dams. The depth distributions of fish were estimated by statistical analyses performed on large 3-D tracking data sets from ~33,500 individual acoustic tagged yearling and subyearling Chinook salmon and juvenile steelhead at the two dams in 2012 and subyearling Chinook salmon at the two dams and the mid-reach reservoir in 2013. This research investigated the correlation between vertical migration behavior and passage routes. The depth distributions of fish within the forebays of the dams were significantly different from fish passing the mid-reach reservoir. Fish residing deeper in the forebay tended to pass the dam using deeper powerhouse routes. This difference in depth distributions indicated that the depth distribution of fish at the mid-reach reservoir was not related to behaviors of fish passing through certain routes of the adjacent dams. For fish that were detected deeper than 17.5 m in the forebays, the probability of powerhouse passage (i.e., turbine) increased significantly. Another important finding was the variation in depth distributions during dam passage associated with the diel period, especially the crepuscular periods.

59 BASIC BIOLOGICAL SCIENCES↗

Extracting the gamma-ray source-count distribution below the Fermi-LAT detection limit with deep learning

We reconstruct the extra-galactic gamma-ray source-count distribution, or dN/dS, of resolved and unresolved sources by adopting machine learning techniques. Specifically, we train a convolutional neural network on synthetic 2-dimensional sky-maps, which are built by varying parameters of underlying source-counts models and incorporate the Fermi-LAT instrumental response functions. The trained neural network is then applied to the Fermi-LAT data, from which we estimate the source count distribution down to flux levels a factor of 50 below the Fermi-LAT threshold. We perform our analysis using 14 years of data collected in the (1,10) GeV energy range. The results we obtain show a source count distribution which, in the resolved regime, is in excellent agreement with the one derived from cataloged sources, and then extends as dN/dS ~ S -2 in the unresolved regime, down to fluxes of 5 · 10 -12 cm -2 s -1 . The neural network architecture and the devised methodology have the flexibility to enable future analyses to study the energy dependence of the source-count distribution.

79 ASTRONOMY AND ASTROPHYSICS↗

Scaled Implementation of Smart Charge Management for Electric Vehicles

Smart charge management (SCM) has become a critical strategy for mitigating potential grid impacts and reducing electricity costs for all customers. This study will evaluate the economics of implementing light-duty EV SCM at scale across the United States. This work will enhance distribution system analysis by estimating SCM implementation costs, exploring viable business cases, and developing a framework for national-level applications. The primary methodology involves leveraging detailed grid modeling from a specific service territory and using spatial extrapolation techniques to generalize findings to other regions. The analysis will develop key metrics to quantify the costs and benefits of SCM, including implementation costs relative to strategy and scale, the cost of distribution system upgrades with and without SCM, and the percentage of peak-load reduction. The objective is to produce a comprehensive report and a parameterized framework that enables utilities to self-assess the value of SCM in their own service territories. This will support the development of cost-effective charging strategies, accelerate the energization of new EV chargers, and facilitate the seamless integration of EVs into the nation's power grid.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Estimating Total Electron Content Using 1,000+ GPS Receivers

A computer program uses data from more than 1,000 Global Positioning System (GPS) receivers in an Internet-accessible global network to generate daily estimates of the global distribution of vertical total electron content (VTEC) of the ionosphere. This program supersedes an older program capable of processing readings from only about 200 GPS receivers. This program downloads the data via the Internet, then processes the data in three stages. In the first stage, raw data from a global subnetwork of about 200 receivers are preprocessed, station by station, in a Kalman-filter-based least-squares estimation scheme that estimates satellite and receiver differential biases for these receivers and for satellites. In the second stage, an observation equation that incorporates the results from the first stage and the raw data from the remaining 800 receivers is solved to obtain the differential biases for these receivers. The only remaining error sources for which an account cannot be given are multipath and receiver noise contributions. The third stage is a postprocessing stage in which all the processed data are combined and used to generate new data products, including receiver differential biases and global and regional VTEC maps and animations.

Komjathy, Attila↗

DSO+T: Valuation Methodology and Economic Metrics (DSO+T Study: Volume 4)

This report summarizes a rigorous valuation analysis methodology used by the Distribution System Operator with Transactive (DSO+T) study to estimate the financial benefits and costs of adopting Transactive Energy coordination of distributed energy resources for key stakeholders (for example distribution system operators and customers). This was achieved by modeling the value exchanges between stakeholders and determining the annualized costs and revenues experienced by stakeholders, enabling the evaluation of overall impact on stakeholder’s annualized cash flow. Extensive work was conducted developing methods to estimate the operating costs of distribution system operators at a level of granularity that would allow the financial impact of implementing a transactive energy approach to be estimated. This work included developing parametric models for labor and software costs, distribution system capital and maintenance costs, growth rates, and factors to determine annualized costs of capital investments. Simulation results were used to calculate wholesale energy costs and revenues from retail sales. Valuation analysis methods were also developed for other stakeholders including customers, the transmission system operator, independent system operator, and generators. The resulting capability allows a complete mapping of the flow of financial value between stakeholders that can be directly integrated with the results of demand flexibility simulations. Example results are provided for a business-as-usual case and compared to a transactive energy case as well as to actual cost data.

24 POWER TRANSMISSION AND DISTRIBUTION↗

Estimating the Contributions to Human Error Probability from the Convolution of the Distribution of Time Available and Time Required

As part of their duties, Human Reliability Analysis must often evaluate if crews in nuclear power plants (NPPs) can complete tasks associated with a human-failure event within time limits. For example, the time required in NPP scenarios is determined by systematic and structured walkthroughs, feasibility studies, recorded times from training exercises, and interviews with experienced operators and experts. Typically, a point estimate is derived for the estimate (mean, maximum, or 95th percentile of time required). Using point-estimate values can mask the risk associated with variability among crews, plant conditions and set-up, environmental conditions, and other impact factors under which these actions are executed. While point estimates for time required and time available have served the industry well, without considering the uncertainty they could lead to biased understanding about the risk. The Integrated Human Event Analysis System - General Methodology (IDHEAS-G) model (developed by the US Nuclear Regulatory Commission, NRC) for human error probability calculates human error probability by summing two probabilities: insufficient time and cognitive error. As such, the model takes a more holistic approach by considering the full distributions for time required and time available to calculate the human error probability because the time available to complete the task is insufficient. In this study, we expand on the work of the NRC and discuss methods for estimating these time considerations. For example, for the time required, the impact of Performance Influencing Factors (PIFs) on the distribution was divided into impacts that are aleatory in nature, such as crew-to-crew variability, and those that are epistemic (i.e., the PIFs). Starting with the factors that introduce aleatory uncertainty, a first-order distribution was developed from a large set of time required (i.e., NPP task completion times) data for the range of operator actions that occur in the NPP control room under simulated accident conditions. The first-order distribution can then be adjusted to account for epistemic uncertainty using research associated with the impact of applicable PIFs on the time required. We also develop guidance for analysts to address the probability distributions for the time available. The guidance we developed on how to estimate time required and time available distributions is based on the identification of pertinent research and data, data analyses, and expert knowledge elicitation.

human error probability, human performance, time e↗