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 109 records · Page 6
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

Combining translational and rotational seismic motions to invert local-scale seismic data for time-variable moment tensors: do rotational motions help for high-frequency seismic data produced by underground explosions?

SUMMARY We present an analysis of combining translational and rotational seismic data in an inversion for the time-variable source time functions corresponding to the components of the seismic moment tensor. We conduct a series of numerical experiments where the data are simulated by a combination of an underground explosion and a co-located double couple shear source and recorded on surface-mounted seismometers within 1–2 km of the source. The experiments are designed to mimic explosion seismology experiments, and thus the data are in the 1–10 Hz frequency range and contain very few surface waves. We use a Monte Carlo method to propagate Earth model uncertainty into the estimates of seismic source parameters. In our experiments, we find that the uncertainty of the estimated seismic source parameters increases when we add rotational seismic motions to the inversion when using a constant number of data channels. In this case, the increased degree of uncertainty in the final results is most likely due to the near-surface Earth model uncertainty that we introduce in our simulations. However, for a fixed number of seismic stations, adding rotational seismic motions to the inversion acts to decrease the uncertainty of the estimated seismic source parameters, most likely due to the increase in the number of data channels used in the inversion.

Poppeliers, Christian (ORCID:0000000159526849)↗
DOE OSTI2016

Model-based analysis of environmental controls over ecosystem primary production in an alpine tundra dry meadow

We investigated several key limiting factors that control alpine tundra productivity by developing an ecosystem biogeochemistry model. The model simulates the coupled cycling of carbon (C), nitrogen (N), and phosphorus (P) and their interactions with gross primary production (GPP). It was parameterized with field observations from an alpine dry meadow ecosystem using a global optimization strategy to estimate the unknown parameters. The model, along with the estimated parameters, was first validated against independent data and then used to examine the environmental controls over plant productivity. Our results show that air temperature is the strongest limiting factor to GPP in the early growing season, N availability becomes important during the middle portion of the growing season, and soil moisture is the strongest limiting factors by late in the growing season. Overall, the controls over GPP during the growing season, from strongest to weakest, are soil moisture content, air temperature, N availability, and P availability. This simulation provides testable predictions of the shifting nature of physical and nutrient limitations on plant growth. The model also indicates that changing environmental conditions in the alpine will likely lead to changes in productivity. For example, warming eliminates the control of P availability on GPP and makes N availability surpass air temperature to become the second strongest limiting factor. In contrast, an increase in atmospheric nutrient deposition eliminates the control of N availability and enhances the importance of P availability. Furthermore, these analyses provide a quantitative and conceptual framework that can be used to test predictions and refine ecological analyses at this long-term ecological research site.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

Simulation of Impedance Changes with Aging in Lithium Titanate-based Cells Using Physics-Based Dimensionless Modeling

Quantifying aging effects in lithium-ion cells with chemistries that have a flat open circuit potential is challenging. We implement a physics-based electrochemical model to track changes in the electrochemical impedance response of lithium titanate-based cells. Frequency domain equations of a pseudo two-dimensional model are made dimensionless, and the corresponding non-dimensional parameters are estimated using a Levenberg-Marquardt routine. The model weighs the relative contributions of changes in diffusion, ionic conduction within the electrolyte phase against solid phase electronic conduction towards cell aging. Solid-phase diffusion, charge transfer resistance and double layer capacitance at the solid-liquid interface are accounted for in the particle impedance. The estimation routine tracks dimensionless parameters using accelerated cycling data from full cells over 1000 cycles. The model can be deployed within a short time for state estimation using physics-based models without requiring prior knowledge of the battery chemistry, format, or capacity.

25 ENERGY STORAGE↗
DOE OSTI2023

Determining circuit model parameters from operation data for PV system degradation analysis: $\mathrm{PVPRO}$

