Chemical Kinetics and the Rössler System. 1 Introduction. 2 The NH 3 - HCl reaction. Dynamics at the Horsetooth Volume 2, 2010.

Save this PDF as:
 WORD  PNG  TXT  JPG

Size: px
Start display at page:

Download "Chemical Kinetics and the Rössler System. 1 Introduction. 2 The NH 3 - HCl reaction. Dynamics at the Horsetooth Volume 2, 2010."

Transcription

1 Dynamics at the Horsetooth Volume 2, 21. Chemical Kinetics and the Rössler System Department of Mathematics Colorado State University Report submitted to Prof. P. Shipman for Math 4, Fall 21 Abstract. Chemical oscillations present questions about both the chemical mechanism behind the oscillations and the potential for mathematical models. The reaction of ammonia with hydrocloric acid, which is know to produce Liesegang ring patterns [1], exhibits oscillation patterns of mathematical interest [6]. This paper describes a mathematical model for this chemical system, as well as a comparison of experimental data to the well-known Rössler System. Keywords: Diffusion Equation, Law of Mass Action, Rössler System, Periodicity, Chaos 1 Introduction During the summer of 21 at Colorado State University, I was awarded a CIMS grant and worked with Dr. Patrick Shipman of the Department of Mathematics and Dr. Stephen Thompson of the Department of Chemistry. We studied the reaction of ammonia with hyrdrocloric acid, a reaction that forms an interesting pattern which falls under the broad classification of Liesegang Rings [1]. My interest in this problem continued into the fall of 21, where I took Math 4, Dynamical Systems with Dr. Patrick Shipman, and I realized that interesting connections could be made from my summer work to concepts learned in class. In this paper, I will summarize some of the results obtained during my summer research which includes a mathematical model of the NH 3 -HCl system. Then I will discuss the Rössler System, a well-known system of equations which exhibits very nice periodic patterns. I will conclude with a comparison of actual experimental data from Tim Lenczycki s work [6] to both the mathematical model and the Rössler System. 2 The NH 3 - HCl reaction For the modeling of this system, we considered a tube with ammonia (NH 3 ) on one end of the tube and hydrocloric acid (HCl) at the other end, as was the case in the data which we were trying to model [6]. We will denote [HCl] = a and [NH 3 ]=b. Since these two chemicals will react while in their gas phase, we certainly know that they will both obey the diffusion equation a t = D 2 a a t 2

2 b t = D 2 b b t 2 where D a and D b are diffusion coefficients. Now we consider the forward reaction of the system, where NH 3 (g) + HCl(g) NH 4 Cl(s) with reaction rate k 1. Denote [NH 4 Cl] = s. Then, by the Law of Mass Action, we have the following set of equations: a t = D 2 a a t 2 k 1ab b t = D 2 b b t 2 k 1ab s t = k 1ab. These equations are the basic building blocks of our model. We now proceed to include more of the details of this system. We know that both heterogeneous and homogeneous nucleation are important in this chemical reaction [6]. Mathematically, we needed to model the large jump in the reaction rate that occurs once ab exceeds a certain amount η. In Figure 1, we see an example of this where η = k 1 H(ab ) ab Figure 1: Jump in the reation rate as a result of nucleation. We modeled this behavior using the following equations a t = D 2 a a t 2 k 1H(ab η)ab b t = D 2 b b t 2 k 1H(ab η)ab s t = k 1H(ab η)ab. Dynamics at the Horsetooth 2 Vol. 2, 21

3 Using these equations, we were able to see oscillations in our solutions as seen in Figures 2, 3, and 4. Figure 2: The reaction front of the NH 3 - HCl system as a function of position and time. Figure 3: A top view of the reaction front; the movement of the reaction front. Figure 4: Diffusion of NH 3 and HCl from either sides of the tube. The oscillations are a promising sign that we are at least heading in the correct direction, Dynamics at the Horsetooth 3 Vol. 2, 21

4 although we should certainly test our model against actual experimental data. Future modifications to this model are discussed in the conclusion of this paper. 3 The Rössler System We now turn to a well-known system of equations, called the Rössler System. The Rössler equations are as follows dx = (y + z) dt dy dt = x + ay dz = b + xz cz dt where a, b, and c are parameters. If we let a = b =.2 and let c vary and we choose initial conditions x() = y() = z() = 1, we will see that we get very interesting patterns. When c =2.3, we can see that we get period-one behavior as shown in Figures and 6. Note that this period-one behavior, which begins after about t = 2, could be represented by a simple oscillatory function such as sin(x). The Rossler System 4 3 z(t) 2 1 y(t) Figure : Period-one limit cycle when c = 2.3. Dynamics at the Horsetooth 4 Vol. 2, 21

5 t Figure 6: Period-one behavior with respect to. When c =3.3, as shown in Figures 7 and 8, we see period-two behavior here. Note that we get two distinct amplitudes here as shown in Figure 8. The Rossler System 8 6 z(t) 4 2 y(t) Figure 7: Period-two limit cycle when c = 3.3. Dynamics at the Horsetooth Vol. 2, 21

6 t Figure 8: Period-two behavior with respect to. In the same way, we can let c =.3 and get period-three behavior, as shown in Figures 9 and 1. The Rossler System 2 1 z(t) 1 1 y(t) Figure 9: Period-three limit cycle when c =.3. Dynamics at the Horsetooth 6 Vol. 2, 21

7 t Figure 1: Period-three behavior with respect to. Finally, if we consider the case where c = 6.3, as in Figure 11, we can see that a slight change in the initial conditions has caused a great change in the solution. In the top part of Figure 11, we have initial conditions x() = y() = z() = 1 while in the bottom part, we have x() = 1.1,y() = z() = 1. Notice that the difference in the two solutions begins right around t =. z(t) z(t) y(t) y(t) The Rossler System The Rossler System t t Figure 11: The system exhibits chaotic behavior when c = 6.3. Dynamics at the Horsetooth 7 Vol. 2, 21

