Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “error vector”

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 145 records · Page 8

Estimation of attitude sensor timetag biases

This paper presents an extended Kalman filter for estimating attitude sensor timing errors. Spacecraft attitude is determined by finding the mean rotation from a set of reference vectors in inertial space to the corresponding observed vectors in the body frame. Any timing errors in the observations can lead to attitude errors if either the spacecraft is rotating or the reference vectors themselves vary with time. The state vector here consists of the attitude quaternion, timetag biases, and, optionally, gyro drift rate biases. The filter models the timetags as random walk processes: their expectation values propagate as constants and white noise contributes to their covariance. Thus, this filter is applicable to cases where the true timing errors are constant or slowly varying. The observability of the state vector is studied first through an examination of the algebraic observability condition and then through several examples with simulated star tracker timing errors. The examples use both simulated and actual flight data from the Extreme Ultraviolet Explorer (EUVE). The flight data come from times when EUVE had a constant rotation rate, while the simulated data feature large angle attitude maneuvers. The tests include cases with timetag errors on one or two sensors, both constant and time-varying, and with and without gyro bias errors. Due to EUVE's sensor geometry, the observability of the state vector is severely limited when the spacecraft rotation rate is constant. In the absence of attitude maneuvers, the state elements are highly correlated, and the state estimate is unreliable. The estimates are particularly sensitive to filter mistuning in this case. The EUVE geometry, though, is a degenerate case having coplanar sensors and rotation vector. Observability is much improved and the filter performs well when the rate is either varying or noncoplanar with the sensors, as during a slew. Even with bad geometry and constant rates, if gyro biases are independently known, the timetag error for a single sensor can be accurately estimated as long as its boresight is not too close to the spacecraft rotation axis.

Sedlak, J.↗

Smart-Divert Powered Descent Guidance to Avoid the Backshell Landing Dispersion Ellipse

A smart-divert capability has been added into the Powered Descent Guidance (PDG) software originally developed for Mars pinpoint and precision landing. The smart-divert algorithm accounts for the landing dispersions of the entry backshell, which separates from the lander vehicle at the end of the parachute descent phase and prior to powered descent. The smart-divert PDG algorithm utilizes the onboard fuel and vehicle thrust vectoring to mitigate landing error in an intelligent way: ensuring that the lander touches down with minimum- fuel usage at the minimum distance from the desired landing location that also avoids impact by the descending backshell. The smart-divert PDG software implements a computationally efficient, convex formulation of the powered-descent guidance problem to provide pinpoint or precision-landing guidance solutions that are fuel-optimal and satisfy physical thrust bound and pointing constraints, as well as position and speed constraints. The initial smart-divert implementation enforced a lateral-divert corridor parallel to the ground velocity vector; this was based on guidance requirements for MSL (Mars Science Laboratory) landings. This initial method was overly conservative since the divert corridor was infinite in the down-range direction despite the backshell landing inside a calculable dispersion ellipse. Basing the divert constraint instead on a local tangent to the backshell dispersion ellipse in the direction of the desired landing site provides a far less conservative constraint. The resulting enhanced smart-divert PDG algorithm avoids impact with the descending backshell and has reduced conservatism.

Carson, John M.↗

The Z3 model of Saturns magnetic field and the Pioneer 11 vector helium magnetometer observations

Magnetic field observations obtained by the Pioneer 11 vector helium magnetometer are compared with the Z(sub 3) model magnetic field. These Pioneer 11 observations, obtained at close-in radial distances, constitute an important and independent test of the Z(sub 3) zonal harmonic model, which was derived from Voyager 1 and Voyager 2 fluxgate magnetometer observations. Differences between the Pioneer 11 magnetometer and the Z(sub 3) model field are found to be small (approximately 1%) and quantitatively consistent with the expected instrumental accuracy. A detailed examination of these differences in spacecraft payload coordinates shows that they are uniquely associated with the instrument frame of reference and operation. A much improved fit to the Pioneer 11 observations is obtained by rotation of the instrument coordinate system about the spacecraft spin axis by 1.4 degree. With this adjustment, possibly associated with an instrumental phase lag or roll attitude error, the Pioneer 11 vector helium magnetometer observations are fully consistent with the Voyager Z(sub 3) model.

Connerney, J. E. P.↗

The Z3 model of Saturn's magnetic field and the Pioneer 11 vector helium magnetometer observations

