Orthogonal collocation is used to analyze the dispersion and chemical conversion of a multicomponent fluid pulse in a laminar isothermal flow in a catalytic tube. Homogeneous and heterogeneous first-order chemical reactions are included, and a full matrix diffusion law is used. The collocation method represents the dynamics well over wide ranges of system parameters and reaction time. Radial averaging approaches prove less convenient, since the dispersion coefficients thus introduced are complicated functions of time, transport properties and chemical kinetics.
IntroductionDispersion effects are widely encountered in chemical reactors and mass transfer systems. In process models, such systems are often treated as one-dimensional by use of transversely averaged mass and energy balances. The local velocity and thermodynamic state in such a model stand for averages over a crosssection of the flow, and the convective dispersion (caused by deviations from these averages) is modelled as enhanced longitudinal diffusion and conduction.The transverse-average approach has been developed extensively for laminar flows in tubes. Taylor (1953) and Aris (1956) derived radially averaged formulas from the long-time asyrnp totes of their two-dimensional problems. Gill and Sankarasubramanian (1 970, 1973, 1974) introduced a set of time-dependent dispersion coefficients Ki(t), and treated first-order reacting systems. Johns (1978, 1979) gave more detailed transient solutions, using a new set of dispersion coefficients & ( t ) based on axial moments of the solute concentration. Multicomponent reacting systems have been studied by DeGance and Johns (1980, 1985, 1987, Aris (1980), andHatton (1981), but only partial results have been obtained.Transversely averaged equations have limited predictive power, since the dispersion coefficients vary with time and with the form of the given problem. The long-time asymptotic forms of these equations are simpler, but also less useful in the presence of fast surface reactions (Aris, 1980). For predictive purposes, an efficient simulation of the process in its full set of dimensions should be preferred. We demonstrate and test the latter approach here, extending the collocation scheme of Wang and Stewart (1983) to reactive multicomponent systems.
Problem StatementA multicomponent stream in steady, developed flow through a tube is initially in chemical equilibrium at a uniform composition wo. At time T = 0 a small mass of miscible fluid is injected. The resulting mass-fraction disturbance vector, Au = w -uo, of dimension N is described by the mass-balance equation
CollocationAn approximate solution for Aw is sought by orthogonal collocation (Villadsen and Stewart, 1967;Stewart, 1984). The vector Aw in Eqs. 1 -4 is approximated by the function gives a modified differential system, expressible in partitioned matrix form as an anEquation 4 then takes the radially discretized form, various cases of which are analyzed below. This completes the restatement of Eqs. 1-4 as a collocation problem. No...