Search NASA⌕ Search

SEARCH · Search NASA

Results for “Estimation method”

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,549 records · Page 86

Do-calculus enables estimation of causal effects in partially observed biomolecular pathways

Abstract Motivation Estimating causal queries, such as changes in protein abundance in response to a perturbation, is a fundamental task in the analysis of biomolecular pathways. The estimation requires experimental measurements on the pathway components. However, in practice many pathway components are left unobserved (latent) because they are either unknown, or difficult to measure. Latent variable models (LVMs) are well-suited for such estimation. Unfortunately, LVM-based estimation of causal queries can be inaccurate when parameters of the latent variables are not uniquely identified, or when the number of latent variables is misspecified. This has limited the use of LVMs for causal inference in biomolecular pathways. Results In this article, we propose a general and practical approach for LVM-based estimation of causal queries. We prove that, despite the challenges above, LVM-based estimators of causal queries are accurate if the queries are identifiable according to Pearl’s do-calculus and describe an algorithm for its estimation. We illustrate the breadth and the practical utility of this approach for estimating causal queries in four synthetic and two experimental case studies, where structures of biomolecular pathways challenge the existing methods for causal query estimation. Availability and implementation The code and the data documenting all the case studies are available at https://github.com/srtaheri/LVMwithDoCalculus. Supplementary information Supplementary data are available at Bioinformatics online.

59 BASIC BIOLOGICAL SCIENCES↗

Uncertainty Quantification of PV Annual Energy Estimates in the System Advisor Model

This work will detail a proposed uncertainty quantification method for PV annual energy estimates. Motivations behind the proposed methodology will be discussed, including the frequent underperformance of installed PV systems to calculated probability of exceedance values being seen in industry. A methodology in which uncertainty factors correlated to different parts of the PV annual energy model chain are applied to annual energy estimates in conjunction with uncertainty from inter-annual solar resource variability will be described. The application of this methodology in the System Advisor Model (SAM) will also be discussed with example cases and graphical outputs of the probability of exceedance data calculated for the PV annual energy estimates.

energy modeling↗

Quantifying Uncertainty in PV Energy Estimates Final Report

Uncertainty in PV energy estimates is "one of the most critical areas of lack of understanding" according to independent engineers, financiers, PV model developers, and other industry stakeholders. The primary problem is a lack of rigorous, transparent, widely accepted methods for quantifying uncertainty in energy production estimates. Uncertainty in energy production estimates arises from variability of the solar resource, inexact PV performance models and their parameters, and system reliability considerations. Uncertainty in annual energy production is frequently calculated for larger projects in order to quantify financial risk. Key statistics for energy, such as the P-values "P50" and "P90" (the annual energy values that are exceeded in future years with 50\% and 90\% probability, respectively) are used by financing institutions to calculate the repayment risk for the project. The current methods to estimate these statistics are typically proprietary, specialized, and involve significant post-processing of commercial performance model results. This black-box approach leads to inconsistent P-value estimates from different parties, which reduces investors' confidence in the results. Since the financial community bases its risk assessment on these estimates, reduced confidence increases perceived project risk, and consequently financing costs. The goal of this project was to establish a set of best practices for quantifying uncertainty in energy production estimates, including identifying what sources of uncertainty must be considered with clear definitions and metrics, determining which sources are the biggest drivers of uncertainty, and providing a computationally efficient framework for combining different sources of uncertainty that is flexible enough to accommodate substitutions of data or methods when better information is available. We engaged a wide set of stakeholders to ensure industry endorsement and adoption, and leveraged complementary projects investigating individual sources of uncertainty in great detail, as well as others' work that started down this path.

ENERGY PLANNING, POLICY, AND ECONOMY,SOLAR ENERGY↗

Estimating basis functions in massive fields under the spatial mixed effects model

Abstract Spatial prediction is commonly achieved under the assumption of a Gaussian random field by obtaining maximum likelihood estimates of parameters, and then using the kriging equations to arrive at predicted values. For massive datasets, fixed rank kriging using the expectation–maximization algorithm for estimation has been proposed as an alternative to the usual but computationally prohibitive kriging method. The method reduces computation cost of estimation by redefining the spatial process as a linear combination of basis functions and spatial random effects. A disadvantage of this method is that it imposes constraints on the relationship between the observed locations and the knots. We develop an alternative method that utilizes the spatial mixed effects model, but allows for additional flexibility by estimating the range of the spatial dependence between the observations and the knots via an alternating expectation conditional maximization algorithm. Experiments show that our methodology improves estimation without sacrificing prediction accuracy while also minimizing the additional computational burden of extra parameter estimation. The methodology is applied to a temperature dataset archived by the United States National Climate Data Center, with improved results over previous methodology.

