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 649 records · Page 36

A stochastic constrained optimization technique and its application to detector array processing.

A stochastic projected gradient algorithm is proposed which can be used for finding a constrained optimum point for a concave or convex objective function subject to nonlinear constraints which form a connected region even when only a noisy estimate of the objective function is available. For a constraint described by a single linear equation, convergence to the constrained optimum value is proved, and the rate of convergence of the algorithm to the constrained optimum value is determined. The algorithm is applied to the nonlinear problem of obtaining automatically an array of detectors which forms a beam in a desired direction in space in the presence of interfering noise so as to maximize the SNR subject to a constraint on the super-gain ratio.

Winkler, L. P.↗

Retrieval of monthly rainfall indices from microwave radiometric measurements using probability distribution functions

An algorithm for the estimation of monthly rain totals for 5 deg cells over the ocean from histograms of SSM/I brightness temperatures has been developed. There are three novel features to this algorithm. First, it uses knowledge of the form of the rainfall intensity probability density function to augment the measurements. Second, a linear combination of the 19.35 and 22.235 GHz channels has been employed to reduce the impact of variability of water vapor. Third, an objective technique has been developed to estimate the rain layer thickness from the 19.35- and 22.235-GHz brightness temperature histograms. Comparison with climatologies and the GATE radar observations suggest that the estimates are reasonable in spite of not having a beam-filling correction. By-products of the retrievals indicate that the SSM/I instrument noise level and calibration stability are quite good.

Wilheit, Thomas T.↗

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↗

The gust-front detection and wind-shift algorithms for the Terminal Doppler Weather Radar system

The Federal Aviation Administration's (FAA) Terminal Doppler Weather Radar (TDWR) system was primarily designed to address the operational needs of pilots in the avoidance of low-altitude wind shears upon takeoff and landing at airports. One of the primary methods of wind-shear detection for the TDWR system is the gust-front detection algorithm. The algorithm is designed to detect gust fronts that produce a wind-shear hazard and/or sustained wind shifts. It serves the hazard warning function by providing an estimate of the wind-speed gain for aircraft penetrating the gust front. The gust-front detection and wind-shift algorithms together serve a planning function by providing forecasted gust-front locations and estimates of the horizontal wind vector behind the front, respectively. This information is used by air traffic managers to determine arrival and departure runway configurations and aircraft movements to minimize the impact of wind shifts on airport capacity. This paper describes the gust-front detection and wind-shift algorithms to be fielded in the initial TDWR systems. Results of a quantitative performance evaluation using Doppler radar data collected during TDWR operational demonstrations at the Denver, Kansas City, and Orlando airports are presented. The algorithms were found to be operationally useful by the FAA airport controllers and supervisors.

Hermes, Laurie G.↗

A model for estimating time-variant rainfall infiltration as a function of antecedent surface moisture and hydrologic soil type

Recent research indicates that the use of remote sensing techniques for the measurement of near surface soil moisture could be practical in the not too distant future. Other research shows that infiltration rates, especially for average or frequent rainfall events, are extremely sensitive to the proper definition and consideration of the role of the soil moisture at the beginning of the rainfall. Thus, it is important that an easy to use, but theoretically sound, rainfall infiltration model be available if the anticipated remotely sensed soil moisture data is to be optimally utilized for hydrologic simulation. A series of numerical experiments with the Richards' equation for an array of conditions anticipated in watershed hydrology were used to develop functional relationships that describe temporal infiltration rates as a function of soil type and initial moisture conditions.

Wilkening, H. A.↗

Least-squares sequential parameter and state estimation for large space structures

This paper presents the formulation of simultaneous state and parameter estimation problems for flexible structures in terms of least-squares minimization problems. The approach combines an on-line order determination algorithm, with least-squares algorithms for finding estimates of modal approximation functions, modal amplitudes, and modal parameters. The approach combines previous results on separable nonlinear least squares estimation with a regression analysis formulation of the state estimation problem. The technique makes use of sequential Householder transformations. This allows for sequential accumulation of matrices required during the identification process. The technique is used to identify the modal prameters of a flexible beam.

