Well, we want to know the temperature distribution

*T*of the ground as a function of time

*t*and position

*x*. So some googlin' or previous knowledge shows that we need to solve the heat equation. For our purposes, we can treat the Earth as flat (I don't plan on digging a basement deep enough to see the curvature of the Earth), so we can assume the temperature only changes with the depth into the ground

*x*:

\[

\frac {\partial T}{\partial t} = a \frac {\partial^2 T} {\partial x^2} \qquad (1)

\]

where

*a*is the thermal diffusivity of the material, in units of square meters per second. It looks like we're going to have to solve some partial differential equations!

Or will we? We can get a very good estimate of how much the temperature changes with depth just by dimensional analysis. Let's measure our time

*t*in terms of a characteristic time of our problem

*w*(it could be 1 year if we were trying to see the change in the ground's temperature from summer to winter, or 1 day if we were looking at the change from day to night). Then we can write:

\[

\frac {\partial T } {\partial t} = \frac 1 w \frac {\partial T} {\partial t/w}

\]

... plugging this in Eq. (1), rearranging, and calling

*l*= sqrt(

*w*a*) gives....

\[

\frac {\partial T}{\partial (t/w)} = \frac {\partial ^2 T} {\partial (x/ l )^2}

\]

Now let's say we didn't know how to or didn't want to solve this equation. (Don't worry, we do & we will). From rearranging this equation, we see right away there is only one "length scale" in the problem,

*l*. So if we had to guess, we could guess that the ground changes temperature a distance

*l*into the ground. A quick look at Wikipedia for thermal diffusivities gives us the following table, for materials we'd find in the ground:

Material | a, mm^2/s |
l (cm), w = 1 day |
l (cm), w = 1 year |
---|---|---|---|

Polycrystalline Silica (glass, sand) | 0.83 | 27 cm | 5.1 meters |

Crystalline Silica (quartz) | 1.4 | 35 | 6.6 |

Sandstone | 1.15 | 32 | 6.0 |

Brick | 0.52 | 21 | 4.0 |

Soil | 0.3-1.25 | 16-33 | 3.1-6.3 |

Just to pat ourselves on the back for our skills at dimensional analysis, a quick check shows that permafrost penetrates 14.6 feet into the ground after 1 year. So our dimensional estimates looks pretty good!

In the next few days I'll solve this equation exactly and throw up a few pretty graphs, and maybe talk a little about PDE's and Fourier series in the process.

?

ReplyDelete