Pazdernik, Karl↗

Assessment of DFT functionals for a minimal nitrogenase [Fe(SH) 4 H]- model employing state-of-the-art ab initio methods

We have designed a [Fe(SH)4H]- model with the fifth proton binding either to Fe or S. We show that the energy difference between these two isomers (ΔE) is hard to estimate with quantum-mechanical (QM) methods. For example, different density functional theory (DFT) methods give ΔE estimates that vary by almost 140 kJ/mol, mainly depending on the amount of exact Hartree–Fock included (0%–54%). The model is so small that it can be treated by many high-level QM methods, including coupled-cluster (CC) and multiconfigurational perturbation theory approaches. With extrapolated CC series (up to fully connected coupled-cluster calculations with singles, doubles, and triples) and semistochastic heat-bath configuration interaction methods, we obtain results that seem to be converged to full configuration interaction results within 5 kJ/mol. Our best result for ΔE is 101 kJ/mol. With this reference, we show that M06 and B3LYP-D3 give the best results among 35 DFT methods tested for this system. Brueckner doubles coupled cluster with perturbaitve triples seems to be the most accurate coupled-cluster approach with approximate triples. CCSD(T) with Kohn–Sham orbitals gives results within 4–11 kJ/mol of the extrapolated CC results, depending on the DFT method. Single-reference CC calculations seem to be reasonably accurate (giving an error of ~5 kJ/mol compared to multireference methods), even if the D1 diagnostic is quite high (0.25) for one of the two isomers.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

When less is more: How increasing the complexity of machine learning strategies for geothermal energy assessments may not lead toward better estimates

Previous moderate- and high-temperature geothermal resource assessments of the western United States utilized data-driven methods and expert decisions to estimate resource favorability. Although expert decisions can add confidence to the modeling process by ensuring reasonable models are employed, expert decisions also introduce human and, thereby, model bias. This bias can present a source of error that reduces the predictive performance of the models and confidence in the resulting resource estimates. Our study aims to develop robust data-driven methods with the goals of reducing bias and improving predictive ability. We present and compare nine favorability maps for geothermal resources in the western United States using data from the U.S. Geological Survey's 2008 geothermal resource assessment. Two favorability maps are created using the expert decision-dependent methods from the 2008 assessment (i.e., weight-of-evidence and logistic regression). With the same data, we then create six different favorability maps using logistic regression (without underlying expert decisions), XGBoost, and support-vector machines paired with two training strategies. The training strategies are customized to address the inherent challenges of applying machine learning to the geothermal training data, which have no negative examples and severe class imbalance. We also create another favorability map using an artificial neural network. We demonstrate that modern machine learning approaches can improve upon systems built with expert decisions. We also find that XGBoost, a non-linear algorithm, produces greater agreement with the 2008 results than linear logistic regression without expert decisions, because the expert decisions in the 2008 assessment rendered the otherwise linear approaches non-linear despite the fact that the 2008 assessment used only linear methods. The F1 scores for all approaches appear low (F1 score < 0.10), do not improve with increasing model complexity, and, therefore, indicate the fundamental limitations of the input features (i.e., training data). Until improved feature data are incorporated into the assessment process, simple non-linear algorithms (e.g., XGBoost) perform equally well or better than more complex methods (e.g., artificial neural networks) and remain easier to interpret.

15 GEOTHERMAL ENERGY↗

Employing Eye Trackers to Reduce Nuisance Alarms

When process operators anticipate an alarm prior to its annunciation, that alarm loses information value and becomes a nuisance. This study investigated using eye trackers to measure and adjust the salience of alarms with three methods of gaze-based acknowledgement (GBA) of alarms that estimate operator anticipation. When these methods detected possible alarm anticipation, the alarm’s audio and visual salience was reduced. A total of 24 engineering students (male = 14, female = 10) aged between 18 and 45 were recruited to predict alarms and control a process parameter in three scenario types (parameter near threshold, trending, or fluctuating). The study evaluated whether behaviors of the monitored parameter affected how frequently the three GBA methods were utilized and whether reducing alarm salience improved control task performance. The results did not show significant task improvement with any GBA methods (F(3,69) = 1.357, p = 0.263, partial η 2 = 0.056). However, the scenario type affected which GBA method was more utilized (X 2 (2, N = 432) = 30.147, p < 0.001). Alarm prediction hits with gaze-based acknowledgements coincided more frequently than alarm prediction hits without gaze-based acknowledgements (X 2 (1, N = 432) = 23.802, p < 0.001, OR = 3.877, 95% CI 2.25–6.68, p < 0.05). Participant ratings indicated an overall preference for the three GBA methods over a standard alarm design (F(3,63) = 3.745, p = 0.015, partial η 2 = 0.151). This study provides empirical evidence for the potential of eye tracking in alarm management but highlights the need for additional research to increase validity for inferring alarm anticipation.

