Modelling reactive distillation

R. Taylor, R. Krishna / Chemical Engineering Science 55 (2000) 5183}5229 5207a function of the binary MS di!usion coe$cients,the mixture composition, and the mass transferrates themselves. However, if, as is usually done, weassume the e!ective di!usivity to be constant, then analyticalsolutions to Eqs. (18) and (19) may be derived.Indeed, the literature on mass transfer with chemicalreaction based on these two equations is vast and suchequations often are used. Results usually are expressed inthe formN "k E (x !x ), (20)where k is the mass transfer coe$cient and E is anenhancement factor that accounts for the in#uence ofchemical reaction on the molar #uxes. Among otherthings, enhancement factors depend quite strongly on thereaction kinetics; thus, there is no universal expressionfor E that can be used in all situations. Indeed, in manycases it is not even possible to obtain analytical expressionsfor the enhancement factor, even for the muchsimpli"edconstitutive relation given by Eq. (19).Danckwerts (1970) and Van Swaaij and Versteeg (1992)review the "eld of mass transfer with chemical reaction.Frank, Kuipers, Krishna and van Swaaij (1995a) citeother studies of non-isothermal mass transfer with simultaneouschemical reaction modelled using the simplerconstitutive relation given by Eq. (19). E!ective di!usivitymodels of this sort have been used with some successin the modelling of amine-based gas treating processes(see, for example, Cornelisse, Beenackers, van Beckum& van Swaaij, 1980; Carey, Hermes & Rochelle, 1991;Altiqi, Sabri, Bouhamra & Alper, 1994 and the extensivelists of references in these papers).A more rigorous approach to modelling mass transferin multicomponent systems is a!orded by the generalizedFick's law: xJ "!c D z , (21)where J is the molar di!usion #ux of speciesi (N "J #x N ) and the D are the Fick di!usion coe$cients for the multicomponent mixture. The Fickdi!usion coe$cients can be related to the binary MSdi!usion coe$cients; see Taylor and Krishna (1993) fordetails of this. Solutions of Eq. (21) and the continuityequation (18) are known for a few cases involving simultaneousdi!usion and reaction (see, for example, Toor,1965; Kenig & Gorak, 1995, 1997). In most cases ofpractical signi"cance Eq. (21) is just as di$cult to solve asthe full MS equations and it is necessary to employnumerical methods here as well.A rigorous approach to modelling heterogeneous systems(Fig. 5(b)) would involve a complete description ofmass transport to the catalyst and di!usion and reactioninside the catalyst particles, if the catalyst is porous, orreaction at the surface if it is not. In either case it isunnecessary to allow for reaction in the vapor or liquid"lms as described above.For all cases involving a heterogeneous reaction wemay model transport through the vapor}liquid interfaceusing the MS equations given above (assuming no reactionin the liquid "lm). We also need to model di!usion ofreactants to and products away from the solid catalystsurface using a separate set of di!usion equations. Forthose cases where the reaction takes place at the catalystsurface the boundary conditions at the liquid}solid interfacewill be determined by the reaction kinetics. However,if the reaction takes place inside a porous catalyst thenwe need to model di!usion and reaction inside the catalystas well as transport from the bulk liquid to/from thecatalyst surface.Modelling multicomponent mass transfer in porousmedia is complicated when the mean free path length ofthe molecules is of the order of magnitude of the porediameter. The di$culties posed by this case may becircumvented by a method originally introduced byMaxwell (1866) and developed further by Mason andco-workers (see Mason & Malinauskas, 1983). Maxwellsuggested that the porous material itself be describedas a supplementary `dusta species, consisting ofvery large molecules that are kept motionless bysome unspeci"ed external force. The Chapman}Enskogkinetic theory is then applied to the new pseudo-gasmixture, in which the interaction between the dust andgas molecules simulates the interaction between the solidmatrix and the gas species. In addition, one is no longerfaced with the problem of #ux and composition variationsacross a pore and problems related to catalystgeometry.The dusty #uid model as developed by Krishna andWesselingh (1997) is a modi"cation of the dusty gasmodel so as to be able to model liquid-phase di!usion inporous media. For a non-ideal mixture we have:x μ ¹ z # x ¹ < p z #x ηB D " x N !x N ! N , (22)c n c D where D is the e!ective Knudsen di!usion coe$cient forspecies i in the porous catalyst.The mass transfer rates are obtained by multiplyingthe #uxes by the interfacial area of the catalyst particles.This is not as straightforward, as it looks, since, dependingon the geometry of the catalyst, the cross-sectionalarea can change along the di!usion path. It is necessaryto take catalyst geometry into account in the solution ofthese equations.The dusty gas model is often used as the basis for thecalculation of a catalyst e!ectiveness factor (Jackson,1977). The extension to non-ideal #uids noted above hasnot been used as often, partly because of the greater