Journal of Combustion

Volume 2013 (2013), Article ID 752068, 8 pages

http://dx.doi.org/10.1155/2013/752068

## Effects of Electric Field on the SHS Flame Propagation of the Si-C System, Examined by the Use of the Heterogeneous Theory

Institute of Aeronautical Technology, Japan Aerospace Exploration Agency, 7-44-1 Jindaiji-Higashi, Chofu, Tokyo, Japan

Received 16 April 2013; Revised 6 August 2013; Accepted 15 August 2013

Academic Editor: Peter F. Nelson

Copyright © 2013 Atsushi Makino. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

#### Abstract

Relevant to the self-propagating high-temperature synthesis (SHS) process, an analytical study has been conducted to investigate the effects of electric field on the combustion behavior because the electric field is indispensable for systems with weak exothermic reactions to sustain flame propagation. In the present study, use has been made of the heterogeneous theory which can satisfactorily account for the premixed mode of the bulk flame propagation supported by the nonpremixed mode of particle consumption. It has been confirmed that, even for the SHS flame propagation under electric field, being well recognized to be facilitated, there exists a limit of flammability, due to heat loss, as is the case for the usual SHS flame propagation. Since the heat loss is closely related to the representative sizes of particles and compacted specimen, this identification provides useful insight into manipulating the SHS flame propagation under electric field, by presenting appropriate combinations of those sizes. A fair degree of agreement has been demonstrated through conducting an experimental comparison, as far as the trend and the approximate magnitude are concerned, suggesting that an essential feature has been captured by the present study.

#### 1. Introduction

Self-propagating high-temperature synthesis (SHS) process, by virtue of a strong exothermic reaction that passes through a compacted mixture of particles, as a flame, has attracted special interests as a rapid and economical way in synthesizing inorganic and/or intermetallic compounds [1–5]. More than hundreds of kinds of materials, including borides, carbides, and silicides, are reported to be synthesized [1–5] by applying this process, while some of them are quite difficult to synthesize in conventional ways. Materials synthesized are now being considered for use as electronic materials, materials resistant to heat, corrosion, and/or wear, and so forth [1–5]. Production of functionally graded materials (FGMs) [6], composed of different components with continuous profiles, is also intended. Note also that near-net-shape fabrication can be anticipated because this is a kind of powder metallurgy.

The SHS process, applicable to various combinations of solid-solid, solid-liquid, and/or solid-gas systems, as pioneered by the group of Merzhanov and Borovinskaya [7], has also been recognized to present diverse phenomena of the flame propagation, such as pulsating, spinning, and/or repeated combustions, as well as steady propagation, as reviewed [1–5, 8, 9]. Among various systems, synthesis of Ti-C system has extensively been examined, because of its simplicity in chemical reaction, as well as its nature of refractory and hardness. It has been confirmed that not the gaseous component but the molten Ti, spreading over carbon particles, plays an important role in the process. Effects of dominant parameters, such as the preheating, the mixture ratio of reactants, the degree of dilution, and the particle size, on the flame propagation speed and/or the range of flammability have been examined. In particular, preheating for rapid synthesis has been examined from the viewpoint of process innovation, to have higher production rates and/or to facilitate melting during the SHS process. Although rapid synthesis can easily be accomplished, with increasing initial temperature, a supplementary heating can sometimes cause an ignition when the compacted mixture is heated over its self- (or spontaneous) ignition temperature [1, 2, 4, 5, 10], as are the cases for usual combustion.

Even for the spontaneous ignition, its utilization has been considered seriously [11–13], for those systems with weak exothermic reactions, because it is not common for those systems that sufficient preheating can be anticipated for sustaining SHS flame propagation. An external heating, by virtue of the Joule heating with applying electric fields [14–20], has also been taken into account, in order to make the materials synthesis by combustion completed.

