0% found this document useful (0 votes)
422 views71 pages

GATE Material. Heat Transfer. Chemical

1. Heat transfer occurs through three modes: conduction, convection, and radiation. Conduction involves the transfer of heat through direct contact of materials. Convection uses fluid movement, and radiation transfers heat without a medium. 2. Fourier's law of heat conduction states that the rate of heat transfer through a material is proportional to the temperature gradient and the area, and inversely proportional to the material thickness. 3. Thermal conductivity is a measure of a material's ability to conduct heat. It depends on factors like the material, density, temperature, and pressure. Materials with higher conductivity are used as heat conductors, while lower conductivity materials are used as insulators.

Uploaded by

suresh
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
0% found this document useful (0 votes)
422 views71 pages

GATE Material. Heat Transfer. Chemical

1. Heat transfer occurs through three modes: conduction, convection, and radiation. Conduction involves the transfer of heat through direct contact of materials. Convection uses fluid movement, and radiation transfers heat without a medium. 2. Fourier's law of heat conduction states that the rate of heat transfer through a material is proportional to the temperature gradient and the area, and inversely proportional to the material thickness. 3. Thermal conductivity is a measure of a material's ability to conduct heat. It depends on factors like the material, density, temperature, and pressure. Materials with higher conductivity are used as heat conductors, while lower conductivity materials are used as insulators.

Uploaded by

suresh
Copyright
© © All Rights Reserved
We take content rights seriously. If you suspect this is your content, claim it here.
Available Formats
Download as PDF, TXT or read online on Scribd
You are on page 1/ 71

Department of Chemical Engineering

GATE Material

Heat transfer

Nature of Heat Flow:

Heat flows from point of high temperature to those of relatively low temperature. The
net flowing is always in the direction of the temperature decreases. Whenever there is an
exchange of heat, heat lost by the hot body is always equal to heat gained by the cold body. This
substance expands on heating.

Three Modes of Heat Transfer


There are three modes of heat transfer: conduction, convection, and radiation. Any energy
exchange between bodies occurs through one of these modes or a combination of them.
Conduction is the transfer of heat through solids or stationery fluids. Convection uses the
movement of fluids to transfer heat. Radiation does not require a medium for transferring heat;
this mode uses the electromagnetic radiation emitted by an object for exchanging heat.
Conduction:

Conduction is at transfer through solids or stationery fluids. When you touch a hot
object, the heat you feel is transferred through your skin by conduction. Two mechanisms
explain how heat is transferred by conduction: lattice vibration and particle collision. Conduction
through solids occurs by a combination of the two mechanisms; heat is conducted through
stationery fluids primarily by molecular collisions. When molecules collide against one another
without moving from a fixed position, such as in a solid, the mechanism of transmitting energy is
called conduction. In solids, the heat is conducted by the following two mechanisms.

I. by Lattice Vibrations:

The faster moving molecules or atoms in hottest part of a body transfer heat by
impacts some of their energy to adjacent molecules

II. by Transport of Free Electrons:


Free electrons provide an energy flux in the direction of decreasing temperature. In
case of gases the mechanism of heat conduction by the kinetic energy of molecule is a function
of temperature.

In case of liquids the mechanism of heat the molecules of more closely spaced
and intra molecular forces come in to play.

Note1: Any rate of flow α (potential (driving) force/resistance).

Flow α (conductance *driving force)

Note2:

Generally the resistance offered by any material is inversely proportional to the conductance.

Resistance α (1/conductance).

Fourier’s law of heat conduction:

The rate of heat flow through a simple homogeneous solid is directly proportional to
area of the cross section at right angles to the direction of heat flow and to change of temperature
with respective to the length of the path of the heat flow.

Mathematically, it can be represented by the eq.

Q α A*(dx/dx)

Where Q is rate of heat flow through a body per unit time –watts;

A is surface area of heat flow –m2.

dT is temperature difference of the bases of thickness dx.

Q=-kA (dT/dx)

Where k is thermal conductivity-watt\m°k


The negative sign is to take care of the decreasing temperature along with the direction of
increasing thickness or the direction of heat flow. The temperature gradient dT/dx always
negative along positive direction and therefore the value of Q becomes positive.

Essential Features for Fourier law:

1. It is applicable to all matter (solid, liquid, gas)

2. It is a vector expression indicating that heat flow rate is in the direction of decreasing

temperature &is normal to an isothermal.

3. It helps to define thermal conductivity.

Thermal Conductivity:

Thermal conductivity Is the quantity of heat passing through a material of unit thickness with
unit heat flow are in unit time when unit temperature difference is maintained across the opposite
faces of material. Units: w/mºc

1. The thermal conductivity is a measure of ability of a substance to conduct heat.

2. Larger the value of k higher will be the amount of heat conducted.

K is the property of the material, depends essential upon the following factors

A. type of material or material structure,

B. moisture content,

C. pressure &temperature,

D. density of material.

A. Type of material:

Thermal conductivity higher for solids compared to liquids and gases.

K S>L>G
1. Thermal conductivity is higher for pure metals compared alloys &nonmetals.

2. Range of thermal conductivity for metal 232 to 420 watt/mºk.

3. Range of K for liquid 0.09 to 0.07 watt/mºk.

4. Range of K for gasses 0.006 to 0.6 watt/mºk.

5. Materials with higher thermal conductivity are usually used as heat conductors.

6. Materials with low thermal conductivity (k<0.21w/mºk) are usually used as insulators.

7. Thermal conductivity is decreases with increasing impurities of the materials.

Density & Temperature:

Thermal conductivity at low temperature is not depends on the material at


moderate &high temperatures.

1) Thermal conductivity of metals decreases with increasing Temperature (Due to decrease in


density and also alluminium&uranium being exception).

2) Thermal conductivity with respective to the temperature for most metals is almost linear.

K=K (1+βT)

Where K is thermal conductivity at 0ºc

β is coefficient thermal conductivity /ºc,it is used for positive for nonmetals, negative for metals.

3) In most of liquid the value of thermal conductivity tends to decreases with temperature (water
being an exception) due to decreases in density with increases in temperature.

4) In case of gases the value of thermal conductivity increases with temperature. Gases with
higher molecules weights have smaller thermal conductivity than with lower molecular weights.

5) This is because the mean molecular path of gas molecule decreases the increasing the
increasing in density and thermal conductivity is directly proportional to the mean prepath of the
molecule.In case of non-metallic solids thermal conductivity is increases density.
Pressure:

In case of solids and liquids thermal conductivity is weekly dependent on pressure. In case
of gases the value of thermal conductivity independent of pressure. Above 10 bars of thermal the
thermal conductivity of gases is affected on pressure.

Note1: Sutherland deduced an equation from the kinetic theory which is applicable to the
variation of the thermal conductivity of gases with temperature.

K=k32 [(492+ck)/ (T+ck)] [T/492]3/2

Where K32 is the thermal conductivity of gases at 32ºF. Ck is Sutherland constant.

T absolute temperature of the gases.

Note2: Generally, therml conductivity of gases is directly proportional to the temperature is


given by

K α T (1/2).

Note3: In case of gases the low molecular weight of molecule having higher thermal
conductivity than with molecular weight of gas molecule.

Note4: Generally, the freezing materials (Ice) having higher thermal conductivity than with l
Heat Conduction through Plane Wall:

Consider that the wall of made a material of thermal conductivity k is of uniform thickness L&
constant cross sectional area A.

Assume that k is independent of temperature and the area of wall is very large in comparison
with the thickness so that the heat losses from the edges are negligible. Hot face is at temperature
T1 and cold face at temperature T2 &both are isothermal surfaces. The direction of heat flow is
perpendicular to the wall & temperature varies in the direction of X-axis.

Q=kA (∆T)/L
Q=∆T/R

Where R is resistance offered by plane wall. R=L/kA, °c/watt.

The thermal conductivity for the plane wall is given by

S=kA/L, watt/ºc.

Heat Flux:

The amount of heat transferred per cross sectional area. Units are watt/m2.

Conduction through composite walls :( compound resistances are in series):

Consider the transmissions of heat through a composite wall consisting of a no of materials A,


B, C. and cross sectional area is same for all materials.

Thermal resistance of A is given by, RA=LA/kA;

Similarly, RB=LB/kB; RC=LC/kC;

Since, the quantity of heat transmitted per unit time through each layer is same.

Q=kAA(T1-T2)/LA =kBA(T2-T3)/LB =kCA(T3-T4)/LC.

Where T2 , T3 are temperatures of interfaces 2, 3.

From the above eqn, the overall temperature drop is ∆T=∆T1+∆T2+∆T3.


T1-T4=QLA/kAA+QLB/kBA+QLC/kCA.

=Q[LA/kA+LB/kB+LC/kC].

=Q[RA+RB+RC]

Q= ∆T/(RA+RB+RC).

Q=∆T/R

Where R=RA+RB+RC

From the above expression,

∆T/R=∆T1/RA=∆T2/RB=∆T3/RC.

We can calculate the interface temperatures by using the above equation.

Note: general equation for heat flow through conduction is given by

(Q/A)Heat Flux= (T1-T(n+1))/∑i =1n Li/ki ( For Unit area).

Mean Area of Heat Transfer:

Since, the cross sectional area of heat transfer may wary along the length of conduction
path, it is necessary to find the mean area for heat transfer rate calculations. If A1, A2 (or) r1, r2
are the areas &radiuses of the inside &outside of the surfaces.

• For a solid of constant cross section, example of large plat plate (A2/A1) <2
The mean area is given by

Am = (A1+A2)/2;

• For a solid whose cross sectional area of heat transfer is proportional to radius,
example of long hallow cylinder, The mean area is given by,

Am=A2-A1/ ln (A2/A1); (or)

The mean radius is given by, Rm = (r2-r1)/ln (r2/r1);


• For a solid whose cross sectional area of heat transfer proportional to the square of radius,
example of hallow sphere,
The mean area is given by, Am=√ (A1*A2); (or)

Rm=4πr1r2; Area of sphere is 4πr1r2 (or) 4πr2;

Heat conduction through cylinders

Case1: Uniform thermal conductivity:

Consider a hallow cylinder made of material having constant thermal conductivity .Let r1,r2
are inner and outer radiuses .T1,T2 both are temperatures of inner and outer surface’s is a
constant thermal conductivity. The temperature distribution for a hallow cylinder is given by,

T-T1/ T2-T1=ln(r/r1)/ln(r2/r1)

From the above eqn:

1) The temperature distribution is logarithmic, not linear as in the case of plain wall.

2) Temperature at any point in the cylinder can be expressed as a function of radius only.

Alternative Method:

Consider an element at radius r and thickness dr for length of hallow cylinder through
which the heat is transmitted. Let dT be the temperature drop over the element. Area through
which heat is transmitted is given by

