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 379 records · Page 21
All sourcesNASA NTRSDOE OSTI
DOE OSTI2017

A new probe concept for internal pipework inspection - 241

The interior visual inspection of nuclear pipework is a critical inspection activity required to ensure the continued safe, reliable operation of plant and thus avoid costly outages. Typically, the video output from a manually deployed probe is viewed by an operator online with the task of identifying and estimating the location of surface defects such as cracks, corrosion and pitting. However, it is very difficult to estimate the nature and spatial extent of defects from the often disorientating small field of view video of a relatively large structure. This work describes a new visual inspection system incorporating a wide field of view machine vision camera and additional sensors designed for inspecting 3 - 6 inch diameter pipes. The output of the system is a photorealistic model of the internal surface of the pipework. The generation of this model relies upon a core component of the system in the form of image feature extraction which estimates camera location. This paper considers the accuracy of this estimation as a function of the number and configuration of the extracted image features. (authors)

42 ENGINEERING↗
DOE OSTI2024

Advancing earth system model calibration: a diffusion-based method

Understanding of climate impact on ecosystems globally requires site-specific model calibration. Here we introduce a novel diffusion-based uncertainty quantification (DBUQ) method for efficient model calibration. DBUQ is a score-based diffusion model that leverages Monte Carlo simulation to estimate the score function and evaluates a simple neural network to quickly generate samples for approximating parameter posterior distributions. DBUQ is stable, efficient, and can effectively calibrate the model given diverse observations, thereby enabling rapid and site-specific model calibration on a global scale. This capability significantly advances Earth system modeling and our understanding of climate impacts on Earth systems. We demonstrate DBUQ's capability in E3SM land model calibration at the Missouri Ozark AmeriFlux forest site. Both synthetic and real-data applications indicate that DBUQ produces accurate parameter posterior distributions similar to those generated by Markov Chain Monte Carlo sampling but with 30X less computing time. This efficiency marks a significant stride in model calibration, paving the way for more effective and timely climate impact analyses.

Liu, Yanfang↗
DOE OSTI2026

How does drought affect residential water demand and price elasticity?

Urban water scarcity is an important social and economic concern, particularly as the intensity, duration, and frequency of droughts is increasing in many regions. We consider whether drought induces changes to water demand and the price elasticity of demand for water that may last beyond a drought’s official end date. If drought shocks prompt long-term changes in water demand behavior, and these changes occur at broad geographic scale, they could have important implications for modeling adaptive responses to water scarcity. We assemble a novel dataset on residential water demand and pricing in the western United States to test empirically for effects of drought on water demand and price elasticity. We perform our analysis with aggregate quantity, price, and drought data, accounting for endogenous prices under increasing-block water tariffs and using both average and marginal water fees in estimating water demand functions. Results are consistent with the hypothesis that households may become less price-sensitive after exposure to drought. However, we find no systematic evidence of long-run, drought-related reductions in water demand, itself.

demand hardening↗
DOE OSTI2022

Weakly Supervised Event Classification Using Imperfect Real-world PMU Data with Scarce Labels

This paper studies event classification using imperfect real-world phasor measurement unit (PMU) data with scarce event types (labels). By investigating the real-world PMU data, it is observed that most real-world PMU data's event type is unknown, which makes it challenging to directly use such dataset to build event classifiers as existing classification techniques require high-quality training data with known event type (i.e., label). To address this challenge, a weakly supervised learning based event classification approach is developed, which can use noisy and low-quality PMU data for the training. First, data quality issues are fixed using data preprocessing techniques and then event features are constructed from the PMU data. Using these features, a series of labeling functions are learnt to generate initial estimates of the labels of large amounts of unlabeled PMU data. As the labeling functions are learnt using the same data with scarce labels, the label estimates from the labeling functions can be correlated, noisy, and bias. To enhance these initial estimates, a generative model is developed to characterize the dependencies among the estimated labels, based on which better labels are obtained for training event classifiers. Numerical experiments using the real-world dataset from the Western Interconnection of the U.S. power transmission grid show that the proposed weakly supervised event classifier trained using the dataset with only 5% labeled data can achieve 78.4% classification accuracy.

