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.
This research program is concerned with the effect of low speed flow on the spreading and extinction processes of flames over solid fuels. Primary attention is given to flame propagation in concurrent flow - the more hazardous situation from the point of view of fire safety.
The interactions between flames spreading over parallel solid sheets of paper are being studied in normal gravity and in microgravity. This geometry is of practical importance since in most heterogeneous combustion systems, the condensed phase is non-continuous and spatially distributed. This spatial distribution can strongly affect burning and/or spread rate. This is due to radiant and diffusive interactions between the surface and the flames above the surfaces. Tests were conducted over a variety of pressures and separation distances to expose the influence of the parallel sheets on oxidizer transport and on radiative feedback.
Our work was primarily theoretical and numerical. We investigated the simplified modeling of heat losses in diffusion flames, then we 'ramped up' the level of complexity in each successive study until the final chapter discussed the general problem of soot/flame interaction. With regard to the specific objective of studying radiative extinction, we conclude that in the steady case a self-extinguishing zero-g flame is unlikely to occur. The soot volume fractions are too small. On the other hand, our work does provide rational means for assessing the mixture of chemical energy release and radiative heat release. It also provides clues for suitable 'tailoring' this balance. Thus heat fluxes to surrounding surfaces can be substantially increased by exploiting and modifying its sooting capability.
Tunable solid-state lasers with direct UV pumping, based on d-f transitions of rare earth ions incorporated in wide band-gap dielectric crystals, are reliable sources of laser radiation that are suitable for excitation of combustion-related free radicals. We have employed such a laser for analytical flame characterization utilizing Laser-Induced Fluorescence (LIF) techniques. LIF spectra of alkane-air flames (used for studying combustion processes under normal and microgravity conditions) excited in the region of the A-X (0,0) OH-absorption band have been recorded and found to be both temperature-sensitive and positionally-sensitive. In addition, also clearly noticeable was the sensitivity of the spectra to the specific wavelength used for data registration. The LiCAF:Ce laser shows good prospects for being able to cover the spectral region between 280 and 340 nm and therefore be used excitation of combustion-intermediates such as the hydroxyl OH, methoxy CH30 and methylthio CH3S radicals.
Papers included address the following topics: (1) Turbulent premixed flame propagation in microgravity; (2) The effect of gravity on turbulent premixed flame propagation - a preliminary cold flow study; and (3) Characteristics of a subgrid model for turbulent premixed combustion.
Studies in the mathematical modeling of the high-speed turbulent combustion has received renewal attention in the recent years. The review of fundamentals, approaches and extensive bibliography was presented by Bray, Libbi and Williams. In order to obtain accurate predictions for turbulent combustible flows, the effects of turbulent fluctuations on the chemical source terms should be taken into account. The averaging of chemical source terms requires to utilize probability density function (PDF) model. There are two main approaches which are dominant in high-speed combustion modeling now. In the first approach, PDF form is assumed based on intuitia of modelliers (see, for example, Spiegler et.al.; Girimaji; Baurle et.al.). The second way is much more elaborate and it is based on the solution of evolution equation for PDF. This approach was proposed by S.Pope for incompressible flames. Recently, it was modified for modeling of compressible flames in studies of Farschi; Hsu; Hsu, Raji, Norris; Eifer, Kollman. But its realization in CFD is extremely expensive in computations due to large multidimensionality of PDF evolution equation (Baurle, Hsu, Hassan).
An acoustic flame-out detection system that renders a large high pressure combustor safe in the event of a flame-out and possible explosive reignition. A dynamic pressure transducer is placed in the fuel and detects the stabilizing fuel pressure oscillations, caused by the combustion process. An electric circuit converts the signal from the combustion vortices, and transmitted to the fuel flow to a series of pulses. A missing pulse detector counts the pulses and continuously resets itself. If three consecutive pulses are missing, the circuit closes the fuel valve. With fuel denied the combustor is shut down or restarted under controlled conditions.
Non-isothermal studies of cool flames and low temperature oxidation reactions in unstirred closed vessels are complicated by the perturbing effects of natural convection at earth gravity. Buoyant convection due to self-heating during the course of slow reaction produces spatio-temporal variations in the thermal and thus specie concentration fields due to the Arrhenius temperature dependence of the reaction rates. Such complexities have never been quantitatively modeled and were the primary impetus for the development of CSTR's (continuously stirred tank reactors) 30 years ago. While CSTR's have been widely adopted since they offer the advantage of spatial uniformity in temperature and concentration, all gradients are necessarily destroyed along with any structure that may otherwise develop. Microgravity offers a unique environment where buoyant convection can be effectively minimized and the need for stirring eliminated. Moreover, eliminating buoyancy and the need for stirring eliminates complications associated with the induced hydrodynamic field whose influence on heat transport and hot spot formation, hence explosion limits, is not fully realized. The objective of this research is to quantitatively determine and understand the fundamental mechanisms that control the onset and evolution of low temperature reactions and cool flames in both static and flow reactors. Microgravity experiments will be conducted to obtain benchmark data on the structure (spatio-temporal temperature, concentration, flow fields), the dynamics of the chemical fronts, and the ignition diagrams (pressure vs. temperature). Ground-based experiments will be conducted to ascertain the role of buoyancy. Numerical simulations including detailed kinetics will be conducted and compared to experiment.
We developed a new experimental technique called: Combined laser sheet particle tracking (LSPT) and laser holographic interferometry (HI), which is capable of measuring the transient behavior of three dimensional structures of temperature and flow both in liquid and gas phases. We applied this technique to a pulsating flame spread over n-butanol. We found a twin vortex flow both on the liquid surface and deep in the liquid a few mm below the surface and a twin vortex flow in the gas phase. The first twin vortex flow at the liquid surface was observed previously by NASA Lewis researchers, while the last two observations are new. These observations revealed that the convective flow structure ahead of the flame leading edge is three dimensional in nature and the pulsating spread is controlled by the convective flow of both liquid and gas.
A principal characteristic of combustion in microgravity is the absence of buoyancy driven flows. In some cases, such as for spherically symmetrical droplet burning, the absence of buoyancy is desirable for matching analytical treatments with experiments. In other cases, however, it can be more valuable to arbitrarily control the flame's convective environment independent of the environmental gravitational condition. To accomplish this, we propose the use of ion generated winds driven by electric fields to control local convection of flames. Such control can produce reduced buoyancy (effectively zero buoyancy) conditions in the laboratory in 1-g facilitating a wide range of laser diagnostics that can probe the system without special packaging required for drop tower or flight tests. In addition, the electric field generated ionic winds allow varying gravitational convection equivalents even if the test occurs in reduced gravity environments.
The interaction of a laminar flame with its self-generated boundary layer in a rectangular channel was numerically simulated using the two-dimensional, reacting, Navier-Stokes equations. A two species chemistry model was implemented which simulates the stoichiometric reaction of acetylene and air. Calculations were performed to investigate the effects of altering the boundary condition of the wall temperature, the Lewis number, the dynamic viscosity, and the ignition method. The purpose of this study was to examine the fundamental physics of the formation of the boundary layer and the interaction of the flame as it propagates into the boundary layer that its own motion has created.
A detailed numerical study was conducted on the dynamics and thermal response of inert, spherical particles in strained, laminar, premixed hydrogen/air flames. The modeling included the solution of the steady conservation equations for both the gas and particle phases along and around the stagnation streamline of an opposed-jet configuration, and the use of detailed descriptions of chemical kinetics and molecular transport, For the gas phase, the equations of mass, momentum, energy, and species are considered, while for the particle phase, the model is based on conservation equations of the particle momentum balance in the axial and radial direction, the particle number density, and the particle thermal energy equation. The particle momentum equation includes the forces as induced by drag, thermophoresis, and gravity. The particle thermal energy equation includes the convective/conductive heat exchange between the two phases, as well as radiation emission and absorption by the particle. A one-point continuation method is also included in the code that allows for the description of turning points, typical of ignition and extinction behavior. As expected, results showed that the particle velocity can be substantially different than the gas phase velocity, especially in the presence of large temperature gradients and large strain rates. Large particles were also found to cross the gas stagnation plane, stagnate, and eventually reverse as a result of the opposing gas phase velocity. It was also shown that the particle number density varies substantially throughout the flowfield, as a result of the straining of the flow and the thermal expansion. Finally, for increased values of the particle number density, substantial flame cooling to extinction states and modification of the gas phase fluid mechanics were observed. As also expected, the effect of gravity was shown to be important for low convective velocities and heavy particles. Under such conditions, simulations indicate that the magnitude and direction of the gravitational force can substantially affect the profiles of the particle velocity, number density, mass flux, and temperature.
We consider the SHS (self-propagating high-temperature synthesis) process for synthesizing materials. In this process a powder mixture of reactants is cold pressed into a sample, which is ignited at one end. A high temperature combustion wave then propagates through the sample converting reactants to the desired product material. In this process, melting of some or all the components is often observed. Therefore, we study combustion waves propagating through a high caloricity inorganic powder mixture whose combustion temperature exceeds the melting temperatures of many components. The solid matrix is thus destroyed by the propagating combustion wave due to melting ahead of the reaction zone, and a liquid bath is formed which contains gaseous bubbles. The waves propagate in the presence of a gravitational field. Due to the effect of gravity, there is relative motion between the rising bubbles and the descending bath, which affects the composition of the medium, its thermophysical properties, the 'liquid flame' structure, and the propagation velocity. To enhance our understanding of phenomena associated with the interaction of the relative motion with the propagating combustion wave we formulate and analyze a relatively simple mathematical model of liquid flames in a gravitational field. We describe the wave structure and combustion characteristics including the combustion velocity. We compare our results to existing experimental observations and suggest new experiments to be performed. We consider the effects of gravity and, in particular, examine both microgravity and large gravity conditions.
Self-propagating high-temperature synthesis (SHS) is a process in which combustion waves, e.g., "solid flames", which are considered here, are employed to synthesize desired materials. Like many other systems, SHS is a pattern forming system. The problem of describing experimentally observed patterns and of predicting new, as yet unobserved, patterns continues to attract the attention of scientists and mathematicians due to the fundamental significance of the phenomena in combustion in particular, and in nonlinear science in general. Here, we analyze the dynamics of solid flame propagation in a 2D region by considering the region to be composed of parallel, identical layers aligned along the direction of propagation and having thermal contact. Each layer is then described by wave propagation in 1D, with the transverse Laplacian replaced by a term describing heat exchange between neighboring layers. This configuration is the simplest model of a 2D system because it accounts, in a simple way, for the principal feature of the problem, i.e., heat exchange between neighbors in the transverse direction. For simplicity, we describe the situation for two layers. Because the layers are identical, uniformly propagating waves in each layer must be identical, independent of the heat exchange rate alpha. When the Zeldovich number Z exceeds a critical value Z(sub c), which depends on alpha, uniformly propagating waves become unstable. The stability diagram for the two coupled layers reproduces that for the full 2D problem after appropriate identification of parameters in the two problems. Depending on parameter values, we determine three different steady-state dynamical behaviors (though additional behaviors are also expected to occur). The three behaviors are: (i) waves in each layer which pulsate in phase as they propagate, so that together they form a single pulsating propagating wave; (ii) the waves in each layer are no longer identical, and antiphase pulsations occur, with the two waves alternately advancing and receding as they propagate. This wave is the analog of the spinning wave on the surface of a circular cylinder; (iii) finally, there is a region of bistability between the in phase and antiphase waves. with each having its own domain of attraction, so that which of the two behaviors occur depends on the condition of initiation of the wave. The results of our computations indicate a qualitative similarity in the behavior of combustion waves in the layer model and in the full 2D model. Specifically, due to the similarity between the small alpha wave behavior in the layer model and the large diameter behavior in the model of waves on the surface of a cylinder, we are able to predict the behavior of the mean velocity for the waves on the cylinder, where computations of the full problem can be rather difficult. We also compute solutions for three or more layers. The results of our computations prompt us to predict that, while planar uniformly propagating waves are unstable, the wave will be quasiplanar, i.e., the resulting spinning waves have very low amplitude hot spots, and travel with the velocity close to that of the uniformly propagating wave. Such waves may be difficult to distinguish from uniformly propagating waves in experiments. We also find that for both the layer model and full 2D problem, steady-state time-dependent waves, e.g., pulsating and spinning wave, have a conserved quantity H which characterizes the excess energy in the wave, just as in the case of uniformly propagating waves. The quantity H, which is generated by dissipation, does not vary in time and is proportional to the diffusivity and caloricity of the system, and inversely proportional to the mean wave velocity.
We report spatially resolved laser-saturated fluorescence measurements of NO concentration in a pre-heated, lean-direct injection (LDI) spray flame at atmospheric pressure. The spray is produced by a hollow-cone, pressure-atomized nozzle supplied with liquid heptane. NO is excited via the Q2(26.5) transition of the gamma(0,0) band. Detection is performed in a 2-nm region centered on the gamma(0,1) band. Because of the relatively close spectral spacing between the excitation (226 nm) and detection wavelengths (236 nm), the gamma(0,1) band of NO cannot be isolated from the spectral wings of the Mie scattering signal produced by the spray. To account for the resulting superposition of the fluorescence and scattering signals, a background subtraction method has been developed that utilizes a nearby non-resonant wavelength. Excitation scans have been performed to locate the optimum off-line wavelength. Detection scans have been performed at problematic locations in the flame to determine possible fluorescence interferences from UHCs and PAHs at both the on-line and off-line excitation wavelengths. Quantitative radial NO profiles are presented and analyzed so as to better understand the operation of lean-direct injectors for gas turbine combustors.
The basic objective of this work is to assess the influence of gravity on "the compositional and the spatial structures" of transitional and turbulent diffusion flames via large eddy simulation (LES), and direct numerical simulation (DNS). The DNS is conducted for appraisal of the various closures employed in LES, and to study the effect of buoyancy on the small scale flow features. The LES is based on our "filtered mass density function"' (FMDF) model. The novelty of the methodology is that it allows for reliable simulations with inclusion of "realistic physics." It also allows for detailed analysis of the unsteady large scale flow evolution and compositional flame structure which is not usually possible via Reynolds averaged simulations.
A detailed numerical study was conducted on the dynamics and thermal response of inert spherical particles in strained, laminar, premixed hydrogen/air flames. The modeling included the solution of the steady conservation equations for both the gas and particle phases along and around the stagnation streamline of an opposed-jet configuration, and the use of detailed descriptions of chemical kinetics and molecular transport. For the gas phase, the equations of mass, momentum, energy, and species are considered, while for the particle phase, the model is based on conservation equations of the particle momentum balance in the axial and radial direction, the particle number density, and the particle thermal energy equation. The particle momentum equation includes the forces as induced by drag, thermophoresis, and gravity. The particle thermal energy equation includes the convective/conductive heat exchange between the two phases, as well as radiation emission and absorption by the particle. A one-point continuation method is also included in the code that allows for the description of turning points, typical of ignition and extinction behavior. As expected, results showed that the particle velocity can be substantially different than the gas phase velocity, especially in the presence of large temperature gradients and large strain rates. Large particles were also found to cross the gas stagnation plane, stagnate, and eventually reverse as a result of the opposing gas phase velocity. It was also shown that the particle number density varies substantially throughout the flowfield, as a result of the straining of the flow and the thermal expansion. Finally, for increased values of the particle number density, substantial flame cooling to extinction states and modification of the gas phase fluid mechanics were observed. As also expected, the effect of gravity was shown to be important for low convective velocities and heavy particles. Under such conditions, simulations indicate that the magnitude and direction of the gravitational force can substantially affect the profiles of the particle velocity, number density, mass flux, and temperature.
Microgravity experiments on flame spread over thermally thick fuels were conducted using foam fuels to obtain low density and thermal conductivity, and thus large spread rate (Sf) compared to dense fuels such as PMMA. This scheme enabled meaningful results to lie obtained even in 2.2 second drop tower experiments. It was found that, in contrast conventional understanding; steady spread can occur over thick fuels in quiescent microgravity environments, especially when a radiatively active diluent gas such as CO2 is employed. This is proposed to be due to radiative transfer from the flame to the fuel surface. Additionally, the transition from thermally thick to thermally thin behavior with decreasing bed thickness is demonstrated.