Thermal conduction

Thermal conduction is the transfer of internal energy by microscopic diffusion and collisions of particles or quasi-particles within a body. The microscopically diffusing and colliding objects include molecules, atoms, and electrons. They transfer disorganized microscopic kinetic and potential energy, which are jointly known as internal energy. Conduction can only take place within an object or material, or between two objects that are in contact with each other. Conduction takes place in all phases of ponderable matter, such as solids, liquids, gases and plasmas, but it is distinctly recognizable only when the matter is undergoing neither chemical reaction nor differential local internal flows of distinct chemical constituents. In the presence of such chemically defined contributory sub-processes, only the flow of internal energy is recognizable, as distinct from thermal conduction. When the processes of conduction yield a net flow of energy across a boundary because of a temperature gradient, the process is characterized as a flow of heat.

Heat spontaneously flows from a hotter to a colder body. In the absence of external drivers, temperature differences decay over time, and the bodies approach thermal equilibrium.

In conduction, the heat flow is within and through the body itself. In contrast, in heat transfer by thermal radiation, the transfer is often between bodies, which may be separated spatially. Also possible is transfer of heat by a combination of conduction and thermal radiation. In convection, internal energy is carried between bodies by a material carrier. In solids, conduction is mediated by the combination of vibrations and collisions of molecules, of propagation and collisions of phonons, and of diffusion and collisions of free electrons. In gases and liquids, conduction is due to the collisions and diffusion of molecules during their random motion. Photons in this context do not collide with one another, and so heat transport by electromagnetic radiation is conceptually distinct from heat conduction by microscopic diffusion and collisions of material particles and phonons. In condensed matter, such as a solid or liquid, the distinction between conduction and radiative transfer of heat is clear in physical concept, but it is often not phenomenologically clear, unless the material is semi-transparent.

In the engineering sciences, heat transfer includes the processes of thermal radiation, convection, and sometimes mass transfer. Usually more than one of these processes occurs in a given situation. The conventional symbol for the material property, thermal conductivity, is $k$ .

Overview

Steady-state conduction

Steady state conduction is the form of conduction that happens when the temperature difference(s) driving the conduction are constant, so that (after an equilibration time), the spatial distribution of temperatures (temperature field) in the conducting object does not change any further. Thus, all partial derivatives of temperature with respect to space may either be zero or have nonzero values, but all derivatives of temperature at any point with respect to time are uniformly zero. In steady state conduction, the amount of heat entering any region of an object is equal to amount of heat coming out (if this were not so, the temperature would be rising or falling, as thermal energy was tapped or trapped in a region).

For example, a bar may be cold at one end and hot at the other, but after a state of steady state conduction is reached, the spatial gradient of temperatures along the bar does not change any further, as time proceeds. Instead, the temperature at any given section of the rod remains constant, and this temperature varies linearly in space, along the direction of heat transfer.

In steady state conduction, all the laws of direct current electrical conduction can be applied to "heat currents". In such cases, it is possible to take "thermal resistances" as the analog to electrical resistances. In such cases, temperature plays the role of voltage, and heat transferred per unit time (heat power) is the analog of electrical current. Steady state systems can be modelled by networks of such thermal resistances in series and in parallel, in exact analogy to electrical networks of resistors. See purely resistive thermal circuits for an example of such a network.

Transient conduction

Main article: Heat equation

In general, during any period in which temperatures change in time at any place within an object, the mode of thermal energy flow is termed transient conduction. Another term is "non steady-state" conduction, referring to time-dependence of temperature fields in an object. Non-steady-state situations appear after an imposed change in temperature at a boundary of an object. They may also occur with temperature changes inside an object, as a result of a new source or sink of heat suddenly introduced within an object, causing temperatures near the source or sink to change in time.

When a new perturbation of temperature of this type happens, temperatures within the system change in time toward a new equilibrium with the new conditions, provided that these do not change. After equilibrium, heat flow into the system once again equals the heat flow out, and temperatures at each point inside the system no longer change. Once this happens, transient conduction is ended, although steady-state conduction may continue if heat flow continues.

If changes in external temperatures or internal heat generation changes are too rapid for equilibrium of temperatures in space to take place, then the system never reaches a state of unchanging temperature distribution in time, and the system remains in a transient state.