99 - GENERAL AND MISCELLANEOUS↗

Robust Online Sequential RVFLNs for Data Modeling of Dynamic Time-Varying Systems with Application of an Ironmaking Blast Furnace

In a world where the increasing complexity of modern industrial processes brings difficulties for accurate mathematical modeling, taking advantage of data has become an efficient solution to complex dynamic process modeling issue. In this paper, we develop a novel robust online sequential version of random vector functional-link networks (RVFLNs) for data-driven modeling of dynamic time-varying system and applied it in a blast furnace (BF) ironmaking process. First, to overcome the time-varying dynamics of process and to enable the RVFLNs to learn online with avoiding data saturation, an improved online sequential version of RVFLNs (OS-RFVLNs) is first presented by online sequential learning with forgetting factor. This improved OS-RVFLNs algorithm is not only suitable for the real-time and large data transfer situation, but also can adjust the sensitivity of the algorithm to different samples with the help of the introduced forgetting factor. Second, since the output weights of the improved OS-RVFLNs as well as other RVFLNs algorithms are obtained by the least squares approach, a robustness problem may occur when the training dataset is contaminated with various outliers. To solve this problem, a Cauchy distribution weighted M-estimator is introduced to improve the robustness of the improved OS- RVFLNs. For this proposed robust OS-RVFLNs (R-OS- RVFLNs), since the weights of different outlier data are properly determined by the Cauchy distribution function, their corresponding contribution on modeling can be properly distinguished. Thus robust and better modeling results can be achieved. Experiments using actual industrial data of BF ironmaking process and comparative studies have demonstrated that the proposed method produces a better estimation accuracy and stronger robustness than other methods.

Blast furnace (BF), Modelling, Dynamic systems↗

A Fully Computable A Posteriori Error Estimate for the Stokes Equations on Polytopal Meshes

Here in this paper, we present a simple a posteriori error estimate for the weak Galerkin finite element method for the Stokes equation. This residual type estimator can be applied to general meshes such as polytopal mesh or meshes with hanging nodes. The reliability and efficiency of the estimator are proved in this paper. Five numerical tests demonstrate the effectiveness and flexibility of the adaptive mesh refinement guided by the designed error estimator.

finite element methods↗

Effects of surface diffusion in electrocatalytic CO 2 reduction on Cu revealed by kinetic Monte Carlo simulations

We report Kinetic Monte Carlo (KMC) methods are frequently used for mechanistic studies of thermally driven heterogeneous catalysis systems but are underused for electrocatalysis. Here, we develop a lattice KMC approach for electrocatalytic CO 2 reduction. The work is motivated by a prior experimental report that performed electroreduction of a mixed feed of 12 CO 2 and 13 CO on Cu; differences in the 13 C content of C2 products ethylene and ethanol (Δ 13 C) were interpreted as evidence of site selectivity. The lattice KMC model considers the effect of surface diffusion on this system. In the limit of infinitely fast diffusion (mean-field approximation), the key intermediates 12 CO* and 13 CO* would be well mixed on the surface and no evidence of site selectivity could have been observed. Using a simple two-site model and adapting a previously reported microkinetic model, we assess the effects of diffusion on the relative isotope fractions in the products using the estimated surface diffusion rate of CO* from literature reports. We find that the size of the active sites and the total surface adsorbate coverage can have a large influence on the values of Δ 13 C that can be observed. Δ 13 C is less sensitive to the CO* diffusion rate as long as it is within the estimated range. We further offer possible methods to estimate surface distribution of intermediates and to predict intrinsic selectivity of active sites based on experimental observations. This work illustrates the importance of considering surface diffusion in the study of electrochemical CO 2 reduction to multi-carbon products. Our approach is entirely based on a freely available open-source code, so will be readily adaptable to other electrocatalytic systems.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

Electromechanical Mode Estimation in the Presence of Periodic Forced Oscillations

