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.
A grid generation program for use with the CAP-TSD transonic small disturbance code is described. The program runs interactively in FORTRAN on the Sun Workstation. A fifth-degree polynomial is used to map the grid index onto the computational coordinate. The grid is plotted to aid in the assessment of its quality and may be saved on file in NAMELIST format.
This document will describe the technology developed at Kennedy Space Center (KSC) to demonstrate automatic tracking, docking, and mating of umbilical systems. Specifically the use of a real time six degree of freedom (6DOF) target tracking vision system, (developed by Adaptive Automation, Inc. under contract to KSC), is discussed in detail. The paper describes the use of mechanical compliance in the docking, mating, and tracking-after-mating operations. The vision system computes six coordinates that define the position and orientation of a three dimensional target using data from a single CCD camera. The camera is mounted on a 6DOF robot arm. After target coordinates are computed, they are transmitted to a supervisory computer which controls the robot motion in real time. Details of the image processing algorithms, image processing hardware, and target configuration used in the vision system are discussed in the paper. The motion of the space vehicle relative to the service structure after mating led us to the development of a compliant system that allows enough displacement of the target relative to the camera so that tracking after mating may continue. This reduces internal stresses between flight and ground hardware.
Under a NASA MSFC SBIR Phase 2 effort an analysis has been developed for liquid flows with vapor bubbles such as those in liquid rocket engine components. The analysis is based on a combined Eulerian-Lagrangian technique, in which Eulerian conservation equations are solved for the liquid phase, while Lagrangian equations of motion are integrated in computational coordinates for the vapor phase. The novel aspect of the Lagrangian analysis developed under this effort is that it combines features of the so-called particle distribution approach with those of the so-called particle trajectory approach and can, in fact, be considered as a generalization of both of those traditional methods. The result of this generalization is a reduction in CPU time and memory requirements. Particle time step (stability) limitations have been eliminated by semi-implicit integration of the particle equations of motion (and, for certain applications, the particle temperature equation), although practical limitations remain in effect for reasons of accuracy. The analysis has been applied to the simulation of cavitating flow through a single-bladed section of a labyrinth seal. Models for the simulation of bubble formation and growth have been included, as well as models for bubble drag and heat transfer. The results indicate that bubble formation is more or less 'explosive'. for a given flow field, the number density of bubble nucleation sites is very sensitive to the vapor properties and the surface tension. The bubble motion, on the other hand, is much less sensitive to the properties, but is affected strongly by the local pressure gradients in the flow field. In situations where either the material properties or the flow field are not known with sufficient accuracy, parametric studies can be carried out rapidly to assess the effect of the important variables. Future work will include application of the analysis to cavitation in inducer flow fields.
Curvilinear coordinate systems were used extensively to solve partial differential equations on arbitrary regions. An analysis of truncation error in the computation of derivatives revealed why numerical results may be erroneous. A more accurate method of computing derivatives is presented.
Cartesian coordinates and partial derivatives of two-body problem as double-precision Fortran 4 program for IBM 7094
Explore the source record for details and available documents.
In connection with the employment of the sigma coordinates introduced by Phillips (1957), problems can arise regarding an accurate finite-difference computation of the pressure gradient force. Over steeply sloped terrain, the calculation of the sigma-coordinate pressure gradient force involves computing the difference between two large terms of opposite sign which results in large truncation error. To reduce the truncation error, several finite-difference methods have been designed and implemented. The present investigation has the objective to provide another method of computing the sigma-coordinate pressure gradient force. Phillips' method is applied for the elimination of a hydrostatic component to a flux formulation. The new technique is compared with four other methods for computing the pressure gradient force. The work is motivated by the desire to use an isentropic and sigma-coordinate hybrid model for experiments designed to study flow near mountainous terrain.
Application of digital computers for deriving motion equations of point mass in curvilinear coordinate systems by use of tensor calculus
The method of Smith and Wiegel was used to generate meshes for mixer lobes and subsonic inlets that are compatible with flow analysis codes requiring a boundary fitted coordinate system. Successful application of this mesh generator required development of procedures to distribute the mesh points along the boundaries, to regulate the dependence of the connecting function to the local boundary slope, to concentrate the mesh into regions of special interest, and to modify the mesh grid so that it possessed a smooth progression of cell metrics and cell volumes in all directions. The method of Smith and Wiegel when coupled with the extensions mentioned above has proven to be easy to use and control for the inlet and mixer lobe geometries investigated.
Preliminary Mariner 9 map products were based on computed aerocentric coordinates at the point where the optic path intercepts a spherical reference figure for Mars. The intercept point on a reference ellipsoid would in general have slightly different values of radius, latitude, and longitude. For near-vertical viewing used for mapping, these differences are on the order of 0.1 in latitude and longitude, which is negligible compared to the uncertainties of the initial uncontrolled mosaics and maps.
A numerical procedure has been developed for the computation of supersonic flows over complex conical geometries. The full potential equation is solved using a finite-volume method with a non-body-fitted rectangular grid. The only mapping done is the transformation of the spherical cross-flow plane to a flat surface using a stereographic projection. A new procedure for very thin fins is described which does not require the resolution of the fin thickness. Applications for simple cones, conical wing-bodies, wave riders and finned geometries compare favorably with existing solutions with body-fitted grids and available experimental data.
An analytical transformation method which can map arbitrary physical coordinate grid distribution into desired computational coordinate with uniform grid distribution is derived. The transformation function and its higher derivatives are differentiable. Salient features include; (1) precise control of grid sizes; (2) more than one location of clustered grids; (3) exact positioning of particular computational nodes in the physical plane; and (4) ensuring several patches of uniformly spaced grids in the physical plane for the higher accuracies (such as at the boundaries), etc., while keeping the variation of grid spacing continuous to avoid numerical instability.
The objective of this study is two-fold: to investigate the absorption of hard x rays along a tangential path through the terrestrial atmosphere, and to locate, other than the use of imaging devices, the flare source of x ray emission in the sun. Observational data were taken with the Hard X Ray Burst Spectrometer (HXRBS) on board the Solar Maximum Mission satellite (SMM). The occultation event chosen for investigation is the flare on 14 October 1981. Night fell on SMM during the decaying phase of the flare. In about fifteen seconds, the x ray flux decreased to zero. The model for computing synthetic light curve, which is directly comparable to the observed curve, consists of two computing programs. One is a modified version of an existing code, which compute the coordinates of the tangent points of the SMM-sun lines. The other program computed the atmospheric densities along the angles (N to S and E to W) measured from the center of the sun. The various model also depends on the destination of the terrestrial atmosphere at geocentric positions of SMM and the sun, the mean atmospheric densities at different heights and latitudes, and the x ray absorption cross of nitrogen, oxygen, and argon. The optical length, at a given time and for a x ray energy, is computed with the two angular parameters. The fit of a particular synthetic light, curve, i.e., the exponentials of the optical lengths, to the observed eclipse curve, yields the desired angular parameters.
Coordinates of points on complicated three-dimensional surface of turbo-machine blade computed by method partly simplifying three-dimensional problem to two-dimensional problem. Method implemented in computer program enabling rapid hydrodynamical evaluation and optimization of blade surface prior to definition of final design of, and analysis of stresses in blade. By incorporating program into procedure for designing turbomachine, one reduces time spent on design iterations.
A computer software system has been developed to automatically generate two-dimensional coordinates from algebraic transformations. For topologically complex regions, a smooth assembly of the transformations can be used to automatically produce a composite mesh where a general gridded format is retained. The algebraic mesh generation system consists of a collection of operator subroutines which are applied to an established data structure and which automatically perform the necessary parts of mesh construction from a sequence of multisurface transformations. The system operators are discussed, taking into account the data base, the order of application, direct surface generators, geometric surface operators, surface generators from existing surfaces, transverse operators, mesh operators, assembly operators, and data visualization operators. Attention is given to applications related to airfoils.
The Framework for Interdisciplinary Design Optimization (FIDO) is a general programming environment for automating the distribution of complex computing tasks over a networked system of heterogeneous computers. For example, instead of manually passing a complex design problem between its diverse specialty disciplines, the FIDO system provides for automatic interactions between the discipline tasks and facilitates their communications. The FIDO system networks all the computers involved into a distributed heterogeneous computing system, so they have access to centralized data and can work on their parts of the total computation simultaneously in parallel whenever possible. Thus, each computational task can be done by the most appropriate computer. Results can be viewed as they are produced and variables changed manually for steering the process. The software is modular in order to ease migration to new problems: different codes can be substituted for each of the current code modules with little or no effect on the others. The potential for commercial use of FIDO rests in the capability it provides for automatically coordinating diverse computations on a networked system of workstations and computers. For example, FIDO could provide the coordination required for the design of vehicles or electronics or for modeling complex systems.
A system of computer programs was developed to model general three dimensional surfaces. Surfaces are modeled as sets of parametric bicubic patches. There are also capabilities to transform coordinates, to compute mesh/surface intersection normals, and to format input data for a transonic potential flow analysis. A graphical display of surface models and intersection normals is available. There are additional capabilities to regulate point spacing on input curves and to compute surface/surface intersection curves. Input and output data formats are described; detailed suggestions are given for user input. Instructions for execution are given, and examples are shown.
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.