Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Parameter estimation using data”

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 55 records · Page 3
All sourcesNASA NTRSDOE OSTI
DOE OSTI2015

A Hubble diagram for quasars

We present a new method to test the ΛCDM cosmological model and to estimate cosmological parameters based on the nonlinear relation between the ultraviolet and X-ray luminosities of quasars. We built a data set of 1138 quasars by merging several samples from the literature with X-ray measurements at 2 keV and SDSS photometry, which was used to estimate the extinction-corrected 2500 Å flux. We obtained three main results: (1) we checked the nonlinear relation between X-ray and UV luminosities in small redshift bins up to z∼6, confirming that the relation holds at all redshifts with the same slope; (2) we built a Hubble diagram for quasars up to z∼6, which is well matched to that of supernovae in the common z = 0–1.4 redshift interval and extends the test of the cosmological model up to z∼6; and (3) we showed that this nonlinear relation is a powerful tool for estimating cosmological parameters. Using the present data and assuming a ΛCDM model, we obtain Ω{sub M} = 0.22{sub −0.08}{sup +0.10} and Ω{sub Λ} = 0.92{sub −0.30}{sup +0.18} (Ω{sub M} = 0.28 ± 0.04 and Ω{sub Λ} = 0.73 ± 0.08 from a joint quasar-SNe fit). Much more precise measurements will be achieved with future surveys. A few thousand SDSS quasars already have serendipitous X-ray observations from Chandra or XMM-Newton, and at least 100,000 quasars with UV and X-ray data will be made available by the extended ROentgen Survey with an Imaging Telescope Array all-sky survey in a few years. The Euclid, Large Synoptic Survey Telescope, and Advanced Telescope for High ENergy Astrophysics surveys will further increase the sample size to at least several hundred thousand. Our simulations show that these samples will provide tight constraints on the cosmological parameters and will allow us to test for possible deviations from the standard model with higher precision than is possible today.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2016

Combining On-Line Characterization Tools with Modern Software Environments for Optimal Operation of Polymerization Processes

Here, this paper discusses the initial steps towards the formulation and implementation of a generic and flexible model centric framework for integrated simulation, estimation, optimization and feedback control of polymerization processes. For the first time it combines the powerful capabilities of the automatic continuous on-line monitoring of polymerization system (ACOMP), with a modern simulation, estimation and optimization software environment towards an integrated scheme for the optimal operation of polymeric processes. An initial validation of the framework was performed for modelling and optimization using literature data, illustrating the flexibility of the method to apply under different systems and conditions. Subsequently, off-line capabilities of the system were fully tested experimentally for model validations, parameter estimation and process optimization using ACOMP data. Experimental results are provided for free radical solution polymerization of methyl methacrylate.

42 ENGINEERING↗
DOE OSTI2017

NONPARAMETRIC EMPIRICAL BAYES ESTIMATION THROUGH DECONVOLUTION IN PROBABILISTIC RISK ANALYSIS

The advent of Markov Chain Monte Carlo (MCMC) simulations and availability of open source software made full hierarchical Bayesian analysis tractable and a popular choice in parameter estimation for Probabilistic Risk Assessment (PRA). However, despite its theoretical attractiveness, the hierarchical Bayesian analysis is not without its problems. Two of the most prominent practical difficulties in applying hierarchical Bayes analysis in practice to analyze source-to-source variability, for example, is its sensitivity to the selection of the first-stage prior and dependence of the rate of convergence on the selection of the first-stage prior. The first-stage prior in hierarchical Bayesian analysis is usually assumed having a parametric form, often conjugate to the likelihood function (aleatory model). To facilitate convergence for hierarchical Bayes models, the first-stage prior is frequently set through empirical Bayes estimate, which can use available historical data to find the parameters of the first-stage prior. The paper discusses a new nonparametric empirical Bayes estimation and compare it to several well-known parametric estimates such as method of moments and maximum likelihood. The new nonparametric technique exploits the prior predictive distribution as integral equation and proceeds to solve it with respect to prior assuming the sampling data distribution and aleatory model are available. The discussion covers topics such as selection of aleatory model for future data, selection of regularization parameter, and density estimation for historical data. For comparison, the paper is using contrived as well as real-world data from utilities. Utilities data represent failure data with Poisson distribution used for aleatory model.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2017