Thau, F. E.↗

Estimation of finite mixtures using the empirical characteristic function

A problem which occurs in analyzing LANDSAT scenes is the problem of separating the components of a finite mixture of several distinct probability distributions. A review of the literature indicates this is a problem which occurs in many disciplines, such as engineering, biology, physiology and economics. Many approaches to this problem have appeared in the literature; however, most are very restrictive in their assumptions or have met with only a limited degree of success when applied to realistic situations. A proceudre is investigated with combines the k-L procedure of (Feurverger and McDunnough, 1981) with the MAICE procedure of (Akaike, 1974). The feasibility of this approach is being investigated numerically via the development of a computer software package enabling a simulation study and comparison with other procedures.

Anderson, C.↗

Satellite Material Type and Phase Function Determination in Support of Orbital Debris Size Estimation

In performing debris surveys of deep-space orbital regions, the considerable volume of the area to be surveyed and the increased orbital altitude suggest optical telescopes as the most efficient survey instruments; but to proceed this way, methodologies for debris object size estimation using only optical tracking and photometric information are needed. Basic photometry theory indicates that size estimation should be possible if satellite albedo and shape are known. One method for estimating albedo is to try to determine the object's material type photometrically, as one can determine the albedos of common satellite materials in the laboratory. Examination of laboratory filter photometry (using Johnson BVRI filters) on a set of satellite material samples indicates that most material types can be separated at the 1-sigma level via B-R versus R-I color differences with a relatively small amount of required resampling, and objects that remain ambiguous can be resolved by B-R versus B-V color differences and solar radiation pressure differences. To estimate shape, a technique advanced by Hall et al. [1], based on phase-brightness density curves and not requiring any a priori knowledge of attitude, has been modified slightly to try to make it more resistant to the specular characteristics of different materials and to reduce the number of samples necessary to make robust shape determinations. Working from a gallery of idealized debris shapes, the modified technique identifies most shapes within this gallery correctly, also with a relatively small amount of resampling. These results are, of course, based on relatively small laboratory investigations and simulated data, and expanded laboratory experimentation and further investigation with in situ survey measurements will be required in order to assess their actual efficacy under survey conditions; but these techniques show sufficient promise to justify this next level of analysis.

Hejduk, M. D.↗

Model-based prognostics for batteries which estimates useful life and uses a probability density function

This invention develops a mathematical model to describe battery behavior during individual discharge cycles as well as over its cycle life. The basis for the form of the model has been linked to the internal processes of the battery and validated using experimental data. Effects of temperature and load current have also been incorporated into the model. Subsequently, the model has been used in a Particle Filtering framework to make predictions of remaining useful life for individual discharge cycles as well as for cycle life. The prediction performance was found to be satisfactory as measured by performance metrics customized for prognostics for a sample case. The work presented here provides initial steps towards a comprehensive health management solution for energy storage devices.

Saha, Bhaskar↗

New Developments in Quaternion Estimation from Vector Observations

This paper contains a critical comparison of estimators minimizing Wahba's loss function. Some new results are presented for the QUaternion ESTimator (QUEST) and EStimators of the Optimal Quaternion (ESOQ and ESOQ2) to avoid the computational burden of sequential rotations in these algorithms. None of these methods is as robust in principle as Davenport's q method or the Singular Value Decomposition (SVD) method, which are significantly slower. Robustness is only an issue for measurements with widely differing accuracies, so the fastest estimators, the modified ESOQ and ESOQ2, are well suited to sensors that track multiple stars with comparable accuracies. More robust forms of ESOQ and ESOQ2 are developed that are intermediate in speed.

Markley, F. Landis↗

