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

Fracture Toughness Evaluation for Spent Nuclear Fuel Clad Systems Using Spiral Notch Torsion Fracture Toughness Test

Many radioactive materials within the nuclear fuel cycle present a significant hazard. Such materials include spent nuclear fuel, high-level waste, legacy waste and other nuclear materials. These materials are often held in long-term storage as an interim stage within their lifecycle, including reuse or disposal. Safety and security of spent nuclear fuel (SNF) interim storage installations and follow-on SNF transportations to its final repository sites are very important, due to a great concentration of fission products, actinides and activation products. Fracture mechanics approach in applying to SNF system reliability investigation, especially for the high burn up (HBU) SNF, during SNF long term dry storage or SNF transportation is necessary due to the inherited flaws and inhomogeneity structures existed in a SNF system from the reactor operation. For instance, such as hydride and oxide formation, surface flaw induced by SNF assembly contact interactions, or internal flaws induced by pellet-clad mechanical interaction (PCMI). However, none of the existing fracture toughness data deal with fuel cladding specific geometry or SNF material conditions, such as cladding structure with the segment fuel pellets configuration and the PCMI mechanism. Consequently, the application of existing fracture toughness data to SNF system reliability investigation can be error prone; thus, the development of an industry-wide consensus fracture mechanics approach is needed. The objective of this research is to develop an in-situ fracture mechanics testing protocol, including the associated analytical procedure, that is suitable for evaluating the SNF fracture toughness. Fracture toughness data was obtained under quasi-static fracture loading using Oak Ridge National Laboratory (ORNL) developed Spiral Notch Torsion Test (SNTT) technology carried out on a biaxial tension/torsion tester. These data will be used to support SNF reliability investigation during long-tern SNF dry storage stage or in the follow-on SNF system transport to the final repository sites. SNTT has been a recent breakthrough in measuring fracture toughness for different materials, including metals, ceramics, concrete and polymers composites. Due to its high geometry constraint and unique loading condition, SNTT can be used to measure the fracture toughness with smaller specimens without concern of size effects. The application of SNTT to brittle materials has been proved to be successful. The micro-cracks induced by original notches in brittle materials could ensure crack growth in SNTT samples. Therefore, no fatigue pre cracks are needed for brittle materials specimens. The application of SNTT to the ductile material to generate valid toughness data will require a test sample with sufficient pre-crack length to increase the sample localized constraint at crack front. Fatigue pre-crack growth techniques with the associated compliance function estimated was developed for estimating the crack penetration depth to monitor the fatigue crack growth evolution. In order to extend SNTT approach to a thin shell cladding structure material, a new testing protocol and the associated analytical procedure for estimating the SNF fracture toughness was developed. The detailed SNTT approach and its estimated fracture toughness for the Zr-4 cladding structure with segment alumina inserts are presented in this report. The SNTT test results indicate that SNTT method is a reliable test approach with good repeatability in applying to Zr-4 cladding material. For the medium and the short crack length, the estimate J Q upon fracture for the baseline Zr-4 cladding is at 285.7 lb/in with 2-sigma uncertainty of 18.59 lb/in., and the associated K Q is at 61.4 Ksi√in. For a long crack length, the crack initialization orientation is apparently deviated from the principle tensile stress contour – 45° spiral crack front (i.e., Mode I + Mode III failure mechanism); and the estimated J Q is at 108 lb/in., the associated K Q is at 37.7 Ksi√in. The estimate Mode I J Q (along the 45° pitch principle tensile stress contour) is at 200 lb/in.; which indicate a significant reduction in the fracture toughness in a tubing structure under mixed mode loading condition, Mode I + Mode III, compared to that of Mode I (tensile stress dominated) alone.

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

Enhanced Dynamic Equivalent Identification Method of Large-Scale Power Systems Using Multiple Event