Physics-based circuit parameters like series and shunt resistance are essential to provide insights into the degradation status of photovoltaic (PV) arrays. However, calculating these parameters typically requires a full current-voltage characteristic (I-V curve), the acquisition of which involves specific measurement devices and costly methods. Thus, I-V curves of the PV system level are often not available. Here this paper proposes a methodology (PVPRO) to estimate these I-V curve parameters using only operation (string-level DC voltage and current) and weather data (irradiance and temperature). PVPRO first performs multi-stage data pre-processing to remove noisy data. Next, the time-series DC data are used to fit an equivalent circuit single-diode model (SDM) to estimate the circuit parameters by minimizing the differences between the measured and estimated values. In this way, the time evolutions of the SDM parameters are obtained. We evaluate PVPRO on synthetic datasets and find an excellent estimation of both SDM and the key I-V parameters (e.g., open-circuit voltage, short-circuit current, maximum power, etc.) with an average relative error of 0.55%. The performance, especially the extracted degradation rate of parameters, is robust to various measurement noises and the presence of faults. In addition, PVPRO is applied to a 271 kW PV field system. The relative error between the real and estimated operation voltage and current is less than 1%, suggesting that degradation trends are well captured. PVPRO represents a promising open-source tool to extract the time-series degradation trends of key PV parameters from routine operation data.

14 SOLAR ENERGY↗
DOE OSTI2013

GRID-BASED EXPLORATION OF COSMOLOGICAL PARAMETER SPACE WITH SNAKE

We present a fully parallelized grid-based parameter estimation algorithm for investigating multidimensional likelihoods called Snake, and apply it to cosmological parameter estimation. The basic idea is to map out the likelihood grid-cell by grid-cell according to decreasing likelihood, and stop when a certain threshold has been reached. This approach improves vastly on the 'curse of dimensionality' problem plaguing standard grid-based parameter estimation simply by disregarding grid cells with negligible likelihood. The main advantages of this method compared to standard Metropolis-Hastings Markov Chain Monte Carlo methods include (1) trivial extraction of arbitrary conditional distributions; (2) direct access to Bayesian evidences; (3) better sampling of the tails of the distribution; and (4) nearly perfect parallelization scaling. The main disadvantage is, as in the case of brute-force grid-based evaluation, a dependency on the number of parameters, N{sub par}. One of the main goals of the present paper is to determine how large N{sub par} can be, while still maintaining reasonable computational efficiency; we find that N{sub par} = 12 is well within the capabilities of the method. The performance of the code is tested by comparing cosmological parameters estimated using Snake and the WMAP-7 data with those obtained using CosmoMC, the current standard code in the field. We find fully consistent results, with similar computational expenses, but shorter wall time due to the perfect parallelization scheme.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Data-Driven Load Modeling and Forecasting of Residential Appliances

The expansion of residential demand side management programs and increased deployment of controllable loads require accurate appliance-level load modeling and forecasting. This paper proposes a conditional hidden semi-Markov model to describe the probabilistic nature of residential appliance demand. Model parameters are estimated directly from power consumption data using scalable statistical learning methods. We also propose an algorithm for short-term load forecasting as a key application for appliance-level load models. Case studies performed using granular sub-metered power measurements from various types of appliances demonstrate the effectiveness of the proposed load model for short-term prediction.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2017

Method for monitoring/managing electrochemical energy device by detecting intercalation stage changes