Data Synthesis and Data Assimilation at Global Change Experiments and Fluxnet Toward Improving Land Process Models

The project was conducted during the period from 7/1/2012 to 6/30/2017 with three major tasks: (1) data synthesis and development of data assimilation (DA) techniques to constrain modeled ecosystem feedback to climate change; (2) applications of DA techniques to improve process models at different scales from ecosystem to regions and the globe; and 3) improvements of modeling soil carbon (C) dynamics by land surface models. During this period, we have synthesized published data from soil incubation experiments (e.g., Chen et al., 2016; Xu et al., 2016; Feng et al., 2016), global change experiments (e.g., Li et al., 2013; Shi et al., 2015, 2016; Liang et al., 2016) and fluxnet (e.g., Niu et al., 2012., Xia et al., 2015; Li et al., 2016). These data have been organized into multiple data products and have been used to identify general mechanisms and estimate parameters for model improvement. We used the data sets that we collected and the DA techniques to improve model performance of both ecosystem models and global land models. The objectives are: 1) to improve model simulations of litter and soil carbon storage (e.g., Schädel et al., 2013; Hararuk and Luo, 2014; Hararuk et al., 2014; Liang et al., 2015); 2) to explore the effects of CO 2 , warming and precipitation on ecosystem processes (e.g., van Groenigen et al., 2014; Shi et al., 2015, 2016; Feng et al., 2017); and 3) to estimate parameters variability in different ecosystems (e.g., Li et al., 2016). We developed a traceability framework, which was based on matrix approaches and decomposed the modeled steady-state terrestrial ecosystem carbon storage capacity into four can trace the difference in ecosystem carbon storage capacity among different biomes to four traceable components: net primary productivity (NPP), baseline C residence times, environmental scalars and climate forcing (Xia et al., 2013). With this framework, we can diagnose the differences in modeled carbon storage across ecosystems, biomes, and models. This framework has been successfully implemented by several global land models, such as CABLE (Xia et al., 2013), LPJ-GUESS (Ahlström et al., 2015), CLM (Hararuk et al., 2014; Huang et al., 2017, submitted; Shi et al., 2017, submitted), and ORCHIDEE (Huang et al., 2017, unpublished). Moreover, we have identified the theoretical foundation of the determinants of transient C storage dynamics by adding another term, C storage potential, to the steady-state traceability framework (Luo et al., 2017). The theoretical foundation of transient C storage dynamics has been applied to develop a transient traceability framework to explore the traceable components of transient C storage dynamics responded to the rising CO 2 and climate change in the two contrasting ecosystem types Duke needleleaved forest and Harvard deciduous broadleaved forest (Jiang et al., 2017, in revision). Overall, with the data synthesis, data assimilation techniques, and the steady-state and transient traceability frameworks, we have greatly improved land process models for predicting responses and feedback of terrestrial C dynamics to global change. The matrix approaches has the potential to be applied in theoretical research on nitrogen and phosphorus cycle, and therefore, the coupling of carbon-nitrogen-phosphorus.

58 GEOSCIENCES↗
DOE OSTI2016

Effects of monohydric alcohols and polyols on the thermal stability of a protein