8 4 Highlights of the Experimental Data Experimental Data was collected for the NH 3 - HCl reaction which measures the position of the reaction front as time progressed [6]. It was noted that the reaction front of this system oscillates and oftentimes in a very periodic manner. We now discuss some specific examples of this. During a specific time segment of the NH 3 - HCl reaction, the reaction front oscillated as shown in Figure 12, which appears to exhibit period-one behavior. This is confirmed by the Fourier Analysis in Figure 13 since we obtain one clear peak.. Figure 12:. Figure 13: In a different time segment of the reaction, we find more periodic behavior. According to the Fourier Analysis, it appears that during this time segment, we have period-two behavior as shown in Figures 14 and 1. Dynamics at the Horsetooth 8 Vol. 2, 21

9 . Figure 14:. Figure 1: Finally, chaotic behavior of the oscillation front is apparent during another time segment of the reaction. This is shown in Figures 16 and 17. Dynamics at the Horsetooth 9 Vol. 2, 21

10 . Figure 16:. Figure 17: Conclusion This project certainly has many avenues of possiblity for future work. First, we would like to add a fourth equation to our NH 3 - HCl model to account for the fact that the NH 4 Cl is, at first, believed to be in the gas phase until enough NH 4 Cl gas molecules form a solid. If we denote [NH 4 Cl (g)] = m (for monomer) and [NH 4 Cl (s)] = s (for solid), these equations should look something like a t = D 2 a a t 2 k 1ab b t = D 2 b b t 2 k 1ab Dynamics at the Horsetooth 1 Vol. 2, 21

11 m t = D 2 m m t 2 + k 1ab k 2 H(m η) s t = k 2H(m η) k 3 s. We would also like to consider the idea of autocatalysis within our model which, simply put, is the idea that the more NH 4 Cl (g) there is, the faster the reaction rate will be. We have also done simple experiment of this NH 3 - HCl Reaction but in a vertical tube. This has resulted in a tornado-like reaction where oscillations are visible but a mathematical model for this would have to take into account things such as fluid flow and the effect that it has on these oscillations. Dynamics at the Horsetooth 11 Vol. 2, 21

12 References [1] Dee, G.T. Patterns Produced by Precipitation at a Moving Reaction Front. Physical Review Letters 7.3 (1986): [2] Hirsch, Morris W., Stephen Smale, and Robert L. Devaney. Differential Equations, Dynamical Systems, and an Introduction to Chaos. Elsevier, 24. [3] Holmes, Mark H. Introduction to the Foundations of Applied Mathematics. New York: Springer, 29. [4] Howard, P. Partial Differential Equations in MATLAB 7.. Spring 2. [] Lebedeva, M.I., D.G. Vlachos, and M. Tsapatsis. Bifurcation Analysis of Liesegang Ring Pattern Formation. Physical Review Letters 92.8 (24): 8831(1)-8831(4). [6] Lenczycki, Timothy T. (Thesis) Oscillations in Gas Phase Periodic Precipitation Systems: The NH 3 -HCl Story. Spring 23. [7] Lynch, Stephen. Dynamical Systems with Applications using MATLAB. Boston: Birkhauser, 24. [8] Strogatz, Steven H. Nonlinear Dynamics and Chaos with Applications to Physics, Biology, Chemistry, and Engineering. Addison-Wesley, [9] Verhulst, Ferdinand. Nonlinear Differential Equations and Dynamical Systems. Second edition. New York: Springer-Verlag Berlin Heidelberg, Dynamics at the Horsetooth 12 Vol. 2, 21

DIFFERENTIAL EQUATIONS, DYNAMICAL SYSTEMS, AND AN INTRODUCTION TO CHAOS

DIFFERENTIAL EQUATIONS, DYNAMICAL SYSTEMS, AND AN INTRODUCTION TO CHAOS DIFFERENTIAL EQUATIONS, DYNAMICAL SYSTEMS, AND AN INTRODUCTION TO CHAOS Morris W. Hirsch University of California, Berkeley Stephen Smale University of California, Berkeley Robert L. Devaney Boston University

More information

Why are Discrete Maps Sufficient?

Why are Discrete Maps Sufficient? Why are Discrete Maps Sufficient? Why do dynamical systems specialists study maps of the form x n+ 1 = f ( xn), (time is discrete) when much of the world around us evolves continuously, and is thus well

More information

Calculating Fractal Dimension of Attracting Sets of the Lorenz System

Calculating Fractal Dimension of Attracting Sets of the Lorenz System Dynamics at the Horsetooth Volume 6, 2014. Calculating Fractal Dimension of Attracting Sets of the Lorenz System Jamie Department of Mathematics Colorado State University Report submitted to Prof. P. Shipman

More information

NONLINEAR DYNAMICS PHYS 471 & PHYS 571

NONLINEAR DYNAMICS PHYS 471 & PHYS 571 NONLINEAR DYNAMICS PHYS 471 & PHYS 571 Prof. R. Gilmore 12-918 X-2779 robert.gilmore@drexel.edu Office hours: 14:00 Quarter: Winter, 2014-2015 Course Schedule: Tuesday, Thursday, 11:00-12:20 Room: 12-919

More information

Lab 5: Nonlinear Systems

Lab 5: Nonlinear Systems Lab 5: Nonlinear Systems Goals In this lab you will use the pplane6 program to study two nonlinear systems by direct numerical simulation. The first model, from population biology, displays interesting

More information

Some Dynamical Behaviors In Lorenz Model

Some Dynamical Behaviors In Lorenz Model International Journal Of Computational Engineering Research (ijceronline.com) Vol. Issue. 7 Some Dynamical Behaviors In Lorenz Model Dr. Nabajyoti Das Assistant Professor, Department of Mathematics, Jawaharlal

More information

Report submitted to Prof. P. Shipman for Math 540, Fall 2009

