Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “error propagation”

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
All sourcesNASA NTRSDOE OSTI
NASA NTRS2004

Adaptive Low Dissipative High Order Filter Methods for Multiscale MHD Flows

Adaptive low-dissipative high order filter finite difference methods for long time wave propagation of shock/turbulence/combustion compressible viscous MHD flows has been constructed. Several variants of the filter approach that cater to different flow types are proposed. These filters provide a natural and efficient way for the minimization of the divergence of the magnetic field [divergence of B] numerical error in the sense that no standard divergence cleaning is required. For certain 2-D MHD test problems, divergence free preservation of the magnetic fields of these filter schemes has been achieved.

Yee, H. C.↗
NASA NTRS2006

Adapting Covariance Propagation to Account for the Presence of Modeled and Unmodeled Maneuvers

This paper explores techniques that can be used to adapt the standard linearized propagation of an orbital covariance matrix to the case where there is a maneuver and an associated execution uncertainty. A Monte Carlo technique is used to construct a final orbital covariance matrix for a 'propagate-burn-propagate' process that takes into account initial state uncertainty and execution uncertainties in the maneuver magnitude. This final orbital covariance matrix is regarded as 'truth' and comparisons between it and three methods using modified linearized covariance propagation are made. The first method accounts for the maneuver by modeling its nominal effect within the state transition matrix but excludes the execution uncertainty by omitting a process noise matrix from the computation. In the second method, the maneuver is not modeled but the uncertainty in its magnitude is accounted for by the inclusion of a process noise matrix. In the third method, which is essentially a hybrid of the first two, the nominal portion of the maneuver is included via the state transition matrix while a process noise matrix is used to account for the magnitude uncertainty. Since this method also correctly accounts for the presence of the maneuver in the nominal orbit, it is the best method for applications involving the computation of times of closest approach and the corresponding probability of collision, Pc. However, applications for the two other methods exist and are briefly discussed. Despite the fact that the process model ('propagate-burn-propagate') that was studied was very simple - point-mass gravitational effects due to the Earth combined with an impulsive delta-V in the velocity direction for the maneuver - generalizations to more complex scenarios, including high fidelity force models, finite duration maneuvers, and maneuver pointing errors, are straightforward and are discussed in the conclusion.

Schiff, Conrad↗
NASA NTRS1994

Optimizations on Designing High-Resolution Finite-Difference Schemes

We describe a general optimization procedure for both maximizing the resolution characteristics of existing finite differencing schemes as well as designing finite difference schemes that will meet the error tolerance requirements of numerical solutions. The procedure is based on an optimization process. This is a generalization of the compact scheme introduced by Lele in which the resolution is improved for single, one-dimensional spatial derivative, whereas in the present approach the complete scheme, after spatial and temporal discretizations, is optimized on a range of parameters of the scheme and the governing equations. The approach is to linearize and Fourier analyze the discretized equations to check the resolving power of the scheme for various wave number ranges in the solution and optimize the resolution to satisfy the requirements of the problem. This represents a constrained nonlinear optimization problem which can be solved to obtain the nodal weights of discretization. An objective function is defined in the parametric space of wave numbers, Courant number, Mach number and other quantities of interest. Typical criterion for defining the objective function include the maximization of the resolution of high wave numbers for acoustic and electromagnetic wave propagations and turbulence calculations. The procedure is being tested on off-design conditions of non-uniform mesh, non-periodic boundary conditions, and non-constant wave speeds for scalar and system of equations. This includes the solution of wave equations and Euler equations using a conventional scheme with and without optimization and the design of an optimum scheme for the specified error tolerance.

Liu, Yen↗
NASA NTRS2020

Detecting Seasonal Ice Dynamics in Satellite Images