The increasing complexity of the interconnected power system makes high-fidelity dynamic simulation models computationally more intensive. To improve computation efficiency, model reduction techniques have been investigated to only preserve the dynamics in a limited area of interest (study area), while deriving equivalent representation for the external area. For this purpose, a measurement-based reduction approach using system identification techniques has been previously proposed. In that case, external areas are represented by dynamic equivalent loads using transfer function estimation. In this paper, in order to enhance the accuracy of the reduced model in preserving dynamics of the study area, multiple “grid” events of different types and at different locations are considered for identifying the parameters of the equivalent loads. Case studies are carried out in the NPCC 140-bus system. Comparison are made between multiple events training and single event training, highlighting the advantages of the proposed method in providing a better representation of the grid dynamics under different operating conditions.

Xue, Yaosuo↗
DOE OSTI2019

Self-consistent redshift estimation using correlation functions without a spectroscopic reference sample

ABSTRACT We present a new method to estimate redshift distributions and galaxy-dark matter bias parameters using correlation functions in a fully data driven and self-consistent manner. Unlike other machine learning, template, or correlation redshift methods, this approach does not require a reference sample with known redshifts. By measuring the projected cross- and auto-correlations of different galaxy sub-samples, e.g. as chosen by simple cells in colour–magnitude space, we are able to estimate the galaxy-dark matter bias model parameters, and the shape of the redshift distributions of each sub-sample. This method fully marginalizes over a flexible parametrization of the redshift distribution and galaxy-dark matter bias parameters of sub-samples of galaxies, and thus provides a general Bayesian framework to incorporate redshift uncertainty into the cosmological analysis in a data-driven, consistent, and reproducible manner. This result is improved by an order of magnitude by including cross-correlations with the cosmic microwave background and with galaxy–galaxy lensing. We showcase how this method could be applied to real galaxies. By using idealized data vectors, in which all galaxy-dark matter model parameters and redshift distributions are known, this method is demonstrated to recover unbiased estimates on important quantities, such as the offset Δz between the mean of the true and estimated redshift distribution and the 68 per cent, 95 per cent, and 99.5 per cent widths of the redshift distribution to an accuracy required by current and future surveys.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

On the impact of the galaxy window function on cosmological parameter estimation

One important source of systematics in galaxy redshift surveys comes from the estimation of the galaxy window function. Up until now, the impact of the uncertainty in estimating the galaxy window function on parameter inference has not been properly studied. In this paper, we show that the uncertainty and the bias in estimating the galaxy window function will be salient for ongoing and next-generation galaxy surveys using a simulation-based approach. With a specific case study of cross-correlating emission-line galaxies from the DESI Legacy Imaging Surveys and the Planck cosmic microwave background lensing map, we show that neural network-based regression approaches to modelling the window function are superior in comparison to linear regression-based models. Here, we additionally show that the definition of the galaxy overdensity estimator can impact the overall signal-to-noise of observed power spectra. Finally, we show that the additive biases coming from the window functions can significantly bias the modes of the inferred parameters and also degrade their precision. Thus, a careful understanding of the window functions will be essential to conduct cosmological experiments.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Robust analytic continuation of Green's functions via projection, pole estimation, and semidefinite relaxation

Green's functions of fermions are described by matrix-valued Herglotz-Nevanlinna functions. Since analytic continuation is fundamentally an ill-posed problem, the causal space described by the matrix-valued Herglotz-Nevanlinna structure can be instrumental in improving the accuracy and in enhancing the robustness with respect to noise. We demonstrate a three-pronged procedure for robust analytic continuation called PES: (1) projection of data to the causal space; (2) estimation of pole locations; and (3) semidefinite relaxation within the causal space. We compare the performance of PES with the recently developed Nevanlinna and Carathéodory continuation methods and find that PES is more robust in the presence of noise and does not require the usage of extended precision arithmetics. We also demonstrate that a causal projection improves the performance of the Nevanlinna and Carathéodory methods. The PES method is generalized to bosonic response functions, for which the Nevanlinna and Carathéodory continuation methods have not yet been developed. Furthermore, it is particularly useful for studying spectra with sharp features, as they occur in the study of molecules and band structures in solids.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2022

