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.
Programs supporting the calculation of topocentric spherical coordinates of satellites from satellite images on negatives are discussed. Particular attention is given to a program for automatic identification of negatives obtained by the AFU-75 camera, and to the minimization of the amount of input information.
Numerical evaluations of heat transfer in a fissioning uranium plasma core reactor cavity, operating with seeded hydrogen propellant, was performed. A two-dimensional analysis is based on an assumed flow pattern and cavity wall heat exchange rate. Various iterative schemes were required by the nature of the radiative field and by the solid seed vaporization. Approximate formulations of the radiative heat flux are generally used, due to the complexity of the solution of a rigorously formulated problem. The present work analyzes the sensitivity of the results with respect to approximations of the radiative field, geometry, seed vaporization coefficients and flow pattern. The results present temperature, heat flux, density and optical depth distributions in the reactor cavity, acceptable simplifying assumptions, and iterative schemes. The present calculations, performed in cartesian and spherical coordinates, are applicable to any most general heat transfer problem.
The objective of the paper is to examine a class of non-force-free fields analytically. Specifically, magnetic fields in magnetohydrostatic equilibrium with a plasma in a gravitational field are treated in an approximation of two independent variables but three vector components. Spherical coordinates are emphasized, although the formal results for cylindrical and Cartesian coordinates are presented in the appendices. Formal solutions for the magnetic-field components are obtained in terms of the plasma variables, and field line equations are derived. A final equation governing the plasma variables is then obtained. Procedures are developed for analyzing this equation and obtaining sets of self-consistent particular solutions to the governing equations; a number of such sets of solutions are presented. As an example, one solution set is examined, illustrating the application of the results to the analysis of solar observational data.
Modifications to explicit finite difference schemes for solving the shallow water equations for meteorological applications by increasing the time step for the fast gravity waves are analyzed. Terms associated with the gravity waves in the shallow water equations are treated on a coarser grid than those associated with the slow Rossby waves, which contain much more of the available energy and must be treated with higher accuracy, enabling a several-fold increase in time step without degrading the accuracy of the solution. The method is presented in Cartesian and spherical coordinates for a rotating earth, using generalized leapfrog, frozen coefficient, and Fourier filtering finite difference schemes. Computational results verify the numerical stability of the approach.
Attention is given to numerical methods and tests of a series of gravitational hydrodynamics computer codes constructed in order to numerically follow the dynamic collapse of spherically symmetric (1-D), axisymmetric (2-D), and non-axisymmetric (3-D) isothermal interstellar clouds. A spherical harmonic expansion is used to solve the Poisson equation for the gravitational potential. It is shown that the use of explicit donor-cell hydrodynamics on a moving spherical coordinate grid ensures mass and momentum conservation and allows the grid to follow the collapse of the fluid. Finally, the performance of the codes is examined.
It is noted that a satellite theory, based on extended phase space and on the true anomaly, was introduced by Scheifele (1970). In the present paper a simple canonical transformation is shown that makes the transition from the classical Delaunay elements to the Scheifele variables. It is stressed that neither spherical coordinates nor Hamilton-Jacobi theory is used. Finally, attention is given to the meaning of the new variables, especially the use of the true anomaly as one of the variables.
A mathematical representation was developed for the electromagnetic force field, the flow field, the temperature field (and for transport controlled kinetics), in a levitation melted metal droplet. The technique of mutual inductances was employed for the calculation of the electromagnetic force field, while the turbulent Navier - Stokes equations and the turbulent convective transport equations were used to represent the fluid flow field, the temperature field and the concentration field. The governing differential equations, written in spherical coordinates, were solved numerically. The computed results were in good agreement with measurements, regarding the lifting force, and the average temperature of the specimen and carburization rates, which were transport controlled.
The equations describing two-dimensional three-component magnetohydrodynamic (MHD) transient flows are formulated for a system of spherical coordinates. With the numerical code based on Implicit Continuous Fluid Eulerian (ICE) scheme, MHD flows resulting from a sudden energy release in a stratified medium are examined. Because of the inclusion of out-of-plane components of velocity and magnetic fields, MHD transverse waves are observed in addition to fast, slow and entropy waves. Numerical results for compressible MHD shocks are found in satisfactory agreement with the theoretical predictions.
In the present investigation, the 'thin layer' Navier-Stokes equations are used to compute the flow about a delta wing and a strake-delta wing configuration. Both configurations possess blunt noses and rounded leading edges. The computational grid about these configurations is generated using a newly developed three-dimensional grid-generation code which solves a set of Poisson equations with spherical coordinate variables using an alternating direction implicit (ADI) scheme. Computational results are obtained for an isolated wing with 60 deg sweep and a strake-wing configuration with 80-60 deg sweep. Computations for angles of attack in the range from 6 to 30 deg are presented. Attention is given to the effect of the angle of attack, the Mach number, and the Reynolds number. Theoretical results are compared with experimental data.
PATRAN-G is an interactive finite element generating and graphics display program. The program is currently executable on the VAX or PRIME computers and is compatible with 19 commercially available graphics terminals. The PATRAN program can be used in two modes: generation of a finite element model, or postprocessing of output data from a finite element analysis program. The finite element generating mode is divided into two phases. Phase I is the geometry construction. During this phase, the geometric boundaries of the model are defined and the material properties are input. Mirror imaging, translation and scaling of lines, surfaces, and volumes, and rotation of regions permit the geometry to be described efficiently. Points may be defined in alternate rectangular, cylindrical, or spherical coordinate systems. Also, PATRAN-G can calculate the intersection of lines or surfaces. This capability allows a user to easily model, for example, the intersection of a cylinder which cuts through a sphere at an angle. Once the geometry of the model has been described in terms of points, lines, surfaces, and volumes. Phase II is used to generate the computer finite element model.
The energy source of stratospheric planetary waves is largely due to vertical propagation from the troposphere, but the stratospheric zonal mean state which changes with time is not stable against small perturbations, whereby transient perturbation waves are generated at the expense of the stratospheric zonal mean available potential energy. This theoretical conclusion is obtained based on the quasi-geostrophic model in spherical coordinates which calculates the characteristic solutions for the basic zonal mean state changing with time.
An algorithm for generating computational grids about arbitrary three-dimensional bodies is developed. The elliptic partial differential equation (PDE) approach developed by Steger and Sorenson and used in the NASA computer program GRAPE is extended from two to three dimensions. Forcing functions which are found automatically by the algorithm give the user the ability to control mesh cell size and skewness at boundary surfaces. This algorithm, as is typical of PDE grid generators, gives smooth grid lines and spacing in the interior of the grid. The method is applied to a rectilinear wind-tunnel case and to two body shapes in spherical coordinates.
In order to determine features of rotating fluid flow that are dependent on the geometry, rotating annulus-type experiments are carried out with a numerical model in spherical coordinates. Rather than constructing and testing a model expressly for this purpose, it is found expedient to modify an existing general circulation model of the atmosphere by removing the model physics and replacing the lower boundary with a uniform surface. A regime diagram derived from these model experiments is presented; its major features are interpreted and contrasted with the major features of rotating annulus regime diagrams. Within the wave regime, a narrow region is found where one or two zonal wave numbers are dominant. The results reveal no upper symmetric regime; wave activity at low rotation rates is thought to be maintained by barotropic rather than baroclinic processes.
The NCOREL computer code was established to compute supersonic flow fields of wings and bodies. The method encompasses an implicit finite difference transonic relaxation method to solve the full potential equation in a spherical coordinate system. Two basic topic to broaden the applicability and usefulness of the present method which is encompassed within the computer code NCOREL for the treatment of supersonic flow problems were studied. The first topic is that of computing efficiency. Accelerated schemes are in use for transonic flow problems. One such scheme is the approximate factorization (AF) method and an AF scheme to the supersonic flow problem is developed. The second topic is the computation of wake flows. The proper modeling of wake flows is important for multicomponent configurations such as wing-body and multiple lifting surfaces where the wake of one lifting surface has a pronounced effect on a downstream body or other lifting surfaces.
An algorithm for implementing a standard least squares procedure on a small computer is presented for use in carrying out a plate reduction to obtain spherical coordinates for comets.
The structure of unstable planetary waves is computed by a quasi-geostrophic model extending from the surface up to 80 km by means of eigenvalue-eigenfunction techniques in spherical coordinates. Three kinds of unstable modes of distinct phase speeds and vertical structures are identified in the winter climate state: (1) the deep Green mode with its maximum amplitude in the stratosphere; (2) the deep Charney mode with its maximum amplitude in the troposphere: and (3) the shallow Charney mode which is largely confined to the troposphere. Both the Green mode and the deep Charney mode are characterized by very slow phase speeds. They are mainly supported by upward wave energy fluxes, but the local baroclinic energy conversion within the stratosphere also contributes in supporting these deep modes. The mesosphere and the troposphere are dynamically independent in the summer season decoupled by the deep stratospheric easterly. The summer mesosphere supports the easterly unstable waves 1-4. Waves 3 and 4 are identified with the observed mesospheric 2-day wave and 1.7-day wave, respectively.
A transformation of the differential equations of motion of the two-body problem in the spherical coordinates to oscillator form is derived. It is shown that the independent variable transformation dt/ds = r squared is a transformation which makes the oscillator form possible.
In the case of an establishment of theoretical models of the hydromagnetic solar wind, the inclusion of the effects of the magnetic field in the solar wind makes it extremely dificult to solve the mathematical problem. This paper has the objective to present a set of particular analytic solutions. The general formulation of Tsinganos (1982) is used to identify a class of analytic solutions to the equations of steady hydromagnetic flows in spherical coordinates. Flow in an open magnetic field are studied, taking into account the problem in dimensionless form, the special case of radial flows with alpha = 0, general radial flows, illustrative examples for flows in which alpha is not equal to 0, a parametric study of nonradial flows in which alpha is not equal to zero, variations in the parameter nu, and variations in the initial speed eta.