PH36010: Numerical Methods - Evaluating the Lorenz Attractor using Runge-Kutta methods Abstract
|
|
- Stella Elliott
- 6 years ago
- Views:
Transcription
1 PH36010: Numerical Methods - Evaluating the Lorenz Attractor using Runge-Kutta methods Mr. Benjamen P. Reed ( ) IMPACS, Aberystwyth University January 31, 2014 Abstract A set of three coupled ordinary differential equations known as the Lorenz equations were evaluated using the Fourth-Order Runge-Kutta method to produce a solution known as the Lorenz Attractor. The parameters of the Lorenz attractor were systematically altered using a FORTRAN program to ascertain their effect on the behaviour of the chaotic system and the possible physical consequences of these changes was discussed. The functionality of the Runge- Kutta method is also considered.
2 Table of Contents 1. Introduction 3 2. Background Theory and Procedure Lorenz Equations Fourth-Order Runge-Kutta Evaluating the Lorenz Equations 5 3. Results Typical Lorenz Attractor Dependence on sigma σ Dependence on beta β Dependence on rho ρ Discussion Conclusion 14 Acknowledgements 14 References 14 Appendixes 15 A. Runge-Kutta Code to Evaluate Lorenz Equations 15 2
3 1. Introduction Of all the dynamic systems that physics, and the physical sciences describe and model, none are more subtle yet potent than chaotic systems; these are systems whose behaviour is highly sensitive to the initial conditions of that system. One example of a chaotic system is the double rod pendulum, which has a full rotation hinge in the middle of the system, allowing the bottom half of the pendulum to swing quasi-independently of the main pendulum. The high initial condition dependence in these systems has been poetically named the Butterfly effect, after Edward Norton Lorenz s (an American mathematician and meteorologist) work on the Lorenz attractor and his 1972 talk, Predictability: Does the Flap of a Butterfly's Wings in Brazil set off a Tornado in Texas? As a meteorologist, Lorenz investigated the functionality and possibility of performing accurate long-term weather forecasts, and his research demonstrated that even the smallest change in a weather system could have significant effects on the system as it evolved in time. He concluded that any model of any weather system would have to take into account the smallest variations in the weather conditions in order to have a high accuracy; today s weather modelling computers are designed to compute thousands of variables in order to predict the short-term forecast. However, despite all the known variables, there are still just as many unknown variables that have non-periodic or chaotic behaviour, and as such they are intrinsically unpredictable in the long-term. Even smaller systems with considerably less variables, still exhibit this long-term unpredictability (Palmer, 1993). One simple system to exhibit chaotic behaviour is the Lorentz equations, a system of three coupled ordinary differential equations that are a simple and idealistic set of twelve hydrodynamical equations that describe turbulent flow in the atmosphere. When the 3- dimensional coordinates of this system are plotted with time as the independent variable, an interesting graph is observed. For a single point moving around in phase space (i.e. x(t), y(t) and z(t) where t is time), it traces out a curve, or an orbit which always converges and then remains on a surface comprised of two wings. For any initial conditions, the point always and eventually links up with this two-winged surface, and hence it is known generally as an attractor. The attractor created with the Lorenz equations, shown in figure 1, is appropriately named the Butterfly attractor for its appearance and also because it visualises solutions that exhibit the butterfly effect (Taylor, 2011). Figure 1 - A visualisation of the Lorenz attractor in phase space; solution to the Lorenz equations as an orbit in phase space. (Source: As previously mentioned, the Lorenz attractor can be constructed using three coupled ordinary differential equations, and as such can be modelled with a relatively simple program and a robust graphing client. In this report, the Lorenz attractor for an arbitrary chaotic system is evaluated using a program written in the script-based programming language FORTRAN. The program employs the use of the Fourth-Order Runge-Kutta method in order to solve the Lorenz equation and thus produce useable data. The parameters and variables of the system were changed systematically in order to ascertain their effect on the shape and profile of the butterfly attractor in phase space. This report will also discuss the effectiveness of using the Fourth-Order Runge- Kutta method with coupled ordinary differential equations. 3
4 2. Background Theory and Procedure 2.1 The Lorenz Equations The Lorenz equations, as written by Edward Lorenz himself are dx dt = σ (y x) (1) dy dt = x(ρ z) y (2) dz dt = xy βz (3) where x, y, and z are variables in phase space and σ, ρ. And β are the constant parameters of the system. It is easy to see in this system, that as one variable of phase space changes, the other two variables in phase space change in response. The parameters in these equations each have physical meanings are not just some mathematical quirk. Sigma σ is called the Prandtl number and is a description of the systems physical characteristics such as viscosity and thermal conductivity in the case of convection flow. Rho ρ represents a control parameter, which is the difference between the extrema driving forces in the system. For the example where the Lorenz attractor is characterising convection in a fluid filled tank, ρ would represent the temperature difference between the top and bottom of the tank. Finally beta β is a value that describes the shape of the box that the Lorenz attractor is contained within. For the convection example, β would be the width-to-height ratio of the convection layer (Stockie & Wong, 2009). 2.2 The Fourth Order Runge-Kutta The Runge-Kutta method is an iterative numerical method for solving systems of coupled ordinary differential equations (CODEs), and is highly considered to be the most accurate numerical method for this purpose for the amount of effort required to code it. It is similar to Euler s method of solving ODEs however whereas Euler s method only evaluates each iterative step at one point, Runge-Kutta evaluates each step at four separate locations as shown by equations 4 to 7. Two evaluations are made at the end points of each interval, and two are made in the middle. F 1 = δtf (y i,t i ) (4)! F 2 = δtf # y i F 1,t i + 1 " 2 δt $ %! F 3 = δtf # y i F,t i " 2 δt $ % & (5) & (6) F 4 = dtf ( y i + F 3,t i +δt) (7) It then combines these solutions into one highly accurate approximation of the function at that point. In essence, this means that any number of time-dependant CODEs can be evaluated by one 4
5 program and hence data can be formulated for data analysis. The final recombination step is characterised by equation 8, which gives the phase space coordinate at that time step (Süli, 2013). y i+1 = y i + 1 ( 6 F + 2F + 2F + F ) (8) (N.B. y in equations 4 8 can represent any value in any axis, i.e. y can be x or z. Cartesian coordinates are evaluated separately for each axis, so Runge-Kutta is executed three times for each point in phase space). 2.3 Evaluating the Lorenz Equations To solve the Lorenz equations and thus produce the Lorenz attractor plot, a program was written in FORTRAN, which used the aforementioned Fourth-Order Runge-Kutta method to evaluate the CODEs hence produce useable data in the form of a comma separating variable file. This data was then sent to gnuplot to produce 3-dimensional phase space plots and 2-dimension time dependent graphs of the variables (Press et al., 2007). The full FORTRAN can be viewed in the appendix. To ascertain the effect of parameters on the Lorenz attractor, a default attractor was plotted and then each parameter was adjusted systematically. The attractor was then re-plotted and compared with the original to ascertain what had changed. 5
6 3. Results 3.1 Typical Lorenz Attractor Figure 2 - A 3D plot of a typical Lorenz attractor in phase space, where σ = 10, β = 8/3, ρ = 28. Figure 3 - A 2D plot of Lorenz attractor phase space coordinates against time, where σ = 10, β = 8/3, ρ = 28. 6
7 3.2 Dependence on sigma σ Figure 4 - A 3D plot of a typical Lorenz attractor in phase space, where σ = 5, β = 8/3, ρ = 28. Figure 5 - A 2D plot of Lorenz attractor phase space coordinates against time, where σ = 5, β = 8/3, ρ = 28. 7
8 Figure 6 - A 3D plot of a typical Lorenz attractor in phase space, where σ = 20, β = 8/3, ρ = 28. Figure 7 - A 2D plot of Lorenz attractor phase space coordinates against time, where σ = 20, β = 8/3, ρ = 28. 8
9 3.3 Dependence on beta β Figure 8 - A 3D plot of a typical Lorenz attractor in phase space, where σ = 10, β =4/3, ρ = 28. Figure 9 - A 2D plot of Lorenz attractor phase space coordinates against time, where σ = 10, β =4/3, ρ = 28. 9
10 Figure 10 - A 3D plot of a typical Lorenz attractor in phase space, where σ = 10, β =16/3, ρ = 28. Figure 11 - A 2D plot of Lorenz attractor phase space coordinates against time, where σ = 10, β =16/3, ρ =
11 3.4 Dependence on rho ρ Figure 12 - A 3D plot of a typical Lorenz attractor in phase space, where σ = 10, β =8/3, ρ = 14. Figure 13 - A 2D plot of Lorenz attractor phase space coordinates against time, where σ = 10, β =8/3, ρ =
12 Figure 14 - A 3D plot of a typical Lorenz attractor in phase space, where σ = 10, β =8/3, ρ = 56. Figure 15 - A 2D plot of Lorenz attractor phase space coordinates against time, where σ = 10, β =8/3, ρ =
13 4. Discussion Figure 2 shows a typical Lorenz attractor with the parameters σ = 10, β = 8/3, ρ = 28. This 3- dimensional plot, along with its 2-dimensional temporal dependence graph in figure 3 was treated as the control, a reference that could be used to determine the changes in the Lorenz attractor shape as the parameters were altered. The 3d representation in figure 2 was typical of the butterfly wings found in the literature and resource material. Hence, it was assumed that these figures were correct and an accurate demonstration of chaotic flow in atmospheric systems. To infer some physical meaning to the results detailed in section 3, it has been assumed that the Lorenz equations in this situation are characterising the convectional flow of a viscous fluid in a heated tank. The first parameter that was altered was sigma σ, which as mentioned in section 2.1, is the Prandtl number and is a description of the systems physical characteristics of flow. It tells information about the fluid s viscosity and thermal conductivity. Figures 4 and 5 show the situation where the parameter sigma has been halved from it original value. Physically, this means that the viscosity of the fluid is higher, so the fluid takes longer to converge on the attractor surface. Furthermore, once on the surface the flow begins to spiral inwards on only one wing. Figure 5 shows that the movement of one unit of fluid decreases over time, which implies that the flow is so viscous that it actually dampens any movement within it. Conversely, in figures 6 and 7 the viscosity is lower and hence the flow quickly converges on the attractor surface, and onto two wings. Once there, it is able to sweep out large arcs of movement in the decreased viscosity hence why there appears to be no inner orbits within the wings. The second parameter to be altered was beta β, which characterises the width-to-height ratio of the convection layer. Intuition infers that a smaller width in the convection zone means that flow will frequently switch between one wing of the attractor surface and the other. Whereas a larger width in the convection zone will mean that flow will switch attractor wings frequently, if at all. The figures in section 3.3 seem to agree with this intuition, with the 2D temporal dependence graphs demonstrating this much more clearly. Figure 9 shows the temporal dependence of the flow when the width-to-height ratio is low; the x(t) and y(t) frequently change phase which implies that the flow is moving through the convection zone frequently and changing which wing it orbits. On the other hand, figure 10 where β is double the control parameter, shows that once the flow is orbiting one wing, it never switches over to other wing, implying that the width of the convection zone is so large, that flow cannot or rarely passes into it. Finally, the last parameter rho ρ details the temperature difference between the top and the bottom of the tank. It is expected that a large difference in temperature between the top and the bottom of the tank will give rise to a large flow rate with frequent changes in attractor orbit. The system would also take an extended period of time to reach equilibrium and so the flow rate would continue to be rapid for a long time. With a small temperature difference, the thermal driving forces would be reduced and the system would reach equilibrium relatively faster. This is certainly what seems to be happening in figures 12 and 13. The low thermodynamic driving force means that the flow hardly changes orbit, or in this case, doesn t at all. Furthermore, once equilibrium is reached, all macroscopic movement stops as shown by figure 13 and the system tends toward equilibrium, a dampening effect is observed much like in figure 5. From this discussion, it is apparent that small changes in the Lorentz attractor system give rise to large macroscopic effects. The Runge-Kutta method has illustrated these changes well and is a very suitable method to use in the modelling of chaotic systems. It does not appear to show signs of divergence, a common issue in the simpler Euler s method, and so confidence is high that there is minimal error in the Lorentz modelling. 13
14 5. Conclusion A model of the Lorenz attractor has been constructed using the Runge-Kutta method to solve the Lorenz equations, a set of three coupled ordinary differential equations that describe chaotic flow in atmospheric or hydrodynamic systems. The modelling proved to be successful and the Lorenz equation parameters were varied to ascertain their effect on the behaviour of the system. It was found that small changes to these parameters could have huge effects on the dynamics of the system, either making it more chaotic, or taming it such that the flow rate of the system tends to zero. In the cases where the flow rate tended to zero, the parameters gave rise to an unstable system where the dynamics could not be maintain for an extended period of time. Conversely, in the situations where the flow became more chaotic, the system remained stable and the orbits of the attractor surface were more pronounced. In respect to the numerical method employed, the Runge-Kutta method performed well and produced an accurate, stable model in which to observe the Lorenz attractor. Acknowledgements The author would like to extend their gratitude to Dr. Martin Wilding for the tuition he provided during the course, and for taking time out of his own schedule to give students on the course the time necessary to provide work of their best effort. The author would also like to thank Mr. Patrick Dixon, Mr. Andrew Leonard and Mr. Tom Knight for helping to run the module and provide lessons in effective coding. Finally, the author would like to thank Miss. Carley Martin, Miss. Fffion Whitehouse and Miss. Rose Cooper for providing much needed discussion and support during the module. References Palmer, T., Exploring Chaos: A Guide to the New Science of Disorder. London: W. W. Norton & Co. Press, W.H., Teukolsky, S.A., Vetterling, W.T. & Flannery, B.P Numerical Recipes: The Art of Scientific Computing. 3rd ed. Cambridge: Cambridge Press. Stockie, J. & Wong, K., A Simple Model of the Unpredictability of Weather: The Lorenz Equations. Clouds, Aerosols and Climate Feedback Publications. Available at: < [Accessed 30 January 2014] Süli, E Numerical Solution of Ordinary Differential Equations. [pdf] Mathematical Institute, University of Oxford. Available at: < [Accessed 30 January 2014] Taylor, R.L.V., Attractors: Nonstrange to Chaotic. SIAM Undergraduate Research Online, 4. 14
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 informationBy Nadha CHAOS THEORY
By Nadha CHAOS THEORY What is Chaos Theory? It is a field of study within applied mathematics It studies the behavior of dynamical systems that are highly sensitive to initial conditions It deals with
More informationIntroduction. Prediction MATH February 2017
21 February 2017 Predicting the future is very difficult, especially if it s about the future. Niels Bohr Can we say what is going to happen: in the next minute? tomorrow? next year? Predicting the future
More informationMaps and differential equations
Maps and differential equations Marc R. Roussel November 8, 2005 Maps are algebraic rules for computing the next state of dynamical systems in discrete time. Differential equations and maps have a number
More informationAPPM 2460 CHAOTIC DYNAMICS
APPM 2460 CHAOTIC DYNAMICS 1. Introduction Today we ll continue our exploration of dynamical systems, focusing in particular upon systems who exhibit a type of strange behavior known as chaos. We will
More informationPhysics 584 Computational Methods
Physics 584 Computational Methods Introduction to Matlab and Numerical Solutions to Ordinary Differential Equations Ryan Ogliore April 18 th, 2016 Lecture Outline Introduction to Matlab Numerical Solutions
More informationEnvironmental Physics Computer Lab #8: The Butterfly Effect
Environmental Physics Computer Lab #8: The Butterfly Effect The butterfly effect is the name coined by Lorenz for the main feature of deterministic chaos: sensitivity to initial conditions. Lorenz has
More informationComparison of Numerical Ordinary Differential Equation Solvers
Adrienne Criss Due: October, 008 Comparison of Numerical Ordinary Differential Equation Solvers Many interesting physical systems can be modeled by ordinary differential equations (ODEs). Since it is always
More informationESTIMATING THE ATTRACTOR DIMENSION OF THE EQUATORIAL WEATHER SYSTEM M. Leok B.T.
This paper was awarded in the I International Competition (99/9) First Step to Nobel Prize in Physics and published in the competition proceedings (Acta Phys. Pol. A 8 Supplement, S- (99)). The paper is
More informationChap. 20: Initial-Value Problems
Chap. 20: Initial-Value Problems Ordinary Differential Equations Goal: to solve differential equations of the form: dy dt f t, y The methods in this chapter are all one-step methods and have the general
More informationTHREE DIMENSIONAL SYSTEMS. Lecture 6: The Lorenz Equations
THREE DIMENSIONAL SYSTEMS Lecture 6: The Lorenz Equations 6. The Lorenz (1963) Equations The Lorenz equations were originally derived by Saltzman (1962) as a minimalist model of thermal convection in a
More informationEen vlinder in de wiskunde: over chaos en structuur
Een vlinder in de wiskunde: over chaos en structuur Bernard J. Geurts Enschede, November 10, 2016 Tuin der Lusten (Garden of Earthly Delights) In all chaos there is a cosmos, in all disorder a secret
More informationSPATIOTEMPORAL CHAOS IN COUPLED MAP LATTICE. Itishree Priyadarshini. Prof. Biplab Ganguli
SPATIOTEMPORAL CHAOS IN COUPLED MAP LATTICE By Itishree Priyadarshini Under the Guidance of Prof. Biplab Ganguli Department of Physics National Institute of Technology, Rourkela CERTIFICATE This is to
More informationChapter 6 - Ordinary Differential Equations
Chapter 6 - Ordinary Differential Equations 7.1 Solving Initial-Value Problems In this chapter, we will be interested in the solution of ordinary differential equations. Ordinary differential equations
More informationMathematical Foundations of Neuroscience - Lecture 7. Bifurcations II.
Mathematical Foundations of Neuroscience - Lecture 7. Bifurcations II. Filip Piękniewski Faculty of Mathematics and Computer Science, Nicolaus Copernicus University, Toruń, Poland Winter 2009/2010 Filip
More informationChapter 6: Ensemble Forecasting and Atmospheric Predictability. Introduction
Chapter 6: Ensemble Forecasting and Atmospheric Predictability Introduction Deterministic Chaos (what!?) In 1951 Charney indicated that forecast skill would break down, but he attributed it to model errors
More informationReal Randomness with Noise and Chaos
Real Randomness with Noise and Chaos by Kevin Fei working with Professor Rajarshi Roy, Professor Tom Murphy, and Joe Hart Random numbers are instrumental to modern computing. They are used by scientists
More informationMATH 415, WEEK 12 & 13: Higher-Dimensional Systems, Lorenz Equations, Chaotic Behavior
MATH 415, WEEK 1 & 13: Higher-Dimensional Systems, Lorenz Equations, Chaotic Behavior 1 Higher-Dimensional Systems Consider the following system of differential equations: dx = x y dt dy dt = xy y dz dt
More informationTHE ESSENCE OF... HAO~ Edward N. Lorenz UNIVERSITY OF WASHINGTON PRESS. Seattle
THE ESSENCE OF... HAO~ Edward N. Lorenz UNIVERSITY OF WASHINGTON PRESS Seattle The Butterfly Effect THE FOLLOWING is the text of a talk that I presented in a session devoted to the Global Atmospheric Research
More informationFORECASTING ECONOMIC GROWTH USING CHAOS THEORY
Article history: Received 22 April 2016; last revision 30 June 2016; accepted 12 September 2016 FORECASTING ECONOMIC GROWTH USING CHAOS THEORY Mihaela Simionescu Institute for Economic Forecasting of the
More informationhttp://www.ibiblio.org/e-notes/mset/logistic.htm On to Fractals Now let s consider Scale It s all about scales and its invariance (not just space though can also time And self-organized similarity
More informationComputational project: Modelling a simple quadrupole mass spectrometer
Computational project: Modelling a simple quadrupole mass spectrometer Martin Duy Tat a, Anders Hagen Jarmund a a Norges Teknisk-Naturvitenskapelige Universitet, Trondheim, Norway. Abstract In this project
More informationEdward Lorenz: Predictability
Edward Lorenz: Predictability Master Literature Seminar, speaker: Josef Schröttle Edward Lorenz in 1994, Northern Hemisphere, Lorenz Attractor I) Lorenz, E.N.: Deterministic Nonperiodic Flow (JAS, 1963)
More informationUsing Matlab to integrate Ordinary Differential Equations (ODEs)
Using Matlab to integrate Ordinary Differential Equations (ODEs) Erica McEvoy (Dated: June 17, 2009) 1. INTRODUCTION Ordinary differential equations tend to arise whenever you need to model changing quantities
More informationEdward Lorenz. Professor of Meteorology at the Massachusetts Institute of Technology
The Lorenz system Edward Lorenz Professor of Meteorology at the Massachusetts Institute of Technology In 1963 derived a three dimensional system in efforts to model long range predictions for the weather
More informationACM/CMS 107 Linear Analysis & Applications Fall 2016 Assignment 4: Linear ODEs and Control Theory Due: 5th December 2016
ACM/CMS 17 Linear Analysis & Applications Fall 216 Assignment 4: Linear ODEs and Control Theory Due: 5th December 216 Introduction Systems of ordinary differential equations (ODEs) can be used to describe
More informationThe 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 informationA Study of the Thermal Properties of a One. Dimensional Lennard-Jones System
A Study of the Thermal Properties of a One Dimensional Lennard-Jones System Abstract In this study, the behavior of a one dimensional (1D) Lennard-Jones (LJ) system is simulated. As part of this research,
More informationlab6 October 25, 2015
lab6 October 25, 2015 Contents 1. Objectives 2. Readings 3. Introduction 3.1 Using the integrator class 4. The Lorenz Equations 4.1 Boundedness of the Solutions 4.2 Steady States 4.3 Linearization of the
More information6.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 informationMAS212 Assignment #2: The damped driven pendulum
MAS Assignment #: The damped driven pendulum Sam Dolan (January 8 Introduction In this assignment we study the motion of a rigid pendulum of length l and mass m, shown in Fig., using both analytical and
More information9. Scientific Computing. Introduction to Computer Science Sedgewick and Wayne Copyright
9. Scientific Computing Introduction to Computer Science Sedgewick and Wayne Copyright 2007 http://www.cs.princeton.edu/introcs Applications of Scientific Computing Science and engineering challenges.
More informationLecture 6. Lorenz equations and Malkus' waterwheel Some properties of the Lorenz Eq.'s Lorenz Map Towards definitions of:
Lecture 6 Chaos Lorenz equations and Malkus' waterwheel Some properties of the Lorenz Eq.'s Lorenz Map Towards definitions of: Chaos, Attractors and strange attractors Transient chaos Lorenz Equations
More informationCorresponding address:
AN ABSTRACT OF THE THESIS OF Solveig T. Mueller for the degree of Honors Baccalaureate of Science in Mathematics presented on May 19, 2011. Title: Numerical Method and Step Size Variation in the Lorenz
More informationReport E-Project Henriette Laabsch Toni Luhdo Steffen Mitzscherling Jens Paasche Thomas Pache
Potsdam, August 006 Report E-Project Henriette Laabsch 7685 Toni Luhdo 7589 Steffen Mitzscherling 7540 Jens Paasche 7575 Thomas Pache 754 Introduction From 7 th February till 3 rd March, we had our laboratory
More informationModule E Unit 2 Lesson 2 Exploration 2 Explaining the accuracy of weather prediction
Module E Unit 2 Lesson 2 Exploration 2 Explaining the accuracy of weather prediction Objective Explore why weather can be predicted Analyze the cause and effect relationships that produce weather Do Now:
More informationDynamical Systems and Chaos Part I: Theoretical Techniques. Lecture 4: Discrete systems + Chaos. Ilya Potapov Mathematics Department, TUT Room TD325
Dynamical Systems and Chaos Part I: Theoretical Techniques Lecture 4: Discrete systems + Chaos Ilya Potapov Mathematics Department, TUT Room TD325 Discrete maps x n+1 = f(x n ) Discrete time steps. x 0
More informationObservability, a Problem in Data Assimilation
Observability, Data Assimilation with the Extended Kalman Filter 1 Observability, a Problem in Data Assimilation Chris Danforth Department of Applied Mathematics and Scientific Computation, UMD March 10,
More informationIn this report, I discuss two different topics of thermodynamics and explain how they have
Thermodynamics of Far-from-Equilibrium Systems: A Shift in Perception of Nature Truong Pham January 31st, 11 AME 36099; Directed Readings Prepared for Professor J. M. Powers 1 Introduction In this report,
More informationSolving PDEs with PGI CUDA Fortran Part 4: Initial value problems for ordinary differential equations
Solving PDEs with PGI CUDA Fortran Part 4: Initial value problems for ordinary differential equations Outline ODEs and initial conditions. Explicit and implicit Euler methods. Runge-Kutta methods. Multistep
More informationLecture 1 Monday, January 14
LECTURE NOTES FOR MATH 7394 ERGODIC THEORY AND THERMODYNAMIC FORMALISM INSTRUCTOR: VAUGHN CLIMENHAGA Lecture 1 Monday, January 14 Scribe: Vaughn Climenhaga 1.1. From deterministic physics to stochastic
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 informationA Classical Approach to the Stark-Effect. Mridul Mehta Advisor: Prof. Enrique J. Galvez Physics Dept., Colgate University
A Classical Approach to the Stark-Effect Mridul Mehta Advisor: Prof. Enrique J. Galvez Physics Dept., Colgate University Abstract The state of an atom in the presence of an external electric field is known
More informationApplications of Hurst Coefficient Analysis to Chaotic Response of ODE Systems: Part 1a, The Original Lorenz System of 1963
Applications of Hurst Coefficient Analysis to Chaotic Response of ODE Systems: Part 1a, The Original Lorenz System of 1963 Dan Hughes December 2008 Abstract I have coded the process for calculating Hurst
More informationWorksheet 8 Sample Solutions
Technische Universität München WS 2016/17 Lehrstuhl für Informatik V Scientific Computing Univ.-Prof. Dr. M. Bader 19.12.2016/21.12.2016 M.Sc. S. Seckler, M.Sc. D. Jarema Worksheet 8 Sample Solutions Ordinary
More informationCHEM-UA 652: Thermodynamics and Kinetics
CHEM-UA 65: Thermodynamics and Kinetics Notes for Lecture I. THE COMPLEXITY OF MULTI-STEP CHEMICAL REACTIONS It should be clear by now that chemical kinetics is governed by the mathematics of systems of
More informationPSTricks. pst-ode. A PSTricks package for solving initial value problems for sets of Ordinary Differential Equations(ODE), v0.12.
PSTricks pst-ode A PSTricks package for solving initial value problems for sets of Ordinary Differential Equations(ODE), v0.12 https://gitlab.com/agrahn/pst-ode 7th June 2018 Package author(s): Alexander
More informationPhysically Based Modeling: Principles and Practice Differential Equation Basics
Physically Based Modeling: Principles and Practice Differential Equation Basics Andrew Witkin and David Baraff Robotics Institute Carnegie Mellon University Please note: This document is 1997 by Andrew
More informationChaos & Recursive. Ehsan Tahami. (Properties, Dynamics, and Applications ) PHD student of biomedical engineering
Chaos & Recursive Equations (Properties, Dynamics, and Applications ) Ehsan Tahami PHD student of biomedical engineering Tahami@mshdiau.a.ir Index What is Chaos theory? History of Chaos Introduction of
More informationATM S 111, Global Warming Climate Models
ATM S 111, Global Warming Climate Models Jennifer Fletcher Day 27: July 29, 2010 Using Climate Models to Build Understanding Often climate models are thought of as forecast tools (what s the climate going
More informationControlling a Novel Chaotic Attractor using Linear Feedback
ISSN 746-7659, England, UK Journal of Information and Computing Science Vol 5, No,, pp 7-4 Controlling a Novel Chaotic Attractor using Linear Feedback Lin Pan,, Daoyun Xu 3, and Wuneng Zhou College of
More informationLesson 9: Predator-Prey and ode45
Lesson 9: Predator-Prey and ode45 9.1 Applied Problem. In this lesson we will allow for more than one population where they depend on each other. One population could be the predator such as a fox, and
More informationPhysically Based Modeling Differential Equation Basics
Physically Based Modeling Differential Equation Basics Andrew Witkin and David Baraff Pixar Animation Studios Please note: This document is 2001 by Andrew Witkin and David Baraff. This chapter may be freely
More informationWhy 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 informationUnit Ten Summary Introduction to Dynamical Systems and Chaos
Unit Ten Summary Introduction to Dynamical Systems Dynamical Systems A dynamical system is a system that evolves in time according to a well-defined, unchanging rule. The study of dynamical systems is
More informationChaos Control for the Lorenz System
Advanced Studies in Theoretical Physics Vol. 12, 2018, no. 4, 181-188 HIKARI Ltd, www.m-hikari.com https://doi.org/10.12988/astp.2018.8413 Chaos Control for the Lorenz System Pedro Pablo Cárdenas Alzate
More informationBoyce/DiPrima/Meade 11 th ed, Ch 1.1: Basic Mathematical Models; Direction Fields
Boyce/DiPrima/Meade 11 th ed, Ch 1.1: Basic Mathematical Models; Direction Fields Elementary Differential Equations and Boundary Value Problems, 11 th edition, by William E. Boyce, Richard C. DiPrima,
More informationThree-Tank Experiment
Three-Tank Experiment Overview The three-tank experiment focuses on application of the mechanical balance equation to a transient flow. Three tanks are interconnected by Schedule 40 pipes of nominal diameter
More informationLesson 4: Population, Taylor and Runge-Kutta Methods
Lesson 4: Population, Taylor and Runge-Kutta Methods 4.1 Applied Problem. Consider a single fish population whose size is given by x(t). The change in the size of the fish population over a given time
More informationNonlinear dynamics & chaos BECS
Nonlinear dynamics & chaos BECS-114.7151 Phase portraits Focus: nonlinear systems in two dimensions General form of a vector field on the phase plane: Vector notation: Phase portraits Solution x(t) describes
More informationComputers, Lies and the Fishing Season
1/47 Computers, Lies and the Fishing Season Liz Arnold May 21, 23 Introduction Computers, lies and the fishing season takes a look at computer software programs. As mathematicians, we depend on computers
More informationChaos in multiplicative systems
Chaos in multiplicative systems Dorota Aniszewska 1 and Marek Rybaczuk 2 1 Institute of Materials Science and Applied Mechanics Wroclaw University of Technology 50-370 Wroclaw, Smoluchowskiego 25 (e-mail:
More informationExploring Nonlinear Oscillator Models for the Auditory Periphery
Exploring Nonlinear Oscillator Models for the Auditory Periphery Andrew Binder Dr. Christopher Bergevin, Supervisor October 31, 2008 1 Introduction Sound-induced vibrations are converted into electrical
More informationThe Nonlinear Pendulum
The Nonlinear Pendulum Evan Sheridan 11367741 Feburary 18th 013 Abstract Both the non-linear linear pendulum are investigated compared using the pendulum.c program that utilizes the trapezoid method for
More informationDiscussion of the Lorenz Equations
Discussion of the Lorenz Equations Leibniz Universität Hannover Proseminar Theoretische Physik SS/2015 July 22, 2015 (LUH) Lorenz Equations July 22, 2015 1 / 42 Outline 1 2 3 4 5 6 7 8 (LUH) Lorenz Equations
More informationNumerical Methods for Ordinary Differential Equations
CHAPTER 1 Numerical Methods for Ordinary Differential Equations In this chapter we discuss numerical method for ODE. We will discuss the two basic methods, Euler s Method and Runge-Kutta Method. 1. Numerical
More information11 Chaos in Continuous Dynamical Systems.
11 CHAOS IN CONTINUOUS DYNAMICAL SYSTEMS. 47 11 Chaos in Continuous Dynamical Systems. Let s consider a system of differential equations given by where x(t) : R R and f : R R. ẋ = f(x), The linearization
More informationIntroduction Knot Theory Nonlinear Dynamics Topology in Chaos Open Questions Summary. Topology in Chaos
Introduction Knot Theory Nonlinear Dynamics Open Questions Summary A tangled tale about knot, link, template, and strange attractor Centre for Chaos & Complex Networks City University of Hong Kong Email:
More informationIntroduction 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 informationFundamentals of Fluid Dynamics: Ideal Flow Theory & Basic Aerodynamics
Fundamentals of Fluid Dynamics: Ideal Flow Theory & Basic Aerodynamics Introductory Course on Multiphysics Modelling TOMASZ G. ZIELIŃSKI (after: D.J. ACHESON s Elementary Fluid Dynamics ) bluebox.ippt.pan.pl/
More informationFORTRAN 77 Lesson 7. Reese Haywood Ayan Paul
FORTRAN 77 Lesson 7 Reese Haywood Ayan Paul 1 Numerical Integration A general all purpose Numerical Integration technique is the Trapezoidal Method. For most functions, this method is fast and accurate.
More informationUse of Differential Equations In Modeling and Simulation of CSTR
Use of Differential Equations In Modeling and Simulation of CSTR JIRI VOJTESEK, PETR DOSTAL Department of Process Control, Faculty of Applied Informatics Tomas Bata University in Zlin nám. T. G. Masaryka
More informationIMHO In My Humble Opinion
IMHO In My Humble Opinion Home Chaos Theory: A Brief Introduction What exactly is chaos? The name "chaos theory" comes from the fact that the systems that the theory describes are apparently disordered,
More informationLorenz Equations. Lab 1. The Lorenz System
Lab 1 Lorenz Equations Chaos: When the present determines the future, but the approximate present does not approximately determine the future. Edward Lorenz Lab Objective: Investigate the behavior of a
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 informationOrdinary Differential Equations. Monday, October 10, 11
Ordinary Differential Equations Monday, October 10, 11 Problems involving ODEs can always be reduced to a set of first order differential equations. For example, By introducing a new variable z, this can
More informationDynamical Systems. 1.0 Ordinary Differential Equations. 2.0 Dynamical Systems
. Ordinary Differential Equations. An ordinary differential equation (ODE, for short) is an equation involves a dependent variable and its derivatives with respect to an independent variable. When we speak
More informationA New Science : Chaos
A New Science : Chaos Li Shi Hai Department of Mathematics National University of Singapore In the new movie Jurassic Park [C], Malcolm, a mathematician specialized in Chaos Theory, explained that Hammond's
More informationChapter 2 Chaos theory and its relationship to complexity
Chapter 2 Chaos theory and its relationship to complexity David Kernick This chapter introduces chaos theory and the concept of non-linearity. It highlights the importance of reiteration and the system
More informationSome tools and methods for determination of dynamics of hydraulic systems
Some tools and methods for determination of dynamics of hydraulic systems A warm welcome to the course in Hydraulic servo-techniques! The purpose of the exercises given in this material is to make you
More informationTheoretical physics. Deterministic chaos in classical physics. Martin Scholtz
Theoretical physics Deterministic chaos in classical physics Martin Scholtz scholtzzz@gmail.com Fundamental physical theories and role of classical mechanics. Intuitive characteristics of chaos. Newton
More informationNumerical Methods - Initial Value Problems for ODEs
Numerical Methods - Initial Value Problems for ODEs Y. K. Goh Universiti Tunku Abdul Rahman 2013 Y. K. Goh (UTAR) Numerical Methods - Initial Value Problems for ODEs 2013 1 / 43 Outline 1 Initial Value
More informationNewton s Second Law of Motion. Isaac Newton (in 1689) The Universe is Regular and Predictable. The Foundation of Science
Key Concepts and Fundamental Theorems Isaac Newton (in 689) of Atmospheric Science This is the earliest portrait of Newton to survive. The artist was Godfrey Kneller, perhaps the greatest portrait painter
More information9. Scientific Computing
9. Scientific Computing Introduction to Computer Science in Java: An Interdisciplinary Approach Robert Sedgewick and Kevin Wayne Copyright 2002 2010 4/16/11 9:11 AM Applications of Scientific Computing
More informationUsing Artificial Neural Networks (ANN) to Control Chaos
Using Artificial Neural Networks (ANN) to Control Chaos Dr. Ibrahim Ighneiwa a *, Salwa Hamidatou a, and Fadia Ben Ismael a a Department of Electrical and Electronics Engineering, Faculty of Engineering,
More informationMath 216 Final Exam 14 December, 2012
Math 216 Final Exam 14 December, 2012 This sample exam is provided to serve as one component of your studying for this exam in this course. Please note that it is not guaranteed to cover the material that
More informationM.Sc. in Meteorology. Numerical Weather Prediction
M.Sc. in Meteorology UCD Numerical Weather Prediction Prof Peter Lynch Meteorology & Climate Centre School of Mathematical Sciences University College Dublin Second Semester, 2005 2006. In this section
More informationNumerical Methods for Ordinary Differential Equations
Numerical Methods for Ordinary Differential Equations By Brian D. Storey 1. Introduction Differential equations can describe nearly all systems undergoing change. They are ubiquitous is science and engineering
More informationSimple conservative, autonomous, second-order chaotic complex variable systems.
Simple conservative, autonomous, second-order chaotic complex variable systems. Delmar Marshall 1 (Physics Department, Amrita Vishwa Vidyapeetham, Clappana P.O., Kollam, Kerala 690-525, India) and J. C.
More informationWright-Fisher Models, Approximations, and Minimum Increments of Evolution
Wright-Fisher Models, Approximations, and Minimum Increments of Evolution William H. Press The University of Texas at Austin January 10, 2011 1 Introduction Wright-Fisher models [1] are idealized models
More informationMath 266: Ordinary Differential Equations
Math 266: Ordinary Differential Equations Long Jin Purdue University, Spring 2018 Basic information Lectures: MWF 8:30-9:20(111)/9:30-10:20(121), UNIV 103 Instructor: Long Jin (long249@purdue.edu) Office
More informationNONLINEAR DYNAMICS AND CHAOS. Numerical integration. Stability analysis
LECTURE 3: FLOWS NONLINEAR DYNAMICS AND CHAOS Patrick E McSharr Sstems Analsis, Modelling & Prediction Group www.eng.o.ac.uk/samp patrick@mcsharr.net Tel: +44 83 74 Numerical integration Stabilit analsis
More informationIntegration of Differential Equations
Integration of Differential Equations Gabriel S. Cabrera August 4, 018 Contents 1 Introduction 1 Theory.1 Linear 1 st Order ODEs................................. 3.1.1 Analytical Solution...............................
More informationTaylor series. Chapter Introduction From geometric series to Taylor polynomials
Chapter 2 Taylor series 2. Introduction The topic of this chapter is find approximations of functions in terms of power series, also called Taylor series. Such series can be described informally as infinite
More information(Refer Slide Time: 0:35)
Fluid Dynamics And Turbo Machines. Professor Dr Shamit Bakshi. Department Of Mechanical Engineering. Indian Institute Of Technology Madras. Part A. Module-1. Lecture-4. Tutorial. (Refer Slide Time: 0:35)
More informationResearch Computing with Python, Lecture 7, Numerical Integration and Solving Ordinary Differential Equations
Research Computing with Python, Lecture 7, Numerical Integration and Solving Ordinary Differential Equations Ramses van Zon SciNet HPC Consortium November 25, 2014 Ramses van Zon (SciNet HPC Consortium)Research
More informationSolving Zhou Chaotic System Using Fourth-Order Runge-Kutta Method
World Applied Sciences Journal 21 (6): 939-944, 2013 ISSN 11-4952 IDOSI Publications, 2013 DOI: 10.529/idosi.wasj.2013.21.6.2915 Solving Zhou Chaotic System Using Fourth-Order Runge-Kutta Method 1 1 3
More informationNumerical Methods - Preliminaries
Numerical Methods - Preliminaries Y. K. Goh Universiti Tunku Abdul Rahman 2013 Y. K. Goh (UTAR) Numerical Methods - Preliminaries 2013 1 / 58 Table of Contents 1 Introduction to Numerical Methods Numerical
More informationMATH 215/255 Solutions to Additional Practice Problems April dy dt
. For the nonlinear system MATH 5/55 Solutions to Additional Practice Problems April 08 dx dt = x( x y, dy dt = y(.5 y x, x 0, y 0, (a Show that if x(0 > 0 and y(0 = 0, then the solution (x(t, y(t of the
More information