Basic Orbits and Mission Analysis

1. Basic Orbits and Mission Analysis

2. Newton's Law of Gravitation To find the velocity increment required for various missions, we must calculate trajectories and orbits. This is done using Newton's Law of Gravitation: Here the lhs is the "radial force" of attraction due to gravitation, between two bodies; the big one of mass M, and the little one of mass m. The universal gravitational constant G is 6.670 E-11 Nm2/kg2.

3. Keplerian Orbits Johannes Kepler (1609: first law, 1619: third law)  Kepler’s First Law: The orbit of each planet is an ellipse, with the sun at one focus.  Second Law: The line joining the planet to the sun sweeps out equal areas in equal times. Third Law: The square of the period of a planet is proportional to the cube of its mean distance from the sun

4. Kepler's Laws applied to satellite of mass m, orbiting a much bigger object of mass M. i.e, m << M): phy-061062.blogspot.com/ 1.The satellite travels in an elliptical path around its center of attraction, which is located at one of focii of the ellipse. The orbit must lie in a plane containing the center of attraction.

5. Second Law: The line joining the planet to the sun sweeps out equal areas in equal times. Implication: The radius vector from the center of attraction sweeps equal areas of the orbit per unit time. As the satellite moves away, its speed decreases. As it nears the center of attraction, its speed increases. Found on the web at a secondary site. Primary developer of this application is unknown

6. Third Law: The square of the period of a planet is proportional to the cube of its mean distance from the sun Implication: The ratio of the squares of the orbital periods of any two satellites about the same body equals the ratio of the cubes of the semi-major axes of the respective orbits.

7. Satellite Equations of Motion Newton’s Second Law Law of Gravitation F=magnitude of force due to gravity (N) G=Universal constant of gravity = M=mass of Earth (kg) M =satellite mass R =distance from center of Earth to the satellite (km) = GM for Earth = 398,600.5 Equation for satellite acceleration (from the above two Newton laws) (2-body equation of motion) The solution to this is the polar equation of a conic section.

8. Conic Energy Semi-major axis a e Circle <0 = radius 0 Ellipse <0 >0 0<e<1 Parabola 0 1 Hyperbola >0 <0 >1 Conic Section Orbits / Trajectories Why Conic Sections are called conic sections: Total specific mechanical Energy of an orbit: (mechanical energy per unit mass) where a is the semi major axis and V is magnitude of velocity. Source: http://www.bartleby.com/61/imagepages/A4conics.html Courtesy Houghton-Miflin Co.

9. A Few Commonly-Used Results on Orbits r in km is radius of orbit Circular orbit around Earth: V in km/s; r in km Escape Velocity from Earth: Speed at any point in an elliptical orbit: Time for one orbit is: Specific angular momentum: is constant for the 2-body problem, per Kepler (no plane change)

10. Orbit Maneuvers where V is the velocity at that point in the corresponding orbit. The assumption made usually in simplifying orbit maneuvers is that we are dealing with a 2-body problem. For transfers between earth orbit and say, Mars orbit, we assume initially an Earth-vehicle problem. When the vehicle leaves the vicinity of Earth, we assume a Sun-vehicle problem, and finally a Mars-vehicle problem Hohmann transfer orbit: Transfer orbit’s ellipse is tangent to initial and final orbits, at the transfer orbit’s pengee and apogee respectively  Most efficient transfer between two circular orbits. Two-burn transfer between circular orbits at altitudes : DVB rB rA DVA

11. Hohmann Transfer DV for a two-burn transfer between circular orbits at altitudes . For Earth, with delta-v in km/s, Example: When the radii of the two orbits are close, the delta-v’s are approximately equal.

12. One-Tangent Burn This is a quicker orbit transfer than the Hohman transfer ellipse, but is more costly in delta-v. DVB DVA

13. Orbit Plane Changes Simple plane change through angle q at an orbital velocity of Vi. Plane changes are often combined with orbit-change burns. Usual practice is to do a little of the plane change at the initial burn (at the lower orbit where velocity is higher) and most of the plane change at the final burn at the higher altitude where the velocity is lower. Plane change combined with Hohman transfer.

14. Launch from Earth Note that the objective is to achieve the required tangential velocity for orbit – the need to gain height above Earth is mainly dictated by the need to get out of the atmosphere and minimize drag losses. As the vehicle trajectory begins to slant over, gravity pulls the trajectory down. At low speeds, this requires a substantial pitch angle between the instantaneous velocity vector and the thrust direction to counter gravity. Gravity Turn 1.Begin ascent with short vertical rise, during which it starts usually and completes a roll from launch-pad azimuth to desired flight azimuth. 2.Pitch-over (pitch because angular change of axis occurs faster than change of velocity vector from vertical) 3.Adjust altitude to keep in trajectory through aerodynamic load region.