A method for determining an operating state (e.g., state-of-charge or state-of-health) and/or generating management (charge/discharge) control information in a system including an electrochemical energy device (EED, e.g., a rechargeable Li-ion battery, supercapacitor or fuel cell) that uses optical sensors to detect the intercalation stage change events occurring in the EED. The externally or internally mounted optical sensors measure operating parameter (e.g., strain and/or temperature) changes of the EED during charge/recharge cycling, and transmit measured parameter data using light signals sent over optical fibers to a detector/converter. A processor then analyzes the measured parameter data, e.g., using a model-based estimation process, to detect intercalation stage changes (i.e., crystalline structure changes caused by migration of guest species, such as Li-ions, between the EED's anode and cathode), and generates the operating state and charge/discharge control information based the analysis.

25 ENERGY STORAGE↗
DOE OSTI2022

Parameter Calibration for Johnson Cook and Preston-Tonks-Wallace Material Strength Models with Uncertainty Quantification

In this study we perform a Bayesian calibration of the parameters in the Johnson Cook (JC) material strength model, with uncertainty, using experiments with for a range of low and medium strain rates. For the parameter calibration, we used a variational Bayesian approach with experimental data from Hopkinson bar and quasi-static tests done on Oxygen Free High Conductivity (OFHC) copper. The estimated parameter values matched well with experimental data for a modest range of strain rates and temperatures. Through this method, we also recovered uncertainty and correlation information for the estimated parameter values. We also compared our results to a calibration of parameters in the Preston-Tonks-Wallace (PTW) material strength model, using the same variational Bayesian method. The parameters estimated for both models provided good agreement with the experimental data.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Characterizing vortex tangle properties in steady-state He II counterflow using particle tracking velocimetry

Historically, there has been little faith in particle tracking velocimetry (PTV) as a tool to make quantitative measurements of thermal counterflow in He II, since tracer particle motion is complicated by influences from the normal fluid, superfluid, and quantized vortex lines, or a combination thereof. Recently, we introduced a scheme for differentiating particles trapped on vortices (G1) from particles entrained by the normal fluid (G2). In this paper, we apply this scheme to demonstrate the utility of PTV for quantitative measurements of vortex dynamics in He II counterflow.We estimate Ι, the mean vortex line spacing, using G2 velocity data, and c 2 , a parameter related to the mean curvature radius of vortices and energy dissipation in quantum turbulence, using G1 velocity data. We find that both estimations show good agreement with existing measurements that were obtained using traditional experimental methods. This is of particular consequence since these parameters likely vary in space, and PTV offers the advantage of spatial resolution. Furthermore, we also show a direct link between power-law tails in transverse particle velocity probability density functions (PDFs) and reconnection of vortex lines on which G1 particles are trapped.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2023

Estimation of hydraulic conductivity in a watershed using sparse multi-source data via Gaussian process regression and Bayesian experimental design

Enhanced water management systems depend on accurate estimation of subsurface hydraulic properties. However, geologic formations can vary significantly, so information from a single source (e.g., widely spaced boreholes) is insufficient in characterizing subsurface aquifer properties. Therefore, multiple sources of information are needed to complement the hydrogeology understanding of a region. Here, this study presents a numerical framework in which information from different measurement sources is combined to characterize the 3D random field in a multi-fidelity prediction model. Coupled with the model, a Bayesian experimental design was used to determine the best future sampling locations. The Upper Sangamon watershed in east-central Illinois was selected as the case study site, where the multi-fidelity Gaussian process model was used to estimate the hydraulic conductivity in the region of interest. Multi-source observation data were obtained from electrical resistivity and borehole pumping tests. The accuracy of the model prediction is dependent on the locations and the distribution of both high- and low-fidelity data. Furthermore, the multi-fidelity model was compared with the single-fidelity model. The uncertainties and confidence in the measurements and parameter estimates were quantified and used to design future cycles of data collection to further improve the confidence intervals.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2015

Assessment of uncertainties in radiation-induced cancer risk predictions at clinically relevant doses

Purpose: Theoretical dose–response models offer the possibility to assess second cancer induction risks after external beam therapy. The parameters used in these models are determined with limited data from epidemiological studies. Risk estimations are thus associated with considerable uncertainties. This study aims at illustrating uncertainties when predicting the risk for organ-specific second cancers in the primary radiation field illustrated by choosing selected treatment plans for brain cancer patients. Methods: A widely used risk model was considered in this study. The uncertainties of the model parameters were estimated with reported data of second cancer incidences for various organs. Standard error propagation was then subsequently applied to assess the uncertainty in the risk model. Next, second cancer risks of five pediatric patients treated for cancer in the head and neck regions were calculated. For each case, treatment plans for proton and photon therapy were designed to estimate the uncertainties (a) in the lifetime attributable risk (LAR) for a given treatment modality and (b) when comparing risks of two different treatment modalities. Results: Uncertainties in excess of 100% of the risk were found for almost all organs considered. When applied to treatment plans, the calculated LAR values have uncertainties of the same magnitude. A comparison between cancer risks of different treatment modalities, however, does allow statistically significant conclusions. In the studied cases, the patient averaged LAR ratio of proton and photon treatments was 0.35, 0.56, and 0.59 for brain carcinoma, brain sarcoma, and bone sarcoma, respectively. Their corresponding uncertainties were estimated to be potentially below 5%, depending on uncertainties in dosimetry. Conclusions: The uncertainty in the dose–response curve in cancer risk models makes it currently impractical to predict the risk for an individual external beam treatment. On the other hand, the ratio of absolute risks between two modalities is less sensitive to the uncertainties in the risk model and can provide statistically significant estimates.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2017

Data-Constrained Projections of Methane Fluxes in a Northern Minnesota Peatland in Response to Elevated CO 2 and Warming

Large uncertainties exist in predicting responses of wetland methane (CH 4 ) fluxes to future climate change. However, sources of the uncertainty have not been clearly identified despite the fact that methane production and emission processes have been extensively explored. In this study, we took advantage of manual CH 4 flux measurements under ambient environment from 2011 to 2014 at the Spruce and Peatland Responses Under Changing Environments (SPRUCE) experimental site and developed a data-informed process-based methane module. The module was incorporated into the Terrestrial ECOsystem (TECO) model before its parameters were constrained with multiple years of methane flux data for forecasting CH 4 emission under five warming and two elevated CO 2 treatments at SPRUCE. We found that 9°C warming treatments significantly increased methane emission by approximately 400%, and elevated CO 2 treatments stimulated methane emission by 10.4%–23.6% in comparison with ambient conditions. The relative contribution of plant-mediated transport to methane emission decreased from 96% at the control to 92% at the 9°C warming, largely to compensate for an increase in ebullition. The uncertainty in plant-mediated transportation and ebullition increased with warming and contributed to the overall changes of emissions uncertainties. At the same time, our modeling results indicated a significant increase in the emitted CH 4 :CO 2 ratio. This result, together with the larger warming potential of CH 4 , will lead to a strong positive feedback from terrestrial ecosystems to climate warming. In conclusion, the model-data fusion approach used in this study enabled parameter estimation and uncertainty quantification for forecasting methane fluxes.

58 GEOSCIENCES↗
DOE OSTI2015

SU-E-T-144: Bayesian Inference of Local Relapse Data Using a Poisson-Based Tumour Control Probability Model

Purpose: The purpose of this project is to develop a robust method of parameter estimation for a Poisson-based TCP model using Bayesian inference. Methods: Bayesian inference was performed using the PyMC3 probabilistic programming framework written in Python. A Poisson-based TCP regression model that accounts for clonogen proliferation was fit to observed rates of local relapse as a function of equivalent dose in 2 Gy fractions for a population of 623 stage-I non-small-cell lung cancer patients. The Slice Markov Chain Monte Carlo sampling algorithm was used to sample the posterior distributions, and was initiated using the maximum of the posterior distributions found by optimization. The calculation of TCP with each sample step required integration over the free parameter α, which was performed using an adaptive 24-point Gauss-Legendre quadrature. Convergence was verified via inspection of the trace plot and posterior distribution for each of the fit parameters, as well as with comparisons of the most probable parameter values with their respective maximum likelihood estimates. Results: Posterior distributions for α, the standard deviation of α (σ), the average tumour cell-doubling time (Td), and the repopulation delay time (Tk), were generated assuming α/β = 10 Gy, and a fixed clonogen density of 10{sup 7} cm−{sup 3}. Posterior predictive plots generated from samples from these posterior distributions are in excellent agreement with the observed rates of local relapse used in the Bayesian inference. The most probable values of the model parameters also agree well with maximum likelihood estimates. Conclusion: A robust method of performing Bayesian inference of TCP data using a complex TCP model has been established.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2025

Updating Nuclear Energy Cost Estimates for Net Zero World Initiative

Energy modeling of decarbonized scenarios in integrated energy systems requires nuclear energy parameters that are critical for forecasting, modeling and cost structure analysis. Using updated real-world data has always been a challenge to estimate current nuclear reactors costs and deployment scenarios. Given this, an updated set of parameters for overnight capital costs and operation and maintenance costs are estimated for the Net Zero World initiative using recent reports that provided a vast set of open sources data inputs. This paper follows the methodology developed in the Net Zero World report and applies the new ranges estimated in the Gateway for Accelerated Innovation in Nuclear report that address many of the current challenges in obtaining accurate cost data for advanced nuclear concepts. The final goal is to provide new estimates of the overnight capital costs and operational costs for different countries. The present paper improves the earlier capital cost estimations, building on recent literature that aims to obtain accurate data for modeling and simulation to enhance energy system evaluations and support decision-making in areas like de-carbonization and capacity expansion. Finally, the paper compares the new cost estimates with the old cost results.

11 - NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2015

SU-E-J-130: Automating Liver Segmentation Via Combined Global and Local Optimization

Purpose: The aim of this work is to develop a robust algorithm for accurate segmentation of liver with special attention paid to the problems with fuzzy edges and tumor. Methods: 200 CT images were collected from radiotherapy treatment planning system. 150 datasets are selected as the panel data for shape dictionary and parameters estimation. The remaining 50 datasets were used as test images. In our study liver segmentation was formulated as optimization process of implicit function. The liver region was optimized via local and global optimization during iterations. Our method consists five steps: 1)The livers from the panel data were segmented manually by physicians, and then We estimated the parameters of GMM (Gaussian mixture model) and MRF (Markov random field). Shape dictionary was built by utilizing the 3D liver shapes. 2)The outlines of chest and abdomen were located according to rib structure in the input images, and the liver region was initialized based on GMM. 3)The liver shape for each 2D slice was adjusted using MRF within the neighborhood of liver edge for local optimization. 4)The 3D liver shape was corrected by employing SSR (sparse shape representation) based on liver shape dictionary for global optimization. Furthermore, H-PSO(Hybrid Particle Swarm Optimization) was employed to solve the SSR equation. 5)The corrected 3D liver was divided into 2D slices as input data of the third step. The iteration was repeated within the local optimization and global optimization until it satisfied the suspension conditions (maximum iterations and changing rate). Results: The experiments indicated that our method performed well even for the CT images with fuzzy edge and tumors. Comparing with physician delineated results, the segmentation accuracy with the 50 test datasets (VOE, volume overlap percentage) was on average 91%–95%. Conclusion: The proposed automatic segmentation method provides a sensible technique for segmentation of CT images. This work is supported by NIH/NIBIB (1R01-EB016777), National Natural Science Foundation of China (No.61471226 and No.61201441), Research funding from Shandong Province (No.BS2012DX038 and No.J12LN23), and Research funding from Jinan City (No.201401221 and No.20120109)

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

