Accretion Energy

I compute here the energy released in the assembly of material from infinity to form a spherical ruble pile. This is also the energy required to disperse it again.

Let ρ be the pile’s constant density. Let m be the mass of the pile which will increase as the pile accretes. When a new mass element, dm, falls in from infinity the attraction is f = G m dm/r2 while the element is still at distance r. The volume of the pile is 4/3πr03 = m/ρ where r0 is the pile’s radius.
It will reach the surface at radius r0 = (4/3πρ/m)–1/3 having released energy = de =
r0f dr = G m dm/r0 = (4/3πρ/m)1/3 G m dm = G dm (4/3πρ)1/3 m2/3.
As the mass increases the total energy is
0Mde = ∫0MG dm (4/3πρ)1/3 m2/3 = G (4/3πρ)1/30M m2/3 dm
= G (4/3πρ)1/3 (3/5)M5/3 = 22/332/35–1 (πρ)1/3G M5/3 = 62/3/5 (πρ)1/3G M5/3
where M is the total mass of the pile.
This agrees almost with this.


Expressing the answer in M and the radius r using M = (4π/3)r3ρ
ρ = M(3/(4π))r–3
E = 62/3/5 (πρ)1/3G M5/3
= 62/3/5 (πM(3/(4π))r–3)1/3G M5/3
= 22/3–2/332/3+1/3/5π0M6/3G r–1
= 3/5M2G/r
= 1.6667 M2G/r
= 1.1124 10–10 M2/r m3 kg–1 s–2

Accretion Energy of Earth:
Mass = 5.98 × 1024 kg.
Radius of Earth = 6.37× 106 m.
E = 6.2449 × 1032 Joules


P.S., Today, 2007 July 24, I fixed a bug that crept in as I expressed E in terms of M and r. These calculations are trivial yet very error prone for me. My old equation had the wrong units. I need a calculator that does units!! (So does NASA!)
(If you plan to disassemble the Earth you should first check my calculations.)