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 595 records · Page 33

Shape Estimation for Elongated Deformable Object using B-spline Chained Multiple Random Matrices Model

In this paper, a B-spline chained multiple random matrix models (RMMs) representation is proposed to model geometric characteristics of an elongated deformable object. The hyper degrees of freedom structure of the elongated deformable object make its shape estimation challenging. Based on the likelihood function of the proposed B-spline chained multiple RMMs, an expectation-maximization (EM) method is derived to estimate the shape of the elongated deformable object. A split and merge method based on the Euclidean minimum spanning tree (EMST) is proposed to provide initialization for the EM algorithm. The proposed algorithm is evaluated for the shape estimation of the elongated deformable objects in scenarios, such as the static rope with various configurations (including configurations with intersection), the continuous manipulation of a rope and a plastic tube, and the assembly of two plastic tubes. The execution time is computed and the accuracy of the shape estimation results is evaluated based on the comparisons between the estimated width values and its ground-truth, and the intersection over union (IoU) metric.

Gang Yao↗

Ion Diffusive Transport Across the Separatrix Between the Low-Latitude Mantle and the Plasma Sheet By Kinetic Alfvén Waves: MMS Observation

To understand the entry of the cool low‐latitude mantle ions into the tail plasma sheet near the flanks under persistent interplanetary magnetic field B y , we evaluate the role of the cross‐field diffusive transport by kinetic Alfvén waves (KAWs) by investigating two events observed by multiscale (MMS) spacecraft. Around the separatrix between the open and closed field‐line regions, a two‐component mixing of hot plasma sheet ions of a few keV with cool mantle ions of a few hundred eV was observed, indicating transport across the separatrix. The waves observed between 0.01 and 10 Hz around the separatrix had characteristics consistent with those of KAWs. The consistency allowed us to estimate the wave vectors as a function of frequency by fitting KAW dispersion to the observations. Using the observed wave powers, plasma moments, and the estimated wave vectors, we computed the cross‐field diffusion rates associated with KAWs. The diffusion rates were found to be comparable to or larger than the Bohm diffusion rates during the intervals when the two‐component mixing was observed, indicating that the KAW diffusive transport can play a role in the entry of low‐latitude mantle ions into the plasma sheet.

Kinetic Alfvén waves↗

The Projection of Space Radiation Environments with a Solar Cycle Statistical Model

A solar cycle statistical model has been developed to project sunspot numbers which represent the variations in the space radiation environment. The resultant projection of sunspot numbers in near future were coupled to space-related quantities of interest in radiation protection, such as the galactic cosmic radiation (GCR) deceleration potential (f) and the mean occurrence frequency of solar particle event (SPE). Future GCR fluxes have been derived from a predictive model, in which GCR temporal dependence represented by f was derived from GCR flux and ground-based Climax neutron monitor rate measurements over the last four decades. Results showed that the point dose equivalent inside a typical spacecraft in interplanetary radiation fields was influenced by solar modulation up to a factor of three. One important characteristic of sporadic SPEs is their mean frequency of occurrence, which is dependent on solar activity. Projections of future mean frequency of SPE occurrence were estimated from a power law function of sunspot number. Furthermore, the cumulative probabilities of SPE during short-period missions were defined with the continuous database of proton fluences of SPE. The analytic representation of energy spectra of SPE was constructed by the Weibull distribution for different event sizes. The representative exposure level at each event size was estimated for the guideline of protection systems for astronauts during future space exploration missions.

Kim, Myung-Hee↗

Overcoming Dynamic Disturbances in Imaging Systems

