Search NASA⌕ Search

SEARCH · Search NASA

Results for “DIRECTION VECTORS”

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 19 records

Autonomous frequency domain identification: Theory and experiment

The analysis, design, and on-orbit tuning of robust controllers require more information about the plant than simply a nominal estimate of the plant transfer function. Information is also required concerning the uncertainty in the nominal estimate, or more generally, the identification of a model set within which the true plant is known to lie. The identification methodology that was developed and experimentally demonstrated makes use of a simple but useful characterization of the model uncertainty based on the output error. This is a characterization of the additive uncertainty in the plant model, which has found considerable use in many robust control analysis and synthesis techniques. The identification process is initiated by a stochastic input u which is applied to the plant p giving rise to the output. Spectral estimation (h = P sub uy/P sub uu) is used as an estimate of p and the model order is estimated using the produce moment matrix (PMM) method. A parametric model unit direction vector p is then determined by curve fitting the spectral estimate to a rational transfer function. The additive uncertainty delta sub m = p - unit direction vector p is then estimated by the cross spectral estimate delta = P sub ue/P sub uu where e = y - unit direction vectory y is the output error, and unit direction vector y = unit direction vector pu is the computed output of the parametric model subjected to the actual input u. The experimental results demonstrate the curve fitting algorithm produces the reduced-order plant model which minimizes the additive uncertainty. The nominal transfer function estimate unit direction vector p and the estimate delta of the additive uncertainty delta sub m are subsequently available to be used for optimization of robust controller performance and stability.

Yam, Yeung↗

ICESat Spacecraft Pointing Support Study Grant

The Geodetic Laser Altimeter System (GLAS) mission is designed to measure changes in the elevations of the polar ice sheets. The ICESat satellite will carry the GLAS altimeter, and will have a nominal orbit altitude of 600 km and orbit inclination of 94deg. The groundtrack repeat period is 182 days and will be maintained to less than 1 km at the equator via routine orbit adjustments. Science requirements for the GLAS mission demand that the laser altimeter be pointed to within 50 meters of a predetermined reference groundtrack. As the actual ICESat groundtrack drifts away from the reference groundtrack, the attitude must be controlled such that the altimeter boresight is pointed, crosstrack, at the reference groundtrack. This orientation may be described by a rotation, theta, about the instantaneous geodetic local horizontal direction vector, which lies in the orbit plane and is oriented in the direction of motion of the satellite. The attitude is further complicated by requirements related to thermal and power considerations for various instruments, spacecraft components, and solar array orientation. In order to keep battery temperatures within the specified operating range, and maintain near normal pointing of the solar array with respect to the sunline direction vector as the orbit precesses relative to the sun, the satellite will be oriented in one of four fixed yaw modes. Each of these yaw modes depends upon the angle between the orbit plane and the sunline direction vector; this angle is designated Beta'. Table 1 shows the satellite yaw angle, Psi, for a given Beta' range. The angle Psi represents a rotation about the satellite z-axis, which points in the geodetic nadir direction; for Psi = 0deg the satellite x-axis points in the direction of motion.

Source record↗

Nested Krylov methods and preserving the orthogonality

Recently the GMRESR inner-outer iteraction scheme for the solution of linear systems of equations was proposed by Van der Vorst and Vuik. Similar methods have been proposed by Axelsson and Vassilevski and Saad (FGMRES). The outer iteration is GCR, which minimizes the residual over a given set of direction vectors. The inner iteration is GMRES, which at each step computes a new direction vector by approximately solving the residual equation. However, the optimality of the approximation over the space of outer search directions is ignored in the inner GMRES iteration. This leads to suboptimal corrections to the solution in the outer iteration, as components of the outer iteration directions may reenter in the inner iteration process. Therefore we propose to preserve the orthogonality relations of GCR in the inner GMRES iteration. This gives optimal corrections; however, it involves working with a singular, non-symmetric operator. We will discuss some important properties, and we will show by experiments that, in terms of matrix vector products, this modification (almost) always leads to better convergence. However, because we do more orthogonalizations, it does not always give an improved performance in CPU-time. Furthermore, we will discuss efficient implementations as well as the truncation possibilities of the outer GCR process. The experimental results indicate that for such methods it is advantageous to preserve the orthogonality in the inner iteration. Of course we can also use iteration schemes other than GMRES as the inner method; methods with short recurrences like GICGSTAB are of interest.

