THE ROLE OF IMPACTS AND MOMENTUM TRANSFER FOR THE EVOLUTION OF ENVISAT S ATTITUDE STATE

Similar documents
Comparison of ENVISAT s attitude simulation and real optical and SLR observations in order to refine the satellite attitude model

Laser de-spin maneuver for an active debris removal mission - a realistic scenario for Envisat

Estimating the Error in Statistical HAMR Object Populations Resulting from Simplified Radiation Pressure Modeling

OPTICAL OBSERVATIONS OF BRIZ-M FRAGMENTS IN GEO

EXTENSIVE LIGHT CURVE DATABASE OF ASTRONOMICAL INSTITUTE OF THE UNIVERSITY OF BERN

An Optical Survey for Space Debris on Highly Eccentric MEO Orbits

An Optical Survey for Space Debris on Highly Eccentric MEO Orbits

EVOLUTION OF SPACECRAFT ORBITAL MOTION DUE TO HYPERVELOCITY IMPACTS WITH DEBRIS AND METEOROIDS

2 INTRODUCTION 3 ORBIT DETERMINATION

IAC-17-F Attitude State Evolution of Space Debris Determined from Optical Light Curve Observations

Effects of Mitigation Measures on the Space Debris Environment

AIUB EFFORTS TO SURVEY, TRACK, AND CHARACTERIZE SMALL-SIZE OBJECTS AT HIGH ALTITUDES

Spacecraft design indicator for space debris

arxiv: v3 [physics.space-ph] 3 Jan 2018

Defunct Satellites, Rotation Rates and the YORP Effect. Antonella A. Albuja and Daniel J. Scheeres University of Colorado - Boulder, Boulder, CO, USA

INITIAL STUDY ON SMALL DEBRIS IMPACT RISK ASSESSMENT DURING ORBIT TRANSFER TO GEO FOR ALL-ELECTRIC SATELLITE

Orbital Environment Simulator

AS3010: Introduction to Space Technology

Phys 270 Final Exam. Figure 1: Question 1

is acting on a body of mass m = 3.0 kg and changes its velocity from an initial

AN ANALYTICAL SOLUTION TO QUICK-RESPONSE COLLISION AVOIDANCE MANEUVERS IN LOW EARTH ORBIT

AP PHYSICS 1 Learning Objectives Arranged Topically

CONCEPTUAL DESIGN FOR EXPERT COORDINATION CENTRES SUPPORTING OPTICAL AND SLR OBSERVATIONS IN A SST SYSTEM

1 Forces. 2 Energy & Work. GS 104, Exam II Review

Satellite Components & Systems. Dr. Ugur GUVEN Aerospace Engineer (P.hD) Nuclear Science & Technology Engineer (M.Sc)

Improved Space Object Orbit Determination Using CMOS Detectors. J. Silha 1,2

1. A force acts on a particle and displaces it through. The value of x for zero work is 1) 0.5 2) 2 4) 6

CHAPTER 5 FUZZY LOGIC FOR ATTITUDE CONTROL

Study Participants: T.E. Sarris, E.R. Talaat, A. Papayannis, P. Dietrich, M. Daly, X. Chu, J. Penson, A. Vouldis, V. Antakis, G.

RISK INDUCED BY THE UNCATALOGUED SPACE DEBRIS POPULATION IN THE PRESENCE OF LARGE CONSTELLATIONS

Mega Constellations and Space Debris - Impact on the Environment and the Constellation Itself -

King Fahd University of Petroleum and Minerals Department of Physics. Final Exam 041. Answer key - First choice is the correct answer

Jitter and Basic Requirements of the Reaction Wheel Assembly in the Attitude Control System

Design of Attitude Determination and Control Subsystem

PhysicsAndMathsTutor.com 1

Passive Orbital Debris Removal Using Special Density Materials

Northwestern Connecticut Community College Course Syllabus

Course Name: AP Physics. Team Names: Jon Collins. Velocity Acceleration Displacement

Solar Activity Space Debris

Passive Magnetic Attitude Control for CubeSat Spacecraft

EXAMPLE: MODELING THE PT326 PROCESS TRAINER

