SEARCH · Engineering Papers

Results for “DIFFUSION”

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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 415 records · Page 23

All sources NASA NTRS DOE OSTI

NASA NTRS1988

Venus' superrotation, mixing length theory and eddy diffusion - A parametric study

The concept of the Hadley mechanism is adopted to describe the axisymmetric circulation of the Venus atmosphere. It is shown that, for the atmosphere of a slowly rotating planet such as Venus, a form of the nonliner 'closure' (self-consistent solution) of the fluid dynamics system which constrains the magnitude of the eddy diffusion coefficients can be postulated. A nonlinear one-layer spectral model of the zonally symmetric circulation was then used to establish the relationship between the heat source, the meridional circulation, and the eddy diffusion coefficients, yielding large zonal velocities. Computer experiments indicated that proportional changes in the heat source and eddy diffusion coefficients do not significantly change the zonal velocities. It was also found that, for large eddy diffusion coefficients, the meridional velocity is virtually constant; below a threshold in the diffusion rate, the meridional velocity decreases; and, for large eddy diffusion and small heating rates, the zonal velocities decrease with decreasing planetary rotation rates.

Mayr, H. G.↗

NASA NTRS1989

Encapsulated diffusion of sulphur into InP

This talk outlines a simple process for the fabrication of n(+)-p solar cells in indium phosphide. Large area cells (greater than 0.25 sq cm) have been made by this process, with a photovoltaic conversion efficiency of 15.21 percent under AM0 conditions of illumination. An ideality factor of 1.1 and a saturation current density of 8 x 10 to the minus 15th power A/sq cm have been observed for these cells. The technique for cell fabrication involves the diffusion of sulfur into InP by an open tube process, and gives highly reproducible results from run to run. A vacuum-deposited layer of gallium sulphide (Ga2S3) was used as the source for sulfur diffusion, with a chemically vapor deposited SiO2 cap layer to prevent decomposition of the InP surface during heat treatment. Diffusions were carried out in a flowing nitrogen ambient at 585 to 708 C, and characterized by their surface carrier concentration and the diffusion constant. The diffusion profile for sulfur in InP is estimated to be of the complementary error function type. The activation energy of the diffusion was estimated to be 1.94 eV. The technique described here is ideally suited for the fabrication of shallow n(+)-p junctions in InP, and has been used for space-borne solar cells.

Parat, K. K.↗

NASA NTRS1992

Estimation of the viscosities of silicate liquids at high pressure from measurements of oxygen diffusivity

The dynamics and evolution of a magma ocean depend significantly on the structure and physical properties of silicate liquids at high pressure. Viscosity is particularly important because of its effect on the dynamics of convection and crystal setting. Traditionally, the viscosities of silicate liquids have been measured at high pressure by falling sphere viscometry but, due to some limitations of this technique, few measurements have been made at pressures exceeding 2.5 GPa. An alternative approach, which has previously been proven up to 2 GPa, is to use high pressure measurements of oxygen self-diffusivity (D) to estimate viscosity (eta) from the Eyring relationship eta = kT/D lambda (where k is the Boltzmann constant, T is absolute temperature, and lambda is the diffusive jump distance). We have performed preliminary experiments on Na2Si4O9 liquid in order to explore the feasibility of this method up to pressures in excess of 10 GPa. Diffusion couples were prepared from glass starting materials with one half of the sample enriched in O-18 and the other half containing the natural abundance (0.2 wt percent). Diffusion experiments were performed at 1600 to 1800 C at pressures in the range of 2.5 to 10 GPa for times up to 6 minutes using a 1200 ton multianvil apparatus. Oxygen diffusivities were derived from the resulting O-18 concentration profiles, which were analyzed using an ion microprobe. At 1800 C, oxygen diffusivities increase continuously from 1 x 10(exp -10) m(exp -2) s(exp -1) at 2.5 Gpa to 5 x 10(exp -10) m(exp -2) s(exp -1) at 10 GPa. From the Eyring relationship, the viscosity of Na2SiO9 liquid is predicted to decrease by approx. 0.7 log units over this pressure range at 1800 C. These trends, which can be related to changes in the structure of the liquid at high pressure, are in agreement with results of molecular dynamics calculations. These results suggest that it should be possible to estimate the viscosities of silicate liquids up to at least 15 GPa and 2200 C by this method.

