Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Function 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,333 records · Page 74
All sourcesNASA NTRSDOE OSTI
DOE OSTI2013

On estimating cosmology-dependent covariance matrices

We describe a statistical model to estimate the covariance matrix of matter tracer two-point correlation functions with cosmological simulations. Assuming a fixed number of cosmological simulation runs, we describe how to build a 'statistical emulator' of the two-point function covariance over a specified range of input cosmological parameters. Because the simulation runs with different cosmological models help to constrain the form of the covariance, we predict that the cosmology-dependent covariance may be estimated with a comparable number of simulations as would be needed to estimate the covariance for fixed cosmology. Our framework is a necessary first step in planning a simulations campaign for analyzing the next generation of cosmological surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

Intensity statistics in the presence of translational noncrystallographic symmetry

In the case of translational noncrystallographic symmetry (tNCS), two or more copies of a component in the asymmetric unit of the crystal are present in a similar orientation. This causes systematic modulations of the reflection intensities in the diffraction pattern, leading to problems with structure determination and refinement methods that assume, either implicitly or explicitly, that the distribution of intensities is a function only of resolution. To characterize the statistical effects of tNCS accurately, it is necessary to determine the translation relating the copies, any small rotational differences in their orientations, and the size of random coordinate differences caused by conformational differences. An algorithm to estimate these parameters and refine their values against a likelihood function is presented, and it is shown that by accounting for the statistical effects of tNCS it is possible to unmask the competing statistical effects of twinning and tNCS and to more robustly assess the crystal for the presence of twinning.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2016

GeoT User’s Guide, A Computer Program for Multicomponent Geothermometry and Geochemical Speciation, Version 2.1

GeoT implements the multicomponent geothermometry method developed by Reed and Spycher (1984, Geochim. Cosmichim. Acta 46 513–528) into a stand-alone computer program, to ease the application of this method and to improve the prediction of geothermal reservoir temperatures using full and integrated chemical analyses of geothermal fluids. Reservoir temperatures are estimated from statistical analyses of mineral saturation indices computed as a function of temperature. The reconstruction of the deep geothermal fluid compositions, and geothermometry computations, are all implemented into the same computer program, allowing unknown or poorly constrained input parameters to be estimated by numerical optimization using existing parameter estimation software, such as iTOUGH2, PEST, or UCODE. This integrated geothermometry approach presents advantages over classical geothermometers for fluids that have not fully equilibrated with reservoir minerals and/or that have been subject to processes such as dilution and gas loss.

58 GEOSCIENCES↗
DOE OSTI2016

iGeoT v1.0: Automatic Parameter Estimation for Multicomponent Geothermometry, User's Guide

GeoT implements the multicomponent geothermometry method developed by Reed and Spycher [1984] into a stand-alone computer program to ease the application of this method and to improve the prediction of geothermal reservoir temperatures using full and integrated chemical analyses of geothermal fluids. Reservoir temperatures are estimated from statistical analyses of mineral saturation indices computed as a function of temperature. The reconstruction of the deep geothermal fluid compositions, and geothermometry computations, are all implemented into the same computer program, allowing unknown or poorly constrained input parameters to be estimated by numerical optimization. This integrated geothermometry approach presents advantages over classical geothermometers for fluids that have not fully equilibrated with reservoir minerals and/or that have been subject to processes such as dilution and gas loss. This manual contains installation instructions for iGeoT, and briefly describes the input formats needed to run iGeoT in Automatic or Expert Mode. An example is also provided to demonstrate the use of iGeoT.

58 GEOSCIENCES↗
DOE OSTI2022

Monitoring Seasonal Shear Wave Velocity Changes in the Top 6 m at Garner Valley in Southern California With Borehole Data