A=2πrL.

Q = 2πkL(T1 − T2 ) / ln r2 / r1

Q=∆T/R;

Where R is resistance offered by the cylinder is given by

R=ln (r2/r1)/2πkL.
Case2:
General Eqn.for conduction heat transfer through hallow cylinders (variable thermal
conductivity)

i =n
Q = 2πL(T1 − Tn +1 ) / ∑ 1 / k i ln(r2 / r1 )
i =1

Where n are no of insulated layers.

Heat conduction through hallow spheres:

Let r1, r2 inner and outer radiuses of sphere.

T1, T2 are temperature of inner and outer surface temperatures, k is thermal conductivity of
sphere. The temperature distribution is given by

(T-T1)/ (T2-T1)= (r2/r)*((r-r1)/ (r2-r1)).

Note: From the above eqn.,it is evident that ,the temperature distribution is associated

with radial conduction through a sphere is represented by a hyperbola.

Alternative Method: Consider a small element of thickness dr at any radius R.

Area through which the heat is transferred is A=4πR2 ,

Q = 4πkr1 r2 ∆T /( r2 − r1 )

From the above Q = ∆T/R R→ Resistance offered by sphere is given by

R=(r2-r1)/4πkr1r2.
Note1: General eqn.for conduction through a composite sphere is given by

i=n
Q = 4π (T1-Tn+1)/ ∑
i =1
(rn+1-rn)/knrnrn+1)

Thermal contact resistance:

In a composite wall the calculations of heat flow are made on the assumptions,
1) The contact between the adjacent layers is perfect.

2) At the interface there is no fall of temperature.

3) At the interface, the temperature is continuous, all though there is discontinuity in temperature
gradient.

In real systems , however due to surface roughness and ,void spaces ,usually, filled with air,
the contact surfaces touches only at discrete locations .Thus, there is not a single plane of contact
which mean that ,the area available for the flow of heat at the interface will be available is less
compared to geometric space area. Due to this reduced area and presence of air voids, the large
resistance to heat flow at the interface occur.

This resistance is known as thermal contact resistance, and it causes the temperature drop
between two materials at the interfaces.

The contact resistance of AB is given by

Rth(AB )= (T2-T3)/(Q/A)

Rth(BC) = (T4-T5)/(Q/A)

Unsteady State Conduction


Transient State is a state of non
non-equilibrium,
equilibrium, when the temperatures are still changing with
time. Some simple examples are the pre
pre-heating
heating time of ovens, when you are waiting for the
temperature inside the oven to reach a specified temperature. However, even in a heated oven,
cooking occurs by transient conduction because the heat is conducted through the food as it
cooks. The internal temperature of the food, therefore, is usually increasing during cooking.

Lumped Body Model The simplest case of transient conduction is the lumped body model.
model This
model considers the heat transfer between a solid and a surrounding fluid. The main assumption
in the model is that the solid can be approximated to have a uniform temperature; that is, there
are no temperature gradients within the solid. Furthermore, the surrounding fluid must be large
enough so that its temperature remains a constant throughout the process. This approximation

can be made when the following condition is satisfied:

where h is the convection coefficient, k is the conductivity of the solid, and Lc is the
characteristic length of the solid, defined as the volume/area ratio. For a wall of thickness 2L,
2 Lc
= L; for a cylinder, Lc = r/2;
/2; for a sphere, it is r/3. The number Bi is a dimensionless parameter
called the Biot number.

For a lumped body, the temperature of the solid which is introduced into the fluid at time t = 0 is
given by:

where A is the surface area of the solid, m is the mass of the solid, cp is the specific heat of the
solid, Ti is the initial temperature of the solid, and T∞ is the temperature of the fluid. This
equation assumes that T∞ remains constant through the process. The newly introduced variable θ/
θi is the "normalized temperature," which measures the temperature in rel
relative
ative terms. That is, the
normalized temperature measures the percent change in temperature rather than exact change in
temperature in °C. Once the normalized temperature is known, the actual temperature of the solid
can be calculated based on initial temp
temperature and the surrounding temperature.

Semi-Infinite Solid
The second approximation often applied to analyze transient conduction is the semi-infinite
approximation. Semi-infinite solids can be visualized as very thick walls with one side exposed
to some fluid. The other side, since the wall is very thick, remains unaffected by the fluid
temperature. This condition is expressed as T(∞, t) = Ti, where Ti is the initial wall temperature.
This is illustrated in Figure. The condition at the exposed side of the wall is called the boundary
condition.

Figure Semi-infinite wall

Constant temperature boundary condition

One possible condition for the wall surface is a constant temperature. In this case, the
temperature inside the wall at time t, at distance x from the surface, is given by:

where Tsurface is the constant wall temperature and α is the thermal diffusivity of the wall.
Thermal diffusivity, is related to the thermal conductivity, and is the measure of how quickly
heat is dissipated in a material. The function erf is called the Gaussian Error Function, and
values for erf(x) are often tabulated or available in graphical form for convenience.
Surface convection boundary condition

Another possible case is to have convection off the surface of the wall (Fig. 4.1b). In this case, h
is assumed to be constant at all times. The temperature in the wall at time t and distance x is:

where erfc is called the Complementary Error Function, equal to 1 - erf(x).

Overall Heat transfer coefficient


Overall Heat transfer coefficient: factors influencing
1. The flow rate
2. Properties of fluid
3. Thickness of material
4. Surface condition of tubes
5. Geometrical configuration of heat exchanger
In case of boiling & condensing process the rate of heat transfer coefficient are high.

• For tubular heat exchangers we must take into account the conduction resistance in the
wall and convection resistances of the fluids at the inner and outer tube surfaces. In
general the overall heat transfer coefficient is given by
1 1 1
= + Rcond +
UA hi Ai ho Ao
ln( Do / Di )
Rcond = 1 1 1
2πkL = =
or UA U i Ai U o Ao

Where,
Ai = πDi L
inner tube surface Ao = πDo L

outer tube surface

Ro=1/hoAo=Out side tube resistance, Ri=1/hiAi=Inside tube resistance

Or

LMTD
Q& = UA(LMTD ) = n
∑ Ri i =1

1 n
1 ln (ro / ri ) 1
= ∑ Ri = + R f ,i + + R f ,o +
UA i =1 2πri hi L 2πk p L 2πro ho L
Individual Overall Heat Transfer coefficients:
1. Based on inside:

1 1 1 1 R"f ,i R "f ,o 1
= = = + + Rcond + +
UA U i Ai U o Ao hi Ai Ai Ao ho Ao

Where, R”fi =Inside fouling resistance, R”fo=Outside fouling resistance of inner tube

2. Based on Outside:

1
Uo = "
Ao AR f ,i 1
+ o + Ao Rcond + R"f ,o +
hi Ai Ai ho
Note:

Uo=Ui=U for thin walled surfaces(Di=Do=D)

1⁄U=1\hi+Xw\k+1\ho ,xw\k represents ratio of plane wall to thermal conductivity


for metal wall resistance is negligible
1\U=1\hi+1\ho ie, U= hiho\hi+ho

Based on outer side surface diameter area, the individual temp drops in two fluids and
metal wall are proportional to the individual resistance
∆T\ (1\Uo) = ∆To\ (1\ho)

Based on the inner surface area


∆T /(1 / U o ) = ∆To /(ro / ro ho )

Fouling resistance

The phenomenon of rust formation & deposited of fluid properties in the tubes of a

heat exchanger, during its normal operation is called fouling.

In a heat exchanger, during the normal operation the tube surface gets covered by deposits

of ash, soot, dirt & scale Due to these surface deposits, the thermal resistance increased

eventually, the performance of heat exchangers lowers.

Since, it is difficult acertile the thickness & thermal conductivity of the scale deposits, the

effect of scale on heat flow is considerable by specifying an equivalent scale heat transfer

coefficient (hs).

• Heat exchanger surfaces are subject to fouling by fluid impurities, rust formation, or
other reactions between the fluid and the wall material. The subsequent deposition of a
film or scale on the surface can greatly increase the resistance to heat transfer between
the fluids.
• An additional thermal resistance can be introduced: The Fouling factor, Rf.
Depends on operating temperature, fluid velocity and length of service of heat
exchanger. It is variable during heat exchanger operation.
If his= heat transfer coefficient of the scale inside surface of the inner tube

hso= heat transfer coefficient of the scale outside surface of the outer tube

1
Rsi= Inside scaling resistance; Rsi =
Ai hsi

1
Rso =Out side scale or fouling resistance of inner tube; Rso =
Ao hso

Total scale resistance is Rf=Rsi+Rso

1. Overall heat transfer coefficient in presence of onside &outside scale resistance


based on inner surface
1
Ui = =Ud
1 1 ri r r
+ + + i + i
hi hsi ro ro ho ro hso
k ln
ri

2. Overall heat transfer coefficient in presence of onside &outside scale resistance


based on outside surface
1
Uo = =Ud
1 1 ro ro ro ro
+ + ln + +
ho hso k ri ri hi ri hsi

Note

Generally the dirt resistance inversely proportional to overall dirt heat transfer coefficient;
1
Rf =
Ud

Uc −Ud
Rf =
U cU d
Where Uc=Clean overall heat transfer coefficient is given by

1 1
Uc = =
1 ri ro r ro r r 1
+ ln + i + o ln o +
hi k ri ro ho ri hi k ri ho

Forced Convection

Forced convection refers to convective transport in flows that are externally driven. The flow
could be external to a surface (flow over a flat plate, flow past a solid object or across a pipe) or
internal (flow between two walls or through a pipe). Generally fluid flows become turbulent as a
result of flow instability, and viscosity resists random motion and therefore tends to make the
flow laminar. Thus the Reynolds number is often an indicator of whether a flow is laminar or
turbulent. Since the length especially and to some extent the velocity are defined in different
ways in different flow situations, the value of Re (unlike that of certain non-dimensional groups
such as Prandtl number, Pr = cµ / k , which are properties of the fluid) cannot be a criterion of
comparison for flows in different flow arrangements.

The Reynolds number (Re) is a fair criterion of transition of a boundary layer from laminar to
turbulent, although the exact condition under which transition occurs depends on boundary
conditions such as surface roughness and the level of free stream turbulence.