Liu, Yunchuan↗
DOE OSTI2024

Estimation of wintertime cloud radiative effects in the Western Arctic, a function of cloud-moisture-coupling and sea ice conditions

Cloud radiative effects are estimated using up- and downward ground-based observations of broadband thermal infrared sensors in cloudy conditions which are provided by the ARM data repository. The ARM Surface Radiation measurement is analyzed. To study the cloud radiative effects, the observed thermal infrared measurements under cloudy conditions are analyzed together with microphysical cloud properties like liquid and ice water path obtained from the CloudNet algorithms applied to the remote-sensing instrumentation suite at the ARM site in Utqiagvik. Cloud observations have been sorted depending on the coupling of the cloud with upwind water vapor transport events that can interact with sea ice upwind. Statistics of surface net thermal fluxes and cloud water/ice properties show a distinguished distribution for coupled and decoupled clouds. It was found that coupled clouds contain more liquid water and are more efficient reduce the cooling of surface LW radiation.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2016

Large covariance matrices: smooth models from the two-point correlation function

We introduce a new method for estimating the covariance matrix for the galaxy correlation function in surveys of large-scale structure. Our method combines simple theoretical results with a realistic characterization of the survey to dramatically reduce noise in the covariance matrix. For example, with an investment of only ≈1000 CPU hours we can produce a model covariance matrix with noise levels that would otherwise require ~35 000 mocks. Non-Gaussian contributions to the model are calibrated against mock catalogues, after which the model covariance is found to be in impressive agreement with the mock covariance matrix. Since calibration of this method requires fewer mocks than brute force approaches, we believe that it could dramatically reduce the number of mocks required to analyse future surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2019

Dwarf spheroidal J-factor likelihoods for generalized NFW profiles

ABSTRACT Indirect detection strategies of particle dark matter (DM) in Dwarf spheroidal satellite galaxies (dSphs) typically entail searching for annihilation signals above the astrophysical background. To robustly compare model predictions with the observed fluxes of product particles, most analyses of astrophysical data – which are generally frequentist – rely on estimating the abundance of DM by calculating the so-called J factor. This quantity is usually inferred from the kinematic properties of the stellar population of a dSph using the Jeans equation, commonly by means of Bayesian techniques that entail the presence (and additional systematic uncertainty) of prior choice. Here, extending earlier work, we develop a scheme to derive the profile likelihood for J factors of dwarf spheroidals for models with five or more free parameters. We validate our method on a publicly available simulation suite, released by the Gaia Challenge, finding satisfactory statistical properties for bias and probability coverage. We present the profile likelihood function and maximum likelihood estimates for the J-factor of 10 dSphs. As an illustration, we apply these profile likelihoods to recently published analyses of γ-ray data with the Fermi Large Area Telescope to derive new, consistent upper limits on the DM annihilation cross-section. We do this for a subset of systems, generally referred to as classical dwarfs. The implications of these findings for DM searches are discussed, together with future improvements and extensions of this technique.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2014

Geological Carbon Sequestration Storage Resource Estimates for the Ordovician St. Peter Sandstone, Illinois and Michigan Basins, USA

