Since the first orbital launch in 1957, the number of artificial objects in Earth orbit has been steadily increasing. This
|
|
- Constance Thornton
- 5 years ago
- Views:
Transcription
1 SpaceOps Conferences 28 May - 1 June 2018, 2018, Marseille, France 2018 SpaceOps Conference / Collision probability through time integration Implementation and operational results Vincent Schaeffer a, Sophie Laurens a, Pierre Seimandi b, and Florian Delmas a a Centre National d Etudes Spatiales (CNES), Toulouse, France b GMV Innovating Solutions, Toulouse, France In 2012, V.T. Coppola proposed a new formulation for the computation of the probability of collision between two satellites based on a time-integration of the probability density flux entering the surface the hard body radius sphere. This formulation takes into account the evolution of the covariance over time, as well as the velocity uncertainty. It should therefore give access to a more realistic probability of collision, especially for low-velocity encounters. This paper describes the work currently done at CNES on the method, which was tested on more than generated cases, in addition to a selection of past collision risks. The paper highlights several issues encountered during the implementation of the method and focuses on a few conjunctions of interest. I. Introduction Since the first orbital launch in 1957, the number of artificial objects in Earth orbit has been steadily increasing. This has led to a corresponding increase in the number of threats to active satellites from hyper-velocity collision, putting in jeopardy crucial services that benefit human society. To mitigate such risks, one of the approaches commonly adopted consists in the planning and the execution of collision avoidance maneuvers for operational satellites. This relies on the analysis of the probability of collision between the satellite of interest and any space object entering a safety volume around it. An accurate estimation of this probability of collision is thus primordial for the management of the spacecraft population. An under-estimated probability could lead to serious risks being ignored, while an over-estimated one could trigger unnecessary maneuvers, leading to mission perturbations and reducing of the operational lifespan of the satellite. To estimate this probability, CNES currently relies on a simple model based on the linearized states of the objects at time of closest approach (TCA) [1, 2]. Assuming that the encounter is of short-term nature, the relative trajectory is locally approximated as rectilinear, the positional uncertainty is considered as time-independent and the velocity uncertainty is neglected. The probability of collision can then be computed by a two-dimensional integration in the encounter plane, making this method easy to implement and fast to compute (an important feature given the number of alerts treated daily). This 2D-model is however not applicable to every encounter. By construction, it is not adapted to low-velocity encounters, where the two objects evolve close to one another for long periods of time. During these encounters, the relative trajectory and the states uncertainties can change significantly: the assumptions described above are simply not valid anymore. More importantly, significant discrepancies have also been observed for several high-velocity alerts. While the 2D-probability computed for these alerts was negligible, the probability returned by a Monte-Carlo (MC) algorithm was sufficiently high to warrant an evasion maneuver. This is problematic operationally since it implies that a serious collision risk was missed. For these reasons, alternative methods for the computation of the probability of collision were investigated. Among them, the method proposed by Coppola [3] in 2012 is of particular interest. This method, referred to here as the 3D-probability method, is based on the integration of the flux of a time-dependent PDF entering the surface the hard body radius (HBR) sphere over a given time-span. Compared to the 2D-probability method currently used operationally, it takes into account the velocity uncertainty and the evolution of the covariance over time. It should therefore give access to a more realistic probability of collision, especially for low-speed encounters, and provide an alternative to expensive Monte-Carlo computations. More precisely, the method defines the probability density flux penetrating the HBR sphere at each time t. The probability of collision is then computed by integrating this flux over the time-span considered. As an example, figure 1 represents the evolution of the flux during a high-velocity encounter. The bell-shaped flux seen Corresponding author: vincent.schaeffer@cnes.fr Corresponding author: sophie.laurens@cnes.fr 1 Copyright 2018 by CNES. Published by the American Institute of Aeronautics and Astronautics, Inc., with permission.
2 here is typical of this kind of encounters, for which the main contribution to the probability of collision occurs over only a few seconds. Fig. 1 Typical evolution of the probability density flux during a high-velocity encounter (TCA is at t = 0) Coppolas formulation relies on the following hypotheses: The uncertainty on the position and velocity of each object remains Gaussian at each time: For our applications, the initial covariance matrices are usually propagated linearly, through a transition matrix. This assumption is thus met by construction. One should however be careful when propagating over long periods of time. Although a covariance matrix could still be computed, there is no guarantee that it would reflect the real uncertainty (which may not even be Gaussian anymore). Here, we will assume that this type of propagation is valid for the time-spans considered; The dynamic model and the probability distribution function of each object are independent: This is a valid assumption for most encounters. However, it would not be appropriate for objects flying in formation, with trajectory correction based upon the location of other objects; Each collision trajectory must cross the combined HBR sphere instantly: This assumption means no collision trajectory must linger on the surface of the HBR sphere for any amount of time. While it is difficult to verify this assumption, we expect that it will be met in most cases; Each collision trajectory crosses the combined HBR sphere only once: Although we expect this assumption to be met for most encounters, it can be disputable for objects staying close to each other for a large period of time, such as objects flying in formation or for some geostationary conjunctions; While these assumptions seem reasonable in most cases, it is difficult to assess beforehand whether they are actually met. Failing that, the method was tested against more than generated collision risks, as well as on a selection of past collision risks. The goal of this first validation was to apply Coppola s method on a large variety of cases in order to identify implementation issues and abnormal behaviors. After a short description of the cases tested in section II, this paper focuses on a few points of interest highlighted by our study. More specifically, section III discusses the formulation for the estimation of the encounter duration proposed by Coppola in a second paper [4]. Section IV then discusses the use of Lebedev s quadrature [5] for the integration over the surface of the HBR sphere and highlight a possible lack of convergence for some configurations. Finally, section V presents two high-velocity cases resulting in a very under-estimated 2D-probability, despite being significant collision risks. II. Test cases description As said in section I, Coppola s method have been tested on more than generated collision risks (28512 exactly). Each case defines the combined HBR, the mean states of the objects at TCA and the corresponding position and velocity uncertainties. These collision risks correspond to 44 unique encounters with respect to the orbital parameters. For each one of these encounters, three combined HBR values (1 m, 10 m or 100 m) are then considered, as well as various initial uncertainties. For the secondary object, the 1-sigma volume of the covariance ranges from 10 m 3 to m 3, every intermediate order of magnitude being covered. The aspect ratio of this volume (i.e. the degree of stretching in a specific direction) is between 1 to 1000 for the most part, but it can reach about for a few cases. The covariance of the primary object is negligible. 2
3 A few key-characteristics on these encounters are listed in table 1. This table shows that, although only 44 unique encounters are considered, a wide-variety of conjunctions are actually represented by these cases. Note that the collision risks tested do not aim to be always realistic (in particular with regard to the covariance used). The main goal of these cases was to test the limits of the method to identify abnormal behaviors or implementation issues. In addition to these generated cases, a selection of past collision risks were also studied. When discussed in this paper, these cases are indicated as "real collision risk". All the results presented here use a basic Keplerian propagation for both the mean orbits and the covariance states (only the central gravitational force is considered). Table 1 Generated cases characteristics Field Values Combined HBR 1 m, 10 m and 100 m Miss distance 1 m, 10 m, 100 m and 2000 m Relative speed 0.79 m s 1 to m s 1 Orbit classification LEO-1/2/3, GEO, HEO Altitude 340 km to km Eccentricity to 0.67 Inclination 0.11 to The probabilities computed for these generated cases are compared figure 2, with the 2D-probability on the x-axis and the 3D-probability on the y-axis. The solid black line indicates a 2D-probability equal to the 3D-probability and the dashed black lines mark a difference of one order of magnitude between the two probabilities. The dotted gray lines represent the 10 8 threshold for both probabilities (below this value, we consider the probability as negligible). High-velocity encounters are plotted in orange and low-velocity encounters in blue, with a threshold of 50 m s 1 to discriminate the two kinds of encounter. In addition to these two graphs, table 2 presents the breakdown of the cases according to the relative error between the 2D-probability and the 3D-probability. As a note, % of these cases are high-velocity encounters and % are the low-velocity encounters. As expected, the two probabilities are in good agreement for most high-velocity encounters. Only 2.08 % yield a difference of more than one order of magnitude and most of these cases return a negligible probability with both methods. However, a few high-velocity cases yield a negligible 2D-probability despite a relatively high 3D-probability. A group in particular returns a 3D-probability above 10 5, which is almost high enough to warrant an evasion maneuver. These cases will be further analyzed in section V, along with a past collision risk exhibiting the same problem. For low-velocity encounters, results are still in relatively good agreement % of the cases result in less than 10% difference and % result in a difference lower than one order of magnitude. However, a notable number of cases yield a negligible 2D-probability and a significant 3D-probability, or vice-versa. These differences can be explained by the nature of these encounters. The 2D-probability model is not adapted and the probability of collision is often poorly evaluated. Table 2 Breakdown of the generated cases according to the relative error on the probability of collision Relative error High-speed Low-speed Less than 1 % % % Between 1 % and 10 % 3.05 % % More than 10 % and less than one order of magnitude 1.79 % % More than one order of magnitude 2.08 % 8.87 % More than one order of magnitude and any P c > % 4.08 % 3
4 (a) All cases (b) Both 2D and 3D probabilities above Fig. 2 Comparison of the 2D and 3D-probability of collision for the generated cases. High-velocity encounters (> 50 m s 1 ) are plotted in orange and low-velocity encounters ( 50 m s 1 ) in blue. III. Encounter duration As described in section I, the formulation proposed by Coppola consists of a time-integration of the probability density flux. To estimate the probability of collision, the time-span to integrate must thus be defined. Coppola addresses this particular problem in a second paper [4], in which he proposes a formulation that characterizes the duration of the encounter using the short-term hypotheses. This formulation yields two parameters, τ 0 and τ 1, that define the duration of the encounter relatively to the TCA. These parameters are: τ 0 = 2α c σ v + q 0 δ max and τ 1 = + 2αc σ v + q 0 δ min, v 0 v 0 where v 0 is the relative speed at TCA and σ v is the relative position uncertainty normal to the encounter plane. The terms q 0, δ min and δ max depend on the relative position in the encounter plan, the combined HBR and the initial uncertainty. The terms δ min and δ max are often negligible compared to σ v and the parameter α c controls the length of the time span returned, while the term q 0 defines its shift relatively to the TCA. Note that the parameter α c must be high enough to ensure a low approximation error under the short-term encounter hypotheses. Coppola indicates in his paper that a value of is enough to yield an error below Here, we chose to use a higher value: n σ = 2α c = 15 (i.e. α c = 10.61). This formulation works well for the vast majority of the cases studied. Figure 1 is representative of the result obtained when it is applied to high-velocity encounters. The time span returned is shifted relatively to the TCA (located at t = 0) and centered on the bell-shaped flux. The shift of the interval computed by the method is however incorrect in a few cases, resulting in an under-evaluated probability induced by a flux either truncated or completely missed. Figure 3 gives an example of a past collision risk impacted by this problem. It is a high-velocity encounter with a relative speed around 14.7 km s 1. While still under investigation, a few things can be considered to palliate the issue. The flux at the beginning and at the end of the integration interval could for example be compared to the maximum value computed in order to check whether the whole flux was actually integrated. If not, the integration interval could be extended and the probability recomputed. Another possible solution could be to automatically include the TCA (i.e. t = 0) in the interval since a lot of cases impacted by the problem (but not all) exhibit a shift too important and the encounter occurs somewhere between the time span computed and the TCA. Although Coppola s formulation is designed for short-term encounters, it can also be useful for low-velocity encounters. Figure 4 shows the probability density flux for a typical low-velocity encounter of a real GEO satellite. During this conjunction, the two objects evolve close to each other for a long period of time with a relative speed of 2 m s 1. The time span returned by Coppola s formulation is still correctly shifted (by about 1800 s) and the flux keeps a similar shape as the one observed for high-velocity encounters, the only difference being the duration of the encounter 4
5 Fig. 3 Truncated probability density flux (real collison risk) which last about 200 s instead of a few seconds. Despite the low-velocity, the relative trajectory near TCA is actually rectilinear, which could explain why Coppola s formulation is still applicable here. The probability density flux for another low-velocity conjunction is given figure 5. In this encounter, the secondary object orbits around the primary object for a long period of time, very closely at first and then slowly drifting away. This results in an important initial risk, followed by periodic spikes in the probability density flux, each one adding to the probability of collision computed. Defining an encounter duration is thus difficult and Coppola s formulation should not be used blindly. Because of this kind of encounter, low-velocity conjunctions should always be analyzed carefully (even if most of the risks studied so far behave as depicted figure 4). Fig. 4 Probability density flux for typical low-velocity encounter (real collision risk) Fig. 5 Probability density flux for a periodic, low-velocity encounter (generated case) 5
6 IV. Integration on the surface of the HBR sphere The formulation proposed by Coppola consists of a time-integration of the probability density flux over a chosen time-span. At each time t, this density probability flux can be computed through a simple two-dimensional integration over the surface of the combined HBR sphere. To compute this surface integral, Coppola suggests to use Lebedev s quadrature [5] with 1454 points (i.e. order 65), noting no significant difference when using higher orders. Since this quadrature uses points evenly distributed on the surface of the sphere, it requires fewer total grid points to achieve an accuracy similar to a standard (θ, φ) discretization and is thus very efficient. For our applications, this quadrature suffers however from several limitations. The determination of the grid weights is achieved by enforcing the scheme to integrate exactly all spherical harmonics up to a given order, which implies solving a set of nonlinear equations. Currently, these weights have been tabulated up to order 131 (5810 points). Although this order is more than enough for most cases, a default of convergence was observed for a few encounters. An example of such an encounter is given figure 6. It depicts the probability density flux obtained using Lebedev s quadrature with 1454 points (figure 6a) or 5810 points (figure 6b) and using a converged (θ, φ) discretization (figure 6c). For this particular case, the covariance in the plane orthogonal to the velocity is very small compared to the combined HBR (100 m) and Lebedev s quadrature is not able to reach convergence, even when using the highest order available. The various probabilities computed for this case are listed in table 3. The differences observed are relatively minor (less than 10 %) and it could be argued that Lebedev s quadrature is accurate enough. A possible default of convergence is however not acceptable operationally, as a poorly evaluated probability of collision could potentially lead to a missed collision risk. The use of an alternate integration scheme able to ensure convergence is thus recommended over Lebedev s quadrature. Even if higher orders were available, ensuring convergence with this quadrature would require either to determine the order to use beforehand (which can be difficult) or to test the available grids incrementally in order to check numerically if the wanted accuracy is reached (which is inefficient, since very few grid points can be reused). Table 3 Probabilities of collision corresponding to the case depicted figure 6 Method Quadrature Probability 2D-probability N.A D-probability (6a) Lebedev (1454 points) D-probability (6b) Lebedev (5810 points) D-probability (6c) Converged (θ, φ) discretization V. Under-estimated 2D-probability of collision Although the 2D and the 3D probabilities of collision are in good agreement for most high-velocity encounters (see section II), that is not systematically the case. In particular, a few conjunctions result in a 2D-probability largely under-estimated compared to the probability of collision returned by Coppola s method. The probability density flux corresponding to two of these cases are presented figures 7 and 8, the first being a generated case and the second being a past collision risk. In addition to these two graphs, table 4 lists the probabilities computed for these conjunctions, including the ones returned by a Monte-Carlo algorithm with 10 billions samples. For both cases, the 3D-probability and the Monte-Carlo probability are very close and sufficiently high to warrant further analysis, while the 2D-probability is under-estimated by at least one order of magnitude. The examples presented here illustrate the two kinds of behavior we observed so far. The probability density flux depicted figure 7 is comprised of two distinct bell-shaped section, located on both side of the TCA and shifted by less than a second. Conversely, the case depicted figure 8 behaves as most other high-velocity cases, characterized by a single bell-shaped flux near TCA lasting about one second. The mechanisms leading to such behaviors are not yet fully understood and must be further investigated. However, it is interesting to note that, for both cases, the covariance of the secondary object is particularly stretched in the direction of the orbital momentum. 6
7 (a) Lebedev s quadrature with 1454 points Fig. 6 (b) Lebedev s quadrature with 5810 points (c) Converged (θ, φ) discretization Probability density flux computed with: (6a) Lebedev s quadrature with 1454 points, (6b) Lebedev s quadrature with 5810 points, (6c) a converged (θ, φ) discretization (generated case) Fig. 7 Dual-bell shaped probability density flux for an encounter yielding a negligible 2D-probability (generated case) 7
8 Table 4 Probabilities of collision computed for the cases depicted figures 7 and 8 Figure 2D-Probablity 3D-Probability MC-Probability Figure Figure Fig. 8 Probability density flux for an encounter yielding negligible 2D-probability (real collision risk) VI. Conclusion & future work The formulation proposed by Coppola [3] for the computation of the probability of collision was tested on a wide-variety of cases, as well as on a few past encounters. As expected, the results for most high-velocity encounters are in good agreement with the 2D-probability of collision currently used operationally. Despite these good results, several high-velocity encounters resulted in important differences the two methods. In particular, some of these encounters yield 2D-probability significantly lower than the 3D-probability and the Monte-Carlo probability, implying that a real collision risk was potentially missed. These cases are thus very problematic operationally and must be further analyzed in order to identify the phenomena causing these differences. For low-velocity encounters, differences between the 2D-probability and the 3D-probability are more pronounced. Although it is not surprising since this kind of encounter is beyond the scope of the 2D-probability model, the results presented stress out the necessity to carefully analyze low-velocity encounters. Even if the 3D-probability seems better suited for this kind of conjunctions, properly defining the duration of the encounter can be difficult, especially for recurring risks. Moreover, the probabilities computed using Coppola s formulation should be considered with caution, as the model relies on hypotheses that might not always be satisfied. Further validations on the method are currently planned in order to increase the confidence in the probability returned. The method will be applied to a large number of recent collision risks alerts and compared to the 2D-probability of collision. The Monte-Carlo probability will then be computed for every case resulting in a strong difference between the two probabilities. A reflection on the application of the method in an operational context is also undergoing, one of the main concerns being the computational time. While relatively quick to compute, the method is still significantly slower than the 2D-probability of collision. Since the computation of the probability of collision done by CNES takes into account a possible dilatation or contraction of the covariance [6], assessing the risk for a single conjunction requires hundreds of individual evaluations. Applying the 3D-probability to every collision risk alert received would thus be inefficient. And while the relative speed at TCA can be used to decide the method to use for a specific case, the results obtained so far indicate that it is not sufficient. Additional criteria must be investigated. References [1] Patera, R., General Method for Calculating Satellite Collision Probability, Journal of Guidance, Control, and Dynamics, Vol. 24, 2001, pp doi: / [2] Alfano, S., Satellite Collision Probability Enhancements, Journal of Guidance Control and Dynamics, Vol. 29, 2006, pp
9 [3] Coppola, V., Including velocity uncertainty in the probability of collision between space objects, AAS/AIAA, Astrodynamics Specialist Conference, AAS Paper [4] Coppola, V., Evaluating the short encounter assumption of the probability of collision formula, AAS/AIAA, Astrodynamics Specialist Conference, AAS Paper [5] Lebedev, V., and Laikov, D., A quadrature formula for the sphere of the 131st algebraic order of accuracy, Russian Academy of Sciences Doklady Mathematics, Vol. 59, No. 3, 1999, pp [6] Laporte, F., JAC software, solving conjunction assessment issues, AMOS Technical Conference, Maui HI,
Conjunction Risk Assessment and Avoidance Maneuver Planning Tools
Conjunction Risk Assessment and Avoidance Maneuver Planning Tools Saika Aida DLR German Space Operations Center (GSOC), Oberpfaffenhofen, 834 Weßling, Germany +49 8153 8-158, saika.aida@dlr.de ASTRACT
More informationDEFINITION OF A REFERENCE ORBIT FOR THE SKYBRIDGE CONSTELLATION SATELLITES
DEFINITION OF A REFERENCE ORBIT FOR THE SKYBRIDGE CONSTELLATION SATELLITES Pierre Rozanès (pierre.rozanes@cnes.fr), Pascal Brousse (pascal.brousse@cnes.fr), Sophie Geffroy (sophie.geffroy@cnes.fr) CNES,
More informationJAC Conjunction Assessment
JAC Conjunction Assessment SSA Operators Workshop Denver, Colorado November 3-5, 2016 François LAPORTE Operational Flight Dynamics CNES Toulouse Francois.Laporte@cnes.fr SUMMARY CA is not an easy task:
More informationA method for calculating probability of collision between space objects
RAA 2014 Vol. 14 No. 5, 601 609 doi: 10.1088/1674 4527/14/5/009 http://www.raa-journal.org http://www.iop.org/journals/raa Research in Astronomy and Astrophysics A method for calculating probability of
More informationThe Applicability Research of Nonlinear Collision Probability
The Applicability Research of Nonlinear Collision Probability XIONG Yongqing, XU Xiaoli Purple Mountain Observatory Chinese Academy of Sciences Nanjing, Jiangsu 210008, China XU Xiaoli Outline Introduction
More informationPROBABILITY OF COLLISION WITH SPECIAL PERTURBATIONS DYNAMICS USING THE MONTE CARLO METHOD
AAS 11-435 PROBABILITY OF COLLISION WITH SPECIAL PERTURBATIONS DYNAMICS USING THE MONTE CARLO METHOD Chris Sabol, Christopher Binz and Alan Segerman, Kevin Roe, and Paul W. Schumacher, Jr. INTRODUCTION
More informationKRATOS: KOLLISION RISK ASSESSMENT TOOL IN ORBITAL ELEMENT SPACES
KRATOS: KOLLISION RISK ASSESSMENT TOOL IN ORBITAL ELEMENT SPACES Joshua T. Horwood, Navraj Singh, and Jeffrey M. Aristoff Numerica Corporation, 5042 Technology Parkway, Suite 100, Fort Collins, CO 80528
More informationAN ANALYTICAL SOLUTION TO QUICK-RESPONSE COLLISION AVOIDANCE MANEUVERS IN LOW EARTH ORBIT
AAS 16-366 AN ANALYTICAL SOLUTION TO QUICK-RESPONSE COLLISION AVOIDANCE MANEUVERS IN LOW EARTH ORBIT Jason A. Reiter * and David B. Spencer INTRODUCTION Collision avoidance maneuvers to prevent orbital
More informationPeriodicity Characterization of Orbital Prediction Error and Poisson Series Fitting
Periodicity Characterization of Orbital Prediction Error and Poisson Series Fitting Bai Xianzong *, Chen Lei, Tang Guoin College of Aerospace and Materials Engineering, National University of Defense Technology,
More informationWhen Does the Uncertainty Become Non-Gaussian. Kyle T. Alfriend 1 Texas A&M University Inkwan Park 2 Texas A&M University
When Does the Uncertainty Become Non-Gaussian Kyle T. Alfriend Texas A&M University Inkwan Park 2 Texas A&M University ABSTRACT The orbit state covariance is used in the conjunction assessment/probability
More informationCOLLISION RISK ASSESSMENT AND MITIGATION STRATEGY FOR THE GSOC GEO SATELLITES
COLLISION RISK ASSESSMENT AND MITIGATION STRATEGY FOR THE GSOC GEO SATELLITES Saika Aida (1), Michael Kirschner (2), Florian Meissner (3) (1) DLR German Space Operations Center (GSOC), Münchner Str.20,
More informationCOVARIANCE DETERMINATION, PROPAGATION AND INTERPOLATION TECHNIQUES FOR SPACE SURVEILLANCE. European Space Surveillance Conference 7-9 June 2011
COVARIANCE DETERMINATION, PROPAGATION AND INTERPOLATION TECHNIQUES FOR SPACE SURVEILLANCE European Space Surveillance Conference 7-9 June 2011 Pablo García (1), Diego Escobar (1), Alberto Águeda (1), Francisco
More informationAccurate numerical orbit propagation using Polynomial Algebra Computational Engine PACE. ISSFD 2015 congress, Germany. Dated: September 14, 2015
Accurate numerical orbit propagation using Polynomial Algebra Computational Engine PACE Emmanuel Bignon (), Pierre mercier (), Vincent Azzopardi (), Romain Pinède () ISSFD 205 congress, Germany Dated:
More informationGEO protected region: ISON capabilities to provide informational support for tasks of spacecraft flight safety and space debris removal
Russian Academy of Sciences Keldysh Institute of Applied Mathematics GEO protected region: ISON capabilities to provide informational support for tasks of spacecraft flight safety and space debris removal
More informationCONJUNCTIONS AND COLLISION AVOIDANCE WITH ELECTRODYNAMIC TETHERS. Eugene M. Levin Electrodynamic Technologies, LLC
CONJUNCTIONS AND COLLISION AVOIDANCE WITH ELECTRODYNAMIC TETHERS Eugene M. Levin Electrodynamic Technologies, LLC Abstract Electrodynamic tethers have very promising applications in debris removal and
More informationIAC-16.A Jason A. Reiter a *, David B. Spencer b
IAC-16.A6.7.5 Trading Spacecraft Propellant Use and Mission Performance to Determine the Optimal Collision Probability in Emergency Collision Avoidance Scenarios Jason A. Reiter a *, David B. Spencer b
More informationATTITUDE CONTROL MECHANIZATION TO DE-ORBIT SATELLITES USING SOLAR SAILS
IAA-AAS-DyCoSS2-14-07-02 ATTITUDE CONTROL MECHANIZATION TO DE-ORBIT SATELLITES USING SOLAR SAILS Ozan Tekinalp, * Omer Atas INTRODUCTION Utilization of solar sails for the de-orbiting of satellites is
More informationSSTD = Standard deviation SMA = Semi Major Axis
- 1 C - EPC-95-212 NOVEL ORBT RASNG STRATEGY MAKES LOW THRUST COMMERCALLY VABLE. ARNON SPTZER* ABSTRACT A new technique for utilizing low thrust Electric Propulsion System (EPS) for orbit raising has been
More informationRELATIVE MISSION ANALYSIS FOR PROBA 3: SAFE ORBITS AND CAM
5TH INTERNATIONAL CONFERENCE ON SPACECRAFT FORMATION FLYING MISSIONS AND TECHNOLOGIES RELATIVE MISSION ANALYSIS FOR PROBA 3: SAFE ORBITS AND CAM Munich, 31 st May 2013 T. V. Peters D. Escorial Presented
More informationMODELLING OF PERTURBATIONS FOR PRECISE ORBIT DETERMINATION
MODELLING OF PERTURBATIONS FOR PRECISE ORBIT DETERMINATION 1 SHEN YU JUN, 2 TAN YAN QUAN, 3 TAN GUOXIAN 1,2,3 Raffles Science Institute, Raffles Institution, 1 Raffles Institution Lane, Singapore E-mail:
More informationLINEARIZED ORBIT COVARIANCE GENERATION AND PROPAGATION ANALYSIS VIA SIMPLE MONTE CARLO SIMULATIONS
LINEARIZED ORBIT COVARIANCE GENERATION AND PROPAGATION ANALYSIS VIA SIMPLE MONTE CARLO SIMULATIONS Chris Sabol, Paul Schumacher AFRL Thomas Sukut USAFA Terry Alfriend Texas A&M Keric Hill PDS Brendan Wright
More informationStatistical methods to address the compliance of GTO with the French Space Operations Act
Statistical methods to address the compliance of GTO with the French Space Operations Act 64 th IAC, 23-27 September 2013, BEIJING, China H.Fraysse and al. Context Space Debris Mitigation is one objective
More informationOperational Implementation of a Pc Uncertainty Construct for Conjunction Assessment Risk Analysis. Lauri K. Newman NASA Goddard Space Flight Center
Operational Implementation of a Pc Uncertainty Construct for Conjunction Assessment Risk Analysis Lauri K. Newman NASA Goddard Space Flight Center Matthew D. Hejduk Astrorum Consulting LLC Lauren C. Johnson
More informationIdentifying Safe Zones for Planetary Satellite Orbiters
AIAA/AAS Astrodynamics Specialist Conference and Exhibit 16-19 August 2004, Providence, Rhode Island AIAA 2004-4862 Identifying Safe Zones for Planetary Satellite Orbiters M.E. Paskowitz and D.J. Scheeres
More informationPlanetary Protection Trajectory Analysis for the Juno Mission
AIAA/AAS Astrodynamics Specialist Conference and Exhibit 18-21 August 2008, Honolulu, Hawaii AIAA 2008-7368 Planetary Protection Trajectory Analysis for the Juno Mission Try Lam 1, Jennie R. Johannesen
More informationLecture D30 - Orbit Transfers
J. Peraire 16.07 Dynamics Fall 004 Version 1.1 Lecture D30 - Orbit Transfers In this lecture, we will consider how to transfer from one orbit, or trajectory, to another. One of the assumptions that we
More informationSPACE DEBRIS CHASER CONSTELLATION
THE VISION SPACE DEBRIS CHASER CONSTELLATION Constellation management Application to SSO debris Inclination between 98 and 99.5 Altitude between 750 km and 850 km Constellation of 38 Cubesats Low-thrust
More informationExhaustive strategy for optical survey of geosynchronous region using TAROT telescopes
Exhaustive strategy for optical survey of geosynchronous region using TAROT telescopes Pascal Richard, Carlos Yanez, Vincent Morand CNES, Toulouse, France Agnès Verzeni CAP GEMINI, Toulouse, France Michel
More informationOn Directional Measurement Representation in Orbit Determination
On Directional Measurement Representation in Orbit Determination Evan M. Ward and Greg Carbott U.S. Naval Research Lab, Washington, DC 20375 Precision Orbit Determination (OD) is often critical for successful
More informationPREDICTING THE ATMOSPHERIC RE-ENTRY OF SPACE DEBRIS THROUGH THE QB50 ENTRYSAT MISSION
PREDICTING THE ATMOSPHERIC RE-ENTRY OF SPACE DEBRIS THROUGH THE QB50 ENTRYSAT MISSION Y. Prevereaud (1), F. Sourgen (2), D. Mimoun (3), A.Gaboriaud (4), J-L. Verant (5), and J-M. Moschetta (6) (1) ONERA
More informationDetection and Conjunction Assessment of Long Encounters
Detection and Conjunction Assessment of Long Encounters F. B a b i k e r V. A b b a s i M D o y o n Detection and Conjunction Assessment of Long Encounters F. B a b i k e r V. A b b a s i M D o y o n Introduction
More informationNumerical Techniques for Computing Orbits
Astrodynamics ASD3: Numerical Techniques for Computing Orbits By Prof. Gregory D. Durgin In this lecture, we present a method for computing the trajectories of satellites around planets or stars with a
More informationSpace Debris Mitigation Activities at ESA
Space Debris Mitigation Activities at ESA Heiner Klinkrad ESA Space Debris Office H. Klinkrad, ESA, Feb 2011, page 1 Space Debris Environment in 2010 4,765 launches and 251 on-orbit break-ups led to 16,200
More informationSpacecraft Orbit Anomaly Representation Using Thrust-Fourier-Coefficients with Orbit Determination Toolbox
Spacecraft Orbit Anomaly Representation Using Thrust-Fourier-Coefficients with Orbit Determination Toolbox Hyun Chul Ko and Daniel J. Scheeres University of Colorado - Boulder, Boulder, CO, USA ABSTRACT
More informationCollision Finding for Particles on Non-linear Paths
Collision Finding for Particles on Non-linear Paths Abstract For many N-body hard-sphere collisional simulations particle trajectories can not be assumed to be linear for the purposes of collision detection.
More informationAccuracy Assessment of SGP4 Orbit Information Conversion into Osculating Elements
Accuracy Assessment of SGP4 Orbit Information Conversion into Osculating Elements Saika Aida (1), Michael Kirschner (2) (1) DLR German Space Operations Center (GSOC), Oberpfaffenhofen, 82234 Weßling, Germany,
More informationProbability of Collision in the Joint Space Operations Center
Probability of Collision in the Joint Space Operations Center For a conjunction 1 between two objects in earth orbit, the Joint Space Operations Center (JSpOC) can compute and report a value commonly referred
More informationYahsat s Approach to GEO Collision Avoidance A Case Stduy with STTW-4
Yahsat s Approach to GEO Collision Avoidance A Case Stduy with STTW-4 John Baker Senior Manager - Flight Dynamics Yahsat Agenda Objective Collision Risk Know your enemy Situational Awareness Yahsat Stationkeeping
More informationAerodynamic Lift and Drag Effects on the Orbital Lifetime Low Earth Orbit (LEO) Satellites
Aerodynamic Lift and Drag Effects on the Orbital Lifetime Low Earth Orbit (LEO) Satellites I. Introduction Carlos L. Pulido Department of Aerospace Engineering Sciences University of Colorado Boulder Abstract
More informationCHANGING INCLINATION OF EARTH SATELLITES USING THE GRAVITY OF THE MOON
CHANGING INCLINATION OF EARTH SATELLITES USING THE GRAVITY OF THE MOON KARLA DE SOUZA TORRES AND A. F. B. A. PRADO Received 3 August 005; Revised 14 April 006; Accepted 18 April 006 We analyze the problem
More informationEutelsat practice & views on long term sustainability. COPUOS/SCST - LTS Workshop Vienna, 14 February 2013
Eutelsat practice & views on long term sustainability COPUOS/SCST - LTS Workshop Vienna, 14 February 2013 Introduction Introduction to Eutelsat Originally set up in 1977 as an intergovernmental organisation
More informationAnalytical Non-Linear Uncertainty Propagation: Theory And Implementation
Analytical Non-Linear Uncertainty Propagation: Theory And Implementation Kohei Fujimoto 1), Daniel J. Scheeres 2), and K. Terry Alfriend 1) May 13, 2013 1) Texas A&M University 2) The University of Colorado
More informationOrekit at the U.S. Naval Research Laboratory. Evan Ward
Orekit at the U.S. Naval Research Laboratory Evan Ward U.S. Naval Research Laboratory Astrodynamics and Navigation Section, Washington DC November 16, 2017 Outline Introduction Geolocation with Orekit
More informationCOMPARISON OF ANGLES ONLY INITIAL ORBIT DETERMINATION ALGORITHMS FOR SPACE DEBRIS CATALOGUING
COMPARISON OF ANGLES ONLY INITIAL ORBIT DETERMINATION ALGORITHMS FOR SPACE DEBRIS CATALOGUING Fran Martinez Fadrique, Alberto Águeda Maté, Joan Jorquera Grau, Jaime Fernández Sánchez, Laura Aivar García
More informationPrevious Lecture. Orbital maneuvers: general framework. Single-impulse maneuver: compatibility conditions
2 / 48 Previous Lecture Orbital maneuvers: general framework Single-impulse maneuver: compatibility conditions closed form expression for the impulsive velocity vector magnitude interpretation coplanar
More informationRAPID GEOSYNCHRONOUS TRANSFER ORBIT ASCENT PLAN GENERATION. Daniel X. Junker (1) Phone: ,
RAPID GEOSYNCHRONOUS TRANSFER ORBIT ASCENT PLAN GENERATION Daniel X. Junker (1) (1) LSE Space GmbH, Argelsrieder Feld 22, 82234 Wessling, Germany, Phone: +49 160 9111 6696, daniel.junker@lsespace.com Abstract:
More informationOptimization of Orbital Transfer of Electrodynamic Tether Satellite by Nonlinear Programming
Optimization of Orbital Transfer of Electrodynamic Tether Satellite by Nonlinear Programming IEPC-2015-299 /ISTS-2015-b-299 Presented at Joint Conference of 30th International Symposium on Space Technology
More informationORBITAL DECAY PREDICTION AND SPACE DEBRIS IMPACT ON NANO-SATELLITES
Journal of Science and Arts Year 16, No. 1(34), pp. 67-76, 2016 ORIGINAL PAPER ORBITAL DECAY PREDICTION AND SPACE DEBRIS IMPACT ON NANO-SATELLITES MOHAMMED CHESSAB MAHDI 1 Manuscript received: 22.02.2016;
More informationAIR FORCE INSTITUTE OF TECHNOLOGY
OPTIMAL ORBITAL COVERAGE OF THEATER OPERATIONS AND TARGETS THESIS Kimberly A. Sugrue, Captain, USAF AFIT/GA/ENY/07-M17 DEPARTMENT OF THE AIR FORCE AIR UNIVERSITY AIR FORCE INSTITUTE OF TECHNOLOGY Wright-Patterson
More informationOrbital Dynamics and Impact Probability Analysis
Orbital Dynamics and Impact Probability Analysis (ISAS/JAXA) 1 Overview This presentation mainly focuses on a following point regarding planetary protection. - How to prove that a mission satisfies the
More informationCHAPTER 3 PERFORMANCE
PERFORMANCE 3.1 Introduction The LM-3A performance figures given in this chapter are based on the following assumptions: Launching from XSLC (Xichang Satellite Launch Center, Sichuan Province, China),
More informationEnd of Life Re-orbiting The Meteosat-5 Experience
End of Life Re-orbiting The Meteosat-5 Experience Milan EUMETSAT, Darmstadt, Germany This article illustrates the orbit maneuver sequence performed during Meteosat- 5 End of Life (EOL) re-orbiting operations
More informationExpanding opportunities for lunar gravity capture
Expanding opportunities for lunar gravity capture Keita Tanaka 1, Mutsuko Morimoto 2, Michihiro Matsumoto 1, Junichiro Kawaguchi 3, 1 The University of Tokyo, Japan, 2 JSPEC/JAXA, Japan, 3 ISAS/JAXA, Japan,
More informationSatellite Constellations for Altimetry
Satellite Constellations for Altimetry By Alain Lamy 1) and Josiane Costeraste 1) 1) CNES - Centre National d'etudes Spatiales, Toulouse, France The framework of the analysis presented in this paper is
More informationUSA Space Debris Environment, Operations, and Modeling Updates
USA Space Debris Environment, Operations, and Modeling Updates Presentation to the 51 st Session of the Scientific and Technical Subcommittee Committee on the Peaceful Uses of Outer Space United Nations
More informationOrbital Anomaly Detection and Application of Space Object
Orbital Anomaly Detection and Application of Space Object Lei CHEN *, Xianzong BAI, Taotao ZHANG College of Aerospace and Materials Engineering, National University of Defense Technology, Changsha 473,China
More informationSEMI-ANALYTICAL COMPUTATION OF PARTIAL DERIVATIVES AND TRANSITION MATRIX USING STELA SOFTWARE
SEMI-ANALYTICAL COMPUTATION OF PARTIAL DERIVATIVES AND TRANSITION MATRIX USING STELA SOFTWARE Vincent Morand, Juan Carlos Dolado-Perez, Hubert Fraysse (1), Florent Deleflie, Jérôme Daquin (2), Cedric Dental
More informationA high order method for orbital conjunctions analysis: Monte Carlo collision probability computation
A high order method for orbital conjunctions analysis: Monte Carlo collision probability computation Alessandro Morselli 1, Dipartimento di Scienze e Tecnologie Aerospaziali, Politecnico di Milano, Via
More informationHybrid (Ion and Chemical) GEO Stationkeeping Maneuver Planning Software
Hybrid (Ion and Chemical) GEO Stationkeeping Maneuver Planning Software J. K. Skipper, D. Racicot, S. Li, R. Provencher and J. Palimaka Telesat Canada, Ottawa, Ontario, Canada. Abstract In the geochronous
More informationCOUPLED OPTIMIZATION OF LAUNCHER AND ALL-ELECTRIC SATELLITE TRAJECTORIES
COUPLED OPTIMIZATION OF LAUNCHER AND ALL-ELECTRIC SATELLITE TRAJECTORIES M. Verlet (1), B. Slama (1), S. Reynaud (1), and M. Cerf (1) (1) Airbus Defence and Space, 66 Route de Verneuil, 78133 Les Mureaux,
More informationA Study of Six Near-Earth Asteroids Summary Introduction
A Study of Six Near-Earth Asteroids John Junkins 1, Puneet Singla 1, Daniele Mortari 1, William Bottke 2, Daniel Durda 2 Summary We consider here 6 Earth-orbit-crossing asteroids (ECAs) as possible targets
More informationRADAR-OPTICAL OBSERVATION MIX
RADAR-OPTICAL OBSERVATION MIX Felix R. Hoots + Deep space satellites, having a period greater than or equal to 225 minutes, can be tracked by either radar or optical sensors. However, in the US Space Surveillance
More informationAS3010: Introduction to Space Technology
AS3010: Introduction to Space Technology L E C T U R E S 8-9 Part B, Lectures 8-9 23 March, 2017 C O N T E N T S In this lecture, we will look at factors that cause an orbit to change over time orbital
More informationDE-ORBITATION STUDIES AND OPERATIONS FOR SPIRALE GTO SATELLITES
DE-ORBITATION STUDIES AND OPERATIONS FOR SPIRALE GTO SATELLITES François BONAVENTURE (1), Slim LOCOCHE (2), Anne-Hélène GICQUEL (3) (1) Tel. (+33) (0)5 62 19 74 27, E-mail. francois.bonaventure@astrium.eads.net
More informationcloseap: GMV S SOLUTION FOR COLLISION RISK ASSESSMENT
closeap: GMV S SOLUTION FOR COLLISION RISK ASSESSMENT D. Escobar Antón (1), A. Pérez Cambriles (1), F. M. Martínez Fadrique (1), A. Águeda Maté (1) (1) GMV S.A., Calle Isaac Newton 11, P.T.M. Tres Cantos,
More informationCHAPTER 3 PERFORMANCE
PERFORMANCE 3.1 Introduction The LM-3B performance figures given in this chapter are based on the following assumptions: Launching from XSLC (Xichang Satellite Launch Center, Sichuan Province, China),
More informationA NONLINEARITY MEASURE FOR ESTIMATION SYSTEMS
AAS 6-135 A NONLINEARITY MEASURE FOR ESTIMATION SYSTEMS Andrew J. Sinclair,JohnE.Hurtado, and John L. Junkins The concept of nonlinearity measures for dynamical systems is extended to estimation systems,
More informationPlasma Physics Prof. V. K. Tripathi Department of Physics Indian Institute of Technology, Delhi
Plasma Physics Prof. V. K. Tripathi Department of Physics Indian Institute of Technology, Delhi Module No. # 01 Lecture No. # 02 Plasma Response to Fields: Fluid Equations Well, friends today I am going
More informationINNOVATIVE STRATEGY FOR Z9 REENTRY
INNOVATIVE STRATEGY FOR Z9 REENTRY Gregor Martens*, Elena Vellutini**, Irene Cruciani* *ELV, Corso Garibaldi, 34 Colleferro (Italy) **Aizoon, Viale Città d Europa 681, 144, Roma (Italy) Abstract Large
More informationAnalytical Method for Space Debris propagation under perturbations in the geostationary ring
Analytical Method for Space Debris propagation under perturbations in the geostationary ring July 21-23, 2016 Berlin, Germany 2nd International Conference and Exhibition on Satellite & Space Missions Daniel
More informationLEAVING CERTIFICATE EXAMINATION 2004 APPLIED MATHEMATICS
Coimisiún na Scrúduithe Stáit State Examinations Commission LEAVING CERTIFICATE EXAMINATION 004 APPLIED MATHEMATICS ORDINARY LEVEL CHIEF EXAMINER S REPORT HIGHER LEVEL CHIEF EXAMINER S REPORT - 1 - CONTENTS
More informationFeedback Control of Spacecraft Rendezvous Maneuvers using Differential Drag
Feedback Control of Spacecraft Rendezvous Maneuvers using Differential Drag D. Pérez 1 and R. Bevilacqua Rensselaer Polytechnic Institute, Troy, New York, 1180 This work presents a feedback control strategy
More informationSingle Lens Lightcurve Modeling
Single Lens Lightcurve Modeling Scott Gaudi The Ohio State University 2011 Sagan Exoplanet Summer Workshop Exploring Exoplanets with Microlensing July 26, 2011 Outline. Basic Equations. Limits. Derivatives.
More informationSatellite Maneuver Detection Using Two-line Element (TLE) Data. Tom Kelecy Boeing LTS, Colorado Springs, CO / Kihei, HI
Satellite Maneuver Detection Using Two-line Element (TLE) Data Tom Kelecy Boeing LTS, Colorado Springs, CO / Kihei, HI Doyle Hall Boeing LTS, Colorado Springs, CO / Kihei, HI Kris Hamada Pacific Defense
More informationOPTIMAL MANEUVERS IN THREE-DIMENSIONAL SWING-BY TRAJECTORIES
OPTIMAL MANEUVERS IN THREE-DIMENSIONAL SWING-BY TRAJECTORIES Gislaine de Felipe and Antonio Fernando Bertachini de Almeida Prado Instituto Nacional de Pesquisas Espaciais - São José dos Campos - SP - 12227-010
More informationExperimental Analysis of Low Earth Orbit Satellites due to Atmospheric Perturbations
Experimental Analysis of Low Earth Orbit Satellites due to Atmospheric Perturbations Aman Saluja #1, Manish Bansal #2, M Raja #3, Mohd Maaz #4 #Aerospace Department, University of Petroleum and Energy
More informationDe Broglie Wavelength and McMahon field theory Abstract Theory: McMahon, C.R. (2015)
Copyright Version: 2 nd February, 2015, updated 24 th September, 2018 Page: 1 of 11 De Broglie Wavelength and McMahon field theory Abstract: Here, I use the concepts of the McMahon field theory to explain
More informationNumerical analysis of the compliance of interplanetary CubeSats with planetary protection requirements
Numerical analysis of the compliance of interplanetary CubeSats with planetary protection requirements Francesca Letizia, Camilla Colombo University of Southampton 5th Interplanetary CubeSat Workshop Session
More informationKeywords: Space Debris, Telescope Network, Geostationary Ring. 1. General Specifications Location
SMARTnet: A SENSOR FOR MONITORING THE GEOSTATIONARY RING Johannes Herzog (1), Hauke Fiedler (2), Martin Weigel (3), Oliver Montenbruck (4), and Thomas Schildknecht (5) (1)(2)(3)(4) DLR German Aerospace
More informationThe CCSDS Navigation Working Group and Space Situational Awareness. GSAW 2013 Session 11C
The CCSDS Navigation Working Group and Space Situational Awareness GSAW 2013 Session 11C David Berry Jet Propulsion Laboratory, California Institute of Technology 20-Mar-2013 Copyright 2013 California
More informationDistributed Coordination and Control of Formation Flying Spacecraft
Distributed Coordination and Control of Formation Flying Spacecraft Michael Tillerson, Louis Breger, and Jonathan P. How MIT Department of Aeronautics and Astronautics {mike t, lbreger, jhow}@mit.edu Abstract
More informationOverview of JAXA Activities on Sustainable Space Development and Space Situational Awareness
Overview of JAXA Activities on Sustainable Space Development and Space Situational Awareness 26 February 2015 Shizuo YAMAMOTO Vice President, JAXA CONTENTS 1. Japan s Space Activities Structure 2. JAXA
More informationIMPLICIT RUNGE-KUTTA METHODS FOR UNCERTAINTY PROPAGATION
IMPLICIT RUNGE-KUTTA METHODS FOR UNCERTAINTY PROPAGATION Jeffrey M. Aristoff, Joshua T. Horwood, and Aubrey B. Poore Numerica Corporation, 485 Hahns Peak Drive, Suite, Loveland CO, 8538 ABSTRACT Accurate
More informationEUROSTAR 3000 INCLINED ORBIT MISSION : LIFETIME OPTIMISATION IN CASE OF INJECTION WITH A LOW INCLINATION
EUROSTAR 3000 INCLINED ORBIT MISSION : LIFETIME OPTIMISATION IN CASE OF INJECTION WITH A LOW INCLINATION Franck Raballand (1), Julie De Lamarzelle (2), François Bonaventure (3), Anne-Hélène Gicquel (4)
More informationAnalytical Non-Linear Uncertainty Propagation: Theory And Implementation
Analytical Non-Linear Uncertainty Propagation: Theory And Implementation Kohei Fujimoto 1 Texas A&M University, College Station, TX, 77843 Daniel J. Scheeres 2 The University of Colorado-Boulder, Boulder,
More informationSession 6: Analytical Approximations for Low Thrust Maneuvers
Session 6: Analytical Approximations for Low Thrust Maneuvers As mentioned in the previous lecture, solving non-keplerian problems in general requires the use of perturbation methods and many are only
More informationAnalysis of Relative Motion of Collocated Geostationary Satellites with Geometric Constraints
www.dlr.de Chart 1 Analysis of Relative Motion of Collocated Geostationary Satellites with Geometric Constraints SFFMT2013, München F. de Bruijn & E. Gill 31 May 2013 www.dlr.de Chart 2 Presentation Outline
More informationDEVIL CHAPTER 6 TEST REVIEW
IB PHYSICS Name: Period: Date: # Marks: 51 DEVIL PHYSICS BADDEST CLASS ON CAMPUS CHAPTER 6 TEST REVIEW 1. A cyclist rides around a circular track at a uniform speed. Which of the following correctly gives
More informationResponsive Imaging Constellations for Support of Geographically Dispersed Theaters
Responsive Imaging Constellations for Support of Geographically Dispersed Theaters Todd J. Mosher Ph.D, 1 Kathryn E. Hamera 2 and Skylar A. Cox 3 MicroSat Systems, Inc., Littleton, Colorado, 80127, University
More informationLong-Term Evolution of High Earth Orbits: Effects of Direct Solar Radiation Pressure and Comparison of Trajectory Propagators
Long-Term Evolution of High Earth Orbits: Effects of Direct Solar Radiation Pressure and Comparison of Trajectory Propagators by L. Anselmo and C. Pardini (Luciano.Anselmo@isti.cnr.it & Carmen.Pardini@isti.cnr.it)
More informationABSTRACT. Nomenclature
ABSTRACT The behavior of two different models of gas-surface interactions is studied using the Direct Simulation Monte Carlo (DSMC) method. The DSMC calculations examine differences in predictions of aerodynamic
More informationAnalysis of Debris from the Collision of the Cosmos 2251 and the Iridium 33 Satellites
Science & Global Security, 18:87 118, 2010 Copyright C Taylor & Francis Group, LLC ISSN: 0892-9882 print / 1547-7800 online DOI: 10.1080/08929882.2010.493078 Analysis of Debris from the Collision of the
More informationSpace Travel on a Shoestring: CubeSat Beyond LEO
Space Travel on a Shoestring: CubeSat Beyond LEO Massimiliano Vasile, Willem van der Weg, Marilena Di Carlo Department of Mechanical and Aerospace Engineering University of Strathclyde, Glasgow 5th Interplanetary
More informationG494. PHYSICS B (ADVANCING PHYSICS) Rise and Fall of the Clockwork Universe ADVANCED GCE. Monday 27 June 2011 Morning. Duration: 1 hour 15 minutes
ADVANCED GCE PHYSICS B (ADVANCING PHYSICS) Rise and Fall of the Clockwork Universe G494 *G410620611* Candidates answer on the question paper. OCR supplied materials: Data, Formulae and Relationships Booklet
More informationEscape Trajectories from Sun Earth Distant Retrograde Orbits
Trans. JSASS Aerospace Tech. Japan Vol. 4, No. ists30, pp. Pd_67-Pd_75, 06 Escape Trajectories from Sun Earth Distant Retrograde Orbits By Yusue OKI ) and Junichiro KAWAGUCHI ) ) Department of Aeronautics
More informationExtending the Patched-Conic Approximation to the Restricted Four-Body Problem
Monografías de la Real Academia de Ciencias de Zaragoza 3, 133 146, (6). Extending the Patched-Conic Approximation to the Restricted Four-Body Problem Thomas R. Reppert Department of Aerospace and Ocean
More informationSpacecraft Collision Probability Estimation for Rendezvous and Proximity Operations
Utah State University DigitalCommons@USU All Graduate Theses and Dissertations Graduate Studies 5-2012 Spacecraft Collision Probability Estimation for Rendezvous and Proximity Operations Michael R. Phillips
More informationADVANCED NAVIGATION STRATEGIES FOR AN ASTEROID SAMPLE RETURN MISSION
AAS 11-499 ADVANCED NAVIGATION STRATEGIES FOR AN ASTEROID SAMPLE RETURN MISSION J. Bauman,* K. Getzandanner, B. Williams,* K. Williams* The proximity operations phases of a sample return mission to an
More informationLecture 1a: Satellite Orbits
Lecture 1a: Satellite Orbits Meteorological Satellite Orbits LEO view GEO view Two main orbits of Met Satellites: 1) Geostationary Orbit (GEO) 1) Low Earth Orbit (LEO) or polar orbits Orbits of meteorological
More informationEVALUATING ORBITS WITH POTENTIAL TO USE SOLAR SAIL FOR STATION-KEEPING MANEUVERS
IAA-AAS-DyCoSS2-14-16-03 EVALUATING ORBITS WITH POTENTIAL TO USE SOLAR SAIL FOR STATION-KEEPING MANEUVERS Thais C. Oliveira, * and Antonio F. B. A. Prado INTRODUCTION The purpose of this paper is to find
More information