Empirical formulae express the Nusselt number, Nu = hL/k (where h is the heat transfer
coefficient, L a characteristic length and k the conductivity of the fluid. The thermal conductivity
of the fluid on the heat transfer coefficient is also contained in the Prandtl number, Pr = ρc /k
(which is the ratio of Kinematic viscosity µ / ρ to the thermal diffusivity α = k/ ρc ). In some
instances, formulae are given for the Stanton number St = h / ρcV (=Nu /Re Pr). Stanton number

is the heat transfer counterpart of the skin friction coefficient (really C f / 2 = τ w / ρV 2 where V is

the characteristic velocity, the free-stream velocity on flow over a plate and bulk mean velocity
in pipe flow). The Peclet number, Pe = Re Pr (= VLk / cµ ) is also used in empirical formulae,
especially in liquid metal heat transfer.

One needs to check whether a formula is for the local value of heat transfer coefficient or the
overall value. For flow over a flat plate, the formulae are often for the local value. In pipe flows,
commonly available formulae are for fully developed flows and where the entry region is short
compared to the entire pipe, thee effect of he entry region can be ignored. In other cases, it needs
to be taken into account.

Prandtl No (Pr)
The dimensionless number associated with momentum and thermal diffusion is called the
“Prandlt No”. Pr has been found to be the parameter which relates the relative thickness of the
velocity and the thermal boundary layer.

µ Dynamic Viscosity (kg/m s)

µ CP Specific Heat (kJ / kg


ViscousDiffusivity υ ρ C µ
Pr = = = = P K)
ThermalDif fusivity α K K
ρC P
K Thermal Conductivity (W / m K)

α Thermal Diffusivity

ρ Density (kg / m3)


Heat transfer by forced convection in laminar flow

Dimensional analysis: Dimensional analysis is a mathematical technique which makes use of


the study of the dimensions for solving several engineering problems. Each physical
phenomenon can be expressed by an equation giving relat
relationship
ionship between different quantities;
such quantities are dimensional and non
non-dimensional.

Applications: 1.To tests the dimensional homogeneity of any equation of fluid motion. 2.
To derive rational formula for a flow phenomenon 3. To derive equatio
equations
ns expressed in terms of
non-dimensional
dimensional parameters to show the relative significance of each parameter.

Dimensional analysis applied to Forced convection heat transfer: Let us assume that the heat
transfer coefficient is in a fully developed forced convec
convection
tion in a tube is a function of the
following variables.

h = f (ρ , d , v, µ , C p , k )

= f 1 (h, ρ , d , v, µ , C p , k )

Note: Generally heat transfer through forced convection the Nusselt num
number
ber is function of Re
and Pr numbers.
The heat transfer coefficient depends on fluid properties (heat capacity, viscosity, thermal
conductivity, and density) and flow properties (pipe diameter, velocity). Dimensional
analysis shows the relation between the variables:

ions that follow are limited to conduit flow without phase change.
The correlations change Different
geometries, boiling, and condensation will be covered in later lectures. Frictional heating
(viscous dissipation) is not included in these correlations. This should not be a problem,
probl since
this phenomena is typically neglected except for highly viscous flows or gases at high mach
numbers.
Unless otherwise specified, fluid properties should be evaluated at the "bulk average"
temperature -- the arithmetic mean of the inlet and outle
outlet temperatures:

Significance of Nusselt number: It is the ratio of heat flow rate by convection process under a
unit temperature gradient to the heat flow rate by conduction process under a unit
temperature gradient through a stationary thickness of L meters.
Case-I: Laminar flow heat transfer to flat plate:

0.332 1 1
For Laminar Flow: Nu = Re 2 Pr 3
xo 34
3 1− ( )
x

Where, xo = distance from the leading edge

X or L =Thickness of the plate. or

The average Nusselt number is given

Where, the average heat transfer coefficient in forced convection is twice that of local heat
transfer coefficient is given by

h = 2hx

Note: For laminar flow in flat plates the Reynolds number range is

Re ≤ 3x10 5

Case: II-Heat transfer by forced convection in tubes or pipes: For laminar flow the Seider-
Tate equation for the calculation of the heat transfer coefficient for the laminar flow of fluids
inside horizontal tubes or pipes is
1 0.14
 d 3 µ 
N u = 1.86  Re Pr   
 L   µw 

hi d i d vρ
Where, N u = , Re = i
k µ

 d
The equation is valid for Re ≤ 2100and  Re Pr  ≥ 100
 L
All the fluid properties are evaluated at the arithmetic average bulk temperature of the
fluid, except µ w which is calculated at the average wall temperature.
From the above equation is related to h & G is

1 1

hi ∝ v orhi ∝ G
3 3

Graetz Number (NGz): It is related only for the heat flow to the fluid flowing through circular
pipe. It is defined as the ratio of heat capacity of fluid flowing through the pipe per unit
length of the pipe to the conductivity of the pipe.

π
Where, m is the mass flow rate is m = ρvA = ρvd 2
4

Therefore the above equation is

π d
Gz = Re Pr
4 L

This value substitute in main equation, the equation is modified as


1 0.14
4 3 µ 
N u = 1.86  G z    , For Gz>100,
π  µ
Two correlations are provided for laminar flow, depending on the magnitude of the
Graetz number. For Gz<100
1
From the above equation, hi = m 3
If the mass flow rate (m) and mass velocity (G) is given, then the Reynolds number is
dG 4m
Re = =
µ πdµ

Peclet number (Pe): The Peclet number is product of Reynolds number and
Prandtl number. Pe = Re Pr

π
Therefore Graetz number is given by, Gz=
Fo

Fully developed flow: A velocity distribution, fluid flow is unchanged with along the
length of the pipe is called fully developed flow.

If the flow is fully developed flow and constant wall temperature (Tw), the heat
transfer coefficient can be calculated by

mC p  T − Ta 
Nu = ln  w
πkL  Tw − Tb 

For constant wall temperature and fully developed flow in laminar flow through
tubes, the equation is

N u = 3.66

For uniform heat flux or constant flux through the entire length of tube, the equation is
given by

N u = 4.36

Analogies between Heat & Momentum transfer

υ ε
1. Reynolds analogy: = pr = M
α εH
By using above conditions in turbulent flow and smooth tube is given by

N u = 0.023( Re ) 0.8 [Where Pr=1]


To divide above equation by product of Re and Pr is given by

Nu − 02
= 0.023Re
Re xPr

−0.2 Nu
S t = 0.023Re [Where S t = ]
Re xPr

f h f
Therefore, S t = or =
2 C pG 2

For smooth tube turbulent flow conditions, the friction factor is given by
−0.2
f = 0.046Re

Note: Generally this analogy used for gases or thermal boundary is equal to hydrodynamic
boundary layer.

2. Chilton-Colburn analogy: For turbulent flow smooth tubes, the equation is given by
0. 8
1 µ 0.14
N u = 0.023Re Pr 3 ( )
µw

From the above equation rearranged is given by

f
JH =
2

Where, JH is called Colburn-JH factor is given by

−0.14
2  µ 
JH = S t Pr 
3 
 µw 

Where µ w is viscosity at wall temperature

3. Prandtl analogy: This is a modification over the Reynolds analogy and the presence of the

Laminar sub layer in the path of transfer is considered.

The equation is given by


f
St = 2
f
1+ (Pr − 1)
2

Note: Prandtl analogy is more accurate analogy compared to other two analogies.

Limitations of analogy equations: All analogy equations connecting f & h have an important
limitation. They apply only to wall or skin friction and must not be used for situations in which form
drag appears.

Heat transfer in transition region:

For this no simple equation is applicable and therefore a graphical method is used.The graph
between JH Vs Re with lines of constant values of L/d.
Heat transfer to liquid metals
1. Liquid metals are used for high temperature heat transfer, especially in nucler reactors.

2. Liquid mercury, sodium and a mixture of sodium and potassium called Nak are

Commonly used as carries of sensible heat.

3. Mercury vapor is also used as a carrier of latent heat.

4. Temperature of 8000C and above is obtained by using such metals.

5. Molten metals have good specific heats, Low viscosities, and high thermal

Conductivities.

The equation is given to calculate heat transfer coefficient is


0.8
 
 
 1.82 Re 
N u = 7 + 0.025 Pe −
 εM  
1.4

   
  υ  

εM
A correction for as a function of the Reynolds number is given in Fig.
υ
Generally the equation is may be written as
0.8
N u = 7 + 0.025Pe

For the critical Peclet number, the Nusselt number in case of liquid metals is given by

Nu = 7

Forced convection in spheres: For heat transfer between a flowing fluid and the surface of
single sphere, the following equation is recommended.
1 1
N u = 2 + 0.6 Re 2 Pr 3

All properties are evaluated at average film temperature Tf.


For stagnant sphere the equation is given by
Nu = 2

Natural Convection
Heat transfer by natural convection occurs when a fluid is in contact with a surface hotter or
colder than itself. As the fluid is heated or cooled it changes its density. This difference in
density causes movement in the fluid that has been heated or cooled and causes the heat transfer
to continue.

General Considerations:

• Free convection refers to fluid motion induced by buoyancy forces.

• Buoyancy forces may arise in fluids for which there are density gradients and a body force that
is proportional to density.

• In heat transfer, density gradients are due to temperature gradients and the body force is
gravitational.

• Stable(Conduction Only, No Fluid Motion) and Unstable(Free Convection Flow and Heat
Transfer)Temperature Gradients.

Flow regimes in natural convection is govern by Grashof number Gr (Like the flow regimes in
Forced convection are governed by Re). Gr represents the ratio of Buoyancy force over viscous
force acting on the fluid.
BuoyancyForces ρ 2 gβ (Tw − T∞ ) L3
Gr = Gr =
ViscousForces µ2

Empirical correlation for natural convection: The Grashof number has a similar significance
for natural convection as the Reynolds number has for forced convection, i.e. it represents a ratio
of buoyancy to viscous forces.

For natural convection process the rate of heat transfer would be expected to follow as relation:

N u = f ( N Gr N Pr )

N u = C ( N Gr N Pr )
n

The fluid properties in the dimensional groups are evaluated at the mean film temperature
Tf.
Tw + T∞
Tf =
2

Coefficient of thermal coefficient (β): It is a property of fluid. It is the fractional

increase in volume, at constant pressure, of the fluid, per degree of temperature

change.

 ∂V 
 
 ∂T  P
β= , / 0C
V

Where V is the specific volume.

• The coefficient of thermal expansion of an ideal gas is equal to the reciprocal of the
absolute temperature.

1 1
β= =
T 273 + T f

For liquids in terms of density:


ρ1 − ρ 2
β=
ρ av (T2 − T1 )

Where ρ1 and ρ 2 are density of fluid at T1 and T2 and ρ av is the average density of fluid.

ρ1 + ρ 2
ρ av =
2
Average Nusselt Number:
1
h.L 4  Gr L  4
3  4 
Nu L = = (Pr) or
k

Nu L = 4 Nu L
3

Rayleigh number:It is the product of the Groshoff’s number and Prandtl numbers.

N Ra = N Gr N Pr

In natural convection the transition from a laminar to turbulent boundary layer occurs
when N Gr N Pr =109.

Laminar conditions in free or natural convection: N Gr N Pr < 10 9 .

Turbulent condition in free or natural convection: N Gr N Pr > 10 9

Many empirical correlations are available for free convection from bodies of different geometries
and flow conditions. The average film heat transfer coefficient in natural convection is generally
given by the correlation of the form:

N Nu = C[N Gr N Pr ]
n

For isothermal surfaces the exponent’s n is ¼ for laminar flow and 1/3 for turbulent flow.
The characteristic dimension L to be used in Nusselt or Grosshoff’s number depends
upon the geometry of the problem. For vertical plate or vertical cylinder it is height of the
plate or cylinder L; for the horizontal cylinder it is the diameter D.The fluid properties
must be evaluated at film temperature.
1. Vertical plates and vertical cylinders:
Ascending flow with the maximum velocity occurring in the boundary layer and zero
velocity at both the surface and outer edge.

1
N Nu = 0.59[N Gr N Pr ]4 For 104< N Gr N Pr < 10 9

1
N Nu = 0.1[N Gr N Pr ]3 For 109< N Gr N Pr < 1012

Characteristic dimensions=Vertical dimension

2. Horizontal cylinders:
1
N Nu = 0.53[N Gr N Pr ]4 For 104< ( N Gr N Pr ) < 10 9

1
N Nu = 0.13[N Gr N Pr ]3 For 10 9 < N Gr N Pr < 1012

1
N Nu = 0.48[N Gr N Pr ]4 For10 4 < N Gr N Pr < 10 7

1
N Nu = 0.125[N Gr N Pr ]3 For10 7 < N Gr N Pr < 1012

Characteristic dimension = Diameter

3. Horizontal plates: The average heat transfer coefficient from horizontal plate depends
upon whether a plate is warmer or cooler compared to the surroundings and also on
whether the plate is facing upward or down ward.
The characteristic dimension L to be used is given by

Surfaceareaoftheplate( A)
L=
Perimeteroftheplate(P )

a. Upper surface of a heated horizontal plate/horizontal plate facing upward


/heated surface facing upward:

1
N Nu = 0.54[N Gr N Pr ]4 For 2.6x104 < N Gr N Pr <8x106

1
N Nu = 0.15[N Gr N Pr ]3 For 8x106< N Gr N Pr <1011
Above correlations are also valid for the lower surface of heated plate.

b .For calculation of natural convection heat transfer from lower surface of

heated plate /horizontal plate facing down ward/heated surface facing down

1
ward: N Nu = 0.27[N Gr N Pr ]4 For 105< N Gr G Pr <1011

Above equation is also valid for the upper surface of a cold vertical plate.

4. Spheres: The following empirical relation for free convection heat transfer from
Spheres of diameter D is:

1
N Nu = 2 + 0.43[N Gr N Pr ]4 For 1< N Gr N Pr <105

1
N Nu = 2 + 0.5[N Gr N Pr ]4 For 3x105< N Gr N Pr <8x108

N Nu = C[N Gr N Pr ]
n

For isothermal surfaces the exponent’s n is ¼ for laminar flow and 1/3 for turbulent flow.
The characteristic dimension L to be used in Nusselt or Grosshoff’s number depends
upon the geometry of the problem. For vertical plate or vertical cylinder it is height of the
plate or cylinder L; for the horizontal cylinder it is the diameter D.The fluid properties
must be evaluated at film temperature.
Heat transfer from condensing vapors

Convection processes associated with a change of phase of a fluid (at a solid-liquid interface).

Latent heat, hfg - this is an important property for phase change processes.

The primary function of a condensation or boiling is to remove latent heat or liberal latent heat
although it is sometimes necessary to remove sensible heat as well. Examples include:
• distillation columns
• reactors
• steam heaters, reboilers
• power plant boilers
• refrigeration systems
High heat transfer rates with small temperature differences & High heat transfer coefficients.
The heat transfer coefficient and rates, due to latent heat associated with phase change, are
generally much higher compared with the normal convection process(without phase change).

Properties affected on Boiling + Condensation:-The phenomenon associated with boiling &


condensations are more complex than the normal convection process due to the following
being very significant.

1. ∆T = Tsat - Ts

2. Buoyancy force arising from phase density difference, g (ρl-ρv)

3. Latent heat, hfg

4. Surface tension, σ

5. Characteristic length, L

6. Thermo physical properties of fluid, ρ , c p , k , µ .

Mechanism of condensation: Condensation of vapor to a liquid involves a change of phase of a


fluid with large heat transfer coefficients.

• Condensation occurs when a saturated vapor such as steam comes in contact with a solid
whose surface temperature is below the saturation temperature to form a liquid such as
water.
∆T = Tsat - Ts (Tsat>Ts)

Where, Tsat= Saturation temperature of vapor; Ts= Surface temperature.

Modes of condensation:

1. Surface condensation: Contact between vapor and surface at Ts < Tsat.Two types surface
condensation mechanisms are commonly observed during condensation of vapors in the
chemical industry. They are
Drop wise condensation: When a saturated vapor comes into contact with the cold solid
surface, the temperature of the solid surface is below the saturation temperature of the vapor,
some of the vapor condenses and drops are formed on the solid surface, if the condensing
liquid does not wet the surface readily. The drops grow larger till their weight overcomes the
surface tension forces holding then to the solid surface. Then the drops fall, cleaning the solid
surface for more condensation. This process is called drop wise condensation.

• Drops form in cracks, pits and cavities. This is facilitated by a non-wetting surface/liquid
combination.
• In general terms, condensate inhibits heat transfer between vapor and surface. Thus, drop
wise condensation is preferred. However, droplets not always stable so we generally
design for film condensation.
• In this type of condensation, a large portion of the area of solid surface is directly
exposed to vapor without an insulating film of condensate liquid, consequently higher
heat transfer rate achieved.
• Drop wise condensation has been observed to occur either on highly polished surfaces or
on surfaces contained with impurities like fatty acids and organic compounds.
• This type of condensation gives coefficient of heat transfer generally 5 to 10 times larger
than with film condensation.
• Although drop wise condensation would be preferred to film wise condensation yet it is
extremely difficult to achieve or maintain.
• This is because most surfaces become wetted after being exposed to condensing vapors
over a period of time.
• Drop wise condensation can be obtained under controlled conditions with the help of
certain additives to the condensate and various surface coatings but its commercial
viability has not yet been approved. For this reason the condensing equipment in use is
desired on the basis of film wise condensation.
• The maximum heat transfer coefficient in drop wise condensation is 110kW/m20K
• Liquid metals usually condense in the drop wise manner.
Note: the mixture of condensed and non-condensable vapors, the process of condensation

is called constant pressure process.


Note: condensation of pure liquid, the process is called constant pressure process.

• After condensation, the liquid flows down the tube surface under the influence of gravity
(unless vapor rates are high enough to produce vapor shear). The flow may be laminar or
turbulent, depending on the fluid, rate of condensation, tube size, etc. The film tends to
thicken as it flows to the bottom of the tube, and the weight of the fluid may cause ripples
to form. These will cause deviations from pure laminar flow.
Film type condensation: When a saturated vapor comes into contact with the cold surface it
condenses and, if condensate wets the surface, it forms a continuous film of condensate through
which heat must be transferred. The liquid flows down the cooling surface under the action of
gravity and the larger continuously grows in thickness because of newly condensing vapors. The
process of condensation is called film type condensation.
• The continuous film offers thermal resistance and checks further transfer of heat between
the vapor and the surface.
• The heat is transferred from the vapor to the condensate formed on the surface by
convection and it is further transferred from the condensate film to the cooling surface by
conduction.
• This combined mode of heat transfer by conduction and convection reduces the rate of
heat transfer(compared with drop wise condensation).That is the reason that transfer rates
of film wise condensation are lower than the drop wise condensation.
• In both cases, nucleation is typically the rate limiting step, rather than heat transfer. Most
industrial applications are based on film mechanisms, since it is tricky and expensive to
build non-wetting surfaces.
Dimensionless parameters: For Boiling and Condensation,

- We normally develop and non-dimensionless the governing equations.


- This is not possible here (two phases, complexity of process).
⇒ Use Buckingham π theorem, i.e. the heat transfer coefficient h may be depends on:

This means there are 10 variables:-

h =f ( ∆T , g ( ρl − ρ v ), h fg , σ , L, ρ , c p , µ , k )
There are 5 fundamental dimensions: m, kg, s, J, K

hL  ρ .g.( ρ − ρ )L3 c p .∆T µc p g ( ρ − ρ )L2 


= f l v l v
These turn out to be , , , 
k  µ2 h fg k σ 

• In calculations, it is common to assume the vapor


vapor-liquid
liquid interface is at thermodynamic
equilibrium at the vapor temperature. Liquid adjacent to the cold surface is assumed to be
at the surface temperature. It is also common to treat condensers as constant pressure
systems, since the total friction losses through an exchanger are usually small.
Correlations for Condensing Heat Transfer:
These notes summarize the correlations recommended for use in this class for condensing heat
transfer. All of these are for condensation of single component vapors. Choice of a correlation
depends on whether you are looking at horizontal or vertical tubes, and whether condensation is
on the inside or outside.
Preliminaries
The condensate loading on a tube is the mass flow of condensate per unit length that must be
traversed by the draining fluid. The length dimension is perpendicular to the direction the
condensate flows; the perimeter for vertical tubes, the length for horizo
horizontal
ntal tubes.

The loading may be subscripted with a b when evaluated at the bottom of a vertical tube. This
can be used to calculate a Reynolds number .

Flow is considered laminar if this Reynolds number is less than 1800.


The driving force for condensation is the temperature difference between the cold wall surface
and the bulk temperature of the saturated vapor

The viscosity and most other properties used in the condensing correlations are evaluated at the
film temperature,, a weighted mean of the cold surface (wall) temperature and the (hot) vapor
saturation temperature.

Wall Temperatures:
As we've seen before, it is often necessary to calculate the wall temperature
by an iterative approach. The summa
summarized procedure is:
1. Assume a film temperature, Tf
2. Evaluate the fluid properties (viscosity, density, etc.) at this temperature
3. Use the properties to calculate a condensing heat transfer coefficient (using the
correlations to be presented)
4. Calculate the wall temperature. The relationship will typically be something like

Although other arrangements may be more appropriate.

5. Use the wall temperature to calculate a film temperature.


6. Compare the calculated film temperature to that from the initial step. If not equal,
reevaluate the properties and repeat.
Laminar film condensation on a vertical plate (Nusselt’s equation):
The basic equations for the rate of heat transfer in film type condensation were first derived by
Nusselt.
• Bulk vapor flow is generally very slow. Gravity controlled flow of the condensate.
- Liquid film originates at the top of the plate.
- Thickness,δ, and condensate mass flow rate, m& , increase with increasing x as a
result of continuous condensation at the liquid-vapor interface, which is at Tsat.

(Ts < Tsat). Heat transfer from liquid-vapor interface → surface.

Assumptions:-
- Laminar flow and const. fluid properties in the liquid film.
- Pure vapor, at uniform temperature of Tsat ⇒ no temperature gradient at liquid-vapor
interface.
- Shear stress at liquid-vapor interface is assumed to be negligible, i.e.
∂u
y =δ =0
∂x

- Momentum + energy transfer by bulk fluid motion (advection) in the


condensate film are assumed to be negligible (very low velocities)

Boundary Conditions:-

∂u
u (o) = 0 ; y =δ = 0
∂y

Velocity Profile in Liquid Film:

g ( ρl − ρv ).δ 2  y y 
2
u( y) =  δ − 12  δ  
µl    

δ = f(x) is unknown.
Condensate Mass Flow Rate per unit width,

.
m( x ) g ⋅ ρl (ρl − ρv )δ 3
Γ( x ) = =
b 3.µl
Considering a differential element:-

Rate of heat transfer, dq, into the film at the liquid-vapor interface = rate of energy release due to
condensation at the interface.

.
⇒ dq = h f .d m
g

No convective heat transfer due to bulk fluid motion

⇒ dq = h fg dm& = qs ″ (b.dx )

.
k .(T − Ts ) h fg .d m
q"s = l sat =
δ b.dx

dΓ k .(T − Ts )
⇒ = l sat
dx δ .h fg

g .ρl .( ρl − ρv )δ 2 dδ
= ⋅
µl dx


1
 4.k .µ .(T − T ).x  4
δ ( x) =  l l sat s 
 g .ρl .( ρl − ρ v )h fg 

Convection heat transfer coefficient, hx:-


kl (Tsat − Ts )
q"s = hx [Tsat − Ts ] =
δ
k
⇒ hx = l
δ
1
 g .ρl .( ρl − ρ v ).kl3.h' fg  4
⇒ hx =  
 4 µl (Tsat − Ts ).x 


Average heat transfer coefficient, h L

h L = 1 ∫ L0 hx dx
L
1
 g .ρ .(ρ − ρ ).k 3.h ' fg  4
⇒ h L = 0.943 l l v l

 µl ..(Tsat − Ts ).L 

and

1
 ρ .g ( ρl − ρ v )h ' fg .L3  4
Nu L = 0.943 l 
 µ l .kl .(Tsat − Ts ) 

Tsat + Ts
Liquid properties at Tf = ; hfg at Tsat
2

Total rate of heat transfer:

q = h L . A.[Tsat − Ts ]

Total condensation rate:-

. h L . A.(Tsat − Ts )
m = q ' =
h fg h' fg

Figure Film condensations on a vertical plate. (a) Condensate rate for plate of

width b. (b) Flow regimes.

Note-1:For inclined surfaces the gravitational force is gsin θ is the angle between the

surface and horizontal.

1
 k L 3 gh fg ρ L (ρ L − ρ V )sin θ  4
hinclined = 1.13 
 µ L L ∆T 

1
Or hinclined = hvertical x(sin θ ) 4
Turbulent Flow outside Vertical Tubes:
When the condensate Reynolds Number is greater than 1800, the recommended

Laminar Flow Outside Horizontal Tubes:


When vapor condenses on the surface of horizontal tubes, the flow is almost always laminar. The
flow path is too short for turbulence to develop. Again, there are two forms of the same
relationship is

changes the constant in the second form from 0.725 to 0.729. The rippling condition (add 20%)
is suggested for condensate Reynolds Numbers greater than 40.

Condenser tubes are typically arranged in banks, so that the condensate which falls off one tube
will typically fall onto a tube below. The bottome tubes in a stack thus have thicker liquid films
and consequently poorer heat transfer. The correlation is adjusted by a factor for the number of
tubes, becoming for the Nth tube in the stack is

Splashing of the falling fluid further reduces heat transfer, so some authors recommend a
different adjustment
Condensation inside Tubes:
Condensate falls from tubes under the influence of gravity or by vapor shear. Outside of tubes,
there is usually enough
ough vapor space to keep velocities low enough that vapor shear is not
important; however, inside tubes the velocities are higher and more likely to effect the heat
transfer coefficient. Thus, for condensation inside tubes it is often necessary to calculate
calculat a second
value of the heat transfer coefficient (assuming vapor shear). The larger of the two values is then
used.
When vapor shear is controlling, the heat transfer coefficient can be determined from

Inside Vertical Tubes:


For condensation inside vertical
rtical tubes, use the gravity controlled "outside vertical" or the shear
controlled value of the heat transfer coefficient, whichever is larger.
Inside Horizontal Tubes:
As condensation occurs inside a horizontal tube, some liquid pools in the bottom. This
Thi reduces
the flow area and changes heat transfer. To account for this, the gravity driven condensation
equation becomes

As with vertical tubes, the designer should use the larger of the gravity or vapor shear controlled
forms.
Radiation between surfaces, radiation shielding

Heat exchange between non-black bodies


The calculation of radiation heat transfer between black surfaces is relatively easy since all the
radiant energy which strikes a surface is absorbed. The main problem is to find out the shape
factors, but once this accomplished, the calculation of heat exchange is very simple.

When non-black bodies are involved, the situation is much more complex and for all the energy
striking a surface will not be absorbed; part will be reflected back to another heat transfer
surface, and part will be reflected back to another heat transfer surface and part may be reflected
out of the system entirely.

To cope with the situation let us assume the surfaces considered in our analysis are diffusive and
uniform in temperature and the reflective and emissive properties are constant over all the
surfaces.

Ex: Two surfaces of area A1 and A2 of emissivities ε1 and ε 2 which exchange heat with each
other and nothing else would be represented by the following circuit diagram.

A2
A1
ε2
ε1

The net heat flow between two large parallel planes is given by

(
Q σ T1 − T2
=
4 4
)
A 1 1
+ −1
ε1 ε2

The net radiant heat energy between two large parallel cylinders is
Q=
( 4
A1σ T1 − T2
4

,
) Where A1 = 2πr1 L , A2 = 2πr2 L
1 A1  1 
+  − 1
ε 1 A2  ε 2 

For spheres A1 = 4πr1 2 , A2 = 4πr2


2

Double pipe and Shell & Tube HE


Heat Exchanger

A heat exchanger may be defined as equipment which transfers the energy


from a hot fluid to a cold fluid, with maximum rate and minimum investment and running
costs.

Ex:-Double heat exchanger, shell and tube heat exchanger.

These exchangers are usually assembled in effective lengths of 3.65, 4.57, 6 m.

Disadvantages of double pipe heat exchanger:

1.Smaller heat transfer surface in large floor space as compared to other type of heat
exchangers.

2. Dismantling requires large time.

3. Maximum leakage points.

Advantages:

1. The double pipe heat exchanger is very alternative, where the total heat transfer surface
required is small, 9.29 m2 to 14 m2 or less.

2. This is simple in construction cheap and easy to clean.

Shell and Tube Heat Exchangers: Introduction


To get substantial heat transfer area from a double pipe exchanger, it must be long. The result
is a high pressure drop, increased pumping costs, and large amounts of (expensive?) metal.
This means we need a more compact arrangement that still simulates countercurrent flow --
the shell and tube exchanger.
Shell and tube (a.k.a. multipass) heat exchangers are the most common industrial application
for liquid/liquid heat exchange.
Shell-Tube Heat Exchanger:-For variety of industrial services where large heat transfer
surface are required, shell –tube heat exchangers are commonly used.
Baffles: baffles (dividers) are installed to direct the flow around the tubes, increase velocity,
and promote cross flow. They also help support the tubes. The baffle cut is the ratio of the
baffle window height to the shell diameter. Typically, baffle cut is about 20 percent. It effects
both heat transfer and pressure drop. Designers also need to specify the baffle spacing;
maximum spacing depends on how much support the tubes need. These are used within the
shell to.

1. Increase the shell-side heat transfer coefficient by increasing turbulence.

2. Control and direct flow patterns of the shell-side liquid and to act as.

3. Dampe. Against vibration and.

4. Structural supports for the tubes.

Baffles are:-

1.Segmented baffles 2.Disk baffles 3.Doughnet-baffles.

Baffle arrangement in heat exchanger is shown above figure.


B—Baffle spacing:

Most common are segmental baffles as shown below.


Segmental baffles are drilld plates with heights which are generally 75% of the inside
diameter of shell. These are known as 25% cut baffles and will be used in process
industries.
Baffle spacing :-(B):-The centre-to-centre distance b/w adjacent baffles are
known as baffle spacing or baffle pitch.

In design considerations:-

The maximum Baffle spacing (B) = inside diameter of shell.


I .dofshell
Minimum baffle spacing = or 2 inches
5
Shell side and tube side passes:-

A pass is when liquid flows all the way across from one end to the other of the exchanger. We
will count shell passes and tube passes. An exchange with one shell pass and two tube passes
is a 1-2 exchanger (shown in the figure). Almost always, the tube passes will be in multiples
of two (1-2, 1-4, 2-4, etc.), since odd numbers of tube passes have more complicated
mechanical stresses, etc. An exception: 1-1 exchangers are sometimes used for vaporizers and
condensers.
With the help of passes we can change the direction of flow in the shell and tubes.

Passes are generally used to obtain higher velocities and longer paths for a fluid to
travel, with out increasing the length of the exchanger, that leads to high heat transfer
rates.
The passes on shell side are single pass, two passes.
The passes on the tube side are one, two, four, and six up to twelve
Tube or pitch arrangement:-

Tube pitch :-(PT):-The shortest centre –to-centre distance the adjacent tubes is called as tube
pitch.

Tube clearance: - (CT):- The shortest distance between two tubes is called as
clearance.

The tubes are commonly arranged.


1. Square pitch. 2. Triangular pitch.

Square pitch:-Tubes arranged such that centers of 4 tubes form a square.

Advantages: - Ease in cleaning, Lower pressure drop, Prepared when the fluid
handled is highly fouling.

Disadvantages: - The heat transfer coefficient will be low.

For the square pitch:-Normally


PT = 1.25 Х Tube outside diameter.

The clearance between tubes (CT) = 0.25 х Tube outside diameter.

The minimum clearance being 4.76mm.


Triangular pitch: - Centre of 3 tubes form an equilateral triangle.

Advantages:-More tubes in a particular shell, higher heat transfer coefficient, higher pressure
drop.

Disadvantages:-Difficult to clean, When the fluid handled is highly fouling this arrangement is not
used. If the fluids are very clean, a triangular pitch is used.

Spacers:-Spacers can be used for support of baffles.

Tie rods:-These can be used; spacers are fluted to the tie rods.

Tube sheet:-It is essentially a flat circular plate with a provision for making gasket-joint
around periphery.

A large number of holes are drilled in the tube sheet according to the pitch
requirements.
Tube sheet thickness ranges from 6mm to 25.4mm for tube outside diameter of
6mm to 4mm.
Note: if the number of shell side passes or tube side passes are increases, the heat transfer
or heat transfer coefficient, pressure drop, frictional are increases.

Note: Compare to 1-1 heat exchanger the 1-4 heat exchanger is

The velocity of 1-4 heat exchanger =4x velocity of 1-1 heat exchanger.
The heat transfer coefficient(h) of 1-4 H.E= 40.8xheat transfer coefficient 1-1 H.E
The pressure drop for 1-4 H.E=41.8xpressure drop of 1-1 H.E.
The frictional loss of 1-4 H.E= 42.8xFrictional loss of 1-1 H.E
Note: Based upon number of passes provide on tube side as well as shell side

,the shell and tube heat exchangers are referred as 1-2 shell and tube, 1-
4 shell and tube heat exchanger, 2-4 shell and tube heat exchanger.

Fluid Allocation
Usually, common sense is a reliable guide when deciding which fluid belongs on the shell side
and which on the tube side:
• Fluids which are prone to fouling belong in tubes, where the higher velocities will
reduce buildup. Mechanical cleaning is also much more practical for tubes than for
shells.
• Corrosive fluids are usually best in tubes, as are very high temperature fluids that
require alloy construction. Tubes are cheaper to fabricate from exotic materials.
• Viscous fluids require a judgment call. Placing the more viscous fluid on the shell side
will improve heat transfer, but using the tube side will lead to lower pressure drop.
• Toxic fluids should go inside the tubes to improve containment.
• Streams with low flow rates typically go on the tube side to increase the velocity and
turbulence.
• High pressure streams usually go inside the tubes, since they are less expensive to
build strong.
• Streams with low allowable pressure drop should usually be placed on the tube side.
General design of heat exchange equipment

Fluids flowing in Annuli-The equivalent diameter: When a fluid flows in a conduit having
other than circular cross-section, Such as an annulus, It is convenient to express heat transfer
coefficients and friction factors by the same type of equations and curves used for pipes and
tubes.

To permit this type of representation for annulus heat transfer it has been found advantageous
to employ an equivalent diameter De.

The equivalent diameter is four times the hydraulic radius (rH).


De = 4 xrH

The hydraulic radius is obtained as the ratio of the flow area to the welted perimeter.
Flowarea
rH =
Wettedperi meter

Flowarea
There fore, De = 4 ×
Wettedperimeter

Case -1: For heat transfer calculations, the wetted perimeter is the outer circumference of the
inner pipe with diameter d1 and De is given by

2 2
d 2 − d1
De =
d1

Case-2: In pressure drop calculations the friction not only results from the resistance of the
outer pipe but is also affected by the outer surface of the inner pipe. The total wetted perimeter
is given by π (d1 + d 2 )

Therefore, the equivalent diameter, in case of pressure drop calculations is given by

' Flowarea
De = 4 ×
Frictionalwettedperimeter

De =
' (
4π d 2 − d 1
2 2
)
4π (d1 + d 2 )

