TP 3:Runge-Kutta Methods-Solar System-The Method of Least Squares
|
|
- Mark Allison
- 5 years ago
- Views:
Transcription
1 TP :Runge-Kutta Methods-Solar System-The Method of Least Squares December 8, Runge-Kutta Method The problem is still trying to solve the first order differential equation dy = f(y, x). (1) dx In the Euler s method we approximate the function y = y(x) in each interval [x n, x n+1 ] by the straight line y n+1 = y n + xf(x n, y n ). (2) The slope f(x n, y n ) of this line is exactly given by the slope of the function y = y(x) at the begining of the inetrval [x n, x n+1 ]. Given the value y n at x n we evaluate the value y n+1 at x n+1 using the method of Runge-Kutta as follows. Firstly the middle of the interval [x n, x n+1 ] which is at the value x n x corresponds to the y value y n+1 calculated using the Euler s method, viz y n+1 = y n k 1 where Secondly the slope at this middle point (x n x, y n k 1) given by k 1 = xf(x n, y n ). () k 2 x = f(x n x, y n k 1) (4) is the value of the slope which will be used to estimate the correct value of y n+1 at x n+1 using again Euler s method,namely y n+1 = y n + k 2. (5) The error in this method is proportional to x. Let us apply this method to the problem of the harmonic oscillator. We have the differential equations Euler s equations read dθ dt = ω dω dt = g θ. (6) l θ n+1 = θ n + tω n ω n+1 = ω n g l θ n t. (7) 1
2 First we consider the function θ = θ(t). The middle point is (t n t, θ n k 1) where k 1 = tω n. For the function ω = ω(t) the middle point is (t n t, ω n k ) where k = g l tθ n. Therefore we have The slope of the function θ(t) at its middle point is The slope of the function ω(t) at its middle point is k 1 = tω n k = g l tθ n. (8) k 2 t = ω n k. (9) The Runge-Kutta solutions are then given by k 4 t = g l (θ n k 1). (10) θ n+1 = θ n + k 2 ω n+1 = ω n + k 4. (11) 2 Kepler s Laws-The Inverse-Square Law Newton s second law: We consider the motion of the earth around the sun. Let r be the distance and M s and M e the masses of the sun and the earth respectively. We neglect the effect of the other planets and the motion of the sun (i.e we assume that M s >> M e ). The goal is to calculate the position of the earth as a function of time. We start from Newton s second law of motion We get the two equations d 2 r M e dt 2 = GM em s r r = GM em s r (x i + y j). (12) d 2 x dt 2 = r x. (1) d 2 y dt 2 = r y. (14) We replace these two second order differential equations by the four first order differential equations dx dt = v x. (15) dv x r x. (16) dy dt = v y. (17) dv y r y. (18) 2
3 Astronomical units: The distance will be measured in astronomical units (AU) whereas time will be measured in years. One astronomical unit of lenght (1 AU) is equal to the average distance between the earth and the sun, viz 1AU = m. The astronomical unit of mass can be found s follows. Assuming a circular orbit we have Equivalently M e v 2 r = M e r 2. (19) = v 2 r. (20) The radius is r = 1AU. The velocity of the earth is v = 2πr/yr = 2πAU/yr. Hence Kepler s laws: These are given by the following three statements: = 4π 2 AU /yr 2. (21) The planets move in elliptical orbits around the sun. The sun resides at one focus. The line joining the sun with any planet sweeps out equal areas in equal times. Given an orbit with a period T and a semimajor axis a the ratio T 2 /a is a constant. The inverse-square law and stability of orbits: Any object with mass generates a gravitational field. In other words gravitational field lines will emanate from the object and radiate outward to infinity. The number of field lines N is proportional to the mass. The density of field lines crossing a sphere of radius r surrounding this object is given by N/4πr 2. This is the origin of the inverse-square law. Therefore any other object placed in this gravitational field will experience a gravitational force proportional to the number of field lines which intersect it. If the distance between this second object and the source is increased the force on it will become weaker because the number of field lines which intersect it will decrease as we are further away from the source. Euler-Cromer Algorithm The time discretization is t t(i) = (i 1) t, i = 1,..., L + 1. (22) The total time interval T = L t. We define x(t) = x(i 1), v x (t) = v x (i 1), y(t) = y(i 1), v y (t) = v y (i 1). Equations (15), (16), (17) and (18) become (with i = 2,,..., L + 1) v x (i) = v x (i 1) x(i 1) t. (2) x(i) = x(i 1) + v x (i 1) t. (24) v y (i) = v y (i 1) y(i 1) t. (25) y(i) = y(i 1) + v y (i 1) t. (26)
4 This is Euler s algorithm. In order to maintain energy conservation we employ Euler-Cromer algorithm. We calculate as in the Euler s algorithm the velocity at time step i by using the position and velocity at time step i 1. However we compute the position at time step i by using the position at time step i 1 and the velocity at time stpe i, viz v x (i) = v x (i 1) x(i 1) t. (27) x(i) = x(i 1) + v x (i) t. (28) v y (i) = v y (i 1) y(i 1) t. (29) y(i) = y(i 1) + v y (i) t. (0) 4 Runge-Kutta Algorithm Let us consider the equations dx dt = v x. (1) dv x r x. (2) dy dt = v y. () dv y r y. (4) First we consider the function x = x(t). The middle point is (t n t, x n k 1) where k 1 = t v xn. For the function v x = v x (t) the middle point is (t n t, v xn k ) where k = GMs t x n. Therefore we have The slope of the function x(t) at the middle point is The slope of the function v x (t) at the middle point is k 4 k 1 = t v xn k = t x n. (5) k 2 t = v xn k. (6) t = (x n k 1). (7) R n Next we consider the function y = y(t). The middle point is (t n t, y n k 1 ) where k 1 = t v yn. For the function v y = v y (t) the middle point is (t n t, v yn k ) where k = GMs t y n. Therefore we have k 1 = t v yn k = t y n. (8) 4
5 The slope of the function y(t) at the middle point is The slope of the function v y (t) at the middle point is k 2 t = v yn k. (9) In the above equations R n = k 4 t = R n (y n k 1 ). (40) (x n k 1) 2 + (y n k 1 )2. (41) The Runge-Kutta solutions are then given by x n+1 = x n + k 2 v x(n+1) = v xn + k 4 y n+1 = y n + k 2 v y(n+1) = v yn + k 4. (42) 5 Solar System-Precession of the Perihelion of Mercury The orbit of mercury is elliptic. The orientation of the axes of the ellipse rotate with time. This is the precession of the perihelion (the point of the orbit nearest to the sun) of mercury. Mercury s perihelion makes one revolution every 2000 years. This is approximately 566 arcseconds per century. The gravitational forces of the other planets (in particular Jupiter) lead to a precession of 52 arcseconds per century. The remaining 4 arcseconds per century are accounted for by general relativity. For objects too close together (like the sun and mercury) the force of gravity predicted by general relativity deviates from the inverse-square law. This force is given by F = M m r 2 (1 + α r 2 ), α = AU 2. (4) Initial Conditions: For the numerical simulations it is important to determine the correct initial conditions. The orbit of mercury is known to be an ellipse with eccentricity e = and radius (semimajor axis) a = 0.9 AU with the sun at one of the foci. The distance between the sun and the center is ea. The first initial condition is x 0 = r 1, y 0 = 0 where r 1 is the maximum distance from mercury to the sun,i.e r 1 = (1 + e)a = 0.47 AU. The second initial condition is the velocity (0, v 1 ) which can be computed using conservation of energy and angular momentum. For example by comparing with the point (0, b) on the orbit where b is the semiminor axis, i.e b = a 1 e 2 the velocity (v 2, 0) there can be obtained in terms of (0, v 1 ) from conservation of angular momentum as follows Next conservation of energy yields M m r 1 r 1 v 1 = bv 2 v 2 = r 1v 1 b. (44) M mv 2 1 = M m r M mv 2 2. (45) In above r 2 = e 2 a 2 + b 2 is the distance between the sun and mercury when at the point (0, b). By substituting the value of v 2 we get an equation for v 1. This is given by GMs 1 e v 1 = = 8.2 AU/yr. (46) a 1 + e 5
6 6 Method of Least Squares We assume that we have N data points (x(i), y(i)). We want to fit this data to some curve say a straight line y fit = mx + b. To this end we define the function N N = (y(i) y fit (i)) 2 = (y(i) mx(i) b) 2. (47) i=1 The goal is to minimize this function with respect to b and m. We have i=1 We get the solution m = 0, = 0. (48) b i b = x(i) j x(j)y(j) i x(i)2 j y(j) ( i x(i))2 N. (49) i x2 i i m = x(i) j y(j) N i x(i)y(i) ( i x(i))2 N. (50) i x2 i 6
Lecture 16. Gravitation
Lecture 16 Gravitation Today s Topics: The Gravitational Force Satellites in Circular Orbits Apparent Weightlessness lliptical Orbits and angular momentum Kepler s Laws of Orbital Motion Gravitational
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 informationVISUAL PHYSICS ONLINE
VISUAL PHYSICS ONLINE PRACTICAL ACTIVITY HOW DO THE PANETS MOVE? One of the most important questions historically in Physics was how the planets move. Many historians consider the field of Physics to date
More informationA SIMULATION OF THE MOTION OF AN EARTH BOUND SATELLITE
DOING PHYSICS WITH MATLAB A SIMULATION OF THE MOTION OF AN EARTH BOUND SATELLITE Download Directory: Matlab mscripts mec_satellite_gui.m The [2D] motion of a satellite around the Earth is computed from
More informationKEPLER S LAWS OF PLANETARY MOTION
KEPLER S LAWS OF PLANETARY MOTION In the early 1600s, Johannes Kepler culminated his analysis of the extensive data taken by Tycho Brahe and published his three laws of planetary motion, which we know
More information5.1. Accelerated Coordinate Systems:
5.1. Accelerated Coordinate Systems: Recall: Uniformly moving reference frames (e.g. those considered at 'rest' or moving with constant velocity in a straight line) are called inertial reference frames.
More informationObservational Astronomy - Lecture 4 Orbits, Motions, Kepler s and Newton s Laws
Observational Astronomy - Lecture 4 Orbits, Motions, Kepler s and Newton s Laws Craig Lage New York University - Department of Physics craig.lage@nyu.edu February 24, 2014 1 / 21 Tycho Brahe s Equatorial
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 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 informationOctober 19, NOTES Solar System Data Table.notebook. Which page in the ESRT???? million km million. average.
Celestial Object: Naturally occurring object that exists in space. NOT spacecraft or man-made satellites Which page in the ESRT???? Mean = average Units = million km How can we find this using the Solar
More informationNewton s Gravitational Law
1 Newton s Gravitational Law Gravity exists because bodies have masses. Newton s Gravitational Law states that the force of attraction between two point masses is directly proportional to the product of
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 information18. Kepler as a young man became the assistant to A) Nicolaus Copernicus. B) Ptolemy. C) Tycho Brahe. D) Sir Isaac Newton.
Name: Date: 1. The word planet is derived from a Greek term meaning A) bright nighttime object. B) astrological sign. C) wanderer. D) nontwinkling star. 2. The planets that were known before the telescope
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 informationLecture 1a: Satellite Orbits
Lecture 1a: Satellite Orbits Meteorological Satellite Orbits LEO view GEO view Two main orbits of Met Satellites: 1) Geostationary Orbit (GEO) 1) Low Earth Orbit (LEO) or polar orbits Orbits of meteorological
More informationDynamical properties of the Solar System. Second Kepler s Law. Dynamics of planetary orbits. ν: true anomaly
First Kepler s Law The secondary body moves in an elliptical orbit, with the primary body at the focus Valid for bound orbits with E < 0 The conservation of the total energy E yields a constant semi-major
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 informationA = 6561 times greater. B. 81 times greater. C. equally strong. D. 1/81 as great. E. (1/81) 2 = 1/6561 as great Pearson Education, Inc.
Q13.1 The mass of the Moon is 1/81 of the mass of the Earth. Compared to the gravitational force that the Earth exerts on the Moon, the gravitational force that the Moon exerts on the Earth is A. 81 2
More informationAy 1 Lecture 2. Starting the Exploration
Ay 1 Lecture 2 Starting the Exploration 2.1 Distances and Scales Some Commonly Used Units Distance: Astronomical unit: the distance from the Earth to the Sun, 1 au = 1.496 10 13 cm ~ 1.5 10 13 cm Light
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 informationChapter 13. Universal Gravitation
Chapter 13 Universal Gravitation Planetary Motion A large amount of data had been collected by 1687. There was no clear understanding of the forces related to these motions. Isaac Newton provided the answer.
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 informationPhysics Lecture 03: FRI 29 AUG
Physics 23 Jonathan Dowling Isaac Newton (642 727) Physics 23 Lecture 03: FRI 29 AUG CH3: Gravitation III Version: 8/28/4 Michael Faraday (79 867) 3.7: Planets and Satellites: Kepler s st Law. THE LAW
More informationName Period Date Earth and Space Science. Solar System Review
Name Period Date Earth and Space Science Solar System Review 1. is the spinning a planetary object on its axis. 2. is the backward motion of planets. 3. The is a unit less number between 0 and 1 that describes
More informationPhys 7221, Fall 2006: Homework # 5
Phys 722, Fal006: Homewor # 5 Gabriela González October, 2006 Prob 3-: Collapse of an orbital system Consier two particles falling into each other ue to gravitational forces, starting from rest at a istance
More informationSatellite meteorology
GPHS 422 Satellite meteorology GPHS 422 Satellite meteorology Lecture 1 6 July 2012 Course outline 2012 2 Course outline 2012 - continued 10:00 to 12:00 3 Course outline 2012 - continued 4 Some reading
More informationThe Law of Ellipses (Kepler s First Law): all planets orbit the sun in a
Team Number Team Members Present Learning Objectives 1. Practice the Engineering Process a series of steps to follow to design a solution to a problem. 2. Practice the Five Dimensions of Being a Good Team
More informationSatellite Communications
Satellite Communications Lecture (3) Chapter 2.1 1 Gravitational Force Newton s 2nd Law: r r F = m a Newton s Law Of Universal Gravitation (assuming point masses or spheres): Putting these together: r
More informationChapter 4 Thrills and Chills +Math +Depth Acceleration of the Moon +Concepts The Moon is 60 times further away from the center of Earth than objects on the surface of Earth, and moves about Earth in an
More informationASTRO 1050 LAB #3: Planetary Orbits and Kepler s Laws
ASTRO 1050 LAB #3: Planetary Orbits and Kepler s Laws ABSTRACT Johannes Kepler (1571-1630), a German mathematician and astronomer, was a man on a quest to discover order and harmony in the solar system.
More informationPhysics 12. Unit 5 Circular Motion and Gravitation Part 2
Physics 12 Unit 5 Circular Motion and Gravitation Part 2 1. Newton s law of gravitation We have seen in Physics 11 that the force acting on an object due to gravity is given by a well known formula: F
More informationLecture 9(+10) Physics 106 Spring 2006
3/22/2006 Andrei Sirenko, NJIT 3 3/22/2006 Andrei Sirenko, NJIT 4 Lecture 9(+10) Physics 106 Spring 2006 Gravitation HW&R http://web.njit.edu/~sirenko/ Gravitation On Earth: the Earth gravitation dominates
More informationComputational Physics (6810): Session 8
Computational Physics (6810): Session 8 Dick Furnstahl Nuclear Theory Group OSU Physics Department February 24, 2014 Differential equation solving Session 7 Preview Session 8 Stuff Solving differential
More informationPHYS 155 Introductory Astronomy
PHYS 155 Introductory Astronomy - observing sessions: Sunday Thursday, 9pm, weather permitting http://www.phys.uconn.edu/observatory - Exam - Tuesday March 20, - Review Monday 6:30-9pm, PB 38 Marek Krasnansky
More informationThe Heliocentric Model of Copernicus
Celestial Mechanics The Heliocentric Model of Copernicus Sun at the center and planets (including Earth) orbiting along circles. inferior planets - planets closer to Sun than Earth - Mercury, Venus superior
More informationGravitation. Kepler s Law. BSc I SEM II (UNIT I)
Gravitation Kepler s Law BSc I SEM II (UNIT I) P a g e 2 Contents 1) Newton s Law of Gravitation 3 Vector representation of Newton s Law of Gravitation 3 Characteristics of Newton s Law of Gravitation
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 informationLesson 9. Luis Anchordoqui. Physics 168. Tuesday, October 24, 17
Lesson 9 Physics 168 1 Static Equilibrium 2 Conditions for Equilibrium An object with forces acting on it but that is not moving is said to be in equilibrium 3 Conditions for Equilibrium (cont d) First
More informationGravity and the Orbits of Planets
Gravity and the Orbits of Planets 1. Gravity Galileo Newton Earth s Gravity Mass v. Weight Einstein and General Relativity Round and irregular shaped objects 2. Orbits and Kepler s Laws ESO Galileo, Gravity,
More informationName. Satellite Motion Lab
Name Satellite Motion Lab Purpose To experiment with satellite motion using an interactive simulation in order to gain an understanding of Kepler s Laws of Planetary Motion and Newton s Law of Universal
More informationGeneral Relativity ASTR 2110 Sarazin. Einstein s Equation
General Relativity ASTR 2110 Sarazin Einstein s Equation Curvature of Spacetime 1. Principle of Equvalence: gravity acceleration locally 2. Acceleration curved path in spacetime In gravitational field,
More information(b) The period T and the angular frequency ω of uniform rotation are related to the cyclic frequency f as. , ω = 2πf =
PHY 302 K. Solutions for problem set #9. Non-textbook problem #1: (a) Rotation frequency of 1 Hz means one revolution per second, or 60 revolutions per minute (RPM). The pre-lp vinyl disks rotated at 78
More informationEarth Science Unit 6: Astronomy Period: Date: Elliptical Orbits
Earth Science Name: Unit 6: Astronomy Period: Date: Lab # 5 Elliptical Orbits Objective: To compare the shape of the earth s orbit (eccentricity) with the orbits of and with a circle. other planets Focus
More informationRadial Acceleration. recall, the direction of the instantaneous velocity vector is tangential to the trajectory
Radial Acceleration recall, the direction of the instantaneous velocity vector is tangential to the trajectory 1 Radial Acceleration recall, the direction of the instantaneous velocity vector is tangential
More informationLecture XIX: Particle motion exterior to a spherical star
Lecture XIX: Particle motion exterior to a spherical star Christopher M. Hirata Caltech M/C 350-7, Pasadena CA 95, USA Dated: January 8, 0 I. OVERVIEW Our next objective is to consider the motion of test
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 informationNovember 16, Henok Tadesse, Electrical Engineer, B.Sc. Ethiopia. or
The outward acceleration of galaxies may be a result of a non uniform and non linear distribution of matter in the universe Non local gravity directed upwards due to higher density outwards! Non elliptical
More informationNotes 10-3: Ellipses
Notes 10-3: Ellipses I. Ellipse- Definition and Vocab An ellipse is the set of points P(x, y) in a plane such that the sum of the distances from any point P on the ellipse to two fixed points F 1 and F
More informationPhys 2101 Gabriela González
Phys 2101 Gabriela González Newton s law : F = Gm 1 m 2 /r 2 Explains why apples fall, why the planets move around the Sun, sciencebulletins.amnh.org And in YouTube! Explains just as well as Newtons why
More informationAstronomy 1143 Quiz 1 Review
Astronomy 1143 Quiz 1 Review Prof. Pradhan September 7, 2017 I What is Science? 1. Explain the difference between astronomy and astrology. Astrology: nonscience using zodiac sign to predict the future/personality
More informationPHYSICS. Chapter 13 Lecture FOR SCIENTISTS AND ENGINEERS A STRATEGIC APPROACH 4/E RANDALL D. KNIGHT Pearson Education, Inc.
PHYSICS FOR SCIENTISTS AND ENGINEERS A STRATEGIC APPROACH 4/E Chapter 13 Lecture RANDALL D. KNIGHT Chapter 13 Newton s Theory of Gravity IN THIS CHAPTER, you will learn to understand the motion of satellites
More informationBlack Holes. Jan Gutowski. King s College London
Black Holes Jan Gutowski King s College London A Very Brief History John Michell and Pierre Simon de Laplace calculated (1784, 1796) that light emitted radially from a sphere of radius R and mass M would
More informationORBITS WRITTEN Q.E. (June 2012) Each of the five problems is valued at 20 points. (Total for exam: 100 points)
ORBITS WRITTEN Q.E. (June 2012) Each of the five problems is valued at 20 points. (Total for exam: 100 points) PROBLEM 1 A) Summarize the content of the three Kepler s Laws. B) Derive any two of the Kepler
More informationPHYS 101 Previous Exam Problems. Gravitation
PHYS 101 Previous Exam Problems CHAPTER 13 Gravitation Newton s law of gravitation Shell theorem Variation of g Potential energy & work Escape speed Conservation of energy Kepler s laws - planets Orbits
More informationPHYS 106 Fall 2151 Homework 3 Due: Thursday, 8 Oct 2015
PHYS 106 Fall 2151 Homework 3 Due: Thursday, 8 Oct 2015 When you do a calculation, show all your steps. Do not just give an answer. You may work with others, but the work you submit should be your own.
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 informationLecture 13. Gravity in the Solar System
Lecture 13 Gravity in the Solar System Guiding Questions 1. How was the heliocentric model established? What are monumental steps in the history of the heliocentric model? 2. How do Kepler s three laws
More informationJohannes Kepler ( ) German Mathematician and Astronomer Passionately convinced of the rightness of the Copernican view. Set out to prove it!
Johannes Kepler (1571-1630) German Mathematician and Astronomer Passionately convinced of the rightness of the Copernican view. Set out to prove it! Kepler s Life Work Kepler sought a unifying principle
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 informationANNEX 1. DEFINITION OF ORBITAL PARAMETERS AND IMPORTANT CONCEPTS OF CELESTIAL MECHANICS
ANNEX 1. DEFINITION OF ORBITAL PARAMETERS AND IMPORTANT CONCEPTS OF CELESTIAL MECHANICS A1.1. Kepler s laws Johannes Kepler (1571-1630) discovered the laws of orbital motion, now called Kepler's laws.
More informationChapter 8 - Gravity Tuesday, March 24 th
Chapter 8 - Gravity Tuesday, March 24 th Newton s law of gravitation Gravitational potential energy Escape velocity Kepler s laws Demonstration, iclicker and example problems We are jumping backwards to
More informationDays of the week: - named after 7 Power (moving) objects in the sky (Sun, Moon, 5 planets) Models of the Universe:
Motions of the Planets ( Wanderers ) Planets move on celestial sphere - change RA, Dec each night - five are visible to naked eye Mercury, Venus, Mars, Jupiter, Saturn Days of the week: - named after 7
More informationToday. Planetary Motion. Tycho Brahe s Observations. Kepler s Laws Laws of Motion. Laws of Motion
Today Planetary Motion Tycho Brahe s Observations Kepler s Laws Laws of Motion Laws of Motion In 1633 the Catholic Church ordered Galileo to recant his claim that Earth orbits the Sun. His book on the
More informationTycho Brahe ( )
Tycho Brahe (1546-1601) Foremost astronomer after the death of Copernicus. King Frederick II of Denmark set him up at Uraniborg, an observatory on the island of Hveen. With new instruments (quadrant),
More informationNotes on Planetary Motion
(1) Te motion is planar Notes on Planetary Motion Use 3-dimensional coordinates wit te sun at te origin. Since F = ma and te gravitational pull is in towards te sun, te acceleration A is parallel to te
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 informationChapter 13 Gravity Pearson Education, Inc. Slide 13-1
Chapter 13 Gravity Slide 13-1 The plan Lab this week: there will be time for exam problems Final exam: sections posted today; some left out Final format: all multiple choice, almost all short problems,
More informationKepler, Newton, and laws of motion
Kepler, Newton, and laws of motion First: A Little History Geocentric vs. heliocentric model for solar system (sec. 2.2-2.4)! The only history in this course is this progression: Aristotle (~350 BC) Ptolemy
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 informationPhysics Mechanics. Lecture 29 Gravitation
1 Physics 170 - Mechanics Lecture 29 Gravitation Newton, following an idea suggested by Robert Hooke, hypothesized that the force of gravity acting on the planets is inversely proportional to their distances
More informationUniform Circular Motion
Circular Motion Uniform Circular Motion Uniform Circular Motion Traveling with a constant speed in a circular path Even though the speed is constant, the acceleration is non-zero The acceleration responsible
More informationChapter 16 The Solar System
Chapter 16 The Solar System Finding the Standard Time and Date at Another Location Example When it is 12 noon in London, what is the standard time in Denver, Colorado (40 N, 105 W)? Section 15.3 Finding
More informationKepler's Laws and Newton's Laws
Kepler's Laws and Newton's Laws Kepler's Laws Johannes Kepler (1571-1630) developed a quantitative description of the motions of the planets in the solar system. The description that he produced is expressed
More informationAPS 1030 Astronomy Lab 79 Kepler's Laws KEPLER'S LAWS
APS 1030 Astronomy Lab 79 Kepler's Laws KEPLER'S LAWS SYNOPSIS: Johannes Kepler formulated three laws that described how the planets orbit around the Sun. His work paved the way for Isaac Newton, who derived
More informationMath 234. What you should know on day one. August 28, You should be able to use general principles like. x = cos t, y = sin t, 0 t π.
Math 234 What you should know on day one August 28, 2001 1 You should be able to use general principles like Length = ds, Area = da, Volume = dv For example the length of the semi circle x = cos t, y =
More informationAstronomy 111 Midterm #1
Astronomy 111 Midterm #1 Prof. Douglass 11 October 2018 Name: You may consult only one page of formulas and constants and a calculator while taking this test. You may not consult any books, digital resources,
More informationResonance In the Solar System
Resonance In the Solar System Steve Bache UNC Wilmington Dept. of Physics and Physical Oceanography Advisor : Dr. Russ Herman Spring 2012 Goal numerically investigate the dynamics of the asteroid belt
More informationToday. Planetary Motion. Tycho Brahe s Observations. Kepler s Laws of Planetary Motion. Laws of Motion. in physics
Planetary Motion Today Tycho Brahe s Observations Kepler s Laws of Planetary Motion Laws of Motion in physics Page from 1640 text in the KSL rare book collection That the Earth may be a Planet the seeming
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 informationGravitation and the Waltz of the Planets
Gravitation and the Waltz of the Planets Chapter Four Guiding Questions 1. How did ancient astronomers explain the motions of the planets? 2. Why did Copernicus think that the Earth and the other planets
More informationGravitation and the Waltz of the Planets. Chapter Four
Gravitation and the Waltz of the Planets Chapter Four Guiding Questions 1. How did ancient astronomers explain the motions of the planets? 2. Why did Copernicus think that the Earth and the other planets
More informationLecture Outline. Chapter 13 Gravity Pearson Education, Inc. Slide 13-1
Lecture Outline Chapter 13 Gravity Slide 13-1 The plan Lab this week: exam problems will put problems on mastering for chapters without HW; will also go over exam 2 Final coverage: now posted; some sections/chapters
More informationAPPENDIX B SUMMARY OF ORBITAL MECHANICS RELEVANT TO REMOTE SENSING
APPENDIX B SUMMARY OF ORBITAL MECHANICS RELEVANT TO REMOTE SENSING Orbit selection and sensor characteristics are closely related to the strategy required to achieve the desired results. Different types
More informationAP Physics C Textbook Problems
AP Physics C Textbook Problems Chapter 13 Pages 412 416 HW-16: 03. A 200-kg object and a 500-kg object are separated by 0.400 m. Find the net gravitational force exerted by these objects on a 50.0-kg object
More informationNewtonian Mechanics. May 16, 2005 Daniel Stump. MSU summer school
May 16, 2005 Daniel Stump CBI course PHY 321 MSU summer school Newtonian Mechanics 1 HISTORY Isaac Newton solved the premier scientific problem of his day, which was to explain the motion of the planets.
More information5. Universal Laws of Motion
5. Universal Laws of Motion If I have seen farther than others, it is because I have stood on the shoulders of giants. Sir Isaac Newton (164 177) Physicist Image courtesy of NASA/JPL Sir Isaac Newton (164-177)
More informationGravitation and the Motion of the Planets
Gravitation and the Motion of the Planets 1 Guiding Questions 1. How did ancient astronomers explain the motions of the planets? 2. Why did Copernicus think that the Earth and the other planets go around
More informationIntroduction to Classical Mechanics
Mechanics 1 Introduction to Classical Mechanics 1 HISTORY Isaac Newton solved the premier scientific problem of his day, which was to explain the motion of the planets. He published his theory in the famous
More informationAM 205 Final Project The N-Body Problem
AM 205 Final Project The N-Body Problem Leah Birch Elizabeth Finn Karen Yu December 14, 2012 Abstract The N-Body Problem can be solved using a variety of numeric integrators. Newton s Law of Universal
More informationUnit: Planetary Science
Orbital Motion Kepler s Laws GETTING AN ACCOUNT: 1) go to www.explorelearning.com 2) click on Enroll in a class (top right hand area of screen). 3) Where it says Enter class Code enter the number: MLTWD2YAZH
More informationGravitation. Makes the World Go Round
Gravitation Makes the World Go Round Gravitational Force The Force of gravity is an attractive force felt between all objects that have mass. G=6.67x10-11 N m 2 /kg 2 Example 1: What is the Force of Gravity
More informationPHYS-2010: General Physics I Course Lecture Notes Section VIII
PHYS-2010: General Physics I Course Lecture Notes Section VIII Dr. Donald G. Luttermoser East Tennessee State University Edition 2.4 Abstract These class notes are designed for use of the instructor and
More informationCelestial Orbits. Adrienne Carter Ottopaskal Rice May 18, 2001
Celestial Orbits Adrienne Carter sillyajc@yahoo.com Ottopaskal Rice ottomanbuski@hotmail.com May 18, 2001 1. Tycho Brache, a Danish astronomer of the late 1500s, had collected large amounts of raw data
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 informationA study upon Eris. I. Describing and characterizing the orbit of Eris around the Sun. I. Breda 1
Astronomy & Astrophysics manuscript no. Eris c ESO 2013 March 27, 2013 A study upon Eris I. Describing and characterizing the orbit of Eris around the Sun I. Breda 1 Faculty of Sciences (FCUP), University
More informationPerihelion Precession in the Solar System
International Journal of Multidisciplinary Current Research ISSN: 2321-3124 Research Article Available at: http://ijmcr.com Sara Kanzi * Hamed Ghasemian! * Currently pursuing Ph.D. degree program in Physics
More informationLecture 2c: Satellite Orbits
Lecture 2c: Satellite Orbits Outline 1. Newton s Laws of Mo3on 2. Newton s Law of Universal Gravita3on 3. Kepler s Laws 4. Pu>ng Newton and Kepler s Laws together and applying them to the Earth-satellite
More informationCopyright 2009, August E. Evrard.
Unless otherwise noted, the content of this course material is licensed under a Creative Commons BY 3.0 License. http://creativecommons.org/licenses/by/3.0/ Copyright 2009, August E. Evrard. You assume
More informationChapter 5 Centripetal Force and Gravity. Copyright 2010 Pearson Education, Inc.
Chapter 5 Centripetal Force and Gravity v Centripetal Acceleration v Velocity is a Vector v It has Magnitude and Direction v If either changes, the velocity vector changes. Tumble Buggy Demo v Centripetal
More informationAS3010: Introduction to Space Technology
AS3010: Introduction to Space Technology L E C T U R E S 8-9 Part B, Lectures 8-9 23 March, 2017 C O N T E N T S In this lecture, we will look at factors that cause an orbit to change over time orbital
More information