The Two -Body Central Force Problem
|
|
- Tracey Norris
- 6 years ago
- Views:
Transcription
1 The Two -Body Central Force Problem Physics W3003 March 6, 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Ṙ m 1Ṙ2 1 + U( R 1 R 2 ) (1) This problem has 6 degrees of freedom (three components of position for each particle). The energy is conserved. The energy is invariant under translation in the three Cartesian directions, thus the the three components of total momentum are conserved. The energy is invariant under rotations about an arbitrary axis, so the three components of angular momentum are conserved. These conservation laws allow us to construct a complete formal solution. These notes take you through the solution, first in the general case (arbitrary potential) and then in the important specific case (gravitational or Coulomb) of U(r) 1/r. 1.2 Separation of motion into center of mass and relative coordinates Define the center of mass coordinate R COM and relative coordinate r via These equations imply R COM = m 1R 1 + m 2 R 2 m 1 + m 2 (2) r = R 1 R 2 (3) m 2 R 1 = R COM + r (4) m 1 + m 2 m 1 R 2 = R COM r (5) m 1 + m 2 1
2 Substituting Eqs. 4,5 into Eq. 1 gives with r = r and the reduced mass µ defined by E = 1 2 (m 1 + m 2 ) Ṙ2 COM µṙ2 + U(r) (6) µ = m 1m 2 m 1 + m 2 (7) Thus the total energy is the sum of a term corresponding to the free motion of a particle with mass equal to the total mass of the system and position equal to the position of the center of mass, and a second term corresponding to the motion of a fictitious second particle which represents the relative motion of the two particles (in other words, the internal degrees of freedom of the system). This fictitious particle feels a force arising from the interparticle potential U and has mass equal to the reduced mass. Conservation of center of mass momentum means that ṘCOM is independent of time so we may write down the solution for the motion of the center of mass as R COM = R init COM + V COM t (8) with initial position R init COM and velocity V COM determined by boundary conditions. From this point on we focus on the relative motion of the two particles 1.3 Simplifying the description of the relative motion Because the center of mass velocity is constant, the energy associated with the center of mass motion (the kinetic energy of the center of mass) is conserved. Therefore the energy of internal motion E int = E 1 2 (m 1 + m 2 ) Ṙ2 COM = 1 2 µṙ2 + U(r) (9) is conserved. The potential U(r) 0 in general so the momentum p = µṙ associated with internal motion is not conserved. However we assume that U(r) depends only on the magnitude of the coordinate r so that it is rotationally invariant. Thus the angular momentum L = r p associated with relative motion is conserved. Conservation of angular momentum has an important consequence. At any given time t the two vectors r and p define a plane and L is perpendicular to this plane. The conservation of L implies that L always remains parallel to its initial direction; thus r and p must also always lie in the same plane (if they did not, there would be a component of L perpendicular to its initial direction. Take the plane defined by r and p to be the x y plane; then the non-zero component of L is the z component L z. Thus the relative motion problem has only two degrees of 2
3 freedom, and has two conserved quantities, E int and L = L = L z, so the problem is integrable (the solution may be written down in terms of definite integrals). It is easiest to write down the solution in polar coordinates. Define Then the chain rule implies and Thus if r(t) is known, θ(t) may be found from implying x = r cos θ (10) y = r sin θ (11) ẋ = ṙ cos θ θ r sin θ (12) ẏ = ṙ sin θ + θ r sin θ (13) L = x(µẏ) y(µẋ) = µ r 2 θ (14) θ = θ(t) θ(t = 0) = L mr 2 (15) t 0 dt L µr 2 (t ) Now, substituting Eqs. 12 and 13 into Eq. 9 for the internal energy gives Finally using Eq. 15 gives (16) E internal = 1 2 µṙ µr2 θ2 + U(r) (17) E internal = 1 2 µṙ2 + L2 + U(r) (18) 2µr2 Thus the time dependence of the radial coordinate is determined by solving the problem of a particle of mass µ moving in the effective potential U eff = L2 + U(r) (19) 2µr2 given by the sum of the physical potential U(r) and the centrifugal potential L 2 /(2µr 2 ). The centrifugal potential blows up as r 0 and thus acts to keep the particle away from r = 0. It expresses the physics that as r decreases, the angular velocity θ must increase rapidly to keep L constant. Thus as r decreases more of the energy goes into rotational motion, leaving less for the radial motion. 3
4 1.4 Formal Solution From Eq. 18 we have (using the definition Eq. 19) 2 ṙ = ± E internal U eff (r) (20) µ Thus t(r) t(r init ) = µ r dr 2 r init Einternal U eff (r ) provides a formal solution for time as a function of position. Inverting this gives r(t); use of Eq. 16 then gives θ(t) completing the solution. Note also that use of Eq. 20 in Eq. 16 gives the explicit relation (21) θ(r) θ(r init ) = L r dr 2µ r init (r ) 2 E internal U eff (r) (22) between position and angle that specifies the orbit. The rest of these notes examine the physics contained in this formal solution. We begin by looking at the effective potential. 2 The effective potential The effective potential U eff (r) is constructed by adding a repulsive term L 2 /(2µr 2 ) to the original interparticle potential U(r). The result depends on the form of U(r). It is useful to distinguish three cases Repulsive U(r). If U(r) 0 for all r then U eff (r) 0 for all r and we simply have a repulsive potential which admits only unbound orbits. An example is sketched in panel a of Fig. 1. Attractive, short-ranged U(r). If U(r) 0 but is short ranged (goes to zero faster than 1/r 2 as r ) and is not divergent (does not increase as fast as 1/r 2 as r 0) then we have the interesting situation sketched in panel (b) of Fig. 1: for small L the potential has a minimum, and thus there can be both bound orbits and unbound orbits, but if L is too large, then there is no minimum and there are only unbound orbits. Gravitational: if U(r) = k/r then for any L at very large r the potential is negative, decreasing as r decreases, but at small r the potential is positive and increasing as r increases, so there is always a minimum and both bound and unbound orbits, as shown in Panel c of Fig. 1 4
5 Figure 1: Left panel: effective potential (solid line) constructed from short ranged repulsive potential U(r) = e r (dashed line) and centrifugal potential computed with µ = 1 and L = 1. Central panel: effective potential (solid line) constructed from short ranged attractive potential U(r) = e r (dashed line) and centrifugal potential computed with µ1 and L = 0.25 (light solid line) and L = 1.1 (heavy line) showing cases with and without minimum. Right panel: effective potential (solid line) constructed from attractive 1/r potential U(r) = 1/r (dashed line) and centrifugal potential computed with µ = 1 and L = (light line) and 0.7 (heavy line). 3 Unbound Orbits The motion in an unbound orbit involves two particles, initially very far from each other. The two particles approach, have their motion affected by the interparticle potential, and then move far away, with the relative motion in a different direction than it was initially. When the particles are far from each other the relative motion is of a free particle r(t) = r init + vt (23) Careful!: these equations only describe the initial state of relative motion, when the particles are far apart. The angular momentum and internal energy are L = µr v = µr init v E internal = µ 2 v2 (24) Let us choose coordinates so that the initial velocity of relative motion is along x, so the initial motion is x(t) = x 0 vt (25) y = y 0 (26) The energy and angular momentum are E internal = µ 2 v2 and L = µy 0 v. If U(r) = 0 the motion would be described by Eqs. 25,26 at all times, and the distance of closest approach r min would be y 0. Let us now describe the motion in radial and angular coordinates. If the potential U(r) = 0 then U ef (r) = L 2 /(2µr 2 ) and the distance of closest approach r min would be 5
6 given by the solution of E internal = U eff (r min ) (27) or (using the expressions for E internal and L in terms of the initial velocity) or so that y Figure 2: Trajectories in xy plane computed as described in text for potentials used in panels (a) and (b) of Fig. 1 with E = 1.5 and L = 1. Upper line corresponds to repulsive interaction, lower to attractive and dashed to no interaction. Dot labels origin of relative coordinates µ 2 v2 = µ2 v 2 y0 2 2µrmin 2 + U(r min ) (28) r 2 min = y U(r min)r 2 min E initial (29) x r 2 min = y U(r min) E initial (30) so that if the physical potential U(r) = 0 then r min = y 0, while if the potential is attractive the minimum distance is smaller than y 0 and if repulsive, the minimum distance is larger than y 0. Thus for an attractive potential the trajectory is bent towards the origin of relative motion and if the potential is repulsive the trajectory is bent away from the origin of relative motion. We can now use Eq. 21 to compute the relation between time and radial position (bearing in mind that we use the negative square root (rdecreases) to describe the approach phase and the positive to describe the separation, and then Eq. 16 to get the angle. From these we can get the x and y positions. I used the software package Mathematica to do the integrals for the short ranged attractive and repulsive potentials shown in panels (a) and (b) of Fig. 1. I took µ = 1 and said that at time t = 0 the particle was at position r = 3 with angular momentum L = 1 and energy E = 1.5. This implies v = 3 so y 0 = 1/ 3 implying that the initial angle is θ init = ArcSin1/(3 3). Using this initial data I performed the integrals. Results are presented in Fig. 2. 6
7 4 Bound orbits 4.1 General properties In a bound orbit the particle is confined between a minimum radius r min and a maximum radius r max which are two solutions of E internal = U eff (r) (31) A very important special case is when E internal is equal to the minimum value of U eff. In this case the minimum and maximum radius are equal: r min = r max = r circ. In this case the centripetal force L 2 /mr 3 arising from the derivative of the centrifugal potential just balances the attractive force du/dr. Because ṙ = 0 the radius is constant in time and we have a circular orbit. For a circular orbit the time dependence of the phase is easy to find: Eq. 16 becomes θ(t) θ(t = 0) = Lt (32) µr 2 circ which is just the statement that for a circular orbit the angular momentum is L = mr 2 theta. If E internal is greater than the minimum value of U eff the radial coordinate oscillates between r min amd r max according to Eq. 21 and the phase advances according to Eq. 16. Panel (a) of Fig. 3 shows the variation with time of the radial coordinate for the attractive potential discussed above. Here I have chosen L = 0.25 and E = 0.2. Notice that r varies much faster near its minimum point because the potential varies much more rapidly at small r, so the force is much greater (cf Fig. 1). Panel (b) of Fig. 3 shows the variation with time of the phase for the same case. We see again that the phase varies more rapidly when r is small, in this case because of conservation of angular momentum. We now ask a simple question: how much time, T r, does it take for the relative coordinate to go from r min to r max and back to r min? From Eq. 21 we find (the 2 is because the total time is twice the time it takes to go from r min to r max. µ T r = 2 2 rmax dr r min Einternal U eff (r ) (33) In the case shown in Fig. 3 the time is about 7. We may now ask another simple question: by how much does the phase advance in the time T r. From Eq. 16 we find Tr θ(t r ) = θ(t = 0) + 0 dt L µr 2 (t ) In the case shown in Fig. 3 the phase advances only by a bit more than π. In general, θ(t r ) 2π; in other words, the system does not come back to the same state after the 7 (34)
8 (a) (c) (b) 1.5 (d) Figure 3: Panel (a) time dependence of relative coordinate computed using attractive potential of Fig 1 with E = 0.2 and L = Panel (b) Time dependence of phase computed for same parameters. (c) and (d) trajectories in xy plane computed from r and θ for a short time (t = 15) and a longer time (t = 50). radial coordinate goes from its minimum to its maximum value and back to its minimum. The location of the point of closest approach changes with time. This is a generic feature of problems with multiple degrees of freedom: if the motion is bounded, the degrees of freedom are periodic in time, but with as many periods as there are degrees of freedom, and generically the periods are incommensurate (not rational multiples of each other). In this case the system explores the entire region energetically accessible to it, as you can see in panels c and d of Fig. 3, which show the trajectories in the xy plane. As we will see in the next section, the gravitational potential is a special case. For the gravitational potential the time needed for the radial coordinate to complete a cycle is the same as the time needed for the phase to advance by 2π so the position of the point of closest approach remains the same and the orbits have a simple geometric structure, in fact an ellipse. 8
9 5 The gravitational potential The gravitational case has two remarkable accidental features. First, for bound orbits the time required for the relative coordinate to move from its minimum to its maximum value is the same as the time required for the phase to advance by π. The orbits therefore have a simple structure. Second, the integrals can be done analytically, so simple expressions can be obtained in terms of known functions. It is interesting to speculate on how different the history of science might have been if this were not the case. In the gravitational case, the potential is k/r with k determined by Newton s gravitational constant and the masses of the two bodies. Thus U eff (r) = L2 2µr 2 k r (35) At very small r, U eff 1/r 2 so is large and positive, and increasing as r decreases. At very large r, U eff 1/r so is negative and gets more negative as r decreases. U eff must therefore have at least one minimum. Setting the r-derivative to zero gives du eff (r) dr implying that the extremum of the function is at = L2 µr 3 + k r 2 = 0 (36) r circ = L2 µk (37) The second derivative is positive so this is indeed a minimum. The value of U eff at the minimum is U circ = µk2 2L 2 (38) If the total energy of radial motion is greater than or equal to U circ but less than 0 the motion is confined between the two solutions of E internal = L2 2µr 2 k r The quantity 1/r appears repeatedly, so it is convenient to use this as a variable. Define η = r circ r Then using Eq. 37 we can write (39) (40) e int = η 2 2η (41) 9
10 with e int = E internal U circ The boundaries of the orbit are then given by (42) η ± = 1 ± 1 + e int (43) with η corresponding to r max and η + corresponding to r min. Instead of evaluating the time dependence we use Eq. 22 to obtain the angle as a function of the radial coordinate. Starting the integration from the minimum value we have θ(r) = θ(r min ) + L r 2µ Expressing this in terms of η and e int gives dr r min (r ) 2 E internal L2 + k 2µr 2 r (44) Now define η dη θ(r) = θ(r min ) η max eint η 2 + 2η (45) η = e int cosψ (46) and note that at r = r min, θ = 0 and at r = r max, θ = π. Eq. 45 becomes θ r circ r 1 θ(r) = dψ = ArcCos (47) θ rmin 1 + eint Note in particular that Eq. 47 implies that θ increases by π when r goes from r min to r max so the particle comes back to exactly the same point after each cycle. Rearranging: or 1 r = 1 ( eint cosθ ) (48) r circ r = r 0 (1 + ɛ cosθ) This is the equation of an ellipse with semimajor axis r 0 = r circ = L 2 /µk and eccentricity ɛ = 1 + 2E internall 2. µk 2 (49) 10
Central 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 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 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 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 informationThe 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 informationPhysics 351, Spring 2015, Homework #7. Due at start of class, Friday, March 3, 2017
Physics 351, Spring 2015, Homework #7. Due at start of class, Friday, March 3, 2017 Course info is at positron.hep.upenn.edu/p351 When you finish this homework, remember to visit the feedback page at positron.hep.upenn.edu/q351
More informationPhysics 351, Spring 2015, Homework #7. Due at start of class, Friday, March 6, 2015
Physics 351, Spring 2015, Homework #7. Due at start of class, Friday, March 6, 2015 Course info is at positron.hep.upenn.edu/p351 When you finish this homework, remember to visit the feedback page at positron.hep.upenn.edu/q351
More informationThe Restricted 3-Body Problem
The Restricted 3-Body Problem John Bremseth and John Grasel 12/10/2010 Abstract Though the 3-body problem is difficult to solve, it can be modeled if one mass is so small that its effect on the other two
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 informationA645/A445: Exercise #1. The Kepler Problem
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
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 informationElastic Collisions. Chapter Center of Mass Frame
Chapter 11 Elastic Collisions 11.1 Center of Mass Frame A collision or scattering event is said to be elastic if it results in no change in the internal state of any of the particles involved. Thus, no
More informationFROM NEWTON TO KEPLER. One simple derivation of Kepler s laws from Newton s ones.
italian journal of pure and applied mathematics n. 3 04 (393 400) 393 FROM NEWTON TO KEPLER. One simple derivation of Kepler s laws from Newton s ones. František Mošna Department of Mathematics Technical
More informationPhysics 115/242 The Kepler Problem
Physics 115/242 The Kepler Problem Peter Young (Dated: April 21, 23) I. INTRODUCTION We consider motion of a planet around the sun, the Kepler problem, see e.g. Garcia, Numerical Methods for Physics, Sec.
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 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 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 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 informationThe first order formalism and the transition to the
Problem 1. Hamiltonian The first order formalism and the transition to the This problem uses the notion of Lagrange multipliers and Legendre transforms to understand the action in the Hamiltonian formalism.
More information3 Space curvilinear motion, motion in non-inertial frames
3 Space curvilinear motion, motion in non-inertial frames 3.1 In-class problem A rocket of initial mass m i is fired vertically up from earth and accelerates until its fuel is exhausted. The residual mass
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 informationSolution Set Two. 1 Problem #1: Projectile Motion Cartesian Coordinates Polar Coordinates... 3
: Solution Set Two Northwestern University, Classical Mechanics Classical Mechanics, Third Ed.- Goldstein October 7, 2015 Contents 1 Problem #1: Projectile Motion. 2 1.1 Cartesian Coordinates....................................
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 informationWe start by noting Kepler's observations (before Newton) regarding planetary motion.
Phys 325 Lecture 8 Thursday 7 February, 2019 Announcements: HW4 is longer than the others, and is due Tuesday 2/19 Midterm I will be in-class on Thursday 2/21. Material covered will be lectures 1-10 and
More informationWe have seen how to calculate forces and potentials from the smoothed density ρ.
9. Orbits in stationary Potentials We have seen how to calculate forces and potentials from the smoothed density ρ. Now we can analyse how stars move in this potential. Because two body interactions can
More informationMechanics Lecture Notes
Mechanics Lecture Notes Lectures 0 and : Motion in a circle. Introduction The important result in this lecture concerns the force required to keep a particle moving on a circular path: if the radius of
More informationMarion and Thornton. Tyler Shendruk October 1, Hamilton s Principle - Lagrangian and Hamiltonian dynamics.
Marion and Thornton Tyler Shendruk October 1, 2010 1 Marion and Thornton Chapter 7 Hamilton s Principle - Lagrangian and Hamiltonian dynamics. 1.1 Problem 6.4 s r z θ Figure 1: Geodesic on circular cylinder
More informationthe EL equation for the x coordinate is easily seen to be (exercise)
Physics 6010, Fall 2016 Relevant Sections in Text: 1.3 1.6 Examples After all this formalism it is a good idea to spend some time developing a number of illustrative examples. These examples represent
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 informationPhysics 106a, Caltech 4 December, Lecture 18: Examples on Rigid Body Dynamics. Rotating rectangle. Heavy symmetric top
Physics 106a, Caltech 4 December, 2018 Lecture 18: Examples on Rigid Body Dynamics I go through a number of examples illustrating the methods of solving rigid body dynamics. In most cases, the problem
More informationClassical Mechanics Comprehensive Exam Solution
Classical Mechanics Comprehensive Exam Solution January 31, 011, 1:00 pm 5:pm Solve the following six problems. In the following problems, e x, e y, and e z are unit vectors in the x, y, and z directions,
More informationMTHE 227 Problem Set 2 Solutions
MTHE 7 Problem Set Solutions 1 (Great Circles). The intersection of a sphere with a plane passing through its center is called a great circle. Let Γ be the great circle that is the intersection of the
More informationPhysics 105 (Fall 2013) Notes: The Two-Body Conservative Central Force Problem
Physics 5 (Fall 3) Notes: The Two-Body Conservative Central Force Problem A.E. Charman Department of Physics, University of California, Berkeley (Dated: /5/3; revised //3; revised /8/3) Here we analyze
More informationCircular motion. Aug. 22, 2017
Circular motion Aug. 22, 2017 Until now, we have been observers to Newtonian physics through inertial reference frames. From our discussion of Newton s laws, these are frames which obey Newton s first
More informationElastic Scattering. R = m 1r 1 + m 2 r 2 m 1 + m 2. is the center of mass which is known to move with a constant velocity (see previous lectures):
Elastic Scattering In this section we will consider a problem of scattering of two particles obeying Newtonian mechanics. The problem of scattering can be viewed as a truncated version of dynamic problem
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 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 informationThe central force problem
1 The central force problem Moment. Das Moment des Eindrucks, den ein Mann auf das gemeine Volk macht, ist ein Produkt aus dem Wert des Rocks in den Titel. Georg Christoph Lichtenberg We start by dealing
More informationPHYSICS I. Lecture 1. Charudatt Kadolkar. Jul-Nov IIT Guwahati
PHYSICS I Lecture 1 Charudatt Kadolkar IIT Guwahati Jul-Nov 2014 Section 1 Introduction to the Course Syllabus Topics Classical Mechanics: Kinetic Energy rest mass energy Syllabus Topics Classical Mechanics:
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 informationProblem Set 5 Solution
Problem Set 5 Solution Friday, 14 October 011 Physics 111 Proble The Setup For a pair of point masses, and m, interacting via a conservative force directed along the line joining the masses, show that
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 informationLecture Notes for PHY 405 Classical Mechanics
Lecture Notes for PHY 405 Classical Mechanics From Thorton & Marion s Classical Mechanics Prepared by Dr. Joseph M. Hahn Saint Mary s University Department of Astronomy & Physics September 1, 2005 Chapter
More informationPhysics 351, Spring 2015, Homework #6. Due at start of class, Friday, February 27, 2015
Physics 351, Spring 2015, Homework #6. Due at start of class, Friday, February 27, 2015 Course info is at positron.hep.upenn.edu/p351 When you finish this homework, remember to visit the feedback page
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 informationPhysics 20 Homework 3 SIMS 2016
Physics 20 Homework 3 SIMS 2016 Due: Thursday, August 25 th Special thanks to Sebastian Fischetti for problems 1, 5, and 6. Edits in red made by Keith Fratus. 1. The ballistic pendulum is a device used
More informationPhysics 103 Homework 5 Summer 2015
Physics 103 Homework 5 Summer 2015 Instructor: Keith Fratus TAs: Michael Rosenthal and Milind Shyani Due Date: Friday, July 31 st, IN LECTURE 1. A particle moving in three dimensions is constrained to
More information21.60 Worksheet 8 - preparation problems - question 1:
Dynamics 190 1.60 Worksheet 8 - preparation problems - question 1: A particle of mass m moes under the influence of a conseratie central force F (r) =g(r)r where r = xˆx + yŷ + zẑ and r = x + y + z. A.
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 informationUniversal Gravitation
Universal Gravitation Newton s Law of Universal Gravitation Every particle in the Universe attracts every other particle with a force that is directly proportional to the product of their masses and inversely
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 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 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 informationCentral Forces. c 2007 Corinne Manogue. March 1999, last revised February 2009
Central Forces Corinne A. Manogue with Tevian Dray, Kenneth S. Krane, Jason Janesky Department of Physics Oregon State University Corvallis, OR 97331, USA corinne@physics.oregonstate.edu c 2007 Corinne
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 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 informationPhysics 5153 Classical Mechanics. Canonical Transformations-1
1 Introduction Physics 5153 Classical Mechanics Canonical Transformations The choice of generalized coordinates used to describe a physical system is completely arbitrary, but the Lagrangian is invariant
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 informationMechanics, Heat, Oscillations and Waves Prof. V. Balakrishnan Department of Physics Indian Institute of Technology, Madras
Mechanics, Heat, Oscillations and Waves Prof. V. Balakrishnan Department of Physics Indian Institute of Technology, Madras Lecture - 21 Central Potential and Central Force Ready now to take up the idea
More informationGravitation. chapter 9
chapter 9 Gravitation Circular orbits (Section 9.3) 1, 2, and 3 are simple exercises to supplement the quantitative calculations of Examples 4, 5, and 6 in Section 9.3. 1. Satellite near Earth s surface
More information1 Reduced Mass Coordinates
Coulomb Potential Radial Wavefunctions R. M. Suter April 4, 205 Reduced Mass Coordinates In classical mechanics (and quantum) problems involving several particles, it is convenient to separate the motion
More informationThe Calculus of Vec- tors
Physics 2460 Electricity and Magnetism I, Fall 2007, Lecture 3 1 The Calculus of Vec- Summary: tors 1. Calculus of Vectors: Limits and Derivatives 2. Parametric representation of Curves r(t) = [x(t), y(t),
More informationThe Principle of Least Action
The Principle of Least Action In their never-ending search for general principles, from which various laws of Physics could be derived, physicists, and most notably theoretical physicists, have often made
More informationPhysical Dynamics (SPA5304) Lecture Plan 2018
Physical Dynamics (SPA5304) Lecture Plan 2018 The numbers on the left margin are approximate lecture numbers. Items in gray are not covered this year 1 Advanced Review of Newtonian Mechanics 1.1 One Particle
More informationPhysical Dynamics (PHY-304)
Physical Dynamics (PHY-304) Gabriele Travaglini March 31, 2012 1 Review of Newtonian Mechanics 1.1 One particle Lectures 1-2. Frame, velocity, acceleration, number of degrees of freedom, generalised coordinates.
More informationOrbits, Integrals, and Chaos
Chapter 7 Orbits, Integrals, and Chaos In n space dimensions, some orbits can be formally decomposed into n independent periodic motions. These are the regular orbits; they may be represented as winding
More informationPHY 5246: Theoretical Dynamics, Fall September 28 th, 2015 Midterm Exam # 1
Name: SOLUTIONS PHY 5246: Theoretical Dynamics, Fall 2015 September 28 th, 2015 Mierm Exam # 1 Always remember to write full work for what you do. This will help your grade in case of incomplete or wrong
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 informationSession 6: Analytical Approximations for Low Thrust Maneuvers
Session 6: Analytical Approximations for Low Thrust Maneuvers As mentioned in the previous lecture, solving non-keplerian problems in general requires the use of perturbation methods and many are only
More informationis conserved, calculating E both at θ = 0 and θ = π/2 we find that this happens for a value ω = ω given by: 2g
UNIVERSITETET I STAVANGER Institutt for matematikk og naturvitenskap Suggested solutions, FYS 500 Classical Mechanics Theory 2016 fall Set 5 for 23. September 2016 Problem 27: A string can only support
More informationGravitation and Central-force motion
Chapter 7 Gravitation and Central-force motion In this chapter we describe motion caused by central forces, especially the orbits of planets, moons, and artificial satellites due to central gravitational
More informationOrbit Characteristics
Orbit Characteristics We have shown that the in the two body problem, the orbit of the satellite about the primary (or vice-versa) is a conic section, with the primary located at the focus of the conic
More informationThe 3 dimensional Schrödinger Equation
Chapter 6 The 3 dimensional Schrödinger Equation 6.1 Angular Momentum To study how angular momentum is represented in quantum mechanics we start by reviewing the classical vector of orbital angular momentum
More information= 0. = q i., q i = E
Summary of the Above Newton s second law: d 2 r dt 2 = Φ( r) Complicated vector arithmetic & coordinate system dependence Lagrangian Formalism: L q i d dt ( L q i ) = 0 n second-order differential equations
More informationConstraints. Noninertial coordinate systems
Chapter 8 Constraints. Noninertial codinate systems 8.1 Constraints Oftentimes we encounter problems with constraints. F example, f a ball rolling on a flo without slipping, there is a constraint linking
More information3 = arccos. A a and b are parallel, B a and b are perpendicular, C a and b are normalized, or D this is always true.
Math 210-101 Test #1 Sept. 16 th, 2016 Name: Answer Key Be sure to show your work! 1. (20 points) Vector Basics: Let v = 1, 2,, w = 1, 2, 2, and u = 2, 1, 1. (a) Find the area of a parallelogram spanned
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 informationCentral Force Motion Challenge Problems
Central Force Motion Challenge Problems Problem 1: Ellitic Orbit A satellite of mass m s is in an ellitical orbit around a lanet of mass m which is located at one focus of the ellise. The satellite has
More informationHW 6 Mathematics 503, Mathematical Modeling, CSUF, June 24, 2007
HW 6 Mathematics 503, Mathematical Modeling, CSUF, June 24, 2007 Nasser M. Abbasi June 15, 2014 Contents 1 Problem 1 (section 3.5,#9, page 197 1 2 Problem 1 (section 3.5,#9, page 197 7 1 Problem 1 (section
More informationMechanics, Heat, Oscillations and Waves Prof. V. Balakrishnan Department of Physics Indian Institute of Technology, Madras
Mechanics, Heat, Oscillations and Waves Prof. V. Balakrishnan Department of Physics Indian Institute of Technology, Madras Lecture 08 Vectors in a Plane, Scalars & Pseudoscalers Let us continue today with
More information/. 3D. a very abstract sort of problem, this research is a step on the way to applicable knowledge.
/. 3D Corwin Atwood-Stone Physics Department, The College of Wooster, Wooster, Ohio 691, USA (Dated: May 7, 009) A simulation of the mutual orbital motion of a point and a stick of equal masses in three
More information06. Lagrangian Mechanics II
University of Rhode Island DigitalCommons@URI Classical Dynamics Physics Course Materials 2015 06. Lagrangian Mechanics II Gerhard Müller University of Rhode Island, gmuller@uri.edu Creative Commons License
More informationMassachusetts Institute of Technology Department of Physics. Final Examination December 17, 2004
Massachusetts Institute of Technology Department of Physics Course: 8.09 Classical Mechanics Term: Fall 004 Final Examination December 17, 004 Instructions Do not start until you are told to do so. Solve
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 informationy(t) = y 0 t! 1 2 gt 2. With y(t final ) = 0, we can solve this for v 0 : v 0 A ĵ. With A! ĵ =!2 and A! = (2) 2 + (!
1. The angle between the vector! A = 3î! 2 ĵ! 5 ˆk and the positive y axis, in degrees, is closest to: A) 19 B) 71 C) 90 D) 109 E) 161 The dot product between the vector! A = 3î! 2 ĵ! 5 ˆk and the unit
More informationMATHEMATICAL MODELLING, MECHANICS AND MOD- ELLING MTHA4004Y
UNIVERSITY OF EAST ANGLIA School of Mathematics Main Series UG Examination 2017 18 MATHEMATICAL MODELLING, MECHANICS AND MOD- ELLING MTHA4004Y Time allowed: 2 Hours Attempt QUESTIONS 1 and 2, and ONE other
More informationDistance travelled time taken and if the particle is a distance s(t) along the x-axis, then its instantaneous speed is:
Chapter 1 Kinematics 1.1 Basic ideas r(t) is the position of a particle; r = r is the distance to the origin. If r = x i + y j + z k = (x, y, z), then r = r = x 2 + y 2 + z 2. v(t) is the velocity; v =
More information5 Home exercise sheet
5 Home exercise sheet 5.1 The central force problem and Scattering Exercise 7.1: Yukawa potential The Yukawa potential is 1. Write down the effective potential. U(r) = αe χr r (1) 2. Sketch the effective
More informationPHYS2330 Intermediate Mechanics Fall Final Exam Tuesday, 21 Dec 2010
Name: PHYS2330 Intermediate Mechanics Fall 2010 Final Exam Tuesday, 21 Dec 2010 This exam has two parts. Part I has 20 multiple choice questions, worth two points each. Part II consists of six relatively
More informationL03: Kepler problem & Hamiltonian dynamics
L03: Kepler problem & Hamiltonian dynamics 18.354 Ptolemy circa.85 (Egypt) -165 (Alexandria) Greek geocentric view of the universe Tycho Brahe 1546 (Denmark) - 1601 (Prague) "geo-heliocentric" system last
More informationRelativistic Dynamics
Chapter 4 Relativistic Dynamics The most important example of a relativistic particle moving in a potential is a charged particle, say an electron, moving in an electromagnetic field, which might be that
More informationRotational & Rigid-Body Mechanics. Lectures 3+4
Rotational & Rigid-Body Mechanics Lectures 3+4 Rotational Motion So far: point objects moving through a trajectory. Next: moving actual dimensional objects and rotating them. 2 Circular Motion - Definitions
More informationPhysics 142 Energy in Mechanics Page 1. Energy in Mechanics
Physics 4 Energy in Mechanics Page Energy in Mechanics This set of notes contains a brief review of the laws and theorems of Newtonian mechanics, and a longer section on energy, because of its central
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 informationStellar Dynamics and Structure of Galaxies
Stellar Dynamics and Structure of Galaxies in a given potential Vasily Belokurov vasily@ast.cam.ac.uk Institute of Astronomy Lent Term 2016 1 / 59 1 Collisions Model requirements 2 in spherical 3 4 Orbital
More information4 Orbits in stationary Potentials (BT 3 to page 107)
-313see http://www.strw.leidenuniv.nl/ franx/college/ mf-sts13-c-1-313see http://www.strw.leidenuniv.nl/ franx/college/ mf-sts13-c Orbits in stationary Potentials (BT 3 to page 17) Now we have seen how
More informationVISUAL PHYSICS ONLINE
VISUAL PHYSICS ONLINE EXCEL SIMULATION MOTION OF SATELLITES DOWNLOAD the MS EXCEL program PA50satellite.xlsx and view the worksheet Display as shown in the figure below. One of the most important questions
More informationTheorem of the Keplerian kinematics
myjournal manuscript No. (will be inserted by the editor) Theorem of the Keplerian kinematics Herve Le Cornec Received: date / Accepted: date Abstract As described in the literature the speed of a Keplerian
More informationPHYSICS 311: Classical Mechanics Final Exam Solution Key (2017)
PHYSICS 311: Classical Mechanics Final Exam Solution Key (017) 1. [5 points] Short Answers (5 points each) (a) In a sentence or two, explain why bicycle wheels are large, with all of the mass at the edge,
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 information