Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Reacting flows”

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.

At least 1,117 records · Page 62

Development and application of the GIM code for the Cyber 203 computer

The GIM computer code for fluid dynamics research was developed. Enhancement of the computer code, implicit algorithm development, turbulence model implementation, chemistry model development, interactive input module coding and wing/body flowfield computation are described. The GIM quasi-parabolic code development was completed, and the code used to compute a number of example cases. Turbulence models, algebraic and differential equations, were added to the basic viscous code. An equilibrium reacting chemistry model and implicit finite difference scheme were also added. Development was completed on the interactive module for generating the input data for GIM. Solutions for inviscid hypersonic flow over a wing/body configuration are also presented.

Stainaker, J. F.↗

Direct numerical simulations of a reacting mixing layer with chemical heat release

In order to study the coupling between chemical heat release and fluid dynamics, direct numerical simulations of a chemically reacting mixing layer with heat release are performed. The fully compressible equations as well as an approximate set of equations that is asymptotically valid for low-Mach-number flows are treated. These latter equations have the computational advantage that high-frequency acoustic waves have been filtered out, allowing much larger time steps to be taken in the numerical solution procedure. A detailed derivation of these equations along with an outline of the numerical solution technique is given. Simulation results indicate that the rate of chemical product formed, the thickness of the mixing layer, and the amount of mass entrained into the layer all decrease with increasing rates of heat release.

Mcmurtry, P. A.↗

Transverse liquid fuel jet breakup, burning, and ignition

An analytical study of the breakup, burning, and ignition of liquid fuels injected transversely into a hot air stream is conducted. The non-reacting liquid jet breakup location is determined by the local sonic point criterion. Two models, one employing analysis of an elliptical jet cross-section and the other employing a two-dimensional blunt body to represent the transverse jet, were used for sonic point calculations. An auxiliary criterion based on surface tension stability is used as a separate means of determining the breakup location. For the reacting liquid jet problem, a diffusion flame supported by a one-step chemical reaction within the gaseous boundary layer is solved along the ellipse surface in subsonic cross flow. Typical flame structures and concentration profiles were calculated for various locations along the jet cross-section as a function of upstream Mach numbers. The integration reaction rate along the jet cross-section is used to predict ignition position, which is found to be situated near the stagnation point. While a multi-step reaction is needed to represent the ignition process more accurately, the present calculation does yield reasonable predictions concerning ignition along a curved surface.

Li, Hsi-Shang↗

Stability of compressible reacting mixing layer

Linear instability of compressible reacting mixing layers is analyzed with emphasis on the effects of heat release and compressibility. Laminar solutions of the compressible boundary-layer equations are used as the base flows. The parameters of this study are the adiabatic flame temperature, the Mach number of the upper stream, frequency, wavenumber, and the direction of propagation of the disturbance wave. Stability characteristics of the flow are presented. Three groups of unstable modes are found when the Mach number and/or heat release are large. Finally, it is shown that the unstable modes are two-dimensional for large heat release even in highly compressible flow.

Shin, D. S.↗

Orion Service Module Reaction Control System Plume Impingement Analysis Using PLIMP/RAMP2

The Orion Crew Exploration Vehicle Service Module Reaction Control System engine plume impingement was computed using the plume impingement program (PLIMP). PLIMP uses the plume solution from RAMP2, which is the refined version of the reacting and multiphase program (RAMP) code. The heating rate and pressure (force and moment) on surfaces or components of the Service Module were computed. The RAMP2 solution of the flow field inside the engine and the plume was compared with those computed using GASP, a computational fluid dynamics code, showing reasonable agreement. The computed heating rate and pressure using PLIMP were compared with the Reaction Control System plume model (RPM) solution and the plume impingement dynamics (PIDYN) solution. RPM uses the GASP-based plume solution, whereas PIDYN uses the SCARF plume solution. Three sets of the heating rate and pressure solutions agree well. Further thermal analysis on the avionic ring of the Service Module showed that thermal protection is necessary because of significant heating from the plume.

Wang, Xiao-Yen J.↗

Orion Service Module Reaction Control System Plume Impingement Analysis Using PLIMP/RAMP2

The Orion Crew Exploration Vehicle Service Module Reaction Control System engine plume impingement was computed using the plume impingement program (PLIMP). PLIMP uses the plume solution from RAMP2, which is the refined version of the reacting and multiphase program (RAMP) code. The heating rate and pressure (force and moment) on surfaces or components of the Service Module were computed. The RAMP2 solution of the flow field inside the engine and the plume was compared with those computed using GASP, a computational fluid dynamics code, showing reasonable agreement. The computed heating rate and pressure using PLIMP were compared with the Reaction Control System plume model (RPM) solution and the plume impingement dynamics (PIDYN) solution. RPM uses the GASP-based plume solution, whereas PIDYN uses the SCARF plume solution. Three sets of the heating rate and pressure solutions agree well. Further thermal analysis on the avionic ring of the Service Module was performed using MSC Patran/Pthermal. The obtained temperature results showed that thermal protection is necessary because of significant heating from the plume.

