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 91 records · Page 5

Development on electromagnetic impedance function modeling and its estimation

Today the Electromagnetic methods such as magnetotellurics (MT) and controlled sources audio MT (CSAMT) is used in a broad variety of applications. Its usefulness in poor seismic areas and its negligible environmental impact are integral parts of effective exploration at minimum cost. As exploration was forced into more difficult areas, the importance of MT and CSAMT, in conjunction with other techniques, has tended to grow continuously. However, there are obviously important and difficult problems remaining to be solved concerning our ability to collect process and interpret MT as well as CSAMT in complex 3D structural environments. This talk aim at reviewing and discussing the recent development on MT as well as CSAMT impedance functions modeling, and also some improvements on estimation procedures for the corresponding impedance functions. In MT impedance modeling, research efforts focus on developing numerical method for computing the impedance functions of three dimensionally (3-D) earth resistivity models. On that reason, 3-D finite elements numerical modeling for the impedances is developed based on edge element method. Whereas, in the CSAMT case, the efforts were focused to accomplish the non-plane wave problem in the corresponding impedance functions. Concerning estimation of MT and CSAMT impedance functions, researches were focused on improving quality of the estimates. On that objective, non-linear regression approach based on the robust M-estimators and the Hilbert transform operating on the causal transfer functions, were used to dealing with outliers (abnormal data) which are frequently superimposed on a normal ambient MT as well as CSAMT noise fields. As validated, the proposed MT impedance modeling method gives acceptable results for standard three dimensional resistivity models. Whilst, the full solution based modeling that accommodate the non-plane wave effect for CSAMT impedances is applied for all measurement zones, including near-, transition-as well as the far-field zones, and consequently the plane wave correction is no longer needed for the impedances. In the resulting robust impedance estimates, outlier contamination is removed and the self consistency between the real and imaginary parts of the impedance estimates is guaranteed. Using synthetic and real MT data, it is shown that the proposed robust estimation methods always yield impedance estimates which are better than the conventional least square (LS) estimation, even under condition of severe noise contamination. A recent development on the constrained robust CSAMT impedance estimation is also discussed. By using synthetic CSAMT data it is demonstrated that the proposed methods can produce usable CSAMT transfer functions for all measurement zones.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

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↗

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↗

Safety and Efficacy of Percutaneous Biliary Covered Stent Placement in Patients with Malignant Biliary Hilar Obstruction; Correlation with Liver Function

PurposeTo estimate the safety and efficacy of percutaneous ePTFE-covered biliary stent placement and the relationship between underlying liver function and stent patency in patients with malignant hilar obstruction.Materials and MethodsFrom March 2012 to June 2015, 41 patients [22 females, 19 males; mean age 69.8 (range 34–94) years] with malignant biliary obstruction underwent percutaneous biliary stent placement (31 patients with unilateral, 10 patients with bilateral side-by-side). Cumulative patient survival and stent patency rate curves were derived using the Kaplan–Meier method. A Cox model was used to explore the relationship between liver function and patient survival, and also biliary stent patency. Pearson correlation coefficient was used to analyze the relationship between patient survival and stent patency.ResultsTechnical success rate was 100 % and clinical success rate was 95 %. During follow-up, four complications occurred (two bilomas and two cases of acute cholecystitis) and were treated successfully with percutaneous drainage. No other complication occurred. Mean serum bilirubin level was 11.34 ± 7.35 mg/dL before drainage and 5.00 ± 4.83 mg/dL 2 weeks after stent placement. The median patent survival duration was 147 days (95 % CI, 69.6–224.4 days). The median stent patency duration was 101 days (95 % CI, 70.0–132.0 days). The cumulative stent patency rates at 1, 3, 6, and 12 months were 97, 57.6, 30.3, and 17.0 %, respectively. Child–Pugh score was correlated significantly with patient survival (P = 0.011) and stent patency (P = 0.007). MELD score was correlated significantly with stent patency (P = 0.044). There was a correlation between patient survival and stent patency (r = 0.778, P < 0.001).ConclusionPercutaneous placement of ePTFE-covered biliary stent was a safe and an effective method for malignant biliary obstruction. Underlying liver function seemed to be one of the important factors affecting patient survival and stent patency, and stent patency showed statistically significant correlation with patient survival.

62 RADIOLOGY AND NUCLEAR MEDICINE↗

Error estimates and specification parameters for functional renormalization

We present a strategy for estimating the error of truncated functional flow equations. While the basic functional renormalization group equation is exact, approximated solutions by means of truncations do not only depend on the choice of the retained information, but also on the precise definition of the truncation. Therefore, results depend on specification parameters that can be used to quantify the error of a given truncation. We demonstrate this for the BCS–BEC crossover in ultracold atoms. Within a simple truncation the precise definition of the frequency dependence of the truncated propagator affects the results, indicating a shortcoming of the choice of a frequency independent cutoff function.

74 ATOMIC AND MOLECULAR PHYSICS↗

Improved estimates of coordinate error for molecular replacement