The Cambro-Ordovician strata of the Midwest of the United States is a primary target for potential geological storage of CO2 in deep saline formations. The objective of this project is to develop a comprehensive evaluation of the Cambro-Ordovician strata in the Illinois and Michigan Basins above the basal Mount Simon Sandstone since the Mount Simon is the subject of other investigations including a demonstration-scale injection at the Illinois Basin Decatur Project. The primary reservoir targets investigated in this study are the middle Ordovician St Peter Sandstone and the late Cambrian to early Ordovician Knox Group carbonates. The topic of this report is a regional-scale evaluation of the geologic storage resource potential of the St Peter Sandstone in both the Illinois and Michigan Basins. Multiple deterministic-based approaches were used in conjunction with the probabilistic-based storage efficiency factors published in the DOE methodology to estimate the carbon storage resource of the formation. Extensive data sets of core analyses and wireline logs were compiled to develop the necessary inputs for volumetric calculations. Results demonstrate how the range in uncertainty of storage resource estimates varies as a function of data availability and quality, and the underlying assumptions used in the different approaches. In the simplest approach, storage resource estimates were calculated from mapping the gross thickness of the formation and applying a single estimate of the effective mean porosity of the formation. Results from this approach led to storage resource estimates ranging from 3.3 to 35.1 Gt in the Michigan Basin, and 1.0 to 11.0 Gt in the Illinois Basin at the P10 and P90 probability level, respectively. The second approach involved consideration of the diagenetic history of the formation throughout the two basins and used depth-dependent functions of porosity to derive a more realistic spatially variable model of porosity rather than applying a single estimate of porosity throughout the entire potential reservoir domains. The second approach resulted in storage resource estimates of 3.0 to 31.6 Gt in the Michigan Basin, and 0.6 to 6.1 Gt in the Illinois Basin. The third approach attempted to account for the local-scale variability in reservoir quality as a function of both porosity and permeability by using core and log analyses to calculate explicitly the net effective porosity at multiple well locations, and interpolate those results throughout the two basins. This approach resulted in storage resource estimates of 10.7 to 34.7 Gt in the Michigan Basin, and 11.2 to 36.4 Gt in the Illinois Basin. A final approach used advanced reservoir characterization as the most sophisticated means to estimating storage resource by defining reservoir properties for multiple facies within the St Peter formation. This approach was limited to the Michigan Basin since the Illinois Basin data set did not have the requisite level of data quality and sampling density to support such an analysis. Results from this approach led to storage resource estimates of 15.4 Gt to 50.1 Gt for the Michigan Basin. The observed variability in results from the four different approaches is evaluated in the context of data and methodological constraints, leading to the conclusion that the storage resource estimates from the first two approaches may be conservative, whereas the net porosity based approaches may over-estimate the resource.

Barnes, David↗
DOE OSTI2021

PSF Estimation in Crowded Astronomical Imagery as a Convolutional Dictionary Learning Problem

In this paper, we present a new algorithm for estimating the Point Spread Function (PSF) in wide-field astronomical images with extreme source crowding. Robust and accurate PSF estimation in crowded astronomical images dramatically improves the fidelity of astrometric and photometric measurements extracted from wide-field sky monitoring imagery. Our radically new approach utilizes convolutional sparse representations to model the continuous functions involved in the image formation. This approach avoids the need to detect and precisely localize individual point sources that is shared by existing methods. In experiments involving simulated astronomical imagery, it significantly outperforms the recent alternative method with which it is compared.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2013

KAPPA DISTRIBUTION MODEL FOR HARD X-RAY CORONAL SOURCES OF SOLAR FLARES

Solar flares produce hard X-ray emission, the photon spectrum of which is often represented by a combination of thermal and power-law distributions. However, the estimates of the number and total energy of non-thermal electrons are sensitive to the determination of the power-law cutoff energy. Here, we revisit an 'above-the-loop' coronal source observed by RHESSI on 2007 December 31 and show that a kappa distribution model can also be used to fit its spectrum. Because the kappa distribution has a Maxwellian-like core in addition to a high-energy power-law tail, the emission measure and temperature of the instantaneous electrons can be derived without assuming the cutoff energy. Moreover, the non-thermal fractions of electron number/energy densities can be uniquely estimated because they are functions of only the power-law index. With the kappa distribution model, we estimated that the total electron density of the coronal source region was {approx}2.4 Multiplication-Sign 10{sup 10} cm{sup -3}. We also estimated without assuming the source volume that a moderate fraction ({approx}20%) of electrons in the source region was non-thermal and carried {approx}52% of the total electron energy. The temperature was 28 MK, and the power-law index {delta} of the electron density distribution was -4.3. These results are compared to the conventional power-law models with and without a thermal core component.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2024

