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

An efficient nonlinear relaxation technique for the three-dimensional, Reynolds-averaged Navier-Stokes equations

An efficient implicit method for the computation of steady, three-dimensional, compressible Navier-Stokes flowfields is presented. A nonlinear iteration strategy based on planar Gauss-Seidel sweeps is used to drive the solution toward a steady state, with approximate factorization errors within a crossflow plane reduced by the application of a quasi-Newton technique. A hybrid discretization approach is employed, with flux-vector splitting utilized in the streamwise direction and central differences with artificial dissipation used for the transverse fluxes. Convergence histories and comparisons with experimental data are presented for several 3-D shock-boundary layer interactions. Both laminar and turbulent cases are considered, with turbulent closure provided by a modification of the Baldwin-Barth one-equation model. For the problems considered (175,000-325,000 mesh points), the algorithm provides steady-state convergence in 900-2000 CPU seconds on a single processor of a Cray Y-MP.

Edwards, Jack R.↗

Clinical benefits of new immobilization system for hypofractionated radiotherapy of intrahepatic hepatocellular carcinoma by helical tomotherapy

Objective: A comprehensive clinical evaluation was conducted, assessing the Body Pro-Lok immobilization and positioning system to facilitate hypofractionated radiotherapy of intrahepatic hepatocellular carcinoma (HCC), using helical tomotherapy to improve treatment precision. Methods: Clinical applications of the Body Pro-Lok system were investigated (as above) in terms of interfractional and intrafractional setup errors and compressive abdominal breath control. To assess interfractional setup errors, a total of 42 patients who were given 5 to 20 fractions of helical tomotherapy for intrahepatic HCC were analyzed. Overall, 15 patients were immobilized using simple vacuum cushion (group A), and the Body Pro-Lok system was used in 27 patients (group B), performing megavoltage computed tomography (MVCT) scans 196 times and 435 times, respectively. Pretreatment MVCT scans were registered to the planning kilovoltage computed tomography (KVCT) for error determination, and group comparisons were made. To establish intrafractional setup errors, 17 patients with intrahepatic HCC were selected at random for immobilization by Body Pro-Lok system, undergoing MVCT scans after helical tomotherapy every week. A total of 46 MVCT re-scans were analyzed for this purpose. In researching breath control, 12 patients, randomly selected, were immobilized by Body Pro-Lok system and subjected to 2-phase 4-dimensional CT (4DCT) scans, with compressive abdominal control or in freely breathing states, respectively. Respiratory-induced liver motion was then compared. Results: Mean interfractional setup errors were as follows: (1) group A: X, 2.97 ± 2.47 mm; Y, 4.85 ± 4.04 mm; and Z, 3.77 ± 3.21 mm; pitch, 0.66 ± 0.62°; roll, 1.09 ± 1.06°; and yaw, 0.85 ± 0.82°; and (2) group B: X, 2.23 ± 1.79 mm; Y, 4.10 ± 3.36 mm; and Z, 1.67 ± 1.91 mm; pitch, 0.45 ± 0.38°; roll, 0.77 ± 0.63°; and yaw, 0.52 ± 0.49°. Between-group differences were statistically significant in 6 directions (p < 0.05). Mean intrafractional setup errors with use of the Body Pro-Lok system were as follows: X, 0.41 ± 0.46 mm; Y, 0.86 ± 0.80 mm; Z, 0.33 ± 0.44 mm; and roll, 0.12 ± 0.19°. Mean liver-induced respiratory motion determinations were as follows: (1) abdominal compression: X, 2.33 ± 1.22 mm; Y, 5.11 ± 2.05 mm; Z, 2.13 ± 1.05 mm; and 3D vector, 6.22 ± 1.94 mm; and (2) free breathing: X, 3.48 ± 1.14 mm; Y, 9.83 ± 3.00 mm; Z, 3.38 ± 1.59 mm; and 3D vector, 11.07 ± 3.16 mm. Between-group differences were statistically different in 4 directions (p < 0.05). Conclusions: The Body Pro-Lok system is capable of improving interfractional and intrafractional setup accuracy and minimizing tumor movement owing to respirations in patients with intrahepatic HCC during hypofractionated helical tomotherapy.

