SEARCH · Engineering Papers
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.
Appendix 4 of the Study of One- and Two-Dimensional Filtering and Deconvolution Algorithms for a Streaming Array Computer discusses coordinate axes, location of origin, and redundancy for the one- and two-dimensional Fourier transform for complex and real data.
A geometrically-based derivation of Kane's dynamical equations is presented. Equations for both holonomic and nonholonomic systems are derived by considering the system's motion in a hypersurface determined from the equations relating Cartesian to generalized coordinates. Using vector space methods, the equations of motion are projected onto the tangent plane to the hypersurface. In the course of this construction, a requirement on the transformation between generalized and Cartesian coordinates is revealed that is often overlooked. This restriction is then shown to be a necessary and sufficient condition for the invertibility of the coefficient matrix of the generalized accelerations appearing in Kane's dynamical equations. Although less succinct than the traditional approach, the present derivation offers some insight into the physics embodied in Kane's equations and appears as a natural generalization of methods used in elementary problems.
The main result of this paper is the derivation of a new expression for the tracer subgrid term in level coordinates S(l) to be employed in O-GCM. The novel feature is the proper account of the random nature of the density field which strongly affects the transformation from isopycnal to level coordinates of the variables of interest, velocity and tracer fields, their correlation functions and ultimately the subgrid terms. In deriving our result we made use of measured properties of vertical ocean turbulence. The major new results are: 1) the new subgrid expression is different from that of the heuristic GM model, 2) u++(tracer)=1/2u+(thickness), where u++ and u+ are the tracer and thickness bolus velocities. In previous models, u++ = u+, 2) the subgrid for a tracer tau is not the same as that for the density rho even when one accounts for the obvious absence of a diffusion term in the latter. The difference stems from a new treatment of the stochastic nature of the density, 3) the mesoscale diffusivity enters both locally and non-locally, as the integral over all z's from the bottom of the ocean to the level z.
A video-rate coordinate remapper includes a memory for storing a plurality of transformations on look-up tables for remapping input images from one coordinate system to another. Such transformations are operator selectable. The remapper includes a collective processor by which certain input pixels of an input image are transformed to a portion of the output image in a many-to-one relationship. The remapper includes an interpolative processor by which the remaining input pixels of the input image are transformed to another portion of the output image in a one-to-many relationship. The invention includes certain specific transforms for creating output images useful for certain defects of visually impaired people. The invention also includes means for shifting input pixels and means for scrolling the output matrix.
This paper explores a VLSI architecture for geometrical mapping address computation. The geometric transformation is discussed in the context of plane projective geometry, which invokes a set of basic transformations to be implemented for the general image processing. The homogeneous and 2-dimensional cartesian coordinates are employed to represent the transformations, each of which is implemented via an augmented CORDIC as a processing element. A specific scheme for a processor, which utilizes full-pipelining at the macro-level and parallel constant-factor-redundant arithmetic and full-pipelining at the micro-level, is assessed to produce a single VLSI chip for HDTV applications using state-of-art MOS technology.
Using 21 days of GPS data from 21 globally distributed receivers operating during early 1991, a 7-parameter transformation between a GPS free-network solution and coordinates of 12 stations listed in the International Terrestrial Reference Frame (ITRF) is solved. Standard errors of GPS coordinates are derived by applying an orthogonal projection operator to the free-network covariance. The weighted rms difference between 33 transformed GPS and ITRF coordinates is 12 mm in the Northern Hemisphere. Best results are obtained by mapping ITRF coordinates to the epoch of this experiment assuming no vertical site motions. Fixing selected sites in the GPS solution to ITRF '90 does not improve the agreement. It is concluded that the use of fiducial constraints is unnecessary for global networks.
Mapping lower stratospheric ozone into potential vorticity (PV)- potential temperature (Theta) coordinates is a common technique employed to analyze sparse data sets. Ozone transformed into a flow-following dynamical coordinate system is insensitive to meteorological variations. Therefore data from a wide range of times/locations can be compared, so long as the measurements were made in the same airmass (as defined by PV). Moreover, once a relationship between ozone and PV/Theta is established, a full 3D ozone field can be estimated from this relationship and the 3D analyzed PV field. However, ozone data mapped in this fashion can be hampered by noisy PV fields, or "mis-matches" in the resolution and/or exact location of the ozone and PV measurements. In this study, we investigate the PV-ozone relationship using output from a recent 50-year run of the Goddard 3D chemical transport model (CTM). Model constituents are transported using off-line dynamics from the finite volume general circulation model (FVGCM). By using the internally consistent model PV and ozone fields, we minimize noise due to mis-matching and resolution issues. We calculate correlations between model ozone and PV throughout the stratosphere, and test the sensitivity of the technique to initial data resolution. To do this we degrade the model data to that of various satellite instruments, then compare the mapped fields derived from the sub-sampled data to the full resolution model data. With these studies we can determine appropriate limits for the PV-theta mapping technique in latitude, altitude, and as a function of original data resolution.
NASA's Human Space Flight program depends heavily on Extra-Vehicular Activities (EVA's) performed by human astronauts. EVA is a high risk environment that requires extensive training and ground support. In collaboration with the Defense Advanced Research Projects Agency (DARPA), NASA is conducting a ground development project to produce a robotic astronaut's assistant, called Robonaut, that could help reduce human EVA time and workload. The project described in this paper designed and implemented a hand-eye calibration scheme for Robonaut, Unit A. The intent of this calibration scheme is to improve hand-eye coordination of the robot. The basic approach is to use kinematic and stereo vision measurements, namely the joint angles self-reported by the right arm and 3-D positions of a calibration fixture as measured by vision, to estimate the transformation from Robonaut's base coordinate system to its hand coordinate system and to its vision coordinate system. Two methods of gathering data sets have been developed, along with software to support each. In the first, the system observes the robotic arm and neck angles as the robot is operated under external control, and measures the 3-D position of a calibration fixture using Robonaut's stereo cameras, and logs these data. In the second, the system drives the arm and neck through a set of pre-recorded configurations, and data are again logged. Two variants of the calibration scheme have been developed. The full calibration scheme is a batch procedure that estimates all relevant kinematic parameters of the arm and neck of the robot The daily calibration scheme estimates only joint offsets for each rotational joint on the arm and neck, which are assumed to change from day to day. The schemes have been designed to be automatic and easy to use so that the robot can be fully recalibrated when needed such as after repair, upgrade, etc, and can be partially recalibrated after each power cycle. The scheme has been implemented on Robonaut Unit A and has been shown to reduce mismatch between kinematically derived positions and visually derived positions from a mean of 13.75cm using the previous calibration to means of 1.85cm using a full calibration and 2.02cm using a suboptimal but faster daily calibration. This improved calibration has already enabled the robot to more accurately reach for and grasp objects that it sees within its workspace. The system has been used to support an autonomous wrench-grasping experiment and significantly improved the workspace positioning of the hand based on visually derived wrench position. estimates.
A time domain subspace iteration technique is presented to compute a set of normal modes from the measured complex modes. By using the proposed method, a large number of physical coordinates are reduced to a smaller number of model or principal coordinates. Subspace free decay time responses are computed using properly scaled complex modal vectors. Companion matrix for the general case of nonproportional damping is then derived in the selected vector subspace. Subspace normal modes are obtained through eigenvalue solution of the (M sub N) sup -1 (K sub N) matrix and transformed back to the physical coordinates to get a set of normal modes. A numerical example is presented to demonstrate the outlined theory.
It was discovered some years ago by Schiff that two equations for fields in vacuum do not carry over without change from an inertial frame to a frame with rotating axes of space coordinates, even for a region with all velocities much lower than the speed of light. However, the belief that all four of the field equations are invariant under such conditions is still prevalent and causes misconceptions in physical applications, including astrophysical and geophysical ones. The purpose of the present paper is therefore to call attention to Schiff's discovery, discussing its basis and its extension to fields in material media, and to interpret the additional terms that must be added to the equations in order to obtain valid transformations to rotating axes of coordinates.
One of the difficulties in determining constituent trends from the ER-2 flight data is the large amount of day to day variability generated by the motion of the polar vortex. To reduce this variability, the observations have been transformed into the conservative (Lagrangian) reference frames consisting of the coordinate pairs, potential temperature (PT) and potential vorticity (PV), or PT and N2O. The requirement of only two independent coordinates rests on the assumption that constituent distributions and their chemical processes are nearly zonal in that coordinate system. Flight data is used everywhere for these transformation except for potential vorticity. Potential vorticity is determined from level flight segments, and NMC PV values during flight dives and takeoffs are combined with flight data in a smooth fashion.
The thinning of the near-Earth current sheet during the growth phase of the Coordinated Data Analysis Workshop (CDAW) 6 magnetospheric substorm is studied. The expansion onset of the substorm occurred at 1054 UT, March 22, 1979. During the growth phase, two spacecraft, International Sun Earth Explorer (ISEE) 1 and ISEE 2, were within the current sheet approximately 13 R(sub E) from the Earth and obtained simultaneous high-resolution magnetic data at two points in the current sheet. Plasma data were also provided by the ISEE spacecraft and solar wind data by IMP 8. To facilitate the analysis, the GSM magnetic field data are transformed to a 'neutral sheet coordinate system' in which the new x axis is parallel to the average magnetic field above and below the neutral sheet and the new y axis lies in the GSM equatorial plane. A model based on the assumption that the current sheet is a time-invariant structure fails to predict neutral sheet crossing times. Consequently, the Harris sheet model, which allows one to remove the restriction of time invariancy, is used instead. It is found that during the growth phase, a model parameter corresponding to the thickness of the current sheet decreased exponentially from about 5 R(sub E) to 1 R(sub E) with a time constant of about 14 min. In addition, the ISEE 1 and ISEE 2 neutral sheet crossings after expansion onset indicate that the neutral sheet was moving upward at 7 km/s relative to the spacecraft. Since both crossings occurred in approximately 80 s, the current sheet thickness is estimated to be about 500 km. These results demonstrate that the near-Earth current sheet undergoes dramatic thinning during the substorm growth phase and expansion onset.
The combination of an orthogonal, curvilinear coordinate generation procedure with a stable forward marching viscous flow solution technique is presently employed in the solution of flow fields for arbitrary, axisymmetric ducts. Coordinate generation is accomplished by means of both potential lines and plane potential flow streamlines. Since the coordinate streamlines approximate actual ones, the equations of motion for viscous compressible flow can be parabolized in order to solve for both the boundary layer and the core flow in a single streamwise pass. The method's versatility is demonstrated by two examples of viscous compressible swirling flow through complex radial gas turbine passages.
Component mode synthesis (CMS) is a method of dynamic analysis, for structures having a large number of degrees of freedom (DOF). These structures often required lengthy computer CPU time and large computer memory resources, if solved directly by the finite-element method (FEM). In CMS, the structure is divided into independent components in which the DOF are defined by a set of generalized coordinates defined by displacement shapes. The number of the generalized coordinates are much less than the original number of physical DOF, in the component. The displacement shapes are used to transform the component property matrices and any applied external loads, to a reduced system of coordinates. Reduced system property matrices are assembled, and any type of dynamic analysis is carried out in the reduced coordinate system. Any obtained results are back transformed to the original component coordinate systems. In all conventional methods of CMS, the mode shapes used for components are dynamic mode shapes, supplemented by static deflected shapes. Historically, all the dynamic mode shapes used in conventional CMS are the natural modes (eigenvectors) of components. A new method of CMS, namely the boundary flexibility vector method of CMS, is presented. The method provides for the incorporation of a set of static Ritz vectors, referred to as boundary flexibility vectors, as a replacement and/or supplement to conventional eigenvectors, as displacement shapes for components. The generation of these vectors does not require the solution of a costly eigenvalue problem, as in the case of natural modes in conventional CMS, and hence a substantial saving in CPU time can be achieved. The boundary flexibility vectors are generated from flexibility (or stiffness) properties of components. The formulation presented is for both free and fixed-interface components, and for both the free and forced vibration problems. Free and forced vibration numerical examples are presented to verify the accuracy of the method and the saving in CPU time. Compared to conventional methods of CMS, the results indicate that by using the new method, more accurate results can be obtained with a substantial saving in CPU time.
A technique for measuring wing twist currently in use at the National Transonic Facility is described. The technique is based upon a single camera photogrammetric determination of two dimensional coordinates with a fixed (and known) third dimensional coordinate. The wing twist is found from a conformal transformation between wind-on and wind-off 2-D coordinates in the plane of rotation. The advantages and limitations of the technique as well as the rationale for selection of this particular technique are discussed. Examples are presented to illustrate run-to-run and test-to-test repeatability of the technique in air mode. Examples of wing twist in cryogenic nitrogen mode are also presented.
The efficient low dissipative high order schemes proposed by Yee et al. is formulated for 3-D curvilinear moving grids. These schemes consists of a high order base schemes combined with nonlinear characteristic filters. The amount of numerical dissipation is minimized by applying the schemes to the entropy splitting form of the inviscid flux derivatives. The analysis is given for a thermally perfect gas. The main difficulty in the extension of higher order schemes that were formulated in Cartesian coordinates to curvilinear moving grids is the higher order transformed metric evaluations. The higher order numerical evaluation of the transformed metric terms to insure freestream preservation is done in a coordinate invariant manner. The formulation is an improvement over existing formulation of high order scheme in curvilinear moving grids.
Gravitational potential of oblate planet, expressed in terms of oblate spheroidal coordinates, generalized by means of transformation of associated Cartesian system
Acoustic propagation in two-dimensional corrugated ducts is analyzed. A transform is made to a coordinate system that matches the duct walls. The resulting wave equation separates in space and time. The spatial differential equation is then solved by expanding in the eigenfunctions of a flat-walled duct. The prescription for this expansion is given by the appropriate application of time-independent perturbation theory. Calculations of first-order corrections for a periodic duct of both finite and infinite length are presented.