Electromechanical modes are inherent to any inter-connected power systems which provide a measure of the small-signal stability margin of the system. A number of algorithms have been developed for the estimation of these modes using synchrophasor measurements. However, most of these algorithms are not designed to operate in the presence of forced oscillations (FO). These FOs are results of periodic rogue input driving the system. When FOs are present, estimates of system modes can be biased depending on the frequency and the amplitude of the FOs. To tackle this problem, a new algorithm is proposed in this paper to estimate system modes in the presence of FOs. In the proposed method, the over-determined modified Yule-Walker method which is used to estimate AR coefficients of an ARMA (Auto-Regressive Moving-Average) signal model is extended to an ARMA with exogenous input (ARMAX) model that incorporates the presence of FOs. Two versions of the proposed method are included in this paper based on the requirement of the information of the duration of FOs in the signal. Furthermore, the results obtained by implementing the proposed algorithm on simulated data and real-world data validate the effectiveness of both versions of the proposed method.

42 ENGINEERING↗

RADIAL FLOW PATTERN OF A SLOW CORONAL MASS EJECTION

Height–time plots of the leading edge of coronal mass ejections (CMEs) have often been used to study CME kinematics. We propose a new method to analyze the CME kinematics in more detail by determining the radial mass transport process throughout the entire CME. Thus, our method is able to estimate not only the speed of the CME front but also the radial flow speed inside the CME. We have applied this method to a slow CME with an average leading edge speed of about 480 km s{sup −1}. In the Lagrangian frame, the speeds of the individual CME mass elements stay almost constant within 2 and 15 R{sub S}, the range over which we analyzed the CME. Hence, we have no evidence of net radial forces acting on parts of the CME in this range or of a pile up of mass ahead of the CME. We find evidence that the leading edge trajectory obtained by tie-pointing may gradually lag behind the Lagrangian front-side trajectories derived from our analysis. Our results also allow a much more precise estimate of the CME energy. Compared with conventional estimates using the CME total mass and leading edge motion, we find that the latter may overestimate the kinetic energy and the gravitational potential energy.

79 ASTRONOMY AND ASTROPHYSICS↗

In Search of a Cost-Effective Approach for Estimating Dam Passage Survival

Abstract Passage of downstream-migrating fishes through hydropower dams poses a management issue in rivers throughout the world. As such, regulations have been enacted at many locales to understand and limit fish losses by requiring dam passage survival estimation. However, limiting estimates of survival to the dam and immediate tailrace presents a challenge to researchers. In this study, two alternative release–recapture methods were used to estimate dam passage survival of juvenile salmonids through Lower Granite Dam on the Snake River, Washington, in 2018. One approach, the virtual/paired-release (ViPRe) model, was extensively used in the Columbia–Snake River basin at federally operated hydroelectric dams during 2010–2014. This existing approach uses three releases of tagged smolts to isolate dam passage survival, defined as survival from the upstream dam face to the tailrace mixing zone 1–2 km downstream of the dam. An alternative approach, the virtual release/dead-fish correction (ViRDCt) model, uses one live-release group paired with a release of dead tagged fish at the dam to estimate the same survival parameter. The alternative estimation approaches were tested on two spring out-migrating stocks, yearling Chinook Salmon Oncorhynchus tshawytscha and steelhead Oncorhynchus mykiss, and on a summer out-migrating stock of subyearling Chinook Salmon. The alternative estimates for these stocks were all within 1 SE. However, the ViRDCt model produced survival estimates with SEs that were 59% smaller than those of the existing ViPRe model and did so using 42% fewer tagged fish and one less downstream acoustic detection array. Because of the reduced sample size, nearly US$540,000 less in tag costs and 4,000 fewer labor-hours were required to implement the ViRDCt model compared to the ViPRe model. As such, the ViRDCt model represents a cost-effective and precise approach to characterizing dam passage survival.

59 BASIC BIOLOGICAL SCIENCES↗

Calculating Probability Densities with Homotopy and Applications to Particle Filters

In this work, we explore a homotopy sampling procedure and its generalization, loosely based on importance sampling, known as annealed importance sampling. The procedure makes use of a known probability distribution to find, via homotopy, the unknown normalization of a target distribution, as well as samples of the target distribution. We propose a reformulation of the method that leads to a rejection sampling alternative. Estimates of the error as a function of homotopy stages and sample averages are derived for the algorithmic version of the method. These estimates are useful in making computational efficiency decisions on how the calculation should proceed, given a computer architecture. Consideration is given to how the procedure can be adapted to Bayesian stationary and nonstationary estimation problems. The connection between homotopy sampling and thermodynamic integration is made. Emphasis is placed on the non-stationary problems, and in particular, on a sequential estimation technique known as particle filtering. It is shown that a modification of the particle filter framework to include the homotopy process can improve the computational robustness of particle filters.

97 MATHEMATICS AND COMPUTING↗

Uncertainty analysis for infrasound waveform inversion: Application to explosion yield estimation

