temperature, heat transfer, heat transfer coefficient, flow conditions, laminar flow, Reynolds number, energy transfer, heat flow, heat transfer coefficients, temperature difference, radiative heat transfer, Heat exchangers, radiation heat exchange, heat exchange, parallel flow, temperature differences, Turbulent Flow, Heat conduction, definition of Joule, temperature scales, temperature T2, Thermodynamics, Btu
Content:
1 BASIC HEAT TRANSFER AND SOME APPLICATIONS IN POLYMER PROCESSING (A version of this was published as a book chapter in Plastics Technician's Toolbox, Volume 2, Pages 2133, SPE 2002) John Vlachopoulos and David Strutt www.polydynamics.com Heat transfer is a branch of engineering science which seeks to determine the rate of energy transfer between bodies as a result of temperature differences. The concept of rate is the basic difference between heat transfer and thermodynamics. Thermodynamics deals with systems in equilibrium and is concerned with the amount of heat required to change a system from one state to another. Thermodynamics does not answer the question of "how fast" a change is accomplished. It is the science of heat transfer that deals with this question.
basic concepts The terms temperature and heat are understood intuitively. Temperature is the
physical property that describes whether a body is "hot or cold". If we touch a hot metal at 120°F (49°C) heat flows from the metal to our hand. If we touch a cold block of ice heat will flow from our hand to the ice. So, the intuitive concept of temperature is really defined by the
heat exchange from one body to another. Heat is a form of energy that flows from one body to another as a result of a temperature difference. The two temperature scales used for measurement purposes are the Fahrenheit (°F) and Celsius (°C) scales. These scales were established by specifying the number of increments between the freezing (32°F, 0°C) and boiling point (212°F, 100°C) of water at standard atmospheric pressure. The absolute Celsius scale is called the Kelvin (K) and the absolute Fahrenheit scale is termed Rankine (R). The following conversion relations apply F = 32 + 9 C 5 C = 5 ( F  32) 9 R= 9 K 5 R = F + 459.69 K = C + 273.16 In other words, 0 degrees R = 459.69°F and 0 degrees K = 273.16°C (absolute zero). Temperature is measured by observing its effect on some easily observable property of a measuring device e.g. expansion of mercury in a glass thermometer.
2 Heat, or energy in general, is usually measured in Btu, cal, kcal and Joule (J). 1 Btu will raise 1 lbm of water 1°F at 68°F, 1 cal will raise 1 g of water 1°C at 20°C, 1 kcal will raise 1 kg of water 1°C at 20°C. The definition of Joule (J) comes from the definition of work done by a force on an object: (Newton meter = N m = J) Work W = F S where F is the force and S the distance travelled. Power is the work done per second and J/s is called Watt (W) Power = F S = F V t where t is time and V the velocity of travel. If we multiply Watts by time, the result is work done. The work done by electricity is usually expressed as kilowatt hours, 1000 Watts times 3600 seconds, or 3.6 106 Joules. The following conversion relations apply: 1 Btu = 1055 J 1 kcal = 4182 J 1 Btu = 252 cal 1 W =1 J = 3.413 Btu/hr s The human metabolism requires about 2500 kcal per day (roughly 400 Btu per hour) according to frequently quoted guidelines for daily food consumption. FUNDAMENTAL THERMODYNAMIC LAWS The First Law of Thermodynamics is the principle of conservation of energy, which states that energy cannot be created nor destroyed. A more technical definition is: The increase of internal energy (E) of a given system is equal to heat (Q) absorbed from the surroundings plus the mechanical work (W) added. Q + W = E If no heat enters or leaves the system the process is referred to as adiabatic thus E = W
3 On the other hand, if no work is done by the surroundings, then E = Q Heat and work are different types of energy. Heat capacity Cp (also called Specific Heat) is the amount of heat required to raise the temperature of a body by one degree. Here are some typical values: Water (at 20°C, 68°F) = 1 kcal/kg°C = 4182 J/kg°C Air (approximately) = 0.239 kcal/kg°C = 1000 J/kg°C Polyethylene (approximately) = 0.550 kcal/°C = 2300 J/kg°C Steel (approximately) = 0.108 kcal/°C = 450 J/kg°C Brick (approximately) = 0.215 kcal/°C = 900 J/kg°C For a melting of a solid heat must be added to shake and demolish the
Crystal Structure present. Heat of fusion (Hf) is the amount of heat required to melt a crystalline solid without raising its temperature. It is equal in magnitude (but opposite in sign) to the heat of crystallization. Here are some typical values: Ice Hf = 333,000 J/kg = (80 kcal/kg) This means that to melt 1 kg of water we need the same amount of energy as that required to raise 1 kg of water by 80°C, and yet with the actual melting, there is no increase in temperature. HDPE Hf = 250,000 J/kg LDPE = Hf = 200,000 J/kg For amorphous polymers like PS, PMMA and PC Hf = 0, since they do not have any crystal structure. Heat of reaction is the amount of heat involved in a
chemical reaction (added or removed). All polymerization reactions (production of polymers) are exothermic, i.e. they involve liberation of heat. The Second Law of Thermodynamics relates to the direction in which energy transfer or conversion may take place. The result is the increase of Entropy, which is a loss in the availability of energy for external purposes. Heat will flow from a high temperature to a low temperature. It is impossible to construct a machine or device which will operate continuously by receiving heat from just a single reservoir and producing work. A heat engine is a device that produces net positive work as a result of heat transfer from a hightemperature body to a lowtemperature body. The thermal efficiency of such engines is defined as the ratio of output (energy sought) to the input (energy that costs). The input or the energy that costs is the heat from the high temperature source (indirectly, the cost for the fuel) We define thermal efficiency as
4 Thermal efficiency = W(energy sought) Q (energy that costs) The efficiency of heat engines is higher if the heat source has higher temperature. However, all such actual devices have low efficiencies (e.g. combustion engines for cars, no more than 35%). In other words, only about 1/3 of the energy in gasoline goes to useful work (motion of the car). The rest is wasted due to thermodynamic implications. This is the reason why there is so much research on fuel cells nowadays, in which fuels react with oxygen to produce electricity, at higher efficiencies (5060%). For more information on the thermodynamic laws and their implications the reader is referred to more specialized textbooks [13]. Example 1 Thirty people gather for a cocktail party in a basement room that can be assumed completely sealed off and insulated. The room dimensions are 24 ft 28 ft with 8 ft ceiling. Calculate the temperature rise in 30 minutes. Solution We apply the first law of thermodynamics Q+W=E Since there is no work being added to the air in the room E = Q We will assume that each person gives off approximately 2500 kcal/day (equal to the average metabolic
energy consumption for light activity, like talking and walking around the room) E = Q = 30 2500 = 75,000 kcal/day = 3125 kcal/hour = 52 kcal/min This amount of energy goes to the air inside the room which is roughly Air volume = 24 ft 28 ft 8 ft = 5376 ft3 = 152.25 m3 We neglect the volume occupied by the people. The density of air is about 1.24 kg/m3. So the total mass of air in the room is m = 1.24 152.25 = 188.79 kg. The internal energy change will be
5
E = m Cp T
where Cp is the heat capacity of air (0.239 kcal/°C).
T= E =
52 kcal/min
=1.15C / min = 2.07 F / min
m Cp 188.79 kg 0.239 kcal/kg C
So, in 30 minutes the temperature would rise by 30 1.15°C/min = 34.5°C (62.1°F)! This means that even a cold room would become quickly very hot, if the assumption of complete insulation is valid. In reality, there would be considerable
heat losses to the surroundings that will slow down the temperature rise.
Example 2
In an injection molding machine 20 kg of LDPE are molded per hour. The melt temperature entering the mold is 180°C and the mold temperature is maintained at 40°C by a cooling
water system. Determine the amount of water required to cool the plastic and keep the mold at 40°C, if the difference in inputoutput temperatures of the water is not to exceed 5°C.
Solution
This is a straightforward application of the first law of thermodynamics, that is the principle of conservation of energy. The heat for coming off the solidifying plastic in the mold must be taken away by the water.
Qplastic = Qwater The heat removed from the plastic is equal to the heat given off as the plastic temperature drops from 180°C to 40°C plus the heat of solidification which is equal to the heat of fusion but opposite in sign. As the plastic solidifies and the molecules stop moving randomly, heat is liberated. Qplastic= m cp T + m Hf m is amount of material molded per hour, Cp its heat capacity and Hf the heat of solidification (200,000 J/kg for LDPE which solidifies around 106°C).
kg hr
J
Qplastic = 20
hr
2300
(180  40C)
3600 s
kg C
+ 20 kg hr 200,000 J/kg hr 3600s = 1789 + 1111 = 2900 J s The heat taken up by the water undergoing a 5°C temperature change is
6
Qwater
=
m
cp
T
=
m
4182
J kg
C
5
C
=
m
20910
Therefore,
m 20910 J = 2900 J
kg
s
m = 0.138 kg/s = 499 kg/hr
For another example of cooling of a plastic in a mold the reader is referred to the last section of this chapter: Special Heat Transfer Problems in Plastics Processing.
HEAT TRANSFER MODES Temperature differences cause the flow of heat from a high temperature to a low temperature. There are three modes of heat transfer: conduction, convection, and radiation. The basic microscopic mechanism of conduction is the motion of molecules and electrons. It can occur in solids, liquids and gases. In nonmetallic solids the transfer of heat energy is due mainly to lattice vibrations. In metallic solids we have both lattice vibrations and random motions of free electrons. Consequently metals are more conductive than nonmetals. In gases, we have mainly random motions of molecules. In liquids we have partly random molecular motions and some sort of vibration of the liquid lattice structure. Convection is associated with the transport of a mass of liquid or gas. It can be forced i.e. when assisted by a pump or fan, or free (also called natural convection) when the motion of a fluid occurs due to density differences. If there is an electrical heating element at the corner of a room and air is blown onto the element by a fan, this is forced convection. In the absence of a fan the air surrounding the heating element will get hotter, its density will decrease and the air will move upwards causing natural circulation within the room, as the hot air is replenished by colder air, which gets hot and rises again. Radiation involves electromagnetic waves which are emitted by a body as a result of its temperature. The electromagnetic radiation has a broad spectrum from radio waves to xrays. Between the two extremes a narrow portion of the radiation spectrum is the visible light and a broader one covers the thermal radiation. The earth is heated by sun's radiation.
CONDUCTION Consider the solid wall shown in Fig. 1. Temperature T1 is higher than temperature T2. Heat flows from the high temperature to the low temperature. If A is the area normal to the direction of heat flow, Fourier's Law states that the amount of heat flow is proportional to the area A, the temperature difference T1 T2 and inversely proportional to the thickness of the wall.
7
Q
=
k
A b
( T1

T 2)
The proportionality constant k is the thermal conductivity measured in W/m °C or Btu/hr ft °F
(see Table 1).
A plastic injected into a mold cavity is cooled by heat conduction through the mold wall. In fact, one of the factors considered in choosing mold materials is their thermal conductivity. Aluminum has roughly 5 times the conductivity of steel and Beryllium
copper alloys have about 3 times higher conductivity than steel. Higher conductivity means faster heat removal (or addition).
It is interesting to note the similarity between Fourier's law of heat conduction and Ohm's law of electricity. We may write:
T1

T2
=
b kA
Q
Fourier's Law
E = R I Ohm's Law
where E is the voltage (corresponds to T1 T2), R the resistance (corresponds to b/kA and I the
electric current (corresponds to Q). The quantity b/kA is sometimes referred to as thermal resistance.
Now, let us consider a wall composed of three different materials with thermal conductivities k1, k2 and k3 as shown in Fig. 2. We will have
8
T1

T2
=
b1 k2 A
Q
T2

T3
=
b2 k2 A
Q
T3

T4
=
b3 k3 A
Q
b1 b2 b3
Figure 2 Heat conduction through a composite wall
Summing up the above expressions, we have
T1

T4
=
b1 k1 A
+
b2 k2 A
+
b3 k3 A
Q
This equation is similar to the expression of three electrical resistors in series
E = (R1 + R2 + R3) I We can then write a general expression to calculate the rate of heat flow through a composite wall for more than three layers as follows:
9
Q=
T1  Tn
b1 + b2 + ... bn
k1 A k2 A
kn A
Table 1
Some typical values of thermal conductivity (k)
W/m°C
Btu/hr ft°F
Copper
380
220
Aluminum
204
118
Carbon Steel
43
25
Glass
0.78
0.45
Polymer
0.2
0.115
Water
0.6
0.347
Air
0.025
0.0144
Table 2
Typical values of convection heat transfer coefficients
Mode
W/m2 °C
Btu/hr ft2 °F
AIR, Free Convection
4  28
0.7  5
AIR, Forced Convection
4  570
0.7  100
WATER, Free Convection
284  1500
50  265
WATER, Forced Convection
284  17,000
50  3,000
WATER, Boiling
2840  57,000
500  10,000
STEAM, Condensing
5680  113,000
1,000  20,000
10 For the case of composite cylinder (as shown in Fig. 3) the above equation takes the form
Q=
T1  T4
2
1 k1
L
n
r2 r1
+
2
1 k2
L
n
r3 r2
+
2
1 k3
L
r4 r3
Figure 3 Heat conduction through a composite cylinder. Unsteady heat conduction involves temperature variations with time. For example, if the surface of a body is suddenly raised to higher temperature how long will it take for the temperature to penetrate inside? Such problems require the solution of
differential equations. However, many solutions of such problems for common geometrical shapes appear in the form of charts in specialized texts [46]. With the help of Figs. 4 and 5, it is possible to determine how long it will take for the temperature, at the midplane of plate or axis of a cylinder, to reach a certain value, if the surface temperature is suddenly raised or lowered. On the vertical axes of Figs. 4 and 5, Tm is the (unknown) temperature at the midplane or axis, Ti is the initial temperature of the material (uniform throughout) and T0 the suddenly imposed surface temperature. On the horizontal axes, the quantity kt Fo = Cp x m 2 is the dimensionless Fourier number, and t represents time and x the distance from the surface to the center.
11 Figure 4 Plot for calculating the temperature Tm at the midplane of a plate as a function of time after the two surfaces are suddenly raised to T0 Figure 5 Plot for calculating the temperature Tm at the axis of a cylinder as a function of time after the surface tempearature is suddenly raised to T0
12 CONVECTION In most heat transfer problems, we are concerned with solid walls separating liquids or gases from each other. In such cases we usually do not know the temperatures on the wall surfaces, but rather the temperatures of the bulk of fluids on both sides. Careful experiments supported also by theoretical considerations, have shown that the greatest temperature drop is confined within a thin fluid layer attached to a solid surface, as shown in Fig. 6.
Figure 6 Heat transfer through a wall separating two fluids.
To explain this observation, we may assume that a thin film, of thickness , adheres to the wall, whereas outside this film all temperature differences vanish as a result of mixing motions. Within the film heat flow takes place by conduction, as in a solid wall. Thus, in general, we may write
Q
=
k
A(T

Tw)
The quantity k/ = h is called the heat transfer coefficient and it is an extremely important concept in heat transfer. This simplified model is very useful for
practical applications, because the calculation of heat transfer can be made in terms of the heat transfer coefficient:
Q = h A (T  Tw)
At this point it suffices to say that the heat transfer coefficient depends on the flow conditions and fluid properties. Typical values are given in Table 2. In the next section we will present some correlations that can be used for the more accurate determination of this coefficient.
For the wall separating the two fluids of Fig. 6, we have
13
Q = h1 A (T1  Tw1)
Q= k b
A
(Tw1  Tw2)
Q = h2 A (Tw2  T2)
or
By summing up
T1

T w1
=
1 h1 A
Q
Tw1 
Tw2
=
b kA
Q
Tw2

T2
=
1 h2 A
Q
T1

T2
=
1 h1 A
+
b kA
+
1 h2 A
Q
or
T1

T2
=
1 A
1 h1
+
b k
+
1 h2
Q
The quantity in the brackets is called the total thermal resistance. A more useful concept, however, is the overall heat transfer coefficient, which is defined as follows
1 = 1 +b+ 1 U h1 k h2 For a composite wall separating two fluids (a and b), we have
Then, in general, we have:
1 = 1 + b1 + b2 + ... + 1
U ha k1 k2
hb
11 T1  T2 = A U Q and Q=U A T
14 The significance of the overall heat transfer coefficient is that it permits the calculation of the rate of heat flow by multiplying this quantity by the heat exchange area (perpendicular to the heat flow direction) and the temperature difference.
Example
Determine the heat loss through an 8ft by 4ft glass window of 4 mm thickness. The inside temperature is assumed to be 24°C (75°F) and the outside temperature is 10°C (14°F). The inside heat transfer coefficient is 5 W/m°C and the outside about 20 W/m°C (due to moderate wind). The thermal conductivity of window glass is 0.78 W/m°C.
Solution
The overall heat transfer coefficient is
1 U
=
1+ h1
b k
+
1 h2
=
5
1 W / m2 C
+
0.004 m 0.78 W / m2 C
+
20
1 W / m2 C
= 0.2 + 0.005 + 0.05 = 0.255 m2 C /W
U = 3.92 W /m2°C The rate of heat flow is Q = U A (T  T2) = 3.92 W / m2 C (2.44 1.22 m2) (24C  (10C)) = 3.92 2.98 34 = 397 W = 1355 Btu / hr
Convection Heat Transfer Coefficient Calculation The most important step in heat convection calculations is the determination of the appropriate heat transfer coefficient. The higher the fluid velocity is, the higher the heat transfer coefficient will be. Numerous correlations have been developed for the calculation of the heat transfer coefficient in terms of dimensionless groups: Nusselt number is the dimensionless heat transfer coefficient defined as Nu = h D k
15 where the heat transfer coefficient h has dimensions of W/m2°C or Btu/hr m2°F, k the thermal conductivity has dimensions of W/m°C or Btu/hr ft°F and D is a characteristic length (m or ft), such as the diameter of a pipe in which a fluid is flowing.
The Reynolds number (see also the chapter on
fluid mechanics), is defined as
Re = V D
where is the density (kg/m3 or lbm/ft3), V the fluid velocity (m/s or ft/s), the viscosity (Pa s or lbfs/ft2) and D a characteristic length (m or ft) such as the diameter of a pipe, in which a fluid is flowing.
When the value of the Reynolds number for pipe flow is less than 2100, the flow is streamlined and regular and is called laminar. Above Re = 2100, the flow is highly chaotic and irregular and is said to be turbulent.
The Prandtl number is defined as
Pr = Cp k where Cp is the heat capacity (J/kg°C or Btu/lbm°F), the viscosity (Pa s or lbfs/ft2) and k the thermal conductivity (W/m°C or Btu/hr ft°F).
For Laminar Flow in tubes, the following correlation applies:
Nu
=
3.66
+
1
0.0668 (D + 0.04 [(D /
/ L) Re Pr L) Re Pr ]2
/
3
For Turbulent Flow in tubes, we have Nu = 0.023 Re0.8 Prn
where n = 0.4 for heating and n = 0.3 for cooling.
Numerous other correlations are available in specialized textbooks [46]. Generally, for laminar flow Nu depends on the 0.3 power of Re, while for turbulent flow on the 0.8 power of Re.
To calculate the heat transfer coefficient we plug the various quantities in the appropriate correlation to calculate the Nusselt Number Nu. Then, from the definition
we can calculate
Nu = h D k
16 h = Nu k D The values given in Table 2 have been calculate from such correlations assuming flow conditions encountered in various practical heat transfer situations (according to reference [4]).
Let us now see how much more effective is turbulent heat transfer over laminar. Assuming D/L = 1/30, Re = 2000 and Pr = 1 (i.e. water at 200°C), the correlation for laminar flow gives
0.0668 1 2000 1
Nu = 3.66 +
30
1 + 0.04
1 30
2000 1 2 / 3
= 6.37
and h = 6.37 k D We repeat the calculation for turbulent flow assuming Re = 3000 the other quantities remaining the same.
Nu = 0.023 (3000 )0.8 10.4 = 13.9 and h =13.9 k D We see that while the Reynolds number (i.e. velocity) was increased by 50% the heat transfer coefficient was increased by 118% i.e. (13.9/6.39) 100 = 218% 100% = 118%. Obviously, to maximize heat transfer rates, we should operate with turbulent flow conditions. For example, by increasing the channel diameter and the coolant velocity, the Reynolds number in the cooling channels of a mold can be increased, beyond the critical value of 2100. The resulting turbulent flow will be much more effective in removing heat from the mold.
Molten polymers are very viscous and the Reynolds number is in the range 0.0001 to 0.01, i.e. polymer melt flows are always laminar. So, for a polymer melt having k = 0.2 W/m°C, Pr 5 106, and flowing through a 8 mm diameter channel the heat transfer coefficient will be roughly
h = 7 0.2 =175 W/ m2 C 0.008
17 RADIATION Many types of thermoforming machines use radiation to heat the plastic sheet. Metal rod heaters, halogen tubes and ceramic plates are used frequently. Heater temperatures can reach 700°C and the corresponding radiative heat transfer coefficients can reach 100 W/m2°C. Radiation is important wherever very
high temperatures are involved. However, radiative effects play a significant role is other less obvious situations. Goose down is known as the most effective insulation for two reasons: reduction of heat conduction by air pockets trapped by clusters of fibers and entrapment of radiation. The goose down fibers have diameter of few microns which is the wavelength of a significant portion of the infrared radiation that escapes from the human body and can go through air pockets and fabrics. Radiation scattering occurs due to equivalence of fiber size and wavelength. Synthetic materials have been manufactured that reportedly match the goose down properties both in reducing conduction and radiation heat loss, again by scattering on micronsized fibers. The radiant heat transfer (which is a form of electromagnetic radiation), depends on the absolute temperature and the nature of the radiating surface. This is stated by the StefanBoltzmann law for a black body (perfect radiator) Q = A T4 where T = absolute temperature of the surface of the body in K or R, A the
surface area (m2 of ft2) and the StefanBoltzmann constant ( = 5.669 08 W/m2K4 or 0.1714 108 Btu/hr ft2 R4). Black bodies are called so, because such materials do not reflect any radiation and appear black to the eye. Thus, a black body absorbs all radiation incident upon it. The "blackness" of a surface to thermal radiation can be quite deceiving insofar as visual observations are concerned [4]. Some visually black surfaces are indeed black to thermal radiation. However, snow and ice appear white and bright, but are essentially "black" in thermal radiation. When a body with surface temperature T1 is placed in a closed environment of temperature T2, the net amount of heat transfer depends on the temperature difference, in the form Q = A1 (T14  T42) For gray bodies, not perfect radiators, we have Q = 1 A1 (T14  T42) where 1 = emissivity (1 = 1 for black bodies, 1 < 1 for gray bodies). For lampblack = 0.96, for oxidized cast iron = 0.7, and for polished steel = 0.1) If surface A1 is not completely enclosed by surface A2 we must introduce an additional factor to account for the relative geometrical orientation of the two radiating surfaces
18
Q = A1 F12 (T14  T42) F12 is called the shape or view factor or angle factor. During the winter in the
Northern Hemisphere, the temperatures are low because the sun's radiation arrives with an unfavorable angle factor.
For practical problems, it is sometimes advisable to define a radiative heat transfer coefficient hr from the following
Q = hr A1 (T1  T2) = A1 F12 (T14  T42)
hr =
F12 (T14  T42) T1  T2
=
F12
(T12 + T22)
(T1 + T2)
The radiative heat transfer coefficient depends strongly on temperature and is less useful as a concept than the convective heat transfer coefficient. However, it is useful for practical problems involving both convection and radiation, for which we can write the rate of heat flow as
Q = (h + hr) Aw (Tw  T)
where Tw is the temperature of a wall and T the temperature of a surrounding medium which completely encloses surface Aw. Due to the fourth power dependence of radiation heat exchange, this mode is more important in very high temperature applications, as in metallurgical operations. In some polymer processing operations, radiation is less important than convection. For example, in blown film extrusion, the film is likely to emerge from the die lips at 200°C (473°K) or so, and a tangentially impinging cooling air jet is likely to have heat transfer coefficient of the order of 100 W/m2°C. Let us calculate the radiative heat transfer coefficient assuming F12 = 1, and T2 = 20°C (293°K)
hr = F12 (T12 + T22) (T1 + T2) = 5.669 108 W / m2 K4 (4732 + 2932) (473 + 293) = 13.44 W / m2 C This is significantly lower than the convective heat transfer coefficient (100 W/m2°C) and may be neglected in some calculations.
heat exchangers Heat exchangers are devices that transfer heat by convection and conduction between two fluids which are separated by a wall. The automobile radiator is a heat exchanger in which convection
19 and conduction take place but no radiative heat transfer. Radiative heat exchangers are used in specialized applications in space vehicles and are beyond the scope of most textbooks. The simplest type of heat exchanger is the double pipe system shown in Fig. 7. It involves tube and annular flow. If the flows are in the same direction, the arrangement is called parallel flow and if the flows are in opposite directions, counterflow. Because of the rather small
Surface areas which are available for heat transfer, double pipe heat exchangers are used for low to moderate heat transfer rates. For high rates of heat, other types of heat exchangers are used that provide large surface areas. In this category belong the shellandtube heat exchangers (that involve several tube passes) and crossflow configurations (with interconnected passageways).
Figure 7 Doublepipe heat exchanger (parallel flow)
For the design and prediction of performance purposes, the concept of the overall heat transfer coefficient is used, which was described earlier in this chapter.
For the double pipe arrangement of Fig. 7, we can write again
Q = U A Tm where U is the overall heat transfer coefficient, A the heat transfer area, and Tm a suitable
average temperature difference between the entering and exiting fluids. It turns out that this average is the socalled logmean temperature difference (LMTD) which is defined as
Tm
=
n
T2 [ T2
/
T1 T1]
where n represents the natural logarithm having base e = 2.71828, while log has base 10.
20 The above relationships are valid for both parallel and counterflow heat exchangers. Injection molding machines involving hot oil and cold water are connected counterflow since this arrangement removes roughly 10% more heat from the oil than parallel flow.
For other types of heat exchangers (multipass shellandtube; cross flow, etc.), a correction factor is needed
Tm = (LMTDfor counter flow) Fcorrectionfactor
F is usually between 0.5 and 1.0 and values of this factor for various configurations of heat exchangers can be found in specialized textbooks [46].
Heat exchangers may become coated with various deposits present in the flow or corroded. Because of fouling a reduction in heat transfer efficiency is observed. The overall effect is represented by the fouling factors
Definition
1
1
Rf =

Ufouled Uclean
tabulated values of Rf can be found in specialized textbooks [46]. For boilers used for long periods of time we may have Rf = 0.0002 m2°C/W, while Uclean = 3000 W/m2°C. Thus 0.0002 = 1  1 Ufouled 3000 Ufouled = 1875 W/m2°C which represents a significant reduction
It is important to keep clean channels and clean circulating water or oil in plastics processing machinery, otherwise fouling will result is significant cooling efficiency reduction.
SPECIAL HEAT TRANSFER PROBLEMS IN PLASTICS PROCESSING There are two properties of plastics that play a very significant role during their processing: their low thermal conductivity and their high viscosity in the molten state. The thermal conductivity of most plastics is around 0.2 W/m°C (0.115 Btu/hr ft°F) which is roughly 200 times smaller than the conductivity of steel and 2000 times smaller than the conductivity of copper. Due to low conductivities, cooling of plastics is slow. In fact, in some processes, such as thick profile extrusion, pipe extrusion and film blowing, cooling might be the output rate limiting step. While the mextruder might be capable of pumping out more product, we might not be able to cool the product fast enough. The viscosity of molten plastics decreases with the rate of shearing, but it always remains very high, roughly of the order of a million times larger than the viscosity of water. Due to high
21 viscosity, polymer melts flowing through channels, dies and process equipment tend to raise their temperature by a mechanism of internal friction. Under certain conditions, the frictional heating might result in high temperatures which can cause degradation of the material.
Cooling of a Plastic Plate in a Mold Let us assume that a HDPE plate of 10 mm thickness initially at 230°C is cooled in a mold and the mold surface is at 30°C. How long will it take for the temperature to reach 90°C at the midplane (xm = 10/2 = 5 mm). We will use Fig. 4 for Ti = 230°C, T0 = 30°C, and Tm = 90°C. The
vertical axis of Fig. 4 is calculated as follows
Tm  T0 = 90  30 = 60 = 0.3 Ti  T0 230  30 200 The corresponding Fourier Number is
kt
Fo =
Cp
x
2 m
= 0.58
Using typical thermal properties for HDPE (see Table 3), we have:
0.25 t
= 0.58
780 2300 0.0052
t =104 seconds
For computer assisted mold cooling calculations and optimization procedures the reader is referred to a thesis [7].
Temperature Rise Due to Frictional Heating The average temperature rise due to frictional heating (also called viscous dissipation) can be calculated assuming that the process is adiabatic (i.e. there is no heat exchange with the surroundings). This is actually not such a bad assumption for flows through channels and dies. Assuming that the mechanical work (due to the pressure pushing the molten plastic) is converted into heat, we can easily prove that the temperature rise will be T= P Cp
22 So, for a die that has a pressure drop P = 4000 psi (27.58 106 Pa) and a molten polymer having = 780 kg/m3 and Cp = 2300 J/kg°C, we get T = 27.58 106 =15.4C (27.7 F) 780 2300 This means that the average temperature of the polymer coming out of the die will be 15.4°C (27.7°F) higher than the temperature coming in.
Power Requirement of a SingleScrew Extruder
This example and the selected thermal properties given in Table 3 are taken from reference [8]. In the extruder, polymer pellets usually coming in at
room temperature (i.e. 20°C) are melted and subsequently pumped at the extrusion temperature through the die (i.e. ~200°C, depending on polymer). Most of the energy comes from turning the screw. Of course some energy is supplied by the heating bands around the barrel. In well running extruders usually net energy input occurs in the first section (near the hopper ) and net output in the second section (near the die) i.e. the heat generated by the viscous dissipation is really heating the barrel.
The power required by the turning screw is needed to: o Raise the temperature from room temperature to extrusion temperature in the die. o Melt the polymer (heat of fusion). o Pump the molten polymer.
The power required by the turning screw is given by the following expression :
Po = Q Cp (Tout  Tin) + Q hf + P Q
Where:
Symbol m m Cp Tin Tout Q P Hf
Name Mass Throughput (kg/hr) Heat Capacity, Average (J/kg °C) Usually room temperature °C Extrusion Temperature in the die °C Density, average (kg/m3) Volume
flow rate (m3/hr) Pressure rise (Pa) Heat of fusion (J/kg)
Typical Values 15003000 J/kg °C 20°C 200°C 300°C 7001200 kg/m3 1050 MPa The heat of fusion can vary from zero (for amorphous polymers) to 300,000 J/kg (for a very crystalline HDPE).
Let us calculate the relative contributions for Q = m =124.22 kg / hr of a typi
Cal Polymer and pressure rise of 30 MPa:
23
Po =112.4 2500 (200  20) +112.4 130,000 + 30 106 112.4 760 Po = 50,580 kJ/hr + 14,612 kJ/hr + 4,437 kJ/hr = 69,629 kJ/hr
Term Raise temperature to 200C Melt the polymer Pump the molten polymer Total Power Requirement
Relative Significance Most Important Somewhat Important Insignificant
Contribution (kJ/hr) 50,580 14,612 4,437 69,629
Let us express the final answer in terms of horsepower:
Po = 69,629 kJ hr =19.34 kW Po =19.34 kW hp = 25.93 hp
hr 3600 s
0.746 kW
So we can size the horsepower of a motor. But we must take into account its efficiency. If we assume 85% efficiency, then
MotorPower = 1 25.93 hp = 30.50 hp = 40.88 kW 0.85
SUMMARY Heat is a form of energy that is transferred from high temperatures to low temperatures by three modes: CONDUCTION, CONVECTION and RADIATION. Conduction through metals e.g. copper, aluminum and beryllium copper alloys is very fast, while through plastics is very slow. The thermal conductivity of aluminum is about 200 W/m°C and of typical plastic about 0.20 W/m°C. This means that heat is conducted one thousand times faster through aluminum than through a plastic. In solid walls the temperature drops linearly from its high value on one side to the low value on the other side. Convection can be natural (free) or forced. Natural convection involves the motion of a gas or liquid due to density differences i.e. hot air is lighter and as it moves upward new air moves near the heating source causing circulation. In forced convection, the movement of a liquid or gas is assisted by a pump or fan. For convective heat transfer calculations it is important to determine the heat transfer coefficient. The heat transfer coefficient increases as the flow velocity increases. Turbulent flow is much more effective than laminar flow for transferring of heat. Radiation is the transfer of heat by electromagnetic waves. It is important whenever high temperatures are involved, such as in plastic sheet heating in thermoforming machines. In plastics processing, heat is generated within flowing highly viscous melts due to internal friction (viscous dissipation). Most of the energy necessary to heat and melt plastic pellets in an extruder comes from the motor turning the screw. This means we have conversion of mechanical energy into heat.
24 It is possible to carry out relatively accurate calculations of temperature rise, drop or distribution and the energy requirements, by applying the principle of conservation of energy and the equations of heat conduction, convection and radiation. For conduction the thermal conductivity of the materials involved must be known. For convection the key step is the determination of the heat transfer coefficients, which are usually available as correlations involving the Reynolds number of the flow. Radiation calculations require the material emissivities and the relative geometrical orientation of the radiating surfaces. REFERENCES 1. Holman, J.P., "Thermodynamics" McGrawHill, New York (1988). 2. Sonntag, R.E. and C. Borgnakke, "Introduction to Engineering Thermodynamics", Wiley (2000). 3. Smith, J.M. and H.C. Van Ness, "Introduction to
Chemical Engineering Thermodynamics", McGrawHill, New York (2001). 4. Holman, J.P., "Heat Transfer", McGrawHill, New York (1990). 5. Bejan, A., "Heat Transfer", Wiley, New York (1993). 6. Pitts, D.R. and L.E. Sisson, "Heat Transfer", (Schaum's Outline Series), McGrawHill, New York (1998). 7. Reynolds, D.W., "Optimization of Cooling Circuits in Injection Molds", M.Eng. Thesis,
mechanical engineering,
McMaster University, Hamilton, ON, Canada (2000). 8. Vlachopoulos, J. and J.R. Wagner, "The SPE Guide on Extrusion Technology and Troubleshooting",
Society of Plastics Engineers, Brookfield, CT (2001). 9. Mark, J.E., "Physical Properties of Polymers Handbook", AIP Press, Woodbury, NY (1996).
25
Table 3 Selected Thermal Properties.
Polymer
Solid Density* (g/cm3)
Glass Transition Tg
Melting Point Tm
Usual Melt Processing Range
Melt Density* (kg/m3)
Thermal Conductivity k (W/m C) (Btu/h ft F)
Heat Capacity Cp (J/kg C) (Btu/lbm F)
Heat of Fusion H (J/kg) (Btu/lb)
HDPE 0.9410.967 130C
130137C 160240C
780
202F
266278F 320464F
0.25 0.145
22002400 210,000300,000
0.520.57
90130
LDPE 0.9150.935 130C
106112C 160240C
760
202F
223234F 320464F
0.20 0.115
22002400 190,000240,000
0.520.57
80100
LLDPE 0.9100.925
130C 202F
125C 257F
160240C
760
320464F
0.20 0.115
22002400 190,000240,000
0.520.57
80100
PP
0.8900.910 20C
4F
165C 329F
180240C
730
356464F
0.18
20002200 210,000260,000
0.10
0.480.52
90110
PVC (Rigid)
1.301.58
80C 176F
175C 347F
165205C 329401F
1250
0.17
10001700 170,000190,000
0.10
0.240.41
7080
PS
1.041.10
100C amorphous** 180240C
1000
212F
356464F
0.15
13002000 amorphous**
0.09
0.310.48
PMMA 1.171.20
105C 221F
amorphous** 180230C 356446F
1050
0.19
14002400 amorphous**
0.11
0.330.57
PET
1.341.39
80C
176F
265C 509F
275290C 527554F
1160
0.18
18002000 120,000140,000
0.10
0.430.48
5060
ABS
1.011.04 105115C amorphous** 200290C
990
221239F
392554F
0.25 0.145
13001700 0.310.41
amorphous**
Nylon66 1.131.15
90C 194F
265C 509F
275290C
980
527554F
0.20 0.115
24002600 190,000205,000
0.570.62
8088
PC
1.2
140C amorphous** 250305C
1050
0.22
13002200 amorphous**
284F
482581F
0.13
0.310.52
* Melt densities have been estimated for roughly the midtemperature of the processing range. See Mark [9] for expressions in the form = ABT CT2.
** Amorphous resin does not possess crystallinity and consequently no melting point or heat of fusion (i.e. heat to break down crystal
structure) can be determined.
26 John Vlachopoulos is Professor of Chemical Engineering and Director of the Centre for Advanced Polymer Processing and Design (CAPPAD) at McMaster University, Hamilton,
Ontario, Canada. He earned his Dipl. Ing. from NTU,
Athens, Greece and an M.S. and D.Sc. from
Washington University,
St. Louis, Mo. He is the author of approximately 200 publications on various aspects of plastics processing, frequently lectures at conferences and seminars, and serves as consultant to industry worldwide. He established and is President of Polydynamics, Inc., a company involved in software research, development and marketing. John was elected an SPE Fellow and received the SPE Education Award in 2001.