Report submitted to Prof. P. Shipman for Math 540, Fall 2009 Dynamics at the Horsetooth Volume 1, 009. Three-Wave Interactions of Spin Waves Aaron Hagerstrom Department of Physics Colorado State University aaronhag@rams.colostate.edu Report submitted to Prof. P.

More information

MATH 312 Section 8.3: Non-homogeneous Systems

MATH 312 Section 8.3: Non-homogeneous Systems MATH 32 Section 8.3: Non-homogeneous Systems Prof. Jonathan Duncan Walla Walla College Spring Quarter, 2007 Outline Undetermined Coefficients 2 Variation of Parameter 3 Conclusions Undetermined Coefficients

More information

An Undergraduate s Guide to the Hartman-Grobman and Poincaré-Bendixon Theorems

An Undergraduate s Guide to the Hartman-Grobman and Poincaré-Bendixon Theorems An Undergraduate s Guide to the Hartman-Grobman and Poincaré-Bendixon Theorems Scott Zimmerman MATH181HM: Dynamical Systems Spring 2008 1 Introduction The Hartman-Grobman and Poincaré-Bendixon Theorems

More information

Math 240: Spring-mass Systems

Math 240: Spring-mass Systems Math 240: Spring-mass Systems Ryan Blair University of Pennsylvania Tuesday March 1, 2011 Ryan Blair (U Penn) Math 240: Spring-mass Systems Tuesday March 1, 2011 1 / 15 Outline 1 Review 2 Today s Goals

More information

Oscillations in Damped Driven Pendulum: A Chaotic System

Oscillations in Damped Driven Pendulum: A Chaotic System International Journal of Scientific and Innovative Mathematical Research (IJSIMR) Volume 3, Issue 10, October 2015, PP 14-27 ISSN 2347-307X (Print) & ISSN 2347-3142 (Online) www.arcjournals.org Oscillations

More information

Unit VI Stoichiometry. Applying Mole Town to Reactions

Unit VI Stoichiometry. Applying Mole Town to Reactions Unit VI Stoichiometry Applying Mole Town to Reactions Learning Goals I can apply mole town to reactions to determine the amount of product based on the amount of a reactant. I can apply mole town to reaction

More information

Math 308 Week 8 Solutions

Math 308 Week 8 Solutions Math 38 Week 8 Solutions There is a solution manual to Chapter 4 online: www.pearsoncustom.com/tamu math/. This online solutions manual contains solutions to some of the suggested problems. Here are solutions

More information

2D-Volterra-Lotka Modeling For 2 Species

2D-Volterra-Lotka Modeling For 2 Species Majalat Al-Ulum Al-Insaniya wat - Tatbiqiya 2D-Volterra-Lotka Modeling For 2 Species Alhashmi Darah 1 University of Almergeb Department of Mathematics Faculty of Science Zliten Libya. Abstract The purpose

More information

Nonlinear systems, chaos and control in Engineering

Nonlinear systems, chaos and control in Engineering Nonlinear systems, chaos and control in Engineering Module 1 block 3 One-dimensional nonlinear systems Cristina Masoller Cristina.masoller@upc.edu http://www.fisica.edu.uy/~cris/ Schedule Flows on the

More information

Introduction to Dynamical Systems Basic Concepts of Dynamics

Introduction to Dynamical Systems Basic Concepts of Dynamics Introduction to Dynamical Systems Basic Concepts of Dynamics A dynamical system: Has a notion of state, which contains all the information upon which the dynamical system acts. A simple set of deterministic

More information

Report submitted to Prof. P. Shipman for Math 641, Spring 2012

Report submitted to Prof. P. Shipman for Math 641, Spring 2012 Dynamics at the Horsetooth Volume 4, 2012. The Weierstrass-Enneper Representations Department of Mathematics Colorado State University mylak@rams.colostate.edu drpackar@rams.colostate.edu Report submitted

More information

States of matter. Particles in a gas are widely spread out and can both vibrate and move around freely. They have the most energy of the three states.

States of matter. Particles in a gas are widely spread out and can both vibrate and move around freely. They have the most energy of the three states. States of matter Particles in a solid are closely packed and can vibrate but cannot move around, they have low energies. Particles in a liquid are still closely packed, but can both vibrate and move around

More information

On new chaotic mappings in symbol space

On new chaotic mappings in symbol space Acta Mech Sin DOI 10.1007/s10409-011-0408-1 RESEARCH PAPER On new chaotic mappings in symbol space Inese Bula Jānis Buls Irita Rumbeniece Received: 24 August 2010 / Revised: 15 September 2010 / Accepted:

More information

Introduction to Partial Differential Equations

Introduction to Partial Differential Equations Introduction to Partial Differential Equations Philippe B. Laval KSU Current Semester Philippe B. Laval (KSU) 1D Heat Equation: Derivation Current Semester 1 / 19 Introduction The derivation of the heat

More information

Lectures on Periodic Orbits

Lectures on Periodic Orbits Lectures on Periodic Orbits 11 February 2009 Most of the contents of these notes can found in any typical text on dynamical systems, most notably Strogatz [1994], Perko [2001] and Verhulst [1996]. Complete

More information

Review Sheet 6 Math and Chemistry

Review Sheet 6 Math and Chemistry Review Sheet 6 Math and Chemistry The following are some points of interest in Math and Chemistry. Use this sheet when answering these questions. Molecular Mass- to find the molecular mass, you must add

More information

MATH 312 Section 4.5: Undetermined Coefficients

MATH 312 Section 4.5: Undetermined Coefficients MATH 312 Section 4.5: Undetermined Coefficients Prof. Jonathan Duncan Walla Walla College Spring Quarter, 2007 Outline 1 Differential Operators 2 Annihilators 3 Undetermined Coefficients 4 Conclusion ODEs

More information

MULTI-LEVEL TECHNIQUES FOR THE SOLUTION OF THE KINETIC EQUATIONS IN CONDENSING FLOWS SIMON GLAZENBORG