We develop and discuss a methodology with the potential to yield a significant reduction in complexity, cost, and risk of space-borne optical systems in the presence of dynamic disturbances. More robust systems almost certainly will be a result as well. Many future space-based and ground-based optical systems will employ optical control systems to enhance imaging performance. The goal of the optical control subsystem is to determine the wavefront aberrations and remove them. Ideally reducing an aberrated image of the object under investigation to a sufficiently clear (usually diffraction-limited) image. Control will likely be distributed over several elements. These elements may include telescope primary segments, telescope secondary, telescope tertiary, deformable mirror(s), fine steering mirror(s), etc. The last two elements, in particular, may have to provide dynamic control. These control subsystems may become elaborate indeed. But robust system performance will require evaluation of the image quality over a substantial range and in a dynamic environment. Candidate systems for improvement in the Earth Sciences Enterprise could include next generation Landsat systems or atmospheric sensors for dynamic imaging of individual, severe storms. The technology developed here could have a substantial impact on the development of new systems in the Space Science Enterprise; such as the Next Generation Space Telescope(NGST) and its follow-on the Next NGST. Large Interferometric Systems of non-zero field, such as Planet Finder and Submillimeter Probe of the Evolution of Cosmic Structure, could benefit. These systems most likely will contain large, flexible optomechanical structures subject to dynamic disturbance. Furthermore, large systems for high resolution imaging of planets or the sun from space may also benefit. Tactical and Strategic Defense systems will need to image very small targets as well and could benefit from the technology developed here. We discuss a novel speckle imaging technique with the potential to separate dynamic aberrations from static aberrations. Post-processing of a set of image data, using an algorithm based on this technique, should work for all but the lowest light levels and highest frequency dynamic environments. This technique may serve to reduce the complexity of the control system and provide for robust, fault-tolerant, reduced risk operation. For a given object, a short exposure image is "frozen" on the focal plane in the presence of the environmental disturbance (turbulence, jitter, etc.). A key factor is that this imaging data exhibits frame-to-frame linear shift invariance. Therefore, although the Point Spread Function is varying from frame to frame, the source is fixed; and each short exposure contains object spectrum data out to the diffraction limit of the imaging system. This novel speckle imaging technique uses the Knox-Thompson method. The magnitude of the complex object spectrum is straightforward to determine by well-established approaches. The phase of the complex object spectrum is decomposed into two parts. One is a single-valued function determined by the divergence of the optical phase gradient. The other is a multi-valued function determined by the circulation of the optical phase gradient-"hidden phase." Finite difference equations are developed for the phase. The novelty of this approach is captured in the inclusion of this "hidden phase." This technique allows the diffraction-limited reconstruction of the object from the ensemble of short exposure frames while simultaneously estimating the phase as a function of time from a set of exposures.

Young, Eric W.↗

Overcoming Dynamic Disturbances in Imaging Systems