Maximum likelihood estimation of difference scaling functions for suprathreshold judgments

Maximum likelihood estimation (MLE) has been used to produce perceptual scales from binary judgments of triads and quadruples. This method relies on Thurstone’s theory of a stochastic perceptual process where the perceived difference of two stimuli is the difference in their perceived strengths. It is possible that the perception of a suprathreshold difference is overestimated when adding smaller differences, a phenomenon referred to as diminishing returns. The current approach to construct a perceptual scale using MLE does not account for this phenomenon. We present a way to model the perception of differences using MLE and Thurstone’s theory, adapted to allow the possibility of diminishing returns. This method is validated using Monte Carlo simulated responses to experimental triads and can correctly model diminishing returns, the absence of diminishing returns, and the opposite of diminishing returns both in the cases when a perceptual scale is known and when the true perceived strengths of the stimuli are unknown. Additionally, this method was applied to empirical data sets to determine its feasibility in investigations of perception. Ultimately, it was found that this analysis allows for more accurate modeling of suprathreshold difference judgments, a more complete understanding of the perceptual processes underlying comparisons, and the evaluation of Thurstone’s theory of difference judgments.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2020

The Effects of Atmospheric Models on the Estimation of Infrasonic Source Functions at the Source Physics Experiment

We invert infrasound signals for an equivalent seismoacoustic source function using different atmospheric models to produce the necessary Green’s functions. The infrasound signals were produced by a series of underground chemical explosions as part of the Source Physics Experiment (SPE). In a previous study, we inverted the infrasound data using so-called predictive atmospheric models, which were based on historic, regional-scaled, publicly available weather observations interpolated onto a 3D grid. For the work presented here, we invert the same infrasound data, but using atmospheric models based on weather data collected in a time window that includes the approximate time of the explosion experiments, which we term postdictive models. We build two versions of the postdictive models for each SPE event: one that is based solely on the regional scaled observations, and one that is based on both regional scaled observations combined with on-site observations obtained by a weather sonde released at the time of the SPE. We then invert the observed data set three times, once for each atmospheric model type. We find that the estimated seismoacoustic source functions are relatively similar in waveform shape regardless of which atmospheric model that we used to construct the Green’s functions. However, we find that the amplitude of the estimated source functions is systematically dependent on the atmospheric model type: using the predictive atmospheric models to invert the data generally yields estimated source functions that are larger in amplitude than those estimated using the postdictive models.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2023

A machine learning approach targeting parameter estimation for plant functional type coexistence modeling using ELM-FATES (v2.0)

Tropical forest dynamics play a crucial role in the global carbon, water, and energy cycles. However, realistically simulating the dynamics of competition and coexistence between different plant functional types (PFTs) in tropical forests remains a significant challenge. This study aims to improve the modeling of PFT coexistence in the Functionally Assembled Terrestrial Ecosystem Simulator (FATES), a vegetation demography model implemented in the Energy Exascale Earth System Model (E3SM) land model (ELM), ELM-FATES. Specifically, we explore (1) whether plant trait relationships established from field measurements can constrain ELM-FATES simulations and (2) whether machine learning (ML)-based surrogate models can emulate the complex ELM-FATES model and optimize parameter selections to improve PFT coexistence modeling. We conducted three ensembles of ELM-FATES experiments at a tropical forest site near Manaus, Brazil. By comparing the ensemble experiments without (Exp-CTR) and with (Exp-OBS) consideration of observed trait relationships, we found that accounting for these relationships slightly improves the simulations of water, energy, and carbon variables when compared to observations but degrades the simulation of PFT coexistence. Using ML-based surrogate models trained on Exp-CTR, we optimized the trait parameters in ELM-FATES and conducted another ensemble of experiments (Exp-ML) with these optimized parameters. The proportion of PFT coexistence experiments significantly increased from 21 % in Exp-CTR to 73 % in Exp-ML. After filtering the experiments that allow for PFT coexistence to agree with observations (within 15 % tolerance), 33 % of the Exp-ML experiments were retained, which is a significant improvement compared to the 1.4 % in Exp-CTR. Exp-ML also accurately reproduces the annual means and seasonal variations in water, energy, and carbon fluxes and the field inventory of aboveground biomass. This study represents a reproducible method that utilizes machine learning to identify parameter values that improve model fidelity against observations and PFT coexistence in vegetation demography models for diverse ecosystems. Our study also suggests the need for new mechanisms to enhance the robust simulation of coexisting plants in ELM-FATES and has significant implications for modeling the response and feedbacks of ecosystem dynamics to climate change.