61 RADIATION PROTECTION AND DOSIMETRY↗

Terrain Relative Navigation for Guided Descent on Titan

Titan’s dense atmosphere, low gravity, and high winds at high altitudes create descent times of >90 minutes with standard entry/descent/landing (EDL) architectures and result in large unguided landing ellipses, with 99% values of 110x110 km and 149x72 km in recent Titan lander proposals. Enabling precision landing on Titan could increase science return for the types of missions proposed to date and make additional types of landing sites accessible, opening up new possibilities for science investigations. Precision landing on Titan has unique challenges, because the hazy atmosphere makes it difficult to see the surface and because it requires guided descent with divert ranges that are one to two orders of magnitude larger than needed for other target bodies, i.e. up to on the order of 100 km. It is conceivable that such a divert capability could be provided economically by a parafoil or other steerable aerodynamic decelerator deployed several 10s of km above the surface. The long descent times lead to large inertial navigation errors, hence a need for terrain relative navigation (TRN). This would require a TRN capability that can operate at such altitudes, despite challenges of seeing the surface sufficiently clearly and of depending on map products that are two orders of magnitude lower in spatial resolution than those for Mars and airless bodies. We then develop algorithms for map matching and feature tracking with descent images and test these with synthetic images created from Cassini/Huygens data sets and our radiative transfer model. We also introduce new possibilities for TRN based on the potential to discriminate some specific types of terrain onboard in descent imagery, such as lake vs adjacent ground and dune vs interdune. We use sensor measurement noise models in simulations of state estimation with an extended Kalman filter that includes coordinates of a set of tracked features in the state vector. Case studies were done for two notional landing sites, one in a site with only dry ground and one in a Titan lake district. In both cases, the filter error model shows 3 position error at touchdown on the order of 2 km. More work is needed to validate these results with higher fidelity camera models and larger data sets, but this is very promising.

Matthies, Larry↗

Neutrino neutral-current elastic scattering on {sup 12}C

The neutral-current elastic scattering of neutrinos on carbon and CH{sub 2} targets is computed using the relativistic distorted-wave impulse approximation with a relativistic optical potential. Results for exclusive and inclusive neutrino reactions on a {sup 12}C target are presented. We show that the nuclear effects on the shape of the four-momentum transferred squared distribution d{sigma}/dQ{sub QE}{sup 2} in neutrino neutral-current and charged-current quasielastic scattering are similar. We also calculate the flux-averaged neutral-current elastic differential cross section d{sigma}/dQ{sub QE}{sup 2} for neutrino scattering from CH{sub 2}, as well as the neutral-current to charged-current cross section ratio as functions of Q{sub QE}{sup 2}. The axial mass M{sub A} is extracted from a fit of the d{sigma}/dQ{sub QE}{sup 2} cross section measured in the MiniBooNE experiment. The extracted value of M{sub A}=1.28{+-}0.05 GeV is consistent within errors with the MiniBooNE result. Additionally, for proton kinetic energies above the Cherenkov threshold, the strange-quark contribution to the neutral-current axial vector form factor at Q{sub QE}{sup 2}=0, {Delta}s, was extracted from a fit to MiniBooNE data for the {nu}p{yields}{nu}p to {nu}N{yields}{nu}N cross section ratio. This value is found to be {Delta}s=-0.11{+-}0.36

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗

Preliminary Theoretical Analysis of Mixed Precision Krylov Subspace Methods (Q3 Report)

The third quarter of the project was spent performing theoretical finite precision analysis of Krylov subspace method variants that use mixed precision. Our focus here is on the Conjugate Gradient (CG) method and the Lanczos method. We have performed an analysis of maximum attainable accuracy for the classical CG method in which 3 precisions are used: a working precision ε, a precision ε IP for the inner product computations, and a precision ε MV for the matrix-vector products. Our results show that performing inner product computations in lower precision does not affect the attainable accuracy. Further, we have performed a complete error analysis of the s-step Lanczos algorithm. In this case, we show that the numerical behavior of the algorithm can be significantly improved by using extra precision in a small part of the computation. We summarize the main theorems in the remainder of the document. Other activities include attending biweekly xSDK meetings. The subsequent quarter will be spent finalizing these results into technical reports and/or manuscripts for submission to journals, as well as identifying opportunities for future work.

