7
$\begingroup$

I recently found the question In a 3D space defined by semimajor axis, eccentricity and inclination, what would be the shape of the space where SDP4 works better than SGP4? and have become interested in trying to systematically explore the idea. However, given that, as explained in @uhoh's question, it would probably be at least a 6-dimensional problem, I was wondering if time (epoch) dimension could be ignored.

The idea would be that, if for any given state vector SV1 at epoch T1 that will result in orbit 1 (let's say the trajectory propagated for 1 period) there is another state vector SV2 at epoch T2 that will result in the same orbit (ignoring that the position of the satellite along the orbit will be different at different propagation times), it might be possible to sample all possible orbits at a given time, instead of also needing to sample along epoch values.

Is this indeed the case? Does it depend on how similar the two resulting trajectories would need to be in order to be considered "the same orbit"?

Intuitively, I'm inclined to think that it might also depend on the applied propagator. For example, if using a high-precision numerical propagator that considers the forces exerted by all other planets in the solar system, it would be almost impossible to find the same situation at two different time points.

$\endgroup$

2 Answers 2

6
$\begingroup$

Partial answer because this is a hard question and because clarifying some of the terms may require more than the exchange of a few comments in this particular case. Kudos for taking on a challenge!

What "same orbit" might mean is hard to pin down

The term "same orbit" is really hard to pin down, different people use it for different levels of "sameness".

Take the A-train for example (well, don't actually Take the A-train, nor this or this one)

enter image description here

Source

These satellites are basically in the "same orbit" in that they will remain fairly close to each other over weeks or longer without (much) station keeping. As the Earth turns beneath them each samples a slightly different part of Earths lumpy gravitational field, but for the most part their travel can be analogous to that of cars in a train, where one goes, so will the others behind it.

Satellite constellations like Planet's Doves, various GNSS systems and all the space-based messaging and internet systems stuff multiple satellites into each of several orbits, sometimes a dozen or more spread out around a full 360°. That's spread out over about 90 minutes, during which the Earth turns by about 22.5° underneath.

We can say satellites share the same orbital plane, but...

We can say they share the same orbital plane, and if someone said they're in the "same orbit" nobody would object. But due to the rotating Earth's lumpy gravity, their orbits won't be "exactly" in the same plane or have "exactly" the same shape or altitudes. Referenced to the center of mass of the Earth in a non-rotating frame, the orbits will differ by several kilometers up to 10 or more km in some cases.

"Same orbit"? It depends on the context of the conversation.

If your time frames T1 vs T2 are days apart or longer, then you have to keep track of where the Sun and Moon are, as their effect is important over long periods of time.

And also remember that inclined orbits slowly precess around Earth. For example the ISS precesses all the way around every 55 or 60 days or so.

About the linked question

Now about the linked question In a 3D space defined by semimajor axis, eccentricity and inclination, what would be the shape of the space where SDP4 works better than SGP4? The difference between "coke classic" SGP41 and SDP4 is that for high orbits SDP4 tries to include the long term effects of the tidal perturbations of the Sun and Moon's gravity on the orbit. It doesn't know exactly where they are at any given moment, it just tries to evolve the orbital elements gradually in a way that these effects might do it.

That said

That said, how can we proceed to discuss what "same orbit" means or doesn't mean in this context?

I think this (and other) answer(s) can evolve as we try to pin this down.


1These days "SGP4" refers to packages that include SDP4 and automatically switch to that propagator when the period of the orbit is longer than 225 minutes.

$\endgroup$
14
  • 1
    $\begingroup$ Thanks a lot for pinning down the main underlying issue, of defining what exactly "same orbit" is. I find a very interesting point that SDP4 does not actually consider the position of Sun/Moon at any point, but instead evolves orbital elements of the target satellite to mimic the effects of the long terms exerted by Sun and Moon. So I guess then, one can indeed sample the possible outcomes of SGP4 and SDP4 without taking into consideration the epoch time (and I guess this makes sense, since SGP4/SDP4 only require a time relative to epoch to propagate, and not the absolute epoch itself). $\endgroup$
    – Rafa
    Oct 18, 2021 at 6:29
  • 1
    $\begingroup$ The problem comes to how to get the reference to compare the outcome of both SGP4 and SDP4 to determine which is more accurate in different regions of the orbital elements multi-dimensional space. For performing a rather exhaustive search, I thought the only possibility would be to also propagate theoretical satellites from the same starting points with a high-precision numerical propagator that considers atmospheric drag, third-body perturbations, solar radiation, etc. However, these do indeed require the time at epoch as an input (to get celestial bodies ephemerides, Earth orientation, etc.) $\endgroup$
    – Rafa
    Oct 18, 2021 at 6:33
  • 1
    $\begingroup$ So I guess one would need to compare the sampled SGP4/SDP4 space of outcomes to the space of outcomes with the accurate propagator at multiple different epochs (and that is without considering what satellite features such as mass, drag and radiation coefficients to choose for propagation... ). It seems like the problem from the original question quickly becomes very highly-dimensional, unless a series of assumptions are made. Maybe sampling along times of the largest period of perturbations accounted for in SDP4 might be a compromise solution? (which I guess might be a lunar month?) $\endgroup$
    – Rafa
    Oct 18, 2021 at 6:38
  • 1
    $\begingroup$ @Rafa quick reply for now: yes I believe that comparing both to numerical integration is the right and only way to answer that question. SGP4/SDP4 are clearly documented but as they are approximations with devils in details I think they are best treated as black boxes. Besides gravitational effects, atmospheric drag is the only other perturbation considered and setting the drag terms in the TLEs to zero (appropriate for the very high orbits where the comparison becomes very interesting) turns those off so it wouldn't have to be numerically simulated. $\endgroup$
    – uhoh
    Oct 18, 2021 at 7:02
  • 1
    $\begingroup$ @Rafa since these are high orbits, Earth's tesseral geopotential terms ("lumpiness" in the longitudinal direction) will be extremely weak so we don't care about Earth's rotation so we don't care about epoch. I'm pretty sure this becomes a 3D parameter space problem depending only on semimajor axis, eccentricity and inclination. Another devil in another detail will be choosing the starting state vector for the numerical part. $\endgroup$
    – uhoh
    Oct 18, 2021 at 7:10
1
$\begingroup$

I feel that I understand what you are trying to solve but there are many « loose ends » in your statement of problem that it is difficult to be affirmative. So, this is not really an answer, but rather a collection of observed « loose ends » and some discussions on how to get rid of them.

  • « orbit » vs « trajectory »

It is not clear whether, for you, they are synonyms. Take this quote from your question

 SV1 at epoch T1 that will result in orbit 1 (let's say the trajectory propagated for 1 period)

This can be interpreted that orbit is a whole, while « trajectory propagated for 1 period » is just a part of that orbit.

Alternatively (and I believe this the right approach), we can say that the term « orbit » conveys the notion of something moving (hence, having a trajectory), but remaining in the vicinity of something else. Therefore, an orbit is a trajectory, but not all trajectories are orbits. I think, for example, we can’t talk about the "Voyagers’ orbits", at least not any more (I am eager to know if somebody can make a rejoinder).

  • Can « orbit » be defined as a closed trajectory ?

An object has a closed trajectory when it returns to the same positions after some time. Do orbits need to be closed trajectories ? Real life tells us there is no such « closed orbits », or at least, they are extremely rare. If the orbital plane changes with time (and it has to because of Earth oblateness, among other things), it is difficult to have systematic closure.

  • What is exactly the « period » of an orbit ?

The same quote above from your question seems to imply that a state vector (position, velocity) at a given epoch determines a closed trajectory in some judicious reference, so that you can define a period as the minimum time needed for an object to come back to its initial position. Hence, the rest of the trajectory is just a repeat of this part, ad-infinitum.

While this notion is correct in a 2-body situation of orbital-mechanic theory, in a judiciously chosen reference frame, it is easy to observe that it breaks down when there are other forces and also when the reference frame is changed.

You may get away with this difficulty by restating that : SV1 at epoch T1 defines a Keplerian orbit, with period P, in an ECI reference. The actual trajectory between T1 and T1+P and the Keplerian orbit are the objects of interest (but distinct).

  • This is my guess of what you are trying to formulate.

Consider an object orbiting a celestial body and its actual trajectory. If we are given exactly a series of state-vectors SV1,SV2,SV3, … at sampling instant T1,T2=T1+P,T3=T2+P, … can we use the Keplerian orbits, which are determined uniquely by SV1, SV2, SV3, … to approximate the actual trajectory, to within a certain tolerance and within a certain finite time duration ? P needs not to be exactly the Kepler period and can be used to control the approximation errors.

$\endgroup$

Your Answer

By clicking “Post Your Answer”, you agree to our terms of service and acknowledge that you have read and understand our privacy policy and code of conduct.

Not the answer you're looking for? Browse other questions tagged or ask your own question.