A645/A445: Exercise #1. The Kepler Problem
|
|
- Rosamund Curtis
- 6 years ago
- Views:
Transcription
1 A645/A445: Exercise #1 The Kepler Problem Due: 2017 Sep 26 1 Numerical solution to the one-degree-of-freedom implicit equation The one-dimensional implicit solution is: t = t o + x x o 2(E U(x)). (1) The turning points are defined as the coordinates for which ẋ = 0: {x : 2(E U(x)) = 0}. There will be two turning points in general: x,x +. It immediately follows that the orbital period found a bound orbit is: P = 2 x+ x 2(E U(x)). (2) To evaluate the integral in equation (2), make the following trigonometric change variables: where θ [ π/2, π/2] and x = a + + a sinθ (3) a ± x + ± x. 2 Substituting equation (3) into equation (2) gives: π/2 dθcosθ P = 2. (4) π/2 2(E U(x)) 1
2 Near a turning point, x to be explicit, the quantity 2(E U(x)) goes to zero linearly, that is: lim 2(E U(x)) = K(x x ) = Ka (1 + cosθ) = 2Ka cos 2 (θ/2) x x where K 2 du. x From this, it is easy to convince yourself from this that the square root divergence is removed from the integrand using this transformation. This works for square root divergences of all types. Similarly, one can solve for the one-dimensional implicit solution: θ(x) dθcosθ t = t o +. (5) θ(x o ) 2(E U(x)) where x and θ are related by equation (3). (a) Use the transformation in equations (3) and (4) to solve for the period of a harmonic oscillator numerically. To perform the integral, use the trapezoidal rule: b a f (x) h N i=1 where h = (b a)/n and x i = a + ih. f (x i 1 ) + f (x i ) 2 (b) Compare with the exact, analytic solution as a function of N. E.g. graph the numerical error as a function of N on a log log plot. Interpret this. 1. Pick a point 2/3 of the way from the lower to upper turning point for your harmonic oscillator and solve equation (??) using equation (6). Repeat your log log error plot for this problem. 2 The full two-body problem In Problem Set 1, you showed that the two-body problem for reduced to a single equation for the radial degreeof-freedom in polar coordinates r,θ is where L = r 2 θ is the angular momentum. (6) r = du dr + L2 r 3 (7) 2017 Sep 10/MDW 2
3 (a) Write down the implicit solution for t = t(r). (b) Using L = const, write down an analogous implicit solution for theta(t) = θ(r(t)). 3 The analytic solution of the gravitational n-body problem The gravitational potential for the gravitational two-body problem is most often written as: U(r) = GM (8) r where G is the gravitational constant and M is the total mass. We will now show that equation (7) can be solved analytically! (a) Change variables equations (7) and (8) to derive: d 2 u dθ 2 + u = GM L 2. This is a linear inhomogeneous 2nd-order differential equation. Solve this to yield: where e and ϖ are constants. u = 1 r = GM [1 + ecos(θ ϖ)] (9) L2 (b) Replace the angular momentum L by another constant of integration a defined by L 2 = GMa(1 e 2 ) to show that the shape of the orbit is given by: r = a(1 e2 ) 1 + ecos f (10) where f = θ ϖ is known as the true anomaly in celestial mechanics. (c) The closest approach of the two bodies occurs at f = 0 or θ = ϖ and ϖ is known, therefore, as the azimuth of periapsis. Using these relations, show that the periapsis value is r p = a(1 e) and that the total energy is E = GM 2a Sep 10/MDW 3
4 (d) Equation (10) is a parametric equation for an ellipse (Kepler s first law)! Show that that major axis (twice the full extent of the long axis) has the value a and that the minor axis (twice the full extent of the short axis) is related to a and e as e 2 = 1 b 2 /a 2. (e) Since θ = L/r 2 we have t = t2 t 1 dt = 1 L f2 f 1 d f r 2. (11) This equation, then implicitly determines θ = ϖ + f (t) and r = r(t) follows using equation (9). Geometrically, (r r d f /2) is the infinitesimal triangle subtended by the arc d f to the origin and therefore the integral above is twice the area of the ellipse if we integrate from f 1 = 0 to f = 2π. Use this to show that P = 2π(a 3 /GM) 1/2, which is Kepler s third law! There are more fancy things we can do with these equations but this is enough for now. 4 Numerical solution of the gravitational n-body problem Now that we know what to expect, let s solve the differential equations numerically. This is also good practice for problems whose characteristics are unknown. (a) To do this, let us return to the general equation of motion for the reduced two-body problem: d 2 x dt 2 = U(r) = GM x x 3. Choose the units G = M = 1. With no loss of generality, we can consider motion on a single plane (why?). Choose the initial conditions x = (0.140,0.078,0.0) and ẋ = ( 1.56,2.78,0.0), t 1 = 0 and solve Newton s equations Cartesian coordinates with an ODE solver to get the trajectory between t 1 and t 2 = 2. Decrease the step size until the solution converges; that is, it gives nearly the same trajectory. (b) Repeat the ODE solution using polar coordinates for the same initial conditions. (c) Finally, compute the solution using the methods describes in the previous problem, and compare with the ODE solutions Sep 10/MDW 4
5 (d) If any of these results are substantially different from each other, try to figure out why... [Hint: This is not a trick question. All methods should give the same trajectory!] 5 Making sense of the Poisson equation For a single point mass, the Poisson equation tells us that: 2 Φ(r) = ( GM r ) r 3 = 4πGMδ 3 (r). On the other hand, away for the origin, we know that ( r ) r 3 = 0. This tension can be explained by the singularity at r = 0 and motivates the definition of the Dirac delta function. This whole thing is a bit easier to swallow if we deal with non-singular potentials. We can remove the singularity by softening the point mass: (a) Show explicitly that where S is a spherical volume. x F ε (x) = (x 2 + ε 2 ) 3/2. lim F ε = 4π ε 0 S (b) Similarly, show that lim F ε 0 ε 0 for x 0. This expression, which is a density, is a representation of the Dirac delta function Sep 10/MDW 5
The two body problem involves a pair of particles with masses m 1 and m 2 described by a Lagrangian of the form:
Physics 3550, Fall 2011 Two Body, Central-Force Problem Relevant Sections in Text: 8.1 8.7 Two Body, Central-Force Problem Introduction. I have already mentioned the two body central force problem several
More informationThe Two -Body Central Force Problem
The Two -Body Central Force Problem Physics W3003 March 6, 2015 1 The setup 1.1 Formulation of problem The two-body central potential problem is defined by the (conserved) total energy E = 1 2 m 1Ṙ2 1
More informationSymmetries 2 - Rotations in Space
Symmetries 2 - Rotations in Space This symmetry is about the isotropy of space, i.e. space is the same in all orientations. Thus, if we continuously rotated an entire system in space, we expect the system
More informationF = ma. G mm r 2. S center
In the early 17 th century, Kepler discovered the following three laws of planetary motion: 1. The planets orbit around the sun in an ellipse with the sun at one focus. 2. As the planets orbit around the
More informationTwo dimensional oscillator and central forces
Two dimensional oscillator and central forces September 4, 04 Hooke s law in two dimensions Consider a radial Hooke s law force in -dimensions, F = kr where the force is along the radial unit vector and
More informationLecture 6, September 1, 2017
Engineering Mathematics Fall 07 Lecture 6, September, 07 Escape Velocity Suppose we have a planet (or any large near to spherical heavenly body) of radius R and acceleration of gravity at the surface of
More information1 Summary of Chapter 2
General Astronomy (9:61) Fall 01 Lecture 7 Notes, September 10, 01 1 Summary of Chapter There are a number of items from Chapter that you should be sure to understand. 1.1 Terminology A number of technical
More informationLecture 22: Gravitational Orbits
Lecture : Gravitational Orbits Astronomers were observing the motion of planets long before Newton s time Some even developed heliocentric models, in which the planets moved around the sun Analysis of
More informationMotion under the Influence of a Central Force
Copyright 004 5 Motion under the Influence of a Central Force The fundamental forces of nature depend only on the distance from the source. All the complex interactions that occur in the real world arise
More informationUse conserved quantities to reduce number of variables and the equation of motion (EOM)
Physics 106a, Caltech 5 October, 018 Lecture 8: Central Forces Bound States Today we discuss the Kepler problem of the orbital motion of planets and other objects in the gravitational field of the sun.
More informationChapter 8. Orbits. 8.1 Conics
Chapter 8 Orbits 8.1 Conics Conic sections first studied in the abstract by the Greeks are the curves formed by the intersection of a plane with a cone. Ignoring degenerate cases (such as a point, or pairs
More informationCentral force motion/kepler problem. 1 Reducing 2-body motion to effective 1-body, that too with 2 d.o.f and 1st order differential equations
Central force motion/kepler problem This short note summarizes our discussion in the lectures of various aspects of the motion under central force, in particular, the Kepler problem of inverse square-law
More informationQuestion 1: Spherical Pendulum
Question 1: Spherical Pendulum Consider a two-dimensional pendulum of length l with mass M at its end. It is easiest to use spherical coordinates centered at the pivot since the magnitude of the position
More informationPhys 7221, Fall 2006: Midterm exam
Phys 7221, Fall 2006: Midterm exam October 20, 2006 Problem 1 (40 pts) Consider a spherical pendulum, a mass m attached to a rod of length l, as a constrained system with r = l, as shown in the figure.
More informationPhysics 202 Laboratory 3. Root-Finding 1. Laboratory 3. Physics 202 Laboratory
Physics 202 Laboratory 3 Root-Finding 1 Laboratory 3 Physics 202 Laboratory The fundamental question answered by this week s lab work will be: Given a function F (x), find some/all of the values {x i }
More information!! r r θ! 2 r!φ 2 sinθ = GM / r 2 r!! θ + 2!r θ! r!φ 2 sinθ cosθ = 0 r!! φ sinθ + 2r θ!!φ cosθ + 2!r!φ sinθ = 0
Phys 325 Lecture 9 Tuesday 12 February, 2019 NB Midterm next Thursday 2/21. Sample exam will be posted soon Formula sheet on sample exam will be same as formula sheet on actual exam You will be permitted
More informationa n b n ) n N is convergent with b n is convergent.
32 Series Let s n be the n-th partial sum of n N and let t n be the n-th partial sum of n N. For k n we then have n n s n s k = a i a i = t n t k. i=k+1 i=k+1 Since t n n N is convergent by assumption,
More informationLecture 17 - Gyroscopes
Lecture 17 - Gyroscopes A Puzzle... We have seen throughout class that the center of mass is a very powerful tool for evaluating systems. However, don t let yourself get carried away with how useful it
More informationM2A2 Problem Sheet 3 - Hamiltonian Mechanics
MA Problem Sheet 3 - Hamiltonian Mechanics. The particle in a cone. A particle slides under gravity, inside a smooth circular cone with a vertical axis, z = k x + y. Write down its Lagrangian in a) Cartesian,
More informationFinal Examination Solutions
Math. 5, Sections 5 53 (Fulling) 7 December Final Examination Solutions Test Forms A and B were the same except for the order of the multiple-choice responses. This key is based on Form A. Name: Section:
More informationPhysics 6010, Fall 2016 Constraints and Lagrange Multipliers. Relevant Sections in Text:
Physics 6010, Fall 2016 Constraints and Lagrange Multipliers. Relevant Sections in Text: 1.3 1.6 Constraints Often times we consider dynamical systems which are defined using some kind of restrictions
More informationPHY321 Homework Set 10
PHY321 Homework Set 10 1. [5 pts] A small block of mass m slides without friction down a wedge-shaped block of mass M and of opening angle α. Thetriangular block itself slides along a horizontal floor,
More informationGeneral Technical Remarks on PDE s and Boundary Conditions Kurt Bryan MA 436
General Technical Remarks on PDE s and Boundary Conditions Kurt Bryan MA 436 1 Introduction You may have noticed that when we analyzed the heat equation on a bar of length 1 and I talked about the equation
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 informationLecture 41: Highlights
Lecture 41: Highlights The goal of this lecture is to remind you of some of the key points that we ve covered this semester Note that this is not the complete set of topics that may appear on the final
More informationPhysics 351, Spring 2017, Homework #4. Due at start of class, Friday, February 10, 2017
Physics 351, Spring 2017, Homework #4. Due at start of class, Friday, February 10, 2017 Course info is at positron.hep.upenn.edu/p351 When you finish this homework, remember to visit the feedback page
More informationTwo-Body Problem. Central Potential. 1D Motion
Two-Body Problem. Central Potential. D Motion The simplest non-trivial dynamical problem is the problem of two particles. The equations of motion read. m r = F 2, () We already know that the center of
More informationPhysics 235 Chapter 7. Chapter 7 Hamilton's Principle - Lagrangian and Hamiltonian Dynamics
Chapter 7 Hamilton's Principle - Lagrangian and Hamiltonian Dynamics Many interesting physics systems describe systems of particles on which many forces are acting. Some of these forces are immediately
More informationSolution Set Five. 2 Problem #2: The Pendulum of Doom Equation of Motion Simple Pendulum Limit Visualization...
: Solution Set Five Northwestern University, Classical Mechanics Classical Mechanics, Third Ed.- Goldstein November 4, 2015 Contents 1 Problem #1: Coupled Mass Oscillator System. 2 1.1 Normal Modes.......................................
More informationPhysics 351, Spring 2015, Homework #3. Due at start of class, Friday, February 6, 2015
Physics 351, Spring 2015, Homework #3. Due at start of class, Friday, February 6, 2015 Course info is at positron.hep.upenn.edu/p351 When you finish this homework, remember to visit the feedback page at
More informationClass Meeting # 1: Introduction to PDEs
MATH 18.152 COURSE NOTES - CLASS MEETING # 1 18.152 Introduction to PDEs, Spring 2017 Professor: Jared Speck Class Meeting # 1: Introduction to PDEs 1. What is a PDE? We will be studying functions u =
More informationPhysics 200 Lecture 4. Integration. Lecture 4. Physics 200 Laboratory
Physics 2 Lecture 4 Integration Lecture 4 Physics 2 Laboratory Monday, February 21st, 211 Integration is the flip-side of differentiation in fact, it is often possible to write a differential equation
More informationPhysics 351 Wednesday, February 14, 2018
Physics 351 Wednesday, February 14, 2018 HW4 due Friday. For HW help, Bill is in DRL 3N6 Wed 4 7pm. Grace is in DRL 2C2 Thu 5:30 8:30pm. Respond at pollev.com/phys351 or text PHYS351 to 37607 once to join,
More informationMATH1231 CALCULUS. Session II Dr John Roberts (based on notes of A./Prof. Bruce Henry) Red Center Room 3065
MATH1231 CALCULUS Session II 2007. Dr John Roberts (based on notes of A./Prof. Bruce Henry) Red Center Room 3065 Jag.Roberts@unsw.edu.au MATH1231 CALCULUS p.1/66 Overview Systematic Integration Techniques
More information2.20 Fall 2018 Math Review
2.20 Fall 2018 Math Review September 10, 2018 These notes are to help you through the math used in this class. This is just a refresher, so if you never learned one of these topics you should look more
More information2 4πε ( ) ( r θ. , symmetric about the x-axis, as shown in Figure What is the electric field E at the origin O?
p E( r, θ) = cosθ 3 ( sinθ ˆi + cosθ ˆj ) + sinθ cosθ ˆi + ( cos θ 1) ˆj r ( ) ( p = cosθ sinθ ˆi + cosθ ˆj + sinθ cosθ ˆi sinθ ˆj 3 r where the trigonometric identit ( θ ) vectors ˆr and cos 1 = sin θ
More informationSolutions 2: Simple Harmonic Oscillator and General Oscillations
Massachusetts Institute of Technology MITES 2017 Physics III Solutions 2: Simple Harmonic Oscillator and General Oscillations Due Wednesday June 21, at 9AM under Rene García s door Preface: This problem
More informationChapter 2 Introduction to Binary Systems
Chapter 2 Introduction to Binary Systems In order to model stars, we must first have a knowledge of their physical properties. In this chapter, we describe how we know the stellar properties that stellar
More informationPHYS 404 Lecture 1: Legendre Functions
PHYS 404 Lecture 1: Legendre Functions Dr. Vasileios Lempesis PHYS 404 - LECTURE 1 DR. V. LEMPESIS 1 Legendre Functions physical justification Legendre functions or Legendre polynomials are the solutions
More informationLecture 15 - Orbit Problems
Lecture 15 - Orbit Problems A Puzzle... The ellipse shown below has one focus at the origin and its major axis lies along the x-axis. The ellipse has a semimajor axis of length a and a semi-minor axis
More informationInvariants and Symmetry
November 007 Terminology Page 1 of 8 Invariants and Symmetry An invariant is some property of an object that does not change even when the object is changed or its description is changed. A semi-invariant
More informationPotential/density pairs and Gauss s law
Potential/density pairs and Gauss s law We showed last time that the motion of a particle in a cluster will evolve gradually, on the relaxation time scale. This time, however, is much longer than the typical
More informationCentral Force Problem
Central Force Problem Consider two bodies of masses, say earth and moon, m E and m M moving under the influence of mutual gravitational force of potential V(r). Now Langangian of the system is where, µ
More information2016 FAMAT Convention Mu Integration 1 = 80 0 = 80. dx 1 + x 2 = arctan x] k2
6 FAMAT Convention Mu Integration. A. 3 3 7 6 6 3 ] 3 6 6 3. B. For quadratic functions, Simpson s Rule is eact. Thus, 3. D.. B. lim 5 3 + ) 3 + ] 5 8 8 cot θ) dθ csc θ ) dθ cot θ θ + C n k n + k n lim
More informationSolutions for the Practice Final - Math 23B, 2016
olutions for the Practice Final - Math B, 6 a. True. The area of a surface is given by the expression d, and since we have a parametrization φ x, y x, y, f x, y with φ, this expands as d T x T y da xy
More informationAST1100 Lecture Notes
AST1100 Lecture Notes 5 The virial theorem 1 The virial theorem We have seen that we can solve the equation of motion for the two-body problem analytically and thus obtain expressions describing the future
More informationCelestial Mechanics II. Orbital energy and angular momentum Elliptic, parabolic and hyperbolic orbits Position in the orbit versus time
Celestial Mechanics II Orbital energy and angular momentum Elliptic, parabolic and hyperbolic orbits Position in the orbit versus time Orbital Energy KINETIC per unit mass POTENTIAL The orbital energy
More informationNote: Final Exam is at 10:45 on Tuesday, 5/3/11 (This is the Final Exam time reserved for our labs). From Practice Test I
MA Practice Final Answers in Red 4/8/ and 4/9/ Name Note: Final Exam is at :45 on Tuesday, 5// (This is the Final Exam time reserved for our labs). From Practice Test I Consider the integral 5 x dx. Sketch
More informationTHE NATIONAL UNIVERSITY OF IRELAND, CORK COLÁISTE NA hollscoile, CORCAIGH UNIVERSITY COLLEGE, CORK. Summer Examination 2009.
OLLSCOIL NA héireann, CORCAIGH THE NATIONAL UNIVERSITY OF IRELAND, CORK COLÁISTE NA hollscoile, CORCAIGH UNIVERSITY COLLEGE, CORK Summer Examination 2009 First Engineering MA008 Calculus and Linear Algebra
More informationTwo- and Three-Dimensional Motion (Symon Chapter Three)
Two- and Three-Dimensional Motion (Symon Chapter Three) Physics A300 Fall 2003 Contents 1 Vector Analysis 2 1.1 Definitions..................................... 2 1.1.1 Notation..................................
More informationSpherical Coordinates and Legendre Functions
Spherical Coordinates and Legendre Functions Spherical coordinates Let s adopt the notation for spherical coordinates that is standard in physics: φ = longitude or azimuth, θ = colatitude ( π 2 latitude)
More informationLecture 21 Gravitational and Central Forces
Lecture 21 Gravitational and Central Forces 21.1 Newton s Law of Universal Gravitation According to Newton s Law of Universal Graviation, the force on a particle i of mass m i exerted by a particle j of
More information4. Complex Oscillations
4. Complex Oscillations The most common use of complex numbers in physics is for analyzing oscillations and waves. We will illustrate this with a simple but crucially important model, the damped harmonic
More information3. On the grid below, sketch and label graphs of the following functions: y = sin x, y = cos x, and y = sin(x π/2). π/2 π 3π/2 2π 5π/2
AP Physics C Calculus C.1 Name Trigonometric Functions 1. Consider the right triangle to the right. In terms of a, b, and c, write the expressions for the following: c a sin θ = cos θ = tan θ =. Using
More informationTwo- and Three-Dimensional Motion (Symon Chapter Three)
Two- and Three-Dimensional Motion (Symon Chapter Three) Physics A300 Spring 2006 Contents 1 Vector Analysis 3 11 Definitions 3 111 Notation 3 12 Basic Vector Arithmetic 3 121 Equality 3 122 Magnitude 4
More informationCurves in the configuration space Q or in the velocity phase space Ω satisfying the Euler-Lagrange (EL) equations,
Physics 6010, Fall 2010 Hamiltonian Formalism: Hamilton s equations. Conservation laws. Reduction. Poisson Brackets. Relevant Sections in Text: 8.1 8.3, 9.5 The Hamiltonian Formalism We now return to formal
More information18.02 Multivariable Calculus Fall 2007
MIT OpenourseWare http://ocw.mit.edu 8.02 Multivariable alculus Fall 2007 For information about citing these materials or our Terms of Use, visit: http://ocw.mit.edu/terms. 8.02 Lecture 8. hange of variables.
More informationChapter 13. Gravitation
Chapter 13 Gravitation 13.2 Newton s Law of Gravitation Here m 1 and m 2 are the masses of the particles, r is the distance between them, and G is the gravitational constant. G =6.67 x10 11 Nm 2 /kg 2
More informationStrauss PDEs 2e: Section Exercise 1 Page 1 of 6
Strauss PDEs 2e: Section 3 - Exercise Page of 6 Exercise Carefully derive the equation of a string in a medium in which the resistance is proportional to the velocity Solution There are two ways (among
More informationFalse. 1 is a number, the other expressions are invalid.
Ma1023 Calculus III A Term, 2013 Pseudo-Final Exam Print Name: Pancho Bosphorus 1. Mark the following T and F for false, and if it cannot be determined from the given information. 1 = 0 0 = 1. False. 1
More informationCalculus III. George Voutsadakis 1. LSSU Math 251. Lake Superior State University. 1 Mathematics and Computer Science
Calculus III George Voutsadakis 1 1 Mathematics and Computer Science Lake Superior State University LSSU Math 251 George Voutsadakis (LSSU) Calculus III January 2016 1 / 76 Outline 1 Parametric Equations,
More informationy = x 3 and y = 2x 2 x. 2x 2 x = x 3 x 3 2x 2 + x = 0 x(x 2 2x + 1) = 0 x(x 1) 2 = 0 x = 0 and x = (x 3 (2x 2 x)) dx
Millersville University Name Answer Key Mathematics Department MATH 2, Calculus II, Final Examination May 4, 2, 8:AM-:AM Please answer the following questions. Your answers will be evaluated on their correctness,
More information1MA6 Partial Differentiation and Multiple Integrals: I
1MA6/1 1MA6 Partial Differentiation and Multiple Integrals: I Dr D W Murray Michaelmas Term 1994 1. Total differential. (a) State the conditions for the expression P (x, y)dx+q(x, y)dy to be the perfect
More informationParticles in Motion; Kepler s Laws
CHAPTER 4 Particles in Motion; Kepler s Laws 4.. Vector Functions Vector notation is well suited to the representation of the motion of a particle. Fix a coordinate system with center O, and let the position
More information1 Exponential Functions Limit Derivative Integral... 5
Contents Eponential Functions 3. Limit................................................. 3. Derivative.............................................. 4.3 Integral................................................
More informationPHY752, Fall 2016, Assigned Problems
PHY752, Fall 26, Assigned Problems For clarification or to point out a typo (or worse! please send email to curtright@miami.edu [] Find the URL for the course webpage and email it to curtright@miami.edu
More informationMATH 1080 Test 2 -Version A-SOLUTIONS Fall a. (8 pts) Find the exact length of the curve on the given interval.
MATH 8 Test -Version A-SOLUTIONS Fall 4. Consider the curve defined by y = ln( sec x), x. a. (8 pts) Find the exact length of the curve on the given interval. sec x tan x = = tan x sec x L = + tan x =
More informationeff (r) which contains the influence of angular momentum. On the left is
1 Fig. 13.1. The radial eigenfunctions R nl (r) of bound states in a square-well potential for three angular-momentum values, l = 0, 1, 2, are shown as continuous lines in the left column. The form V (r)
More informationCelestial Mechanics Lecture 10
Celestial Mechanics Lecture 10 ˆ This is the first of two topics which I have added to the curriculum for this term. ˆ We have a surprizing amount of firepower at our disposal to analyze some basic problems
More informationThe answers below are not comprehensive and are meant to indicate the correct way to solve the problem. sin
Math : Practice Final Answer Key Name: The answers below are not comprehensive and are meant to indicate the correct way to solve the problem. Problem : Consider the definite integral I = 5 sin ( ) d.
More informationIn this chapter we study elliptical PDEs. That is, PDEs of the form. 2 u = lots,
Chapter 8 Elliptic PDEs In this chapter we study elliptical PDEs. That is, PDEs of the form 2 u = lots, where lots means lower-order terms (u x, u y,..., u, f). Here are some ways to think about the physical
More informationGraphical Analysis and Errors MBL
Graphical Analysis and Errors MBL I Graphical Analysis Graphs are vital tools for analyzing and displaying data Graphs allow us to explore the relationship between two quantities -- an independent variable
More informationThe first term involves the cross product of two parallel vectors and so it vanishes. We then get
Physics 3550 Angular Momentum. Relevant Sections in Text: 3.4, 3.5 Angular Momentum You have certainly encountered angular momentum in a previous class. The importance of angular momentum lies principally
More information3.4 Conic sections. Such type of curves are called conics, because they arise from different slices through a cone
3.4 Conic sections Next we consider the objects resulting from ax 2 + bxy + cy 2 + + ey + f = 0. Such type of curves are called conics, because they arise from different slices through a cone Circles belong
More informationSpace flying vehicles orbital motion control system synthesis: power invariants
CHAOS 010, 3rd Chaotic Modeling and Simulation 1-4 June 010 Space flying vehicles orbital motion control system synthesis: power invariants Alexander A. Kolesnikov (PhD) Taganrog Institute of Technology
More informationP321(b), Assignement 1
P31(b), Assignement 1 1 Exercise 3.1 (Fetter and Walecka) a) The problem is that of a point mass rotating along a circle of radius a, rotating with a constant angular velocity Ω. Generally, 3 coordinates
More informationPhysics 6010, Fall Relevant Sections in Text: Introduction
Physics 6010, Fall 2016 Introduction. Configuration space. Equations of Motion. Velocity Phase Space. Relevant Sections in Text: 1.1 1.4 Introduction This course principally deals with the variational
More informationGo over the illustrated examples in each section.
Math 242 Fall 2009 Please before you start the practice problems, go over the illustrated examples in each section first. Cover up the solutions and try to work out the answers on your own. Practice Problems
More informationMore Final Practice Problems
8.0 Calculus Jason Starr Final Exam at 9:00am sharp Fall 005 Tuesday, December 0, 005 More 8.0 Final Practice Problems Here are some further practice problems with solutions for the 8.0 Final Exam. Many
More informationMA 114 Worksheet #01: Integration by parts
Fall 8 MA 4 Worksheet Thursday, 3 August 8 MA 4 Worksheet #: Integration by parts. For each of the following integrals, determine if it is best evaluated by integration by parts or by substitution. If
More informationAssignment 2 - Complex Analysis
Assignment 2 - Complex Analysis MATH 440/508 M.P. Lamoureux Sketch of solutions. Γ z dz = Γ (x iy)(dx + idy) = (xdx + ydy) + i Γ Γ ( ydx + xdy) = (/2)(x 2 + y 2 ) endpoints + i [( / y) y ( / x)x]dxdy interiorγ
More informationMcGill University April 16, Advanced Calculus for Engineers
McGill University April 16, 2014 Faculty of cience Final examination Advanced Calculus for Engineers Math 264 April 16, 2014 Time: 6PM-9PM Examiner: Prof. R. Choksi Associate Examiner: Prof. A. Hundemer
More informationChapter 13. Gravitation
Chapter 13 Gravitation e = c/a A note about eccentricity For a circle c = 0 à e = 0 a Orbit Examples Mercury has the highest eccentricity of any planet (a) e Mercury = 0.21 Halley s comet has an orbit
More informationfor changing independent variables. Most simply for a function f(x) the Legendre transformation f(x) B(s) takes the form B(s) = xs f(x) with s = df
Physics 106a, Caltech 1 November, 2018 Lecture 10: Hamiltonian Mechanics I The Hamiltonian In the Hamiltonian formulation of dynamics each second order ODE given by the Euler- Lagrange equation in terms
More informationGravity 3. Gravity 3. Gravitational Potential and the Geoid. Chuck Connor, Laura Connor. Potential Fields Geophysics: Week 2.
Gravitational Potential and the Geoid Chuck Connor, Laura Connor Potential Fields Geophysics: Week 2 Objectives for Week 1 Gravity as a vector Gravitational Potential The Geoid Gravity as a vector We can
More informationChaos in the Hénon-Heiles system
Chaos in the Hénon-Heiles system University of Karlstad Christian Emanuelsson Analytical Mechanics FYGC04 Abstract This paper briefly describes how the Hénon-Helies system exhibits chaos. First some subjects
More information21 Laplace s Equation and Harmonic Functions
2 Laplace s Equation and Harmonic Functions 2. Introductory Remarks on the Laplacian operator Given a domain Ω R d, then 2 u = div(grad u) = in Ω () is Laplace s equation defined in Ω. If d = 2, in cartesian
More informationConstrained motion and generalized coordinates
Constrained motion and generalized coordinates based on FW-13 Often, the motion of particles is restricted by constraints, and we want to: work only with independent degrees of freedom (coordinates) k
More informationMath 1272 Solutions for Fall 2005 Final Exam
Math 7 Solutions for Fall 5 Final Exam ) This fraction appears in Problem 5 of the undated-? exam; a solution can be found in that solution set. (E) ) This integral appears in Problem 6 of the Fall 4 exam;
More informationISIMA lectures on celestial mechanics. 1
ISIMA lectures on celestial mechanics. 1 Scott Tremaine, Institute for Advanced Study July 2014 The roots of solar system dynamics can be traced to two fundamental discoveries by Isaac Newton: first, that
More informationLecture 6. Angular Momentum and Roataion Invariance of L. x = xcosφ + ysinφ ' y = ycosφ xsinφ for φ << 1, sinφ φ, cosφ 1 x ' x + yφ y ' = y xφ,
Lecture 6 Conservation of Angular Momentum (talk about LL p. 9.3) In the previous lecture we saw that time invariance of the L leads to Energy Conservation, and translation invariance of L leads to Momentum
More informationL = 1 2 a(q) q2 V (q).
Physics 3550, Fall 2011 Motion near equilibrium - Small Oscillations Relevant Sections in Text: 5.1 5.6 Motion near equilibrium 1 degree of freedom One of the most important situations in physics is motion
More informationRotational motion of a rigid body spinning around a rotational axis ˆn;
Physics 106a, Caltech 15 November, 2018 Lecture 14: Rotations The motion of solid bodies So far, we have been studying the motion of point particles, which are essentially just translational. Bodies with
More informationKepler s laws. M. Plakhotnyk. March 2, 2017
Kepler s laws arxiv:1702.06537v2 [math.ho] 1 Mar 2017 M. Plakhotnyk March 2, 2017 Kepler s laws with introduction to differential calculus (book for scholars, who are interested in physics and mathematics)
More informationIdeas from Vector Calculus Kurt Bryan
Ideas from Vector Calculus Kurt Bryan Most of the facts I state below are for functions of two or three variables, but with noted exceptions all are true for functions of n variables..1 Tangent Line Approximation
More informationMITOCW watch?v=fxlzy2l1-4w
MITOCW watch?v=fxlzy2l1-4w PROFESSOR: We spoke about the hydrogen atom. And in the hydrogen atom, we drew the spectrum, so the table, the data of spectrum of a quantum system. So this is a question that
More information4.10 Dirichlet problem in the circle and the Poisson kernel
220 CHAPTER 4. FOURIER SERIES AND PDES 4.10 Dirichlet problem in the circle and the Poisson kernel Note: 2 lectures,, 9.7 in [EP], 10.8 in [BD] 4.10.1 Laplace in polar coordinates Perhaps a more natural
More informationFact: Every matrix transformation is a linear transformation, and vice versa.
Linear Transformations Definition: A transformation (or mapping) T is linear if: (i) T (u + v) = T (u) + T (v) for all u, v in the domain of T ; (ii) T (cu) = ct (u) for all scalars c and all u in the
More informationAPPM 1360 Final Exam Spring 2016
APPM 36 Final Eam Spring 6. 8 points) State whether each of the following quantities converge or diverge. Eplain your reasoning. a) The sequence a, a, a 3,... where a n ln8n) lnn + ) n!) b) ln d c) arctan
More informationLagrangian for Central Potentials
Physics 411 Lecture 2 Lagrangian for Central Potentials Lecture 2 Physics 411 Classical Mechanics II August 29th 2007 Here we will review the Lagrange formulation in preparation for the study of the central
More information