An example of a new source of heat "turning on" within an object, causing transient conduction, is an engine starting in an automobile. In this case the transient thermal conduction phase for the entire machine is over, and the steady state phase appears, as soon as the engine reaches steady-state operating temperature. In this state of steady-state equilibrium, temperatures vary greatly from the engine cylinders to other parts of the automobile, but at no point in space within the automobile does temperature increase or decrease. After establishing this state, the transient conduction phase of heat transfer is over.

New external conditions also cause this process: for example the copper bar in the example steady-state conduction experiences transient conduction as soon as one end is subjected to a different temperature from the other. Over time, the field of temperatures inside the bar reach a new steady-state, in which a constant temperature gradient along the bar is finally set up, and this gradient then stays constant in space. Typically, such a new steady state gradient is approached exponentially with time after a new temperature-or-heat source or sink, has been introduced. When a "transient conduction" phase is over, heat flow may still continue at high power, so long as temperatures do not change.

An example of transient conduction that does not end with steady-state conduction, but rather no conduction, occurs when a hot copper ball is dropped into oil at a low temperature. Here the temperature field within the object begins to change as a function of time, as the heat is removed from the metal, and the interest lies in analyzing this spatial change of temperature within the object over time, until all gradients disappear entirely (the ball has reached the same temperature as the oil). Mathematically, this condition is also approached exponentially; in theory it takes infinite time, but in practice it is over, for all intents and purposes, in a much shorter period. At the end of this process with no heat sink but the internal parts of the ball (which are finite), there is no steady state heat conduction to reach. Such a state never occurs in this situation, but rather the end of the process is when there is no heat conduction at all.

Analysis of non steady-state conduction systems is more complex than steady-state systems. If the conducting body has a simple shape then exact analytical mathematical expressions and solutions may be possible (see heat equation for the analytical approach).^[1] However, most often, because of complicated shapes with varying thermal conductivities within the shape (i.e., most complex objects, mechanisms or machines in engineering) often the application of approximate theories is required, and/or numerical analysis by computer. One popular graphical method involves the use of Heisler Charts.

Occasionally transient conduction problems may be considerably simplified if regions of the object being heated or cooled can be identified, for which thermal conductivity is very much greater than that for heat paths leading into the region. In this case, the region with high conductivity can often be treated in the lumped capacitance model, as a "lump" of material with a simple thermal capacitance consisting of its aggregate heat capacity. Such regions warm or cool, but in the process, show no significant temperature variation across their extent, during the process (as compared to the rest of the system). This is due to their far higher conductance. During transient conduction, therefore, the temperature across their conductive regions changes uniformly in space, and as a simple exponential in time. An example of such systems are those that follow Newton's law of cooling during transient cooling (or the reverse during heating). The equivalent thermal circuit consists of a simple capacitor in series with a resistor. In such cases, the remainder of the system with high thermal resistance (comparatively low conductivity) plays the role of the resistor in the circuit.

Relativistic conduction

The theory of relativistic heat conduction is a model that is compatible with the theory of special relativity. For most of the last century, it was recognized that Fourier equation is in contradiction with the theory of relativity because it admits an infinite speed of propagation of heat signals. For example, according to Fourier equation, a pulse of heat at the origin would be felt at infinity instantaneously. The speed of information propagation is faster than the speed of light in vacuum, which is physically inadmissible within the framework of relativity. Alterations to the Fourier model provided for a relativistic model of heat conduction, avoiding this problem.

Quantum conduction

Second sound is a quantum mechanical phenomenon in which heat transfer occurs by wave-like motion, rather than by the more usual mechanism of diffusion. Heat takes the place of pressure in normal sound waves. This leads to a very high thermal conductivity. It is known as "second sound" because the wave motion of heat is similar to the propagation of sound in air.

Fourier's law

The law of heat conduction, also known as Fourier's law, states that the time rate of heat transfer through a material is proportional to the negative gradient in the temperature and to the area, at right angles to that gradient, through which the heat flows. We can state this law in two equivalent forms: the integral form, in which we look at the amount of energy flowing into or out of a body as a whole, and the differential form, in which we look at the flow rates or fluxes of energy locally.

Newton's law of cooling is a discrete analog of Fourier's law, while Ohm's law is the electrical analogue of Fourier's law.