15. Gravity turn gldt glCos g dt V(t) -dg V(t+dt) g Local horizontal If the angle of the trajectory with respect to the local horizontal is gamma (g) then, this angle changes as (rad/s) (rad.) or

16. Controlled Angle of Attack – Linear Tangent Steering Angle of attack or the tangent of the thrust angle with respect to the local horizontal, are controlled as functions of time to minimize losses. Angle of attack can be increased as dynamic pressure decreases beyond the “max q” point, subject to limitations on vehicle bending moments. Tangent of the thrust angle usually starts at a high value and is reduced according to a linear or other specified schedule, i.e., the slope B is usually negative. .

17. Launch Trajectories

19. Acceleration During Launch

20. Interplanetary Trajectories Earth-Mars Hohmann scheme • Escape velocity from Earth – to an orbit around the Sun, with zero velocity with • respect to Earth. • Impulsive burn to get into transfer orbit. • Impulsive burn at Mars orbit to match velocity with Mars. DUMars Change in orbital energy per unit mass (w.r.t. Earth) DUT

21. Data on Solar System Destinations Mass of Earth = 5.975 E24 kg. Mean Diameter = 12742. 46km

22. Earth-Mars Oberth Scheme Escape from Earth while carrying enough propellant to produce the transfer velocity is wasteful in work against gravity. Can be avoided by achieving the right velocity magnitude and direction in the first maneuver. The launch velocity is larger at Earth, directly placing the vehicle in the Hohman transfer ellipse. Thus Complications: 1) Atmospheric drag 2) Orbit desired @ Mars

23. Faster Trajectories for Human-Carrying Missions • Four impulsive burns: • Escape from Earth towards Mars with an excess over the Hohman transfer velocity. • Increment to get captured into Mars orbit: Most trajectories put the vehicle in Mars orbit, moving slower around the Sun than Mars – need a burn to get captured. • Opposite of the above to escape Mars into a transfer trajectory towards Earth • Burn to decelerate to a safe Earth re-entry velocity.

24. Values for Sample Missions: Garrison & Stuhler 1 Astronautical Unit = mean distance from Earth to Sun = 149,558,000km Impulsive burns assumed in all the above.

25. Spiral Transfer: Low, Continuous Thrust Delta v required is greater than that for Hohman transfer, but a small engine can be used. Also, possibility of using solar energy as the heat source. Thrust applied along the circumferential direction is then given by or with initial conditions at where Solution shows that the delta v required with a very small thrust-to-EarthWeight ratio is about 2.4 times the impulsive-burn delta v. As thrust-to-weight approaches 1, this comes down close to 1.

26. Spiral Transfer: Earth-Moon “SMART-1 to continue alternating between long and short thrust arcs until 15 Nov lunar capture: orbital period 4 days, 19 hrs; apogee altitude has increased from 35,880 km to 200,000 km; lunar commissioning, science operations begin in mid-Jan 2005” http://sci.esa.int/

27. Gravity Assist Maneuvers http://science.howstuffworks.com/framed.htm? parent=question102.htm&url=http://www.jpl.nasa.gov /basics/bsf4-1.html Vehicle approaches another planet with velocity greater than escape velocity for that planet. Planet’s gravity adds momentum to the vehicle in desired direction. With respect to the planet, vehicle gains speed on approach and decelerates as it moves away. But when the planet’s velocity is considered, there is net gain in angular momentum with respect to the Sun.

29. Orbital Elements a: semi-major axis. e: eccentricity . i: inclination of the angular momentum vector from the Z-axis. (Z is usually Earth’s axis of rotation. W: right ascension of the ascending node: Angle from the vernal equinox to the “ascending node.” “node” : one of two points where the satellite orbit intersects the X-Y (equatorial) plane. Ascending node is the point where the satellite passes through the equatorial plane moving from south to north. Right ascension is measured in degrees of right-handed rotation about the pole (Z axis). w: argument of perigee. The angle from the ascending node to the eccentricity vector, measured in degrees in the direction of the satellite’s position. The eccentricity vector points from the center of the Earth to the perigee with magnitude equal to the eccentricity of the orbit, in degrees. • : “true anomaly”. Angle from the eccentricity vector to satellite’s position vector, measured in the direction of the satellite motion. Can also use the time since perigee passage, in degrees.