Heat Transfer Notes Ebook
Heat Transfer Notes Ebook
This temperature difference is thought of as a driving force that causes heat to flow. The concepts of heat transfer and temperature, the key words in the discipline of heat transfer, are 2 of the most basic concepts of thermodynamics.
drivingfffff ffffffffff fffffforcef fffffffff Rate of transport process= fffffffffff or rate = coefficient Bdriving force resistance
System: a region in space containing a quantity of matter which is separated from its surroundings by a boundary. Closed system (no flow system): no exchange of matter with the surroundings, only heat and work cross the boundary. Open system (flow system): there is matter exchange with the surroundings in addition to heat and work. Work (W): is a transient quantity (energy) which only appears at the boundary when a system changes its state due to the movement of a part of the boundary under the action of a force. Sign convention: system + the work is done by the system on the surroundings: the work exits the
(Although there cannot be said to be any work in a system either before or after the change has taken place, work may be said to flow or be transferred across the boundary.) Heat (Q): is something (energy transfer), which only appears at the boundary when a system changes its state due to a difference in temperature between the system and its surrounding. Heat, like work, is a transient quantity, which only appears at a boundary while a change is taking place within the system. (Although there cannot be said to be any heat in a system before and after a change of state, loosely speaking heat may be said to flow or be transferred across the boundary. Strictly speaking it is energy which is transferred, but to say heat is transferred is a shorthand way of saying the energy transferred by virtue of a temperature difference.) Sigh convention: + if heat flows into a system from the surroundings: If heat flows from the system to the surroundings.
1.1 First Law of Thermodynamics It is the principle of conservation of energy. It is an axiom. The first law of thermodynamics says that there exists a property of a closed system (U) such that a change in its value is equal to the difference between the heat supplied and the work done during any change of state:
1
X Q @ W =U 2 @U 1
1
Writing Q and W for the quantities of heat and work crossing the boundary during the change of non heat energy equation state: Q-W=U2-U1 By words: any quantity of heat supplied to a closed system must equal the increase of internal energy plus the work done by the system. The internal energy of a closed system remains unchanged. For isolated systems Q=0, W=0 therefore U = 0 For irreversible non-flow processes the energy equation can only be applied in integrated form:
Q @ W = U
For reversible processes, the energy equation can only be applied in differential form:
V
d Q @ W = d U Q @ pdV = dU W = pdV
The potential energy term is either zero or negligible compared with the other terms. The 1st law of thermodynamics does not make any distinction between heat transfer and work transfer: to it they are both energy interactions (non-properties) that must be distinguished from the energy change (property) 1.2 Second Law of Thermodynamics It is an axiom. It says that it is impossible to construct a system which will operate in a cycle, extract heat from a reservoir, and does an equivalent amount of work on the surroundings. The 1st law says: the net work can never be greater than the heat supplied.
2
The 2nd law says: it must always be less. Therefore if a system is to undergo a cycle and produce work, it must operate between at least 2 reservoirs of different temperature. As a consequence, work is a more valuable form of energy transfer than heat: heat can never be transformed continuously and completely into heat.
dQf fff fff ff f There exists a property of a closed system (s) such that a change in its value is equal to Z T
1
for any reversible process undergone by the system between state 1 and state 2
Z
1 2f
Mathematically:
= S 2 @ S 1 where S = entropy, J / K
f
rev
or in differential for:
= dS
rev
The entropy of a reversible isolated closed system remains constant. (For reversible adiabatic a dQf fff ff f process dQ = 0 but dS = fff d S = 0 [ S = constant ) [ T The entropy of an irreversible isolated closed system increases:
2f 2f
g g
=0
A,irreversible
since dQ = 0
Z
1
= S 2 @ S1
1f
B,irreversible
dQf dQf fff ff f fff ff f The cycle as a wjhole is irreversible, and E fff Z fff = T T
2
<0
B,irreversible
Therefore S1-S2<0
and
S2-S1
The proper distinction between heat transfer and work transfer is made by the 2nd law of thermodynamics. There is always an entropy transfer associated with the heat transfer. A work transfer carries zero entropy. The 2nd law of only provides definite quantitative statement about reversible processes. Only for such processes is it possible to predict the work and heat transfers crossing the boundary of a system. For irreversible processes the law merely provides statements of trend (that the entropy of an isolated irreversible system must increase), and quantitative prediction of energy transfer
3
cannot be made. Therefore empirical relations (such as Newtons law of viscosity, Fouriers law of heat conduction, and Ficks law of diffusion) are used. They all take the form of proportionality between a quantity transferred in an irreversible process and the property gradient. These laws are approximate empirical relations, whose accuracy is manifested by the fact that the proportionality constants are not in fact constants but vary with the conditions of the experiment (temperature, pressure, or their gradient,)
1.3 TEMPERATURE It is the system property that determines whether the system is in thermal equilibrium with another system. In thermal equilibrium of system A and B: TA=TB. The temperature of a system is measured by placing the system in contact with a special system (a test system) called thermometer. Temperature scales have 2 reference points. Celsius scale 0oC Ice-point of water (ice and air-saturated water in equilibrium at standard atmosphere pressure) 100 oC: Boiling point of water (liquid water in equilibrium with its own vapor at standard atmospheric pressure) Fahrenheit scale 32 oF: Ice-point of water 212 oF: Boiling point of water ToC =5/9 [T(oF)-32] Kelvin scale (thermodynamic scale or absolute scale) 273.16K: Tripe-point of water (it is only slightly above that of the ice-point 273.15K) 0 K: Absolute zero The unit of thermodynamic temperature is 1/273.16 = 1K This odd unit makes Tbp-Ticepoint = 100K and T C = K T The temperature difference in Kelvin is the same as it is in oC 0 oC=273.15K and T(K)= T(oC)+273.15
4
` a ` a
1.4 The Modes of Heat Transfer There are 3 basic modes of heat transport: 1. Conduction 2. Convection 3. Radiation Conduction (Fouriers law) Heat is transported on a molecular scale with no movement of macroscopic portions of matter relative to one another. Conduction can take place through solids, liquids, and gases. In fluids: The kinetic energy of the molecules is associated with the property we call temperature. In high temperature regions, molecules have higher kinetic energies than those in a low temperature region. The molecules in the high- temperature region transfer their energy, through collisions, to molecules in the low-temperature region. In solids: Conduction is due to motion of free electrons in metals, lattice waves (vibrations of the crystal lattice structure) in non metals, magnetic excitations and electromagnetic radiation. Convection (sensible heat; and Newtons law) It is heat transport due to the bulk motion of the fluid. The process occurs through the movement of macroscopic volume elements of the fluid in space from a region of one temperature to that of another. The energy transported in this way is called sensible heat.
Q = mc p T
b
fl in
@T
fl out
,J
Convection also involves the energy exchange between a solid and a fluid (interface transport). Convection is only possible in a fluid medium 2 types of convection are distinguished (a) free or natural convection The fluid moves because of the density difference resulting from the temperature difference in the fluid (b) forced convection The fluid is forced to flow past a solid surface
Radiation It is the transport of energy by electromagnetic radiation having a defined range of wavelength. The wavelength range of interest in thermal radiation extends from 0.2-20m. (The wavelength range of visible radiation is 0.3-0.8m, and that of infrared (IR) radiation is 0.8-400m). All substances emit radiant heat but the net flow of heat is from the high to low temperature region. So the cooler substance will absorb more radiant heat energy than it emits. The radiation from all bodies depends on temperature (increases with temperature). No physical medium is needed for the propagation of radiation. Often all 3 forms of radiation are involved simultaneously. It is then usual to calculate the heat transport by each mode separately and adding the separate effects to provide an estimate of the total rate of heat transport. In a number of cases, 1 mode of heat transport is dominant. 1.5 Fouriers law of heat conduction Fouriers law of heat conduction is based on the empirical observation of one- dimensional steady heat flow through a solid. (Steady flow means that the temperature at any point does not vary with time; one-dimensional means that the temperature is uniform over surfaces perpendicular to the direction of heat flow.)
The slab is of Area A Heat is conducted through the slab in only 1 direction (x)
A= area normal to the heat flow direction, m2 dt/dx= temperature gradient or slope of the temperature curve k= thermal conductivity, w/mk Qx/A= Heat flux (heat flow rate per unit area) W/m2 As x increases, the temperature decreases, so the right-hand side of the equation is a positive quantity. The sign ensures that the equation properly indicates that the heat flow is in the direction of temperature fall. Fouriers law of heat conduction is an experimentally observed law and serves as a definition of the property of substances called the thermal conductivity. Solving Fouriers equation for the slab:
qff f f x Z ff = Z @ kdT At x=0 T= T1, at x=L T=T2, constant heat flux, constant k dx A
0 T1
T2
qff ff fZ f x dx = @ k Z dT A
0 T1
T2
k varies over a wide range. It is highest for metals and lowest for finely powdered materials. The thermal conductivities of most liquids are rather small, except for metallic liquids. A few examples of thermal conductivities Material Metals Alloys Water Air Insulators K, W/mK 50-400 10-120 0.598 (at 20oC) 0.0251(at 20 oC) 0.04-0.2
While the rate of heat is transported in a body (q) is dependent on the thermal conductivity (k) and the temperature gradient (T), the rise in temperature that this heat will produce will vary with the specific heat (cp) and the density () of the body.
1.6 Interphase transport (convective heat transfer) - Newtons law of cooling Convective heat transfer occurs when a fluid acts as a carrier or conveyor belt for the energy that it draws from (or delivers to) a solid wall. The characteristics of the flow affect greatly the heat transfer rate between the wall and the stream. The heat flow may be related to the temperature difference between the temperature at the interphase (Tw) and that in the fluid (Tfl)
q = hA T w @T
b
fl
A heated plate with a uniform wall temperature of Tw is immersed in a uniform fluid flow of velocity v V 1 and temperature T 1 . At any location the velocity at the wall is 0 and increases with increasing y to V 1 free stream velocity at some vertical distance away, known as the hydrodynamic boundary layer. The temperature distribution T(y) is also drawn. The temperature decreases from Tw at the wall to T 1 free stream temperature at some distance from the wall. This distance is called thermal boundary layer thickness and is highly dependent on the flow velocity. Heat is transferred from the wall to the fluid.
q = hA T w @ T 1 Where h= average heat transfer coefficient, that is customarily assumed over
b c
the length of plate. In many cases, empirical correlations are available to predict h, since it often cannot be predicted theoretically.
Bulk temperature is also called cup mixing temperature or flow average temperature When a fluid is being heated or cooled, the temperature will vary throughout the cross section of the stream. The bulk temperature is the temperature that would be attained if all the fluid flowing across the section in question were withdrawn and mixed adiabatically to a uniform temperature
c b
Fluid and condition Air in natural convection Superheated steam or air in forced convection Oil in forced convection Water in forced convection Boiling of water
h, W/m2K 5-25
30-300
Condensation of steam
6000-120000
Stephan-Boltzmann equation
Real materials do not emit electromagnetic radiation ideally as black bodies, they have an emissivity
< 1.0 and are called grey bodies. Their emissive power is reduced by
qr = AT4 (<1.0) Because of the 4th power of the absolute temperature, the contribution of radiation to heat transfer is not significant at low temperatures, but becomes predominant at the temperature levels encountered in pyrometallurgical processing.
11
Emissivity
The surface of a substance highly influences its radiation characteristics and therefore the amount of radiative heat the surface can absorb, transmit, reflect and emit. Emissivity is a property that describes how radiant energy interacts with the surface of the material. The emissivity is determined entirely by the properties of the surface of the material and its temperature. Most metals have emissivities that are quite high approaching the black body limit of unity. Some normal emissivity values of various materials at room temperature
Material Aluminum polished Aluminum rough Asphalt Concrete Glass, window Paint black Paint white Red brick skin Soil Water Wood
0.04 0.07 0.85-0.93 0.88-0.93 0.90-0.95 0.98 0.89-0.97 0.93-0.96 0.95 0.93-0.96 0.96 0.80-0.90
Clean, highly polished metallic surfaces have very low emissivities. The emissivity of metallic surfaces increases with increasing temperature. Non metallic materials are characterized by high emissivities in the infrared region. Smooth and rough surfaces have almost the same emissivity. This however holds only for low temperature radiation. The maximum intensity of radiation shifts to greater wavelengths with decreasing temperature. Thus, radiation at low temperature is mainly long-wave radiation. If short-wave radiation, like that of the sun strikes a white surface or ice, the absorptivity is much smaller than when it hits a black surface. The good reflection of sunlight by ice or white fabrics is well known. It is less known that the reflection of the same bodies is very small for long-wave radiation.
Kirchhoffs law states that at the same temperature T1 1 = 1 of a given surface. Net heat transfer between 2 radiating surfaces
In an enclosure all surfaces will simultaneously emit and absorb radiant energy from the other surfaces. The corresponding rates of emission and absorption will depend on the temperature of
12
each surface. When 2 surfaces are at a given distance from each other, the resultant rate of heat transfer by radiation from one to the other may be determined. Two long concentric cylinders and concentric spheres Cylinder (1) with a surface A1 and the temperature T1 encloses cylinder (2) with a surface area A2 and the temperature T2 (T1>T2). The net radiation heat transfer from cylinder (1) to cylinder (2)
b c 1 4 4 ffffffffffff ffffffffffff ffffffffffff f q = ffffffffffffA1 T 1 @ T 2 g 1f Aff 1f ff f1f ff f f f ff f f f f + f f@ 1 1 A 2 2
If plate (1) is at higher temperature than plate (2), then the net heat flux from (1) to (2) is
b c 1 4 4 ffffffffff ffffffffff fffffffff q. . = ffffffffff T 1 @ T 2 f g 1f 1f ff f f f f f f + ff 1 @ 1 2
Lets examine a small gray object (area A1, temperature T1) in a large enclosure at a higher temperature ofT2 so there is a net radiation to the small object. The small object emits an amount of radiation of A1 1 T 1 It absorbs energy from the surroundings at T2 by A1 1 T 2 The net heat of absorption:
4 4
1 at T 2 1 at T 1
qr = A1 1 T 2 @ T 1
4
qr = A1 1 F1 @ 2 T 2 @ T 1
4
qrad = A1 hr T 2 @ T 1
When radiation heat transfer occurs from a surface it is usually accompanied by convective heat transfer unless the surface is in vacuum: qconv = A1 hc T 2 @ T 1 The total heat transfer is the sum of the two:
q = qconv + qrad = hc + hr A1 T 2 @ T 1
Note: The convective heat transfer is not strongly dependent on temperature; however the radiation heat transfer coefficient is a very strong function of temperature.
14
This equation assumes that the surroundings are either very large or black. This equation shows that the temperature indicated by the thermometer is not the true gas temperature but some radiation convection equilibrium temperature. Vey large errors can result temperature measurements if this energy balance is not taken into account. An estimate of the errors in temperature measurement should be made in each installation. The heat transfer coefficient h for the above heat balance can be estimated from: Nu=0.3 Re0.57 where Nu=hd/k= Nusselt number Re=vd/= Reynolds number d= outer diameter of thermocouple well (house), m
Radiation shields are frequently employed to alleviate this difficulty. Radiation shields are made of metals which are highly reflective outside the blackened inside. They do not deliver or remove any heat from the overall system; they only place another resistance in the heat flow path so that the overall heat transfer is reduced. This insulation power depends upon the emissivity of the shield surface. The radiation shield must be designed so that it protects the temperature sensing element of the thermometer from the radiation coming from the walls but does not inhibit convective heat transfer at the thermometer surface. A short, hollow cylinder open on both ends is a desirable configuration. Radiation from the sun
The sun radiates very nearly like a black circular disk with a temperature of 5600K. As a consequence of the high temperature, the maximum radiation intensity is found at 0.5m wavelength and approximately 5 % of the radiation occurs in the ultravoilet (UV) range, 40% in the visible wavelength range and 55% in the infrared (IR) up to ~3m. In a yearly average the earth absorbs ~43% of the radiation coming from the sun (27% directly and 16% directly as diffuse sky radiation), 42% is reflected or refracted back to space from clouds and the air and reflected from the earths surface, and 15% is absorbed in the temperature. The amount of the solar radiation, which a surface absorbs, depends on its absorptivity. Solar radiation has a short wavelength and the absorptivity for such radiation may be considerably different from the absorptivity for long wave radiation .For instance, a white surface has considerably smaller absorptivity than the Al surface for solar radiation.
15
f f f f
5 is a vector, it cannot stand alone but must operate on a scalar, vector or tensor function. It has dimensions of 1/L (1/m)
5=
F ff ff f G ff ff ff ff ff ff ff ff ff
, , x y z
f f f f f f f f f
+e z
rsin
Fouriers law then in the different coordinate systems: Rectangular Cylindrical Spherical
Tf fff ff f q. . x = @ k fff , x
q. . y = @ k Tf fff fff ff f y
Tf fff ff f q. . z = @ k fff z
Thermal energy may be produced by: 1. The degradation of electrical energy (as a result of electric current passing through the material) 2. The conversion of chemical energy into heat (reaction heat of chemical reactions) 3. The Degradation of mechanical energy (viscous dissipation) 4. Nuclear reaction Thermal energy in or out: 1. by heat conduction- consider only this now 2. by virtue of overall fluid motion (sensible heat)
17
Rate of thermal energy in by conduction across surface at x across surface at y across surface at z
yz q. . xLx
L
Rate of thermal energy out by conduction Across surface at x+ x Across the surface at y + y Across the surface at z + z Rate of thermal energy production: Rate of thermal energy accumulation:
L xy q. . zLz
L yz q. . yL y
L xz q. . yLy + y
yz q. . xLx + x
L xy q. . zLz + z
x yz q/
x~ }~ B { ~ ~ ~yz~~y ~~ ~ ~ ~
mass
uf ff ff ff f w t
Substitute Fouriers law in all direction and from thermodynamics: du=dh=cpdT [u=h-pv, differentiating gives du=dh- pdv- vdp. Since in solids p and are constants therefore v=constant so du=dh. From the definition of specific heat dh=cpdT]
18
c p
Tf ff Tf ff ff Tf f Tf fff fff ff f ff ff ff fff fff ff f ff ff fff fff ff ff f ff ff fff fff ff f =@ @k @ @k + @k + q. . . t x x y y z z Tf fff fff ff f =k t
h i T T Tf ffff ffff ffff ffff ffff fffk fff fff fff j fff fff fff
2 2 2
c p
x 2
y2
z2
+ q. . .
This is the general 3-dimensional unsteady-state conduction equation in rectangular coordinate system. In vector notation:
Tf 2 fff ff f c p fff k T + q. . . = t
Where
=
2
x 2
q.fff .f 1fTf 2 fffff f fff f ff f ffff ff. f = T + ffff which is another form of the general 3-dimension unsteady-state t k conduction equation
Simplified forms of the heat conduction equation
No heat source
Tf 2 fff ff f c p fff k T = t
2
Tf 2 fff ff f or fff T = t
T =0
Laplacian equation
4 kinds of conditions are commonly used: 1. The temperature at the surface may be specified e.g. T=Tw 2. The heat flux at a surface may be specified e.g.
q . = qc .
4. At solid-solid interfaces the continuity of temperature and heat flux may be specified. Initial conditions are commonly of 2 types: 1. An initial time condition: the temperature or heat flux may be specified at the start. 2. Time goes to infinity ( for a number of problems when
t =1
20
Tf 2 2 fff ff f c p fff K T + q / [ T = 0 = t
21
In x @ direction:
Rate of flow= driving force / resistance Rate of flow: Driving force: Resistance: qx heat flow rate, W temperature difference, o C thermal resistance, K/W
T= T1-T2
R=/kA
A composite wall consist o 3 materials placed next to each other For each layer we can apply the solution of the slab. qx=T/R
22
T 1 = T 0 @ T 1
T 2 = T 1 @ T 2 T 3 = T 2 @ T 3
R1 = 1 / k 1 A R2 = 2 / k 2 A R3 = 3 / k 3 A
b b c c
Since all the heat passes through the first layer must pass through the second and the third one,
qx1 = qx2 = qx3 = qx = constant T = T 1 + T 2 + T 3 = qx R1 + R 2 + R 3
b c
T difference ffffffffff overallfffffffffffffffffff fffffffff ffffffffffffffffffffffff fffffffff ffffftemperature fffffffff qx = ffffffffff fffffffffffffffffffffffff = fffffffffffffff R1 + R 2 + R 3 thermal resistances P
T ffff3 fffffffffffff ffffo @fffffff ffffff T fffff qx = ffffffffffffff ff fff fff ff fff fff f1ff fff fff fff fff f2ff d 3ff f f f + fff+ fff
k1 A k2 A k3 A
23
Frequently it is required to calculate the heat flow through a wall separating 2 fluids of known temperature, when the surface temperatures of the walls are unknown. The fluid temperatures an be the bulk temperatures, when the fluids flow through ducts; or free-stream temperatures, when the fluid flows are external. Most of the temperature drop in the fluids occurs very near the wall in a relatively stagnant boundary layer (laminar sublayer) which adheres to the wall. In practice the entire resistance to the heat transfer is regarded as being in the laminar sublayer . Despite the fact that the heat flow across this layer by conduction, the mechanism is called convective heat transfer. Newtons law of cooling applies. hhot ,hcold and k are known. h cold and k are known
b c 1 ffff ffff ffff q = hhot A T hot @ T 1 [ T 1 = T hot @ T 1 = q fffff hhot A
24
T hot @ T cold = T = T 1 + T 2 + T 3 = q
fffffffffffffff ffTfffffffffff fffff@T ffffff f fffffff f Therefore q = ffffhotffffcoldf1fff 1 f ffff fffff fff fffff ffff ffff ffff ff ffff ff fffff f ffff + ff+ fffff
hhot A
kA
hcold A
h hot A
kA
h cold A
It is convenient to express the heat transfer rate in terms of an overall heat transfer coefficient for both convection and conduction resistances:
q = UA T hot @ T cold
b c
Where U =overall heat transfer coefficient based on an area A, W/m2K Comparing the last 2 equations:
1 ffffffffffff fffffffffffff fffffffffff fffffffffffff fffffffffff fffff1fffffff fffffffffff ffffffffffff f c for a plane wall U= 1 = b fff f f1ff fff f fff f ff ff f fff f fff f ff fff f A R1 + R2 + R3 + +
h hot k h cold
Steady state, No internal heat generation K= constant, and T=T(r) the temperature varies only with r (radial direction)
25
1fTf 1ff Tf 1f fff fff q/ff T T fffff f ff f fff fff f ff f ff f = ff r fff ffffff ffff fff + fff ffffff ffff fff + 2f fff+ fff+ ff 2 t r r r r k z2
2 2
The temperature distribution is a logarithmic function. To determine the integration constants: B.C 1 B.C 2 at r=R1 at r=R2 T=T1 T=T2
26
T 1 = C 1 ln R1 + C 2 [ C 2 = T 1 @ C 1 ln R1 T 2 = C 1 ln R 2 + C 2 [ C 2 = T 2 @ C 1 ln R1 T 1 @ C 1 ln R1 = T 2 @ C1 ln R2 Tffffff Tffffff ffffff fffff ffffff f@ T f f 1 @T f 1 f ff C 1 = fffff2f C 2 = T 1 @ C1 ln R1 = T 1 @ fffff2fln R1 [ Rff Rff ff ff f f ff f f 1 ln ln f1f
R2 R2
ln b c ln R @T1 ffffff Tffffff ffffff fffff fffff f 1 ffffff fffff fR 1 T = T 1 + T 1 @ T 2 ffffff ffffff= ffffff [ ffffff fffff g g T 1 @T 2 Rff Rff ff f f ff f f ln f1f ln f1f
R2 R2
rff ff ff f f
rf ff ff ff f
2kl
27
A composite cylindrical wall consists of 3 materials placed in series. We can apply the results obtained for cylindrical wall. Using the temperature distributions and thermal resistances given in the next table, the heat flow rate in radial direction in the composite cylindrical wall can be determined
Temperature distribution
Location
d rf ff ff ff f e
Thermal resistance
g Rff ff ff f f 2 R1 ffffff ffffff fffff fffff f
Tffffff ffffff @T 1 R ffffff ffffff ffffff fffff ffffff fffff = f 1g T 2 @T1 Rff ff f f ln f2f
R1
ln
R1 r R2
RT1 =
ln
2k 1 L
f
Tffffff ffffff @T 2 R ffffff ffffff ffffff fffff ffffff fffff = f 2g T 3 @T1 Rff ff f f ln f3f
R2
ln
rf ff ff ff f
R2 r R 3
RT2 =
2k 2 L
f
Tffffff ffffff @T 3 R ffffff ffffff ffffff fffff ffffff fffff = f 3g T 4 @T 3 Rff ff f f ln f4f
R3
ln
rf ff ff ff f
R3 r R 4
RT3 =
ln
2k 3 L
2L T f T 4 @ Tffffff ffffffffff1fffffffffff @T 4 fffffff ffffffffffffffffffff fffffff ffffffffffffffffffff ffffff ffffffffffffffffffff d e d e qr = 1 = d e f2 f ff f f ff ff ff f4 f ff f f 3 X RT i ln Rff ln Rff ln Rff Rf ffffff ffffff ffffff ffffff ffffff ffffff ffffff ffffff ffffff fff1ff fff f ffR 2ff fff f ff ffR 3ff fff f ff i + +
k1 k2 k3
28
In many problems fluids with different temperatures are flowing inside and around tubes and we need to determine the heat transfer rate knowing the bulk fluid temperatures. L= length of tube
= Ro @ Ri = wall thickness =
Rhot = hi Ai
ln
d
Ai = 2LRi
R wall = R cold
T hot @ T cold 2 L f hot @ T ff Tff fffffffffffffffff fffffTfffffffffff fffff fffffffffffffffff ffffffffffffcoldfff ffff ff f d e d qr = fffff ffffffffffffffffff ffffffffffefffffff = fffffffffffffffff ffffffffffffffff = Rff Rff fof ff f f fof ff f f XR ln ln 1 1 ffff ffff ffffff ffff f1 ff ffffff f1ff ffff ffffff ffff ffff ffffff ffff ffff ffR iff f f f ff f ffff ffffff ffff ffff ffR iff fff fffff f fff + fff ff+ ffff ffff+ ffffff ffff + fff
h o Ao 2k L hi Ai ho Ro k h i Ri
29
Often it is desirable to express this equation in terms of an overall heat transfer coefficient U: qr = UA T hot @ T cold
b c
The question is which area to select for A. It is customary to select the outside surface area (Ao) of the cylinder (pipe), but sometimes the inside surface area (Ai) is chosen. Comparing the last 2 equations, the overall heat transfer coefficient Uo based on Ao is
Uo =
R o ln
h i Ri
Uo =
R i ln
hi
To= outside tube wall temperature Ro= outside tube radius T=ambient fluid temperature The overall resistance is the sum o the cylindrical shell resistance of the insulation and the external convective heat transfer resistance:
30
e Rfff fff ff ff f insf ln Rf 1 ffffff ffffffffff ffffff ffffffffff fffo ff fffffffff c R = fffffff bffffffffff +
2kL
2Rins L ho
b
2 T O @ T 1 Tfffffff qf fffffffffffff ffffff f@ffff ff T 1 f ffffffffffff d e q = f0ffffff f= fffffffffffff [ f fffffffffffff f ff ff ff R L ln Rffff ins 1 fffffff fffff ffffff ffffff ffffff fffff ffR o ff fffff ff ffff + ffff
k h o R ins
Rins
=0
Rins
=0
Rff fff1fffff fffff 1ff 1 ffff fffffffff fffffffff ffff fffffffff fffffffff fff ff ffffffff f fffff ff f =` a fff@ h 2 L 2 1ff ff ff Rins 2k L fff o Rins R ins
Solving for the critical insulation radius for the cylinder Rins crit
f g 1 1 fffffffff 1f fffffff fffffffff f fffffff fff f1fff f fffffff fff ffff f fffffff f
2 L Rins crit k
kf ff f Rins crit = ff ho
ho Rins crit
The critical insulation radius corresponds to maximum heat transferred. If Rins< Rins crit , then adding insulation increases the heat transferred. In case of a thick bare cylinder: Ro> Rins crit , the addition of insulation always increases the insulation effect. In case of a thin bare cylinder: Ro< Rins crit , the addition of the first layer of insulation decreses the overall thermal resistance R, so enhances heat transfer, not insulating until Rins>Rins crit. (Eg. A wire insulated with a covering of not too low thermal conductivity may run cooler, for a given current, than the same wire if bare. A practical application is the problem of insulating electrical wires where the objective would be the provision of adequate electrical insulation at the same time providing the maximum wire cooling.) The critical insulation radius for a sphere is Rins crit= 2(k/ho)
31
In 2-dimension:
For relatively simple geometries the above differential equation can be solved, but for more complicated systems we need numerical or graphical techniques. There are 2 types of resistances to the heat transport:
We can identify 3 types of problems that arise for a transient problem. For some problems the internal resistance is negligible compared to the surface resistance. For other problems the surface resistance is negligible. Still, for others both are important.
The temperature profile within the object is uniform and everywhere equal to the surface temperature, At any instant in time there are no temperature gradients within the object, only a temperature difference between the surface and the free stream The heat transferred between the object and the fluid is controlled by the surface resistance (ie the value of the heat transfer coefficient h)
33
The heat transfer process is called Newtonian cooling (or heating). Because the temperature does not change (significantly) withi the solid, we can formulate a model by performing an energy balance on the object.
This equation describes the time-temperature history of the solid object. It is called the lumped capacitance equation.
cp V -
lumped capacitance
cpV/hA - the characteristic time of decay, or thermal time constant, or decay time constant.
It takes a longer time for the body to reach equilibrium with the surroundings fluid when its lumped capacitance is large and / or its product hA is small. Examining the exponent of the lump capacitance equation: The ratio of volume to surface area for an object is often used as a characteristic length Lch
Where
Fourier number is a dimensionless time parameter. It presents the ratio of heat transfer by conduction to the energy storage rate within the material. Thus the lump capacitance equation becomes:
35
The assumption of negligible internal resistance is reasonably accurate when Bi < 0.1. This is the criterion for the lump- capacitance equation to apply.
At any time, the instantaneous rate of heat transfer q(t) in Watt from the body is
Where T= instantaneous temperature. The instantaneous temperature from the lump capacitance equation is
This implies:
The amout of heat Q in Joules transferred from the solid from time t=0 to t= t is
36
Consider an infinite plate of thickness 2L at an initial uniform temperature Ti. At some time zero, the left and right surfaces of the plate are suddenly cooled to T1 by a fluid at a temperature T. The thickness of the plate is 2L. The aim is to predict the temperature profile in the plate with time, an the heat transfer rate. Fourier solved this problem by developing the method of Fourier series.) If the heat transfer coefficient h is very large then the film resistance (1/hA) is negligible so the surface temperature, once exposed, will assume the same temperature as the fluid: T1=T which holds when Bi>1. The descriptive differential equation:
37
=(Bi, Fo, x/L) the temperature distribution is a function of Fourier and Biot numbers and x/L.
The analytical solution got by the separation of variables method:
38
This equation is for the case where the left and right faces of an infinite wall are maintained at the same temperature Ti. The centerline of the plate is a line of symmetry and therefore analogous to an insulated or adiabatic surface. The temperature history at the midplane (x=L) is of practical interest and is a function of Fourier number:
For values of Fo much greater than 0.1, only the first term in the series need to be used in the solution (except at points very close to the boundaries).This makes possible a graphical presentation . On a chart (Fig 6.6), knowing Fo number the dimensionless temperature for the center plane can be read for different shapes. An infinite plate is a plate whose thickness is small in relation to the other dimensions. An infinite cylinder is a cylinder where the diameter is small compared to the length.
The instantaneous heat flow rate (q) and the amout of heat transferred (q) for the infinite plate are :
39
On a chart (fig 6.7) q (in watt) and Q (in Joule) can be read as a function of Fo
A similar analysis can be formulated for an infinitely long cylinder , a semi-infinite solid, a sphere and other types of common geometries. The solutions are conveniently summarized graphically on the above charts.
Heisler charts
(1) The dimensionless temperature history at the center of an infinite plate of thickness.
40
This chart can be used to calculate the midplane temperature from Fourier number, which in turn can be used in the next chart to determine the temperature at any other x location within the plate.
(2) The dimensionless temperature distribution in an infinite plate at any x location within the plate is given in Fig 3.9
Q= the total amount of energy (J) that has passed through the wall up to any time t
41
Qo = the initial internal energy of the wall relative to the fluid temperature T( it is equal to the maximum internal energy change that can occur) Qo = cpV(Ti-T) Similar charts for infinite cylinder (L/R > or = 10 at least; Lch=R) are also available. Note: the characteristic length (Lch) used in the lumped-capacitance method is different from the characteristic length used in the Heisler and Grober charts. For the negligible surface resistance case the heat transfer coefficient is large (ie. H ). For this case 1/Bi=k/hLch=0. This is also represented on the charts.
A solid body occupying the space from x=0 to x= is initially at temperature Ti. At time T=0, the surface at x=0 is suddenly raised to temperature Ts and maintained at the temperature fot t>0 constant wall temperature. The solution to this problem appeared in fig 6.7 for the case of negligible surface resistance ( h is high Ts =T). Lets examine the analytical solution to the same problem.
42
Assumptions:
heat conduction occurs only in x direction, no heat generation, constant physical properties.
This integral has no closed form. The right hand side is called Gauss error function or probability integral and denoted by error function of lamda . This integral is a function of its upper limit. Its values are available in tables. Some properties of the error function:
Therefore:
43
The error function reaches a value of 0.99 when the argument is about 2, therefore we can define a thermal penetration thickness T as that didtance x for which the temperature dropped to a value of 0.001(Ts-Ti)
For distances x> T , the initial temperature is changed by less than 1 % of the difference ( TsTi)
44
When T is small with respect to the plate thickness L, then the semi-infinite solid solution can be used. When V is of the order of magnitude of plate thickness L, then the finite slab solution (rectangular parallelopipedon) has to be used. The heat flux is
The total amount of heat (Q) transferred from the solid from t=0 to t=t is
The problem can be redefined and the surface convection effect or constant heat flux at the surface could be accounted for. The next figure compares the solutions the 3 types of possible boundary conditions.
45
2(y,t)
For the purpose of simplifying the notation, we shall use the following symbols:
46
The following figures show how the solutions to 2 one dimensional problems are combined to obtain the solution to a 2-dimensional problem. Other combinations can be made. In each case the solid exists at an initial temperature Ti, and at the time zero it is subjected to a fluid temperature T with a heat transfer coefficient h.
Finite cylinder
47
Semi-infinite plate
48
This body can be viewed as the intersection of 3 infinite plates that are perpendicular.
5. CONVECTION
In most transfer processes heat transfer in fluids is accompanied by some fluid motion so that the heat transfer does not occur by conduction. There are 2 types of heat transport involving fluids:
1. Forced convection
The e fluid flow is established by an external driving force (pressure difference provided by a pump or gravitational force). The method of solution: First the velocity profile is found, th en it is used to derive the temperature profile.
49
The velocity profile and temperature profile are intimately related and therefore determined together. Whether the heat transfer mechanism is natural or forced convection, the fluid motion can be described by the equations of fluid mechanics. At low velocities the flow is laminar throughout the system; at high velocities it is laminar near the heating surface and turbulent at some distance away. All problems of convective heat transfer can be expressed in terms o f differential mass, energy, and momentum balance (equations of change). However, the mathematical difficulties with the integration of these simultaneous non-linear partial differential equations are such that analytical solutions exist only for relatively simple problems. In general, we have 2 types of flow situation: 1. External flows- they are not enclosed but have contact with a solid surface. 2. Internal flows- they are enclosed completely by solid surfaces (flow in a pipe or duct)
50
It is more convenient to have the thermal energy equation in terms of fluid properties ( T cp). If , , and k are constants, then the thermal energy equation (Table 7.3) becomes
51
Where v= viscous dissipation function . It is usually negligible; it is needed only in special situations. Simplification of the equation of thermal energy: For ideal gas
52
For very simple geometry, these equations can be solved, but for most of the situations we cannot. The refore alternative solution methods have been devised such as boundary layer theory, dimensional analysis
These equations were considerably simplified by Prandtl. He recognized in 1904 that the effect of the solid boundary on the flow is confined too a very thin layer of fluid immediately adjacent to the solid wall and can be considered negligible beyond d it. This layer is called the boundary layer.
53
When a fluid flows over a surface the elements in contact with the surface will be brought to rest and the adjacent layers layers retarded by the viscous drag of the fluid. Thus, the velocity in the neighborhood of the surface will change with distance perpendicular to the flow. Suppose a fluid approaches the plate with uniform v (or free stream velocity) and T temperature when the fluid reaches the surface, a velocity gradient is set up because of the viscous forces acting within the liquid. The fluid in contact with the surface will be brought to rest and will gradually approach the free stream velocity v at some distance from the surface. We divide the flow regions into 2 parts: (1) A non viscous region away from any solid surfaces (ideal fluid: =0, =constant) (2) A boundary layer (the fluid adheres to the surface due to viscous effects)
The thickness of the layer in which the fluid is retarded becomes greater with distance in the direction of flow. This layer was termed boundary layer by Prandtl. Anywhere along the plate the velocity profile changes from the uniform free- stream profile. The velocity increases from zero at the surface to nearly the free stream value at some distance away . The free stream velocity is approached asymptotically, and therefore the boundary layer strictly has no precise outer limit. However, it is convenient to define the boundary layer thickness such that the velocity at its outer edge equals 99% of the free stream velocity. That is the hydrodynamic boundary layer thickness is where Vx is approximately 0.99V. Within the boundary layer the velocity in x direction varies only with y: Vx=Vx(y)
Similar behavior occurs with temperature. The thermal boundary layer thickness T where the temperature approaches T: T is approximately T. The boundary layer grows from the leading edge with distance x. Typically , and t are a few thousandths of a millimeter thick. Despite their size, both are important in determining drag forces and heat transfer rates.
54
The mathematical simplifications for the equations of change: 1. The boundary layer is thin compared to the distance measured from the leading edge. 2. As a result of (1), d2Vx/dx2 is much smaller than d2Vx/dy2 3. vy<<vx (so all terms in the Navier stokes equation y-component which involve Vy may be neglected). 4. The pressure across the boundary layer is assumed to be constant: dp/dy is approx 0 5. There is no pressure gradient or gravity force in x direction 6. There is negligible heat conduction in the direction of flow: d2T/dx2<< d2T/dy2 The equations for the boundary layer become:
Solutions for the temperature profiles depend on the boundary conditions specified (heating could begin at the leading edge of the plate or somewhere downstream, qw or Tw can be constant)
There are 2 strategies for solving the boundary momentum equation for the velocity profile: 1. The exact solution by Blasius (1908),
55
2. The approximate solution from the integral method developed by von Karman and Pohlhausen (1921)
where
The Reynolds number based on the downstream distance x. That is , the thickness of the boundary layer is proportional to the square root of the kinematic viscosity () and the square root of the downstream distance x. Some literature gives the constant 5.0 instead of 4.96.
these much simpler equations do not reveal anything new about the temperature and velocity profiles, but they do give accurate explicit equations for w and qw. The solution leads to the boundary layer thickness
The result is in close agreement with that obtained in the exact solution.
The total skin friction drag obtained by integrating the wall shear stress over the plate area ( with W as the plate width):
The physical significance of the Prandtl number of a fluid is that it represents the ratio of the transport coefficient for viscous momentum transport to the transport coefficient for heat conduction.
57
When the Prandtl is large, as in the case of slags, the ability of the fluid to transport momentum is greater than that to transport thermal energy. Under these conditions, the hydrodynamic boundary layer is greater than the thermal boundary layer. >T. For Prandtl numbers much smaller than unity, as in the case of liquid metals, the hydrodynamic boundary layer is relatively small and we neglect the effect of the solid surface on the velocity distribution when we formulate the thermal energy balance; that is we may ignore any velocity gradients and assume that all the fluid moves at the free stream velocity. When Prandtl is close to unity, as in the case with most gases, the hydrodynamic and thermal boundary layers practically coincide. =T
(a) For Pr=1 The equation of energy is similar to the equation of motion. The boundary conditions are identical for dimensionless velocity and dimensionless temperature. If =, then the differential equations are identical, =T, so each problem must have the same solution( i.e. the dimensionless velocity profile solution is the same as the dimensionless temperature profile solution and for any point in the flow system, the dimensionless velocity is equal to the dimensionless temperature):
In the case we can immediately calculate the heat transfer coefficient from the heat flux at the wall from Blasius exact solution:
58
It is customary to express the heat transfer coefficient in terms of dimensionless group Nux.
(b) For Pr not equal to 1 From the integral method of solution Pohlhausen found that
So the heat flux at the wall, local heat transfer coefficient and the local Nusselt number is:
59
These equations give the local heat transfer coefficient (hx) and the dimensionless local heat transfer coefficient (Nux) at a distance x from the leading edge of the plate. The average heat transfer coefficient h (or h-) for the entire surface is
Since the analysis was based on constant fluid properties, which, in real situations are not, they depend on the local temperature, fluid temperature; fluid properties have to be evaluated at the average temperature of the fluid in the boundary layer:
60
As the boundary layer thickens, at distances farther from the leading edge, a point is reached where turbulence appears. The onset of turbulence is characterized by a sudden rapid increase in the thickness of the boundary layer. In turbulent region, however, there is still a laminar sublayer near the surface. Transition usually occurs over a range 2 x 10 5 < or = Rex < or =3 x 106 and not a single point. For purposes of calculation it is customarily assumed that transition occurs at Rex= 5 x 105. The velocity profiles in laminar and turbulent flows are different, the turbulent profile is flatter over a greater portion of the boundary layer due to mixing effect. It is not possible to solve the boundary layer equations to obtain an exact solution. However, the integral method gives very good results. The hydrodynamic boundary layer thickness is
61
Although the solution for the differential energy balance in the laminar flow regime is available, engineers prefer to work with heat transfer coefficients. Ass a result, it is customary to express heat transfer performance, even for laminar flow, in terms of the heat transfer coefficient. The heat transfer coefficients for the laminar flow have a strong dependence on position. This is not usually the case for the heat transfer with turbulent flow. There are 2 objectives 1. To relate the heat transferred to the temperature change experienced by the fluid, 2. To relate the heat transferred to the fluid temperature between the wall a nad the fluid
62
Therefore for Pr> 1 the flow develops hydrodynamically more rapidly than it does thermally. The local heat transfer coefficient hz is
Tw-Tb is the temperature difference between the wall and the fluid. It is an important quantity. Both Tw and Tb can vary in the z-direction. If we integrate the local coefficient hz over a definite area (such as the surface area As), the integrated result would be h- average heat transfer coefficient.
6.2 Heat transfer in laminar flow- constant heat flux at the wall
The heat flux at the wall qw = q/A = constant. Experimentally this can be achieved by passing electric current through a metal pipe by wrapping the duct with a material through which an electric current is passed. The fluid temperature at the inlet is T and varies within the tube. At any axial location z, T varies from Tw to some value at the center line. The temperature profile reaches a non-changing shape (thermally developed flow), but with the constant heat input that exists, T-T increases with increasing axial distance z. The solution to the problem:
63
The is the dimensionless heat transfer coefficient based on the local heat transfer coefficient is:
64
The following graph summarizes the results obtained from the solution of the differential equation of energy. The temperature profile is thermally developed when 1/Gz> 0.05
65
It applies to steady laminar flow of Newtonian fluid in a tube having constant wall temperature under the following conditions:
Heat transfer coefficients are higher with turbulent flow than with laminar flow. It is not possible to obtain a closed form solution for the velocity profile in turbulent flow. Dimensional analysis predicts that the Nusselt number depends on Reynolds and Prandtl numbers. Combining with experimental results, the following is found to apply:
Sieder-Tate equation:
This correlation is applicable in case of moderate temperature differences between the wall and the fluid. This equation should be used for water. The use of this equation may be trial and error, since value of hL must be known to evaluate Tw, and hence w, at the wall temperature. Since the rate of heat transfer is greater I turbulent flow than in laminar flow, most equipment are operated in the turbulent range. In place of this type of correlation another is often used, which consists of using a dimensionless group j called the Colburn factor:
67
These correlations for the turbulent regime predict the local Nusselt number. The average Nusselt number based on a heat transfer coefficient for the entire tube is not available. When dealing with non-circular conduits, the hydraulic diameter or equivalent diameter is used in the heat transfer correlations whenever D diameter appears.
68
7. HEAT EXCHANGERS
A heat exchanger is a device designed for exchanging heat between two fluids.
the amount of heat transferred the outlet temperature of one or both fluid streams , the required the heat transfer area , depending on what is known.
It is usually assumed, that all the heat lost by the hot fluid is transferred to the cold fluid (that is there are no heat losses).
69
In the counter-flow arrangement t2 ( cold fluid outlet temperature) can be greater than T2 (hot fluid outlet temperature). In parallel flow, however, the upper limit of t2 is T2
70
One fluid flows through the center pipe while another flows in the annulus. Such an exchanger may consist of several passes arranged in a vertical stack. Effective lengths: -12ft (3.65m)
It is useful, when not more, than 9-14m2 of heat transfer surface is required It is preferred for small capacity, high pressure and countercurrent operation For low flowrates,
71
If the heat duty (q) is less than 500kW, If for the shell and tube heat exchanger the number of tubes are less than about 30(3/4in=19mm OD tubes) and the diameter of the shell is less than 200mm.
Under these conditions a shell-and-tube heat exchanger is uneconomic; a large number of small diameter shells in series are required to provide adequate velocities and near-countercurrent flow. Disadvantages -Each double-pipe exchanger introduces no fewer than 14 points at which leakage might occur - Considerable time and expense required for dismantling and periodically cleaning it compared with the other types of equipment. -when a large number are connected, they require considerable space.
72
The number of tubes in the tube bundle ranges from about 20 to over 1000. The tubes are attached to perforated flat plates, termed tube sheets, at each end. The tubes pass through a number of flat plates, called baffles, along their length, which serve to support them and direct the fluid flow in the shell in such a way that the heat transfer is enhanced. The assembly of tubes and baffles, called the tie bundle is held together by a system of tie rods and spacer tubes.
73
A mixture of condensable and non-condensable substances (e.g. steam and air) A mixture of 2 or more condensable vapors.
The condensing temperature of a single pure substance depends on ly on the pressure therefore the process is isothermal. Mixed vapors condensing at constant pressure condense over a temperature range and yield a condensate of variable composition until the entire vapor stream is condensed. There are 2 types of condensation -dropwise condensation -film-type condensation
Dropwise condensation
The condensate begins to form at microscopic nucleation sites (eg tiny pits, scratches, dust, specks ). The drops grow and coalesce with their neighbors to form visible fine drops. The fine drops, in turn, coalesce into small streams, which flow down the surface under the force of gravity, sweep away condensate, and clear the surface for more droplets. During drop wise condensation, large portion of the surface area is covered with an extremely thin film of liquid of negligible thermal resistance, consequently the heat transfer coefficient at these bare areas is very high. The average coefficient for dropwise condensation may be 5-10 times that for the film-type condensation and can be as high as 110 000W/m2K. Dropwise condensation occurs when the condensate does not wet the surface (e.g. ethylene glycol, glycerin, nitrobenzene, isoheptane, steam). This type of condensation is so unstable and is so difficult to maintain (because the surfaces become wetted after prolonged exposure to a condensing vapor) that the method is not uncommon. For normal design purposes film-type condensation is assumed.
74
liquid between the surface and vapour that forms the main resistance to heat transfer. The thickness of the film increase rapidly in the first few centimeters and then more and more slowly. Condensing film coefficients are much greater than those in forced convection and are of the order of magnitude of several thousand of W/m2K. The film may flow in the laminar or in the turbulent regime, depending upon: -the rate of condensation -fluid properties -length of path of the growing condensate film, and -geometry The usual design for condensers employs tubes, with the condensate formed inside or outside of the tubes. Pure vapors are usually condensed on the outside of tubes (shell side). These same basic design equations used for heat exchangers are valid, but the nature of the coefficient correlations is somewhat different, depending upon whether the condensation process takes place on vertical or horizontal tubes. Friction losses in a condenser are normally small, so that condensation is a constant pressure process. The basic equations for the rate of heat transfer in film-type condensation were first developed by Nusselt (1916)
75
On a vertical flat surface of width W a vapor is condensing in a film-type manner. The film is 2dimensional and has a thickness of at any z location [=(z)]. Nusselt Assumptions: 1. Pure vapor is at its Tsat saturation temperature 2. The condensate film flows in laminar flow 3. Heat is transferred through the film solely by conduction ( 4. The temperature distribution through the film is linear. 5. The temperature of the condensing surface Tw is constant. 6. The physical properties of the condensate are constant and evaluated at a mean film temperature Tf hz = k/)
Velocity distribution
The velocity distribution in the falling film can be determined from the Navier-Stokes equation as
77
78
79
80
81
82
PART 2 1. INTTRODUCTION
A large number of the unit operations used in the recovery and purification of metals are concerned with the problem of changing the composition of solutions and mixtures through methods not involving chemical reactions. These operations separate mixtures into its component parts. Such operations may be entirely mechanical (e.g. filtration, screening). If the operations involve changes in composition of solutions, they are called mass transfer operations. The mass transfer operations are characterized by transfer of a substance through another on a molecular scale. The mass transfer is a result of concentration difference, or gradient (which is the driving force), the diffusing substance moving from a place of high to one of low concentration. Mass transfer can occur in the same phase orfrom one phase to the other. Mass transfer operations include distillation, gas absorption, adsorption, dehumidification, liquid extraction, leaching, crystallization, drying, (dialysis, permeation, electrodialysis, etc but these are non-equilibrium operations). The operations may be conducted in either batchwise or continuous manner. The function of the operations: Distillation: is to separate, by vaporization, a liquid of miscible and volatile substances into individual components or into groups of components (e.g mixture of alcohol and water, liquid air into N2, O2 Ar, crude oil into LPG, petrol, diesel, kerosene, fuel oil,) Gas adsorption: a soluble vapor is absorbed by means of a liquid in which the solute gas is more or less soluble, from its mixture with an inert gas (e.g. washing of NH3 from a mixture of NH3 air by liquid water). The solute is then recovered from the liquid by distillation. Adsorption: A solute is removed from either a liquid or gas through contact with a solid adsorbent, the surface of which has a special affinity for the solute. Dehumidification: a pure liquid is partially removed from an inert or carrier gas by condensation (e.g. H2O vapor of a stream of N2 by condensation). Liquid extraction (solvent extraction): a mixture of 2 components is treated by a solvent that dissolves 1 or more of the components in the mixture Leaching (extraction of solids): a soluble material is dissolved from its mixture with an inert solid by means of a liquid solvent.
83
Auxilliary operations -heat addition orremoval ( to change temperature or phase conditions) -work addition or removal ( to change pressure) - mixing or dividing, -Solids separation by size Mass transfer problems can be solved by 2 distinctly different methods: -based on diffusional rate processes, -utilizing the concept of equilibrium stages. The choice of method depends on the kind of equipment in which the operation is carried out. All mass transfer calculations, however, involve knowledge of the equilibrium relationships between phases. You have already seen that analogy between momentum and energy transfer, we will extend the analogy to include mass transfer. If the mixture is stagnant, the transfer occurs through stagnant layers of solid or fluid by molecular diffusion. This mechanism is analogous to the heat transfer by conduction. If there is a bulk mixing of phases or layers of fluids by mechanical stirring or because of a density gradient, mass transfer occurs by the mechanism of forced or natural convection. This is called eddy diffusion. Sometimes the diffusion process is accompanied by bulk flow of the mixture in a direction parallel to the direction of diffusion and it is often associated with heat flow. We will consider only binary mixtures
84
2. MOLECULAR DIFFUSION 2.1 Molecular Diffusion is the movement of individual molecules through a substance by virtue of their thermal energy.
The kinetic theory of gases provides a means of visualizing what occurs. A molecule is imagined to travel in a straight line at a uniform velocity until it collides with another molecule, whereupon its velocity changes n both magnitude and direction. The average distance the molecule travels between collisions is its mean free path, and the average velocity depends on the temperature. The molecule thus travels in a zig-zag path, the net distance in one direction which it moves in a given time, the rate of diffusion, being only a small fraction of the length of its actual path. For this reason, the diffusion rate is very slow. (To demonstrate the importance of the barrier molecular collision presents to diffusive movement: the rate of evaporation of water at 25oc into complete vacuum is approximately 3.3kg/s per 1m2 of water surface. But placing a 0.1 mm thick layer of stagnant air at 1 atm pressure above the water surface reduces the rate by a factor of approximately 600) The same general mechanism prevails also for the liquid state, but because of the considerably higher molecular concentration, the diffusion rates are even slower than in gases. (e.g. A 1.5m diameter tank is filled with a salt solution to a depth of 0.75m. A 0.75m deep layer of pure water is carefully placed over the brine without disturbing the brine in any way. If the contents of the tank are left completely undisturbed, the salt will, by molecular diffusion, completely penetrate the liquid, ultimately coming everywhere to one-half its concentration in the original brine. But the process is very slow. It can be calculated that the salt concentration at the top surface will still be only 87.5% of its final value after 10 years and will reach 99% of is final value after 28 years. On the other hand, a simple paddle agitator rotating in the tank at 22rpm will bring about complete uniformity in approximately 60s.) Diffusion also takes place in fluids by physical mixing and by eddies of turbulent flow. This method of solute transfer is called eddy or turbulent diffusion. Within each eddy, no matter how small, uniformity is achieved only by molecular diffusion, which is the ultimate process. The analogy between heat and mass transfer is so straightforward that equations developed for the former are often found to apply to the latter by simply changing the meaning of the symbols. A concentration gradient tends to move the compound in such direction as to equalize concentrations and destroy the gradient. When the gradient is maintained by constantly supplying the diffusing component to the high-concentration end of the gradient and removing it at the low-concentration end, the flow of the diffusing component is continuous. In a two-phase system not at equilibrium, sum as layer of NH3-air as a gas solution in contact with a layer of liquid water, spontaneous alteration through molecular diffusion also occurs,
85
ultimately bringing the entire system to a state of equilibrium whereupon alteration stops. At the end, the concentration of any constituent is the same throughout a phase, but it will not be the same in both phases. The chemical potential ( or activity ) of the components, which is differently dependent on concentration in the 2 phases, will be uniform everywhere throughout the system at equilibrium and this brings the diffusive process to a halt. In all mass transfer operations, diffusion occurs in at least 1 phase and often in both phases (e.g. In gas absorption, solute diffuses through the gas phase to the interface between the phases and through the liquid phase from the interface.
86
I In a diffusing mixture the various chemical species are moving at different velocities. Let vi denote the velocity of species I with respect to stationary coordinate axes. Then, for a mixture the local mass average velocity v is defined as
87
This is the local velocity one would measure by a Pitot tube and identical to the bulk flow velocity v as used in for pure fluids.
When the whole fluid is moving in convective flow and component A is diffusing in the same direction as the flow, its diffusion velocity is measured relative to the moving fluid. The velocity of a particular species relative to the molar average velocity (v*) is called the diffusion velocity (vAd ). To a stationary observer, A is moving faster than the bulk of the phase, since the diffusion velocity vAd is added to that of the bulk phase molar average velocity v*. Expressed mathematically, the velocity of A relative to the stationary point VA is then the sum of the diffusion vAd and the molar average velocity v*.
88
The mass or molar flux of a species is a vector quantity denoting the mass or moles of the species that passes through a unit surface area per unit time. The motion may be referred to as stationary coordinates or to the local molar average velocity v*
In any of the equations derived for the integral or differential balances, which apply to mixtures, the velocity used has been the mass average velocity. Therefore, the mass flux is not equal to the molar flux N times the average molecular mass M, because the mass average velocity and the molar average velocity are different. ( In the flow of a uniform mixture through a pipe, there is no interdiffusion effect , vA and vB are identical, and the mas flux is then equal to he molar flux N multiplied by the average molecular mass M. However, in cases foe which there is a concentration gradient in the direction of flow, vA and vB differ with effect noted above.)
89
When the mole fraction xA or mass fraction A are small, the flux with respect to stationary coordinates (NA, nA) is approximately equal to the flux relative to the average velocity of the fluid due to diffusion (JA, jA)
90
Ficks first law states that the species A diffuses (moves relative to the mixture) in a direction of decreasing concentration of A, just as heat flows by conduction in the direction of decresing temperature. Ficks law is the definition of diffusivity DAB . For steady state conditions the net flux is NA+NB=N It is usually convenient to describe the mass flux by diffusion (diffusion flux) relative to a stationary coordinate system, rather than the average molar velocity of the fluid. This can be done by the following consideration: -The movement of A carried by the convective bulk fluid flow, -The movement of a carried by molecular diffusion. NA= NxA+JA since xA=cA/c and N =cvz
All the above has considered diffusion in only 1 direction , but in general, concentration gradients , velocities, and diffusional fluxes exist in all direction.
91
2.4 Diffusivity
Diffusivities are best established by experimental measurements. Often the desired values are available and they must be estimated from published correlations.
Gases
DAB values fall in the range of 0.1- 1.0 cm2/s. The diffusivities can be considered independent of concentration at pressures below ~5atm. Diffusivities in gases can be predicted with considerable accuracy from kinetic theory. The diffusivity of gases varies inversely with pressure:
Liquids
The theory of diffusion in liquids is not advanced nor the experimental data as plentiful as gas diffusion. DAB ranges 10-5 to 10-6cm2/s and it varies appreciably with the concentration and temperature. Diffusivities for dilute solutions can be calculated approximately from formulas presented in literature. The diffusivity in concentrated solutions differs from that in dilute solutions because of changes in viscosity with concentration and also because of changes in the degree of non-ideality of the solution. Unlike the case for binary gas mixtures the diffusivity for a dilute solution of A in B is not the same as for a dilute solution of B in A. For intermediate concentrations, an approximate value of D is sometimes obtained, but this method can lead to large errors for non-ideal solutions. For strong electrolytes dissolved in water, the diffusion rates are those of the individual ions, which move rapidly than the large, undissociated molecules. Solids
92
One component in a solid will diffuse through the other at a measurable rate if there is a suitable concentration gradient and the temperature is high enough. The effect of diffusion in solids is very important in metallurgy and ceramic industry. The rate of reaction in some chemical processes is determined by solid diffusion. The atoms in a crystalline solid are localized at certain positions in the crystal lattice. An atom vibrates about it equilibrium position, and occasionally a vibration becomes strong enough to cause the atom to escape its original position and jump to a new place. Its energy is low at the lattice position but higher in between; the atom needs certain activation energy to jump over this energy barrier. The diffusion process in solids may be strongly nonisotropic. DAB values are given in reference books. It is of the order of 10-8 to 10-10 cm2/s. For concentrated systems the diffusivity is a strong function of composition and temperature.
93
If there are no chemical reactions (steady stse, since generation of A is zero), and if the bulk flow through the volume element is negligible, the component balance equation for A results in a second-order differential equation in terms of the concentration:
3.2 Diffusion through a stagnant fluid (one-component mass transfer or one-way diffusion)
In several important processes, 1 component in a gaseos mixture will be transported relative to a fixed plane (say, a liquid interface) and the other will undergo no net movement. In gas absorption a soluble gas A will be ranferred to a liquid surface where it dissolves, whereas the insoluble gas B wll undergo no net movement with respect to the interface. In evaporation from a free surface, the vapor will move away from the surface but the air will have kno movement.
Consider liquid A (water) is evaporating into gas B (air), and we maintain the liquid level at z=z1. The vapor is diffusing at a constant rate (steady-state0 from the liquid surface up through the layer of stagnant air in the tube. This implies that NB =0
94
95
96
97
kc may apply to forced or natural convection. Its value is a function of the geometry of the system, the location along the surface, the fluid velocity , the properties of the fluid , and concentration. For unsteady state situation, the mass transfer coefficient may also depend on time. NA is defined a the flux at the phase interface, or boundary, wher substance leaves or enters the phase for which kc is the mass transfer coefficient.
98
Since the convective flow term has been ignored in the solution, these mass-transfer coeeficients are useful only for low mass transfer rates. Comparing the rate equation with the solution of Transfer of A through stationary B
In a few limited situations mass transfer coefficients can be deduced from theoretical principles. In th eegreat majority of cases however, we depend on direct measurement under known conditions, for use later in design.
99
When a fluid flows turbulently past a solid surface, with mass transfer occurring from the surface to the fluid, the concentration distance relation ia sa shown by the full curve. The shape of the curve is controlled by ED eddy diffusivity to DAB. Instead of attempting to solve the equation of the fluid flow, we use film theory. The film theory assumes that the concentration will follow the broken curve of the figure, such that the entire concentration difference (cA1-cA2) is attributed to molecular diffusion within an effective film of thickness . Thus, the the concentration at the outer edge of this film is cA2.That is the resistance to diffusion in the entire flowing fluid is considered equivalent to that in a stagnant film of thickness. (This is a fictitious film; its thickness cannot be determined by any experiment, since it does not exist.) The problem is now reduced to steady-state diffusion through a stagnant film of thickness , because of a concentration difference of (cA1-cA2). The problem was solved to give the molar flux of A as
100
That is the k-type mass transfer of coefficients for different solutes are directly proportional to the DABs for the solutes. Unfortunately this last equation does not allow a direct way to evaluate kc, since the thickness of this hypothetical film , is a function of the same parameters as the mass transfer coefficient itself: velocity and properties of the fluid, geometry, concentration levels etc. The film theory is in conflict with our knowledge of turbulent flow, but does well in case of high mass transfer flux, the effect of mass transfer on heat transfer, and in predicting the effect o reaction rate on mass transfer. (Although there exists a laminar sublayer next to the wall in turbulent flow, this is not equivalent to a laminar-flowing film, since only part of the concentration gradient is in the laminar sublayer and all of it is postulated to occur across a laminar film. So this fictitious thickness of laminar-flowing fluid offers the same resistance to transfer as actually exists in the entire turbulently flowing fluid. The thickness of this fictitious film extends beyond the laminar sublayer to allow for the change in concentration in the buffer layer and turbulent core. Obviously can never be measured, since it does not exist.) Penetration theory In many situations the time of exposure of a fluid to mass transfer is short, so that the concentration gradient of the film theory, characteristic of steady state, would not have time to develop. Higbie (1935) described the contact of 2 fluids with this theory.
101
102
This independence on DAB is typical of short exposure times, where the depth of solute penetration is small relative to the depth of absorbing pool. Difficulty: tT is not usually known. The depth of penetration, defined as the distance at which the concentration change is 1% of the final value is
Surface renewal theory The Higbie theory assumed a constant time of the eddies of fluid at the surface. Danckwerts gave more realistic picture, where the eddies are exposed for varying length of time. The liquid-gas interface is then a mosaic of surface elements of different exposure-time histories. Since the rate of solute penetration depends on exposure time, the average rate for a unit surface area must be determined by summing up the individual values.
A fixed amount of water and NH3-air mixture is placed in a closed container and maintained at constant temperature and pressure. Since NH3is very soluble in water, some NH3 molecules will instantly transfer from the gas into the liquid, crossing the interfacial surface separating the 2 phases. A portion of the NH3 escapes back into the gas, at a rate proportional to their concentration in the liquid. As more NH3 enters the liquid, the concentration increases in the liquid and the rate at which NH3 returns to the gas increases, until eventually the rate at which it enters the liquid exactly equals that at which it leaves. At the same time, through the mechanism of diffusion, the concentration throughout each phase becomes uniform. A dynamic equilibrium now exists, and while NH3 molecules continue to transfer back and forth from phase to the other, the net transfer falls to zero. The concentration within each phases no longer change. To the observer who cannot see the individual molecules the diffusion has apparently stopped.
Equilibrium relations
There are many useful equilibrium equations in books on physical chemistry. A few simple equations relating to vapor-liquid equilibrium and gas absorption: Daltons law According to this law the partial pressure of a gas within a gas mixture is equal to the product of the total pressure (P) and its mole fraction (yi) of the gas phase. This law applies whether the gas is ideal or nonideal
104
Roults law When the liquid phase can be considered ideal, it is possible to predict the equilibrium partial pressure (pi*)of an ideal gas/vapor from the solution concentration and the physical properties of the pure components.
In reality there are no ideal solutions. Ideality would require that the mplcules of the constituents be similar in size, structure and chemical nature. Practically however, many solutions are so nearly ideal that for engineering purposes they can be so considered. Roults law matches experimental data best for components present in high concentration.
Henrys law
For liquid solutions, which are not ideal, Roults law gives highly incorrect results. If the solution is sufficiently diluted, th eproportionality constant is no longer the vapor pressure of the pure componenet Pi but Hi
105
When Henrys law holds, the solubility is defined by giving the Henry law constant and the temperature. Most gases follow Henrys law to equilibrium pressure to ~5 std atm.
The bulk concentrations cs and cm are notequilibrium values, since diffusion of the solute A occurs. But these bulk concentrations cannot be used directly with a mass transfer coefficient since the 2 concentrations are differently related to the chemical potential, which is real driving force. To get around the problem, Lewis and Whitman (1923) assumed:
106
The only diffusional resistances are those residing in the fluids themselves. There is no resistance to solute transfer across the interface separating the phases. This implies that The concentrations cs,I and cm,I are equilibrium values==>
(Experiments verify the validity of the assumption) The concentrations are shown graphically in the next figure. Point P: represents the 2 bulk-phase concentrations Point M: represents the concentrations at the interface
107
For steady-state mass transfer, the rate at which A reaches the interface from Liquid 1 must equal that at which it diffuses to the bulk liquid II , so that no accumulation or depletion of A at the interface occur
The overall resistance to mass transfer is the sum of individual resistances (hence tworesistance theory) The films in the model need not be stagnant layers of a certain thickness in order for the 2-film theory to apply. Mass transfer in either film may be by diffusion through a laminar boundary layer or by unsteady-state diffusion, as in the penetration theory.
108
6. MASS TRANSFER WITH LAMINAR AND TURBULENT FLOW. 6.1 Mass transfer with laminar flow
In laminar flow molecular diffusion prevails and therelationships developed for steady-state diffusion of A through stagnant B can be used to calculate mass transfer rates. Laminar flow occurs in only a few systems, because most mass transfer systems contain more than one fluid phase, so that stable laminar boundary layers cannot build up. However, laminar boundary layer exists when mass transfer occurs between a solid and a fluid and in a falling liquid film.
Dimensionless groups
109
110