Fully understanding how glaciers respond to environmental change will require new methods to help us identify the onset of ice acceleration events and observe how dynamic signals propagate within glaciers. In particular, observations of ice dynamics on seasonal timescales may offer insights into how a glacier interacts with various forcing mechanisms throughout the year. The task of generating continuous ice velocity time series that resolve seasonal variability is made more difficult by a spotty satellite record that contains no optical observations throughout the dark, polar winters. Furthermore, velocities obtained by feature tracking are marked by high noise when image pairs are separated by short time intervals and contain no direct insights into variability that occurs between images separated by long time intervals. In this paper, we describe a method of analyzing optical or SAR-derived feature-tracked velocities to characterize the magnitude and timing of seasonal ice dynamic variability. Our method is agnostic to data gaps and is able to recover climatological average winter velocities regardless of the availability of direct observations during winter. Using characteristic image acquisition times and error distributions from Antarctic image pairs in the ITS_LIVE dataset, we generate synthetic ice velocity time series, then apply our method to recover imposed magnitudes of seasonal variability within ±1.4 m yr−1. We then validate the techniques by comparing our results to GPS data collected on Russell Glacier in Greenland. The methods presented here may be applied to better understand how ice dynamic signals propagate on seasonal timescales, and what mechanisms control the flow of the world’s ice.

Chad A. Greene↗
NASA NTRS1984

Numerical simulations of the Ross Sea tides

Tidal currents below the floating Ross Ice shelf are reconstructed by using a numerical tidal model. They are predominantly diurnal, achieve maximum strength in regions near where the ice shelf runs aground, and are significantly enhanced by topographic Rossby wave propagation along the ice front. A comparison with observations of the vertical motion of the ice shelf surface indicates that the model reproduces the diurnal tidal characteristics within 20 percent. Similar agreement for the relatively weak semidiurnal tides was not obtained, and this calls attention to possible errors of the open boundary forcing obtained from global-ocean tidal simulations and to possible errors in mapping zones of ice shelf grounding. Air-sea contact below the ice shelf is eliminated by the thick ice cover. The dominant sub-ice-shelf circulation may thus be tidally induced. A preliminary assessment of sub-ice-shelf conditions based on the numerical tidal simulations suggests that (1) strong barotropic circulation is driven along the ice front and (2) tidal fronts may form in the sub-ice-shelf cavity where the water column is thin and where the buoyancy input is weak.

Macayeal, D. R.↗
NASA NTRS1994

An alternative approach to estimating rainfall rate by radar using propagation differential phase shift

In this work it is shown that for frequencies from 3 to 13 GHz, the ratio of the specific propagation differential phase shift phi(sub DP) to the rainfall rate can be specified essentially independently of the form of the drop size distribution by a function only of the mass-weighted mean drop size D(sub m). This significantly reduces one source of substantial bias errors common to most other techniques for measuring rain by radar. For frequencies 9 GHz and greater, the coefficient can be well estimated from the ratio of the specific differential attenuation to phi(sub DP), while at nonattenuating frequencies such as 3 GHz, the coefficient can be well estimated using the differential reflectivity. In practice it appears that this approach yields better estimates of the rainfall rate than any other current technique. The best results are most likely at 13.80 GHz, followed by those at 2.80 GHz. An optimum radar system for measuring rain should probably include components at a both frequencies so that when signals at 13.8 GHz are lost because of attenuation, good measurements are still possible at the lower frequency.

Jameson, A. R.↗
NASA NTRS2003

First Year Wilkinson Microwave Anisotropy Probe(WMAP) Observations: Data Processing Methods and Systematic Errors Limits