Part Two: Earlier Material

Sub-millimeter size debris monitoring system with IDEA OSG 1

MULTIPLE CHOICE. Choose the one alternative that best completes the statement or answers the question.

Center of Gravity Pearson Education, Inc.

Physics 180A Test Points

Chapter 9: Momentum Tuesday, September 17, :00 PM

Physics 2A Chapter 8 - Impulse & Momentum Fall 2017

DEVIL PHYSICS THE BADDEST CLASS ON CAMPUS AP PHYSICS

Net Capture Mechanism for Debris Removal Demonstration Mission

GUIDANCE, NAVIGATION, AND CONTROL TECHNIQUES AND TECHNOLOGIES FOR ACTIVE DEBRIS REMOVAL

A DETAILED IMPACT RISK ASSESSMENT OF POSSIBLE PROTECTION ENHANCEMENTS TO TWO LEO SPACECRAFT

spacecraft mass = kg xenon ions speed = m s 1 Fig. 2.1 Calculate the mass of one xenon ion. molar mass of xenon = 0.

Chapter 9: Momentum and Conservation. Newton s Laws applied

Development of Laser Measurement to Space Debris at Shanghai SLR Station

Third Body Perturbation

PHYSICS I RESOURCE SHEET

Scalar product Work Kinetic energy Work energy theorem Potential energy Conservation of energy Power Collisions

AP practice ch 7-8 Multiple Choice

What Happened to BLITS? An Analysis of the 2013 Jan 22 Event

Analysis of Debris from the Collision of the Cosmos 2251 and the Iridium 33 Satellites

Graz in Space Graz SLR System. Daniel Kucharski. IWF / SatGeo

1. (a) Describe the difference between over-expanded, under-expanded and ideallyexpanded

Lecture Module 5: Introduction to Attitude Stabilization and Control

Feedback Control of Spacecraft Rendezvous Maneuvers using Differential Drag

Quaternion-Based Tracking Control Law Design For Tracking Mode

ATTITUDE CONTROL MECHANIZATION TO DE-ORBIT SATELLITES USING SOLAR SAILS

Phys101 Second Major-173 Zero Version Coordinator: Dr. M. Al-Kuhaili Thursday, August 02, 2018 Page: 1. = 159 kw

Plane Motion of Rigid Bodies: Momentum Methods

Optical Surveys for Space Debris in MEO

Northwestern CT Community College Course Syllabus. Course Title: CALCULUS-BASED PHYSICS I with Lab Course #: PHY 221

MAE 142 Homework #5 Due Friday, March 13, 2009

PHYSICS 8A, Lecture 2 Spring 2017 Midterm 2, C. Bordel Thursday, April 6 th, 7pm-9pm

Physics A - PHY 2048C

The basic principle to be used in mechanical systems to derive a mathematical model is Newton s law,

Attitude Control Strategy for HAUSAT-2 with Pitch Bias Momentum System

CHALLENGES RELATED TO DISCOVERY, FOLLOW-UP, AND STUDY OF SMALL HIGH AREA-TO-MASS RATIO OBJECTS AT GEO

Angular Momentum. Objectives CONSERVATION OF ANGULAR MOMENTUM

Swiss Contributions to a Better Understanding of the Space Debris Environment

3. ROTATIONAL MOTION

CONCLUSIONS FROM ESA SPACE DEBRIS TELESCOPE OBSERVATIONS ON SPACE DEBRIS ENVIRONMENT MODELLING

The Torque Rudder: A Novel Semi-Passive Actuator for Small Spacecraft Attitude Control

Pointing Control for Low Altitude Triple Cubesat Space Darts

Tailoring the Observation Scenarios and Data Processing Techniques for Supporting Conjunction Event Assessments

COLLEGE OF FOUNDATION AND GENERAL STUDIES PUTRAJAYA CAMPUS FINAL EXAMINATION TRIMESTER I 2012/2013

NAME. (2) Choose the graph below that represents the velocity vs. time for constant, nonzero acceleration in one dimension.

This Week. 9/5/2018 Physics 214 Fall

