Donald Fussell. October 28, Computer Science Department The University of Texas at Austin. Point Masses and Force Fields.
|
|
- Randall Lamb
- 5 years ago
- Views:
Transcription
1 s Vector Moving s and Coputer Science Departent The University of Texas at Austin October 28, 2014
2 s Vector Moving s Siple classical dynaics - point asses oved by forces Point asses can odel particles They can be grouped to odel acroscopic objects Forces can ove particles Forces can constrain point asses to be part of larger objects Used in particle systes, cloth, rigid bodies, etc.
3 s At tie t, a particle has s Vector Moving s A position A velocity x(t) = v(t) = ẋ(t) = An acceleration [ ] x(t) y(t) a(t) = ẍ(t) = d v(t) dt [ dx(t) ] dt dy(t) dt = d2 x(t) dt 2 = d 2 x(t) dt 2 d 2 y(t) dt 2 y x And a ass
4 Vector s Vector Moving s Ignore acceleration and ass for now and let the velocity be dictated by a function g of position and tie, so ẋ(t) = g( x(t), t). At any tie t, the function g defines a vector field over x. This vector field specifies the velocity of a particle located at position x(t) at tie t. y x
5 s Vector Moving s The equation ẋ(t) = g( x(t), t) is a first-order ordinary differential equation (ODE). We need to use it to copute the path that a particle takes through the vector field if it starts at a given initial position. This path is called an integral curve. The proble is called an initial value proble, a special case of a boundary value proble. y x
6 Moving s s Vector Moving s Given ẋ(t) = g( x(t), t), we want to solve for x(t), the particle s trajectory. We can do this nuerically. Euler s Method Choose a discrete tiestep t. Iterate, taking linear steps in the direction of the velocity at the beginning of each step. x(t + t) = x(t) + ẋ(t) t = x(t) + g( x(t), t) t In sipler tiestep iteration notation x i+1 = x i + v i t This is the siplest (but not the best) way to nuerically integrate a first-order ODE.
7 s Vector Moving s Proble: discrete steps in instantaneously valid directions Bigger steps, bigger errors Tiestep error O( t 2 ), global error O( t) Need to take pretty sall steps, so aybe not efficient Step size needed depends on rate of velocity change ethods allow bigger tiesteps, we ll see soe later y x
8 s Vector Moving s For particle dynaics, we won t have a g( x(t), t) to specify velocities directly, rather we ll copute velocities fro the forces acting on particles Since particles obey Newton s Law: f = a = ẍ And since in general force (and thus acceleration) can depend on the position and velocity of a particle We have a(t) = ẍ(t) = f( x(t), ẋ(t), t)
9 s Vector Moving s Equations with a second derivative in the (like ẍ) are called second-order ODEs. Euler s Method and related techniques handle first-order ODEs. Turning higher-order ODEs into first-order ODEs We have ẍ(t) = f( x(t), ẋ(t), t) If we introduce the velocity variable v(t), we can rewrite this as a pair of coupled first-order ODEs [ ] ẋ(t) = v(t) v(t) = f( x(t), v(t),t) We can treat this as a first-order ODE in phase space [ẋ(t) ] [ ] = v(t) v(t) f( x(t), v(t),t)
10 s Vector Moving s We can concatenate vectors to for vectors in higher-diensional space, called phase space. If x and v are each 3-diensional, we get a 6-diensional phase space. [ ] x(t) v(t)
11 s Vector Moving s We can concatenate vectors to for vectors in higher-diensional space, called phase space. If x and v are each 3-diensional, we get a 6-diensional phase space. We can take tie derivatives in phase space. [ ] x(t) v(t) [ẋ(t) ] v(t)
12 s Vector Moving s We can concatenate vectors to for vectors in higher-diensional space, called phase space. If x and v are each 3-diensional, we get a 6-diensional phase space. We can take tie derivatives in phase space. We can have vector fields in phase space, like our force field. (Yes, I know it sounds like Star Trek.) [ ] x(t) v(t) [ẋ(t) ] v(t) f( x(t), v(t), t)
13 s Vector Moving s We can concatenate vectors to for vectors in higher-diensional space, called phase space. If x and v are each 3-diensional, we get a 6-diensional phase space. We can take tie derivatives in phase space. We can have vector fields in phase space, like our force field. (Yes, I know it sounds like Star Trek.) We can have ODEs in phase space too. [ ] x(t) v(t) [ẋ(t) ] v(t) f( x(t), v(t), t) [ẋ(t) ] [ = v(t) v(t) f( x(t), v(t),t) ]
14 s Vector Moving s [ẋ(t) ] [ We can solve = v(t) v(t) f( x(t), v(t),t) Euler s Method in ] using Euler s Method. Choose a discrete tiestep t and iterate, as before. With tiestep indices x(t + t) = x(t) + v(t) t v(t + t) = v(t) + f( x(t), v(t), t) x i+1 = x i + v i t v i+1 = v i + f i( x i, v i) t t This is just what we did before except we re updating v iteratively along with x.
15 s Vector Moving s Euler s Method is first-order, so The error per tiestep (local error) is proportional to t 2. The cuulative (global) error is proportional to t. This isn t very good, sall tiesteps are often needed to keep the coputation fro diverging. This is a particular proble for stiff systes. For our purposes, these are roughly systes with really stiff springs in the, or soething exerting very strong forces that can change rapidly. That s not a definition, which is uch harder. A related proble is that Euler s Method is not syplectic. It does not conserve energy (or coe close to doing so) as the siulation proceeds.
16 s Vector Moving s We said Euler s Method is local order O( t 2 ) and global order O( t). What do we ean? We siulate fro t 0 to T by taking n tiesteps of size t. The total error at tie T, e( t, T ) = x(t ) x n, the difference between the correct position at T and the siulated position after the nth tiestep. e( t, T ) is a function of both the tiestep size t and the total siulation tie T. The order of an integration ethod describes how e( t, T ) changes as a function of t, holding T constant.
17 s Vector Moving s Local order O( t 2 ) eans changing t changes the worst-case truncation error proportionally to t 2. Global order O( t) eans a change in t changes the worst case cuulative truncation error of the entire siulation proportionally to t. This doesn t tell you uch about how the siulation error grows as the overall siulation tie T grows. This depends on any factors, what ODEs you re integrating, the initial conditions, etc. A ethod with a higher global order is better than a lower-order ethod because it gets worse ore slowly as you increase the tiestep size, not because you have any idea how far off the right answer it will get in the actual siulation.
18 s Vector Moving s Assue that f is only a function of position and not of velocity, i.e. it is f( x(t), t). We don t odel wind resistance, for exaple, but we can still handle a lot of interesting forces. Now we can easily fix Euler s Method to ake it syplectic. Sei-iplicit Euler Integration We can use v i+1 instead of v i to copute x i+1 v i+1 = v i + f i( x i) t x i+1 = x i + v i+1 t Alternatively, we can use x i+1 instead of x i to copute v i+1 x i+1 = x i + v i t v i+1 = v i + f i( x i+1) t
19 s Vector Moving s Euler s Method can be obtained fro the Taylor series for x(t + t) and v(t + t). x(t + t) = x(t) + ẋ(t) t + O( t 2 ) v(t + t) = v(t) + v(t) t + O( t 2 ) We just throw away the quadratic and higher order ters, which is why the local error is proportional to t 2. We can use siilar derivations to get better ethods. For instance, the techniques we ve seen so far aren t tie-reversible, which proper physics should be. Let s try for that property.
20 s Vector Moving s We can do Taylor series in both tie directions x(t + t) = x(t) + ẋ(t) t + 1/2ẍ(t) t 2 + 1/6... x(t) t 3 + O( t 4 ) x(t t) = x(t) ẋ(t) t + 1/2ẍ(t) t 2 1/6... x(t) t 3 + O( t 4 ) Integration Add the and apply soe algebra to get x(t + t) = 2 x(t) x(t t) + ẍ(t) t 2 + O( t 4 ) Throw away the O( t 4 ) local error ter, write as an iteration, put in the forcing function for the acceleration, and we have x i+1 = 2 x i x i 1 + f i( x i) t2 We re given initial x 0 and v 0, but we also need x 1, so x 1 = x 0 + v 0 t fi( x 0) t2
21 s Vector Moving s advantages It is tie-reversible, syplectic, local O(t 4 ) and global O(t 2 ), or second-order. Easy to use for constraints between particles (e.g. infinitely stiff springs), especially for rigid bodies, cloth, etc. See Jakobsen s faous GaaSutra article guide/ /jacobson pfv.ht disadvantages Multistep ethod, needs adjustent for non-constant t. x i+1 = x i + ( x i x i 1 ) t i + f i ( x i ) t i 1 t i + t i 1 t i 2 Doesn t calculate velocities.
22 s Vector Moving s If we need velocities, we can just calculate the fro the point positions. v i = xi+1 xi 1 + O( t 2 xi+1 xi ) or v i+1 = + O( t) t i + t i 1 t i More coonly, they are incorporated into the iteration as the Method Expand x(t + t) = x(t) + v(t) t a t2 + O( t 3 ) Expand forward v(t + t) = v(t) + a(t) t + O( t 2 ) and reverse v(t) = v(t + t) a(t + t) t + O( t 2 ) Add the latter two to get v(t + t) = v(t) + 1 ( a(t) + a(t + t)) t 2 So we have x i+1 = x i + v i t fi( x i) t2 v i+1 = v i + f i( x i) + f i+1( x i+1) t 2
Physically Based Modeling CS Notes Spring 1997 Particle Collision and Contact
Physically Based Modeling CS 15-863 Notes Spring 1997 Particle Collision and Contact 1 Collisions with Springs Suppose we wanted to ipleent a particle siulator with a floor : a solid horizontal plane which
More informationPh 20.3 Numerical Solution of Ordinary Differential Equations
Ph 20.3 Nuerical Solution of Ordinary Differential Equations Due: Week 5 -v20170314- This Assignent So far, your assignents have tried to failiarize you with the hardware and software in the Physics Coputing
More informationKinematics and dynamics, a computational approach
Kineatics and dynaics, a coputational approach We begin the discussion of nuerical approaches to echanics with the definition for the velocity r r ( t t) r ( t) v( t) li li or r( t t) r( t) v( t) t for
More informationma x = -bv x + F rod.
Notes on Dynaical Systes Dynaics is the study of change. The priary ingredients of a dynaical syste are its state and its rule of change (also soeties called the dynaic). Dynaical systes can be continuous
More informationPHY307F/407F - Computational Physics Background Material for Expt. 3 - Heat Equation David Harrison
INTRODUCTION PHY37F/47F - Coputational Physics Background Material for Expt 3 - Heat Equation David Harrison In the Pendulu Experient, we studied the Runge-Kutta algorith for solving ordinary differential
More informationHORIZONTAL MOTION WITH RESISTANCE
DOING PHYSICS WITH MATLAB MECHANICS HORIZONTAL MOTION WITH RESISTANCE Ian Cooper School of Physics, Uniersity of Sydney ian.cooper@sydney.edu.au DOWNLOAD DIRECTORY FOR MATLAB SCRIPTS ec_fr_b. This script
More informationNote-A-Rific: Mechanical
Note-A-Rific: Mechanical Kinetic You ve probably heard of inetic energy in previous courses using the following definition and forula Any object that is oving has inetic energy. E ½ v 2 E inetic energy
More informationSome Perspective. Forces and Newton s Laws
Soe Perspective The language of Kineatics provides us with an efficient ethod for describing the otion of aterial objects, and we ll continue to ake refineents to it as we introduce additional types of
More informationlecture 36: Linear Multistep Mehods: Zero Stability
95 lecture 36: Linear Multistep Mehods: Zero Stability 5.6 Linear ultistep ethods: zero stability Does consistency iply convergence for linear ultistep ethods? This is always the case for one-step ethods,
More informationNB1140: Physics 1A - Classical mechanics and Thermodynamics Problem set 2 - Forces and energy Week 2: November 2016
NB1140: Physics 1A - Classical echanics and Therodynaics Proble set 2 - Forces and energy Week 2: 21-25 Noveber 2016 Proble 1. Why force is transitted uniforly through a assless string, a assless spring,
More informationLesson 24: Newton's Second Law (Motion)
Lesson 24: Newton's Second Law (Motion) To really appreciate Newton s Laws, it soeties helps to see how they build on each other. The First Law describes what will happen if there is no net force. The
More informationPhysics 201, Lecture 15
Physics 0, Lecture 5 Today s Topics q More on Linear Moentu And Collisions Elastic and Perfect Inelastic Collision (D) Two Diensional Elastic Collisions Exercise: Billiards Board Explosion q Multi-Particle
More informationForce and dynamics with a spring, analytic approach
Force and dynaics with a spring, analytic approach It ay strie you as strange that the first force we will discuss will be that of a spring. It is not one of the four Universal forces and we don t use
More information26 Impulse and Momentum
6 Ipulse and Moentu First, a Few More Words on Work and Energy, for Coparison Purposes Iagine a gigantic air hockey table with a whole bunch of pucks of various asses, none of which experiences any friction
More informationChapter 1: Basics of Vibrations for Simple Mechanical Systems
Chapter 1: Basics of Vibrations for Siple Mechanical Systes Introduction: The fundaentals of Sound and Vibrations are part of the broader field of echanics, with strong connections to classical echanics,
More information12 Towards hydrodynamic equations J Nonlinear Dynamics II: Continuum Systems Lecture 12 Spring 2015
18.354J Nonlinear Dynaics II: Continuu Systes Lecture 12 Spring 2015 12 Towards hydrodynaic equations The previous classes focussed on the continuu description of static (tie-independent) elastic systes.
More informationPY /005 Practice Test 1, 2004 Feb. 10
PY 205-004/005 Practice Test 1, 2004 Feb. 10 Print nae Lab section I have neither given nor received unauthorized aid on this test. Sign ature: When you turn in the test (including forula page) you ust
More informationFrame with 6 DOFs. v m. determining stiffness, k k = F / water tower deflected water tower dynamic response model
CE 533, Fall 2014 Undaped SDOF Oscillator 1 / 6 What is a Single Degree of Freedo Oscillator? The siplest representation of the dynaic response of a civil engineering structure is the single degree of
More informationToday s s topics are: Collisions and Momentum Conservation. Momentum Conservation
Today s s topics are: Collisions and P (&E) Conservation Ipulsive Force Energy Conservation How can we treat such an ipulsive force? Energy Conservation Ipulsive Force and Ipulse [Exaple] an ipulsive force
More informationI. Understand get a conceptual grasp of the problem
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Departent o Physics Physics 81T Fall Ter 4 Class Proble 1: Solution Proble 1 A car is driving at a constant but unknown velocity,, on a straightaway A otorcycle is
More informationlecture 37: Linear Multistep Methods: Absolute Stability, Part I lecture 38: Linear Multistep Methods: Absolute Stability, Part II
lecture 37: Linear Multistep Methods: Absolute Stability, Part I lecture 3: Linear Multistep Methods: Absolute Stability, Part II 5.7 Linear ultistep ethods: absolute stability At this point, it ay well
More informationFinite fields. and we ve used it in various examples and homework problems. In these notes I will introduce more finite fields
Finite fields I talked in class about the field with two eleents F 2 = {, } and we ve used it in various eaples and hoework probles. In these notes I will introduce ore finite fields F p = {,,...,p } for
More informationSupport Vector Machines MIT Course Notes Cynthia Rudin
Support Vector Machines MIT 5.097 Course Notes Cynthia Rudin Credit: Ng, Hastie, Tibshirani, Friedan Thanks: Şeyda Ertekin Let s start with soe intuition about argins. The argin of an exaple x i = distance
More informationLecture 21. Interior Point Methods Setup and Algorithm
Lecture 21 Interior Point Methods In 1984, Kararkar introduced a new weakly polynoial tie algorith for solving LPs [Kar84a], [Kar84b]. His algorith was theoretically faster than the ellipsoid ethod and
More informationOne Dimensional Collisions
One Diensional Collisions These notes will discuss a few different cases of collisions in one diension, arying the relatie ass of the objects and considering particular cases of who s oing. Along the way,
More informationCHAPTER 15: Vibratory Motion
CHAPTER 15: Vibratory Motion courtesy of Richard White courtesy of Richard White 2.) 1.) Two glaring observations can be ade fro the graphic on the previous slide: 1.) The PROJECTION of a point on a circle
More informationOSCILLATIONS AND WAVES
OSCILLATIONS AND WAVES OSCILLATION IS AN EXAMPLE OF PERIODIC MOTION No stories this tie, we are going to get straight to the topic. We say that an event is Periodic in nature when it repeats itself in
More informationProjectile Motion with Air Resistance (Numerical Modeling, Euler s Method)
Projectile Motion with Air Resistance (Nuerical Modeling, Euler s Method) Theory Euler s ethod is a siple way to approxiate the solution of ordinary differential equations (ode s) nuerically. Specifically,
More informationCHAPTER 1 MOTION & MOMENTUM
CHAPTER 1 MOTION & MOMENTUM SECTION 1 WHAT IS MOTION? All atter is constantly in MOTION Motion involves a CHANGE in position. An object changes position relative to a REFERENCE POINT. DISTANCE is the total
More informationSIMPLE HARMONIC MOTION: NEWTON S LAW
SIMPLE HARMONIC MOTION: NEWTON S LAW siple not siple PRIOR READING: Main 1.1, 2.1 Taylor 5.1, 5.2 http://www.yoops.org/twocw/it/nr/rdonlyres/physics/8-012fall-2005/7cce46ac-405d-4652-a724-64f831e70388/0/chp_physi_pndul.jpg
More information2.003 Engineering Dynamics Problem Set 2 Solutions
.003 Engineering Dynaics Proble Set Solutions This proble set is priarily eant to give the student practice in describing otion. This is the subject of kineatics. It is strongly recoended that you study
More informationComparison of Stability of Selected Numerical Methods for Solving Stiff Semi- Linear Differential Equations
International Journal of Applied Science and Technology Vol. 7, No. 3, Septeber 217 Coparison of Stability of Selected Nuerical Methods for Solving Stiff Sei- Linear Differential Equations Kwaku Darkwah
More informationVector Spaces in Physics 8/6/2015. Chapter 4. Practical Examples.
Vector Spaces in Physics 8/6/15 Chapter 4. Practical Exaples. In this chapter we will discuss solutions to two physics probles where we ae use of techniques discussed in this boo. In both cases there are
More informationXI PHYSICS M. AFFAN KHAN LECTURER PHYSICS, AKHSS, K. https://promotephysics.wordpress.com
XI PHYSICS M. AFFAN KHAN LECTURER PHYSICS, AKHSS, K affan_414@live.co https://prootephysics.wordpress.co [MOTION] CHAPTER NO. 3 In this chapter we are going to discuss otion in one diension in which we
More informationModel Fitting. CURM Background Material, Fall 2014 Dr. Doreen De Leon
Model Fitting CURM Background Material, Fall 014 Dr. Doreen De Leon 1 Introduction Given a set of data points, we often want to fit a selected odel or type to the data (e.g., we suspect an exponential
More informationDimensions and Units
Civil Engineering Hydraulics Mechanics of Fluids and Modeling Diensions and Units You already know how iportant using the correct diensions can be in the analysis of a proble in fluid echanics If you don
More informationAnalyzing Simulation Results
Analyzing Siulation Results Dr. John Mellor-Cruey Departent of Coputer Science Rice University johnc@cs.rice.edu COMP 528 Lecture 20 31 March 2005 Topics for Today Model verification Model validation Transient
More information2 Q 10. Likewise, in case of multiple particles, the corresponding density in 2 must be averaged over all
Lecture 6 Introduction to kinetic theory of plasa waves Introduction to kinetic theory So far we have been odeling plasa dynaics using fluid equations. The assuption has been that the pressure can be either
More informationChapter VI: Motion in the 2-D Plane
Chapter VI: Motion in the -D Plane Now that we have developed and refined our vector calculus concepts, we can ove on to specific application of otion in the plane. In this regard, we will deal with: projectile
More informationNUMERICAL MODELLING OF THE TYRE/ROAD CONTACT
NUMERICAL MODELLING OF THE TYRE/ROAD CONTACT PACS REFERENCE: 43.5.LJ Krister Larsson Departent of Applied Acoustics Chalers University of Technology SE-412 96 Sweden Tel: +46 ()31 772 22 Fax: +46 ()31
More informationUSEFUL HINTS FOR SOLVING PHYSICS OLYMPIAD PROBLEMS. By: Ian Blokland, Augustana Campus, University of Alberta
1 USEFUL HINTS FOR SOLVING PHYSICS OLYMPIAD PROBLEMS By: Ian Bloland, Augustana Capus, University of Alberta For: Physics Olypiad Weeend, April 6, 008, UofA Introduction: Physicists often attept to solve
More informationFigure 1: Equivalent electric (RC) circuit of a neurons membrane
Exercise: Leaky integrate and fire odel of neural spike generation This exercise investigates a siplified odel of how neurons spike in response to current inputs, one of the ost fundaental properties of
More informationAn Approximate Model for the Theoretical Prediction of the Velocity Increase in the Intermediate Ballistics Period
An Approxiate Model for the Theoretical Prediction of the Velocity... 77 Central European Journal of Energetic Materials, 205, 2(), 77-88 ISSN 2353-843 An Approxiate Model for the Theoretical Prediction
More informationNational 5 Summary Notes
North Berwick High School Departent of Physics National 5 Suary Notes Unit 3 Energy National 5 Physics: Electricity and Energy 1 Throughout the Course, appropriate attention should be given to units, prefixes
More informationHonors Lab 4.5 Freefall, Apparent Weight, and Friction
Nae School Date Honors Lab 4.5 Freefall, Apparent Weight, and Friction Purpose To investigate the vector nature of forces To practice the use free-body diagras (FBDs) To learn to apply Newton s Second
More informationNow multiply the left-hand-side by ω and the right-hand side by dδ/dt (recall ω= dδ/dt) to get:
Equal Area Criterion.0 Developent of equal area criterion As in previous notes, all powers are in per-unit. I want to show you the equal area criterion a little differently than the book does it. Let s
More informationEnergy and Momentum: The Ballistic Pendulum
Physics Departent Handout -10 Energy and Moentu: The Ballistic Pendulu The ballistic pendulu, first described in the id-eighteenth century, applies principles of echanics to the proble of easuring the
More informationSolving initial value problems by residual power series method
Theoretical Matheatics & Applications, vol.3, no.1, 13, 199-1 ISSN: 179-9687 (print), 179-979 (online) Scienpress Ltd, 13 Solving initial value probles by residual power series ethod Mohaed H. Al-Sadi
More informationlecture 35: Linear Multistep Mehods: Truncation Error
88 lecture 5: Linear Multistep Meods: Truncation Error 5.5 Linear ultistep etods One-step etods construct an approxiate solution x k+ x(t k+ ) using only one previous approxiation, x k. Tis approac enoys
More informationName Period. What force did your partner s exert on yours? Write your answer in the blank below:
Nae Period Lesson 7: Newton s Third Law and Passive Forces 7.1 Experient: Newton s 3 rd Law Forces of Interaction (a) Tea up with a partner to hook two spring scales together to perfor the next experient:
More information(a) Why cannot the Carnot cycle be applied in the real world? Because it would have to run infinitely slowly, which is not useful.
PHSX 446 FINAL EXAM Spring 25 First, soe basic knowledge questions You need not show work here; just give the answer More than one answer ight apply Don t waste tie transcribing answers; just write on
More informationDefinition of Work, The basics
Physics 07 Lecture 16 Lecture 16 Chapter 11 (Work) v Eploy conservative and non-conservative forces v Relate force to potential energy v Use the concept of power (i.e., energy per tie) Chapter 1 v Define
More informationPhysics 221A: HW3 solutions
Physics 22A: HW3 solutions October 22, 202. a) It will help to start things off by doing soe gaussian integrals. Let x be a real vector of length, and let s copute dxe 2 xt Ax, where A is soe real atrix.
More informationEE5900 Spring Lecture 4 IC interconnect modeling methods Zhuo Feng
EE59 Spring Parallel LSI AD Algoriths Lecture I interconnect odeling ethods Zhuo Feng. Z. Feng MTU EE59 So far we ve considered only tie doain analyses We ll soon see that it is soeties preferable to odel
More informationEN40: Dynamics and Vibrations. Midterm Examination Tuesday March
EN4: Dynaics and Vibrations Midter Exaination Tuesday March 8 16 School of Engineering Brown University NAME: General Instructions No collaboration of any kind is peritted on this exaination. You ay bring
More information1 Brownian motion and the Langevin equation
Figure 1: The robust appearance of Robert Brown (1773 1858) 1 Brownian otion and the Langevin equation In 1827, while exaining pollen grains and the spores of osses suspended in water under a icroscope,
More informationBALLISTIC PENDULUM. EXPERIMENT: Measuring the Projectile Speed Consider a steel ball of mass
BALLISTIC PENDULUM INTRODUCTION: In this experient you will use the principles of conservation of oentu and energy to deterine the speed of a horizontally projected ball and use this speed to predict the
More informationLecture 23: Central Force Motion
Lectue 3: Cental Foce Motion Many of the foces we encounte in natue act between two paticles along the line connecting the Gavity, electicity, and the stong nuclea foce ae exaples These types of foces
More informationOn the role of quadrature rules and system dimensions in variational multirate integrators
The 3 rd Joint International Conference on Multibody Syste Dynaics The 7 th Asian Conference on Multibody Dynaics June 30-July 3, 04, BEXCO, Busan, Korea On the role of quadrature rules and syste diensions
More informationAnalysis of Impulsive Natural Phenomena through Finite Difference Methods A MATLAB Computational Project-Based Learning
Analysis of Ipulsive Natural Phenoena through Finite Difference Methods A MATLAB Coputational Project-Based Learning Nicholas Kuia, Christopher Chariah, Mechatronics Engineering, Vaughn College of Aeronautics
More informationUnderstanding the coefficient of restitution (COR) using mass/spring systems
Understanding the coefficient of restitution (COR) using ass/spring systes Dr. David Kagan Departent of Physics California State University, Chico Chico, CA 9599-00 dkagan@csuchico.edu The coefficient
More informationANALYSIS OF HALL-EFFECT THRUSTERS AND ION ENGINES FOR EARTH-TO-MOON TRANSFER
IEPC 003-0034 ANALYSIS OF HALL-EFFECT THRUSTERS AND ION ENGINES FOR EARTH-TO-MOON TRANSFER A. Bober, M. Guelan Asher Space Research Institute, Technion-Israel Institute of Technology, 3000 Haifa, Israel
More informationEULER EQUATIONS. We start by considering how time derivatives are effected by rotation. Consider a vector defined in the two systems by
EULER EQUATIONS We now consider another approach to rigid body probles based on looking at the change needed in Newton s Laws if an accelerated coordinate syste is used. We start by considering how tie
More informationProblem Set 2. Chapter 1 Numerical:
Chapter 1 Nuerical: roble Set 16. The atoic radius of xenon is 18 p. Is that consistent with its b paraeter of 5.15 1 - L/ol? Hint: what is the volue of a ole of xenon atos and how does that copare to
More informationMomentum. February 15, Table of Contents. Momentum Defined. Momentum Defined. p =mv. SI Unit for Momentum. Momentum is a Vector Quantity.
Table of Contents Click on the topic to go to that section Moentu Ipulse-Moentu Equation The Moentu of a Syste of Objects Conservation of Moentu Types of Collisions Collisions in Two Diensions Moentu Return
More informationU V. r In Uniform Field the Potential Difference is V Ed
SPHI/W nit 7.8 Electric Potential Page of 5 Notes Physics Tool box Electric Potential Energy the electric potential energy stored in a syste k of two charges and is E r k Coulobs Constant is N C 9 9. E
More informationThe path integral approach in the frame work of causal interpretation
Annales de la Fondation Louis de Broglie, Volue 28 no 1, 2003 1 The path integral approach in the frae work of causal interpretation M. Abolhasani 1,2 and M. Golshani 1,2 1 Institute for Studies in Theoretical
More informationWhat is the instantaneous acceleration (2nd derivative of time) of the field? Sol. The Euler-Lagrange equations quickly yield:
PHYSICS 75: The Standard Model Midter Exa Solution Key. [3 points] Short Answer (6 points each (a In words, explain how to deterine the nuber of ediator particles are generated by a particular local gauge
More informationCh 12: Variations on Backpropagation
Ch 2: Variations on Backpropagation The basic backpropagation algorith is too slow for ost practical applications. It ay take days or weeks of coputer tie. We deonstrate why the backpropagation algorith
More informationPage 1. Physics 131: Lecture 22. Today s Agenda. SHM and Circles. Position
Physics 3: ecture Today s genda Siple haronic otion Deinition Period and requency Position, velocity, and acceleration Period o a ass on a spring Vertical spring Energy and siple haronic otion Energy o
More informationChaotic Coupled Map Lattices
Chaotic Coupled Map Lattices Author: Dustin Keys Advisors: Dr. Robert Indik, Dr. Kevin Lin 1 Introduction When a syste of chaotic aps is coupled in a way that allows the to share inforation about each
More information8.1 Force Laws Hooke s Law
8.1 Force Laws There are forces that don't change appreciably fro one instant to another, which we refer to as constant in tie, and forces that don't change appreciably fro one point to another, which
More information27 Oscillations: Introduction, Mass on a Spring
Chapter 7 Oscillations: Introduction, Mass on a Spring 7 Oscillations: Introduction, Mass on a Spring If a siple haronic oscillation proble does not involve the tie, you should probably be using conservation
More informationPhysics 221B: Solution to HW # 6. 1) Born-Oppenheimer for Coupled Harmonic Oscillators
Physics B: Solution to HW # 6 ) Born-Oppenheier for Coupled Haronic Oscillators This proble is eant to convince you of the validity of the Born-Oppenheier BO) Approxiation through a toy odel of coupled
More informationA Self-Organizing Model for Logical Regression Jerry Farlow 1 University of Maine. (1900 words)
1 A Self-Organizing Model for Logical Regression Jerry Farlow 1 University of Maine (1900 words) Contact: Jerry Farlow Dept of Matheatics Univeristy of Maine Orono, ME 04469 Tel (07) 866-3540 Eail: farlow@ath.uaine.edu
More informationPhysicsAndMathsTutor.com
. A raindrop falls vertically under gravity through a cloud. In a odel of the otion the raindrop is assued to be spherical at all ties and the cloud is assued to consist of stationary water particles.
More informationTactics Box 2.1 Interpreting Position-versus-Time Graphs
1D kineatic Retake Assignent Due: 4:32p on Friday, October 31, 2014 You will receive no credit for ites you coplete after the assignent is due. Grading Policy Tactics Box 2.1 Interpreting Position-versus-Tie
More informationTutorial Exercises: Incorporating constraints
Tutorial Exercises: Incorporating constraints 1. A siple pendulu of length l ass is suspended fro a pivot of ass M that is free to slide on a frictionless wire frae in the shape of a parabola y = ax. The
More informationNumerical Studies of a Nonlinear Heat Equation with Square Root Reaction Term
Nuerical Studies of a Nonlinear Heat Equation with Square Root Reaction Ter Ron Bucire, 1 Karl McMurtry, 1 Ronald E. Micens 2 1 Matheatics Departent, Occidental College, Los Angeles, California 90041 2
More informationCHECKLIST. r r. Newton s Second Law. natural frequency ω o (rad.s -1 ) (Eq ) a03/p1/waves/waves doc 9:19 AM 29/03/05 1
PHYS12 Physics 1 FUNDAMENTALS Module 3 OSCILLATIONS & WAVES Text Physics by Hecht Chapter 1 OSCILLATIONS Sections: 1.5 1.6 Exaples: 1.6 1.7 1.8 1.9 CHECKLIST Haronic otion, periodic otion, siple haronic
More informationPhysics Circular Motion: Energy and Momentum Conservation. Science and Mathematics Education Research Group
F FA ACULTY C U L T Y OF O F EDUCATION E D U C A T I O N Departent of Curriculu and Pedagogy Physics Circular Motion: Energy and Moentu Conservation Science and Matheatics Education Research Group Supported
More informationSimple Schemes of Multi anchored Flexible Walls Dynamic Behavior
6 th International Conference on Earthquake Geotechnical Engineering -4 Noveber 05 Christchurch, New Zealand Siple Schees of Multi anchored Flexible Walls Dynaic Behavior A. D. Garini ABSTRACT Siple schees
More informationMassachusetts Institute of Technology Quantum Mechanics I (8.04) Spring 2005 Solutions to Problem Set 4
Massachusetts Institute of Technology Quantu Mechanics I (8.04) Spring 2005 Solutions to Proble Set 4 By Kit Matan 1. X-ray production. (5 points) Calculate the short-wavelength liit for X-rays produced
More information16.333: Lecture # 7. Approximate Longitudinal Dynamics Models
16.333: Lecture # 7 Approxiate Longitudinal Dynaics Models A couple ore stability derivatives Given ode shapes found identify sipler odels that capture the ain responses Fall 24 16.333 6 1 More Stability
More informationFeedforward Networks
Feedforward Networks Gradient Descent Learning and Backpropagation Christian Jacob CPSC 433 Christian Jacob Dept.of Coputer Science,University of Calgary CPSC 433 - Feedforward Networks 2 Adaptive "Prograing"
More informationBoosting with log-loss
Boosting with log-loss Marco Cusuano-Towner Septeber 2, 202 The proble Suppose we have data exaples {x i, y i ) i =... } for a two-class proble with y i {, }. Let F x) be the predictor function with the
More informationPrinciples of Optimal Control Spring 2008
MIT OpenCourseWare http://ocw.it.edu 16.323 Principles of Optial Control Spring 2008 For inforation about citing these aterials or our Ters of Use, visit: http://ocw.it.edu/ters. 16.323 Lecture 10 Singular
More informationFeedforward Networks
Feedforward Neural Networks - Backpropagation Feedforward Networks Gradient Descent Learning and Backpropagation CPSC 533 Fall 2003 Christian Jacob Dept.of Coputer Science,University of Calgary Feedforward
More informationCOS 424: Interacting with Data. Written Exercises
COS 424: Interacting with Data Hoework #4 Spring 2007 Regression Due: Wednesday, April 18 Written Exercises See the course website for iportant inforation about collaboration and late policies, as well
More informationBayesian inference for stochastic differential mixed effects models - initial steps
Bayesian inference for stochastic differential ixed effects odels - initial steps Gavin Whitaker 2nd May 2012 Supervisors: RJB and AG Outline Mixed Effects Stochastic Differential Equations (SDEs) Bayesian
More information1 Generalization bounds based on Rademacher complexity
COS 5: Theoretical Machine Learning Lecturer: Rob Schapire Lecture #0 Scribe: Suqi Liu March 07, 08 Last tie we started proving this very general result about how quickly the epirical average converges
More informationPhysics 2107 Oscillations using Springs Experiment 2
PY07 Oscillations using Springs Experient Physics 07 Oscillations using Springs Experient Prelab Read the following bacground/setup and ensure you are failiar with the concepts and theory required for
More informationFeedforward Networks. Gradient Descent Learning and Backpropagation. Christian Jacob. CPSC 533 Winter 2004
Feedforward Networks Gradient Descent Learning and Backpropagation Christian Jacob CPSC 533 Winter 2004 Christian Jacob Dept.of Coputer Science,University of Calgary 2 05-2-Backprop-print.nb Adaptive "Prograing"
More informationm potential kinetic forms of energy.
Spring, Chapter : A. near the surface of the earth. The forces of gravity and an ideal spring are conservative forces. With only the forces of an ideal spring and gravity acting on a ass, energy F F will
More informationOn Lotka-Volterra Evolution Law
Advanced Studies in Biology, Vol. 3, 0, no. 4, 6 67 On Lota-Volterra Evolution Law Farruh Muhaedov Faculty of Science, International Islaic University Malaysia P.O. Box, 4, 570, Kuantan, Pahang, Malaysia
More informationTHERMAL ENDURANCE OF UNREINFORCED UNSATURATED POLYESTERS AND VINYL ESTER RESINS = (1) ln = COMPOSITES & POLYCON 2009
Aerican Coposites Manufacturers Association January 15-17, 29 Tapa, FL USA Abstract THERMAL ENDURANCE OF UNREINFORCED UNSATURATED POLYESTERS AND VINYL ESTER RESINS by Thore M. Klaveness, Reichhold AS In
More informationPhysics 2210 Fall smartphysics 20 Conservation of Angular Momentum 21 Simple Harmonic Motion 11/23/2015
Physics 2210 Fall 2015 sartphysics 20 Conservation of Angular Moentu 21 Siple Haronic Motion 11/23/2015 Exa 4: sartphysics units 14-20 Midter Exa 2: Day: Fri Dec. 04, 2015 Tie: regular class tie Section
More informationFeature Extraction Techniques
Feature Extraction Techniques Unsupervised Learning II Feature Extraction Unsupervised ethods can also be used to find features which can be useful for categorization. There are unsupervised ethods that
More informationWe last left off by talking about how the area under a force vs. time curve is impulse.
Lecture 11 Ipulse and Moentu We last left off by talking about how the area under a force vs. tie curve is ipulse. Recall that for our golf ball we had a strongly peaked force curve: F F avg t You have
More informationPY241 Solutions Set 9 (Dated: November 7, 2002)
PY241 Solutions Set 9 (Dated: Noveber 7, 2002) 9-9 At what displaceent of an object undergoing siple haronic otion is the agnitude greatest for the... (a) velocity? The velocity is greatest at x = 0, the
More information