Rubie, D. C.↗

NASA NTRS1992

Cu-Zn binary phase diagram and diffusion couples

The objectives of this paper are to learn: (1) what information a binary phase diagram can yield; (2) how to construct and heat treat a simple diffusion couple; (3) how to prepare a metallographic sample; (4) how to operate a metallograph; (5) how to correlate phases found in the diffusion couple with phases predicted by the phase diagram; (6) how diffusion couples held at various temperatures could be used to construct a phase diagram; (7) the relation between the thickness of an intermetallic phase layer and the diffusion time; and (8) the effect of one species of atoms diffusing faster than another species in a diffusion couple.

Mccoy, Robert A.↗

NASA NTRS1992

Diffusion mechanisms in chemical vapor-deposited iridium coated on chemical vapor-deposited rhenium

Radiation-cooled rocket thruster chambers have been developed which use CVD Re coated with CVD Ir on the interior surface that is exposed to hot combustion gases. The Ir serves as an oxidation barrier which protects the structural integrity-maintaining Re at elevated temperatures. The diffusion kinetics of CVD materials at elevated temperatures is presently studied with a view to the prediction and extension of these thrusters' performance limits. Line scans for Ir and Re were fit on the basis of a diffusion model, in order to extract relevant diffusion constants; the fastest diffusion process is grain-boundary diffusion, where Re diffuses down grain boundaries in the Ir overlayer.

Hamilton, J. C.↗

NASA NTRS1992

Self-diffusion of magnesium in spinel and in equilibrium melts - Constraints on flash heating of silicates

An isotopic tracer is used to measure Mg self-diffusion in spinel and coexisting melt at bulk chemical equilibrium. The diffusion coefficients were calculated from the measured isotope profiles using a model that includes the complementary diffusion of Mg-24, Mg-25, and Mg-26 in both phases with the constraint that the Mg content of each phase is constant. The activation energy and preexponential factor for Mg self-diffusion in spinel are, respectively, 384 +/- 7 kJ and 74.6 +/- 1.1 sq cm/s. These data indicate Mg diffusion in spinel is much slower than previous estimates. The activation energy for Mg self-diffusion in coexisting melt is 343 +/- 25 kJ and the preexponential factor is 7791.9 +/- 1.3 sq cm/s. These results are used to evaluate cooling rates of plagioclase-olivine inclusions (POIs) in the Allende meteorite. Given a maximum melting temperature for POIs of about 1500 C, these results show that a 1-micron radius spinel would equilibrate isotopically with a melt within about 60 min.

Sheng, Y. J.↗

NASA NTRS1993

Effects of Buoyancy on Laminar, Transitional, and Turbulent Gas Jet Diffusion Flames

Gas jet diffusion flames have been a subject of research for many years. However, a better understanding of the physical and chemical phenomena occurring in these flames is still needed, and, while the effects of gravity on the burning process have been observed, the basic mechanisms responsible for these changes have yet to be determined. The fundamental mechanisms that control the combustion process are in general coupled and quite complicated. These include mixing, radiation, kinetics, soot formation and disposition, inertia, diffusion, and viscous effects. In order to understand the mechanisms controlling a fire, laboratory-scale laminar and turbulent gas-jet diffusion flames have been extensively studied, which have provided important information in relation to the physico-chemical processes occurring in flames. However, turbulent flames are not fully understood and their understanding requires more fundamental studies of laminar diffusion flames in which the interplay of transport phenomena and chemical kinetics is more tractable. But even this basic, relatively simple flame is not completely characterized in relation to soot formation, radiation, diffusion, and kinetics. Therefore, gaining an understanding of laminar flames is essential to the understanding of turbulent flames, and particularly fires, in which the same basic phenomena occur. In order to improve and verify the theoretical models essential to the interpretation of data, the complexity and degree of coupling of the controlling mechanisms must be reduced. If gravity is isolated, the complication of buoyancy-induced convection would be removed from the problem. In addition, buoyant convection in normal gravity masks the effects of other controlling parameters on the flame. Therefore, the combination of normal-gravity and microgravity data would provide the information, both theoretical and experimental, to improve our understanding of diffusion flames in general, and the effects of gravity on the burning process in particular.