The thermal stability of a protein is lowered by the addition of a monohydric alcohol, and this effect becomes larger as the size of hydrophobic group in an alcohol molecule increases. By contrast, it is enhanced by the addition of a polyol possessing two or more hydroxyl groups per molecule, and this effect becomes larger as the number of hydroxyl groups increases. Here, we show that all of these experimental observations can be reproduced even in a quantitative sense by rigid-body models focused on the entropic effect originating from the translational displacement of solvent molecules. The solvent is either pure water or water-cosolvent solution. Three monohydric alcohols and five polyols are considered as cosolvents. In the rigid-body models, a protein is a fused hard spheres accounting for the polyatomic structure in the atomic detail, and the solvent is formed by hard spheres or a binary mixture of hard spheres with different diameters. The effective diameter of cosolvent molecules and the packing fractions of water and cosolvent, which are crucially important parameters, are carefully estimated using the experimental data of properties such as the density of solid crystal of cosolvent, parameters in the pertinent cosolvent-cosolvent interaction potential, and density of water-cosolvent solution. We employ the morphometric approach combined with the integral equation theory, which is best suited to the physical interpretation of the calculation result. It is argued that the degree of solvent crowding in the bulk is the key factor. When it is made more serious by the cosolvent addition, the solvent-entropy gain upon protein folding is magnified, leading to the enhanced thermal stability. When it is made less serious, the opposite is true. The mechanism of the effects of monohydric alcohols and polyols is physically the same as that of sugars. However, when the rigid-body models are employed for the effect of urea, its addition is predicted to enhance the thermal stability, which conflicts with the experimental fact. We then propose, as two essential factors, not only the solvent-entropy gain but also the loss of protein-solvent interaction energy upon protein folding. The competition of changes in these two factors induced by the cosolvent addition determines the thermal-stability change.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2017

Leveraging 35 years of Pinus taeda research in the southeastern US to constrain forest carbon cycle predictions: regional data assimilation using ecosystem experiments

Predicting how forest carbon cycling will change in response to climate change and management depends on the collective knowledge from measurements across environmental gradients, ecosystem manipulations of global change factors, and mathematical models. Formally integrating these sources of knowledge through data assimilation, or model–data fusion, allows the use of past observations to constrain model parameters and estimate prediction uncertainty. Data assimilation (DA) focused on the regional scale has the opportunity to integrate data from both environmental gradients and experimental studies to constrain model parameters. Here, we introduce a hierarchical Bayesian DA approach (Data Assimilation to Predict Productivity for Ecosystems and Regions, DAPPER) that uses observations of carbon stocks, carbon fluxes, water fluxes, and vegetation dynamics from loblolly pine plantation ecosystems across the southeastern US to constrain parameters in a modified version of the Physiological Principles Predicting Growth (3-PG) forest growth model. The observations included major experiments that manipulated atmospheric carbon dioxide (CO 2 ) concentration, water, and nutrients, along with nonexperimental surveys that spanned environmental gradients across an 8.6 × 10 5 km 2 region. We optimized regionally representative posterior distributions for model parameters, which dependably predicted data from plots withheld from the data assimilation. While the mean bias in predictions of nutrient fertilization experiments, irrigation experiments, and CO 2 enrichment experiments was low, future work needs to focus modifications to model structures that decrease the bias in predictions of drought experiments. Predictions of how growth responded to elevated CO 2 strongly depended on whether ecosystem experiments were assimilated and whether the assimilated field plots in the CO 2 study were allowed to have different mortality parameters than the other field plots in the region. We present predictions of stem biomass productivity under elevated CO 2 , decreased precipitation, and increased nutrient availability that include estimates of uncertainty for the southeastern US. Overall, we (1) demonstrated how three decades of research in southeastern US planted pine forests can be used to develop DA techniques that use multiple locations, multiple data streams, and multiple ecosystem experiment types to optimize parameters and (2) developed a tool for the development of future predictions of forest productivity for natural resource managers that leverage a rich dataset of integrated ecosystem observations across a region.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Using ASCEM Modeling and Visualization to Inform Stakeholders of Contaminant Plume Evolution and Remediation Efficacy at F-Basin Savannah River, SC – 15156

