Non-Dimensionalization of Differential Equations. September 9, 2014
|
|
- Charla Watkins
- 5 years ago
- Views:
Transcription
1 Non-Dimensionalization of Differential Equations September 9, Non-dimensionalization The most elementary idea of dimensional analysis is this: in any equation, the dimensions of the right hand side and the left hand side must match up. For example, an equation describing the velocity of a ball v = f(v 0, g, t) must have units of velocity, say, meters per second, on both sides or it does not make any sense. If we are measuring velocity in meters per second, the equation will have some constants in it that are related to our units of measure, meters and seconds. The specific form of the velocity equation is v = v 0 + gt, where v 0 is the original velocity, g is the gravitational constant of the Earth and t is time. Including the units, the equation becomes v m/sec = v 0 m/sec 9.8 m/sec 2 t sec, and we can see that the dimensions of the equations are m/sec = m/sec + (m/sec 2 ) sec, which matches up. The problem with this equation is that if we gave it to someone who measured velocities in feet per second, she would get the wrong answer if she didn t convert the measurements on her instruments to meters per second. In fact, the multi million dollar Mars Climate Probe was destroyed because of just this kind of error. Engineers working on different parts of the problem were using different units and one group forgot to convert from English units to SI units before sending their findings to another group. The solution to this problem, it seems, is to express it in non-dimensional form v = v 0 + gt. As long as v is measured in feet per minute, or meters per second, or even decicubits per fortnight, with the other physical quantities being similarly measured and the constant g given in that system, the solution will be correct. We naturally do dimensional matching almost unconsciously (if we have been at this long enough!) For instance, suppose we wish to model the concentration of salt in a lake. More specifically, suppose the lake has cubic meters of water before the salt started polluting it, and that it has a single inlet and a single outlet, both of which flow at a rate of 3 cubic meters per second, and that at some time t 0 a manufacturing plant upriver begins dumping salt in the river, so that when the river flows into the lake, it has a concentration of 1 kilogram of salt per cubic meter. To develop a differential equations describing the situation, we form a difference that describes the change in salt over an interval of time t. Represent the salt at time t by S(t), and the salt at t later by S(t + t); then the change in salt over the period t is S(t + t) S(t),. This goes on the left hand side and has units of kilograms. On the right hand side, we describe the physical situation, change in mass of salt per unit time times the time t: this is salt per unit time in minus salt per unit time out, times t. Specifically, change in salt is (1kg/m 3 3m 3 /sec) (S m 3 /sec), and the change 1
2 in time t sec is ( S) kg t. The last step is to divide both sides by t and take the limit as t goes to zero (now ignoring the units of kilograms): S(t + t) S(t) lim = ds t 0 t dt = S. We can solve this ODE and have a formula for the total amount of salt in the lake at any time, and we can use the volume of the lake to determine the concentration. Now, while we were able to cancel out the SI unit names specifically, and we ignored the remaining mass unit, we still have units involved in this equation implicitly in the form of the constants 10 6 and 3, both of which arose from the units used to measure the volume of the lake, the concentration of salt and the rate of flow of the rivers. How can we make an equation that automatically accounts for differing units of measure, like we did for the velocity equation? And how do we find equations like the velocity equation and much more complicated modeling tasks, for instance, when we don t know for sure what the physics or the chemistry and so on actually are, or when we don t know what variables might be the most important and which might be ignored? The answer lies in a more formal and exacting application of dimensional analysis, sometimes called non-dimensionalizing equations or finding dimensionally homogeneous equations. 1.1 Example from physics Before we get into the full machinery of the method, let s do an example from physics in the way it is usually done in a physics class, and analyze the results. We begin by assembling the differential equation for the unforced frictionless pendulum using Newtons second law F = ma, force equals mass times acceleration, and with the assumption that all of the mass of the pendulum is concentrated in the center of the bob. On the left, the mathematical expression for the mass times angular acceleration m r θ, where m is the mass of the bob, r its length, θ is measured counterclockwise from the straight down position, and we have used the standard physicists notation for the second derivative with respect to time, θ d2 θ dt. On the right will be the specific physical forces 2 in the θ direction, g m r sin θ, so we have m r θ = g m sin θ. Now divide through by mr, which is never zero, to get θ = g r sin θ. Here is a trick worth remembering. There is no reason to choose to measure time in seconds. Let s choose a new unit of time measurement τ = k t for some k which we haven t yet determined. Now suppose we are measuring time in these units. Then θ = θ(τ) = θ(τ(t)). Using the chain rule, we find θ = d dt ( ) dθ = d ( dθ dt dt dτ ) dτ = d dt dt ( ) dθ k = d2 θ dτ dτ dτ 2 dt k = d2 θ dτ 2 k2 = k 2 θ, where the primes in θ indicate differentiation with respect to the new time variable τ. This simply says that the change in a property with respect to one measure of time is a constant times the change in that property with respect to another measure of time. This makes good sense. Now we have k 2 θ = g r sin θ, 2
3 Figure 1: The physical pendulum. which we divide by k to get θ = g 1 sin θ. r k2 If we let k = g r, then we see that k has units of L becomes θ + sin θ = 0. T 2 L = 1 T and the differential equation In terms of the actual constants of gravity and length, we are measuring time as a function of the gravitational constant and the length of the pendulum rod. There is something much more interesting happening here, as we shall soon see. Now, to see what k = g r really means, let s go back to the original equation θ = g r sin θ and solve it using the usual methods from differential equations. Because it is nonlinear, we will use a trick to approximate the solution by approximating the sine function with the linear term of its power series (Taylor Series) expansion sin θ = θ 1 3 θ The linear approximation sin θ θ is good to within about 1% for angles θ < 15, so if we solve θ + g r θ = 0, the solution will be very close to correct for small excursions of the pendulum. 3
4 Figure 2: The sine function and the function θ. Using the usual methods for second order linear ODEs, we find g g θ = K 1 sin r t + K 2 cos r t. We see that the scaling we chose for time in the full model of the pendulum is the frequency of the solution to the linearized problem. So the natural time scale to use when modeling a pendulum is the time it takes to complete a small amplitude swing, actually the time it takes for a small amplitude swing in the limit as the amplitude goes to zero (Why?). We also see that it is the magnitude of the gravitational acceleration and the length of the shaft that determine the period of the pendulum, and the mass does not figure into it at all. Actually, this is not quite correct, because we assumed the pendulum shaft was rigid and had no mass (incompatible assumptions!) and that the mass of the bob was concentrated at the center (more realistic). In a real pendulum, the distributed mass of the shaft displaces the actual center of gravity of the pendulum (see any calculus book for moment of inertia formulas for various shapes). Now, why should this be? Why should a physical object like a pendulum prefer one particular measure of time over another? Does it have a preferred length scale too? A preferred gravitational acceleration or velocity? The answer is that a system prefers a system of measurements that makes it dimensionless. A dimensionless system is the simplest possible version of the system, and the best, because it is scalable, in other words, it models the phenomena over a range of sizes. In particular, if we have two pendula one a meter long and one a foot long, measuring time in units of r g with r in meters and g as 9.8 meters per second squared in the first, and r in feet and g as 16 feet per second squared in the second will result in a dimensionally homogeneous equation. Thus the same solution is valid for either case, and any variation of the gravity or pendulum length will only scale the time axis. Adding friction, the original equation would become ml θ = γ θ mg sin θ. 4
5 Applying the same technique as before, we find resulting in, after a little algebra, k 2 θ = γ m kθ g sin θ, r θ + γ m r g θ + sin θ = 0. What should we do about this last constant? The answer will come later. 1.2 Example from ecology In a closed ecological system reside a population of rabbits who have all the clover they need to grow without bound, and a population of foxes who depend entirely on the rabbits for their sustenance. Because the population of rabbits will grow exponentially, while the population of foxes will die out exponentially, a first pass at a system, without interaction, is f = k 1 f ṙ = k 2 r. Adding an interactioon term, with the assumption that the interaction between foxes and rabbits (fox catches and eats rabbit) can be modeled as a proportionality depending on the joint probability of a rabbit and fox occupying some small capture region, we obtain 5
6 f = k 1 f + k 3 fr ṙ = k 2 r k 4 fr. In this formulation, we have four free parameters, k 1, k 2, k 3, k 4, some of which are probably unnecessary. Because we are free to choose the units we use to measure time (days, weeks, years) and number of creatures (number of individuals, number of average family sizes, percentage of total population, etc.) we can choose our units to eliminate some of these parameters, by combining them in dimensionless groups. With that in mind, let us propose new variables proportional to our current ones, say τ = at, F = bf, R = cr. We will adjust a, b, c to give the simplest possible model. Then substitute f = 1 b F, r = 1 c R, f df dt = a b df dτ a b F, ṙ dr dt = a c dr dτ a c R to get or a b F 1 = k 1 b F + k 1 3 bc F R (1) a 1 c R = k 2 c R k 1 4 F R, bc (2) F 1 = k 1 a F + k 1 3 ac F R (3) R 1 = k 2 a R k 1 4 F R. ba (4) Let s eliminate constant multipliers one by one. It is clear that let a = k 1 or a = k 2. Which one we choose is arbitrary we are either choosing the time rescaling a to be the growth rate of the rabbits or the decay rate of the foxes. Choosing a = k 1, we obtain F = F + k 3 k 1 c F R (5) R = k 2 k 1 R k 4 bk 1 F R. (6) Now, if we choose c = k3 k 1 and b = k4 k 1, and rename the ratio k2 k 1 α, we have eliminated all but one parameter. The single parameter α, the ratio of the growth of rabbits to the decay of foxes apparently has all the important information in it. F = F + F R (7) R = αr F R. (8) By choosing to rescale our units for the fox and rabbit populations and our measurement of time, we have simplified the problem immensely. We can now run numerical experiments and analytical procedures on our model over a realistic rance of a single parameter. It is important to realize that we have done noting but change the density of tics on the three axes. Several plots of solutions to this system can be seen in Figure(1.2). 6
7 2 Stability Analysis For linear systems (and also for linearization of non-linear systems, as we shall see) the eigenvalues and eigenvectors completely characterize the local qualitative dynamics of most equilibrium points. Suppose we have a linear system ẋ = Ax, or in more detail, [ ẋ ẏ ] = [ a b c d ] [ x y If b = c = 0 then the system is decoupled and the solutions have the form x = x 0 e at, y = y 0 e dt. Otherwise solutions will be (generically) of the form x = x 0 e λ1t v 1, ȳ = ȳ 0 e λ2t v 2 for some constants λ 1, λ 2 and vectors v 1, v 2. Substituting x for x, we get Av 1 = λ 1 v 1 or (A λ 1 I)v 1 = 0, which, for v 1 0, has a solution for λ 1 which can be found by solving det(a λ 1 I) = 0, and likewise for λ 2, v 2. Upon finding λ 1 and λ 2, these can be substituted into (A λi)v = 0 to find v 1, v 2. With the eigenvectors v 1, v 2 and eigenvalues λ 1, λ 2 in hand, we can decompose the matrix A in the problem ẋ = Ax as A = P DP 1, where P is a matrix whose columns are the eigenvectors v 1, v 2, D has the eigenvalues λ 1, λ 2 on its diagonal and zeroes elsewhere, to obtain ẋ = P DP 1 x or P 1 ẋ = DP 1 x. Now two vectors are identical whether written in one basis or another, thus Ix = P x represent the same vector, so the components of x, the coefficients of the column vectors v 1, v 2 making up P are given by P 1 x = x. In other words, P 1 rewrites x in the basis spanned by the columns of P, and the differential equation is now in the form x = D x, which is uncoupled in the new variables. Now we solve x = D x to get x = x 0 e λ1t v 1 and ȳ = ȳ 0 e λ2t v 2. Linear combinations of these two solutions (with coefficients x 1, ȳ 0 give solutions beginning anywhere in the phase space. We can then convert this solution back to a solution in the standard basis by writing v 1, v 2 as linear combinations [ ] of the[ standard ] basis vectors, [ ] and[ taking ] linear combinations of these. Write v 1 = p 0 + q 1 and v 2 = r 0 + s 1. Then [ e λ1t v 1 = e λ1t 1 (p 0 ]. ] [ 0 + q 1 ] ) and so [ ] [ ] e λ2t v 2 = e λ2t (r + s ), x = pe λ1t + re λ2t 7
8 y = qe λ1t + se λ2t. Obviously λ i can be any complex conjugate pair. Depending on the nature of these numbers, the solutions can have very different natures. The stability type of the equilibrium point [ is ] determined by the form of the eigenvalues. The eigenvalues for a matrix A = a b c d are found by solving det(a λi) = [ ] [ ] [ ] a b 1 0 c d λ 0 1 = a λ b c d λ = 0 or (a λ)(d λ) bc = λ 2 (a + d)λ + (ad bc) = 0, which has solutions λ 1,2 = (a + d) ± (a + d) 2 4(ad bc) 2 = τ 2 ± ( τ 2 )2, where we have renamed the trace τ = a+d and the determinant = ad bc of the matrix A. Then, depending on whether τ or the quantity under the radical is negative, positive or zero, we will have either stable, unstable or neutral dynamics, which can be subdivided into several topological types shown in the chart of Figure 3. Figure 3: This diagram encodes the equilibrium type in terms of the sign of the trace and the relative sizes of the trace and determinant. In addition to the equilibria shown, there are other types of limiting behavior, most notably limit cycles and strange attractors. A limit cycle is a solution z(t) with the property that z(t + T ) = z(t) for some finite T. A limit cycle is a topological circle, and may be knotted if it lives in three dimensions. A strange attractor is the limit set for a system that is neither an equilibrium nor a limit cycle. It is a much more complicated object that we will not discuss here, except to note that ordinary differential equation models for such simple systems as a forced oscillator, pendulum or epidemic model can exhibit the chaotic behavior that leads to a strange attractor. 8
9 Figure 4: The various types of equilibria determined by the eigenvalues of the local linear manifolds. A simple example of a differential equation that can lead to a limit cycle is dy dx = y 3 +x 2 y x y x 3 +xy 2 x+y. This can be written more simply in polar coordinates as θ = C, ṙ = ar(1 r), and its solution curves are shown in 5. By plotting ṙ against r in the second equation of this system we can see that except for the equilibria r = 0, r = 1, any initial conditions will approach the circle r = 1 by spiraling counterclockwise in or out. Figure 5: An attracting cycle with equation θ = C, ṙ = ar(1 r) in polar coordinates. The Lotka-Volterra predator-prey model ẋ = ax bxy, ẏ = cy + dxy exhibits a limit cycle, as does the van der Pol oscillator ÿ ɛ(1 y 2 )ẏ + y = 0. The graphs for solutions to these equations are shown below in Figure 6. 9
10 Figure 6: Limit cycles of the Lotka-Volterra system (left) and the van der Pol oscillator (right). A linear system has an equilibrium point at (0, 0) and nowhere else. Nonlinear systems can have any number of equilibrium points. We can perform the same kind of analysis near each of these equilibrium points if we linearize the local dynamics. The next section elaborates on this idea. 3 Local Linearization If a system is linear, we can determine the stability type of the equilibrium point by the method above for 2 2 systems, and can, at least in principle, find an algebraic solution for the eigenvalues and eigenvectors of linear systems up to 4 4. Beyond this, we are solving a fifth degree or higher polynomial to get the eigenvalues, and there is no general solution to these in terms of radicals and powers. More often than not, our model is a nonlinear system of ODEs, and in addition there are equilibrium points other than the origin. In these cases, we would like to build a locally linear system and do the same analysis we did for linear systems near each equilibrium point. For most types of equilibria, the eigenvalues of the linear system are sufficient to determine entirely the qualitative dynamics near the equilibrium point. Of all the types of equilibria, only the true nonlinear dynamics of centers of the linearized system remain a mystery. We will introduce a method to deal with problem cases like these shortly. Now we look at the method of local linearization for a two dimensional system. The method generalizes easily to higher dimensional systems. Suppose we have an equilibrium point at (x, y ). Then near that point we can set up a coordinate system (u, v) = (x x, y = [ y ),] and[ nearby, we ] may expand the right hand[ side of ] the [ nonlinear differential ] ẋ f(x, y) f(u + x equations ẏ = g(x, y) in a Taylor series. We get u v =, v + y ) g(u + x, v + y, ) where f(u + x, v + y ) = 1 0! f(x, y ) + 1 1! (u1 f x (x, y ) + v 1 f y (x, y )) ! (u2 f xx (x, y ) + 2uvf xy (x, y ) + v 2 f yy (x, y ))
11 Figure 7: Nonlinear stable and unstable manifolds W s and W u may be well approximated by stable and unstable linear manifolds E s and E u near the equilibrium. and g(u + x, v + y ) = 1 0! (u0 + v 0 )g(x, y ) + 1 1! (u1 g x (x, y ) + v 1 g y (x, y )) ! (u2 g xx (x, y ) + 2uvg xy (x, y ) + v 2 g yy (x, y )) Now because (x, y ) is an equilibrium point, f(x, y ) = 0 and g(x, y ) = 0. Ignoring quadratic and higher terms, we obtain for f and g or in matrix form, f(x, y) = uf x (x, y ) + vf y (x, y ), g(x, y) = ug x (x, y ) + vg y (x, y ), [ u v ] = [ ] fx f y g x g y (x,y ) [ ] fx f The matrix J = y g x g is known as the Jacobian matrix, and we can find its eigenvectors and eigenvalues at each equilibrium point (x, y ) i to determine the local y dynamics. [ u v ]. 11
Fundamentals of Dynamical Systems / Discrete-Time Models. Dr. Dylan McNamara people.uncw.edu/ mcnamarad
Fundamentals of Dynamical Systems / Discrete-Time Models Dr. Dylan McNamara people.uncw.edu/ mcnamarad Dynamical systems theory Considers how systems autonomously change along time Ranges from Newtonian
More information28. Pendulum phase portrait Draw the phase portrait for the pendulum (supported by an inextensible rod)
28. Pendulum phase portrait Draw the phase portrait for the pendulum (supported by an inextensible rod) θ + ω 2 sin θ = 0. Indicate the stable equilibrium points as well as the unstable equilibrium points.
More informationPhysics: spring-mass system, planet motion, pendulum. Biology: ecology problem, neural conduction, epidemics
Applications of nonlinear ODE systems: Physics: spring-mass system, planet motion, pendulum Chemistry: mixing problems, chemical reactions Biology: ecology problem, neural conduction, epidemics Economy:
More informationLinearization of Differential Equation Models
Linearization of Differential Equation Models 1 Motivation We cannot solve most nonlinear models, so we often instead try to get an overall feel for the way the model behaves: we sometimes talk about looking
More informationENGI 9420 Lecture Notes 4 - Stability Analysis Page Stability Analysis for Non-linear Ordinary Differential Equations
ENGI 940 Lecture Notes 4 - Stability Analysis Page 4.01 4. Stability Analysis for Non-linear Ordinary Differential Equations A pair of simultaneous first order homogeneous linear ordinary differential
More informationLotka Volterra Predator-Prey Model with a Predating Scavenger
Lotka Volterra Predator-Prey Model with a Predating Scavenger Monica Pescitelli Georgia College December 13, 2013 Abstract The classic Lotka Volterra equations are used to model the population dynamics
More informationAutonomous Systems and Stability
LECTURE 8 Autonomous Systems and Stability An autonomous system is a system of ordinary differential equations of the form 1 1 ( 1 ) 2 2 ( 1 ). ( 1 ) or, in vector notation, x 0 F (x) That is to say, an
More informationProblem set 7 Math 207A, Fall 2011 Solutions
Problem set 7 Math 207A, Fall 2011 s 1. Classify the equilibrium (x, y) = (0, 0) of the system x t = x, y t = y + x 2. Is the equilibrium hyperbolic? Find an equation for the trajectories in (x, y)- phase
More informationSystems of Linear ODEs
P a g e 1 Systems of Linear ODEs Systems of ordinary differential equations can be solved in much the same way as discrete dynamical systems if the differential equations are linear. We will focus here
More informationSolutions to Dynamical Systems 2010 exam. Each question is worth 25 marks.
Solutions to Dynamical Systems exam Each question is worth marks [Unseen] Consider the following st order differential equation: dy dt Xy yy 4 a Find and classify all the fixed points of Hence draw the
More informationSystems of Ordinary Differential Equations
Systems of Ordinary Differential Equations Scott A. McKinley October 22, 2013 In these notes, which replace the material in your textbook, we will learn a modern view of analyzing systems of differential
More information1 The pendulum equation
Math 270 Honors ODE I Fall, 2008 Class notes # 5 A longer than usual homework assignment is at the end. The pendulum equation We now come to a particularly important example, the equation for an oscillating
More information7 Planar systems of linear ODE
7 Planar systems of linear ODE Here I restrict my attention to a very special class of autonomous ODE: linear ODE with constant coefficients This is arguably the only class of ODE for which explicit solution
More informationNonlinear Autonomous Systems of Differential
Chapter 4 Nonlinear Autonomous Systems of Differential Equations 4.0 The Phase Plane: Linear Systems 4.0.1 Introduction Consider a system of the form x = A(x), (4.0.1) where A is independent of t. Such
More informationCoordinate Curves for Trajectories
43 The material on linearizations and Jacobian matrices developed in the last chapter certainly expanded our ability to deal with nonlinear systems of differential equations Unfortunately, those tools
More informationSystems of Differential Equations: General Introduction and Basics
36 Systems of Differential Equations: General Introduction and Basics Thus far, we have been dealing with individual differential equations But there are many applications that lead to sets of differential
More informationChapter 6. Second order differential equations
Chapter 6. Second order differential equations A second order differential equation is of the form y = f(t, y, y ) where y = y(t). We shall often think of t as parametrizing time, y position. In this case
More informationMAT 22B - Lecture Notes
MAT 22B - Lecture Notes 4 September 205 Solving Systems of ODE Last time we talked a bit about how systems of ODE arise and why they are nice for visualization. Now we'll talk about the basics of how to
More informationSTABILITY. Phase portraits and local stability
MAS271 Methods for differential equations Dr. R. Jain STABILITY Phase portraits and local stability We are interested in system of ordinary differential equations of the form ẋ = f(x, y), ẏ = g(x, y),
More informationSolving Differential Equations: First Steps
30 ORDINARY DIFFERENTIAL EQUATIONS 3 Solving Differential Equations Solving Differential Equations: First Steps Now we start answering the question which is the theme of this book given a differential
More informationContinuous time population models
Continuous time population models Jaap van der Meer jaap.van.der.meer@nioz.nl Abstract Many simple theoretical population models in continuous time relate the rate of change of the size of two populations
More informationSolutions for B8b (Nonlinear Systems) Fake Past Exam (TT 10)
Solutions for B8b (Nonlinear Systems) Fake Past Exam (TT 10) Mason A. Porter 15/05/2010 1 Question 1 i. (6 points) Define a saddle-node bifurcation and show that the first order system dx dt = r x e x
More informationComplex Dynamic Systems: Qualitative vs Quantitative analysis
Complex Dynamic Systems: Qualitative vs Quantitative analysis Complex Dynamic Systems Chiara Mocenni Department of Information Engineering and Mathematics University of Siena (mocenni@diism.unisi.it) Dynamic
More informationChapter 9 Global Nonlinear Techniques
Chapter 9 Global Nonlinear Techniques Consider nonlinear dynamical system 0 Nullcline X 0 = F (X) = B @ f 1 (X) f 2 (X). f n (X) x j nullcline = fx : f j (X) = 0g equilibrium solutions = intersection of
More informationChapter 4. Systems of ODEs. Phase Plane. Qualitative Methods
Chapter 4 Systems of ODEs. Phase Plane. Qualitative Methods Contents 4.0 Basics of Matrices and Vectors 4.1 Systems of ODEs as Models 4.2 Basic Theory of Systems of ODEs 4.3 Constant-Coefficient Systems.
More informationthen the substitution z = ax + by + c reduces this equation to the separable one.
7 Substitutions II Some of the topics in this lecture are optional and will not be tested at the exams. However, for a curious student it should be useful to learn a few extra things about ordinary differential
More informationMATH 320, WEEK 11: Eigenvalues and Eigenvectors
MATH 30, WEEK : Eigenvalues and Eigenvectors Eigenvalues and Eigenvectors We have learned about several vector spaces which naturally arise from matrix operations In particular, we have learned about the
More informationPeriodic Motion. Periodic motion is motion of an object that. regularly repeats
Periodic Motion Periodic motion is motion of an object that regularly repeats The object returns to a given position after a fixed time interval A special kind of periodic motion occurs in mechanical systems
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 informationA plane autonomous system is a pair of simultaneous first-order differential equations,
Chapter 11 Phase-Plane Techniques 11.1 Plane Autonomous Systems A plane autonomous system is a pair of simultaneous first-order differential equations, ẋ = f(x, y), ẏ = g(x, y). This system has an equilibrium
More informationCHALMERS, GÖTEBORGS UNIVERSITET. EXAM for DYNAMICAL SYSTEMS. COURSE CODES: TIF 155, FIM770GU, PhD
CHALMERS, GÖTEBORGS UNIVERSITET EXAM for DYNAMICAL SYSTEMS COURSE CODES: TIF 155, FIM770GU, PhD Time: Place: Teachers: Allowed material: Not allowed: April 06, 2018, at 14 00 18 00 Johanneberg Kristian
More informationOld Math 330 Exams. David M. McClendon. Department of Mathematics Ferris State University
Old Math 330 Exams David M. McClendon Department of Mathematics Ferris State University Last updated to include exams from Fall 07 Contents Contents General information about these exams 3 Exams from Fall
More informationMath Ordinary Differential Equations
Math 411 - Ordinary Differential Equations Review Notes - 1 1 - Basic Theory A first order ordinary differential equation has the form x = f(t, x) (11) Here x = dx/dt Given an initial data x(t 0 ) = x
More informationLS.2 Homogeneous Linear Systems with Constant Coefficients
LS2 Homogeneous Linear Systems with Constant Coefficients Using matrices to solve linear systems The naive way to solve a linear system of ODE s with constant coefficients is by eliminating variables,
More informationMath 266: Phase Plane Portrait
Math 266: Phase Plane Portrait Long Jin Purdue, Spring 2018 Review: Phase line for an autonomous equation For a single autonomous equation y = f (y) we used a phase line to illustrate the equilibrium solutions
More information2D-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 informationReview: control, feedback, etc. Today s topic: state-space models of systems; linearization
Plan of the Lecture Review: control, feedback, etc Today s topic: state-space models of systems; linearization Goal: a general framework that encompasses all examples of interest Once we have mastered
More informationDifferential Equations and Modeling
Differential Equations and Modeling Preliminary Lecture Notes Adolfo J. Rumbos c Draft date: March 22, 2018 March 22, 2018 2 Contents 1 Preface 5 2 Introduction to Modeling 7 2.1 Constructing Models.........................
More informationRotational Kinetic Energy
Lecture 17, Chapter 10: Rotational Energy and Angular Momentum 1 Rotational Kinetic Energy Consider a rigid body rotating with an angular velocity ω about an axis. Clearly every point in the rigid body
More informationMATH 415, WEEKS 7 & 8: Conservative and Hamiltonian Systems, Non-linear Pendulum
MATH 415, WEEKS 7 & 8: Conservative and Hamiltonian Systems, Non-linear Pendulum Reconsider the following example from last week: dx dt = x y dy dt = x2 y. We were able to determine many qualitative features
More information!T = 2# T = 2! " The velocity and acceleration of the object are found by taking the first and second derivative of the position:
A pendulum swinging back and forth or a mass oscillating on a spring are two examples of (SHM.) SHM occurs any time the position of an object as a function of time can be represented by a sine wave. We
More informationLecture 9: Eigenvalues and Eigenvectors in Classical Mechanics (See Section 3.12 in Boas)
Lecture 9: Eigenvalues and Eigenvectors in Classical Mechanics (See Section 3 in Boas) As suggested in Lecture 8 the formalism of eigenvalues/eigenvectors has many applications in physics, especially in
More informationOutline. Calculus for the Life Sciences. What is a Differential Equation? Introduction. Lecture Notes Introduction to Differential Equa
Outline Calculus for the Life Sciences Lecture Notes to Differential Equations Joseph M. Mahaffy, jmahaffy@mail.sdsu.edu 1 Department of Mathematics and Statistics Dynamical Systems Group Computational
More information4 Insect outbreak model
4 Insect outbreak model In this lecture I will put to a good use all the mathematical machinery we discussed so far. Consider an insect population, which is subject to predation by birds. It is a very
More information1. < 0: the eigenvalues are real and have opposite signs; the fixed point is a saddle point
Solving a Linear System τ = trace(a) = a + d = λ 1 + λ 2 λ 1,2 = τ± = det(a) = ad bc = λ 1 λ 2 Classification of Fixed Points τ 2 4 1. < 0: the eigenvalues are real and have opposite signs; the fixed point
More informationQuestion: Total. Points:
MATH 308 May 23, 2011 Final Exam Name: ID: Question: 1 2 3 4 5 6 7 8 9 Total Points: 0 20 20 20 20 20 20 20 20 160 Score: There are 9 problems on 9 pages in this exam (not counting the cover sheet). Make
More informationPhysics 5A Final Review Solutions
Physics A Final Review Solutions Eric Reichwein Department of Physics University of California, Santa Cruz November 6, 0. A stone is dropped into the water from a tower 44.m above the ground. Another stone
More informationPHYS2330 Intermediate Mechanics Fall Final Exam Tuesday, 21 Dec 2010
Name: PHYS2330 Intermediate Mechanics Fall 2010 Final Exam Tuesday, 21 Dec 2010 This exam has two parts. Part I has 20 multiple choice questions, worth two points each. Part II consists of six relatively
More informationClassical Mechanics Comprehensive Exam Solution
Classical Mechanics Comprehensive Exam Solution January 31, 011, 1:00 pm 5:pm Solve the following six problems. In the following problems, e x, e y, and e z are unit vectors in the x, y, and z directions,
More informationVectors, matrices, eigenvalues and eigenvectors
Vectors, matrices, eigenvalues and eigenvectors 1 ( ) ( ) ( ) Scaling a vector: 0.5V 2 0.5 2 1 = 0.5 = = 1 0.5 1 0.5 ( ) ( ) ( ) ( ) Adding two vectors: V + W 2 1 2 + 1 3 = + = = 1 3 1 + 3 4 ( ) ( ) a
More informationApplications of Linear Higher-Order DEs
Week #4 : Applications of Linear Higher-Order DEs Goals: Solving Homogeneous DEs with Constant Coefficients - Second and Higher Order Applications - Damped Spring/Mass system Applications - Pendulum 1
More informationEigenvalues and Eigenvectors: An Introduction
Eigenvalues and Eigenvectors: An Introduction The eigenvalue problem is a problem of considerable theoretical interest and wide-ranging application. For example, this problem is crucial in solving systems
More informationLinear Algebra and ODEs review
Linear Algebra and ODEs review Ania A Baetica September 9, 015 1 Linear Algebra 11 Eigenvalues and eigenvectors Consider the square matrix A R n n (v, λ are an (eigenvector, eigenvalue pair of matrix A
More informationAIMS Exercise Set # 1
AIMS Exercise Set #. Determine the form of the single precision floating point arithmetic used in the computers at AIMS. What is the largest number that can be accurately represented? What is the smallest
More informationConsider a particle in 1D at position x(t), subject to a force F (x), so that mẍ = F (x). Define the kinetic energy to be.
Chapter 4 Energy and Stability 4.1 Energy in 1D Consider a particle in 1D at position x(t), subject to a force F (x), so that mẍ = F (x). Define the kinetic energy to be T = 1 2 mẋ2 and the potential energy
More information18. Linearization: the phugoid equation as example
79 18. Linearization: the phugoid equation as example Linearization is one of the most important and widely used mathematical terms in applications to Science and Engineering. In the context of Differential
More informationOscillatory Motion. Simple pendulum: linear Hooke s Law restoring force for small angular deviations. small angle approximation. Oscillatory solution
Oscillatory Motion Simple pendulum: linear Hooke s Law restoring force for small angular deviations d 2 θ dt 2 = g l θ small angle approximation θ l Oscillatory solution θ(t) =θ 0 sin(ωt + φ) F with characteristic
More informationMathematical Models. MATH 365 Ordinary Differential Equations. J. Robert Buchanan. Spring Department of Mathematics
Mathematical Models MATH 365 Ordinary Differential Equations J. Robert Buchanan Department of Mathematics Spring 2018 Ordinary Differential Equations The topic of ordinary differential equations (ODEs)
More informationMathematical Models. MATH 365 Ordinary Differential Equations. J. Robert Buchanan. Fall Department of Mathematics
Mathematical Models MATH 365 Ordinary Differential Equations J. Robert Buchanan Department of Mathematics Fall 2018 Ordinary Differential Equations The topic of ordinary differential equations (ODEs) is
More information6 EIGENVALUES AND EIGENVECTORS
6 EIGENVALUES AND EIGENVECTORS INTRODUCTION TO EIGENVALUES 61 Linear equations Ax = b come from steady state problems Eigenvalues have their greatest importance in dynamic problems The solution of du/dt
More informationOscillatory Motion. Simple pendulum: linear Hooke s Law restoring force for small angular deviations. Oscillatory solution
Oscillatory Motion Simple pendulum: linear Hooke s Law restoring force for small angular deviations d 2 θ dt 2 = g l θ θ l Oscillatory solution θ(t) =θ 0 sin(ωt + φ) F with characteristic angular frequency
More informationOrdinary Differential Equations
Ordinary Differential Equations Michael H. F. Wilkinson Institute for Mathematics and Computing Science University of Groningen The Netherlands December 2005 Overview What are Ordinary Differential Equations
More informationUnit 7: Oscillations
Text: Chapter 15 Unit 7: Oscillations NAME: Problems (p. 405-412) #1: 1, 7, 13, 17, 24, 26, 28, 32, 35 (simple harmonic motion, springs) #2: 45, 46, 49, 51, 75 (pendulums) Vocabulary: simple harmonic motion,
More informationChaotic motion. Phys 420/580 Lecture 10
Chaotic motion Phys 420/580 Lecture 10 Finite-difference equations Finite difference equation approximates a differential equation as an iterative map (x n+1,v n+1 )=M[(x n,v n )] Evolution from time t
More informationDifferential Equations 2280 Sample Midterm Exam 3 with Solutions Exam Date: 24 April 2015 at 12:50pm
Differential Equations 228 Sample Midterm Exam 3 with Solutions Exam Date: 24 April 25 at 2:5pm Instructions: This in-class exam is 5 minutes. No calculators, notes, tables or books. No answer check is
More informationClassification of Phase Portraits at Equilibria for u (t) = f( u(t))
Classification of Phase Portraits at Equilibria for u t = f ut Transfer of Local Linearized Phase Portrait Transfer of Local Linearized Stability How to Classify Linear Equilibria Justification of the
More informationPH 120 Project # 2: Pendulum and chaos
PH 120 Project # 2: Pendulum and chaos Due: Friday, January 16, 2004 In PH109, you studied a simple pendulum, which is an effectively massless rod of length l that is fixed at one end with a small mass
More informationAddendum: Lyapunov Exponent Calculation
Addendum: Lyapunov Exponent Calculation Experiment CP-Ly Equations of Motion The system phase u is written in vector notation as θ u = (1) φ The equation of motion can be expressed compactly using the
More informationPhysics 1A Lecture 10B
Physics 1A Lecture 10B "Sometimes the world puts a spin on life. When our equilibrium returns to us, we understand more because we've seen the whole picture. --Davis Barton Cross Products Another way to
More informationBasic Theory of Dynamical Systems
1 Basic Theory of Dynamical Systems Page 1 1.1 Introduction and Basic Examples Dynamical systems is concerned with both quantitative and qualitative properties of evolution equations, which are often ordinary
More information8 Example 1: The van der Pol oscillator (Strogatz Chapter 7)
8 Example 1: The van der Pol oscillator (Strogatz Chapter 7) So far we have seen some different possibilities of what can happen in two-dimensional systems (local and global attractors and bifurcations)
More informationChapter 2 Overview: Introduction to Limits and Derivatives
Chapter 2 Overview: Introduction to Limits and Derivatives In a later chapter, maximum and minimum points of a curve will be found both by calculator and algebraically. While the algebra of this process
More informationThe Liapunov Method for Determining Stability (DRAFT)
44 The Liapunov Method for Determining Stability (DRAFT) 44.1 The Liapunov Method, Naively Developed In the last chapter, we discussed describing trajectories of a 2 2 autonomous system x = F(x) as level
More informationLecture 20/Lab 21: Systems of Nonlinear ODEs
Lecture 20/Lab 21: Systems of Nonlinear ODEs MAR514 Geoffrey Cowles Department of Fisheries Oceanography School for Marine Science and Technology University of Massachusetts-Dartmouth Coupled ODEs: Species
More informationPHY411 Lecture notes Part 5
PHY411 Lecture notes Part 5 Alice Quillen January 27, 2016 Contents 0.1 Introduction.................................... 1 1 Symbolic Dynamics 2 1.1 The Shift map.................................. 3 1.2
More informationChapter 11 Vibrations and Waves
Chapter 11 Vibrations and Waves If an object vibrates or oscillates back and forth over the same path, each cycle taking the same amount of time, the motion is called periodic. The mass and spring system
More information7 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 informationImaginary. Axis. Real. Axis
Name ME6 Final. I certify that I upheld the Stanford Honor code during this exam Monday December 2, 25 3:3-6:3 p.m. ffl Print your name and sign the honor code statement ffl You may use your course notes,
More informationREVIEW OF DIFFERENTIAL CALCULUS
REVIEW OF DIFFERENTIAL CALCULUS DONU ARAPURA 1. Limits and continuity To simplify the statements, we will often stick to two variables, but everything holds with any number of variables. Let f(x, y) be
More informationChapter 6 Nonlinear Systems and Phenomena. Friday, November 2, 12
Chapter 6 Nonlinear Systems and Phenomena 6.1 Stability and the Phase Plane We now move to nonlinear systems Begin with the first-order system for x(t) d dt x = f(x,t), x(0) = x 0 In particular, consider
More informationr CM = ir im i i m i m i v i (2) P = i
Physics 121 Test 3 study guide Thisisintendedtobeastudyguideforyourthirdtest, whichcoverschapters 9, 10, 12, and 13. Note that chapter 10 was also covered in test 2 without section 10.7 (elastic collisions),
More informationFigure 5.16 Compound pendulum: (a) At rest in equilibrium, (b) General position with coordinate θ, Freebody
Lecture 27. THE COMPOUND PENDULUM Figure 5.16 Compound pendulum: (a) At rest in equilibrium, (b) General position with coordinate θ, Freebody diagram The term compound is used to distinguish the present
More informationElementary Differential Equations
Elementary Differential Equations George Voutsadakis 1 1 Mathematics and Computer Science Lake Superior State University LSSU Math 310 George Voutsadakis (LSSU) Differential Equations January 2014 1 /
More informationOscillatory Motion SHM
Chapter 15 Oscillatory Motion SHM Dr. Armen Kocharian Periodic Motion Periodic motion is motion of an object that regularly repeats The object returns to a given position after a fixed time interval A
More informationSecond-Order Linear ODEs
Chap. 2 Second-Order Linear ODEs Sec. 2.1 Homogeneous Linear ODEs of Second Order On pp. 45-46 we extend concepts defined in Chap. 1, notably solution and homogeneous and nonhomogeneous, to second-order
More informationFIRST-ORDER SYSTEMS OF ORDINARY DIFFERENTIAL EQUATIONS III: Autonomous Planar Systems David Levermore Department of Mathematics University of Maryland
FIRST-ORDER SYSTEMS OF ORDINARY DIFFERENTIAL EQUATIONS III: Autonomous Planar Systems David Levermore Department of Mathematics University of Maryland 4 May 2012 Because the presentation of this material
More information154 Chapter 9 Hints, Answers, and Solutions The particular trajectories are highlighted in the phase portraits below.
54 Chapter 9 Hints, Answers, and Solutions 9. The Phase Plane 9.. 4. The particular trajectories are highlighted in the phase portraits below... 3. 4. 9..5. Shown below is one possibility with x(t) and
More informationFor a rigid body that is constrained to rotate about a fixed axis, the gravitational torque about the axis is
Experiment 14 The Physical Pendulum The period of oscillation of a physical pendulum is found to a high degree of accuracy by two methods: theory and experiment. The values are then compared. Theory For
More informationLecture 7 - Separable Equations
Lecture 7 - Separable Equations Separable equations is a very special type of differential equations where you can separate the terms involving only y on one side of the equation and terms involving only
More informationUpdated 2013 (Mathematica Version) M1.1. Lab M1: The Simple Pendulum
Updated 2013 (Mathematica Version) M1.1 Introduction. Lab M1: The Simple Pendulum The simple pendulum is a favorite introductory exercise because Galileo's experiments on pendulums in the early 1600s are
More informationChaotic motion. Phys 750 Lecture 9
Chaotic motion Phys 750 Lecture 9 Finite-difference equations Finite difference equation approximates a differential equation as an iterative map (x n+1,v n+1 )=M[(x n,v n )] Evolution from time t =0to
More informationAP Physics. Harmonic Motion. Multiple Choice. Test E
AP Physics Harmonic Motion Multiple Choice Test E A 0.10-Kg block is attached to a spring, initially unstretched, of force constant k = 40 N m as shown below. The block is released from rest at t = 0 sec.
More informationAPPM 2360: Final Exam 10:30am 1:00pm, May 6, 2015.
APPM 23: Final Exam :3am :pm, May, 25. ON THE FRONT OF YOUR BLUEBOOK write: ) your name, 2) your student ID number, 3) lecture section, 4) your instructor s name, and 5) a grading table for eight questions.
More informationAPPM 2360 Lab #3: The Predator Prey Model
APPM 2360 Lab #3: The Predator Prey Model 1 Instructions Labs may be done in groups of 3 or less. One report must be turned in for each group and must be in PDF format. Labs must include each student s:
More informationLecture Notes. W. Tang. School of Mathematical & Statistical Sciences Arizona State University
Applied Math / Dynamical Systems Workshop Lecture Notes W. Tang School of Mathematical & Statistical Sciences Arizona State University May. 23 rd -25 th, 2011 Tang (ASU) AMDS Workshop Notes May. 23 rd
More informationT = 1 2 mv2. For a system of N particles, the kinetic energy of the system is obtained by adding up the individual kinetic energies: T = 2 m ivi 2.
Physics 3550, Fall 2011 Energy. Relevant Sections in Text: Chapter 4 Kinetic Energy. I probably don t need to impress upon you the importance of the notion of energy in analyzing/understanding physical
More informationChapter 5: Limits and Derivatives
Chapter 5: Limits and Derivatives Chapter 5 Overview: Introduction to Limits and Derivatives In a later chapter, maximum and minimum points of a curve will be found both by calculator and algebraically.
More informationMeasuring the Universal Gravitational Constant, G
Measuring the Universal Gravitational Constant, G Introduction: The universal law of gravitation states that everything in the universe is attracted to everything else. It seems reasonable that everything
More informationPhysics 351, Spring 2017, Homework #3. Due at start of class, Friday, February 3, 2017
Physics 351, Spring 2017, Homework #3. Due at start of class, Friday, February 3, 2017 Course info is at positron.hep.upenn.edu/p351 When you finish this homework, remember to visit the feedback page at
More information1. The accumulated net change function or area-so-far function
Name: Section: Names of collaborators: Main Points: 1. The accumulated net change function ( area-so-far function) 2. Connection to antiderivative functions: the Fundamental Theorem of Calculus 3. Evaluating
More information