Search NASA⌕ Search

SEARCH · Search NASA

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 271 records · Page 15

Spatially Explicit Modeling of Criteria Air Pollutants from Agricultural and Forestry Feedstock Production

The Feedstock Production Emissions to Air Model (FPEAM) is a spatially explicit supply chain model that quantifies criteria air pollutants and precursors generated by biomass production, harvest and transportation to biorefineries. FPEAM was most recently used for the 2016 Billion Ton Study (BTS) published by the U.S. Department of Energy Bioenergy Technologies Office. This presentation describes several use cases for the model, and provides an overview of recent restructuring that increases FPEAM's usability and functionality. FPEAM currently estimates emissions of volatile organic compounds, particulate matter, nitrogen oxides, sulfur oxides, and carbon monoxide from on-farm equipment operation, on-farm fugitive dust, on-farm chemical (fertilizer, herbicide and insecticide) application, and biomass transportation. These calculations are implemented within FPEAM as modules, which independently perform calculations for each emissions category. This structure allows users to include or exclude emissions categories according to the goal and scope of their analysis. Each module's input data defines the activities taking place with information including fuel combusted, chemicals applied, and equipment used. Users can replace or augment the default input data to quantify additional pollutant types, emissions categories, farming practices and biomass types. Emissions from biomass transportation and on-farm equipment operation are by default calculated by the U.S. Environmental Protection Agency's Motor Vehicle Emission Simulator (MOVES) and NONROAD model, respectively. These models - which are computationally intensive but provide detailed results - can optionally be replaced with user-provided emissions factors for fuel combustion and vehicle operation. This approach reduces FPEAM's runtime and allows users to run sensitivity analyses around biomass production and transportation scenarios. Other enhancements currently under development are the use of farm-to-biorefinery route information to increase the spatial resolution of transportation emissions, and a module to connect FPEAM and the Intervention Model for Air Pollution (InMAP), which will enable the quantification of human health impacts.

09 BIOMASS FUELS↗

Compressive sensing optical design and simulation tool

A system or method for designing optical systems includes generating a mathematical model of an imaging system; tracing a plurality of light rays from a source through the model of the imaging system in a forward direction and a backward direction incident upon an arbitrarily defined detector location; calculating the total light rays incident upon the arbitrarily defined detector location in response to the plurality of light rays; tracing a series of calibration images first through the computational or compressive imager mathematical model; measuring an output set of parameters based on a plurality of detected rays calibration images at the defined detector location; storing the output set of parameters in software, and determining an estimated system transfer function.

42 ENGINEERING↗

Beyond Generalized Eigenvalues in Lattice Quantum Field Theory

Two analysis techniques, the generalized eigenvalue method (GEM) or Prony's (or related) method (PM), are commonly used to analyze statistical estimates of correlation functions produced in lattice quantum field theory calculations. GEM takes full advantage of the matrix structure of correlation functions but only considers individual pairs of time separations when much more data exists. PM can be applied to many time separations and many individual matrix elements simultaneously but does not fully exploit the matrix structure of the correlation function. We combine both these methods into a single framework based on matrix polynomials. As these algebraic methods are well known for producing extensive spectral information about statistically-noisy data, the method should be paired with some information criteria, like the recently proposed Bayesean model averaging.

Fleming, George T.↗

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↗

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↗

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↗

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↗

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↗

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↗

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↗

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 ↗

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↗

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↗

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↗

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↗

Accurate shear estimation with fourth-order moments

ABSTRACT As imaging surveys progress in exploring the large-scale structure of the Universe through the use of weak gravitational lensing, achieving sub-per cent accuracy in estimating shape distortions caused by lensing, or shear, is imperative for precision cosmology. In this paper, we extend the Fourier power function shapelets (FPFS) shear estimator using fourth-order shapelet moments and combine it with the original second-order shear estimator to reduce galaxy shape noise. We calibrate this novel shear estimator analytically to a sub-per cent level-accuracy using the AnaCal framework. This higher order shear estimator is tested with realistic image simulations, and after analytical correction for the detection/selection bias and noise bias, the multiplicative shear bias $|m|$ is below $3\times 10^{-3}$ (99.7 per cent confidence interval) for both isolated and blended galaxies. Once combined with the second-order FPFS shear estimator, the shape noise is reduced by $\sim 35~{{\ \rm per\ cent}}$ for isolated galaxies in simulations with Hyper Suprime-Cam and Vera C. Rubin Observatory Legacy Survey of Space and Time observational conditions. However, for blended galaxies, the effective number density does not significantly improve with the combination of the two estimators. Based on these results, we recommend exploration of how this framework can further reduce the systematic uncertainties in shear due to point spread function leakage and modelling error, and potentially provide improved precision in shear inference in high-resolution space-based images.

79 ASTRONOMY AND ASTROPHYSICS↗

A complete non-perturbative renormalization prescription for quasi-PDFs

In this work we present, for the first time, the non-perturbative renormalization for the unpolarized, helicity and transversity quasi-PDFs, in an RI' scheme. The proposed prescription addresses simultaneously all aspects of renormalization: logarithmic divergences, finite renormalization as well as the linear divergence which is present in the matrix elements of fermion operators with Wilson lines. Furthermore, for the case of the unpolarized quasi-PDF, we describe how to eliminate the unwanted mixing with the twist-3 scalar operator. We utilize perturbation theory for the one-loop conversion factor that brings the renormalization functions to the $\overline{MS}$-scheme at a scale of 2 GeV. We also explain how to improve the estimates on the renormalization functions by eliminating lattice artifacts. The latter can be computed in one-loop perturbation theory and to all orders in the lattice spacing. We apply the methodology for the renormalization to an ensemble of twisted mass fermions with $N_f = 2+1+1$ dynamical quarks, and a pion mass of around 375 MeV.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