MULTI-LEVEL TECHNIQUES FOR THE SOLUTION OF THE KINETIC EQUATIONS IN CONDENSING FLOWS SIMON GLAZENBORG MULTI-LEVEL TECHNIQUES FOR THE SOLUTION OF THE KINETIC EQUATIONS IN CONDENSING FLOWS SIMON GLAZENBORG CONTENTS Introduction Theory Test case: Nucleation pulse Conclusions & recommendations 2 WHAT IS CONDENSATION

More information

STOICHIOMETRY via ChemLog

STOICHIOMETRY via ChemLog STOICHIOMETRY via ChemLog 3 Mg + N (g) Mg 3 N by Dr. Stephen Thompson Mr. Joe Staley Ms. Mary Peacock The contents of this module were developed under grant award # P116B-001338 from the Fund for the Improvement

More information

1. (a) The volume of a piece of cake, with radius r, height h and angle θ, is given by the formula: [Yes! It s a piece of cake.]

1. (a) The volume of a piece of cake, with radius r, height h and angle θ, is given by the formula: [Yes! It s a piece of cake.] 1. (a The volume of a piece of cake, with radius r, height h and angle θ, is given by the formula: / 3 V (r,h,θ = 1 r θh. Calculate V r, V h and V θ. [Yes! It s a piece of cake.] V r = 1 r θh = rθh V h

More information

CHEM 116 Collision Theory and Reaction Mechanisms

CHEM 116 Collision Theory and Reaction Mechanisms CHEM 116 Collision Theory and Reaction Mechanisms Lecture 13 Prof. Sevian Note: If there is anything we do not finish about reaction mechanisms today, that is where we will start on Tuesday with Lecture

More information

[A] a [B] b [C] c [D] d

[A] a [B] b [C] c [D] d Unit 4, Lesson #4: Calculations using The quilibrium Constant For the general homogeneous, reversible chemical reaction: aa + bb cc + dd At equilibrium, when the rate of the forward and reverse reactions

More information

Nonlinear Control Lecture 1: Introduction

Nonlinear Control Lecture 1: Introduction Nonlinear Control Lecture 1: Introduction Farzaneh Abdollahi Department of Electrical Engineering Amirkabir University of Technology Fall 2011 Farzaneh Abdollahi Nonlinear Control Lecture 1 1/15 Motivation

More information

ANALYTICAL PROCESS FOR DETERMINATION OF HOPF BIFURCATIONS

ANALYTICAL PROCESS FOR DETERMINATION OF HOPF BIFURCATIONS Volume 1, No. 8, August 2013 Journal of Global Research in Mathematical Archives RESEARCH PAPER Available online at http://www.jgrma.info ANALYTICAL PROCESS FOR DETERMINATION OF HOPF BIFURCATIONS Tarini

More information

B. The rate will increase by a factor of twelve. C. The rate will increase by a factor of twenty-four. D. The rate will decrease by a factor of six.

B. The rate will increase by a factor of twelve. C. The rate will increase by a factor of twenty-four. D. The rate will decrease by a factor of six. 1. If O 2 (g) disappears at a rate of 0.250 M/s at a particular moment in the reaction below, what is the rate of appearance of H 2 O(g) at the same time? C 3 H 8 (g) + 5 O 2 (g) 3 CO 2 (g) + 4 H 2 O(g)

More information

1. Which of the following reactions is the slowest at room temperature?

1. Which of the following reactions is the slowest at room temperature? Kinetics Quiz #2 Factors that Change the Reaction Rate 1. Which of the following reactions is the slowest at room temperature? A. Zn (s) + S (s) ZnS (s) B. Ba 2+ 2- (aq) + SO 4 (aq) BaSO 4(s) C. NH 3(g)

More information

y(n) Time Series Data

y(n) Time Series Data Recurrent SOM with Local Linear Models in Time Series Prediction Timo Koskela, Markus Varsta, Jukka Heikkonen, and Kimmo Kaski Helsinki University of Technology Laboratory of Computational Engineering

More information

UPPER AND LOWER SOLUTIONS FOR A HOMOGENEOUS DIRICHLET PROBLEM WITH NONLINEAR DIFFUSION AND THE PRINCIPLE OF LINEARIZED STABILITY

UPPER AND LOWER SOLUTIONS FOR A HOMOGENEOUS DIRICHLET PROBLEM WITH NONLINEAR DIFFUSION AND THE PRINCIPLE OF LINEARIZED STABILITY ROCKY MOUNTAIN JOURNAL OF MATHEMATICS Volume 30, Number 4, Winter 2000 UPPER AND LOWER SOLUTIONS FOR A HOMOGENEOUS DIRICHLET PROBLEM WITH NONLINEAR DIFFUSION AND THE PRINCIPLE OF LINEARIZED STABILITY ROBERT

More information

Tangent Space and Derivative Mapping on Time Scale

Tangent Space and Derivative Mapping on Time Scale International J.Math. Combin. Vol.2 (2009, 01-10 Tangent Space and Derivative Mapping on Time Scale Emin ÖZYILMAZ (Department of Mathematics,University of Ege, 35100, İzmir, Turkey E-mail: emin.ozyilmaz@ege.edu.tr

More information

Qualitative Analysis of Tumor-Immune ODE System

Qualitative Analysis of Tumor-Immune ODE System of Tumor-Immune ODE System L.G. de Pillis and A.E. Radunskaya August 15, 2002 This work was supported in part by a grant from the W.M. Keck Foundation 0-0 QUALITATIVE ANALYSIS Overview 1. Simplified System

More information

The Volume of a Hypersphere

The Volume of a Hypersphere The hypersphere has the equation The Volume of a Hypersphere x 2 y 2 x 2 w 2 = 2 if centered at the origin (,,,) and has a radius of in four dimensional space. We approach the project of determining its

More information

Solving Differential Equations Using MATLAB

Solving Differential Equations Using MATLAB Solving Differential Equations Using MATLAB Abraham Asfaw aasfaw.student@manhattan.edu November 28, 2011 1 Introduction In this lecture, we will follow up on lecture 2 with a discussion of solutions to

More information

I/O monotone dynamical systems. Germán A. Enciso University of California, Irvine Eduardo Sontag, Rutgers University May 25 rd, 2011

I/O monotone dynamical systems. Germán A. Enciso University of California, Irvine Eduardo Sontag, Rutgers University May 25 rd, 2011 I/O monotone dynamical systems Germán A. Enciso University of California, Irvine Eduardo Sontag, Rutgers University May 25 rd, 2011 BEFORE: Santa Barbara, January 2003 Having handed to me a photocopied

More information

Energy Barriers and Rates - Transition State Theory for Physicists

Energy Barriers and Rates - Transition State Theory for Physicists Energy Barriers and Rates - Transition State Theory for Physicists Daniel C. Elton October 12, 2013 Useful relations 1 cal = 4.184 J 1 kcal mole 1 = 0.0434 ev per particle 1 kj mole 1 = 0.0104 ev per particle

More information

Phase Synchronization of Coupled Rossler Oscillators: Amplitude Effect

Phase Synchronization of Coupled Rossler Oscillators: Amplitude Effect Commun. Theor. Phys. (Beijing, China) 47 (2007) pp. 265 269 c International Academic Publishers Vol. 47, No. 2, February 15, 2007 Phase Synchronization of Coupled Rossler Oscillators: Amplitude Effect

More information

Position, Speed and Velocity Position is a variable that gives your location relative to an origin. The origin is the place where position equals 0.

Position, Speed and Velocity Position is a variable that gives your location relative to an origin. The origin is the place where position equals 0. Position, Speed and Velocity Position is a variable that gives your location relative to an origin. The origin is the place where position equals 0. The position of this car at 50 cm describes where the

More information

Waves in a Shock Tube

Waves in a Shock Tube Waves in a Shock Tube Ivan Christov c February 5, 005 Abstract. This paper discusses linear-wave solutions and simple-wave solutions to the Navier Stokes equations for an inviscid and compressible fluid

More information

12. Lines in R 3. through P 0 and is parallel to v is written parametrically as a function of t: Using vector notation, the same line is written

12. Lines in R 3. through P 0 and is parallel to v is written parametrically as a function of t: Using vector notation, the same line is written 12. Lines in R 3 Given a point P 0 = (x 0, y 0, z 0 ) and a direction vector v 1 = a, b, c in R 3, a line L that passes through P 0 and is parallel to v is written parametrically as a function of t: x(t)

More information

Phase Response Curves, Delays and Synchronization in Matlab

Phase Response Curves, Delays and Synchronization in Matlab Phase Response Curves, Delays and Synchronization in Matlab W. Govaerts and B. Sautois Department of Applied Mathematics and Computer Science, Ghent University, Krijgslaan 281-S9, B-9000 Ghent, Belgium

More information

Notes for Expansions/Series and Differential Equations

Notes for Expansions/Series and Differential Equations Notes for Expansions/Series and Differential Equations In the last discussion, we considered perturbation methods for constructing solutions/roots of algebraic equations. Three types of problems were illustrated

More information

Lecture 22: The Arrhenius Equation and reaction mechanisms. As we wrap up kinetics we will:

Lecture 22: The Arrhenius Equation and reaction mechanisms. As we wrap up kinetics we will: As we wrap up kinetics we will: Lecture 22: The Arrhenius Equation and reaction mechanisms. Briefly summarize the differential and integrated rate law equations for 0, 1 and 2 order reaction Learn how

More information

7 Pendulum. Part II: More complicated situations

7 Pendulum. Part II: More complicated situations MATH 35, by T. Lakoba, University of Vermont 60 7 Pendulum. Part II: More complicated situations In this Lecture, we will pursue two main goals. First, we will take a glimpse at a method of Classical Mechanics

More information

Lecture 3. The Kinetic Molecular Theory of Gases

Lecture 3. The Kinetic Molecular Theory of Gases Lecture 3. The Kinetic Molecular Theory of Gases THE IDEAL GAS LAW: A purely empirical law solely the consequence of experimental observations Explains the behavior of gases over a limited range of conditions

More information

Factors that Affect Reaction Rates

Factors that Affect Reaction Rates Factors that Affect Reaction Rates Preface: There are 2 kinds of reactions: Homogeneous reactions - all reactants are in the same phase (don't consider products) eg.) 3H 2(g) + N 2(g) 2NH 3(g) Ag + (aq)

More information

Ordinal Analysis of Time Series

Ordinal Analysis of Time Series Ordinal Analysis of Time Series K. Keller, M. Sinn Mathematical Institute, Wallstraße 0, 2552 Lübeck Abstract In order to develop fast and robust methods for extracting qualitative information from non-linear

More information

2.3 Damping, phases and all that

2.3 Damping, phases and all that 2.3. DAMPING, PHASES AND ALL THAT 107 2.3 Damping, phases and all that If we imagine taking our idealized mass on a spring and dunking it in water or, more dramatically, in molasses), then there will be

More information

Numerical Methods for Partial Differential Equations. CAAM 452 Spring 2005 Lecture 2 Instructor: Tim Warburton

Numerical Methods for Partial Differential Equations. CAAM 452 Spring 2005 Lecture 2 Instructor: Tim Warburton Numerical Methods for Partial Differential Equations CAAM 452 Spring 2005 Lecture 2 Instructor: Tim Warburton Note on textbook for finite difference methods Due to the difficulty some students have experienced

More information

The Conley Index and Rigorous Numerics of Attracting Periodic Orbits

The Conley Index and Rigorous Numerics of Attracting Periodic Orbits The Conley Index and Rigorous Numerics of Attracting Periodic Orbits Marian Mrozek Pawe l Pilarczyk Conference on Variational and Topological Methods in the Study of Nonlinear Phenomena (Pisa, 2000) 1

More information

Molecular Physics. Attraction between the ions causes the chemical bond.

Molecular Physics. Attraction between the ions causes the chemical bond. Molecular Physics A molecule is a stable configuration of electron(s) and more than one nucleus. Two types of bonds: covalent and ionic (two extremes of same process) Covalent Bond Electron is in a molecular

More information

Laplace s Method for Ordinary Differential Equations

Laplace s Method for Ordinary Differential Equations Physics 2400 Spring 2016 Laplace s Method for Ordinary Differential Equations Lecture notes by M. G. Rozman Last modified: March 31, 2016 We can [1, pp. 124-8], [2, h. VIII], [3, Apps. a-b], [4, h. 18],

More information

Research Article Hamilton-Poisson Realizations for the Lü System

Research Article Hamilton-Poisson Realizations for the Lü System Mathematical Problems in Engineering Volume 011, Article ID 8435, 13 pages doi:10.1155/011/8435 Research Article Hamilton-Poisson Realizations for the Lü System Camelia Pop, 1 Camelia Petrişor, 1 and Dumitru

More information

SAMPLE CHAPTERS UNESCO EOLSS WAVES IN THE OCEANS. Wolfgang Fennel Institut für Ostseeforschung Warnemünde (IOW) an der Universität Rostock,Germany

SAMPLE CHAPTERS UNESCO EOLSS WAVES IN THE OCEANS. Wolfgang Fennel Institut für Ostseeforschung Warnemünde (IOW) an der Universität Rostock,Germany WAVES IN THE OCEANS Wolfgang Fennel Institut für Ostseeforschung Warnemünde (IOW) an der Universität Rostock,Germany Keywords: Wind waves, dispersion, internal waves, inertial oscillations, inertial waves,

More information

A Measure Theory Tutorial (Measure Theory for Dummies) Maya R. Gupta

A Measure Theory Tutorial (Measure Theory for Dummies) Maya R. Gupta A Measure Theory Tutorial (Measure Theory for Dummies) Maya R. Gupta {gupta}@ee.washington.edu Dept of EE, University of Washington Seattle WA, 98195-2500 UW Electrical Engineering UWEE Technical Report

More information

Period-doubling cascades of a Silnikov equation

Period-doubling cascades of a Silnikov equation Period-doubling cascades of a Silnikov equation Keying Guan and Beiye Feng Science College, Beijing Jiaotong University, Email: keying.guan@gmail.com Institute of Applied Mathematics, Academia Sinica,

More information

The first order quasi-linear PDEs

The first order quasi-linear PDEs Chapter 2 The first order quasi-linear PDEs The first order quasi-linear PDEs have the following general form: F (x, u, Du) = 0, (2.1) where x = (x 1, x 2,, x 3 ) R n, u = u(x), Du is the gradient of u.

More information

The Lost Cousin of the Fundamental Theorem of Algebra

The Lost Cousin of the Fundamental Theorem of Algebra 290 MATHEMATICS MAGAZINE The Lost Cousin of the Fundamental Theorem of Algebra TIMO TOSSAVAINEN Department of Teacher Education Joensuu University Savonlinna, Finland timo.tossavainen@joensuu.fi Have you

More information

which has an equilibrium constant of Which of the following diagrams represents a mixture of the reaction at equilibrium?

which has an equilibrium constant of Which of the following diagrams represents a mixture of the reaction at equilibrium? Chapter 9 Quiz: Chemical Equilibria 1. Which of the following statements is true regarding chemical equilibrium? I. The concentrations of reactants and products at equilibrium are constant, which means

More information

Double replacement reactions

Double replacement reactions 1. Learn to predict Double replacement reaction If and when a reaction occurs, what are the products? 2. Learn to write Double replacement reaction: (i) Balanced chemical reaction (ii) Net ionicreaction

More information

PH36010: Numerical Methods - Evaluating the Lorenz Attractor using Runge-Kutta methods Abstract

PH36010: Numerical Methods - Evaluating the Lorenz Attractor using Runge-Kutta methods Abstract PH36010: Numerical Methods - Evaluating the Lorenz Attractor using Runge-Kutta methods Mr. Benjamen P. Reed (110108461) IMPACS, Aberystwyth University January 31, 2014 Abstract A set of three coupled ordinary

More information

Transient Heat Conduction in a Circular Cylinder

Transient Heat Conduction in a Circular Cylinder Transient Heat Conduction in a Circular Cylinder The purely radial 2-D heat equation will be solved in cylindrical coordinates using variation of parameters. Assuming radial symmetry the solution is represented

More information

6.2 Brief review of fundamental concepts about chaotic systems

6.2 Brief review of fundamental concepts about chaotic systems 6.2 Brief review of fundamental concepts about chaotic systems Lorenz (1963) introduced a 3-variable model that is a prototypical example of chaos theory. These equations were derived as a simplification

More information

Math 304, Differential Forms, Lecture 1 May 26, 2011

Math 304, Differential Forms, Lecture 1 May 26, 2011 Math 304, Differential Forms, Lecture 1 May 26, 2011 Differential forms were developed and applied to vector fields in a decisive way by the French Mathematician Élie artan (1869 1951) in the first quarter

More information

Proceedings of Neural, Parallel, and Scientific Computations 4 (2010) xx-xx PHASE OSCILLATOR NETWORK WITH PIECEWISE-LINEAR DYNAMICS

Proceedings of Neural, Parallel, and Scientific Computations 4 (2010) xx-xx PHASE OSCILLATOR NETWORK WITH PIECEWISE-LINEAR DYNAMICS Proceedings of Neural, Parallel, and Scientific Computations 4 (2010) xx-xx PHASE OSCILLATOR NETWORK WITH PIECEWISE-LINEAR DYNAMICS WALTER GALL, YING ZHOU, AND JOSEPH SALISBURY Department of Mathematics

More information

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

Chaotic and non-chaotic mixed oscillations in a logistic systems with delay Chaotic and non-chaotic mixed oscillations in a logistic systems with delay Marek Berezowski a, b, Artur Grabski b a Institute of Chemical Engineering and Physical Chemistry, Cracow University of Technology,

More information

Modelling biological oscillations

Modelling biological oscillations Modelling biological oscillations Shan He School for Computational Science University of Birmingham Module 06-23836: Computational Modelling with MATLAB Outline Outline of Topics Van der Pol equation Van

More information

Complex Numbers: A Brief Introduction. By: Neal Dempsey. History of Mathematics. Prof. Jennifer McCarthy. July 16, 2010

Complex Numbers: A Brief Introduction. By: Neal Dempsey. History of Mathematics. Prof. Jennifer McCarthy. July 16, 2010 1 Complex Numbers: A Brief Introduction. By: Neal Dempsey History of Mathematics Prof. Jennifer McCarthy July 16, 2010 2 Abstract Complex numbers, although confusing at times, are one of the most elegant

More information

SELF-ORGANIZATION IN NONRECURRENT COMPLEX SYSTEMS

SELF-ORGANIZATION IN NONRECURRENT COMPLEX SYSTEMS Letters International Journal of Bifurcation and Chaos, Vol. 10, No. 5 (2000) 1115 1125 c World Scientific Publishing Company SELF-ORGANIZATION IN NONRECURRENT COMPLEX SYSTEMS PAOLO ARENA, RICCARDO CAPONETTO,

More information

Introducing the Normal Distribution

Introducing the Normal Distribution Department of Mathematics Ma 3/103 KC Border Introduction to Probability and Statistics Winter 2017 Lecture 10: Introducing the Normal Distribution Relevant textbook passages: Pitman [5]: Sections 1.2,

More information

Opleiding Wiskunde & Informatica

Opleiding Wiskunde & Informatica Opleiding Wiskunde & Informatica Numerics and continuation for Reaction-Diffusion equations Renzo Baasdam s1524054 Supervisors: Martina Chirilus-Bruckner (MI) & Michael Emmerich (LIACS) BACHELOR THESIS

More information

Introduction to Vectors

Introduction to Vectors Introduction to Vectors K. Behrend January 31, 008 Abstract An introduction to vectors in R and R 3. Lines and planes in R 3. Linear dependence. 1 Contents Introduction 3 1 Vectors 4 1.1 Plane vectors...............................

More information

A CAUCHY PROBLEM OF SINE-GORDON EQUATIONS WITH NON-HOMOGENEOUS DIRICHLET BOUNDARY CONDITIONS 1. INTRODUCTION

A CAUCHY PROBLEM OF SINE-GORDON EQUATIONS WITH NON-HOMOGENEOUS DIRICHLET BOUNDARY CONDITIONS 1. INTRODUCTION Trends in Mathematics Information Center for Mathematical Sciences Volume 4, Number 2, December 21, Pages 39 44 A CAUCHY PROBLEM OF SINE-GORDON EQUATIONS WITH NON-HOMOGENEOUS DIRICHLET BOUNDARY CONDITIONS

More information

CALCULATION OF NONLINEAR VIBRATIONS OF PIECEWISE-LINEAR SYSTEMS USING THE SHOOTING METHOD

CALCULATION OF NONLINEAR VIBRATIONS OF PIECEWISE-LINEAR SYSTEMS USING THE SHOOTING METHOD Vietnam Journal of Mechanics, VAST, Vol. 34, No. 3 (2012), pp. 157 167 CALCULATION OF NONLINEAR VIBRATIONS OF PIECEWISE-LINEAR SYSTEMS USING THE SHOOTING METHOD Nguyen Van Khang, Hoang Manh Cuong, Nguyen

More information

USING COUPLED OSCILLATORS TO MODEL THE SINO-ATRIAL NODE IN THE HEART

USING COUPLED OSCILLATORS TO MODEL THE SINO-ATRIAL NODE IN THE HEART USING COUPLED OSCILLATORS TO MODEL THE SINO-ATRIAL NODE IN THE HEART A Thesis Presented to the Faculty of the Graduate School of Cornell University in Partial Fulfillment of the Requirements for the Degree

More information

ELEMENTS, COMPOUNDS AND MIXTURES AND HOW THEY ARE REPRESENTED. Jan 12-13, 2014

ELEMENTS, COMPOUNDS AND MIXTURES AND HOW THEY ARE REPRESENTED. Jan 12-13, 2014 ELEMENTS, COMPOUNDS AND MIXTURES AND HOW THEY ARE REPRESENTED Jan 12-13, 2014 WHAT ARE ELEMENTS? Elements are pure substances Made of only one kind of material Has definite properties, and Is the same

More information

St Robert of Newminster Catholic School and Sixth Form College

St Robert of Newminster Catholic School and Sixth Form College St Robert of Newminster Catholic School and Sixth Form College Year 12 Pre-Course Tasks: CHEMISTRY Exercise Mark Grade Atomic structure Chemical bonding Chemical equations Maths for chemists Moles Name:

More information

The Spike Response Model: A Framework to Predict Neuronal Spike Trains

The Spike Response Model: A Framework to Predict Neuronal Spike Trains The Spike Response Model: A Framework to Predict Neuronal Spike Trains Renaud Jolivet, Timothy J. Lewis 2, and Wulfram Gerstner Laboratory of Computational Neuroscience, Swiss Federal Institute of Technology

More information

Section 6.4 DEs with Discontinuous Forcing Functions

Section 6.4 DEs with Discontinuous Forcing Functions Section 6.4 DEs with Discontinuous Forcing Functions Key terms/ideas: Discontinuous forcing function in nd order linear IVPs Application of Laplace transforms Comparison to viewing the problem s solution

More information

Polymer Dynamics and Rheology

Polymer Dynamics and Rheology Polymer Dynamics and Rheology 1 Polymer Dynamics and Rheology Brownian motion Harmonic Oscillator Damped harmonic oscillator Elastic dumbbell model Boltzmann superposition principle Rubber elasticity and

More information

1 Random walks and data

1 Random walks and data Inference, Models and Simulation for Complex Systems CSCI 7-1 Lecture 7 15 September 11 Prof. Aaron Clauset 1 Random walks and data Supposeyou have some time-series data x 1,x,x 3,...,x T and you want

More information

March 9, :18 Int J. Bifurcation and Chaos/INSTRUCTION FILE Morfu2v2 EFFECT OF NOISE AND STRUCTURAL INHOMOGENEITIES IN REACTION DIFFUSION MEDIA

March 9, :18 Int J. Bifurcation and Chaos/INSTRUCTION FILE Morfu2v2 EFFECT OF NOISE AND STRUCTURAL INHOMOGENEITIES IN REACTION DIFFUSION MEDIA March 9, 2007 10:18 Int J. Bifurcation and Chaos/INSTRUCTION FILE Int J. Bifurcation and Chaos Submission Style EFFECT OF NOISE AND STRUCTURAL INHOMOGENEITIES IN REACTION DIFFUSION MEDIA S. Morfu Laboratoire

More information

Chapter Test B. Chapter: Chemical Equilibrium. following equilibrium system? 2CO(g) O 2 (g) ^ 2CO 2 (g)

Chapter Test B. Chapter: Chemical Equilibrium. following equilibrium system? 2CO(g) O 2 (g) ^ 2CO 2 (g) Assessment Chapter Test B Chapter: Chemical Equilibrium PART I In the space provided, write the letter of the term or phrase that best completes each statement or best answers each question. 1. What is

More information

First and Second Order Differential Equations Lecture 4

First and Second Order Differential Equations Lecture 4 First and Second Order Differential Equations Lecture 4 Dibyajyoti Deb 4.1. Outline of Lecture The Existence and the Uniqueness Theorem Homogeneous Equations with Constant Coefficients 4.2. The Existence

More information

3 Stability and Lyapunov Functions

3 Stability and Lyapunov Functions CDS140a Nonlinear Systems: Local Theory 02/01/2011 3 Stability and Lyapunov Functions 3.1 Lyapunov Stability Denition: An equilibrium point x 0 of (1) is stable if for all ɛ > 0, there exists a δ > 0 such

More information

Linear Second Order ODEs

Linear Second Order ODEs Chapter 3 Linear Second Order ODEs In this chapter we study ODEs of the form (3.1) y + p(t)y + q(t)y = f(t), where p, q, and f are given functions. Since there are two derivatives, we might expect that

More information

Work sheet / Things to know. Chapter 3

Work sheet / Things to know. Chapter 3 MATH 251 Work sheet / Things to know 1. Second order linear differential equation Standard form: Chapter 3 What makes it homogeneous? We will, for the most part, work with equations with constant coefficients

More information

Stability of the Lagrange Points, L 4 and L 5

Stability of the Lagrange Points, L 4 and L 5 Stability of the Lagrange Points, L 4 and L 5 Thomas Greenspan January 7, 014 Abstract A proof of the stability of the non collinear Lagrange Points, L 4 and L 5. We will start by covering the basics of

More information

Equilibrium and LeChatelier s Principle

Equilibrium and LeChatelier s Principle 1 Equilibrium and LeChatelier s Principle Purpose: To examine LeChatelier s Principle by studying disturbances applied to several equilibrium systems. Introduction Many chemical reactions reach a state

More information

Nick Pritchard 1 School of Physics, Georgia Institute of Technology, Atlanta, Georgia 30332, USA

Nick Pritchard 1 School of Physics, Georgia Institute of Technology, Atlanta, Georgia 30332, USA An Investigation of the Chaotic Dripping Faucet Nick Pritchard 1 School of Physics, Georgia Institute of Technology, Atlanta, Georgia 30332, USA (Dated: 16 December 2011) The dripping faucet experiment

More information

Mathematics for Chemists 2 Lecture 14: Fourier analysis. Fourier series, Fourier transform, DFT/FFT

Mathematics for Chemists 2 Lecture 14: Fourier analysis. Fourier series, Fourier transform, DFT/FFT Mathematics for Chemists 2 Lecture 14: Fourier analysis Fourier series, Fourier transform, DFT/FFT Johannes Kepler University Summer semester 2012 Lecturer: David Sevilla Fourier analysis 1/25 Remembering

More information

MATH 1910 Limits Numerically and Graphically Introduction to Limits does not exist DNE DOES does not Finding Limits Numerically

MATH 1910 Limits Numerically and Graphically Introduction to Limits does not exist DNE DOES does not Finding Limits Numerically MATH 90 - Limits Numerically and Graphically Introduction to Limits The concept of a limit is our doorway to calculus. This lecture will explain what the limit of a function is and how we can find such

More information

Take home Exam due Wednesday, Aug 26. In class Exam will be the that morning class multiple choice questions.

Take home Exam due Wednesday, Aug 26. In class Exam will be the that morning class multiple choice questions. Announcements Take home Exam due Wednesday, Aug 26. In class Exam will be the that morning class. 15-20 multiple choice questions. Updated projects Aug 28: answer what lab chemistry needs to get done to

More information

Study Guide for Module 13 An Introduction to Equilibrium

Study Guide for Module 13 An Introduction to Equilibrium Chemistry 1020, Module 13 Name Study Guide for Module 13 An Introduction to Equilibrium Reading Assignment: Section 12.1 and Chapter 13 of Chemistry, 6th Edition by Zumdahl. Guide for Your Lecturer: 1.

More information

The Nonlinear Pendulum

The Nonlinear Pendulum The Nonlinear Pendulum - Pádraig Ó Conbhuí - 08531749 TP Monday 1. Abstract This experiment was performed to examine the effects that linearizing equations has on the accuracy of results and to find ways

More information

Solutions: Section 2.5

Solutions: Section 2.5 Solutions: Section 2.5. Problem : Given dy = ay + by2 = y(a + by) with a, b > 0. For the more general case, we will let y 0 be any real number. Always look for the equilibria first! In this case, y(a +

More information