Differential form

The differential form of Fourier's Law of thermal conduction shows that the local heat flux density, $\overrightarrow{q}$ , is equal to the product of thermal conductivity, $k$ , and the negative local temperature gradient, $-\nabla T$ . The heat flux density is the amount of energy that flows through a unit area per unit time.

\overrightarrow{q} = - k {\nabla} T

where (including the SI units)

\overrightarrow{q}

is the local heat flux density, W·m⁻²

\big.k\big.

is the material's conductivity, W·m⁻¹·K⁻¹,

\big.\nabla T\big.

is the temperature gradient, K·m⁻¹.

The thermal conductivity, $k$ , is often treated as a constant, though this is not always true. While the thermal conductivity of a material generally varies with temperature, the variation can be small over a significant range of temperatures for some common materials. In anisotropic materials, the thermal conductivity typically varies with orientation; in this case $k$ is represented by a second-order tensor. In non uniform materials, $k$ varies with spatial location.

For many simple applications, Fourier's law is used in its one-dimensional form. In the x-direction,

q_x = - k \frac{d T}{d x}

Integral form

By integrating the differential form over the material's total surface $S$ , we arrive at the integral form of Fourier's law:

\frac{\partial Q}{\partial t} = -k \oint_S{{\nabla} T \cdot \,\overrightarrow{dA}}

where (including the SI units):

$\big. \frac{\partial Q}{\partial t}\big.$ is the amount of heat transferred per unit time (in W), and
$\overrightarrow{dA}$ is an oriented surface area element (in m²)

The above differential equation, when integrated for a homogeneous material of 1-D geometry between two endpoints at constant temperature, gives the heat flow rate as:

\big. \frac{\Delta Q}{\Delta t} = -k A \frac{\Delta T}{\Delta x}

where

A is the cross-sectional surface area,

\Delta T

is the temperature difference between the ends,

\Delta x

is the distance between the ends.

This law forms the basis for the derivation of the heat equation.

Conductance

Writing

\big. U = \frac{k}{\Delta x}, \quad

where U is the conductance, in W/(m² K).

Fourier's law can also be stated as:

\big. \frac{\Delta Q}{\Delta t} = U A\, (-\Delta T).

The reciprocal of conductance is resistance, R, given by:

\big. R = \frac{1}{U} = \frac{\Delta x}{k} = \frac{A\, (-\Delta T)}{\frac{\Delta Q}{\Delta t}}.

Resistance is additive when several conducting layers lie between the hot and cool regions, because A and Q are the same for all layers. In a multilayer partition, the total conductance is related to the conductance of its layers by:

\big. \frac{1}{U} = \frac{1}{U_1} + \frac{1}{U_2} + \frac{1}{U_3}+ \cdots

So, when dealing with a multilayer partition, the following formula is usually used:

\big. \frac{\Delta Q}{\Delta t} = \frac{A\,(-\Delta T)}{\frac{\Delta x_1}{k_1} + \frac{\Delta x_2}{k_2} + \frac{\Delta x_3}{k_3}+ \cdots}.

For heat conduction from one fluid to another through a barrier, it is sometimes important to consider the conductance of the thin film of fluid that remains stationary next to the barrier. This thin film of fluid is difficult to quantify because its characteristics depend upon complex conditions of turbulence and viscosity—but when dealing with thin high-conductance barriers it can sometimes be quite significant.

Intensive-property representation

The previous conductance equations, written in terms of extensive properties, can be reformulated in terms of intensive properties.

Ideally, the formulae for conductance should produce a quantity with dimensions independent of distance, like Ohm's Law for electrical resistance: $R = V/I\,\!$ , and conductance: $G = I/V \,\!$ .

From the electrical formula: $R = \rho x / A \,\!$ , where ρ is resistivity, x is length, and A is cross-sectional area, we have $G = k A / x \,\!$ , where G is conductance, k is conductivity, x is length, and A is cross-sectional area.

For Heat,

\big. U = \frac{k A} {\Delta x}, \quad

where U is the conductance.

Fourier's law can also be stated as:

\big. \dot{Q} = U \, \Delta T \quad

analogous to Ohm's law: $I = V/R \,\!$ or $I = V G. \,\!$

