Diffusion/second-order kinetics, numerical

Notice that the molar density of key-limiting reactant A on the external surface of the catalytic pellet is always used as the characteristic quantity to make the molar density of component i dimensionless in all the component mass balances. This chapter focuses on explicit numerical calculations for the effective diffusion coefficient of species i within the internal pores of a catalytic pellet. This information is required before one can evaluate the intrapellet Damkohler number and calculate a numerical value for the effectiveness factor. Hence, 50, effective is called the effective intrapellet diffusion coefficient for species i. When 50, effective appears in the denominator of Ajj, the dimensionless scaling factor is called the intrapellet Damkohler number for species i in reaction j. When the reactor design focuses on the entire packed catalytic tubular reactor in Chapter 22, it will be necessary to calcnlate interpellet axial dispersion coefficients and interpellet Damkohler nnmbers. When there is only one chemical reaction that is characterized by nth-order irreversible kinetics and subscript j is not required, the rate constant in the nnmerator of equation (21-2) is written as instead of kj, which signifies that k has nnits of (volume/mole)"" per time for pseudo-volumetric kinetics. Recall from equation (19-6) on page 493 that second-order kinetic rate constants for a volnmetric rate law based on molar densities in the gas phase adjacent to the internal catalytic surface can be written as... 

Fig. 15.7 Calculation of the equilibrium between oxygen consumption occurring during the oxidation of sedimentary organic matter near the sediment surface and resupply of oxygen by means of diffusion. Here, oxygen consumption is assumed to follow first-order kinetics. Application of Pick s Second Law in an explicit numeric solution permits a reliable calculation of the times required for the adjustment of a stationary condition.

TABLE 18-1 Numerical Solution of the Mass Transfer Equation for One-Dimensional Diffusion and Second-Order Irreversible Chemical Kinetics in Porous Catalysts ... 

The heterogeneous rate law in (22-57) is dimensionalized with pseudo-volumetric nth-order kinetic rate constant k that has units of (volume/mol)" per time. k is typically obtained from equation (22-9) via surface science studies on porous catalysts that are not necessarily packed in a reactor with void space given by interpellet. Obviously, when axial dispersion (i.e., diffusion) is included in the mass balance, one must solve a second-order ODE instead of a first-order differential equation. Second-order chemical kinetics are responsible for the fact that the mass balance is nonlinear. To complicate matters further from the viewpoint of obtaining a numerical solution, one must solve a second-order ODE with split boundary conditions. By definition at the inlet to the plug-flow reactor, I a = 1 at = 0 via equation (22-58). The second boundary condition is d I A/df 0 as 1. This is known classically as the Danckwerts boundary condition in the exit stream (Danckwerts, 1953). For a closed-closed tubular reactor with no axial dispersion or radial variations in molar density upstream and downstream from the packed section of catalytic pellets, Bischoff (1961) has proved rigorously that the Danckwerts boundary condition at the reactor inlet is... 

This is a system of stiff, second-order partial differential equations which can be solved numerically to yield both transient, and steady state concentration profiles within the layer. Comparison of the experimental calibration curves and of the time response curves with the calculated ones provides the verification of the proposed model from which it is possible to determine the optimum thickness of the enzyme layer. Because the Thiele modulus is the controlling parameter in the diffusion-reaction equation it is obvious from Eq.6 that the optimum thickness L will depend on the other constants and functions included in the Thiele modulus. Because of this the optimum thickness will vary from one kinetic scheme to another. 

In Chapter 14, we discussed the case of a single-component band. In practice, there are almost always several components present simultaneously, and they have different mass transfer properties. As seen in Chapter 4, the equilibrium isotherms of the different components of a mixture depend on the concentrations of all the components. Thus, as seen in Chapters 11 to 13, the mass balances of the different components are coupled, which makes more complex the solution of the multicomponent kinetic models. Because of the complexity of these models, approximate analytical solutions can be obtained only under the assumption of constant pattern conditions. In all other cases, only numerical solutions are possible. The problem is further complicated because the diffusion coefficients and the rate constants depend on the concentrations of the corresponding components and of all the other feed components. However, there are still relatively few papers that discuss this second form of coupling between component band profiles in great detail. In most cases, the investigations of mass transfer kinetics and the use of the kinetic models of chromatography in the literature assume that the rate constants and the diffusion coefficients are concentration independent. This seems to be an acceptable first-order approximation in many cases, albeit separation problems in which more sophisticated theoretical approaches are needed begin to appear as the accuracy of measru ments improve and more interest is paid to complex... 

---