In this study, silicon carbide (SiC) has been chosen for the research subject since it has attracted special attention in the field of materials science as the light-weight materials in high-temperature structural applications, relevant to the advanced aviation gas turbines. Focus is put on the SHS flame propagation that proceeds in the compacted mixture of Si-C system, from the viewpoint of combustion engineering, which has been rare, although examinations from the viewpoint of materials science have been conducted extensively [14–20]. Since self-propagation of the SHS flame cannot be anticipated without preheating and/or applying electric field, especially for the Si-C system, here, it is intended to elucidate effects of electric field on the burning velocity and/or the limit of flammability, by the use of the heterogeneous theory [9] for the SHS process, which can satisfactorily account for the premixed mode of the bulk flame propagation supported by the nonpremixed mode of particle consumption. Experimental comparisons are also conducted, which are found to be satisfactory, as far as the trend and approximate magnitude are concerned.

#### 2. Formulation

##### 2.1. Model Definition

The problem of interest, as shown in Figure 1, is the one-dimensional, planner, heterogeneous flame propagation under electric field, being applied perpendicular to the direction of the propagation, in accordance with the experiment [19]. It is assumed in the heterogeneous theory [9] that the flame front propagates in the doubly infinite domain of a condensed medium, originally consisting of a mixture of particles of nonmetal (or higher melting-point metal) *N*, lower melting-point metal , and an inert that can also be the product *P* of the reaction between and according to
where is the stoichiometric coefficient of the reaction. For simplicity, it is assumed that the nonmetal particle-size distribution is monodisperse, with an initial radius and number density *n*_{0}. It is also assumed that there is no reaction between and until the mixture has been heated to the melting point of the metal , at which all the metal particles melt instantaneously. The reaction between the solid and molten is then assumed to take place at the particle surface at a finite rate and proceed until all of the particles or species are consumed. The enthalpy of phase change is neglected because it is much smaller than the heat of combustion. The lateral heat loss to the ambience, occurring throughout the entire flame structure, is also taken into account for the flame propagation considered.

##### 2.2. Governing Equations and Boundary Conditions

By the use of the heterogeneous theory [9] for the SHS flame propagation, which can fairly formulate the situation involving a premixed mode of propagation for the bulk flame, supported by the nonpremixed reaction of the dispersed nonmetal particles in the melt, the governing equations are expressed as follows.

Continuity:

energy conservation:

*M*-conservation:

*N*-consumption:
where is the velocity along the coordinate, the mass burning rate (per unit area) of the flame, the temperature, the specific heat, the heat of reaction per unit mass of the mass fraction of fluid the fluid density, the mass of solid per unit spatial volume, the total density the thermal conductivity, the mass fraction of in the fluid, the stoichiometric mass ratio the molecular weight, and the subscript 0 designates the unburned state.

The surface regression rate , defined by , is given [9] as where , , , and are, respectively, the reduced Damköhler number, the Damköhler number, the frequency factor, and the activation temperature of the surface reaction. Note here that the nondimensional surface regression rate in (6) is presented by the use of the transfer number in terms of the natural logarithmic term, just like that for the droplet combustion, and that the transfer number here is expressed as that similar to the well-known expression for the combustion rate of solid for the first-order kinetics, suggesting that roles of diffusion and chemical kinetics are both included in the present formulation. In addition, is the Arrhenius mass diffusivity with an activation temperature , being anticipated to increase markedly over a relatively thin, high temperature combustion zone, under the condition that the diffusion-limited situation prevails in the chemical kinetics at high temperatures.

A set of these governing equations is basically the same as that in the previous works [9], except for the term that appears in the energy conservation equation, in order to take account of an effect of electric field. While the heat-loss rate , as expressed in the previous works, is given as
where is the emissivity and is the Stefan-Boltzmann constant, the energy supply rate *P *due to the electric field is expressed as
where is the electric field (V/m) and is the electric conductivity , with the Joule heating being taken into account. In accordance with experimental observations [17, 19, 20] that the Joule heating due to external energy supply has mainly been conducted in the combustion zone, here, it is a priori assumed that it can be made by virtue of the conductible, nonmetal particles in the melt within the combustion zone, so that the mass ratio appears in (8). It is also assumed that no energy supplies before arrival of the SHS flame can be expected, because of the extremely high resistance of contact among particles, prior to the appearance of melt. After completion of the materials synthesis, we do not expect any energy supplies, because of the high electric resistance of SiC, which is too high to allow for the electric current.

