Nonlinear dynamics of three-way catalyst with microkinetics and internal diffusion

Similar documents
Overview of LNT Modeling Approaches

Kinetic Parameter Identification for a DOC Catalyst Using SGB test and Advanced Optimization Algorithms

Multiscale Modeling. a. Ab initio quantum chemical calculations

Prediction of Catalytic Converter Efficiency in Natural Gas Engine for Nitrogen Monoxide and Carbon Monoxide Emission Control

Notes on reaction-diffusion cases with effectiveness factors greater than one! Richard K. Herz,

AUTOMOTIVE EXHAUST AFTERTREATMENT

INTRODUCTION TO CATALYTIC COMBUSTION

Simulation of Selective Catalytic Reduction using DARS 1D Tool

2011 DOE Crosscut Workshop on Lean Emissions Reduction Simulation April 2011 Dearborn, MI

DARS Digital Analysis of Reactive Systems

UC San Diego UC San Diego Electronic Theses and Dissertations

Correlation between LNT NH 3 and N 2 O selectivities under fast cycling conditions

Part A: Operando FT-IR Studies of heterogeneous catalytic reactions: pitfalls and benefits.

A First Course on Kinetics and Reaction Engineering Unit 12. Performing Kinetics Experiments

Mercury oxidation across SCR catalyst. Karin Madsen on April 14th 2011 At McIlvaine Company Hot Topic Hour

Mathematical Investigation and CFD Simulation of Monolith Reactors: Catalytic Combustion of Methane

Multi-channel monolith reactors as dynamical systems

Supporting Information

Identification of catalytic converter kinetic model using a genetic algorithm approach

Kinetic Parameters Estimation using Vehicle Data for Exhaust Aftertreatment Devices

Department of Chemical and Petroleum Engineering Sharif University of Technology, Iran

DETAILED MODELLING OF SHORT-CONTACT-TIME REACTORS

Hydrogen addition to the Andrussow process for HCN synthesis

Numerical simulation of fluid flow in a monolithic exchanger related to high temperature and high pressure operating conditions

Modeling, Simulation and Optimization of a Cross Flow Molten Carbonate Fuel Cell

CFD study of gas mixing efficiency and comparisons with experimental data

Nonlinear Behaviour of a Low-Density Polyethylene Tubular Reactor-Separator-Recycle System

Chemical Engineering Journal

Emissions Catalyst Design using GT-SUITE. Dr. Chaitanya Sampara Viridis Chemicals Private Limited

Emissions Catalyst Design using GT-SUITE. Dr. Chaitanya Sampara Mr. Dominik Artukovic

Understanding NH 3 emissions over a three-way catalyst in lean/rich conditions.

CFD Simulation of Catalytic Combustion of Benzene

Semi-analytical solutions for cubic autocatalytic reaction-diffusion equations; the effect of a precursor chemical

Predicting the Effect of Catalyst Axial Active Site Distributions on a Diesel Oxidation Catalyst Performance

Dynamics and Control of Reactor - Feed Effluent Heat Exchanger Networks

Providing an Entry Length in Heterogeneous Catalytic Reactors with Fast Diffusion

Chaotic and non-chaotic mixed oscillations in a logistic systems with delay

A NOTE ON NON-ISOTHERMAL DIFFUSION-REACTION PROCESSES

DESIGN OF MICRO-FLUIDIC BIO-REACTORS USING TOPOLOGY OPTIMIZATION

Zeolite HC Traps: Experiments and Simulations

Fundamentals of Combustion

ONE-DIMENSIONAL HIGH-FIDELITY AND REDUCED-ORDER MODELS FOR THREE- WAY CATALYTIC CONVERTER

The Seeding of Methane Oxidation

PROCEEDINGS of the 5 th International Conference on Chemical Technology 5 th International Conference on Chemical Technology

Two Stage-CLC: a novel reactor configuration for packed bed CLC with syngas

Topical Workshop MOF Catalysis. Microkinetics in Heterogeneous Catalysis. DFG Priority Program 1362

Dehydrogenation of propane with selective hydrogen combustion: A mechanistic study by transient analysis of products

Numerical Modelling of the Original and Advanced Version of the TEMKIN-Reactor for Catalysis Experiments in Laboratory Scale