Wang, Xiao-Yen↗

CFD Analysis of Modular Thrusters Performance

The effective performance of modular thrusters in an aerospike configuration is difficult to determine. Standard analytical tools are applicable to conventional nozzle shapes, but are limited when applied to an aerospike nozzle (An aerospike nozzle is an altitude compensating external nozzle). Three baseline nozzle shapes are derived using standard analytical procedures. The baseline nozzle sizes are restricted to fill a volume envelope. The three shapes are an axi-symmetric round nozzle, a two dimensional planar square exit nozzle, and a super elliptic round to nearly square nozzle. The integrated (thruster/aerospike) performance of the three nozzles is determined through the use of three dimensional viscous computational fluid dynamic (CFD) calculation where complex features of the flow field can be accurately captured. The resulting installed performance is then used to evaluate the efficiency of these nozzle shapes for aerospike applications. The determination of effective performance of a thruster nozzle integrated into an aerospike nozzle require the solution of the three dimensional turbulent Navier-Stokes equations. The model used in this study consisted of two zones; one of the upstream thruster cowl surface so freestream conditions can be accurately predicted, and two, the aerospike surface beginning with the thruster outflow and extending to the end of the aerospike surface. The numerical grid consisted of over 120,000 nodes and used symmetry on the thruster centerline and edge. A two species non-reacting chemistry model was used to capture the variation of fluid properties between the hot plume base and freestream air. From the results of the three baseline nozzle aerospike calculations, the effictive performance of the nozzle was determined. The flow fields of these calculations do show some variation between the cases. Recirculation zones on the cowl surface is predicted for the two dimensional planar nozzle and a smaller one for the super elliptic nozzle. The recirculation is caused by the strong pressure gradient between the plume and freestream flows. The axi-symmetric nozzle results indicate recirculation zones on the thruster face. These recirculation zones smooth the pressure gradient between the plume and freestream flow limiting the formation of recirculation on the cowl surface. Thruster to thruster interaction is evident for the axi-symmetric and supper elliptic calculation while the two dimensional planar nozzle did not have any lateral expansion in the nozzle, so thruster to thruster interaction is limited. The integrated performance results, at the altitude choosen, show very little variation between the three thruster shapes. This result allows for nozzle shape determination based on additional considerations (thermal, structural, weight) besides performance.

Ungewitter, Ronald J.↗

Atomic oxygen testing with thermal atom systems - A critical evaluation

The use of thermal atom (kinetic energy near 0.04 eV) test methods as a materials selection and screening technique for LEO spacecraft is critically evaluated in this paper. The physics and chemistry of the thermal atom environments are shown to produce specific mass loss rates (mg/sq cm per min) and reaction efficiencies (Re) radically different from those produced in the LEO environment. A response surface study shows that specific mass loss rates change rapidly with plasma-asher parameters and seldom agree with flight data. FEP Teflon is shown to react by a different mechanism than Kapton, polyethylene, or graphite. The Re (Re = volume of material removed/oxygen atom) of Kapton, polyethylene, Mylar, Tedlar, FEP Teflon, and graphite measured in a flowing afterglow apparatus are 0.001 to 0.0001 those measured with high-energy atoms (kinetic energy 1.5 eV or greater) in beam systems or in LEO. The effect of sample temperature and atom impact energy on Re is discussed. A simple kinetic model describing the reaction of atomic oxygen with polymer surfaces is developed. Guidelines and recommendations for thermal atom testing and interpretation of test results are presented.

Koontz, Steven L.↗

Development of CFD code evaluation criteria and a procedure for assessing predictive capability and performance

Careful validation of computational fluid dynamic (CFD) codes is essential if they are to be used as engineering design tools. Validation must be carried out in a systematic manner to ensure that all code aspects as they apply to the application of interest are understood and, to the greatest extent possible, quantified. A study is being conducted in which a general code validation procedure is defined and demonstrated. A four phase validation procedure is defined in which a series of validation test cases are computed and compared with available analytical solutions and test data. The procedure is demonstrated using the REACT CFD code to compute validation cases for each of the four phases. For phase four, the application of interest, the Space Shuttle Main Engine high pressure fuel turbopump impeller flow field, is computed.

Lin, S. J.↗

Oxidation Mechanisms of Toluene and Benzene