Bahadori, M. Yousef↗

NASA NTRS1953

Performance Characteristics of Plane-Wall Two-Dimensional Diffusers

Experiments have been made at Stanford University to determine the performance characteristics of plane-wall, two-dimensional diffusers which were so proportioned as to insure reasonable approximation of two-dimensional flow. All of the diffusers had identical entrance cross sections and discharged directly into a large plenum chamber; the test program included wide variations of divergence angle and length. During all tests a dynamic pressure of 60 pounds per square foOt was maintained at the diffuser entrance and the boundary layer there was thin and fully turbulent. The most interesting flow characteristics observed were the occasional appearance of steady, unseparated, asymmetric flow - which was correlated with the boundary-layer coalescence - and the rapid deterioration of flow steadiness - which occurred as soon as the divergence angle for maximum static pressure recovery was exceeded. Pressure efficiency was found to be controlled almost exclusively by divergence angle, whereas static pressure recovery was markedly influenced by area ratio (or length) as well as divergence angle. Volumetric efficiency. diminished as area ratio increased, and at a greater rate with small lengths than with large ones. Large values of the static-pressure-recovery coefficient were attained only with long diffusers of large area ratio; under these conditions pressure efficiency was high and. volumetric efficiency low. Auxiliary tests with asymmetric diffusers demonstrated that longitudinal pressure gradient, rather than wall divergence angle, controlled flow separation. Others showed that the addition of even a short exit duct of uniform section augmented pressure recovery. Finally, it was found that the installation of a thin, central, longitudinal partition suppressed flow separation in short diffusers and thereby improved pressure recovery

Reid, Elliott G↗

NASA NTRS1994

Technical report series on global modeling and data assimilation. Volume 2: Direct solution of the implicit formulation of fourth order horizontal diffusion for gridpoint models on the sphere

High order horizontal diffusion of the form K Delta(exp 2m) is widely used in spectral models as a means of preventing energy accumulation at the shortest resolved scales. In the spectral context, an implicit formation of such diffusion is trivial to implement. The present note describes an efficient method of implementing implicit high order diffusion in global finite difference models. The method expresses the high order diffusion equation as a sequence of equations involving Delta(exp 2). The solution is obtained by combining fast Fourier transforms in longitude with a finite difference solver for the second order ordinary differential equation in latitude. The implicit diffusion routine is suitable for use in any finite difference global model that uses a regular latitude/longitude grid. The absence of a restriction on the timestep makes it particularly suitable for use in semi-Lagrangian models. The scale selectivity of the high order diffusion gives it an advantage over the uncentering method that has been used to control computational noise in two-time-level semi-Lagrangian models.

Max J. Suarez↗

NASA NTRS1994

Design of a variable area diffuser for a 15-inch Mach 6 open-jet tunnel

