Conduction
Conduction
Conduction
A T1
T2 T3
qx
T
LA x qx T1 LA kA A T2 LB kBA T3 LC kC A T4 LB LC
2005 The Board of Trustees of the University of Illinois All Rights Reserved
II. Background
A. Fourier's Law r r The fundamental relation between the heat flux, q" , and the temperature gradient, T , for conduction heat transfer in an isotropic medium* is
r r q" = kT
(1)
where k (W/(m.K)) is an important property of the material called the thermal conductivity. This relationship is based on observations, and as such we refer to it as phenomenological; it is named Fouriers Law. Thoughtful inspection of Fouriers Law can tell you a lot about conduction heat transfer. Namely, the conductive heat flux is a vectorit has a magnitude and direction. The magnitude of the conductive heat flux is proportional to the temperature gradient, and its direction is coincident with the temperature gradient. Since the temperature gradient is in the direction of maximum temperature change, the heat flux is in the direction of maximum temperature change; i.e., it is perpendicular to the isotherms. The
*
The words isotropic medium mean a material in which the properties have no directional dependence.
(2)
Cylindrical:
(3)
B. The Heat Diffusion Equation The partial differential equation governing the temperature distribution in a stationary medium is called the heat diffusion equation (HDE). You can derive the HDE by performing an energy balance on an appropriate differential control volume (see ref. [1]). In general vector form, the result may be written as
r r & = c p T q"+q t r or substituting for q" from Fourier's law r r & = c p T k T + q t
(4)
(5)
& is the energy generation rate per unit volume (W/m3) in the medium. Physically, this equation where q reflects an energy balance at any arbitrary point within the medium. It states that the net rate of conduction per unit volume to the point, plus the volumetric energy generation rate at the point, must equal the time rate of change of thermal energy per unit volume at that point.
(6)
Cylindrical:
1 T 1 T T & T k + k kr + 2 + q = c p r r r r z z t
(7)
x=L
Assuming negligible heat loss from the cylinder periphery and that the x-coordinate is aligned with the cylinder axis, the heat transfer and temperature gradient will occur in the x-direction only (i.e., heat transfer is by one-dimensional conduction). For steady state conditions (T/t = 0) with no heat generation, the heat diffusion equation reduces to
dq x " =0 dx
(8)
or d dT k =0 dx dx (9)
The first form states that under these conditions the heat flux, qx , is a constant, independent of x. If the thermal conductivity can be assumed to be constant, the second form reduces to
What is important here is that the area perpendicular to the conductive flux is constant. Although a circular cylinder has been shown, the cross-sectional area could be any shapesquare, triangular, rectangular, etc. We used a cylindrical geometry in the sketch because it looks like the specimens you will use in the laboratory. You may consider the system as Cartesian if the area for conduction is constant in the direction of heat transfer.
(10)
which can be integrated twice to obtain the general form of the temperature distribution
T (x ) = C1 x + C 2
(11)
Applying the boundary condition T(0) = T1at the left end gives
C 2 = T1
(12)
Likewise, applying the boundary condition at the right end T(L) = T2yields
C1 =
(T
T1 ) L
(13)
Substituting into the general solution for T(x): x x T(x ) = (T2 T1 ) + T1 = T1 (T1 T2 ) L L or in dimensionless form
T(x ) T1 x = T2 T1 L
(14)
(15)
This result demonstrates the well-known finding that for steady, one-dimensional heat conduction through a plane wall with constant thermal conductivity and no volumetric heat generation, the temperature varies linearly with x. Now, using Fourier's law to relate the heat flux to the temperature gradient
qx " = qx dT k (T2 T1 ) k (T1 T2 ) = k = = A dx L L
(16)
qx is the heat transfer rate in the x-direction, and A is the cross-sectional area normal to the x-direction. Rearranging for qx
qx = kA (T1 T2 ) L
(17)
Note that by solving this equation for (T1- T2)/L and substituting into the solution for T(x), the temperature distribution may also be expressed in the form
(18)
In analogy with Ohm's law for electrical resistance, Relec= E/I, the thermal resistance is defined as the ratio of the driving potential T to the corresponding transfer rate qx. For steady heat transfer through a plane wall, the thermal resistance for conduction is therefore given by
R cond =
(T
T2 ) L = qx kA
(19)
D. Composite Walls One-dimensional conduction through more complex systems, such as composite walls, is easily analyzed using the thermal resistance concept. For example, consider the three-layer plane wall shown below.
T1
T2 T3
qx
T
LA
x
LB T2 T3 LB kBA
LC T4 LC kC A
qx
T1 LA kA A
From the definition of thermal resistance, the heat transfer rate qx can be expressed in terms of the overall temperature difference T = (T1- T4)as
qx T (T1 T4 ) , = R R
(20)
(21)
or
qx =
(T1 T4 ) LA LC LB k A + k A + k A B C A
(22)
E. Contact Resistance Up to this point, we have neglected the finite resistance that occurs at the interface between materials. This contact resistance is primarily a result of surface roughness, in that regions of contact between two materials are interspersed with gaps, which are usually filled with air.
Gap T(x) TA q
x
B T
TB
Contact Region
The heat transfer that occurs at the interface may be thought of as the parallel contributions of conduction across the contact areas and conduction, convection, and/or radiation across the interstitial gaps. Especially for rough surfaces, the contact area at the interface is small and the major resistance is due to the gaps. In analogy to the definition of the conduction thermal resistance given above, the contact resistance for a unit area of the interface is defined as the ratio of the temperature difference at the interface (driving potential) to the heat flux (transfer rate)
R cont T (TA TB ) = qx qx
(23)
r2 qr r1 T2 T T1 T1
T2
r1
r2
For steady state conditions with no internal heat generation, and assuming a circumferentially symmetric temperature distribution (T/ = 0), with no heat conduction in the axial direction (T/z = 0), the heat transfer will occur only radially through the cylindrical walls. The heat diffusion equation, written in cylindrical coordinates, reduces to d dT kr =0 dr dr (24)
Note that for radial heat conduction, the area perpendicular to the radial heat flux, 2rL, changes with radial location, r. Radial conduction cannot be analyzed without considering the cylindrical geometry.
(25)
where qr is the radial heat transfer rate, q r is the radial heat flux, and A = 2rL is the cross-sectional area normal to the heat flux at any radius r. Multiplying the first equation above by the constant (-2L) and moving this quantity inside the first derivative, we recognize that it states dqr =0 dr
(26)
or, for steady state conditions with no volumetric heat generation, the radial heat transfer rate qr, (but not This result should be obvious from the heat flux q r ) is constant, independent of the radial location. energy balance considerations. Returning to the simplified heat diffusion equation above, and now assuming constant thermal conductivity, it reduces to
d dT r =0 dr dr
(27)
(28)
Applying the boundary conditions T(r1) = T1 and T(r2) = T2, solving simultaneously for the constants C1 and C2, and substituting into the expression above, the particular solution for the radial temperature profile is given by
T (r ) =
(T T ) ln r +T ln(r r ) r
2 1 2 1 1
= T1
(T
) ln r 1 T2 ln(r2 r1 ) r 1
(29)
or in dimensionless form
T(r ) T1 ln(r r1 ) = T2 T1 ln(r2 r1 )
(30)
This result shows that, under these assumptions, the temperature varies logarithmically for radial conduction through cylindrical walls, rather than linearly as for plane walls. Applying Fourier's Law to relate the heat-transfer rate to the temperature difference across the cylindrical walls
(31)
or
qr =
(32)
If this expression is solved for (T1 - T2)/ln(r2/r1) and the result substituted into the equation for T(r), the temperature distribution can also be written as
T (r ) = T1 r qr ln 2Lk r1
(33)
(34)
we see that the conduction resistance for radial heat transfer through cylindrical walls is given by
R cond = ln(r2 r1 ) 2Lk
(35)
III. Experiments
Objectives The objectives of this experiment are as follows: (1) To determine the thermal conductivity of several metallic and insulating materials; (2) To investigate conduction heat transfer in a composite wall system; (3) To evaluate the contact resistance at the interface between materials; (4) To investigate radial conduction through a cylindrical-walled system; (5) To compare the experimental measurements to handbook or theoretical results and to discuss and explain any differences.
Brass (70% Cu, 30% Zn): 304 stainless steel: 316 stainless steel: Teflon:
For the materials listed above, the value of thermal conductivity depends on the specific composition and temperaturegiven above are nominal values. Unfortunately, the equipment manufacturer did not specify the thermal conductivity for the samples supplied with the apparatus you will use in your experiments.
Apparatus Measurements of conduction heat transfer will be obtained using an experimental setup consisting of two components: a "linear module" for the study of steady, one-dimensional conduction through plane walls; and a "radial module" for investigation of steady, one-dimensional radial conduction through cylindrical walls. Linear Module The linear module is sketched below. This apparatus has a control unit, shown on the left-hand side of the sketch, and a thermal unit, shown on the right-hand side of the sketch. Within the thermal unit, the conduction sample is clamped between heater and cooler sections as shown. In the sketch, the heater section is at the left end of the thermal unit; it consists of an electrical resistance heater with a 25mm-diameter brass rod immediately to its right. Three temperature sensors (1, 2, and 3) are installed at 10-mm intervals along the length of this brass rod. The cooler section is a 25-mm-diameter brass rod with chilled water running over its right end (see schematic). Three temperature sensors (7, 8, and 9) are located along this rod at 10-mm intervals. Several samples are available for study as listed below.
10
54 C
60.4
Cooling water "On/Off" switch Thermocouple selector switch Heater power control
Nominal Dimensions 25 mm dia. x 30 mm long 25 mm dia. x 30 mm long 25 mm dia. heat transfer area
The metallic samples are instrumented with three temperature sensors (4, 5, and 6) located at 10-mm intervals; the Teflon sample is not instrumented. The ends of the brass heater and cooler sections and the metallic samples have all been given a fine surface finish to promote good thermal contact; a heat-conducting compound may also be applied at these surfaces to reduce thermal resistance. In order to minimize energy losses from the periphery of the linear module, the heater and cooler sections and the metallic conduction samples are enclosed in an insulating case that includes an air gap. Electrical power is supplied to the heater section by means of the control unit, shown in the left portion of the sketch above. Heater power is controlled over the range from 0 - 100 Watts, by a variable transformer and can be read by means of the wattmeter digital display located on the front panel of the console. A nine-position selector switch permits the temperature to be read, via a digital display, at any one of the sensor locations on the linear module.
Radial Module The radial module is shown in the schematic below. Like the linear module, the radial module has a control unit and a thermal unit. The thermal unit consists of an annular brass disk of nominal thickness 3 mm and inner and outer radii of r1 = 4 mm and r2 = 55 mm, respectively. An electrical heater is located in the center of the disk (r < 4 mm), and the outer periphery of the disk is cooled by means of water circulating in a copper tube. Six temperature sensors are located along the disk radius, starting at the center, r = 0, and then in 10-mm intervals. Again, energy loss is minimized by preserving an air gap around the disk with an insulating case. Electrical power is controlled and supplied to the heater by
11
6 5 4 3 2 1
Disk
54 C
60.4
Experiment 1: Thermal Conductivity Measurements In this experiment the linear module is used to measure the thermal conductivity of several typical metals and insulators: brass, stainless steel, and Teflon. In order to maximize the resolution of the temperature measurements, the electrical power supplied to the heater should be adjusted so that the maximum temperature in the heater section is in the range T1 = 80 to 100C. Under no circumstances should you exceed T1 = 100Cexceeding this temperature will damage the experimental apparatus. A data sheet will be supplied for your use during the laboratory exercise; it will include suggested heater powers. In order to obtain accurate measurements, steady state conditions should be achieved as closely as possible. You should monitor system temperatures to determine when steady state is attained. For the thermal conductivity measurements, the thermal grease (a heat-conductive compound) supplied with the apparatus should be applied at each interface, metal-to-metal or metal-to-Teflon. This compound will minimize contact resistance at these joints. For each conductivity determination, the conduction sample is clamped between the heating and cooling sections, the appropriate heater power is supplied, the temperature readings are monitored until steady state conditions are reached, and these temperature measurements are then recorded. To measure the thermal conductivity of brass, use the 25-mm-diameter x 30-mm long instrumented brass sample. With this specimen, the brass heater, sample, and cooler sections will provide an 80-mm long continuous brass sample with 9 temperature sensors (T1 - T9) located every 10 mm along its length. A least-squares line can be fit to the temperature readings to determine the temperature gradient along the brass rod.
12
0 10 20 30 40 50 60 70 80
x (mm)
The thermal conductivity of stainless steel is determined by measuring the temperature distribution in the 25-mm-diameter x 30-mm-long SS sample at steady-state conditions (sensors T4 - T6). A least-squares linear fit to these three temperature readings can then be used to determine the temperature gradient in the stainless steel. Note the sensors in this sample are spaced at 10-mm intervals with T4 and T6 assumed to be 5 mm from the adjacent ends.
Stainless Steel Sample T4
Heater
Cooler T6
0 10 20 30 40 50 60 70 80
x (mm)
The same general experimental method is used for determining k for the Teflon sample. However, since this sample is not instrumented with temperature sensors, a different procedure must be used for determining the temperature gradient. In this case, the three temperature readings in the heating section (T1 - T3) are least-squares fit and the resulting linear equation is used to extrapolate for the heater/sample interface temperature TA. A similar procedure is used for the three cooling-section temperature readings (T7 - T9) in order to extrapolate for the sample/cooler interface temperature TB Assuming a linear temperature profile through the sample then allows determination of its temperature gradient. Note that we will assume that sensors T3 and T7 are located exactly 5 mm from the adjacent interface. The thickness of the Teflon sample should be measured with a micrometer. The procedure just described could also have been used to determine the thermal conductivity of brass and stainless steel, but direct measurement of the temperature gradient using the instrumented samples, as outlined above, is preferable.
13
Once the temperature gradient in each sample is determined, its thermal conductivity can be found from Fourier's Law
k= q x A dT dx
(36)
The values determined from the measurements should be compared to handbook values for each material.
Experiment 2: Composite Wall Heat Conduction In this experiment, steady heat conduction through a composite plane wall will be studied using the linear module. The composite wall to be investigated consists of the following system:
Teflon sample (25 mm diameter heat transfer area) Brass heater (25 mm diameter) T1 Brass sample (25 mm diameter) Brass cooler (25 mm diameter) T9 Teflon thickness 25 mm R1 30 mm R2 R3 25 mm R4 T9
qx
qx
T1
Conductive paste should be applied at all interfaces before the samples are clamped into the linear module. The procedure is similar to that discussed previously. An appropriate heater power is supplied such that T1 is 80 to 100C, and the system is allowed to come to steady state. In this case, however, only temperatures T1 and T9 are recorded. Using the conductivity values determined in the earlier experiments, estimate the individual resistances R1, R2, R3, R4 indicated in the sketch above and the total resistance R. From the total thermal resistance and the measured temperatures T1 and T9, compute the expected heat transfer rate
14
(T
T9 ) R
(39)
Heater T1 TB
Sample Cooler TA TC T9 TD
0 10 20 30 40 50 60 70 80
x (mm)
Note that the interfaces are assumed to be at x = 25 and 55 mm, respectively. Using the temperature difference and measured heat flux, qx = qx/A, determine and compare the contact resistance for a unit area of each interface based on the definition
R cont T qx "
(40)
15
(T
T6 ) r ln ln(r6 r1 ) r1
1
(41)
T (r ) = T1
r qr ln 2Lk r1
(42)
In the first expression use the measured values of T1 and T6, and in the second, use the measured values of T1 and qr and the value of thermal conductivity for brass determined in your experiments. Note that the disk thickness, L, is nominally 3 mm. Discuss the agreement and/or discrepancies between the data and the two theoretical temperature profiles.
16
[1] Incropera, F. P., and DeWitt, D. P., Fundamentals of Heat and Mass Transfer, 4th edition, Wiley, New York, 1996.
17
Apparatus Number________
Insert
Brass (25mm x 30mm) Stainless Steel (25 x 30mm) Teflon 25mm dia. Teflon thickness
Apparatus Number________
Insert
T4
T5
T6
T7
T8
T9
Apparatus Number________
Notes: 1. One interface should be clean and dry; the other should have conductive compound applied. 2. Set T1 to be 90oC.
Insert Heater Power T1 T2 T3 T4 T5 T6 T7 T8 T9
Apparatus Number________
1. The "inner radius" is 4mm; T1 can be assumed to equal the temperature at r=4mm. 2. Set T1 to be 90oC.
Heater Power T1 T2 T3 T4 T5 T6
Radial Setup