By introducing nondimensional variables and parameters as the governing equations are as follows: In the above, and are, respectively, the mass burning rate and the mass burning rate eigenvalue in the adiabatic condition. The nondimensional temperature is defined as . Note here that , called the heat-loss parameter and defined as a ratio of the heat-loss and heat-release rates, depends not only on the physicochemical parameters but also on the particle diameter and representative size of a sample specimen [9]. As for the parameter introduced here, we shall call it a heat-input parameter, hereafter.

Equations (9) through (11), together with the mass burning rate , are to be solved, subject to the boundary conditions Note that the cold boundary difficulty is eliminated by specifying the reaction to be initiated at the melting point of . Thus, the present problem is reduced to be an eigenvalue problem for obtaining as an eigenvalue, under a nonadiabatic condition in which there exists not only heat loss but also external energy supply.

##### 2.3. Burning Velocity under Adiabatic Condition

When there is neither heat loss nor external energy supply, the burning velocity is obtained as [9] once the specific value of the eigenvalue is determined by conducting a numerical calculation. It is seen that the burning velocity is inversely proportional to the particle diameter . It may be informative to note that the product of and has even been called the SHS rate-constant [9]. The mass ratio of the lower melting-point metal , before and after the combustion, is expressed as Here, is the mixture ratio, defined as the initial molar ratio of to divided by the corresponding stoichiometric molar ratio; is the degree of dilution, defined as the initial mass fraction of diluent.

##### 2.4. Parameters

Values of physicochemical parameters employed are those of the Si-C system: , , , , , and . Although some of the chemical kinetic data are reported in [17], it is not the chemical kinetics but the mass diffusivity, at high temperatures, as explained in Section 2.2. Therefore, use has been made of the mass diffusivity, [21], so that the representative Lewis number is set to be . The total density before and after the combustion is assumed to be equal. Note that the thermophysical properties used here implicitly account for effects of compact density, gases in void spaces, gas evolution, and so forth. Other thermophysical properties, indispensable for examining the limit of flammability, related to the heat loss, are for the thermometric conductivity [22] and for the emissivity [23]. In addition, use has been made of the electric conductivity of carbon given by being obtained by the use of a curve fitting of data in the literature [24], as a function of the representative temperature , defined as an arithmetic mean of the melting point of Si and the maximum flame temperature.

#### 3. Results

##### 3.1. Effects of Heat-Input Parameter in the Adiabatic Condition

First, effects of heat-input parameter on the mass burning rate are examined in the adiabatic condition. Here, mass burning rate under electric field, , has been normalized by the use of that without electric field, . Figure 2 shows the normalized mass burning rate for the stoichiometric Si-C system, as a function of the heat-input parameter , with the initial temperature taken as a parameter. With increasing heat-input parameter , the normalized mass burning rate increases first gradually and then rapidly, up to the limiting value determined in (28), to be mentioned later. It is seen that external energy supply is effective in enhancing mass burning rate and that its contribution to becomes large with increasing initial temperature .

##### 3.2. Effects of Heat-Input Parameter under Heat Loss Condition

Effects of heat loss on the mass burning rate are then examined. In this examination, mass burning rate is normalized by the adiabatic mass burning rate. Figure 3 shows the normalized mass burning rate , as a function of the heat-loss parameter , with the heat-input parameter taken as a parameter. When there is no electric field (), the well-known, characteristic extinction turning-point behavior is exhibited, with the upper branch being the stable solution and the turning point designating the state of extinction [9]. In contrast, the trend becomes quite different from this, with increasing heat-input parameter . When , although there exists the turning-point behavior, the solution ceases to exist at certain and , because of other restrictions. When or more, there appears the limit of flammability, designated by an open circle, before reaching the turning point. It is seen that, with increasing , the range of for the flame propagation contracts, while the range of expands.

