3

THE ORBITAL ELEMENTS 

In this chapter, we will introduce, describe, and present mathematical equations to find the Classical Orbital Elements, COEs, which are used to locate a satellite in orbit. So far, we have learned to describe where satellites are in orbit by it position, [latex]\vec{R}[/latex] and velocity, [latex]\vec{V}[/latex] vectors. However, these are constantly changing as a satellite moves around its orbit and it is difficult to actually “visualize” where it is.

The provided image illustrates the geometry associated with the analytical solution for the two-body equation of motion in terms of polar coordinates. The equation expresses the position vector (R) as a function of the semimajor axis (a), eccentricity (e), true anomaly (ν), and the cosine of the true anomaly. The formula, R = a(1 - e^2) / (1 + e * cos ν), is applicable when the eccentricity (e) falls within the range of 0 to 1. The depiction highlights the geometrical relationships involved in describing the orbit of a satellite around a central body using these polar coordinates.
Figure 1: Associated Geometry for the Two-Body Equation of Motion (Source: Siegenthaler and Saylor, 2017).

In this chapter we will reference these vectors to an approximate inertial coordinate frame. We will then learn how these can be used to understand the COEs six parameters that completely describe an orbit and a satellite’s location in it.

1. Its size, a, or semi-major axis, which was introduced in chapter 2.

2. Its shape, e, or eccentricity, which was also introduced in chapter 2.

3.  Its orbital tilt, or how it is oriented relative to the earth’s equator.

4.  Its orientation, or swivel, locating the orbital plane in space.

5.  The orientation of the orbit within its orbital plane.

6.  The location of the satellite within its orbit, or true anomaly, which was also introduced in chapter 2

We will first introduce the coordinate frame used for most earth-orbiting satellites, the geocentric equatorial coordinate system, then introduce each of the Classical Orbital Elements relative to that fixed coordinate frame.

THE GEOCENTRIC EQUATORIAL COORDINATE SYSTEM

As we discussed in chapter 2, when the two-body equation of motion was derived, we found that it was based on Newton’s Second Law. To be valid, Newton’s laws must be based on an inertial, or non-accelerating, coordinate frame.

A reference frame is a collection of unit vectors, or vectors with a magnitude of one, at right angles to each other. In the figure below, the black axis would be the coordinate frame, with the unit vectors, [latex]\hat{I}[/latex], [latex]\hat{J}[/latex], and [latex]\widehat{K}[/latex] each 90 degrees apart.

In the provided image, a reference frame is depicted as a coordinate frame with three axes: (i) (along the x-axis in red), (j) (along the y-axis in green), and (k) (along the z-axis in blue). These unit vectors are at right angles to each other, forming a three-dimensional Cartesian coordinate system. The image specifically illustrates a right-handed coordinate frame, where the (k) axis is determined by the cross product of the (i) and (j) vectors. The right-hand rule is applied to determine the direction of the (k) vector, emphasizing the conventions used in physics and engineering for establishing a consistent and universally accepted coordinate system.
Figure 2: A Right-Handed Coordinate Frame (Source: Nav Ganti, 2020).

A right-handed coordinate frame is one in which the third axis, [latex]\widehat{K}[/latex] in our case, is equal to the cross product of the other two vectors; i.e.

[latex]\hat{I} \times \hat{J} = \widehat{K}[/latex]

The geocentric-equatorial coordinate system is used for earth-orbiting satellites. It has its origin at the center of the earth. The fundamental plane is the earth’s equatorial plane and the principle axis, [latex]\hat{I}[/latex], points in the vernal equinox direction, [latex]\Upsilon[/latex].  The z-axis, [latex]\widehat{K}[/latex], points in the direction of the North Pole. The [latex]\hat{J}[/latex] axis completes the right-handed coordinate frame.  The coordinate system is shown in the figure below.

In the provided image, the geocentric-equatorial coordinate system is illustrated. This coordinate system is employed for referencing earth-orbiting satellites. The origin is situated at the center of the Earth. The fundamental plane is aligned with the Earth's equatorial plane, denoted by the (i) axis. The principal axis, (j), points in the vernal equinox direction, and the (k) axis points towards the North Pole, thus completing the right-handed coordinate frame. This system serves as a crucial reference for specifying orbital elements, providing a standardized way to describe the positions and movements of satellites in Earth's orbit.
Figure 3: The Geocentric Equatorial Coordinate Frame (Source: FAA, 2015).

So what is the Vernal Equinox direction [latex]\Upsilon[/latex]? Just as the Greenwich meridian has been arbitrarily chosen as the zero point for measuring longitude on the surface of the Earth, the first point of Aries has been chosen as the zero point in the celestial sphere. It is called the first point of Aries because it was originally in the Zodiac constellation, Aries, as shown in the figure below. This is where the star was located in 150 B.C. when Ptolemy first mapped the constellations. Due to precession of the equinoxes it moved from Aries into Pisces in about 70 B.C., and is now approaching Aquarius (Martin, 2009).

The first point of Aries, serving as the zero point in the celestial sphere, was initially located in the Zodiac constellation Aries around 150 B.C., as depicted in the figure provided. This point, similar to the Greenwich meridian's role in measuring longitude on Earth's surface, serves as a reference for celestial coordinates. The choice of this point as the starting reference allows astronomers to specify the positions of celestial objects and map the celestial sphere. Over time, due to the precession of the Earth's axis, the first point of Aries has shifted relative to the background stars, but its historical significance in celestial coordinate systems persists.
Figure 4: The Aries Constellation (Source: Jess Waid, 2013).

It is also found by drawing a line from the earth through the Sun on the first day of spring, as shown in the figure below (Kwantus, 2004).

The image depicts the celestial coordinate system with the first point of Aries as the reference, marked by the intersection of the celestial equator and the ecliptic. This point is a fundamental reference in astronomy and serves as the starting point for measuring celestial longitude. Over centuries, the precession of the Earth's axis has caused a shift in the location of the first point of Aries relative to the background stars, but it remains a crucial reference for celestial navigation and mapping.
Figure 5: First Point of Aries on the Vernal Equinox (Source: Kwantus, 2004).

The reason we use this star as our reference is that it is very, very far away (over 160 light-years) from our solar system. As the earth travels around the sun throughout the year, the vernal equinox direction will remain essentially fixed in the same direction, allowing astrodynamics to use this as an inertial coordinate frame.

INTRODUCTION TO THE CLASSICAL ORBITAL ELEMENTS

For any physical object we generally need six parameters to describe its location in three-dimensional space. There are three degrees of freedom, x, y, and z as shown in the figure below, and three roll parameters about each axis.

The image illustrates the six degrees of freedom of a quadcopter, representing its ability to move in different directions. The six degrees of freedom are often described as follows: (a) A quadcopter can move longitudinally (forward and backward), vertically (upward and downward), and laterally (right and left). (b) It can also move rotationally among each axis to produce roll, pitch, and yaw movements. These degrees of freedom provide the quadcopter with maneuverability and control in various directions, allowing it to navigate through three-dimensional space.
Figure 6: A Quadcopter Moving (a) Longitudinally, Vertically, and Laterally. Among Each Axis, It Can Produce (b) Roll, Pitch, and Yaw Movements (Source: ResearchGate, 2017).

For example, if you are give a position and velocity vector for a spacecraft, it will have six components in the [latex]\hat{I}[/latex], [latex]\hat{J}[/latex], [latex]\widehat{K}[/latex] coordinate frame:

[latex]\vec{R} = R_{x} \hat{I}+ R_{y} \hat{J} + R_{z} \widehat{K}[/latex]

[latex]\vec{V} = V_{X} \hat{I} + V_{Y} \hat{J} + V_{Z}[/latex]

However, this does not help much in visualizing the spacecraft in its orbit. You cannot tell much about the orbit’s size or shape, nor where the satellite is located in it. Kepler once again came to the rescue and defined a set of six classical orbital elements, sometimes called the Keplerian elements, that were introduced at the beginning of this chapter. We will further investigate these elements one-by-one and see what they tell use about a satellite’s orbit and the satellite’s location in it.

