Dynamics of space particles and spacecrafts passing by the atmosphere of the Earth.
The close approach between two or more celestial bodies is very much studied in the literature. This phenomenon is involved in astronomical problems like the close approaches of three bodies [1-3], encounters of single stars and hard binaries , crossing between planets and a rotating star , or encounters between a star and a massive black hole . The study of the dynamics of comets during close approaches is also available in the literature [7, 8]. This phenomenon also has many applications in the astronautical field. Several missions used this technique to gain or lose energy, saving fuel in their maneuvers. One of the most famous applications was the Voyager mission, which visited the exterior planets of the Solar System based on the use of close approaches [9, 10]. Other applications of this maneuver are available in the literature, like the use of swing-bys to send a spacecraft to the giant planets [11-14] or even to the Sun 15], the use of Venus in a trip to Mars [16,17], the studies to make a three-dimensional close approach to Jupiter to change the orbital plane of the spacecraft , the use of the passages by the Moon to increase the energy of the spacecraft , and the use of multiple passages by the secondary body to find trajectories linking the primaries or the Lagrangian points [20, 21].
A good description of this maneuver is available in . There are many other researches available in the literature considering this problem, several of them optimizing parameters to obtain some desired results, like combining this approach with impulsive maneuvers . The eccentricity of the primaries was also included in some researches, like .
The present research studies the effects of the atmosphere of the Earth on the trajectory of a particle or a spacecraft that makes this close approach with the Earth [25, 26]. The mathematical model considers the problem as a circular restricted planar three-body problem  plus the atmospheric drag force. Some other papers also consider the behavior of a cloud of particles making the close approach [28, 29].
The equations of motion are numerically integrated starting at the periapsis in both directions of time, until the particle or spacecraft reaches distances far enough from the Earth and the system Sun-particle can be considered a two-body system. At those two points, one for each direction of time, it is possible to compute the two-body energy Sun-particle before and after the passage by the Earth for each particle. If the atmosphere was not present, this type of maneuver would be the standard swing-by described in the references cited previously. Now, with the inclusion of the atmosphere of the Earth, there is a drag force changing the trajectory of the particle or the spacecraft. This drag force is obtained by assuming dependence with the square of the velocity multiplied by the density of the atmosphere, which is calculated based on an exponential relation with the altitude. The objective of this research is to map the energy before and after this passage for different geometries of the approach, velocity, and altitude of the periapsis.
2. Definition of the Problem and Mathematical Model
It is assumed that the spacecraft comes from an orbit around the Sun (that can be elliptic or hyperbolic), outside the sphere of influence of the Earth. When approaching the Earth, its trajectory causes the entrance of the particle or spacecraft in the atmosphere and then its motion is governed by the gravity of the planet, assumed to be a point of mass, and the atmospheric drag. The hypothesis for the motion of the Earth is that it is travelling in a circular orbit around the center of mass of the system Sun-Earth. The spacecraft starts its motion in an orbit that is also around the center of mass of the Sun-Earth system. Figure 1 shows the basic geometry The particle or spacecraft approaches the Earth from point A, makes a close approach with the Earth, including the passage by the atmosphere, and then goes to point B. The perigee of the trajectory is marked by point P. The atmosphere of the Earth is represented by the circle around its center. At points A and B, the particle or spacecraft is assumed to be far from the Earth and its motion can be modeled by a two-body problem with the Sun. The equations of motion written in the rotating frame, which is a reference system that rotates together with the Sun-Earth system for each particle or spacecraft, are shown below (see (1) to (4)). The canonical system of units is also used, means that the Sun-Earth distance is the unit of distances, the total mass of the Sun and the Earth is the unit of mass, and the unit of time is defined in such a way that the period of the motion of the Sun-Earth system is 2n. So, the equations of motion are 
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (1)
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (2)
where [OMEGA] is the potential, given by
[OMEGA] = [1/2]([x.sup.2] + [y.sup.2]) + (1 - [mu])/[r.sub.1] + [mu]/[r.sub.2], (3)
where x and y represent the coordinates of the position of the particle/spacecraft; [r.sub.1] and [r.sub.2] are the distances from the particle/spacecraft to the Sun and the Earth, respectively; [mu] is the gravitational parameter of the Earth; [F.sub.x] and [F.sub.y] represent the components of the drag force ([??]), which is given by
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII], (4)
where [C.sub.B] = [C.sub.D]A/2m is the ballistic coefficient; [C.sub.D] is the drag coefficient, which is a parameter that considers the form of the particle or spacecraft; A is the cross section area of the particle or the spacecraft; [??] is the velocity of the spacecraft or the particles with respect to the atmosphere; m is its mass; [rho] is the density of the atmosphere, which can be estimated by an exponential function representing a barometric equilibrium:
[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII]. (5)
In (5), [[rho].sub.0] is the density of the atmosphere at an altitude [h.sub.0]; h is the altitude of the spacecraft; H is a constant that specifies the decay velocity of the density with the altitude. This is a usual approximation for the density of the atmosphere of the Earth, and the values of the constants involved are obtained by observations.
3. Identifying One Trajectory
The set of variables used to define each trajectory is the following:
(a) [psi]: the angle of approach, that is, the angle between the line of the periapsis of the trajectory of the spacecraft around the Earth and the line that connects the Sun and the Earth;
(b) [v.sub.p]: the velocity of the particle or spacecraft at the periapsis;
(c) [h.sub.p]: the altitude of the periapsis of the trajectory of the particle or spacecraft around the Earth.
So, the sequence of steps shown as follows can show the importance of the drag in the trajectories.
(i) A numerical integration is performed in backward time, with the particle or spacecraft starting at the periapsis until it reaches point A (Figure 1). This point is far enough from the Earth, so the two-body celestial mechanics can be used to obtain the two-body energy Sun-particle or spacecraft before the close approach.
(ii) Then, the numerical integration is made again, with the particle or spacecraft starting at the periapsis again, but this time in forward time. This process is done until the particle or spacecraft collides with the planet or reaches point B (see Figure 1), which is assumed to be far from the Earth such that the motion is again Keplerian and the two-body energy Sun-particle or spacecraft after the close approach is obtained.
So, by repeating this strategy for a large range of values of the ballistic coefficient, it is possible to understand the importance of the drag force in the process. The numerical integration is made using a Fourth order Runge-Kutta with step-size control and a numerical precision of [10.sup.-8].
The results are shown in the figures that display the ballistic coefficient (in [m.sup.2]/kg) in the horizontal axis and the energy before and after the close approach (in canonical units) in the vertical axis. Several values were considered for the velocity at periapsis [v.sub.p]: 0.3, 0.5, and 0.7 canonical units, where 1.0 canonical unit is equal to 29.78 km/s, that is, the velocity of the Earth in its motion around the Sun. Those values were chosen after making several simulations with different values. It is important to choose values that are not too low, because that causes too many captures by the Earth, and also those which are not too high, which causes the Swing-Bys to have very few effects on the trajectory of the spacecraft. Four values were used for the angle of approach [psi]: 0[degrees], 90[degrees], 180[degrees], and 270[degrees]. The idea was to consider values that cover the particular cases like maximum gains and losses of energy and regions with no effects from the Swing-Bys. Figures 2, 3, 4, 5, 6, 7, and 8 show the results, where [E.sub.-] and [E.sub.+] represent the energy before and after the Swing-By, respectively. The values for CB go from zero (a satellite with [C.sub.D] = 0), that represents a situation where the effects of the drag are not present, to 10 [m.sup.2]/kg, that would represent a satellite with a large area/mass ratio, in the order of 1000 [m.sup.2] for a 100 kg satellite, assuming [C.sub.D] = 2. It implies the use of a panel in the satellite. The value of [C.sub.D] varies from 2 to 2.3, depending on the interaction of the atmospheric constituents with the satellite surface .
There are several facts to be observed from those results. First, it is noticed that the presence of the atmospheric drag increases the variation of energy, always removing energy from the particle in both senses of time. The relation between the [C.sub.B] and the energy is almost linear for all the situations simulated here. The inclination of those lines increases with the velocity at the periapsis, because the drag force is proportional to the square of the velocity of the particle or spacecraft. This means that a higher value for the velocity at the periapsis causes higher values for the velocities during the whole trajectory inside the atmosphere, so the loss of energy is larger.
Considering, in some detail, the specific cases, it is clear that the maneuvers where [psi] =0[degrees] and [psi] = 180[degrees] have no change in the energy due to the Swing-By itself. The lines for the energy before and after the passages start at the same points, which means that, for [C.sub.B] = 0 (situation with no drag) the variation of energy is zero. This is in agreement with the literature . It is also clear that the orbits are all elliptic because the energy is negative before and after the close approach. Another interesting characteristic is that there is an inversion in the order of the lines. For [psi] = 0[degrees] the upper lines represent the higher velocities and for [psi] = 180[degrees] the upper lines represent the lower velocities. This fact can be explained by the geometry of the problem. The energy is related to the two-body system Sun-particle or spacecraft and the velocities are measured with respect to the Earth. So, to obtain the energy, it is necessary to make the vectorial addition of the velocities of the particle or spacecraft with respect to the Earth with the velocity of the Earth with respect to the Sun. The velocity of the Earth with respect to the Sun is always pointing towards the vertical axis, in the positive direction. The velocity of the particle or spacecraft with respect to the Earth points towards the positive direction of the vertical axis when [psi] = 0[degrees], so both velocities have the same sense and their magnitudes added together cause an increase in the energy with the increase of the velocity at the periapsis. In the situation where [psi] = 180[degrees], the velocity of the particle with respect to the Earth is in the negative direction of the vertical axis, so the velocities have opposite senses and their magnitudes have to be subtracted from each other to get its value with respect to the Sun, causing a decrease in the energy with the increase of the velocity at the periapsis.
It is also noted that the situation where [psi] = 0[degrees] shows the presence of hyperbolic and elliptic orbits, depending on the velocity at the periapsis.
The maneuvers with [psi] = 90[degrees] represent the conditions that give the maximum loss of energy due to the Swingby. This loss can be noted by the gap in the initial points of both energies before and after the passage. Note that this gap increases when the velocity at the periapsis decreases, due to the fact that lower velocities allow longer interactions between the bodies, and so the variation of energy is stronger.
On the opposite side, maneuvers with [psi] = 270[degrees] represent the situations that give the maximum gain of energy due to the Swing-by. Once again the gaps in the initial points show that fact. These gaps also increase when the velocity at periapsis increases, for the same reasons already explained. Another characteristic of this geometry is the crossing of the lines of energy before and after the close approach. It happens because the case with [psi] = 270[degrees] and no drag ([C.sub.B] = 0) has an increase in the energy, so the line of the energy before the maneuver starts at a point higher than the line of the energy after the maneuver. When considering the presence of the atmospheric drag, the loss in the energy starts to occur, and this loss increases with the ballistic coefficient. So, there are values for the [C.sub.B], marked by the crossing of the lines of the energy before and the energy after the maneuver, where the loss due to the atmospheric drag is equal to the gain due to the Swing-By and the net result is a maneuver with zero variation of energy. This point can even be used for a spacecraft that needs to realize a close approach with the Earth to complete the goals of the mission but does not want to change its orbit. After this point, the loss of energy is larger and the net result is a modification of the orbit that reduces its energy.
After that, a more detailed study is made for near-parabolic orbits. Those are the orbits with two-body energy Sun-particle or spacecraft near zero in the situation without drag. In the situation where [psi] = 90[degrees], that causes loss of energy due to the Swing-By, the initial orbit is hyperbolic (positive energy) and the final orbit is elliptic (negative energy). This is a situation where a capture occurred. The difference in those energies increases with the ballistic coefficient, so the presence of the atmospheric drag works in favor of those captures and orbits with larger energy can be captured with the help of the atmosphere. At the same time, those orbits become even more closed, with lower values of the energy.
On the opposite side, the situation where [psi] = 270[degrees] generates gains of energy due to the Swing-By. Considering the case with no drag ([C.sub.B] = 0), the initial orbit is elliptic (negative energy) and the final orbit is hyperbolic (positive energy). This is a situation where an escape occurred. This also is in agreement with the literature . When adding the presence of drag, the mechanism of removal of the energy by the atmosphere starts. Once again there is a point where the variation of energy is zero, as a consequence of the equilibrium between the gain due to the swing-by and the loss due to the atmosphere. The difference among those energies increases with the ballistic coefficient, as expected. So, the presence of the atmosphere causes captures, even in this geometry.
After that, a study is made to take into account the effects of the periapsis altitude in the passage of the particles or the spacecraft by the atmosphere. The simulations are made for the velocity of the periapsis of 1.0 canonical units and angle of approach [psi] = 270[degrees]. Other situations were simulated with similar results, so they are not shown here. The values used for the periapsis altitude are 130,140, and 150 km above the surface of the Earth. Figure 8 shows the results. As expected, the effects of the passage are much stronger for lower values of the periapsis altitude because the density of the atmosphere is larger. In that sense, the results quantify these effects and show that the variation in energy is about four times larger when the periapsis is 130 km when compared to the situation where the periapsis is 140 km and fifteen times larger when compared to the value of periapsis altitude of 150 km. Another characteristic shown in the results is that the relation between energy and ballistic coefficient is near linear in most of the cases, but it starts to become quadratic when the periapsis distance decreases. This is a consequence of the exponential model for the density, since the distances are becoming larger and linear approximation of the density is no longer valid.
A numerical sequence of steps to measure the effects of the atmosphere in a close approach between a particle or a spacecraft with a planet that has an atmosphere is developed and implemented. This approach is then applied to a particle or a spacecraft passing by the atmosphere of the Earth. The results show regions where escapes and captures occur and also quantify the effects of the atmosphere in the trajectories. The atmosphere reduces the energy of the particles or spacecraft, helping captures to occur, even for geometries where the energy should increase due to the close approach. The effects of the velocity at the periapsis are measured in all the situations, in order to quantify its effects by the present algorithm. Near-parabolic orbits are studied in more detail, to emphasize the study of the mechanism of captures of the particles or spacecraft by the atmosphere of the Earth. The study of the variation of the altitude of the periapsis is also made, showing the limits where the relation between energy and ballistic coefficients is linear. On the astronautical side, the present research can be used to help to plan missions with the goal of capturing a spacecraft by the atmosphere of the Earth, giving numerical results that help to define values for the ballistic coefficient to complete a capture. In terms of astronomical problems, the algorithm developed here can help to predict the percentage of captured particles as a function on the initial conditions at the periapsis or before the close approach, by using the values obtained after the backward integration.
Conflict of Interests
The authors declare that there is no conflict of interests regarding the publication of this paper and that they do not have a direct financial relation that might lead to one.
The authors wish to express their appreciation for the support provided by Grant Nos. 473387/2012-3 and 304700/2009-6 from the National Council for Scientific and Technological Development (CNPq); Grants Nos. 2011/09310-7, 2011/081713, and 2011/13101-4 from the Sao Paulo Research Foundation (FAPESP), and the financial support from the National Council for the Improvement of Higher Education (CAPES).
 V. Szebehely, "Triple close approaches in the problem of three bodies," Celestial Mechanics, vol. 8, no. 2, pp. 163-167, 1973.
 J. Laskar and C. Marchal, "Triple close approach in the three-body problem: a limit for the bounded orbits," Celestial Mechanics, vol. 32, no. 1, pp. 15-28, 1984.
 C. Marchal, "n-Body problem: selective effects of triple close approaches," Acta Astronautica, vol. 7, no. 1, pp. 123-126, 1980.
 J. G. Hills and L. W. Fullerton, "Computer Simulations of close encounters between single stars and hard binaries," The Astronomical Journal, vol. 85, pp. 1281-1291, 1980.
 P. B. Ivanov and J. C. B. Papaloizou, "Close encounters of a rotating star with planets in parabolic orbits of varying inclination and the formation of hot Jupiters," Celestial Mechanics and Dynamical Astronomy, vol. 111, no. 1-2, pp. 51-82, 2011.
 A. Khokhlov, I. D. Novikov, and C. J. Pethick, "Strong effects during close encounters of a star with a massive black hole," The Astrophysical Journal Letters, vol. 418, no. 1, pp. 181-186, 1993.
 J. L. Zhou, Y. Y. Sun, J. Q. Zheng, and M. J. Valtonen, "The transfer of comets from near-parabolic to short-period orbits: map approach," Astronomy and Astrophysics, vol. 364, no. 2, pp. 887-893, 2000.
 E. Everhat, "Close encounters of comets and planets," The Astronomical Journal, vol. 74, no. 5, p. 735, 1969.
 M. A. Minovich, "A method for determining interplanetary free-fall reconnaissance trajectories," JPL Technical Memo 312130, JPL, Pasadena, Calif, USA, 1961.
 G. A. Flandro, "Fast reconnaissance missions to the outer solar system utilizing energy derived from the gravitational field of Jupiter," Astronautica Acta, vol. 12, no. 4, pp. 329-337, 1966.
 L. A. DAmario, D. V. Byrnes, and R. H. Stanford, "Interplanetary trajectory optimization with application to Galileo," Journal of Guidance, Control, and Dynamics, vol. 5, no. 5, pp. 465-471, 1982.
 L. A. D'Amario, D. V Byrnes, andR. H. Stanford, "A new method for optimizing multiple-flyby trajectories," Journal of guidance and control, vol. 4, no. 6, pp. 591-596, 1981.
 N. J. Strange and J. M. Longuski, "Graphical method for gravity-assist trajectory design," Journal of Spacecraft and Rockets, vol. 39, no. 1, pp. 9-16, 2002.
 T. T. McConaghy, T. J. Debban, A. E. Petropoulos, and J. M. Longuski, "Design and optimization of low-thrust trajectories with gravity assists," Journal of Spacecraft and Rockets, vol. 40, no. 3, pp. 380-387, 2003.
 A. A. Sukhanov, "Close approach to sun using gravity assists of the inner planets," Acta Astronautica, vol. 45, no. 4-9, pp. 177185, 1999.
 W. M. Hollister and J. E. Prussing, "Optimum transfer to mars via venus," Astronautica Acta, vol. 12, no. 2, pp. 169-179, 1966.
 S. A. Striepe and R. D. Braun, "Effects of a venus swingby periapsis burn during an Earth-Mars trajectory," Journal of the Astronautical Sciences, vol. 39, no. 3, pp. 299-312, 1991.
 R. Carvell, "Ulysses--the sun from above and below," Space, vol. 1, pp. 18-55, 1986.
 D. Dunham and S. Davis, "Optimization of a multiple lunar-swing-by trajectory sequence," Journal of Astronautical Sciences, vol. 33, no. 3, pp. 275-288, 1985.
 A. F. B. A. Prado and R. A. Broucke, "Transfer orbits in restricted problem," Journal of Guidance Control and Dynamics, vol. 18, no. 3, pp. 593-598, 1995.
 A. F. B. A. Prado and R. A. Broucke, "Transfer orbits in the Earth-moon system using a regularized model," Journal of Guidance, Control, and Dynamics, vol. 19, no. 4, pp. 929-933, 1996.
 R. A. Broucke, "The celestial mechanics of gravity assist," in Proceedings of the AIAA/AAS Astrodynamics Conference, AIAA paper 88-4220, Minneapolis, Minn, USA, August 1988.
 A. E B. A. Prado, "Powered swing-by," Journal of Guidance, Control, and Dynamics, vol. 19, no. 5, pp. 1142-1147, 1996.
 A. F. B. A. Prado, "Close-approach trajectories in the elliptic restricted problem," Journal of Guidance, Control, and Dynamics, vol. 20, no. 4, pp. 797-802, 1997
 A. F. B. A. Prado and R. A. Broucke, "Effects of atmospheric drag in swing-by trajectory," Acta Astronautica, vol. 36, no. 6, pp. 285-290, 1995.
 W. Schulz, A. F. B. A. Prado, and R. V de Moraes, "An analytical and numerical study of plane change maneuvers using aerodynamic force," Journal of the Astronautical Sciences, vol. 50, no. 3, pp. 289-303, 2003.
 V. G. Szebehely, Theory of Orbits, Academic Press, New York, NY, USA, 1967.
 V. M. Gomes and A. F. B. A. Prado, "Swing-by maneuvers for a cloud of particles with planets of the solar system," WSEAS Transactions on Applied and Theoretical Mechanics, vol. 3, no. 11, pp. 869-878, 2008.
 V. M. Gomes and A. F. B. A. Prado, "A study of the impact of the initial energy in a close approach of a cloud of particles," WSEAS Transactions on Mathematics, vol. 9, no. 10, pp. 811-820, 2010.
 O. Montenbruck and E. Gill, Satellite Orbits--Models, Methods, and Applications, Springer, Berlin, Germany, 2000.