Desturler, Eric↗

Whistler mode turbulence generated by electron beams in earth's bow shock

The Landau and cyclotron growth rates of whistler mode waves in the earth's bow shock are calculated by using electron distribution functions obtained with the fast plasma experiment on ISEE 2. Three electron distribution functions measured within the transition region of the shock are analyzed. These functions spontaneously generate whistler mode waves with plasma rest frame frequencies between about 0.1 and 100 Hz. The wave normal angles of the generated whistlers range from 0 deg to the resonance cone angle. Electromagnetic Landau resonance and/or cyclotron resonance contribute to wave growth over the range of observed velocity distributions. Waves generated by the normal cyclotron resonance have wave vectors directed toward the solar wind, while those generated by the Landau and the anomalous cyclotron resonances have vectors directed toward the magnetosheath. The preditions of the study are in qualitative agreement with observations of whistler mode waves near the earth's bow shock.

Tokar, R. L.↗

Analysis of three-component aeromagnetic data

The numerical method of obtaining three field components from total field measurements, using double Fourier series expansion, is presented. The expressions for moments of the anomalous field components over a finite area are given. The magnitude and direction of the magnetization vector indicate that the vertical component of the magnetic field calculated from total field observations is more accurate at higher geomagnetic latitudes than at lower latitudes. The opposite is true for the horizontal components. The error in determining the magnetization vector directions are significantly large over most of the range of variation of declination and inclination of the vector, demonstrating the practical limitations of computing field components from total field data even under the best of conditions.

Bhattacharyya, B. K.↗

A comparison of two multi-variable integrator windup protection schemes

Two methods are examined for limit and integrator wind-up protection for multi-input, multi-output linear controllers subject to actuator constraints. The methods begin with an existing linear controller that satisfies the specifications for the nominal, small perturbation, linear model of the plant. The controllers are formulated to include an additional contribution to the state derivative calculations. The first method to be examined is the multi-variable version of the single-input, single-output, high gain, Conventional Anti-Windup (CAW) scheme. Except for the actuator limits, the CAW scheme is linear. The second scheme to be examined, denoted the Modified Anti-Windup (MAW) scheme, uses a scalar to modify the magnitude of the controller output vector while maintaining the vector direction. The calculation of the scalar modifier is a nonlinear function of the controller outputs and the actuator limits. In both cases the constrained actuator is tracked. These two integrator windup protection methods are demonstrated on a turbofan engine control system with five measurements, four control variables, and four actuators. The closed-loop responses of the two schemes are compared and contrasted during limit operation. The issue of maintaining the direction of the controller output vector using the Modified Anti-Windup scheme is discussed and the advantages and disadvantages of both of the IWP methods are presented.

Mattern, Duane↗

First geodetic measurement of convergence across the Java Trench

Convergence across the Java Trench has been estimated for the first time, from annual Global Positioning System (GPS) measurements commencing in 1989. The directions of motion of Christmas and Cocos Island are within 1 deg of that predicted by the No-Net Rotation (NNR) NUVEL-1 plate motion model for the Australian plate although their rates are 25% and 37% less than predcited, respectively. The motion of West Java differs significantly from the NNR NUVEL-1 prediction for the Eurasian plate with a 1 deg difference in direction and a 40% increase in rate. We infer that either West Java moves with a distinct Southeast Asian plate or this region experiences plate margin deformation. The convergence of Christmas Island with respect to West Java is 67 +/- mm/yr in a direction N11 deg E +/- 4 deg which is orthogonal to the trench. The magnitude of convergence agrees well with rescaled NUVEL-1 relative plate model which predicts a value of 71 mm/yr between Australia and Eurasia. The direction of motion matches the direction inferred from earthquake slip vectors at the trench but may be more northerly than the N20 deg E +/- 3 deg predicted by NUVEL-1. On June 2, 1994, almost a year after the last GPS survey, an M(sub W) = 7.5 earthquake with slip vector direction N5 deg occurred south of central Java.