De ' = d 2 − d1

Design calculations and procedure for Double pipe heat exchanger: The prime objective
in the design of an exchanger is to determine the surface area, required for the specified duty
using the temperature differences available.

Calculations: The calculations procedure consists of computing ‘ho’and’hi’ to obtain clean


overall heat transfer coefficient (Uc).A value of design heat transfer coefficient (Ud) is
obtained by taking into account fouling resistance knowing Rd ,the surface area can be found
with the help of equation.

Q = U D A∆Tm
In the step wise procedure of calculation given below.
Step 1: Process Conditions Required:

1. Hot fluid:- Th1 , Th 2 , m h , C ph , µ , k , ∆p, Rsi , Rso

Cold Fluid: Tc1 , Tc 2 , mc , C pc , µ , k , ∆p, Rsi , Rso

2. Check the heat balance, of the heat exchanger, of the temperature is unknown by using
energy balance equation to calculate unknown temperature. The following equation is
given by
Q = mh C ph (Th1 − Th 2 ) = mc C pc (Tc 2 − Tc1 )

Parallel Flow
3. Calculate log mean temperature difference (LMTD) for co or current flow.

∆T1 − ∆T2
LMTD =
 ∆T 
ln 1  ∆T1 = Th ,i − Tc ,i
 ∆T2  ∆T2 = Th ,o − Tc ,o
