Search NASA⌕ Search

Engineering topics

Matkowsky, B. J.

Publications and source records attributed to Matkowsky, B. J..

6 records
All sourcesNASA NTRSDOE OSTI
NASA NTRS2000

Diffusion Driven Combustion Waves in Porous Media

Filtration of gas containing oxidizer, to the reaction zone in a porous medium, due, e.g., to a buoyancy force or to an external pressure gradient, leads to the propagation of Filtration combustion (FC) waves. The exothermic reaction occurs between the fuel component of the solid matrix and the oxidizer. In this paper, we analyze the ability of a reaction wave to propagate in a porous medium without the aid of filtration. We find that one possible mechanism of propagation is that the wave is driven by diffusion of oxidizer from the environment. The solution of the combustion problem describing diffusion driven waves is similar to the solution of the Stefan problem describing the propagation of phase transition waves, in that the temperature on the interface between the burned and unburned regions is constant, the combustion wave is described by a similarity solution which is a function of the similarity variable x/square root of(t) and the wave velocity decays as 1/square root of(t). The difference between the two problems is that in the combustion problem the temperature is not prescribed, but rather, is determined as part of the solution. We will show that the length of samples in which such self-sustained combustion waves can occur, must exceed a critical value which strongly depends on the combustion temperature T(sub b). Smaller values of T(sub b) require longer sample lengths for diffusion driven combustion waves to exist. Because of their relatively small velocity, diffusion driven waves are considered to be relevant for the case of low heat losses, which occur for large diameter samples or in microgravity conditions, Another possible mechanism of porous medium combustion describes waves which propagate by consuming the oxidizer initially stored in the pores of the sample. This occurs for abnormally high pressure and gas density. In this case, uniformly propagating planar waves, which are kinetically controlled, can propagate, Diffusion of oxidizer decreases the wave velocity. In addition to the reaction and diffusion layers, the uniformly propagating wave structure includes a layer with a pressure gradient, where the gas motion is induced by the production or consumption of the gas in the reaction as well as by thermal expansion of the gas. The width of this zone determines the scale of the combustion wave in the porous medium.

Aldushin, A. P.↗
NASA NTRS2000

Gas-Phase Influence on Quasisteady "Liquid Flames" in Gravitational Fields

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.

Shkadinsky, K. G.↗
NASA NTRS2000

Slowly and Rapidly Propagating "Liquid Flames" in Gravitational Fields

We consider the combustion, in a gravitational field, of a heterogeneous powder mixture compressed into a solid sample, in which the high temperature ahead of the reaction zone destroys the solid, due, e.g. to melting of some of components of the mixture. Thus, a suspension is formed, consisting of a liquid bath containing solid or liquid particles. Processes such as heat and mass transfer as well as chemical reactions in the suspension determine the structure of the combustion wave and its propagation velocity. Under the influence of gravitational forces there is the possibility of relative motion of the liquid and solid. Previous theoretical analyses considered the rate of beat transfer between the solid and liquid phases to be sufficiently large that their two distinct temperatures rapidly equilibrated to a single temperature. In addition to this case, we also consider the case when the rate of heat transfer is not so large and the model involves the separate temperatures of the solid and liquid phases. We find that multiplicity of traveling wave structures is possible. In particular, in addition to a low velocity structure, which is essentially the same as that obtained from the one temperature description, we find a high velocity structure, which does not exist in the one temperature description, but rather depends on the fact that the solid and fluid temperatures differ from each other. Both structures can exist for the same parameter values in a given range. We describe the dependence of the combustion characteristics of the two structures on gravitational forces and other factors. In particular, we compare the characteristics in gravity and microgravity environments.

Shkadinsky, K. G.↗
NASA NTRS2000

Dynamics in Layer Models of Solid Flame Propagation

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.

Aldushin, A. P.↗
NASA NTRS1989

An adaptive pseudo-spectral method for reaction diffusion problems

The spectral interpolation error was considered for both the Chebyshev pseudo-spectral and Galerkin approximations. A family of functionals I sub r (u), with the property that the maximum norm of the error is bounded by I sub r (u)/J sub r, where r is an integer and J is the degree of the polynomial approximation, was developed. These functionals are used in the adaptive procedure whereby the problem is dynamically transformed to minimize I sub r (u). The number of collocation points is then chosen to maintain a prescribed error bound. The method is illustrated by various examples from combustion problems in one and two dimensions.

Bayliss, A.↗
NASA NTRS1987

An adaptive pseudo-spectral method for reaction diffusion problems

The spectral interpolation error was considered for both the Chebyshev pseudo-spectral and Galerkin approximations. A family of functionals I sub r (u), with the property that the maximum norm of the error is bounded by I sub r (u)/J sub r, where r is an integer and J is the degree of the polynomial approximation, was developed. These functionals are used in the adaptive procedure whereby the problem is dynamically transformed to minimize I sub r (u). The number of collocation points is then chosen to maintain a prescribed error bound. The method is illustrated by various examples from combustion problems in one and two dimensions.

Bayliss, A.↗