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.
The transonic aspect of helicopter flow analysis is addressed. The equations of motion and their implementations are examined, and the computation of real rotor flows is considered. Nonlifting rotor flows, high-speed hover, high advance ratio lifting rotor flows, and strong blade/vortex interaction computations are discussed.
A finite-element model of the Navier-Stokes equations is used for numerical simulation of flow past two normal flat plates arranged side by side at Reynolds number 80 and 160. The results from this simulation indicate that when the gap between the plates is twice the width of a single plate, the individual wakes of the plates behave independently, with the antiphase vortex shedding being dominant. At smaller gap sizes, the in-phase vortex shedding, with strong wake interaction, is favored. The gap flow in those cases becomes biased, with one of the wakes engulfing the other. The direction of the biased flow was found to be switching at irregular intervals, with the time histories of the indicative flow parameters and their power spectra resembling those of a chaotic system.
A 3-D axisymmetric hypersonic engine inlet was investigated using PARC2D, an ideal gas Computational Fluid Dynamics code. The code was used to predict the results of tests conducted in the Rensselaer Polytechnic Institute Hypersonic Shock Tunnel which measured surface and pitot pressures, and shock positions (through Schlieren photography) at freestream Mach numbers of 10, 13, and 15. A strong viscous/shock interaction was observed in both the experiment and the CFD results, due to the model's parabolic compression ramp. Good agreement was found between the experimental results and the CFD solution both for surface pressures and shock positions. Agreement between pitot pressures was less reliable.
The lunar telescope project was an excellent model for the CSTI study because a telescope is a very sensitive instrument, and thermal expansion or mechanical vibration of the mirror assemblies will rapidly degrade the resolution of the device. Consequently, the interactions are strongly coupled. The lunar surface experiences very large temperature variations that range from approximately -180 C to over 100 C. Although the optical assemblies of the telescopes will be well insulated, the temperature of the mirrors will inevitably fluctuate in a similar cycle, but of much smaller magnitude. In order to obtain images of high quality and clarity, allowable thermal deformations of any point on a mirror must be less than 1 micron. Initial estimates indicate that this corresponds to a temperature variation of much less than 1 deg through the thickness of the mirror. Therefore, a lunar telescope design will most probably include active thermal control, a means of controlling the shape of the mirrors, or a combination of both systems. Historically, the design of a complex vehicle was primarily a sequential process in which the basic structure was defined without concurrent detailed analyses or other subsystems. The basic configuration was then passed to the different teams responsible for each subsystem, and their task was to produce a workable solution without requiring major alterations to any principal components or subsystems. Consequently, the final design of the vehicle was not always the most efficient, owing to the fact that each subsystem design was partially constrained by the previous work. This procedure was necessary at the time because the analysis process was extremely time-consuming and had to be started over with each significant alteration of the vehicle. With recent advances in the power and capacity of small computers, and the parallel development of powerful software in structural, thermal, and control system analysis, it is now possible to produce very detailed analyses of intermediate designs in a much shorter period of time. The subsystems can thus be designed concurrently, and alterations in the overall design can be quickly adopted into each analysis; the design becomes an iterative process in which it is much easier to experiment with new ideas, configurations, and components. Concurrent engineering has the potential to produce efficient, highly capable designs because the effect of one subystem on another can be assessed in much more detail at a very early point in the program. The research program consisted of several tasks: scale a prototype telescope assembly to a 1 m aperture, develop a model of the telescope assembly by using finite element (FEM) codes that are available on site, determine structural deflections of the mirror surfaces due to the temperature variations, develop a prototype control system to maintain the proper shape of the optical elements, and most important of all, demonstrate the concurrent engineering approach with this example. In addition, the software used for the finite element models and thermal analysis was relatively new within the Program Development Office and had yet to be applied to systems this large or complex; understanding the software and modifying it for use with this project was also required. The I-DEAS software by Structural Dynamics Research Corporation (SDRC) was used to build the finite element models, and TMG developed by Maya Heat Transfer Technologies, Ltd. (which runs as an I-DEAS module) was used for the thermal model calculations. All control system development was accomplished with MATRIX(sub X) by Integrated Systems, Inc.
Numerical flowfields around a scramjet inlet model are simulated and analyzed. The present inlet flowfield is characterized by thick boundary-layer ingestion and strong viscous/inviscid interaction because of a combined effect of high hypersonic freestream Mach and low Reynolds numbers. Shock-induced separation further enlarges regions of viscous flows which occupy most of the inlet flowfield. Results obtained from the computations with the PARC code developed for ideal gas are presented for several 2D cases at various hypersonic Mach numbers ranging from 10 to 25, and two 3D simulations at Mach numbers of 12 and 19 are also discussed. Comparison between computation and experiment is made in terms of pressure distributions at the wall center line. Large discrepancy is observed and may be partially attributed to the lack of real gas and/or 3D effects in the simulation as well as to the uncertainty of the experiment.
The study of ecosystem function and the contributions of particular ecosystems to global elemental cycling will require data collection from orbital platforms. Presently used ground sampling techniques are not adequate for large area vegetation assessments. The research was directed toward an understanding of the relationship between spectral radiance and plant canopy biomass. Changes in the architecture of the canopy alters the observed spectral radiance patterns. The green components of a plant canopy interacts most strongly with incident radiation and usually reflect edaphic conditions as they relate to plant vigor. Plant biomass serves as an excellent integrator of the environmental factors encountered by the plant and can be used as an indicator for predicting soil redox, nutrient status, microbial activity, etc. The use was examined of TM spectral data to discriminate biomass classes and, in turn, to relate the biomass classes to the chemical status of the soil. Several preliminary studies have indicated that solar angle can be used to improve biomass predictions. This important observation also suggests that latitudinal extrapolation of the biomass predicting models will be possible if selected characteristics of the incident radiation are known.
Turbulent thermal convection is of considerable importance in fluid dynamical transport phenomena occurring, for example, in the planetary boundary layer of the Earth, the interiors of stars, and accretion disks. In particular, during a significant portion of the evolutionary phase of many stars having convectively unstable cores or outer envelopes, a substantial fraction of energy is transported from the central layers to the outer layers by thermal convection. Moreover, as much of the interior of a star is in highly turbulent motion, a complete theory of stellar structure and evolution requires the explicit consideration of turbulence in order to have expressions for the turbulent quantities arising in the stellar structure equations, and particularly, the turbulent fluxes that appear in the total flux conservation equation, such as the convective flux, kinetic energy flux, etc. A reliable quantification of these fluxes continues to present a challenge in astrophysical fluid dynamics, primarily because astrophysical turbulence is almost always fully-developed and nearly inviscid, and therefore governed by strong nonlinear interactions that distribute the energy among a very wide spectrum of eddies with scales ranging from the characteristic dimension of the flow to those sufficiently small to be affected by viscous dissipation. Furthermore, astrophysical flows are invariably compressible, anisotropic, and inhomogeneous, which requires the consideration of the dynamics of longitudinal modes and their interaction with the transverse modes, as well as complicated boundary conditions. In order to reach a compromise between analytical and numerical tractability and the basic physics of turbulent convection, we have constructed a model of stationary turbulent convection that yields various turbulence statistics, including the convective flux, that are required in stellar evolution models.
This study examines the influence of ambient vertical wind shear on microburst intensity, asymmetry, and translation. Results show that microburst asymmetry is influenced by the magnitude of the low-level ambient vertical shear. The microburst outflow elongates in the direction of the shear vector (which is not necessarily in the direction of translation), and generates the greatest hazard (for commercial jet transports) along paths orthogonal to the shear vector. The model results also show that the asymmetry increases with increasing shear magnitude. One implication of these results concerns the detection of a microburst by a ground-based doppler systems. These systems may underestimate the hazard for landing and departing aircraft that are on trajectories orthogonal to both the sensor beam and shear vector, especially if the magnitude of the shear is large. Another implication is that microburst are more likely to be asymmetrical in regions (seasons) where there is climatologically a significant low-level shear. The model results also show that the rotor microbursts and severe wind damage can be a product of the microburst interaction with strong ambient wind shear.
The high resolution (HiRes) algorithm has been used to analyze the far infrared emission of shocked gas and dust in supernova remnants. In the case of supernova remnant IC 443, we find a very good match between the resolved features in the deconvolved images and the emissions of shocked gas mapped in other wavelengths (lines of H2, CO, HCO+, and HI). Dust emission is also found to be surrounding hot bubbles of supernova remnants which are seen in soft X-ray maps. Optical spectroscopy on the emission of the shocked gas suggests a close correlation between the FIR color and local shock speed, which is a strong function of the ambient (preshock) gas density. These provide a potentially effective way to identify regions of strong shock interaction, and thus facilitate studies of kinematics and energetics in the interstellar medium.
Recent studies of low-Z molecular materials including hydrogen to multimegabar pressures (less than 300 GPa) have uncovered a range of phenomena relevant to understanding the nature of the interiors of the outer planets and their satellites. Synchrotron x ray diffraction measurements (to 42 GPa) have been used to determine the crystal structure of the solid (hexagonal-close packed) and equation of state. Sound velocities in fluid and solid hydrogen (to 24 GPa) have been inverted to obtain elastic constants and aggregate bulk and shear moduli. In addition, an improved intermolecular potential has been determined which fits both static and shock-wave data. Use of the new potential for the molecular envelope of Jupiter suggests the need for major revisions of existing Jovian models or a reanalysis of reported free oscillations for the planet. Studies at higher pressures (greater than 100 GPa) reveal a sequence of pressure-induced symmetry-breaking transitions in molecular hydrogen, giving rise to three high-pressure phases (1, 2, and 3). Phase 1 is the rotationally disordered hcp phase which persists from low pressure to well above 100 GPa at high temperature (e.g., 300 K). Phase 2 is a low-temperature, high-pressure phase (transition at 100 GPa and 77 K in H2) with spectral features indicative of partial rotational ordering and crystallographic distortion. The transition to Phase 3 at 150 GPa is accompanied by a weakening of the molecular bond, gradual changes in orientational ordering, strong enhancement of the infrared intramolecular vibrational absorption, and strong intermolecular interactions similar to those of ambient-pressure network solids. Studies of the phase diagram reveal a triple point near 130 K and 160 GPa. Higher pressure measurements of vibrational spectra place a lower bound of approximately 250 GPa on the predicted transition pressure for dissociation of molecular hydrogen to form a monatomic metal.
As a compact binary inspirals due to the emission of gravitational waves, its orbital period decreases continuously down to approximately 1 ms, its value at coalescence. During the last part of the inspiral, the two stars are close together, and their tidal interactions become strong. Neutron stars have many normal modes (core g-modes, crustal discontinuity modes, shear modes, etc.) whose periods lie in the range (approximately several ms) swept by the orbital period. Some of these modes are resonantly excited by the tidal force. The amount of energy a mode absorbs is proportional to the square of the overlap integral between its displacement field and the tidal force field. For all modes of interest, this overlap is poor, resulting in relatively weak excitation. For the best case, the absorbed energy is only a small fraction (approximately 10(exp -6)) of the orbital energy, so the orbital phase shift is too weak to be detected by observations of the gravitational wave signal emitted by the inspiraling binary. However, with displacement amplitudes of excited quadrupole modes ranging up to 0.5% of the stellar radius, the possibility of a detectable electromagnetic signature cannot be dismissed. Both the periods of the modes and the energy they absorb depend quite strongly on the internal structure of the star. Their observation could shed light on the correct high-density equation of state.
Analytical and numerical studies are to be carried out to examine time-averaged thermal effects which are induced by the interaction of strong acoustic fields with a rigid boundary (thermoacoustic streaming). Also of interest is the significance of a second-order thermal expansion coefficient that emerges from this analysis. The model problem to be considered is that of a sphere that is acoustically levitated such that it is effectively isolated in a high-intensity standing acoustic field. The solution technique involves matched asymptotic analysis along with numerical solution of the boundary layer equations. The objective of this study is to predict the thermoacoustic streaming behavior and fully understand the role of the associated second-order thermodynamic modulus.
In spite of significant advancement in the field of high speed computing, flow calculations involving complex geometries and/or flow behavior still require large amounts of CPU time and memory. In order to predict such flows without sacrificing grid convergence and accuracy, adaptive gridding techniques that provide optimal resolution are highly desirable. The present work combines multigrid techniques and domain decomposition concepts to provide local, solution adaptive, grid refinement. Several viscous compressible and incompressible, two and three-dimensional, flows with strong inviscid interaction and/or axial flow reversal, are considered with a segmented multigrid domain decomposition (SMGDD) procedure for which uniform meshes result in each domain. A pressure-based form of flux-vector splitting is applied to the Navier-Stokes equations, which are represented by an implicit lowest-order reduced Navier-Stokes (RNS) system and a purely diffusive, higher-order, deferred-corrector. A trapezoidal or box-like form of discretization insures that all mass conservation properties are satisfied at interfacial and outflow boundaries, even for this primitive-variable non-staggered grid computation. The SMGDD technique presented herein has previously been applied for incompressible two dimensional flows. The present work offers improvement in the gridding strategy, by allowing for disjoint subdomains that provide optimal resolution of disparate flow features. It also extends the SMGDD technique to three dimensional compressible flows. Laminar and turbulent flow in a backward facing step channel is considered; although the procedure is applicable to more severe geometries. The standard K-epsilon model is applied for turbulence closure. For Re greater than 400, differences between two-dimensional theory and experiment are resolved through a three dimensional simulation, which confirms the experimentally observed three dimensionality of the recirculation patterns on the upper and lower surfaces.
Year-round satellite records of sea ice distribution now extend over more than two decades, providing a valuable tool to investigate related characteristics and circulations in the Southern Ocean. We have studied a variety of features indicative of oceanic and atmospheric interactions with Antarctic sea ice. In the Amundsen & Bellingshausen Seas, sea ice extent was found to have decreased by approximately 20% from 1973 through the early 1990's. This change coincided with and probably contributed to recently warmer surface conditions on the west side of the Antarctic Peninsula, where air temperatures have increased by approximately 0.5 C/decade since the mid-1940's. The sea ice decline included multiyear cycles of several years in length superimposed on high interannual variability. The retreat was strongest in summer, and would have lowered the regional mean ice thickness, with attendant impacts upon vertical heat flux and the formation of snow ice and brine. The cause of the regional warming and loss of sea ice is believed to be linked to large-scale circulation changes in the atmosphere and ocean. At the eastern end of the Weddell Gyre, the Cosmonaut Polyna revealed greater activity since 1986, a recurrence pattern during recent winters and two possible modes of formation. Persistence in polynya location was noted off Cape Ann, where the coastal current can interact more strongly with the Antarctic Circumpolar Current. As a result of vorticity conservation, locally enhanced upwelling brings warmer deep water into the mixed layer, causing divergence and melting. In the Ross Sea, ice extent fluctuates over periods of several years, with summer minima and winter maxima roughly in phase. This leads to large interannual cycles of sea ice range, which correlate positively with meridinal winds, regional air temperatures and subsequent shelf water salinities. Deep shelf waters display considerable interannual variability, but have freshened by approximately 0.03/decade since the late 1950's. That could have slowed the thermohaline circulation beneath the Ross Ice Shelf and the properties or volume of local bottom water production.
Lead iodide (PbI2) clusters were synthesized from the chemical reaction of NaI (or KI) with Pb(NO3)2 in H2O, D2O, CH3OH, and C3H7OH media. The observation of the absorption features above 350 nm with the help of integrating sphere accessory strongly suggests the quantum dot formation of PbI2 in solution. Spectral comparison between the synthesized PbI2 clusters in solution and PbI2 nanophase by impregnation of PbI2 in four different pore-sized porous silica indicates that the PbI2 cluster size in solution is less than 2.5 nm in lateral dimension. Atomic force microscopy (AFM) measurements show that the PbL clusters deposited onto three different molecularly flat surfaces are single-layered. The measured height is 1.0 - 0.1 nm. The swollen layer thickness can be attributed to the intralayer contraction from the strong lateral interaction among PbI2 molecules, which is supported by ab initio calculation. Raman scattering measurement of LO and TO modes of PbI2 in bulk and in the confined state were also conducted in 50-150 cu cm region. The observed three bands at 74, %, 106 1/cm are assigned to TO2, LO2, and LO, mode, respectively. The relatively small red-shift in LO modes may be caused by the surface phonon polaritons of PbI2 nanophase in the porous silica.
Transient disturbances in the solar wind initiated by coronal eruptions have been modeled for many years, beginning with the self-similar analytical models of Parker and Simon and Axford. The first numerical computer code (one-dimensional, gas dynamic) to study disturbance propagation in the solar wind was developed in the late 1960s, and a variety of other codes ranging from simple one-dimensional gas dynamic codes through three-dimensional gas dynamic and magnetohydrodynamic codes have been developed in subsequent years. For the most part, these codes have been applied to the problem of disturbances driven by fast CMEs propagating into a structureless solar wind. Pizzo provided an excellent summary of the level of understanding achieved from such simulation studies through about 1984, and other reviews have subsequently become available. More recently, some attention has been focused on disturbances generated by slow CMEs, on disturbances driven by CMEs having high internal pressures, and disturbance propagation effects associated with a structured ambient solar wind. Our purpose here is to provide a brief tutorial on fluid aspects of solar wind disturbances derived from numerical gas dynamic simulations. For the most part we illustrate disturbance evolution by propagating idealized perturbations, mimicking different types of CMEs, into a structureless solar wind using a simple one-dimensional, adiabatic (except at shocks), gas dynamic code. The simulations begin outside the critical point where the solar wind becomes supersonic and thus do not address questions of how the CMEs themselves are initiated. Limited to one dimension (the radial direction), the simulation code predicts too strong an interaction between newly ejected solar material and the ambient wind because it neglects azimuthal and meridional motions of the plasma that help relieve pressure stresses. Moreover, the code ignores magnetic forces and thus also underestimates the speed with which pressure disturbances propagate in the wind.
An initial investigation was conducted to determine the effect of platinum addition on the activities of aluminum and nickel in beta-NiAl(Pt) over the temperature range 1354 to 1692 K. These measurements were made with a multiple effusion-cell configured mass spectrometer (multi-cell KEMS). The results of this study show that Pt additions act to decreased alpha(Al) and increased the alpha(Ni) in beta-NiAl(Pt) for constant X(sub Ni)/X(sub Al) approx. = 1.13, while at constant X(sub Al) the affect of Pt on Al is greatly reduced. The measured partial enthalpies of mixing indicate Al-atoms have a strong self interaction while Ni- and Pt-atoms in have similar interactions with Al-atoms. Conversely the binding of Ni-atoms in beta-NiAl decreases with Pt addition independent of Al concentration. These initial results prove the technique can be applied to the Ni-Al-Pt system but more activity measurements are required to fully understand the thermodynamics of this system and how Pt additions improved the scaling behavior of nickel-based superalloys. In addition, with the choice of a suitable oxide material for the effusion-cell, the "closed" isothermal nature of the effusion-cell allows the direct investigation of an alloy-oxide equilibrium which resembles the "local-equilibrium" description of the metal-scale interface observed during high temperature oxidation. It is proposed that with an Al(l) + Al2O3(s) experimental reference state together with the route measurement of the relative partial-pressures of Al(g) and Al2O(g) allows the activities of O and Al2O3 to be determined along with the activities of Ni and Al. These measurements provide a direct method of investigating the thermodynamics of the metal-scale interface of a TGO-scale.
We report data from the Cassini radio and plasma wave instrument during the approach and first orbit at Saturn. During the approach, radio emissions from Saturn showed that the radio rotation period is now 10 hours 45 minutes 45 k 36 seconds, about 6 minutes longer than measured by Voyager in 1980 to 1981. In addition, many intense impulsive radio signals were detected from Saturn lightning during the approach and first orbit. Some of these have been linked to storm systems observed by the Cassini imaging instrument. Within the magnetosphere, whistler-mode auroral hiss emissions were observed near the rings, suggesting that a strong electrodynamic interaction is occurring in or near the rings.