SEARCH · Search NASA
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.
Nonlinear equations solved with linear approximations applicable during small increments of time. Report presents theoretical study of selected aspects of dynamics and control of robotic manipulator that includes multiple flexible links. Addresses problem of control primarily from perspective of inverse-dynamics problem.
Conventional ultrasonic C scan used to detect flaws in graphite/epoxy composite panels modified to enhance sensitivity to transverse cracks. Modification involves reorientation and repositioning of ultrasonic transmitter and receiver to take advantage of mode-conversion phenomena at interfaces between liquid couplant and panel to be probed. Angle of incidence chosen so only shear waves propagate obliquely through panel. Although mode-conversion equations solved to take full advantage of technique somewhat complicated, they are straightforward, and technique implemented easily in practice.
A parallel program to analyze transient finite element problems was written and implemented on a system of transputer processors. The program uses the explicit time integration algorithm which eliminates the need for equation solving, making it more suitable for parallel computations. An interprocessor communication scheme was developed for arbitrary two dimensional grid processor configurations. Several 3-D problems were analyzed on a system with a small number of processors.
Backscatter lidars are very useful tools in gaining information about the atmosphere. Inversion of backscatter signals can be used to retrieve the spatial distribution of clouds and aerosols, e.g., vertical extensions of clouds and aerosol layers, detection of 'subvisible' cirrus, etc., very accurately compared to passive radiometers. To retrieve extinction coefficient profiles, in general, the so-called lidar equation is solved. This equation describes a laser pulse which propagates through an atmosphere with transmission exp(- integral of sigma(z')dz'), which is scattered singly under 180 degrees towards a receiver, and which is attenuated on its way back again by exp(- integral of sigma(z')dz'). For inversion, further information on optical properties of the scatterers are required (extinction to backscatter ratio). However, it is evident that, at least for optically thick clouds and/or large receiver's field-of-views, this approach may be erroneous because of multiply scattered photons. Several studies were performed to simulate multiple scattering in case of collimated beams or pulses including Monte Carlo, transport, and diffusion theory. Unfortunately, most of the approaches require an enormous amount of computer time. Thus, it is meaningful to elaborate approximate methods, and one such method is discussed. Our method is based on Bissonnett (1988). We have improved the numerical simulation of limited field-of-view detectors, and compared the results to Monte Carlo calculations and an exact model of Altmann (1989).
We present a method for solving plasma magnetohydrodynamic (MHD) problems arising from the interaction of plasmas with magnetic fields in stellar atmospheres. Our approach, in contrast to previous methods, is not based on solving equations for the magnetic field and plasma velocity but rather studies the evolution of the electric current and density (and the related gas pressure). We have applied the method to several studies involving linearized departures from static, current-free equilibria. The applications show explicit solutions for cases found in astrophysics and to problems encountered with earlier studies where the gas pressure was neglected. The method is particularly well suited for studying situations which involve a transition between high and low plasma-beta regions. It shows precisely how electric currents, and magnetic free-energy, build up in the plasma as a result of the slow stressing of a potential magnetic field configuration. The method also demonstrates how transverse-current waves, a mix of Alfven and magneto-acoustic modes, propagate in a low-beta plasma for any density stratification and background field geometry.
Several forms of trapped liquid found within nakhlite meteorites have been examined, including interstitial melt and magmatic inclusions within the cores of large olivine grains. Differences in the mineralogy and texture between two types of trapped melt inclusions, and between these inclusions and the mesostasis, indicate that vitrophyric inclusions are most appropriate for estimating the composition of a nakhlite parental magma in equilibrium with early-forming olivine and augite. Parent liquids were calculated from the mineralogy of large inclusions in Nakhla and Governador Valadares, using a system of mass-balance equations solved by linear regression methods. The chosen parental liquids were cosaturated in olivine and augite and had Mg/Fe values consistent with measured augite/liquid Kds. These parental magma compositions are similar to other published compositions for Nakhla, Chassigny, and Shergotty parental melts, and may correspond to a significant magma type on Mars.
One-equation model of turbulence based on standard equations of k-epsilon model of turbulence, where k is turbulent energy and e is rate of dissipation of k. Derivation of one-equation model motivated partly by inaccuracies of flows computed by some Navier-Stokes-equations-solving algorithms incorporating algebraic models of turbulence. Satisfies need to avoid having to determine algebraic length scales.
Thermal Radiation Analyzer System (TRASYS) computer program is software program having generalized capability to solve equations of radiation-related aspects of thermal-analysis problems. Computes total thermal-radiation environment for spacecraft in orbit. Software calculates internode-radiation-interchange data as well as data on rates of incidence and absorption of heat originating from environmental radiant sources. Provides data of both types in format directly usable by such thermal-analyzer programs as SINDA '85/FLUINT (available from COSMIC, program number MSC-21528). CRAY version of TRASYS (P25) written in FORTRAN 77. Other versions available upon request.
Two related computer programs written for use in analyzing leakage through cracks. Leakage flow laminar or turbulent. One program used to determine dimensions of crack under given flow conditions and given measured rate of leakage. Other used to determine rate of leakage of gas through crack of given dimensions under given flow conditions. Programs, written in BASIC language, accelerate and facilitate iterative calculations and parametric analyses. Solve equations of Fanno flow. Enables rapid solution of leakage problem.
Submicron size monodisperse particles are of interest in many industrial and scientific applications. These include the manufacture of ceramic parts using fine ceramic particles, the production of thin films by deposition of ionized clusters, monodisperse seed particles for laser anemometry, and the study of size dependence of cluster chemical and physical properties. An inexpensive and relatively easy way to generate such particles is by utilizing the phenomenon of spontaneous condensation. The phenomenon occurs when the vapor or a mixture of a vapor and a noncondensing gas is expanded at a high expansion rate. The saturation line is crossed with the supercooled vapor behaving like a gas, until all of a sudden at the so called Wilson point, condensation occurs, resulting in a large number of relatively monodisperse droplets. The droplet size is a function of the expansion rate, inlet conditions, mass fraction of vapor, gas properties, etc. Spontaneous condensation of steam and water vapor and air mixture in a one dimensional nozzle was modeled and the resulting equations solved numerically. The droplet size distribution at the exit of various one dimensional nozzles and the flow characteristics such as pressure ratio, mean droplet radius, vapor and droplet temperatures, nucleation flux, supercooling, wetness, etc., along the axial distance were obtained. The numerical results compared very well with the available experimental data. The effect of inlet conditions, nozzle expansion rates, and vapor mass fractions on droplet mean radius, droplet size distribution, and pressure ratio were examined.
A procedure is presented for performing three-dimensional stress analysis of spiral bevel gears in mesh using the finite element method. The procedure involves generating a finite element model by solving equations that identify tooth surface coordinates. Coordinate transformations are used to orientate the gear and pinion for gear meshing. Contact boundary conditions are simulated with gap elements. A solution technique for correct orientation of the gap elements is given. Example models and results are presented.
The characteristics of nonlinear two-dimensional horizontally periodic double-diffusive fingering convection in a saturated porous medium is investigated, using the Darcy equation including Brinkman and Forchheimer terms for the momentum equation. To solve the equations and the corresponding initial and boundary conditions, a Galerkin method is applied in the horizontal direction, and a hybrid finite difference method is used in the vertical direction. The developed computer code was used to compute the thermal convection case, and the results were found to be in good agreement with existing results.
A structural model is developed for the linear elastic response of structures consisting of multiple layers of varying thickness such as laminated composites containing internal ply drop-offs. The assumption of generalized plane deformation is used to reduce the solution domain to two dimensions while still allowing some out-of-plane deformation. The Hellinger-Reissner variational principle is applied to a layerwise assumed stress distribution with the resulting governing equations solved using finite differences.
A procedure is presented for performing three-dimensional stress analysis of spiral bevel gears in mesh using the finite element method. The procedure involves generating a finite element model by solving equations that identify tooth surface coordinates. Coordinate transformations are used to orientate the gear and pinion for gear meshing. Contact boundary conditions are simulated with gap elements. A solution technique for correct orientation of the gap elements is given. Example models and results are presented.
TDIGG is fast and versatile computer program for generating two-dimensional computational grids for use in programs solving equations of flow by finite-difference methods. Both algebraic and elliptic grid-generation systems included. Enables user to view results of each iteration. Written in FORTRAN 77.
A Workstation has been designed and constructed for rapidly simulating motions of rigid and elastic multibody systems. We examine the Workstation from the point of view of analysts who use the machine in an industrial setting. Two aspects of the device distinguish it from other simulation programs. First, one uses a series of windows and menus on a computer terminal, together with a keyboard and mouse, to provide a mathematical and geometrical description of the system under consideration. The second hallmark is a facility for animating simulation results. An assessment of the amount of effort required to numerically describe a system to the Workstation is made by comparing the process to that used with other multibody software. The apparatus for displaying results as a motion picture is critiqued as well. In an effort to establish confidence in the algorithms that derive, encode, and solve equations of motion, simulation results from the Workstation are compared to answers obtained with other multibody programs. Our study includes measurements of computational speed.
Three dimensional stress analysis of spiral bevel gears in mesh using the finite element method is presented. A finite element model is generated by solving equations that identify tooth surface coordinates. Contact is simulated by the automatic generation of nonpenetration constraints. This method is compared to a finite element contact analysis conducted with gap elements.
A new and very efficient numerical method for solving equations of the Helmholtz type is specialized for problems having axisymmetric geometry. It is then demonstrated by application to the classical problem of acoustic radiation from a vibrating piston set in a stationary infinite plane. The method utilizes 'Green's Function Discretization', to obtain an accurate resolution of the waves using only 2-3 points per wave. Locally valid free space Green's functions, used in the discretization step, are obtained by quadrature. Results are computed for a range of grid spacing/piston radius ratios at a frequency parameter, omega R/c(sub 0), of 2 pi. In this case, the minimum required grid resolution appears to be fixed by the need to resolve a step boundary condition at the piston edge rather than by the length scale imposed by the wave length of the acoustic radiation. It is also demonstrated that a local near-field radiation boundary procedure allows the domain to be truncated very near the radiating source with little effect on the solution.