We describe the calibration and data processing methods used to generate full-sky maps of the cosmic microwave background (CMB) from the first year of Wilkinson Microwave Anisotropy Probe (WMAP) observations. Detailed limits on residual systematic errors are assigned based largely on analyses of the flight data supplemented, where necessary, with results from ground tests. The data are calibrated in flight using the dipole modulation of the CMB due to the observatory's motion around the Sun. This constitutes a full-beam calibration source. An iterative algorithm simultaneously fits the time-ordered data to obtain calibration parameters and pixelized sky map temperatures. The noise properties are determined by analyzing the time-ordered data with this sky signal estimate subtracted. Based on this, we apply a pre-whitening filter to the time-ordered data to remove a low level of l/f noise. We infer and correct for a small (approx. 1 %) transmission imbalance between the two sky inputs to each differential radiometer, and we subtract a small sidelobe correction from the 23 GHz (K band) map prior to further analysis. No other systematic error corrections are applied to the data. Calibration and baseline artifacts, including the response to environmental perturbations, are negligible. Systematic uncertainties are comparable to statistical uncertainties in the characterization of the beam response. Both are accounted for in the covariance matrix of the window function and are propagated to uncertainties in the final power spectrum. We characterize the combined upper limits to residual systematic uncertainties through the pixel covariance matrix.

Hinshaw, G.↗
NASA NTRS2009

Error-Rate Bounds for Coded PPM on a Poisson Channel

Equations for computing tight bounds on error rates for coded pulse-position modulation (PPM) on a Poisson channel at high signal-to-noise ratio have been derived. These equations and elements of the underlying theory are expected to be especially useful in designing codes for PPM optical communication systems. The equations and the underlying theory apply, more specifically, to a case in which a) At the transmitter, a linear outer code is concatenated with an inner code that includes an accumulator and a bit-to-PPM-symbol mapping (see figure) [this concatenation is known in the art as "accumulate-PPM" (abbreviated "APPM")]; b) The transmitted signal propagates on a memoryless binary-input Poisson channel; and c) At the receiver, near-maximum-likelihood (ML) decoding is effected through an iterative process. Such a coding/modulation/decoding scheme is a variation on the concept of turbo codes, which have complex structures, such that an exact analytical expression for the performance of a particular code is intractable. However, techniques for accurately estimating the performances of turbo codes have been developed. The performance of a typical turbo code includes (1) a "waterfall" region consisting of a steep decrease of error rate with increasing signal-to-noise ratio (SNR) at low to moderate SNR, and (2) an "error floor" region with a less steep decrease of error rate with increasing SNR at moderate to high SNR. The techniques used heretofore for estimating performance in the waterfall region have differed from those used for estimating performance in the error-floor region. For coded PPM, prior to the present derivations, equations for accurate prediction of the performance of coded PPM at high SNR did not exist, so that it was necessary to resort to time-consuming simulations in order to make such predictions. The present derivation makes it unnecessary to perform such time-consuming simulations.

Moision, Bruce↗
NASA NTRS1991

Phase retrieval for the Hubble Space Telescope using iterative propagation algorithms

Phase retrieval algorithms, including the iterative transform algorithm and gradient search algorithms, were generalized to include the effects of propagation through a complicated optical system and to discount the effects of bad CCD pixels. For the gradient search algorithms, analytic gradients were derived that greatly speed up the computation over finite difference methods. For the Hubble Space Telescope (HST), the aperture function was reconstructed and the phase errors were retrieved. This information is useful to design correction optics for the telescope and for the deconvolution of blurred images from the HST.

Fienup, J. R.↗
NASA NTRS1994

A new approach to estimating rainwater content by radar using propagation differential phase shift

As microwaves propagate through rain, the rate of phase change with increasing distance is different depending upon whether the transmissions are polarized horizontally or vertically. This rate of change is the so-called specific propagation differential phase shift phi(sub DP). This paper demonstrates that at several frequencies and over a wide domain the ratio of phi(sub DP) to the rainwater content W is nearly linearly related to D(sub m), the mass-weighted mean drop size. An investigation of errors indicates that this new approach is likely to yield more accurate estimates of W than the other classical reflectivity factor Z, attenuation, or polarization techniques. The most accurate estimates of W are most likely at the highest frequency considered, 13.80 GHz. In lieu of such high-frequency measurements, these somewhat esoteric results are made more concrete through an analysis of 3-GHz radar measurements collected during the Convection and Precipitation Experiment in a tropical rainstorm in Florida. Among the principal advantages of using phi(sub DP) to measure rain are that an absolute calibration of the radar is no longer required and the estimates are decoupled from measurements of the radar reflectivity factor. Consequently, temporal and spatial structures of rain estimates do not simply mimic those of the reflectivity factor, as happens for classical estimation techniques using Z.