An expanded and improved version of a previously published benzene oxidation mechanism is presented and shown to model published experimental data fairly successfully. This benzene submodel is coupled to a modified version of a toluene oxidation submodel from the recent literature. This complete mechanism is shown to successfully model published experimental toluene oxidation data for a highly mixed flow reactor and for higher temperature ignition delay times in a shock tube. A comprehensive sensitivity analysis showing the most important reactions is presented for both the benzene and toluene reacting systems. The NASA Lewis toluene mechanism's modeling capability is found to be equivalent to that of the previously published mechanism which contains a somewhat different benzene submodel.

Bittker, David A.↗

Ceramic distribution members for solid state electrolyte cells and method of producing

A solid state electrolyte cells apparatus and method of producing is disclosed. The apparatus can be used for separating oxygen from an oxygen-containing feedstock or as a fuel cell for reacting fluids. Cells can be stacked so that fluids can be introduced and removed from the apparatus through ceramic distribution members having ports designed for distributing the fluids in parallel flow to and from each cell. The distribution members can also serve as electrodes to membranes or as membrane members between electrodes. The distribution member design does not contain any horizontal internal ports which allows the member to be thin. A method of tape casting in combination with an embossing method allows intricate radial ribs and bosses to be formed on each distribution member. The bosses serve as seals for the ports and allow the distribution members to be made without any horizontal internal ports.

Clark, Douglas J.↗

Method of producing ceramic distribution members for solid state electrolyte cells

A solid state electrolyte cells apparatus and method of producing is disclosed. The apparatus can be used for separating oxygen from an oxygen-containing feedstock or as a fuel cell for reacting fluids. Cells can be stacked so that fluids can be introduced and removed from the apparatus through ceramic distribution members having ports designed for distributing the fluids in parallel flow to and from each cell. The distribution members can also serve as electrodes to membranes or as membrane members between electrodes, The distribution member design does not contain any horizontal internal ports which allows the member to be thin. A method of tape casting in combination with an embossing method allows intricate radial ribs and bosses to be formed on each distribution member. The bosses serve as seals for the ports and allow the distribution members to be made without any horizontal internal ports.

Clark, Douglas J.↗

Gas-solid reaction pathway for chlorination of rare earth and actinide metals using hydrogen and chlorine gas

Hydriding and chlorination of buttons of metallic cerium have been tested as a head-end step to separating metal impurities from actinides. The objective here was to achieve complete conversion to chlorides while maintaining high surface area via small particles and open porosity. Hydriding by reaction with 100% H 2 gas at 573 K was successful at reducing cm-scale metal buttons into small particles/fragments. This step was followed by reaction with flowing 99.5% anhydrous Cl 2 gas at temperatures ranging from 423 to 573 K for 60 min. In experiments reported here, up to 92% conversion to cerium chloride was achieved from reacting with the anhydrous chlorine gas. The final chloride product has high porosity coupled with small particle sizes (microns to 1600 μm), which makes it suitable for subsequent purification processes. Higher chlorine concentrations in the gas favor higher conversion to chloride. The optimal temperature to achieve high chlorination is 523 K. Oxygen is a minor impurity in the chlorine gas which appears to prevent complete chlorination from being achieved.

38 RADIATION CHEMISTRY, RADIOCHEMISTRY, AND NUCLEA↗

Active Aerodynamic Load Control for Wind Turbines