Abstract Subsurface structures play important roles in seismic ground motion, crustal hydrology, stability of the built environment, and more. Constraining temporal changes of subsurface shear wave velocity ( V S ) can provide useful information to all these topics and the growing field of hydrological monitoring with seismic velocity. Using borehole records at Garner Valley, CA, we estimate seasonal subsurface V S variations from impulse response functions (IRFs) of earthquake data (2005–2018) along with IRFs and cross‐correlation of cross‐hole experiment data (2015–2018). The inferred V S variations are up to ∼25% in the top 6 m and ∼10% at 2–5 m in depth. The V S variations correlate strongly with the water table depth changes, suggesting that the changes are mostly due to fluctuations of pore pressure in the shallow material. The shallow velocity changes alter the near‐surface conditions, can affect seismic hazard estimation, and may be improperly attributed to deeper processes without careful analysis.

58 GEOSCIENCES↗
DOE OSTI2014

A convex pseudolikelihood framework for high dimensional partial correlation estimation with convergence guarantees

Sparse high dimensional graphical model selection is a topic of much interest in modern day statistics. A popular approach is to apply l1-penalties to either parametric likelihoods, or regularized regression/pseudolikelihoods, with the latter having the distinct advantage that they do not explicitly assume Gaussianity. As none of the popular methods proposed for solving pseudolikelihood-based objective functions have provable convergence guarantees, it is not clear whether corresponding estimators exist or are even computable, or if they actually yield correct partial correlation graphs. We propose a new pseudolikelihood-based graphical model selection method that aims to overcome some of the shortcomings of current methods, but at the same time retain all their respective strengths. In particular, we introduce a novel framework that leads to a convex formulation of the partial covariance regression graph problem, resulting in an objective function comprised of quadratic forms. The objective is then optimized via a coordinate wise approach. The specific functional form of the objective function facilitates rigorous convergence analysis leading to convergence guarantees; an important property that cannot be established by using standard results, when the dimension is larger than the sample size, as is often the case in high dimensional applications. These convergence guarantees ensure that estimators are well defined under very general conditions and are always computable. In addition, the approach yields estimators that have good large sample properties and also respect symmetry. Furthermore, application to simulated and real data, timing comparisons and numerical convergence is demonstrated. We also present a novel unifying framework that places all graphical pseudolikelihood methods as special cases of a more general formulation, leading to important insights.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2024

A Nonlinear Least Squares Phasor Estimation Algorithm with a Trust Metric

The paper presents a separable nonlinear least squares (NLLS)-based approach for estimation of fundamental-frequency phasors from sampled point-on-wave measurements. An analytical connection is established between the NLLS cost function and the discrete Fourier transform (DFT)-based periodogram of the input signal. This periodogram-based interpretation of the cost function offers an intuitive and easy-to-implement solution for the frequency estimate. Using the residual error of the NLLS-fit, the paper also presents an insightful measure for ascertaining the quality of the phasor estimates and their validity, especially for data windows containing signal transients.

Chatterjee, Kaustav↗
DOE OSTI2018

Elastic functional principal component regression

We study regression using functional predictors in situations where these functions contains both phase and amplitude variability. In other words, the functions are misaligned due to errors in time measurements, and these errors can significantly degrade both model estimation and prediction performance. The current techniques either ignore the phase variability, or handle it via preprocessing, that is, use an off–the–shelf technique for functional alignment and phase removal. We develop a functional principal component regression model which has a comprehensive approach in handling phase and amplitude variability. The model utilizes a mathematical representation of the data known as the square–root slope function. These functions preserve the L 2 norm under warping and are ideally suited for simultaneous estimation of regression and warping parameters. Furthermore, using both simulated and real–world data sets, we demonstrate our approach and evaluate its prediction performance relative to current models. In addition, we propose an extension to functional logistic and multinomial logistic regression.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2012

Cost estimation for solid waste management in industrialising regions - Precedents, problems and prospects

