Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “error localization”

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 181 records · Page 10

Long Term Mean Local Time of the Ascending Node Prediction

Significant error has been observed in the long term prediction of the Mean Local Time of the Ascending Node on the Aqua spacecraft. This error of approximately 90 seconds over a two year prediction is a complication in planning and timing of maneuvers for all members of the Earth Observing System Afternoon Constellation, which use Aqua's MLTAN as the reference for their inclination maneuvers. It was determined that the source of the prediction error was the lack of a solid Earth tide model in the operational force models. The Love Model of the solid Earth tide potential was used to derive analytic corrections to the inclination and right ascension of the ascending node of Aqua's Sun-synchronous orbit. Additionally, it was determined that the resonance between the Sun and orbit plane of the Sun-synchronous orbit is the primary driver of this error. The analytic corrections have been added to the operational force models for the Aqua spacecraft reducing the two-year 90-second error to less than 7 seconds.

McKinley, David P.↗

Development of local hybrid density functionals to treat self-interaction error and many-electron effects

The goal of this project was to improve density functional theory (DFT) for systems where conventional semilocal approximations are limited by self-interaction error (SIE) and by near-degeneracy or strong many-electron effects. While DFT remains the only broadly practical first-principles framework for large-scale materials simulations, its predictive accuracy is often challenged in situations involving stretched bonds, charge transfer, transition-metal chemistry, magnetic couplings, band gaps, and correlated electronic states. This project addressed these limitations by developing physically grounded and numerically robust exchange–correlation functionals that retain the efficiency of modern DFT while extending its predictive scope.

36 MATERIALS SCIENCE↗

Slope-induced errors in radar altimetry over continental ice sheets

Altimeter measurements of ranges to modelled irregular surfaces are simulated and two correction schemes are used to reconstruct the modelled surfaces from the simulated data. When the cross-track slope is negligible, the problem can be treated as two-dimensional and the relocation method is preferred to the slope correction method. In a two-dimensional test case, 85 percent of the slope-induced rms error was removed by the relocation method. The mean error along profiles of about 75 km or longer is usually reduced more than the rms error. An alternative slope correction scheme, which uses the local slope to calculate the expected error, is less effective in the two-dimensional case. Over a simulated three-dimensional surface, where groundtracks are widely spaced and cross-track slopes are significant, the slope-correction method must be used in at least the cross-track direction.

Brenner, A. C.↗

A second catalog of gamma ray bursts: 1978 - 1980 localizations from the interplanetary network

Eighty-two gamma ray bursts were detected between 1978 September 14 and 1980 February 13 by the experiments of the interplanetary network (Prognoz 7, Venera 11 and 12 SIGNE experiments, Pioneer Venus Orbiter, International Sun-Earth Explorer 3, Helios 2, and Vela). Sixty-five of these events have been localized to annuli or error boxes by the method of arrival time analysis. The distribution of sources is consistent with isotropy, and there is no statistically convincing evidence for the detection of more than one burst from any source position. The localizations are compared with those of two previous catalogs.

Atteia, J. L.↗

Local Discontinuous Galerkin Methods for Partial Differential Equations with Higher Order Derivatives

In this paper we review the existing and develop new continuous Galerkin methods for solving time dependent partial differential equations with higher order derivatives in one and multiple space dimensions. We review local discontinuous Galerkin methods for convection diffusion equations involving second derivatives and for KdV type equations involving third derivatives. We then develop new local discontinuous Galerkin methods for the time dependent bi-harmonic type equations involving fourth derivatives, and partial differential equations involving fifth derivatives. For these new methods we present correct interface numerical fluxes and prove L(exp 2) stability for general nonlinear problems. Preliminary numerical examples are shown to illustrate these methods. Finally, we present new results on a post-processing technique, originally designed for methods with good negative-order error estimates, on the local discontinuous Galerkin methods applied to equations with higher derivatives. Numerical experiments show that this technique works as well for the new higher derivative cases, in effectively doubling the rate of convergence with negligible additional computational cost, for linear as well as some nonlinear problems, with a local uniform mesh.

Yan, Jue↗