Tregoning, P.↗

The Effect of Companions on the SIM Reference Frame

The Space Interferometry Mission (SIM) is a 10-m Michelson space-based optical interferometer designed for precision astrometry (4 microarcseconds, 3 microarcseconds/year) with better accuracy than before over a narrow field of view. One of the primary objectives of the SIM instrument is to determine accurately the directions to a grid of stars, together with their proper motions and parallax, improving a priori knowledge by nearly three orders of magnitude over Hipparcos and one order of magnitude over FAME's planned accuracy (Johnston, 2000). The instrument does not measure directly the angular separation between stars, but rather it measures the projection of each star's direction vector onto the interferometer baseline vector by measuring the pathlength delay of starlight as it passes through the two arms of the interferometer. The accuracy and stability of SIM's celestial reference frame is subject to degradation over the 5-year mission from the reflex motion induced by massive companions of the objects used to construct the celestial reference frame. The authors present the results of simulations that show the sensitivity of reference frame accuracy to companions as a function of mass and period. They assume that pre-launch ground surveys will eliminate all objects with RMS radial velocity greater than 10 m/s. They further assume that the standard astrometric parameters of position, parallax, and proper motion plus acceleration terms in right ascension and declination will be allowed to absorb reflex motion.

ASTROMETRY↗

The use of direct solvers for compressible flow computations

The feasibility of using a vectorized direct solver for the compressible Euler and Navier-Stokes equations is examined for both single-grid and multi-grid strategies. A procedure is developed for comparing the computational effort required for the direct method with that of the vertical line Gauss-Seidel iteration scheme in order to provide a criteria for choosing between the two techniques. The direct method is shown to have a relatively wide range of application on a vector processor with large memory. Results for both inviscid and viscous test problems in two-dimensions are examined.

Riggins, David↗

A vectorization of the Jameson-Caughey NYU transonic swept-wing computer program FLO-22-V1 for the STAR-100 computer

The computer program FLO-22 for analyzing inviscid transonic flow past 3-D swept-wing configurations was modified to use vector operations and run on the STAR-100 computer. The vectorized version described herein was called FLO-22-V1. Vector operations were incorporated into Successive Line Over-Relaxation in the transformed horizontal direction. Vector relational operations and control vectors were used to implement upwind differencing at supersonic points. A high speed of computation and extended grid domain were characteristics of FLO-22-V1. The new program was not the optimal vectorization of Successive Line Over-Relaxation applied to transonic flow; however, it proved that vector operations can readily be implemented to increase the computation rate of the algorithm.

Smith, R. E.↗

Experiments on near-wall structure of three-dimensional boundary layers

Investigations of three-dimensional turbulent boundary layers have shown basic differences between two- and three-dimensional flows. These differences can significantly impact the modeling of three-dimensional flows since many flow models are based on results from two-dimensional boundary layers. In many cases the shear stress vector direction has been shown to lag relative to the direction of the mean velocity gradient as the cross flow grows downstream. Coincidence of these vectors is necessary for a scalar eddy viscosity modeling assumption. A second effect is a reduction in magnitude of the shear stress and/or the shear stress to turbulence energy ratio, a(sub 1). This reduction has been observed in several experiments. Recent numerical simulations also indicate wall-layer structural differences between two- and three-dimensional boundary layers. The differences in structure between two- and three-dimensional boundary layers was also addressed in the experiment of Littell & Eaton. The experiment used two-point correlations to investigate the vortical structures in a three dimensional boundary layer on a spinning disk. It was found that each sign of longitudinal vortex is equally likely to exist, but one sign of vorticity is associated with a structure which is better at producing ejections. The goal of the current investigation is to study the structure of the inner layers. Among other questions, the differences between the effects deduced from the three-dimensional flow simulations and the effects seen in experiments can be examined. The research concentrates on the structure of the wall-layer through flow visualization and direct turbulence measurements down to y(+) = 5.

Flack, Karen A.↗

Parafrase restructuring of FORTRAN code for parallel processing