The Langley 15-inch Mach 6 High Temperature Tunnel was recently converted from a Mach 10 Hypersonic Flow Apparatus. This conversion was effected to improve the capability of testing in Mach 6 air at relatively high reservoir temperatures not previously possible at Langley. Elevated temperatures allow the matching of the Mach numbers, Reynolds numbers, and ratio of wall-to-adiabatic-wall temperatures (TW/Taw) between this and the Langley 20-inch Mach 6 CF4 Tunnel. This ratio is also matched for Langley's 31-inch Mach 10 Tunnel and is an important parameter useful in the simulation of slender bodies such as National Aerospace Plane (NASP) configurations currently being studied. Having established the nozzle's operating characteristics, the decision was made to install another test section to provide model injection capability. This test section is an open-jet type, with an injection system capable of injecting a model from retracted position to nozzle centerline between 0.5 and 2 seconds. Preliminary calibrations with the new test section resulted in Tunnel blockage. This blockage phenomenon was eliminated when the conical center body in the diffuser was replaced. The issue then, is to provide a new and more efficient variable area diffuser configuration with the capability to withstand testing of larger models without sending the Tunnel into an unstart condition. Use of the 1-dimensional steady flow equation with due regard to friction and heat transfer was employed to estimate the required area ratios (exit area / throat area) in a variable area diffuser. Correlations between diffuser exit Mach number and area ratios, relative to the stagnation pressure ratios and diffuser inlet Mach number were derived. From these correlations, one can set upper and lower operating pressures and temperatures for a given diffuser throat area. In addition, they will provide appropriate input conditions for the full 3-dimensional computational fluid dynamics (CFD) code for further simulation studies.

Loney, Norman W.↗

NASA NTRS1994

Radial diffusion models of energetic electrons and Jupiter's synchrotron radiation. 2: Time variability

We used a radial diffusion code for energetic electrons in Jupiter's magnetosphere to investigate variations in Jupiter's radio emission due to changes in the electron phase space density at L shells between 6 and 50, and due to changes in the radial diffusion parameters. We suggest that the observed variations in Jupiter's radio emission are likely caused by changes in the electron phase space density at some boundary L(sub 1) is greater than 6, if the primary mode of transport of energetic electrons is radial diffusion driven by fluctuating electric and/or magnetic fields induced by upper atmospheric turbulence. We noticed an excellent empirical correlation, both in phase and relative amplitude, between changes in the solar wind ram pressure and Jupiter's synchrotron radiation if the electron phase space density at the boundary L(sub 1) (L(sub 1) is approximately equal to 20-50) varies linearly with the square root of the solar wind ram pressure, f is approximately (N(sub s)nu(exp 2 sub s))(exp 1/2). The calculations were carried out with a diffusion coefficient D(sub LL) = D(sub n)L(exp n) with n = 3. The diffusion coefficient which best fit the observed variations in Jupiter's synchrotron radiation D(sub 3) = 1.3 +/- 0.2 x 10(exp -9)/s is approximately 0.041/yr, which corresponds to a lagtime of approximately 2 years. We further show that the observed short term (days-weeks) variations in Jupiter's radio emission cannot be explained adequately when radial diffusion is taken into account.

De Pater, Imke↗

NASA NTRS1995

Unsteady planar diffusion flames: Ignition, travel, burnout

In microgravity, a thin planar diffusion flame is created and thenceforth travels so that the flame is situated at all times at an interface at which the hydrogen and oxygen meet in stoichiometric proportion. If the initial amount of hydrogen is deficient relative to the initial amount of oxygen, then the planar flame will travel further and further into the half volume initially containing hydrogen, until the hydrogen is (virtually) fully depleted. Of course, when the amount of residual hydrogen becomes small, the diffusion flame is neither vigorous nor thin; in practice, the flame is extinguished before the hydrogen is fully depleted, owing to the finite rate of the actual chemical-kinetic mechanism. The rate of travel of the hydrogen-air diffusion flame is much slower than the rate of laminar flame propagation through a hydrogen-air mixture. This slow travel facilitates diagnostic detection of the flame position as a function of time, but the slow travel also means that the time to burnout (extinction) probably far exceeds the testing time (typically, a few seconds) available in earth-sited facilities for microgravity-environment experiments. We undertake an analysis to predict (1) the position and temperature of the diffusion flame as a function of time, (2) the time at which extinction of the diffusion flame occurs, and (3) the thickness of quench layers formed on side walls (i.e., on lateral boundaries, with normal vectors parallel to the diffusion-flame plane), and whether, prior to extinction, water vapor formed by burning will condense on these cold walls.