The Interplanetary Network Supplement to the Fermi Gbm Catalog of Cosmic Gamma-Ray Bursts

We present Interplanetary Network (IPN) data for the gamma-ray bursts in the first Fermi Gamma-Ray BurstMonitor (GBM) catalog. Of the 491 bursts in that catalog, covering 2008 July 12 to 2010 July 11, 427 wereobserved by at least one other instrument in the nine-spacecraft IPN. Of the 427, the localizations of 149 could beimproved by arrival time analysis (or triangulation). For any given burst observed by the GBM and one otherdistant spacecraft, triangulation gives an annulus of possible arrival directions whose half-width varies betweenabout 0. 4 and 32, depending on the intensity, time history, and arrival direction of the burst, as well as the distancebetween the spacecraft. We find that the IPN localizations intersect the 1 GBM error circles in only 52 of thecases, if no systematic uncertainty is assumed for the latter. If a 6 systematic uncertainty is assumed and added inquadrature, the two localization samples agree about 87 of the time, as would be expected. If we then multiply theresulting error radii by a factor of three, the two samples agree in slightly over 98 of the cases, providing a goodestimate of the GBM 3 error radius. The IPN 3 error boxes have areas between about 1 arcmin2 and 110 deg2,and are, on the average, a factor of 180 smaller than the corresponding GBM localizations. We identify two burstsin the IPNGBM sample that did not appear in the GBM catalog. In one case, the GBM triggered on a terrestrialgamma flash, and in the other, its origin was given as uncertain. We also discuss the sensitivity and calibration ofthe IPN.

INTERPLANET↗

Control of nonlinear systems with applications to constrained robots and spacecraft attitude stabilization

This thesis is organized in two parts. In Part 1, control systems described by a class of nonlinear differential and algebraic equations are introduced. A procedure for local stabilization based on a local state realization is developed. An alternative approach to local stabilization is developed based on a classical linearization of the nonlinear differential-algebraic equations. A theoretical framework is established for solving a tracking problem associated with the differential-algebraic system. First, a simple procedure is developed for the design of a feedback control law which ensures, at least locally, that the tracking error in the closed loop system lies within any given bound if the reference inputs are sufficiently slowly varying. Next, by imposing additional assumptions, a procedure is developed for the design of a feedback control law which ensures that the tracking error in the closed loop system approaches zero exponentially for reference inputs which are not necessarily slowly varying. The control design methodologies are used for simultaneous force and position control in constrained robot systems. The differential-algebraic equations are shown to characterize the slow dynamics of a certain nonlinear control system in nonstandard singularly perturbed form. In Part 2, the attitude stabilization (reorientation) of a rigid spacecraft using only two control torques is considered. First, the case of momentum wheel actuators is considered. The complete spacecraft dynamics are not controllable. However, the spacecraft dynamics are small time locally controllable in a reduced sense. The reduced spacecraft dynamics cannot be asymptotically stabilized using continuous feedback, but a discontinuous feedback control strategy is constructed. Next, the case of gas jet actuators is considered. If the uncontrolled principal axis is not an axis of symmetry, the complete spacecraft dynamics are small time locally controllable. However, the spacecraft attitude cannot be asymptotically stabilized using continuous feedback, but a discontinuous stabilizing feedback control strategy is constructed. If the uncontrolled principal axis is an axis of symmetry, the complete spacecraft dynamics cannot be stabilized. However, the spacecraft dynamics are small time locally controllable in a reduced sense. The reduced spacecraft dynamics cannot be asymptotically stabilized using continuous feedback, but again a discontinuous feedback control strategy is constructed.

Krishnan, Hariharan↗

Psychophysical evaluation of three-dimensional auditory displays

Work during this reporting period included the completion of our research on the use of principal components analysis (PCA) to model the acoustical head related transfer functions (HRTFs) that are used to synthesize virtual sources for three dimensional auditory displays. In addition, a series of studies was initiated on the perceptual errors made by listeners when localizing free-field and virtual sources. Previous research has revealed that under certain conditions these perceptual errors, often called 'confusions' or 'reversals', are both large and frequent, thus seriously comprising the utility of a 3-D virtual auditory display. The long-range goal of our work in this area is to elucidate the sources of the confusions and to develop signal-processing strategies to reduce or eliminate them.

