Rocket Performance Parallel staging Modular staging Standard atmospheres Orbital decay due to drag Straight-line (no gravity) entry based on atmospheric density 1 2014 David L. Akin - All rights reserved http://spacecraft.ssl.umd.edu
Parallel Staging Multiple dissimilar engines burning simultaneously Frequently a result of upgrades to operational systems General case requires brute force numerical performance analysis 2
Parallel-Staging Rocket Equation Momentum at time t: M = mv Momentum at time t+δt: (subscript b =boosters; c =core vehicle) M = (m m b m c )(v + v) + m b (v V e,b ) + m c (v V e,c ) Assume thrust (and mass flow rates) constant 3
Parallel-Staging Rocket Equation Rocket equation during booster burn V = V e ln mfinal m initial = V min,b +m e ln in,c + m pr,c +m 0,2 m in,b +m pr,b +m in,c +m pr,c +m 0,2 where χ= fraction of core propellant remaining after booster burnout, and where V e = V e,bṁ b +V e,c ṁ c ṁ b +ṁ c = V e,bm pr,b +V e,c (1 χ)m pr,c m pr,b +(1 χ)m pr,c 4
Analyzing Parallel-Staging Performance Parallel stages break down into pseudo-serial stages: Stage 0 (boosters and core) V 0 = V e ln Stage 1 (core alone) V 1 = Subsequent stages are as before ( ) min,b +m in,c +χm pr,c +m 0,2 m in,b +m pr,b +m in,c +m pr,c +m 0,2 V e,c ln m in,c +m 0,2 m in,c + m pr,c +m 0,2 5
Parallel Staging Example: Space Shuttle 2 x solid rocket boosters (data below for single SRB) Gross mass 589,670 kg Empty mass 86,183 kg Ve 2636 m/sec Burn time 124 sec External tank (space shuttle main engines) Gross mass 750,975 kg Empty mass 29,930 kg Ve 4459 m/sec Burn time 480 sec Payload (orbiter + P/L) 125,000 kg 6
Shuttle Parallel Staging Example V e = V e,b = 2636 m sec χ = 480 124 480 V e,c =4459 m sec =0.7417 2636(1, 007, 000) + 4459(721, 000)(1.7417) 1, 007, 000 + 721, 000(1.7417) =2921 m sec V 0 = 2921 ln V 1 = 4459 ln 862, 000 3, 062, 000 =3702m sec 154, 900 689, 700 =6659m sec V tot =10, 360 m sec 7
Modular Staging Use identical modules to form multiple stages Have to cluster modules on lower stages to make up for nonideal ΔV distributions Advantageous from production and development cost standpoints 8
Module Analysis All modules have the same inert mass and propellant mass Because δ varies with payload mass, not all modules have the same δ! Use module-oriented parameters ε Conversions m in m in + m pr ε = δ 1 λ σ m in m pr σ = δ 1 δ λ 9
Rocket Equation for Modular Boosters Assuming n modules in stage 1, r 1 = n(m in )+m o2 n(m in + m pr )+m o2 = nε + mo2 m mod n + m o2 m mod If all 3 stages use same modules, n j for stage j, where r 1 = n 1ε + n 2 + n 3 + ρ pl n 1 + n 2 + n 3 + ρ pl ρ pl m pl ; m mod = m in + m pr m mod 10
Example: Conestoga 1620 (EER) Small launch vehicle (1 flight, 1 failure) Payload 900 kg Module gross mass 11,400 kg Module empty mass 1,400 kg Exhaust velocity 2754 m/sec Staging pattern 1st stage - 4 modules 2nd stage - 2 modules 3rd stage - 1 module 4th stage - Star 48V (gross mass 2200 kg, empty mass 140 kg, V e 2842 m/sec) 11
Conestoga 1620 Performance 4th stage V V 4 = V e4 ln m f4 = 2842 ln m o4 Treat like three-stage modular vehicle; M pl =3100 kg = m in m mod = 1400 11400 =0.1228 pl = 900 + 140 900 + 2200 = 3104 m sec m pl = 3100 m mod 11400 =0.2719 n 1 = 4; n 2 = 2; n 3 =1 12
Constellation 1620 Performance (cont.) r 1 = n 1 + n 2 + n 3 + pl = n 1 + n 2 + n 3 + pl r 2 = n 2 + n 3 + pl = n 2 + n 3 + pl r 3 = n 3 + pl = n 3 + pl 13 4 0.1228 + 2 + 1 + 0.2719 4+2+1+0.2719 2 0.1228 + 1 + 0.2719 2+1+0.2719 1 0.1228 + 0.2719 1+0.2719 V 1 = 1814 m sec ; V 2 = 2116 m sec V 3 = 3223 m sec ; V 4 = 3104 m V total = 10, 257 m sec sec =0.5175 =0.4638 =0.3103
Discussion about Modular Vehicles Modularity has several advantages Saves money (smaller modules cost less to develop) Saves money (larger production run = lower cost/ module) Allows resizing launch vehicles to match payloads Trick is to optimize number of stages, number of modules/stage to minimize total number of modules Generally close to optimum by doubling number of modules at each lower stage Have to worry about packing factors, complexity 14
OTRAG - 1977-1983 15
Modular Example Let s build a launch vehicle out of seven Space Shuttle Solid Rocket Boosters M in =86,180 kg M pr =503,500 kg ε m in m in + m pr =0.1461 σ m in m pr =0.1711 Look at possible approaches to sequential firing 16
Modular Sequencing - SRB Example Assume no payload All seven firing at once - ΔV tot =5138 m/sec 3-3-1 sequence - ΔV tot =9087 m/sec 4-2-1 sequence - ΔV tot =9175 m/sec 2-2-2-1 sequence - ΔV tot =9250 m/sec 2-1-1-1-1-1 sequence - ΔV tot =9408 m/sec 1-1-1-1-1-1-1 sequence - ΔV tot =9418 m/sec Sequence limited by need to balance thrust laterally 17
Atmospheric Density with Altitude Ref: V. L. Pisacane and R. C. Moore, Fundamentals of Space Systems Oxford University Press, 1994 18
Energy Loss Due to Atmospheric Drag Drag D 1 2 v2 Ac D Drag acceleration a d = D m = m c D A v2 2 Ac D m <== Ballistic Coefficient a d = v2 2 orbital energy E = µ 2a de dt = µ da 2a 2 dt 19
Energy Loss Due to Atmospheric Drag Since drag is highest at perigee, the first effect of atmospheric drag is to circularize the orbit (high perigee drag lowers apogee) de drag dt = a d v v 2 circ = µ a de drag dt = v2 2 µ a de drag dt = µ a 2 µ a = µ 3 2 a 2 20
Derivation of Orbital Decay Due to Drag Set orbital energy variation equal to energy lost by drag = o e µ da 2a 2 dt = 2 µ a da dt = µa 3 2 h hs a = h + r E = da dt = dh dt dh dt = µ (h + r E ) o e h hs 21
Derivation of Orbital Decay (2) This is a separable differential equation... 1 re + h e h hs dh = µ o dt h re + h e h hs dh = h o 1 µ o t t o dt Assume 1 re h r E + h h o e h hs dh = re for r E h µ o (t t o ) 22
Derivation of Orbital Decay (3) h s e h hs e h o hs re = µ o (t t o ) e h hs h(t) = h s ln e h o hs Note that some variables typically use km, and others are in meters - you have to make sure unit conversions are done properly to make this work out correctly! e h o hs = 23 µre o (t t o ) h s µre o (t t o ) h s
Orbit Decay from Atmospheric Drag 250 200 Altitude (km) 150 100 50 β=500 β=1500 β=5000 0 0 10000 20000 30000 40000 Time (sec) 24
Time Until Orbital Decay e h hs e h o hs = µre h s o (t t o ) To find the time remaining (t o =0) until the orbit reaches any given critical altitude, some algebra gives t(h crit )= h s p e h o hs µre o h crit e hs t(h crit ) / 25
Decay Time to r=120 km 600 500 Altitude (km) 400 300 200 100 β=500 β=1500 β=5000 0 0 20 40 60 80 100 Decay Time (yrs) 26
(no lift) s = distance along the flight path dv dt = g sin D m dv dt = dv ds ds dt = V dv ds = 1 2 d(v 2 ) ds D mg horizontal v, s 1 2 d(v 2 ) ds = g sin D m Drag D 1 2 v2 Ac D 1 2 sin 2 d(v 2 ) ds = g sin v 2 2m Ac D d(v 2 ) dh = g sin v 2 2m Ac D ds dh ds = dh sin 27
(2) Exponential atmosphere = o e h hs d o = e h hs dh h s = o e o h hs dh h s = o dh h s dh = hs d sin 2 sin 2 d(v 2 ) d d(v 2 ) dh = g sin v 2 2m Ac D v 2 = g sin 2 h s Ac D m d(v 2 ) d = 2gh s 28 + h sv 2 sin Ac D m
(3) Let m c D A Ballistic Coe cient d(v 2 ) d h s sin v2 = 2gh s Assume mg D to get homogeneous ODE Use d(v 2 ) d h s sin v2 = 0 v 2 as integration variable d(v 2 ) v 2 = h s sin d v v e d(v 2 ) v 2 = h s sin 0 d v e = velocity at entry 29
(4) Note that the effect of ignoring gravity is that there is no force perpendicular to velocity vector constant flight path angle γ straight line trajectories ln v2 v 2 e = 2 ln v v e = h s sin v v e = exp h s 2 sin v v e = exp h s o 2 sin o Check units: m kg m 3 kg m 2 30
Earth Entry, γ=-60 v/v e 1 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0 0 0.2 0.4 0.6 0.8 1 / o Beta=100 kg/m^3 300 1000 3000 10000 31
What About Peak Deceleration? To find n max, set d dt d 2 v dt 2 = 1 2 d 2 v dt 2 = 1 2 n dv dt = dv dt v2 2 = d2 v dt 2 = 0 2v dv dt 2 2 v 3 2 + v2 d dt + v 2 d dt = 0 = 0 2 v 3 = v 2 d dt 2 v = d dt 32
Peak Deceleration (2) From exponential atmosphere, d dt = From geometry, dh dt = v sin d dt = o h s e v h s sin 2 v = h hs dh dt = dh h s dt 2 v = v sin h s d dt Remember that this refers to the conditions at max deceleration nmax = hs sin 33
Critical β for Deceleration Before Impact At surface, = o Value of at which vehicle hits crit = oh s sin ground at point of maximum deceleration How large is maximum deceleration? dv dt = v2 2 dv = v2 dt max 2 dv = n v2 max dt max 2 sin h s = 1 2 v 2 h s sin Note that this value of v is actually v nmax 34
Peak Deceleration (3) From page 14, v nmax v e dv dt max = 1 2 v v e = exp 35 h s 2 sin h s = exp sin 2 sin h s 2 v e e 1 2 = e 1 2 h s sin = v2 e sin 2h s e Note that the velocity at which maximum deceleration occurs is always a fixed fraction of the entry velocity - it doesn t depend on ballistic coefficient, flight path angle, or anything else! Also, the magnitude of the maximum deceleration is not a function of ballistic coefficient - it is dependent on the entry trajectory (v e and γ) but not spacecraft parameters (i.e., ballistic coefficient).
Terminal Velocity Full form of ODE - d v 2 d h s sin v2 = 2gh s At terminal velocity, v = constant v T h s sin v2 T = 2gh s v 2 T = 2g sin 36
Cannon Ball γ=-90 6.75 diameter sphere, c D =0.2, V E =6000 m/sec Iron Aluminum Balsa Wood Weight 40 lb 15.6 lb 14.5 oz β (kg/m 2 ) 3938 1532 89 ρmd (kg/m 3 ) 0.555 0.216 0.0125 hmd (m) 5600 12,300 32,500 Vimpact (m/s) 1998 355 0* Vterm (m/sec) 251 156 38 *Artifact of assumption that D 37 mg
Atmospheric Density with Altitude Pressure=the integral of the atmospheric density in the column above the reference area = f(h) P o = Z 1 o gdh = o g Z 1 o e h hs dh = o gh s he h hs = o gh s [0 1] i 1 o Earth: o = 1.226 kg m 3 ; h s = 7524m; P o = o gh s P o (calc) = 90, 400 Pa; P o (act) = 101, 300 Pa o, P o 38
Nondimensional Ballistic Coefficient v v e = exp h s o 2 sin o Po =exp 2 g sin Let b = g (Nondimensional form of ballistic coe cient) o h s P o Note that we are using the estimated value of P o = o gh s, not the actual surface pressure. v 1 = exp v e 2 sin o o crit = oh s sin crit = 1 sin 39
Entry Velocity Trends, γ=-90 Density Ratio 0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 Velocity Ratio 0 0.2 0.4 0.6 0.8 1 0.03 0.1 0.3 1 3 40