Physics for Scientist and Engineers third edition Rotational Motion About a Fixed Axis Problems

Hybrid spacecraft attitude control system

III. Work and Energy

Critical Density of Spacecraft in Low Earth Orbit: Using Fragmentation Data to Evaluate the Stability of the Orbital Debris Environment

PHY218 SPRING 2016 Review for Final Exam: Week 14 Final Review: Chapters 1-11, 13-14

The Evaluation Methods of Collision. Between Space Debris and Spacecraft

CHAPTER 4 CONVENTIONAL CONTROL FOR SATELLITE ATTITUDE

Spacecraft Attitude Dynamics for Undergraduates

SELENE TRANSLUNAR TRAJECTORY AND LUNAR ORBIT INJECTION

Study Guide Solutions

ISS Intergovernmental Agreement

FORMATION FLYING GUIDANCE FOR SPACE DEBRIS OBSERVATION, MANIPULATION AND CAPTURE

The Cosmic Perspective Seventh Edition. Making Sense of the Universe: Understanding Motion, Energy, and Gravity. Chapter 4 Lecture

Transcription:

2017 AMOS Technical Conference, 19-22 September, 2017, Maui, Hawaii, USA THE ROLE OF IMPACTS AND MOMENTUM TRANSFER FOR THE EVOLUTION OF ENVISAT S ATTITUDE STATE Thomas Schildknecht, Jiri Silha, Astronomical Institute, University of Bern, Sidlerstr. 5, CH-3012 Bern, Switzerland Holger Krag Space Debris Office, ESA/ESOC, Germany ABSTRACT The currently proposed space debris remediation measures include the active removal of large objects and just in time collision avoidance by deviating the objects using, e.g., ground-based lasers. These techniques require precise knowledge of the attitude state and state changes of the target objects. In the former case, e.g. to devise methods to capture the target with a tug spacecraft, in the latter, to precisely propagate the orbits of potential collision partners, as disturbing forces like air drag and solar radiation pressure depend on the attitude of the objects. The long-term evolution of the attitude motion is, among many other causes, depending on the effects of possible impacts of debris and meteoroid, while momentum transfer from reaction wheels or other moving internal components may contribute to the root cause of the initial attitude motion. Impacts of small particles like meteoroids and space debris pieces on compact space objects are unavoidable events, which were already observed several times, e.g., on International Space Station, or rather recently on the Sentinel-1A on August 23, 2016. This paper will discuss a detailed analysis of the effects of momentum transfer from the reaction wheels and of debris and meteoroid impacts for the particular case of Envisat. Based on the physical model of Envisat and the MASTER environment model, the likelihood to have an impact-related attitude rate increase in ten years larger than selected threshold rates was determined. 1. INTRODUCTION Recently the determination of attitude states, and in particular rotation rates, of space objects became a topic of interest in the space debris community. This is to be seen in the context of the multitude of techniques which are currently proposed to remove space debris from orbit or to re-orbit them into disposal orbits. The majority of the techniques to remove large objects, which are driving the evolution of the space debris population on the long term, require capturing the target with a robotic arm, a net, a harpoon, or another mechanism. The attitude motion is in all these cases a critical parameter and the maximum tolerable target rotation rate is limited. The attitude motion of an inactive Earth orbiting object may change either due to dissipation of stored internal energy or internal momentum transfer, or because of external forces. Internal energy may be released, e.g., in form of leaking attitude control systems or pressurizes tanks, any type of breakup events, and momentum can be transferred from reaction wheels. External torques may results from the interaction with the gravity field gradients and the magnetic field of the Earth or from impulse transfer when colliding with other objects. The aim of this study is to assess the maximum angular momentum change to be expected for the defunct Envisat spacecraft, on one hand due to angular momentum transfer from its reaction wheels, and, on the other hand, from impacts with space debris objects or micrometeorites. The latter was analyzed once for a single impact event and for a statistical particle flux over a then year period. The study was in particular triggered by the sudden acceleration of Envisat s rotation observed about two months after its failure when Envisat s attitude changed from an almost stable NADIR fixed attitude to a rotation with 0.547 rpm [1], [2], [3] (inertial rotation period ~ 110s; see Figure 1).