Water-Gas Shift Kinetics over Iron Oxide Catalysts at Membrane Reactor Conditions

Accuracy of Mathematical Model with Regard to Safety Analysis of Chemical Reactors*

Chemical Kinetics and Reaction Engineering

1D Simulation Modeling of SCR Catalyst at Steady State Condition

Modeling of Packed Bed Reactors: Hydrogen Production by the Steam Reforming of Methane and Glycerol

The Effect of Grid Resolution and Oxygen Storage in a One-Dimensional Monolithic Three-Way Catalyst Model

Analysis of homogeneous combustion in Monolithic structures

chromatographic reactors

Chemical Engineering Science 55 (2000) 2827}2833

Catalysis by supported metal clusters and model systems

Ugur Pasaogullari, Chao-Yang Wang Electrochemical Engine Center The Pennsylvania State University University Park, PA, 16802

USE OF DETAILED KINETIC MODELS FOR MULTISCALE PROCESS SIMULATIONS OF SULFUR RECOVERY UNITS

Overview of Reacting Flow

Advanced Transport Phenomena Module 6 - Lecture 28

Model-based analysis of TWCcoated filters performance aspects

APPLICATION OF CHEMICAL KINETICS IN THE HETEROGENEOUS CATALYSIS STUDIES

Mathematical Investigation and CFD Simulation of Monolith Reactors: Catalytic Combustion of Methane

Global kinetic modelling of the reaction of soot with NO x and O 2 on Fe 2 O 3 catalyst

Supplementary Information. The role of copper particle size in low pressure methanol synthesis via CO 2 hydrogenation over Cu/ZnO catalysts

NUMERICAL ASPECTS OF A BACTERIA GROWTH MODEL

Adsorption and conversion of various hydrocarbons on monolithic hydrocarbon adsorber

PART I Fundamentals COPYRIGHTED MATERIAL

Three Way Catalyst and Lean NO x Trap Modeling for a Lean Burn SIDI Gasoline engine

Lifetime and Performance Prediction of SOFC Anodes Operated with Trace Amounts of Hydrogen Sulfide

ChE 344 Winter 2011 Final Exam. Open Book, Notes, and Web

Extinction Limits of Premixed Combustion Assisted by Catalytic Reaction in a Stagnation-Point Flow

Numerical and experimental aspects of Temporal Analysis of Products: a tool for understanding and designing heterogeneous catalysts

Boreskov Institute of Catalysis. From the SelectedWorks of Andrey N Zagoruiko. Andrey N Zagoruiko. October, 2006

A Machine Learning Approach to Modeling and Identification of Automotive Three-Way Catalytic Converters

Diffusion and Reaction in Fe-Based Catalyst for Fischer- Tropsch Synthesis Using Micro Kinetic Rate Expressions

Interactions between oxygen permeation and homogeneous-phase fuel conversion on the sweep side of an ion transport membrane

In-Situ Differential Scanning Calorimetry in the Studies of Oxide Catalysts and Mechanisms in Oxidative Heterogeneous Catalysis

Nonlinear Operability of a Membrane Reactor for Direct Methane Aromatization

A First Course on Kinetics and Reaction Engineering Example 1.2

MODELING AND SIMULATION OF AN AUTOTHERMAL REFORMER

A Multi-Level Mathematical Model of the CO Catalytic Conversion Process

Shadow system for adsorbate-induced phase transition model

TRAP A ZERO DIMENSIONAL COMPUTER CODE FOR THE PRELIMINARY ANALYSIS OF SOOT REGENERATION PROCESS IN DIESEL PARTICULATE FILTER

Parametric Sensitivity in Chemical Systems

International Journal of Scientific Research and Innovative Technology ISSN: Vol. 4 No. 4; April CO Oxidation.

Numerical Investigation of Wall Heat Conduction Effects in Catalytic Combustion Systems

Investigation of CNT Growth Regimes in a Tubular CVD Reactor Considering Growth Temperature

SCR-Catalyst Materials for Exhaust Gas Detection D. Schönauer-Kamin, R. Moos

Non-steady-state approach t o steady-state kinetics: cas e study of H2S oxidation by oxygen

Chemical Reactions and Kinetics of the Carbon Monoxide Coupling in the Presence of Hydrogen

