Tag: Orbital Mechanics

Earth-Mars Free Return

Let’s try using the Fortran Astrodynamics Toolkit and Pikaia to solve a real-world orbital mechanics problem. In this case, computing the optimal Earth-Mars free return trajectory in the year 2018. This is a trajectory that departs Earth, and then with

Tagged with: , , , , ,

IAU Rotation Models

The IAU Working Group on Cartographic Coordinates and Rotational Elements (WGCCRE) is the keeper of official models that describe the cartographic coordinates and rotational elements of planetary bodies (such as the Earth, satellites, minor planets, and comets).  Periodically, they release a report containing the coefficients

Tagged with: , , , ,

Nonsingular Geopotential Models

The gravitational potential of a non-homogeneous celestial body at radius , latitude , and longitude can be represented in spherical coordinates as: Where is the radius of the body, is the gravitational parameter of the body, and  are spherical harmonic coefficients,

Tagged with: , , , ,

Speeding up SPICE

The SPICE Toolkit software is an excellent package of very well-written and well-documented routines for a variety of astrodynamics applications.  It is produced by NASA’s Navigation and Ancillary Information Facility (NAIF).  Versions are available for Fortran 77, C, IDL, and Matlab. To speed up

Tagged with: , ,

Rocket Equation

The rocket equation describes the basic principles of a rocket in the absence of external forces.  Various forms of this equation relate the following fundamental parameters: engine thrust (), the engine specific impulse (), the engine exhaust velocity (), the initial mass

Tagged with: , ,

Porkchop Plots

A porkchop plot is a data visualization tool used in interplanetary mission design which displays contours of various quantities as a function of departure and arrival date.  Example pork chop plots for 2016 Earth-Mars transfers are shown here.  The x-axis is the Earth departure

Tagged with: , , ,

Lambert’s Problem

Lambert’s problem is to solve for the orbit transfer that connects two position vectors in a given time of flight.  It is one of the basic problems of orbital mechanics, and was solved by Swiss mathematician Johann Heinrich Lambert. A standard

Tagged with: , ,

Object Oriented Runge-Kutta Module

Fortran 2003 was a significant update to Fortran 95 (probably something on the order of the update from C to C++).  It brought Fortran into the modern computing world with the addition of object oriented programming features such as type extension

Tagged with: , ,

SPICE N65 Released

JPL just released a new version of the SPICE Toolkit (N65).  According to the What’s New page, changes include: support for some new environments and termination of some old environments new Geometry Finder (GF) interfaces — illumination angle search, phase angle search, user-defined

Tagged with: , ,

Fortran Astrodynamics Toolkit

I’m starting a new project on GitHub: the Fortran Astrodynamics Toolkit. Hardly anyone is developing open source orbital mechanics software for modern Fortran, so the time has come. Most of the code from this blog will eventually find its way

Tagged with: , , ,