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.
Temperature gradients for analyzing surface distortions of magnesium solar collector in ecliptic earth orbit determined by analog computer simulation
The source of audible tones occurring with a coaxial nozzle in a range of low Mach numbers is explored experimentally as well as computationally. The hardware is comprised of an inner and an outer nozzle, without a center-body, that are held together by a set of four struts. With increasing jet Mach number (M(sub j)), first a tone occurred at about 2550 Hz around M(sub j)=0.06. At higher values of M(sub j) a tone at 5200 Hz dominated the noise spectra. The corresponding nondimensional frequency, based on effective thickness of the inner nozzle lip and jet exit velocity, turned out to be about 0.2, a value characteristic of Karmann vortex shedding. Thus, vortex shedding from the inner nozzle lip could be linked to the tones. From a comparison of the acoustic wavelengths and the nozzle dimensions, it was inferred that the vortex shedding excited a one-quarter-wave resonance within the divergent section of either the inner nozzle or the outer nozzle. This led to the generation of the sharp tones.
A fully three-dimensional numerical simulation of vortex breakdown using the unsteady, incompressible Navier-Stokes equations has been performed. Solutions to four distinct types of breakdown are identified and compared with experimental results. The computed solutions include weak helical, double helix, spiral, and bubble-type breakdowns. The topological structure of the various breakdowns as well as their interrelationship are studied. The data reveal that the asymmetric modes of breakdown may be subject to additional breakdowns as the vortex core evolves in the streamwise direction. The solutions also show that the freestream axial velocity distribution has a significant effect on the position and type of vortex breakdown.
The e(N) method is employed with Navier-Stokes calculations to study 3D supersonic boundary layers at the leading-edge region of swept wings with attention given to the conditions for linear stability. Mean-flow profiles are computed as solutions to the Navier-Stokes equations, and the e(N) method is used to identify the onset of the transition. Specific treatment is given to the identification of parameters which affect the transition such as free-stream Reynolds number, wall cooling, and boundary-layer suction levels. The boundary layer can be stabilized to below N=10 levels in all locations by utilizing distributed suction, and enhanced stabilization can also be achieved by employing wall cooling and reduced free-stream Reynolds numbers. This study presents key calculations for corresponding experimental investigations by shedding light on requirements for test conditions and measurement requirements.
Design and development of input/output processor and voter-comparator switch concept of guidance and control computer for space stations - Vol. 4
The effect of free-stream turbulence level on stall conditions in plane diffusers is investigated with a finite-volume numerical method which utilizes a k-epsilon turbulence model or a differential Reynolds stress model for closure. Computed results show qualitative agreement with measured data. The ability to eliminate stall and thereby increase pressure recovery in diffusers at large included angles by utilizing high free-stream turbulence is correctly predicted.
The unsteady, compressible, thin-layer Navier-Stokes equations are used to obtain three-dimensional, asymmetric, vortex-flow solutions around cones and cone-cylinder configurations. The equations are solved using an implicit, upwind, flux-difference splitting, finite-volume scheme. The computational applications cover asymmetric flows around a 5 semi-apex angle cone of unit length at various Reynolds number. Next, a cylindrical afterbody of various length is added to the conical forebody to study the effect of the length of cylindrical afterbody on the flow asymmetry. All the asymmetric flow solutions are obtained by using a short-duration side-slip disturbance.
In the fuel preburner of the Space Shuttle Main Engine, face plate, injector, and baffle erosion have been observed. The observed patterns of erosion suggest that flame attachment to the walls is a contributing factor. To better understand the physical phenomena involved, a portion of the preburner was modeled computationally. The simulated 'preburner' had three two-dimensional jets entering a cavity adjacent to a baffle. The computational model employed the Patankar Spalding algorithm with upwind differencing. The turbulence model was a standard k-epsilon model with wall functions. The effects of incoming boundary conditions on turbulent kinetic energy and dissipation, k and epsilon, was studied. The results indicate a very strong sensitivity to these boundary conditions over certain ranges of values.
The fringe-imaging skin friction (FISF) technique, which was originally developed by D. J. Monson and G. G. Mateer at Ames Research Center and recently extended to 3-D flows, is the most accurate skin friction measurement technique currently available. The principle of this technique is that the skin friction at a point on an aerodynamic surface can be determined by measuring the time-rate-of-change of the thickness of an oil drop placed on the surface under the influence of the external air boundary layer. Lubrication theory is used to relate the oil-patch thickness variation to shear stress. The uncertainty of FISF measurements is estimated to be as low as 4 percent, yet little is known about the effects of surface tension and wall adhesion forces on the measured results. A modified version of the free-surface Navier-Stokes solver RIPPLE, developed at Los Alamos National Laboratories, was used to compute the time development of an oil drop on a surface under a simulated air boundary layer. RIPPLE uses the volume of fluid method to track the surface and the continuum surface force approach to model surface tension and wall adhesion effects. The development of an oil drop, over a time period of approximately 4 seconds, was studied. Under the influence of shear imposed by an air boundary layer, the computed profile of the drop rapidly changes from its initial circular-arc shape to a wedge-like shape. Comparison of the time-varying oil-thickness distributions computed using RIPPLE and also computed using a greatly simplified numerical model of an oil drop equation which does not include surface tension and wall adhesion effects) was used to evaluate the effects of surface tension on FISF measurement results. The effects of surface tension were found to be small but not necessarily negligible in some cases.
A hybrid plasma simulation model is described and applied to the study of electrostatic wave propagation in a one-dimensional Maxwellian plasma with periodic boundary conditions. The model employs a cloud-in-cell scheme which can drastically reduce the fluctuations in particle simulation models and greatly ease the computational difficulties of the Vlasov equation approach. A grid in velocity space is introduced and the particles are represented by points in the x-v phase space. The model is tested first in the absence of an applied signal and then in the presence of a small-amplitude perturbation. The method is also used to study propagation of an essentially monochromatic plane wave. Results on amplitude oscillation and nonlinear frequency shift are compared with available theories.
Outline: Introduction; Launch Ascent and Vehicle Aerodynamics (LAVA) Framework; 1st AIAA Sonic Boom Prediction Workshop Oblique Shock/Plume Interaction; Summary.
This paper presents results from a study conducted to investigate the effect of incident oblique shocks on the effectiveness of a coolant film at Mach numbers, typical of those expected in a scramjet combustor at Mach 15 to 20 flight. Computations with a parabolic code are in good agreement with the measured pressures and heat fluxes, after accounting for the influence of the shock upstream of its point of impingement on the plate, and the expansion from the trailing edge of the shock generator. The test data shows that, for the blowing rates tested, the film is rendered largely ineffective by the shock. Computations show that coolant blowing rates five to ten times those tested are required to protect against shock-induced heating. The implications of the results to scramjet combustor design are discussed.
Axisymmetric and three-dimensional, multi-nozzle plume flows around generic rocket geometries are investigated with a three-dimensional Navier-Stokes solver to study the interactive effects between hard body and the plume. Time-asymptotic, laminar, ideal-gas solutions obtained with a two-factor, flux-split scheme and a diagonal, upwind scheme are presented. Computed solutions to three-dimensional, multi-nozzle problems and single-nozzle, axisymmetric problems demonstrate flow field features including three-dimensionality and hard-body effects. Geometry and three-dimensional effects are shown to be significant in multi-nozzle flows.
A symmetric total variation diminishing (STVD) algorithm has been applied to the solution of the three-dimensional hypersonic flowfield surrounding the Aeroassist Flight Experiment (AFE) vehicle. Both perfect-gas and chemical nonequilibrium models have been used. The perfect-gas flows were computed at two different Reynolds numbers, including a flight trajectory point at maximum dynamic pressure, and on two different grids. Procedures for coupling the solution of the species continuity equations with the Navier-Stokes equations in the presence of chemical nonequilibrium are reviewed and tested on the forebody of the AFE and on the complete flowfield assuming noncatalytic wall and no species diffusion. Problems with the STVD algorithm unique to flows with variable thermodynamic properties (real gas) are identified and algorithm modifications are suggested. A potential heating problem caused by strong flow impingement on the nozzle lip in the near wake at 0-deg angle of attack has been identified.
The present investigation is a continuation of a previous joint project between the Florida State University and the NASA Langley Research Center Liner Physics Team. In the previous project, a study of acoustic liners, in two dimensions, inside a normal incidence impedance tube was carried out. The study consisted of two parts. The NASA team was responsible for the experimental part of the project. This involved performing measurements in an impedance tube with a large aspect ratio slit resonator. The FSU team was responsible for the computation part of the project. This involved performing direct numerical simulation (DNS) of the NASA experiment in two dimensions using CAA methodology. It was agreed that upon completion of numerical simulation, the computed values of the liner impedance were to be sent to NASA for validation with experimental results. On following this procedure good agreements were found between numerical results and experimental measurements over a wide range of frequencies and sound-pressure-level. Broadband incident sound waves were also simulated numerically and measured experimentally. Overall, good agreements were also found.
A viscous-inviscid interaction technique is advocated as both an efficient and accurate means of predicting the performance of two-dimensional thrust augmenting ejectors. The flow field is subdivided into a viscous region that contains the turbulent jet and an inviscid region that contains the ambient fluid drawn into the device. The inviscid region is computed with a higher-order panel method, while an integral method is used for the description of the viscous part. The strong viscous-inviscid interaction present within the ejector is simulated in an iterative process where the two regions influence each other en route to a converged solution. The model is applied to a variety of parametric and optimization studies involving ejectors having either one or two primary jets. The effects of nozzle placement, inlet and diffuser shape, free stream speed, and ejector length are investigated. The inlet shape for single jet ejectors is optimized for various free stream speeds and Reynolds numbers. Optimal nozzle tilt and location are identified for various dual-ejector configurations.
Erosion of carbon due to oxidation does not occur uniformly but through the formation of localized etch pits because of active surface sites. These active sites are formed due to the presence of atomic defects on the carbon surface, and have much higher reactivity compared to average non-defective sites. Thus, these active sites are first to react during ablation, resulting in their removal. This causes all the neighboring atoms to be defective and increase their reactivity, thus leading to the localized carbon removal around these “active” sites. In this manner, these highly reactive defective sites serve as nucleation sites for the formation and growth of etch pits with potentially detrimental effects on the structural integrity. In order to understand the influence of these etch pits on the material properties of carbon fiber microstructures, we have developed a new capability within direct simulation Monte Carlo (DSMC) to capture the etch pit formation process. This capability is developed within the DSMC code SPARTA (Stochastic PArallel Rarefied-gas Time-accurate Analyzer) and can model the material removal in presence of active sites leading to the formation of etch pits. The focus of the current work will be to study the effect of etch pits on the material properties of FiberForm, a commonly used base material within many thermal protection system materials (TPS). The microstructure of virgin FiberForm obtained directly from X-ray microtomography experiments is used within SPARTA to obtain the ablated geometries with etch pits. These pitted microstructures are then imported within the Porous Microstructure Analysis (PuMA) software and various material properties such as elasticity, thermal conductivity, and permeability are computed. The variation of these properties as a result of the complex evolution of the surface topology due to etch pit formation is studied and analyzed. Furthermore, the effect of pitting is compared to the case of shrinking fibers, which has been the standard for modelling ablation of carbon structures; and significant differences are observed. Thus, such a physically realistic modeling of material removal through the formation of etch pits will be helpful in predicting the degradation of carbon-based TPS more accurately during oxidation; as well as other mechanisms such as spallation, which involves the removal of chunks of material into the flow due to etch pit growth. This will ultimately improve our understanding of the failure modes in these materials due to ablation.
Erosion of carbon due to oxidation does not occur uniformly but through the formation of localized etch pits because of active surface sites. These active sites are formed due to the presence of atomic defects on the carbon surface, and have much higher reactivity compared to average non-defective sites. Thus, these active sites are first to react during ablation, resulting in their removal. This causes all the neighboring atoms to be defective and increase their reactivity, thus leading to the localized carbon removal around these “active” sites. In this manner, these highly reactive defective sites serve as nucleation sites for the formation and growth of etch pits with potentially detrimental effects on the structural integrity. In order to understand the influence of these etch pits on the material properties of carbon fiber microstructures, we have developed a new capability within direct simulation Monte Carlo (DSMC) to capture the etch pit formation process. This capability is developed within the DSMC code SPARTA (Stochastic PArallel Rarefied-gas Time-accurate Analyzer) and can model the material removal in presence of active sites leading to the formation of etch pits. The focus of the current work will be to study the effect of etch pits on the material properties of FiberForm, a commonly used base material within many thermal protection system materials (TPS). The microstructure of virgin FiberForm obtained directly from X-ray microtomography experiments is used within SPARTA to obtain the ablated geometries with etch pits. These pitted microstructures are then imported within the Porous Microstructure Analysis (PuMA) software and various material properties such as elasticity, thermal conductivity, and permeability are computed. The variation of these properties as a result of the complex evolution of the surface topology due to etch pit formation is studied and analyzed. Furthermore, the effect of pitting is compared to the case of shrinking fibers, which has been the standard for modelling ablation of carbon structures; and significant differences are observed. Thus, such a physically realistic modeling of material removal through the formation of etch pits will be helpful in predicting the degradation of carbon-based TPS more accurately during oxidation; as well as other mechanisms such as spallation, which involves the removal of chunks of material into the flow due to etch pit growth. This will ultimately improve our understanding of the failure modes in these materials due to ablation.