54 ENVIRONMENTAL SCIENCES↗
DOE OSTI2025

A note on the reliability of goal-oriented error estimates for Galerkin finite element methods with nonlinear functionals

Here, we consider estimating the discretization error in a nonlinear functional J (u) in the setting of an abstract variational problem: find u ϵ $\mathscr{V}$ such that B (u, φ) = L (φ) ∀φ ϵ $\mathscr{V}$, as approximated by a Galerkin finite element method. Here, $\mathscr{V}$ is a Hilbert space, B (. , .) is a bilinear form, and L (∙) is a linear functional. We consider well-known error estimates η of the form J (u) - J (u h ) ≈ η = L (z) - B (u h , z), where u h denotes a finite element approximation to u, and z denotes the solution to an auxiliary adjoint variational problem. We show that there exist nonlinear functionals for which error estimates of this form are not reliable, even in the presence of an exact adjoint solution z. An estimate η is said to be reliable if there exists a constant C ϵ $\mathbb{R}$ >0 independent of u h such that |J (u) - J (u h )| ≤ C|η|. We present several example pairs of bilinear forms and nonlinear functionals where reliability of η is not achieved.

A posteriori↗
DOE OSTI2022

Statistical Estimation of Strain Using Spatial Correlation Functions

Ex-situ estimation of strains from deformed micrographs is not possible as there are no persistent features which can be tracked. Two point spatial statistics enable the rigorous quantification of spatial patterns in heterogeneous media. In this paper, we propose a novel method for estimating strains directly from dissimilar micrographs using a continuum mechanics approach. Rather than operating directly on images from sequential frames, as is done in digital image correlation, we operate on different microstructure realizations. This is made possible by comparing the spatial autocorrelation maps of deformed and undeformed micrographs rather than direct comparison of images. Additionally, a Bayesian framework is proposed for quantifying uncertainty. We first illustrate the efficacy of this method on speckle pattern images from digital image correlation experiments. Then, we demonstrate that the method is capable of operating on dissimilar micrographs using deformed synthetic binary microstructures. Finally, we present a case study on polycrystalline additively manufactured 316L deformed via tension. The proposed method works well and we discuss implications and limitations of the presented work.

36 MATERIALS SCIENCE↗
DOE OSTI2024

An Empirical Quantile Estimation Approach for Chance-Constrained Nonlinear Optimization Problems

We investigate an empirical quantile estimation approach to solve chance-constrained nonlinear optimization problems. Our approach is based on the reformulation of the chance constraint as an equivalent quantile constraint to provide stronger signals on the gradient. In this approach, the value of the quantile function is estimated empirically from samples drawn from the random parameters, and the gradient of the quantile function is estimated via a finite-difference approximation on top of the quantile-function-value estimation. We establish a convergence theory of this approach within the framework of an augmented Lagrangian method for solving general nonlinear constrained optimization problems. The foundation of the convergence analysis is a concentration property of the empirical quantile process, and the analysis is divided based on whether or not the quantile function is differentiable. In contrast to the sampling-and-smoothing approach used in the literature, the method developed in this paper does not involve any smoothing function and hence the quantile-function gradient approximation is easier to implement and there are less accuracy-control parameters to tune. Furthermore, we demonstrate the effectiveness of this approach and compare it with a smoothing method for the quantile-gradient estimation. Numerical investigation shows that the two approaches are competitive for certain problem instances.