SEMI-MAJOR AXIS, a

The first orbital element we will consider describes the orbit’s size.

Recall the semi-major axis measures half the distance across the major (long) axis of the orbit, as shown in the figure below. It is a constant parameter in the absence of forces other than gravity since the orbit’s specific mechanical energy does not change.

The presented image guides the process of integrating the two-body equation of motion for an elliptical orbit, specifically focusing on the orbit's semi-major axis (v) and acceleration vector (R). The image encourages recalling the relevant parameters associated with an elliptical orbit, emphasizing the need to establish a coordinate system. The goal is to numerically integrate the second-order nonlinear vector differential equation to unveil the shape of the elliptical orbit. The depicted graphic serves as a visual aid to understand the steps involved in this analytical process.
Figure 7: Coordinate System Within an Elliptical Orbit.

In chapter 2, we related the size of the orbit to its specific mechanical energy using the relationship:

[latex]\varepsilon = - \frac{\mu}{2a}[/latex]

Where:

[latex]\varepsilon[/latex] = specific mechanical energy (km2/s2)

[latex]\mu[/latex] = gravitational parameter of central body (km3/s2)

a = semi-major axis (km)

Given a position and velocity vector for a satellite, we can then find specific mechanical energy of the orbit using the relationship:

[latex]\varepsilon = \frac{V^{2}}{2} - \frac{\mu}{R}[/latex]

Where:

V = magnitude of the velocity vector, [latex]\vec{V}[/latex], in km/s

[latex]\mu[/latex] = gravitational parameter of central body (km3/s2)

R = magnitude of the position vector, [latex]\vec{R}[/latex], in km

Recall the range of values for each orbit type:

Orbit type Semi-major axis value Comments
Circular a = R >(6378 + 200) km R radius
Elliptical a such that Rp > 6378 km Rp = radius of perigee
Parabolic a = [latex]\infty[/latex]
Hyperbolic a < 0

ECCENTRITY, e

Eccentricity describes the shape of an orbit. It is the ratio of the distance between the two foci and the length of the major axis.

[latex]e = \frac{c}{a}[/latex]

Thus, it is a unitless quantity. As shown in the figure below, the eccentricity tells you how “out of roundness” the conic section is:

The provided image illustrates the geometry associated with the analytical solution for the two-body equation of motion in terms of polar coordinates. The equation expresses the position vector (R) as a function of the semimajor axis (a), eccentricity (e), true anomaly (ν), and the cosine of the true anomaly. The formula, R = a(1 - e^2) / (1 + e * cos ν), is applicable when the eccentricity (e) falls within the range of 0 to 1. The symbol Beta () refers to the angle between the position and velocity vectors. The depiction highlights the geometrical relationships involved in describing the orbit of a satellite around a central body using these polar coordinates.
Figure 8: Associated Geometry for the Two-Body Equation of Motion With Emphasis on Angle Beta (Source: Siegenthaler and Saylor, 2017).

However, this relationship does not provide a convenient way to get the eccentricity from the position and velocity vectors.

The eccentricity of an orbit is also a vector, [latex]\vec{e}[/latex], which points from the center of the earth to perigee and whose magnitude is equal to the eccentricity, e. It is related to the position and velocity vectors using one of two equations, which are equivalent:

[latex]\vec{e} = \frac{1}{\mu}[(V^{2} - \frac{\mu}{R})\vec{R} - (\vec{R} \times \vec{V})\vec{V}][/latex]

Or

[latex]\vec{e} = \frac{1}{\mu}(\vec{V} \times \vec{h}) - \frac{\vec{R}}{\left| \vec{R} \right|}[/latex]

Where:

V = magnitude of the velocity vector, [latex]\vec{V}[/latex], in km/s

[latex]\mu[/latex] = gravitational parameter of central body in km3/s2

R = magnitude of the position vector, [latex]\vec{R}[/latex], in km

[latex]\vec{h} = \vec{R} \times \vec {V}[/latex]= specific angular momentum vector in km2/s

The eccentricity is then the magnitude of the vector, or

[latex]e = \left| \vec{e} \right|[/latex]

Alternately, if the Radius of Perigee and semi-major axis, a, are known, you can also find the magnitude of the eccentricity vector using the relationship:

[latex]R_{p} = a(1-e)[/latex]

Or if Radius of Apogee and semi-major axis, a, are known, you can find the magnitude of the eccentricity vector using the relationship:

[latex]R_{a} = a(1+e)[/latex]

Recall the range of values for each orbit type:

Orbit Type Eccentricity Value(s)
Circular  e = 0
Elliptical 0 < e < 1
Parabolic e = 1
Hyperbolic e > 1

 

INCLINATION, i

The inclination indicates the tilt of an orbit. It is the angle measured between the [latex]\widehat{K}[/latex] axis and the angular momentum vector, [latex]\vec{h}[/latex], as shown in the figure below:

The image depicts the concept of inclination in orbital mechanics. Inclination represents the tilt of an orbit and is measured as the angle between the axis and the angular momentum vector (h). The angular momentum vector is a vector quantity pointing along the axis of rotation and is a fundamental parameter in describing the orientation of an orbit. The inclination angle provides information about how much an orbit deviates from a reference plane, typically the celestial equator or the ecliptic plane. The figure helps visualize how the inclination angle is determined in relation to the angular momentum vector and the reference axis.
Figure 9: The Geocentric Equatorial Coordinate Frame Showcasing Orbit Inclination (Source: FAA, 2018).

You might be wondering why the angle between the equatorial plane and the orbital plane is not used to describe inclination. Recall from the mathematics review in chapter 2 that the angle between two vectors [latex]\vec{A}[/latex] and [latex]\vec{B}[/latex] is given by:

[latex]\cos\theta = \frac{\vec{A} \times \vec{B}}{AB}[/latex]

The image illustrates the concept of measuring the inclination of an orbit using the angle between two vectors. Instead of directly using the angle between the equatorial plane and the orbital plane, the inclination is determined by the angle between two specific vectors: the orbital angular momentum vector and the unit vector along the axis perpendicular to the equatorial plane. As mentioned in the mathematics review, the angle between two vectors can be calculated using the dot product. This approach offers a convenient and mathematically sound way to quantify the tilt of an orbit without directly relying on the angles between planes.
Figure 10: An Angle Between Vectors A and B.

Thus this provides a convenient way to measure the inclination as it is the angle between the two known vectors. The angular moment vector, [latex]\vec{h}[/latex] is given by:

[latex]\vec{h} = \vec{R} \times \vec{V}[/latex]

The [latex]\widehat{K}[/latex] vector is the unit vector that points from the center of the earth to the north pole, and it is give by [latex]\widehat{K} = 0 \hat{I} + 0\hat{J} + 1\widehat{K}[/latex]. Hence, the inclination is given by:

 

[latex]cos i = \frac{\vec{h} \times \widehat{K}}{hK}[/latex]

Since,

 

[latex]\vec{h} \times \widehat{K} = h_{I}(0) + h_{J}(0) + h_{k}(1)[/latex]

And K is the magnitude of the [latex]\widehat{K}[/latex] vector, or 1, so this expression becomes:

[latex]i = \cos^{-1}(\frac{h_{k}}{h})[/latex]

Where:

hk = the kth component of the angular momentum vector, [latex]\vec{h}[/latex]

h = the magnitude of the angular momentum vector, [latex]\vec{h}[/latex]

The inclination, i, only takes values between 0° and 180°. An explanation of the types of orbits is found in the figure below:

The image illustrates different types of orbital inclinations in the context of orbital mechanics. Inclination refers to the tilt of an orbit relative to a reference plane. The first type shown is a low inclination, where the orbit is nearly aligned with the equatorial plane. The second type is a moderate inclination, indicating a moderate tilt relative to the equatorial plane. The third type is a high inclination, representing a polar orbit, where the orbit is nearly perpendicular to the equatorial plane, passing over the Earth's poles. Understanding these variations in inclination is essential for planning satellite trajectories for different purposes.
Figure 11: Types of Orbits and Their Inclination Table (Source: FAA, 2015).