Figure 1: Inertial spin rate of Envisat measured by Satellite Laser Ranging observations (from [3]). 2. MOMENTUM ANDD IMPULSE TRANSFER ANALYSIS The interaction between compact spacee objects, such as spacecraft and upper stages on a geocentric orbit, and smalll particles like meteoroids and space debris pieces (hereafter collectively calledd MMOD) is unavoidable and has been observed many times. In case that the impact particle will not cause a catastrophic collision, meaning that t the targett body is not completely disrupted, the impulse of thee impacting particle will bee partially transferred to thee angular mo- mentum of the impacted body and partially to its linear momentum. Such an event e will cause, among other effects, a sudden change of the attitude state of the impacted body. The size and direction of the resulting angular momentumm change will strongly depend on the impact particle mass, the impact velocity, and on the impact geometry, i.e. the loca- tion and impact angle with respect to the target body center of mass. In our investigation we used a simple model for the Envisat satellite consisting of a bus (a rectangular cube with dimensions 4m x 4m x 10m, 7512 kg) and solar panel (a thin plate with dimensionss 14m x 5m x 0m, 388 kg) ) (see Figure 2). By using the parallel axis theorem, we determinedd the momenta of inertia for this model. The values obtained for,, and for the given principal axis of inertia are listed in the Table 1. [kg.m 2 ] [kg.m 2 ] [kg.m 2 ] 135284..9 21163.7 134476.5 Table 1: Momenta of inertia for the Envisat model consisting of a solid rectangular cube (satellite bus) and a thin plate (solar panel). 3. SINGLE IMPACT EVENT In this case we investigated the probability for a collision between Envisat and a particle withh a mass within the interval 0.0016 g to 1 kg. In addition we derivedd the resultingg angular momenta that particle would transfer to the satellite model during the collision by assuming given scenario. A representative synthetic population p has been generated by using ESA MASTER tool [5] and the orbital elements of Envisat from October 2014. For every generated particle which was defined by the mass, impact velocity, and impact angle we calculated the resulting received angular momentum, and consequently for a given case also the angular a velocity change and thee resulting spin period P. Three different cases a, b or c, corresponding to a rotation around the principal axes, or, respectively, r were analyzed. These casess are shown in Figure 2. The impacting particle hits thee target at four different locations (greenn points defined by position vector displayed as blue arrow) to achieve that the target would start rotating either around the principal axis (Figure 2 a), (Figure 2 b) or (Figuree 2 c and Figure 2 d). 1,000 synthetic particles per mass interval were generated and the angular momentum that these particles would transfer to Envisat was computed. The resulting probability for

different angular momentum changes around a the three principal axes over a period of ten years are given in Figure 3. The corresponding angular velocities are provided in Figure 4. Figure 2: Visualization of four different collision scenarios to determine the transferred angular momentumm and the re- with targett in the position marked by the green point and blue position vector.. sulting angular velocity change of the target after a collision. The impactor is showed as a redd point colliding These resultss clearly show that there is i an very small probability for an impact which would result in a substantial change of the angular momentum and thus t in a substantial change of the spin rate around any axis of Envisat. There is, e.g., only a 0.9% probability that the bus would be hit by a particle of 6.31 g which in the worst case (casee (a) Figure 2) would change the angular velocity by 0.15 deg/s, corresponding to spin period of o ~40 min. An impact of an a object with a mass of ~100 g is required to spin up Envisat fromm a stable attitude to about 0.4 rpm, which would be comparable c to the angular momentum change observed two monthss after the contact with the spacecraft s hass been lost (Figure 1). Such an event is predicted with a probability of 0.24% over a time span of 10 yearss or with a probability of 0.004% over a period of two months. As a conclusionn we may clearly excludedd an impact event for the explanation of the observed spin rate change of Envisat as a collision with an object of this mass would have damaged the spacecraft severely but remote sensing observation show no evidence of any damage on Envisat. Figure 3: Probability of collision during ten years of exposure (vertical axis) versus received angular momentum assum- ing a perfect inelastic collision for three different cases causing maximum rotation around the principal axes (red), (blue) and (green). The figures show the result of f 1,000 simulations using the MASTER MMOD population (mass and impact velocity distributions) generated by the MASTER program assumingg Envisat orbit and 10 yearss of exposure.