Relevant to the limit of flammability under electric field, an attempt has been made to examine a range for the existence of SHS flame, which yields a relation in (29), to be mentioned later. In Figure 3, those limits are also shown by dotted curves, below which the steady propagation ceases to exist. It should be noted that, in determining the position of the SHS flame, theoretical consideration in the next section is indispensable.

#### 4. Theoretical Consideration

##### 4.1. Flame Position and Temperature

In the limit of large Zeldovich number, melting, diffusion, and consumption/convection are confined to a thin layer in the neighborhood of the SHS flame. Outside this layer, the diffusion and consumption/convection terms are exponentially small. If we consider the situation that the flame locates at and that the electric field is supplied in the region , temperature profile outside the flame is expressed as
where
Relations among constants , , , and *C* are determined by the use of the continuity in temperature and heat flux at the position , the continuity in temperature at the flame position (), and the jump condition at the flame:
Then, we have

If we further note that the temperature at is the melting point , we can determine the flame position as by the use of (21). The flame temperature in the adiabatic condition (), which is indispensable for nondimensionalization, is obtained from (24) as We see that the adiabatic flame temperature is increased due to the application of the electric field.

##### 4.2. Range of Electric Field

Even though the electric field is effective in sustaining SHS flame propagation, there exists an indispensable requirement for the flame position to be positive (), which further requires that the logarithmic term in (27) should be positive. From this requirement, we can determine the range for the valid heat-input parameter as In the same manner, when there exists heat loss, we have from (25) the following relation between the heat-loss parameter and the heat-input parameter :

#### 5. Experimental Comparison

Although an introduction of the heat input parameter has turned out to be useful in examining the limit of flammability, effects of electric field on the mass burning rate are to be examined directly because use has been made of in defining , as shown in (12). Figure 4 shows the SHS rate-constant as a function of the electric field , with the initial temperature taken as a parameter. The SHS rate-constant expressed in (15) is used, instead of the normalized mass burning rate used in Figures 2 and 3. It is seen that the SHS rate-constant in the adiabatic condition, shown by solid curves, increases gradually with an increase in the electric field , because of the external energy supply. It is also seen that, at a certain electric field, the SHS rate-constant becomes high, with increasing initial temperature , because of the preheating effect.

A dashed curve in Figure 4 shows the limit of flammability. As the representative size of the cross-sectional area in the compacted specimen, indispensable in determining heat-loss parameter in (11), use has been made of 9.175 mm for , determined by use of the relation for a rectangular cross-section as [9] with mm and mm, in accordance with the experiment by Feng and Munir [19]. As for the particle size of carbon, use has been made of , reported. Although the limit of flammability, shown in Figure 4, is only that for K, it has been confirmed in a preparatory study that it is nearly independent of the initial temperature , within the range examined.

Data points in Figure 4 are experimental in the literature [19]; a symbol designates a result for relative density 0.64 and for relative density 0.52. It is seen that there is no remarkable effect of the relative density and that a fair degree of agreement exists between experimental and theoretical results, as far as the trend and approximate magnitude are concerned. It is also seen that the SHS flame propagation under low electric fields, say, less than about 1000 V/m, proceeds close to the limit of flammability at K. To the contrary, at high electric fields, the SHS flame propagation proceeds under the condition, free from heat loss, because of the external energy supply.

#### 6. Concluding Remarks

In the present study, flame propagation in the SHS process of the Si-C system has been examined by the use of the heterogeneous theory, with taking account of an effect of electric field. It has been shown that the burning velocity is enhanced by applying an electric field, due to an increase in the external energy supply, as reported in the literature. Furthermore, it has been confirmed that there exists the limit of flammability, even for the SHS flame propagation under electric field, as is the case for usual SHS flame propagation, although its dependence on the heat-loss rate is quite different from that without electric field. A fair degree of agreement between experimental and theoretical results suggests that the present study has captured the essential feature of the SHS flame propagation under electric field, providing useful insight into manipulating the SHS flame propagation, by the use of representative sizes of particles and compacted specimen, which has not been recognized in the previous studies.