We develop and discuss a methodology with the potential to yield a significant reduction in complexity, cost, and risk of space-borne optical systems in the presence of dynamic disturbances. More robust systems almost certainly will be a result as well. Many future space-based and ground-based optical systems will employ optical control systems to enhance imaging performance. The goal of the optical control subsystem is to determine the wavefront aberrations and remove them. Ideally reducing an aberrated image of the object under investigation to a sufficiently clear (usually diffraction-limited) image. Control will likely be distributed over several elements. These elements may include telescope primary segments, telescope secondary, telescope tertiary, deformable mirror(s), fine steering mirror(s), etc. The last two elements, in particular, may have to provide dynamic control. These control subsystems may become elaborate indeed. But robust system performance will require evaluation of the image quality over a substantial range and in a dynamic environment. Candidate systems for improvement in the Earth Sciences Enterprise could include next generation Landsat systems or atmospheric sensors for dynamic imaging of individual, severe storms. The technology developed here could have a substantial impact on the development of new systems in the Space Science Enterprise; such as the Next Generation Space Telescope(NGST) and its follow-on the Next NGST. Large Interferometric Systems of non-zero field, such as Planet Finder and Submillimeter Probe of the Evolution of Cosmic Structure, could benefit. These systems most likely will contain large, flexible optormechanical structures subject to dynamic disturbance. Furthermore, large systems for high resolution imaging of planets or the sun from space may also benefit. Tactical and Strategic Defense systems will need to image very small targets as well and could benefit from the technology developed here. We discuss a novel speckle imaging technique with the potential to separate dynamic aberrations from static aberrations. Post-processing of a set of image data, using an algorithm based on this technique, should work for all but the lowest light levels and highest frequency dynamic environments. This technique may serve to reduce the complexity of the control system and provide for robust, fault-tolerant, reduced risk operation. For a given object, a short exposure image is "frozen" on the focal plane in the presence of the environmental disturbance (turbulence, jitter, etc.). A key factor is that this imaging data exhibits frame-to-frame linear shift invariance. Therefore, although the Point Spread Function is varying from frame to frame, the source is fixed; and each short exposure contains object spectrum data out to the diffraction limit of the imaging system. This novel speckle imaging technique uses the Knox-Thompson method. The magnitude of the complex object spectrum is straightforward to determine by well-established approaches. The phase of the complex object spectrum is decomposed into two parts. One is a single-valued function determined by the divergence of the optical phase gradient. The other is a multi-valued function determined by, the circulation of the optical phase gradient-"hidden phase." Finite difference equations are developed for the phase. The novelty of this approach is captured in the inclusion of this "hidden phase." This technique allows the diffraction-limited reconstruction of the object from the ensemble of short exposure frames while simultaneously estimating the phase as a function of time from a set of exposures.

Young, Eric W.↗

Sampling and Reconstruction of the Pupil and Electric Field for Phase Retrieval

This technology is based on sampling considerations for a band-limited function, which has application to optical estimation generally, and to phase retrieval specifically. The analysis begins with the observation that the Fourier transform of an optical aperture function (pupil) can be implemented with minimal aliasing for Q values down to Q = 1. The sampling ratio, Q, is defined as the ratio of the sampling frequency to the band-limited cut-off frequency. The analytical results are given using a 1-d aperture function, and with the electric field defined by the band-limited sinc(x) function. Perfect reconstruction of the Fourier transform (electric field) is derived using the Whittaker-Shannon sampling theorem for 1<Q<2. The Fourier transform is constructed by periodic extension, i.e., by spacing copies of the transform in a definite way, recognizing that no aliasing occurs for values of the sampling ratio such that 1<Q<2, which can be used to advantage in the application of phase retrieval estimation. A method was developed for propagating the electro magnetic field with no aliasing, which has been extended to 2-d optical apertures.

Dean, Bruce↗

Fabrication and properties of functionally graded NiAl/Al2O3 composites

A modified sedimentation process was used in the production of a functionally gradient material (FGM), NiAl/Al2O3. A simple finite element model was used to guide our design and fabrication efforts by estimating residual stress states as a function of composite structure. This approach could lead to tailored designs that enhance or avoid specific residual stress states. Thermal cycling tests were factored into the model to predict time dependent or steady-state internal temperature and stress profiles. Four-point bend tests were conducted to establish the mechanical load-displacement behavior of a single interlayer FGM at room temperature, 800 and 1000 K. Room temperature bend strength of the FGM was 3-4 times that of the base NiAl. At elevated temperatures, composite fracture occurred in a gradual, noncatastrophic mode involving NiAl retardation of a succession of cracks originating in the alumina face.

Miller, D. P.↗

Cloverleaf microgyroscope with electrostatic alignment and tuning

A micro-gyroscope (10) having closed loop output operation by a control voltage (V.sub.ty), that is demodulated by a drive axis (x-axis) signal V.sub.thx of the sense electrodes (S1, S2), providing Coriolis torque rebalance to prevent displacement of the micro-gyroscope (10) on the output axis (y-axis) V.sub.thy.about.0. Closed loop drive axis torque, V.sub.tx maintains a constant drive axis amplitude signal, V.sub.thx. The present invention provides independent alignment and tuning of the micro-gyroscope by using separate electrodes and electrostatic bias voltages to adjust alignment and tuning. A quadrature amplitude signal, or cross-axis transfer function peak amplitude is used to detect misalignment that is corrected to zero by an electrostatic bias voltage adjustment. The cross-axis transfer function is either V.sub.thy/V.sub.ty or V.sub.tnx/V.sub.tx. A quadrature signal noise level, or difference in natural frequencies estimated from measurements of the transfer functions is used to detect residual mistuning, that is corrected to zero by a second electrostatic bias voltage adjustment.