An orbit with an inclination of 0° is called an equatorial orbit because satellites in these orbits travel around the earth’s equator. A Satellite at a geostationary altitude, as calculated in example 2 in chapter 2, is in an equatorial orbit. This allows it to remain above the same spot on the earth as its eastward velocity is approximately equal to the the velocity of the ground directly beneath it. An orbit with an inclination of 180° could theoretically be used, but it would be travelling in the opposite direction of the earth’s rotation and would be highly inefficient to launch into. A satellite in a polar orbit has an inclination of 90° and hence travels over the north and south geographical poles. This type of orbit allows for maximum viewing of all of the latitudes of the earth.

Satellites can also be categorized as prograde or retrograde. Prograde orbits have inclinations between 0° and 90° and travel in the direction of the earth’s rotation. These are, by far, the most prevalent types of orbits. They are more efficient to launch into because they already have the benefit of the earth’s rotation at the launch site to give the satellite an eastward boost. Retrograde orbits have inclinations greater than 90°. Retrograde orbits with greater than about 110° inclinations are rare as satellite in them travel against the earth’s rotation. There is, however, a well-used orbital inclination at about 98° to 100°, called a sun-synchronous orbit. We will learn later how sun-synchronous orbits are designed and how they are very valuable for earth-observations.

Inclination applies to all orbit types and conic sections. Here is a summary of what we know so far about inclinations:

Inclination

0

Type

Equatorial

 

 

Comments

Travels above Earth’s equator

If at geo altitude, is geostationary

 

180 Equatorial Not normally used
90 Polar Covers all latitudes

 

0 – 90 Prograde Satellite travels in the same direction the earth’s rotates

 

90 – 180 Retrograde Satellite travels in the opposite direction the earth’s rotates

Seldom use orbital inclinations greater than 110 degrees

 

RIGHT ASCENSION OF ASCENDING NODE, RAAN, [latex]\Omega[/latex]

The RAAN indicates the swivel or orientation of an orbit. It is the angle measured between the [latex]\hat{I}[/latex] axis and the ascending node vector, or the vector pointing to the location where the satellite passes from the southern to northern hemisphere. First, though, we must define the ascending node vector, [latex]\vec{n}[/latex], sometimes called the line of nodes, as shown in the figure below:

[latex]\vec{n} = \widehat{K} \times \vec{h}[/latex]

The image depicts the concept of the ascending node in orbital mechanics, where the ascending node is a critical element in determining the orientation of an orbit. In the illustration, the gray plane represents the orbital plane of a celestial body, while the equatorial plane serves as the reference plane for measuring orbital elements. The ascending node (AN) is identified as the point where the orbit intersects from below to above the reference plane, crucial in defining the orbit's orientation. Additionally, the ascending node vector, depicted by the arrow pointing upward from the ascending node, signifies the direction perpendicular to both the orbital and reference planes. This vector is essential for understanding and specifying the orbit's position relative to the equatorial plane of the celestial body.
Figure 12: RAAN Omega With Ascending Node Vector (Source: FAA, 2015).

The ascending node vector lies along the intersection of two planes – the orbital plane and the equatorial plane.

RAAN, or [latex]\Omega[/latex], is then the angle measured from the vernal equinox direction, [latex]\hat{I}[/latex], to the ascending node vector, [latex]\vec{n}[/latex]. Hence, RAAN is given by:

[latex]\cos\Omega = \frac{\hat{I} \times \vec{n}}{In}[/latex]

And since the [latex]\hat{I}[/latex] vector is the unit vector that points from the center of the earth to the vernal equinox, hence is give by [latex]\hat{I} = 1\hat{I} + 0\hat{J} + 0\widehat{K}[/latex]. Hence, RAAN is given by:

 

[latex]\Omega = \cos^{-1}(\frac{n_{i}}{n})[/latex]

Where:

ni = the ith component of the ascending node vector, [latex]\vec{n}[/latex]

n = the magnitude of the ascending node vector, [latex]\vec{n}[/latex]

RAAN can take any value between 0° and 360° so a half plane check is needed. As can be seen in the figure above, if the jth component of the ascending node vector is negative, then RAAN is in the second half plane. This can be written mathematically as:

If nj < 0 then 180 < [latex]\Omega[/latex] < 360 degrees.

As described in the Math review section, remember that for a cosine term, to find the angle in the second half plane, use 360° minus the angle in the first half plane, or:

[latex]\Omega[/latex] second half plane = 360 degrees – [latex]\Omega[/latex] first half plane

If the ascending node vector points along the vernal equinox direction, then [latex]\Omega[/latex] = 0.The RAAN in the figure shown above is about 45°. In the figure below, it is about 90° and in the second figure below, it is about 270°. Our next COE will be referenced from the ascending node vector.

The image illustrates the concept of the Right Ascension of the Ascending Node (RAAN) in orbital mechanics. The gray plane represents the orbital plane of a celestial body, while the equatorial plane serves as the reference plane. The Ascending Node (AN) is the point where the orbital plane intersects the reference plane, transitioning from below to above. The RAAN vector, depicted by the arrow, signifies the direction of RAAN, an angle measured along the celestial equator from a reference point to the ascending node. In this specific configuration, RAAN equals 90 degrees, indicating that the ascending node is positioned at a right angle from the reference point along the celestial equator. This orientation is crucial for defining the tilt and positioning of an orbit in the celestial sphere. Understanding RAAN helps specify the alignment of the orbital plane within the broader astronomical context.
Figure 13: RAAN Omega With Ascending Node Vector at 90 Degrees (Source: Orbits and Trajectories, 2005).
The image depicts key elements in orbital mechanics, showcasing the concept of orbital inclination. The gray plane represents the orbital plane of a celestial body, with the equatorial plane serving as the reference plane. Inclination (i), the angle between the orbital and reference planes, is measured at the ascending node (AN), the point where the orbit crosses from below to above the reference plane. The yellow plane represents the ecliptic plane of Earth's orbit around the Sun. The inclination vector (i Vector) illustrates the orbit's tilt. Specifically, the image indicates that the Right Ascension of the Ascending Node (RAAN) is approximately 270°, signifying the angle at which the orbit intersects the reference plane. Understanding these elements is crucial for describing the orientation and tilt of an orbit in celestial mechanics.
Figure 14: RAAN Omega With Ascending Node Vector at 270 Degrees (Source: ResearchGate, 2007).

Note [latex]\Omega[/latex] is not defined if there is no ascending node, in other words if the orbit is equatorial with an inclination of 0° or 180°.

 

ARGUMENT OF PERIGEE, [latex]\omega[/latex]

The argument of perigee, [latex]\omega[/latex] indicates theorientation of the orbit, or where its perigee is located. It is the angle measured between the ascending node vector, [latex]\vec{n}[/latex], and the location of perigee, at which the eccentricity vector, [latex]\vec{e}[/latex], points.

The image presents an illustration related to celestial mechanics, specifically focusing on the argument of perigee (arg perigee) in orbital dynamics. The diagram features an orbital plane, and the key parameter, argument of perigee, is highlighted. The argument of perigee is the angle measured in the orbital plane between the ascending node and the perigee of the orbit. This parameter is essential for characterizing the orientation of an orbit, providing information about the positioning of the closest approach (perigee) to the central body. Understanding the argument of perigee contributes to a comprehensive description of orbital elements and behavior in space.
Figure 15: An Angle, Omega, Between the Ascending Node and Perigee (Source: FAA, 2015).

Hence, argument of perigee is given by:

[latex]\cos\omega = \frac{\vec{n} \times \vec{e}}{ne}[/latex]

Or

[latex]\omega = \cos^{-1} (\frac{\vec{n} \times \vec{e}}{ne})[/latex]

 