Applied Probability↗
DOE OSTI2019

Estimation of split renal function using different volumetric methods: inter- and intraindividual comparison between MRI and CT

Purpose: This study aims to determine whether contrast-enhanced (CE)-magnetic resonance imaging (MRI) is comparable to CE-computed tomography (CT) for estimation of split renal function (SRF). For this purpose, two different kidney volumetry methods, the renal cortex volumetry (RCV) and modified ellipsoid volume (MELV), are compared for both acquisition types (CT vs. MRI) with regard to accuracy and reliability, subsequently referred to as RCV{sub CT}/RCV{sub MRI} and MELV{sub CT}/MELV{sub MRI}. Methods: This retrospective study included 29 patients (18 men and 11 women; mean age 62.8 ± 12.4 years) who underwent CE-MRI and CE-CT of the abdomen within a period of 3 months. Two independent readers (R1/R2) performed RCV and MELV in all datasets with corresponding semiautomated software tools. RCV was performed with datasets in the arterial phase and MELV in the venous phase. Statistics were calculated using one-way ANOVA, two-tailed Student’s t test, Pearson´s correlation, and Bland–Altman plots with p ≤ 0.05 being considered statistically significant. Results: In all datasets, SRF was almost identical for both volumetry methods with a mean difference of < 1%. Bland–Altman analysis comparing RCV in CT and MRI showed very good agreement for R1/R2. Interreader agreement was strong for RCV{sub CT} and good for RCV{sub MRI} (r = 0.89; r = 0.69). MELV{sub CT/MRI} interreader agreement was only moderate (r = 0.54; r = 0.50) with a high range of values. Intrareader agreement was excellent for all measurements, except MELV{sub MRI} which showed a high mean bias and range of values (RCV{sub CT}: r = 0.93, RCV{sub MRI}: r = 0.98, MELV{sub CT}: r = 0.89, MELV{sub MRI}: r = 0.54). Conclusion: Renal volumetric estimates of SRF are almost as accurate and reliable with CE-MRI as with CE-CT using RCV method. In distinction, the calculation of SRF using MELV was inferior to RCV with respect to accuracy and reliability. Thus, RCV method is recommended to estimate SRF, primarily using CT datasets. However, RCV with MRI datasets for kidney volumetry allows for comparable accuracy and reliability while sparing patients and healthy donors of unnecessary radiation exposure.

62 RADIOLOGY AND NUCLEAR MEDICINE↗
DOE OSTI2018

Generators for Spaces of Entire Functions with a System of Weighted Estimates

We consider spaces of entire functions with systems of weighted estimates. The case of twoterm weight sequences consisting of radial and nonradial components is studied. Under some assumptions on the weight sequence, we obtain a complete description of generators in these spaces. We apply this result to the problem of normal solvability of systems of convolution equations in the Roumieu spaces of ultradifferentiable functions and, as a particular case, in Gevrey classes.

97 MATHEMATICS AND COMPUTING↗
DOE OSTI2023

Tangent functional connectomes uncover more unique phenotypic traits

Functional connectomes (FCs) containing pairwise estimations of functional couplings between pairs of brain regions are commonly represented by correlation matrices. As symmetric positive definite matrices, FCs can be transformed via tangent space projections, resulting into tangent-FCs. Tangent-FCs have led to more accurate models predicting brain conditions or aging. Motivated by the fact that tangent-FCs seem to be better biomarkers than FCs, we hypothesized that tangent-FCs have also a higher fingerprint. We explored the effects of six factors: fMRI condition, scan length, parcellation granularity, reference matrix, main-diagonal regularization, and distance metric. Our results showed that identification rates are systematically higher when using tangent-FCs across the “fingerprint gradient” (here including test-retest, monozygotic and dizygotic twins). Highest identification rates were achieved when minimally (0.01) regularizing FCs while performing tangent space projection using Riemann reference matrix and using correlation distance to compare the resulting tangent-FCs. Such configuration was validated in a second dataset (resting-state).

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2022