Communication with stakeholders, regulatory agencies, and the public is an essential part of implementing different remediation and monitoring activities, and developing site closure strategies at contaminated sites. Modeling of contaminant plume evolution plays a critical role in estimating the benefit, cost, and risk of particular options. At the same time, effective visualization of monitoring data and modeling results are particularly important for conveying the significance of the results and observations. In this paper, we present the results of the Advanced Simulation Capability for Environmental Management (ASCEM) project, including the discussion of the capabilities of newly developed ASCEM software package, along with its application to the F-Area Seepage Basins located in the U.S. Department of Energy Savannah River Site (SRS). ASCEM software includes state-of-the-art numerical methods for simulating complex flow and reactive transport, as well as various toolsets such as a graphical user interface (GUI), visualization, data management, uncertainty quantification, and parameter estimation. Using this software, we have developed an advanced visualization of tritium plume migration coupled with a data management system, and simulated a three-dimensional model of flow and plume evolution on a high-performance computing platform. We evaluated the effect of engineered flow barriers on a nonreactive tritium plume, through advanced plume visualization and modeling of tritium plume migration. In addition, we developed a geochemical reaction network to describe complex geochemical processes at the site, and evaluated the impact of coupled hydrological and geochemical heterogeneity. These results are expected to support SRS’s monitoring activities and operational decisions.

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2015

Using ASCEM modeling and visualization to inform stakeholders of contaminant plume evolution and remediation efficacy at F-Basin Savannah river, SC - 15156

Communication with stakeholders, regulatory agencies, and the public is an essential part of implementing different remediation and monitoring activities, and developing site closure strategies at contaminated sites. Modeling of contaminant plume evolution plays a critical role in estimating the benefit, cost, and risk of particular options. At the same time, effective visualization of monitoring data and modeling results are particularly important for conveying the significance of the results and observations. In this paper, we present the results of the Advanced Simulation Capability for Environmental Management (ASCEM) project, including the discussion of the capabilities of newly developed ASCEM software package, along the its application to the F-Area Seepage Basins located in the U.S. Department of Energy Savannah River Site (SRS). ASCEM software includes state-of-the-art numerical methods for simulating complex flow and reactive transport, as well as various tool-sets such as a graphical user interface (GUI), visualization, data management, uncertainty quantification, and parameter estimation. Using this software, we have developed an advanced visualization of tritium plume migration coupled with a data management system, and simulated a three-dimensional model of flow and plume evolution on a high-performance computing platform. We evaluated the effect of engineered flow barriers on a nonreactive tritium plume, through advanced plume visualization and modeling of tritium plume migration. In addition, we developed a geochemical reaction network to describe complex geochemical processes at the site, and evaluated the impact of coupled hydrological and geochemical heterogeneity. These results are expected to support SRS's monitoring activities and operational decisions. (authors)

12 MANAGEMENT OF RADIOACTIVE AND NON-RADIOACTIVE W↗
DOE OSTI2014

Risk assessment of turbine rotor failure using probabilistic ultrasonic non-destructive evaluations

The study presents a method and application of risk assessment methodology for turbine rotor fatigue failure using probabilistic ultrasonic nondestructive evaluations. A rigorous probabilistic modeling for ultrasonic flaw sizing is developed by incorporating the model-assisted probability of detection, and the probability density function (PDF) of the actual flaw size is derived. Two general scenarios, namely the ultrasonic inspection with an identified flaw indication and the ultrasonic inspection without flaw indication, are considered in the derivation. To perform estimations for fatigue reliability and remaining useful life, uncertainties from ultrasonic flaw sizing and fatigue model parameters are systematically included and quantified. The model parameter PDF is estimated using Bayesian parameter estimation and actual fatigue testing data. The overall method is demonstrated using a realistic application of steam turbine rotor, and the risk analysis under given safety criteria is provided to support maintenance planning.

42 ENGINEERING↗
DOE OSTI2021

Does Bayesian model averaging improve polynomial extrapolations? Two toy problems as tests

We assess the accuracy of Bayesian polynomial extrapolations from small parameter values, x, to large values of x. We consider a set of polynomials of fixed order, intended as a proxy for a fixed-order effective field theory (EFT) description of data. We employ Bayesian model averaging (BMA) to combine results from different order polynomials (EFT orders). Our study considers two 'toy problems' where the underlying function used to generate data sets is known. We use Bayesian parameter estimation to extract the polynomial coefficients that describe these data at low x. A 'naturalness' prior is imposed on the coefficients, so that they are $\mathcal{O}(1)$. We BMA different polynomial degrees by weighting each according to its Bayesian evidence and compare the predictive performance of this BMA with that of the individual polynomials. In conclusion, the credibility intervals on the BMA forecast have the stated coverage properties more consistently than does the highest evidence polynomial, though BMA does not necessarily outperform every polynomial.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
DOE OSTI2014