Argument of perigee can take any value between 0° and 360° so a half plane check is needed. As can be seen in the figure above, if the kth component of the eccentricity vector is negative, then argument of perigee is in the second half plane. This can be written mathematically as:

If ek < 0 then 180 < [latex]\omega[/latex] < 360 degrees.

 

If perigee is at the ascending node, then [latex]\omega = 0[/latex]. In the figure above, argument of perigee is about 90°. In the figure below, [latex]\omega[/latex] is about 270°.

The image depicts the concept of the argument of perigee in celestial mechanics, specifically with a focus on a scenario where the argument of perigee is approximately 270 degrees. In orbital dynamics, the argument of perigee represents the angle within the orbital plane between the ascending node and the perigee of the orbit. In this illustration, the orbital plane is portrayed, and the angle measuring 270 degrees in the argument of perigee is highlighted. This information is crucial for understanding the orientation and geometry of an orbit, providing insights into the position of the orbit's closest approach (perigee) relative to the central body.
Figure 16: An Argument of Perigee at 270 Degrees (Source: GlobalSecurity.org, 2003).

Note [latex]\omega[/latex] is not defined if there is no ascending node, in other words if the orbit is equatorial with an inclination of 0° or 180°. It is also not defined if the orbit is circular, or has an eccentricity of 0, since there is no perigee in a circular orbit.

 

TRUE ANOMALY, [latex]\nu[/latex]

The last COE is the true anomaly, [latex]\nu[/latex], which indicates the location of the satellite in the orbit. It is the angle measured between the eccentricity vector, [latex]\vec{e}[/latex], and the position vector to the satellite, [latex]\vec{R}[/latex] , as shown in the figure below:

The presented image guides the process of integrating the two-body equation of motion for an elliptical orbit, specifically focusing on the orbit's semi-major axis (v) and acceleration vector (R). The image encourages recalling the relevant parameters associated with an elliptical orbit, emphasizing the need to establish a coordinate system. The goal is to numerically integrate the second-order nonlinear vector differential equation to unveil the shape of the elliptical orbit. The depicted graphic serves as a visual aid to understand the steps involved in this analytical process.
Figure 17: Coordinate System Within an Elliptical Orbit.

Hence, true anomaly is given by:

[latex]\cos\nu = \frac{\vec{e} \times \vec{R}}{eR}[/latex]

 

Or

[latex]\nu = \cos^{-1} (\frac{\vec{e} \times \vec{R}}{eR})[/latex]

 

True anomaly can take any value between 0° and 360° so a half plane check is needed. This check is not quite as easy to visualize as the others. First, recall that the angle between two vectors is given by:

[latex]\cos\theta = \frac{\vec{A} \times \vec{B}}{AB}[/latex]

 

In the case of the [latex]\vec{R}[/latex] and [latex]\vec{V}[/latex] vectors, this relationship becomes:

 

[latex]\vec{R} \times \vec{V} = RV\cos\theta[/latex]

where [latex]\theta[/latex]is the angle between the [latex]\vec{R}[/latex] and [latex]\vec{V}[/latex] vectors, which can be related to the flight path angle by [latex]\phi = 90 - \theta[/latex].  Hence:

 

[latex]\vec{R} \times \vec{V} = RV\sin\phi[/latex]

If the satellite is travelling from perigee to apogee, the angle between the [latex]\vec{R}[/latex] and [latex]\vec{V}[/latex] vectors is between 0° and 90° and [latex]\phi[/latex] is positive, hence sin [latex]\phi[/latex] is also positive. If the satellite is travelling from apogee to perigee, the angle between the [latex]\vec{R}[/latex] and [latex]\vec{V}[/latex] vectors is between 90° and 180° and [latex]\phi[/latex] is negative, hence sin [latex]\phi[/latex] is also negative. Thus the sign of the dot product, [latex]\vec{R} \times \vec{V}[/latex] provides a convenient way of determining which half plane the satellite is in.

The illustration presents the concept of the dot product in celestial mechanics, specifically concerning the relative position of a satellite in its orbit. The dot product involves the vectors R and V, where R represents the radial velocity vector and V is the velocity vector of the satellite. The angle between these vectors varies as the satellite moves from perigee to apogee or vice versa. When traveling from perigee to apogee, the angle is between 0° and 90°, resulting in a positive sine value, indicating a positive dot product. Conversely, when moving from apogee to perigee, the angle is between 90° and 180°, leading to a negative sine value and, consequently, a negative dot product. This dot product's sign serves as a convenient means to determine the satellite's position within the orbital half-plane.
Figure 18: Half Plane Check for True Anomaly.

The half plane check for true anomaly thus becomes:

If [latex]\vec{R} \times \vec{V}[/latex] < 0 then 180 < [latex]\nu[/latex] < 360 degrees.

If the satellite is at perigee, then the true anomaly is 0 and if the satellite is at apogee, then the true anomaly is 180°.

Note [latex]\nu[/latex] is not defined if the orbit is circular, or has an eccentricity of 0, since there is no perigee in a circular orbit.

The Classical Orbital Elements are summarized in the table below:

Element Name Description Range of values Undefined
a Semi-major Axis Orbital Size

[latex]a = \frac{\mu}{2\epsilon}[/latex]

Depends on Orbit Type Never
e Eccentricity      Orbit Shape

[latex]\small\vec{e} =\frac{1}{\mu}[(V^{2} - \frac{\mu}{R})\vec{R} - (\vec{R} \times \vec{V})\vec{V}][/latex]

[latex]\vec{e} =\frac{1}{\mu}(\vec{V} \times \vec{h})- \frac{\vec{R}}{\left| \vec{R} \right|}[/latex]

Depends on Orbit Type Never
i Inclination Orbital Tilt

[latex]i = \cos^{-1}(\frac{h_{k}}{h})[/latex]

[latex]0 \le i \le 180[/latex] Never
[latex]\Omega[/latex] Right Ascension of Ascending Node (RAAN)

Orbital Swivel

[latex]\Omega = \cos^{-1}(\frac{n_{i}}{n})[/latex]

If [latex]n_{j} \lt 0[/latex] then

[latex]180 \lt \Omega \lt 360[/latex]

 

[latex]0 \le \Omega \lt 360[/latex] When i = 0 or 180 degrees (equatorial orbit)

[latex]\Omega[/latex]

Argument of Perigee

Orbit Orientation

[latex]\omega = cos^{-1}(\frac{\vec{n} \times \vec{e}}{ne})[/latex]

If [latex]e_{k} \lt 0[/latex] then

[latex]180 \lt \omega \lt 360[/latex]

[latex]0 \le \omega \lt 360[/latex] When i = 0 or 180 degrees (equatorial orbit) or   e = 0 (circular orbit)
[latex]\nu[/latex] True Anomaly

Location of the Satellite

[latex]\nu = \cos^{-1}(\frac{\vec{e} \times \vec{R}}{eR})[/latex]

If [latex]\vec{R} \times \vec{V} \lt 0[/latex] then

[latex]180 \lt \nu \lt 360[/latex]

[latex]0 \le \nu \lt 360[/latex] When e = 0 (circular orbit)

In the next section, we will address what to use when one or more of the Classical Orbital Elements are not defied. But first, work through this example problem to see how well you understand the concepts from this section. The answers and solutions are at the end of this chapter.

EXAMPLE 1

[latex]\vec{R} = 0\hat{I} + 0\hat{J} + 10,000\widehat{K}[/latex]km

[latex]\vec{V} = 6\hat{I} + 0\hat{J} + 0\widehat{K}[/latex] km/s

Answer the following questions:

  1. Where is the satellite currently?
  2. What is the flight path angle at the satellite’s current position?
  3. What is the orbit’s specific angular momentum?
  4. What is the orbit’s semi-major axis?
  5. What is the orbit’s eccentricity?
  6. What is its inclination, i?
  7. What is its RAAN, [latex]\Omega[/latex]?
  8. What is its argument of perigee, [latex]\omega[/latex]?
  9. What is the satellite’s true anomaly, [latex]\nu[/latex]?

