16.522, Space Propulsion Prof. Manuel martinez-Sanchez Lecture 4: Re-positioning in Orbits Suppose we want to move a satellite in a circular orbit to a position A9 apart in the same orbit in a time At (assumed to be several orbital times at least). The general approach is to transfer to a lower(for positive A9)or higher(for 49<0)nearby orbit then drift in this faster (or slower) orbit for a certain time then return to the original orbit The analysis is similar for low and high thrust, because we have radius ratios very close to 1, so that in either case the satellite is nearly in orbit"even during thrusting periods, and Av's for orbit transfer amount (in magnitude) to the difference of the beginning and ending orbital speeds. In detail of course, if done at high thrust the maneuver involves a two-impulse Hohmann transfer to the drift orbit and one other two-impulse Hohmann transfer back to the original orbit For the low-thrust case, continuo us thrusting is used during both legs, with some guidance required to remove the very slight radial component of v picked up during spiral flight (although ignored here) We will do the analysis for the low-thrust case only then adapt the result for high- et 89 be the advance angle relative to a hypothetical satellite remaining in the original orbit and left undisturbed the general shape of the maneuver is sketched Forward △e Thrust △6-6 coasting se Δt-t1 Backwards Thrust Since the orbital angular velocity is its variation with orbit radius is 16.522, Space Propulsion Lecture 4 Prof. Manuel martinez-Sanchez Page 1 of 10
16.522, Space Propulsion Prof. Manuel Martinez-Sanchez Lecture 4: Re-positioning in Orbits Suppose we want to move a satellite in a circular orbit to a position ∆ϑ apart in the same orbit, in a time ∆t (assumed to be several orbital times at least). The general approach is to transfer to a lower (for positive ∆ϑ ) or higher (for ∆ϑ < 0 ) nearby orbit, then drift in this faster (or slower) orbit for a certain time, then return to the original orbit. The analysis is similar for low and high thrust, because we have radius ratios very close to 1, so that, in either case the satellite is nearly “in orbit” even during thrusting periods, and ∆V's for orbit transfer amount (in magnitude) to the difference of the beginning and ending orbital speeds. In detail, of course, if done at high thrust the maneuver involves a two-impulse Hohmann transfer to the drift orbit and one other two-impulse Hohmann transfer back to the original orbit. For the low-thrust case, continuous thrusting is used during both legs, with some guidance required to JG remove the very slight radial component of v picked up during spiral flight (although ignored here). We will do the analysis for the low-thrust case only, then adapt the result for highthrust. Let δϑ be the advance angle relative to a hypothetical satellite remaining in the original orbit and left undisturbed. The general shape of the maneuver is sketched below: Since the orbital angular velocity is Ω = µ 3 , r its variation with orbit radius is 16.522, Space Propulsion Lecture 4 Prof. Manuel Martinez-Sanchez Page 1 of 10
During thrusting, Sr is varying according to d v≡a F 日= d 1dr 2ar1/2 (3) and so d(a2) d(89)3 2ar1/2 3 (4) For integration, we will regard r=roas a constant(small variations) d(89) 3a =--t+ constant Starting from t=0, 69=00(89)20 e obtain t69=3at(t<t,) After(t=ti), we continue to drift at a constant rate dt and, since we start from (t,)2To the 59 during the coasting phase is 3 a t2-t1(t-t1)=- t At the end of coasting(t=At-t), we have then 16.522, Space P rtinez- sanchez Lecture 4 Prof. Manuel martinez
3 Ω δ r d (δϑ) δΩ = - = (1) 2 r dt During thrusting, δ r is varying according to d ⎛ µ ⎞ F r ⎝ a= M ⎠ ⎟ (2) dt ⎜ ⎝ - 2r ⎟ ⎠ ≅ M v ≅ a µ ⎛ F ⎞ ⎜ or µ dr ≅ a µ → 1dr = 2ar 12 (3) 2r2 dt r rdt µ and so d (δΩ) d2 (δϑ) 3 2ar12 3a = = - Ω = − (4) dt dt2 2 µ r For integration, we will regard r ≅ r0 as a constant (small variations): d (δϑ) 3a = - t constant + (5) dt r0 Starting from t= 0, = 0, δϑ d (δϑ) = 0, dt we obtain d (δϑ) = - 3a t ; δϑ = - 3 a t2 (t < t1 ) (6) dt r0 2 r0 After (t=t1), we continue to drift at a constant rate d (δϑ) 3a = - t1 dt r0 and, since we start from δϑ ( ) 3 a 2 t= - t1 , 1 2 r0 the δϑ during the coasting phase is t1 ⎞ δϑ = - 3 a t1 2 - 3a t1 (t - t1 ) = - 3at1 ⎜ ⎛ t - ⎟ (7) 2 r0 r r0 0 ⎝ 2 ⎠ At the end of coasting (t= ∆t t - ) , we have then 1 16.522, Space Propulsion Lecture 4 Prof. Manuel Martinez-Sanchez Page 2 of 10
Bat 6(△t+t)=-2 (8) and after the period t of reversed thrust we return to the initial orbit with O, and with 89 as in( 8), plus a further 89(t,) This is then the total 49 3at At -t (9) Clearly, the mission(given A9 and At )can be accomplished with different choices of thrusting time ti(but notice that t, At/ 2 in any case). The required thrust/mass ratio(a)and AV=2 at depend on this choice (10 3t1(△t-t1) 4v=2 3 4『o△6 2 We can see here that low thrust ends up as a penalty on Av, so that the thrusting time should be selected as short as possible within the available on-board power. 16.522, Space P pessan Lecture 4 Prof. Manuel martinez
3 ⎞ δϑ (∆tt ) = - 3at1 ⎜ ⎛ ∆t - t1 ⎟ - (8) 1 r0 ⎝ 2 ⎠ and, after the period t1 of reversed thrust, we return to the initial orbit with 0 d (δϑ ) =, and with δϑ as in (8), plus a further δϑ (t1 ) . dt This is then the total ∆ϑ : 3at1 ⎜ ⎛ ∆t - 3 ⎞ 3 a 2 ∆ϑ = - r0 ⎝ 2 t1 ⎠ ⎟ - 2 r0 t1 ∆ϑ = - 3at1 (∆t - t1 ) (9) r0 Clearly, the mission (given ∆ϑ and ∆t ) can be accomplished with different choices of thrusting time t1 (but notice that t1 < ∆t 2 in any case). The required thrust/mass ratio (a) and V=2 ∆ t a 1depend on this choice: 1 r0∆ϑ -a = 3 t1 (∆t - t1 ) (10) 2 r0∆ϑ ∆V= 3 (∆t t1 ) (11) - We can see here that low thrust ends up as a penalty on ∆V , so that the thrusting time should be selected as short as possible within the available on-board power. 16.522, Space Propulsion Lecture 4 Prof. Manuel Martinez-Sanchez Page 3 of 10
In the limit of impulsive thrust we realize that t, cannot really be any less than the Hohmann transfer time(1/2 orbit, or t,=-). A more detailed analysis of this case confirms that, for the high thrust case, Equations(10)and(11)are indeed valid witht,=r/Q2 The power per unit mass required is: P Fcac M 2n M P 1 040C M6nt1(△t-t1) It is useful to express results in terms of the coasting time At t=At-t 3 At +t t1(△t-t1) We then have △V 4。40 3△t+t P2c△6 M3n△t2-t coasting reduces Av, but increases P/M (not much if te /At is small) Sun-Moon Effects on a Geosynchronous Orbit- The Ns Drift Geosynchronous orbits are in the Equatorial plane, which is inclined 23.44 to the plane of the Ecliptic (path of Earth around the Sun). Because of this the tidal or gravity-gradient forces exerted by the Sun on a geosynchronous spacecraft will sometimes exert a" torque"on the orbit. This happens primarily at Solstice ( January and July) when the orbit dips the most out of the Ecliptic 16.522, Space P pessan Lecture 4 Prof. Manuel martinez
In the limit of impulsive thrust, we realize that t1 cannot really be any less than the π Hohmann transfer time (1 2 orbit, or t = ). A more detailed analysis of this case 1 Ω confirms that, for the high thrust case, Equations (10) and (11) are indeed valid with t1 = π Ω . The power per unit mass required is: P 1 F c a c = = M 2η M 2η P = 1 r0∆θ c M 6η t - 1 (∆t t 1 ) It is useful to express results in terms of the coasting time, t = ∆t - 2t c 1 ∆t t - c ∆t t + c t - 1 = 2 ⇒ ∆t t 1 = 2 , 2 c t - 1 (∆t t ) = ∆t2 + t , 1 4 We then have 4 r0∆θ ∆V = 3 ∆t t + c P = 2c r0∆θ 2 - 2 M 3η ∆t t c coasting reduces ∆V , but increases P M (not much if tc ∆t is small). Sun-Moon Effects on a Geosynchronous Orbit – The N/S Drift Geosynchronous orbits are in the Equatorial plane, which is inclined 23.44D to the plane of the Ecliptic (path of Earth around the Sun). Because of this, the tidal or gravity-gradient forces exerted by the Sun on a geosynchronous spacecraft will sometimes exert a “torque” on the orbit. This happens primarily at Solstice (January and July) when the orbit dips the most out of the Ecliptic: 16.522, Space Propulsion Lecture 4 Prof. Manuel Martinez-Sanchez Page 4 of 10
SUMMER GEOGRAPHIC SIDEREA SOLSTICE 23° S/C IN GEO s TIDAL ong vernal Equinox FORCES Ecliptic, solsticial line z: North sidereal EQUIN If r is the relative position vector of the spacecraft and s is the unit vector from the Sun, the Tidal Force( Gravity gradient) per unit mass is R (1) (This the imbalance between solar attraction and centrifugal force, which cancels on the earth, r=0). Of this, the r part gives no torque about the Earth's center. To compute the effect on the orbit, we smear"the torque r x f about the geosynchronous orbit, and calculate an orbit-average torque(per unit mass) (3) We project onto a set of axis(x, y, z)as defined in the sketch( fixed axis) r=y(=raising cosy here RG is the geosynchronous orbit radius(42, 200 Km)and y=23. 44. is the angle between the equatorial plane and the ecliptic 16.522, Space Propulsion Lecture 4 Prof. Manuel martinez-Sanchez Page 5 of 10
θ r r r Z f R ψ S Fi l x y o ~ G.N. TIDAL TORQUE, SUMMER SOLSTICE VERNAL EQUINOX TIDAL FORCES ES GEOGRAPHICAL NORTH SIDEREAL NORTH EARTH xed Axes: x: A ong Verna y: In Ec z: North s SUN 23 /C IN GEO s l Equinox liptic, solsticial line idereal G G If r is the relative position vector of the spacecraft, and s is the unit vector from the Sun, the Tidal Force (Gravity gradient) per unit mass is JG ( GG G G f = µs ⎡3 r.s s − r⎤ ) (1) 3 R ⎣ Es ⎦ (This the imbalance between solar attraction and centrifugal force, which cancels on G G the Earth, r = 0 ). Of this, the r part gives no torque about the Earth’s center. To G G compute the effect on the orbit, we “smear” the torque r f × about the geosynchronous orbit, and calculate an orbit-average torque (per unit mass): JG 1 2π G G q = × 2π ∫0 r f d θ (2) JG 3µ 1 2π G G G G q = 3 s (r.s ) (r × s) dθ (3) REs 2π ∫0 We project onto a set of axis (x, y, z) as defined in the sketch (fixed axis): ⎧ ⎫ x ⎧cos θ ⎫ G ⎪⎪ ⎪ ⎪ r= ⎨ ⎬y = RG ⎨sin θ cos γ⎬ (4) ⎪ ⎪ ⎪ ⎩ ⎭ ⎩sin θ sin γ ⎪ z ⎭ where RG is the geosynchronous orbit radius (42, 200 Km) and γ = 23.44 D is the angle between the equatorial plane and the ecliptic. 16.522, Space Propulsion Lecture 4 Prof. Manuel Martinez-Sanchez Page 5 of 10