3.4.22 · D3 · Physics › Rocket Flight Mechanics › Thermal protection systems — ablators (PICA, SLA), metallic
Intuition Yeh page kyun hai
Parent note ne tumhe teen formulas diye. Lekin ek formula jo sirf "beech mein" chalti ho, usska koi fayda nahi — reentry extreme inputs leke aati hai: bahut chota nose radius, near-vacuum densities, itni speed jahan radiation suddenly dominate karne lagti hai, aur re-radiation terms jo recession rate ko zero kar dete hain. Yeh page har corner case walk karta hai taaki koi bhi scenario tumhe unprepared na pakde. Hum sirf wahi teen tools reuse karenge jo parent mein hain, aur unhe neeche dobara define kar rahe hain taaki kuch bhi assume na ho.
Koi bhi number se pehle, teen tools ko simple words mein dobara state karte hain taaki har symbol earn kiya gaya ho .
k ki hidden units (taaki answer sach mein W/m² mein aaye)
Bracketed pieces units carry karte hain: ρ / R n ki units hain ( kg m − 3 ) / m = kg 1/2 m − 2 , aur v 3 ki units hain m 3 s − 3 . In dono ko multiply karne par milta hai kg 1/2 m s − 3 . Ise W/m 2 = kg s − 3 mein convert karne ke liye, constant k khud units carry karta hai
[ k ] = kg 1/2 m s − 3 kg s − 3 = kg 1/2 m − 1 .
Toh k ke "SI units" ka matlab hai k = 1.7 × 1 0 − 4 kg 1/2 m − 1 . Yahi reason hai ki k ek empirical constant hai: yeh silently woh half-power of mass aur inverse metre absorb kar leta hai jo physics scaling khud fix nahi kar sakti.
Tool 2 — Re-radiation (Stefan–Boltzmann). Ek garam surface heat wapas bahar phenk deti hai:
q ˙ rad,out = ε σ T w 4
ε = emissivity, ek pure number 0 se (shiny, poorly radiate karta hai) lekar 1 tak (perfect black radiator).
σ = 5.67 × 1 0 − 8 W m − 2 K − 4 = Stefan–Boltzmann constant.
T w = wall temperature kelvin (K) mein. Dekho Radiative Heat Transfer and Stefan–Boltzmann Law .
Tool 3 — Ablator energy balance. Jo aata hai, minus jo re-radiate hota hai, woh material khaane mein jaata hai:
s ˙ = ρ abl Q ∗ q ˙ conv − ε σ T w 4
s ˙ = recession rate : surface kitni tezi se andar khaati jaati hai, m/s mein.
ρ abl = ablator ki density (kg/m³).
Q ∗ = effective heat of ablation (J/kg) — ek kilogram shield hatane mein kitni energy lagti hai.
Intuition Jab Tool 1 khud break ho jaata hai (radiation-dominated entry)
Sutton–Graves (Tool 1) ek convective correlation hai: yeh assume karta hai ki heat wall tak garm gas ke touch karne se pahunchi. Extreme speeds par (roughly v > 10 km/s , Mach numbers ~25 ke baad — lunar aur Mars return) shock gas itna garam ho jaata hai ki woh glow karta hai, aur ek doosra channel khulta hai: gas heat seedha wall par radiate karta hai. Woh radiative input q ˙ rad,in Tool 1 mein bilkul nahi hai, aur yeh roughly ρ a R n b v c ke saath scale hota hai jisme v ka power bahut steep hota hai (often c ≈ 8 –12 , convective v 3 se kaafi upar). Toh bahut high speed par, convective correlations real load ko under-predict karti hain aur tumhe radiative term add karna padta hai — yeh exactly woh regime hai jiske edge par Ex 4 hai, aur isliye PICA (ek high-Q ∗ ablator) wahan mandatory hai. Dekho Radiative Heat Transfer and Stefan–Boltzmann Law aur Reentry Aerothermodynamics .
Is page ka har problem neeche diye table ka ek cell hai. Agar koi cell cover nahi hua, toh ek din tum use unprepared milte — isliye hum sab cover kar rahe hain.
#
Cell class
Kya extreme / tricky hai
Example
C1
Baseline convective heating
ordinary LEO capsule
Ex 1
C2
Chota R n limit (sharp nose)
R n → small → q ˙ → huge
Ex 2
C3
Bada R n / density-zero limit
R n → ∞ ya ρ → 0 → q ˙ → 0
Ex 3
C4
Speed scaling (v 3 sensitivity, radiation onset)
v double → ×8 heat; v > 10 km/s Tool 1 break karta hai
Ex 4
C5
Radiative equilibrium wall temp
tile survival, T 4 invert karo
Ex 5
C6
Ablator recession (saare terms live)
re-radiation subtract karo, s ˙ nikalo
Ex 6
C7
Degenerate ablation (s ˙ → 0 )
re-radiation ≥ incoming → koi recession nahi
Ex 7
C8
Real-world word problem (mass budget)
Mars entry ke liye total shield mass
Ex 8
C9
Exam twist (kaun sa TPS + kyun)
T w melting se zyada → ablator force karta hai
Ex 9
Poore page mein use hone wale constants (inhe yaad karo):
k = 1.7 × 1 0 − 4 kg 1/2 m − 1 , σ = 5.67 × 1 0 − 8 W m − 2 K − 4 .
Intuition Matrix figure ko kaise padhein
Yeh map Tool 1 ko picture ke roop mein dikhata hai. Horizontal axis = nose radius R n (right ki taraf blunter); vertical axis = speed v km/s mein; colour = log 10 q ˙ (colour bar exponent hai, toh + 1 ka jump matlab das guna zyada heat). Kyunki q ˙ ∝ 1/ R n , left slide karna (sharper) quickly bright kar deta hai — yeh hai Ex 2. Kyunki q ˙ ∝ v 3 , upar slide karna (faster) aur bhi tezi se bright karta hai — yeh hai Ex 4. Chaar coloured pins hamare examples hain: white Ex1 baseline map ke beech mein; red Ex2 sharp brightest zone mein left wall se chipka; green Ex3 blunt neeche-right ke dark cool corner mein chhupta; yellow Ex4 fast speed axis par upar ride karta. Har worked example neeche bas "is pin par jao aur heat padho" hai.
Worked example Ex 1 — C1 Baseline: ek LEO capsule
Ek capsule v = 7.5 km/s par reenter karta hai jahan air density ρ = 1.0 × 1 0 − 3 kg/m 3 hai. Nose radius R n = 1.0 m , k = 1.7 × 1 0 − 4 . q ˙ conv nikalo.
Forecast: Pehle order of magnitude guess karo — millions of watts per square metre? Tens of thousands? Aage padhne se pehle apna guess likho.
Speed ko SI mein convert karo: v = 7500 m/s .
Yeh step kyun? Formula SI mein hai; km/s silently answer ko 100 0 3 se multiply kar deta.
Density/radius factor compute karo: ρ / R n = 1 0 − 3 /1.0 = 1 0 − 3 = 0.03162 .
Yeh step kyun? Yeh "kitne molecules, kitne spread out" wala part hai — ise isolate karo taaki v ka cube alag rahe.
Speed cube karo: v 3 = ( 7500 ) 3 = 4.219 × 1 0 11 .
Yeh step kyun? Speed sabse dominant lever hai; ise sahi karna sabse zyada important hai.
Multiply karo: q ˙ conv = 1.7 × 1 0 − 4 × 0.03162 × 4.219 × 1 0 11 = 2.27 × 1 0 6 W/m 2 .
Answer: q ˙ conv ≈ 2.27 MW/m 2 .
Verify: Units: [ k ] ⋅ kg m − 3 / m ⋅ ( m/s ) 3 = kg 1/2 m − 1 ⋅ kg 1/2 m − 2 ⋅ m 3 s − 3 = kg s − 3 = W/m 2 . ✓ Sanity: kuch MW/m 2 LEO return ke liye textbook figure hai — dekho Reentry Aerothermodynamics . Tumhara "millions" ka forecast sahi tha.
Worked example Ex 2 — C2 Sharp-nose limit: capsules round kyun hote hain
Ex 1 lo lekin nose ko R n = 0.01 m (1 cm sharp edge, jaise Shuttle wing leading edge ) tak shrink karo. Same ρ , same v . Naya q ˙ aur Ex 1 se ratio nikalo.
Forecast: R n ×100 gira. Kyunki q ˙ ∝ 1/ R n , compute karne se pehle heating multiplier guess karo.
Naya factor: ρ / R n = 1 0 − 3 /0.01 = 0.1 = 0.3162 .
Yeh step kyun? Sirf R n badla; Ex 1 ki baaki sab cheez reusable hai.
Ex 1 se ratio: 1.0/0.01 = 100 = 10 .
Yeh step kyun? q ˙ ∝ 1/ R n , toh ×100 chota radius 100 = 10 × heat deta hai — yeh limit dekhne layak hai.
Naya heating: q ˙ = 10 × 2.27 MW/m 2 = 2.27 × 1 0 7 W/m 2 = 22.7 MW/m 2 .
Answer: 22.7 MW/m 2 — blunt case ka das guna .
Verify: R n → 0 push karo (truly sharp): 1/ R n → ∞ , toh q ˙ → ∞ . Yeh mathematical proof hai ki razor nose roast ho jaata hai — exactly isliye leading edges ko RCC chahiye, silica tiles nahi.
Worked example Ex 3 — C3 Vanishing limits: high altitude aur bada shield
Do degenerate checks. (a) Vehicle abhi bhi near-vacuum mein hai, ρ = 1.0 × 1 0 − 8 kg/m 3 (R n = 1 m, v = 7500 m/s). (b) Ek giant inflatable shield, R n = 100 m (wapas ρ = 1 0 − 3 ).
Forecast: Dono mein heating choti honi chahiye. 1 0 5 ke factor ke change ke liye — density zyada tezi se giregi ya radius?
(a) Density factor: 1 0 − 8 /1 = 1 0 − 4 . Phir q ˙ = 1.7 × 1 0 − 4 × 1 0 − 4 × 4.219 × 1 0 11 = 7.17 × 1 0 3 W/m 2 .
Yeh step kyun? Top-of-atmosphere limit dikhata hai: almost koi molecule nahi, almost koi heat nahi. Jaise ρ → 0 , q ˙ → 0 smoothly.
(b) Radius factor: 1 0 − 3 /100 = 1 0 − 5 = 3.162 × 1 0 − 3 . Phir q ˙ = 1.7 × 1 0 − 4 × 3.162 × 1 0 − 3 × 4.219 × 1 0 11 = 2.27 × 1 0 5 W/m 2 .
Yeh step kyun? Jaise R n → ∞ , q ˙ → 0 — mathematical reason "blunter is cooler."
Answer: (a) ≈ 7.2 kW/m 2 ; (b) ≈ 0.227 MW/m 2 .
Verify: Density root ke neeche linearly enter hoti hai (ρ ); ρ mein 1 0 − 5 ki drop 1 0 − 5 ≈ 3 × 1 0 − 3 deti hai — same factor jo 1/ R n mein 1 0 − 5 ki drop deta. Toh symmetry hold karti hai: jab bhi koi variable gentle extreme ki taraf jaata hai, dono channels shut off ho jaate hain. Case (a) case (b) se 1000 × chota hai kyunki ρ 1 0 5 gira jabki R n sirf 1 0 2 badha.
Worked example Ex 4 — C4 Speed scaling:
v 3 lever (aur jahan Tool 1 break hota hai)
Ek lunar-return capsule LEO ke 7.5 km/s ki jagah v = 11 km/s hit karta hai (same ρ = 1 0 − 3 , R n = 1 ). Convective heating kitne factor se jump karta hai, aur uski value kya hai? Phir caveat note karo.
Forecast: Speed 11/7.5 ≈ 1.47 badhi. v 3 ki wajah se, heat multiplier guess karo.
Speed ratio: 11/7.5 = 1.4667 .
Yeh step kyun? Sirf woh isolate karo jo badla.
Cube karo: 1.466 7 3 = 3.155 .
Yeh step kyun? q ˙ ∝ v 3 , toh heating speed ratio ke cube ke saath scale hoti hai — reentry mein sabse important sensitivity.
Naya convective heating: q ˙ conv = 3.155 × 2.27 MW/m 2 = 7.16 MW/m 2 .
Answer: ≈ 3.16 × zyada, ≈ 7.16 MW/m 2 — sirf convective .
Verify: Sirf 47% tezi se triple heat ho gai. Caveat (radiation edge case): 11 km/s par hum 10 km/s ke baad hain, toh glowing-gas radiative input q ˙ rad,in jo Tool 1 omit karta hai ab matter karta hai aur convective value ke barabar ya zyada ho sakta hai. 7.16 MW/m 2 isliye ek lower bound hai; actual load zyada hai, aur exactly isliye deep-space returns (dekho Specific Impulse and Energy Budgets ) high-Q ∗ PICA ablator maangti hain. Direct check: ( 11000 ) 3 / ( 7500 ) 3 = 3.155 . ✓
Worked example Ex 5 — C5 Radiative equilibrium: kya silica tile survive karegi?
Ek reusable silica tile Shuttle belly par q ˙ conv = 3.0 × 1 0 5 W/m 2 face karti hai. Uski coating ki emissivity ε = 0.85 hai. Agar saari incoming heat re-radiate ho, toh steady wall temperature T w nikalo. (Tiles ~1500 K ke upar fail hoti hain.)
Forecast: Balance kehta hai q ˙ conv = ε σ T w 4 . Hum fourth power invert karenge — kya T w 1500 K ke neeche ya upar aayega?
Balance likho: incoming = outgoing radiation, q ˙ conv = ε σ T w 4 .
Yeh step kyun? "Steady" ka matlab hai tile tab heating band karti hai jab jo radiate hota hai woh arrive hone ke barabar ho — energy in = energy out.
T w ke liye solve karo: T w = ( ε σ q ˙ conv ) 1/4 .
Fourth root kyun? Radiation T 4 ke saath badhta hai; undo karne ke liye fourth root lete hain — inverse operation, bilkul jaise square undo karta hai.
Plug in karo: 0.85 × 5.67 × 1 0 − 8 3.0 × 1 0 5 = 6.226 × 1 0 12 . Fourth root: T w = ( 6.226 × 1 0 12 ) 1/4 = 1580 K .
Yeh step kyun? Inverted formula mein numbers substitute karo aur fourth root actually lo — yeh abstract balance ko ek concrete kelvin value mein convert karta hai jo hum 1500 K survival limit se compare kar saken.
Answer: T w ≈ 1580 K .
Verify: 1500 K limit se thoda upar — yeh belly location tile ke survival edge ke kaafi paas hai, jo exactly real design margin ki kahani hai. Fourth root par sanity: ε ko 0.85 se 1 ki taraf badhao toh T w girta hai (better radiator cooler chalta hai), jo high-ε black coating choice confirm karta hai.
Worked example Ex 6 — C6 Full ablator: re-radiation live ke saath recession
Ek PICA shield (ρ abl = 270 kg/m 3 , Q ∗ = 1.2 × 1 0 7 J/kg , ε = 0.9 ) char temperature T w = 3000 K par q ˙ conv = 8.0 × 1 0 6 W/m 2 ke neeche hai. Recession rate s ˙ nikalo.
Forecast: Pehle kuch heat re-radiate hoti hai; baaki material khaata hai. Guess karo ki 8 MW/m² mein se re-radiation bada ya chota slice hai.
Re-radiated flux: ε σ T w 4 = 0.9 × 5.67 × 1 0 − 8 × ( 3000 ) 4 . Kyunki 300 0 4 = 8.1 × 1 0 13 , yeh hai 4.135 × 1 0 6 W/m 2 .
Yeh step kyun? Char white-hot glow karta hai; incoming heat ka bada chunk light ke roop mein nikal jaata hai, mass ablate hone se pehle hi.
Ablation ke liye net: q ˙ conv − ε σ T w 4 = 8.0 × 1 0 6 − 4.135 × 1 0 6 = 3.865 × 1 0 6 W/m 2 .
Yeh step kyun? Sirf woh heat jo re-radiate nahi hoti, woh dying material ke zariye carry-off hoti hai.
Recession: s ˙ = 270 × 1.2 × 1 0 7 3.865 × 1 0 6 = 3.24 × 1 0 9 3.865 × 1 0 6 = 1.193 × 1 0 − 3 m/s .
Yeh step kyun? Net heat ko (shield ki energy per m³) = ρ abl Q ∗ se divide karo taaki metres per second mein recession mile.
Answer: s ˙ ≈ 1.19 × 1 0 − 3 m/s ≈ 1.19 mm/s .
Verify: Re-radiation ne incoming heat ka 4.1/8.0 ≈ 52% khaya — aadhe se zyada — isliye hot char ek genuine third defence hai, footnote nahi. Units: ( W/m 2 ) / [( kg/m 3 ) ( J/kg )] = ( J s − 1 m − 2 ) / ( J m − 3 ) = m/s . ✓
Worked example Ex 7 — C7 Degenerate ablation: recession zero par ruk jaata hai
Same PICA jaise Ex 6, lekin reentry ke end mein heating q ˙ conv = 4.0 × 1 0 6 W/m 2 par drop ho gayi hai jabki char abhi bhi T w = 3000 K par glow kar raha hai. s ˙ nikalo.
Forecast: Ex 6 se re-radiation term 4.135 MW/m 2 hai. Incoming 4.0 MW/m 2 hai. Kya hota hai jab shield zyada radiate kare jितना arrive ho raha hai?
Ablation ke liye net: 4.0 × 1 0 6 − 4.135 × 1 0 6 = − 1.35 × 1 0 5 W/m 2 (negative).
Yeh step kyun? Tool 3 mein bracket negative ja sakta hai — yeh woh degenerate case hai jo parent ne kabhi nahi dikhaya.
Negative ko interpret karo: negative "recession" unphysical hai — mass un-ablate nahi ho sakta. Physical outcome hai s ˙ = 0 ; surface receding band ho jaati hai aur cool hone lagti hai (actual T w girta hai jab tak radiation reduced input ko balance na kare).
Yeh step kyun? Hum zero par clamp karte hain kyunki formula ne steady ablation assume ki thi, jo tab valid nahi hoti jab re-radiation alone saari incoming heat dump kar sake.
Answer: s ˙ = 0 (ablation ruk jaati hai; shield ab re-radiating tile ki tarah behave karta hai).
Verify: Yeh q ˙ conv ≤ ε σ T w 4 ka physical meaning hai: us threshold ke neeche koi material consume nahi hota. Yeh parent ka steel-man mistake "reusable tiles ablators se hamesha better hain" bhi reinforce karta hai — parent ka fix tha ki high flux par koi bhi reusable solid survive nahi kar sakta, isliye ablation required hai; yeh example flip side dikhata hai, ki low flux par ablator bhi ablate karna band kar deta hai aur reusable-jaisa behave karta hai. Dono milke prove karte hain ki TPS choice flux peak se govern hoti hai, average se nahi.
Worked example Ex 8 — C8 Real-world word problem: Mars shield mass
Ek Mars entry vehicle (dekho Mars Entry Descent and Landing ) SLA-561V shield use karta hai, ρ abl = 260 kg/m 3 . Engineers ko 4.5 m diameter ke circular heatshield par thickness 5.0 cm chahiye. Ablator ki mass kitni hogi?
Forecast: 4.5 m disc, 5 cm thick, roughly quarter-water density par. Guess karo: tens of kg? Hundreds?
Shield frontal area: A = π r 2 = π ( 2.25 ) 2 = 15.90 m 2 .
Yeh step kyun? Ablator front face par ek layer hai; uska volume = area × thickness.
Volume: V = A × t = 15.90 × 0.05 = 0.7952 m 3 .
Yeh step kyun? 5 cm ko pehle 0.05 m mein convert karo taaki saari lengths SI mein rahein.
Mass: m = ρ abl V = 260 × 0.7952 = 206.8 kg .
Yeh step kyun? Mass = density × volume — density ki definition rearranged.
Answer: ≈ 207 kg ablator.
Verify: Kuch tonnes ke lander ke liye, ~200 kg heatshield ek believable single-digit-percent mass fraction hai — real EDL budgets se match karta hai. Units: ( kg/m 3 ) ( m 2 ) ( m ) = kg . ✓
Worked example Ex 9 — C9 Exam twist: kaun sa TPS, aur numbers se prove karo
Ek wing leading edge q ˙ conv = 1.5 × 1 0 6 W/m 2 dekh raha hai. Ek candidate reusable ceramic ε = 0.8 ke saath radiate karta hai lekin 2000 K par melt ho jaata hai. Pure re-radiation use karte hue, decide karo ki ceramic survive karega ya nahi; agar nahi, toh TPS class ka naam batao jo ki isski jagah use hogi aur kyun.
Forecast: Radiative-equilibrium T w compute karo aur 2000 K se compare karo. Upar → material melt ho jaata hai.
Equilibrium temp: T w = ( ε σ q ˙ conv ) 1/4 = ( 0.8 × 5.67 × 1 0 − 8 1.5 × 1 0 6 ) 1/4 .
Yeh step kyun? Same balance jaise Ex 5 — ek reusable material tab survive kar sakta hai jab uska equilibrium temperature uske melting point ke neeche rahe.
Bracket ke andar: 4.536 × 1 0 − 8 1.5 × 1 0 6 = 3.307 × 1 0 13 . Fourth root: T w = 2398 K .
Yeh step kyun? T 4 ko fourth root se invert karo, phir numbers substitute karke concrete kelvin value lo jise 2000 K melting point se compare kar sakein.
Compare karo: 2398 K > 2000 K → ceramic melt ho jaata hai. Isliye ek reusable radiator is spot par survive nahi kar sakta. Ya toh RCC (carbon–carbon, >1900 K rated — Shuttle nose/leading-edge material) use karo ya ek ablator , jiska m ˙ Q ∗ blocking term woh energy remove karta hai jo koi solid re-radiate nahi kar sakta.
Yeh step kyun? Yeh parent ka steel-man mistake quantitative banata hai: high flux par, ( q ˙ / ε σ ) 1/4 har melting point se zyada ho sakta hai, jo ablation ya highest-temperature reusable material force karta hai.
Answer: T w ≈ 2398 K , 2000 K se upar → ceramic fail; RCC ya ablator choose karo.
Verify: Yeh exactly Space Shuttle Columbia Accident ka lesson hai — leading edges RCC lete hain precisely isliye kyunki ordinary tiles ka equilibrium temperature unki limit se zyada ho jaata hai. 2398 K > 2000 K number bina kisi hand-waving ke decision settle kar deta hai.
Recall Matrix ko self-test karo
Kaun se cell mein fourth-root chahiye? ::: C5 aur C9 — ε σ T w 4 ko invert karke T w nikalna.
Kaun sa cell q ˙ → ∞ karta hai? ::: C2, jaise R n → 0 (sharp nose).
Kaun sa cell s ˙ ko zero par clamp karta hai? ::: C7, jab re-radiation ≥ incoming heat ho.
Speed double karne par heating kitne se multiply hoti hai? ::: 2 3 = 8 (v 3 law, cell C4).
Kis speed ke upar Tool 1 akela enough nahi rehta? ::: ~10 km/s (Mach >25), jahan glowing-gas radiative input oopar se add ho jaata hai.
Mnemonic Koi bhi TPS problem padhne ka tarika
"Cube the speed, root the rest, radiate before you ablate." — v cubed enter hoti hai, ρ / R n root ke neeche enter hota hai, aur recession compute karne se pehle hamesha ε σ T w 4 subtract karo.