Challoner, A. Dorian↗

Separating Direct and Indirect Turbofan Engine Combustion Noise While Estimating Post-Combustion (Post-Flame) Residence Time Using the Correlation Function

A previous investigation on the presence of direct and indirect combustion noise for a full-scale turbofan engine using a far-field microphone at 130 is extended by also examining signals obtained at two additional downstream directions using far-field microphones at 110 deg and 160 deg. A generalized cross-correlation function technique is used to study the change in propagation time to the far field of the combined direct and indirect combustion noise signal as a sequence of low-pass filters are applied. The filtering procedure used produces no phase distortion. As the low-pass filter frequency is decreased, the travel time increases because the relative amount of direct combustion noise is reduced. The indirect combustion noise signal travels more slowly because in the combustor entropy fluctuations move with the flow velocity, which is slow compared to the local speed of sound. The indirect combustion noise signal travels at acoustic velocities after reaching the turbine and being converted into an acoustic signal. The direct combustion noise is always propagating at acoustic velocities. The results show that the estimated indirect combustion noise time delay values (post-combustion residence times) measured at each angle are fairly consistent with one another for a relevant range of operating conditions and demonstrate source separation of a mixture of direct and indirect combustion noise. The results may lead to a better idea about the acoustics in the combustor and may help develop and validate improved reduced-order physics-based methods for predicting turbofan engine core noise.

Miles, Jeffrey Hilton↗

Gaussian windows: A tool for exploring multivariate data

Presented here is a method for interactively exploring a large set of quantitative multivariate data, in order to estimate the shape of the underlying density function. It is assumed that the density function is more or less smooth, but no other specific assumptions are made concerning its structure. The local structure of the data in a given region may be examined by viewing the data through a Gaussian window, whose location and shape are chosen by the user. A Gaussian window is defined by giving each data point a weight based on a multivariate Gaussian function. The weighted sample mean and sample covariance matrix are then computed, using the weights attached to the data points. These quantities are used to compute an estimate of the shape of the density function in the window region. The local structure of the data is described by a method similar to the method of principal components. By taking many such local views of the data, we can form an idea of the structure of the data set. The method is applicable in any number of dimensions. The method can be used to find and describe simple structural features such as peaks, valleys, and saddle points in the density function, and also extended structures in higher dimensions. With some practice, we can apply our geometrical intuition to these structural features in any number of dimensions, so that we can think about and describe the structure of the data. Since the computations involved are relatively simple, the method can easily be implemented on a small computer.

Jaeckel, Louis A.↗

Analysis of condensates formed at the Viking 2 lander site - The first winter

Relative surface albedo, spectral reflectance estimates and a limited photometric function are reduced from Viking 2 lander data obtained during a 249 Mars day period, in the lander's first year, when a light ground covering appeared on the surface. During the deposition, surface broadband albedo more than doubled and blue reflectance increased by a factor of 4.0. Comparison of lander data with earlier laboratory measurements of CO2 and H2O frosts and snows shows that reflectance estimates do not resemble either of those solids. The condensate reflectance resembles that of the surface after the covering disappeared. The covering may have been colored by dust which fell before it, by dust mixed with it, or by dust on top of it; but the data strongly support a mixture of dust with H2O and CO2 solids. The covering thickness is estimated to be between 0.5 and a few millimeters.

Wall, S. D.↗

Estimating three-demensional energy transfer in isotropic turbulence