1. Evaluate the physical properties of hot and cold fluid at the arithmetic mean of Th1 and
µ 0.14
Th2 and Tc1 and Tc2.For non-viscous fluids ( ) may be taken as 1 as Assumed
µw
below.
Step 2: Lanner pipe Calculations:-

1. Calculate inside flow area:


π 2
Ai = di
4

Step-3: Tube side calculations:

d i vρ diG
1. Reynolds number is given by N Re = =
µ µ
m
Where’d=Tube inner diameter, m, u=Inner flow velocity, , ρ= Density of inner
sec
fluid, kg/m3, µ=viscosity of inner fluid, kg/m.sec, G=Mass velocity of fluid, kg/m2.sec.

m
G=
Ai
Th1 + Th2 Tc1 + Tc2
2. From C p , µ , k all obtained at Thm = , Tcm = , The value of
2 2
Prandtl number is given by

Cpµ
N Pr =
k

3. Obtain the value of inside heat transfer coefficient hi using Seider-Tate or Dittus-
0.14
 µ 
Boelter equation taking   =1
 µw 
hi d i
k
= 0.023 N Re ( ) (N )0.8
Pr
n

Where, n=0.4 for heating, n=0.3 for cooling.

Step-3: Annulus or outside pipe calculations:

1. Flow area = Aa =
(
π d 2 2 − d1 2 )
4
Where, d1 = Inside diameter of inner pipe, d2 = Inside diameter of outer pipe.

