🌱 Tony's Memex Tony's Intertwingled Memex
Digital Garden | OER Curriculum
Home / Digital Garden / Atmospheric Density and Aerodynamic Drag (Aerospace Engineering) Solution

Atmospheric Density and Aerodynamic Drag (Aerospace Engineering) Solution

πŸ–¨οΈ Printable Reference / PDF

Problem Overview

We calculate the dynamic time-rate of change of aerodynamic drag force (dFddt\frac{dF_d}{dt}) on an ascending rocket casing through a three-layer nested system:

Fd(ρ)=450,000ρ,ρ(z)=1.2eβˆ’0.15z,z(t)=0.05t2+2tF_d(\rho) = 450,000\rho, \qquad \rho(z) = 1.2e^{-0.15z}, \qquad z(t) = 0.05t^2 + 2t

Step-by-Step Mathematical Solution

Step 1: Deconstruct the Three-Layer Composition

Because drag force FdF_d depends on air density ρ\rho, density ρ\rho depends on altitude zz, and altitude zz depends on flight time tt, the Chain Rule links all three rates:

dFddt=dFddρ⋅dρdzβ‹…dzdt\frac{dF_d}{dt} = \frac{dF_d}{d\rho} \cdot \frac{d\rho}{dz} \cdot \frac{dz}{dt}

Step 2: Differentiate Each Structural Layer Independently

  1. Outer Layer (Force vs. Density):
    dFddρ=ddρ[450,000ρ]=450,000 N/(kg/m3)\frac{dF_d}{d\rho} = \frac{d}{d\rho}[450,000\rho] = 450,000\text{ N}/(\text{kg/m}^3)
  2. Middle Layer (Density vs. Altitude): Using the Exponential Chain Rule with u=βˆ’0.15zu = -0.15z:
    dρdz=1.2(βˆ’0.15)eβˆ’0.15z=βˆ’0.18eβˆ’0.15zΒ (kg/m3)/km\frac{d\rho}{dz} = 1.2(-0.15)e^{-0.15z} = -0.18e^{-0.15z}\text{ (kg/m}^3)/\text{km}
  3. Inner Layer (Altitude vs. Time / Vertical Velocity):
    dzdt=ddt[0.05t2+2t]=0.1t+2Β km/s\frac{dz}{dt} = \frac{d}{dt}[0.05t^2 + 2t] = 0.1t + 2\text{ km/s}

Step 3: Evaluate Localized Kinematics at t=10Β secondst = 10\text{ seconds}

  1. Altitude z(10)z(10):
    z(10)=0.05(10)2+2(10)=0.05(100)+20=5+20=25Β kmz(10) = 0.05(10)^2 + 2(10) = 0.05(100) + 20 = 5 + 20 = 25\text{ km}
  2. Vertical Velocity dzdt\frac{dz}{dt} at t=10t = 10:
    dzdt∣t=10=0.1(10)+2=1+2=3 km/s\frac{dz}{dt}\Big|_{t=10} = 0.1(10) + 2 = 1 + 2 = 3\text{ km/s}

Step 4: Assemble and Compute the Chain Rule Product

Substitute each evaluated component into the primary three-layer Chain Rule:

dFddt=(450,000)β‹…(βˆ’0.18eβˆ’0.15(25))β‹…(3)\frac{dF_d}{dt} = (450,000) \cdot \left(-0.18e^{-0.15(25)}\right) \cdot (3)
dFddt=(450,000)β‹…(βˆ’0.18eβˆ’3.75)β‹…(3)=βˆ’243,000eβˆ’3.75Β N/s\frac{dF_d}{dt} = (450,000) \cdot (-0.18e^{-3.75}) \cdot (3) = -243,000 e^{-3.75}\text{ N/s}
Numerically:
βˆ’243,000eβˆ’3.75β‰ˆβˆ’243,000(0.0235177)β‰ˆβˆ’5,714.8Β N/s-243,000 e^{-3.75} \approx -243,000(0.0235177) \approx -5,714.8\text{ N/s}

Engineering Conclusion

At t=10Β secondst = 10\text{ seconds}, the sounding rocket is climbing through 25Β km25\text{ km} at 3Β km/s3\text{ km/s}, and aerodynamic drag is decreasing at approximately 5,715Β N/s5,715\text{ N/s}. The negative sign demonstrates that the exponential thinning of atmospheric density over altitude dominates the velocity increase, reducing structural skin-friction drag. Guidance, Navigation, and Control (GNC) engineers track dFddt\frac{dF_d}{dt} to identify when the vehicle has safely cleared "Max-Q" (maximum dynamic pressure, where dFddt=0\frac{dF_d}{dt} = 0), authorizing throttled-up engine thrust.