Wightman, Frederic L.↗

Image Correlation Pattern Optimization for Micro-Scale In-Situ Strain Measurements

The accuracy and precision of digital image correlation (DIC) is a function of three primary ingredients: image acquisition, image analysis, and the subject of the image. Development of the first two (i.e. image acquisition techniques and image correlation algorithms) has led to widespread use of DIC; however, fewer developments have been focused on the third ingredient. Typically, subjects of DIC images are mechanical specimens with either a natural surface pattern or a pattern applied to the surface. Research in the area of DIC patterns has primarily been aimed at identifying which surface patterns are best suited for DIC, by comparing patterns to each other. Because the easiest and most widespread methods of applying patterns have a high degree of randomness associated with them (e.g., airbrush, spray paint, particle decoration, etc.), less effort has been spent on exact construction of ideal patterns. With the development of patterning techniques such as microstamping and lithography, patterns can be applied to a specimen pixel by pixel from a patterned image. In these cases, especially because the patterns are reused many times, an optimal pattern is sought such that error introduced into DIC from the pattern is minimized. DIC consists of tracking the motion of an array of nodes from a reference image to a deformed image. Every pixel in the images has an associated intensity (grayscale) value, with discretization depending on the bit depth of the image. Because individual pixel matching by intensity value yields a non-unique scale-dependent problem, subsets around each node are used for identification. A correlation criteria is used to find the best match of a particular subset of a reference image within a deformed image. The reader is referred to references for enumerations of typical correlation criteria. As illustrated by Schreier and Sutton and Lu and Cary systematic errors can be introduced by representing the underlying deformation with under-matched shape functions. An important implication, as discussed by Sutton et al., is that in the presence of highly localized deformations (e.g., crack fronts), error can be reduced by minimizing the subset size. In other words, smaller subsets allow the more accurate resolution of localized deformations. Contrarily, the choice of optimal subset size has been widely studied and a general consensus is that larger subsets with more information content are less prone to random error. Thus, an optimal subset size balances the systematic error from under matched deformations with random error from measurement noise. The alternative approach pursued in the current work is to choose a small subset size and optimize the information content within (i.e., optimizing an applied DIC pattern), rather than finding an optimal subset size. In the literature, many pattern quality metrics have been proposed, e.g., sum of square intensity gradient (SSSIG), mean subset fluctuation, gray level co-occurrence, autocorrelation-based metrics, and speckle-based metrics. The majority of these metrics were developed to quantify the quality of common pseudo-random patterns after they have been applied, and were not created with the intent of pattern generation. As such, it is found that none of the metrics examined in this study are fit to be the objective function of a pattern generation optimization. In some cases, such as with speckle-based metrics, application to pixel by pixel patterns is ill-conditioned and requires somewhat arbitrary extensions. In other cases, such as with the SSSIG, it is shown that trivial solutions exist for the optimum of the metric which are ill-suited for DIC (such as a checkerboard pattern). In the current work, a multi-metric optimization method is proposed whereby quality is viewed as a combination of individual quality metrics. Specifically, SSSIG and two auto-correlation metrics are used which have generally competitive objectives. Thus, each metric could be viewed as a constraint imposed upon the others, thereby precluding the achievement of their trivial solutions. In this way, optimization produces a pattern which balances the benefits of multiple quality metrics. The resulting pattern, along with randomly generated patterns, is subjected to numerical deformations and analyzed with DIC software. The optimal pattern is shown to outperform randomly generated patterns.

Bomarito, G. F.↗

Efficient Implementation of an Optimal Interpolator for Large Spatial Data Sets