97 MATHEMATICS AND COMPUTING↗

Magnetic field of Jupiter: A generalized inverse approach

The estimation of planetary magnetic fields from observations of the magnetic field gathered along a spacecraft flyby trajectory is examined with the aid of generalized inverse techniques, with application to the internal magnetic field of Jupiter. Model nonuniqueness resulting from the limited spatial extent of the observations and noise on the data is explored and quantitative estimates of the model parameter resolution are found. The presence of a substantial magnetic field of external origin due to the currents flowing in the Jovian magnetodisc is found to be an important source of error in estimates of the internal Jovian field, and new models explicitly incorporating these currents are proposed. New internal field models are derived using the vector helium magnetometer observations and the high field fluxgate observations of Pioneer 11, and knowledge of the external current system gained from the Pioneer 10 and Voyagers 1 and 2 encounters.

Connerney, J. E. P.↗

The magnetic field of Jupiter - A generalized inverse approach

The estimation of planetary magnetic fields from observations of the magnetic field gathered along a spacecraft flyby trajectory is examined with the aid of generalized inverse techniques, with application to the internal magnetic field of Jupiter. Model non-uniqueness resulting from the limited spatial extent of the observations and noise on the data is explored and quantitative estimates of the model parameter resolution are found. The presence of a substantial magnetic field of external origin due to the currents flowing in the Jovian magnetodisc is found to be an important source of error in estimates of the internal Jovian field, and new models explicitly incorporating these currents are proposed. New internal field models are derived using the vector helium magnetometer observations and the high field fluxgate observations of Pioneer 11, and knowledge of the external current system gained from the Pioneer 10 and Voyagers 1 and 2 encounters.

Connerney, J. E. P.↗

MAGSAT and aeromagnetic data in the North American continent

Problems were encountered in deriving a proper reference field to be subtracted from the aeromagnetic data obtained from Project MAGNET. Field models tried thus far do not seem to eliminate properly the main field. The MAGSAT data in the North American continent for the period November 1 to December 22, 1979 are being compiled and compared with MAGNET data. Efforts are being made to eliminate the orbital bias errors. A computer program was developed and successfully tested which computes a topographic profile of the Curie depth isotherm which fits best to the observed vector or scalar field magnetic data.

Source record↗

Laser Doppler anemometer measurements using nonorthogonal velocity components - Error estimates

Laser Doppler anemometers (LDAs) that are arranged to measure nonorthogonal velocity components (from which orthogonal components are computed through transformation equations) are more susceptible to calibration and sampling errors than are systems with uncoupled channels. In this paper uncertainty methods and estimation theory are used to evaluate, respectively, the systematic and statistical errors that are present when such devices are applied to the measurement of mean velocities in turbulent flows. Statistical errors are estimated for two-channel LDA data that are either correlated or uncorrelated. For uncorrelated data the directional uncertainty of the measured velocity vector is considered for applications where mean streamline patterns are desired.

Orloff, K. L.↗

Reduction of flow-measurement uncertainties in laser velocimeters with nonorthogonal channels

An analysis of certain geometrical limitations inherent in the application of laser velocimeters with nonorthogonal channels has led to the development of advanced-LDA-calibration and data-acquisition techniques that minimize systematic and statistical errors, respectively. The data-acquisition technique optimizes the number of velocity samples collected from three velocimeter channels as a function of local turbulence intensity, vector direction, and prescribed confidence interval. Linear velocity surveys and streamline traces measured in a turbulent flow field with a three-dimensional laser velocimeter are presented and the validity and accuracy of the theoretical analysis are discussed.

Snyder, P. K.↗

State estimation Kalman filter using optical processings Noise statistics known