Obtain equivalent diameter from following relation (based on the heat transfer) is given
by

2 2
d − d1
De = 2 ,m
d1

2. The mass velocity in annulus is given by


m
Ga = , kg / m 2 sec
Aa

De G a
3. The Reynolds number is given by N Re =
µ
4. Obtain C p , µ , k at Thm and Tcm (take mean temperature) , Prandtl number is given

by
Cpµ
N Pr =
k

5. Obtain ho either by using Dittus-Boelter or Seider-Tate equation taking

0.14
 µ 
  =1
 µw 

ho De
= 0.023( N Re ) ( N Pr )
0.8 n

5. Clean overall heat transfer coefficient is given by


1 1 r
= + o
U c ho ri hi

6. The design or dirty overall heat transfer coefficient (UD) is given by


1 1
= + Rd
UD Uc

Where Rd is the overall Dirt factor is given by Rd = R fi + R f o

7. Finally,

1
UD = or
1 r r r R f ro
+ R f o + o ln o + o + i
ho k ri hi ri ri

1
UD =
1 1 ro ro r r
+ + ln + o + o
ho hso k ri hi ri hsi ri

7. The surface area of the double pipe heat exchanger is given by


Q
A=
U D ∆Tm

A
8 .The length of heat exchanger is L=
πd o
9. For a given fouling period, the total Dirt resistance is given by

Uc −UD
Rd =
U cU D

Note: When Rd(deposited) >Rd(Allowed), as after a period of service , the Apparatus no


longer delivers a quantity of heat equal to the process requirements and must be cleaned.

Step-4: Pressure drop calculations (inner pipe): The pressure drop allowance in an
exchange is the static fluid pressure which may be expended to drive the fluid through the
exchanger. The pump selected for the circulation of a process fluid is one which develops
sufficient head at the desired capacity to overcome the frictional losses caused by connecting
piping, fitting, and control regulators.

1. The pressure drop due to friction is calculated by using fanning equation is


L v2
∆p f = 4 f
d 2g

Where, f is fanning friction factor, depends on the type of flow. For laminar flow,
smooth pipe is

16
f =
N Re

0.125
For smooth tube and turbulent flow is f = 0.0014 + 0.32
N Re

2. The pressure drop exited due to velocity head of fluid enters or leaves exchanger is
given by
v2
∆pv =
2g

3. The total pressure drop in tube side is ∆p = ∆p f + ∆p v

Note: If the Total pressure drop>Maximum allowable pressure drop, the design of
double pipe heat exchanger is not satisfied.

Step-5: Pressure drop in annulus: By using same equations, but in presence

of diameter, use equivalent diameter.

From the above equations


∆p ∝ G 2

hi ∝ G 0.8

The best use of available pressure is to increase the mass velocity which also increases
hi and less the size and cost of the apparatus.
Shell and Tube Heat Exchangers: Design Calculations
Step-1: Thermo physical properties of hot and cold fluids are given below.
Hot fluid: Th1 , Th 2 , mh , C ph , s, µ , k , Rd , ∆p

Cold fluid: Tc1 , Tc 2 , mc , C pc , s, µ , k , Rd , ∆p .

Step-2: for the shell-tube heat exchanger the following data must be known.
Shell-side: ID of shell, Baffle spacing, passes on shell side.
Tube side: BWG number, Length, Outside diameter, pitch and passes.
Step-3: Heat balance of exchanger is
Q = mh C ph (Th1 − Th 2 ) = mc C pc (Tc 2 − Tc1 )

Step-4: Temperature correction factor (F): In double pipe or 1-1 heat exchangers, when the
hot and cold fluids in heat exchanger are in true counter current flow or in parallel flow, the
log mean temperature difference should be used. In the cases where a multipass heat
exchanger is involved, the mean temperature is depending on the arrangement of the shell and
tube passes. Consider 1-2 heat exchanger, the cold fluid in the first tube pass is in counter flow
with the hot fluid. In the second tube pass the cold fluid is in parallel flow with the hot fluid.
So, in case of multi pass heat exchanger to use temperature correction factor (F) to calculate
true temperature drop. The basic design calculation for any heat exchanger is the
determination of heat transfer area. Most generally, this is done using
although in practice it is more common to assume that fluid properties can be treated as
constant at the bulk average values, and approximate the design equation with:

To use this equation, it is necessary to determine the heat transfer coefficient and the
temperature difference. We have seen that for a double pipe heat exchanger the required
average temperature difference is the log mean temperature difference (LMTD).
Unfortunately, the flow patterns in shell and tube exchangers are such that the LMTD by itself
is no longer adequate. It must first be adjusted by means of a correction factor. The second
parameter that must be calculated for a typical process design is the pressure drop in the fluids
moving through the exchanger.