For each impact geometry (related also to given moment of inertia, see Error! Reference source not found.) and total angular momentum we can calculate thee obtained rotation velocity of the model. In Figure 4 can be seen received angu- lar velocity (lower horizontal axis) and resulting spinn period (upper horizontal axis) a as a function of probability of event (left vertical axis). Figure 4: Probability of collision during ten years of exposure (left vertical axis) ) versus received angular velocity (lower horizontal axis) and resulting spin period (upper horizontal axis) assuming a perfect inelastic collision for three t different cases causing maximum rotation around the principal axes (red), (blue) and (green). Thee figures show the result of 1,000 simulations using the MASTER MMOD population (masss and impact velocity v distributions) generated by the MASTER program assuming Envisat orbit and 10 years of exposure. 4. STATISTICA AL PARTICLE FLUX IMPACTS To investigate the effect that would be caused by a flux of particles impacting Envisat E duringg ten years of exposure, the ESA MASTER model MMOD population of particles heavier than 10-15 kg was used. The properties of this population are given in Figure 5. The MASTER tool was then used to predict impacts of this MMOD population with Envisat over a 10 year period. In total about 5 10 6 particles were impacting Envisat. We mayy now again consider a worst case situa- tion and assume that all impacts occurred at one of the positions (a), (b), or (c) defined in Figure 2. The total amount of angular momentum received from this flux of particles over 10 years is relatively small, namely 29.4 kg m 2 s -1 for case (a), 12.2 kg m 2 s -1 for case (b), and 97..6 kg m 2 s -1 for case (c), respectively. Some characteristics of these impacts are given in Figure 6.The results of these cumulative c momentum transfers are equivalent to an impact of a single s particle with a mass equal to ~0.45 g, ~0.45g and ~0.55g, respectively, resulting in angular a velocity changes off 0.013 deg/s, 0.033 deg/s, and 0.042 deg/s for the three cases. Figure 5: Cumulative distributions forr the mass (onn the left) and size (on the right) of the MMOD particles. The flux was generated by MASTER for the Envisat orbit (reference to 21102014). The simulation s covered 10 years from 01-11-- 2014 to 01-11-2024 using the MASTERR MMOD population in thee mass range 10-15 g - 1 kg.

Figure 6: Impact velocity distribution (left) and impact azimuth and elevation distribution d (right) of the MMOD particovered 10 cles. The flux was generated by MASTER for the Envisat orbit (reference to 21102014). 2 The simulation years from 01-11-2014 to 01-11-2024 using the MASTER MMOD population in the mass range 10-15 g - 1 kg. 5. MOMENTUM TRANSFER FROM REACTION WHEELSS There are in total five reaction wheels installed on Envisat [1]. The configuration of the reaction wheels is shown in Figure 7. One reaction wheel is aligned along the axis (ESA reference system) or alongg the axis (reference sys- tem used in this paper). Two reaction wheels are aligned along the axis (along the axis) and two reaction wheels are aligned along the (along the axis) a Figure 7: Distribution and direction off the reaction wheels installed on Envisat in the satellite body coordinate system and in the ESA satellite body coordinate system. Each of the five reaction wheels has a capacity c of 400 N m s (40 kg m -2 /s) with a torque range of ±0.45 Nm. The maxi- mum angular momentum which can be transferred from the wheels too the body amounts thus to t = = 80 N m s and for it is = 40 N m s. For this analysis we used slightly different values forr moment of inertia of the principal p axes of inertia of Envisat whichh were determined by ESA (Table 2). [kg m 2 ] [kg m 2 ] [kg.m 2 ] 17023.3 124825.7 129112.2 Table 2: Momenta of inertia for the Envisat spacecraft determined by ESA. The maximum change in angular velocity for the three axes, if the momentum from fully loaded reaction wheel would be transferred to the body, is given in Table 3.

