We here collect a few numerical tests, in order to put into evidence the potentialities of HBVMs [4, 6, 7].
|
|
- Jacob Blake
- 5 years ago
- Views:
Transcription
1 Chaper Numerical Tess We here collec a few numerical ess, in order o pu ino evidence he poenialiies of HBVMs [4, 6, 7]. Tes problem Le us consider he problem characerized by he polynomial Hamilonian (4.) in [9], H(p, q) = p3 3 p + q6 3 + q4 4 q , (.) havingdegree ν = 6, saringa he iniial poin y (q(), p()) T = (, ) T, so ha H(y ) =. For such a problem, in [9] i has been experienced a numerical drif in he discree Hamilonian, when usinghe fourh-order Lobao IIIA mehod wih sepsize h =.6, as confirmed by he plo in Figure.. When using he fourh-order Gauss-Legendre mehod he drif disappears, even hough he Hamilonian is no exacly preserved alonghe discree soluion, as is confirmed by he plo in Figure.. On he oher hand, by usinghe fourh-order HBVM(6,) wih he same sepsize, he Hamilonian urns ou o be preserved up o machine precision, as shown in Figure.3, since such mehod exacly preserves polynomial Hamilonians of degree up o 6. In such a case, accordingo he las iem in Remark7, he numerical soluions obained by usinghe Lobao nodes {c =, c,..., c 6 = } or he Gauss-Legendre nodes {c,..., c 6 } are he same. The fourh-order convergence of he mehod is numerically verified by he resuls lised in Table.. 5
2 6 CHAPTER. NUMERICAL TESTS x 6 H Figure.:Fourh-order Lobao IIIA mehod, h =.6, problem (.):drif in he Hamilonian. 3.5 x H Figure.:Fourh-order Gauss-Legendre mehod, h =.6, problem (.):H 6. x 6 H Figure.3:Fourh-order HBVM(6,) mehod, h =.6, problem (.):H 6.
3 7 Tes problem The second es problem, havinga highly oscillaingsoluion, is he Fermi-Pasa- Ulam problem (see [, Secion I.5.]), modellinga chain of m mass poins conneced wih alernaingsof nonlinear and siff linear springs, and fixed a he end poins. The variables q,..., q m sand for he displacemens of he mass poins, and p i = q i for heir velociies. The correspondinghamilonian, represeninghe oal energy, is H(p, q) = m ( p i + pi) ω + i= 4 m m (q i q i ) + (q i+ q i ) 4, (.) i= i= wih q = q m+ =. In our simulaion we have used he followingvalues: m = 3, ω = 5, and saringvecor p i =, q i = (i )/, i =,..., 6. In such a case, he Hamilonian funcion is a polynomial of degree 4, so ha he fourh-order HBVM(4,) mehod, eiher when usinghe Lobao nodes or he Gauss-Legendre nodes, is able o exacly preserve he Hamilonian, as confirmed by he plo in Figure.6, obained wih sepsize h =.5. Conversely, by using he same sepsize, boh he fourh-order Lobao IIIA and Gauss-Legendre mehods provide only an approximae conservaion of he Hamilonian, as shown in he plos in Figures.4 and.5, respecively. The fourh-order convergence of he HBVM(4,) mehod is numerically verified by he resuls lised in Table..
4 8 CHAPTER. NUMERICAL TESTS x H H Figure.4:Fourh-order Lobao IIIA mehod, h =.5, problem (.): H H 3. 4 x 4 H H Figure.5:Fourh-order Gauss-Legendre mehod, h =.5, problem (.): H H 3..5 x 4.5 H H Figure.6:Fourh-order HBVM(4,) mehod, h =.5, problem (.): H H 4.
5 9 Tes problem 3(non-polynomial Hamilonian) In he previous examples, he Hamilonian funcion was a polynomial. Neverheless, as observed above, also in his case HBVM(k,s) are expeced o produce a pracical conservaion of he energy when applied o sysems defined by a nonpolynomial Hamilonian funcion ha can be locally well approximaed by a polynomial. As an example, we consider he moion of a charged paricle in a magneic field wih Bio-Savar poenial. I is defined by he Hamilonian [6] H(x, y, z, ẋ, ẏ, ż) = (.3) [ ( ẋ α x ) + (ẏ α yϱ ) ] + (ż + α log(ϱ)), m ϱ wih ϱ = x + y, α = e B, m is he paricle mass, e is is charge, and B is he magneic field inensiy. We have used he values wih saringpoin m =, e =, B =, x =.5, y =, z =, ẋ =., ẏ =.3, ż =. By usinghe fourh-order Lobao IIIA mehod, wih sepsize h =., a drif is again experienced in he numerical soluion, as is shown in Figure.7. By usinghe fourh-order Gauss-Legendre mehod wih he same sepsize, he drif disappears even hough, as shown in Figure.8, he value of he Hamilonian is preserved wihin an error of he order of 3. On he oher hand, when using he HBVM(6,) mehod wih he same sepsize, he error in he Hamilonian decreases o an order of 5 (see Figure.9), hus givinga pracical conservaion. Finally, in Table.4 we lis he maximum absolue difference beween he numerical soluions over 3 inegraion seps, compued by he HBVM(k, ) mehods based on Lobao abscissae and on Gauss-Legendre abscissae, as k grows, wih sepsize h =.. We observe ha he difference ends o, as k increases. Finally, also in his case, one verifies a fourh-order convergence, as he resuls lised in Table.3show. This kind of moion causes he well known phenomenon of aurora borealis.
6 CHAPTER. NUMERICAL TESTS 6 x 3 4 H H Figure.7:Fourh-order Lobao IIIA mehod, h =., problem (.3):drif in he Hamilonian. 5 x 3 H H Figure.8:Fourh-order Gauss-Legendre mehod, h =., problem (.3): H H 3. x 5 H H Figure.9:Fourh-order HBVM(6,) mehod, h =., problem (.3): H H 5.
7 Table.:Numerical order of convergence for he HBVM(6,) mehod, problem (.). h error order Table.:Numerical order of convergence for he HBVM(4,) mehod, problem (.). h error order Table.3:Numerical order of convergence for he HBVM(6,) mehod, problem (.3). h error order Table.4:Maximum difference beween he numerical soluions obained hrough he fourh-order HBVM(k, ) mehods based on Lobao abscissae and Gauss- Legendre abscissae for increasingvalues of k, problem (.3), 3 seps wih sepsize h =.. k h =
8 CHAPTER. NUMERICAL TESTS Tes problem 4(Sinikovproblem) The main problem in Celesial Mechanics is he so called N-body problem, i.e. o describe he moion of N poin paricles of posiive mass movingunder Newon s law of graviaion when we know heir posiions and velociies a a given ime. This problem is described by he Hamilonian funcion: H(q, p) = N i= p i m i G N i= i m j m i, (.4) q i q j j= where q i is he posiion of he ih paricle, wih mass m i, and p i is is momenum. The Sinikov problem is a paricular configuraion of he 3-body dynamics (see, e.g., [3]). In his problem wo bodies of equal mass (primaries) revolve abou heir cener of mass, here assumed a he origin, in ellipic orbis in he xyplane. A hird, and much smaller body (planeoid), is placed on he z-axis wih iniial velociy parallel o his axis as well. The hird body is small enough ha he wo body dynamics of he primaries is no desroyed. Then, he moion of he hird body will be resriced o he z-axis and oscillaingaround he origin bu no necessarily periodic. In fac his problem has been shown o exhibi a chaoic behavior when he eccenriciy of he orbis of he primaries exceeds a criical value ha, for he daa se we have used, is ē.75 (see Figure.). We have solved he problem defined by he Hamilonian funcion (.4) by he Gauss mehod of order 4 (i.e., HBVM(,) a Gaussian nodes) and by HBVM(8,) a 8 Gaussian nodes (order 4, fundamenal and 6 silen sages), wih he followingse of parameers in (.4): N G m m m 3 e d h max where e is he eccenriciy, d is he disance of he apocenres of he primaries (poins a which he wo bodies are he furhes), h is he used ime-sep, and [, max] is he ime inegraion inerval. The eccenriciy e and he disance d may be used o define he iniial condiion [q, p ] (see [3] for he deails): q = [ 5,,, 5,,,,, 9 ] T, p = [,,,,,,,, ]T. Firs of all, we consider he wo picures in Figure. reporinghe relaive errors in he Hamilonian funcion and in he angular momenum evaluaed alonghe numerical soluions compued by he wo mehods. We know ha he
9 HBVM(8,) precisely conserves Hamilonian polynomial funcions of degree a mos 8. This accuracy is high enough o guaranee ha he nonlinear Hamilonian funcion (.4) is as well conserved up o he machine precision (see he upper picure):from a geomerical poin of view his means ha a local approximaion of he level curves of (.4) by a polynomial of degree 8 leads o a negligible error. The Gauss mehod exhibis a cerain error in he Hamilonian funcion while, beinghis formula symplecic, i precisely conserves he angular momenum, as is confirmed by lookinga he down picure of Figure.. The error in he numerical angular momenum associaed wih he HBVM(8,) undergoes some bounded periodic-like oscillaions. Figures. and.3show he numerical soluion compued by he Gauss mehod and HBVM(8,), respecively. Since he mehods leave he xy-plane invarian for he moion of he primaries and he z-axis invarian for he moion of he planeoid, we have jus repored he moion of he primaries in he xy-phase plane (upper picures) and he space-ime diagram of he planeoid (down picure). We observe ha, for he Gauss mehod, he orbis of he primaries are irregular in characer so ha he hird body, afer performingsome oscillaions around he origin, will evenually escape he sysem (see he down picure of Figure.). On he conrary (see he upper picure of Figure.3), he HBVM(8,) mehod generaes a quie regular phase porrai. Due o he large sepsize h used, a sham roaion of he xy-plane appears which, however, does no desroy he global symmery of he dynamics, as esified by he bounded oscillaions of he planeoid (down picure of Figure.3) which lookvery similar o he reference ones in Figure.. This aspec is also confirmed by he picures in Figure.4 displayinghe disance of he primaries as a funcion of he ime. We see ha he disance of he apocenres (correspondingo he maxima in he plos), as he wo bodies wheel around he origin, are preserved by he HBVM(8,) (down picure) while he same is no rue for he Gauss mehod (upper picure). 3
10 4 CHAPTER. NUMERICAL TESTS Figure.:The upper picure displays he configuraion of 3-bodies in he Sinikov problem. To an eccenriciy of he orbis of he primaries e =.75, here correspond bounded chaoic oscillaions of he planeoid as is argued by looking a he space-ime diagram in he down picure.
11 5 5 HBVM(8,) Gauss HBVM(8,) Gauss Figure.:Upper picure:relaive error H(y n ) H(y ) / H(y ) of he Hamilonian funcion evaluaed alonghe numerical soluion of he HBVM(8,) and he Gauss mehod. Down picure:relaive error M(y n ) M(y ) / M(y ) of he angular momenum evaluaed alonghe numerical soluion of he HBVM(8,) and he Gauss mehod.
12 6 CHAPTER. NUMERICAL TESTS.5 Trajecory of m Trajecory of m Figure.:The Sinikov problem solved by he Gauss mehod of order 4, wih sepsize h =.5, in he ime inerval [, 5]. The rajecories of he primaries in he xy-plane (upper picure) exhibi a very irregular behavior which causes he planeoid o evenually escape he sysem, as illusraed by he space-ime diagram in he down picure.
13 7.5 Trajecory of m Trajecory of m Figure.3:The Sinikov problem solved by he HBVM(8,) mehod (order 4), wih sepsize h =.5, in he ime inerval [, 5]. Upper picure:he rajecories of he primaries are ellipse shape. The discreizaion inroduces a ficiious uniform roaion of he xy-plane which however does no aler he global symmery of he sysem. Down picure:he space-ime diagram of he planeoid on he z-axis displayed (for clearness) on he ime inerval [, 35] shows ha, alhough a large value of he sepsize h has been used, he overall behavior of he dynamics is well reproduced (compare wih he down picure in Figure.).
14 8 CHAPTER. NUMERICAL TESTS Figure.4:Disance beween he wo primaries as a funcion of he ime, relaed o he numerical soluions generaed by he Gauss mehod (upper picure) and HBVM(8,) (down picure). The maxima correspond o he disance of apocenres. These are conserved by HBVM(8,) while he Gauss mehod inroduces pachy oscillaions ha desroy he overall symmery of he sysem.
Physics 235 Chapter 2. Chapter 2 Newtonian Mechanics Single Particle
Chaper 2 Newonian Mechanics Single Paricle In his Chaper we will review wha Newon s laws of mechanics ell us abou he moion of a single paricle. Newon s laws are only valid in suiable reference frames,
More informationMath Week 14 April 16-20: sections first order systems of linear differential equations; 7.4 mass-spring systems.
Mah 2250-004 Week 4 April 6-20 secions 7.-7.3 firs order sysems of linear differenial equaions; 7.4 mass-spring sysems. Mon Apr 6 7.-7.2 Sysems of differenial equaions (7.), and he vecor Calculus we need
More informationSome Basic Information about M-S-D Systems
Some Basic Informaion abou M-S-D Sysems 1 Inroducion We wan o give some summary of he facs concerning unforced (homogeneous) and forced (non-homogeneous) models for linear oscillaors governed by second-order,
More informationMEI Mechanics 1 General motion. Section 1: Using calculus
Soluions o Exercise MEI Mechanics General moion Secion : Using calculus. s 4 v a 6 4 4 When =, v 4 a 6 4 6. (i) When = 0, s = -, so he iniial displacemen = - m. s v 4 When = 0, v = so he iniial velociy
More informationLecture 2-1 Kinematics in One Dimension Displacement, Velocity and Acceleration Everything in the world is moving. Nothing stays still.
Lecure - Kinemaics in One Dimension Displacemen, Velociy and Acceleraion Everyhing in he world is moving. Nohing says sill. Moion occurs a all scales of he universe, saring from he moion of elecrons in
More informationVehicle Arrival Models : Headway
Chaper 12 Vehicle Arrival Models : Headway 12.1 Inroducion Modelling arrival of vehicle a secion of road is an imporan sep in raffic flow modelling. I has imporan applicaion in raffic flow simulaion where
More informationWEEK-3 Recitation PHYS 131. of the projectile s velocity remains constant throughout the motion, since the acceleration a x
WEEK-3 Reciaion PHYS 131 Ch. 3: FOC 1, 3, 4, 6, 14. Problems 9, 37, 41 & 71 and Ch. 4: FOC 1, 3, 5, 8. Problems 3, 5 & 16. Feb 8, 018 Ch. 3: FOC 1, 3, 4, 6, 14. 1. (a) The horizonal componen of he projecile
More informationEXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE
Version April 30, 2004.Submied o CTU Repors. EXPLICIT TIME INTEGRATORS FOR NONLINEAR DYNAMICS DERIVED FROM THE MIDPOINT RULE Per Krysl Universiy of California, San Diego La Jolla, California 92093-0085,
More informationDifferential Equations
Mah 21 (Fall 29) Differenial Equaions Soluion #3 1. Find he paricular soluion of he following differenial equaion by variaion of parameer (a) y + y = csc (b) 2 y + y y = ln, > Soluion: (a) The corresponding
More informationRobotics I. April 11, The kinematics of a 3R spatial robot is specified by the Denavit-Hartenberg parameters in Tab. 1.
Roboics I April 11, 017 Exercise 1 he kinemaics of a 3R spaial robo is specified by he Denavi-Harenberg parameers in ab 1 i α i d i a i θ i 1 π/ L 1 0 1 0 0 L 3 0 0 L 3 3 able 1: able of DH parameers of
More informationSolutions from Chapter 9.1 and 9.2
Soluions from Chaper 9 and 92 Secion 9 Problem # This basically boils down o an exercise in he chain rule from calculus We are looking for soluions of he form: u( x) = f( k x c) where k x R 3 and k is
More information4.5 Constant Acceleration
4.5 Consan Acceleraion v() v() = v 0 + a a() a a() = a v 0 Area = a (a) (b) Figure 4.8 Consan acceleraion: (a) velociy, (b) acceleraion When he x -componen of he velociy is a linear funcion (Figure 4.8(a)),
More informationThe Contradiction within Equations of Motion with Constant Acceleration
The Conradicion wihin Equaions of Moion wih Consan Acceleraion Louai Hassan Elzein Basheir (Daed: July 7, 0 This paper is prepared o demonsrae he violaion of rules of mahemaics in he algebraic derivaion
More informationWall. x(t) f(t) x(t = 0) = x 0, t=0. which describes the motion of the mass in absence of any external forcing.
MECHANICS APPLICATIONS OF SECOND-ORDER ODES 7 Mechanics applicaions of second-order ODEs Second-order linear ODEs wih consan coefficiens arise in many physical applicaions. One physical sysems whose behaviour
More informationPredator - Prey Model Trajectories and the nonlinear conservation law
Predaor - Prey Model Trajecories and he nonlinear conservaion law James K. Peerson Deparmen of Biological Sciences and Deparmen of Mahemaical Sciences Clemson Universiy Ocober 28, 213 Ouline Drawing Trajecories
More information!!"#"$%&#'()!"#&'(*%)+,&',-)./0)1-*23)
"#"$%&#'()"#&'(*%)+,&',-)./)1-*) #$%&'()*+,&',-.%,/)*+,-&1*#$)()5*6$+$%*,7&*-'-&1*(,-&*6&,7.$%$+*&%'(*8$&',-,%'-&1*(,-&*6&,79*(&,%: ;..,*&1$&$.$%&'()*1$$.,'&',-9*(&,%)?%*,('&5
More informationQ.1 Define work and its unit?
CHP # 6 ORK AND ENERGY Q.1 Define work and is uni? A. ORK I can be define as when we applied a force on a body and he body covers a disance in he direcion of force, hen we say ha work is done. I is a scalar
More information1. VELOCITY AND ACCELERATION
1. VELOCITY AND ACCELERATION 1.1 Kinemaics Equaions s = u + 1 a and s = v 1 a s = 1 (u + v) v = u + as 1. Displacemen-Time Graph Gradien = speed 1.3 Velociy-Time Graph Gradien = acceleraion Area under
More informationHamiltonian Boundary Value Methods (Energy Preserving Discrete Line Integral Methods) 1 2
European Sociey of Compuaional Mehods in Sciences and Engineering (ESCMSE) Journal of Numerical Analysis, Indusrial and Applied Mahemaics (JNAIAM) vol. 5, no. -2, 2, pp. 7-37 ISSN 79 84 Hamilonian Boundary
More informationTwo Coupled Oscillators / Normal Modes
Lecure 3 Phys 3750 Two Coupled Oscillaors / Normal Modes Overview and Moivaion: Today we ake a small, bu significan, sep owards wave moion. We will no ye observe waves, bu his sep is imporan in is own
More informationTraveling Waves. Chapter Introduction
Chaper 4 Traveling Waves 4.1 Inroducion To dae, we have considered oscillaions, i.e., periodic, ofen harmonic, variaions of a physical characerisic of a sysem. The sysem a one ime is indisinguishable from
More informationModule 2 F c i k c s la l w a s o s f dif di fusi s o i n
Module Fick s laws of diffusion Fick s laws of diffusion and hin film soluion Adolf Fick (1855) proposed: d J α d d d J (mole/m s) flu (m /s) diffusion coefficien and (mole/m 3 ) concenraion of ions, aoms
More informationElectrical Circuits. 1. Circuit Laws. Tools Used in Lab 13 Series Circuits Damped Vibrations: Energy Van der Pol Circuit
V() R L C 513 Elecrical Circuis Tools Used in Lab 13 Series Circuis Damped Vibraions: Energy Van der Pol Circui A series circui wih an inducor, resisor, and capacior can be represened by Lq + Rq + 1, a
More informationwhere the coordinate X (t) describes the system motion. X has its origin at the system static equilibrium position (SEP).
Appendix A: Conservaion of Mechanical Energy = Conservaion of Linear Momenum Consider he moion of a nd order mechanical sysem comprised of he fundamenal mechanical elemens: ineria or mass (M), siffness
More informationk 1 k 2 x (1) x 2 = k 1 x 1 = k 2 k 1 +k 2 x (2) x k series x (3) k 2 x 2 = k 1 k 2 = k 1+k 2 = 1 k k 2 k series
Final Review A Puzzle... Consider wo massless springs wih spring consans k 1 and k and he same equilibrium lengh. 1. If hese springs ac on a mass m in parallel, hey would be equivalen o a single spring
More informationWeek 1 Lecture 2 Problems 2, 5. What if something oscillates with no obvious spring? What is ω? (problem set problem)
Week 1 Lecure Problems, 5 Wha if somehing oscillaes wih no obvious spring? Wha is ω? (problem se problem) Sar wih Try and ge o SHM form E. Full beer can in lake, oscillaing F = m & = ge rearrange: F =
More information18 Biological models with discrete time
8 Biological models wih discree ime The mos imporan applicaions, however, may be pedagogical. The elegan body of mahemaical heory peraining o linear sysems (Fourier analysis, orhogonal funcions, and so
More informationSingle and Double Pendulum Models
Single and Double Pendulum Models Mah 596 Projec Summary Spring 2016 Jarod Har 1 Overview Differen ypes of pendulums are used o model many phenomena in various disciplines. In paricular, single and double
More informationModule 3: The Damped Oscillator-II Lecture 3: The Damped Oscillator-II
Module 3: The Damped Oscillaor-II Lecure 3: The Damped Oscillaor-II 3. Over-damped Oscillaions. This refers o he siuaion where β > ω (3.) The wo roos are and α = β + α 2 = β β 2 ω 2 = (3.2) β 2 ω 2 = 2
More informationKinematics and kinematic functions
Kinemaics and kinemaic funcions Kinemaics deals wih he sudy of four funcions (called kinemaic funcions or KFs) ha mahemaically ransform join variables ino caresian variables and vice versa Direc Posiion
More informationTHE MYSTERY OF STOCHASTIC MECHANICS. Edward Nelson Department of Mathematics Princeton University
THE MYSTERY OF STOCHASTIC MECHANICS Edward Nelson Deparmen of Mahemaics Princeon Universiy 1 Classical Hamilon-Jacobi heory N paricles of various masses on a Euclidean space. Incorporae he masses in he
More informationHOMEWORK # 2: MATH 211, SPRING Note: This is the last solution set where I will describe the MATLAB I used to make my pictures.
HOMEWORK # 2: MATH 2, SPRING 25 TJ HITCHMAN Noe: This is he las soluion se where I will describe he MATLAB I used o make my picures.. Exercises from he ex.. Chaper 2.. Problem 6. We are o show ha y() =
More informationNon-uniform circular motion *
OpenSax-CNX module: m14020 1 Non-uniform circular moion * Sunil Kumar Singh This work is produced by OpenSax-CNX and licensed under he Creaive Commons Aribuion License 2.0 Wha do we mean by non-uniform
More informationKINEMATICS IN ONE DIMENSION
KINEMATICS IN ONE DIMENSION PREVIEW Kinemaics is he sudy of how hings move how far (disance and displacemen), how fas (speed and velociy), and how fas ha how fas changes (acceleraion). We say ha an objec
More informationLecture 4 Kinetics of a particle Part 3: Impulse and Momentum
MEE Engineering Mechanics II Lecure 4 Lecure 4 Kineics of a paricle Par 3: Impulse and Momenum Linear impulse and momenum Saring from he equaion of moion for a paricle of mass m which is subjeced o an
More informationSecond-Order Differential Equations
WWW Problems and Soluions 3.1 Chaper 3 Second-Order Differenial Equaions Secion 3.1 Springs: Linear and Nonlinear Models www m Problem 3. (NonlinearSprings). A bod of mass m is aached o a wall b means
More informationON THE BEAT PHENOMENON IN COUPLED SYSTEMS
8 h ASCE Specialy Conference on Probabilisic Mechanics and Srucural Reliabiliy PMC-38 ON THE BEAT PHENOMENON IN COUPLED SYSTEMS S. K. Yalla, Suden Member ASCE and A. Kareem, M. ASCE NaHaz Modeling Laboraory,
More informationAnnouncements: Warm-up Exercise:
Fri Apr 13 7.1 Sysems of differenial equaions - o model muli-componen sysems via comparmenal analysis hp//en.wikipedia.org/wiki/muli-comparmen_model Announcemens Warm-up Exercise Here's a relaively simple
More information3.6 Derivatives as Rates of Change
3.6 Derivaives as Raes of Change Problem 1 John is walking along a sraigh pah. His posiion a he ime >0 is given by s = f(). He sars a =0from his house (f(0) = 0) and he graph of f is given below. (a) Describe
More informationLinear Surface Gravity Waves 3., Dispersion, Group Velocity, and Energy Propagation
Chaper 4 Linear Surface Graviy Waves 3., Dispersion, Group Velociy, and Energy Propagaion 4. Descripion In many aspecs of wave evoluion, he concep of group velociy plays a cenral role. Mos people now i
More information15. Vector Valued Functions
1. Vecor Valued Funcions Up o his poin, we have presened vecors wih consan componens, for example, 1, and,,4. However, we can allow he componens of a vecor o be funcions of a common variable. For example,
More informationInterpretation of special relativity as applied to earth-centered locally inertial
Inerpreaion of special relaiviy as applied o earh-cenered locally inerial coordinae sysems in lobal osiioning Sysem saellie experimens Masanori Sao Honda Elecronics Co., Ld., Oyamazuka, Oiwa-cho, Toyohashi,
More informationCH.7. PLANE LINEAR ELASTICITY. Continuum Mechanics Course (MMC) - ETSECCPB - UPC
CH.7. PLANE LINEAR ELASTICITY Coninuum Mechanics Course (MMC) - ETSECCPB - UPC Overview Plane Linear Elasici Theor Plane Sress Simplifing Hpohesis Srain Field Consiuive Equaion Displacemen Field The Linear
More informationdt = C exp (3 ln t 4 ). t 4 W = C exp ( ln(4 t) 3) = C(4 t) 3.
Mah Rahman Exam Review Soluions () Consider he IVP: ( 4)y 3y + 4y = ; y(3) = 0, y (3) =. (a) Please deermine he longes inerval for which he IVP is guaraneed o have a unique soluion. Soluion: The disconinuiies
More informationStructural Dynamics and Earthquake Engineering
Srucural Dynamics and Earhquae Engineering Course 1 Inroducion. Single degree of freedom sysems: Equaions of moion, problem saemen, soluion mehods. Course noes are available for download a hp://www.c.up.ro/users/aurelsraan/
More informationTHE 2-BODY PROBLEM. FIGURE 1. A pair of ellipses sharing a common focus. (c,b) c+a ROBERT J. VANDERBEI
THE 2-BODY PROBLEM ROBERT J. VANDERBEI ABSTRACT. In his shor noe, we show ha a pair of ellipses wih a common focus is a soluion o he 2-body problem. INTRODUCTION. Solving he 2-body problem from scrach
More informationClass Meeting # 10: Introduction to the Wave Equation
MATH 8.5 COURSE NOTES - CLASS MEETING # 0 8.5 Inroducion o PDEs, Fall 0 Professor: Jared Speck Class Meeing # 0: Inroducion o he Wave Equaion. Wha is he wave equaion? The sandard wave equaion for a funcion
More informationFinal Spring 2007
.615 Final Spring 7 Overview The purpose of he final exam is o calculae he MHD β limi in a high-bea oroidal okamak agains he dangerous n = 1 exernal ballooning-kink mode. Effecively, his corresponds o
More informationFrom Particles to Rigid Bodies
Rigid Body Dynamics From Paricles o Rigid Bodies Paricles No roaions Linear velociy v only Rigid bodies Body roaions Linear velociy v Angular velociy ω Rigid Bodies Rigid bodies have boh a posiion and
More informationBasilio Bona ROBOTICA 03CFIOR 1
Indusrial Robos Kinemaics 1 Kinemaics and kinemaic funcions Kinemaics deals wih he sudy of four funcions (called kinemaic funcions or KFs) ha mahemaically ransform join variables ino caresian variables
More informationLinear Response Theory: The connection between QFT and experiments
Phys540.nb 39 3 Linear Response Theory: The connecion beween QFT and experimens 3.1. Basic conceps and ideas Q: How do we measure he conduciviy of a meal? A: we firs inroduce a weak elecric field E, and
More informationSPH3U: Projectiles. Recorder: Manager: Speaker:
SPH3U: Projeciles Now i s ime o use our new skills o analyze he moion of a golf ball ha was ossed hrough he air. Le s find ou wha is special abou he moion of a projecile. Recorder: Manager: Speaker: 0
More informationGuest Lectures for Dr. MacFarlane s EE3350 Part Deux
Gues Lecures for Dr. MacFarlane s EE3350 Par Deux Michael Plane Mon., 08-30-2010 Wrie name in corner. Poin ou his is a review, so I will go faser. Remind hem o go lisen o online lecure abou geing an A
More informationMotion along a Straight Line
chaper 2 Moion along a Sraigh Line verage speed and average velociy (Secion 2.2) 1. Velociy versus speed Cone in he ebook: fer Eample 2. Insananeous velociy and insananeous acceleraion (Secions 2.3, 2.4)
More informationNotes on Kalman Filtering
Noes on Kalman Filering Brian Borchers and Rick Aser November 7, Inroducion Daa Assimilaion is he problem of merging model predicions wih acual measuremens of a sysem o produce an opimal esimae of he curren
More informationOrdinary dierential equations
Chaper 5 Ordinary dierenial equaions Conens 5.1 Iniial value problem........................... 31 5. Forward Euler's mehod......................... 3 5.3 Runge-Kua mehods.......................... 36
More informationSection 7.4 Modeling Changing Amplitude and Midline
488 Chaper 7 Secion 7.4 Modeling Changing Ampliude and Midline While sinusoidal funcions can model a variey of behaviors, i is ofen necessary o combine sinusoidal funcions wih linear and exponenial curves
More informationChapter 3 Boundary Value Problem
Chaper 3 Boundary Value Problem A boundary value problem (BVP) is a problem, ypically an ODE or a PDE, which has values assigned on he physical boundary of he domain in which he problem is specified. Le
More informationSection 3.8, Mechanical and Electrical Vibrations
Secion 3.8, Mechanical and Elecrical Vibraions Mechanical Unis in he U.S. Cusomary and Meric Sysems Disance Mass Time Force g (Earh) Uni U.S. Cusomary MKS Sysem CGS Sysem fee f slugs seconds sec pounds
More informationMatlab and Python programming: how to get started
Malab and Pyhon programming: how o ge sared Equipping readers he skills o wrie programs o explore complex sysems and discover ineresing paerns from big daa is one of he main goals of his book. In his chaper,
More informationSolutions for homework 12
y Soluions for homework Secion Nonlinear sysems: The linearizaion of a nonlinear sysem Consider he sysem y y y y y (i) Skech he nullclines Use a disincive marking for each nullcline so hey can be disinguished
More informationPhysics 221 Fall 2008 Homework #2 Solutions Ch. 2 Due Tues, Sept 9, 2008
Physics 221 Fall 28 Homework #2 Soluions Ch. 2 Due Tues, Sep 9, 28 2.1 A paricle moving along he x-axis moves direcly from posiion x =. m a ime =. s o posiion x = 1. m by ime = 1. s, and hen moves direcly
More informationRC, RL and RLC circuits
Name Dae Time o Complee h m Parner Course/ Secion / Grade RC, RL and RLC circuis Inroducion In his experimen we will invesigae he behavior of circuis conaining combinaions of resisors, capaciors, and inducors.
More informationGravity and SHM Review Questions
Graviy an SHM Review Quesions 1. The mass of Plane X is one-enh ha of he Earh, an is iameer is one-half ha of he Earh. The acceleraion ue o raviy a he surface of Plane X is mos nearly m/s (B) 4m/s (C)
More informationMath Final Exam Solutions
Mah 246 - Final Exam Soluions Friday, July h, 204 () Find explici soluions and give he inerval of definiion o he following iniial value problems (a) ( + 2 )y + 2y = e, y(0) = 0 Soluion: In normal form,
More information10.6 Parametric Equations
0_006.qd /8/05 9:05 AM Page 77 Secion 0.6 77 Parameric Equaions 0.6 Parameric Equaions Wha ou should learn Evaluae ses of parameric equaions for given values of he parameer. Skech curves ha are represened
More informationToday in Physics 218: radiation reaction
Today in Physics 18: radiaion reacion Radiaion reacion The Abraham-Lorenz formula; radiaion reacion force The pah of he elecron in oday s firs example (radial decay grealy exaggeraed) 6 March 004 Physics
More informationPhysics 1402: Lecture 22 Today s Agenda
Physics 142: ecure 22 Today s Agenda Announcemens: R - RV - R circuis Homework 6: due nex Wednesday Inducion / A curren Inducion Self-Inducance, R ircuis X X X X X X X X X long solenoid Energy and energy
More informationMost Probable Phase Portraits of Stochastic Differential Equations and Its Numerical Simulation
Mos Probable Phase Porrais of Sochasic Differenial Equaions and Is Numerical Simulaion Bing Yang, Zhu Zeng and Ling Wang 3 School of Mahemaics and Saisics, Huazhong Universiy of Science and Technology,
More informationLab 10: RC, RL, and RLC Circuits
Lab 10: RC, RL, and RLC Circuis In his experimen, we will invesigae he behavior of circuis conaining combinaions of resisors, capaciors, and inducors. We will sudy he way volages and currens change in
More informationOscillations. Periodic Motion. Sinusoidal Motion. PHY oscillations - J. Hedberg
Oscillaions PHY 207 - oscillaions - J. Hedberg - 2017 1. Periodic Moion 2. Sinusoidal Moion 3. How do we ge his kind of moion? 4. Posiion - Velociy - cceleraion 5. spring wih vecors 6. he reference circle
More informationChapter 7: Solving Trig Equations
Haberman MTH Secion I: The Trigonomeric Funcions Chaper 7: Solving Trig Equaions Le s sar by solving a couple of equaions ha involve he sine funcion EXAMPLE a: Solve he equaion sin( ) The inverse funcions
More informationVariable acceleration, Mixed Exercise 11
Variable acceleraion, Mixed Exercise 11 1 a v 1 P is a res when v 0. 0 1 b s 0 0 v d (1 ) 1 0 1 0 7. The disance ravelled by P is 7. m. 1 a v 6+ a d v 6 + When, a 6+ 0 The acceleraion of P when is 0 m
More informationUCLA Math 135, Winter 2015 Ordinary Differential Equations Review of First-Order Equations from Math 33b
UCLA Mah 35, Winer 05 Ordinary Differenial Equaions Review of Firs-Order Equaions from Mah 33b 7. Numerical Mehods C. David Levermore Deparmen of Mahemaics Universiy of Maryland January 9, 05 7. Numerical
More informationASTR415: Problem Set #5
ASTR45: Problem Se #5 Curran D. Muhlberger Universi of Marland (Daed: April 25, 27) Three ssems of coupled differenial equaions were sudied using inegraors based on Euler s mehod, a fourh-order Runge-Kua
More informationToroidal Field Ripple Induced Excursion of Banana Orbit
Toroidal Field Ripple Induced Excursion of Banana Orbi in Tokamak Plasmas Gao Qingdi Key words TF ripple Banana orbi Local magneic well Thermonuclear plasma ions wihin a okamak mus be achieved wih a finie
More information72 Calculus and Structures
72 Calculus and Srucures CHAPTER 5 DISTANCE AND ACCUMULATED CHANGE Calculus and Srucures 73 Copyrigh Chaper 5 DISTANCE AND ACCUMULATED CHANGE 5. DISTANCE a. Consan velociy Le s ake anoher look a Mary s
More informationLab #2: Kinematics in 1-Dimension
Reading Assignmen: Chaper 2, Secions 2-1 hrough 2-8 Lab #2: Kinemaics in 1-Dimension Inroducion: The sudy of moion is broken ino wo main areas of sudy kinemaics and dynamics. Kinemaics is he descripion
More informationMulti-scale 2D acoustic full waveform inversion with high frequency impulsive source
Muli-scale D acousic full waveform inversion wih high frequency impulsive source Vladimir N Zubov*, Universiy of Calgary, Calgary AB vzubov@ucalgaryca and Michael P Lamoureux, Universiy of Calgary, Calgary
More informationChapter 12: Velocity, acceleration, and forces
To Feel a Force Chaper Spring, Chaper : A. Saes of moion For moion on or near he surface of he earh, i is naural o measure moion wih respec o objecs fixed o he earh. The 4 hr. roaion of he earh has a measurable
More informationOnline Appendix to Solution Methods for Models with Rare Disasters
Online Appendix o Soluion Mehods for Models wih Rare Disasers Jesús Fernández-Villaverde and Oren Levinal In his Online Appendix, we presen he Euler condiions of he model, we develop he pricing Calvo block,
More information2.3 SCHRÖDINGER AND HEISENBERG REPRESENTATIONS
Andrei Tokmakoff, MIT Deparmen of Chemisry, 2/22/2007 2-17 2.3 SCHRÖDINGER AND HEISENBERG REPRESENTATIONS The mahemaical formulaion of he dynamics of a quanum sysem is no unique. So far we have described
More informationBiol. 356 Lab 8. Mortality, Recruitment, and Migration Rates
Biol. 356 Lab 8. Moraliy, Recruimen, and Migraion Raes (modified from Cox, 00, General Ecology Lab Manual, McGraw Hill) Las week we esimaed populaion size hrough several mehods. One assumpion of all hese
More informationSolution: b All the terms must have the dimension of acceleration. We see that, indeed, each term has the units of acceleration
PHYS 54 Tes Pracice Soluions Spring 8 Q: [4] Knowing ha in he ne epression a is acceleraion, v is speed, is posiion and is ime, from a dimensional v poin of view, he equaion a is a) incorrec b) correc
More information2002 November 14 Exam III Physics 191
November 4 Exam III Physics 9 Physical onsans: Earh s free-fall acceleraion = g = 9.8 m/s ircle he leer of he single bes answer. quesion is worh poin Each 3. Four differen objecs wih masses: m = kg, m
More informationPhysics for Scientists and Engineers I
Physics for Scieniss and Engineers I PHY 48, Secion 4 Dr. Beariz Roldán Cuenya Universiy of Cenral Florida, Physics Deparmen, Orlando, FL Chaper - Inroducion I. General II. Inernaional Sysem of Unis III.
More informationSection 3.5 Nonhomogeneous Equations; Method of Undetermined Coefficients
Secion 3.5 Nonhomogeneous Equaions; Mehod of Undeermined Coefficiens Key Terms/Ideas: Linear Differenial operaor Nonlinear operaor Second order homogeneous DE Second order nonhomogeneous DE Soluion o homogeneous
More informationDamped mechanical oscillator: Experiment and detailed energy analysis
1 Damped mechanical oscillaor: Experimen and deailed energy analysis Tommaso Corridoni, DFA, Locarno, Swizerland Michele D Anna, Liceo canonale, Locarno, Swizerland Hans Fuchs, Zurich Universiy of Applied
More informationand v y . The changes occur, respectively, because of the acceleration components a x and a y
Week 3 Reciaion: Chaper3 : Problems: 1, 16, 9, 37, 41, 71. 1. A spacecraf is raveling wih a veloci of v0 = 5480 m/s along he + direcion. Two engines are urned on for a ime of 84 s. One engine gives he
More informationOrdinary Differential Equations
Ordinary Differenial Equaions 5. Examples of linear differenial equaions and heir applicaions We consider some examples of sysems of linear differenial equaions wih consan coefficiens y = a y +... + a
More informationChapter 11 VIBRATORY MOTION
Ch. 11--Vibraory Moion Chaper 11 VIBRATORY MOTION Noe: There are wo areas of ineres when discussing oscillaory moion: he mahemaical characerizaion of vibraing srucures ha generae waves and he ineracion
More informationSimulation-Solving Dynamic Models ABE 5646 Week 2, Spring 2010
Simulaion-Solving Dynamic Models ABE 5646 Week 2, Spring 2010 Week Descripion Reading Maerial 2 Compuer Simulaion of Dynamic Models Finie Difference, coninuous saes, discree ime Simple Mehods Euler Trapezoid
More informationMOMENTUM CONSERVATION LAW
1 AAST/AEDT AP PHYSICS B: Impulse and Momenum Le us run an experimen: The ball is moving wih a velociy of V o and a force of F is applied on i for he ime inerval of. As he resul he ball s velociy changes
More informationLet us start with a two dimensional case. We consider a vector ( x,
Roaion marices We consider now roaion marices in wo and hree dimensions. We sar wih wo dimensions since wo dimensions are easier han hree o undersand, and one dimension is a lile oo simple. However, our
More information(1) (2) Differentiation of (1) and then substitution of (3) leads to. Therefore, we will simply consider the second-order linear system given by (4)
Phase Plane Analysis of Linear Sysems Adaped from Applied Nonlinear Conrol by Sloine and Li The general form of a linear second-order sysem is a c b d From and b bc d a Differeniaion of and hen subsiuion
More informationSTUDY ON A METHOD TO IDENTIFY EXTERNAL FORCES AND DAMPING PARAMETERS OF STRUCTURE FROM RANDOM RESPONSE VALUES
STUDY ON A METHOD TO IDENTIFY ETERNAL FORCES AND DAMPING PARAMETERS OF STRUCTURE FROM RANDOM RESPONSE VALUES 114 M KANDA 1, E MARUTA, D INABA 3 And T MATSUYAMA 4 SUMMARY Proposed in his paper is a mehod
More informationA VARIATIONAL APPROACH TO MULTIRATE INTEGRATION FOR CONSTRAINED SYSTEMS
MULTIBODY DYNAMICS, ECCOMAS Themaic Conference J.C. Samin, P. Fisee eds. Brussels, Belgium, 4-7 July A VARIATIONAL APPROACH TO MULTIRATE INTEGRATION FOR CONSTRAINED SYSTEMS Sigrid Leyendecer, Sina Ober-Blöbaum
More information4.6 One Dimensional Kinematics and Integration
4.6 One Dimensional Kinemaics and Inegraion When he acceleraion a( of an objec is a non-consan funcion of ime, we would like o deermine he ime dependence of he posiion funcion x( and he x -componen of
More information10.1 EXERCISES. y 2 t 2. y 1 t y t 3. y e
66 CHAPTER PARAMETRIC EQUATINS AND PLAR CRDINATES SLUTIN We use a graphing device o produce he graphs for he cases a,,.5,.,,.5,, and shown in Figure 7. Noice ha all of hese curves (ecep he case a ) have
More informationChapter 7 Response of First-order RL and RC Circuits
Chaper 7 Response of Firs-order RL and RC Circuis 7.- The Naural Response of RL and RC Circuis 7.3 The Sep Response of RL and RC Circuis 7.4 A General Soluion for Sep and Naural Responses 7.5 Sequenial
More information