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 =
∫r0∞f 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/3∫0M 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
I need a calculator that does units!!
(So does NASA!)