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.
A three-dimensional numerical method has been developed to analyze the complex flow field inside a turbine volute. Comparisons are made between solutions with different boundary conditions.
Traditionally, Solid Rocket Motor (SRM) internal ballistic performance has been analyzed and predicted with either zero-dimensional (volume filling) codes or one-dimensional ballistics codes. One dimensional simulation of SRM performance is only necessary for ignition modeling, or for motors that have large length to port diameter ratios which exhibit an axial "pressure drop" during the early burn times. This type of prediction works quite well for many types of motors, however, when motor aspect ratios get large, and port to throat ratios get closer to one, two dimensional effects can become significant. The initial propellant grain configuration for the Space Shuttle Reusable Solid Rocket Motor (RSRM) was analyzed with 2-D, steady, axi-symmetric computational fluid dynamics (CFD). The results of the CFD analysis show that the steady-state performance prediction at the initial burn geometry, in general, agrees well with 1-D transient prediction results at an early time, however, significant features of the 2-D flow are captured with the CFD results that would otherwise go unnoticed. Capturing these subtle differences gives a greater confidence to modeling accuracy, and additional insight with which to model secondary internal flow effects like erosive burning. Detailed analysis of the 2-D flowfield has led to the discovery of its hidden 1-D isentropic behavior, and provided the means for a thorough and simplified understanding of internal solid rocket motor flow. Performance parameters such as nozzle stagnation pressure, static pressure drop, characteristic velocity, thrust and specific impulse are discussed in detail and compared for different modeling and prediction methods. The predicted performance using both the 1-D codes and the CFD results are compared with measured data obtained from static tests of the RSRM. The differences and limitations of predictions using ID and 2-D flow fields are discussed and some suggestions for the design of large L/D motors and more critically, motors with port to throat ratios near one, are covered.
An efficient grid-interfacing zonal algorithm was developed for computing the three-dimensional transonic flow field about wing/nacelle configurations. the algorithm uses the full-potential formulation and the AF2 approximate factorization scheme. The flow field solution is computed using a component-adaptive grid approach in which separate grids are employed for the individual components in the multi-component configuration, where each component grid is optimized for a particular geometry such as the wing or nacelle. The wing and nacelle component grids are allowed to overlap, and flow field information is transmitted from one grid to another through the overlap region using trivariate interpolation. This report represents a discussion of the computational methods used to generate both the wing and nacelle component grids, the technique used to interface the component grids, and the method used to obtain the inviscid flow solution. Computed results and correlations with experiment are presented. also presented are discussions on the organization of the wing grid generation (GRGEN3) and nacelle grid generation (NGRIDA) computer programs, the grid interface (LK) computer program, and the wing/nacelle flow solution (TWN) computer program. Descriptions of the respective subroutines, definitions of the required input parameters, a discussion on interpretation of the output, and the sample cases illustrating application of the analysis are provided for each of the four computer programs.
An analytical investigation based on the Mobility Power Flow method is presented for the determination of the vibrational response and power flow for two coupled flat plate structures in an L-shaped configuration, subjected to acoustical excitation. The principle of the mobility power flow method consists of dividing the global structure into a series of subsystems coupled together using mobility functions. Each separate subsystem is analyzed independently to determine the structural mobility functions for the junction and excitation locations. The mobility functions, together with the characteristics of the junction between the subsystems, are then used to determine the response of the global structure and the power flow. In the coupled plate structure considered here, mobility power flow expressions are derived for excitation by an incident acoustic plane wave. In this case, the forces (acoustic pressures) acting on the structure are dependent on the response of the structure because of the scattered pressure component. The interaction between the structure and the fluid leads to the derivation of a corrected mode shape for the plates' normal surface velocity and also for the structure mobility functions. The determination of the scattered pressure components in the expressions for the power flow represents an additional component in the power flow balance for the source plate and the receiver plate. This component represents the radiated acoustical power from the plate structure.
An analytical investigation based in the Mobility Power Flow (MPF) method is presented for the determination of the vibrational response and power flow for two coupled flat plate structures in an L-shaped configuration, subjected to distributed excitation. The principle of the MPF method consists of dividing the global structure into a series of subsystems coupled together using mobility functions. Each separate subsystem is analyzed independently to determine the structural mobility functions for the junction and excitation locations. The mobility functions, together with the characteristics of the junction between the subsystems, are then used to determine the response of the global structure and the MPF. In the considered coupled plate structure, MPF expressions are derived for distributed mechanical excitation which is independent of the structure response. However using a similar approach with some modifications excitation by an acoustic plane wave can be considered. Some modifications are required to deal with the latter case are necessary because the forces (acoustic pressure) acting on the structure are dependent on the response of the structure due to the presence of the scattered pressure.
The Fluid Mechanics and Acoustics Laboratory at Hampton University (HU/FM&AL) jointly with the NASA Glenn Research Center has conducted four connected subprojects under the reporting project. Basically, the HU/FM&AL Team has been involved in joint research with the purpose of theoretical explanation of experimental facts and creation of accurate numerical simulation techniques and prediction theory for solution of current problems in propulsion systems of interest to the NAVY and NASA agencies. This work is also supported by joint research between the NASA GRC and the Institute of Mechanics at Moscow State University (IM/MSU) in Russia under a CRDF grant. The research is focused on a wide regime of problems in the propulsion field as well as in experimental testing and theoretical and numerical simulation analyses for advanced aircraft and rocket engines. The FM&AL Team uses analytical methods, numerical simulations and possible experimental tests at the Hampton University campus. The fundamental idea uniting these subprojects is to use nontraditional 3D corrugated and composite nozzle and inlet designs and additional methods for exhaust jet noise reduction without essential thrust loss and even with thrust augmentation. These subprojects are: (1) Aeroperformance and acoustics of Bluebell-shaped and Telescope-shaped designs; (2) An analysis of sharp-edged nozzle exit designs for effective fuel injection into the flow stream in air-breathing engines: triangular-round, diamond-round and other nozzles; (3) Measurement technique improvement for the HU Low Speed Wind Tunnel; a new course in the field of aerodynamics, teaching and training of HU students; experimental tests of Mobius-shaped screws: research and training; (4) Supersonic inlet shape optimization. The main outcomes during this reporting period are: (l) Publications: The AIAA Paper #00-3170 was presented at the 36th AIAA/ASME/SAE/ASEE Joint Propulsion Conference, 17-19 June, 2000, Huntsville, AL. The AIAA Paper #01-1893 has been accepted for the AIAA/NAL-NASDA-ISAS 10th International Space Planes and Hypersonic Systems and Technologies Conference, 24-27 April 2001, Kyoto, Japan. The AIAA Paper #01 -3204 has been accepted for presentation at the 37th AIAA/ASME/SAE/ASEE Joint Propulsion Conference, being held on 08-11 July, in Salt Lake City, UT; (2) A U.S. patent #6,082,635 was granted on July 4, 2000; (3) Grants and proposals: The H U/ FM&AL was awarded the NASA grant NAG-3-2495 in October 2000 and the laboratory is a primary U.S. research team in a joint project under the CRDF award granted to the NASA GRC and IM/MSU (Russia) in July 2000; (4) Theory and numerical simulations: Analytical theory, numerical simulation, comparison of theoretical with experimental results, and modification of theoretical approaches, models, grids, etc., have been conducted for several complicated 2D and 3D nozzle and inlet designs using NASA, ICASE, and IM/MSU codes based on full Euler and Navier-Stokes solvers: CFL3D, FLUENT, and GODUNOV, and others; (5) Experimental Tests: (a) A new course: "Advanced Aerodynamics and Aircraft Performance" presented in spring semester, 2001; training and experimental test research using the HU LSWT. (b) Small-scale M6bius-shaped screws were tested in different conditions and their application has shown essential benefits by comparison with traditional designs; (6) Installation in the FM&AL computer system: second software TECPLOT 8.0 for the UNIX SGI workstation and free TECPLOT 7.5 for the PC Dell computer, and 2D and 3D GRIDGEN (version 9) for the UNIX SGI as well as installation of two free NASA codes, 3D MAG and VULCAN; (7) Student Research Activity: Involvement of two undergraduate students as research assistants in the current research project.
The research is focused on a wide regime of problems in the propulsion field as well as in experimental testing and theoretical and numerical simulation analyses for advanced aircraft and rocket engines. Results obtained are based on analytical methods, numerical simulations and experimental tests at the NASA LaRC and Hampton University computer complexes and experimental facilities. The main objective of this research is injection, mixing and combustion enhancement in propulsion systems. The sub-projects in the reporting period are: (A) Aero-performance and acoustics of Telescope-shaped designs. The work included a pylon set application for SCRAMJET. (B) An analysis of sharp-edged nozzle exit designs for effective fuel injection into the flow stream in air-breathing engines: triangular-round and diamond-round nozzles. (C) Measurement technique improvements for the HU Low Speed Wind Tunnel (HU LSWT) including an automatic data acquisition system and a two component (drag-lift) balance system. In addition, a course in the field of aerodynamics was developed for the teaching and training of HU students.
A general purpose, three-dimensional computational fluid dynamics code named PHOENICS, developed by CHAM Inc., is used to model the flow in the aft-platform seal cavity in the high pressure fuel pump of the space shuttle main engine. The model is used to predict the temperatures, velocities, and pressures in the cavity for six different sets of boundary conditions. The results are presented as input for further analysis of two known problems in the region, specifically: erratic pressures and temperatures in the adjacent coolant liner cavity and cracks in the blade shanks near the outer diameter of the aft-platform seal.
The heat transfer and pressure drop through turbine vanes made of a sintered, porous metal coated with a thin layer of ceramic and convection cooled by spanwise flow of cooling air were analyzed. The analysis was made to determine the feasibility of using this concept for cooling very small turbines, primarily for short duration applications such as in missile engines. The analysis was made for gas conditions of approximately 10 and 40 atm and 1644 K and with turbine vanes made of felt type porous metals with relative densities from 0.2 to 0.6 and ceramic coating thicknesses of 0.076 to 0.254 mm.
A key component of the small unmanned aircraft systems traffic management ecosystem is the design of scalable algorithms for strategic deconfliction of drones prior to takeoff. In this work, we focus on efficient flow management of drones on a network of intersecting edges subject to two kinds of spacing constraints: 1) between any two adjacent vehicles on an edge and 2) between any two vehicles on two different edges arriving one after the other at an intersection. The spacing is designed to enable non-intersection of operational volumes corresponding to two different vehicles thereby properly separating the vehicles inside each volume. For simplicity, we assume a constant ground speed for the drones and fixed dimensions for the operational volume blocks. The deconfliction is managed by adjusting the takeoff time of the drones, thereby regulating their arrival time at various crossing waypoints in the network. This framework allows us to study the maximum flow (throughput) of vehicles on a network of edges connecting depots to drop off sites subject to the temporal spacing constraints. The departure scheduling of individual drones results in a combinatorial optimization problem. To alleviate this, we solve a max-flow formulation and use queueing theory to simplify the analysis and provide upper bounds to the underlying optimization problem for individual drone departure scheduling. Our results indicate that throughput drops rapidly after the density of drones in the network passes the max-flow limits.
An unsteady linearized formulation based on Oswatitsch-Keune's parabolic method is developed to analyze transonic flow past oscillating slender bodies. In contrast to the widely used integral transform method, it is shown that all solutions can be derived by a simpler method directly in the physical plane. By various expansion procedures, low-frequency solutions then are derived according to two clearly defined frequency ranges. Adams-Sears' iteration is employed to account for the second-order effects. Stability derivatives are compared with available theories and data. It is found that the derivatives depend more sensitively on thickness than on the reduced frequency. Finally, a critical assessment of the present method is given.
A hybrid Lax-Wendroff/Method of Characteristics computer code has been developed for numerical simulation of flow transients associated with the operation of MHD generator facilities. The code employs the shock-fitting method, with an Eulerian formulation of the basic conservation equations and explicit tracking of shock waves. Pressure, temperature, and velocity are used as primary integration variables to simplify interfacing of the code with real-gas thermodynamic and transport property tables. Application of the code to the simulation of selected transients for the AEDC/HPDE MHD generator produced results that are in good agreement with experimental observations.
Three-dimensional unsteady viscous effects can significantly influence the performance of fixed and rotary wing aircraft. These effects are important in both flows about helicopter rotors in forward flight and flows about 3-D (swept and tapered) supercritical wings. A computational procedure for calculating such flow field is developed, and therefore would be of great value in the design process as well as in understanding the corresponding flow phenomena. The procedure is based upon an alternating direction technique employing the Linearized Block Implicit method for solving 3-D viscous flow problems. In order to demonstrate the viability of this method, 2-D and 3-D problems are computed. These include the flow over a 2-D NACA 0012 airfoil under steady and oscillating conditions, and the steady, skewed, 3-D flow on a flat plate. Although actual 3-D flows over wings were not obtained, the ground work was laid for considering such flows. The description of the computational procedure and results are given.
The Level 4 Mission Sequence Test (MST) was studied to develop strategies and recommendations to facilitate information flow. Recommendations developed as a result of this study include revised format of the Test and Assembly Procedure (TAP) document and a conceptualized software based system to assist in the management of information flow during the MST.
A hypersonic-wind-tunnel nozzle concept which incorporates a hot-core flow surrounded by an annular flow of cold air offers a promising technique for maximizing the model size while minimizing the power required to heat the test core. This capability becomes especially important when providing the true-temperature duplication needed for hypersonic propulsion testing. Several two-dimensional wind-tunnel nozzle configurations that are designed according to this concept are analyzed by using recently developed analytical techniques for prediction of the boundary-layer growth and the mixing between the hot and cold coaxial supersonic airflows. The analyses indicate that introduction of the cold annular flow near the throat results in an unacceptable test core for the nozzle size and stagnation conditions considered because of both mixing and condensation effects. Use of a half-nozzle with a ramp on the flat portion does not appear promising because of the thick boundary layer associated with the extra length. However, the analyses indicate that if the cold annular flow is introduced at the exit of a full two-dimensional nozzle, an acceptable test core will be produced. Predictions of the mixing between the hot and cold supersonic streams for this configuration show that mixing effects from the cold flow do not appreciably penetrate into the hot core for the large downstream distances of interest.
A new method is presented for calculating the quasi-steady transonic flow over a lifting or non-lifting rotor blade in both hover and forward flight by using Euler equations. The approach is to solve Euler equations in a rotor-fixed frame of reference using a finite volume method. A computer program was developed and was then verified by comparison with wind-tunnel data. In all cases considered, good agreement was found with published experimental data.
This report summarizes the methods developed for the aerodynamic analysis and the shape optimization of the nozzle-afterbody section of a hypersonic vehicle. Initially, exhaust gases were assumed to be air. Internal-external flows around a single scramjet module were analyzed by solving the 3D Navier-Stokes equations. Then, exhaust gases were simulated by a cold mixture of Freon and Ar. Two different models were used to compute these multispecies flows as they mixed with the hypersonic airflow. Surface and off-surface properties were successfully compared with the experimental data. The Aerodynamic Design Optimization with Sensitivity analysis was then developed. Pre- and postoptimization sensitivity coefficients were derived and used in this quasi-analytical method. These coefficients were also used to predict inexpensively the flow field around a changed shape when the flow field of an unchanged shape was given. Starting with totally arbitrary initial afterbody shapes, independent computations were converged to the same optimum shape, which rendered the maximum axial thrust.
Turbine-based propulsion systems for access to space have been investigated at NASA Glenn Research center. A ground demonstrator engine for validation testing has been developed as a part of the program. The demonstrator, the Revolutionary Turbine Accelerator (RTA-1), is a variable cycle turbofan ramjet designed to transition from an augmented turbofan to a ramjet that produces the thrust required to accelerate the vehicle to Mach 4. The RTA-1 is designed to accommodate a large variation in bypass ratio from sea level static to Mach 4 flight condition. A key component of this engine is a new fan stage that accommodates these large variations in bypass ratio and flow ranges. In the present study, unsteady flow behavior in the fan of the RTA-1 is studied in detail with large eddy simulation (LES) and the numerical results are compared with measured data. During the experimental study of the fan stage, humming sound was detected at 100 % speed near stall operation. The main purpose of the study is to investigate details of the unsteady flow behavior at near stall operation and to identify a possible cause of the hum. The large eddy simulation of the current flow field reproduces main features of the measured flow very well. The LES simulation indicates that non-synchronous flow instability develops as the fan operates toward the stall limit. The FFT analysis of the calculated wall pressure shows that the rotating flow instability has the characteristic frequency that is about 50% of the blade passing frequency.