Jameson, A. R.↗
NASA NTRS1990

Inflight estimation of gyro noise

A method is described and demonstrated for estimating single-axis gyro noise levels in terms of the Farrenkopf model parameters. This is accomplished for the Cosmic Background Explorer (COBE) by comparing gyro-propagated attitudes with less accurate single-frame solutions and fitting the squared differences to a third-order polynomial in time. Initial results are consistent with the gyro specifications, and these results are used to determine limits on the duration of batches used to determine attitude. Sources of error are discussed, and guidelines for a more elegant implementation, as part of a batch estimator or filter, are included for future work.

Filla, O. H.↗
NASA NTRS1993

A ground track control algorithm for the Topographic Mapping Laser Altimeter (TMLA)

The results of an analysis of an algorithm that will provide autonomous onboard orbit control using orbits determined with Global Positioning System (GPS) data. The algorithm uses the GPS data to (1) compute the ground track error relative to a fixed longitude grid, and (2) determine the altitude adjustment required to correct the longitude error. A program was written on a personal computer (PC) to test the concept for numerous altitudes and values of solar flux using a simplified orbit model including only the J sub 2 zonal harmonic and simple orbit decay computations. The algorithm was then implemented in a precision orbit propagation program having a full range of perturbations. The analysis showed that, even with all perturbations (including actual time histories of solar flux variation), the algorithm could effectively control the spacecraft ground track and yield more than 99 percent Earth coverage in the time required to complete one coverage cycle on the fixed grid (220 to 230 days depending on altitude and overlap allowance).

Blaes, V.↗
NASA NTRS2000

Neural Network and Response Surface Methodology for Rocket Engine Component Optimization

The goal of this work is to compare the performance of response surface methodology (RSM) and two types of neural networks (NN) to aid preliminary design of two rocket engine components. A data set of 45 training points and 20 test points obtained from a semi-empirical model based on three design variables is used for a shear coaxial injector element. Data for supersonic turbine design is based on six design variables, 76 training, data and 18 test data obtained from simplified aerodynamic analysis. Several RS and NN are first constructed using the training data. The test data are then employed to select the best RS or NN. Quadratic and cubic response surfaces. radial basis neural network (RBNN) and back-propagation neural network (BPNN) are compared. Two-layered RBNN are generated using two different training algorithms, namely solverbe and solverb. A two layered BPNN is generated with Tan-Sigmoid transfer function. Various issues related to the training of the neural networks are addressed including number of neurons, error goals, spread constants and the accuracy of different models in representing the design space. A search for the optimum design is carried out using a standard gradient-based optimization algorithm over the response surfaces represented by the polynomials and trained neural networks. Usually a cubic polynominal performs better than the quadratic polynomial but exceptions have been noticed. Among the NN choices, the RBNN designed using solverb yields more consistent performance for both engine components considered. The training of RBNN is easier as it requires linear regression. This coupled with the consistency in performance promise the possibility of it being used as an optimization strategy for engineering design problems.

Vaidyanathan, Rajkumar↗
NASA NTRS2000

Seven Years of ACTS Technology Verification Experiments Reviewed

The Advanced Communications Technology Satellite (ACTS) was designed to achieve a 99.5-percent system availability rate and signals with less than one error in 10(exp 7) bits throughout the continental United States. To accomplish such a high rate of system availability, ACTS uses multiple narrow hopping beams and very small aperture terminal (VSAT) technology. In addition, ACTS uses an adaptive rain fade compensation protocol to reduce the negative effects of propagation on the system. To enhance knowledge on how propagation and system variances affect system availability, researchers at the NASA Glenn Research Center at Lewis Field performed technology verification experiments over a 7-yr period (from September 1993 to the present). These experiments include T1VSAT System Availability, Statistical Rain Fade Compensation Characterization, Statistical Characterization of Ka-Band Propagation Effects on Communication Link Performance, and Multibeam Antenna Performance.