Scalable deep learning for watershed model calibration

Watershed models such as the Soil and Water Assessment Tool (SWAT) consist of high-dimensional physical and empirical parameters. These parameters often need to be estimated/calibrated through inverse modeling to produce reliable predictions on hydrological fluxes and states. Existing parameter estimation methods can be time consuming, inefficient, and computationally expensive for high-dimensional problems. In this paper, we present an accurate and robust method to calibrate the SWAT model (i.e., 20 parameters) using scalable deep learning (DL). We developed inverse models based on convolutional neural networks (CNN) to assimilate observed streamflow data and estimate the SWAT model parameters. Scalable hyperparameter tuning is performed using high-performance computing resources to identify the top 50 optimal neural network architectures. We used ensemble SWAT simulations to train, validate, and test the CNN models. We estimated the parameters of the SWAT model using observed streamflow data and assessed the impact of measurement errors on SWAT model calibration. We tested and validated the proposed scalable DL methodology on the American River Watershed, located in the Pacific Northwest-based Yakima River basin. Our results show that the CNN-based calibration is better than two popular parameter estimation methods (i.e., the generalized likelihood uncertainty estimation [GLUE] and the dynamically dimensioned search [DDS], which is a global optimization algorithm). For the set of parameters that are sensitive to the observations, our proposed method yields narrower ranges than the GLUE method but broader ranges than values produced using the DDS method within the sampling range even under high relative observational errors. The SWAT model calibration performance using the CNNs, GLUE, and DDS methods are compared using R 2 and a set of efficiency metrics, including Nash-Sutcliffe, logarithmic Nash-Sutcliffe, Kling-Gupta, modified Kling-Gupta, and non-parametric Kling-Gupta scores, computed on the observed and simulated watershed responses. The best CNN-based calibrated set has scores of 0.71, 0.75, 0.85, 0.85, 0.86, and 0.91. The best DDS-based calibrated set has scores of 0.62, 0.69, 0.8, 0.77, 0.79, and 0.82. The best GLUE-based calibrated set has scores of 0.56, 0.58, 0.71, 0.7, 0.71, and 0.8. The scores above show that the CNN-based calibration leads to more accurate low and high streamflow predictions than the GLUE and DDS sets. Our research demonstrates that the proposed method has high potential to improve our current practice in calibrating large-scale integrated hydrologic models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2013