Kinetic Modeling of C3H6 Inhibition on NO Oxidation over Pt Catalyst

Engineering and. Tapio Salmi Abo Akademi Abo-Turku, Finland. Jyri-Pekka Mikkola. Umea University, Umea, Sweden. Johan Warna.

In situ molecular beam mass spectrometry for high temperature catalysis research: Heterogeneous Homogeneous Catalytic Oxidations

Thermodynamically Coupled Transport in Simple Catalytic Reactions

Broensted Acidity of Fibreglass Materials

GPC 2001 GLOBAL POWERTRAIN CONGRESS:

Calculation of Excess Adiabatic Compressibility from Ternary Liquid Experimental Data

Transcription:

Nonlinear dynamics of three-way catalyst P. Kočí, V. Nevoral, M. Kubíček, M. Marek Center for Nonlinear Dynamics of Chemical and Biological Systems Prague Institute of Chemical Technology Technická 5, CZ-166 28 Praha 6, Czech Republic petr.koci@vscht.cz, milos.marek@vscht.cz Tel: +420-22435 3104, Fax: +420-23333 7335

Three-way catalytic monolith reactor (TWC) converted gas CO 2, N 2, H 2 O exhaust gas CO, HC, NO X monolith channels (~1 mm) support (wall) porous washcoat γ-al 2 O 3, CeO 2 (~10 µm) catalytic metals Pt, Rh, Pd (~10 nm)

CO oxidation on Pt/CeO 2 /γ-al 2 O 3 catalyst Harmsen J.M., Schouten J.C. & Hoebink J.H., 2000, 2001

C 2 H 2 oxidation on Pt/CeO 2 /γ-al 2 O 3 catalyst Harmsen J.M., Schouten J.C. & Hoebink J.H., 2000, 2001

C 2 H 4 oxidation on Pt/CeO 2 /γ-al 2 O 3 catalyst Harmsen J.M., Schouten J.C. & Hoebink J.H., 2000, 2001

NO X transformation and reduction on Pt/CeO 2 /γ-al 2 O 3 catalyst Harmsen J.M., Schouten J.C. & Hoebink J.H., 2000, 2001

Model mass balances Bulk gas Washcoat (Pt) (Ce) Pt/CeO 2 /γ-al 2 O 3 catalyst Surface sites (Al)

Model enthalpy balances Bulk gas Washcoat Model boundary conditions

Composition of model gas at the converter inlet Component CO C 2 H 2 C 2 H 4 NO O 2 CO 2 N 2 Mol. fraction 1.22 % 280 ppm 380 ppm 1130 ppm 0.5-0.9 % 12.2 % balance

CO oxidation evolution diagram in oxygen concentration stoichiometric Outlet CO concentration in dependence on inlet O 2 concentration. CO oxidation on Pt/γ-Al 2 O 3, δ=20 µm, D eff =6x10-7 m 2.s -1, L NM =80 mol.m -3, T=630 K (isothermal). Cf. Kočí, Marek & Kubíček, Proc. of ESCAPE 13, Lappeenranta, Finland, 2003

CO oxidation solution diagram obtained by continuation stoichiometric Outlet CO concentration in dependence on inlet O 2 concentration. CO oxidation on Pt/γ-Al 2 O 3, model with no diffusional resistance within the washcoat. δ=20 µm, L NM =80 mol.m -3, T=630 K (isothermal). Only stable periodic solutions are shown. Computations performed with continuation software CONT. Cf. Kočí, Marek & Kubíček, Proc. of ESCAPE 13, Lappeenranta, Finland, 2003 Cf. also poster of Nevoral, Kočí, Schreiber, Marek & Kubíček, ESF Reactor meeting, Budapest, 2003

TWC operation coexistence of stable steady state and oscillations Outlet HC concentrations - results for different initial conditions. TWC operation of Pt/CeO 2 /γ-al 2 O 3 catalyst, δ=10 µm, D eff =2x10-6 m 2.s -1, L NM =80 mol.m -3, L OSC =1000 mol.m -3, inlet y O2 =0.74 % (small excess), T=630 K (isothermal).

