M-RICh v0.5 MATLAB Rate Integrator for Chemical equations
|
|
- Eleanore Reed
- 6 years ago
- Views:
Transcription
1 M-RICh v0.5 MATLAB Rate Integrator for Chemical equations Abstract: Accurately simulating the chemical conditions within a reactor tests both the stability and accuracy of a numerical scheme. M-RICh is a MATLAB-based implementation of an array of implicit and explicit finite-differencing techniques with the specific aim of simulating arbitrary chemical systems relevant to atmospheric chemistry in a user-friendly fashion. We cover the aims and development of M-RICh before testing the program with a small number of known-outcome problems, demonstrating the stability and accuracy features of the implemented solvers while demonstrating the inherent stiffness limit on chemical systems. Author: Sebastian D. Eastham MIT Department of Aeronautics and Astronautics Dates: Project start: Project report: Project presentation:
2 18.086: M-RICh Sebastian D. Eastham Contents 1 Introduction Background Mathematical basis Calculating k First, second and third order reactions Termolecular reactions Photolysis Other chemistry Consequences for numerical integration Running M-RICh Demonstration of ability Toy problem Explicit methods Implicit methods Rosenbrock Analysis Conclusions Closing remarks References 11 2
3 Sebastian D. Eastham : M-RICh 1 Introduction 1.1 Background As the scientific debate regarding pollution and climate change has grown, the need to accurately model the atmosphere has increased at an exponential rate. A typical modern atmospheric model consists of a 3-D grid, on which chemistry and transport are solved independently by operator splitting. Each part of the model presents challenges both in terms of what is known about the underlying science and in terms of solving the specific numerical problem. For example, the tendency for much of the Earth s data to be defined in terms of latitude and longitude has resulted in a prevalence of models which grid by latitude and longitude. However, this results in cells with small CFL numbers as they approach the poles, with a polar singularity which necessitates the use of complex semi-lagrangian transport schemes. However, the focus of this project is the atmospheric chemistry specifically, what are the difficulties associated with modeling the chemical processes occurring in a typical atmospheric grid box? As a general note, this report will not address the computational specifics of implementation; interested readers are directed to inspect the code itself, which is included as an appendix and can be freely acquired by request to seastham@mit.edu 1.2 Mathematical basis Chemical reactions can be represented mathematically as follows. For an example reaction aa + bb xx + yy (1) we can see that the rate of loss of each of the reactants, A and B, wil be proportional to the rate of increase of the products X and Y. Representating the number density of A in molec cm 3 as [A], we find x d[x] = y d[y] = a d[a] dt dt dt Through standard chemical kinetic theory, we then find = b d[b] dt d[x] = 1 dt x k[a]a [B] b cm 3 molec 1 s 1 where k is a rate constant independent of chemical concentrations (see section 2). Assuming that a = b = 1 as is usually the case, the example given here is a second-order reaction and therefore not linear. For n reactions involving m species we can produces one non-linear ODE for each species of the form d[x 1 ] = f i (k, [X]) cm 3 molec 1 s 1 dt Integrating these equations is then the aim of any atmospheric chemistry module. 2 Calculating k The reaction rates k are determined by a number of parameterizations. The parameters for these equations such as A and B in equation 2 or σ in equation 3 can usually be found in either recent studies or in the JPL atmospheric chemistry data publications First, second and third order reactions The basic mechanism by which molecules interact require a discrete number of reacting molecules. For reaction i involving n reactants, we can expess the rate constant k i as follows ( ) Bi k i (T ) = A i exp (cm 3 molec -1 ) n s -1 (2) T which gives use the rate of change a product [P] with respect to reaction i involving n reactants [R] j as d[p] n dt = k i (T ) [R] j molec cm 3 s 1 i j=1 3
4 18.086: M-RICh Sebastian D. Eastham The Arrhenius parameter A and activation energy B are constants which can be taken from the literature. 2.2 Termolecular reactions Some reactions are pressure-dependent; this pressure dependency is represented by including an extra reactant, [M] which represents the number density of all molecules in the local environment, resulting in a termolecular reaction. The conventional approach above would be to use a third-order reaction involving [M] as a species; however, this does not adequately take into account the mechanisms by which termolecular reactions take place. Given the very low concentrations of atmospheric species relative to conventional lab conditions, simultaneous interaction of three molecules is exceedingly rare, such that the reaction X + Y + M Z + M would almost never happen. Instead, the reaction X + Y Z occurs, with Z corresponding to an unstable, high-energy form of Z. A series of collisions with M can then remove the excess energy from Z to yield Z. Representing these reactions requires special consideration, with low-pressure and high-pressure limits calculated as ( ) n T k 0 (T ) = k0 300 cm 6 molec 2 s ( ) m T k (T ) = k 300 cm 3 molec 1 s where the 300 K limiting rates k0 300 and k 300 are given parameters, as are n and m. For cases with pure low-pressure or high-pressure limits, these can be used in isolation; however, the general formula for the termolecular reaction rate constant k f is k f ([M], T ) = k 0 (T )[M] 1 + k 0(T )[M] k (T ) 0.6 [ ( k0 (T )[M] 1 + log 10 k (T ) )] 1 cm 6 molec 2 s 1 with a slightly adjusted version used for termolecular reactions involving bound intermediates. [M] is then removed as a reactant, and the rate of the reaction R 1 + R 2 + M P would be calculated as d[p] dt = k f,i ([M], T )[R] 1 [R] 2 molec cm 3 s 1 i Further details, such as reactions with only one limiting case, and explanations of the individual terms can be found in the JPL handbook Photolysis Another relevant set of reactions are photolytic decompositions. These reactions are crucial to almost every aspect of atmospheric chemistry, since photolytic decomposition of O 2 and O 3 is responsible for the presence of OH in the atmosphere, which in turn drives processes such as oxidation of CO to CO 2. An example decomposition such as R P 1 + P 2 can be represented in the usual format as d[p] 1 dt = j i (T, p)[r] molec cm 3 s 1 i Calculation of j involves is achieved as j i (T, p) = 0 σ(λ, T, p)dλ (3) where σ(λ, T ) is the species photolytic cross section as a function of wavelength λ, temperature T and pressure p. However, calculating this online requires estimation of UV penetration to each cell, which is beyond the scope of this project. Since temperature, pressure and actinic flux can all be approximated as a function of altitude and local time, we instead use literature values for noontime j-rates at a range of altitudes, multiplied by the cosine of the solar zenith angle as an analogy for the intensity of the local UV flux. 4
5 Sebastian D. Eastham : M-RICh 2.4 Other chemistry Although the bulk of atmospheric chemistry is captured by the above mechanisms included in M-RICh, some mechanisms such as heterogeneous chemistry on and within aerosols are not included. 2.5 Consequences for numerical integration Clearly, the above problems are non-linear. For even the basic equation NO + O 3 NO 2 + O 2, the representative ODE is non-linear, requiring the operation [NO][O 3 ]. This means that any implicit method will require not only calculation and inversion of a Jacobian matrix, but also repeated iteration using Newton s method or something similar, incurring huge computational cost. Explicit methods, however, are almost never practical in the solution of such problems due to the vast differences in timescales exhibited by atmospheric chemistry. Investigation of these problems is the basis of the rest of this report. 3 Running M-RICh 3.1 Demonstration of ability M-RICh is capable of accepting any chemical mechanism with an arbitrary number of species, either live or dead, and an arbitrary number of reactions of the form given in section 2. Specifically, all of the solvers are set up to require no recoding for new reaction mechanisms; Jacobian matrix generation code, for example, is determined at solver initialization based on the spcfile.dat and rxnfile.dat specification files. Figure 1 shows the results of M-RICh integrating the set of reactions shown in table 2. This particular reaction set is a partial implementation of the mechanism in Crutzen s seminal paper investigating the mechanisms underlying the ozone layer 2, but is shown here simply to demonstrate the broad capability of M-RICh. O 3 Reaction O 2 + O( 1 D) O 3 O 2 + O( 3 P) O 2 2 O( 3 P) N 2 O N 2 + O( 1 D) Type Photolysis Photolysis Photolysis Photolysis NO 2 NO + O( 3 P) Photolysis NO 3 NO 2 + O( 3 P) Photolysis NO 3 NO + O 2 Photolysis N 2 O 5 NO 2 + NO 3 Photolysis O 3 + O( 3 P) 2 O 2 Bimolecular O( 1 D) + O 2 O( 3 P) + O 2 Bimolecular O( 1 D) + N 2 O( 3 P) + N 2 Bimolecular O( 1 D) + N 2 O 2 NO Bimolecular O( 1 D) + N 2 O N 2 + O 2 Bimolecular NO + O 3 NO 2 + O 2 Bimolecular NO 2 + O 3 NO 3 + O 2 Bimolecular O( 3 P) + O 2 + M O 3 + M Termolecular NO 3 + NO 2 + M N 2 O 5 + M Termolecular N 2 O 5 + M NO 3 + NO 2 + M Termolecular Table 1: Reactions used for the demonstration integration in figure Toy problem However, for the purpose of this investigation, we will focus on a simpler problem: the Chapman mechanism 3. This mechanism was proposed in 1930 to explain the presence of the ozone layer in the lower stratosphere. For the purposes of this report, I will be borrowing heavily from the analysis given in 5
6 18.086: M-RICh Sebastian D. Eastham Abundance (v/v) O 3 O( 3 P) O( 1 D) NO NO 2 NO 3 N 2 O Time (days) Figure 1: Integration of the system shown in table 2 over 60 days based on standard conditions at 25 km altitude (extended internation standard atmosphere). Solution found using backward Euler method with a 30 minute timestep as described in section??. Atmospheric Chemistry and Physics by Seinfeld and Pandis 4. Reaction numbering therefore follows their conventions. The relevant reactions are as follows. Using a fixed concentration of background oxygen, atomic oxygen is produced either at the ground state or in the excited form by photolysis of oxygen and ozone; it is assumed that excited atomic oxygen, O( 1 D), immediately returns to the ground state O( 3 P) for the purposes of this exercise. These photolytic decompositions proceed at the rate j O2 and j O3 respectively. O( 3 P) then reacts with oxygen through the termolecular reaction O( 3 P) + O 2 + M O 3 + M. This termolecular reaction rate constant we shall call k 2. Finally, O( 3 P) and O 3 can be returned to oxygen by the bimolecular reaction O( 3 P) + O 3 2 O 2, the rate constant for which we call k 4. From this point onwards I will refer to O( 3 P) as simply O, and will ignore O( 1 D), assuming that all ozone photolysis results in O( 3 P) either directly or almost immediately. This reaction set is shown in table??. Reaction Type Rate constant O 3 O 2 + O Photolysis j O3 O 2 2 O Photolysis j O2 O + O 2 + M O 3 + M Termolecular k 2 O 3 + O 2 O 2 Bimolecular k 4 Table 2: The Chapman mechanism We therefore have the following set of coupled ODEs: d[o] dt d[o 3 ] dt = k 2 [O][O 2 ][M] j O3 [O 3 ] k 4 [O 3 ][O] = 2j O2 [O 2 ] k 2 [O][O 2 ][M] + j O3 [O 3 ] k 4 [O 3 ][O] For a given temperature T in K and pressure p in Pa, the number density of the surrounding air can be estimated using the ideal gas law. Multiplying by Avogadro s number A and a factor of 10 6 yields [M] = pa 10 6 molec cm 3 RT where R J K 1 mol 1 is the universal gas constant. We further define [O 2 ] = 0.21 [M] on the basis that oxygen consistently makes up approximately 21% of the atmosphere by volume throughout 6
7 Sebastian D. Eastham : M-RICh the troposphere and stratosphere. If we assume that conversion between O 3 and O is rapid compared to production from photolysis of oxygen and loss by reaction 4, we can further simplify to find d[o 3 ] dt 2j O2 [O 2 ] 2k 4j O3 [O 3 ] 2 k 2 [O 2 ][M] This equation, of the form dy dt = β αy2, has an analytical solution if time variance in reaction rates and conditions is ignored specifically, if it is assumed that there is always sunlight present. Substituting for [O 2 ] as explained, and setting [O 3 ](t = 0) = 0, this yields [O 3 ](t) α = 2k 4j O3 [O 3 ] k 2 [M] 2 β = j O2 [M] ( ) 0.5 β 1 exp [ 2(ab) 0.5 t ] α 1 + exp [ 2(αβ) 0.5 t] This will be used as the reference solution. Using A variety of different integration methods have been coded into M-RICh, with a range of different accuracy and stability characteristics. For each of the coded solvers, I will perform the following analysis: Find the longest integer timestep (in seconds) which results in a stable and accurate solution to the simplified Chapman mechanism with the conditions below: Altitude: 40 km Pressure: 2.78 hpa Temperature: K Latitude: 0 Longitude: 180 Storage timestep: 30 minutes End time: 3 days Accuracy is assessed by achieving within 10% of the reference solution at all times and 2% after 3 simulated days. The timestep length, average accuracy and runtime are summarized in table 3. The interface developed for this project is shown in figure 2. Figure 2: The M-RICh GUI Please note that I have assumed reader familiarity with all of the explored methods. Although some other methods were implemented, notably Adams-Bashforth and Adams-Moulton variable order methods along with adaptive-timestep RK4, I have not tested these methods due to simple lack of time. 7
8 18.086: M-RICh Sebastian D. Eastham 3.3 Explicit methods Forward Euler (FE): FE could not converge even for a timestep of 1 second. This is a result of the very low stability of the forward Euler method, which must contend with the extreme stiffness of the problem at hand. Runge-Kutta 4th-Order (RK4): An implementation of Runge-Kutta 4 was found to remain stable for the problem when a timestep of 2 seconds was chosen. This is, in some ways, highly desirable; the increased stability of RK4 relative to forward Euler allows it to successfully integrate the problem, while the explicit nature of RK4 is easy to implement and can respond rapidly to changes in rate (see discussion in Backward Euler regarding the damping problem). However, the extremely long solution time 206 seconds for a 3 day simulation compared to for the same simulation with Backward Euler to a similar level of accuracy mean that RK4 is impractical. The general outcome from the attempt to implement explicit methods is that the problem is simply too stiff to be reasonably solved by a purely explicit approach. Although such methods have the advantage of high per-iteration computational efficiency, they are not practical for these purposes. 3.4 Implicit methods Backward Euler (BE): Backwards Euler converges easily, with a long timestep yielding good results. At a coarse testing resolution, a timestep of 2,600 seconds was found to be the marginal case, yielding an initial error around 10% which rapidly fell, converging to within 0.3% of the steady state solution. An interesting problem arises, however; implicit methods, when used with exponential problems, will inherently damp the solution. This is because they implicitly use the gradient at t + t for time t, which, in this case, is lower than the true gradient. Backward Difference variable-order (BDFn): For n = 1, backward difference schemes are identical to the backward Euler scheme. However, the order of the solver increases as more steps are included. M-RICh allows up to 5 steps to be taken to maximize accuracy, and the solution time does not appreciably increase as the accuracy is increased. Unfortunately, the increasing number of previous steps becomes a problem when dealing with initial conditions. For this project, a naive approach was taken, simply replicating the initial conditions for the phantom previous steps; however, this results in a significant damping of the solution and large maximum errors. As a result, BDF experiences larger and larger initial errors as the order increases. This resulted in a very small timestep 50 s needed to bring the peak error with 10%. Note the results mentioned previously; although RK4 produced the most accurate solution, it did so at massive computational expense. 3.5 Rosenbrock A future goal is to implement the Rosenbrock solver, a mixed implicit-explicit solver known to be effective for numerically stiff problems and used widely in the atmospheric chemistry modeling community. The Rosenbrock solver should be relatively simple to implement. 3.6 Analysis Figure 3 shows the error as a function of solution time for all of the successful solvers. This figure, along with table 3, can be used to highlight some important issues. Firstly, the problem of initial conditions is shown clearly for the BDF2 solver; the large initial error takes almost half a day to be reduced. Furthermore, the three different solvers converge on the same solution, suggesting that it may be the reference solution that is erroneous rather than the solvers themselves. A future investigation may wish to use, for example, the high-fidelity RK4 simulation rather than the approximated analytical solution as a reference. This conclusion is further reinforced by the lowest result on table 3, which shows that a small-timestep Backward Euler simulation still observes the same difference from the reference solution as the RK4 solution. The middle table, which shows the results of varying the order of the BDF solver while holding the timestep constant, is yet another reinforcement; increasing the order does not close the gap between the reference solution and the integrated solution. Rather, it simply increases the initial error, poisoning the solution for a greater period after initialization. 8
9 Sebastian D. Eastham : M-RICh Error (%) RK4 ( t = 2 s) Backward Euler ( t = 2600 s) BDF2 ( t = 50 s) Time (days) Figure 3: Error as a function of time using each solver Solver Order t (s) Max Error (%) Final Error (%) Runtime (s) Forward Euler Runge-Kutta 4 th order Backward Euler BDF Backward Differences Backward Differences Backward Euler Table 3: Project results 9
10 18.086: M-RICh Sebastian D. Eastham 4 Conclusions A chemical solver capable of taking an arbitrary chemical mechanism was implemented in MATLAB. This solver, M-RICh v0.5, is verified using an analytical solution to the simplified Chapman mechanism from atmospheric chemistry. Implicit methods are found to be hugely preferable to the nearly-unusable explicit methods due to their ability to solve stiff problems, although the damping effect which would be neutralized by use of an explicit method is found to be a problem. Future work includes implementation of a more complex chemical mechanism along with some of the unfulfilled stretch goals such as implementation of simplified transport. Future analysis would include comparison against a more reliable reference solution and implementation of the Rosenbrock solution scheme. 4.1 Closing remarks Thank you for reading this report, and I hope you have found it interesting or at least enlightening. The Chapman mechanism described above was revolutionary but significantly overestimated the ozone number density; the NO x mechanism described above was one of the first steps towards understanding why this was the case. Figure 4 shows how the two mechanisms are represented by M-RICh. I look forward to further exploring these results as part of my research. 18 x [O3] (molec/cm 3 ) Chapman NO x Time (days) Figure 4: Ozone at 40 km under sunlit conditions under 2 different mechanisms one without NO x, and the other with 10
11 Sebastian D. Eastham : M-RICh References [1] S. P. Sander, R. R. Friedl, J. R. Barker, D. M. Golden, M. J. Kurylo, G. E. Sciences, P. H. Wine, J. P. D. Abbatt, J. B. Burkholder, C. E. Kolb, G. K. Moortgat, R. E. Huie, and V. L. Orkin, Chemical Kinetics and Photochemical Data for Use in Atmospheric Studies Evaluation Number 17 NASA Panel for Data Evaluation :, no. 17, [2] P. Crutzen, Ozone production rates in an oxygen hydrogen nitrogen oxide atmosphere, Journal of Geophysical Research, vol. 76, [3] S. Chapman, XXXV. On ozone and atomic oxygen in the upper atmosphere, The London, Edinburgh, and Dublin Philosophical Magazine and Journal of Science, vol. 10, no. 26, [4] J. H. Seinfeld and S. N. Pandis, Atmospheric Chemistry and Physics. Wiley, 2 ed.,
The Chemical Kinetics Time Step a detailed lecture. Andrew Conley ACOM Division
The Chemical Kinetics Time Step a detailed lecture Andrew Conley ACOM Division Simulation Time Step Deep convection Shallow convection Stratiform tend (sedimentation, detrain, cloud fraction, microphysics)
More information0. Introduction 0.1 Concept The air / environment (geosphere): Is it a reactor? It s a matter of reactions and transports and mixing!
0. Introduction 0.1 Concept The air / environment (geosphere): Is it a reactor? It s a matter of reactions and transports and mixing! mixing times: vertically lower few kilometers (boundary layer) 1h-1d,
More informationOrdinary differential equations - Initial value problems
Education has produced a vast population able to read but unable to distinguish what is worth reading. G.M. TREVELYAN Chapter 6 Ordinary differential equations - Initial value problems In this chapter
More informationSupplement of SOA formation from the photooxidation of α-pinene: systematic exploration of the simulation of chamber data
Supplement of Atmos. Chem. Phys., 16, 278 282, 216 http://www.atmos-chem-phys.net/16/278/216/ doi:1.194/acp-16-278-216-supplement Author(s) 216. CC Attribution 3. License. Supplement of SA formation from
More informationModeling & Simulation 2018 Lecture 12. Simulations
Modeling & Simulation 2018 Lecture 12. Simulations Claudio Altafini Automatic Control, ISY Linköping University, Sweden Summary of lecture 7-11 1 / 32 Models of complex systems physical interconnections,
More informationA SPARC Success Story: The Role of Halogen Chemistry in Polar Stratospheric Ozone Depletion
A SPARC Success Story: The Role of Halogen Chemistry in Polar Stratospheric Ozone Depletion An Update on the Initiative Sponsored by the Stratospheric Processes and their Role in Climate (SPARC) Project
More informationCHM 5423 Atmospheric Chemistry Notes on kinetics (Chapter 4)
CHM 5423 Atmospheric Chemistry Notes on kinetics (Chapter 4) Introduction A mechanism is one or a series of elementary reactions that convert reactants into products or otherwise model the chemistry of
More informationLinear Multistep Methods I: Adams and BDF Methods
Linear Multistep Methods I: Adams and BDF Methods Varun Shankar January 1, 016 1 Introduction In our review of 5610 material, we have discussed polynomial interpolation and its application to generating
More information2 Reaction kinetics in gases
2 Reaction kinetics in gases October 8, 2014 In a reaction between two species, for example a fuel and an oxidizer, bonds are broken up and new are established in the collision between the species. In
More informationChapter 14 Chemical Kinetics
7/10/003 Chapter 14 Chemical Kinetics 14-1 Rates of Chemical Reactions 14- Reaction Rates and Concentrations 14-3 The Dependence of Concentrations on Time 14-4 Reaction Mechanisms 14-5 Reaction Mechanism
More informationOn Stationary state, also called steady state. Lifetimes and spatial scales of variability
On sources and sinks ATOC 3500/CHEM 3151 Week 5-6 Additional Notes February 16/18, 2016 On lifetimes, variability, and models On Stationary state, also called steady state Lifetimes and spatial scales
More informationChapter 14 Chemical Kinetics
4//004 Chapter 4 Chemical Kinetics 4- Rates of Chemical Reactions 4- Reaction Rates and Concentrations 4-3 The Dependence of Concentrations on Time 4-4 Reaction Mechanisms 4-5 Reaction Mechanism and Rate
More informationIntroduction to Chemical Kinetics AOSC 433/633 & CHEM 433/633 Ross Salawitch
Introduction to Chemical Kinetics AOSC 433/633 & CHEM 433/633 Ross Salawitch Class Web Site: http://www.atmos.umd.edu/~rjs/class/spr2013 Goals for today: Overview of Chemical Kinetics in the context of
More informationIntroduction to Chemical Kinetics AOSC 433/633 & CHEM 433 Ross Salawitch
Introduction to Chemical Kinetics AOSC 433/633 & CHEM 433 Ross Salawitch Class Web Site: http://www.atmos.umd.edu/~rjs/class/spr2017 Goals for today: Loose ends from last lecture Overview of Chemical Kinetics
More informationAnalysis Methods in Atmospheric and Oceanic Science
Analysis Methods in Atmospheric and Oceanic Science AOSC 652 Ordinary Differential Equations Week 12, Day 1 1 Differential Equations are central to Atmospheric and Ocean Sciences They provide quantitative
More informationLecture Notes to Accompany. Scientific Computing An Introductory Survey. by Michael T. Heath. Chapter 9
Lecture Notes to Accompany Scientific Computing An Introductory Survey Second Edition by Michael T. Heath Chapter 9 Initial Value Problems for Ordinary Differential Equations Copyright c 2001. Reproduction
More informationCHAPTER 1. MEASURES OF ATMOSPHERIC COMPOSITION
1 CHAPTER 1. MEASURES OF ATMOSPHERIC COMPOSITION The objective of atmospheric chemistry is to understand the factors that control the concentrations of chemical species in the atmosphere. In this book
More informationScientific Computing: An Introductory Survey
Scientific Computing: An Introductory Survey Chapter 9 Initial Value Problems for Ordinary Differential Equations Prof. Michael T. Heath Department of Computer Science University of Illinois at Urbana-Champaign
More informationAP Chemistry - Notes - Chapter 12 - Kinetics Page 1 of 7 Chapter 12 outline : Chemical kinetics
AP Chemistry - Notes - Chapter 12 - Kinetics Page 1 of 7 Chapter 12 outline : Chemical kinetics A. Chemical Kinetics - chemistry of reaction rates 1. Reaction Rates a. Reaction rate- the change in concentration
More informationAdvanced Physical Chemistry CHAPTER 18 ELEMENTARY CHEMICAL KINETICS
Experimental Kinetics and Gas Phase Reactions Advanced Physical Chemistry CHAPTER 18 ELEMENTARY CHEMICAL KINETICS Professor Angelo R. Rossi http://homepages.uconn.edu/rossi Department of Chemistry, Room
More informationCHEM Chemical Kinetics. Reaction Mechanisms
Chemical Kinetics Deri ed Rate La s from Derived Rate Laws from Reaction Mechanisms Reaction Mechanism Determine the rate law by experiment Devise a reaction mechanism If the predicted and experimental
More informationTananyag fejlesztés idegen nyelven
Tananyag fejlesztés idegen nyelven Prevention of the atmosphere KÖRNYEZETGAZDÁLKODÁSI AGRÁRMÉRNÖKI MSC (MSc IN AGRO-ENVIRONMENTAL STUDIES) Fundamentals to atmospheric chemical reactions. The stratospheric
More informationX i t react. ~min i max i. R ij smallest. X j. Physical processes by characteristic timescale. largest. t diff ~ L2 D. t sound. ~ L a. t flow.
Physical processes by characteristic timescale Diffusive timescale t diff ~ L2 D largest Sound crossing timescale t sound ~ L a Flow timescale t flow ~ L u Free fall timescale Cooling timescale Reaction
More informationCHAPTER 10 CHEMICAL KINETICS
CHAPTER 10 CHEMICAL KINETICS Introduction To this point in our study of chemistry, we have been concerned only with the composition of the equilibrium mixture, not the length of time required to obtain
More informationECE257 Numerical Methods and Scientific Computing. Ordinary Differential Equations
ECE257 Numerical Methods and Scientific Computing Ordinary Differential Equations Today s s class: Stiffness Multistep Methods Stiff Equations Stiffness occurs in a problem where two or more independent
More informationReview of Fitting Kinetic Data
L6-1 Review of Fitting Kinetic Data True or false: The goal of fitting kinetic data is to find the true rate expression. What are the two general methods used to fit kinetic data? L6-2 Advantages and Drawbacks
More informationBindel, Fall 2011 Intro to Scientific Computing (CS 3220) Week 12: Monday, Apr 18. HW 7 is posted, and will be due in class on 4/25.
Logistics Week 12: Monday, Apr 18 HW 6 is due at 11:59 tonight. HW 7 is posted, and will be due in class on 4/25. The prelim is graded. An analysis and rubric are on CMS. Problem du jour For implicit methods
More informationFrancis X. Giraldo,
1 Time-Integrators Francis X. Giraldo, giraldo@nrlmry.navy.mil, www.nrlmry.navy.mil/~giraldo/projects/nseam.html 1.1 Introduction Roughly speaking, there are 2 classes of TIs: 1. EulerianMethods(fixed-frame-e.g.,arockatthebottomofaflowing
More informationChapter 13 Lecture Lecture Presentation. Chapter 13. Chemical Kinetics. Sherril Soman Grand Valley State University Pearson Education, Inc.
Chapter 13 Lecture Lecture Presentation Chapter 13 Chemical Kinetics Sherril Soman Grand Valley State University Ectotherms Lizards, and other cold-blooded creatures, are ectotherms animals whose body
More information1. Introduction to Chemical Kinetics
1. Introduction to Chemical Kinetics objectives of chemical kinetics 1) Determine empirical rate laws H 2 + I 2 2HI How does the concentration of H 2, I 2, and HI change with time? 2) Determine the mechanism
More informationHow fast reactants turn into products. Usually measured in Molarity per second units. Kinetics
How fast reactants turn into products. Usually measured in Molarity per second units. Kinetics Reaction rated are fractions of a second for fireworks to explode. Reaction Rates takes years for a metal
More informationChapter 12. Chemical Kinetics
Chapter 12 Chemical Kinetics Section 12.1 Reaction Rates Reaction Rate Change in concentration of a reactant or product per unit time. Rate = concentration of A at time t t 2 1 2 1 concentration of A at
More informationFinite Difference Solution of the Heat Equation
Finite Difference Solution of the Heat Equation Adam Powell 22.091 March 13 15, 2002 In example 4.3 (p. 10) of his lecture notes for March 11, Rodolfo Rosales gives the constant-density heat equation as:
More informationExamples of fast and slow reactions
1 of 10 After completing this chapter, you should, at a minimum, be able to do the following. This information can be found in my lecture notes for this and other chapters and also in your text. Correctly
More informationOutline. Chemical lifetime. Photochemistry. Ozone chemistry Chapman model Catalytic cycles Ozone hole. Institute of Applied Physics University of Bern
Institute of Applied Physics University of Bern Outline Introduction Chemical reactions between stable molecules are quite slow in planetary s Absorption of solar UV-radiation leads to the production of
More information2. Sketch a plot of R vs. z. Comment on the shape. Explain physically why R(z) has a maximum in the atmospheric column.
190 PROBLEMS 10. 1 Shape of the ozone layer Consider a beam of solar radiation of wavelength λ propagating downward in the vertical direction with an actinic flux I at the top of the atmosphere. Assume
More informationT(K) k(cm 3 /molecule s) 7.37 x x x x x 10-12
CHM 5423 Atmospheric Chemistry Problem Set 3 Due date: Tuesday, February 19 th. The first hour exam is on Thursday, February 21 st. It will cover material from the first four handouts for the class. Do
More informationJim Lambers MAT 772 Fall Semester Lecture 21 Notes
Jim Lambers MAT 772 Fall Semester 21-11 Lecture 21 Notes These notes correspond to Sections 12.6, 12.7 and 12.8 in the text. Multistep Methods All of the numerical methods that we have developed for solving
More informationCS 450 Numerical Analysis. Chapter 9: Initial Value Problems for Ordinary Differential Equations
Lecture slides based on the textbook Scientific Computing: An Introductory Survey by Michael T. Heath, copyright c 2018 by the Society for Industrial and Applied Mathematics. http://www.siam.org/books/cl80
More informationAppendix C: Recapitulation of Numerical schemes
Appendix C: Recapitulation of Numerical schemes August 31, 2009) SUMMARY: Certain numerical schemes of general use are regrouped here in order to facilitate implementations of simple models C1 The tridiagonal
More informationChemical Kinetics. Chapter 13. Copyright The McGraw-Hill Companies, Inc. Permission required for reproduction or display.
Chemical Kinetics Chapter 13 Copyright The McGraw-Hill Companies, Inc. Permission required for reproduction or display. Chemical Kinetics Thermodynamics does a reaction take place? Kinetics how fast does
More informationChemical Kinetics I: The Dry Lab. Up until this point in our study of physical chemistry we have been interested in
Chemical Kinetics I: The Dry Lab Up until this point in our study of physical chemistry we have been interested in equilibrium properties; now we will begin to investigate non-equilibrium properties and
More information1.2 Derivation. d p f = d p f(x(p)) = x fd p x (= f x x p ). (1) Second, g x x p + g p = 0. d p f = f x g 1. The expression f x gx
PDE-constrained optimization and the adjoint method Andrew M. Bradley November 16, 21 PDE-constrained optimization and the adjoint method for solving these and related problems appear in a wide range of
More informationPHYSICS OF THE SPACE ENVIRONMENT
PHYSICS OF THE SPACE ENVIRONMENT PHYS/EATS 380 Winter 006 Notes Set 6 Ionospheric Electron Densities The D, E, F1 and F Layers With the advent of radio communication in the early part of the last century
More informationCHEM Chemical Kinetics. & Transition State Theory
Chemical Kinetics Collision Theory Collision Theory & Transition State Theory The rate of reaction is markedly affected by temperature. k versus T Ae E a k RT Two theories were developed to explain the
More informationMath Numerical Analysis Homework #4 Due End of term. y = 2y t 3y2 t 3, 1 t 2, y(1) = 1. n=(b-a)/h+1; % the number of steps we need to take
Math 32 - Numerical Analysis Homework #4 Due End of term Note: In the following y i is approximation of y(t i ) and f i is f(t i,y i ).. Consider the initial value problem, y = 2y t 3y2 t 3, t 2, y() =.
More informationChapter 12. Chemical Kinetics
Chapter 12 Chemical Kinetics Section 12.1 Reaction Rates Section 12.1 Reaction Rates Section 12.1 Reaction Rates Section 12.1 Reaction Rates Section 12.1 Reaction Rates Section 12.1 Reaction Rates Section
More informationELEMENTARY CHEMICAL KINETICS
ELEMENTARY CHEMICAL KINETICS EDR Chapter 25... a knowledge of the rate, or time dependence, of chemical change is of critical importance for the successful synthesis of new materials and for the utilization
More informationName AP CHEM / / Chapter 12 Outline Chemical Kinetics
Name AP CHEM / / Chapter 12 Outline Chemical Kinetics The area of chemistry that deals with the rate at which reactions occur is called chemical kinetics. One of the goals of chemical kinetics is to understand
More informationBrown et al, Chemistry, 2nd ed (AUS), Ch. 12:
Kinetics: Contents Brown et al, Chemistry, 2 nd ed (AUS), Ch. 12: Why kinetics? What is kinetics? Factors that Affect Reaction Rates Reaction Rates Concentration and Reaction Rate The Change of Concentration
More informationPreliminary report: Analyses of tcfp s potential impact on atmospheric ozone
Preliminary report: Analyses of tcfp s potential impact on atmospheric ozone Dong Wang, Seth Olsen, and Donald Wuebbles Department of Atmospheric Sciences University of Illinois, Urbana, IL 61801 Abstract
More informationIndex. higher order methods, 52 nonlinear, 36 with variable coefficients, 34 Burgers equation, 234 BVP, see boundary value problems
Index A-conjugate directions, 83 A-stability, 171 A( )-stability, 171 absolute error, 243 absolute stability, 149 for systems of equations, 154 absorbing boundary conditions, 228 Adams Bashforth methods,
More informationElementary Reactions
Elementary Reactions Elementary reactions occur in a single encounter Unimolecular: A Rate = k[a] Bimolecular: A + B Rate = k[a][b] Termolecular: A + B + C Rate = k[a][b][c] Termolecular reactions are
More informationCFC: chlorofluorocarbons
The rate of reaction is markedly affected by temperature. Chemical Kinetics & k versus T Two theories were developed to explain the temperature effects. 1. 2. 2 UV radiation strikes a CFC molecule causing
More informationFourth Order RK-Method
Fourth Order RK-Method The most commonly used method is Runge-Kutta fourth order method. The fourth order RK-method is y i+1 = y i + 1 6 (k 1 + 2k 2 + 2k 3 + k 4 ), Ordinary Differential Equations (ODE)
More informationQuantitative Understanding in Biology Module IV: ODEs Lecture I: Introduction to ODEs
Quantitative Understanding in Biology Module IV: ODEs Lecture I: Introduction to ODEs Ordinary Differential Equations We began our exploration of dynamic systems with a look at linear difference equations.
More informationAppendix B: Aqueous chemistry and gas-phase halogen chemistry
1 Appendix B: Aqueous chemistry and gasphase halogen chemistry SANFORD SILLMAN, FRANK MARSIK, KHALID I. ALWALI, GERALD J. KEELER AND MATTHEW S. LANDIS* Department of Atmospheric, Oceanic and Space Sciences
More informationChapter 14 Chemical Kinetics
Chapter 14 14.1 Factors that Affect Reaction Rates 14.2 Reaction Rates 14.3 Concentration and Rate Laws 14.4 The Change of Concentration with Time 14.5 Temperature and Rate 14.6 Reaction Mechanisms 14.7
More informationATM 507 Lecture 5. Text reading Chapter 4 Problem Set #2 due Sept. 20 Today s topics Photochemistry and Photostationary State Relation
ATM 507 Lecture 5 Text reading Chapter 4 Problem Set #2 due Sept. 20 Today s topics Photochemistry and Photostationary State Relation Beer-Lambert Law (for the absorption of light) Used to describe the
More informationKINETICS STUDY GUIDE- Written INTRODUCTION
Written Kinetics KINETICS STUDY GUIDE- Written Section: What follows is a comprehensive guide to the written component of the Chemistry 12 Provincial exam for the Unit. The questions below are from previous
More informationModule 6 : Reaction Kinetics and Dynamics Lecture 28 : Elementary Reactions and Reaction Mechanisms
Module 6 : Reaction Kinetics and Dynamics Lecture 28 : Elementary Reactions and Reaction Mechanisms Objectives In this Lecture you will learn to do the following Define what is an elementary reaction.
More informationFactors Affecting Reaction Rate
Factors Affecting Reaction Rate Outcomes: Formulate an operational definition of reaction rate. State the collision theory. Perform a lab to identify factors that affect reaction rate. Describe, qualitatively,
More informationFactors That Affect Rates. Factors That Affect Rates. Factors That Affect Rates. Factors That Affect Rates
KINETICS Kinetics Study of the speed or rate of a reaction under various conditions Thermodynamically favorable reactions DO NOT mean fast reactions Some reactions take fraction of a second (explosion)
More informationMTH 452/552 Homework 3
MTH 452/552 Homework 3 Do either 1 or 2. 1. (40 points) [Use of ode113 and ode45] This problem can be solved by a modifying the m-files odesample.m and odesampletest.m available from the author s webpage.
More informationCHEM/ENVS 380 S14, Midterm Exam ANSWERS 1 Apr 2014
PART- A. Multiple Choice Questions (5 points each): Each question may have more than one correct answer. You must select ALL correct answers, and correct answers only, to receive full credit. 1. Which
More informationInitial-Value Problems for ODEs. Introduction to Linear Multistep Methods
Initial-Value Problems for ODEs Introduction to Linear Multistep Methods Numerical Analysis (9th Edition) R L Burden & J D Faires Beamer Presentation Slides prepared by John Carroll Dublin City University
More informationENVIRONMENTAL STRUCTURE AND FUNCTION: EARTH SYSTEM - Chemistry Of The Atmosphere - I.L. Karol and A.A. Kiselev
CHEMISTRY OF THE ATMOSPHERE I.L. Karol and A.A. Main Geophysical Observatory, St. Petersburg, Russia Keywords: Atmospheric composition, gas phase reactions, heterogeneous reactions, catalytic cycles, lifetime
More informationReaction Rate. Rate = Conc. of A at t 2 -Conc. of A at t 1. t 2 -t 1. Rate = Δ[A] Δt
Kinetics The study of reaction rates. Spontaneous reactions are reactions that will happen - but we can t tell how fast. Diamond will spontaneously turn to graphite eventually. Reaction mechanism- the
More informationLecture V: The game-engine loop & Time Integration
Lecture V: The game-engine loop & Time Integration The Basic Game-Engine Loop Previous state: " #, %(#) ( #, )(#) Forces -(#) Integrate velocities and positions Resolve Interpenetrations Per-body change
More informationCONTENTS 1 MEASURES OF ATMOSPHERIC COMPOSITION
i CONTENTS 1 MEASURES OF ATMOSPHERIC COMPOSITION 1 1.1 MIXING RATIO 1 1.2 NUMBER DENSITY 2 1.3 PARTIAL PRESSURE 6 PROBLEMS 10 1.1 Fog formation 10 1.2 Phase partitioning of water in cloud 10 1.3 The ozone
More informationChemical Kinetics. Rate = [B] t. Rate = [A] t. Chapter 12. Reaction Rates 01. Reaction Rates 02. Reaction Rates 03
Chapter Chemical Kinetics Reaction Rates 0 Reaction Rate: The change in the concentration of a reactant or a product with time (M/s). Reactant Products aa bb Rate = [A] t Rate = [B] t Reaction Rates 0
More informationApplied Math for Engineers
Applied Math for Engineers Ming Zhong Lecture 15 March 28, 2018 Ming Zhong (JHU) AMS Spring 2018 1 / 28 Recap Table of Contents 1 Recap 2 Numerical ODEs: Single Step Methods 3 Multistep Methods 4 Method
More informationChapter 14: Chemical Kinetics II. Chem 102 Dr. Eloranta
Chapter 14: Chemical Kinetics II Chem 102 Dr. Eloranta Rate Laws If you are familiar with calculus Experiments would allow you to determine the reaction order and rate constant, but what if you wanted
More informationChemical Kinetics. What quantities do we study regarding chemical reactions? 15 Chemical Kinetics
Chemical Kinetics Chemical kinetics: the study of reaction rate, a quantity conditions affecting it, the molecular events during a chemical reaction (mechanism), and presence of other components (catalysis).
More informationReview for Exam 2 Ben Wang and Mark Styczynski
Review for Exam Ben Wang and Mark Styczynski This is a rough approximation of what we went over in the review session. This is actually more detailed in portions than what we went over. Also, please note
More informationCHEMICAL KINETICS. LECTURE Introduction
LECTURE-2. 2. Introduction CHEMICAL KINETICS 09//03 We have shown in the previous chapter that upon listing of the plausible reaction stoichiometries we can calculate the composition of the system in its
More information, but bursts into flames in pure oxygen.
Chemical Kinetics Chemical kinetics is concerned with the speeds, or rates of chemical reactions Chemical kinetics is a subject of broad importance. How quickly a medicine can work The balance of ozone
More informationAstronomy 8824: Numerical Methods Notes 2 Ordinary Differential Equations
Astronomy 8824: Numerical Methods Notes 2 Ordinary Differential Equations Reading: Numerical Recipes, chapter on Integration of Ordinary Differential Equations (which is ch. 15, 16, or 17 depending on
More informationMeridional Transport in the Stratosphere of Jupiter
Submitted to Astrophysical Journal Letters Meridional Transport in the Stratosphere of Jupiter Mao-Chang Liang 1, Run-Lie Shia 1, Anthony Y.-T. Lee 1, Mark Allen 1,2,A.James Friedson 2, and Yuk L. Yung
More informationCHM 5423 Atmospheric Chemistry Notes on reactions of organics in the troposphere (Chapter 5)
CHM 5423 Atmospheric Chemistry Notes on reactions of organics in the troposphere (Chapter 5) 5.1 Introduction In general, the lifetime of a molecule in the troposphere is governed by a variet of processes.
More informationChapter 12. Chemical Kinetics
Chapter 12 Chemical Kinetics Chapter 12 Table of Contents 12.1 Reaction Rates 12.2 Rate Laws: An Introduction 12.3 Determining the Form of the Rate Law 12.4 The Integrated Rate Law 12.5 Reaction Mechanisms
More informationTheoretical Models for Chemical Kinetics
Theoretical Models for Chemical Kinetics Thus far we have calculated rate laws, rate constants, reaction orders, etc. based on observations of macroscopic properties, but what is happening at the molecular
More informationMOX EXPONENTIAL INTEGRATORS FOR MULTIPLE TIME SCALE PROBLEMS OF ENVIRONMENTAL FLUID DYNAMICS. Innsbruck Workshop October
Innsbruck Workshop October 29 21 EXPONENTIAL INTEGRATORS FOR MULTIPLE TIME SCALE PROBLEMS OF ENVIRONMENTAL FLUID DYNAMICS Luca Bonaventura - Modellistica e Calcolo Scientifico Dipartimento di Matematica
More information14.4 Reaction Mechanism
14.4 Reaction Mechanism Steps of a Reaction Fred Omega Garces Chemistry 201 Miramar College 1 Reaction Mechanism The Ozone Layer Ozone is most important in the stratosphere, at this level in the atmosphere,
More informationAn Overly Simplified and Brief Review of Differential Equation Solution Methods. 1. Some Common Exact Solution Methods for Differential Equations
An Overly Simplified and Brief Review of Differential Equation Solution Methods We will be dealing with initial or boundary value problems. A typical initial value problem has the form y y 0 y(0) 1 A typical
More informationNumerical solution of ODEs
Numerical solution of ODEs Arne Morten Kvarving Department of Mathematical Sciences Norwegian University of Science and Technology November 5 2007 Problem and solution strategy We want to find an approximation
More informationOrdinary Differential Equations
Ordinary Differential Equations We call Ordinary Differential Equation (ODE) of nth order in the variable x, a relation of the kind: where L is an operator. If it is a linear operator, we call the equation
More informationCh 13 Chemical Kinetics. Modified by Dr. Cheng-Yu Lai
Ch 13 Chemical Kinetics Modified by Dr. Cheng-Yu Lai Outline 1. Meaning of reaction rate 2. Reaction rate and concentration 3. Writing a Rate Law 4. Reactant concentration and time 5. Reaction rate and
More informationNumerical Methods for the Solution of Differential Equations
Numerical Methods for the Solution of Differential Equations Markus Grasmair Vienna, winter term 2011 2012 Analytical Solutions of Ordinary Differential Equations 1. Find the general solution of the differential
More informationChemical Kinetics. Topic 7
Chemical Kinetics Topic 7 Corrosion of Titanic wrec Casón shipwrec 2Fe(s) + 3/2O 2 (g) + H 2 O --> Fe 2 O 3.H 2 O(s) 2Na(s) + 2H 2 O --> 2NaOH(aq) + H 2 (g) Two examples of the time needed for a chemical
More informationChapter 5 Exercises. (a) Determine the best possible Lipschitz constant for this function over 2 u <. u (t) = log(u(t)), u(0) = 2.
Chapter 5 Exercises From: Finite Difference Methods for Ordinary and Partial Differential Equations by R. J. LeVeque, SIAM, 2007. http://www.amath.washington.edu/ rjl/fdmbook Exercise 5. (Uniqueness for
More informationElevated atmospheric escape of atomic hydrogen from Mars induced by high-altitude water
In the format provided by the authors and unedited. SUPPLEMENTARY INFORMATION DOI: 10.1038/NGEO2887 Elevated atmospheric escape of atomic hydrogen from Mars induced by high-altitude water M. S. Chaffin,
More informationChapter: Chemical Kinetics
Chapter: Chemical Kinetics Rate of Chemical Reaction Question 1 Nitrogen pentaoxide decomposes according to equation: This first order reaction was allowed to proceed at 40 o C and the data below were
More information11/9/2012 CHEMICAL REACTIONS. 1. Will the reaction occur? 2. How far will the reaction proceed? 3. How fast will the reaction occur?
CHEMICAL REACTIONS LECTURE 11: CHEMICAL KINETICS 1. Will the reaction occur? 2. How far will the reaction proceed? 3. How fast will the reaction occur? CHEMICAL REACTIONS C(s, diamond) C(s, graphite) G
More informationKinetics - Chapter 14. reactions are reactions that will happen - but we can t tell how fast. - the steps by which a reaction takes place.
The study of. Kinetics - Chapter 14 reactions are reactions that will happen - but we can t tell how fast. - the steps by which a reaction takes place. Factors that Affect Rx Rates 1. The more readily
More information2. Outline of the MRI-EPS
2. Outline of the MRI-EPS The MRI-EPS includes BGM cycle system running on the MRI supercomputer system, which is developed by using the operational one-month forecasting system by the Climate Prediction
More informationBasic Aspects of Discretization
Basic Aspects of Discretization Solution Methods Singularity Methods Panel method and VLM Simple, very powerful, can be used on PC Nonlinear flow effects were excluded Direct numerical Methods (Field Methods)
More informationChapter 14 Chemical Kinetics
How fast do chemical processes occur? There is an enormous range of time scales. Chapter 14 Chemical Kinetics Kinetics also sheds light on the reaction mechanism (exactly how the reaction occurs). Why
More informationChemical Reaction Engineering Prof. Jayant Modak Department of Chemical Engineering Indian Institute of Science, Bangalore
Chemical Reaction Engineering Prof. Jayant Modak Department of Chemical Engineering Indian Institute of Science, Bangalore Lecture No. #40 Problem solving: Reactor Design Friends, this is our last session
More informationRate of a chemical reaction = Change in concentration Change in time
1) 2) 1) The nature of reactants and products 2) The concentration of reacting species 3) Temperature 4) Catalyst [A] Rate of a chemical reaction = Change in concentration Change in time [B] Rate of disappearance
More information