2D orbit increasing over time

Question

I have a method for orbiting a sprite (fleet) around another sprite (planet). The problem is that, over time, the fleet moves further away (very slowly, almost unnoticeable unless you let the orbit run for a few minutes).

public void Orbit()
{
    // Planet's position - fleet position
    Vector2 distanceToDestination = currentLocation.Position - position;

    // Calculate fleet's rotation ( - right angle so fleet rotates clockwise)
    rotation = MathHelper.WrapAngle((float)Math.Atan2(distanceToDestination.Y, distanceToDestination.X)) - Helpers.RightAngle;

    // Get the direction that the fleet is facing
    direction = new Vector2((float)Math.Cos(rotation), (float)Math.Sin(rotation));

    // Fleet position, thrust is set to 0.4f
    position = Vector2.Add(position, direction * thrust);

    // Fleet rectangle
    rect = new Rectangle((int)position.X, (int)position.Y, fleetTexture.Width, fleetTexture.Height);
}

public static class Helpers
{
    public const float RightAngle = (float)Math.PI / 2;
}

I would be grateful if the community could point out why my fleet is not maintaining a consistent orbit, which is what I want to achieve! Many thanks in advance.

Thank you both for your excellent answers. @dbc - I selected your answer as I used a variation of the method you posted. Thanks again. — user3256944 salutes Monica, Jul 14 '14 at 07:41

score 2 · Accepted Answer · edited May 23 '17 at 12:13

If I understand what you are doing, you are trying to simulate orbital physics using a model sort of like that of a car driving around a circle (i.e. you're not modeling the acceleration of gravity perpendicular to the current velocity and numerically integrating it). Instead, you're doing this:

Move perpendicular to the center of the circle of radiusr by a small amount d
Reorient so the "car" is still facing perpendicular to the center of the circle.
Repeat.

Now, after each step, how far is your object from the center of the circle? Well, by the Pythagorean theorem, it's sqrt(r*r + d*d), which is going to be slightly larger than r. Using a very large value for d for emphasis:

enter image description here

That should account for your drift.

Incidentally, you don't need to use trig functions to construct a perpendicular to the distanceToDestination. It is much easier to simply cross it with the ±Z vector, the result will be (in 2d) ±(-distanceToDestination.Y, distanceToDestination.X).

Update: since you're not really simulating orbital physics, why not solve the problem analytically?

    public static Vector2 Orbit(Vector2 center, Vector2 startPos, float speed, float time)
    {
        // positive speed means CCW around the center, negative means CW.
        var radiusVec = (startPos - center);
        var radius = radiusVec.Length();
        var angularVelocity = speed / radius; // Add check for divide by zero
        Vector2 perpendicularVec = new Vector2(-radiusVec.Y, radiusVec.X);
        return center + (float)Math.Cos(time * angularVelocity) * radiusVec + (float)Math.Sin(time * angularVelocity) * perpendicularVec;
    }

Blau · Answer 2 · 2014-07-13T18:08:48.943

This is what you are calculating:

enter image description here

To get the next point in the orbit you have to keep the vector lenght.

So your calculationss should be:

// Planet's position - fleet position
Vector2 distanceToDestination = currentLocation.Position - position;

// Calculate fleet's rotation ( - right angle so fleet rotates clockwise)
rotation = MathHelper.WrapAngle((float)Math.Atan2(distanceToDestination.Y, distanceToDestination.X));
rotation -= Helpers.RightAngle;

// Get the direction that the fleet is facing
direction = new Vector2((float)Math.Cos(rotation), (float)Math.Sin(rotation));
direction *= distanceToDestination.Length;

position = currentLocation.Position - direction;

2D orbit increasing over time

2 Answers2