Spatial and Temporal Patterns of Global Onshore Wind Speed Distribution

Wind power, a renewable energy source, can play an important role in electrical energy generation. Information regarding wind energy potential is important both for energy related modeling and for decision-making in the policy community. While wind speed datasets with high spatial and temporal resolution are often ultimately used for detailed planning, simpler assumptions are often used in analysis work. An accurate representation of the wind speed frequency distribution is needed in order to properly characterize wind energy potential. Using a power density method, this study estimated global variation in wind parameters as fitted to a Weibull density function using NCEP/CFSR reanalysis data. The estimated Weibull distribution performs well in fitting the time series wind speed data at the global level according to R2, root mean square error, and power density error. The spatial, decadal, and seasonal patterns of wind speed distribution were then evaluated. We also analyzed the potential error in wind power estimation when a commonly assumed Rayleigh distribution (Weibull k = 2) is used. We find that the assumption of the same Weibull parameter across large regions can result in substantial errors. While large-scale wind speed data is often presented in the form of average wind speeds, these results highlight the need to also provide information on the wind speed distribution.

Zhou, Yuyu↗
DOE OSTI2025

Data-Efficient Methods for Determining Flory–Huggins χ Parameters in Multicomponent Polymer Formulations

Polymer formulations are essential in diverse applications including personal care products, coatings, paints, adhesives, and plastic materials. Designing these formulations requires navigating large, complex design spaces, where phase and self-assembly behavior critically impact performance. The Flory–Huggins χ parameter, which quantifies segmental miscibility, is widely used to parametrize the excess free energy of mixing in formulation models. In this work, we introduce two data-efficient, top-down methods for estimating χ parameters using the Random Phase Approximation (RPA): (i) Boundary Nonlinear Regression (Boundary-NLR), which fits theoretical spinodal boundaries to experimental phase boundaries, and (ii) Surrogate Model Inverse Parameter Estimation (SMIPE), which uses a Gaussian Process Classifier to fit sparse phase maps via a surrogate model. Both methods allow rapid parametrization of polymer field-theoretic models without the need for additional experiments. We evaluate these approaches on data sets involving polymer–solvent–nonsolvent ternary mixtures and block copolymer–solvent systems, demonstrating their robustness to experimental noise and their relevance for real-world formulation design.

copolymers↗