Temperature and heat

A preliminary caution

Before we even get started on our modeling project it is worthwhile considering the meaning of the terms, temperature (T) and heat (Q), to see if we really understand what we are modeling. One of the problems of the science of thermodynamics that it is difficult to give a satisfactory definition of temperature without implying the concept of heat, and that it is equally difficult to give a satisfactory definition of heat without bringing in the concept of temperature. Our course is not a course in thermodynamics, and it is not appropriate to launch ourselves into the subtleties of thermodynamic argument. It will be sufficient to note that it is usual to start with the idea of thermal equilibrium. Two bodies are said to be at the same temperature if no heat (energy) flows (or they remain in thermal equilibrium) when they are brought into contact with each other.  Temperatures, or at least scales of temperature, are then defined in terms of certain measurable physical properties such as the volumetric expansion of mercury. It is shown in textbooks of thermodynamics that one particular scale of temperature, the perfect-gas scale, has special significance because it is identical with an absolute thermodynamic scale of temperature based on the second law of thermodynamics.  (Unfortunately there is no perfect gas, and if you delve deeply into the temperature literature, you discover that temperature is only defined with respect to physical standards, and that therefore they are arbitrary)

Heat is said to pass from one system to another if the two systems are initially at different temperatures and are then placed in contact with each other (note that this is actually a circular definition with the above paragraph). The quantity of heat passing to a system is measured in terms of the product of mass (m), specific heat (cp) and temperature rise. Note that heat is only one form of energy in a body, therefore it is difficult to define an absolute heat quantity in a body. Usually we can really only talk about heat transfer and changes in heat quantity.  (Again, if you delve into the heat literature, you find that the quantity heat is just as poorly defined as temperature.)

The above comments are, in practice, largely irrelevant, except in experiments where you are worried about temperatures to 3 decimal places.  However, the concept that the underlying equations or physics of a process are poorly know or understood at some level should be a warning to us as modelers. At this early stage in our thermal modeling, we really don't fully understand what we are modeling. For thermal problems, this is not much of a problem (we understand well enough for the questions we are asking), but for modeling in general this is a serious warning, since we often do not fully understand the processes. No amount of modeling will improve our understanding of the basic processes (although it will hopefully highlight errors to a perceptive modeler).

Heat Transport

The movement of heat energy from one location to another (when there is no radiation or movement of material) is given by the heat flux (q, energy per area per time). The heat flux is the flow of heat energy into a unit volume of material and therefore by heat conservation,

dQ/dt = -dq/dx           (A)

(in one dimension, per unit volume, also to be strictly correct the derivatives need to partials [we will discuss this in class]). Note this is just a statement of conservation, i.e. the amount of heat (Q) at a point varies in time solely on difference between the heat flux (q) coming in one side and going out the other (dq/dx). Heat flux is controlled by another important thermal quantity, the thermal conductivity K. Fourier's Law of Heat Conduction establishes the linear proportionality between heat flow qi and temperature gradient in a thermally isotropic medium

qi= -K dT/dxi   (B)

where the minus sign denotes that heat flows down the temperature gradient. As heat flow is a flux of energy through a unit area per unit time, K has dimensions [mlt-3T-1], and is measured in W m-1 K-1 in the SI system. Its value for silicates is in the range 1-5 W m-1 K-1. (The subscripts in the above equation indicate that the Law holds for all directions simultaneously. For anisotropic media the thermal conductivity is a tensor Kij, which is symmetric because heat flow from an isolated point occurs along straightlines. Isotropy may be safely assumed for polycrystals.)

The change of heat in a region causes a temperature change

m cp D T = D Q          (C)      

which actually is used to define the specific heat (cp). In materials with large thermal expansions there is a significant difference between the temperature change if the system is kept at constant pressure or constant volume during heating. Fortunately, this difference is small in earth materials and we can often get away with using a simple value of specific heat (in geo-problems usually cp, constant pressure, is appropriate).

The conductive properties of a material are also conveniently expressed by its thermal diffusivity k, defined as

k = K / r cp

where the dimensions of k are [l2t-1], and its value is of the order of 10-6 m2s-1 in silicates.

A Thermal Model

At this point we have a basic knowledge that heat flows down temperature gradients, and that the temperature of a region rises proportionally to the quantity of heat added. This gives us sufficient information to do back of the envelope calculations and have a general view of thermal conduction problems. To proceed with a model (for whatever our ultimate purpose [to be discussed in class]) we need to chose a basic approach. Thermal problems are simple enough that it is typical to approach them in a direct mathematical simulation. So we proceed to develop a mathematical description on which we can base our model.