A semi-supervised learning method to produce explainable radioisotope proportion estimates for NaI-based synthetic and measured gamma spectra

Quantifying the radioactive sources present in gamma spectra is an ever-present and growing national security mission and a time-consuming process for human analysts. While machine learning models exist that are trained to estimate radioisotope proportions in gamma spectra, few address the eventual need to provide explanatory outputs beyond the estimation task. In this work, we develop two machine learning models for a NaI detector measurements: one to perform the estimation task, and the other to characterize the first model’s ability to provide reasonable estimates. To ensure the first model exhibits a behavior that can be characterized by the second model, the first model is trained using a custom, semi-supervised loss function which constrains proportion estimates to be explainable in terms of a spectral reconstruction. The second auxiliary model is an out-of-distribution detection function (a type of meta-model) leveraging the proportion estimates of the first model to identify when a spectrum is sufficiently unique from the training domain and thus is out-of-scope for the model. In demonstrating the efficacy of this approach, we encourage the use of meta-models to better explain ML outputs used in radiation detection and increase trust.

46 INSTRUMENTATION RELATED TO NUCLEAR SCIENCE AND ↗
DOE OSTI2015

Standard Model thermodynamics across the electroweak crossover

Even though the Standard Model with a Higgs mass m{sub \tiny H}=125 GeV possesses no bulk phase transition, its thermodynamics still experiences a “soft point” at temperatures around T=160 GeV, with a deviation from ideal gas thermodynamics. Such a deviation may have an effect on precision computations of weakly interacting dark matter relic abundances if their mass is in the few TeV range, or on leptogenesis scenarios operating in this temperature range. By making use of results from lattice simulations based on a dimensionally reduced effective field theory, we estimate the relevant thermodynamic functions across the crossover. The results are tabulated in a numerical form permitting for their insertion as a background equation of state into cosmological particle production/decoupling codes. We find that Higgs dynamics induces a non-trivial “structure” visible e.g. in the heat capacity, but that in general the largest radiative corrections originate from QCD effects, reducing the energy density by a couple of percent from the free value even at T>160 GeV.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2014

Model-based PSF and MTF estimation and validation from skeletal clinical CT images

Purpose: A method was developed to correct for systematic errors in estimating the thickness of thin bones due to image blurring in CT images using bone interfaces to estimate the point-spread-function (PSF). This study validates the accuracy of the PSFs estimated using said method from various clinical CT images featuring cortical bones. Methods: Gaussian PSFs, characterized by a different extent in the z (scan) direction than in the x and y directions were obtained using our method from 11 clinical CT scans of a cadaveric craniofacial skeleton. These PSFs were estimated for multiple combinations of scanning parameters and reconstruction methods. The actual PSF for each scan setting was measured using the slanted-slit technique within the image slice plane and the longitudinal axis. The Gaussian PSF and the corresponding modulation transfer function (MTF) are compared against the actual PSF and MTF for validation. Results: The differences (errors) between the actual and estimated full-width half-max (FWHM) of the PSFs were 0.09 ± 0.05 and 0.14 ± 0.11 mm for the xy and z axes, respectively. The overall errors in the predicted frequencies measured at 75%, 50%, 25%, 10%, and 5% MTF levels were 0.06 ± 0.07 and 0.06 ± 0.04 cycles/mm for the xy and z axes, respectively. The accuracy of the estimates was dependent on whether they were reconstructed with a standard kernel (Toshiba's FC68, mean error of 0.06 ± 0.05 mm, MTF mean error 0.02 ± 0.02 cycles/mm) or a high resolution bone kernel (Toshiba's FC81, PSF FWHM error 0.12 ± 0.03 mm, MTF mean error 0.09 ± 0.08 cycles/mm). Conclusions: The method is accurate in 3D for an image reconstructed using a standard reconstruction kernel, which conforms to the Gaussian PSF assumption but less accurate when using a high resolution bone kernel. The method is a practical and self-contained means of estimating the PSF in clinical CT images featuring cortical bones, without the need phantoms or any prior knowledge about the scanner-specific parameters.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2013