A function for estimating the effective root-mean-square deviation in coordinates between two proteins has been developed that depends on both the sequence identity and the size of the protein and is optimized for use with molecular replacement in Phaser. A top peak translation-function Z-score of over 8 is found to be a reliable metric of when molecular replacement has succeeded. The estimate of the root-mean-square deviation (r.m.s.d.) in coordinates between the model and the target is an essential parameter for calibrating likelihood functions for molecular replacement (MR). Good estimates of the r.m.s.d. lead to good estimates of the variance term in the likelihood functions, which increases signal to noise and hence success rates in the MR search. Phaser has hitherto used an estimate of the r.m.s.d. that only depends on the sequence identity between the model and target and which was not optimized for the MR likelihood functions. Variance-refinement functionality was added to Phaser to enable determination of the effective r.m.s.d. that optimized the log-likelihood gain (LLG) for a correct MR solution. Variance refinement was subsequently performed on a database of over 21 000 MR problems that sampled a range of sequence identities, protein sizes and protein fold classes. Success was monitored using the translation-function Z-score (TFZ), where a TFZ of 8 or over for the top peak was found to be a reliable indicator that MR had succeeded for these cases with one molecule in the asymmetric unit. Good estimates of the r.m.s.d. are correlated with the sequence identity and the protein size. A new estimate of the r.m.s.d. that uses these two parameters in a function optimized to fit the mean of the refined variance is implemented in Phaser and improves MR outcomes. Perturbing the initial estimate of the r.m.s.d. from the mean of the distribution in steps of standard deviations of the distribution further increases MR success rates.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗

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↗

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↗

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↗

A tight distance-dependent estimator for screening three-center Coulomb integrals over Gaussian basis functions

A new estimator for three-center two-particle Coulomb integrals is presented. Our estimator is exact for some classes of integrals and is much more efficient than the standard Schwartz counterpart due to the proper account of distance decay. Although it is not a rigorous upper bound, the maximum degree of underestimation can be controlled by two adjustable parameters. We also give numerical evidence of the excellent tightness of the estimator. The use of the estimator will lead to increased efficiency in reduced-scaling one- and many-body electronic structure theories.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗

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↗

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↗

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↗

An Empirical Orthogonal Function-Based Algorithm for Estimating Terrestrial Latent Heat Flux from Eddy Covariance, Meteorological and Satellite Observations

Accurate estimation of latent heat flux (LE) based on remote sensing data is critical in characterizing terrestrial ecosystems and modeling land surface processes. Many LE products were released during the past few decades, but their quality might not meet the requirements in terms of data consistency and estimation accuracy. Merging multiple algorithms could be an effective way to improve the quality of existing LE products. In this paper, we present a data integration method based on modified empirical orthogonal function (EOF) analysis to integrate the Moderate Resolution Imaging Spectroradiometer (MODIS) LE product (MOD16) and the Priestley-Taylor LE algorithm of Jet Propulsion Laboratory (PT-JPL) estimate. Twenty-two eddy covariance (EC) sites with LE observation were chosen to evaluate our algorithm, showing that the proposed EOF fusion method was capable of integrating the two satellite data sets with improved consistency and reduced uncertainties. Further efforts were needed to evaluate and improve the proposed algorithm at larger spatial scales and time periods, and over different land cover types.

54 ENVIRONMENTAL SCIENCES↗

RECONSTRUCTING THE SHAPE OF THE CORRELATION FUNCTION

We develop an estimator for the correlation function which, in the ensemble average, returns the shape of the correlation function, even for signals that have significant correlations on the scale of the survey region. Our estimator is general and works in any number of dimensions. We develop versions of the estimator for both diffuse and discrete signals. As an application, we apply them to Monte Carlo simulations of X-ray background measurements. These include a realistic, spatially inhomogeneous population of spurious detector events. We discuss applying the estimator to the averaging of correlation functions evaluated on several small fields, and to other cosmological applications.

79 ASTRONOMY AND ASTROPHYSICS↗

A probability density function model describing height estimation uncertainty due to image pixel intensity noise in digital fringe projection measurements

Digital fringe projection is a surface-profiling technique used for highly accurate non-contact measurements. As with any measurement technique, a variety of sources degrade to the measurement accuracy of the method. Here, this paper presents an analytically-derived probability density function that explicitly models the surface height measurement error due to inevitable phase measurement error, and it includes the specific case of pixel noise inducing the phase measurement error that ultimately leads to the height estimation error. The accuracy of the model was validated through Monte-Carlo simulations of resultant height distributions subject to arbitrarily correlated pixel intensity noise and experimental digital fringe projection measurements where the pixel-by-pixel height uncertainty estimations were compared to the predictions of the derived model.

42 ENGINEERING↗

Series of semi-Markov processes to model infrastructure resilience under multihazards

Civil infrastructure systems are subjected to multiple hazards, including natural and anthropogenic, that disrupt their function or the level of service offered. Estimating the function recovery of these systems (or how soon normalcy of operations will be restored) when subjected to repeated hazard events by considering the inter-event dependencies is an important problem in multihazard infrastructure resilience. However, this problem has been less addressed in the field. This paper proposes a series of semi-Markov processes model to capture the inter-event dependencies in infrastructure recovery when subjected to successive hazard events. Recovery after each new hazard event is represented by a unique semi-Markov process that models the reduced recovery rates and the increased recovery times caused by the system’s incomplete recovery from the preceding event. Two novel formulations of the inter-event dependency modeling, namely Maximal Effects Dependency (considers the worst impact of two successive hazard events) and Cumulative Effects Dependency (considers the aggregated impacts of two successive hazard events), are proposed and discussed. The model is demonstrated by considering the following applications: Three-state system subjected to deterministic and random occurrences of identical hazard events; and Multihazard resilience of a building in Charleston, SC, considering earthquake and hurricane hazards. Results indicate that considering inter-event dependencies in recovery modeling can lead to lesser-predicted resilience, thereby affecting resilience-based decision-making.

42 ENGINEERING↗