If the heat flux is steady in both space and time, then Fourier's law of heat conduction can be combined with the concept of heat conservation (or continuity of heat assuming no internal heat production), to allow direct calculations of steady state temperature gradients or heat flows. For example, the heat flux near the surface of the earth is generally close to one heat flow unit (hfu), which is defined as 10-6 calories per square centimetre per second (41.84 mW m-2). This is a very small quantity compared to the solar heat flux, so it has no effect on the climate of the earth. It is insufficient to prevent the development of permafrost near the poles, though it does produce a measurable warming of the water in the lowest meter or two of the ocean basins. However, it is simple to directly calculate the approximate background temperature gradient that this flux implies in the near surface crust (order 30 K/km).

Geoscientists are often more interested in the unsteady case. For example, how will the temperature vary with both time and distance from the intrusive contact when a thick sill of hot magma is emplaced in a flat-lying set of strata? This problem is an example of a one-dimensional heat flow problem. In its simplest form, suppose that one side of a uniform layer of material (of density r ) is suddenly raised in temperature by an increment D T. Heat starts to diffuse into the region of lower temperature, but at first the temperature gradient is much larger close to the interface than further away from it. So long as a variable temperature gradient persists, more heat flows into a volume element from one side than leaves from the other: this heat is used to increase the temperature of the element, and the amount used per unit volume depends on the density and specific heat cp. (cp is defined as the heat necessary to raise the temperature of one kilogram by one kelvin). For the moment, we need not be concerned whether this is the specific heat at constant volume, cv, or at constant pressure, cp).

By considering the conservation of energy (heat) in equation A, we can see that a space (or time) varying temperature gradient in equation B leads to a time varying variation in the local heat quantity (Q) and therefore a time variation in the local temperature (T) by equation C. Taking equation A and substituting B (actually it's space derivative) into the right hand side and C (actually it's per unit time per unit volume form) into the left hand side gives

r cp T/ t = dQ/dt = -dq/dx = K 2T/ x2

(again in one dimension, x, and using partial derivative symbols in B and C). This is usually expressed as

T/ t = k 2T/ x2    (D)

where k is the thermal diffusivity. This equation, called the one-dimensional heat flow equation states that the rate at which the temperature at a point is changing with time is proportional to the rate at which the temperature gradient at that point is changing in the direction of heat flow. So, if the heat flow is steady (i.e., if the left-hand side is equal to zero) then the temperature gradient must be a constant: the temperature must vary linearly in the direction of heat flow.

The internal temperature gradient will be constant, if the 1D problem is steady state in time, unless there are internal heat sources.  Internal heat sources or sinks are easily included by adding a source term to the right-hand side

T/ t = k 2T/ x2 + q

Where q is the heat source or sink in terms of Watts per cubic meter.

If the temperature gradient is not initially linear (as in the example of the intrusion of a sill), then as heat diffuses into the overlying layer of thickness d, the temperature gradient will eventually become almost linear and the heat flow almost uniform from one side of the layer to the other (assuming the magma doesn't cool and the surface stays at a constant temperature).

It is very useful when dealing with problems to be able to get an idea of the space and time scales involved. An example would involve answering the question: how long do we need to model, seconds, years or millions of years. In thermal problems this question of scaling can often be approached with the thermal relaxation time. This is a characteristic time for heat to diffuse through a layer and for the layer to approach a linear temperature gradient. It may be obtained by dividing d2 by k:

Tr = d2/k

For a more complete discussion of heat flow problems of geological interest see Turcotte and Schubert (1982, Chapter 4).

More Than One Dimension

So far we have mainly talked about problems with one spatial dimension. For a homogeneous and isotropic material it is relatively straightforward to write down the multi-dimensional form of equation D (although there are several logical and theoretical difficulties, most of which stem from the heterogeneity of any real system). The multi-dimensional form of D just takes the sum of the spatial temperature curvatures in orthogonal directions,

T/ t = k Ñ 2T           (E)

(which I will explain in class), where the Ñ 2 operator stands for taking the second spatial derivative in each direction and summing the result.

Towards a Model

Equation E represents the distillation of many of our ideas about heat flow, and will form the basis of our thermal model. Remember, however, that there are many things that we have already thrown out in constructing E (and we haven't even started modeling). We demand isotropy and homogeneity. We disallow convection, radiation or any heat sources or sinks. We haven't talked about it, but in fact we have disallowed phase changes or crystal habit changes. We have even disallowed dimensional changes, expansions, and an almost infinite list of other possible variations (many of which we may not even know about!).

Well, at least equation E represents an approximate mathematical description of the process that we want to model. We will talk in class about how to go about turning E into a model.

A final note, or problem

Equation E is often used as the starting point of thermal modeling.  In sloppy work, equation E is used even when k varies in space, and modeling is carried out by varying k in different regions of the model, using equation E.  Note however that you actually get a different equation than E if you decide at the outset to let k be a variable.  Can you derive the correct equation for variable k?