Magnetic field observations obtained by the Pioneer 11 vector helium magnetometer are compared with the Z(sub 3) model magnetic field. These Pioneer 11 observations, obtained at close-in radial distances, constitute an important and independent test of the Z(sub 3) zonal harmonic model, which was derived from Voyager 1 and Voyager 2 fluxgate magnetometer observations. Differences between the Pioneer 11 magnetometer and the Z(sub 3) model field are found to be small (approximately 1 percent) and quantitatively consistent with the expected instrumental accuracy. A detailed examination of these differences in spacecraft payload coordinates shows that they are uniquely associated with the instrument frame of reference and operation. A much improved fit to the Pioneer 11 observations is obtained by rotation of the instrument coordinate system about the spacecraft spin axis by 1.4 degree. With this adjustment, possibly associated with an instrumental phase lag or roll attitude error, the Pioneer 11 vector helium magnetometer observations are fully consistent with the Voyager Z(sub 3) model.

Connerney, J. E. P.↗

Biweekly Maps of Wind Stress for the North Pacific from the ERS-1 Scatterometer

The European Remote-sensing Satellite (ERS-1) was launched in July 1991 and contained several instruments for observing the Earth's ocean including a wind scatterometer. The scatterometer measurements were processed by the European Space Agency (ESA) and the Jet Propulsion Laboratory (JPL). JPL reprocessed (Freilich and Dunbar, 1992) the ERS-1 backscatter measurements to produced a 'value added' data set that contained the ESA wind vector as well as a set of up to four ambiguities. These ambiguities were further processed using a maximum-likelihood estimation (MLE) and a median filter to produce a 'selected vector.' This report describes a technique developed to produce time-averaged wind field estimates with their expected errors using only scatterometer wind vectors. The processing described in this report involved extracting regions of interest from the data tapes, checking the quality and creating the wind field estimate. This analysis also includes the derivation of biweekly average wind vectors over the North Pacific Ocean at a resolution of 0.50 x 0.50. This was done with an optimal average algorithm temporally and an over-determined biharmonic spline spatially. There have been other attempts at creating gridded wind files from ERS-1 winds, e.g., kriging techniques (Bentamy et al., 1996) and successive corrections schemes (Tang and Liu, 1996). There are several inherent problems with the ERS-1 scatterometer. Since this is a multidisciplinary mission, the satellite is flown in different orbits optimized for each phase of the mission. The scatterometer also shares several sub-systems with the Synthetic Aperture Radar (SAR) and cannot be operated while the SAR is in operation. The scatterometer is also a single-sided instrument and only measures backscatter along the right side of the satellite. The processing described here generates biweekly wind maps during the wktwo years analysis period regardless of the satellite orbit or missing data.

SYNTHETIC APERTURE RADAR↗

$B \to \pi \ell \nu$ and $B_s \to K \ell \nu$ form factors and $|V_{ub}|$ from 2+1-flavor lattice QCD with domain-wall light quarks and relativistic heavy quarks

We calculate the form factors for B → πℓν and B s → Kℓν decay in dynamical lattice quantum chromodynamics (QCD) using domain-wall light quarks and relativistic b-quarks. We use the (2+1)-flavor gauge-field ensembles generated by the RBC and UKQCD collaborations with the domain-wall fermion action and Iwasaki gauge action. For the b-quarks we use the anisotropic clover action with a relativistic heavy-quark interpretation. We analyze data at two lattice spacings of a ≈ 0.11, 0.086 fm with unitary pion masses as light as M π ≈ 290 MeV. We simultaneously extrapolate our numerical results to the physical light-quark masses and to the continuum and interpolate in the pion/kaon energy using SU(2) “hard-pion” chiral perturbation theory for heavy-light meson form factors. We provide complete systematic error budgets for the vector and scalar form factors f + (q 2 ) and f 0 (q 2 ) for both B → πℓν and B s → Kℓν at three momenta that span the q 2 range accessible in our numerical simulations. Next we extrapolate these results to q 2 = 0 using a model-independent z-parametrization based on analyticity and unitarity. We present our final results for f + (q 2 ) and f 0 (q 2 )as the coefficients of the series in z and the matrix of correlations between them; this provides a parametrization of the form factors valid over the entire allowed kinematic range. Our results agree with other three-flavor lattice-QCD determinations using staggered light quarks, and have comparable precision, thereby providing important independent cross-checks. Both B → πℓν and B s → Kℓν decays enable determinations of the Cabibbo-Kobayashi-Maskawa matrix element |V ub |. Furthermore, we perform a combined z-fit of our numerical B → πℓν form-factor data with the experimental measurements of the branching fraction from BABAR and Belle leaving the relative normalization as a free parameter; we obtain |V ub | = 3.61(32)×10 -3 , where the error includes statistical and all systematic uncertainties. The same approach can be applied to the decay B s → Kℓν to provide an alternative determination of |V ub | once the process has been measured experimentally. In anticipation of future experimental measurements, we make predictions for B → πℓν and B s → Kℓν differential branching fractions and forward-backward asymmetries in the Standard Model.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗

Structure Factors of Neutron Matter at Finite Temperature

