Visual walkthrough — Atmospheric drag — exponential atmosphere model, orbit decay

3.2.34 · D2 · Physics › Orbital Mechanics & Astrodynamics › Atmospheric drag — exponential atmosphere model, orbit decay

---

Step 0 — Ek arrow (vector) ka matlab kya hai, aur do arrows ko "dot" kaise karte hain

KYA HAI: Ek vector ek aisi quantity hai jisme size aur direction dono hote hain — hum ise ek arrow ki tarah draw karte hain. Hum letter ke upar ek chhota sa hat lagate hain, jaise $\vec{v}$, yeh kehne ke liye ki "yeh ek arrow hai, sirf ek number nahi". Plain letter $v$ (bina arrow ke) sirf iski length — yaani speed — ek single positive number ko batata hai.

YAHAN PEHLE KYO INTRODUCE KAREIN? Steps 3 aur 5 mein hum drag ko ek arrow ki tarah handle karte hain jo peeche ki taraf point karta hai, aur hum poochte hain "ek arrow ka kitna hissa doosre ke saath align hai". Uske liye do ideas chahiye: arrow notation $\vec{v}$, aur dot product.

PICTURE: Ek arrow $\vec{v}$ jiske saath iski length $v$ label ki gayi hai; unit arrow $\hat{v}$ (length bilkul 1) jo pure direction dikhata hai; aur do arrows ek angle par milte hue, jinका dot product ek shadow ki tarah draw kiya gaya hai.

Do facts jinpar hum rely karenge: ek arrow ko apne aap se dot karo toh length squared milti hai, $\vec{v}\cdot \vec{v}=v^2$ (angle $0$, $\cos 0=1$); aur do opposite arrows ek negative dot product dete hain ($\cos 180^{\circ }=-1$). Yahi negative sign exactly woh reason hai kyun drag energy drain karta hai.

---

Step 1 — Satellite actually kya hit kar raha hai?

KYA HAI: Ek satellite mass $m$ ka ek solid body hai jo bahut patli hawa mein move karta hai jiska density $\rho$ hai (density = kitne kilograms gas har cubic metre mein baithe hain). Yeh speed $v$ se move karta hai. Iske aage ki face ka area $A$ hai — woh shadow jo yeh dalega agar tum iske travel direction ke saath light chalao.

KYO: Kisi bhi force se pehle, hum count karte hain ki satellite har second kitni gas se milta hai. Woh gas mass hi woh hai jo push back karta hai.

PICTURE: Ek chhote time $dt$ (ek second ka bahut chhota sa hissa) mein, satellite aage $vdt$ distance move karta hai. Yeh ek tube sweep karta hai: ek cylinder jiska cross-section $A$ aur length $vdt$ hai. Neela tube dekho.

Us tube ke andar pakdi gayi gas mass hai density $\times$ volume:

$d{m}_{air}=\underset{{\text{kg per m}}^3}{\underbrace{\rho }}\times \underset{{\text{m}}^2}{\underbrace{A}}\times \underset{\text{m (tube length)}}{\underbrace{vdt}}$

Har piece apni jagah earn karta hai: $\rho$ volume ko mass mein badalta hai; $A\cdot vdt$ hi tube ka volume hai. $d{m}_{air}$ (gas ka ek chhota parcel) aur $m$ (satellite ki fixed mass) mein fark note karo — dono ki zaroorat hogi.

---

Step 2 — "Gas hit" ko ek force mein badalna

KYA HAI: Force hai momentum given away per second. Momentum hai mass $\times$ velocity. Satellite us swept gas ko roughly apni speed $v$ tak push karta hai.

YAHAn force kyun, energy nahi? Hum woh push-back chahte hain jo satellite feel karta hai — woh ek force hai. Moving mass se force tak pahunchne ka sabse clean raasta Newton ka original form hai: $F=\frac{d(\text{momentum})}{dt}$. Yahi exactly woh tool hai "cheez aa rahi hai aur accelerate ho rahi hai" ke liye.