The reciprocal of conductance is resistance, R, given by:

\big. R = \frac{\, \Delta T}{\dot{Q}}, \quad

analogous to Ohm's law: $R = V/I. \,\!$

The rules for combining resistances and conductances (in series and in parallel) are the same for both heat flow and electric current.

Cylindrical shells

Conduction through cylindrical shells (e.g. pipes) can be calculated from the internal radius, $r_1$ , the external radius, $r_2$ , the length, $\ell$ , and the temperature difference between the inner and outer wall, $T_2 - T_1$ .

The surface area of the cylinder is $A_r= 2 \pi r \ell$

When Fourier's equation is applied:

\dot{Q} = -k A_r \frac{\mathrm{d}T}{\mathrm{d}r} = -2 k \pi r \ell \frac{\mathrm{d}T}{\mathrm{d}r}

and rearranged:

\dot{Q} \int_{r_1}^{r_2} \frac{1}{r} \mathrm{d}r = -2 k \pi \ell \int_{T_1}^{T_2} \mathrm{d}T

then the rate of heat transfer is:

\dot{Q} = 2 k \pi \ell \frac{T_1 - T_2}{\ln (r_2 /r_1)}

the thermal resistance is:

R_c = \frac{\Delta T}{\dot{Q}}= \frac{\ln (r_2 /r_1)}{2 \pi k \ell}

and $\dot{Q} = 2 \pi k \ell r_m \frac{T_1-T_2}{r_2-r_1}$ , where $r_m = \frac{r_2-r_1}{\ln (r_2 /r_1)}$ . It is important to note that this is the log-mean radius.

Spherical

The conduction through a spherical shell with internal radius, $r_1$ , and external radius, $r_2$ , can be calculated in a similar manner as for a cylindrical shell.

The surface area of the sphere is: $\!A = 4\pi r^2.$

Solving in a similar manner as for a cylindrical shell (see above) produces: $\dot{Q} = 4 k \pi \frac{T_1 - T_2}{1/{r_1}-1/{r_2}} = 4 k \pi \frac{(T_1 - T_2) r_1 r_2}{r_2-r_1}$

Transient thermal conduction

Main article: heat equation

Interface heat transfer

The heat transfer at an interface is considered a transient heat flow. To analyze this problem, the Biot number is important to understand how the system behaves. The Biot number is determined by: $\textit{Bi}= \frac{hL}{k}$ The heat transfer coefficient, h, is introduced in this formula, and is measured in $\frac{J}{m^{2}sK}$ . If the system has a Biot number of less than 0.1, the material behaves according to Newtonian cooling, i.e. with negligible temperature gradient within the body. If the Biot number is greater than 0.1, the system behaves as a series solution. The temperature profile in terms of time can be determined by the function can be derived from the equation:

$q=-h \Delta T$

Which becomes:

$\frac{T-T_{f}}{T_{i} - T_{f}} = \operatorname{exp} \left [ \frac{-hAt}{\rho C_p V} \right ]$

The heat transfer coefficient, h, is measured in $\mathrm{\frac{W}{m^{2}K}}$ , and represents the transfer of heat at an interface between two materials. This value is different at every interface, and is an important concept in understanding heat flow at an interface.

The Series Solution can be analyzed with a nomogram. A nomogram has relative temperature as the y coordinate and the Fourier number, which is calculated by: $\textit{Fo}= \frac{\alpha t}{L^2}$ The Biot number increases as the Fourier number decreases. There are 5 steps to determine a temperature profile in terms of time.

Calculate the Biot number
Determine which relative depth matters, either x or L.
Convert time to the Fourier number.
Convert $T_i$ to relative temperature with the boundary conditions.
Compared required point to trace specified Biot number on the nomogram.

Thermal conduction applications

Splat cooling

Splat cooling is a method for quenching small droplets of molten materials by rapid contact with a cold surface. The particles undergo a characteristic cooling process, with the heat profile at $t=0$ for initial temperature as the maximum at $x=0$ and $T = 0$ at $x = -\infin$ and $x = \infin$ , and the heat profile at $t=\infin$ for $-\infin \le x \le \infin$ as the boundary conditions. Splat cooling rapidly ends in a steady state temperature, and is similar in form to the Gaussian diffusion equation. The temperature profile, with respect to the position and time of this type of cooling, varies with:

$T(x,t) - T_i = \frac{T_i \Delta X}{2\sqrt{\pi \alpha t}} \operatorname{exp} \left ( -\frac{x^2}{4 \alpha t} \right )$

Splat cooling is a fundamental concept that has been adapted for practical use in the form of thermal spraying. The thermal diffusivity coefficient, represented as $\alpha$ , can be written as $\alpha =\frac{k}{\rho C_p}$ . This varies according to the material.^[2]^[3]

Metal quenching

Metal quenching is a transient heat transfer process in terms of the Time Temperature Transformation (TTT). It is possible to manipulate the cooling process to adjust the phase of a suitable material. For example, appropriate quenching of steel can convert a desirable proportion of its content of austenite to martensite, creating a very tough product. To achieve this it is necessary to quench at the "nose" (or eutectic) of the TTT diagram. Since materials differ in their Biot numbers, the time it takes for the material to quench, or the Fourier number, varies in practice.^[4] In steel, the quenching temperature range is generally from 600 °C to 200 °C. To control the quenching time and to select suitable quenching media, it is necessary to determine the Fourier number from the desired quenching time, the relative temperature drop, and the relevant Biot number. Usually the correct figures are read from a standard nomogram. By calculating the heat transfer coefficient from this Biot number, one can find a liquid medium suitable for the application.^[5]

Zeroth law of thermodynamics

One statement of the so-called zeroth law of thermodynamics is directly focused on the idea of conduction of heat. Bailyn (1994) writes that "... the zeroth law may be stated:

All diathermal walls are equivalent."^[6]

A diathermal wall is a physical connection between two bodies that allows the passage of heat between them. Bailyn is referring to diathermal walls that exclusively connect two bodies, especially conductive walls.

This statement of the 'zeroth law' belongs to an idealized theoretical discourse, and actual physical walls may have peculiarities that do not conform to its generality.

For example, the material of the wall must not undergo a phase transition, such as evaporation or fusion, at the temperature at which it must conduct heat. But when only thermal equilibrium is considered and time is not urgent, so that the conductivity of the material does not matter too much, one suitable heat conductor is as good as another. Conversely, another aspect of the zeroth law is that, subject again to suitable restrictions, a given diathermal wall is indifferent to the nature of the heat bath to which it is connected. For example the glass bulb of a thermometer acts as a diathermal wall whether exposed to a gas or to a liquid, provided they do not corrode or melt it.

These indifferences are amongst the defining characteristics of heat transfer. In a sense they are symmetries of heat transfer.

References

↑ The Exact Analytical Conduction Toolbox contains a variety of transient expressions for heat conduction, along with algorithms and computer code for obtaining precise numerical values.
↑ Sam Zhang; Dongliang Zhao (19 November 2012). Aeronautical and Aerospace Materials Handbook. CRC Press. pp. 304–. ISBN 978-1-4398-7329-8. Retrieved 7 May 2013.
↑ Martin Eein (2002). Drop-Surface Interactions. Springer. pp. 174–. ISBN 978-3-211-83692-7. Retrieved 7 May 2013.
↑ Rajiv Asthana; Ashok Kumar; Narendra B. Dahotre (9 January 2006). Materials Processing and Manufacturing Science. Butterworth-Heinemann. pp. 158–. ISBN 978-0-08-046488-6. Retrieved 7 May 2013.
↑ George E. Totten (2002). Handbook of Residual Stress and Deformation of Steel. ASM International. pp. 322–. ISBN 978-1-61503-227-3. Retrieved 7 May 2013.
↑ Bailyn, M. (1994). A Survey of Thermodynamics, American Institute of Physics, New York, ISBN 0-88318-797-3, page 23.

Dehghani, F 2007, CHNG2801 – Conservation and Transport Processes: Course Notes, University of Sydney, Sydney
John H Lienhard IV and John H Lienhard V, 'A Heat Transfer Textbook', Third Edition, Phlogyston Press, Cambridge Massachusetts

External links

Heat conduction - Thermal-FluidsPedia
Newton's Law of Cooling by Jeff Bryant based on a program by Stephen Wolfram, Wolfram Demonstrations Project.