◀ Back to Note

Tony's Intertwingled Memex • Reference Manual

Variable-Density Flywheel Rotational Inertia (Aerospace & Mechanical Engineering)

Formulas & Definitions

An aerospace guidance, navigation, and control (GNC) engineer at the UNLV Howard R. Hughes College of Engineering is designing an energy-storage reaction flywheel for a low-Earth-orbit (LEO) small satellite.

The flywheel is constructed as a flat circular disk of outer radius R=20 cmR = 20\text{ cm}R=20 cm. To maximize rotational kinetic energy while minimizing center-of-mass balancing jitter, the disk is fabricated using a radially graded carbon nanotube matrix whose planar mass density increases quadratically toward the outer rim. In Cartesian coordinates, the area mass density distribution is modeled as:

ρ(x,y)=3(x2+y2)[grams per square centimeter, g/cm2]\rho(x, y) = 3(x^2 + y^2) \quad [\text{grams per square centimeter, g/cm}^2]ρ(x,y)=3(x2+y2)[grams per square centimeter, g/cm2]
over the disk x2+y2≤400x^2 + y^2 \le 400x2+y2≤400.

The satellite's attitude control agility and kinetic storage capacity depend directly on the flywheel's mass moment of inertia (IzI_zIz​) about its central perpendicular spin axis:

Iz=∬D(x2+y2) ρ(x,y) dAI_z = \iint_D (x^2 + y^2) \, \rho(x, y) \, dAIz​=∬D​(x2+y2)ρ(x,y)dA

  1. Express the circular boundary of the flywheel disk in polar coordinates rrr and θ\thetaθ.
  2. Transform the moment of inertia integrand into polar form, incorporating the mass density function ρ(r)\rho(r)ρ(r) and the Jacobian differential area element dA=r dr dθdA = r \, dr \, d\thetadA=rdrdθ.
  3. Use Fubini's product separability theorem to evaluate the double integral, computing the exact moment of inertia in g⋅cm2\text{g}\cdot\text{cm}^2g⋅cm2 and converting the result to standard SI units (kg⋅m2\text{kg}\cdot\text{m}^2kg⋅m2).
  4. Explain how concentrating mass toward the outer perimeter of a reaction flywheel optimizes angular momentum storage per unit total mass in space satellite attitude control systems.