The goal of this project was to develop and demonstrate an advanced dielectric barrier discharge (DBD) plasma actuator technology. We set out to demonstrate the efficacy and impact of the new actuator technology as a key component of an active load control system for wind turbines. A DBD plasma actuator consists of a thin layer of dielectric material separating a pair of offset electrodes. One of those electrodes is embedded between the dielectric layer and a non-conductive substrate, while the other is exposed to air. When driven by an appropriate high voltage waveform, the device ionizes the air near the surface of the dielectric and adjacent to the exposed electrode. Collisions between the ions in the plasma and neutral air molecules result in a wall-jet – a region of induced air velocity that can be used to modify the flow around a lifting body. Use of the device near the trailing edge of a wind turbine blade, designed in such a way as to amplify the effects of the flow-modifying device, can result in large changes to the global forces experienced by the blade. Because of the fast response time of the device, it can allow the turbine to react in real time to changes in the wind associated with turbulence, wind shear, gusts, and the like, when paired with appropriate sensors and control algorithms. The goal of this project was to increase the capacity of the device to induce velocity on its surface at the levels required by large, utility-scale wind turbines. The first technical goal was to increase the induced velocity from the current industry-best of about 4 m/s to 10 m/s by modifying the electrical waveform used to drive the device and by introducing a semi-conductive surface coating to control electrostatic charge build-up on the surface. The second goal was to use the device to modify the lift on a representative airfoil in a wind tunnel, demonstrating a reduction in lift coefficient of 0.2 or better. The third goal was to produce an actuator-induced velocity of 20 m/s. Alongside our partners at the University of Texas at Dallas, we also applied modern, advanced design methods to optimize the impact of the technology on the design of wind turbines. We also had planned to install a segmented, active load control system on a test turbine to demonstrate the ability to reduce unsteady aerodynamic forces on the turbine associated with changes in the wind. The highest induced velocity achieved during this research was 11 m/s. However, practical design constraints limited the change in lift coefficient to about 0.12 for a representative airfoil in the wind tunnel at a Reynolds number of 400,000. The primary conclusion is that the plasma actuator control authority remains insufficient for practical purposes when extrapolated to Reynolds numbers over 1 million. The potential impact of the active lift control concept was evaluated through detailed simulations for three wind turbine sizes: a 3.4MW onshore turbine, a 10MW offshore turbine, and a 15MW offshore turbine. A feedback control system was designed for each turbine within two scenarios: one where the active lift control is used as a retrofit capability on the baseline design, and the other where the designers were permitted to “upscale” the turbines in order increase annual energy production. The levelized cost of energy was then evaluated for the range of turbine sizes and design configuration. It was found that the LCOE reduction associated with active lift control fell in the range of 0.7% to 7.2%, with the highest reduction associated with upscaling the 3.4MW turbine.

17 WIND ENERGY↗

Thermospheric wind effects on the global distribution of helium in the earth's upper atmosphere

The momentum and continuity equations for a minor gas are combined with the momentum equation for the major constituents to obtain the time dependent continuity equation for the minor species reflecting a wind field in the background gas. This equation is used to study the distributions of helium and argon at times of low, medium, and high solar activity for a variety of latitudinal-seasonal wind cells. For helium, the exospheric return flow at the higher thermospheric temperatures dominates the distribution to the extent that much larger latitudinal gradients can be maintained during periods of low solar activity than during periods of high activity. By comparison to the exospheric flow, the smoothing effect of horizontal diffusion is almost negligible. The latitudinal variation of helium observed by satellite mass spectrometers can be reproduced by the effect of a wind system of air rising in the summer hemisphere, flowing across the equator with speeds on the order of 100 to 200 m/sec, and descending in the winter hemisphere. Argon, being heavier than the mean mass in the lower thermosphere, reacts oppositely to helium in that it is enhanced in the summer hemisphere and depleted in the winter.

Reber, C. A.↗

Agglomeration-Free Distributor for Fluidized Beds

New gas distributor for fluidized beds prevents hot particles from reacting on it and forming hard crust. In reduction of iron ore in fluidized bed, ore particles do not sinter on distributor and perhaps clog it or otherwise interfere with gas flow. Distributor also relatively cool. In fluidized-bed production of silicon, inflowing silane does not decompose until within bed of hot silicon particles and deposits on them. Plates of spiral distributor arranged to direct incoming gas into spiral flow. Turbulence in flow reduces frequency of contact between fluidized-bed particles and distributor.

Ouyang, F.↗

Transition control of instability waves over a flexible surface in the presence of an acoustic field

Experimental results are presented which demonstrate the coupling of a laminar boundary layer flow with a typical flexible aircraft panel. It is shown that the boundary layer induces plate oscillations which, in turn, perturb the flow at the same frequencies. This feedback mechanism is an inherent property of laminar boundary layer flow passing over a flexible plate. As a result, the flexibility of the plate becomes a source of early transition. The laminar boundary layer at the leading edge of the plate reacts to small, upstream, unsteady disturbances due to a streamwise pressure gradient. The experiments demonstrate that a nominal sound pressure incident at the leading edge triggers early transition. It is shown that transition can be delayed by activating a heat source at the leading edge of the plate which results in downstream cooling.

Maestrello, L.↗

Three-dimensional rarefied flow simulations for the Aeroassist Flight Experiment vehicle

Hypersonic rarefied flow about the Aeroassist Flight Experiment vehicle has been investigated using a three-dimensional direct simulation Monte Carlo method. Calculations are performed for the transitional flows encountered during the vehicle's atmospheric entry for altitudes of 110 and 100 km with an entry velocity of 9.9 km/s. The simulations are performed using a five-species reacting gas model that account for rotational and vibrational internal energies. The solutions indicate that dissociation is important at altitudes of 110 km and below. Results are presented for surface pressures, convective heating, flowfield structure, and aerodynamic coefficient variations with altitude.

Celenligil, M. Cevdet↗