Ensemble-Based Parameter Estimation in a Coupled General Circulation Model

Parameter estimation provides a potentially powerful approach to reduce model bias for complex climate models. Here, in a twin experiment framework, the authors perform the first parameter estimation in a fully coupled ocean–atmosphere general circulation model using an ensemble coupled data assimilation system facilitated with parameter estimation. The authors first perform single-parameter estimation and then multiple-parameter estimation. In the case of the single-parameter estimation, the error of the parameter [solar penetration depth (SPD)] is reduced by over 90% after ~40 years of assimilation of the conventional observations of monthly sea surface temperature (SST) and salinity (SSS). The results of multiple-parameter estimation are less reliable than those of single-parameter estimation when only the monthly SST and SSS are assimilated. Assimilating additional observations of atmospheric data of temperature and wind improves the reliability of multiple-parameter estimation. The errors of the parameters are reduced by 90% in ~8 years of assimilation. Finally, the improved parameters also improve the model climatology. With the optimized parameters, the bias of the climatology of SST is reduced by ~90%. Altogether, this study suggests the feasibility of ensemble-based parameter estimation in a fully coupled general circulation model.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2024

Increasing the power of weak lensing data with multipole-based intrinsic alignment estimators

ABSTRACT It has long been known that galaxy shapes align coherently with the large-scale density field. Characterizing this effect is essential to interpreting measurements of weak gravitational lensing, the deflection of light from distant galaxies by matter overdensities along the line of sight, as it also produces coherent galaxy alignments that we wish to interpret in terms of a cosmological model. Existing direct measurements of intrinsic alignments using galaxy samples with high-quality shape and redshift measurements typically use well-understood but sub-optimal projected estimators, which do not make good use of the information in the data when comparing those estimators to theoretical models. We demonstrate a more optimal estimator, based on a multipole expansion of the correlation functions or power spectra, for direct measurements of galaxy intrinsic alignments. We show that even using the lowest order multipole alone increases the significance of inferred model parameters using simulated and real data, without any additional modelling complexity. We apply this estimator to measurements of parameters of the non-linear alignment model using data from the Sloan Digital Sky survey, demonstrating consistent results with a factor of ∼2 greater precision in parameter fits to intrinsic alignments models. This result is functionally equivalent to quadrupling the survey area, but without the attendant costs – thereby demonstrating the value in using this new estimator in current and future intrinsic alignments measurements using spectroscopic galaxy samples.

Astronomy & Astrophysics↗
DOE OSTI2021

Santa Barbara Desalination Digital Twin Technical Report

A digital twin is a computer model that has sufficient breadth and fidelity to enable diagnostic and predictive analyses for outputs of interest. We have constructed a digital twin of the SB desal plant using the ProteusLib software developed for NAWI under task 4.2.2, and specifically the ProteusLib models of pumps, filters, mixers, splitters, and reverse osmosis units and a seawater property model. The steps needed to configure this digital twin with the real conditions of the SB desal plant are: (1) Data reconciliation, which is gathering the data and plant connectivity and ensure that all relevant values are present and in the correct units and validating the data using material and energy balances, and correct for errors, gaps, or deviations; (2) Parameter estimation, which is matching model parameters with data from previous step, filling missing parameters with initial estimates, then iteratively refining estimated values. Once we completed the digital twin, we used the optimization capabilities of ProteusLib to compute outputs of interest over a chosen range of inputs. The primary output of interest to the owners of the plant is its specific energy consumption (kWh/acre-foot). The remainder of this report will provide more detail on the data reconciliation, parameter estimation, and analyses we performed.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2016

Distribution system model calibration with big data from AMI and PV inverters