Reference is made to a study by Casasent et al. (1983), which gave a description of a frequency-multiplexed acoustooptic processor and showed how it was capable of performing all the individual operations required in Kalman filtering. The data flow and organization of all required operations however, were not detailed in that study. Consideration is given here to a simpler Kalman filter state estimation problem. Equally spaced time-sampled intervals (k times T sub s, with k the iterative time index) are assumed. It is further assumed that the system noise vector w and the measurement noise vector v are uncorrelated and Gaussian distributed and that the noise statistics (Q and R) and the system model (Phi, Gamma, H) are known. The error covariance matrix P and the extrapolated error covariance matrix M can thus be precomputed and the Kalman gain matrix K sub k can be precomputed and stored for each input time sample.

Jackson, J.↗

Modal cost analysis for simple continua

The most popular finite element codes are based upon appealing theories of convergence of modal frequencies. For example, the popularity of cubic elements for beam-like structures is due to the rapid convergence of modal frequencies and stiffness properties. However, for those problems in which the primary consideration is the accuracy of response of the structure at specified locations, it is more important to obtain accuracy in the modal costs than in the modal frequencies. The modal cost represents the contribution of a mode in the norm of the response vector. This paper provides a complete modal cost analysis for simple continua such as beam-like structures. Upper bounds are developed for mode truncation errors in the model reduction process and modal cost analysis dictates which modes to retain in order to reduce the model for control design purposes.

Hu, A.↗

Tropical synoptic composites in VAS, OLR and TOVS radiance data

The horizontal structure of specific tropical synoptic features are compared with the features of a zonal mean ITCZ using four different forms of satellite radiance data: daily gridded OLR data, GOES VAS water vapor imagery, brightness temperatures from TOVS water vapor channels (6.7 and 7.3 microns), and brightness temperature vectors for a complete set of TOVS channels. All of the data are composited over many events to minimize the effects of missing data, sampling errors, and system variability. The major differences between the observing systems are examined, focusing on resolution and clouds. It is found that the ITCZ exhibits a zonally symmetric convective phase and an enhanced regional convective phase (plumes), both of which can be differentiated from climatological mean behavior.

Mcguirk, James P.↗

Multilevel Algorithm for Atmospheric Data Assimilation

A multiscale algorithm for the problem of optimal statistical interpolation of observed data has been developed. This problem includes the calculation of the vector of the 'analyzed' (best estimated) atmosphere flow field w(sup a) by the formula: w(sup a) = w(sup f) + P(sup f) H(sup T) y, where the quantity y is defined by the equation (H P(sup f) H(sup T) + R)y = w(sup o) - H w(sup f), using the given model forecast first guess w(sup f) and the vector of observations w(sup o); H is an interpolation operator from the regular grid to the observation network, P(sup f) is the forecast error covariance matrix, and R is the observation error covariance matrix. At this initial stage the case of univariate analysis of single level radiosonde height data is considered. The matrix R is assumed to be diagonal, and the matrix P(sup f) is assumed to be given by the formula P(sub ij)(sup f) = sigma(sub i)(sup f) mu(sub ij) sigma(sub j)(sub f), where mu(sub ij) is a smooth, decreasing function of the distance between the i-th and the j-th points. In this paper we describe a multiscale iterative process based on a multiresolution, simultaneous displacement technique and a localized variational calculation of iteration parameters.

Brandt, Achi↗

Improved Modeling in a Matlab-Based Navigation System

An innovative approach to autonomous navigation is available for low earth orbit satellites. The system is developed in Matlab and utilizes an Extended Kalman Filter (EKF) to estimate the attitude and trajectory based on spacecraft magnetometer and gyro data. Preliminary tests of the system with real spacecraft data from the Rossi X-Ray Timing Explorer Satellite (RXTE) indicate the existence of unmodeled errors in the magnetometer data. Incorporating into the EKF a statistical model that describes the colored component of the effective measurement of the magnetic field vector could improve the accuracy of the trajectory and attitude estimates and also improve the convergence time. This model is identified as a first order Markov process. With the addition of the model, the EKF attempts to identify the non-white components of the noise allowing for more accurate estimation of the original state vector, i.e. the orbital elements and the attitude. Working in Matlab allows for easy incorporation of new models into the EKF and the resulting navigation system is generic and can easily be applied to future missions resulting in an alternative in onboard or ground-based navigation.