Fendell, F.↗

NASA NTRS1996

Analysis of Subgrid Eddy Diffusivity in Turbulence Using Direct Numerical Simulation Databases

Direct numerical simulation data in isotropic turbulence at Taylor-scale Reynolds number 90 are used to compute subgrid transfer and spectral eddy diffusivities for passive scalars at Schmidt numbers 1 and 1/8. The results are consistent with spectral transfer processes dominated by a local forward cascade in wavenumber space which is a consequence of moderately nonlocal interactions in the scalar fields. Data decompositions show that resolvable-subgrid interactions make (in general) substantially stronger contributions than fully-subgrid interactions to the eddy diffusivity, which has a non-negligible inverse cascade, or backscatter, component. The eddy diffusivity has a strong cusp at wavenumbers close to the cutoff value separating resolvable and subgrid ranges. At low wavenumbers the eddy diffusivity is generally small, but there the details depend on Schmidt number. The behavior of the subgrid turbulent Schmidt number (the ratio of effective eddy viscosity to eddy diffusivity) is observed to depend strongly on Schmidt number. In interpreting the results at a given spectral cutoff it is important to note that the spectrum of less diffusive scalars (of higher Schmidt number) has more high-wavenumber content.

Yeung, P. K.↗

NASA NTRS1996

Theoretical and Numerical Investigation of Radiative Extinction of Diffusion Flames

The influence of soot radiation on diffusion flames was investigated using both analytical and numerical techniques. Soot generated in diffusion flames dominate the flame radiation over gaseous combustion products and can significantly lower the temperature of the flame. In low gravity situations there can be significant accumulation of soot and combustion products in the vicinity of the primary reaction zone owing to the absence of any convective buoyant flow. Such situations may result in substantial suppression of chemical activities in a flame, and the possibility of a radiative extinction may also be anticipated. The purpose of this work was to not only investigate the possibility of radiative extinction of a diffusion flame but also to qualitatively and quantitatively analyze the influence of soot radiation on a diffusion flame. In this study, first a hypothetical radiative loss profile of the form of a sech(sup 2) was assumed to influence a pure diffusion flame. It was observed that the reaction zone can, under certain circumstances, move through the radiative loss zone and locate itself on the fuel side of the loss zone contrary to our initial postulate. On increasing the intensity and/or width of the loss zone it was possible to extinguish the flame, and extinction plots were generated. In the presence of a convective flow, however, the movement of the temperature and reaction rate peaks indicated that the flame behavior is more complicated compared to a pure diffusional flame. A comprehensive model of soot formation, oxidation and radiation was used in a more involved analysis. The soot model of Syed, Stewart and Moss was used for soot nucleation and growth and the model of Nagle and Strickland-Constable was used for soot oxidation. The soot radiation was considered in the optically thin limit. An analysis of the flame structure revealed that the radiative loss term is countered both by the reaction term and the diffusion term. The essential balance for the soot volume fraction was found to be between the processes of soot convection and soot growth. Such a balance yielded to analytical treatment and the soot volume fraction could be expressed in the form of an integral. The integral was evaluated using two approximate methods and the results agreed very well with the numerical solutions for all cases examined.

Ray, Anjan↗

NASA NTRS1998

NASA/DOD Aerospace Knowledge Diffusion Research Project. Paper 67: Maximizing the Results of Federally-Funded Research and Development Through Knowledge Management: A Strategic Imperative for Improving US Competitiveness

