9.3. Diffusional Kinetics
In the case of heterogeneous surface burning of a particle, consideration must be given to the question of whether diffusion rates or surface kinetic reaction rates are controlling the overall burning rate of the material. In many cases, it cannot be assumed that the surface oxidation kinetic rate is fast compared to the rate of diffusion of oxygen to the surface. The surface temperature determines the rate of oxidation, and this temperature is not always known a priori. Thus, in surface combustion the assumption that chemical kinetic rates are much faster than diffusion rates cannot be made.
Consider, for example, a carbon surface burning in a concentration of oxygen in the free stream specified by
ρmo∞. The burning is at a steady mass rate. Then the concentration of oxygen at the surface is some value
mos. If the surface oxidation rate follows first-order kinetics, as Frank-Kamenetskii
[21] assumed,
Gox=Gfi=ksρmos
(9.9)
where G is the flux in grams per second per square centimeter; ks the heterogeneous specific reaction rate constant for surface oxidation in units reflecting a volume to surface area ratio, that is, centimeters per second; and i the mass stoichiometric index. The problem is that mos is unknown. But one knows that the consumption rate of oxygen must be equal to the rate of diffusion of oxygen to the surface. Thus, if hDρ is designated as the overall convective mass transfer coefficient (conductance), one can write
Gox=ksρmos=hDρ(mo∞−mos)
(9.10)
What is sought is the mass burning rate in terms of mo∞. It follows that
hDmos=hDmo∞−ksmos
(9.11)
ksmos+hDmos=hDmo∞
(9.12)
mos=(hDks+hD)mo∞
(9.13)
Gox=(ρkshDks+hD)mo∞=ρKmo∞
(9.14)
K≡kshDks+hD
(9.15)
1K=ks+hDkshD=1hD+1ks
(9.16)
When the kinetic rates are large compared to the diffusion rates,
K =
hD; when the diffusion rates are large compared to the kinetic rates,
K =
ks. When
k <<
hD,
mos ≅
mo∞ from Eqn
(9.13); thus
Gox=ksρmo∞
(9.17)
When
ks >>
hD, Eqn
(9.13) gives
mos=(hDks)mo∞
(9.18)
But since
ks >>
hD, it follows from Eqn
(9.18) that
mos<<mo∞
(9.19)
This result permits one to write Eqn
(9.10) as
Gox=hDρ(mo∞−mos)≅hDρmo∞
(9.20)
Consider the case of rapid kinetics,
ks >>
hD, further. In terms of Eqn
(9.14), or examining Eqn
(9.20) in light of
K,
K=hD
(9.21)
Of course, Eqn
(9.20) also gives one the mass burning rate of the fuel
Gfi=Gox=hDρmo∞
(9.22)
where
hD is the convective mass transfer coefficient for an unspecified geometry. For a given geometry,
hD would contain the appropriate boundary layer thickness, or it would have to be determined by independent measurements giving correlations that permit
hD to be found from other parameters of the system. More interestingly, Eqn
(9.22) should be compared to Eqn
(6.179) in
Chapter 6, which can be written as
Gox≅λcpδmo∞HLv=λcpρρδmo∞HLv=αδHLvρmo∞=DδHLvρmo∞
(9.23)
Thus one notes that the development of Eqn
(9.22) is for a small
B number, in which case
hD=DδHLv
(9.24)
where the symbols are as defined in
Chapter 6. (
H/
Lv) is a simplified form of the
B number. Nevertheless, the approach leading to Eqn
(9.22) gives simple physical insight into surface oxidation phenomena where the kinetic and diffusion rates are competitive.