ALTERNATE ORBITAL ELEMENTS

As indicated in the table above, there are situations where one or more of the COEs are not defined. What do we do in those situations? In those cases, we use one of the alternate orbital elements.

Argument of latitude, u

When a satellite is in a circular orbit, its argument of perigee, [latex]\omega[/latex], and true anomaly, [latex]\mu[/latex], are not defined. In this case, we use the argument of latitude, u, which indicates the location of the satellite in the orbit. It is the angle measured between the ascending node vector, [latex]\vec{n}[/latex], and the position vector to the satellite, [latex]\vec{R}[/latex] , as shown in the figure below.

The diagram depicts the argument of latitude in the context of celestial mechanics. The argument of latitude is an orbital parameter describing the position of a satellite in its elliptical orbit around a celestial body. In the figure, the angle (u) is formed between the satellite's position vector, denoted by the arrow pointing from the celestial body to the satellite, and the satellite's velocity vector. The argument of latitude provides insight into the satellite's location within its orbit, offering a measure of its angular distance from the point of perigee. Understanding this parameter is crucial for predicting and analyzing the satellite's orbital behavior and dynamics.
Figure 19: An Argument of Latitude (u).

Hence, argument of latitude is given by:

[latex]\cos u = \frac{\vec{n} \times \vec{R}}{nR}[/latex]

 

Or

[latex]u = \cos^{-1} (\frac{\vec{n} \times \vec{R}}{nR})[/latex]

 

Argument of latitude can take any value between 0° and 360° so a half plane check is needed. As can be seen in the figure above, if the kth component of the position vector is negative, then argument of perigee is in the second half plane. This can be written mathematically as:

If Rk < 0 then 180 < u < 360 degrees.

In the figure above, argument of latitude is about 25°.

 

Longitude of Perigee

When a satellite is in an equatorial orbit, it has no ascending node vector, [latex]\vec{n}[/latex], and hence RAAN, [latex]\Omega[/latex] and argument of perigee, [latex]\omega[/latex] are undefined. In this case, we use the longitude of perigee, [latex]\Pi[/latex] which indicates the orientation of the orbit, or where its perigee is located. It is the angle measured between the vernal equinox direction, [latex]\hat{I}[/latex] and the eccentricity vector, [latex]\vec{e}[/latex], which points to perigee.  The longitude of perigee is show in the figure below.

The image illustrates the longitude of perigee in the context of celestial mechanics. The longitude of perigee is an orbital parameter that defines the angular distance along the celestial body's equator to the point of perigee in the satellite's orbit. In the figure, the celestial body is depicted at the center, and the orbital plane is represented by the horizontal plane. The line extending from the celestial body's center to the perigee point defines the longitude of perigee. This parameter is significant in characterizing the orientation and geometry of the satellite's elliptical orbit, providing valuable information for orbital analysis and prediction.
Figure 20: A Longitude of Perigee.

Hence, longitude of perigee is given by:

[latex]\cos\Pi = \frac{\hat{I} \times \vec{e}}{e}[/latex]

Or

[latex]\Pi = \cos^{-1}(\frac{e_{i}}{e})[/latex]

Longitude of perigee can take any value between 0° and 360° so a half plane check is needed. As can be seen in the figure above, if the jth component of the eccentricity vector is negative, then longitue of perigee is in the second half plane. This can be written mathematically as:

If ej < 0 then 180 < [latex]\Pi[/latex] < 360 degrees.

In the figure above, longitude of perigee is about 90°.

True Longitude

When a satellite is in an equatorial orbit and is circular, it has no ascending node vector, [latex]\vec{n}[/latex], and hence RAAN, [latex]\Omega[/latex], argument of perigee, [latex]\omega[/latex], and true anomaly, [latex]\nu[/latex], are undefined. In this case, we use the true longitude, l, which indicates the location of the satellite in its orbit. Thus, It is the angle measured between the vernal equinox direction, [latex]\hat{I}[/latex] and the position vector, [latex]\vec{R}[/latex]. The true longitude is shown in the figure below:

The image portrays the concept of true longitude in the context of celestial mechanics. True longitude is an orbital parameter that specifies the angular distance of a celestial body along its orbit measured from a reference point, often the vernal equinox. In the figure, the celestial body is situated at the center, and the orbital plane is represented as the horizontal plane. The line extending from the celestial body's center to the point on the orbit indicates the true longitude. This parameter is crucial for determining the celestial body's position in its orbit and is a key element in orbital calculations and predictions.
Figure 21: A True Longitude of Perigee.

Hence, true longitude is given by:

[latex]\cos l = \frac{\hat{I} \times \vec{R}}{R}[/latex]

Or

[latex]l =\cos^{-1} (\frac{R_{i}}{R})[/latex]

True longitude can take any value between 0° and 360° so a half plane check is needed. As can be seen in the figure above, if the jth component of the position vector is negative, then true longitude is in the second half plane. This can be written mathematically as:

If Rj < 0 then 180 < [latex]l[/latex]  < 360 degrees.

In the figure above, true longitude is about 115°.

The Alternate Orbital Elements are summarized in the table below.

Element Name Description Range of Values When To Use
u Argument of Latitude  

[latex]u = \cos^{-1}(\frac{\vec{n} \times \vec{R}}{nR})[/latex]

If Rk<0 then 180 < u < 360 deg

[latex]0 \le u \lt 360[/latex] [latex]e = 0[/latex] and [latex]i \not= 0[/latex]
[latex]\Pi[/latex] Longitude of Perigee  

[latex]\Pi = \cos^{-1}(\frac{e_{i}}{e})[/latex]

 

If ej < 0 then 180 < [latex]\Pi[/latex] < 360 deg

 

[latex]0 \le \Pi \lt 360[/latex]

 

 

 

[latex]i = 0[/latex] and [latex]e \not= 0[/latex]
[latex]l[/latex] True longitude  

[latex]l = \cos^{-1}(\frac{R_{i}}{R})[/latex]

If Rj < 0 then 180 < [latex]l[/latex] < 360 deg

           [latex]0 \le l \lt 360[/latex]

 

[latex]i = 0[/latex] and [latex]e = 0[/latex]

To see if you understand the alternate COEs, complete examples 2 and 3. The answers and solutions are at the end of this chapter.

EXAMPLE 2

[latex]\vec{R} = 10,000 \hat{I} + 0 \hat{J} + 0\widehat{K}[/latex] km

[latex]\vec{V} = 0 \hat{I} + 4.464 \hat{J} - 4.464\widehat{K}[/latex]  km/s

Answer the following questions:

  1.  Where is the satellite currently?
  2. What is the flight path angle at the satellite’s current position?
  3. What is the orbit’s specific angular momentum?
  4. What are the orbit’s orbital elements?

EXAMPLE 3

[latex]\vec{R} = 0 \hat{I} - 7,000 \hat{J} + 0\widehat{K}[/latex] km

[latex]\vec{V} = 9\hat{I} + 0 \hat{J} + 0\widehat{K}[/latex] km/s

Answer the following questions:

  1.  Where is the satellite currently?
  2. What is the flight path angle at the satellite’s current position?
  3. What is the orbit’s specific angular momentum?
  4. What are the orbit’s orbital elements?

Complete this quick quiz in order to see how well you understand the concepts from this section.

In the next section, we will continue considering the “big picture” of how to determine where a satellite is located in orbit given earth-based observations. This will introduce us to two more coordinate systems, PQW and SEZ.

OPTIONAL MATERIAL

ORBITAL VISUALIZATION TOOLS

As part of a first astrodynamics class, the author often uses 3-D orbit visualization tools. A custom-designed and built orbital demonstrate at use at the University of Colorado, Colorado Springs is shown below.

The provided image displays an orbital demonstrator utilized in a first astrodynamics class at the University of Colorado, Colorado Springs. This custom-designed and built demonstrator serves as a valuable educational tool for visualizing three-dimensional orbits. It likely incorporates various elements to represent key astrodynamics concepts, offering students a hands-on experience in understanding the complexities of orbital mechanics. The use of such demonstrators enhances the learning experience by providing a tangible and interactive way to explore celestial motion principles.
Figure 22: Orbital Demonstrator in UCCS Classroom.

