Targeting Mars. 19th June University of Ljubljana. Faculty of Mathematic and Physics. Seminar: 4. year-university program
|
|
- Ellen Flowers
- 5 years ago
- Views:
Transcription
1 University of Ljubljana Faculty of Mathematic and Physics Targeting Mars Seminar: 4. year-university program 19th June 016 Author: Jo²t Vrabi Koren Advisor: Prof. Dr. Tomaº Zwitter Ljubljana, June 016 ABSTRACT Spacetravel to Mars is an actual topic. I start with a two body problem, followed by a N-body problem. For a solution an integrator is needed. I describe the IAS15 integrator which was used for all calculations of the orbits. I then present an approximation for determining initial conditions for a rocket destined to travel to Mars and reach an areocentric orbit around Mars or even land.
2 Contents 1 Introduction Two body problem 3 Orbital elements 4 4 N-body problem and integrator 5 5 Targeting Mars 7 6 Results 10 7 Conclusion 10 1 Introduction The rst theoretical proposal of space travel using rockets was published in 1861 by Scottish astronomer and mathematician William Leitch. In the middle of the 0th century rst satellites were launched into space. Since then we have launched several thousand spacecrafts. Some of them have landed on other planets and some are orbiting other panets in our solar system. Typical journey of a spacecraft destined to other planets is consisting of launch from the Earth, a geocentric orbit, transfer to a heliocentric transfer orbit and transfer to a areocentric orbit (orbit around a planet). One of the main tasks of planning a mission is determining ight parameters (initial velocity, time of launch, end velocity,...) with n-body simulations using powerful n-body integrators. Mars is the next human landing destination. That is why it is extremely popular. It has several spacecraft already in its orbit, observing its surface and characteristics. Landing of a sonde on Mars was also already achieved, several times (Viking, Curiosity, Pathnder,...). Two body problem Figure 1: Particles with masses m 1 and m, distance r apart. Let us start with a conservative central force, where the potential is V(r), a function of r only, which means that the force is always along r. We take a particle with reduced mass m, moving about a xed center of force, which is taken as the origin of the coordinate system. Reduced mass is m = m1m m 1+m, [1] where m 1 and m are masses of our bodies (Figure 1). Since potential energy is dependent on radial distance only, the problem has spherical symmetry. Total angular momentum vector
3 L = r p, [] is then conserved. It follows that that r is always perpendicular to L, which has a xed direction in space. This can only be true if r lies in a plane with a normal parallel to L. Central force motion is thus always restricted to a plane. Motion of a particle in space is described by three coordinates. In spherical coordinate system, these are azimuth angle θ, the zenizh angle ψ and radial distance r. We choose a system with a polar axis in the direction of L so the motion is restricted to the plane perpendicular to the polar axis. ψ then has the constant value π / and can be dropped from the following sections. Conservation of total energy can be written as: Here we dene gravitational potential V (r) as: It can be shown that the equation of the orbit for our problem is: 1 l mṙ + 1 mr + V = konst. = E. [3] V (r) = k r. [4] 1 r = mk ( El² cos(θ l² mk² θ )). [5] where θ is one of the angles where the orbit turns. Eccentricity is dened as: e = 1 + El mk. [6] Orbits can be classied according to eccentricity e: e > 1 E > 0 hyperbola e = 1 E = 0 parabola e < 1 E < 0 ellipse e = 0 E = mk l circle Let us take a more detailed look at ellipse, as this is the orbit that is the most important for this seminar. r 1 in r are two apsidal distances (Figure ). Figure : Ellipse with apsidal distances r 1 and r and semimajor axis a. These two points corespond to angle θ where distance is the smallest. The radial velocity at these points is zero. It can be shown that the major axis depends solely upon the energy: 3
4 We can now rewrite the elliptical orbit equation [6] as a = r1+r = k E. [7] r = a(1 e²) 1+e cos(θ θ ). [8] Let us describe the motion of a particle on an elliptic orbit in time. The relation between the radial distance and the time is given by the enrgy equation [3], which we rewrite as ˆ m r t = dr k. [9] r l² +E This integral is most conveniently calculated with a new variable ψ. It is dened by the relation r 0 mr² r = a(1 e cos ψ). [10] By comparing this equation to equation [5], we see that ψ covers the interval 0 to π, as θ goes trough a complete revolution. Perihelion occurs at ψ = 0 (also θ = 0) and aphelion at ψ = π = θ. With this new variable (plus some algebra and substitutions), we can now rewrite equation [9] as t = ma³ k ψ (1 e cos ψ)dψ. [11] 0 If we integrate this equation over full range of π, we get the expression for the period of an elliptic movement τ = πa 3/ m k. [1] Let us remember that a particle moving around the Sun is a two body problem, so we have to take both masses into account. We substitute m with the reduced mass from equation [1], where mass m 1 is the particle mass and mass m Solar mass. Further, the gravitational law of attraction is Equation [1] then becomes f = G m1m r². [13] πa τ = 3/ πa3/ G(m1+m ) Gm, [14] if we neglect the mass of the particle compared to the Sun. In this seminar, as we are dealing with spacecrafts, this can be done. But planetary mass is not always negligible compared to the Sun - mass of Jupiter, for example, is approximately 0,1% of the mass of the Sun. 3 Orbital elements We need 6 parameters to determine an orbit. In cartesian coordinates, this are radial vector r and velocity vector v in a chosen system. In our case that is the Solar system with Sun in the center, with z coordinate equal to zero in Earths ecliptic plane and x axis pointing towards vernal equinox. In astronomy, an elliptical orbit is often described with orbital elements a, e, i, ω, Ω and T : ˆ Semimajor axis (a). 4
5 ˆ ˆ ˆ ˆ Eccentricity (e). Inclination (i) - vertical tilt (in degrees) of the ellipse with respect to the reference plane, measured at the ascending node (where the orbit passes upward through the reference plane). Longitude of the ascending node (Ω) - horizontally orients the ascending node of the ellipse (where the orbit passes upward through the reference plane) with respect to the reference frame's vernal point. Argument of periapsis (ω) - denes the orientation of the ellipse in the orbital plane, as an angle measured from the ascending node to the periapsis. ˆ Time at epoch (T) - denes the position of the orbiting body along the ellipse at a specic time. Figure 3: Orbital elements. Note that true anomaly in our case is θ. True anomaly at epoch is related to Time at epoch (T). Source: 4 N-body problem and integrator The n-body problem considers N point masses m i, i=1,...n in an inertial reference frame in three dimensional space moving under the inuence of mutual gravitational attraction. Position vector q i is appointed to each mass m i. Gravitational force on mass m i by mass m j is given as and sum over al masses gives us the complete force on particle i: F i = F ij = Gmimj(qj qi) q j q i 3, [15] N j=1,j i Gm im j(q j q i) q j q i 3. (16] 5
6 This system cannot be solved analiticaly, so we have to use numerical integration. In calculations presented in later chapters, I used Rebound package for python that uses IAS15 integrator (Ref.[4]) for N-body problems. Here I will describe this integrator. Let us start with the equation y = F [y, y, t], (17) with y being the acceleration of a particle and F the specic force, depending on the velocity of the particle y, position of the particle y, and time t. Let us expand equation (17) into series: y [t] y 0 + a 0 t + a 1 t²+...+a 6 t 7. (18) The constant term is the force at the beginning of a timestep, y 0 = y [0] = F [t = 0]. By introducing the step size dt, h = t/dt and b k = a k dt k+1, the series becomes y [h] y 0 + b 0 h + b 1 h² b 6 h 7. (19) Coecients b k have dimensions of acceleration, since h is dimensionless. By integrating equation (19) y [h] y 0 + hdt(y 0 + h (b 0 + h 3 (b 1+...))), (0) we get an estimate of the velocity y during and at the end of the timstep. By integrationg once more y[h] y 0 + y 0hdt + h²dt² (y 0 + h 3 (b 0 + h (b 1+...))), (1) we get an estimate of the positions. The rst two terms are the position and the velocity at the beginning of the timestep. The approximation of this integral is very accurate as we choose Gauss- Radau spacing. It is closely related to Gaussian quadrature for approximating an integral and uses the starting point at h=0. We use a quadrature with 8 function evaluations to construct a 15th-order scheme. Now, we have to nd good estimates for the coecients b k. We can resolve this with a predictor corrector scheme. We get an estimate for the position and velocities to calculate the forces which is the predictor. In the rst iteration, we set all b k = 0, which corresponds to a particle moving along a path of constant acceleration. We then use these forces to calculate better estimates, which is the corrector. This process is iterated until the positions and velocities converge below machine precision. The number of predictor-corrector iterations is determined dynamically. The iterations converges with iterations in most cases. 1 iterations are set as the upper limit to prevent innite loops. If it has not converged by then, the timestep is too large. The timestep can be set to be choosen automatically. Authors claim this integrator is superior to all other integrators (symplectin and non-symplectic) and one of the main reasons is adaptive timestep (Ref. [4]). This is perticularly usefull when dealing with close encounters. 6
7 5 Targeting Mars Let us now imagine we have a rocket orbiting Earth on a circular ecliptic orbit. We want to increase its velocity at the right time, so we can escape Earth's gravitation and begin our journey to space. We want to choose the correct direction and magnitude of the velocity, so we can reach Mars. For this seminar, I simplied the problem. As we have seen in previous chapters, all orbits in a gravitational eld are conic sections in 3D space. The planets have almost elliptical orbits. First I asumed the orbits of Mars and Earth are circular. Eccentricites of Earth's and Mars's orbits are e Earth = and e Mars = We see that Earth's orbit is in fact almost circular (Fig [4]). Figure 4: Left: Comparison of real and circular approximation orbits of Earth and Mars in the same plane (i = 0 ). Right: Two possible starting positions and velocities of our rocket. Then I asumed the orbits are coplanar. Inclinations, in degrees, of the Earth's orbit is 0, and of the Mars's orbit i Mars = Orbital period of Mars is P Mars = days and of Earth P Earth = days. From equations [14], we see that semimajor axis is a = 3 GM Sun P ² 4π, () where M Sun is the mass of Sun. In the case of a circular orbit, a is radius of rotation r 0. By our asumption we get r 0Mars = a.e. and r 0Earth = a.e.. Our rocket is circling the Earth. We turn on the engines, when its velocity is maximum. This occurs when the rocket velocity vector is in the same direction as the Earth velocity vector v r v Earth. We have two dierent solutions (Fig [4]), depending on the rotation direction. Choosing magnitude is not so simple. We rewrite the energy equation [3] as ε = v rel² GM Mars r rel, [3] where ε is now energy per unit mass, v rel velocity of our rocket relative to Mars, r rel distance of the rocket from Mars and M Mars mass of Mars. We want the rocket to get close to Mars and dene close as F Mars F Sun > 10, [4] 7
8 where F Mars is the gravitational force of Mars magnitude and F Sun of the Sun. This way we can neglect the eects of the Sun when we are close and the equation [3] is then correct. The aproach orbit turns out to be nearly hyperbolic (with respect to Mars) and hyperbolic orbits have positive energy. We can then lower rocket's velocity, so the energy becomes negative and we get a closed orbit around Mars. So our desire is that the energy is as low as possible as we come close to Mars. This way we are closest to a closed orbit. The velocity magnitude of a particle on an elliptical orbit is greatest at perihelium and lowest at aphelium. That means we want to reach Mars at aphelium, so the rst term in equation [3] is minimal. Figure 5: Orbits representing velocities v a, v b and v n and Mars. Now, we choose two initial velocities. One so that aphelium is inside the Mars orbit and the other so that it is outside (Figure 4). The rst one is chosen to be the sum of Earth's escape velocity and Earth's velocity: GM v a = v escape + v Earth = Earth r Earth + v Earth, [5] where M Earth is Earth mass and r Earth distance of our rocket from the center of Earth. The second velocity is chosen to be slightly smaller than Sun escape velocity: GM v b < Sun r Sun, [6] where r Sun is the distance of our rocket to the center of Sun. Velocities v a and v b are chosen so they have a physical meaning. We could have ofcourse chosen two dierent velocities, as long as one would correspond to an orbit with an aphelion inside the Mars orbit and the other with aphelion outside. The following procedure is some kind of bisection. We appoint a new velocity v N = va+v b, [7] and send our spacecraft toward Mars. We nd the aphelium position and check if its distance from the Sun is greater or smaller than r 0Mars. If it is greater we assign v b = v N and if it is smaller v a = v N. We continue untill we achieve the desired accuracy for the initial velocity v 0 (Figure 5). This is all 8
9 carried out in a system with the Sun in the origin and consisting only of Earth and our spacecraft (no Mars). Figure 6: System without Mars showing paths of rockets, starting in a circular orbit around the Earth at 8000 km. Left one starts closer to the Sun and right one further from the Sun than Earth. Figure 7: System with Mars showing a rocket with initial condition v 0 and t, calculated by bisection. Rocket was red when it was further from the Sun than Earth. Let us now introduce Mars into our system. We start by putting it at the point of orbit intersection and then take it an angle ωt back on its orbit. t is the travel time the spacecraft. We have come close to Mars, but not close enough (Figure 6). To come closer, we have to solve a two-dimensional problem as we have two parameters now, initial velocity v 0 and travel time t, which also denes the angle of the position of Mars. It was shown before that we are searching for the energy minimum. We choose suciently small v 0 and t and start searching for the minimum. Eventually we get there (Figure 7). 9
10 Figure 8: Sheme of an energy minimum. We start with initial conditions v 0 and t and make our way to minimum. Note that we rst have to get close to Mars (Eq.[4]), before we start searching for the minimum. That means we rst have to get close to the gravitational minimum. This is done in the same fashion as our search for energy minimum. 6 Results I used Nasa's Horizons database to nd coordinates in time of Mars Atmosphere and Volatile EvolutioN Mission (MAVEN) and Mars Orbiter Mission (MOM). This are two Mars orbiting sondes, launched from Earth in November 013. They were rst in heliocentric orbits and then accelerated to reach v 0 and begin their journey toward Mars. When reaching v 0, MOM was approximately 7900km from the center of the Earth and MAVEN approximately 8000 km. Our rocket is at 8000 km. The main parameter I compared here is the change of velocity v final near Mars, needed to put the selected spacecraft on a closed orbit around Mars. This change denes fuel consumption at the end of the travel. MOM MAVEN Rocket under Earth Rocket above Earth time of travel 98 days 308 days 7 days 81 days v 0 40,1 km/s 40, km/s 40,5 km/s 40,6 km/s v final 0,31 km/s 0,39 km/s 0,48 km/s 0,46 km/s Table 1: Comparison of dierent spacecraft. v final was calculated via equation [3]. 7 Conclusion We have shown that by making several approximations to our problem, we still get comparable results. This was however expected, as the orbits of Earth and Mars are relatively circular and lying almost in the same plane. It also indicates that studying an approximated problem can lead us to initial conditions, which can be then used on real orbits in real time. 10
11 References [1] H. Goldstein, C. Poole and J. Safko, 00, CLASSICAL MECHANICS, third edition, Pearson Education, San Francisco, 638 p. [] [3] H. Rein and S. Liu, 011, REBOUND:An open-source multi-purpose N-body code for collisional dynamics, Astronomy & Astrophysics, Volume 537, id.a18, 10 pp. [4] H. Rein and D.S. Spiegel, 015, IAS15: A fast, adaptive, high order integrator for gravitational dynamics, accurate to machine precision over billion orbits, Monthly Notices of the Royal Astronomical Society, Volume 446, Issue, p [5] [6] [7]
ANNEX 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 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 informationorbits Moon, Planets Spacecrafts Calculating the and by Dr. Shiu-Sing TONG
A Science Enrichment Programme for Secondary 3-4 Students : Teaching and Learning Resources the and Spacecrafts orbits Moon, Planets Calculating the 171 of by Dr. Shiu-Sing TONG 172 Calculating the orbits
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 informationOrbits in Geographic Context. Instantaneous Time Solutions Orbit Fixing in Geographic Frame Classical Orbital Elements
Orbits in Geographic Context Instantaneous Time Solutions Orbit Fixing in Geographic Frame Classical Orbital Elements Instantaneous Time Solutions Solution of central force motion, described through two
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 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 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 informationAST111, Lecture 1b. Measurements of bodies in the solar system (overview continued) Orbital elements
AST111, Lecture 1b Measurements of bodies in the solar system (overview continued) Orbital elements Planetary properties (continued): Measuring Mass The orbital period of a moon about a planet depends
More informationExperimental Analysis of Low Earth Orbit Satellites due to Atmospheric Perturbations
Experimental Analysis of Low Earth Orbit Satellites due to Atmospheric Perturbations Aman Saluja #1, Manish Bansal #2, M Raja #3, Mohd Maaz #4 #Aerospace Department, University of Petroleum and Energy
More informationCelestial Mechanics and Satellite Orbits
Celestial Mechanics and Satellite Orbits Introduction to Space 2017 Slides: Jaan Praks, Hannu Koskinen, Zainab Saleem Lecture: Jaan Praks Assignment Draw Earth, and a satellite orbiting the Earth. Draw
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 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 informationGravitational Potential Energy and Total Energy *
OpenStax-CNX module: m58347 Gravitational Potential Energy and Total Energy * OpenStax This work is produced by OpenStax-CNX and licensed under the Creative Commons Attribution License 4.0 By the end of
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 informationAstromechanics. 6. Changing Orbits
Astromechanics 6. Changing Orbits Once an orbit is established in the two body problem, it will remain the same size (semi major axis) and shape (eccentricity) in the original orbit plane. In order to
More informationFundamentals of Satellite technology
Fundamentals of Satellite technology Prepared by A.Kaviyarasu Assistant Professor Department of Aerospace Engineering Madras Institute Of Technology Chromepet, Chennai Orbital Plane All of the planets,
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 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 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 informationASE 366K Spacecraft Dynamics
ASE 366K Spacecraft Dynamics Homework 2 Solutions 50 Points Total: 10 points each for 1.16, 1.19, 2.6, 2.7, and 10 points for completing the rest. 1.13 Show that the position vector is a min or max at
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 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 informationCelestial Mechanics Lecture 10
Celestial Mechanics Lecture 10 ˆ This is the first of two topics which I have added to the curriculum for this term. ˆ We have a surprizing amount of firepower at our disposal to analyze some basic problems
More 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 informationKeplerian Elements Tutorial
Keplerian Elements Tutorial This tutorial is based on the documentation provided with InstantTrack, written by Franklin Antonio, N6NKF. Satellite Orbital Elements are numbers that tell us the orbit of
More informationORBITAL CHARACTERISTICS DUE TO THE THREE DIMENSIONAL SWING-BY IN THE SUN-JUPITER SYSTEM
ORBITAL CHARACTERISTICS DUE TO THE THREE DIMENSIONAL SWING-BY IN THE SUN-JUPITER SYSTEM JORGE K. S. FORMIGA 1,2 and ANTONIO F B A PRADO 2 National Institute for Space Research -INPE 1 Technology Faculty-FATEC-SJC
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 informationChapter 4. Motion and gravity
Chapter 4. Motion and gravity Announcements Labs open this week to finish. You may go to any lab section this week (most people done). Lab exercise 2 starts Oct 2. It's the long one!! Midterm exam likely
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 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 informationChapter 2: Orbits and Launching Methods
9/20/ Chapter 2: Orbits and Launching Methods Prepared by Dr. Mohammed Taha El Astal EELE 6335 Telecom. System Part I: Satellite Communic ations Winter Content Kepler s First, Second, and Third Law Definitions
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 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 informationarxiv: v1 [physics.pop-ph] 29 May 2009
Gravity assist in 3D like in Ulysses mission Herbert Morales Escuela de Física, Universidad de Costa Rica, San José, Costa Rica arxiv:0905.4788v1 [physics.pop-ph] 29 May 2009 (Dated: October 22, 2018)
More informationGravitation. Luis Anchordoqui
Gravitation Kepler's law and Newton's Synthesis The nighttime sky with its myriad stars and shinning planets has always fascinated people on Earth. Towards the end of the XVI century the astronomer Tycho
More informationAP Physics QUIZ Gravitation
AP Physics QUIZ Gravitation Name: 1. If F1 is the magnitude of the force exerted by the Earth on a satellite in orbit about the Earth and F2 is the magnitude of the force exerted by the satellite on the
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 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 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 informationMAE 180A: Spacecraft Guidance I, Summer 2009 Homework 2 Due Tuesday, July 14, in class.
MAE 180A: Spacecraft Guidance I, Summer 2009 Homework 2 Due Tuesday, July 14, in class. Guidelines: Please turn in a neat and clean homework that gives all the formulae that you have used as well as details
More informationLecture D30 - Orbit Transfers
J. Peraire 16.07 Dynamics Fall 004 Version 1.1 Lecture D30 - Orbit Transfers In this lecture, we will consider how to transfer from one orbit, or trajectory, to another. One of the assumptions that we
More informationPW-Sat two years on orbit.
13th of February 2014 is the second anniversary of launch of the first polish student-made satellite PW-Sat. Currently Students' Space Association on Warsaw University of Technology is working on another
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 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 informationThe Astrodynamics and Mechanics of Orbital Spaceflight
The Astrodynamics and Mechanics of Orbital Spaceflight Vedant Chandra 11-S1, TSRS Moulsari 1 1 Introduction to Rocketry Before getting into the details of orbital mechanics, we must understand the fundamentals
More informationOrbital Mechanics! Space System Design, MAE 342, Princeton University! Robert Stengel
Orbital Mechanics Space System Design, MAE 342, Princeton University Robert Stengel Conic section orbits Equations of motion Momentum and energy Kepler s Equation Position and velocity in orbit Copyright
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 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 informationExtending the Patched-Conic Approximation to the Restricted Four-Body Problem
Monografías de la Real Academia de Ciencias de Zaragoza 3, 133 146, (6). Extending the Patched-Conic Approximation to the Restricted Four-Body Problem Thomas R. Reppert Department of Aerospace and Ocean
More informationAstronomy Section 2 Solar System Test
is really cool! 1. The diagram below shows one model of a portion of the universe. Astronomy Section 2 Solar System Test 4. Which arrangement of the Sun, the Moon, and Earth results in the highest high
More informationAstrodynamics (AERO0024)
Astrodynamics (AERO0024) 10. Interplanetary Trajectories Gaëtan Kerschen Space Structures & Systems Lab (S3L) Motivation 2 6. Interplanetary Trajectories 6.1 Patched conic method 6.2 Lambert s problem
More informationSatellite Orbital Maneuvers and Transfers. Dr Ugur GUVEN
Satellite Orbital Maneuvers and Transfers Dr Ugur GUVEN Orbit Maneuvers At some point during the lifetime of most space vehicles or satellites, we must change one or more of the orbital elements. For example,
More information1 Summary of Chapter 2
General Astronomy (9:61) Fall 01 Lecture 7 Notes, September 10, 01 1 Summary of Chapter There are a number of items from Chapter that you should be sure to understand. 1.1 Terminology A number of technical
More informationPrevious Lecture. Orbital maneuvers: general framework. Single-impulse maneuver: compatibility conditions
2 / 48 Previous Lecture Orbital maneuvers: general framework Single-impulse maneuver: compatibility conditions closed form expression for the impulsive velocity vector magnitude interpretation coplanar
More informationProjectile Motion. Conceptual Physics 11 th Edition. Projectile Motion. Projectile Motion. Projectile Motion. This lecture will help you understand:
Conceptual Physics 11 th Edition Projectile motion is a combination of a horizontal component, and Chapter 10: PROJECTILE AND SATELLITE MOTION a vertical component. This lecture will help you understand:
More informationConceptual Physics 11 th Edition
Conceptual Physics 11 th Edition Chapter 10: PROJECTILE AND SATELLITE MOTION This lecture will help you understand: Projectile Motion Fast-Moving Projectiles Satellites Circular Satellite Orbits Elliptical
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 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 informationAP Physics Multiple Choice Practice Gravitation
AP Physics Multiple Choice Practice Gravitation 1. Each of five satellites makes a circular orbit about an object that is much more massive than any of the satellites. The mass and orbital radius of each
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 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 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 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 informationCelestial Mechanics III. Time and reference frames Orbital elements Calculation of ephemerides Orbit determination
Celestial Mechanics III Time and reference frames Orbital elements Calculation of ephemerides Orbit determination Orbital position versus time: The choice of units Gravitational constant: SI units ([m],[kg],[s])
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 informationCelestial Mechanics I. Introduction Kepler s Laws
Celestial Mechanics I Introduction Kepler s Laws Goals of the Course The student will be able to provide a detailed account of fundamental celestial mechanics The student will learn to perform detailed
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 informationChapter 13: universal gravitation
Chapter 13: universal gravitation Newton s Law of Gravitation Weight Gravitational Potential Energy The Motion of Satellites Kepler s Laws and the Motion of Planets Spherical Mass Distributions Apparent
More informationSpacecraft Dynamics and Control
Spacecraft Dynamics and Control Matthew M. Peet Arizona State University Lecture 5: Hyperbolic Orbits Introduction In this Lecture, you will learn: Hyperbolic orbits Hyperbolic Anomaly Kepler s Equation,
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 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 informationNewton s Laws of Motion and Gravity ASTR 2110 Sarazin. Space Shuttle
Newton s Laws of Motion and Gravity ASTR 2110 Sarazin Space Shuttle Discussion Session This Week Friday, September 8, 3-4 pm Shorter Discussion Session (end 3:40), followed by: Intro to Astronomy Department
More informationHyperbolic-Type Orbits in the Schwarzschild Metric
Hyperbolic-Type Orbits in the Schwarzschild Metric F.T. Hioe* and David Kuebel Department of Physics, St. John Fisher College, Rochester, NY 468 and Department of Physics & Astronomy, University of Rochester,
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 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 informationAstronomy 111, Fall October 2011
Astronomy 111, Fall 011 4 October 011 Today in Astronomy 111: asteroids, perturbations and orbital resonances Leftovers: proofs of Kepler s second and third laws Elliptical orbits and center of mass More
More informationGRAVITATION. F = GmM R 2
GRAVITATION Name: Partner: Section: Date: PURPOSE: To explore the gravitational force and Kepler s Laws of Planetary motion. INTRODUCTION: Newton s law of Universal Gravitation tells us that the gravitational
More informationName: Earth 110 Exploration of the Solar System Assignment 1: Celestial Motions and Forces Due on Tuesday, Jan. 19, 2016
Name: Earth 110 Exploration of the Solar System Assignment 1: Celestial Motions and Forces Due on Tuesday, Jan. 19, 2016 Why are celestial motions and forces important? They explain the world around us.
More informationPHYS 160 Astronomy Test #1 Fall 2017 Version B
PHYS 160 Astronomy Test #1 Fall 2017 Version B 1 I. True/False (1 point each) Circle the T if the statement is true, or F if the statement is false on your answer sheet. 1. An object has the same weight,
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 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 informationOutline for Today: Newton s Law of Universal Gravitation The Gravitational Field Orbital Motion Gravitational Potential Energy. Hello!
PHY131H1F - Class 13 Outline for Today: Newton s Law of Universal Gravitation The Gravitational Field Orbital Motion Gravitational Potential Energy Under the Flower of Kent apple tree in the Woolsthorpe
More informationCHAPTER 8 PLANETARY MOTIONS
1 CHAPTER 8 PLANETARY MOTIONS 8.1 Introduction The word planet means wanderer (πλάνητες αστέρες wandering stars); in contrast to the fixed stars, the planets wander around on the celestial sphere, sometimes
More informationlightyears observable universe astronomical unit po- laris perihelion Milky Way
1 Chapter 1 Astronomical distances are so large we typically measure distances in lightyears: the distance light can travel in one year, or 9.46 10 12 km or 9, 600, 000, 000, 000 km. Looking into the sky
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 informationPSI AP Physics 1 Gravitation
PSI AP Physics 1 Gravitation Multiple Choice 1. Two objects attract each other gravitationally. If the distance between their centers is cut in half, the gravitational force A) is cut to one fourth. B)
More information1 The Solar System. 1.1 a journey into our galaxy
1 The Solar System Though Pluto, and the far-flung depths of the Solar System, is the focus of this book, it is essential that Pluto is placed in the context of the planetary system that it inhabits our
More informationPHYSICS 1030 Homework #9
PHYSICS 1030 Homework #9 (Due Dec. 6, 2017) Find the position of the planet Mars at time t D December 6, 2017, 5:00 am EST. You will do this by following the steps shown below. (a) Convert the time t to
More informationPhysical Science 1 Chapter 16 INTRODUCTION. Astronomy is the study of the universe, which includes all matter, energy, space and time.
INTRODUCTION Astronomy is the study of the universe, which includes all matter, energy, space and time. Although the universe is vast and almost beyond imagination, much is known about its make-up and
More informationGravitation Part I. Ptolemy, Copernicus, Galileo, and Kepler
Gravitation Part I. Ptolemy, Copernicus, Galileo, and Kepler Celestial motions The stars: Uniform daily motion about the celestial poles (rising and setting). The Sun: Daily motion around the celestial
More informationLecture Module 2: Spherical Geometry, Various Axes Systems
1 Lecture Module 2: Spherical Geometry, Various Axes Systems Satellites in space need inertial frame of reference for attitude determination. In a true sense, all bodies in universe are in motion and inertial
More informationASTR-1010: Astronomy I Course Notes Section III
ASTR-1010: Astronomy I Course Notes Section III Dr. Donald G. Luttermoser Department of Physics and Astronomy East Tennessee State University Edition 2.0 Abstract These class notes are designed for use
More informationAnnouncements. Topics To Be Covered in this Lecture
Announcements! Tonight s observing session is cancelled (due to clouds)! the next one will be one week from now, weather permitting! The 2 nd LearningCurve activity was due earlier today! Assignment 2
More informationNewton s Legacy. 1- accelerate to break free of Earth. Rocket Science: How to send a spacecraft to Mars
Reading: today: web-based reading on satellite orbits; Chap. 3 Sec. 5 Chap. 7, Sect. 1, 2 (for next week) Exam 1: Tuesday, September 26, 6:45-8:00. Room assignments on course website ESSAY QUESTION Homework
More informationThe Three Body Problem
The Three Body Problem Joakim Hirvonen Grützelius Karlstad University December 26, 2004 Department of Engineeringsciences, Physics and Mathematics 5p Examinator: Prof Jürgen Füchs Abstract The main topic
More informationDynamics of the Earth
Time Dynamics of the Earth Historically, a day is a time interval between successive upper transits of a given celestial reference point. upper transit the passage of a body across the celestial meridian
More informationNAVIGATION & MISSION DESIGN BRANCH
c o d e 5 9 5 National Aeronautics and Space Administration Michael Mesarch Michael.A.Mesarch@nasa.gov NAVIGATION & MISSION DESIGN BRANCH www.nasa.gov Outline Orbital Elements Orbital Precession Differential
More information2 so that the (time-dependent) torque, as measured in space-xed coordinates, is ~N = ma2! 2 (^x cos!t +^y sin!t) 10 = a p 5 2 ^z ~ F 2 ~F 1 = a p 5^z
1 PHYS 321 Homework Assignment #9 Due: Friday, 22 November 2002 (4 probs) 1. Problem 7-14, B&O p. 280. The situation is as shown above. We take the body-xed axes ^x 0 and ^y 0 to lie in the plane of the
More information