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.
Numerical calculation of inviscid adiabatic flow field around blunt bodies at hypersonic speeds
Quantitative analyses were developed to the point where numerical results can be easily obtained in the optimal transfer trajectory problem. Analyses of deployment dynamics and control of the paired satellite concept are progressing. Automatic control in space, optimal retrieval transfer trajectories, deployment attitude dynamics and control, and automatic control systems for retrieval are also considered.
The work of several graduate students which relate to the NASA grant is briefly summarized. One student has worked on a detailed analysis of the so-called ijk forms of Gaussian elemination and Cholesky factorization on concurrent processors. Another student has worked on the vectorization of the incomplete Cholesky conjugate method on the CYBER 205. Two more students implemented various versions of Gaussian elimination and Cholesky factorization on the FLEX/32.
Category 1, Problems 1 and 2, Category 2, Problem 2, and Category 3, Problem 2 are solved computationally using the Dispersion-Relation-Preserving (DRP) scheme. All these problems are governed by the linearized Euler equations. The resolution requirements of the DRP scheme for maintaining low numerical dispersion and dissipation as well as accurate wave speeds in solving the linearized Euler equations are now well understood. As long as 8 or more mesh points per wavelength is employed in the numerical computation, high quality results are assured. For the first three categories of benchmark problems, therefore, the real challenge is to develop high quality numerical boundary conditions. For Category 1, Problems 1 and 2, it is the curved wall boundary conditions. For Category 2, Problem 2, it is the internal radiation boundary conditions inside the duct. For Category 3, Problem 2, they are the inflow and outflow boundary conditions upstream and downstream of the blade row. These are the foci of the present investigation. Special nonhomogeneous radiation boundary conditions that generate the incoming disturbances and at the same time allow the outgoing reflected or scattered acoustic disturbances to leave the computation domain without significant reflection are developed. Numerical results based on these boundary conditions are provided.
The author has identified the following significant results. Numerous fracture traces were detected on both the color transparencies and black and white spectral bands. Fracture traces of value to mining hazards analysis were noted on the EREP imagery which could not be detected on either the ERTS-1 or high altitude aircraft color infrared photography. Several areas of mine subsidence occurring in the Busseron Creek area near Sullivan, Indiana were successfully identified using color photography. Skylab photography affords an increase over comparable scale ERTS-1 imagery in level of information obtained in mined lands inventory and reclamation analysis. A review of EREP color photography permitted the identification of a substantial number of non-fuel mines within the Southern Indiana test area. A new mine was detected on the EREP photography without prior data. EREP has definite value for estimating areal changes in active mines and for detecting new non-fuel mines. Gob piles and slurry ponds of several acres could be detected on the S-190B color photography when observed in association with large scale mining operations. Apparent degradation of water quality resulting from acid mine drainage and/or siltation was noted in several ponds or small lakes and appear to be related to intensive mining activity near Sullivan, Indiana.
The critical ionization velocity (CIV) phenomenon involves the transfer of energy from neutrals moving relative to a plasma to electrons. Three mechanisms are discussed and illustrated with numerical simulations. If the elastic electron neutral collision frequency is large (comparable to the electron cyclotron frequency) the electrons can gain enough energy from the neutrals to ionize them. Newly created ions form beams in the plasma which are unstable to lower hybrid waves propagating obliquely to the magnetic field. These waves have an electric field component parallel to the background magnetic field which can accelerate electrons. However, the efficiency of this heating is significantly reduced by electromagnetic effects which occur when the relative velocity exceeds the Alfven speed or when the plasma beta is large. A third mechanism involves the E x B drift of electrons in the secondary electric field due to charge separation at the front of a neutral gas cloud. This E field is opposite to the neutral gas velocity and can be much larger than the polarization field. The electron drift excites the modified two-stream instability which causes rapid electron heating and thus ionization in the front.
Three of Uranus' satellites are involved in an interaction similar to the Laplace relation among the Galilean satellites of Jupiter. However, there are no pair-wise commensurabilities in the Uranus system as there are among the Galilean satellites. Thus analysis by Fourier expansion of disturbing functions becomes too unwieldy for practical calculation. Instead, the unexpanded perturbations are numerically averaged over a time scale short compared to the important variations in the Laplace relation. Perturbations of Miranda's orbit are estimated to be of observable magnitude for reasonable masses of Ariel and Umbriel.
A simulation of the flow about the integrated space shuttle vehicle in ascent mode has been undertaken for various flight conditions using the Chimera composite grid discretization approach. Overset body-conforming grids were used to represent each geometric component, and an implicit approximately factored finite-difference procedure was used to solve the three-dimensional thin-layer Navier-Stokes equations. The computational results have been compared with both wind tunnel and flight test data. Although relatively good agreement is obtained with the experimental data, further refinement and evaluation of numerical error is under way.
The first Tethered-Satellite-System (TSS-1), scheduled for a flight in late 1992, is expected to provide relevant information related to the concept of generating an emf in a 20-km-long (or longer) conducting wire. This paper presents numerical simulations of the electrodynamic interactions between the TSS system and space plasma, using a 2D and 3D models of the system. The 2D case code simulates the motion of a long cylinder past a plasma, which is composed of electrons and H(+) ions. The system is solved by allowing the plasma to flow past the cylinder with an imposed magnetic field. The more complex 3D case is considered to study the dynamics in great detail. Results of 2D simulation show that the interaction of a satellite with plasma flowing perpendicularly to the magnetic field results in an enhancement in the current collection.
It is pointed out that an increase in the understanding of weather-related aircraft accidents can save hundreds of human lives and million of dollars. A better understanding regarding the interaction between aircraft operation and severe weather conditions can be obtained with the aid of flight simulator facilities. It is shown that numerical weather modeling is one of the most precise and cost-effective inputs for flight simulators in the long run. A comprehensive weather modeling system is being developed for the simulation of different scales of atmospheric phenomena. The modeling system utilizes two numerical weather models, including the Mesoscale Atmospheric Simulation system, and the Terminal Area Simulation System.
Three types of CVD tungsten (fluoride-produced, chloride-produced and the combination of the two which is termed duplex) were evaluated to determine their weldability, high temperature strength and structural stability during 5000 hour exposure to temperatures of 1540 C and 1700 C. Each type of CVD tungsten could be successfully electron beam welded but the results for the chloride product were not as satisfactory as those of the other two materials. The high temperature strength behavior of the three materials did not differ greatly. However a large difference was noted for the grain growth behavior of the two basic CVD tungsten materials. Fluoride tungsten was found to be relatively stable while for the most part the grain size of chloride tungsten increased appreciably. The examination of freshly fractured surfaces with a scanning electron microscope revealed numerous bubbles in the fluoride material following its exposure to 1700 C for 5000 hours. Less severe thermal treatments produced relatively few bubbles in this material. Only at certain locations within the chloride material associated with the interruption of tungsten were bubbles noted.
Topographic profiles and surface characteristics of small (5 - 25 km diameter) plains-style shield volcanoes on the eastern Snake River Plain (ESRP) provide a method to evaluate eruptive processes and magmatic evolution on Martian volcanic plains. The ESRP is an ideal place to observe Mars-like volcanic features where hundreds of small monogenetic basaltic shields dominate the volcanic-sedimentary depositional sequence, and numerous planetary analogues are evident: coalescent mafic shields, hydromagmatic explosive eruptions, the interaction of lava flows with surficial water and glacial ice, and abundant eolian sand and loess. Single flows cannot be correlated over great distances, and are spatially restricted. These relations are useful for planetary exploration when inferring volcanic evolutionary patterns in lava plains represented by numerous eruptive vents. High spatial resolution imagery and digital topographic data for Mars from MOC, MOLA, and THEMIS is allowing for improvements in the level of detail of stratigraphic mapping of fields of small (< 25 km in diameter) volcanoes as well as studies of the morphological characteristics of individual volcanoes. In order to compare Mars and Earth volcanic features, elevation data from U.S.G.S. 10-meter digital elevation models (DEMs) and high-precision GPS field measurements are used in this study to generate approx. 20m spacing topographic profiles from which slope and surface morphology can be extracted. Average ESRP flank and crater slopes are calculated using 100 - 200 m spacing for optimum comparison to MOLA data, and to reduce the effects of surface irregularities.
The authors numerically studied tidal effects of a massive perturber on a satellite galaxy. The model consists of a spherical satellite galaxy and a point mass perturber and the encounter is non-penetrating. A wide range of density ratios and eccentricities of the relative orbits have been considered. The disruption of the satellite galaxy has been observed when the numerical value of the fractional change in the energy is greater than two. The changes in the energy and angular momentum show smooth variation in the case of unbound orbits and irregular variation in the bound orbit cases. It is shown that for a constant pericentral distance, increasing the density ratio decreases the tidal effects; and for a given density ratio an increase in the eccentricity decreases the tidal effects.
The relationships between pilot fatigue, pilot deviations, reported incidents, and time of day are examined. A sample of 200 Aviation Safety Reporting System (ASRS) reports were analyzed from 1985 and 200 reports from 1987, plus 100 reports from late 1987 and early 1988 that were selected because of possible association with fatigue. The FAA pilot deviation data and incident data were analyzed in relation to denominator data that summarized the hourly operations (landings and takeoffs of scheduled flights) at major U.S. airports. Using as numerators FAA data on pilot deviations and incidents reported to the FAA, the rates by time of day were calculated. Pilot age was also analyzed in relation to the time of day, phase of flight, and type of incident.
In several papers, the solar cycle (SC) effect in the lower atmosphere has been linked to the Quasi-biennial Oscillation (QBO), which is generated primarily by small-scale gravity waves. Salby and Callaghan (2000) analyzed the observed zonal winds of the QBO over more than 40 years and found that it contains a relatively large SC signature at 20 km. Following up on an earlier 2D study with our global-scale Numerical Spectral Model (NSM), we discuss here a 3D study with the QBO under the influence of the SC. For a SC period of 10 years, the amplitude of the relative variations of radiative forcing is taken to vary for simplicity from 0.2% at the surface to 2% at 50 km to 20% at 100 km and above. Covering a limited time span of 40 years, this model produces in the lower stratosphere a relatively large modulation of the QBO, which appears to be related to the SC and is in qualitative agreement with the observations. Some of the energy in the QBO, confined to low latitudes primarily, is redistributed globally by the meridional circulation and planetary waves presumably, so that a measurable SC modulation is generated in the tropospheric temperatures of the polar regions. Further studies are needed, (1) to determine whether the effect is real and prevails in more extensive simulations and whether the results are robust when shorter integration steps are employed, and (2) to explore the mechanism(s) that may ample the apparent SC influence of the UV radiation extending into the lower atmosphere. Quasi-decadal oscillations, generated internally by the QBO interacting with the seasonal cycles, may interfere with or aid the SC effect.
A numerical model for the prediction of local properties of statistically stationary spray-combusting flows is evaluated by comparison with experimental data. To appraise the relative performance of turbulence models, computations were carried out by the k-epsilon model and the algebraic stress model. The present numerical results show the qualitative agreement with experimental data. In terms of overall local flow properties, the algebraic stress model improves a degree of conformity to the experimental data due to its ability to introduce the nonisotropic turbulence effects. Two swirl numbers are considered to investigate the influence of swirl on the droplet evaporation and trajectories, and the effects of droplet/turbulence interactions in flow properties. It is found that the large swirl produces a higher evaporationn rate, and more intensive turbulent mixing and burning. The discrepancies observed in the results are attributed mainly to uncertainties in the initial spray size and velocity distributions, the droplet/wall impingement interaction, the combustion model with the fast chemistry and the turbulence models dealing with the strong streamline curvature and complex interactions between the dispersed droplets and the continuous gas-phase flows.
The main objective of this investigation was the characterization of soil moisture using imaging radars. In order to accomplish this task, a number of intermediate steps had to be undertaken. In this proposal, the theoretical, numerical, and experimental aspects of electromagnetic scattering from natural surfaces was considered with emphasis on remote sensing of soil moisture. In the general case, the microwave backscatter from natural surfaces is mainly influenced by three major factors: (1) the roughness statistics of the soil surface, (2) soil moisture content, and (3) soil surface cover. First the scattering problem from bare-soil surfaces was considered and a hybrid model that relates the radar backscattering coefficient to soil moisture and surface roughness was developed. This model is based on extensive experimental measurements of the radar polarimetric backscatter response of bare soil surfaces at microwave frequencies over a wide range of moisture conditions and roughness scales in conjunction with existing theoretical surface scattering models in limiting cases (small perturbation, physical optics, and geometrical optics models). Also a simple inversion algorithm capable of providing accurate estimates of soil moisture content and surface rms height from single-frequency multi-polarization radar observations was developed. The accuracy of the model and its inversion algorithm is demonstrated using independent data sets. Next the hybrid model for bare-soil surfaces is made fully polarimetric by incorporating the parameters of the co- and cross-polarized phase difference into the model. Experimental data in conjunction with numerical simulations are used to relate the soil moisture content and surface roughness to the phase difference statistics. For this purpose, a novel numerical scattering simulation for inhomogeneous dielectric random surfaces was developed. Finally the scattering problem of short vegetation cover above a rough soil surface was considered. A general scattering model for grass-blades of arbitrary cross section was developed and incorporated in a first order random media model. The vegetation model and the bare-soil model are combined and the accuracy of the combined model is evaluated against experimental observations from a wheat field over the entire growing season. A complete set of ground-truth data and polarimetric backscatter data were collected. Also an inversion algorithm for estimating soil moisture and surface roughness from multi-polarized multi-frequency observations of vegetation-covered ground is developed.
Output from numerical weather models is often used to simulate observations from remote sensing instruments, for purposes ranging from data assimilation, synthetic retrievals of geophysical quantities, and optimization of observing systems. However, when hydrometeors are present, the level of detail provided by the weather model is generally insufficient to fully constrain the input to the radiative transfer model (RTM), and further assumptions must be made by the RTM user in order to produce synthetic observations. Using a hierarchy of models including cloud-resolving, double-moment, bin microphysical, and ice-habit predicting models, along with scattering properties from the OpenSSP, Atmospheric Radiative Transfer Simulator (ARTS) databases, as well as relatively simple geometries (e.g., cylindrical plates and columns), we demonstrate the spread in synthetic observation output and the extent to which it is reduced when microphysics is more strongly constrained by the model. As an intermediate step, an error budget for the RTM simulations was derived and from that we developed and will describe best practices for observation simulation (e.g., optimal number of hydrometeor size bins, truncation of the particle size distribution, angular resolution of scattering phase function). Some statistical comparisons with observations will also be presented.