TWC operation 2-periodic oscillations gas washcoat layer wall Spatiotemporal concentration profile of atomic oxygen deposited on Pt sites (O*) in the washcoat layer. TWC operation of Pt/CeO 2 /γ-al 2 O 3 catalyst. δ=10 µm, D eff =2x10-6 m 2.s -1, L NM =80 mol.m -3, L OSC =1000 mol.m -3, inlet y O2 =0.74 % (small excess), T=630 K (isothermal).

TWC operation complex oscillations Outlet CO concentrations - complex oscillatory behaviour. TWC operation of Pt/CeO 2 /γ-al 2 O 3 catalyst, δ=20 µm, D eff =2x10-6 m 2.s -1, L NM =80 mol.m -3, L OSC =100 mol.m -3, inlet y O2 =0.74 % (small excess), T=630 K (isothermal).

TWC operation forced inlet O 2 concentration X=38% X=75% X=77% Stabilized outlet NO X concentrations for periodic variation of inlet O2 concentration. Oscillation frequency f=1 Hz, amplitude A=10%. Pt/CeO 2 /γ-al 2 O 3 catalyst, δ=30 µm, D eff =1x10-6 m 2.s -1, L NM =30 mol.m -3, inlet <y O2 >=0.7375 % (stoichiometric), T in =573 K (adiabatic). Cf. Kočí, Kubíček & Marek, Proc. of IcheaP6, Pisa, Italy, 2003

TWC operation forced inlet O 2 concentration N 2 formation, NO X reduction CO 2 formation, CO and HC oxidation gas gas wall wall Spatiotemporal concentration profile of species deposited on Pt sites (N* and OCO*) in the washcoat layer. Inlet O 2 oscillation frequency f=1 Hz, amplitude A=10%. Pt/CeO 2 /γ-al 2 O 3 catalyst, D eff =1x10-6 m 2.s -1, L NM =30 mol.m -3, inlet <y O2 >=0.7375 % (stoichiometric). Cf. Kočí, Kubíček & Marek, Proc. of IcheaP6, Pisa, Italy, 2003

Conclusions Mathematical model of three-way catalytic converter (TWC) has been developed, including diffusion in the porous structure of catalytic washcoat layer and simultaneous reactions described by a complex microkinetic scheme. Multiple steady and oscillatory states (simple and doubly-periodic) and complex spatiotemporal patterns have been found for a certain range of operation parameters. The above effects arise from the used kinetic scheme, but they are influenced and modified by transport processes within the washcoat layer. Parametric continuation study of the respective system with no diffusional resistance in the washcoat has been performed. Thorough bifurcation analysis of the full model is desirable, but computationally complex (~ 10 3 stiff ODEs). TWC operation with periodic variation of inlet oxygen concentration has been simulated. It has been shown that such operation can significantly improve time-averaged outlet conversions. Acknowledgements This work has been supported by the Project MSM 223400007 of Czech Ministry of Education and ESF Reactor.

References Kočí P., Marek M. & Kubíček M., 2003, Oscillatory behaviour in mathematical model of TWC with microkinetics and internal diffusion, Proc. of 13 th Europ. Symp. on Comp. Aided Proc. Engng. (ESCAPE13), Lappeenranta, Finland, June 1-4 2003. Kočí P., Kubíček M. & Marek M., 2003, Dynamics of three-way catalyst with diffusion in the monolith washcoat, Proc. of 6 th Italian Conf. on Chem. And Proc. Engng. (ICHEAP6), Pisa, Italy, June 8-11 2003. Harmsen J., Hoebink J. & Schouten J., 2001, Acetylene and carbon monoxide oxidation over a Pt/Rh/CeO 2 /γ-al 2 O 3 automotive exhaust gas catalyst: kinetic modelling of transient experiments, Chem. Eng. Sci. 56, 2019-2035. Harmsen J., Hoebink J. & Schouten J., 2001, NO reduction by CO over automotive exhaust gas catalyst in the presence of O 2, Catal. Letters 71, 81-90. Harmsen J., Hoebink J. & Schouten J., 2001, Transient kinetic modeling of the ethylene and carbon monoxide oxidation over a commercial automotive exhaust gas catalyst, Ind. Eng. Chem. Res. 39, 599-609. Marek M., 1985, Periodic and aperiodic regimes in forced chemical oscillations, in Temporal Order, Eds. L. Rensing & N. Jaeger, Springer-Verlag, Berlin, p.105.