Scattered data interpolation is a problem of interest in numerous areas such as electronic imaging, smooth surface modeling, and computational geometry. Our motivation arises from applications in geology and mining, which often involve large scattered data sets and a demand for high accuracy. The method of choice is ordinary kriging. This is because it is a best unbiased estimator. Unfortunately, this interpolant is computationally very expensive to compute exactly. For n scattered data points, computing the value of a single interpolant involves solving a dense linear system of size roughly n x n. This is infeasible for large n. In practice, kriging is solved approximately by local approaches that are based on considering only a relatively small'number of points that lie close to the query point. There are many problems with this local approach, however. The first is that determining the proper neighborhood size is tricky, and is usually solved by ad hoc methods such as selecting a fixed number of nearest neighbors or all the points lying within a fixed radius. Such fixed neighborhood sizes may not work well for all query points, depending on local density of the point distribution. Local methods also suffer from the problem that the resulting interpolant is not continuous. Meyer showed that while kriging produces smooth continues surfaces, it has zero order continuity along its borders. Thus, at interface boundaries where the neighborhood changes, the interpolant behaves discontinuously. Therefore, it is important to consider and solve the global system for each interpolant. However, solving such large dense systems for each query point is impractical. Recently a more principled approach to approximating kriging has been proposed based on a technique called covariance tapering. The problems arise from the fact that the covariance functions that are used in kriging have global support. Our implementations combine, utilize, and enhance a number of different approaches that have been introduced in literature for solving large linear systems for interpolation of scattered data points. For very large systems, exact methods such as Gaussian elimination are impractical since they require 0(n(exp 3)) time and 0(n(exp 2)) storage. As Billings et al. suggested, we use an iterative approach. In particular, we use the SYMMLQ method, for solving the large but sparse ordinary kriging systems that result from tapering. The main technical issue that need to be overcome in our algorithmic solution is that the points' covariance matrix for kriging should be symmetric positive definite. The goal of tapering is to obtain a sparse approximate representation of the covariance matrix while maintaining its positive definiteness. Furrer et al. used tapering to obtain a sparse linear system of the form Ax = b, where A is the tapered symmetric positive definite covariance matrix. Thus, Cholesky factorization could be used to solve their linear systems. They implemented an efficient sparse Cholesky decomposition method. They also showed if these tapers are used for a limited class of covariance models, the solution of the system converges to the solution of the original system. Matrix A in the ordinary kriging system, while symmetric, is not positive definite. Thus, their approach is not applicable to the ordinary kriging system. Therefore, we use tapering only to obtain a sparse linear system. Then, we use SYMMLQ to solve the ordinary kriging system. We show that solving large kriging systems becomes practical via tapering and iterative methods, and results in lower estimation errors compared to traditional local approaches, and significant memory savings compared to the original global system. We also developed a more efficient variant of the sparse SYMMLQ method for large ordinary kriging systems. This approach adaptively finds the correct local neighborhood for each query point in the interpolation process.

Memarsadeghi, Nargess↗

Periodic variations stratospheric temperature from 20-65 km at 80 deg N to 30 deg S

Results for a seasonally varying diurnal tide in temperature at Churchill are presented, and possible significant aliasing of longer period waves by this tide is discussed. A diurnal tide whose amplitude and phase are coherent throughout the year is found to have little effect on periodic amplitudes other than the long-term mean, because most rocketsonde observations are taken near the same local time each day. Errors in periodic components arising from lack of solar radiation corrections are found to be largest for the long-term mean with a small influence noted in the annual wave's amplitude. Spatial variations of the amplitudes and phases of long-period waves are examined through the use of height-latitude sections, 20-65 km, at 80 deg N to 30 deg S. The quasi-biennial oscillation and semiannual waves have tropical maxima of 2 and 3C near 30 and 40 km respectively. The annual wave's maximum is over 22C near 45 km at 70 deg N and the terannual wave's maximum is over 6C near 55 km at 80 deg N. The semiannual wave has to polar maxima: 7C near 75 deg N at 32 km and 3C above 60 km north of 35 deg N.

Nastrom, G. D.↗

Mach and Reynolds number effects on turbulent skin friction reduction by injection

The investigation reported is concerned with questions regarding a possible Mach number influence on skin friction reduction caused by injection. The investigation shows that data considered by Danberg (1967) for the no-blowing skin friction coefficient are in error. Accurate profiles and local skin friction coefficient values are obtained when the influence of low Reynolds number amplification in the outer region of the boundary layer is included in a calculation method.

