THE SOLAR SYSTEM. We begin with an inertial system and locate the planet and the sun with respect to it. Then. F m. Then

Similar documents
Inverse Kinematics 1 1/21/2018

Where Standard Physics Runs into Infinite Challenges, Atomism Predicts Exact Limits

Einstein's Energy Formula Must Be Revised

, an inverse square law.

DISCHARGE MEASUREMENT IN TRAPEZOIDAL LINED CANALS UTILIZING HORIZONTAL AND VERTICAL TRANSITIONS

Chapter 6 Control Systems Design by Root-Locus Method. Lag-Lead Compensation. Lag lead Compensation Techniques Based on the Root-Locus Approach.

@(; t) p(;,b t) +; t), (; t)) (( whih lat line follow from denition partial derivative. in relation quoted in leture. Th derive wave equation for ound

To determine the biasing conditions needed to obtain a specific gain each stage must be considered.

UVa Course on Physics of Particle Accelerators

Chapter 4. Simulations. 4.1 Introduction

Thermochemistry and Calorimetry

Practice Exam 2 Solutions

Laplace Transformation

The Laws of Electromagnetism Maxwell s Equations Displacement Current Electromagnetic Radiation

Energy-Work Connection Integration Scheme for Nonholonomic Hamiltonian Systems

ES 247 Fracture Mechanics Zhigang Suo. Applications of Fracture Mechanics

COMM 602: Digital Signal Processing. Lecture 8. Digital Filter Design

Moment of Inertia of an Equilateral Triangle with Pivot at one Vertex

Final Comprehensive Exam Physical Mechanics Friday December 15, Total 100 Points Time to complete the test: 120 minutes

Bogoliubov Transformation in Classical Mechanics

Wave Phenomena Physics 15c

Sound Propagation through Circular Ducts with Spiral Element Inside

q-expansions of vector-valued modular forms of negative weight

V V The circumflex (^) tells us this is a unit vector

Physics 2. Angular Momentum. Prepared by Vince Zaccone For Campus Learning Assistance Services at UCSB

5.2.6 COMPARISON OF QUALITY CONTROL AND VERIFICATION TESTS

Physics 741 Graduate Quantum Mechanics 1 Solutions to Final Exam, Fall 2014

DIFFERENTIAL EQUATIONS

DIFFERENTIAL EQUATIONS Laplace Transforms. Paul Dawkins

Pulsed Magnet Crimping

Period #8: Axial Load/Deformation in Indeterminate Members

Mechanics Physics 151

Correction for Simple System Example and Notes on Laplace Transforms / Deviation Variables ECHE 550 Fall 2002

Analysis of Feedback Control Systems

Physics 6A. Angular Momentum. Prepared by Vince Zaccone For Campus Learning Assistance Services at UCSB

TMA4125 Matematikk 4N Spring 2016

Chapter 26 Lecture Notes

Lecture 17: Analytic Functions and Integrals (See Chapter 14 in Boas)

Factor Analysis with Poisson Output

into a discrete time function. Recall that the table of Laplace/z-transforms is constructed by (i) selecting to get

Heat transfer and absorption of SO 2 of wet flue gas in a tube cooled

EE 333 Electricity and Magnetism, Fall 2009 Homework #11 solution

online learning Unit Workbook 4 RLC Transients

Effect of Anisotropy Asymmetry on the Switching. Behavior of Exchange Biased Bilayers

Econ 455 Answers - Problem Set 4. where. ch ch ch ch ch ch ( ) ( ) us us ch ch us ch. (world price). Combining the above two equations implies: 40P

Green s function for the wave equation

Deepak Rajput

arxiv:gr-qc/ v2 6 Feb 2004

A Queueing Model for Call Blending in Call Centers

ELECTROMAGNETIC NORMAL MODES AND DISPERSION FORCES.

231 Outline Solutions Tutorial Sheet 7, 8 and January Which of the following vector fields are conservative?

Halliday/Resnick/Walker 7e Chapter 6

Singular perturbation theory

The elementary spike produced by a pure e + e pair-electromagnetic pulse from a Black Hole: The PEM Pulse

Critical Percolation Probabilities for the Next-Nearest-Neighboring Site Problems on Sierpinski Carpets

Notes on Phase Space Fall 2007, Physics 233B, Hitoshi Murayama

MATEMATIK Datum: Tid: eftermiddag. A.Heintz Telefonvakt: Anders Martinsson Tel.:

The elementary spike produced by a pure e + e pair-electromagnetic pulse from a Black Hole: The PEM Pulse

The Prime Number Theorem

Tutorial 2 Euler Lagrange ( ) ( ) In one sentence: d dx

Chapter 4. The Laplace Transform Method

Given the following circuit with unknown initial capacitor voltage v(0): X(s) Immediately, we know that the transfer function H(s) is

PHYS 110B - HW #6 Spring 2004, Solutions by David Pace Any referenced equations are from Griffiths Problem statements are paraphrased

Dr G. I. Ogilvie Lent Term 2005

INDIVIDUAL OVERTOPPING EVENTS AT DIKES

Design Manual to EC2. LinkStudPSR. Version 3.1 BS EN : Specialists in Punching Shear Reinforcement.

On the Stationary Convection of Thermohaline Problems of Veronis and Stern Types

p. (The electron is a point particle with radius r = 0.)

Nonlinear Dynamics of Single Bunch Instability in Accelerators *

Astr 5465 Mar. 29, 2018 Galactic Dynamics I: Disks

Lecture 15 - Current. A Puzzle... Advanced Section: Image Charge for Spheres. Image Charge for a Grounded Spherical Shell

Bernoulli s equation may be developed as a special form of the momentum or energy equation.

Name Solutions to Test 1 September 23, 2016

PHYSICS 211 MIDTERM II 12 May 2004

Hylleraas wavefunction for He. dv 2. ,! r 2. )dv 1. in the trial function. A simple trial function that does include r 12. is ) f (r 12.

Strauss PDEs 2e: Section Exercise 3 Page 1 of 13. u tt c 2 u xx = cos x. ( 2 t c 2 2 x)u = cos x. v = ( t c x )u

ME 375 EXAM #1 Tuesday February 21, 2006

Solving Radical Equations

Green s function for the wave equation

Circular Motion Problem Solving

Reading assignment: In this chapter we will cover Sections Definition and the Laplace transform of simple functions

Practice Problems - Week #7 Laplace - Step Functions, DE Solutions Solutions

Math Skills. Scientific Notation. Uncertainty in Measurements. Appendix A5 SKILLS HANDBOOK

The machines in the exercise work as follows:

Nonlinear Single-Particle Dynamics in High Energy Accelerators

13.Prandtl-Meyer Expansion Flow

Advanced Computational Fluid Dynamics AA215A Lecture 4

Universal Gravitation

V = 4 3 πr3. d dt V = d ( 4 dv dt. = 4 3 π d dt r3 dv π 3r2 dv. dt = 4πr 2 dr

Reading assignment: In this chapter we will cover Sections Definition and the Laplace transform of simple functions

The gravitational phenomena without the curved spacetime

Constant Force: Projectile Motion

Impulse. calculate the impulse given to an object calculate the change in momentum as the result of an impulse

Relativity in Classical Physics

Solving Differential Equations by the Laplace Transform and by Numerical Methods

Control Systems Analysis and Design by the Root-Locus Method

The Hand of God, Building the Universe and Multiverse

Sample Problems. Lecture Notes Related Rates page 1

Linear Momentum. calculate the momentum of an object solve problems involving the conservation of momentum. Labs, Activities & Demonstrations:

V V V V. Vectors. Mike Bailey. Vectors have Direction and Magnitude. Magnitude: x y z. Computer Graphics.

Transcription:

THE SOLAR SYSTEM We now want to apply what we have learned to the olar ytem. Hitorially thi wa the great teting ground for mehani and provided ome of it greatet triumph, uh a the diovery of the outer planet. We will ue a implified model of the ytem baed on the fat that the ma of the un i muh greater than that of any of the planet or in fat of all of them ombined. We will alo neglet the effet of the planet on eah other, although it i preiely thi whih led to the diovery of the outer planet. Hene we will uppoe that the only objet in the univere are the un and the planet other objet we are onidering. The firt tep i to hooe a atifatory oordinate ytem. It mut be inertial if we want to avoid the ue of the fititiou fore we have diued earlier. Ideally we would hooe a oordinate ytem either entered on the earth (Pre Coperniu) or on the un (Pot Coperniu). At firt glane neither would be atifatory beaue of Newton third law. Eah exert a fore on the other and thu neither i unaelerated. However thi problem an be overome, a we will now ee. We begin with an inertial ytem and loate the planet and the un with repet to it. Then Rp R r Ap A a F F m A m A a p p p p p F F Fp ma A m p Then Fp Fp mp a m

mp Fp 1 mpa m and mm p Fp a a m m p where μ (the redued ma) = m m p /(m +m p). Thu we an hooe the enter of the un a the origin of our oordinate ytem if we merely replae the ma of the planet by it redued ma. In the ae of the olar ytem all mae are mall ompared to that of the un, and hene the orretion i uually negligible. However if we were to onider the moon going around the earth thi would not be the ae. The phyial reaon for thi effet i eay to undertand. An oberver on either objet ee both her aeleration and the aeleration of the other objet toward her, a diued in la. We now have a very imple fore diagram for planet. There i only one fore the gravitational attration of the un. Sine it i direted through the axi it produe no torque: RFp 0 Then beaue of Newton eond law written in term of angular quantitie: We mut have the angular momentum ontant. But ine angular momentum i a vetor it diretion mut be ontant. Sine: dr dt J rp rmv and both r and v are in the plane of the orbit we ee that J i perpendiular to the plane of the orbit. Then ine that diretion i fixed we get our firt onluion the planet mut move in a plane. We will find the orbit exatly in a moment, but firt we onider a peial ae in whih the orbit i irular. It turn out that mot of the planet move in nearly irular orbit o thi i a good tarting point. In thi ae the problem i very imple. The aeleration toward the enter of the irle mut be provided by the gravitational attration:

mp vp FG m p R R FM Rv p But vp R P where P i the period of the planet the time to go one around the un. Then or 3 4 R P 3 R p 4 (the ame for all objet in the olar ytem). Thi reult i Kepler third law, and provide a mean of meauring the ma of the un. We now turn to the general problem where the orbit are not irular. Beaue the orbit may repeat it will be muh more onvenient to ue polar rather than Carteian oordinate (note that we are in a plane o pherial oordinate are not needed). We now need to alulate the aeleration in polar oordinate. We do thi a before. To get ˆ dr / dt onider the keth r rrˆ dr dr drˆ v rˆ r dt dt dt

Then drˆ d ˆ dt dt Hene dr d v rˆ r ˆ dt dt Of oure thi i obviou. dr/dt i jut the omponent of veloity along the radiu d r r dt i jut the tangential veloity along the irumferene. Next dv d r dr drˆ dr d d d dˆ a rˆ ˆ r ˆ r dt dt dt dt dt dt dt dt dt Hene we need d ˆ /dt. We find thi from the keth Thu dˆ d ˆr dt dt

Then d r d ˆ d dr d a rˆ r r dt dt dt dt dt Again mot of thee term are obviou. d r/dt i jut the aeleration along r. i jut the familiar entripetal aeleration. d v v r r r dt r r d d r r r dt dt i jut the tangential aeleration along the urve. The lat term dr d dt dt i new beaue we have not previouly onidered aeleration in both the radial and tangential diretion at the ame time. Sine we know the fore we an now write Newton eond law in our hoen oordinate ytem. m F r p rˆ d r dr dt dt r (1) d dr d r 0 dt dt dt () Thi look a bit intimidating until we ue our knowledge of the phyi of the ituation. We know that angular momentum i ontant. We have already ued the fat that the diretion i ontant to ee that the motion i in a plane. We now ue the fat that the magnitude i ontant to olve ().

Hene dr d ˆ d J mr v mrrˆ rˆ r mr zˆ dt dt dt d r ontant dt Thu d d dr d d d dr d r 0 r r r r dt dt dt dt dt dt dt dt Thu d r ont A dt i the olution of (). We now ubtitute thi reult in (1) to get: A r dr dt r r (3) To olve (3) it i ueful to make a hange of variable. Let: 1 r u Then dr dt 1 du u dt dr du 1 du 3 dt dt u u dt and du 1 d u 3 Au dt u u dt 3 u (4)

Sine we are intereted in the hape of the orbit we really want u(θ) rather than u(t). We therefore make another hange of variable. We have: We now ue the hain rule to get: Subtituting thi in (4) give: or u(t) = u(θ(t)) du du d du A du Au dt d dt d r d d u d du d u du Au Au Au Au Au d d dt d d du 1 d u 1du 3 Au Au Au A u 3 u u d u d ud du 3 Au Au u d d u u d A But thi i an equation we know how to olve. We do it a the um of the general olution of the homogeneou equation plu one olution of the whole equation. The olution of the homogeneou equation i obviouly: uh Bo where B and φ are arbitrary ontant. A olution of the whole equation i: Thu u p A u uk up Bo A

To fix the arbitrary ontant we hooe θ to be zero at the point of loet approah. At that point we alo know the angular momentum i: But Hene Chooe φ = 0. Then rmv J A rv m 1 1 1 r r u Bo rv Bo rv r Br o 1 rv 1 r Br 1 rv r rv 1 r r rv B 1 r r r rv 1 o 1 o rv rv 1 where 1 rv rv 1 rv

But thi i jut the equation of a oni etion with eentriity ε. There are now everal ae to onider depending on the value of ε. We begin by aking whether or not the planet will return, i.e. i it in a loed orbit or doe it ultimately leave the olar ytem. To anwer thi we need only look at the energy. Thi i given at the point of loet approah by: m 1 1 E mv m v r r If the planet i to reah infinite ditane it will have to have a veloity greater than or equal to zero at infinite ditane. Hene E mut be greater than or equal to zero. 1 1 r v v 0 v 1 r r If thi ondition i not met the planet will return and alo we an predit how far from the un it will get. To do o we note that both the angular momentum and the energy are onerved. Thu: or rv r v m m 1 1 v vm r r 1 r v p m m E m r r E 1 rm rm rv 0 m p r m E r v E m p mp 1/ Sine E/m p < 0 we an ue either ign. One give r, the other r m (the maximum ditane).

The ae to be onidered are then: Cae 1: < 1 Cae : ε = 1 Cae 3: ε = 0 Cae 4: ε > 1 CASE 3 Thi i the implet ituation. We then have: r v r r v rv or r 4 P Thi i jut the ae of irular orbit diued above. CASE Thi i a bit more ompliated, but not muh. In thi ae we have: rv rv 11 1 Thi i the ituation in whih the planet would jut reah infinite ditane with zero veloity. We an readily find the equation of the orbit a follow. Swithing to Carteian oordinate we have rv r 1 o

rv 1/ x y x 1 1/ x y 1/ r v rv rv rv x y x x y x xx But thi i jut the equation of a parabola. y r r x y x r r CASE 1 Thi i the ae for whih the objet will return, and hene the one that applie for planet in the olar ytem. Thi i more ompliated to analyze, o ome patiene i required. Now we have: rv rv r r 1/ r x y 1o x 1 1/ x y 1/ rv x y x r rv rv rv x y r x r r x x rv rv x 1 r xy r

We now omplete the quare on the x term x rv r y rv r x 1 1 1 rv rv rv y r r r x 1 1 1 1 rv r rv 1 r 1 1 1 Sine < 1, thi an equation of the form x y 1 a b whih i the equation of an ellipe entered at x = -α with emi-major axi a and emi-minor axi b. Note that But rv rv rv rv rv rv r r 1 r 1 r 1 1 rv rv rv rv 1 1

r o r rm E 1 v m p r i the enter of the ellipe. Then rv rv o 1 1 v rv v r r r 1 r r r Thu we have the equation of an ellipe with origin at the enter of the ellipe. CASE 4 Now all that hange i the ign of (1 ε ). Thi hange the equation to: x y 1 a b whih i the equation of a hyperbola. We now turn to the quetion of the motion a a funtion of time. To do thi we note that: d A r d dt dt r A 3 r A d t d A 1 o o o 3 1 tan A in 1 1 tan 3/ 1/ o 1 1o 1 1 For example we an find the period of the orbit of a planet a:

P 3 1 tan 3 A 1 A tan 3/ 1/ 3/ 1 1 1 0 But 1/ r r v A A a 1 1 a 1 3 3/ P A 1 a A a 3/ 1/ 4 3 P a Thi i Kepler third law. In other word the planet have the ame period a irular orbit of radiu equal to the emi-major axi would have. It i now a trivial problem to find the period of any planet if we know it ditane and veloity relative to the un at loet approah.