Highlights: Black-Right-Pointing-Pointer We review cost estimation approaches for solid waste management. Black-Right-Pointing-Pointer Unit cost method and benchmarking techniques used in industrialising regions (IR). Black-Right-Pointing-Pointer Variety in scope, quality and stakeholders makes cost estimation challenging in IR. Black-Right-Pointing-Pointer Integrate waste flow and cost models using cost functions to improve cost planning. - Abstract: The importance of cost planning for solid waste management (SWM) in industrialising regions (IR) is not well recognised. The approaches used to estimate costs of SWM can broadly be classified into three categories - the unit cost method, benchmarking techniques and developing cost models using sub-approaches such as cost and production function analysis. These methods have been developed into computer programmes with varying functionality and utility. IR mostly use the unit cost and benchmarking approach to estimate their SWM costs. The models for cost estimation, on the other hand, are used at times in industrialised countries, but not in IR. Taken together, these approaches could be viewed as precedents that can be modified appropriately to suit waste management systems in IR. The main challenges (or problems) one might face while attempting to do so are a lack of cost data, and a lack of quality for what data do exist. There are practical benefits to planners in IR where solid waste problems are critical and budgets are limited.

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

Optical properties of alkali halide crystals from all-electron hybrid TD-DFT calculations

We present a study of the electronic and optical properties of a series of alkali halide crystals AX, with A = Li, Na, K, Rb and X = F, Cl, Br based on a recent implementation of hybrid-exchange time-dependent density functional theory (TD-DFT) (TD-B3LYP) in the all-electron Gaussian basis set code CRYSTAL. We examine, in particular, the impact of basis set size and quality on the prediction of the optical gap and exciton binding energy. The formation of bound excitons by photoexcitation is observed in all the studied systems and this is shown to be correlated to specific features of the Hartree-Fock exchange component of the TD-DFT response kernel. All computed optical gaps and exciton binding energies are however markedly below estimated experimental and, where available, 2-particle Green’s function (GW-Bethe-Salpeter equation, GW-BSE) values. We attribute this reduced exciton binding to the incorrect asymptotics of the B3LYP exchange correlation ground state functional and of the TD-B3LYP response kernel, which lead to a large underestimation of the Coulomb interaction between the excited electron and hole wavefunctions. Considering LiF as an example, we correlate the asymptotic behaviour of the TD-B3LYP kernel to the fraction of Fock exchange admixed in the ground state functional c{sub HF} and show that there exists one value of c{sub HF} (∼0.32) that reproduces at least semi-quantitatively the optical gap of this material.

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

Electric train energy consumption modeling

For this paper we develop an electric train energy consumption modeling framework considering instantaneous regenerative braking efficiency in support of a rail simulation system. The model is calibrated with data from Portland, Oregon using an unconstrained non-linear optimization procedure, and validated using data from Chicago, Illinois by comparing model predictions against the National Transit Database (NTD) estimates. The results demonstrate that regenerative braking efficiency varies as an exponential function of the deceleration level, rather than an average constant as assumed in previous studies. The model predictions are demonstrated to be consistent with the NTD estimates, producing a predicted error of 1.87% and -2.31%. The paper demonstrates that energy recovery reduces the overall power consumption by 20% for the tested Chicago route. Furthermore, the paper demonstrates that the proposed modeling approach is able to capture energy consumption differences associated with train, route and operational parameters, and thus is applicable for project-level analysis. The model can be easily implemented in traffic simulation software, used in smartphone applications and eco-transit programs given its fast execution time and easy integration in complex frameworks.

32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
DOE OSTI2014

SU-E-T-333: Towards Customizable Radiotherapy Enhancement (CuRE) for Prostate Cancer Using Cisplatin Nanoparticles

Purpose: Replacing routinely used brachytherapy spacers with multifunctional ones loaded with cisplatin nanoparticles (CNP), which can be released into the tumor after implantation, could enable customizable radiation boosting to the prostate tumor in addition to chemotherapy effect. This study investigates the feasibility of customizing the intra-tumor biodistribution and corresponding dose enhancement (DEF) over time for the released CNP as a function of nanoparticle size. Methods: Dose enhancement factors (DEF) due to photon-induced emission of photo-/Auger electrons from CNPs were calculated as a function of concentration using previously published analytical calculation method. An experimentally determined diffusion coefficient (D) for 10 nm nanoparticles in mouse tumor model was employed to estimate D for other sizes using the Stoke- Einstein equation. The error function diffusion model in the experimental study was applied to generate the intra-tumor concentration profile for a burst release of CNPs from the spacer over time. The corresponding DEF profiles were then determined for brachytherapy using Pd-103 and I-125 sources. Results: As expected, the generated profiles showed greater DEF over time for smaller CNP sizes at sample distances from the spacer. For example, for a centrally located spacer, clinically significant DEF (> 20%) could be achieved near the tumor periphery (ca. 0.85 cm distance from the spacer for average PCa tumor size) after 20, and 100 days, respectively for CNPs sizes of 2 nm, and 10 nm, using I-125. Meanwhile for Pd-103, clinically significant DEF could be achieved at the same position after 22 and 108 days, respectively, for same size particles. Conclusion: Our preliminary results demonstrate the feasibility of customizing dose enhancement to prostate tumors as a function of spacer location, brachytherapy source type or size of CNPs released from multifunctional spacers. Such an approach could enable customizable radiation boosting to tumor sub-volumes, while minimizing dose to healthy tissues.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2010

