MODELING CHAR OXIDATION AS A FUNCTION OF PRESSURE ...
MODELING CHAR OXIDATION AS A FUNCTION OF PRESSURE ... MODELING CHAR OXIDATION AS A FUNCTION OF PRESSURE ...
Boundary Layer Diffusion The molar flux of oxygen, N O2, in the bulk phase can be related to the surface mole fraction (Bird et al., 1960): N O2 − x s (N O2 + N CO + N CO 2 ) = k xm (x ∞ − x s ) (2.6) where N denotes the molar flux of a substance, x is the oxygen mole fraction in the bulk stream, x s is the oxygen mole fraction at the external surface of the particle, and k xm is the mass transfer coefficient and can be obtained from the Sherwood number correlation for spheres in a convective flow (Bird et al., 1960; Field et al., 1967; Mulcahy and Smith, 1969) k xm d p C f D ABf = Sh = 2.0 + 0.60Re 1/2 Sc 1/3 where Sh is the Sherwood number, C f is the total gas concentration at the film 10 (2.7) temperature, D ABf is the molecular diffusivity at the film temperature, Re is the Reynolds number, and Sc is the Schmidt number. Note that in Eq. 2.6, the positive flux direction is designated as the direction from the bulk phase to the particle. As a result, N CO and N CO2 take negative values. Using the stoichiometric relations in Eq. (2.3), N CO and N CO2 in Eq. (2.4) can be expressed in terms of N O2. Eq. 2.6 can be re-written as where = NO2 − x s NO 2 = C f DABf Sh (x∞ − x s ) (2.8) d p − 1 + 1 The oxygen molar flux can be converted to carbon consumption rate q diff (gC/cm 2 /sec) by multiplying the molecular weight of carbon and the reciprocal of the stoichiometric coefficient of oxygen (2.9)
q diff (1 − x s ) = M C It is convenient to define a new parameter k D as o o c f DABf Sh (x∞ − xs ) (2.10) d p k D = M c C f DABf Sh 1 P = MC d p Eq. 2.10 can thus be written as o D ABf Sh d p 11 1 RT f (2.11) qdiff (1 − Ps P ) = kD (P ∞ − Ps ) (2.12) The net mass diffusion rate (N O2 + N CO + N CO2) is often neglected (equivalent to assuming equil-molar counter diffusion), and the above equation is further simplified to q diff = k D (P ∞ − P s ) (2.13) Despite the widespread use of this simplified equation (Smith, 1982; Essenhigh, 1988), the more accurate form (Eq. 2.12) is recommended. At high temperatures, surface reaction is so fast that the surface oxygen partial pressure approaches zero, and the overall reaction rate approaches the maximum value allowed by boundary layer diffusion: q max = q diff Ps = 0 = k D P ∞ (2.14) In this case the overall reaction rate is solely controlled by boundary layer diffusion. This situation is also called Zone III combustion. In the char combustion literature, the factor is often used to determine the importance of boundary layer diffusion effects. The factor is defined as the observed reaction rate (g/sec/cm 2 ) over the maximum reaction rate allowed by boundary layer diffusion (g/sec/cm 2 ):
- Page 1 and 2: MODELING CHAR OXIDATION AS A FUNCTI
- Page 3 and 4: BRIGHAM YOUNG UNIVERSITY As chair o
- Page 5 and 6: CBK model uses: 1) an intrinsic Lan
- Page 7 and 8: Table of Contents List of Figures..
- Page 9: Appendices.........................
- Page 12 and 13: Figure A.2. Mass releases of the Ko
- Page 14 and 15: Table 7.6. Parameters Used in Model
- Page 16 and 17: Ed activation energy of desorption,
- Page 18 and 19: vc the volume of combustible materi
- Page 21 and 22: Background 1. Introduction The rate
- Page 23: the CBK model developed at Brown Un
- Page 26 and 27: Zone III rate ∝ C og E obs → 0
- Page 28 and 29: coal-general kinetic rate constants
- Page 32 and 33: = q obs q max The factor can be use
- Page 34 and 35: where k 1 and K are two kinetic par
- Page 36 and 37: particle can therefore be convenien
- Page 38 and 39: This is the first time that the gen
- Page 40 and 41: Data of Mathias Mathias (1996) perf
- Page 42 and 43: urn with shrinking diameters, and t
- Page 45 and 46: 3. Objectives and Approach The obje
- Page 47 and 48: Introduction 4. Analytical Solution
- Page 49 and 50: Task and Methodology Task One of th
- Page 51 and 52: 2 [ (i +1) − (i − 1)] i b = −
- Page 53 and 54: Table 4.1. The Relative Error * (%)
- Page 55 and 56: The resulting observations regardin
- Page 57 and 58: correction. The values of functions
- Page 59 and 60: Table 4.6. The Relative Error* (%)
- Page 61 and 62: Table 4.8. The Relative Error* (%)
- Page 63 and 64: general asymptotic solution. An arc
- Page 65 and 66: 5. Theoretical Developments The int
- Page 67 and 68: order of a reaction is usually dete
- Page 69 and 70: nobs = 1 (KCs ) 2 2 1 [KCs − ln(1
- Page 71 and 72: The observed reaction order in Zone
- Page 73 and 74: Bulk Diffusion vs. Knudsen Diffusio
- Page 75 and 76: where D K is in cm 2 /sec, r p is t
- Page 77 and 78: where T p is in K, P is in atm. The
- Page 79 and 80: Both of these assumptions are argua
Boundary Layer Diffusion<br />
The molar flux of oxygen, N O2, in the bulk phase can be related to the surface mole<br />
fraction (Bird et al., 1960):<br />
N O2 − x s (N O2 + N CO + N CO 2 ) = k xm (x ∞ − x s ) (2.6)<br />
where N denotes the molar flux of a substance, x is the oxygen mole fraction in the bulk<br />
stream, x s is the oxygen mole fraction at the external surface of the particle, and k xm is the<br />
mass transfer coefficient and can be obtained from the Sherwood number correlation for<br />
spheres in a convective flow (Bird et al., 1960; Field et al., 1967; Mulcahy and Smith,<br />
1969)<br />
k xm d p<br />
C f D ABf<br />
= Sh = 2.0 + 0.60Re 1/2 Sc 1/3<br />
where Sh is the Sherwood number, C f is the total gas concentration at the film<br />
10<br />
(2.7)<br />
temperature, D ABf is the molecular diffusivity at the film temperature, Re is the Reynolds<br />
number, and Sc is the Schmidt number. Note that in Eq. 2.6, the positive flux direction is<br />
designated as the direction from the bulk phase to the particle. As a result, N CO and N CO2<br />
take negative values. Using the stoichiometric relations in Eq. (2.3), N CO and N CO2 in Eq.<br />
(2.4) can be expressed in terms of N O2. Eq. 2.6 can be re-written as<br />
where =<br />
NO2 − x s NO 2 = C f DABf Sh<br />
(x∞ − x s ) (2.8)<br />
d p<br />
− 1<br />
+ 1<br />
The oxygen molar flux can be converted to carbon consumption rate q diff<br />
(gC/cm 2 /sec) by multiplying the molecular weight of carbon and the reciprocal of the<br />
stoichiometric coefficient of oxygen<br />
(2.9)