A Review of Kinetic Modeling Methodologies for Complex Processes Luís Pereira De Oliveira, Damien Hudebine, Denis Guillaume, Jan Verstraete

Total Page:16

File Type:pdf, Size:1020Kb

A Review of Kinetic Modeling Methodologies for Complex Processes Luís Pereira De Oliveira, Damien Hudebine, Denis Guillaume, Jan Verstraete A Review of Kinetic Modeling Methodologies for Complex Processes Luís Pereira de Oliveira, Damien Hudebine, Denis Guillaume, Jan Verstraete To cite this version: Luís Pereira de Oliveira, Damien Hudebine, Denis Guillaume, Jan Verstraete. A Review of Kinetic Modeling Methodologies for Complex Processes. Oil & Gas Science and Technology - Revue d’IFP Energies nouvelles, Institut Français du Pétrole, 2016, 71 (3), pp.45. 10.2516/ogst/2016011. hal- 01395195 HAL Id: hal-01395195 https://hal.archives-ouvertes.fr/hal-01395195 Submitted on 10 Nov 2016 HAL is a multi-disciplinary open access L’archive ouverte pluridisciplinaire HAL, est archive for the deposit and dissemination of sci- destinée au dépôt et à la diffusion de documents entific research documents, whether they are pub- scientifiques de niveau recherche, publiés ou non, lished or not. The documents may come from émanant des établissements d’enseignement et de teaching and research institutions in France or recherche français ou étrangers, des laboratoires abroad, or from public or private research centers. publics ou privés. Oil & Gas Science and Technology – Rev. IFP Energies nouvelles (2016) 71,45 Ó L.P. de Oliveira et al., published by IFP Energies nouvelles, 2016 DOI: 10.2516/ogst/2016011 Dossier Methodology for Process Development at IFP Energies nouvelles Méthodologies pour le développement de procédés à IFP Energies nouvelles A Review of Kinetic Modeling Methodologies for Complex Processes Luís P. de Oliveira, Damien Hudebine, Denis Guillaume and Jan J. Verstraete* IFP Energies nouvelles, Rond-point de l’échangeur de Solaize, BP 3, 69360 Solaize - France e-mail: [email protected] * Corresponding author Abstract — In this paper, kinetic modeling techniques for complex chemical processes are reviewed. After a brief historical overview of chemical kinetics, an overview is given of the theoretical background of kinetic modeling of elementary steps and of multistep reactions. Classic lumping techniques are introduced and analyzed. Two examples of lumped kinetic models (atmospheric gasoil hydrotreating and residue hydroprocessing) developed at IFP Energies nouvelles (IFPEN) are presented. The largest part of this review describes advanced kinetic modeling strategies, in which the molecular detail is retained, i.e. the reactions are represented between molecules or even subdivided into elementary steps. To be able to retain this molecular level throughout the kinetic model and the reactor simulations, several hurdles have to be cleared first: (i) the feedstock needs to be described in terms of molecules, (ii) large reaction networks need to be automatically generated, and (iii) a large number of rate equations with their rate parameters need to be derived. For these three obstacles, molecular reconstruction techniques, deterministic or stochastic network generation programs, and single-event micro-kinetics and/or linear free energy relationships have been applied at IFPEN, as illustrated by several examples of kinetic models for industrial refining processes. Résumé — Une revue de méthodes de modélisation cinétique pour des procédés complexes — Dans cet article, les techniques de modélisation cinétique des processus chimiques complexes sont examinées. Après un bref aperçu historique de la cinétique chimique, un aperçu des bases théoriques de la modélisation cinétique d’étapes élémentaires et de réactions globales est présenté. Les techniques classiques de regroupement (lumping) sont ensuite présentées et analysées. Deux exemples de modèles cinétiques regroupés (pour l’hydrotraitement de gazole atmosphérique et pour l’hydrotraitement de résidus) développés à IFP Energies nouvelles (IFPEN) sont présentés. La plus grande partie de cette revue décrit des stratégies avancées de modélisation cinétique, dans lesquelles le détail moléculaire est retenu : les réactions entre les molécules sont représentées ou même subdivisées en étapes élémentaires. Pour être en mesure de conserver ce niveau moléculaire à la fois dans le modèle cinétique et dans les simulations de réacteurs, plusieurs obstacles doivent d’abord être éliminés : (i) la charge doit être décrite en termes de molécules, (ii) les grands réseaux réactionnels doivent être générés automatiquement et (iii) un grand nombre d’équations de vitesse avec leurs paramètres de vitesse doit être dérivé. Pour ces trois obstacles, des techniques de reconstruction moléculaire, des programmes de génération de réseaux déterministes ou stochastiques, et des modèles microcinétiques basés sur des événements constitutifs (single events) et/ou des This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited. Page 2 of 49 Oil & Gas Science and Technology – Rev. IFP Energies nouvelles (2016) 71,45 relations linéaires d’énergie libre (Linear Free Energy Relationships) ont été utilisés à IFPEN, comme illustré par plusieurs exemples de modèles cinétiques pour des procédés de raffinage industriels. LIST OF ABBREVIATIONS NIST National Institute of Standards and Technology NMR Nuclear Magnetic Resonance AEBP Atmospheric Equivalent Boiling Point ODE Ordinary Differential Equation AED Atomic Emission Detector PCP Protonated CycloPropane AGO Atmospheric Gas Oil PCB Protonated CycloButane ARO Aromatics PDE Partial Differential Equation ASP Asphaltenes PDF Probability Distribution Function ASTM American Society for Testing Materials PES Potential Energy Surface BT BenzoThiophene(s) Pr Protonation CARB Carbazoles QFER Quadratic Free-Energy Relationship CME Chemical Master Equation QSAR Quantitative Structure-Activity Relationships DBT DiBenzoThiophene(s) QSRC Quantitative Structure-Reactivity Correlations DFRN Deterministic Full Reaction Network QSS Quasi-Steady State DFT Density Functional Theory RDS Rate Determining Step DH Dehydrogenation REM Reconstruction by Entropy Maximization DI Di-aromatics RES Resins DIPPR Design Institute for Physical PRoperties RRK Rice–Ramsperger–Kassel DLRN Deterministic Limited Reaction Network SARA Saturates/Aromatics/Resins/Asphaltenes FCC Fluid Catalytic Cracking SAT Saturates FID Flame Ionization Detector SCD Sulfur Chemiluminescence Detector FT-ICR-MS Fourier Transform – Ion Cyclotron Resonance SEC Size Exclusion Chromatography – Mass Spectrometry SEK Single Event Kinetics GC Gas Chromatography SLRN Stochastically Limited Reaction Network HDS HydroDeSulfurization SOL Structure Oriented Lumping HPLC High-Performance Liquid Chromatography SR Stochastic Reconstruction HS Hydride Shift SSA Stochastic Simulation Algorithm IUPAC International Union of Pure and Applied TRC Thermodynamics Research Center Chemistry TRI Tri-aromatics kMC Kinetic Monte-Carlo TS Transition State LC Liquid Chromatography TST Transition State Theory LCO Light Cycle Oil VR Vacuum Residue LFER Linear Free-Energy Relationships LHSV Liquid Hourly Space Velocity LSODE Livermore Solver for Ordinary Differential LIST OF SYMBOLS Equations MARI Most Abundant Reaction Intermediate MASI Most Abundant Surface Intermediate Ai Pre-exponential factor for the rate constant of MC Monte-Carlo reaction i MONO Mono-aromatics c0 Unit concentration MS Mass Spectrometry ci Stochastic rate constant for reaction i MTHS Molecular Type and Homologous Series Cj Concentration of species j approach E Shannon’s information entropy criterion NCD Nitrogen Chemiluminescence Detector Ea Activation energy NES Nearly Empty Surface Ei Activation energy for reaction i Oil & Gas Science and Technology – Rev. IFP Energies nouvelles (2016) 71,45 Page 3 of 49 Fk Mixing rule for analytical property k Dn Difference in molecularity between the reverse h Planck’s constant (6.626068 9 10À34 m2 kg sÀ1) and the forward reaction D 6¼ – Ik Intermediate k n 1 molecularity of the reaction that forms the k Rate constant activated complex (TST) ~ D k Intrinsic single-event rate coefficient of the SR Entropy of reaction elementary step DS6¼ Activation entropy (Transition State Theory) k ki Rate constant for reaction i Reorganization energy (Quadratic Free-Energy kB Boltzmann’s constant Relationship) À23 2 À2 À1 (1.3806503 9 10 m kg s K ) mj Stoichiometric coefficient for component j K Number of analytical properties or constraints mij Stoichiometric coefficient for component j in Ka Ionization constant reaction i q ’ Kc Equilibrium constant (concentration basis) H Hammett s sensitivity constant r Keq Overall equilibrium constant i Stoichiometric number for reaction/elementary step i Ki Equilibrium constant for reaction / elementary r step i AB Reaction cross-section (collision theory) r ’ Mj Molecular weight of component j H Hammett s substitution constant rr MAB Reduced molecular weight of components A and Symmetry number for the reactant B (collision theory) r6¼ Symmetry number for the activated complex h ne Number of single events i Reactivity index of molecule i N Number of components hij Reactivity index j of molecule i NA Avogadro constant (6.0221408577 9 1023 molÀ1) NAlkyl Number of alkyl substituents in b position SUBSCRIPTS of a sulfur atom c Concentration basis NAR Number of aromatic rings eq At equilibrium conditions NSR Number of saturated rings glob For global symmetry and chirality NTR Number of thiophene rings i For reaction i P Pressure int Intrinsic Pk Value of analytical property k j For component j r Reaction rate overall For the overall reaction ri Reaction rate of reaction
Recommended publications
  • Rate Constants and Kinetic Isotope Effects for Methoxy Radical Reacting with NO2 and O2 Jiajue Chai, Hongyi Hu, Theodore
    Rate Constants and Kinetic Isotope Effects for Methoxy Radical Reacting with NO2 and O2 Jiajue Chai, Hongyi Hu, Theodore. S. Dibble* Department of Chemistry, College of Environmental Science and Forestry, State University of New York, Syracuse, New York 13210, United States Geoffrey. S. Tyndall,* and John. J. Orlando Atmospheric Chemistry Division, National Center for Atmospheric Research, Boulder, Colorado 80307, United States Abstract Relative rate studies were carried out to determine the temperature dependent rate constant ratio k1/k2a: CH3O• + O2 → HCHO + HO2• (R1) CH3O• + NO2 (+M) → CH3ONO2 (+ M) (R2a) over the temperature range 250-333 K in an environmental chamber at 700 Torr using FTIR detection. Absolute rate constants k2 were determined using laser flash photolysis/laser induced fluorescence under the same conditions. The analogous experiments were carried out for the reactions of perdeuterated methoxy radical (CD3O•). Absolute rate constants k2 were in excellent agreement with the recommendations of the JPL Data Evaluation panel. The combined data (i.e., 3 -1 k1/k2 and k2) allow determination of k1 as cm s , corresponding to 1.4 × 10-15 cm3 s-1 at 298 K. The rate constant at 298 K is in excellent agreement with previous work, but the observed temperature dependence is less than previously reported. The deuterium isotope effect, kH/kD, can be expressed in Arrhenius form as The deuterium isotope effect does not appear to be greatly influenced by tunneling, which is consistent with previous theoretical work by Hu and Dibble (Hu, H.; Dibble, T. S., J. Phys. Chem. A 2013, 117, 14230–14242). Keywords: methylnitrite, formaldehyde, methylnitrate, alkoxy radicals, kinetic isotope effect * [email protected], [email protected] 1 I.
    [Show full text]
  • Tables of Rate Constants for Gas Phase Chemical Reactions of Sulfur Compounds (1971-1979)
    NBSIR 80-2118 Tables of Rate Constants for Gas Phase Chemical Reactions of Sulfur Compounds (1971-1979) Francis Westley Chemical Kinetics Division Center for Thermodynamics and Molecular Science National Measurement Laboratory National Bureau of Standards U.S. Department of Commerce Washington. D.C. 20234 July 1980 Interim Report Prepared for Morgantown Energy Technology Center Department of Energy Morgantown, West Virginia 26505 Fice of Standard Reference Data itional Bureau of Standards ashington, D.C. 20234 c, 2 Jl 1 m * rational bureau 07 STANDARDS LIBRARY NBSIR 80-2118 wov 2 4 1980 TABLES OF RATE CONSTANTS FOR GAS f)Oi Clcc -Ore PHASE CHEMICAL REACTIONS OF SULFUR COMPOUNDS (1971-1979) 1 5 (o lo-z im Francis Westley Chemical Kinetics Division Center for Thermodynamics and Molecular Science National Measurement Laboratory National Bureau of Standards U.S. Department of Commerce Washington, D.C. 20234 July 1980 Interim Report Prepared for Morgantown Energy Technology Center Department of Energy Morgantown, West Virginia 26505 and Office of Standard Reference Data National Bureau of Standards Washington, D.C. 20234 U.S. DEPARTMENT OF COMMERCE, Philip M. Klutznick, Secretary Luther H. Hodges, Jr., Deputy Secretary Jordan J. Baruch, Assistant Secretary for Productivity, Technology, and Innovation NATIONAL BUREAU OF STANDARDS. Ernest Ambler, Director . i:i : a 1* V'.OTf in. te- ?') t aa.‘ .1 J Table of Contents Abstract 1 Introduction 2 Guidelines for the User 4 Table of Arrhenius Parameters for Chemical Reactions of Sulfur Compounds 12 Reactions of: S(Sulfur atom) 12 $2(Sulfur dimer) 20 S0(Sulfur monoxide) 21 S02(Sulfur dioxide 23 S0 (Sulfur trioxide) 28 3 S20(Sul fur oxide^O) ) 29 SH (Mercapto free radical) 29 H2S (Hydroqen sulfide 31 CS (Carbon monosulfide free radical) 32 CS2 (Carbon disulfide) 33 COS (Carbon oxide sulfide) 35 CH^S* ( Methyl thio free radical) 36 CH^SH (Methanethiol ) 36 cy-CF^CF^S (Thiirane) 37 CH^SCF^* (Methyl, (methyl thio)-, free radical) ...
    [Show full text]
  • Classical and Modern Methods in Reaction Rate Theory
    J. Phys. Chem. 1988, 92, 3711-3725 3711 loss associated with the transition between states is offset by a For the Re( MQ+)-based complexes an additional complication gain in energy in the intramolecular vibrations (eq 3). to intramolecular electron transfer exists arising from a "flattening" On the other hand, if (AGe,,2 + A0,2) > (AG,,l + A,,,) (Figure of the relative orientations of the two rings of the MQ+ ligand 4B) and nuclear tunneling effects in the librational modes are when it accepts an However, the analysis presented unimportant, intramolecular electron transfer can not occur in here does provide a general explanation for the role of free energy a frozen environment. The situation is the same as in Figure 3A. change in environments where dipole motions are restricted. For Even direct Re(1) - MQ' excitation (process a in Figure 4B) a directed, intramolecular electron transfer such as reaction 1, must either lead to the ground state by emission or to the Re- the change in electronic distribution between states will, in general, (11)-bpy'--based MLCT state by reverse electron transfer: lead to a difference between Ao,l and A0,2. This difference creates a solvent dipole barrier to intramolecular electron transfer. It [ (bpy)Re"(CO) 3( MQ')] 2+* -+ (bpy'-)Re"( CO) 3(MQ')] 2'* follows from eq 4 that if light-induced electron transfer is to occur (5) in a frozen environment, the difference in solvent reorganizational The inversion in the sense of the electron transfer in the frozen energy, A& = bS2- A,,, must be compensated for by a favorable environment is induced energetically by the greater solvent re- free energy change, A(AG) = - AGes,l,as organizational energy for the lower excited state, A0,2 > Ao,l - Aces,~- AGes.2 > Ao-2 - &,I (AGes,2 - AGes,l)* or In a macroscopic sample, a distribution of solvent dipole en- vironments exists around the assembly of ground states.
    [Show full text]
  • Eyring Activation Energy Analysis of Acetic Anhydride Hydrolysis in Acetonitrile Cosolvent Systems Nathan Mitchell East Tennessee State University
    East Tennessee State University Digital Commons @ East Tennessee State University Electronic Theses and Dissertations Student Works 5-2018 Eyring Activation Energy Analysis of Acetic Anhydride Hydrolysis in Acetonitrile Cosolvent Systems Nathan Mitchell East Tennessee State University Follow this and additional works at: https://dc.etsu.edu/etd Part of the Analytical Chemistry Commons Recommended Citation Mitchell, Nathan, "Eyring Activation Energy Analysis of Acetic Anhydride Hydrolysis in Acetonitrile Cosolvent Systems" (2018). Electronic Theses and Dissertations. Paper 3430. https://dc.etsu.edu/etd/3430 This Thesis - Open Access is brought to you for free and open access by the Student Works at Digital Commons @ East Tennessee State University. It has been accepted for inclusion in Electronic Theses and Dissertations by an authorized administrator of Digital Commons @ East Tennessee State University. For more information, please contact [email protected]. Eyring Activation Energy Analysis of Acetic Anhydride Hydrolysis in Acetonitrile Cosolvent Systems ________________________ A thesis presented to the faculty of the Department of Chemistry East Tennessee State University In partial fulfillment of the requirements for the degree Master of Science in Chemistry ______________________ by Nathan Mitchell May 2018 _____________________ Dr. Dane Scott, Chair Dr. Greg Bishop Dr. Marina Roginskaya Keywords: Thermodynamic Analysis, Hydrolysis, Linear Solvent Energy Relationships, Cosolvent Systems, Acetonitrile ABSTRACT Eyring Activation Energy Analysis of Acetic Anhydride Hydrolysis in Acetonitrile Cosolvent Systems by Nathan Mitchell Acetic anhydride hydrolysis in water is considered a standard reaction for investigating activation energy parameters using cosolvents. Hydrolysis in water/acetonitrile cosolvent is monitored by measuring pH vs. time at temperatures from 15.0 to 40.0 °C and mole fraction of water from 1 to 0.750.
    [Show full text]
  • Empirical Chemical Kinetics
    Empirical Chemical Kinetics time dependence of reactant and product concentrations e.g. for A + 2B → 3C + D d[A] 1 d[B] 1 d[C] d[D] rate =− =− = = d2d3ddtttt =ν In general, a chemical equation is 0∑ ii R i nt()− n (0 ) The extent of a reaction (the advancement) = ii ν i For an infinitesmal advancement dξ each reactant/product =ν ξ concentration changes by d[Ri ]i d . d1d[R]ξ By definition, the rate == i ν ddtti Reaction rates usually depend on reactant concentrations, e.g., rate= k [A]xy [B] order in B total order = x+y rate constant In elementary reaction steps the orders are always integral, but they may not be so in multi-step reactions. The molecularity is the number of molecules in a reaction step. Rate Laws d[A] Zero order: −=k dt 0 [A] = [A]00 -[A]t kt [A]0 t1 = 2 2k 0 0 t −=d[A] First order: k1[A] [A] ln[A] dt − = kt1 [A]t [A]0 e ln 2 0 t 0 t t1 = 2 k1 d[A] 2 Second order: −=2[A]k [A] dt 2 [A] [A] = 0 t + 12[A]kt20 0 t -1=+ -1 [A]t [A]02 2kt 1 t1 = 2 2[A]k20 [A]-1 0 t Second-order Kinetics: Two Reactants A + B k → products da rate=− =kab a = [A], b =[B] dt =− − ka()()00 x b x ddxa Since aa=− x, =− 0 ddtt dx So =−ka()() x b − x dt 00 x dx kt = ∫ ()()−− 0 axbx00 −111x =−dx ()−−−∫ ()() ab000 ax 0 bx 0 −1 ab = ln 0 ()− ab00 ab 0 or more usefully, aa=+−0 ln ln ()abkt00 bb0 Data Analysis “Classical” methods of data analysis are often useful to explore the order of reactions, or to display the results (e.g.
    [Show full text]
  • Kinetics and Atmospheric Chemistry
    12/9/2017 Kinetics and Atmospheric Chemistry Edward Dunlea, Jose-Luis Jimenez Atmospheric Chemistry CHEM-5151/ATOC-5151 Required reading: Finlayson-Pitts and Pitts Chapter 5 Recommended reading: Jacob, Chapter 9 Other reading: Seinfeld and Pandis 3.5 General Outline of Next 3 Lectures • Intro = General introduction – Quick review of thermodynamics • Finlayson-Pitts & Pitts, Chapter 5 A. Fundamental Principles of Gas-Phase Kinetics B. Laboratory Techniques for Determining Absolute Rate Constants for Gas-Phase Reactions C. Laboratory Techniques for Determining Relative Rate Constants for Gas-Phase Reactions D. Reactions in Solution E. Laboratory Techniques for Studying Heterogeneous Reactions F. Compilations of Kinetic Data for Atmospheric Reactions 1 12/9/2017 Kinetics and Atmospheric Chemistry • What we’re doing here… – Photochemistry already covered – We will cover gas phase kinetics and heterogeneous reactions – Introductions to a few techniques used for measuring kinetic parameters • What kind of information do we hope to get out of “atmospheric kinetics”? – Predictive ability over species emitted into atmosphere • Which reactions will actually proceed to products? • Which products will they form? • How long will emitted species remain before they react? Competition with photolysis, wash out, etc. • Pare down list of thousands of possible reactions to the ones that really matter – Aiming towards idea practical predictive abilities • Use look up tables to decide if reaction is likely to proceed and determine an “atmospheric lifetime”
    [Show full text]
  • On the Determination of the Reaction Rate Constant and Selectivity in Gas and Liquid-Phase Organic Reactions: Temperature and Solvent Effects
    On the Determination of the Reaction Rate Constant and Selectivity in Gas and Liquid-Phase Organic Reactions: Temperature and Solvent Effects Aikaterini Diamanti A thesis submitted for the Doctor of Philosophy degree of Imperial College London and the Diploma of Imperial College. Centre for Process Systems Engineering Department of Chemical Engineering Imperial College London London SW7 2AZ United Kingdom 9th December 2016 Declaration of Originality I, Aikaterini Diamanti, certify that this thesis, and the research to which it refers, are the product of my own work, and that any ideas or quotations from the work of other people, published or otherwise, are fully acknowledged in accordance with the standard referencing practices of the discipline. Copyright The copyright of this thesis rests with the author and is made available under a Creative Commons Attribution Non-Commercial No Derivatives licence. Researchers are free to copy, distribute or transmit the thesis on the condition that they attribute it, that they do not use it for commercial purposes and that they do not alter, transform or build upon it. For any reuse or distribution, researchers must make clear to others the licence terms of this work. 2 Abstract Chemical reactions occur abundantly in nature and the rates at which they proceed are critically influenced by factors such as the temperature and the solvent medium in which they take place. These two factors do not only affect the rate of the reaction but may also have a significant impact on other metrics, such as the selectivity and the catalytic activity, as well as on the overall performance of a process application.
    [Show full text]
  • Electrochemical Hydrogen Evolution and Kinetic Isotope Effect
    Title ELECTROCHEMICAL HYDROGEN EVOLUTION AND KINETIC ISOTOPE EFFECT Author(s) TSIONSKY, V. M.; KRISHTALIK, L. I. Citation JOURNAL OF THE RESEARCH INSTITUTE FOR CATALYSIS HOKKAIDO UNIVERSITY, 28(3), 371-380 Issue Date 1981-06 Doc URL http://hdl.handle.net/2115/25106 Type bulletin (article) File Information 28(3)_P371-380.pdf Instructions for use Hokkaido University Collection of Scholarly and Academic Papers : HUSCAP .T. Res. In st. Catalysis, Hokkaido Unic·., 1'01. 28, No.8, j)p. 871 to 880 (1980) ELECTROCHEMICAL HYDROGEN EVOLUTION AND KINETIC ISOTOPE EFFECT By V. M. TSIONSKY,*) and L. 1. KRISHTALIK*) (Received September 8, 1980) Abstract It has been considered the relationship between the isotope effect (separation factor) of electrochemical hydrogen evolution and the set of reaction steps. On that basis has been shown the possibility of the isotope effect analysis in solving of many problems of hydrogen evolution kinetics: elementary act of a reaction, stepwise mechanism of molecular hydrogen formation, influence of ~~l-potential on the kinetics, problems associated with the nature of intermediates and proton donors. Professor Horiuti was the first to use hydrogen isotopes for investigation of the reaction mechanism at the hydrogen electrode. In his lecture at the lapan-Soviet Seminar in Sapporo in October 1974 he gave a vivid picture of the initiation of this direction of research early in the thirties. v Since then Horituti and his School have carried out numerous studies in which a number of problems of the mechanism of hydrogen evolution reaction were investigated by the isotope exchange and the kinetic isotope effect methods (for literature, see e.
    [Show full text]
  • Chapter 13. Chemical Kinetics
    Chapter 13. Chemical Kinetics What we will learn: • The rate of a reaction • The rate law • The relation between reactant concentration and time • Activation energy • Reaction energy • Reaction mechanism • Catalysis GCh13-1 Chemical kinetics: • The area of chemistry investigating the speed (rate) of chemical reactions Reactants Products Reaction rate: The change in the concentration of a reactant or product with time (M/s) M-molarity During a chemical reaction the concentration of reactants decreases and the concentration of products increases GCh13-2 D[A] The change in the reactant D[A] concentration (molarity) rate = Dt Dt The change in time GCh13-3 Reaction between Br2 and HCOOH: Br2(aq) + HCOOH(aq) 2HBr(aq) + CO2(g) D [Br2] average rate = - Dt [Br2]final - [Br2]initial = - tfinal - tinitial Average rate = the reaction rate averaged at a certain period of time Instantaneous rate = the reaction rate at a specific instant in time GCh13-4 k1,k2 the instantaneous rates GCh13-5 The rate constant: The rate constant is a constant of proportionality between the reaction rate and the concentration of a reactant rate = k[A] k - rate constant The reaction rate is proportional to the The rate constant “k” is concentration of the reactant. The more specific for each reaction reactants, the bigger probability that the reaction occurs GCh13-6 Reaction rates and stoichiometry 2A B Two moles of the reactant A disappear, and one mole of the product B is formed at the same time. The reaction rate can be expressed as: GCh13-7 General formula Example Consider the reaction 4NO2(g) + O2(g) g 2N2O5(g) At a particular moment of the reaction, O2 is reacting at the rate 0.024 M/s.
    [Show full text]
  • Pang, Genny (2012). Experimental Determination of Rate Constants For
    EXPERIMENTAL DETERMINATION OF RATE CONSTANTS FOR REACTIONS OF THE HYDROXYL RADICAL WITH ALKANES AND ALCOHOLS A DISSERTATION SUBMITTED TO THE DEPARTMENT OF MECHANICAL ENGINEERING AND THE COMMITTEE ON GRADUATE STUDIES OF STANFORD UNIVERSITY IN PARTIAL FULFILLMENT OF THE REQUIREMENTS FOR THE DEGREE OF DOCTOR OF PHILOSOPHY Genny Anne Pang August 2012 © 2012 by Genny Anne Pang. All Rights Reserved. Re-distributed by Stanford University under license with the author. This work is licensed under a Creative Commons Attribution- Noncommercial 3.0 United States License. http://creativecommons.org/licenses/by-nc/3.0/us/ This dissertation is online at: http://purl.stanford.edu/rh393tq1232 ii I certify that I have read this dissertation and that, in my opinion, it is fully adequate in scope and quality as a dissertation for the degree of Doctor of Philosophy. Ronald Hanson, Primary Adviser I certify that I have read this dissertation and that, in my opinion, it is fully adequate in scope and quality as a dissertation for the degree of Doctor of Philosophy. Craig Bowman I certify that I have read this dissertation and that, in my opinion, it is fully adequate in scope and quality as a dissertation for the degree of Doctor of Philosophy. David Golden Approved for the Stanford University Committee on Graduate Studies. Patricia J. Gumport, Vice Provost Graduate Education This signature page was generated electronically upon submission of this dissertation in electronic format. An original signed hard copy of the signature page is on file in University Archives. iii iv Abstract Over one quarter of the energy usage in the United States currently occurs in the transportation sector.
    [Show full text]
  • New Method of Determining Kinetic Parameters for Decomposition of Hydrogen Peroxide by Catalase
    catalysts Article New Method of Determining Kinetic Parameters for Decomposition of Hydrogen Peroxide by Catalase Ilona Trawczy´nska Department of Chemical and Bioprocess Engineering, Faculty of Chemical Technology and Engineering, University of Science and Technology in Bydgoszcz, Seminaryjna 3, 85-326 Bydgoszcz, Poland; [email protected] Received: 19 February 2020; Accepted: 9 March 2020; Published: 12 March 2020 Abstract: The presented study investigates the kinetic properties of catalase during hydrogen peroxide decomposition reaction. A novel and simple method is hereby proposed for the determination of the enzyme deactivation rate constant (kd) and the decomposition of H2O2 reaction rate constant (kr). Available methods allow the kd constant to be determined only based on previously experimentally determined kr. The presented method differs from the conventional procedure. Known initial and final concentrations of hydrogen peroxide enable determination of both constants at the same time based on data from only one experiment. The correctness of the new method proposed here in determining the reaction rate constant was checked by comparing the obtained constant values with the calculated values according to the commonly used Aebi method. The method was used to analyze in detail the effect of pH (3–10) and temperature (10–45 ◦C) of the reaction medium on kinetic constants. The value of the constant kd increases together with the value of pH and temperature. In addition, the activation energy for decomposition reaction and deactivation reaction was found to be 1 1 Er = 14 kJ mol− and Ed = 56.8 kJ mol− respectively. Keywords: catalase; kinetic parameters; enzyme deactivation; enzyme activity 1.
    [Show full text]
  • Reaction Rate Coefficient K20 and Temperature Coefficient Θ in Organic Waste Thermal Disintegration
    ITM Web of Conferences 23, 00026 (2018) https://doi.org/10.1051/itmconf/20182300026 XLVIII Seminar of Applied Mathematics Reaction rate coefficient k20 and temperature coefficient Θ in organic waste thermal disintegration Sylwia Myszograj1,* 1Univeristy of Zielona Góra, Faculty of Civil Engineering, Architecture and Environmental Engineering, Poland Abstract. It was described the test of sewage sludge and organic fraction of municipal mixed solid waste thermal disintegration process. The waste activated sludge used during the tests was collected from the secondary settlement tank in a mechanical-biological wastewater treatment plant. The biowaste used in the studies was collected from an area of new buildings. It was noticed from means values of Soluble Chemical Oxygen Demand (SCOD) plot that both heating temperature and time, influence the amount of dissolved COD. The observations indicate that changes of SCOD can be described by an increasing, differentiable function of time and the rate of change of the soluble COD in the hydrolysates, in time is proportional to the difference of the maximum values of SCOD and its value in time, which leads to the relationship of the first-order ordinary differential equation. The process effectiveness depending on the temperature was described with the mathematical model including Van't Hoff-Arrhenius equation. Inspection of the data and some preliminary fits indicates, that for the description of changes in SCOD terms of time and temperature were adopted the form of nonlinear mixed model. Values of k20 indicator and Θ parameter depend on the substrate type. For waste activated sludge thermal disintegration, value of reaction speed indicator k20 was 0.028 h-1 (0,67 d-1), and value of temperature indicator equalled Θ = 1.024.
    [Show full text]