Here, we compute continuum and infinite volume limit extrapolations of the structure factors of neutron matter at finite temperature and density. Using a lattice formulation of leading-order pionless effective field theory, we compute the momentum dependence of the structure factors at finite temperature and at densities beyond the reach of the virial expansion. The Tan contact parameter is computed and the result agrees with the high momentum tail of the vector structure factor. All errors, statistical and systematic, are controlled for. This calculation is a first step towards a model-independent understanding of the linear response of neutron matter at finite temperature.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

SU-G-TeP4-12: Individual Beam QA for a Robotic Radiosurgery System Using a Scintillator Cone

Purpose: The targeting accuracy of the Cyberknife system is measured by end-to-end tests delivering multiple isocentric beams to a point in space. While the targeting accuracy of two representative beams can be determined by a Winston-Lutz-type test, no test is available today to determine the targeting accuracy of each clinical beam. We used a scintillator cone to measure the accuracy of each individual beam. Methods: The XRV-124 from Logos Systems Int’l is a scintillator cone with an imaging system that is able to measure individual beam vectors and a resulting error between planned and measured beam coordinates. We measured the targeting accuracy of isocentric and non-isocentric beams for a number of test cases using the Iris and the fixed collimator. The average difference between plan and measured beam position was 0.8–1.2mm across the collimator sizes and plans considered here. The max error for a single beam was 2.5mm for the isocentric plans, and 1.67mm for the non-isocentric plans. The standard deviation of the differences was 0.5mm or less. Conclusion: The CyberKnife System is specified to have an overall targeting accuracy for static targets of less than 0.95mm. In E2E tests using the XRV124 system we measure average beam accuracy between 0.8 to 1.23mm, with maximum of 2.5mm. We plan to investigate correlations between beam position error and robot position, and to quantify the effect of beam position errors on patient specific plans. Martina Descovich has received research support and speaker honoraria from Accuray.

60 APPLIED LIFE SCIENCES↗

Adaptive Algebraic Multigrid for Finite Element Elliptic Equations with Random Coefficients

This thesis presents a two-grid algorithm based on Smoothed Aggregation Spectral Element Agglomeration Algebraic Multigrid (SA-ρAMGe) combined with adaptation. The aim is to build an efficient solver for the linear systems arising from discretization of second-order elliptic partial differential equations (PDEs) with stochastic coefficients. Examples include PDEs that model subsurface flow with random permeability field. During a Markov Chain Monte Carlo (MCMC) simulation process, that draws PDE coefficient samples from a certain distribution, the PDE coefficients change, hence the resulting linear systems to be solved change. At every such step the system (discretized PDE) needs to be solved and the computed solution used to evaluate some functional(s) of interest that then determine if the coefficient sample is acceptable or not. The MCMC process is hence computationally intensive and requires the solvers used to be efficient and fast. This fact that at every step of MCMC the resulting linear system changes, makes an already existing solver built for the old problem perhaps not as efficient for the problem corresponding to the new sampled coefficient. This motivates the main goal of our study, namely, to adapt an already existing solver to handle the problem (with changed coefficient) with the objective to achieve this goal to be faster and more efficient than building a completely new solver from scratch. Our approach utilizes the local element matrices (for the problem with changed coefficients) to build local problems associated with constructed by the method agglomerated elements (a set of subdomains that cover the given computational domain). We solve a generalized eigenproblem for each set in a subspace spanned by the previous local coarse space (used for the old solver) and a vector, component of the error, that the old solver cannot handle. A portion of the spectrum of these local eigen-problems (corresponding to eigenvalues close to zero) form the coarse basis used to define the new two-level method of our interest. We illustrate the performance of this adaptive two-level procedure with a large set of numerical experiments that demonstrate its efficiency over building the solvers from scratch.

97 MATHEMATICS AND COMPUTING↗

ML Enhancements via the Calculation of Rigid Body Modes (RBMs) for Mechanics Problems Implemented within the Albany Code Base

The algebraic multigrid approach known as smoothed aggregation is very efficient at solving systems that arise from elasticity problems. In order to construct an efficient algebraic multilevel method, a multigrid solver should be provided with a small set of vectors that represent the error components that are difficult to resolve. It is well-known that for linear elasticity problems, these components correspond to the the so-called rigid body modes (RBMs). The present document summarizes some new development within the Albany code base that has enabled the application of algebraic multigrid preconditioners from the ML package of Trilinos to mechanics problems implemented within Albany via a new function that calculates the RBMs using information about the problem's underlying mesh. The performance of these preconditioners is evaluated on four problems: a 3D static elasticity problem, a 3D non-linear elasticity problem, a 3D thermo-elasticity problem, and a 3D thermo-poro-plasticity problem. The tests reveal the superiority of the ML preconditioners over ILU preconditioners from the Trilinos Ifpack package for mechanics problems in Albany.