Parafrase transforms a FORTRAN code, subroutine by subroutine, into a parallel code for a vector and/or shared-memory multiprocessor system. Parafrase is not a compiler; it transforms a code and provides information for a vector or concurrent process. Parafrase uses a data dependency to reveal parallelism among instructions. The data dependency test distinguishes between recurrences and statements that can be directly vectorized or parallelized. A number of transformations are required to build a data dependency graph.

Wadhwa, Atul↗

An efficient approach for computing the geometrical optics field reflected from a numerically specified surface

An approach for computing the geometrical optic fields reflected from a numerically specified surface is presented. The approach includes the step of deriving a specular point and begins with computing the reflected rays off the surface at the points where their coordinates, as well as the partial derivatives (or equivalently, the direction of the normal), are numerically specified. Then, a cluster of three adjacent rays are chosen to define a 'mean ray' and the divergence factor associated with this mean ray. Finally, the ampilitude, phase, and vector direction of the reflected field at a given observation point are derived by associating this point with the nearest mean ray and determining its position relative to such a ray.

Mittra, R.↗

The auxiliary problem for feasible directions

The present work compares the computational efficiency of two means of determining a direction vector S at constraint boundaries and demonstrates the advantages of using the alternative that formulates the S-determination subproblem with a single quadratic constraint, which in turn can be cast as a special linear problem and solved directly.

Gwin, L. B.↗

Ponderomotive Force in the Presence of Electric Fields

This paper presents averaged equations of particle motion in an electromagnetic wave of arbitrary frequency with its wave vector directed along the ambient magnetic field. The particle is also subjected to an E cross B drift and a background electric field slowly changing in space and acting along the magnetic field line. The fields, wave amplitude, and the wave vector depend on the coordinate along the magnetic field line. The derivations of the ponderomotive forces are done by assuming that the drift velocity in the ambient magnetic field is comparable to the particle velocity. Such a scenario leads to new ponderomotive forces, dependent on the wave magnetic field intensity, and, as a result, to the additional energy exchange between the wave and the plasma particles. It is found that the parallel electric field can lead to the change of the particle-wave energy exchange rate comparable to that produced by the previously discussed ponderomotive forces.

Plasmas↗

Buoyancy Effects upon Vapor Flame and Explosion Processes

The objective of this microgravity project is to develop an experimental and theoretical analyses critical to the understanding of the coupling of buoyancy and turbulence generation and its effect on fuel-air mixing, flame intensity and flame propagation in jet diffusion flames. The experiment is designed to examine certain effects of buoyancy acting on a diffusion flame in which the flame is directed either upward or downward. This change from negative to positive g is observed to significantly alter the flame shape although all other operating conditions are the same for both configurations. However, to perform this experiment a significant coaxial secondary air flow is needed in order to prevent flow reversal when the flame is inverted. The theoretical analysis that has been developed handles the secondary air flow and the extreme change in gravity vector direction. Thus the data will provide a measure of credibility of the analysis which will then be used to assist in the design of the actual zero-g experiment.

Edelman, R. B.↗

The collective gyration of a heavy ion cloud in a magnetized plasma

In both the ionospheric barium injection experiments CRIT 1 and CRIT 2, a long duration oscillation was seen with a frequency close to the gyro frequency of barium and a time duration of about one second. A model for the phenomena which was proposed for the CRIT 1 experiment is compared to the results from CRIT 2 which made a much more complete set of measurements. The model follows the motion of a low Beta ion cloud through a larger ambient plasma. The internal field of the model is close to antiparallel to the injection direction v sub i but slightly tilted towards the self polarization direction E sub p = -V sub i by B. As the ions move across the magnetic field, the space charge is continuously neutralized by magnetic field aligned electron currents from the ambient ionosphere, drawn by the divergence in the perpendicular electric field. These currents give a perturbation of the magnetic field related to the electric field perturbation by Delta E/Delta B approximately equal to V sub A. The model predictions agree quite well with the observed vector directions, field strengths, and decay times of the electric and magnetic fields in CRIT 2. The possibility to extend the model to the active region, where the ions are produces in this type of self-ionizing injection experiments, is discussed.

Brenning, N.↗