|
Houstoun in [1] in 1929 estimates the age of the Earth using heat conduction and Earth temperature. Let's see how close the estimate comes to current day estimate of
 years [2].
The general formula for linear flow in a semi-infinite solid is
![$\displaystyle v = \frac{c}{2\sqrt{\pi \kappa t}} \left[ \int_{0}^{\infty} f(\xi... ... - \int_{0}^{\infty} f(\xi) e^{-\frac{(\xi + x)^2}{4 \kappa t}} \,d \xi \right]$](https://images.physicslibrary.org/cache/objects/936/l2h/img2.png "$\displaystyle v = \frac{c}{2\sqrt{\pi \kappa t}} \left[ \int_{0}^{\infty} f(\xi... ... - \int_{0}^{\infty} f(\xi) e^{-\frac{(\xi + x)^2}{4 \kappa t}} \,d \xi \right]$") |
(1) |
For this problem let  be constant and equal to  ,
 .
Then
![$\displaystyle v = \frac{c}{2\sqrt{\pi \kappa t}} \left[ \int_{0}^{\infty} e^{-\... ...\,d \xi - \int_{0}^{\infty} e^{-\frac{(\xi + x)^2}{4 \kappa t}} \,d \xi \right]$](https://images.physicslibrary.org/cache/objects/936/l2h/img6.png "$\displaystyle v = \frac{c}{2\sqrt{\pi \kappa t}} \left[ \int_{0}^{\infty} e^{-\... ...\,d \xi - \int_{0}^{\infty} e^{-\frac{(\xi + x)^2}{4 \kappa t}} \,d \xi \right]$") |
(2) |
sub in
 |
(3) |
since
 is an even function of  . Write
 ; then the result takes the simpler form
 |
(4) |
Tables of values of this integral have been drawn up, and hence  can be determined as a function of the upper limit.
If we descend into the earth we find that after we pass the points where the diurnal and annual variation cease to be appreciable, the temperature begins to increase. The rate of increase varies from place to place, but may be taken roughly as
 for every 50 feet of descent for depths up to about 1 mile. This increase of temperature is easily explained on the assumption that the center of the earth is at a high temperature and that heat is flowing outwards.
If we assume that the earth was originally at a uniform temperature and that its surface has been always at a constant temperature zero, we can use the above result to find how long it has taken to cool. We neglect the convexity of the earth's surface.
We find from (3) that
Kelvin found by the method indicated in paragraph 81 (need to reference new entry) that  for the material of the earth's surface has the value 400, the units of length and time being the foot and year. Assume that the earth was initially at the temperature of molten rock, i.e. about
 . Insert the value for the gradient at the surface, namely,
 for every 50 feet. Then, writing  in the exponential, we obtain
i.e.
If we write  miles,
At that depth the gradient is only
 of its surface value after  years. We see, therefore, from the first result, that according to our assumptions years have elapsed since the earth was at a temperature of
, and we see from the second result that it is permissible to neglect the convexity of the earth's surface. The assumption throughout all the temperature change of constant conductivity, specific heat and density is, of course, open to questions. Also, the earth may have taken much longer to cool, owing to the liberation of heat due to the radioactive disintegration of some of its material.
So vs
 is not such a bad estimate.
This article is a derivative work of the public domain in [1].
- [1] Houstoun, R.A., "An Introduction to Mathematical Physics" Longmans, Green and CO., London, 1929.
[2] Wikipedia contributors, "Age of Earth," Wikipedia, The Free Encyclopedia (accessed February 20, 2025).
|
