Diffusion Effectiveness factor

To have a quantitative idea of the problem of intraparticle diffusion, effectiveness factors for the two catalysts were calculated from the observed second order rate constants (based on surface area) using the "triangle method" suggested by Saterfield (4). The effectiveness factors for Monolith and Nalcomo 474 catalysts on Synthoil liquid at 371°C (700 F) were calculated to be 0.94 and 0.216, respectively. In applying the relationship between the "Thiele Modulus," 4>> and the "effectiveness factor," n> the following simplifying assumptions were made ... 

The differential material balances contain a large number of physical parameters describing the structure of the porous medium, the physical properties of the gaseous mixture diffusing through it, the kinetics of the chemical reaction and the composition and pressure of the reactant mixture outside the pellet. In such circumstances it Is always valuable to assemble the physical parameters into a smaller number of Independent dimensionless groups, and this Is best done by writing the balance equations themselves in dimensionless form. The relevant equations are (11.20), (11.21), (11.22), (11.23), (11.16) and the expression (11.27) for the effectiveness factor. 

Table II.1 which depends on the pellet size, so the familiar plot of effectiveness factor versus Thiele modulus shows how t varies with pellet radius. A slightly more interesting case arises if it is desired to exhibit the variation of the effectiveness factor with pressure as the mechanism of diffusion changes from Knudsen streaming to bulk diffusion control [66,...

At the limit of Knudsen diffusion control it is not reasonable to expect that any of the proposed approximation methods will perform well since, as we know, percentage variations in pressure are quite large. Nevertheless it is interesting to examine their results, which are shown in Figure 11 4 At this limit it is easy to check algebraically that equations (11.54) and (11.55) become the same, while (11.60) differs from the other two. Correspondingly the values of the effectiveness factor calculated using the approximation of Kehoe and Aris coincide with the results of Apecetche et al., and with the exact solution, ile Hite and Jackson s effectiveness factors differ substantially. 

Thermal transpiration and thermal diffusion will not be considered here, but it would be incorrect to assume that their influence is negligible, or even small in all circumstances. Recent results of Wong et al. [843 indi cate that they may Influence computed values of the effectiveness factor iby as much as 30. An account of thermal transpiration and thermal diffu-Ision is given in Appendix I. 

Catalyst Effectiveness. Even at steady-state, isothermal conditions, consideration must be given to the possible loss in catalyst activity resulting from gradients. The loss is usually calculated based on the effectiveness factor, which is the diffusion-limited reaction rate within catalyst pores divided by the reaction rate at catalyst surface conditions (50). The effectiveness factor E, in turn, is related to the Thiele modulus,

first-order rate constant, a the internal surface area, and the effective diffusivity. It is desirable for E to be as close as possible to its maximum value of unity. Various formulas have been developed for E, which are particularly usehil for analyzing reactors that are potentially subject to thermal instabilities, such as hot spots and temperature mnaways (1,48,51). 

Figure 10 shows that Tj is a unique function of the Thiele modulus. When the modulus ( ) is small (- SdSl), the effectiveness factor is unity, which means that there is no effect of mass transport on the rate of the catalytic reaction. When ( ) is greater than about 1, the effectiveness factor is less than unity and the reaction rate is influenced by mass transport in the pores. When the modulus is large (- 10), the effectiveness factor is inversely proportional to the modulus, and the reaction rate (eq. 19) is proportional to k ( ), which, from the definition of ( ), implies that the rate and the observed reaction rate constant are proportional to (1 /R)(f9This result shows that both the rate constant, ie, a measure of the intrinsic activity of the catalyst, and the effective diffusion coefficient, ie, a measure of the resistance to transport of the reactant offered by the pore stmcture, influence the rate. It is not appropriate to say that the reaction is diffusion controlled it depends on both the diffusion and the chemical kinetics. In contrast, as shown by equation 3, a reaction in solution can be diffusion controlled, depending on D but not on k. 

The mass transport influence is easy to diagnose experimentally. One measures the rate at various values of the Thiele modulus the modulus is easily changed by variation of R, the particle size. Cmshing and sieving the particles provide catalyst samples for the experiments. If the rate is independent of the particle size, the effectiveness factor is unity for all of them. If the rate is inversely proportional to particle size, the effectiveness factor is less than unity and

experimental points allow triangulation on the curve of Figure 10 and estimation of Tj and ( ). It is also possible to estimate the effective diffusion coefficient and thereby to estimate Tj and ( ) from a single measurement of the rate (48). 

Treatment of thermal conductivity inside the catalyst can be done similarly to that for pore diffusion. The major difference is that while diffusion can occur in the pore volume only, heat can be conducted in both the fluid and solid phases. For strongly exothermic reactions and catalysts with poor heat conductivity, the internal overheating of the catalyst is a possibility. This can result in an effectiveness factor larger than unity. 

Diffusion effects can be expected in reactions that are very rapid. A great deal of effort has been made to shorten the diffusion path, which increases the efficiency of the catalysts. Pellets are made with all the active ingredients concentrated on a thin peripheral shell and monoliths are made with very thin washcoats containing the noble metals. In order to convert 90% of the CO from the inlet stream at a residence time of no more than 0.01 sec, one needs a first-order kinetic rate constant of about 230 sec-1. When the catalytic activity is distributed uniformly through a porous pellet of 0.15 cm radius with a diffusion coefficient of 0.01 cm2/sec, one obtains a Thiele modulus y> = 22.7. This would yield an effectiveness factor of 0.132 for a spherical geometry, and an apparent kinetic rate constant of 30.3 sec-1 (106). 

The ratio of the overall rate of reaction to that which would be achieved in the absence of a mass transfer resistance is referred to as the effectiveness factor rj. SCOTT and Dullion(29) describe an apparatus incorporating a diffusion cell in which the effective diffusivity De of a gas in a porous medium may be measured. This approach allows for the combined effects of molecular and Knudsen diffusion, and takes into account the effect of the complex structure of the porous solid, and the influence of tortuosity which affects the path length to be traversed by the molecules. 

The effectiveness factor depends, not only on the reaction rate constant and the effective diffusivity, but also on the size and shape of the catalyst pellets. In the following analysis detailed consideration is given to particles of two regular shapes ... 

A hydrocarbon is cracked using a silica-alumina catalyst in the form of spherical pellets of mean diameter 2.0 mm. When the reactant concentration is 0.011 kmol/m3, the reaction rate is 8.2 x 10"2 kmol/(m3 catalyst) s. If the reaction is of first-order and the effective diffusivity De is 7.5 x 10 s m2/s, calculate the value of the effectiveness factor r). It may be assumed that the effect of mass transfer resistance in the. fluid external Lo the particles may be neglected. 

Few fixed-bed reactors operate in a region where the intrinsic kinetics are applicable. The particles are usually large to minimize pressure drop, and this means that diffusion within the pores. Steps 3 and 7, can limit the reaction rate. Also, the superficial fluid velocity may be low enough that the external film resistances of Steps 2 and 8 become important. A method is needed to estimate actual reaction rates given the intrinsic kinetics and operating conditions within the reactor. The usual approach is to define the effectiveness factor as... 

Example 10.6 A commercial process for the dehydrogenation of ethylbenzene uses 3-mm spherical catalyst particles. The rate constant is 15s , and the diffusivity of ethylbenzene in steam is 4x 10 m /s under reaction conditions. Assume that the pore diameter is large enough that this bulk diffusivity applies. Determine a likely lower bound for the isothermal effectiveness factor. 

Many theoretical embellishments have been made to the basic model of pore diffusion as presented here. Effectiveness factors have been derived for reaction orders other than first and for Hougen and Watson kinetics. These require a numerical solution of Equation (10.3). Shape and tortuosity factors have been introduced to treat pores that have geometries other than the idealized cylinders considered here. The Knudsen diffusivity or a combination of Knudsen and bulk diffusivities has been used for very small pores. While these studies have theoretical importance and may help explain some observations, they are not yet developed well enough for predictive use. Our knowledge of the internal structure of a porous catalyst is still rather rudimentary and imposes a basic limitation on theoretical predictions. We will give a brief account of Knudsen diffusion. 

The concentration of gas over the active catalyst surface at location / in a pore is ai [). The pore diffusion model of Section 10.4.1 linked concentrations within the pore to the concentration at the pore mouth, a. The film resistance between the external surface of the catalyst (i.e., at the mouths of the pore) and the concentration in the bulk gas phase is frequently small. Thus, a, and the effectiveness factor depends only on diffusion within the particle. However, situations exist where the film resistance also makes a contribution to rj so that Steps 2 and 8 must be considered. This contribution can be determined using the principle of equal rates i.e., the overall reaction rate equals the rate of mass transfer across the stagnant film at the external surface of the particle. Assume A is consumed by a first-order reaction. The results of the previous section give the overall reaction rate as a function of the concentration at the external surface, a. ... 

Equation (10.29) is the appropriate reaction rate to use in global models such as Equation (10.1). The reaction rate would be —ka if there were no mass transfer resistance. The effectiveness factor rj accounts for pore diffusion and him resistance so that the effective rate is — 7ka. 

Suppose that catalyst pellets in the shape of right-circular cylinders have a measured effectiveness factor of r] when used in a packed-bed reactor for a first-order reaction. In an effort to increase catalyst activity, it is proposed to use a pellet with a central hole of radius i /, < Rp. Determine the best value for RhjRp based on an effective diffusivity model similar to Equation (10.33). Assume isothermal operation ignore any diffusion limitations in the central hole, and assume that the ends of the cylinder are sealed to diffusion. You may assume that k, Rp, and eff are known. 

If pore diffusion is controlUng, we repeat the effectiveness factor calculations in Chapter 10. Equation (10.29) has the form of Equation (11.48), and it includes both film resistance and pore diffusion. 

Figure 2.1 Dependence of the effectiveness factor on the Thiele modulus for a first-order irreversible reaction. Steady-state diffusion and reaction, slab model, and isothermal conditions are assumed.

We define an effectiveness factor e as the ratio of the diffusion limited rate to the rate in the absence of diffusion limitation ... 

Figure S.3S. Effectiveness factor e plotted as a function of the Thiele diffusion modulus Og. The effective factor is well approximated by 3/Og for Og > 10.

In this exercise we shall estimate the influence of transport limitations when testing an ammonia catalyst such as that described in Exercise 5.1 by estimating the effectiveness factor e. We are aware that the radius of the catalyst particles is essential so the fused and reduced catalyst is crushed into small particles. A fraction with a narrow distribution of = 0.2 mm is used for the experiment. We shall assume that the particles are ideally spherical. The effective diffusion constant is not easily accessible but we assume that it is approximately a factor of 100 lower than the free diffusion, which is in the proximity of 0.4 cm s . A test is then made with a stoichiometric mixture of N2/H2 at 4 bar under the assumption that the process is far from equilibrium and first order in nitrogen. The reaction is planned to run at 600 K, and from fundamental studies on a single crystal the TOP is roughly 0.05 per iron atom in the surface. From Exercise 5.1 we utilize that 1 g of reduced catalyst has a volume of 0.2 cm g , that the pore volume constitutes 0.1 cm g and that the total surface area, which we will assume is the pore area, is 29 m g , and that of this is the 18 m g- is the pure iron Fe(lOO) surface. Note that there is some dispute as to which are the active sites on iron (a dispute that we disregard here). 

The size of the cataly.st particle influences the observed rate of reaction the smaller the particle, the less time required for the reactants to move to the active catalyst sites and for the products to diffuse out of the particle. Furthermore, with relatively fast reactions in large particles the reactants may never reach the interior of the particle, thus decreasing the catalyst utilization. Catalyst utilization is expressed as the internal effectiveness factor //,. This factor is defined as follows ... 

In order for diffusional limitations to be negligible, the effectiveness factor must be close to 1, i.e. nearly complete catalyst utilization, which requires that the Thiele modulus is suffieiently small (< ca. 0.5), see Figure 3.32. Therefore, the surface-over-volume ratio must be as large as possible (particle size as small as possible) from a diffusion (and heat-transfer) point of view. There are many different catalyst shapes that have different SA/V ratios for a given size. 

If diffusion of reactants to the active sites in pores is slower than the chemical reaction, internal mass transfer is at least partly limiting and the reactant concentration decreases along the pores. This reduces the reaction rate compared to the rate at external surface conditions. A measure of the reaction rate decrease is the effectiveness factor, r, which has been defined as ... 

Inspection of Fig. 15.3 reveals that while for jo 0.1 nAcm , the effectiveness factor is expected to be close to 1, for a faster reaction with Jo 1 p,A cm , it will drop to about 0.2. This is the case of internal diffusion limitation, well known in heterogeneous catalysis, when the reagent concentration at the outer surface of the catalyst grains is equal to its volume concentration, but drops sharply inside the pores of the catalyst. In this context, it should be pointed out that when the pore size is decreased below about 50 nm, the predominant mechanism of mass transport is Knudsen diffusion [Malek and Coppens, 2003], with the diffusion coefficient being less than the Pick diffusion coefficient and dependent on the porosity and pore stmcture. Moreover, the discrete distribution of the catalytic particles in the CL may also affect the measured current owing to overlap of diffusion zones around closely positioned particles [Antoine et ah, 1998]. 

---