Error handling strategies in multiphase inverse modeling

Parameter estimation by inverse modeling involves the repeated evaluation of a function of residuals. These residuals represent both errors in the model and errors in the data. In practical applications of inverse modeling of multiphase flow and transport, the error structure of the final residuals often significantly deviates from the statistical assumptions that underlie standard maximum likelihood estimation using the least-squares method. Large random or systematic errors are likely to lead to convergence problems, biased parameter estimates, misleading uncertainty measures, or poor predictive capabilities of the calibrated model. The multiphase inverse modeling code iTOUGH2 supports strategies that identify and mitigate the impact of systematic or non-normal error structures. We discuss these approaches and provide an overview of the error handling features implemented in iTOUGH2.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2018

An atomistic investigation of the interaction of dislocations with Guinier-Preston zones in Al-Cu alloys

The interaction between edge dislocations and Guinier-Preston zones in an Al-Cu alloy was analyzed by means of atomistic simulations. The different thermodynamic functions that determine the features of these obstacles for the dislocation glide were computed using molecular statics, molecular dynamics and the nudged elastic band method. It was found that Guinier-Preston zones are sheared by dislocations and the rate at which dislocations overcome the precipitate is controlled by the activation energy, Δ U , in agreement with the postulates of the harmonic transition state theory. Moreover, the entropic contribution to the Helmholtz activation free energy was in the range 1.3–1.8 k b , which can be associated with the typical vibrational entropy of solids. Finally, an estimation of the initial shear flow stress as a function of temperature was carried out from the thermodynamic data provided by the atomistic simulations. Furthermore, comparison with experimental results showed that the effect of the random precipitate distribution and of the dislocation character and dislocation/precipitation orientation has to be taken into account in the simulations to better reproduce experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2019

Impedance-Based Characterization of Power System Frequency Response: Preprint

This paper presents an impedance-based noninvasive approach for the characterization of power system frequency response in real time in the absence of a transient event. It measures transfer function response from the injected active power to the frequency at the point of interconnection for the estimation of inertia and primary frequency response of the system. The so-called frequency response function also shows the effects of the speed of primary frequency control. The paper also derives the relationship between the frequency response function and the network impedance seen from the point of interconnection. This relationship highlights how impedance captures the frequency response behavior of a network. The proposed methodology is demonstrated on a modified IEEE 9-bus system with 25% penetration of wind and PV generation; a 5-MW battery energy storage system is used for the injection of active power perturbations for the measurement of the frequency response function.

24 POWER TRANSMISSION AND DISTRIBUTION↗
DOE OSTI2012

Relationship between x-ray illumination field size and flat field intensity and its impacts on x-ray imaging