The LMTD was developed for a model restricted to parallel and countercurrent flow patterns.
In shell and tube exchangers, the flow pattern is a mixture of concurrent, countercurrent, and
cross flow, so the LMTD does no
not directly apply. Instead, a corrected LMTD must be used.
Correction factors (F correction factors) are tabulated and used to adjust:

The correction factor charts are available from many sources. They aare
re based on two
parameters:
Th1 − Th 2 T − Tc1
R= And P = c 2
Tc 2 − Tc1 Th1 − Tc1
These parameters are cross-referenced on the appropriate chart to find the F factor. A different
chart is needed for each exchanger layout (1- 2, 1-4, etc.). Note that the charts provided in
MSH are written using hot and cold subscripts and assume that the hot fluid is on the shell
side of the exchanger. Be careful and do what is right for your problem.
F factor curves drop off rapidly below 0.8. Consequently, if your design is indicating
an F less than 0.8, you probably need to redesign (add tube passes, increase
temperature differences, etc.) to get a better approximation of countercurrent flow and
thus higher F values.
Similar correction factors are available for exchanger types other than shell and tube,
including cross flow and compact exchangers.
Tube Side Heat Transfer Coefficients
Tube side heat transfer coefficients are easy to determine, since the Seider-Tate equation (or
equivalent) applies. The most common slipup made at this stage is to incorrectly divide the
stream flow among the tubes. Be sure to adjust all flow rates and velocities for the number of
tubes and the number of tube passes. If an exchanger has 200 tubes in 2 passes, the total flow
will be moving through 100 tubes at a time; if there are 4 passes, it will go through 50 tubes.
The mass velocity in tube side is given by
Massflowrate int ubeside × n
Gt =
Flowareaofeachtube × N
Where, n = is the number of tube passes, N= is the number of tubes.
0.8
From the Dittus-boilter equation, hi ∝ Gt , or hi ∝ n 0.8
Step-5: Shell side calculations: The shell side heat transfer coefficient (ho) can
be calculated by using following relation
ho d e
= 0.36( N Re ) (N pr )
0.55 0.33

k
Where, de=Equivalent diameter of shell side heat transfer is depends
upon the pitch arrangement.
 2 π 2
4  PT − d o 
de =  
4
For square pitch arrangement
πd o

 2 π 2
4 0.86 PT − d o 
de =  
4
For triangular pitch
πd o

d e Gs Cpµ
Where, N Re = , N Pr =
µ k
massflowrateinshellside
Gs = Shell side mass velocity =
shellsideflowarea
ID × B × CT
Shell side flow area (As) is given by As =
PT
Where ID=inner dia of shell
B = Baffle spacing
CT =Tube clearance = PT − Outsidediaoftube
PT = Tube pitch
0.55 1
From the above equation hoαGs or ho ∝
B 0.55
Step-6: Calculation of clean overall heat transfer coefficient:
1
U o orU c =
ro 1 r r
+ + o ln o
hi ri ho k ri
Step-7: Overall design heat transfer coefficient:

1
UD =
1 r R fi ro ro ro
+ R fo + o + + ln
ho hi ri ri k ri
Step-7: Overall dirt factor or scale factor (Rd):

Uc −UD
Rd =
U cU D
Note: If Rd (deposited)>Rd (allowed), the heat exchanger cannot be operate.
Step-8: Number of tubes: In shell-tube heat exchanger, the overall rate of heat transfer is
Q = FU D Ao ∆Tlm

Where Ao = πd o nNL

Step-9: Tube side pressure drop:


Tube side pressure drop is calculated using the same pipe flow factors developed in your fluid
mechanics class. Again, be careful to divide the total flow among the correct number of tubes.
The isothermal friction factor can be obtained from the Moody/Stanton charts or an
appropriate correlation.
This friction factor must be corrected for the effect of temperature on viscosity. MSH6 pp.
115-16 explains correction factors, but you can probably get good enough results using the
Seider-Tate viscosity correction:
This can then be used with the mechanical energy balance to get the pressure drop:

The pressure drop due to friction is given by

2
fGt Ln N
∆P f = , Where, s=Specific gravity,
d e sφ t m2
φt = Tube side viscosity correction factor.
The pressure drop due to velocity gradient is given by
4n v 2
∆Pv =
s 2g
Total pressure drop in tube side is
∆PT = ∆Pf + ∆Pv
2
From the above correlations ∆P ∝ Gt Ln or ∆P ∝ n 3
Step-10: Shell side pressure drop:
The shell side flow path is considerably more complex than for a tube, so the calculation of
shell side pressure drop is significantly more complicated as well. For our purposes, we will
use a correlation presented by Kern. It is somewhat dated, but has the advantage of mirroring
the tube side calculation.Begin by determining a shell side equivalent diameter, de.
To transform the friction factor to a pressure drop, we need a count of how many times the
fluid crosses the tube bundle. It crosses between the baffles, so the "crosses" will be one more
than the number of baffles, NB. If the number of baffles isn't already known, it can be
determined using the baffle spacing:
Tubelength L
Number of crosses = N B + 1 = =
Bafflespacing B
The pressure drop is then determined using the equivalent diameter, cross flow velocity,
friction factor, number of crosses, and fluid properties:

fG s d s (N B + 1) N
2

∆Ps = ,
2 gρ d e φ s m2

Note: In shell-tube heat exchanger, there should always be an odd number of crosses (Even
number of baffles), if both shell nozzles are on opposite sides of the shell and an even number,
if both shell nozzles are on the same side of the shell.
From the above equation,
1
∆Ps ∝ Gs ( N B + 1) Or ∆Ps ∝
2

B3
NTU-Method: 1.The LMTD method may be applied to design problems for which the fluid
flow rates and inlet temperatures, as well as a desired out let temperature are prescribed. For a
specified Heat exchanger type, the required sizes, as well as the other out let temperature, are
readily determined.

2. If the LMTD is used in performance calculations for which both out let temperatures must
be determined from the knowledge of the inlet temperatures, the solution procedure is
iterative.

3. For both design and performance calculations, the effectiveness-NTU method may be used
without iteration.

4. NTU-method is used when the temperatures of the fluids leaving the exchanger are not
known.

5. This method is based on the effectiveness of a heat exchanger in transferring a given


amount of heat.

Definition of NTU: is defined as the change in temperature of one stream divided by the
average temperature difference or driving force in the exchanger.

NTU can be based on either stream, but generally the stream with the lower capacity
(mCp) is chosen.
The temperature change is always taken as a positive number, when the Cp and U are
UA
constant. NTU =
(mC p )min
Capacity ratio(R): It is the ratio of lower capacity stream to higher-capacity stream.

mC p (lowercapacitystrem )
R=
mC p (highercapasitystream )

If m h C ph > mc C pc (Hot fluid capacity > Cold fluid capacity),the capacity ratio is

mc C pc
R=
mh C ph

Heat exchanger effectiveness ( ε ): is defined as the ratio of actual heat transfer to the
maximum possible heat transfer. Thus
Actualheattransferhotorcoldfluid
ε=
Maximumpossibleheattransfer

mh C ph (Th1 − Th 2 ) Th1 − Th 2
Case-1: If m h C ph < mc C pc , ε= =
mh C ph (Th1 − Th 2 ) Th1 − Tc1

mc C pc (Tc 2 − Tc1 )
Case-2: If m h C ph > mc C pc , ε=
m h C ph (Th1 − Tc1 )

For the above equations, generally less capacity will take on numerator.
Case-3: For counter flow heat exchanger (mh C ph < mc C pc ) , the effective ness is given by

1 − Exp[− NTU (1 − R )]
ε=
1 − RExp[− NTU (1 − R )]

Case-4: For parallel flow heat exchanger (mh C ph < mc C pc ) , the effectiveness is given by

1 − Exp[− NTU (1 + R )]
ε=
1+ R

Case-5: For steam condensers and evaporators, by using above case, we arrive at the
following common expression for parallel flow as well as counter flow heat

exchangers. Where R=0.

ε = 1 − Exp(− NTU )

Case-6: For regenerators, when R=1, in case of parallel flow, the effectiveness is

1 − Exp (− 2 NTU )
ε=
2

NTU
For counter flow regenerators the effectiveness is ε =
1 + NTU

Evaporation
There are three main measures of evaporator performance:
1. Capacity (kg vaporized / time)
2. Economy (kg vaporized / kg steam input)
3. Steam Consumption (kg / hr)
Capacity: Capacity of an evaporator is defined as the number of kilograms of water
evaporated per hour.

The rate of heat transfer Q through the heating surface of evaporator is the product of area of
heat transfer surface A, the overall heat transfer coefficient U and the overall temperature drop
∆T

Q = UA∆T

Where, ∆T = Ts − TL

Ts= saturation temperature of steam, TL= Boiling point of solution.

If the feed solution is at the boiling temperature corresponding to the pressure in vapor
space of evaporator, then all the heat transferred through heating surface is available
for evaporation and capacity is proportional to rate of heat transfer.
Steam economy

Evaporator performance is rated on the basis of steam economy. Steam economy is defined as
the kilograms of the solvent evaporated per kilogram of steam used. The greatest increase in
steam economy is achieved in a multiple effect evaporator, by using the vapor from one effect
as the heating medium for the next effect, in which boiling takes place at a lower temperature
and pressure. The steam economy of a multiple effect evaporator will increase in proportion to
the number of effects.

Economy calculations are determined using enthalpy balances.


The key factor in determining the economy of an evaporator is the number of effects. The
economy of a single effect evaporator is always less than 1.0. Multiple effect evaporators have
higher economy but lower capacity than single effect. The thermal condition of the evaporator
feed has an important impact on economy and performance. If the feed is not already at its
boiling point, heat effects must be considered. If the feed is cold (below boiling) some of the
heat going into the evaporator must be used to raise the feed to boiling before evaporation can
begin; this reduces the capacity. If the feed is above the boiling point, some flash evaporation
occurs on entry.
Steam Consumption(S) = Capacity/Economy.
Multiple effect evaporators

Single effect evaporators: When a single evaporator is used, the vapor from the boiling liquid
is condensed and discarded the method is known as single-effect evaporation.

The economy of single-effect evaporator is always less than one.


Generally for evaporation of 1 kg of water from a solution, 1.3 kg of steam is required.
The simplest mode of evaporation is to use a single stage, where steam is fed into the steam
chest, concentrate and vapor are removed, and the vapor is condensed into hot water.

However, the vapors produced are still steam, and thus can be used to provide the heat for
evaporation in a subsequent stage. Therefore, steam can be used many times to provide
evaporation in a series of operations.