To obtain an estimate of the spectral transfer function that indicates the rate of decay of energy, an x-wire probe was set at a fixed position, and two single wire probes were set at a number of locations in the same plane perpendicular to the mean flow in the wind tunnel. The locations of the single wire probes are determined by pseudo-random numbers (Monte Carlo). Second order spectra and cross spectra are estimated. The assumption of isotropy relative to second order spectra is examined. Third order spectra are also estimated corresponding to the positions specified. A Monte Carlo Fourier transformation of the downstream bispectra corresponding to integration across the plane perpendicular to the flow is carried out assuming isotropy. Further integration is carried out over spherical energy shells.

Li, K. S.↗

Evaluation of Preliminary Buffet Forcing Function Development based on Wind-Tunnel Tests of Geometrically-similar Models

The current state-of-the-art for launch vehicle buffet load estimation involves wind tunnel testing of highly instrumented rigid models and data processing of experimental unsteady pressure measurements to estimate the fluctuating aerodynamic loads. The model and instrumentation requirements to accomplish this type of testing are significant. The rigid model, in particular, must be built with relatively high geometric fidelity to approximate similar aerodynamic phenomena that will occur during the full-scale vehicle flight. Due to this requirement, several iterative wind tunnel tests are often required during the design maturation of the vehicle. For preliminary buffet load estimation, however, data from previous launch vehicle wind tunnel tests of geometrically similar models can be used to provide initial buffet environment estimates before a high geometric fidelity wind tunnel test has been conducted. In this paper, two components of this preliminary buffet estimation process are described: (1) assembly of buffet forcing functions from multiple data sources and (2) scaling of these buffet forcing functions to geometrically-similar components of a target vehicle configuration. The impacts of these estimation methods on the full-vehicle buffet forcing functions, including their longitudinal distribution along the vehicle, frequency-content at individual stations, and station-to-station coherence, are investigated. Considerations for using these methods to extrapolate wind tunnel data to geometrically similar models are also discussed.

Patrick S Heaney↗

Adjoint-Based, Three-Dimensional Error Prediction and Grid Adaptation

Engineering computational fluid dynamics (CFD) analysis and design applications focus on output functions (e.g., lift, drag). Errors in these output functions are generally unknown and conservatively accurate solutions may be computed. Computable error estimates can offer the possibility to minimize computational work for a prescribed error tolerance. Such an estimate can be computed by solving the flow equations and the linear adjoint problem for the functional of interest. The computational mesh can be modified to minimize the uncertainty of a computed error estimate. This robust mesh-adaptation procedure automatically terminates when the simulation is within a user specified error tolerance. This procedure for estimating and adapting to error in a functional is demonstrated for three-dimensional Euler problems. An adaptive mesh procedure that links to a Computer Aided Design (CAD) surface representation is demonstrated for wing, wing-body, and extruded high lift airfoil configurations. The error estimation and adaptation procedure yielded corrected functions that are as accurate as functions calculated on uniformly refined grids with ten times as many grid points.

Park, Michael A.↗

Modeling, Monitoring and Fault Diagnosis of Spacecraft Air Contaminants

Control of air contaminants is a crucial factor in the safety considerations of crewed space flight. Indoor air quality needs to be closely monitored during long range missions such as a Mars mission, and also on large complex space structures such as the International Space Station. This work mainly pertains to the detection and simulation of air contaminants in the space station, though much of the work is easily extended to buildings, and issues of ventilation systems. Here we propose a method with which to track the presence of contaminants using an accurate physical model, and also develop a robust procedure that would raise alarms when certain tolerance levels are exceeded. A part of this research concerns the modeling of air flow inside a spacecraft, and the consequent dispersal pattern of contaminants. Our objective is to also monitor the contaminants on-line, so we develop a state estimation procedure that makes use of the measurements from a sensor system and determines an optimal estimate of the contamination in the system as a function of time and space. The real-time optimal estimates in turn are used to detect faults in the system and also offer diagnoses as to their sources. This work is concerned with the monitoring of air contaminants aboard future generation spacecraft and seeks to satisfy NASA's requirements as outlined in their Strategic Plan document (Technology Development Requirements, 1996).