However, since it is not possible for every student to have a large orbital demonstrator, we provide a hand-held orbital element model, or “whiz wheel”, to help students with the visualization of orbits. A 3-D printed “whiz wheel” is shown below:

The image depicts a hand-held orbital element model, often referred to as a "whiz wheel," designed to aid students in visualizing orbits. This 3D-printed wheel serves as a portable educational tool, allowing students to interactively explore orbital elements and better understand the dynamics of celestial motion. While not as large as a full-scale orbital demonstrator, the whiz wheel offers a practical and accessible means for students to grasp fundamental astrodynamics concepts. Its design likely includes various components representing key orbital parameters, providing a hands-on learning experience in a more compact form.
Figure 23: A 3-D “Whiz Wheel” Orbital Demonstrator.

This device has the capability of demonstrating the following:
1. The physical characteristics of the COEs including how they are measured relative to the principle axes and fundamental plane of the geocentric-equatorial coordinate system, and in which direction the COEs are measured.
2. The effect that change in the values of the COEs will have on the orbit.
3. The resulting ground tracks caused by the earth rotating on its axis under the orbit fixed in space.

Take this quick quiz to see if you have a good understanding of orbital visualization.

When the author taught at the US Air Force Academy in the Department of Astronautics, cardboard paper versions of whiz wheels were used as an orbital visualization tool. When 3-D printing technology came along, it was a perfect candidate for development of a more durable tool, while keeping it small enough to still be handheld and easily carried in a backpack.

For instructions on how to build your own whiz wheel out of cardboard paper and some examples on how to use it, see the following attachments:

10 New whiz_wheel (ivory print at 74-76 percent)

Whiz Wheel Examples

Email the author if you would like the files for the 3-D version.

 

PROGRAMMING THE ORBITAL ELEMENTS

  1. Given position and velocity vectors for any earth satellite, write a program to calculate the COEs. Your output should be formatted to clearly indicate the orbit type (circular, elliptical, parabolic, hyperbolic, and/or equatorial) as well as the COEs or alternate COEs.

For practical purposes, consider orbits:

  1. circular if eccentricity < .001
  2. parabolic if < .001
  3. equatorial if inclination < .001°
  1. Run your program for the following cases. Include an accompanying set of hand calculations for cases 2 and 5. All units are in kms and seconds.

Case 1

[latex]\vec{R} = -424.0961\hat{I} - 369.963 \hat{J} + 7757.78\widehat{K}[/latex]

[latex]\vec{V} = -1.364721 \hat{I} + 7.9109 \hat{J} + 2.86777\widehat{K}[/latex]

 

Case 2

[latex]\vec{R} = 10000\hat{I} + 0 \hat{J} + 0\widehat{K}[/latex]

[latex]\vec{V} = 0 \hat{I} + 4.464 \hat{J} - 4.464\widehat{K}[/latex]

 

Case 3

[latex]\vec{R} = -12208 \hat{I} - 25698 \hat{J} - 8680\widehat{K}[/latex]

[latex]\vec{V} = 4\hat{I} + 0 \hat{J} - 6\widehat{K}[/latex]

 

Case 4

[latex]\vec{R} = 19455 \hat{I} + 8305 \hat{J} + 0\widehat{K}[/latex]

[latex]\vec{V} = 3\hat{I} + 3 \hat{J} + 0\widehat{K}[/latex]

 

Case 5

[latex]\vec{R} = 24912.16\hat{I} + 0 \hat{J} + 0\widehat{K}[/latex]

[latex]\vec{V} = 0 \hat{I} + 4 \hat{J} + 0\widehat{K}[/latex]

 

Case 6

[latex]\vec{R} = 7199\hat{I} + 9700\hat{J} + 15940\widehat{K}[/latex]

[latex]\vec{V} = 4.464 \hat{I} + 4.464 \hat{J} + 0\widehat{K}[/latex]

 

Solutions to Cases 1 and 5 are provided.

Program Case 1

Program Case 5

 

SOLUTIONS TO EXAMPLES:

***EXAMPLE 1 SOLUTION***

[latex]\vec{R} = 0 \hat{I} + 0 \hat{J} + 10,000\widehat{K}[/latex] km

[latex]\vec{V} = 6 \hat{I} + 0 \hat{J} + 0\widehat{K}[/latex] km/s

Answer the following questions:

  1. Where is the satellite currently?  Over the north pole. You know this because the position vector has only a positive [latex]\widehat{K}[/latex] component.
The image illustrates a scenario where a satellite is currently positioned over the North Pole. This conclusion is drawn from the statement that the position vector has only a positive component. The depiction likely represents a simplified 3D coordinate system where the positive direction aligns with the North Pole, and the absence of negative components in the position vector indicates the satellite's location in the northern hemisphere. This example serves to highlight the relationship between position vectors and the satellite's geographical position in the context of orbital dynamics.
Figure 24: A Positive Vector With a Positive (K) Component.

2.  Flight path angle, [latex]\phi[/latex] = 0. We know this because the position and velocity vectors are perpendicular (their dot product is zero).

3. Specific angular momentum vector, [latex]\vec{h}[/latex]

[latex]\vec{h} = \vec{R} \times \vec{V}[/latex]= [latex]\begin{vmatrix} \hat{I} & \hat{J} & \widehat{K}\\ 0 & 0 & 10,000 \\ 6 & 0 & 0 \end{vmatrix}[/latex] = [latex]60,000 \hat{J}[/latex]km2/s

[latex]\vec{h}[/latex]= 60,000 [latex]\hat{J}[/latex] km2/s

4.  Semi-major axis, a

[latex]\varepsilon = \frac{V^{2}}{2} - \frac{\mu}{R} = \frac{6^{2}}{2} - \frac{389600.5}{10000} = -21.86[/latex] km2/s2

 

[latex]a = -\frac{\mu}{2\varepsilon} = -\frac{398600.5}{2(-21.68)} = 9117[/latex] km

        a = 9117 km

      5.  Eccentricity, e

In this example, we will look at three methods to find e:

a.  Method 1. Use the long form of the eccentricity vector equation; however, in this case it simplified a lot since [latex]\vec{R} \times \vec{V}[/latex].

 

[latex]\vec{e} = \frac{1}{\mu}[(V^{2} - \frac{\mu}{R}) \vec{R} - (\vec{R} \times \vec{V})\vec{V}] = \frac{1}{389600.5}[(6^{2} - \frac{398600.5}{10000})10000\widehat{K} -0] = -0.0968\widehat{K}[/latex]

[latex]e = \vec{e} = 0.0968[/latex]

       e = 0.0968

b.  Method 2. Use the dot product form of the eccentricity vector equation.

 

[latex]\vec{e} = \frac{1}{\mu}(\vec{V} \times \vec{h}) - \frac{\vec{R}}{\left| \vec{R} \right|}[/latex]

 

[latex]\vec{V} \times \hat{h} = \begin{vmatrix} \hat{I} & \hat{J} & \widehat{K}\\ 6 & 0 & 0 \\ 0 & 60000 & 0 \end{vmatrix} = 360,000\widehat{K}[/latex]

 

[latex]\hat{e} = \frac{1}{398600.5}(360,000\widehat{K})- \frac{10000\widehat{K}}{10000} = -0.0968\widehat{K}[/latex]

       e = 0.0968

c.   Method 3. Use radius of apogee to determine the magnitude of e.  This takes a little thought, but once you have done these a few times it becomes more and more obvious.   Follow this logic:

i.  Since [latex]\vec{R} \times \vec{V}[/latex], the two vectors are perpendicular. Hence, the flight path angle will be zero, which only occurs at two points in an elliptical orbit: apogee and perigee.

ii. Since the semi-major axis, a = 9117 km and the position vector is currently at a magnitude of 10,000 km, the satellite must be at apogee as can be seen in the figure below:

