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 current status of computational methods for unsteady aerodynamics is reviewed. The need to match the fluid dynamic flow equation level to the complexity of the type of unsteady flow under consideration is discussed. Comparisons of computational predictions with experimental unsteady pressures and flutter boundaries are presented. The treatment of complex aircraft geometries is also described.
NASA technical memorandum compares results of upwind and symmetric shock-capturing methods in numerical simulation of gas-dynamic flows. Both methods find shocks as sharp variations in fluid properties over few grid points. Methods differ in type of artificial viscosity introduced to stabilize computations. Symmetric method shown as accurate as upwind method, but with fewer and simpler time steps in transient case or fewer iteration steps in steady-state case.
The heavy rain effects program is presented in the form of the view-graphs. The following topics are covered: rain effects on airfoil performance; two-phase flow dynamics; wind tunnel test results; issues; large-scale results; and summary.
In many previous studies, it has been assumed that the streamlines of the solar wind are in a purely radial direction beyond a few solar radii. Therefore, the nonradial modulations by the magnetic force and gas pressure is neglected at large heliocentric distances. The work reported in this paper includes a two-dimensional magnetohydrodynamic (MHD) study of the latitudinal structure of solar wind flows. The dynamical effects of the interplanetary magnetic field (IMF) are investigated by MHD solutions and by comparing them with both one-dimensional and two-dimensional hydrodynamic (HD) solutions. The corotational effect of the sun and the influence of transient events in the solar atmosphere are neglected. With such a simplification, the latitudinal structure of the steady solar wind is controlled by the dynamic balance of the magnetic force and the gas pressure. Several different patterns of the evolution of the solar wind latitudinal structure are exhibited in the calculation. The results indicate the existence of a proton number density maximum at the magnetic neutral line, whether or not there is a density maximum or minimum at the inner boundary. The drift motion of magnetic field lines toward the magnetic neutral line enhances the magnetic field strength around the neutral sheet, which may provide a possible explanation of the discrepancy between the measured IMF at 1 AU and that extrapolated from the photospheric magnetic field by the current source-surface modeling.
Dynamic measurements of fluctuating static pressure levels were made using flush mounted high frequency response pressure transducers at eleven locations in the circuit of the National Transonic Facility (NTF) over the complete operating range of this wind tunnel. Measurements were made at test section Mach numbers from 0.2 to 1.2, at pressure from 1 to 8.6 atmospheres and at temperatures from ambient to -250 F, resulting in dynamic flow disturbance measurements at the highest Reynolds numbers available in a transonic ground test facility. Tests were also made independently at variable Mach number, variable Reynolds number, and variable drivepower, each time keeping the other two variables constant thus allowing for the first time, a distinct separation of these three important variables. A description of the NTF emphasizing its flow quality features, details on the calibration of the instrumentation, results of measurements with the test section slots covered, downstream choke, effects of liquid nitrogen injection and gaseous nitrogen venting, comparisons between air and nitrogen, isolation of the effects of Mach number, Reynolds number, and fan drive power, and identification of the sources of significant flow disturbances is included. The results indicate that primary sources of flow disturbance in the NTF may be edge-tones generated by test section sidewall re-entry flaps and the venting of nitrogen gas from the return leg of the tunnel circuit between turns 3 and 4 in the cryogenic mode of operation. The tests to isolate the effects of Mach number, Reynolds number, and drive power indicate that Mach number effects predominate. A comparison with other transonic wind tunnels shows that the NTF has low levels of test section fluctuating static pressure especially in the high subsonic Mach number range from 0.7 to 0.9.
The structure of slow magnetosonic shocks in the low beta regime is analyzed with attention given to ion heating and the effects of waves upstream of the electromagnetic ion/ion cyclotron (EMIIC) instability. Shock formation is assessed by means of three methods - a relaxation method and two based on dynamic flow interactions - to determine the effects of initialization and boundary conditions on the formation. Good solutions are found with the piston method and the similar flow-flow method in which the plasma is injected from two boundaries to form two slow shocks. Plasma parameters and shock normal angle are found to be the key variables dictating the structure of the magnetosonic shocks. Four unique classes of resultant shock structures are described in which classical, steady, or nonsteady behavior is found. The analysis also yields insight into the relationship between EMIIC instability and ion dissipation.
Two-dimensional hybrid simulations with particle ions and fluid electrons are used to calculate the kinetic evolution of the Kelvin-Helmholtz instability for a magnetopauselike configuration. The unidirectional magnetic field is essentially transverse to the plasma flow velocity, which is the most unstable case according to linear theory and models the flow dynamics in the subsolar region of the magnetopause for northward interplanetary magnetic field. We recover effects analogous to those found in MHD simulations, including a mode cascade to longer wavelengths. The boundary layer consists of coherent structure and is not well described by a diffusive process. Isolated structures on the order of the ion gyroradius are formed which can cross the boundary in either direction. We describe how the time evolution of these structures represents transport across boundary layers, and we consider the possible connection of these entities to flux transfer events and other structure seen in the low-latitude boundary layer at the Earth's magnetopause as well as to flux ropes commonly observed near the ionopause of Venus. We also discuss the relation of the hybrid calculations to previous MHD simulations and to observations.
In order to simulate the effects of the impingement of hot exhaust jets of High Performance Aircraft on landing surfaces a multi-disciplinary computation coupling flow dynamics to heat conduction in the runway needs to be carried out. Such simulations, which are essentially unsteady, require very large computational power in order to be completed within a reasonable time frame of the order of an hour. Such power can be furnished by the latest generation of massively parallel computers. These remove the bottleneck of ever more congested data paths to one or a few highly specialized central processing units (CPU's) by having many off-the-shelf CPU's work independently on their own data, and exchange information only when needed. During the past year the first phase of this project was completed, in which the optimal strategy for mapping an ADI-algorithm for the three dimensional unsteady heat equation to a MIMD parallel computer was identified. This was done by implementing and comparing three different domain decomposition techniques that define the tasks for the CPU's in the parallel machine. These implementations were done for a Cartesian grid and Dirichlet boundary conditions. The most promising technique was then used to implement the heat equation solver on a general curvilinear grid with a suite of nontrivial boundary conditions. Finally, this technique was also used to implement the Scalar Penta-diagonal (SP) benchmark, which was taken from the NAS Parallel Benchmarks report. All implementations were done in the programming language C on the Intel iPSC/860 computer.
There is no definitive knowledge of which of several concurrent processes ultimately results in unstable combustion within liquid rocket chambers employing shear coaxial injectors. Possible explanations are a detrimental change in the atomization characteristics due to a decrease in the gas-to-liquid velocity ratio, a change in the gas side injector pressure drop allowing acoustic coupling to the propellant feed system or the disappearance of a stabilizing recirculation region at the base of the LOX post. The aim of this research effort is to investigate these proposed mechanisms under conditions comparable to actual engine operation. Spray characterization was accomplished with flash photography and planar laser imaging to examine the overall spray morphology and liquid jet breakup processes and with a PDPA to quantify the spatial distribution of droplet size and mean axial velocity. A simplified stability model based on the Rayleigh criterion was constructed for the flow dynamics occurring within the chamber and injector to evaluate the potential coupling between the chamber and injector acoustic modes and was supported by high frequency measurements of chamber and injector pressure oscillations. To examine recirculation within the LOX post recess, velocity measurements were performed in the recess region by means of LDV. Present experiments were performed under noncombusting conditions using LOX/GH2 stimulants at pressures up to 4 MPa.
A large set of bow shock crossings (i.e., 1392) observed by 17 spacecraft has been used to explore the three-dimensional shape and location of the Earth's bow shock and its dependence on solar wind and interplanetary magnetic field (IMF) conditions. This study investigates deviations from gas dynamic flow models associated with the magnetic terms in the magnetohydrodynamic (MHD) equations. Empirical models predicting the statistical position and shape of the bow shock for arbitrary values of the solar wind pressure, IMF, and Alfvenic Mach number (M(sub A)) have been derived. The resulting data set has been used to fit three-dimensional bow shock surfaces and to explore the variations in these surfaces with sonic (M(sub S)), Alfvenic (M(sub A)) and magnetosonic (M(sub MS)) Mach numbers. Analysis reveals that among the three Mach numbers, M(sub A) provides the best ordering of the least square bow shock curves. The subsolar shock is observed to move Earthward while the flanks flare outward in response to decreasing M(sub A); the net change represents a 6-10% effect. Variations due to changes in the IMF orientation were investigated by rotating the crossings into geocentric interplanetary medium coordinates. Past studies have suggested that the north-south extent of the bow shock surface exceeds the east-west dimension due to asymmetries in the fast mode Mach cone. This study confirms such a north-south versus east-west asymmetry and quantifies its variation with M(sub S), M(sub A), M(sub MS), and IMF orientation. A 2-7% effect is measured, with the asymmetry being more pronounced at low Mach numbers. Combining the bow shock models with the magnetopause model of Roelof and Sibeck (1993), variations in the magnetosheath thickness at different local times are explored. The ratio of the bow shock size to the magnetopause size at the subpolar point is found to be 1.46; at dawn and dusk, the ratios are found to be 1.89 and 1.93, respectively. The subsolar magnetosheath thickness is used to derive the polytropic index gamma according to the empirical relation of Spreiter et al. (1966). The resulting gamma = 2.3 suggests the empirical formula is inadequate to describe the MHD interaction between the solar wind and the magnetosphere.
Dynamic measurements of fluctuating static pressure levels were taken with flush-mounted, high-frequency response pressure transducers at 11 locations in the circuit of the National Transonic Facility (NTF) across the complete operating range of this wind tunnel. Measurements were taken at test-section Mach numbers from 0.1 to 1.2, at pressures from 1 to 8.6 atm, and at temperatures from ambient to -250 F, which resulted in dynamic flow disturbance measurements at the highest Reynolds numbers available in a transonic ground test facility. Tests were also made by independent variation of the Mach number, the Reynolds number, or the fan drive power while the other two parameters were held constant, which for the first time resulted in a distinct separation of the effects of these three important parameters.
In a effort to address current needs for efficient, air propulsion systems, we have developed some new analytical predictive tools for understanding and alleviating aircraft engine instabilities which have led to accelerated high cycle fatigue and catastrophic failures of these machines during flight. A frequent cause of failure in Jets engines is excessive resonant vibrations and stall flutter instabilities. The likelihood of these phenomena is reduced when designers employ the analytical models we have developed. These prediction models will ultimately increase the nation's competitiveness in producing high performance Jets engines with enhanced operability, energy economy, and safety. The objectives of our current threads of research in the final year are directed along two lines. First, we want to improve the current state of blade stress and aeromechanical reduced-ordered modeling of high bypass engine fans, Specifically, a new reduced-order iterative redesign tool for passively controlling the mechanical authority of shroudless, wide chord, laminated composite transonic bypass engine fans has been developed. Second, we aim to advance current understanding of aeromechanical feedback control of dynamic flow instabilities in axial flow compressors. A systematic theoretical evaluation of several approaches to aeromechanical feedback control of rotating stall in axial compressors has been conducted. Attached are abstracts of two .papers under preparation for the 1998 ASME Turbo Expo in Stockholm, Sweden sponsored under Grant No. NAG3-1571. Our goals during the final year under Grant No. NAG3-1571 is to enhance NASA's capabilities of forced response of turbomachines (such as NASA FREPS). We with continue our development of the reduced-ordered, three-dimensional component synthesis models for aeromechanical evaluation of integrated bladeddisk assemblies (i.e., the disk, non-identical bladeing etc.). We will complete our development of component systems design optimization strategies for specified vibratory stresses and increased fatigue life prediction of assembly components, and for specified frequency margins on the Campbell diagrams of turbomachines. Finally, we will integrate the developed codes with NASA's turbomachinery aeromechanics prediction capability (such as NASA FREPS).
Analysis of the cratering records on the Moon, Mercury, and Mars have shown similar size-frequency distributions of craters produced during the late heavy bombardment of the inner solar system (Strom, 1988). Venus provides a valuable data base of information on the impacting population for more recent time. Because of resurfacing events, the Venusian cratering record has been estimated to be only about 500 million years old, and for the most part is in pristine condition, thereby producing an opportunity to discover the properties (size/velocity distribution) of the objects that recently impacted its surface. The Magellan IR mapping of the Venusian surface has produced an extremely high-quality set of crater topographies. The observed deficit of small craters is qualitatively explained by atmospheric effects on impactor breakup and the retardation effects of pressure on crater formation. Information about resurfacing history and impactor flux population can only be conjectured using arrant approximations for atmospheric effects on crater size scaling, such as assuming the absence of pressure effects or using other ad hoc approximations for this dependence. The recent work by Ivanov et al.(1986; 1992); Phillips et al.(199 1; 1992); Schaber et al. (1992) and others support the notion that atmospheric effects may have strongly influenced the Venusian cratering record. The work reported here looks at the potential synergism of aerodynamic entry and the gas dynamic flow fields that govern during the time scale and in the vicinity of crater formation.
A post-test analysis of a set of inside-diameter/outside-diameter (ID/OD) bidirectional brush seals used in three-port wave rotor tests was undertaken to determine brush bristle and configuration wear, pullout, and rotor coating wear. The results suggest that sharp changes in the pressure profiles were not well reflected in bristle tip configuration patterns or wear. Also, positive-to-negative changes in axial pressure gradients appeared to have little effect on the backing plates. Although the brushes had similar porosities, they had very different unpacked arrays. This difference could explain the departure of experimental data from computational fluid dynamics flow predictions for well-packed arrays at higher pressure drops. The rotor wear led to "car-track" scars (upper and lower wear bands) with a whipped surface between the bands. Those bands may have resulted from bristle stiffening at the fence and gap plates during alternate portions of the rotor cycle. Within the bristle response range the wear surface reflected the pressure distribution effect on bristle motion. No sacrificial metallurgical data were taken. The bristles did wear, with correspondingly more wear on the ID brush configurations than on the OD configurations; the complexity in constructing the ID brush was a factor.
Version 3.1 of the NPARC computational fluid dynamics flow solver introduces a capability to solve unsteady flow on moving multi-block, structured grids with nominally second-order time accuracy. The grid motion is due to segments of the boundary grid that translate and rotate in a rigid-body manner or deform. The grid is regenerated at each time step to accommodate the boundary grid motion. The flow equations and computational models sense the moving grid through the grid velocities, which are computed from a time-difference of the grids at two consecutive time levels. For three-dimensional flow domains, it is assumed that the grid retains a planar character with respect to one coordinate. The application and accuracy of NPARC v3.1 is demonstrated for flow about a flying wedge, rotating flap, a collapsing bump in a duct, and the upstart / restart flow in a variable-geometry inlet. The results compare well with analytic and experimental results.
The north polar residual ice deposits of mars are thought to be relatively young, based on the reported lack of any fresh impact craters in Viking Orbiter images. A handful of possible impact features were identified, but available data were inconclusive. Determining the number and current topographic characteristics of any craters on the polar residual ice surface is important in constraining the surface age, relative importance of polar cap resurfacing processes and ice flow dynamics, and the role of the cap deposits in the global volatile and climate cycles. Subsequent image processing advances and new Mars Orbiter Laser Altimeter (MOLA) data in the north polar region are now a considerable aid in impact feature identification. This study reconsiders the abundant Viking high-resolution polar images along with the new altimetry data. We examine possible impact features, compare their topography with known mars high-latitude impact features, and use their morphology as a constraint on crater modification processes.
In this paper we present an evaluation of critical requirements of organometallic chemical vapor deposition (OMCVD) at elevated pressure for a channel flow reactor in a microgravity environment. The objective of using high pressure is to maintain single-phase surface composition for materials that have high thermal decomposition pressure at their optimum growth temperature. Access to microgravity is needed to maintain conditions of laminar flow, which is essential for process analysis. Based on ground based observations we present an optimized reactor design for OMCVD at high pressure and reduced gravity. Also, we discuss non-intrusive real-time optical monitoring of flow dynamics coupled to homogeneous gas phase reactions, transport and surface processes. While suborbital flights may suffice for studies of initial stages of heteroepitaxy experiments in space are essential for a complete evaluation of steady-state growth.
The purpose of this paper is to review modeling and real-time monitoring by robust methods of reflectance spectroscopy of organometallic chemical vapor deposition (OMCVD) processes in extreme regimes of pressure. The merits of p-polarized reflectance spectroscopy under the conditions of chemical beam epitaxy (CBE) and of internal transmission spectroscopy and principal angle spectroscopy at high pressure are assessed. In order to extend OMCVD to materials that exhibit large thermal decomposition pressure at their optimum growth temperature we have designed and built a differentially-pressure-controlled (DCP) OMCVD reactor for use at pressures greater than or equal to 6 atm. We also describe a compact hard-shell (CHS) reactor for extending the pressure range to 100 atm. At such very high pressure the decomposition of source vapors occurs in the vapor phase, and is coupled to flow dynamics and transport. Rate constants for homogeneous gas phase reactions can be predicted based on a combination of first principles and semi-empirical calculations. The pressure dependence of unimolecular rate constants is described by RRKM theory, but requires variational and anharmonicity corrections not included in presently available calculations with the exception of ammonia decomposition. Commercial codes that include chemical reactions and transport exist, but do not adequately cover at present the kinetics of heteroepitaxial crystal growth.