Meta-Study of Particulate Detection Losses on Radioactive Air Sample Filters

Several mathematical relationships between air sample filter mass loading and the correlated analytical self-absorption factor were developed using data from other published research in this meta-study. Gross-alpha and -beta applications are addressed for this research. As filter media becomes loaded with particulate matter, there is potential for measurement losses due to self-absorption by mass loading. Components contributing to absorption include particulate dust, radioactive particulates, and filter material. Standards indicate a correction factor should be used when the penetration of radioactive material into the collection media or self-absorption of radiation by the material collected would reduce the detection rate by more than 5%. Previously, losses due to self-absorption have been reported up to 100% over a range up to ~10 mg∙cm -2 mass loading. These absorption losses then can be used to determine a correction factor for sample results. For low mass loadings (e.g., ≤0.1 mg∙cm -2 ) corrections factors in the 0.85 - 1 range have been recommended and used, while at higher mass loadings nearer to 10 mg∙cm -2 correction factors closer to 0 (representing near 100% losses) are used. Based on data from published studies, the different methods for relating percent loss due to self-absorption to mass loading include linear, exponential, quadratic, and trinomial derived functions. Where applicable, both forced zero and non-forced zero results were evaluated. From the derived functions evaluated, the trinomial function provided the best fit. Once the sample filter mass loading is known, the trinomial function can be applied to estimate losses and the corresponding self-absorption factor. When applied to routine operating conditions for radiological facility stacks monitored at the Pacific Northwest National Laboratory for an average sample filter mass loading of 0.09 ± 0.12 (2σ) mg∙cm -2 (excluding negative values and outliers) and a range from 0 - 0.24 mg∙cm -2 , the estimated trinomial function nominal self-absorption losses are less than 5% at 0.09 mg∙cm -2 and less than 10% at 0.24 mg∙cm -2 . The trinomial function is one method that may be used to adjust the activity results of an air sample when the sample-specific mass loading is determined. The application of no correction factor when the ANSI/HPS N13.1-2021 guidance of a 5% threshold for loss is not reached with typical stack sample mass loadings may be reasonable in high-efficiency particulate air filtered systems. For simplicity, it would be conservative in assigning the self-absorption correction factor at the 5% threshold (i.e., 0.95) for general uses but in cases of heavy mass loading to calculate the factor.

air sampling↗
DOE OSTI2025

Neural Posterior Estimation for Cataloging Astronomical Images with Spatially Varying Backgrounds and Point Spread Functions

Neural posterior estimation (NPE), a type of amortized variational inference, is a computationally efficient means of constructing probabilistic catalogs of light sources from astronomical images. To date, NPE has not been used to perform inference in models with spatially varying covariates. However, ground-based astronomical images exhibit spatially varying sky backgrounds and point spread functions (PSFs), and accounting for this variation is essential for constructing accurate catalogs of imaged light sources. In this work, we introduce a novel NPE-based cataloging method that trains an inference network with semisynthetic astronomical images generated using PSFs and backgrounds sampled from the Sloan Digital Sky Survey. In experiments with semisynthetic images, we evaluate the method on key cataloging tasks: light source detection, star/galaxy separation, and flux measurement. A “generalist” inference network—trained with diverse PSFs and backgrounds—performs as well as a “specialist” network even when both are evaluated on the specialist’s particular PSF/background combination. This result suggests that a single NPE network can generalize across spatial variations, eliminating the need for retraining on each observational condition.

astronomy image processing↗