Deutschmann, Julie↗

The Influence of Antenna Pattern on Faraday Rotation in Remote Sensing at L-band

Faraday rotation is a change in the polarization vector of electromagnetic radiation that occurs as the waves propagate from the Earth surface through the ionosphere to a spaceborne sensor. This change can cause errors in monitoring parameters at the surface such as soil moisture and sea surface salinity and it is an important consideration for radiometers on future missions in space such as NASA's Aquarius mission and ESA's SMOS mission. Two prominent strategies for compensating for Faraday rotation are using a sum of the signal at two polarizations and using the correlation between the signals at the two polarizations. These strategies work for an idealized antenna. This paper evaluates the strategies in the context of realistic antennas such as will be built for the Aquarius radiometer. Realistic antennas will make small differences that need to be included in planning for retrieval algorithms in future missions.

LeVine, David M.↗

Satellite Data Assimilation in the GEOS Atmospheric Data Assimilation System

The Goddard Earth Observing System (GEOS) is used for weather, climate, and air quality forecasts and producing reanalysis datasets. Its atmospheric data assimilation system (ADAS) is a hybrid–4DEnVar system. Various satellite radiance, atmospheric motion vector (winds), aircraft, and convectional data are assimilated by this system to produce global atmospheric states. The GEOS ADAS is enhanced by considering correlated observational errors among infrared radiance observations and assimilating new microwave radiance data in all-sky conditions. Historical microwave radiance data are tested within the all-sky microwave radiance assimilation framework and show beneficial effect in the preparation of the next GEOS reanalysis product. In this seminar, I will present the procedures of to assimilate microwave radiance data in all-sky conditions in the GEOS-ADAS. I will also discuss our efforts to update GEOS ADAS by assimilating more satellite data such as geostationary infrared observations and new observations made by the Advanced Technology Microwave Sounder (ATMS) aboard NOAA 21 polar orbit satellite.

Jianjun Jin↗

TH-CD-207A-05: Lung Surface Deformation Vector Fields Prediction by Monitoring Respiratory Surrogate Signals

Purpose: In this study, we developed and evaluated a method for predicting lung surface deformation vector fields (SDVFs) based on surrogate signals such as chest and abdomen motion at selected locations and spirometry measurements. Methods: A Patient-specific 3D triangular surface mesh of the lung region at end-expiration (EE) phase was obtained by threshold-based segmentation method. For each patient, a spirometer recorded the flow volume changes of the lungs; and 192 selected points at a regular spacing of 2cm X 2cm matrix points over a total area of 34cm X 24cm on the surface of chest and abdomen was used to detect chest wall motions. Preprocessing techniques such as QR factorization with column pivoting (QRCP) were employed to remove redundant observations of the chest and abdominal area. To create a statistical model between the lung surface and the corresponding surrogate signals, we developed a predictive model based on canonical ridge regression (CRR). Two unique weighting vectors were selected for each vertex on the surface of the lung, and they were optimized during the training process using the all other phases of 4D-CT except the end-inspiration (EI) phase. These parameters were employed to predict the vertices locations of a testing data set, which was the EI phase of 4D-CT. Results: For ten lung cancer patients, the deformation vector field of each vertex of lung surface mesh was estimated from the external motion at selected positions on the chest wall surface plus spirometry measurements. The average estimation of 98th percentile of error was less than 1 mm (AP= 0.85, RL= 0.61, and SI= 0.82). Conclusion: The developed predictive model provides a non-invasive approach to derive lung boundary condition. Together with personalized biomechanical respiration modelling, the proposed model can be used to derive the lung tumor motion during radiation therapy accurately from non-invasive measurements.

60 APPLIED LIFE SCIENCES↗