Federally-funded research and development (R&D) represents a significant annual investment (approximately $79 billion in fiscal year 1996) on the part of U.S. taxpayers. Based on the results of a 10-year study of knowledge diffusion in U.S. aerospace industry, the authors take the position that U.S. competitiveness will be enhanced if knowledge management strategies, employed within a capability-enhancing U.S. technology policy framework, are applied to diffusing the results of federally-funded R&D. In making their case, the authors stress the importance of knowledge as the source of competitive advantage in today's global economy. Next, they offer a practice-based definition of knowledge management and discuss three current approaches to knowledge management implementation-mechanistic, "the learning organization," and systemic. The authors then examine three weaknesses in existing U.S. public policy and policy implementation-the dominance of knowledge creation, the need for diffusion-oriented technology policy, and the prevalence of a dissemination model- that affect diffusion of the results of federally-funded R&D. To address these shortcomings, they propose the development of a knowledge management framework for diffusing the results of federally-funded R&D. The article closes with a discussion of some issues and challenges associated with implementing a knowledge management framework for diffusing the results of federally-funded R&D.

Pinelli, Thomas E.↗

NASA NTRS1999

Radiant Extinction of Gaseous Diffusion Flames

The absence of buoyancy-induced flows in microgravity (mu-g) and the resulting increase in the reactant residence time significantly alters the fundamentals of many combustion processes. Substantial differences between normal gravity (ng) and mu-g flames have been reported in experiments on candle flames, flame spread over solids, droplet combustion, and others. These differences are more basic than just in the visible flame shape. Longer residence times and higher concentration of combustion products in the flame zone create a thermochemical environment that changes the flame chemistry and the heat and mass transfer processes. Processes such as flame radiation, that are often ignored in ng, become very important and sometimes even controlling. Furthermore, microgravity conditions considerably enhance flame radiation by: (1) the build-up of combustion products in the high-temperature reaction zone which increases the gas radiation; and (2) longer residence times make conditions appropriate for substantial amounts of soot to form which is also responsible for radiative heat loss. Thus, it is anticipated that radiative heat loss may eventually extinguish the "weak" (low burning rate per unit flame area) mu-g diffusion flame. Yet, space shuttle experiments on candle flames show that in an infinite ambient atmosphere, the hemispherical candle flame in mu-g will burn indefinitely. This may be because of the coupling between the fuel production rate and the flame via the heat-feedback mechanism for candle flames, flames over solids and fuel droplet flames. Thus, to focus only on the gas-phase phenomena leading to radiative extinction, aerodynamically stabilized gaseous diffusion flames are examined. This enables independent control of the fuel flow rate to help identify conditions under which radiative extinction occurs. Also, spherical geometry is chosen for the mu-g experiments and modeling because: (1) It reduces the complexity by making the problem one-dimensional; (2) The spherical diffusion flame completely encloses the soot which is formed on the fuel rich side of the reaction zone. This increases the importance of flame radiation because now both soot and gaseous combustion products co-exist inside the high temperature spherical diffusion flame; (3) For small fuel injection velocities, as is usually the case for a pyrolyzing solid, the diffusion flame in mu-g around the solid naturally develops spherical symmetry. Thus, spherical diffusion flames are of interest to fires in mu-g and identifying conditions that lead to radiation-induced extinction is important for spacecraft fire safety.

Berhan, Sean↗

NASA NTRS1999

Flame Shapes of Luminous NonBuoyant Laminar Coflowing Jet Diffusion Flames

Laminar diffusion flames are of interest as model flame systems that are more tractable for analysis and experiments than practical turbulent diffusion flames. Certainly understanding laminar flames must precede understanding more complex turbulent flames while man'y laminar diffusion flame properties are directly relevant to turbulent diffusion flames using laminar flamelet concepts. Laminar diffusion flame shapes have been of interest since the classical study of Burke and Schumann because they involve a simple nonintrusive measurement that is convenient for evaluating flame structure predictions. Motivated by these observations, the shapes of laminar flames were considered during the present investigation. The present study was limited to nonbuoyant flames because most practical flames are not buoyant. Effects of buoyancy were minimized by observing flames having large flow velocities at small pressures. Present methods were based on the study of the shapes of nonbu,3yant round laminar jet diffusion flames of Lin et al. where it was found that a simple analysis due to Spalding yielded good predictions of the flame shapes reported by Urban et al. and Sunderland et al.

Lin, K.-C.↗

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.↗