Calculate Orbit Velocity Vector in Eci Frame with Position Vector

Calculating the velocity vector in the Earth-Centered Inertial (ECI) frame from a position vector is a fundamental task in orbital mechanics. This process involves determining the velocity components based on the position vector and other orbital parameters. The ECI frame is a non-rotating coordinate system centered on Earth, making it ideal for tracking satellites and other space objects.

Introduction

The Earth-Centered Inertial (ECI) frame is a coordinate system that is fixed relative to the stars and centered on Earth. In this frame, the x-axis points towards the vernal equinox, the z-axis points towards the North Pole, and the y-axis completes the right-handed system. Calculating the velocity vector in the ECI frame from a position vector involves using the orbital mechanics principles to determine the velocity components based on the position vector and other orbital parameters.

This calculation is essential for satellite tracking, orbital maneuver planning, and other space-related applications. The velocity vector in the ECI frame provides the necessary information to predict the future position of an object in space.

Formula

The velocity vector in the ECI frame can be calculated using the following formula:

Velocity Vector Formula

v = (r × h) / (μr) + (h × r) / (μr)

Where:

v is the velocity vector in the ECI frame
r is the position vector
h is the specific angular momentum vector
μ is the standard gravitational parameter of Earth

This formula is derived from the vis-viva equation and the conservation of angular momentum in orbital mechanics. The specific angular momentum vector h is calculated as the cross product of the position vector r and the velocity vector v.

Calculation

To calculate the velocity vector in the ECI frame, follow these steps:

Determine the position vector r in the ECI frame.
Calculate the specific angular momentum vector h using the cross product of r and v.
Use the formula v = (r × h) / (μr) + (h × r) / (μr) to determine the velocity vector.
Verify the results using orbital mechanics principles and ensure consistency with the given position vector.

Assumptions

The calculation assumes a two-body system (Earth and the satellite) and neglects perturbations from other celestial bodies and non-spherical gravity effects. The standard gravitational parameter μ is approximately 3.986004418 × 10¹⁴ m³/s² for Earth.

Example

Consider a satellite with a position vector r = [7000, 0, 0] km in the ECI frame. The specific angular momentum vector h is [0, 42164, 0] km²/s. Using the formula:

Example Calculation

v = (r × h) / (μr) + (h × r) / (μr)

Substituting the values:

v = ([7000, 0, 0] × [0, 42164, 0]) / (3.986004418 × 10¹⁴ × 7000) + ([0, 42164, 0] × [7000, 0, 0]) / (3.986004418 × 10¹⁴ × 7000)

Calculating the cross products:

r × h = [0, 0, -295148000]

h × r = [0, 0, 295148000]

Substituting back:

v = [-295148000 / (3.986004418 × 10¹⁴ × 7000), 0, 0] + [295148000 / (3.986004418 × 10¹⁴ × 7000), 0, 0]

Simplifying:

v = [0, 0, 0] km/s

This example demonstrates the calculation of the velocity vector in the ECI frame using the given position vector and specific angular momentum vector. The result shows the velocity components in the ECI frame.

FAQ

What is the Earth-Centered Inertial (ECI) frame?

The Earth-Centered Inertial (ECI) frame is a coordinate system that is fixed relative to the stars and centered on Earth. It is used for tracking satellites and other space objects because it provides a stable reference frame.

How is the velocity vector calculated in the ECI frame?

The velocity vector in the ECI frame is calculated using the formula v = (r × h) / (μr) + (h × r) / (μr), where r is the position vector, h is the specific angular momentum vector, and μ is the standard gravitational parameter of Earth.