Purpose: X-ray cone-beam CT (CBCT) is being increasingly used for various clinical applications, while its performance is still hindered by image artifacts. This work investigates a new source of reconstruction error, which is often overlooked in the current CBCT imaging. The authors find that the x-ray flat field intensity (I{sub 0}) varies significantly as the illumination volume size changes at different collimator settings. A wrong I{sub 0} value leads to inaccurate CT numbers of reconstructed images as well as wrong scatter measurements in the CBCT research. Methods: The authors argue that the finite size of x-ray focal spot together with the detector glare effect cause the I{sub 0} variation at different illumination sizes. Although the focal spot of commercial x-ray tubes typically has a nominal size of less than 1 mm, the off-focal-spot radiation covers an area of several millimeters on the tungsten target. Due to the large magnification factor from the field collimator to the detector, the penumbra effects of the collimator blades result in different I{sub 0} values for different illumination field sizes. Detector glare further increases the variation, since one pencil beam of incident x-ray is scattered into an area of several centimeters on the detector. In this paper, the authors study these two effects by measuring the focal spot distribution with a pinhole assembly and the detector point spread function (PSF) with an edge-spread function method. The authors then derive a formula to estimate the I{sub 0} value for different illumination field sizes, using the measured focal spot distribution and the detector PSF. Phantom studies are carried out to investigate the accuracy of scatter measurements and CT images with and without considering the I{sub 0} variation effects. Results: On our tabletop system with a Varian Paxscan 4030CB flat-panel detector and a Varian RAD-94 x-ray tube as used on a clinical CBCT system, the focal spot distribution has a measured full-width-at-half-maximum (FWHM) of around 0.4 mm, while non-negligible off-focal-spot radiation is observed at a distance of over 2 mm from the center. The measured detector PSF has an FWHM of 0.510 mm, with a shape close to Gaussian. From these two distributions, the author calculate the estimated I{sub 0} values at different collimator settings. The I{sub 0} variation mainly comes from the focal spot effect. The estimation matches well with the measurements at different collimator widths in both horizontal and vertical directions, with an average error of less than 3%. Our method improves the accuracy of conventional scatter measurements, where the scatter is measured as the difference between fan-beam and cone-beam projections. On a uniform water cylinder phantom, more accurate I{sub 0} suppresses the unfaithful high-frequency signals at the object boundaries of the measured scatter, and the SPR estimation error is reduced from 0.158 to 0.014. The proposed I{sub 0} estimation also reduces the reconstruction error from about 20 HU on the Catphan Copyright-Sign 600 phantom in the selected regions of interest to less than 4 HU. Conclusions: The I{sub 0} variation is identified as one additional error source in x-ray imaging. By measuring the focal-spot distribution and detector PSF, the authors propose an accurate method of estimating the I{sub 0} value for different illumination field sizes. The method obtains more accurate scatter measurements and therefore facilitates scatter correction algorithm designs. As correction methods for other CBCT artifacts become more successful, our research is significant in further improving the CBCT imaging accuracy.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2021

Evaluating Gaussian process metamodels and sequential designs for noisy level set estimation

Abstract We consider the problem of learning the level set for which a noisy black-box function exceeds a given threshold. To efficiently reconstruct the level set, we investigate Gaussian process (GP) metamodels. Our focus is on strongly stochastic simulators, in particular with heavy-tailed simulation noise and low signal-to-noise ratio. To guard against noise misspecification, we assess the performance of three variants: (i) GPs with Student- t observations; (ii) Student- t processes (TPs); and (iii) classification GPs modeling the sign of the response. In conjunction with these metamodels, we analyze several acquisition functions for guiding the sequential experimental designs, extending existing stepwise uncertainty reduction criteria to the stochastic contour-finding context. This also motivates our development of (approximate) updating formulas to efficiently compute such acquisition functions. Our schemes are benchmarked by using a variety of synthetic experiments in 1–6 dimensions. We also consider an application of level set estimation for determining the optimal exercise policy of Bermudan options in finance.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2016

On a stochastic approach to a code performance estimation

The main goal of an efficient profiling of software is to minimize the runtime overhead under certain constraints and requirements. The traces built by a profiler during the work, affect the performance of the system itself. One of important aspect of an overhead arises from the randomness of variability in the context in which the application is embedded, e.g., due to possible cache misses, etc. Such uncertainty needs to be taken into account in the design phase. In order to overcome these difficulties we propose to investigate this issue through the analysis of the probability distribution of the difference between profiler’s times for the same code. The approximating model is based on the finite normal mixtures within the framework of the method of moving separation of mixtures. We demonstrate some results for the MATLAB profiler using plotting of 3D surfaces by the function surf. The idea can be used for an estimating of a program efficiency.

97 MATHEMATICS AND COMPUTING↗