While the acoustic waveform inversion method is increasingly used in geophysical acoustics to constrain source parameters, the inversion results are often provided without any uncertainty analysis. Here, this study presents a probabilistic representation for acoustic waveform inversion and method to evaluate the inversion uncertainty using ground-truth data. A posteriori probability distribution of source estimate is described by a priori waveform misfit covariance and the variance of acoustic source model. The probabilistic framework is applied to local explosion infrasound to estimate the yields of explosions and uncertainty. Estimated yields showed overall good agreement with the true yields (less than 25% errors). The uncertainty of the estimated yield is represented by the sum of the waveform inversion uncertainty and source model uncertainty. It is shown that the yield uncertainty attributed to local infrasound inversion (within 10 km) is as small as the uncertainty caused by 10% prediction errors in the acoustic source model. Finally, these results indicate that the acoustic source model uncertainty should also be considered for accurate yield estimation and that local infrasound can be a valuable tool to understand the magnitude of the source uncertainty.

58 GEOSCIENCES↗

CO2-SCREEN

CO2 Storage prospeCtive Resource Estimation Excel aNalysis (CO2-SCREEN) is a tool which was developed by the Department of Energy’s National Energy Technology Laboratory (DOE-NETL) to screen saline formations for prospective CO2 storage resources. CO2-SCREEN applies U.S. DOE methods and equations for estimating prospective CO2 storage resources for saline formations. This provides a dependable method for calculating prospective CO2 storage resources which allows for consistent comparison of results between different research efforts. CO2-SCREEN latest version (v5.0) is a python application that utilizes a Java front end user interface. The 3.0 version consists of an Excel spreadsheet containing geologic inputs and outputs, linked to a GoldSim Player model that calculates prospective CO2 storage resources via Monte Carlo simulation. https://www.youtube.com/watch?v=lhakk-HYfOI

CO2 Storage Resource,Prospective Storage,Screening↗

SU-G-IeP3-04: Effective Dose Measurements in Fast Kvp Switch Dual Energy Computed Tomography

Purpose: The objective of this study was two-fold: (a) to test a new approach to approximating organ dose by using the effective energy of the combined 80kV/140kV beam in dual-energy (DE) computed tomography (CT), and (b) to derive the effective dose (ED) in the abdomen-pelvis protocol in DECT. Methods: A commercial dual energy CT scanner was employed using a fast-kV switch abdomen/pelvis protocol alternating between 80 kV and 140 kV. MOSFET detectors were used for organ dose measurements. First, an experimental validation of the dose equivalency between MOSFET and ion chamber (as a gold standard) was performed using a CTDI phantom. Second, the ED of DECT scans was measured using MOSFET detectors and an anthropomorphic phantom. For ED calculations, an abdomen/pelvis scan was used using ICRP 103 tissue weighting factors; ED was also computed using the AAPM Dose Length Product (DLP) method and compared to the MOSFET value. Results: The effective energy was determined as 42.9 kV under the combined beam from half-value layer (HVL) measurement. ED for the dual-energy scan was calculated as 16.49 ± 0.04 mSv by the MOSFET method and 14.62 mSv by the DLP method. Conclusion: Tissue dose in the center of the CTDI body phantom was 1.71 ± 0.01 cGy (ion chamber) and 1.71 ± 0.06 (MOSFET) respectively; this validated the use of effective energy method for organ dose estimation. ED from the abdomen-pelvis scan was calculated as 16.49 ± 0.04 mSv by MOSFET and 14.62 mSv by the DLP method; this suggests that the DLP method provides a reasonable approximation to the ED.

60 APPLIED LIFE SCIENCES↗

Comparing Capacity Credit Calculations for Wind: A Case Study in Texas

The degree to which wind energy can contribute to the capacity needed to meet resource adequacy requirements, also known as capacity credit (CC), varies regionally with wind resource and correlation to net load. CC is an important metric widely used for resource planning and resource adequacy assessments. However, there are multiple methods for computing and estimating CC, depending on specific needs, access to data, and computational burden. It is unclear the extent to which the CC computation method may influence the result. To address this, we use a probabilistic resource adequacy tool and multiple approximation methods to systematically assess the CC of wind for near-term wind deployment under a case study in Texas. We find that proper consideration of transmission constraints is important; some approximation methods may overestimate the CC of wind due to a lack of consideration of transmission constraints, while other approximation methods may underestimate the CC by not capturing the ability of wind to be shipped to neighboring regions. In this case study, we find that several approximation methods do come close to the CC calculated by more robust probabilistic methods. However, the best approximation method may vary on a case-by-case basis, depending on system-specific considerations.

17 WIND ENERGY↗