This conclusion is drawn from the given parameters, specifically the semi-major axis (a) being 9117 km and the magnitude of the position vector currently at 10,000 km. In an elliptical orbit, apogee represents the point farthest from the central body (in this case, Earth). The figure visually conveys the relationship between the semi-major axis, position vector magnitude, and the satellite's position in the orbit, emphasizing the concept of apogee in orbital dynamics.
Figure 25: A Position Vector With a Magnitude of 10,000 KM.

Since Ra = a(1 + e)

10000 = 9117 (1 + e)

e =  10000/9117 – 1

e = 0.0968

     6. What is its inclination, i?

Since the position vector has only a [latex]\widehat{K}[/latex] component,

i = 90 degrees

Alternately, you may calculate this using the equation:

 

[latex]i = \cos^{-1}(\frac{h_{k}}{h})[/latex]

From 3,

[latex]\vec{h} = 60,000 \hat{J}[/latex] km2/s

The k component of [latex]\vec{h}[/latex] is 0, hence,

i = cos-1(0) = 90 degrees

i = 90 degrees

  7.   What is its RAAN, [latex]\Omega[/latex]?

Since the velocity vector has only an [latex]\hat{I}[/latex] component, by inspection (see Figure “Example 1 Orbit” under #1),

[latex]\Omega[/latex] = 180 degrees

Alternately, you may calculate RAAN using the equation:

[latex]\Omega = \cos^{-1}(\frac{n^{i}}{n})[/latex]

First, you will need to find the ascending node vector =  [latex]\vec{n} = \widehat{K} \times \vec{h}[/latex]

[latex]\vec{K} \times \hat{h} = \begin{vmatrix} \hat{I} & \hat{J} & \widehat{K}\\ 0 & 0 & 1 \\ 0 & 60000 & 0 \end{vmatrix} = - 60,000\widehat{I}[/latex]

 

[latex]\omega[/latex] = cos-1(-60000/60000) = cos-1(-1)

[latex]\omega[/latex] = 180 degrees

8.  What is its argument of perigee, [latex]\omega[/latex]?

Since the ascending node vector, [latex]\vec{n}[/latex], only has a negative[latex]\hat{I}[/latex] component and the eccentricity vector points through the south pole,

[latex]\omega[/latex] = 270 degrees

Alternately, you may calculate argument of perigee using the equation:

 

[latex]\omega = \cos^{-1}(\frac{\vec {n} \times \vec{e}}{ne})[/latex]

since [latex]\vec{n} = -60000 \hat{I}[/latex]

and [latex]\vec{e} = -0.0968 \hat{K}[/latex]

Their dot product is equal to zero.  Hence

[latex]\omega[/latex] = cos-1(0) = 90 degrees or 270 degrees

However, since, ek<0,

[latex]\omega[/latex] = 270 degrees

9.  What is the satellite’s true anomaly, [latex]\nu[/latex]?

Since the satellite is currently at apogee,

[latex]\nu[/latex] = 180 degrees

Alternately, you may calculate true anomaly using the equation:

 

[latex]\nu = \cos^{-1}(\frac{(-0.0968)(10000)}{(0.0968)(10000)})[/latex]

[latex]\nu[/latex] = cos-1(-1)

[latex]\nu[/latex] = 180 degrees

 

***EXAMPLE 2 SOLUTION***

[latex]\vec{R}= 10,000 \hat{I} + 0\hat{J} + 0\widehat{K}[/latex] km

[latex]\vec{V} = 0\hat{I} + 4.464\hat{J} - 4.464\widehat{K}[/latex] km/s [latex]\vec{V} =  0\hat{I} + 4.464\hat{J} - 4.464\widehat{K}[/latex] km/s

1. Where is the satellite currently?    Over the equator at the descending node.

2. What is the flight path angle at the satellite’s current position?  0 (position and velocity vectors are perpendicular)

3. What is the orbit’s specific angular momentum?

[latex]\vec{h} = \vec{R} \times \vec{V} = \begin{vmatrix} \hat{I} & \hat{J} & \widehat{K}\\ 10000 & 0 & 0 \\ 0 & 4.464 & 4.464 \end{vmatrix} = - 44,640\hat{J} + 44,640\widehat{K}[/latex] km2/s

[latex]\vec{h} = 0\hat{I} + 44,640 \hat{J} + 44,640 \widehat{K}[/latex] km2/s

4. What are the orbit’s orbital elements?

Semi-major axis, a

 

[latex]V = \left| \vec{V} \right| = \sqrt{4.464^{2} + 4.464^{2}} \frac{km}{s}[/latex]

V = 6.313 km/s

 

[latex]R = \left| \vec{R} \right| = 10,000[/latex] km

so [latex]\varepsilon = \frac{V^{2}}{2} - \frac{\mu}{R} = \frac{6.313^{2}}{2} - \frac{398600.5}{10000} = -19.93 \frac{km}{s}[/latex]

 

and [latex]a = -\frac{\mu}{2\varepsilon} = -\frac{398600.5}{2 \times(-19.93)}[/latex]

a = 10000 km

Thus, this is a circular orbit.

 

eccentricity, e

Since this is a circular orbit, e = 0

 

inclination, i

By inspection using a drawing or visualization tool (introduced later), it can be deduced that:

i = 45 degrees

Alternately, it can be calculated by [latex]\cos i = \frac{\vec{h} \times \widehat{K}}{hK}[/latex]

From part 3, [latex]h = \left| \vec{h} \right| = \sqrt{44640^{2} + 44640^{2}} \frac{km^{2}}{s}[/latex]

h = 63130 km2/s

So [latex]\cos i = \frac{44640}{63130} = \frac{1}{\sqrt{2}}[/latex]

i = 45 degrees

Right Ascension of Ascending Node (RAAN), [latex]\Omega[/latex]

Again, by inspection it can be determined that:

[latex]\Omega[/latex] = 180 degrees

Alternatively, it can be calculated by [latex]\cos \Omega = \frac{\hat{I} \times \widehat{n}}{In}[/latex]

Where [latex]\vec{n} = \widehat{K} \times \vec{h}[/latex]

[latex]\vec{n} = \begin{vmatrix} \hat{I} & \hat{J} & \widehat{K}\\ 0 & 0 & 1 \\ 0 & 44640 & 44640 \end{vmatrix}[/latex]

[latex]\vec{n} = -44640\hat{I} + 0\hat{J} + 0\widehat{K}[/latex] km2/s

so [latex]\cos \Omega = \frac{-44640}{44640} = -1[/latex]

[latex]\Omega[/latex] = 180 degrees

Now since this is a circular orbit, there is no argument of perigee, so we must use an alternate COE.

 

Argument of latitude, u

u is measured from the ascending node vector, [latex]\vec{n}[/latex] to the position vector, [latex]\vec{R}[/latex]

Again, by a good drawing or visualization tool, it can be seen by inspection that:

u = 180 degrees

Alternatively, it can be calculated by [latex]\cos u = \frac{\vec{n} \times \vec{R}}{nR}[/latex]

From before, we found:

[latex]\vec{n} = -44640\hat{I} + 0\hat{J} + 0\widehat{K}[/latex] km2/s

and since [latex]\vec{R} = 10,000\hat{I} + 0\hat{J} + 0\widehat{K}[/latex] km

so [latex]\cos u = \frac{(-44640)(10000)}{(44640)(10000)}[/latex]

u = 180 degrees

 

***EXAMPLE 3 SOLUTION***

[latex]\vec{R} = 0\hat{I} - 7000\hat{J} + 0\widehat{K}[/latex] km

[latex]\vec{V} = 9\hat{I} + 0\hat{J} + 0\widehat{K}[/latex] km/s

Answer the following questions:

  1.  Where is the satellite currently?

About all you can tell at this point is that the satellite is over the equator at either apogee, perigee, or in a circular orbit.  There are no [latex]\widehat{K}[/latex] components in either the velocity or position vectors, so this must be an equatorial orbit. 

2. What is the flight path angle at the satellite’s current position?  Zero

We observe that the [latex]\vec{R}[/latex] and [latex]\vec{V}[/latex] are perpendicular, so at that instant the flight path angle is zero, hence at apogee, perigee (in an elliptical orbit) or in a circular orbit.

3.  What is the orbit’s specific angular momentum?

[latex]\vec{h} = \vec{R} \times \vec{V}[/latex]

[latex]\vec{h} = \begin{vmatrix} \hat{I} & \hat{J} & \widehat{K}\\ 0 & -7000 & 0 \\ 9 & 0 & 0 \end{vmatrix}[/latex]

[latex]\vec{h} = 63000\widehat{K}[/latex] km2/s

4. What are the orbit’s orbital elements?

Semi-major axis, a

[latex]\varepsilon = \frac{V^{2}}{2} - \frac{\mu}{R} = -\frac{\mu}{2a}[/latex]

so

[latex]\varepsilon = \frac{9^{2}}{2} - \frac{398600.5}{7000} = -16.44[/latex] km2/s2

[latex]a = -\frac{\mu}{2\varepsilon}[/latex] so [latex]a = -\frac{398600.5}{7000}[/latex]

a = 12120 km

Eccentricity, e

In this example, we will look at three methods to find e.

a.  Method 1. Use the long form of the eccentricity vector equation; however, in this case it simplified a lot since [latex]\vec{R} \times \vec{V} = 0[/latex].

[latex]\vec{e} = \frac{1}{\mu}[(V^{2}-\frac{\mu}{R})\vec{R} - (\vec{R} \times \vec{V})\vec{V}] = -.422\hat{J}[/latex]

[latex]e = \left| \vec{e} \right| = 0.422[/latex]

       e = 0.422

 

b.  Method 2. Use the dot product form of the eccentricity vector equation.

 

[latex]\vec{e} = \frac{1}{\mu}(\vec{V} \times \vec{h}) - \frac{\vec{R}}{\left| \vec{R} \right|}[/latex]

where

 

[latex]\vec{V} \times \vec{h} = \begin{vmatrix} \hat{I} & \hat{J} & \widehat{K}\\ 9 & 0 & 0 \\ 0 & 0 & 63000 \end{vmatrix} = -567000\hat{J}[/latex]

so

 

[latex]\vec{e} = \frac{1}{398600.5}(-567000\hat{J}) - \frac{(-7000\hat{J})}{7000}[/latex]

 

[latex]\vec{e} = -1.4225\hat{J} + 1\hat{J}[/latex]

[latex]\vec{e} = 0.4225\hat{J}[/latex]

e = 0.422, which is the same result we had previously.

 

c.   Method 3. Use radius of apogee to determine the magnitude of e.  This takes a little thought, but once you have done these a few times it becomes more and more obvious. Follow this logic:

i.  Since [latex]\vec{R} \times \vec{V} = 0[/latex], the two vectors are perpendicular. Hence, the flight path angle will be zero, which only occurs at two points in an elliptical orbit:  apogee and perigee.

ii. Since the semi-major axis, a = 12120 km and the position vector is currently at a magnitude of 7000 km, the satellite must be at perigee as can be seen in the figure below.

The presented image illustrates a scenario where the satellite is at perigee in its orbit. This determination is based on the provided parameters, specifically the semi-major axis (a) being 12120 km and the magnitude of the position vector currently at 7,000 km. In an elliptical orbit, perigee represents the point closest to the central body (in this case, Earth). The figure visually emphasizes the relationship between the semi-major axis, position vector magnitude, and the satellite's position in the orbit, highlighting the concept of perigee in orbital dynamics.
Figure 26: A Position Vector With a Magnitude of 7,000 KM.

Since Ra = a(1 – e)

7000 = 12120 (1 – e)

e =  1 – 7000/12120

e = 0.422

 

inclination, i

By inspection using a drawing or visualization tool and since there are no [latex]\widehat{K}[/latex] components in either the velocity or position vectors, so this must be an equatorial orbit, it can be deduced that:

i = 0 degrees

Alternately, you could calculate it by using:

[latex]\cos i = \frac{\vec{h} \times \widehat{K}}{hK}\large[/latex]

Since from part 3)  [latex]\vec{h}[/latex] = 63000 [latex]\widehat{K}[/latex] km2/s