Hybrid diversity method utilizing adaptive diversity function for recovering unknown aberrations in an optical system

A method of recovering unknown aberrations in an optical system includes collecting intensity data produced by the optical system, generating an initial estimate of a phase of the optical system, iteratively performing a phase retrieval on the intensity data to generate a phase estimate using an initial diversity function corresponding to the intensity data, generating a phase map from the phase retrieval phase estimate, decomposing the phase map to generate a decomposition vector, generating an updated diversity function by combining the initial diversity function with the decomposition vector, generating an updated estimate of the phase of the optical system by removing the initial diversity function from the phase map. The method may further include repeating the process beginning with iteratively performing a phase retrieval on the intensity data using the updated estimate of the phase of the optical system in place of the initial estimate of the phase of the optical system, and using the updated diversity function in place of the initial diversity function, until a predetermined convergence is achieved.

Dean, Bruce H.↗

Kalman filter estimation of attitude and gyro bias with the QUEST observation model

The loss function for the Wahba attitude estimation problem employs unit vector observations with scalar weights. It is usually associated with the QUEST observation model, where the actual sensor noise is assumed the same for all components of the observed vector, regardless of the position in the field of view. The QUEST model has the great advantage of being sensor-independent except for the scalar parameter that characterizes the sensor errors. Although efficient algorithms for solving the Wahba problem exist, extending these algorithms to estimate gyro biases or sensor alignments has had mixed success. However, it is straightforward to estimate bias and alignment parameters with a Kalman filter. This paper investigates the use of an extended Kalman filter for the attitude and gyro bias that incorporates the QUEST observation model, to be referred to as the Unit Vector Filter (UVF). The UVF results are compared with those from a more conventional filter, the Real-Time Sequential Filter (RTSF), for which the residual is the two-dimensional projection of the unit vector onto the plane perpendicular to the sensor boresight. The RTSF is similar in design to that used by Multimission Modular Spacecraft for onboard attitude determination. An apparent obstacle to the use of unit vectors as measurements is their singular noise covariance matrix. Shuster has shown that this problem should not affect filter performance. The UVF and RTSF are tested using actual flight data from the Extreme Ultraviolet Explorer (EUVE). It is found that these filters generate nearly identical attitude and gyro bias estimates, thus validating the use of unit vectors and the QUEST noise model.

Sedlak, J.↗

A study of planetary meteorology

Preliminary estimates are prepared of the weighting functions of the infrared bands of methane, ammonia, and hydrogen for the Jovian atmosphere. Together with the infrared spectral emission observations of Jupiter, these weighting functions are used to obtain a rough estimate of the vertical temperature profile in the Jovian atmosphere.

Ohring, G.↗

A Fast Monte Carlo Method for Model-Based Prognostics Based on Stochastic Calculus

This work proposes a fast Monte Carlo method to solve differential equations utilized in model-based prognostics. The methodology is derived from the theory of stochastic calculus, and the goal of such a method is to speed up the estimation of the probability density functions describing the independent variable evolution over time. In the prognostic scenarios presented in this paper, the stochastic differential equations describe variables directly or indirectly related to the degradation of a monitored system. The method allows the estimation of the probability density functions by solving the deterministic equation and approximating the stochastic integrals using samples of the model noise. By so doing, the prognostic problem is solved without the Monte Carlo simulation based on Euler's forward method, which is typically the most time consuming task of the prediction stage. Three different prognostic scenarios are presented as proof of concept: (i) life prediction of electrolytic capacitors, (ii) remaining time to discharge of Lithium-ion batteries, and (iii) prognostic of cracked structures under fatigue loading. The paper shows how the method produces probability density functions that are statistically indistinguishable from the distributions estimated with Euler's forward Monte Carlo simulation. However, the proposed solution is orders of magnitude faster when computing the time-to-failure distribution of the monitored system. The approach may enable complex real-time prognostics and health management solutions with limited computing power.

Corbetta, M.↗