FAR-ULTRAVIOLET OBSERVATIONS OF THE TAURUS-PERSEUS-AURIGA COMPLEX

We have constructed a far-ultraviolet (FUV) continuum map of the Taurus-Auriga-Perseus complex, one of the largest local associations of dark clouds, by merging the two data sets of Galaxy Evolution Explorer and FUV Imaging Spectrograph, which made observations at similar wavelengths. The FUV intensity varies significantly across the whole region, but the diffuse FUV continuum is dominated by dust scattering of stellar photons. A diffuse FUV background of {approx}1000 CU is observed, part of which may be attributable to the scattered photons of foreground FUV light, located in front of the thick clouds. The fluorescent emission of molecular hydrogen constitutes {approx}10% of the total FUV intensity throughout the region, generally proportional to the local continuum level. We have developed a Monte Carlo radiative transfer code and applied it to the present clouds complex to obtain the optical properties of dust grains and the geometrical structures of the clouds. The albedo and the phase function asymmetry factor were estimated to be 0.42{sup +0.05}{sub -0.05}, and 0.47{sup +0.11}{sub -0.27}, respectively, in accordance with theoretical estimations as well as recent observations. The distance and thickness of the four prominent clouds in this complex were estimated using a single-slab model applied individually to each cloud. The results obtained were in good agreement with those from other observations of the Taurus cloud, as its geometrical structure is rather simple. For other clouds that were observed to have multiple components, the results gave distances and thicknesses encompassing all of the components of each cloud. The distance and thickness estimations were not crucially sensitive to the exact values of the albedo and the phase function asymmetry factor, while the locations of the bright field stars relative to the clouds as initial photon sources seem to be the most important factor in the process of fitting.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Joint Bayesian Estimation of Quasar Continua and the Lyα Forest Flux Probability Distribution Function

We present a new Bayesian algorithm making use of Markov Chain Monte Carlo sampling that allows us to simultaneously estimate the unknown continuum level of each quasar in an ensemble of high-resolution spectra, as well as their common probability distribution function (PDF) for the transmitted Lyα forest flux. This fully automated PDF regulated continuum fitting method models the unknown quasar continuum with a linear principal component analysis (PCA) basis, with the PCA coefficients treated as nuisance parameters. The method allows one to estimate parameters governing the thermal state of the intergalactic medium (IGM), such as the slope of the temperature–density relation γ−1, while marginalizing out continuum uncertainties in a fully Bayesian way. Using realistic mock quasar spectra created from a simplified semi-numerical model of the IGM, we show that this method recovers the underlying quasar continua to a precision of ≃7% and ≃10% at z = 3 and z = 5, respectively. Given the number of principal component spectra, this is comparable to the underlying accuracy of the PCA model itself. Most importantly, we show that we can achieve a nearly unbiased estimate of the slope γ−1 of the IGM temperature–density relation with a precision of ±8.6% at z = 3 and ±6.1% at z = 5, for an ensemble of ten mock high-resolution quasar spectra. Applying this method to real quasar spectra and comparing to a more realistic IGM model from hydrodynamical simulations would enable precise measurements of the thermal and cosmological parameters governing the IGM, albeit with somewhat larger uncertainties, given the increased flexibility of the model.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2017

Diffusion coefficients of transition metals in fcc cobalt