Acosta, Roberto J.↗
NASA NTRS1983

LANDSAT 4 band 6 data evaluation

Satellite data collected over Lake Ontario were processed to observed surface temperature values. This involved computing apparent radiance values for each point where surface temperatures were known from averaged digital count values. These radiance values were then converted by using the LOWTRAN 5A atmospheric propagation model. This model was modified by incorporating a spectral response function for the LANDSAT band 6 sensors. A downwelled radiance term derived from LOWTRAN was included to account for reflected sky radiance. A blackbody equivalent source radiance was computed. Measured temperatures were plotted against the predicted temperature. The RMS error between the data sets is 0.51K.

Source record↗
NASA NTRS1986

Performance analysis of replication ALOHA for fading mobile communications channels

This paper describes an ALOHA random access protocol for fading communications channels. A two-state Markov model is used for the channel error process to account for the channel fading memory. The ALOHA protocol is modified to send multiple contiguous copies of a message at each transmission attempt. Both pure and slotted ALOHA channels are considered. The analysis is applicable to fading environments where the channel memory is short compared to the propagation delay. It is shown that smaller delay may be achieved using replications and, in noisy conditions, can also improve throughput.

Yan, Tsun-Yee↗
NASA NTRS1991

A proposed time transfer experiment between the USA and the South Pacific

This paper describes the concept, architecture and preliminary details of an experiment directc<l towards providing continuous Ultra High Precision (UHP) time transfer between Washington, DC; Salisbury, SA Australia; Orroral Valley, ACT Australia; arid Lower Hutt, New Zealand. It further describes a proposed method of distributing UTC(USNO) at a high level of precision to passive users over a broad area of the South Pacific. The concept is based on active two-way satellite time transfer from the United States Naval Observatory (USNO) to the proposed USNO Master Clock West (MCW) in Wahiwa, HI USA at the 1 nanosecond level using active satellite two-way time transfer augmented by Precise Positioning Service (PPS) of the Global Positioning System (GPS). MCW would act as an intermediate transfer/reference station, again linked to Salisbury at the 1 nanosecond level using active satellite two-way time transfer augmented by PPS GPS. From this point, time would bc distributed within the region by two methods. The first is an existing TV line sync system using an Australian communications satellite (AUSSAT K1) which is useful to the 20 nanosecond level. The second approach is RF ranging and multilateration between Salisbury, Orroral Observatory, Lower Hutt and the AUSSAT B1 and B2 to be launched in 1992. Orroral Observatory will provide precise laser ranging to the AUSSAT B1/B2 retro reflectors which will reduce ephemeris related time transfer errors to below 1 nanosecond. The corrected position will be transmitted by both the time transfer modem and the existing TV line sync dissemination process. Multilateration has the advantage of being an all weather approach and when used with the laser ranging technique will provide a precise measurement of the propagation path delays. This will result in time transfer performance levels on the order of 10 nanoseconds to passive users in both Australia and New Zealand.

John Luck↗
NASA NTRS1990

Asymptotic boundary conditions for dissipative waves: General theory

An outstanding issue in the computational analysis of time dependent problems is the imposition of appropriate radiation boundary conditions at artificial boundaries. Accurate conditions are developed which are based on the asymptotic analysis of wave propagation over long ranges. Employing the method of steepest descents, dominant wave groups are identified and simple approximations to the dispersion relation are considered in order to derive local boundary operators. The existence of a small number of dominant wave groups may be expected for systems with dissipation. Estimates of the error as a function of domain size are derived under general hypotheses, leading to convergence results. Some practical aspects of the numerical construction of the asymptotic boundary operators are also discussed.

Hagstrom, Thomas↗