PICTURE: Red arrow = satellite par push-back, motion ke opposite pointing. Har second, ek fresh slug of gas $d{m}_{air}/dt$ speed $\sim v$ tak pheka jaata hai.

$F\sim \frac{d{m}_{air}}{dt}\cdot v=\left(\rho Av\right)\cdot v=\rho A{v}^2$

Term by term: $\frac{d{m}_{air}}{dt}=\rho Av$ gas-mass-per-second hai (Step 1 se, $dt$ se divide karke); $v$ se multiply karna (woh speed jो hum ise dete hain) mass-rate ko momentum-rate = force mein badal deta hai.

Real bodies apni saari speed gas ko nahi dete — kuch sides ke around slip ho jaata hai. Hum us fudge ko ek dimensionless number $C_D$ (the drag coefficient) mein fold karte hain aur conventional factor $\frac{1}{2}$:

Acceleration paane ke liye hum force ko satellite ki apni mass $m$ se divide karte hain (Newton's $a=F/m$), aur satellite ke fixed properties ko ek symbol $B\equiv \frac{C_{D}A}{m}$ mein group karte hain:

$a_{drag}=\frac{1}{2}B\rho v^2,B\equiv \frac{C_{D}A}{m}\left(\text{units }\frac{{\text{m}}^2}{\text{kg}}\right).$

---

Step 3 — Direction: kyun $v\vec{v}$ hota hai aur kabhi ${\vec{v}}^2$ nahi

KYA HAI: Drag hamesha path ke saath backwards point karta hai. Step 0 ki arrow notation use karte hue, hum ek aisa vector likhna chahte hain jiska length $v^2$ ki tarah badhe lekin jiska arrow $-\vec{v}$ ki direction mein point kare (minus arrow ko opposite direction mein flip karta hai).

KYO: "${\vec{v}}^2$" meaningless hai — tum ek arrow ko square nahi kar sakte aur direction nahi rakh sakte. Trick: $v^2$ ko (ek scalar length $v$) $\times$ (the arrow $\vec{v}$) mein split karo. Step 0 se yaad karo $\vec{v}=v\hat{v}$, isliye $v\vec{v}=v^2\hat{v}$ — size $v^2$, direction $\hat{v}$.

PICTURE: Velocity $\vec{v}$ green mein; drag acceleration red mein, bilkul anti-parallel, iski length $\propto v^2$.

---

Step 4 — Do orbit facts jo hum import karne chahiye (assumed nahi, built kiye)

KYA HAI: Orbit ke badle jaane ki baat karne ke liye, hum pure gravity problem (abhi tak koi drag nahi) se do relations chahiye. Radius/semi-major-axis $a$ ki near-circular orbit ke liye jo ek planet ke around hai jisme $\mu =GM$ hai:

$v=\sqrt{\frac{\mu }{a}}\text{aur}\epsilon =-\frac{\mu }{2a}.$

YEH DUE KYO? Drag speed $v$ par act karta hai; lekin decay ko shrinking size $a$ ke roop mein measure kiya jaata hai. Hume $v$, energy $\epsilon$, aur $a$ ke beech ek bridge chahiye. Yeh seedhe Two-Body Problem se aate hain aur Vis-Viva Equation ki jaan hain.

PICTURE: Ek green circular orbit. Jaise $a$ (radius) shrink hota hai, speed arrow (yellow) lamba hota jaata hai — ek paradox ki pehli jhalak.

$\epsilon =-\mu /2a$ padho: energy negative hai (bound orbit) aur jaise $a$ shrink hota hai aur zyada negative hoti jaati hai — energy khona orbit ko andar ki taraf kheenchta hai. $v=\sqrt{\mu /a}$ padho: chhota $a$ ⇒ bada $v$.

---

Step 5 — Drag kitni tez energy drain karta hai? (Power)

KYA HAI: Power hai energy changed per second. Ek force ke liye, power $=$ force $\cdot$ velocity — aur woh "$\cdot$" Step 0 ka dot product hai. Kyunki drag $\vec{v}$ ko oppose karta hai, woh dot product negative hai — energy bahar nikalti jaati hai.

DOT PRODUCT KYO? Dot product $\vec{a}\cdot \vec{v}$ sirf force ka woh hissa pick out karta hai jo motion ke saath hai — sirf wahi hissa jo work karta hai. Drag fully anti-parallel hai (angle $180^{\circ }$, $\cos 180^{\circ }=-1$), isliye woh maximum negative work karta hai.

PICTURE: Drag arrow (red) aur velocity arrow (green) opposite directions mein point karte hain; opposite arrows ka dot product negative hota hai — ek downhill energy slide ke roop mein dikhaya gaya hai.

$\frac{d\epsilon }{dt}={\vec{a}}_{drag}\cdot \vec{v}=(-\frac{1}{2}B\rho v\vec{v})\cdot \vec{v}=-\frac{1}{2}B\rho v(\vec{v}\cdot \vec{v})=-\frac{1}{2}B\rho v^3.$

Term by term: $\vec{v}\cdot \vec{v}=v^2$ (ek arrow ko apne aap se dot karo toh length squared milti hai — Step 0); bache hue $v$ se multiply karne par $v^3$ milta hai; minus sign kehta hai energy sirf neche hi jaati hai.

---

Step 6 — "Energy lost" ko "orbit shrinks" mein convert karna

KYA HAI: Hamare paas $d\epsilon /dt$ hai. Hum $da/dt$ chahte hain. $\epsilon =-\mu /2a$ use karo aur differentiate karo taaki dono rates ko link kiya ja sake.

DIFFERENTIATE KYO? $\epsilon$ depend karta hai $a$ par. Chain rule hume batata hai: agar $\epsilon$ itna per second drop hota hai, toh $a$ kitna move karta hai? Differentiation woh exact tool hai jo "ek cheez ki rate doosri ki rate ke terms mein" ka jawaab deta hai.

PICTURE: $\epsilon$-vs-$a$ curve. $\epsilon$ mein ek chhota sa neecha slide (red) slope ke through $a$ mein ek chhote andar slide (blue) mein map hota hai.

$\epsilon =-\frac{\mu }{2a}$ ko time ke respect se differentiate karo:

$\frac{d\epsilon }{dt}=\frac{d}{dt}\left(-\frac{\mu }{2a}\right)=\frac{\mu }{2a^2}\frac{da}{dt}.$

Yahan $\frac{d}{da}\left(-\frac{\mu }{2a}\right)=+\frac{\mu }{2a^2}$ ($1/a$ ka derivative $-1/a^2$ hai; do minus signs cancel ho jaate hain), aur extra $\frac{da}{dt}$ chain rule hai.

Ab $d\epsilon /dt$ ke dono expressions ko equal set karo (Step 5 = Step 6):

$\frac{\mu }{2a^2}\frac{da}{dt}=-\frac{1}{2}B\rho v^3.$

---

Step 7 — Orbit facts substitute karo aur clean up karo

KYA HAI: $v^3$ ko $v^2=\mu /a$ use karke replace karo, isliye $v^3=(\mu /a{)}^{3/2}$, aur $da/dt$ ke liye solve karo.

KYO: Hum answer sirf $a$ ke terms mein chahte hain (aur constants), kyunki jaise orbit decay karta hai hum $a$ ko track karte hain.

PICTURE: Algebra ek cancellation cascade ki tarah — $\mu$'s aur $a$'s milke ek single clean $\sqrt{\mu a}$ mein aa jaate hain.

$\frac{\mu }{2a^2}\frac{da}{dt}=-\frac{1}{2}B\rho {\left(\frac{\mu }{a}\right)}^{3/2}.$

Dono sides ko $\frac{2a^2}{\mu }$ se multiply karo:

$\frac{da}{dt}=-B\rho \frac{a^2}{\mu }{\left(\frac{\mu }{a}\right)}^{3/2}=-B\rho \frac{a^2}{\mu }\cdot \frac{{\mu }^{3/2}}{a^{3/2}}.$

Powers collect karo: ${\mu }^{3/2}/\mu ={\mu }^{1/2}$ aur $a^2/a^{3/2}=a^{1/2}$:

---

Step 8 — Runaway edge case (kyun end achanak hota hai)

KYA HAI: $\rho$ constant nahi hai — exponential model se, $\rho ={\rho }_0{e}^{-(h-h_0)/H}$. Jaise $a$ shrink karta hai, altitude $h$ drop hoti hai, isliye $\rho$ exponentially climb karta hai. Use wapas plug in karo: $|da/dt|$ mein $\rho$ hai, isliye jaise orbit girta hai, iska khud ka shrink-rate accelerate hota hai.

YEH CASE KYO DIKHAO? Ek reader jo constant $\rho$ assume karta hai woh ek gentle linear glide predict karta hai. Reality ek nonlinear plunge hai — model ko yeh hide nahi karna chahiye.

PICTURE: Altitude vs time. Saalon tak almost flat, phir end mein ek cliff — "catastrophic reentry".

---

Ek picture mein summary

YEH FIGURE KYA DIKHATA HAI (top-left se bottom-left padho): har box derivation ka ek step hai, aur har arrow ek se doosre tak logical move hai. Box 1 — satellite gas $dm=\rho Avdt$ ka ek tube sweep karta hai (Step 1). Arrow to Box 2 — woh swept mass, speed $v$ tak phenka gaya, drag force $\frac{1}{2}B\rho v^2$ ban jaata hai (Step 2). Arrow to Box 3 — force ko velocity ke saath dot karne se power drain $-\frac{1}{2}B\rho v^3$ milta hai (Step 5). Arrow down to Box 4 — $\epsilon =-\mu /2a$ aur chain rule use karna us energy loss ko $a$ se link karta hai (Step 6). Arrow to Box 5 — $v=\sqrt{\mu /a}$ substitute karo (Step 7). Arrow to Box 6 — decay law $\frac{da}{dt}=-B\rho \sqrt{\mu a}$ bahar aa jaata hai. Final arrow down — exponential $\rho$ feed karna ise ek runaway plunge bana deta hai (Step 8).

Recall Feynman retelling — plain words mein poora walk

Ek satellite har second hawa ka ek patla tube scoop karta hai (Step 1). Us hawa ko side mein shoving karne ki cost ek backward push hai — ek drag force jo speed ke square ke saath badha hai (Steps 2–3). Woh push hamesha motion ke against point karta hai, isliye yeh orbit ki energy ko second by second steadily drain karta rehta hai (Step 5). Ab yahan twist hai: ek circle ke liye, kam energy matlab ek chhota circle, aur chhota circle matlab ek tez satellite (Step 4). Hum "energy baar baar bahar nikalti hai" ko "circle shrink ho raha hai" mein translate karte hain thodi si calculus use karke (Step 6), algebra theek karte hain (Step 7), aur ek clean rule bahar aata hai: orbit ek rate $B\rho \sqrt{\mu a}$ se shrink karta hai. Finally, kyunki jitna neecha girta hai utni mooti hawa (altitude $h=a-R$ girta hai, density blast hoti hai), poori cheez snowball ho jaati hai — saalon ki gentle drift, phir ek achanak fiery plunge (Step 8).

Recall One-line self-test

Decay rate $\rho$ ke proportional kyun hai, aur yeh reentry ko achanak kyun banata hai? ::: $\frac{da}{dt}=-B\rho \sqrt{\mu a}$; $\rho$ exponentially badha hai jaise altitude girta hai, isliye shrink-rate khud ko ek runaway mein feed karta hai.

---

Connections

• Parent topic — poori theory aur worked numbers.
• Two-Body Problem · Vis-Viva Equation — $v=\sqrt{\mu /a}$ ka source.
• Orbital Energy & Specific Mechanical Energy — woh $\epsilon =-\mu /2a$ jo humne differentiate kiya.
• Hydrostatic Equilibrium (Atmospheres & Stars) — jahan se $\rho (h)$ ka exponential aata hai.
• Perturbations in Orbital Mechanics · Reentry Aerodynamics & Ballistic Coefficient — endgame.