Using first-principles density functional theory (DFT), we calculate the diffusivities of 32 different solute elements—all transition metals, together with Al and Si—in fcc cobalt within the formalism of the five-frequency model. For self-diffusion in fcc cobalt, we compare the accuracy of various approximations to the exchange-correlation energy functional of DFT in estimating the activation energy, and find that only the Perdew-Burke-Ernzerhof (PBE) approximation agrees well with experimental reports and all other functionals largely overestimate it. Our calculations also show that an accurate estimation of the self-diffusion coefficient requires explicit calculation of the effective jump frequency and vacancy formation entropy via phonons. Using accurate self-diffusion data and scaling all solute-related attempt frequencies with respect to the attempt frequency for self-diffusion using a simple relation involving the atomic mass and melting temperature of the solute yields solute diffusivities in excellent agreement with experiments, where such data is available. Additionally, we find that large solutes spontaneously relax toward the nearest neighbor vacancy to relieve the misfit strain, and the extent of this relaxation correlates negatively with the migration energy. Thus, in general, larger solutes have lower migration energies and diffuse faster than smaller solutes in fcc cobalt. However, extremely large solutes, e.g., group III elements Sc, Y, Lu, tend to be trapped in an energy valley located halfway toward the vacancy, and monovacancy mediated diffusion may no longer be valid in such cases. Finally, for all the solutes considered, we systematically tabulate the diffusion-related quantities calculated—diffusion prefactors, migration and activation energies—constructing an extensive and accurate first-principles database for solute diffusion in fcc cobalt.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Electrodynamic Shaker Capability Estimation Through Experimental Dynamic Substructuring [Thesis]

Electrodynamic shaker systems are an essential tool in shock and vibration testing of dynamic environments. However, the specific performance capability of these systems is difficult to characterize. The dynamics of the shaker itself, the device under test and the specific test configuration used all couple to create a dynamic response unique to each test. Poorly predicted limitations in shaker capability affect the ability to achieve test specifications, delay testing schedules, and create difficulties for choosing test equipment. To predict shaker capability for a specific test configuration prior to setup, a lumped parameter model of the shaker system and a modal model of a device under test was developed. These models were then analytically coupled using LaGrange multiplier frequency based substructuring to estimate their coupled frequency response functions. The coupled frequency response functions were used to predict electrical inputs required to meet a given test specification. These input requirements were finally compared to a validation test using the specification and setup. Input requirements estimated using the substructuring estimated frequency response functions showed significant error. However, results using an ideal frequency response function showed very little error. These results indicate that with a better method of experimental dynamic substructuring employed it would be possible to accurately predict shaker capability for a given test configuration prior to setup.

42 ENGINEERING↗
DOE OSTI2018

The Effects of Stochastic Velocity Variations on Estimating Time Dependent Seismic Moment Tensors: Applications to the Blue Mountain Well Perforation Data

Resolving the time dependent terms in the seismic moment tensor provides important information that can be used to interpret the source process of an explosion, including the separation of isotropic explosion terms from shear forces and potentially isolated force couples. In this report, we detail our method of inverting three component seismic data for the seismic moment tensor. We review possible seismic source models from the simplest isotropic explosion type source to those incorporating the six independent moment tensor terms. The inversion we describe is formulated in the frequency domain, and results in estimates of time dependent moment tensor components. The inversion relies on an accurate estimate of the Green's functions of the Earth. However, given the complexity of the Earth, we explore the effects of inaccuracies in the presumed Earth model used to estimate the Green's functions needed for the inversion. Specifically, we explore the effects of stochastic variations in the Earth models on the inversion results. These tests are synthetic throughout, and show that adding stochastic density/velocity heterogeneity in the presumed Earth model results in reduced amplitude seismic moment tensor estimates, as well as degrading the data misfit. We suggest two mitigation strategies. First, produce a suite of Green's functions using different realizations of the stochastic field within the Earth Model. Secondly, perform the in- version in the power spectral domain, eliminating all phase information. Finally, we analyze actual seismic data collected in winter 2017/2018. The seismic data was collected at in active geothermal well site outside of Winnimucca, NV, and was produced during well stimulation operations. In general, the inversion results were poor, with a high degree of data misfit. We hypothesize that the poor results are a function of a poorly constrained Earth model as well as noisy, high-frequency data being used in the inversion.

58 GEOSCIENCES↗