Efficient management and coordination of distributed energy resources with advanced automation schemes requires accurate distribution system modeling and monitoring. Big data from smart meters and photovoltaic (PV) micro-inverters can be leveraged to calibrate existing utility models. This paper presents computationally efficient distribution system parameter estimation algorithms to improve the accuracy of existing utility feeder radial secondary circuit model parameters. The method is demonstrated using a real utility feeder model with advanced metering infrastructure (AMI) and PV micro-inverters, along with alternative parameter estimation approaches that can be used to improve secondary circuit models when limited measurement data is available. Lastly, the parameter estimation accuracy is demonstrated for both a three-phase test circuit with typical secondary circuit topologies and single-phase secondary circuits in a real mixed-phase test system.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2017

Bayesian calibration of terrestrial ecosystem models: A study of advanced Markov chain Monte Carlo methods

Calibration of terrestrial ecosystem models is important but challenging. Bayesian inference implemented by Markov chain Monte Carlo (MCMC) sampling provides a comprehensive framework to estimate model parameters and associated uncertainties using their posterior distributions. The effectiveness and efficiency of the method strongly depend on the MCMC algorithm used. In this study, a Differential Evolution Adaptive Metropolis (DREAM) algorithm was used to estimate posterior distributions of 21 parameters for the data assimilation linked ecosystem carbon (DALEC) model using 14 years of daily net ecosystem exchange data collected at the Harvard Forest Environmental Measurement Site eddy-flux tower. The DREAM is a multi-chain method and uses differential evolution technique for chain movement, allowing it to be efficiently applied to high-dimensional problems, and can reliably estimate heavy-tailed and multimodal distributions that are difficult for single-chain schemes using a Gaussian proposal distribution. The results were evaluated against the popular Adaptive Metropolis (AM) scheme. DREAM indicated that two parameters controlling autumn phenology have multiple modes in their posterior distributions while AM only identified one mode. The calibration of DREAM resulted in a better model fit and predictive performance compared to the AM. DREAM provides means for a good exploration of the posterior distributions of model parameters. Lastly, it reduces the risk of false convergence to a local optimum and potentially improves the predictive performance of the calibrated model.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2017

Modelling bivariate extreme precipitation distribution for data-scarce regions using Gumbel-Hougaard copula with maximum entropy estimation

A new method of parameter estimation in data scarce regions is valuable for bivariate hydrological extreme frequency analysis. Here, this paper proposes a new method of parameter estimation (maximum entropy estimation, MEE) for both Gumbel and Gumbel–Hougaard copula in situations when insufficient data are available. MEE requires only the lower and upper bounds of two hydrological variables. To test our new method, two experiments to model the joint distribution of the maximum daily precipitation at two pairs of stations on the tributaries of Heihe and Jinghe River, respectively, were performed and compared with the method of moments, correlation index estimation, and maximum likelihood estimation, which require a large amount of data. Both experiments show that for the Ye Niugou and Qilian stations, the performance of MEE is nearly identical to those of the conventional methods. For the Xifeng and Huanxian stations, MEE can capture information indicating that the maximum daily precipitation at the Xifeng and Huanxian stations has an upper tail dependence, whereas the results generated by correlation index estimation and maximum likelihood estimation are unreasonable. Moreover, MEE is proved tobe generally reliable and robust by many simulations under three different situations. TheGumbel–Hougaard copula with MEE can also be applied to the bivariate frequency analysis ofother extreme events in data-scarce regions.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2020

Calibrating Synchronous-Generator-Interfaced DG Models in Microgrids Using Multiple Event Data

Microgrids, consisting of distributed generators (DGs), loads, and energy storage systems, play an important role in future smart grids. In order to evaluate how microgrids operate and their dynamic response, it is necessary to develop accurate dynamic models of DGs in microgrids. This paper proposes a systematic method to calibrate synchronous-generator-interfaced DG models. Challenges associated with model validation in microgrids are analyzed. Underlying parameters are categorized into two groups, i.e., the steady-state parameters and time constants, which are estimated in two successive stages using multiple event data. This methodology ensures that the validated model with the estimated parameters are applicable for all events. A modified unscented Kalman filter (UKF) is proposed to deal with the cases in which measurement of power angle is unavailable. The effectiveness of the proposed method is validated using field test data of a 2.6 MVA diesel generator in a real microgrid.

Wang, Zhiwen↗