#### Nomenclature

: | Constant for temperature profile |

: | Length of a rectangular sample specimen [m] |

: | Constant for temperature profile |

: | Width of a rectangular sample specimen [m] |

: | Constant for temperature profile |

: | Specific heat [J/(kg·K)] |

: | Diffusivity [m^{2}/s] |

: | Electric field [V/m] |

: | Stoichiometric mass ratio |

: | Heat input parameter |

: | Heat loss [J/(m^{3}·s)] |

: | Lewis number |

: | Representative size of the rectangular sample specimen [m] |

: | Mass burning rate [kg/(m^{2}·s)] |

: | Energy supply due to electric field [J/(m^{3}·s)] |

: | Heat of combustion per unit mass of [J/kg] |

: | Particle radius [m] |

: | Radius of sample specimen [m] |

: | Temperature [K] |

: | Activation temperature of mass diffusivity [K] |

: | Burning velocity [m/s] |

: | Molecular weight |

: | Mass fraction |

: | Mass ratio of species . |

*Greek*

: | Emissivity |

: | Nondimensional temperature |

: | Degree of dilution |

: | Eigenvalue of the mass burning rate |

: | Thermal conductivity [W/(m·K)] or exponent for temperature profile |

: | Mixture ratio |

: | Stoichiometric coefficient |

: | Normalized mass ratio of |

: | Nondimensional length |

: | Electric conductivity [] |

: | Stefan-Boltzmann constant [] |

: | Density [kg/m^{3}] |

: | Surface regression rate [m^{2}/s] |

: | Heat loss parameter. |

*Subscript*

: | Adiabatic condition |

: | Downstream |

: | Flame |

f: | Fluid or melt |

: | Electric field |

: | Hot zone |

: | Lower melting-point metal |

: | Melting point |

: | Nonmetal or higher melting-point metal |

: | Combustion product or diluent |

: | Solid or surface |

: | Total |

: | Upstream |

0: | Initial or unburned state. |

*Superscript*

Nondimensional or stoichiometrically mass weighted. |

#### References

