Light-off criterion and transient analysis of catalytic monoliths

被引:113
作者
Ramanathan, K
Balakotaiah, V
West, DH
机构
[1] Univ Houston, Dept Chem Engn, Houston, TX 77204 USA
[2] Dow Chem Co USA, Freeport, TX 77541 USA
关键词
ignition; light-off; monolith reactor; catalytic converter; transient time; cumulative emissions;
D O I
10.1016/S0009-2509(02)00679-6
中图分类号
TQ [化学工业];
学科分类号
0817 ;
摘要
A one-dimensional two-phase model is used to derive an analytical light-off criterion for a straight channeled catalytic monolith with washcoat, in which the flow is laminar. For the case of uniform catalyst loading and a first order reaction, the light-off criterion is given by [GRAPHICS] Here, T-f,T-in is the inlet fluid temperature, DeltaT(ad) is the adiabatic temperature rise, R-Omega is one-half the channel hydraulic radius (R-Omega =A(Omega)/P-Omega, A(Omega), P-Omega = channel cross-section area, perimeter), L is the channel length, (u) over bar is the fluid velocity, D-e is the reactant effective diffusivity in the washcoat, delta(c) is the effective washcoat thickness, k(f) is the fluid thermal conductivity and k(p)(T-f,T-in) is the first order rate constant per unit washcoat volume at the inlet fluid temperature. Nu(H,infinity) is the asymptotic Nusselt number in the channel. The function f accounts for diffusional. limitations in the washcoat and is given by f(phi) = I for 9 < 0.5 and f(phi) = 2phi for p > 0.5. The factor g(Pe(h)) depends on the solid conductivity, or more precisely, the heat Peclet number, Pe(h) = (u) over barLrho(f)c(pf)R(Omega)/delta(w)k(w), where delta(w)(k(w)) is the effective wall thickness (thermal conductivity). The function g(Pe(h)) decreases monotonically from 2.718 for Pe(h) = 0 to unity for Pe(h) = infinity. We also show that if the second term is negligible and the first exceeds unity, then ignition occurs at the back-end. If the second term exceeds unity then ignition occurs at the front-end. If the sum exceeds unity with the second term less than unity and not negligible compared to the first term then ignition occurs in the middle of the channel. This analytical ignition criterion is verified by numerical simulations using an accurate transient model that uses position dependent heat and mass transfer coefficients. We show that the plot of exit concentration versus time consists of two regions: kinetically controlled transient region and the mass transfer controlled steady-state asymptote. For the case of high solid conductivity, we present an analytical expression for the transient time at which the monolith shifts from the kinetically controlled to the mass transfer controlled regime. We also determine the influence of various parameters such as the washcoat thickness, channel dimensions, catalyst loading and initial solid temperature on this transient time and the cumulative emissions. Examination of the influence of solid conduction and channel geometry on cumulative emissions showed that designs that are optimum for steady-state operation lead to higher transient emissions and vice versa. Finally, we discuss the transient and steady-state behavior of the monolith for the special case of Le(f) < 1 (hydrogen oxidation). (C) 2003 Elsevier Science Ltd. All rights reserved.
引用
收藏
页码:1381 / 1405
页数:25
相关论文
共 26 条
[21]   LIGHT-OFF BEHAVIOR OF CATALYTIC-CONVERTERS [J].
PLEASE, CP ;
HAGAN, PS ;
SCHWENDEMAN, DW .
SIAM JOURNAL ON APPLIED MATHEMATICS, 1994, 54 (01) :72-92
[22]  
SHAH RK, 1978, LAMINAR FORCED CONNE
[23]   ADEQUACY OF LUMPED PARAMETER MODELS FOR SCR REACTORS WITH MONOLITH STRUCTURE [J].
TRONCONI, E ;
FORZATTI, P .
AICHE JOURNAL, 1992, 38 (02) :201-210
[24]  
WEI J, 1975, CATALYTIC MUFFLER, V148, P1
[25]   MATHEMATICAL-MODELS OF MONOLITH CATALYTIC-CONVERTER .1. DEVELOPMENT OF MODEL AND APPLICATION OF ORTHOGONAL COLLOCATION [J].
YOUNG, LC ;
FINLAYSON, BA .
AICHE JOURNAL, 1976, 22 (02) :331-343
[26]  
ZYGOURAKIS K, 1993, COMPUT AIDED DESIGN, P297