[latex]\cos i = \frac{63000}{63000} = 1 \large[/latex]

so i = 0 degrees

 

Longitude of Perigee, [latex]\Pi[/latex]

Since there is no ascending node, there is no right ascension of ascending node [latex]\Omega[/latex]. However, this is an eccentric orbit so it has a perigee and since longitude of perigee is measured from the [latex]\hat{I}[/latex]vector to perigee, by inspection, we can deduce:

[latex]\Pi[/latex] = 270 degrees

Alternatively you could calculate it using:

 

[latex]\cos \Pi = \frac{\hat{I} \times \hat{e}}{e}[/latex]

From above, [latex]\vec{e}[/latex] = -0.4225 [latex]\hat{J}[/latex]

cos[latex]\Pi[/latex] = 0

so [latex]\Pi[/latex] = 90 degrees or 270 degrees

However, our quadrant check, if ej < 0 then 180 < [latex]\Pi[/latex] < 360 deg, tells us [latex]\Pi[/latex] must be in the second half plane, hence,

[latex]\Pi[/latex]270 degrees

True Anomaly, [latex]\nu[/latex]

Since the satellite is currently at perigee,

[latex]\nu[/latex] = 0 degrees

 

Alternately, you may calculate true anomaly using the equation:

[latex]\cos \nu = \frac{\vec{e} \times \vec{R}}{eR}[/latex]

[latex]\cos i = \frac{63000}{63000} = 1[/latex]

[latex]\nu[/latex] = cos-1(1)

[latex]\nu[/latex] = 0 degrees

 

 

REFERENCES

Bate, R. R., Mueller, D. D., White, J. E., & Saylor, W. W. (1971). Fundamentals of astrodynamics. Mineola (N.Y.): Dover Publications.

Exo Cruiser (2016, March).  Where is the Solar System Zero? dodlithr.blogspot.com

Martin, Reeve. (2009, June 28). First point of Aries. https://stargazerslounge.com/topic/44684-first-point-of-aries/

Rune, Signa. (2014, July 18). The two-body equation of motionhttps://www.slideserve.com/signa/the-two-body-equation-of-motion

Sanad, Ibrihim. (2013, January).  (Pdf) Design of Remote Sensing Satellite Orbit. (n.d.).  https://www.researchgate.net/publication/311486304_Design_of_Remote_Sensing_Satellite_Orbit

Strimel, Greg & Bartholomew, Scott & Kim, Eunhye. (2017). Engaging Children in Engineering Design Through the World of Quadcopters. Children’s Technology and Engineering Journal. 21. 7-11.

Sellers, Jerry Jon.  (2005).  Understanding Space:  An Introduction to Astronautics, 3rd edition.  McGraw-Hill.

Siegenthaler, Kenneth, David Barnart, Jerry Sellers, James White and Tim White (2004).  Educational Methods/Aids Used in the
Astronautics Core Course at the United States Air Force Academy.  Proceedings of the 2004 American Society for Engineering
Education Annual Conference & Exposition, American Society for Engineering Education.”

Taghinivarami, Khodayar (2016, July).  (Pdf) Rotational Evolution of Low-Mass Stars in Close Binaries.  https://www.researchgate.net/publication/309727000_Rotational_Evolution_of_Low-Mass_stars_in_close_binaries

United States Air Force Academy Department of Astronautical Engineering. (n.d.)  Astro 310 Whiz Wheel.

United States Air Force Academy Department of Astronautical Engineering (n.d.)  Whiz Wheel Examples.

License

Icon for the Creative Commons Attribution-NonCommercial 4.0 International License

Introduction to Orbital Mechanics Copyright © by lgeorge2 is licensed under a Creative Commons Attribution-NonCommercial 4.0 International License, except where otherwise noted.

Share This Book