HybridSail. Hybrid Solar Sails for Active Debris Removal Final Report


 Bryce Porter
 6 years ago
 Views:
Transcription
1 HybridSail Hybrid Solar Sails for Active Debris Removal Final Report Authors: Lourens Visagie (1), Theodoros Theodorou (1) Affiliation: 1. Surrey Space Centre  University of Surrey ACT Researchers: Leopold Summerer Date: 27 June 211 Contacts: Vaios Lappas Tel: +44 () Fax: +44 () Leopold Summerer (Technical Officer) Tel: +31 () Fax: +31 () Available on the ACT website Ariadna ID: b Ariadna study type: Standard Contract Number: /1/NL/CBi
2 Abstract The historical practice of abandoning spacecraft and upper stages at the end of mission life has resulted in a polluted environment in some earth orbits. The amount of objects orbiting the Earth poses a threat to safe operations in space. Studies have shown that in order to have a sustainable environment in low Earth orbit, commonly adopted mitigation guidelines should be followed (the Inter Agency Space Debris Coordination Committee has proposed a set of debris mitigation guidelines and these have since been endorsed by the United Nations) as well as Active Debris Removal (ADR). HybridSail is a proposed concept for a scalable de orbiting spacecraft that makes use of a deployable drag sail membrane and deployable electrostatic tethers to accelerate orbital decay. The HybridSail concept consists of deployable sail and tethers, stowed into a nano satellite package. The nano satellite, deployed from a mothership or from a launch vehicle will home in towards the selected piece of space debris using a small thruster propulsion firing and magnetic attitude control system to dock on the debris. The docking process relies on either a net that envelops the debris piece or a mechanical grabber that can grip the target object. Once docked, the magnetic attitude control system on the HybridSail craft will de tumble the debris object so that it is in a stable roll, pitch and yaw orientation. At this point the sail and tethers will deploy and the combined drag force and electrostatic tether force will proceed to de orbit the joined satellites. This report gives results on the de orbiting capabilities of such a concept, together with an analysis of the effects that influence the de orbit times. We also report on the conceptual design of the satellite and docking system and give results from attitude control simulations. This study was performed within the ESA Ariadna scheme.
3 Contents 1 Background Orbital debris Mitigation Measures Active Debris Removal (ADR) Hybrid Solar Sail for Active Debris Removal CubeSail ADR using HybridSail Theory Aerodynamic drag Electrostatic tether Solar Radiation Pressure Satellite orbit propagation Models SGP Cowell Integration Semi analytic Satellite Theory (SST) Validation of SST & Cowell methods Attitude Coordinate system definitions Attitude representation Attitude propagation Satellite de orbit analysis Drag sail de orbiting in LEO Solar activity Sail sizing Orbit Eccentricity Electrostatic de orbiting in LEO Solar activity Tether Parameters Combined sail tether de orbiting in LEO Collision probability Solar sailing orbit manoeuvres... 33
4 4.5.1 Solar sail de orbiting in GEO Conclusion HybridSail design Sail design Booms Sail membrane and folding Sail shape Electrodynamic tethers Conclusions Docking system design Nets Robotic arms Tethers Unlimited Inc (TUI) The Grasp technology Grabber Comparison/Conclusion Attitude Control De tumbling using magnetic control De tumbling controller Torquer sizing Results Passive stabilisation due to sail shape Conclusions Future work References... 57
5 1 Background 1.1 Orbital debris The historical practice of abandoning spacecraft and upper stages at the end of mission life has resulted in a polluted environment in low earth orbit. The amount of objects orbiting the Earth poses a threat to safe operations in space. A collision in space does not only destroy the colliding objects, but also creates a cloud of additional debris pieces that can further harm other intact satellites (as demonstrated by the collision between Russia's Cosmos 2251 and a commercial Iridium satellite in 29). If additional debris is generated during such a break up at higher altitudes, it is not naturally removed by drag and can cause the on set of a cascading effect during which the particles would engage in chain reactions until they are ground to a limiting size. Altitudes with such a critical particle concentration are already suspected to exist near 8, 1 and 15 km [1] The majority of debris objects occur in low Earth orbit (LEO), at altitudes below 2km. The figures below shows the distribution of objects by orbit and object type, taken from the 26 catalogue (this was prior to the Cosmos/Iridium collision). Highly MEO eccentric 4% orbits 11% GEO 1% Other 8% LEO 67% Fragments 4% Operafonal payloads 6% Non operafonal payloads 26% Mission related objects 1% Rocket bodies 18% Figure 1 Distribution of Earth orbiting objects by orbit type (left) and object type (right) [2] It can be seen that only a small portion of the objects orbiting the Earth are operational, and that most objects are concentrated around LEO. The trend is set to continue, as most of the planned launches for the current decade (21 to 219) are destined for LEO 41 of the 122 planned payloads will be launched to LEO. This number excludes rocket bodies, and satellites with a mass below 5kg. HybridSail : Hybrid Solar Sails for Active Debris Removal Page 1
6 Other (GTO, MEO, HEO, deep space) 36% LEO 34% GEO (communica fons) 17% LEO/MEO (communica fon) 13% Figure 2 Satellites to be launched between 21 and 219 by orbit type [3] The NASA orbital debris evolutionary model, LEGEND, predicts moderate linear increases in the MEO and GEO ranges, using a repetition of past launch cycles and assuming no post mission disposal was carried out. But for LEO ranges there is a rapid non linear increase in the orbital population as a result of new launches and collisions [4]. Their analysis has shown that 2 key strategies have to be implemented to stabilize the future LEO environment [5]: 1. A good implementation of the commonly adopted mitigation measures 2. Active removal of existing debris at a rate of at least 5 significant objects per year 1.2 Mitigation Measures The current accepted mitigation guidelines are based on the Inter Agency Space Debris Coordination Committee (IADC) Space Debris Mitigation Guidelines [6]. The IADC is an international forum for the coordination of activities related to space debris and is comprised of the space agencies of 1 countries including ESA. In February 27, the Scientific and Technical Subcommittee (STSC) of the United Nations' Committee on the Peaceful Uses of Outer Space (COPUOS) adopted a consensus set of space debris mitigation guidelines based on that of the IADC. The guidelines were accepted by the COPUOS and endorsed by the United Nations in January 28. The guidelines relevant to LEO satellites and to this study can be summarized as [6]: 1. Limit debris released during normal operations 2. Minimise the potential for on orbit break up a. Minimise the potential for post mission break ups resulting from stored energy b. Minimise the potential for break ups during operational phase 3. Post mission disposal: A spacecraft or orbital stage should be left in an orbit in which atmospheric drag will limit the orbital lifetime after completion of operation. 25 years have been found to be a reasonable and appropriate limit. 1.3 Active Debris Removal (ADR) Active Debris Removal is the action of identifying and disposing of existing orbital debris by moving it into a disposal orbit (one where atmospheric drag will limit the lifetime or a graveyard orbit where it will not interfere with other objects). The problem is quite complex because it requires either experimental contactless technologies, or complicated manoeuvring and docking operations, in addition to efficient propulsion with re usability requirements. HybridSail : Hybrid Solar Sails for Active Debris Removal Page 2
7 Perhaps the easier part of ADR is the identification of priority pieces of debris. Some debris objects pose a greater collision risk than others. These objects are better candidates for ADR than others. The NASA LEGEND model suggests a selection criterion based on the product of the object mass and the collision probability, where the latter is a function of the object cross section area and the object density of the orbit in which it occurs. The US Space Surveillance Network (SSN) maintains a publicly available catalogue of tracked objects in Earth orbit. The information on each object is in the form of orbital elements (two line element sets) that describes the orbit. Information was extracted from the catalogue (on 12 Dec 21) to obtain the distribution plots below. 8 AlPtude DistribuPon in LEO AlPtude/InclinaPon DistribuPon in LEO No. of objects No. of objects AlPtude (km) 15 1 AlPtude (km) 5 14,4 28,8 43,2 57, ,4 1,8 115,2 129, ,4 172,8 InclinaPon (deg) Figure 3 Distribution of objects in low Earth orbit It can be seen that peak object densities occur around 8km and 15km altitude, and also at 97.2 inclinations (at an altitude of 86 km) and 73.8 (72km). Unfortunately the SSN catalogue does not contain information about the mass or cross section area of objects, and no publicly available source exists for this information. Without the satellite mass and area knowledge, the best that the selection criterion can achieve is to identify objects in the higher density orbits. For the purposes of this study, the ESA Advanced Concepts Team has provided a list of priority debris objects. All the objects in the list have perigee altitude above 7km, apogee altitude below 9km, eccentricity below.1 and mass exceeding 5kg. Average cross section area (m 2 ) Mass (kg) Table 11 Priority list of debris objects for active debris removal Semi major axis Eccentricity Inclination (deg) Right ascension of the ascending node (deg) Argument of perigee (deg) HybridSail : Hybrid Solar Sails for Active Debris Removal Page 3
8 E E E HybridSail : Hybrid Solar Sails for Active Debris Removal Page 4
9 Hybrid Solar Sail for Active Debris Removal HybridSail is a concept that combines a large deployable reflective sail with embedded tethers and docking capability to yield a satellite that can be used to de orbit debris pieces using a combination of aerodynamic drag, solar radiation pressure and electrostatic braking. In drag sail mode the deployed membranes will be used to increase the area of the spacecraft that will interact with the atmospheric particles, causing an increased drag and a faster de orbiting. The size of the sail required to successfully de orbit a satellite will depend on the mass of the spacecraft and its initial orbit. The embedded tethers will operate on the principle of an electric sail [16]. The de orbiting is accelerated by making use of momentum exchange with the ionospheric plasma stream. Electrostatic braking is a result of maintaining a negative charge in the tethers. The negative charge in the tether will deflect electrons in the ionosphere plasma, thereby exchanging momentum with it. The momentum loss in the satellite will result in a lower velocity and consequently a faster de orbiting time. The concept will build on an existing solar sail/drag sail project currently underway at the Surrey Space Centre, CubeSail, as well as the research carried out on plasma brake de orbiting [16]. 2.1 CubeSail The objectives of the CubeSail project is to develop and flight test a deployable sail to demonstrate solar radiation pressure and drag de orbiting. This project is a pathfinder mission that hopes to create options for de orbiting of satellites and paves the way for solar sail propulsion. The stowed satellite will have a CubeSat 3U form factor (1x1x34cm) and the square deployable sail has a size of 5m x 5m. The figure below shows this concept in its stowed configuration ready to be launched and in the deployed configuration. HybridSail : Hybrid Solar Sails for Active Debris Removal Page 5
10 Figure 4 CubeSail stowed and deployed configurations A scaled version of the deployment system has already been constructed at the Surrey Space Centre. This is a 1.7x1.7 meter breadboard model used to test different deployment configurations and folding patterns. Figure 5 CubeSail 1.7m x 1.7m prototype Launch of this spacecraft is planned for late 211 and the project is funded by EADS/Astrium. 2.2 ADR using HybridSail The HybridSail concept consists of deployable sail and tethers, stowed into a nano satellite package. The nano satellite, deployed from a mothership or from a launch vehicle will home in towards the selected piece of space debris. The nano satellite will use a small thruster propulsion firing and use its magnetic attitude control system to dock on the debris. The docking process relies on either a net that envelops the debris piece or a mechanical grabber that can grip the target object. The docking method is detailing in chapter 6. HybridSail : Hybrid Solar Sails for Active Debris Removal Page 6
11 Once docked, the magnetic attitude control system on the HybridSail craft will de tumble the debris object so that it is in a stable roll, pitch and yaw orientation. At this point the sail and tethers will deploy and the combined drag force and electrostatic tether force will proceed to de orbit the joined debris object HybridSail. Figure 6 HybridSail docking and deployment 3 Theory 3.1 Aerodynamic drag The presence of gas molecules in low Earth orbit results in both a disturbance force and a disturbance torque. The first will perturb the orbit of the satellite and cause the semi major axis to decay over time, and the latter will influence the attitude of the satellite. This study considers two different models for the aerodynamic drag force. The first is the widely accepted formula for atmospheric drag "#$ = 1 2 " "#. "# "# "# 1 HybridSail : Hybrid Solar Sails for Active Debris Removal Page 7
12 Where is the density of the atmosphere, is the area presented to the approaching molecules, is the coefficient of drag (typically varies between 2. and 3. for a satellite) and "# is the velocity of the satellite relative to the atmosphere. The above formula is an approximation that is often used when estimating orbit lifetime. It assumes hypervelocity, continuum flow [11]. It is also common to use the average cross section area for a rotating satellite. For this study we will only consider the area of the sail membrane because the satellite body will be negligible compared to it. We also assume that the upper atmosphere rotates along with the Earth to give "# = + "# = + 2 Where is the velocity of the satellite, is the position vector and = is the angular velocity vector of the Earth. The second model for atmospheric drag will be used when studying the attitude of the satellite. The disturbing torque due to aerodynamic drag is given by the formula [12]: "#$ = "# "# + "# + 2 cos 3 Where "# = "# "# is the local atmospheric unit vector is the centre of pressure vector relative to the satellite centre of mass is the normal vector of the sail (see explanation below) and are the normal and tangential accommodation coefficients "# is the ratio of molecular exit velocity to local atmospheric velocity cos = "# is the cosine of the angle between the local velocity vector and the sail normal vector The projected area,, is calculated from Where is the actual area of the sail. = cos 4 To allow for situations where the sail flips over, the sail normal vector is chosen such that cos = "# is greater or equal to zero. If the dot product is negative, the negative of the sail normal is taken to ensure that the angle will always be 9 degrees or smaller. We will use the assumptions from [13]: = =.8 and "# =.5 The atmospheric density is an important factor when evaluating drag. Various models exist for calculating atmospheric density with varying consideration of contributing factors. The simplest HybridSail : Hybrid Solar Sails for Active Debris Removal Page 8
13 static model assumes that the density varies exponentially with altitude. More sophisticated models make use of indices of solar and geomagnetic activity (F 1.7, a p and K p ) to determine the time dependant density. It was found that the input data to these models affects the accuracy just as much as the choice of the actual model [11]. It is therefore important to sufficiently estimate the solar and geomagnetic activity in order to model the atmosphere. Because propagations are carried out into the future and real data is not available in this case, the problem becomes estimating what the most likely trend will be. The atmospheric model that was used in this study is the US Naval Research Laboratory NRLMSISE model. It is particularly suited for space applications. The inputs to the model are: Current time Position F 1.7 solar flux index a p geomagnetic index Historic values for the F 1.7 and a p indices are available, but for future dates these values have to be estimated. The estimation method that was followed in this study is described in [14]. The historic solar data was used to construct a pattern of mean values over one cycle for both indices. These mean cycles are then used for future dates, with a regression method to estimate the remainder of the current cycle. F1.7 Solar Flux  Mean Solar Cycle Ap geomagnetic index  Mean Solar Cycle time (years) time (years) Figure 7 Mean solar flux and geomagnetic indices for use in atmospheric density estimation 3.2 Electrostatic tether Tethers have proven to be a popular research topic for satellite de orbiting. Different physical phenomena are exploited by the various tether types. Electrodynamic tethers make use of electromagnetic principles and can act either as electric power generators or as thrusters [15]. It relies on the motion of the satellite through a magnetic field. The tether principle exploited in this study relies on the momentum exchange between a charged tether and the ionospheric plasma. This principle is described in [16]. The braking force that the tether will experience is given by: HybridSail : Hybrid Solar Sails for Active Debris Removal Page 9
14 "#$%&'$($)# = 1.72 "# "# "# "# 5 Where "# = "# is the flow dynamic pressure is the tether length is the ionosphere electron number density is the ionospheric mean ion mass "# is the velocity of the spacecraft relative to the plasma is the tether voltage is the permittivity of vacuum (8.854 x 112 A 2 s 4 kg  1 m  3 ) is the charge of an electron (1.69 x 119 C) The force calculation requires a model of the ionosphere. The International Reference Ionosphere 27 (IRI 27) model was used in this study to obtain the ion and electron densities. 3.3 Solar Radiation Pressure Solar Radiation Pressure (SRP) is normally considered a disturbing force on the orbit of a satellite. This force can also work to our advantage when using the sail as a solar sail to perform orbit manoeuvres. The solar radiation pressure forces are due to the momentum exchange of photons with the sail. A portion of the photons that impact the sail are absorbed, and for the reflected photons, a portion is specularly reflected and the rest is diffusely reflected. Specular reflection is mirror like reflection where the reflected photon always travels on a direction that is determined by the incoming direction. Diffuse reflection is reflection from a rough surface where the incident photon is seemingly reflected at a number of angles. The amount of photons that are absorbed, specularly and diffusely reflected are expressed as fractions. + + = 1 6 From [17], the SRP force acting on a flat Lambertian surface at a distance of = 1 AU from the sun is given by = " Where = x 16 N/m 2 is the nominal solar radiation pressure constant 1 AU from the sun, is the surface area of the sail, is the normal vector to the sail surface and is a unit vector pointing from the sun to the surface. HybridSail : Hybrid Solar Sails for Active Debris Removal Page 1
15 The SRP force for a distance other than 1 AU is obtained by correcting the solar radiation pressure constant,, for the new solar distance,. This results in a SRP force of "#$% = " We will be using equation (7) throughout this study because the orbits considered here all have an approximate sun distance of 1AU. Figure 8 Solar Radiation Pressure normal and tangential components Typical values for a reflective sail membrane are (from the optical properties described in [18]) =.17 =.83 = 9 The disturbance torque caused by solar radiation pressure on the sail is given by "#$% = "#$% 1 Where is the position vector from the satellite centre of mass to the centre of pressure of the sail. 3.4 Satellite orbit propagation In order to examine the effectiveness of Hybridsail as a de orbiting device, a means is required to accurately propagate the orbit of the satellite taking all the disturbing forces into account. The equation that describes the motion of a satellite in an inertial coordinate frame is [7] Where = + + V 11 is the position vector from the centre of mass of the central body to the satellite = is the acceleration vector, = " is the gravitational parameter for the central body, is the acceleration due to non conservative perturbing forces such as aerodynamic drag, solar radiation pressure (SRP) and electromagnetic forces and is a potential function that describes the acceleration due to conservative perturbing forces (third body point mass and central body spherical harmonics). HybridSail : Hybrid Solar Sails for Active Debris Removal Page 11
16 Orbit propagation is carried out in the Earth centred Inertial (ECI) coordinate frame, referenced to the J2 epoch. Although the ECI frame is strictly not an inertial frame (the Earth is accelerating as it orbits around the Sun) equation (11) is still valid, as long as the third body potential function includes the acceleration of the ECI frame (the second term in equation (14)). Orbit propagation methods fall into three categories [8]: 1. Special perturbations orbit propagators based on numerical integration of the equations of motion 2. General perturbations orbit propagators that use analytic theories 3. Semi analytic Satellite Theory (SST) propagators that use numerical integration to propagate mean element values and analytical models for the short term periodic variations The first method yields the most accurate results but requires small integration steps (typically 1 steps per revolution). Adaptive step sizes are usually employed when performing the numeric integration with the step size determined by the error bounds on the state and state derivatives. General perturbations propagators like the SGP4 orbit model can give quick results but suffers from inaccuracies due to inherent approximations [8]. The results also decay quickly when propagated far into the future relative to the time when the defining element sets were generated, although this is as a result of inaccuracies contained in the element set, not the propagation technique. Semi analytic propagators can approach the accuracy of special perturbations propagators with the advantage of integrating over much larger time steps (typically 1 day or longer). Orbit propagation requires mathematical models for the perturbing effects, such as the central body gravitational potential, third body effects and non conservative forces like drag and SRP. Two different orbit propagators were implemented for this study. The first one performs numerical integration on the motion equations in Cartesian form (Cowell integration). The second uses Semi analytical Satellite Theory to integrate the mean orbital elements. A reference implementation of the SGP4 orbit propagator was used to validate both methods Models Earth gravity potential The Earth gravity potential is modelled using spherical harmonics [9]. "#$ = " 1 + " sin " cos " + " sin " 12 Where is the gravitational potential function (m 2 /s 2 ), GM is the gravitational parameter of the Earth, is the distance from the satellite to the Earth s centre of mass, is the semi major axis of the WGS84 ellipsoid, " and " are the normalized gravitational constants and,, is the position of the satellite in spherical coordinate representation. = = atan2, 13 = tan HybridSail : Hybrid Solar Sails for Active Debris Removal Page 12
17 The EGM96 (Earth Gravity Model 1996) coefficients were used in the model, up to order and degree Sun and Moon (third body) gravity potential The sun and moon are modelled as point masses. In this case, the disturbing potential becomes [7] "# "#$ =. 14 With equal to the third body gravitational constant, the position vector from the Earth centre of mass to the satellite and the length of this vector, the vector from the Earth centre of mass to the third body and the length of this vector. The position of the sun is calculated by assuming an unperturbed elliptical motion of the Earth around the Sun, as described in [1]. The position of the moon relative to the Earth centre of mass is calculated using periodic terms, also described in [1] Nonconservative forces The non conservative forces are described by the equations in section 3.1 to SGP4 SGP4 collectively refers to 5 mathematical models (SGP, SGP4, SDP4, SGP8 and SDP8) for predicting the perturbation effects on satellites in order to produce analytic solutions for satellite position and velocity vectors [19]. They are compatible with the two line element sets produced by NORAD for all tracked Earth orbiting objects. The Simplified General Perturbations (SGP) models apply to satellites with an orbital period of 225 minutes or less. The position error of the SGP4 model is in the order of 1km at epoch and it grows by 1 to 3km per day. NORAD releases updated two line elements for all tracked objects for this reason. Simplified Deep space Perturbations (SDP) models are used for satellites with higher orbits (period greater than 225 minutes) and assume a simplified drag model and add third body (sun and moon) gravity effects. Because of the availability of up to date two line elements, the SGP4 models have become the standard for coarse satellite position determination. The theory for the SGP4 models is contained in reference [19]. For this study, the reference code implementation from the Center for Space Standards & Innovation (CSSI) was used Cowell Integration The state of the satellite is given by = 15 with = the position vector of the satellite in inertial Cartesian coordinates and = the velocity. The equation of motion then takes on the form of an initial value problem: HybridSail : Hybrid Solar Sails for Active Debris Removal Page 13
18 = 16 "# (, ) =, = "#$ + "# + ""# + " Where is the mass of the satellite, and "#, = "#$% "#$"%$&' "#$$%"# + "#$ + "#$%&'()*+,# 17 is the sum of all the forces acting on the satellite. Their values are calculated from the models described in sections 3.1 to 3.3 The gravity acceleration due to the Earth gravity potential is found by calculating the gradient "#$ = "#$ " "#$ " "#$ " 18 And similarly for the third bodies (Sun and moon) = = " " " 19 The initial value problem in equation (16) can then be solved by any numeric ordinary difference equation method. For this implementation, the classical 4th order Runge Kutta (RK) method was used with a fixed time step, ℎ. ℎ , ℎ ℎ +, ℎ ℎ +, ℎ, + ℎ = + = = = = 2 This choice was made because of the ease of implementation. An adaptive step size RK method will yield a more optimal selection of the integration time step because then each time step is chosen to satisfy error bounds Semi analytic Satellite Theory (SST) Semi analytic Satellite Theory, as used in this study, is described completely by reference [8]. This section is a short summary of the theory. Semi analytic Satellite Theory represents the state of the satellite as an equinoctial element set. The set consists of 6 elements,,..,. It relies on the equinoctial reference frame, defined by the 3 basis vectors, and. These vectors have the following properties: HybridSail : Hybrid Solar Sails for Active Debris Removal Page 14
19 1. Vectors and lie in the orbital plane 2. is parallel to the angular momentum vector of the satellite 3. The angle between and the ascending node is equal to the longitude of the ascending node There are 2 choices for and and this choice will determine whether the direct or retrograde equinoctial element set is being used. The 6 elements of the equinoctial element set are =  the Keplerian semi major axis, = h,  the components of the eccentricity vector along and respectively, =,  the components of the ascending node vector along and respectively =  the mean longitude The retrograde factor,, determines whether the direct or retrograde equinoctial element set is being used. For the direct set, = 1, and for the retrograde set = 1. This choice is necessary to avoid singularities that exist at equatorial and polar orbits (inclination of and 9 ). The first 5 elements are slowly varying with time, while the mean longitude is fast varying. SST separates the osculating equinoctial elements into mean values and a remainder which varies periodically with the fast variable,. = +,..,, 21 Where represents the osculating element and the mean element. The effect of orbital perturbations on the mean elements can be described by ordinary differential equations " = +,..,, 22 Where " = 1 if = and otherwise (Kronecker delta). is the mean motion of the satellite. The short period variations are described by Fourier series expressions =,..,, cos " +,..,, sin " 23 Because the effect of perturbations on the mean elements are described by differential equations, the mean elements can be integrated numerically. Because the mean elements vary slowly, the integration time step can be quite long (typically a day or longer) much longer than Cowell integration. The same 4 th order Runge Kutta scheme as in equation (2) was used for the implemented SST propagator. Because we are only interested in satellite orbit lifetime for this study (and not in accurate position at every time step) the short period variations were omitted. HybridSail : Hybrid Solar Sails for Active Debris Removal Page 15
20 Reference [8] gives analytic expressions for all the contributions to the mean element differential equations for central body gravity potential, third body gravitational potential, aerodynamic drag and solar radiation pressure. The effect of the electrostatic force generated by the HybridSail tether is modelled in the same way as atmospheric drag in SST. That is, " = h "#$%&'$($)# " 24 Where the integration spans over the true longitude,. The perturbing acceleration due to electrostatic force is "#$%&'$($)# = "#$%&'$($)# 25 Where is the mass of the satellite, and "#$%&'$($)# is the electrostatic force calculated from (5) Validation of SST & Cowell methods The following test orbit was used to validate the developed orbit propagation methods. Satellite name Orbit Altitude Sumbandila LEO with drag, circular sun synchronous 5 km Eccentricity.15 Mass 8 kg Drag surface area.5 m 2 Two line elements U 949F Cowell orbit propagator The Cowell orbit propagator was initialized using position and velocity points from the SGP4 propagator. The reference SGP4 implementation was used to generate position and velocity at 1s intervals, for 1 orbit. The position and velocity at time t was used as a first guess for the Cowell propagator state, and the Cowell integration was performed for one orbit. The position error was used to refine the initial condition until a best fit was achieved. The graphs below shows the comparison between the SGP4 output and Cowell propagator output for 4 orbits. HybridSail : Hybrid Solar Sails for Active Debris Removal Page 16
21 SGP4 orbit posipon Cowell orbit posipon 1 X Y Z 1 X Y Z PosiPon (km) 55 PosiPon (km) Time (s) Time (s) Argument of perigee (deg) SGP4 orbit  argument of perigee Cowell orbit  argument of perigee Time (s) Time (s) Figure 9 SGP4 Cowell orbit comparison (short term) The position error for the first 4 orbits is plotted below. Argument of perigee (deg) 3 Cowell orbit propagator  posipon error 25 PosiPon error (m) Time (s) Figure 1 Cowell orbit propagator position error If the propagation is carried on past this point, the position error begins to increase as expected due to numeric errors and the simplified models of SGP4. The trend for slowly varying parameters should still remain the same for both methods. The plots below show the semi major axis and right ascension of the ascending node for 1 month. HybridSail : Hybrid Solar Sails for Active Debris Removal Page 17
22 689 SGP4 Semi major axis 689 Cowell Semi major axis Semi major axis (km) Semi major axis (km) Time (days) Time (days) 2 SGP4 Right ascension of the ascending node Cowell Right ascension of the ascending 2 node RAAN (deg) 1 RAAN (deg) Time (days) Figure 11 SGP4 Cowell orbit comparison (long term) Time (days) Semianalytic Satellite Theory propagator The same initial SGP4 orbit as before was used to initialize the SST propagator. 5 data points spread out over 3 orbits was used to perform a regression fit to the initial mean elements. Because the SST propagator only operates on the mean elements, a short term comparison with the SGP4 points is not possible. It is however possible to compare the slowly varying elements, taking into account that the SST orbital elements are mean elements with the osculating part removed. HybridSail : Hybrid Solar Sails for Active Debris Removal Page 18
23 689 Semi major axis Semi major axis (km) SGP4 SST Time (days),4,3,2,1 ,1 ,2 ,3 ,4 H component of eccentricity vector SGP4 SST 1 Time 2 (days) 3 4,4,3,2,1 ,1 ,2 ,3 ,4 K component of eccentricity vector 1 Time 2 (days) 3 4 SGP4 SST Figure 12 SGP4 SST orbit comparison 3.5 Attitude Coordinate system definitions The satellite coordinate system is defined as follows: the satellite body x axis is aligned with the normal vector of the sail and the y and z axes are aligned in the sail plane. The satellite attitude is modelled in the orbit coordinate system. The orbit coordinate system rotates with the orbit so that the z reference axis points towards nadir, the x reference axis points towards the velocity vector and the y axis points along the orbit anti normal. The satellite coordinate system will nominally be aligned with the orbit coordinate system when roll pitch and yaw equal. The Earth centered inertial coordinate frame (ECI) is used to propagate the position of the satellite and third bodies. A transformation is thus required to convert from ECI to orbit coordinate frames. This transformation can be obtained from the satellite position,, and velocity,, vector. / = 26 HybridSail : Hybrid Solar Sails for Active Debris Removal Page 19
24 3.5.2 Attitude representation The spacecraft attitude is the transformation between the orbit coordinate frame and the satellite body frame. This transformation can be expressed as a quaternion, direction cosine matrix (DCM) or Euler angles. When propagating the attitude we will make use of the quaternion representation to avoid singularities. / = 27 The quaternion attitude representation can be converted to a DCM using the following: / = We will also express the attitude as pitch,, roll,, and yaw,, Euler angles, defined as successive rotations. We will use a rotation, so the pitch rotation takes place first. The Euler angles can be computed from the DCM and vice versa. = atan2 ", = sin " = atan2 ", 29 / "# + "# = "#$ + "#" "#$ "#$ "# " "#$ + "#$% "# + "#$% "# 3 Where C is the cosine function and S the sine function Attitude propagation " " ", (the rotation rate of the The orbit referenced body rate vector, = spacecraft body frame relative to the orbit reference frame) is used to propagate the quaternion attitude. 1 " = 2 " " " " " " " " " " " 31 The relationship between the orbit referenced body rates and inertially referenced body rates,, is given by Where = = + / 32 is the orbit angular rate. The attitude dynamics equation gives the rate of change of the inertially referenced body rates. HybridSail : Hybrid Solar Sails for Active Debris Removal Page 2
25 = 33 Where is the sum of all the disturbance and control torques on the satellite and is the moment of inertia tensor. The disturbance torques that will be taken into consideration are the aerodynamic disturbance torque from equation (3), solar radiation pressure disturbance torque (equation 1) and gravity gradient torque. Rotations take place about the centre of mass of the satellite. The aerodynamic and solar radiation pressure torques are moments caused by the solar radiation force and aerodynamic force. These forces are applied at the centre of pressure of the sail, and the lever arm (connecting the centre of mass to the point where the force is applied) is the position vector from the centre of mass to the centre of pressure. The gravity gradient torqure,, is given by: The vector is the orbit nadir vector in body coordinates. = 3 34 = / 1 35 The attitude state can be propagated using a numerical integration method, using the differential equations of (31) and (33) and some initial state for the attitude and orbit referenced body rates. A modified Euler integration scheme (pyramid integration) was used to propagate the attitude for this study. 4 Satellite deorbit analysis This section will investigate the effectiveness of HybridSail as a de orbiting device, as well as the applications where it can be used and the parameters that influence the effectiveness. For this purpose the time to de orbit is evaluated and compared for the criteria under investigation. The de orbit time is the time it takes for a satellite to reach a post mission disposal orbit. For this study we consider only LEO satellites that are affected by drag and Geosynchronous Earth Orbit (GEO) satellites. The latter can be de orbited using the sail as a solar sail, and in the first case the satellite can be de orbited using a combination of aerodynamic drag and electrostatic braking. For a de orbit strategy to be considered effective, the de orbit time has to be below 25 years from the spacecraft end of life, to comply with [6]. The first sub section (4.1) will evaluate only aerodynamic drag augmentation. The next (4.2) will investigate the capabilities of electrostatic tethers as a de orbiting means. Section 4.3 considers the combined effect of drag de orbiting augmented with tethers. Section 4.4 looks at the collision probability of a sail in its de orbiting phase and section 4.5 looks at the option of using the sail as a solar sail to propel the host satellite to a higher graveyard orbit. HybridSail : Hybrid Solar Sails for Active Debris Removal Page 21
26 Only in section 4.5 is the solar force used constructively to de orbit. In all other cases the sail is used purely as a drag augmentation device although the disturbing effect of solar radiation pressure is still considered throughout. 4.1 Drag sail deorbiting in LEO Aerodynamic drag augmentation is achieved by maximizing the surface area perpendicular to the satellite velocity. This is depicted in the image below. Figure 13 Aerodynamic drag sail de orbiting In this case the attitude of the sail remains constant. Attitude is the orientation of the satellite body relative to the orbit coordinate system. For the attitude to remain constant the satellite has to rotate once per obit relative to an inertial reference frame Solar activity Solar activity plays a significant role in atmospheric models. It has been observed that there is a day night variation in the density of the upper atmosphere, sometimes referred to as the diurnal bulge. The bulge has a peak around 2p.m. local solar time and a minimum around 4a.m. [2]. But is has also been found that the 11 year cyclic variation in solar activity influences the density of the atmosphere. The daily variations in atmospheric density is such that the effect on aerodynamic drag over a long time period can be approximated by an average density, but the 11 year solar cycle can significantly influence de orbit times. HybridSail : Hybrid Solar Sails for Active Debris Removal Page 22
27 Atmospheric density at 5km alptude 4,5E 13 4,E 13 Atmosphere desnsity (kg/m3) 3,5E 13 3,E 13 2,5E 13 2,E 13 1,5E 13 1,E 13 5,E 14,E LaPtude (deg) Local solar Pme (hours) 22 Figure 14 Atmospheric density at 5km altitude showing the diurnal bulge Most atmospheric models take as input indices that reflect the current solar activity. The most widely used indices are the F 1.7 solar flux and a p or K p geomagnetic indices. These values are recorded by observatories over the world and historic data is readily available, but for future predicted solar activity, an estimation algorithm is required, such as the one described in section 3.1. The difference in density as a result of solar activity is illustrated by plotting the atmospheric density per altitude at the maximum and minimum epoch. The last minimum and maximum epochs were used. Table 41 Previous solar minimum and maximum Date a p index F 1.7 index Solar min December Solar max February HybridSail : Hybrid Solar Sails for Active Debris Removal Page 23
28 1,1,1 Atmospheric density comparison Solar Min Solar Max Density (kg/m3),1,1 1E 9 1E 11 1E 13 1E AlPtude (km) Figure 15 Atmospheric density as a function of altitude at solar max and min ( deg longitude, deg latitude, 12p.m. local solar time) The density at solar maximum can sometimes be more than 2 times that at corresponding conditions at solar minimum. The effect can further be demonstrated by comparing de orbit times under differing solar conditions. The following table shows de orbit times for the same satellite, starting to de orbit at solar maximum and solar minimum. It also shows the ratio of the two de orbiting times. Table 42 Comparison of de orbit times starting at solar maximum and minimum epochs Sail size (m 2 ) De orbit time starting at solar max (days) De orbit time starting at solar min (days) Ratio of solar min de orbit time vs. solar max The parameters for the simulation in the above table are: Physical parameters Satellite mass 1 kg Initial orbit Altitude 6 km Eccentricity.1 Inclination deg The plot below shows how the ratio of the de orbiting times change as a function of the solar max de orbit time. It can be seen that only when the time to de orbit is in the order of two times the solar cycle duration (>2 years) will the effect of solar activity even out and the starting point of de orbitation not play a role anymore. HybridSail : Hybrid Solar Sails for Active Debris Removal Page 24
29 RaPo: de orbit Pme starpng at solar min / de orbit Pme starpng at solar max Figure 16 Ratio of solar maximum and minimum de orbit times Sail sizing The parameters that influences the de orbit time of a drag sail is the physical satellite parameters (mass and sail area) and the initial orbit (altitude, eccentricity, inclination). To evaluate the effect of the physical satellite parameters, we note that the acceleration due to the drag force is a function of sail area over satellite mass. "#$ = "#$ = " 2 "#. "# "# "# 36 The ballistics coefficient, " =, is often used to describe this relation. Because the sail area scales together with the mass, it is only necessary to find the desired ballistics coefficient for a certain scenario and the sail area can be selected given the mass of the satellite. The table below shows the ballistics coefficient for a few sail sizes and satellite masses. A drag coefficient,, of 2.2 is assumed throughout. Table 43 Ballistics coefficient as a function of mass and sail area Ballistics Coefficient (kg/m 2 ) 5 1m 2 25m 2 1m 2 5m 2 RaPo of de orbit Pmes 2m 2 1m Time to de orbit starpng at solar max (years) Satellite mass (kg) Sail area (m 2 ) To simplify the analysis, we consider circular orbits at various altitudes. Also, the initial inclination for the orbit is chosen so that the orbit is sun synchronous. That is, it will have a nodal precession rate equal to the rate at which the earth rotates around the sun. Consequently the satellite will have the same lighting conditions every orbit, and thus the same atmospheric density profile. The graphs below show the de orbit times as a function of initial altitude and ballistics coefficient. The plotted de orbit times are the averages taken from a solar maximum and solar minimum starting point. HybridSail : Hybrid Solar Sails for Active Debris Removal Page 25
30 De orbit Pmes using aerodynamic drag De orbit Pme (years) BallisPc coefficient.6395 (kg/m2) IniPal alptude (km) De orbit Pmes (years) using aerodynamic drag BallisPc coefficient (kg/m2) IniPal alptude (km) Figure 17 De orbit times using aerodynamic drag From the above plot, it is possible to select the required sail size for a given satellite mass and initial altitude. For instance, to de orbit within 25 years from an initial 9km circular orbit, a ballistics coefficient of 1.6 is needed. If the mass of the satellite is know, the sail area can be calculated from HybridSail : Hybrid Solar Sails for Active Debris Removal Page 26
Orbital Mechanics and Space Geometry
Orbital Mechanics and Space Geometry AERO4701 Space Engineering 3 Week 2 Overview First Hour Coordinate Systems and Frames of Reference (Review) Kepler s equations, Orbital Elements Second Hour Orbit
More informationPenn State University Physics 211 ORBITAL MECHANICS 1
ORBITAL MECHANICS 1 PURPOSE The purpose of this laboratory project is to calculate, verify and then simulate various satellite orbit scenarios for an artificial satellite orbiting the earth. First, there
More informationSection 4: The Basics of Satellite Orbits
Section 4: The Basics of Satellite Orbits MOTION IN SPACE VS. MOTION IN THE ATMOSPHERE The motion of objects in the atmosphere differs in three important ways from the motion of objects in space. First,
More informationCHAPTER 2 ORBITAL DYNAMICS
14 CHAPTER 2 ORBITAL DYNAMICS 2.1 INTRODUCTION This chapter presents definitions of coordinate systems that are used in the satellite, brief description about satellite equations of motion and relative
More informationPhysics 9e/Cutnell. correlated to the. College Board AP Physics 1 Course Objectives
Physics 9e/Cutnell correlated to the College Board AP Physics 1 Course Objectives Big Idea 1: Objects and systems have properties such as mass and charge. Systems may have internal structure. Enduring
More informationAttitude and Orbit Dynamics of High AreatoMass Ratio (HAMR) Objects and
Attitude and Orbit Dynamics of High AreatoMass Ratio (HAMR) Objects and Carolin Früh National Research Council Postdoctoral Fellow, AFRL, cfrueh@unm.edu Orbital Evolution of Space Debris Objects Main
More informationOrbital Mechanics. Angular Momentum
Orbital Mechanics The objects that orbit earth have only a few forces acting on them, the largest being the gravitational pull from the earth. The trajectories that satellites or rockets follow are largely
More informationSittiporn Channumsin Coauthors
28 Oct 2014 Space Glasgow Research Conference Sittiporn Channumsin Sittiporn Channumsin Coauthors S. Channumsin Outline Background Objective The model Simulation Results Conclusion and Future work 2 Space
More information2. Orbits. FERZagreb, Satellite communication systems 2011/12
2. Orbits Topics Orbit types Kepler and Newton laws Coverage area Influence of Earth 1 Orbit types According to inclination angle Equatorial Polar Inclinational orbit According to shape Circular orbit
More informationFlight and Orbital Mechanics
Flight and Orbital Mechanics Lecture slides Challenge the future 1 Material for exam: this presentation (i.e., no material from text book). Sunsynchronous orbit: used for a variety of earthobserving
More informationChapter 2. Mission Analysis. 2.1 Mission Geometry
Chapter 2 Mission Analysis As noted in Chapter 1, orbital and attitude dynamics must be considered as coupled. That is to say, the orbital motion of a spacecraft affects the attitude motion, and the attitude
More informationSUNSYNCHRONOUS ORBIT SLOT ARCHITECTURE ANALYSIS AND DEVELOPMENT. A Thesis. Presented to. the Faculty of California Polytechnic State University
SUNSYNCHRONOUS ORBIT SLOT ARCHITECTURE ANALYSIS AND DEVELOPMENT A Thesis Presented to the Faculty of California Polytechnic State University San Luis Obispo In Partial Fulfillment of the Requirements
More informationCan Hubble be Moved to the International Space Station? 1
Can Hubble be Moved to the International Space Station? 1 On January 16, NASA Administrator Sean O Keefe informed scientists and engineers at the Goddard Space Flight Center (GSFC) that plans to service
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 informationEVOLUTION OF THE DEBRIS CLOUD GENERATED BY THE FENGYUN1C FRAGMENTATION EVENT
EVOLUTION OF THE DEBRIS CLOUD GENERATED BY THE FENGYUN1C FRAGMENTATION EVENT Carmen Pardini and Luciano Anselmo Space Flight Dynamics Laboratory Istituto di Scienza e Tecnologie dell Informazione Alessandro
More informationDynamics of Iain M. Banks Orbitals. Richard Kennaway. 12 October 2005
Dynamics of Iain M. Banks Orbitals Richard Kennaway 12 October 2005 Note This is a draft in progress, and as such may contain errors. Please do not cite this without permission. 1 The problem An Orbital
More informationUS ACTIVE DEBRIS REMOVAL (ADR) EFFORTS
US ACTIVE DEBRIS REMOVAL (ADR) EFFORTS DR. DARREN MCKNIGHT VIENNA UN COPUOS FEBRUARY 2013 1 DISCLAIMER This presentation solely represents the opinions of the author and should not be construed as being
More informationExamination Space Missions and Applications I (AE2103) Faculty of Aerospace Engineering Delft University of Technology SAMPLE EXAM
Examination Space Missions and Applications I AE2103 Faculty of Aerospace Engineering Delft University of Technology SAMPLE EXAM Please read these instructions first: This are a series of multiplechoice
More informationSatellite Mission Analysis
CARLETON UNIVERSITY SPACECRAFT DESIGN PROJECT 2004 FINAL DESIGN REPORT Satellite Mission Analysis FDR Reference Code: FDRSAT20043.2.A Team/Group: Satellite Mission Analysis Date of Submission: April
More informationOrbital Dynamics: Formulary
Orbital Dynamics: Formulary 1 Introduction Prof. Dr. D. Stoffer Department of Mathematics, ETH Zurich Newton s law of motion: The net force on an object is equal to the mass of the object multiplied by
More informationRS platforms. Fabio Dell Acqua  Gruppo di Telerilevamento
RS platforms Platform vs. instrument Sensor Platform Instrument The remote sensor can be ideally represented as an instrument carried by a platform Platforms Remote Sensing: Groundbased airborne spaceborne
More informationSatellite Posi+oning. Lecture 5: Satellite Orbits. Jan Johansson jan.johansson@chalmers.se Chalmers University of Technology, 2013
Lecture 5: Satellite Orbits Jan Johansson jan.johansson@chalmers.se Chalmers University of Technology, 2013 Geometry Satellite Plasma Posi+oning physics Antenna theory Geophysics Time and Frequency GNSS
More informationCoverage Characteristics of Earth Satellites
Coverage Characteristics of Earth Satellites This document describes two MATLAB scripts that can be used to determine coverage characteristics of single satellites, and Walker and userdefined satellite
More informationDevelopment of a Sun Synchronous. Conjunctions
Development of a Sun Synchronous Zoning Architecture to Minimize Conjunctions Kevin Shortt Brian Weeden Secure World Foundation www.secureworldfoundation.org Overview Current Situation in Sun synchronous
More informationSpacecraft Dynamics and Control. An Introduction
Brochure More information from http://www.researchandmarkets.com/reports/2328050/ Spacecraft Dynamics and Control. An Introduction Description: Provides the basics of spacecraft orbital dynamics plus attitude
More informationAnalysis on the Longterm Orbital Evolution and Maintenance of KOMPSAT2
Analysis on the Longterm Orbital Evolution and Maintenance of KOMPSAT2 OkChul Jung 1 Korea Aerospace Research Institute (KARI), 45 Eoeundong, Daejeon, South Korea, 305333 JungHoon Shin 2 Korea Advanced
More informationGravity Field and Dynamics of the Earth
Milan Bursa Karel Pec Gravity Field and Dynamics of the Earth With 89 Figures SpringerVerlag Berlin Heidelberg New York London Paris Tokyo HongKong Barcelona Budapest Preface v Introduction 1 1 Fundamentals
More informationDevelopment of an automated satellite network management system
Development of an automated satellite network management system Iasonas Kytros Christos Porios Nikitas Terzoudis Varvara Chatzipavlou Coach: Sitsanlis Ilias February 2013 Abstract In this paper we present
More informationDEOS. Deutsche Orbitale Servicing Mission. The Inflight Technology Demonstration of Germany s Robotics Approach to Service Satellites
DEOS Deutsche Orbitale Servicing Mission The Inflight Technology Demonstration of Germany s Robotics Approach to Service Satellites B. Sommer, K. Landzettel, T. Wolf, D. Reintsema, German Aerospace Center
More informationHalliday, Resnick & Walker Chapter 13. Gravitation. Physics 1A PHYS1121 Professor Michael Burton
Halliday, Resnick & Walker Chapter 13 Gravitation Physics 1A PHYS1121 Professor Michael Burton II_A2: Planetary Orbits in the Solar System + Galaxy Interactions (You Tube) 21 seconds 131 Newton's Law
More informationInteraction of Energy and Matter Gravity Measurement: Using Doppler Shifts to Measure Mass Concentration TEACHER GUIDE
Interaction of Energy and Matter Gravity Measurement: Using Doppler Shifts to Measure Mass Concentration TEACHER GUIDE EMR and the Dawn Mission Electromagnetic radiation (EMR) will play a major role in
More informationArtificial Satellites Earth & Sky
Artificial Satellites Earth & Sky Name: Introduction In this lab, you will have the opportunity to find out when satellites may be visible from the RPI campus, and if any are visible during the activity,
More informationLecture L17  Orbit Transfers and Interplanetary Trajectories
S. Widnall, J. Peraire 16.07 Dynamics Fall 008 Version.0 Lecture L17  Orbit Transfers and Interplanetary Trajectories In this lecture, we will consider how to transfer from one orbit, to another or to
More informationG U I D E T O A P P L I E D O R B I T A L M E C H A N I C S F O R K E R B A L S P A C E P R O G R A M
G U I D E T O A P P L I E D O R B I T A L M E C H A N I C S F O R K E R B A L S P A C E P R O G R A M CONTENTS Foreword... 2 Forces... 3 Circular Orbits... 8 Energy... 10 Angular Momentum... 13 FOREWORD
More informationPresentation of problem T1 (9 points): The Maribo Meteorite
Presentation of problem T1 (9 points): The Maribo Meteorite Definitions Meteoroid. A small particle (typically smaller than 1 m) from a comet or an asteroid. Meteorite: A meteoroid that impacts the ground
More informationDEVELOPMENT OF AN ARCHITECTURE OF SUNSYNCHRONOUS ORBITAL SLOTS TO MINIMIZE CONJUNCTIONS. Brian Weeden Secure World Foundation
DEVELOPMENT OF AN ARCHITECTURE OF SUNSYNCHRONOUS ORBITAL SLOTS TO MINIMIZE CONJUNCTIONS Brian Weeden Secure World Foundation Sunsynchronous orbit (SSO) satellites serve many important functions, primarily
More informationCalculation of Azimuth, Elevation and Polarization for nonhorizontal aligned Antennas
Calculation of Azimuth, Elevation and Polarization for nonhorizontal aligned Antennas Algorithm Description Technical Document TD1205a Version 1.1 23.10.2012 In Cooperation with 1 Objective Many SatCom
More informationThe µtorque MomentumExchange Tether Experiment
The µtorque MomentumExchange Tether Experiment Robert P. Hoyt Tethers Unlimited, Inc., 19011 36 h Ave. W., Suite F, Lynnwood, WA 980365752 (425) 7440400 TU@tethers.com Abstract. Long, highstrength
More informationSection 2. Satellite Orbits
Section 2. Satellite Orbits References Kidder and Vonder Haar: chapter 2 Stephens: chapter 1, pp. 2530 Rees: chapter 9, pp. 174192 In order to understand satellites and the remote sounding data obtained
More informationPhysics 2A, Sec B00: Mechanics  Winter 2011 Instructor: B. Grinstein Final Exam
Physics 2A, Sec B00: Mechanics  Winter 2011 Instructor: B. Grinstein Final Exam INSTRUCTIONS: Use a pencil #2 to fill your scantron. Write your code number and bubble it in under "EXAM NUMBER;" an entry
More informationState of Stress at Point
State of Stress at Point Einstein Notation The basic idea of Einstein notation is that a covector and a vector can form a scalar: This is typically written as an explicit sum: According to this convention,
More informationUnderstanding Orbital Mechanics Through a StepbyStep Examination of the SpaceBased Infrared System (SBIRS)
Understanding Orbital Mechanics Through a StepbyStep Examination of the SpaceBased Infrared System (SBIRS) Denny Sissom Elmco, Inc. May 2003 Pg 1 of 27 SSMD1102366 [1] The GroundBased Midcourse Defense
More informationAttitude Control and Dynamics of Solar Sails
Attitude Control and Dynamics of Solar Sails Benjamin L. Diedrich A thesis submitted in partial fulfillment of the requirements for the degree of Master of Science in Aeronautics & Astronautics University
More informationHalliday, Resnick & Walker Chapter 13. Gravitation. Physics 1A PHYS1121 Professor Michael Burton
Halliday, Resnick & Walker Chapter 13 Gravitation Physics 1A PHYS1121 Professor Michael Burton II_A2: Planetary Orbits in the Solar System + Galaxy Interactions (You Tube) 21 seconds 131 Newton's Law
More informationMechanics lecture 7 Moment of a force, torque, equilibrium of a body
G.1 EE1.el3 (EEE1023): Electronics III Mechanics lecture 7 Moment of a force, torque, equilibrium of a body Dr Philip Jackson http://www.ee.surrey.ac.uk/teaching/courses/ee1.el3/ G.2 Moments, torque and
More informationUSING MS EXCEL FOR DATA ANALYSIS AND SIMULATION
USING MS EXCEL FOR DATA ANALYSIS AND SIMULATION Ian Cooper School of Physics The University of Sydney i.cooper@physics.usyd.edu.au Introduction The numerical calculations performed by scientists and engineers
More informationState Newton's second law of motion for a particle, defining carefully each term used.
5 Question 1. [Marks 28] An unmarked police car P is, travelling at the legal speed limit, v P, on a straight section of highway. At time t = 0, the police car is overtaken by a car C, which is speeding
More informationExperiment #1, Analyze Data using Excel, Calculator and Graphs.
Physics 182  Fall 2014  Experiment #1 1 Experiment #1, Analyze Data using Excel, Calculator and Graphs. 1 Purpose (5 Points, Including Title. Points apply to your lab report.) Before we start measuring
More informationSample Questions for the AP Physics 1 Exam
Sample Questions for the AP Physics 1 Exam Sample Questions for the AP Physics 1 Exam Multiplechoice Questions Note: To simplify calculations, you may use g 5 10 m/s 2 in all problems. Directions: Each
More informationNewton s Law of Gravity
Gravitational Potential Energy On Earth, depends on: object s mass (m) strength of gravity (g) distance object could potentially fall Gravitational Potential Energy In space, an object or gas cloud has
More informationNewton s Law of Universal Gravitation
Newton s Law of Universal Gravitation The greatest moments in science are when two phenomena that were considered completely separate suddenly are seen as just two different versions of the same thing.
More informationRotation: Moment of Inertia and Torque
Rotation: Moment of Inertia and Torque Every time we push a door open or tighten a bolt using a wrench, we apply a force that results in a rotational motion about a fixed axis. Through experience we learn
More informationThe Sun. Solar radiation (Sun EarthRelationships) The Sun. The Sun. Our Sun
The Sun Solar Factoids (I) The sun, a mediumsize star in the milky way galaxy, consisting of about 300 billion stars. (Sun EarthRelationships) A gaseous sphere of radius about 695 500 km (about 109 times
More informationAstrodynamics (AERO0024)
Astrodynamics (AERO0024) 6. Interplanetary Trajectories Gaëtan Kerschen Space Structures & Systems Lab (S3L) Course Outline THEMATIC UNIT 1: ORBITAL DYNAMICS Lecture 02: The TwoBody Problem Lecture 03:
More informationThis paper is also taken for the relevant Examination for the Associateship. For Second Year Physics Students Wednesday, 4th June 2008: 14:00 to 16:00
Imperial College London BSc/MSci EXAMINATION June 2008 This paper is also taken for the relevant Examination for the Associateship SUN, STARS, PLANETS For Second Year Physics Students Wednesday, 4th June
More informationSatellite Orbital Decay Calculations
IPS RADIO AND SPACE SERVICES Level 15, Tower C 300 Elizabeth Street Sydney, NSW 2000 PO Box 1386 Haymarket, NSW 1240 Australia Tel: +61 2 92138000 Faz: +61 2 92138060 office@ips.gov.au www.ips.gov.au Satellite
More information11.1. Objectives. Component Form of a Vector. Component Form of a Vector. Component Form of a Vector. Vectors and the Geometry of Space
11 Vectors and the Geometry of Space 11.1 Vectors in the Plane Copyright Cengage Learning. All rights reserved. Copyright Cengage Learning. All rights reserved. 2 Objectives! Write the component form of
More informationPhysics Notes Class 11 CHAPTER 6 WORK, ENERGY AND POWER
1 P a g e Work Physics Notes Class 11 CHAPTER 6 WORK, ENERGY AND POWER When a force acts on an object and the object actually moves in the direction of force, then the work is said to be done by the force.
More informationAnalysis of the 2007 Chinese ASAT Test and the Impact of its Debris on the Space Environment. T.S. Kelso Center for Space Standards & Innovation
Analysis of the 2007 Chinese ASAT Test and the Impact of its Debris on the Space Environment T.S. Kelso Center for Space Standards & Innovation ABSTRACT On 2007 January 11, the People s Republic of China
More informationSolar System. 1. The diagram below represents a simple geocentric model. Which object is represented by the letter X?
Solar System 1. The diagram below represents a simple geocentric model. Which object is represented by the letter X? A) Earth B) Sun C) Moon D) Polaris 2. Which object orbits Earth in both the Earthcentered
More informationAbstract: Formation flying is described as harnessing and exploring the
FORMATION FLYING MISSION WITH SIMPLIFIED GUIDANCE M.P.RAMACHANDRAN Scientist/Engineer, Flight Dynamics Group, ISRO Satellite Centre, Airport Road, Bangalore 560 017, India, 91802508 4409, mprama@isac.gov.in
More informationEN4 Dynamics and Vibrations. Design Project. Orbital Design for a Lunar Impact Mission. Synopsis
EN4 Dynamics and Vibrations Design Project Orbital Design for a Lunar Impact Mission Synopsis NASA has identified a need for a lowcost mission to launch a satellite that will impact the moon. You will
More informationBasic Principles of Inertial Navigation. Seminar on inertial navigation systems Tampere University of Technology
Basic Principles of Inertial Navigation Seminar on inertial navigation systems Tampere University of Technology 1 The five basic forms of navigation Pilotage, which essentially relies on recognizing landmarks
More information1. Fluids Mechanics and Fluid Properties. 1.1 Objectives of this section. 1.2 Fluids
1. Fluids Mechanics and Fluid Properties What is fluid mechanics? As its name suggests it is the branch of applied mechanics concerned with the statics and dynamics of fluids  both liquids and gases.
More informationPHY121 #8 Midterm I 3.06.2013
PHY11 #8 Midterm I 3.06.013 AP Physics Newton s Laws AP Exam Multiple Choice Questions #1 #4 1. When the frictionless system shown above is accelerated by an applied force of magnitude F, the tension
More informationACTIVE DEBRIS REMOVAL MISSION WITH SMALL SATELLITE
ACTIVE DEBRIS REMOVAL MISSION WITH SMALL SATELLITE Marcello Valdatta Niccolò Bellini Davide Rastelli University of Bologna, ITALY 2nd Mission Idea Contest 4th Nanosatellite Symposium on October 10th 
More informationName: Earth 110 Exploration of the Solar System Assignment 1: Celestial Motions and Forces Due in class Tuesday, Jan. 20, 2015
Name: Earth 110 Exploration of the Solar System Assignment 1: Celestial Motions and Forces Due in class Tuesday, Jan. 20, 2015 Why are celestial motions and forces important? They explain the world around
More informationDoes currently available technology have the capacity to facilitate a manned mission to Mars?
Furze Platt Senior School Does currently available technology have the capacity to facilitate a manned mission to Mars? Daniel Messias Date: 8/03/2015 Candidate Number: 7158 Centre Number: 51519 Contents
More informationSearching for space debris elements with the Pi of the Sky system
Searching for space debris elements with the Pi of the Sky system Marcin Sokołowski msok@fuw.edu.pl Soltan Institute for Nuclear Studies ( IPJ ) Warsaw, Poland 7th Integral / BART Workshop ( IBWS), 1418
More informationSpacecraft Power Systems
Spacecraft Power Systems AOE 4065 Space Design Refs: SMAD Chap 11.4, G&F Chap 10, F&S Chap 11, P&M Chap 6 Electrical Power Subsystem (EPS) Functions Supply electrical power to spacecraft loads Control
More informationMobile Communications Chapter 5: Satellite Systems
Mobile Communications Chapter 5: Satellite Systems History Basics Localization Handover Routing Systems History of satellite communication 1945 Arthur C. Clarke publishes an essay about Extra Terrestrial
More informationEDMONDS COMMUNITY COLLEGE ASTRONOMY 100 Winter Quarter 2007 Sample Test # 1
Instructor: L. M. Khandro EDMONDS COMMUNITY COLLEGE ASTRONOMY 100 Winter Quarter 2007 Sample Test # 1 1. An arc second is a measure of a. time interval between oscillations of a standard clock b. time
More informationSatellite Communications
Satellite Communications Department of Electrical Engineering Faculty of Engineering Chiangmai University Origin of Satellite Communications Arthur C. Clark (1945) British Science fiction writer propose
More informationIndiana's Academic Standards 2010 ICP Indiana's Academic Standards 2016 ICP. map) that describe the relationship acceleration, velocity and distance.
.1.1 Measure the motion of objects to understand.1.1 Develop graphical, the relationships among distance, velocity and mathematical, and pictorial acceleration. Develop deeper understanding through representations
More informationINFLIGHT CALIBRATION OF THE MICROSCOPE SPACE MISSION INSTRUMENT: DEVELOPMENT OF THE SIMULATOR
SF2A 2011 G. Alecian, K. Belkacem, R. Samadi and D. VallsGabaud (eds) INFLIGHT CALIBRATION OF THE MICROSCOPE SPACE MISSION INSTRUMENT: DEVELOPMENT OF THE SIMULATOR E. Hardy 1, A. Levy 1, G. Métris 2,
More informationDIRECT ORBITAL DYNAMICS: USING INDEPENDENT ORBITAL TERMS TO TREAT BODIES AS ORBITING EACH OTHER DIRECTLY WHILE IN MOTION
1 DIRECT ORBITAL DYNAMICS: USING INDEPENDENT ORBITAL TERMS TO TREAT BODIES AS ORBITING EACH OTHER DIRECTLY WHILE IN MOTION Daniel S. Orton email: dsorton1@gmail.com Abstract: There are many longstanding
More informationMobile Communications: Satellite Systems
Mobile Communications: Satellite Systems Mobile Communication: Satellite Systems  Jochen Schiller http://www.jochenschiller.de 1 History of satellite communication 1945 Arthur C. Clarke publishes an essay
More informationQuasiSynchronous Orbits
QuasiSynchronous Orbits and Preliminary Mission Analysis for Phobos Observation and Access Orbits Paulo J. S. Gil Instituto Superior Técnico Simpósio Espaço 50 anos do 1º Voo Espacial Tripulado 12 de
More informationThe Threat of Orbital Debris and Protecting NASA Space Assets from Satellite Collisions 28 April 2009
The Threat of Orbital Debris and Protecting NASA Space Assets from Satellite Collisions 28 April 2009 Executive Summary Collision avoidance processes are in place for all NASA human space flight missions
More informationSTEREO Guidance & Control
STEREO Guidance & Control J. Courtney Ray J.C.Ray@jhuapl.edu J. C. Ray 98/11/19 1 STEREO G&C Requirements Baseline System Software Some Analysis J. C. Ray 98/11/19 2 G&C Requirements  Drivers Spacecraft
More informationIonosphere Properties and Behaviors  Part 2 By Marcel H. De Canck, ON5AU
Ionosphere Properties and Behaviors  Part 2 By Marcel H. De Canck, ON5AU I n the previous issue I explained that gyrofrequency depends on the earth s magnetic field and mentioned that this magnetic field
More information5. Satellite Systems. History of Satellite Communications
5. Satellite Systems History and Orbits Routing, Localization, and Handover Systems 2005 Burkhard Stiller and Jochen Schiller FU Berlin M5 1 History of Satellite Communications 1945 Arthur C. Clarke about
More informationSIX DEGREEOFFREEDOM MODELING OF AN UNINHABITED AERIAL VEHICLE. A thesis presented to. the faculty of
SIX DEGREEOFFREEDOM MODELING OF AN UNINHABITED AERIAL VEHICLE A thesis presented to the faculty of the Russ College of Engineering and Technology of Ohio University In partial fulfillment of the requirement
More informationName Class Date. true
Exercises 131 The Falling Apple (page 233) 1 Describe the legend of Newton s discovery that gravity extends throughout the universe According to legend, Newton saw an apple fall from a tree and realized
More informationGünter Seeber. Satellite Geodesy 2nd completely revised and extended edition
Günter Seeber Satellite Geodesy 2nd completely revised and extended edition Walter de Gruyter Berlin New York 2003 Contents Preface Abbreviations vii xvii 1 Introduction 1 1.1 Subject of Satellite Geodesy...
More informationThe TwoBody Problem
The TwoBody Problem Abstract In my short essay on Kepler s laws of planetary motion and Newton s law of universal gravitation, the trajectory of one massive object near another was shown to be a conic
More informationCRANFIELD UNIVERSITY CONTROL OF COLOCATED GEOSTATIONARY SATELLITES
CRANFIELD UNIVERSITY S. HARDACRE CONTROL OF COLOCATED GEOSTATIONARY SATELLITES COLLEGE OF AERONAUTICS PhD THESIS CRANFIELD UNIVERSITY COLLEGE OF AERONAUTICS PhD THESIS Academic Year 19956 S. HARDACRE
More informationAstronomy 110 Homework #04 Assigned: 02/06/2007 Due: 02/13/2007. Name:
Astronomy 110 Homework #04 Assigned: 02/06/2007 Due: 02/13/2007 Name: Directions: Listed below are twenty (20) multiplechoice questions based on the material covered by the lectures this past week. Choose
More informationORBITAL DISPERSION SIMULATION OF NEAREARTH OBJECT DEFLECTION/FRAGMENTATION BY NUCLEAR EXPLOSIONS
IAC09C1.10.2 ORBITAL DISPERSION SIMULATION OF NEAREARTH OBJECT DEFLECTION/FRAGMENTATION BY NUCLEAR EXPLOSIONS B. Kaplinger Iowa State University, United States bdkaplin@iastate.edu B. Wie Iowa State
More informationThe Four Seasons. A Warm Up Exercise. A Warm Up Exercise. A Warm Up Exercise. The Moon s Phases
The Four Seasons A Warm Up Exercise What fraction of the Moon s surface is illuminated by the Sun (except during a lunar eclipse)? a) Between zero and onehalf b) The whole surface c) Always half d) Depends
More informationScience Standard 4 Earth in Space Grade Level Expectations
Science Standard 4 Earth in Space Grade Level Expectations Science Standard 4 Earth in Space Our Solar System is a collection of gravitationally interacting bodies that include Earth and the Moon. Universal
More informationDepartment of Aeronautics and Astronautics School of Engineering Massachusetts Institute of Technology. Graduate Program (S.M., Ph.D., Sc.D.
Department of Aeronautics and Astronautics School of Engineering Massachusetts Institute of Technology Graduate Program (S.M., Ph.D., Sc.D.) Field: Space Propulsion Date: October 15, 2013 1. Introduction
More informationForces on the Rocket. Rocket Dynamics. Equation of Motion: F = Ma
Rocket Dynamics orces on the Rockets  Drag Rocket Stability Rocket Equation Specific Impulse Rocket otors Thrust orces on the Rocket Equation of otion: = a orces at through the Center of ass Center of
More informationAP1 Oscillations. 1. Which of the following statements about a springblock oscillator in simple harmonic motion about its equilibrium point is false?
1. Which of the following statements about a springblock oscillator in simple harmonic motion about its equilibrium point is false? (A) The displacement is directly related to the acceleration. (B) The
More informationSound. References: L.D. Landau & E.M. Lifshitz: Fluid Mechanics, Chapter VIII F. Shu: The Physics of Astrophysics, Vol. 2, Gas Dynamics, Chapter 8
References: Sound L.D. Landau & E.M. Lifshitz: Fluid Mechanics, Chapter VIII F. Shu: The Physics of Astrophysics, Vol., Gas Dynamics, Chapter 8 1 Speed of sound The phenomenon of sound waves is one that
More informationNotes: Most of the material in this chapter is taken from Young and Freedman, Chap. 13.
Chapter 5. Gravitation Notes: Most of the material in this chapter is taken from Young and Freedman, Chap. 13. 5.1 Newton s Law of Gravitation We have already studied the effects of gravity through the
More informationRenewable Energy. Solar Power. Courseware Sample 86352F0
Renewable Energy Solar Power Courseware Sample 86352F0 A RENEWABLE ENERGY SOLAR POWER Courseware Sample by the staff of LabVolt Ltd. Copyright 2009 LabVolt Ltd. All rights reserved. No part of this
More informationMobile Computing. Chapter 5: Satellite Systems
Mobile Computing Chapter 5: Satellite Systems Prof. SangJo Yoo History of satellite communication 1945 Arthur C. Clarke publishes an essay about Extra Terrestrial Relays 1957 First satellite SPUTNIK by
More informationCBERS Program Update Jacie 2011. Frederico dos Santos Liporace AMS Kepler liporace@amskepler.com
CBERS Program Update Jacie 2011 Frederico dos Santos Liporace AMS Kepler liporace@amskepler.com Overview CBERS 3 and 4 characteristics Differences from previous CBERS satellites (CBERS 1/2/2B) Geometric
More informationLearn From The Proven Best!
Applied Technology Institute (ATIcourses.com) Stay Current In Your Field Broaden Your Knowledge Increase Productivity 349 Berkshire Drive Riva, Maryland 1140 888501100 4109568805 Website: www.aticourses.com
More information