97 MATHEMATICS AND COMPUTING↗

Solar sensor with autocollimator

Monolithic structure, consisting of fused quartz fabricated solar sensor and autocollimator, measures angle between reference surface on spacecraft and solar vector. Device compensates for errors caused by changes in the alignment of its mountings or component parts. Integrated concept is advantageous whenever sensor mounting errors constitute problems.

Fain, D. L.↗

Stability of matrix equations arising in electromagnetics

Many problems in electromagnetics can be formulated as a Fredholm equation of the first kind. In these problems it is desired to determine some function f(x) from its integral transform g(y). Except for special cases for which an explicit inverse operator is known, the equation is usually discretized by the method of moments into a matrix equation Af = g, which is then solved on a computer. A method of calculating the sensitivity of the solution of Af = g due to errors in either the vector g or the matrix A is developed. The usefulness of this method is demonstrated with problems from antenna theory, scattering theory, and remote sensing.

Klein, C.↗

Use of sensitivity analysis to predict pilot performance as a function of different displays

A technique for objectively evaluating different displays by sensitivity analysis is described. First, the mathematical model used to analyze static displays is developed. The technique is based on formulating functional relationships between the state variables and the variables observable in the display. The matrix of the partial derivatives of the display variables with respect to the state variables, together with the observer's acuity function, is used to calculate expected errors in the state vector estimation. The technique is expanded by the use of Kalman filtering to process a time series of observation vectors. This provides a tool for analyzing displays of dynamic processes by means of a dynamic display evaluation computer program. Results are reported using this program to simulate an Instrument Landing System approach.

Mout, M. L.↗

Adaptive antenna arrays for weak interfering signals

The interference protection provided by adaptive antenna arrays to an Earth station or satellite receive antenna system is studied. The case where the interference is caused by the transmission from adjacent satellites or Earth stations whose signals inadverently enter the receiving system and interfere with the communication link is considered. Thus, the interfering signals are very weak. To increase the interference suppression, one can either decrease the thermal noise in the feedback loops or increase the gain of the auxiliary antennas in the interfering signal direction. Both methods are examined. It is shown that one may have to reduce the noise correlation to impractically low values and if directive auxiliary antennas are used, the auxiliary antenna size may have to be too large. One can, however, combine the two methods to achieve the specified interference suppression with reasonable requirements of noise decorrelation and auxiliary antenna size. Effects of the errors in the steering vector on the adaptive array performance are studied.

Gupta, I. J.↗

A study of optimal abstract jamming strategies vs. noncoherent MFSK

The present investigation is concerned with the performance of uncoded MFSK modulation in the presence of arbitrary additive jamming, taking into account the objective to devise robust antijamming strategies. An abstract model is considered, giving attention to the signal strength as a nonnegative real number X, the employment of X as a random variable, its distribution function G(x), the transmitter's strategy G, the jamming noise as an M-dimensional random vector Z, and the error probability. A summary of previous work on the considered problem is provided, and the results of the current study are presented.

Mceliece, R. J.↗

Spacecraft attitude determination using a second-order nonlinear filter

The stringent attitude determination accuracy and faster slew maneuver requirements demanded by present-day spacecraft control systems motivate the development of recursive nonlinear filters for attitude estimation. This paper presents the second-order filter development for the estimation of attitude quaternion using three-axis gyro and star tracker measurement data. Performance comparisons have been made by computer simulation of system models and filter mechanization. It is shown that the second-order filter consistently performs better than the extended Kalman filter when the performance index of the root sum square estimation error of the quaternion vector is compared. The second-order filter identifies the gyro drift rates faster than the extended Kalman filter. The uniqueness of this algorithm is the online generation of the time-varying process and measurement noise covariance matrices, derived as a function or the process and measurement nonlinearity, respectively.

Vathsal, S.↗

Principled halftoning based on human vision models

When models of human vision adequately measure the relative quality of candidate halftonings of an image, the problem of halftoning the image becomes equivalent to the search problem of finding a halftone that optimizes the quality metric. Because of the vast number of possible halftones, and the complexity of image quality measures, this principled approach has usually been put aside in favor of fast algorithms that seem to perform well. We find that the principled approach can lead to a range of useful halftoning algorithms, as we trade off speed for quality by varying the complexity of the quality measure and the thoroughness of the search. High quality halftones can be obtained reasonably quickly, for example, by using as a measure the vector length of the error image filtered by a contrast sensitivity function, and, as the search procedure, the sequential adjustment of individual pixels to improve the quality measure. If computational resources permit, simulated annealing can find nearly optimal solutions.

Mulligan, Jeffrey B.↗