Ramirez, W. Fred↗

A Ground Footprint Eccentricity Model For Asteroid Airbursts

Uncertainties in early observations of potentially hazardous asteroids result in preliminary impact corridors that can stretch across large portions of the Earth’s surface. At this early stage of detection, the corridor width and potential for damage are typically estimated using techniques from nuclear weapons research. These estimates often employ spherical blast assumptions resulting in a constant width impact corridor. In actuality, however, the ground damage footprint of obliquely entering asteroids is generally roughly elliptical or “butterfly” shaped, with the major axis extending in the cross range direction and the minor axis aligned with ground-track of the meteoroid. Since actual ground footprints for oblique entries may have aspect ratios greater than two or three, the assumption of a circular blast may significantly underestimate the area of the impact swath and the at-risk population. This work develops an engineering model that can be used to quickly estimate the eccentricity of the ground footprint as a function of local impact parameters. This yields vastly improved local estimates of the corridor width and can significantly enhance the accuracy of risk analysis.

Aftosmis, Michael J.↗

Exploring the Connection Between Sampling Problems in Bayesian Inference and Statistical Mechanics

The Bayesian and statistical mechanical communities often share the same objective in their work - estimating and integrating probability distribution functions (pdfs) describing stochastic systems, models or processes. Frequently, these pdfs are complex functions of random variables exhibiting multiple, well separated local minima. Conventional strategies for sampling such pdfs are inefficient, sometimes leading to an apparent non-ergodic behavior. Several recently developed techniques for handling this problem have been successfully applied in statistical mechanics. In the multicanonical and Wang-Landau Monte Carlo (MC) methods, the correct pdfs are recovered from uniform sampling of the parameter space by iteratively establishing proper weighting factors connecting these distributions. Trivial generalizations allow for sampling from any chosen pdf. The closely related transition matrix method relies on estimating transition probabilities between different states. All these methods proved to generate estimates of pdfs with high statistical accuracy. In another MC technique, parallel tempering, several random walks, each corresponding to a different value of a parameter (e.g. "temperature"), are generated and occasionally exchanged using the Metropolis criterion. This method can be considered as a statistically correct version of simulated annealing. An alternative approach is to represent the set of independent variables as a Hamiltonian system. Considerab!e progress has been made in understanding how to ensure that the system obeys the equipartition theorem or, equivalently, that coupling between the variables is correctly described. Then a host of techniques developed for dynamical systems can be used. Among them, probably the most powerful is the Adaptive Biasing Force method, in which thermodynamic integration and biased sampling are combined to yield very efficient estimates of pdfs. The third class of methods deals with transitions between states described by rate constants. These problems are isomorphic with chemical kinetics problems. Recently, several efficient techniques for this purpose have been developed based on the approach originally proposed by Gillespie. Although the utility of the techniques mentioned above for Bayesian problems has not been determined, further research along these lines is warranted

Pohorille, Andrew↗

How to Estimate Attitude from Vector Observations

The most robust estimators minimizing Wahba's loss function are Davenport's q method and the Singular Value Decomposition (SVD) method. The q method is faster than the SVD method with three or more measurements. The other algorithms are less robust since they solve the characteristic polynomial equation to find the maximum eigenvalue of Davenport's K matrix. They are only preferable when speed or processor power is an important consideration. Of these, Fast Optimal Attitude Matrix (FOAM) is the most robust and faster than the q method. Robustness is only an issue for measurements with widely differing accuracies, so the fastest algorithms, Quaternion ESTimator (QUEST), EStimator of the Optimal Quaternion (ESOQ), and ESOQ2, are well suited to star sensor applications.

Markley, F. Landis↗