Vivian Martins Gomes, (1) Antonio Fernando Bertachini de Almeida Prado, (2) and Justyna Golebiewska (3)
(1) Grupo de Dinamica Orbital e Planetologia, Universidade Estadual Paulista, FEG/UNESP, 12516-410 Guaratingueta, SP, Brazil
(2) Instituto Nacional de Pesquisas Espaciais (INPE), 12227-010 Sao Jose dos Campos, SP, Brazil
(3) Astronomical Observatory Institute, Faculty of Physics, Adam Mickiewicz University, 60-286 Poznan, Poland
Correspondence should be addressed to Antonio Fernando Bertachini de Almeida Prado; firstname.lastname@example.org
Received 11 September 2013; Accepted 7 November 2013
Academic Editors: G. Pascazio and C. Zhang
|Printer friendly Cite/link Email Feedback|
|Title Annotation:||Research Article|
|Author:||Gomes, Vivian Martins; de Almeida Prado, Antonio Fernando Bertachini; Golebiewska, Justyna|
|Publication:||The Scientific World Journal|
|Date:||Jan 1, 2013|
|Previous Article:||Establishing an efficient way to utilize the drought resistance germplasm population in wheat.|
|Next Article:||Effects of resveratrol on methotrexate-induced testicular damage in rats.|