- W. L. Frankhouser, K. W. Brendley, M. C. Kieszek, and S. T. Sullivan,
*Gasless Combustion Synthesis of Refractory Compounds*, Noyes, Park Ridge, NJ, USA, 1985. - Z. A. Munir and U. Anselmi-Tamburini, “Self-propagating exothermic reactions: the synthesis of high-temperature materials by combustion,”
*Materials Science Reports*, vol. 3, no. 7-8, pp. 277–365, 1989. View at Google Scholar · View at Scopus - A. G. Merzhanov,
*Self-Propagating High-Temperature Synthesis: Twenty Years of Search and Findings*, Combustion and Plasma Synthesis of High-Temperature Materials, VCH Publishers, New York, NY, USA, 1990, Edited by Z. A. Munir and J. B. Holt. - A. Varma and J.-P. Lebrat, “Combustion synthesis of advanced materials,”
*Chemical Engineering Science*, vol. 47, no. 9–11, pp. 2179–2194, 1992. View at Google Scholar · View at Scopus - A. Varma, A. S. Rogachev, A. S. Mukasyan, and S. Hwang, “Combustion synthesis of advanced materials: principles and applications,”
*Advances in Chemical Engineering*, vol. 24, no. C, pp. 79–226, 1998. View at Publisher · View at Google Scholar · View at Scopus - W. G. J. Bunk,
*Gradient Materials For Structural and Functional Applications*, Advanced Materials '93, III/B, Elsevier, Amsterdam, The Netherlands, 1994. - A. G. Merzhanov and I. P. Borovinskaya, “A new class of combustion processes,”
*Combustion Science and Technology*, vol. 10, no. 5-6, pp. 195–201, 1975. View at Google Scholar · View at Scopus - S. B. Margolis, “The transition to nonsteady deflagration in gasless combustion,”
*Progress in Energy and Combustion Science*, vol. 17, no. 2, pp. 135–162, 1991. View at Google Scholar · View at Scopus - A. Makino, “Fundamental aspects of the heterogeneous flame in the self-propagating high-temperature synthesis (SHS) process,”
*Progress in Energy and Combustion Science*, vol. 27, no. 1, pp. 1–74, 2001. View at Publisher · View at Google Scholar · View at Scopus - A. G. Merzhanov, “Nonisothermal methods in chemical kinetics,”
*Combustion, Explosion, and Shock Waves*, vol. 9, no. 1, pp. 3–28, 1975. View at Publisher · View at Google Scholar · View at Scopus - A. S. Rogachev, A. S. Mukasyan, and A. Varma, “Volume combustion modes in heterogeneous reaction systems,”
*Journal of Materials Synthesis and Processing*, vol. 10, no. 1, pp. 31–36, 2002. View at Publisher · View at Google Scholar · View at Scopus - J. P. Lebrat, A. Varma, and A. E. Miller, “Combustion synthesis of Ni
_{3}Al and Ni_{3}Al-matrix composites,”*Metallurgical Transactions A*, vol. 23, no. 1, pp. 69–76, 1992. View at Publisher · View at Google Scholar · View at Scopus - L. Thiers, A. S. Mukasyan, and A. Varma, “Thermal explosion in Ni-Al system: influence of reaction medium microstructure,”
*Combustion and Flame*, vol. 131, no. 1-2, pp. 198–209, 2002. View at Publisher · View at Google Scholar · View at Scopus - V. A. Knyazik, A. G. Merzhanov, V. B. Solomonov, and A. S. Shteinberg, “Macrokinetics of high-temperature titanium interaction with carbon under electrothermal explosion conditions,”
*Combustion, Explosion, and Shock Waves*, vol. 21, no. 3, pp. 333–337, 1985. View at Publisher · View at Google Scholar · View at Scopus - O. Yamada, Y. Miyamota, and M. Koizumi, “Self-propagating high-temperature synthesis of the SiC,”
*Journal of Materials Research*, vol. 1, no. 2, pp. 275–279, 1986. View at Publisher · View at Google Scholar - V. A. Knyazik, A. G. Merzhanov, and A. S. Shteinberg, “About combustion mechanism of the Ti-C system,”
*Doklady Akademii Nauk SSSR*, vol. 301, no. 3, pp. 899–902, 1989. View at Google Scholar - V. I. Gorovenko, V. A. Knyazik, and A. S. Shteinberg, “High-temperature interaction between silicon and carbon,”
*Ceramics International*, vol. 19, no. 2, pp. 129–132, 1993. View at Google Scholar · View at Scopus - A. Feng and Z. A. Munir, “The effect of an electric field on self-sustaining combustion synthesis: part I. modeling studies,”
*Metallurgical and Materials Transactions B*, vol. 26, no. 3, pp. 581–586, 1995. View at Publisher · View at Google Scholar · View at Scopus - A. Feng and Z. A. Munir, “The effect of an electric field on self-sustaining combustion synthesis: part II. field-assisted synthesis of
*β*-SiC,”*Metallurgical and Materials Transactions B*, vol. 26, no. 3, pp. 587–593, 1995. View at Publisher · View at Google Scholar · View at Scopus - A. Feng and Z. A. Munir, “Relationship between field direction and wave propagation in activated combustion synthesis,”
*Journal of the American Ceramic Society*, vol. 79, no. 8, pp. 2049–2058, 1996. View at Google Scholar · View at Scopus - Thermophysical Properties Handbook,
*Japan Society of Thermophysical Properties*, Yokendo, Tokyo, Japan, 1990. - Thermophysical Properties Handbook,
*Japan Society of Thermophysical Properties*, Yokendo, Tokyo, Japan, 1990. - S. Hanzawa, “Refractories of furnaces to reduce environmental impact,”
*Annual Report of Ceramics Research Laboratory*, vol. 9, pp. 33–42, 2009 (Japanese). View at Google Scholar *Thermophysical Properties of High Temperature Solid Materials*, vol. 1, MacMillan, New York, NY, USA, 1965, Edited by: Y. S. Touloukian.