Bushnell, D. M.↗

Exospheric density measurements from the drag-free satellite Triad

The Triad satellite is maintained within 1 mm of its proof mass (a small metal ball) in a purely gravitational orbit, using a drag-free control system, named Discos, as the disturbance compensation system. The Discos proof mass is surrounded by an outer shell which holds fuel tanks and cold gas jets. The shell shields the proof mass from such nongravitational forces as radiation pressure, atmospheric drag, and micrometeorite impact. Whenever these forces displace the outer shell relative to the proof mass, an error signal is generated by a capacitive bridge sensor, and the satellite is propelled by gas jets to remain centered on the proof mass when the error reaches 1 mm. Local atmospheric densities near 800 km were measured, using ball position data, and the observed accelerations were corrected for the effects of solar radiation pressure. The measured densities reveal a greater dependence on latitude than that indicated by balloon satellites, and also a dependence on longitude. Only a small dependence on Kp, however, was observed in the region where the measurements were made.

Moe, K.↗

Digital image correlation techniques applied to LANDSAT multispectral imagery

The author has identified the following significant results. Automatic image registration and resampling techniques applied to LANDSAT data achieved accuracies, resulting in mean radial displacement errors of less than 0.2 pixel. The process method utilized recursive computational techniques and line-by-line updating on the basis of feedback error signals. Goodness of local feature matching was evaluated through the implementation of a correlation algorithm. An automatic restart allowed the system to derive control point coordinates over a portion of the image and to restart the process, utilizing this new control point information as initial estimates.

Bonrud, L. O.↗

A preliminary estimate of geoid-induced variations in repeat orbit satellite altimeter observations

Altimeter satellites are often maintained in a repeating orbit to facilitate the separation of sea-height variations from the geoid. However, atmospheric drag and solar radiation pressure cause a satellite orbit to drift. For Geosat this drift causes the ground track to vary by + or - 1 km about the nominal repeat path. This misalignment leads to an error in the estimates of sea surface height variations because of the local slope in the geoid. This error has been estimated globally for the Geosat Exact Repeat Mission using a mean sea surface constructed from Geos 3 and Seasat altimeter data. Over most of the ocean the geoid gradient is small, and the repeat-track misalignment leads to errors of only 1 to 2 cm. However, in the vicinity of trenches, continental shelves, islands, and seamounts, errors can exceed 20 cm. The estimated error is compared with direct estimates from Geosat altimetry, and a strong correlation is found in the vicinity of the Tonga and Aleutian trenches. This correlation increases as the orbit error is reduced because of the increased signal-to-noise ratio.

Brenner, Anita C.↗

On High-Order Radiation Boundary Conditions

In this paper we develop the theory of high-order radiation boundary conditions for wave propagation problems. In particular, we study the convergence of sequences of time-local approximate conditions to the exact boundary condition, and subsequently estimate the error in the solutions obtained using these approximations. We show that for finite times the Pade approximants proposed by Engquist and Majda lead to exponential convergence if the solution is smooth, but that good long-time error estimates cannot hold for spatially local conditions. Applications in fluid dynamics are also discussed.

Hagstrom, Thomas↗

Eulerian-Lagrangian Simulations of Transonic Flutter Instabilities

This paper presents an overview of recent applications of Eulerian-Lagrangian computational schemes in simulating transonic flutter instabilities. This approach, the fluid-structure system is treated as a single continuum dynamics problem, by switching from an Eulerian to a Lagrangian formulation at the fluid-structure boundary. This computational approach effectively eliminates the phase integration errors associated with previous methods, where the fluid and structure are integrated sequentially using different schemes. The formulation is based on Hamilton's Principle in mixed coordinates, and both finite volume and finite element discretization schemes are considered. Results from numerical simulations of transonic flutter instabilities are presented for isolated wings, thin panels, and turbomachinery blades. The results suggest that the method is capable of reproducing the energy exchange between the fluid and the structure with significantly less error than existing methods. Localized flutter modes and panel flutter modes involving traveling waves can also be simulated effectively with no a priori knowledge of the type of instability involved.

Bendiksen, Oddvar O.↗