2017 AMOS Technical Conference, 19-22 September, 2017, Maui, Hawaii, USA Principal axis No. of reaction wheels aligned with axis Reaction wheel capacity [kg m 2 s 1 ] ENVISAT moment of inertia [kg m 2 ] Maximum angular velocity change [deg s 1 ] Corresponding rotational period [hour] 1 40 17023.3 0.1346 0.743 2 40 + 40 124825.7 0.0367 2.723 2 40 + 40 129112.2 0.0355 2.817 Table 3: The maximum angular velocities and corresponding rotational periods of Envisat which can be obtained from reaction wheel(s) assuming different principal axis of inertia. These results show that a momentum transfer from the reaction wheels cannot explain the observed spin-up of Envisat as the maximum change in angular velocity from this transfer is almost two orders of magnitudes smaller than the observed one (period of ~120s). 6. SUMMARY Effects of debris and meteoroid impacts and of momentum transfer from the reaction wheels for the particular case of Envisat were analyzed. Based on the physical model of Envisat and the MASTER environment model, the likelihood to have an impact-related spin rate change in ten years was determined. There is an very small probability for a single impact which would result in a substantial change of the angular momentum and thus in a substantial change of the spin rate around any axis of Envisat over a ten year time period. A particle of ~100 g is required to spin up Envisat from a stable attitude to about 0.4 rpm, which would be comparable to the angular velocity change observed two months after the contact with the spacecraft has been lost. The probability for such an event is less than 0.25%. Simulations showed that the accumulated momentum transfer from the entire MASTER MMOD population is also minor and would in with worst case assumption result in angular velocity change of 0.04 deg/s. The maximum momentum which could be transfer from the Envisat s reaction wheels to its body, assuming that the wheels were fully loaded, would result in an angular velocity change of 0.14 deg/s, 0.037 deg/s, and 0.036 deg/s, respectively, for the three main axes of inertia. Again, this is by far not sufficient to explain the observed spin-up of Envisat after its failure. 7. ACKNOWLEDGEMENTS This work was performed in the context of the ESA study Debris Attitude Motion Measurements and Modelling, (ESA/ESOC 4000112447/14/D/SR). 8. REFERENCES [1] Sommer, S., Rosebrock, J., Cerutti-Maori, D., & Leushacke, L. (2017). Temporal analysis of ENVISAT s rotational motion. In 7th European Conference on Space Debris, ESA/ESOC, Darmstadt/Germany, 2017. [2] D. Kucharski, G. Kirchner, F. Koidl, C. Fan, R. Carman, C. Moore, A. Dmytrotsa, M. Ploner, G. Bianco, M. Medvedskij, A. Makeyev, G. Appleby, M. Suzuki, J. M. Torre, Z. Zhongping, L. Grunwaldt, and Q. Feng, Attitude and spin period of space debris envisat measured by satellite laser ranging, IEEE T. Geosci. Remote, vol. 52, no. 12, pp. 7651 7657, 2014. [3] Kirchner, G., Steindorfer, M., Wang, P., Koidl, F., Kucharski, D., Silha, J., Schildknecht, T., Krag, H. & Flohrer, T. Determination of Attitude and Attitude Motion of Space Debris, Using Laser Ranging and Single-Photon Light Curve Data. In 7th European Conference on Space Debris, ESA/ESOC, Darmstadt/Germany, 2017. [4] Flegel, S., Gelhaus, J., Möckel, M., Wiedemann, C., Kempf, D., Krag, H. and Vörsmann, P. Maintenance of the

2017 AMOS Technical Conference, 19-22 September, 2017, Maui, Hawaii, USA ESA MASTER Model Final Report. Technical Report ESA Contract Number: 21705/08/D/HK, Institute of Aerospace Systems (ILR), June 2011. [5] Bargellini, P.; Garcia Matatoros, M. A.; Ventimiglia, L.; Suen, D., Envisat Attitude and Orbit Control In-Orbit Performance: AN Operational View, 6th International ESA Conference on Guidance, Navigation and Control Systems, held 17-20 October 2005 in Loutraki, Greece. Edited by D. Danesy. ESA SP-606. European Space Agency, 2006. Published on CDROM., id.52.1