1. Introduction
Uranus and its satellites are one of the least explored systems in our solar system. The only mission sent to visit Uranus took place in 1986 with Voyager 2. A mission to explore the Uranus system would be extremely important to understanding the formation and evolution of the solar system.
The main information about the Uranus system comes from observations made from Earth and data sent by space telescopes, but some important questions still remain unanswered. The distances of the ice giant planets are still an enigma, as some models suggest they should be closer to the Sun. Their high obliquity and low luminosity are intriguing. In addition to all these issues, information about Uranus’ satellites would also be of great importance, as they have irregular surfaces and there is evidence of oceanic formation [
1,
2]. According to [
3], NASA is working in partnership with ESA to plan possible missions to the ice giants during the period 2024–2037. These possible missions can help to understand the origin and evolution of this system.
Titania, Uranus’ largest natural satellite, has a surface with flaws that could indicate past or still active seismic activity. In addition, there are indications that its surface has a mixture of water and ice with carbon dioxide [
4,
5,
6]. An orbiter around Titania could investigate, for example, whether this satellite is or once was an ocean world, or even determine recent endogenous geological activity [
2].
The best way to study this natural satellite is to send a probe around it. Therefore, it is necessary to look for the best orbits around Titania. In the literature, we still do not have references to the study of orbits around Titania that include perturbations due to its non-sphericity. Using maps that show the orbital duration, we present in this work a study to find the best conditions to obtain longer lifetimes of a probe around Titania. Here, Titania is at the center of the system and is perturbed by Uranus. Due to the asymmetry caused by Titania’s ellipticity, the coefficient is analyzed, as well as , which generates a vertical symmetry in the gravitational field around the satellite.
Lifetime mapping is a technique used to find the best initial conditions for orbits as a function of their duration. The map provides the best values of the orbital elements for which orbits last the longest. With this, it is possible to extract significant values of the semi-major axis (a), eccentricity (e), inclination (I), argument of periapsis (), longitude of ascending node (), and mean anomaly (M) to analyze the trajectory of a space vehicle.
This technique is used in the work of [
7] to study a hypothetical system, with masses and distances similar to the Earth-Moon system. In this work, the authors looked for the longevity of a probe around a natural satellite with high eccentricity. The study was done through a second-order expansion of the potential of the perturbing planet, assuming that it is in an elliptical orbit. The results indicated that the high variation of the eccentricity and inclination significantly reduce the useful life of the probe. The authors also showed that this high eccentricity variation is responsible for a 7–8% decrease in the probe lifetime.
These maps were used in [
8] to analyze the best orbits for a probe around Callisto. In this study, the authors used the double-averaged method to investigate the importance of hypothetical eccentricity values for the perturber, Jupiter. Callisto’s non-sphericity was considered and an analysis of the effect of its gravity coefficients was presented. Lifetime maps were used to find the best initial conditions for long-life orbits. In addition, the authors highlighted the importance of adopting specific values for
and
in order to increase the useful life of the probe. The results obtained by them showed that for this system, the contributions due to Callisto gravitational coefficients
and
are small.
The double-averaging method was also used to study another Galilean satellite in [
9]. The dynamics of orbits around Europa were analyzed considering the perturbation of Jupiter in an elliptical orbit and the gravity coefficients
,
, and
of Europa. The study used two processes to find long-duration orbits: the double-average and single-average methods. The authors concluded that even with not eliminating the coefficient
, the simple-average process was more realistic and can be used to find orbits that reach 80 days of life. Orbits up to 200 days of duration have also been found using this method, however, the eccentricity increases slightly in this case. Orbits around Europa were studied again in [
10]. In this work, the authors found a set of long-lived frozen orbits with initial conditions given by
km,
,
,
, and
.
Frozen orbits are of great interest for projecting future space missions. These orbits are characterised by keeping the variation of eccentricity, inclination, and argument of periapsis constant. For these orbits to be found, the argument of periapsis must be fixed at or and the short period elements eliminated. These elements are responsible for causing large variations in eccentricity and inclination due to the dependence on the term .
In [
11], frozen orbits around the Moon were investigated considering the Earth perturbation and the effects of the coefficients
and
. The authors performed a closed-form transformation between osculating elements and averaged elements in order to find initial osculating elements for frozen orbits. This transformation was done only for circular orbits. The results showed that it was necessary to transform the average elements into osculating elements in order to add the
coefficient, as this term caused large variations in high-altitude orbits and can totally change the behaviour of the orbits.
Orbits that keep their eccentricity, inclination, and argument of periapsis constant over a period resemble frozen orbits. In [
12] three families of periodic lunar orbits were found. The results showed that these orbits are very sensitive to the initial conditions and that these conditions have small differences with respect to the frozen orbits. Families of frozen orbits around the Moon were also found in the work of [
13]. The authors eliminated the short period terms and considered the oblateness of the Moon to order 4 (
). For eccentric frozen orbits, all families had a critical inclination around
. In the case of very low quasi-circular lunar orbits, several families with different critical inclinations were obtained, although they quickly disappeared due to small changes in the semi-major axis and eccentricity.
Despite causing large variations in some orbital elements, such as inclination, eccentricity, and argument of periapsis, the term
has an important role in the dynamics of orbits around natural satellites. The combination of this term with the zonal term
is able to attenuate the effects caused by the perturbation due to the third body acting as a “protection mechanism” to prolong/increase the lifetime of the desired orbit [
14,
15]. In [
16] the combination of this coefficient with the
term was also investigated in orbits around the Moon. The results showed that the values of
and
strongly affected the near-critical inclination and the libration amplitude of the argument of periapsis. This work also showed that these effects were greater when added to the contribution of the rotation rate.
In several works found in the literature, the perturbation due to the third body was considered when it comes to orbits around natural satellites. The effects of this perturbation were mainly observed in the eccentricity of the orbit around the observed satellite. These high eccentricity values achieved by the probe’s orbit are related to the effect of the Kozai–Lidov mechanism. This mechanism was studied for the first time in [
17,
18]. This phenomenon causes the argument of periapsis to oscillate at a constant value. Whereas the argument of periapsis oscillates around a specific value, the eccentricity and inclination suffer periodic oscillations (e.g., see [
19,
20] and the references therein).
The structure of this work is organised as follows: in
Section 2 the mathematical model adopted to analyze the system is described, as well as the methodology. In
Section 3 we discuss the results obtained. In
Section 4 we show the isolated effect of each perturbation through the analysis of the orbital elements of each specific case. In this same section, we also present the mechanism for obtaining long-time orbits through the evolution of the orbital elements. In
Section 5 we present an estimate of the probe’s lifetime when affected by possible errors in the values of Titania’s gravitational coefficients. Finally, in
Section 6, we present our conclusions.
2. Mathematical Model
The system addressed in this work consists of a central body (Titania), a space probe, and a perturber (Uranus). We assume that the perturber body is in a Keplerian orbit around Titania with a radius of
km and mass
kg, (
https://ssd.jpl.nasa.gov/), accessed on 5 April 2021. The orbital parameters are given in
Table 1. To make an accurate investigation of the problem, we consider the main gravity coefficients (
and
) of Titania, regarding the oblateness and ellipticity of the body. These terms are the most important in the gravitational field, after the Keplerian term. Previous work [
11,
16,
21,
22] shows that these coefficients are the most relevant to the pertubative effects caused on a probe by a satellite and, in most cases, the most significant to model the irregular shape of a body. It is important to note that these two coefficients are the only coefficients available in the literature. As our focus is to analyze orbits very close to the surface of this satellite, these will be the main terms to be considered in this work. Titania has a mass of
kg and a radius of
km, (
https://ssd.jpl.nasa.gov/), accessed on 5 April 2021. The other parameters of Titania are shown in
Table 2. The equations of motion are described according to [
22]:
where
and
are the expansions of the gravitational potential up to second order for Titania and up to fourth order for Uranus, respectively;
,
, and
are components of the
vector and
,
, and
the components of the
vector;
m is the mass of the space probe,
the mass of Titania, and
the mass of Uranus. The terms
and
are the radius vector of the probe and Uranus, respectively.
The investigation of orbits around natural satellites with a non-uniform mass distribution has been the subject of several studies in recent years. Studies [
8,
9,
14,
22,
23,
24] emphasise the importance of considering the oblateness and ellipticity of the body to be orbited by the probe. However, it is understood that analyzing the influence of higher oblateness coefficients of the disturbing body can provide a much more accurate study of the orbits.
Uranus is very far from Earth, approximately 3 billion km, and the time to send signals from Earth is on the order of 2.8 h; therefore it takes approximately 5.6 h to send a command and receive a signal back. Because of this time lag, orbital maneuvering is a logistical problem. An important reason to search for long-term natural orbits in which to place a space probe is because they do not require frequent orbital maneuvers, thereby saving fuel and simplifying the orbit control. For the sake of comparison, we did a single test considering the highest order gravity coefficients
and
of the disturbing body. This implementation was made through an adaptation of the Mercury package [
25], which was used in our numerical simulations.
Table 1.
Parameters of Titania with respect to Uranus.
Table 1.
Parameters of Titania with respect to Uranus.
Parameter | Value |
---|
Semi-major axis (km) | |
Eccentricity | |
Inclination () | |
Argument of periapsis () | |
Longitude of ascending node () | |
Mean anomaly () | |
Table 2.
Gravity coefficients of Titania and Uranus.
Table 2.
Gravity coefficients of Titania and Uranus.
Body | | | |
---|
Uranus | | | - |
Titania | | - | |
Our main goal is to study circular orbits with high inclinations. Circular orbits are better for observations as eccentric orbits have a lower perigee and are more susceptible to collisions. Highly inclined orbits take advantage of the central body’s rotation to observe the entire natural satellite. However, we understand that analyzing other values of eccentricities can contribute to a more complete approach to the problem, giving more choices for the mission designers. Therefore, we considered the following eccentricity values for the probe’s orbit: , , , , and . In order to analyze those orbits very close to the surface of Titania, we adopted a semi-major axis interval (a) ranging from 810 km to 2500 km for each value of the eccentricity. The inclination is in the range between and .
We considered two cases regarding the argument of periapsis and the longitude of the ascending node ( and , respectively). First, we consider and, after that, we choose specific values in regions with longer lifetimes, to explore best results.
All cases were simulated for 1000 days, equivalent to approximately orbital periods for the closest orbit to Titania’s surface and revolutions for the farthest orbit. For our purposes, this time is adequate, sufficient to analyze the orbital evolution of a space probe and collect and send data for analysis. The lifetime of the orbit is defined by the time that the space probe remains in orbit around Titania without colliding with the central body. Therefore, the maps of versus show the lifetime of the probes numerically simulated up to 1000 days. Our main goal is looking for longer lived orbits around Titania under the effects of its gravitational coefficients and Uranus.
In order to analyze the influence of the gravity coefficients of Titania and the third body, we constructed life maps considering two cases. Point mass models are considered without taking into account the gravity coefficients of both bodies, with gravity coefficients for both bodies, Titania and Uranus, and including
and
of Titania. “Difference maps” (third body—
,
from Titania) are also presented. The case in which we analyzed the
and
coefficients showed insignificant changes in lifetimes compared to the case where these coefficients are not included. Therefore, we will adopt Uranus as a point of mass. Our purpose with these maps is to find regions where the orbits are affected by the gravity coefficients of Titania and the third body. These effects are in opposite directions, and the orbits can last longer without collision or escape. The time evolution of the orbital elements
a,
e,
I,
, and
of the orbits with a long lifetime are shown
Section 4.
3. Analyzing the Results
The results show diagrams
considering: (i) only the contribution of the disturbing body (Uranus), modelled as a point of mass; (ii) the gravity coefficients of Titania (
and
), and (iii) the more complete model, where both effects are considered. In the figures, the color bar represents the lifetimes of the orbits. We constructed “difference maps” by subtracting the lifetime of the orbits obtained when Uranus is considered. The lifetime of the orbits when considering both effects, the gravitational effects of Uranus, and the gravity coefficients of Titania, are presented in
Section 3.
It should be noted that, in these figures showing the differences in the lifetimes, a positive sign indicates those values where only the presence of Uranus prolonged the probe’s life, whereas negative values indicate that the gravity coefficients of Titania helps to naturally control the third-body perturbations Uranus, thus increasing the lifetime of the orbits.
Figure 1 presents a map where we considered only the third body (a) and the third-body and the gravity coefficients of Titania (b). For the case where only the gravity coefficients of Titania are considered,
Figure 1b shows an increase of 50 days in the life of the probe for orbits closer to the surface of Titania. Orbits with lifetimes up to 450 days are observed in the same altitudes and inclinations. Lifetimes in the range 100–300 days for all values of inclinations occur for the semi-major axis in the range 1000–2200 km (
Figure 1a), and
a = 1000–1800 km (
Figure 1b) (shown in light green, green, and blue squares).
We show lifetime “difference maps” in
Figure 1c. The maximum value of the semi-major axis is
km; after that, there are no significant changes in the orbits’ lifetime because the space probe is too far from Titania. Its irregular shape is unable to balance the third-body perturbations from Uranus.
Figure 1c shows that the largest differences in lifetime appear for values of
a = 810–1000 km and
I = 75–83° (blue region), where the equilibrium between the gravity coefficients of Titania and Uranus occurs. We found increases in the lifetime on the order of 60 days, which is not negligible. In addition, those orbits are closer to Titania, which helps in the observation of this body. It is noted that the most negative numbers occur for lower values of inclinations and semi-major axes. In these regions, the third body perturbations are smaller, and the gravity coefficients dominate the dynamics and reduce the duration of the orbits.
For
, the results present similar behaviour compared to
. The simulations were performed using the same initial conditions used in
Figure 1, except for the eccentricity, which now equals
.
We found an equal maximum lifetime when the third body and the and coefficients of Titania are considered. However, when considering only Uranus as a perturber, the region with the longest orbit life ranges in the interval a = 810–900 km and I = 75–80°. In this region, the longest lifetime is approximately 450 days. When including the gravity coefficients of Titania, this region increases to a = 810–1000 km and I = 75–84°. Another fact is that the number of orbits with a lifetime between 300–350 days decreases when the gravity coefficients of Titania are considered.
The results showed that the orbits in the range a = 810–1050 km and I = 75–84° have increased by 60 and 20 days. Orbits with a = 810–900 km and I = 85–90° have their lifetime extended by approximately 40–80 days, thanks to the effect caused by the equilibrium of the perturbations.
For
, shown in
Figure 2, our results, given in the diagram
when only the gravitational effect of Uranus is considered, are similar to those observed for
and
. As already noted, the terms
and
of Titania are more important to increasing the lifetime of the probe, confirming that the effects caused by the gravity coefficients of Titania compensate for the effects caused by Uranus’ gravitational attraction, as expected.
The lifetime differences maps presented in
Figure 2c have two nearby regions where the changes in the lifetimes are expressive (regions in purple, red, and orange squares). In the region with purple squares, the lifetime for those orbits only considering the third body is up to 250 days. At the red dots, the probe’s lifetime can reach up to 200 days when
coefficients of Titania are considered.
Figure 2b, including the gravity coefficients of Titania, presents a more extensive region with lifetimes of 350–500 days. This region is located very close to the surface of Titania (
a = 810–1000 km), where the terms
and
cancel the effects of the third body.
Figure 2a shows the numerical simulations where the gravitational effects of Uranus are considered. Orbits with a lifetime close to these values appear only on a small island, for all inclinations, and a semi-major axis closer to 1400 km.
In the regions where the probe remains around Titania for longer times, there is a balance between the perturbation caused by the third body and the perturbation due to the gravitational coefficients of Titania, causing the probe’s lifetime to be extended. This balance is described in previous work [
9,
14], where it is shown that, when combined with
and
of the central body, this resultant perturbation acts as a “protection mechanism” responsible for softening the effects caused by the third body perturbation on the variation of the eccentricity.
The results presented in
Figure 3 are obtained for
. Note that there are no large differences between the results obtained when considering only the third body perturbation (a) and those when the
and
from Titania are added in the system (b). There is only a decrease in the region with lifetimes between 160 and 180 days, covering all the inclinations, when
and
of Titania are included. This region is located in the range
a = 900–1050 km and
I = 75–80°.
For
, the orbital duration maps show that the results are very close when all the perturbations are considered, and the behaviour is similar to the one obtained for
. The difference is that for
, the probe’s maximum lifetime is much shorter, 60 days, but the regions of greater and lesser duration are very similar. Previous studies [
7,
28] found similar results, where more eccentric orbits have shorter lifetimes, as expected, as these orbits have smaller periapsis and, therefore, collide with the central body in less time. In addition, they have larger apoapsis, which increases the effect of the third body.
Previous work [
8,
29,
30] shows the relevance of the contribution of the argument of periapsis (
) and longitude of the ascending node (
) in the duration of highly inclined orbits. In [
8], the initial conditions capable of increasing the lifetime of the probe are
and
. In [
29], the best values found for these angles are
= 140.5–148°,
= 321–327°, and
. Given the above, we analyzed the regions with longer lifetimes, shown in
Figure 1,
Figure 2 and
Figure 3, in a range of values of the periapsis argument and the node longitude. Within these regions, we chose the best values of
a and
I and built maps as a function of
and
with the four values of eccentricities adopted earlier. For this analysis, we consider the perturbation due to the third body and also the gravity coefficients of Titania.
The first case to explore is for
. The results provided by the simulations showed that the longest orbits occurred for
a = 900–1058 km and
I = 70–80°. The result of this analysis can be seen in
Figure 4.
We fixed a semi-major axis value, 900 km, and varied the inclination to
(
Figure 4a) and
(
Figure 4b). We then fixed the inclination at
and vary the semi-major axis to 1018 km, (
Figure 4c) and 1058 km (
Figure 4d). The angles
and
range from
to
.
Figure 4a,b show similar regions for longer life orbits,
= 0– 150°,
= 220–360° and
= 0–100°,
= 200–300°. In
Figure 4a, where the considered inclination is
, the longer duration orbits have a lifetime of up to 900 days. As we increase the inclination by
, the lifetime increases by 100 days (
Figure 4b). This is due to the action of the gravitational coefficients
and
, which are responsible for attenuating the effects caused by the third body for more inclined orbits.
The regions with longer duration orbits present in
Figure 4c,d, have values of
= 25–150°,
= 50–100°, and
= 150–200°. Orbits with a semi-major axis of 1058 km have a maximum lifetime of up to 650 days, whereas, in
Figure 4d, this time is 800 days for
km.
We chose the initial conditions of a specific point within these longer life regions to redo the simulations for
. The chosen point is highlighted in
Figure 4d (black circle), and has values of
,
,
km, and
.
Figure 4a,b present orbits with longer lifetimes compared to the point chosen in
Figure 4d. However, our main objective is to show that the adoption of these angles can increase the lifetime of the orbit even for regions with shorter orbits. In order to analyze the importance of these angles in the lifetimes, we made a new simulation with values of
,
for
.
The insertion of the values and in the new simulation for increase the probe’s life by approximately 78%, a very large value. We extended the numerical integration for the two cases analyzed: third body and and from Titania. However, the new results show distinct regions with more extended life for each perturbation considered. In the case of the gravitational effects of Uranus, a small island of orbits with a lifetime of 600–800 days appears with a semi-major axis ranging from km to 900 km and inclinations from to . If we include the gravity coefficients of Titania, the maximum lifetime is approximately 600 days. It is located in the region a = 810–1000 km and can be seen for all values of I, increasing from I = 78° to 90°. Orbits with a lifetime of 800 days are found for a value of a close to 900 km and I = 78–90°.
The differences indicate that after the inclusion of the angles and , the gravity coefficients of Titania are more relevant for inclined orbits, I = 80–90°.
For a = 875–950 km, the orbits reach lifetimes between 200 and 400 days. This lifetime is also experienced for orbits with a semi-major axis ranging from 950 km to 1050 km and I = 75–80°, shown in the smaller island. From 810 km to 1100 km, the gravity coefficients of Titania can still extend the duration of the orbits up to 100 days for inclinations in the range of 80–90°. Before the inclusion of and , the orbits that were most affected by the gravity coefficients of Titania had a maximum lifetime up to 60 days. Considering these specific values of those angles, this time has been extended to 400 days, a very large difference.
In the case of Uranus, the differences point out that, when and angles are not zero, there is a deviation in the location of the orbits where their contribution was most significant. When , the most inclined orbits were the most affected by the third body. After that there are orbits with inclinations between 75 and 76° with lifetimes of 100–300 days.
For the case
, the best initial conditions presented in
Figure 2 point to a semi-major axis in the range of 834–1265 km and an inclination of 70–80°. In
Figure 5 we present a diagram
for some specific cases. In
Figure 5a,
km and
, whereas in
Figure 5b, the initial conditions are
km and
.
Figure 5a shows large islands with orbits that reach 450 days. These islands are located at
= 100–200°,
= 300–360°,
, and
= 250–360°. In the case of
Figure 5b, where the inclination is
, the “protection mechanism” caused by the combination of the terms
and
increases the orbital duration to 550 days. However, these islands are much smaller and can be seen at
= 125–180°,
and
and
.
Figure 5c,d have the same inclination
, but different semi-major axes. In
Figure 5c, with
km, the maximum lifetime is almost 3 years. However, the regions are smaller compared to those seen in
Figure 5a,b. They are found with
around
and
close to
. More distant orbits, with a semi-major axis equal to 1000 km, (
Figure 5d), it is noted that the lifetime decreases, and the values of
and
for orbits with the longest duration are
and
, respectively.
An initial condition (a point) with values of
and
is highlighted (black circle) in
Figure 5b. The values of these angles were added in the new simulations for the initial conditions of
Figure 2, in order to investigate how they would affect the lifetime of the probe. The new investigation was carried out by considering only the perturbation of the third body and with the complete system, third body plus the
,
coefficients of Titania. Shown in
Figure 6, it is observed that the use of the non-zero values for these angles in
increases by 100% of the lifetime of the probe in both cases. When both effects are included, the gravitational attraction of Uranus and the gravity coefficients of Titania, the use of non-zero values for
and
also changed those regions where the lifetimes are longer. Orbits with a lifetime up to 600 days are found for inclinations ranging from
to
and
km to
km (third-body effects).
When considering the gravity coefficients of Titania, orbits with this lifetime appear for an interval of semi-major axes equal to a = 810–1050 km and I in the range 75–90° (light green squares). The location of the orbits with a lifetime of approximately 400 days is practically the same for both cases, a = 810–1000 km and I = 75–90° (green squares). The blue dot region has lifetimes in the range 200–300 days for all inclinations and with a in the range 1000–1600 km.
Figure 6c shows that Uranus affects the lifetime of the more inclined orbits. This effect appears in a small orange island that goes from
I = 84–90° for
a = 810–1300 km. In this region, the duration of the orbits is approximately 400 days, considering the lifetime of the probe affected by the gravity coefficients of Titania and Uranus. The effect of
and
of Titania and the third-body acts on orbits with inclinations of 75–90° and semi-major axes of 810–1000 km. The lifetimes of these orbits range from 100 to 600 days compared to the case where only the effect of Uranus is considered; there are no significant changes for
a larger than 1400 km.
For
, the longest-lived orbits are located in the range 900–1160 km and
I = 70–80°. Therefore, the initial conditions used to build the
graphs are within this range of values. In
Figure 7 we present four scenarios for analysis of these angles. In
Figure 7a, the semi-major axis considered is 1000 km and the inclination is
, whereas in
Figure 7b,
a also assumed the value 1000 km and
.
There is a similarity between the regions of orbits with longer orbital duration shown in
Figure 7a,b. In both cases, the islands with the longest lifetime are located for values of
close to
, whereas for
the best values are around
and
. Despite the high value of the eccentricity, with the combination of Titania’s gravity terms with the third body perturbation, long-lived orbits can be found. In
Figure 7a the maximum lifetime is 550 days, whereas in
Figure 7b, where the inclination is greater, the duration of the orbits is up to 650 days.
The best orbits for
(
Figure 7c,d), are for
and for
close to
and
. Due to the proximity to the surface of Titania, the lifetime for the cases analyzed in
Figure 7d is longer, 700 days.
New numerical simulations for
were made according to the initial condition highlighted in
Figure 7c (black circle), where the orbital elements are equal to
a = 1162 km,
,
, and
. These values for the argument of periapsis and the longitude of the ascending node were used in the new simulations with the same values of
a,
e, and
I used in
Figure 3. The results of the inclusion of these new values of
and
are shown in
Figure 8.
When considering and non-zero, the increase of the lifetime of the probe is 2.5 times greater in comparison with . The lifetimes increased from 180 to 450 days.
In
Figure 8a, the orbits closer to the surface of Titania have lifetimes of 250–300 days with a semi-major axis in the range 810–1000 km. It occurs for all inclinations when we considered the perturbation of Uranus. When combining the effects of Uranus and the gravity coefficients of Titania,
Figure 8b, orbits with this lifetime occupy a slightly larger region,
a = 810–1350 km, also for all analyzed values of
I. The islands with lifetimes up to 450 days are seen in both cases: when only the perturbation of Uranus is considered, and when the gravity coefficients of Titania are also taken into account. However, these orbits change locations according to each perturbation investigated. In
Figure 8a, these orbits can be seen from
I = 75–90° and
a = 1300–1500 km, whereas in
Figure 8b, this lifetime is observed for orbits in the range
a = 1100–1400 km and
I =75–76° (orange and red dots).
The orbits with
a = 900–1450 km and
I =75–82° have longer lifetimes when subjected to the effect of the gravity coefficients of Titania. The difference in the lifetimes is up to 50 days, in comparison with the situation where only the influence of the third body is considered. These orbits can be seen in
Figure 8c (green dots).
Numerical integration for showed that the best values for a and I would be 900 km and , respectively. We use these values and plot the lifetimes as a function of varying these angles over the range 0–360°. The results of this analysis are not shown here in this work, as they are similar to the case . However, it is important to note that the adoption of specific values for and also increased the orbit lifetime.
The results shown in
Figure 6b point out that the angles
and
contribute to the increase in the lifetimes for orbits with high eccentricities, as in the case
. With the values found for these angles, the duration of the orbits around Titania increased by two times when considering only the third body perturbation and more than eight times when the gravity coefficients of Titania were also included. In the analysis involving the contribution from Uranus gravitational attraction, the regions of greater and smaller duration remained similar, increasing only the maximum lifetime.
However, exploring the scenario where and of Titania were added to the model, the orbits with zero and non-zero values of and had very different maximum lifetimes as well as larger regions with longer lifetimes. Without the inclusion of specific values for these angles, the orbits with the longest lifetime were located at I = 75–90° and a = 900–1200 km. In this interval, the maximum times reached by these orbits were 60 days. With the adoption of better values for and in the new simulations, these orbits started to have lifetimes from 200 to 500 days, for a = 900–1000 km and I = 75–90°.
The analysis carried out in this section shows the importance of a different approach, first including variations in and , and then mapping for the best values of and . In all cases investigated, the best values for can be found in two regions, in the ranges 100–150° and 27–360°. In the case of the ascending node longitude, the best values are close to 50–100° and 250–360°. The orbits that had their lifetimes extended were those orbits closer to Titania’s surface and with the greatest inclination. The reason for this is associated with the fact that lower and more inclined orbits are subject to the action of the zonal oblateness coefficient , which, added to the term and the perturbation caused by the third body, promotes an equilibrium capable of increasing the duration of the orbit.
4. Long-Duration Orbits
Previous studies of systems similar to the one presented here often use the double averaging method to find frozen orbits. According to the work of [
9], in order to find a frozen orbit, the relations
,
,
must be satisfied. For such relationships to be met, it is common in the literature to define
or
, because, according to the Kozai–Lidov mechanism, under the presence of the third body, the argument of periapsis librate around these values as the eccentricity and inclination oscillate [
17,
18]. However, when these values are fixed,
is non-zero due to the term
, which causes the inclination to increase significantly and hence the eccentricity. In other work, where the double mean model is used, it is possible to eliminate the short period terms that contain the
coefficient and then obtain the necessary conditions to find the frozen orbits.
In this work, we are not interested in finding frozen orbits, whose noted derivatives are equal to zero. Despite showing similarity to a frozen orbit, orbits that last longer are our main interest, independently of the value of these derivatives. This is because our goal is to analyze more realistic models and investigate the isolated effect of the term on the desired orbits.
In this section, we present a study on the isolated effect of each perturbation on the evolution of some orbital elements of a probe for some specific cases. The perturbations considered here are: the gravitational effect of Uranus, the gravity coefficient
from Titania, the ellipticity of Titania (
), Uranus +
from Titania, Uranus +
of Titania, and Uranus +
+
of Titania. This analysis is presented in
Figure 9 and
Figure 10, where the temporal evolution of the orbital elements
e,
I, and
are shown. In all cases analyzed, the numerical integration stops when the probe collides with the surface of Titania.
Figure 9a shows the temporal evolution of the eccentricity under the effect of these perturbations. It can be observed that the effect of the
and
gravity coefficients of Titania are small compared to the effect caused by the third body. This is seen even when the other perturbations are added. When
and
are included in the system, the lifetime increases, as already noted; however, the eccentricity continues to increase due to the presence of the third body.
In the circular restricted three body problem, the Kozai-Lidov effect is common in some cases. This phenomenon is directly related to the presence of the third body, which, as shown in
Figure 9, causes the eccentricity of the probe to reach high values in a short period of time. The Kozai–Lidov mechanism is still capable of making the argument of periapsis librate around a constant value. These values are usually
or
. However, in cases where the eccentricity grows rapidly, as is the case of the present work, a region of libration is not found for the diagram
(e.g., see [
9]). When
and
of Titania are included, there is a balance between the effects caused by the gravity coefficients of Titania and those caused by the Kozai–Lidov effect. This balance causes a slow increase in the eccentricity and, consequently, the lifetime of the orbit increases.
Regarding the evolution of the inclination, shown in
Figure 9b, the third body and the
coefficient do not cause large variations in the probe inclination. The variation caused by
is small and does not exceed
, although the lifetime is 1000 days. When Uranus is considered, the variation is small and the lifetime is just 300 days. When only Uranus and
are taken into account, it is noted that the presence of the third body reduces lifetime compared to the case when only the
coefficient was considered.
The effect responsible for causing large variations in the inclination is the tesseral coefficient
. When we analyze the dynamics only assuming
, the inclination increases from
to
degrees and the lifetime is 1000 days. In all cases in which
is considered, the variation of
I is large. This can be seen in the cases Uranus
and Uranus
, in which the variation of
I can reach
. This effect caused by
is expected, as shown in previous work [
9,
15,
16,
21].
In the analysis of the argument of periapsis, presented in
Figure 10, it is possible to note that, as well as in the evolution of the eccentricity shown in
Figure 9a, the third body is responsible for causing a greater perturbation in this element. The behaviour of
is similar in all cases where Uranus is added. The argument undergoes a great variation in the first days, with oscillations ranging from
to
. After a certain time, the argument of periapsis begins to librate around ∼50° for a short period of time. An expected behaviour for
would be its libration around
or
, as described by Kozai–Lidov [
17,
18]. However, as the short period terms that contain the
coefficient were not eliminated, the inclination grows very quickly and, as the eccentricity depends on the slope, it also reaches high values in a few days.
When only the effects of of Titania are assumed, the argument of periapsis reaches and then decays to approximately until a collision occurs in 1000 days. When only the perturbation due to is considered, the lifetime of the argument of periapsis is also 1000 days; however, the argument of periapsis increases to approximately before the collision.
The relevance of the
and
angles have already been evident in the lifetime maps shown in
Figure 1,
Figure 2,
Figure 3,
Figure 4,
Figure 5,
Figure 6,
Figure 7 and
Figure 8. In
Figure 11, an analysis of the influence of these angles on the orbital elements of two orbits is presented. The orbital elements analyzed in
Figure 11 have different values of
. The first adopted value of
and the other values of the orbital elements were taken from the point highlighted in
Figure 7 (black circle).
The second orbit has the same initial conditions as the highlighted initial condition, except . The value of for the highlighted point is . The use of this value is able to significantly increase the lifetimes of the orbits. For the sake of comparison, a value of just below the circled point was chosen, .
In all cases shown in
Figure 11, the value of
causes an increase in the lifetimes of the probe. The variation of the semi-major axis, shown in
Figure 11a, is very small in both cases. For
and
, the variations are no more than 1 km. For the eccentricity, shown in
Figure 11b, the value of
attenuates its rapid growth. It remains almost constant for 300 days and only then begins to grow, reaching a maximum of 0.3 in 430 days. Using
, the eccentricity increases much faster, and in just 135 days, the collision with the surface of Titania occurs.
The evolution of the inclination, shown in
Figure 11c, presents a variation of
when
, but has a long lifetime, 430 days. Using
, the variation is only
, but the collision happens in a few days.
In the case of the evolution of the argument of periapsis,
Figure 11d, the behaviour is similar to the one shown in
Figure 11 in the cases where the third body was considered. For both
and
, the argument of periapsis oscillates around a constant value for some time. When we consider
, the argument decays and librates around
for approximately 80 days. For
, the argument of periapsis circulates for about 200 days and then proceeds to librate around
for almost 230 days. This behaviour is expected due to the aforementioned Kozai–Lidov mechanism.
In
Figure 11e, the evolution of the ascending node is presented for
and
. For
, the ascending node decays almost
in 140 days, whereas for
this decay is almost
during the first 250 days of integration. After that, the node starts to describe a periodic behaviour until the collision occurs.
In order to gather all the best orbits investigated, a summary of the best lifetimes for each initial condition and for each case analyzed is presented in
Table 3. This table presents the best initial conditions for the five values of the eccentricity studied in this work. In addition to the results for
, the results considering
and
different from
are also presented. The table shows the best initial conditions obtained for long-duration orbits, with these conditions given by:
km,
,
,
, and
considering only the third body perturbation. For the system, third-body
of Titania, we have
km,
,
,
, and
. In both cases, the orbits have a lifetime of 1000 days.
In
Figure 12 we present different projections of one of the orbits with the best initial conditions found in this work. The trajectory was numerically simulated for 991 days, considering the complete system, the gravitational effects of the third body (Uranus), and Titania’s gravitational coefficients (
and
). The initial conditions are:
km,
,
,
, and
.
6. Final Comments
In this work we study the dynamics of orbits around the natural satellite Titania, the largest moon in the Uranus system. The gravitational attraction due to Uranus is considered, as well as the gravity coefficients and of Titania. Titania is at the center of the system, and Uranus is in an elliptical orbit and acts as the third body to perturb the motion of the space probe around Titania.
Through a set of numerical simulatons of the equations of motion, we build lifetime maps for different values of eccentricities of the orbit of the space probe. The construction of these maps allowed us to find values of a, e, I, , and capable of increasing the lifetime of the probe. We also constructed lifetime “difference maps” in order to analyze the regions where the effect of the third body and the gravity coefficients of Titania were in equilibrium, allowing longer life for the space probe.
We also present an analysis of the importance of the angles and in the lifetime of the probe. In addition, we studied the response of each effect on the orbital element.
Our results showed that, for the altitudes adopted in this work, the oblateness coefficients and of Uranus do not affect the orbits of the space probe. We also showed that the orbits with a longer lifetime had a semi-major axis close to 900 km, whereas the best eccentricity is between and when considering . The results also point out that the inclusion of specific values of and are extremely important in increasing the lifetime of the probe. In all cases, the increase in the lifetime was more than 50% after adopting specific values for these two angles. The best values for are around 145–165° and for around 55–205°. In some cases, including these values resulted in the probe lifetime increasing up to eight times compared to cases with .
Our main objective is to find long-duration orbits around Titania. After analyzing the results including the values of and , we found ideal initial conditions for an orbit lasting up to 1000 days. This orbit has initial orbital elements given by: km, , , , and .
Analyzing the effect of each perturbation, we show that the third body is responsible for causing a greater perturbation in the eccentricity and argument of periapsis, whereas the ellipticity coefficient of Titania causes great oscillations in the inclination. The zonal coefficient , when acting alone, has little influence on the orbital elements; however, when added to the other perturbations, it helps to promote a balance between third body and gravity coefficients, and then we have the so-called “protection mechanism” prolonging the lifetimes of the orbits.
The gravitational coefficients of Titania have not yet been used to investigate orbits close to its surface. Thus, we present an original study on the sensitivity of the lifetime as a function of some possible errors attributed to the values of
and
. Our results show that the hypothetical errors in the values of Titania’s gravitational coefficients can greatly increase or decrease the lifetime of the orbit, especially when the eccentricity is equal to
. To try to predict the lifetime as a function of these errors, we present multiple regression models for the lifetime as a function of the coefficients
and
. We found a function with a value of
equals to
, capable of providing a good approximation of the lifetime for values of
and
between
and
of the values given in [
27].
The results presented in this paper can help plan missions to Titania, which can bring important scientific data, in particular considering that the Uranus system has been studied mostly from Earth.