Evaporation by heating with steam in a single vessel is called evaporation in a single effect.
They are used when the required capacity is small, steam is cheap and are operated in batch,
semi-batch or continuous modes. This method of operation is rarely used, since it requires that
the body be large enough to hold the entire charge of feed.

Multiple effect evaporators: Most of the evaporators use low pressure steam for heating
purpose. Due to addition of heat to solution in evaporator by condensation of steam, the
solutions in the evaporator boil.

Vapors leaving the evaporator are fed to same form of condenser then heat associated
with vapor will be lost and system is said to make poor use of system.
Vapor coming out of evaporator can be used as heating media another evaporator
which will be operating at pressure lower than the pressure in the evaporator from
which vapors are issuing so as to provide sufficient temperature gradient for heat
transfer in the evaporator.
Method of increasing the evaporation per kilo gram of steam by employing series of
evaporator’s between steam supply and condenser is known as multiple – effect –
evaporation.
It is one way to increase the economy of evaporator system.
Multiple effect evaporation system is commonly used in large scale operation. In such
system, evaporators are arranged in series so that the vapor from one evaporator is
used as heat source for the next one that is operating under lower pressure and
temperature each unit in series is called as an effect.
In case of triple effect evaporator , If the first effect is operating at atmospheric
pressure , then the second and third effect operate under vacuum , steam is fed the first
effect and vapor from third effect is condensed in condenser connected to vacuum
pump.
The economy of multiple effect evaporation is grater than one.
With multiple effect evaporation system, it is possible theoretically to evaporate, ‘N’
kgs of water from 1kg of initial steam fed, where ‘N’ is the number of effects.
Note: - In multiple effect evaporation, the temperatures and pressure gradually decrease from
effect to effect.

If the vapor from the first evaporator is taken for heating in a second evaporator, we have a
double effect. The heat in the original steam is reused in the second effect. The evaporation
achieved by a unit mass of steam fed to the first effect is approximately doubled. Additional
effects can be added in the same manner.

Thus, with three vessels, we have a triple effect and so on like a quadruple, quintuple or a
sextuple effect. There are thus as many evaporator vessels as the effects. The most common
evaporators encountered in practice are the triple, quadruple and quintuple effects.

Evaporators are classified by the number of effects. In a single-effect evaporator, steam


provides energy for vaporization and the vapor product is condensed and removed from the
system. In a double-effect evaporator, the vapor product off the first effect is used to provide
energy for a second vaporization unit. Triple- effect evaporator problems are familiar to
generations of engineering students. This cascading of effects can continue for many stages.
Multiple-effect evaporators can remove much larger amounts of solvent than is possible in a
single effect.
In a multiple effect arrangement, the latent heat of the vapor product off of an effect is used to
heat the following effect. Effects are thus numbered beginning with the one heated by steam.
It will have the highest pressure. Vapor from Effect I will be used to heat Effect II, which
consequently will operate at lower pressure. This continues through the train: pressure drops
through the sequence so that the hot vapor will travel from one effect to the next.
Normally, all effects in an evaporator will be physically the same in terms of size,
construction, and heat transfer area. Unless thermal losses are significant, they will all have
the same capacity as well.
Evaporator trains may receive their feed in several different ways. The feed order is NOT
related to the numbering of effects. Effects are always numbered according to decreasing
pressure (steam flow).
Methods of increasing the economy:-

1. Use of multiple effect evaporation system

2. Vapor recompression.

Use of multiple –effect-evaporation system:- In the multiple effect evaporation system ,


vapor produced in first effect is fed to the steam first of second effect as heating medium in
which boiling takes at low pressure and temperature and so on .Thus in triple –effect
evaporator, 1kg of steam fed to first effect evaporates approximately 2.5kgs of water.

Boiling Point Elevation


Since evaporators dealing with boiling solutions, and in particular with solutions with non-
volatile solutes, any calculations must account for the effect of boiling point elevation.
The vapor pressure of an aqueous solution is less than that of pure water at the same
temperature; so the boiling point of the solution will be higher than that of the water. This is
called Boiling Point Elevation (BPE) or vapor pressure lowering.
The boiling point of a solution is a colligative property -- it depends on the concentration of
solute in the solution, but not on what the solute and solvent are.
When working problems involving heat transfer to or from boiling solutions, it is necessary to
adjust the temperature difference driving force for the boiling point elevation.
Note that the equilibrium vapor rising from a solution exhibiting boiling point elevation will
exist at a temperature and pressure such that it is superheated with respect to pure vapor. The
vapor rises at the solution boiling point, elevated with respect to the pure component boiling
point. The vapor, however, is solute free, so it won't condense until the extra heat
corresponding to the elevation is removed, thus it is superheated

Determination for strong solutions, one can take advantage of Duhring's Rule
The boiling point of a given solution is a linear function of the boiling point of water at the
same temperature.

This lets us plot TBP solution against TBP water and get a straight line for each concentration.
Another way of thinking of these plots -- they plot the temperature where the vapor pressure
of the solution is equal to some fixed value against the temperature where the vapor pressure
of water equals the same value.

For Duhring Plots to be valid the range of boiling points must be relatively narrow and the
solution must obey Raoult's Law.
A Duhring Plot for the NaOH/Water system can be found in McCabe.
To use a Duhring plot:
1. For a particular system pressure, determine the boiling temperature of pure water. This
can be done from a vapor pressure equation or steam table.
2. Enter the plot from the bottom (the water boiling point), trace up to the diagonal line
representing the NaOH fraction, and then trace left to read the solution boiling point
from the vertical axis.
3. The boiling point elevation is the difference between the two temperatures.
Liquid head ( ∆Th ): The pressure exerted by the column of boiling liquid above the heating
surface enhances the boiling point of solution, which is called the hydrostatic head effect in
boiling point elevation.

The temperature drop due to hydrostatic effect is given by


0.87v R TR
∆Th = ∆P
λs

Where, ∆Th = Hydrostatic elevation in boiling point 0C. v R = Specific volume of water, m3/kg,
TR= Solution boiling point,0K, λ s = Latent heat at saturation pressure, Kcal/kg.

1
∆P = × Liquidlevelabovetheheatingsurface, m
2
The true temperature drop in evaporator is given by
∆TTrue = ∆T App − BPE − ∆Th

Where, ∆T App = Ts − TL

BPE and Hydrostatic effect on multiple- effect evaporation: 1. for same terminal
conditions (input and out put temperature), and BPE is negligible, the capacity of single effect
is equal to the multiple effect evaporators.

2. When the BPE is negligible, the effective overall ∆T is equals the sum of the ∆T ' s

In each effect, and the amount of water evaporated per unit area of surface in an N-effect
multiple-effect evaporator is approximately (1/N)th that in the single effect.

3. If the BPE is considered, the capacity of multiple effect evaporator is lees than the single
effect evaporation..

4. The economy of a multiple-effect evaporator is not influenced by BPE. The economy of


multiple-effect evaporator depends on heat balance considerations and not on the rate of heat
transfer.

5. The capacity of the double-effect evaporator concentrating a solution with a BPE is


generally less than half the capacity of two single-effects, each operating with the same overall
temperature drop. The capacity of a triple effect is generally less than one third of three single-
effects with same terminal temperatures.

6. The liquid head and the BPE influence the capacity of multiple-effect evaporator even more
than does that of a single effect.

7. The liquid head reduces the temperature drop available in each effect of a multiple-effect
evaporator, as does the BPE.
Methods of feeding
Methods of Feeding:-Most common methods used in chemical industry for feeding the
multiple evaporation system are.

1. Forward feeding 2. Backward feed 3. Mixed feeding 4. Parallel feeding

Forward feed arrangement: - In this arrangement, the liquid feed flows in the same direction
as the vapor flows.

The thin liquid is fed to the first effect by means of a pump and flows through

The remaining effects in the same direction as the vapor.


This type of operation is used when the feed is hot and product is damaged at high
Temperature.
Since the flow is in the direction of decreasing pressure, this operation does not
require pumps for transferring concentrate material from effect to effect.
Back ward feeding:-In this arrangement the feed solution and vapor flow in opposite
direction.

Feed is introduced in the effect and the concentrated product is taken out from
the first effect.
The feed flows in the direction of increasing pressure and temperature, which is
opposite to that of the vapor flow.
Here liquid pumps are used in each effect
This type of operation is used when the feed is cold, since it saves a huge
quantity of live steam and makes for batter steam economy.
The operation is also used when the concentrated product is highly viscous the
high temperature in the early effects reduces the viscosity, thereby increasing
the heat transfer coefficients.
Mixed Feeding:- In this feed arrangement, steam is admitted to steam chest of first effect and
vapors leaving first effect evaporator are fed to the steam chest of first of second effect and so
on.

Feed solution is admitted to an intermediate effect and flows to the first effect from
where it is fed to last effect for final concentration.
This arrangement is combination of forward and backward feed adopted for best
overall performance.
Parallel Feeding:-In the parallel operation the feed is introduced directly in to each effect and
no material is transferred from effect to effect.

The concentrated liquid is withdrawn from each effect.


This type of operation is used when the feed is a saturated solution and solid crystals
are the product, as in the manufacture of common salt.
This type of operation used in crystallization.
Comparison of forward feed and back ward feed arrangements:-

1. In forward feed, the flow of solution to be concentrated is parallel to steam flow.

In backward feed the flow of the solution to be concentrated is in opposite direction


to that of steam flow.

2. Forward feed arrangement does not need pump for moving the solution from effect to
effect as vacuum is maintained on the last effect.

Back ward feed arrangement needs pumps for moving the solution from effect
to effect as transfer of solution is to be done from evaporator operating at high pressure.

3. In forward feed as heating of cold feed solution is done in first effect, less vapor is

Produced per kilograms of steam fed, resulting into lower economy

In backward feed, the solution is heated in each effect which usually

Results in better economy than that with forward flow.

4. With forward feed most concentrated liqueur is in the lost effect where in the

Temperature is lowest and viscosity is highest. These conditions lead to reduction

in capacity of the system.

In case of backward feed, the concentrated liquor is in first effect where


in the temperature is highest as steam is admitted to that effect and the viscosity is

lowest , thus the overall coefficient can be moderate high in spite of high viscosity.

5. In case of forward feed, the maintenance charging and power cost is less.

In case of backward feed, the maintenance charging and power cost is more for same
duty.

7. Forward feed is more economical in steam, at low values of feed temperature,


back ward gives higher economy.

You might also like

pFad - Phonifier reborn

Pfad - The Proxy pFad of © 2024 Garber Painting. All rights reserved.

Note: This service is not intended for secure transactions such as banking, social media, email, or purchasing. Use at your own risk. We assume no liability whatsoever for broken pages.


Alternative Proxies:

Alternative Proxy

pFad Proxy

pFad v3 Proxy

pFad v4 Proxy