r. shanthini 19 may 2010 1 content of lectures 10 to 12: heat transfer: mathematical problems on...

42
R. Shanthini 19 May 1 Content of Lectures 10 to 12: Heat transfer: Mathematical problems on heat transfer PM3125: Lectures 10 to 12

Upload: rhoda-fletcher

Post on 11-Jan-2016

220 views

Category:

Documents


3 download

TRANSCRIPT

Page 1: R. Shanthini 19 May 2010 1 Content of Lectures 10 to 12: Heat transfer: Mathematical problems on heat transfer PM3125: Lectures 10 to 12

R. Shanthini 19 May 2010

1

Content of Lectures 10 to 12:

Heat transfer: • Mathematical problems on heat transfer

PM3125: Lectures 10 to 12

Page 2: R. Shanthini 19 May 2010 1 Content of Lectures 10 to 12: Heat transfer: Mathematical problems on heat transfer PM3125: Lectures 10 to 12

R. Shanthini 19 May 2010

2

Heat Exchanger in the industry

Shell-and-Tube Heat Exchanger

Page 3: R. Shanthini 19 May 2010 1 Content of Lectures 10 to 12: Heat transfer: Mathematical problems on heat transfer PM3125: Lectures 10 to 12

R. Shanthini 19 May 2010

3

Heat Exchanger in the industry

Shell-and-Tube Heat Exchanger

Page 4: R. Shanthini 19 May 2010 1 Content of Lectures 10 to 12: Heat transfer: Mathematical problems on heat transfer PM3125: Lectures 10 to 12

R. Shanthini 19 May 2010

4

Heat Exchanger in the industry

Shell-and-Tube Heat Exchanger

Page 5: R. Shanthini 19 May 2010 1 Content of Lectures 10 to 12: Heat transfer: Mathematical problems on heat transfer PM3125: Lectures 10 to 12

R. Shanthini 19 May 2010

5

Heat Exchanger in the industry

Shell

Tubes

Baffle

Shell-and-Tube Heat Exchanger

Page 6: R. Shanthini 19 May 2010 1 Content of Lectures 10 to 12: Heat transfer: Mathematical problems on heat transfer PM3125: Lectures 10 to 12

R. Shanthini 19 May 2010

6

Heat Exchanger in the industry

Shell-and-Tube Heat Exchanger

Page 7: R. Shanthini 19 May 2010 1 Content of Lectures 10 to 12: Heat transfer: Mathematical problems on heat transfer PM3125: Lectures 10 to 12

R. Shanthini 19 May 2010

7

Heat Exchanger Analysis

Tc,in

Tc,out

Th,in

Th,out

. .Q = mc cc (Tc,out – Tc,in) = mh ch (Th,in – Th,out)

.

Page 8: R. Shanthini 19 May 2010 1 Content of Lectures 10 to 12: Heat transfer: Mathematical problems on heat transfer PM3125: Lectures 10 to 12

R. Shanthini 19 May 2010

8

Tc,in

Th,out

Th,in

Tc,out

Tc,in

Tc,out

Th,in

Th,out

Parallel-flow heat exchanger

high heat transfer

low heat transfer

Heat Exchanger Analysis

Page 9: R. Shanthini 19 May 2010 1 Content of Lectures 10 to 12: Heat transfer: Mathematical problems on heat transfer PM3125: Lectures 10 to 12

R. Shanthini 19 May 2010

9

Tc,in

Th,out

Th,in

Tc,out

a b

Q = U A ΔTLM .

where ΔTLM = ΔTa - ΔTb

ln(ΔTa / ΔTb)is the log mean temperature difference (LMTD)

ΔTa ΔTb

Parallel-flow heat exchanger

Heat Exchanger Analysis – LMTD method

Page 10: R. Shanthini 19 May 2010 1 Content of Lectures 10 to 12: Heat transfer: Mathematical problems on heat transfer PM3125: Lectures 10 to 12

R. Shanthini 19 May 2010

10

Tc,out

Tc,in

Th,in

Th,out

. .Q = mc cc (Tc,out – Tc,in) = mh ch (Th,in – Th,out)

.

Counter-flow heat exchanger

Heat Exchanger Analysis

Page 11: R. Shanthini 19 May 2010 1 Content of Lectures 10 to 12: Heat transfer: Mathematical problems on heat transfer PM3125: Lectures 10 to 12

R. Shanthini 19 May 2010

11

Tc,in

Tc,out

Th,in

Th,out

Th,in

Th,out

Counter-flow heat exchanger

Tc,in

Tc,out

Heat Exchanger Analysis

Page 12: R. Shanthini 19 May 2010 1 Content of Lectures 10 to 12: Heat transfer: Mathematical problems on heat transfer PM3125: Lectures 10 to 12

R. Shanthini 19 May 2010

12

a b

Q = U A ΔTLM .

where ΔTLM = ΔTa - ΔTb

ln(ΔTa / ΔTb)is the log mean temperature difference (LMTD)

ΔTa

ΔTb

Counter-flow heat exchanger

Tc,in

Tc,out

Th,in

Th,out

Heat Exchanger Analysis – LMTD method

Page 13: R. Shanthini 19 May 2010 1 Content of Lectures 10 to 12: Heat transfer: Mathematical problems on heat transfer PM3125: Lectures 10 to 12

R. Shanthini 19 May 2010

13

An exhaust pipe, 75 mm outside diameter, is cooled by surrounding it by an annular space containing water.

The hot gases enters the exhaust pipe at 350oC, gas flow rate being 200 kg/h, mean specific heat capacity at constant pressure 1.13 kJ/kg K, and comes out at 100oC.

Water enters from the mains at 25oC, flow rate 1400 kg/h, mean specific heat capacity 4.19 kJ/kg K.

The heat transfer coefficient for gases and water may be taken as 0.3 and 1.5 kW/m2 K and pipe thickness may be taken as negligible.

Calculate the required pipe length for (i) parallel flow, and for (ii) counter flow.

Heat Exchanger Design – LMTD method

Page 14: R. Shanthini 19 May 2010 1 Content of Lectures 10 to 12: Heat transfer: Mathematical problems on heat transfer PM3125: Lectures 10 to 12

R. Shanthini 19 May 2010

14

Solution:

. .Q = mc cc (Tc,out – Tc,in) = mh ch (Th,in – Th,out)

.

(1400 kg/hr) (4.19 kJ/kg K) (Tc,out – 25)oC

= (200 kg/hr) (1.13 kJ/kg K) (350 – 100)oC

The temperature of water at the outlet = Tc,out = 34.63oC.

Heat Exchanger Design – LMTD method

Page 15: R. Shanthini 19 May 2010 1 Content of Lectures 10 to 12: Heat transfer: Mathematical problems on heat transfer PM3125: Lectures 10 to 12

R. Shanthini 19 May 2010

15

(i) Parallel flow:

ΔTa = 350 – 25 = 325oC

ΔTb = 100 – 34.63 = 65.37oC

ΔTLM = ΔTa - ΔTb

ln(ΔTa / ΔTb)

325 – 65.37

ln(325 / 65.37)= = 162oC

Q = U A ΔTLM .

= (UA) 162oC

What is UA?

Heat Exchanger Design – LMTD method

Solution continued:

Page 16: R. Shanthini 19 May 2010 1 Content of Lectures 10 to 12: Heat transfer: Mathematical problems on heat transfer PM3125: Lectures 10 to 12

R. Shanthini 19 May 2010

16

1/U = 1/hwater + 1/hgases

= 1/1.5 + 1/0.3 = 4 (kW/m2 K)-1

Therefore, U = 0.25 kW/m2 K

A = π (outer diameter) (L) = π (0.075 m) (L m)

Q .

= (UA) 162oC

What is Q?

= (0.25) π (0.075) L (162) kW

.

Heat Exchanger Design – LMTD method

Solution continued:

Page 17: R. Shanthini 19 May 2010 1 Content of Lectures 10 to 12: Heat transfer: Mathematical problems on heat transfer PM3125: Lectures 10 to 12

R. Shanthini 19 May 2010

17

Q .

= (UA) 162oC = (0.25) π (0.075) L (162) kW

. .Q = mc cc (Tc,out – Tc,in) = mh ch (Th,in – Th,out)

.

= (200 kg/h) (1.13 kJ/kg K) (350 – 100)oC = 15.69 kW

Substituting the above in

we getL = 1.64 m

Heat Exchanger Design – LMTD method

Solution continued:

Page 18: R. Shanthini 19 May 2010 1 Content of Lectures 10 to 12: Heat transfer: Mathematical problems on heat transfer PM3125: Lectures 10 to 12

R. Shanthini 19 May 2010

18

(ii) Counter flow:

ΔTa = 350 – 34.63 = 315.37oC

ΔTb = 100 – 25 = 75oC

ΔTa - ΔTb

ln(ΔTa / ΔTb)

315.37 – 75

ln(315.37 / 75)= = 167.35oC

.= (UA) 167.35oC

Q = 15.69 kW; U = 0.25 kW/m2 K ;.

A = π (0.075) L m2

Therefore, L = 1.59 m

Heat Exchanger Design – LMTD method

Q = U A ΔTLM

ΔTLM =

Solution continued:

Page 19: R. Shanthini 19 May 2010 1 Content of Lectures 10 to 12: Heat transfer: Mathematical problems on heat transfer PM3125: Lectures 10 to 12

R. Shanthini 19 May 2010

19

Other Heat Exchanger Types

Cross-flow heat exchanger with both fluids unmixed

The direction of fluids are perpendicular to each other.

In calculating the required surface area for this heat exchanger using LMTD method, a correction factor is applied on the LMTD.

It lies between the required surface area for counter-flow and parallel-flow heat exchangers.

Page 20: R. Shanthini 19 May 2010 1 Content of Lectures 10 to 12: Heat transfer: Mathematical problems on heat transfer PM3125: Lectures 10 to 12

R. Shanthini 19 May 2010

20

Cross-flow Heat Exchanger Analysis

Q = U A ΔTLM .

ΔTLM = F ΔTLM, counter flow

Page 21: R. Shanthini 19 May 2010 1 Content of Lectures 10 to 12: Heat transfer: Mathematical problems on heat transfer PM3125: Lectures 10 to 12

R. Shanthini 19 May 2010

21

In calculating the required surface area for this heat exchanger using LMTD method, a correction factor is applied on the LMTD.

Other Heat Exchanger Types

One shell pass and two tube passes

T1

T2

t2

t1

t2

t1

T1

T2

Page 22: R. Shanthini 19 May 2010 1 Content of Lectures 10 to 12: Heat transfer: Mathematical problems on heat transfer PM3125: Lectures 10 to 12

R. Shanthini 19 May 2010

22

Two shell passes and two tube passes

Th,inTc,in

Th,out

Tc,out

Other Heat Exchanger Types

In calculating the required surface area for this heat exchanger using LMTD method, a correction factor is applied on the LMTD.

Page 23: R. Shanthini 19 May 2010 1 Content of Lectures 10 to 12: Heat transfer: Mathematical problems on heat transfer PM3125: Lectures 10 to 12

R. Shanthini 19 May 2010

23

Multi-pass Heat Exchanger Analysis

See lecture notes provided separately.

Page 24: R. Shanthini 19 May 2010 1 Content of Lectures 10 to 12: Heat transfer: Mathematical problems on heat transfer PM3125: Lectures 10 to 12

R. Shanthini 19 May 2010

24

Batch Sterilization (method of heating):

Steamheating

Electricalheating

Direct steam sparging

Page 25: R. Shanthini 19 May 2010 1 Content of Lectures 10 to 12: Heat transfer: Mathematical problems on heat transfer PM3125: Lectures 10 to 12

R. Shanthini 19 May 2010

25

For batch heating with constant rate heat flow:

M c (T - T0) q t =

M - mass of the medium

T0 - initial temperature of the medium

T - final temperature of the medium

c - specific heat of the medium

q - rate of heat transfer from the electrical coil to the medium

t - duration of electrical heating

Electricalheating

Total heat lost by the coil to the medium

= heat gained by the medium

.

.

Page 26: R. Shanthini 19 May 2010 1 Content of Lectures 10 to 12: Heat transfer: Mathematical problems on heat transfer PM3125: Lectures 10 to 12

R. Shanthini 19 May 2010

26

Direct steam sparging = (M + ms t) c (T – T0)

M - initial mass of the raw medium

T0 - initial temperature of the raw medium

ms - steam mass flow rate

t - duration of steam sparging

H - enthalpy of steam relative to the enthalpy at the initial temperature of the raw medium (T0)

T - final temperature of the mixture

c - specific heat of medium and water

For batch heating by direct steam sparging:

.

(ms t) (H + cT0).

+ M c T0 = (M + mst) c T.

ms t H. .

Page 27: R. Shanthini 19 May 2010 1 Content of Lectures 10 to 12: Heat transfer: Mathematical problems on heat transfer PM3125: Lectures 10 to 12

R. Shanthini 19 May 2010

27

For batch heating with isothermal heat source:

M - mass of the medium

T0 - initial temperature of the medium

TH - temperature of heat source (steam)

T - final temperature of the medium

c - specific heat of the medium

t - duration of steam heating

U - overall heat transfer coefficient

A - heat transfer area

Steamheating

U A t = M c ln

Could you prove the above?

T0 - TH ( ) T - TH

Page 28: R. Shanthini 19 May 2010 1 Content of Lectures 10 to 12: Heat transfer: Mathematical problems on heat transfer PM3125: Lectures 10 to 12

R. Shanthini 19 May 2010

28

For batch heating with isothermal heat source:

Steamheating

T = TH + (T0 - TH) exp - U A tc M ( )

U A t = M c lnT0 - TH ( ) T - TH

Page 29: R. Shanthini 19 May 2010 1 Content of Lectures 10 to 12: Heat transfer: Mathematical problems on heat transfer PM3125: Lectures 10 to 12

R. Shanthini 19 May 2010

29

Example of batch heating by direct steam sparging:

A fermentor containing 40 m3 medium at 25oC is going to be sterilized by direct injection of saturated steam. The steam at 350 kPa absolute pressure is injected with a flow rate of 5000 kg/hr, which will be stopped when the medium temperature reaches 122oC. Determine the time taken to heat the medium.

Additional data required:

Enthalpy of saturated steam at 350 kPa = ??

Enthalpy of water at 25oC = ??

The heat capacity of the medium 4.187 kJ/kg.K

The density of the medium are 4.187 kJ/kg.K and 1000 kg/m3, respectively.)

Page 30: R. Shanthini 19 May 2010 1 Content of Lectures 10 to 12: Heat transfer: Mathematical problems on heat transfer PM3125: Lectures 10 to 12

R. Shanthini 19 May 2010

30

Example of batch heating by direct steam sparging:

A fermentor containing 40 m3 medium at 25oC is going to be sterilized by direct injection of saturated steam. The steam at 350 kPa absolute pressure is injected with a flow rate of 5000 kg/hr, which will be stopped when the medium temperature reaches 122oC. Determine the time taken to heat the medium.

Additional data: The enthalpy of saturated steam at 350 kPa and water at 25oC are 2732 and 105 kJ/kg, respectively. The heat capacity and density of the medium are 4.187 kJ/kg.K and 1000 kg/m3, respectively.

Solution:

= (M + ms t) c (T – T0)ms t H. .

Use the equation below:

Page 31: R. Shanthini 19 May 2010 1 Content of Lectures 10 to 12: Heat transfer: Mathematical problems on heat transfer PM3125: Lectures 10 to 12

R. Shanthini 19 May 2010

31

(5000 kg/hr) (th) (2732-105) kJ/kg

= [(40 m3)(1000 kg/m3) + (5000 kg/hr)(th)](4.187 kJ/kg.K)(122-25)K

(5000 th) (2627) kJ = [40000 + 5000 t](4.187)(97)kJ

Taking the heating time (th) to be in hr, we get

(5000 th) [2627 – 4.187 x 97] = 40000 x 4.187 x 97

th = 1.463 hr

= (M + ms t) c (T – T0)ms t H. .

Therefore, the time taken to heat the medium is 1.463 hours.

Page 32: R. Shanthini 19 May 2010 1 Content of Lectures 10 to 12: Heat transfer: Mathematical problems on heat transfer PM3125: Lectures 10 to 12

R. Shanthini 19 May 2010

32

Example of batch heating with isothermal heat source:

A fermentor containing 40 m3 medium at 25oC is going to be sterilized by an isothermal heat source, which is saturated steam at 350 kPa absolute pressure. Heating will be stopped when the medium temperature reaches 122oC. Determine the time taken to heat the medium.

Additional data: The saturated temperature of steam at 350 kPa is 138.9oC. The heat capacity and density of the medium are 4.187 kJ/kg.K and 1000 kg/m3, respectively. U = 2500 kJ/hr.m2.K and A = 40 m2

Solution:

Use the equation below:

U A t = M c lnT0 - TH ( ) T - TH

Page 33: R. Shanthini 19 May 2010 1 Content of Lectures 10 to 12: Heat transfer: Mathematical problems on heat transfer PM3125: Lectures 10 to 12

R. Shanthini 19 May 2010

33

U A t = M c lnT0 - TH ( ) T - TH

Taking the heating time (th) to be in hr, we get

(2500 kJ/hr.m2.K) (40 m2) (tc)

= (40 m3) (1000 kg/m3) (4.187 kJ/kg.K) ln[(25-138.9)/(122-138.9)]

(2500 kJ/K) (40) (th) = (40) (1000) (4.187 kJ/K) ln[113.9/16.9]

(2500 kJ/K) (40) (th) = (40) (1000) (4.187 kJ/K) (1.908)

th = 3.1955 hr

Therefore, the time taken to heat the medium is 3.1955 hours.

Page 34: R. Shanthini 19 May 2010 1 Content of Lectures 10 to 12: Heat transfer: Mathematical problems on heat transfer PM3125: Lectures 10 to 12

R. Shanthini 19 May 2010

34

Explain why heating with isothermal heat source takes twice the time taken by heating with steam sparging, even though we used the same steam.

Page 35: R. Shanthini 19 May 2010 1 Content of Lectures 10 to 12: Heat transfer: Mathematical problems on heat transfer PM3125: Lectures 10 to 12

R. Shanthini 19 May 2010

35

2) Briefly discuss the applications of the following laws.

(a) Fourier’s law [2.5 marks](b) Newton’s law of cooling [2.5

marks](c) Stefan–Boltzmann Law [2.5

marks](d) Fick’s law [2.5

marks]

Questions from PM3125 / Feb 2010 past paper

1)(a) What are the major differences between diffusion and

convection? [5 marks](b) What are the major differences among conduction, convection

and radiation? [5 marks]

Page 36: R. Shanthini 19 May 2010 1 Content of Lectures 10 to 12: Heat transfer: Mathematical problems on heat transfer PM3125: Lectures 10 to 12

R. Shanthini 19 May 2010

36

8) A steel pipeline (inside diameter = 52.50 mm; outside diameter = 60.32 mm) contains saturated steam at 121.1oC. The line is insulated with 25.4 mm of asbestos. Assume that the inside surface temperature of the metal wall is at 121.1oC and the outer surface of the insulation is at 26.7oC. Taking the average value of ksteel as 45 W/m.K and that of

kasbestos as 0.182 W/m.K, calculate the following:

(a) Heat loss for 30.5 m of pipe length. [10 marks]

(b) Mass (in kg) of steam condensed per hour in the pipe due to the heat loss. [10 marks]

Additional data given on the next slide:

Questions from PM3125 / Feb 2010 past paper

Page 37: R. Shanthini 19 May 2010 1 Content of Lectures 10 to 12: Heat transfer: Mathematical problems on heat transfer PM3125: Lectures 10 to 12

R. Shanthini 19 May 2010

37

Additional Data:

i) Heat transfer rate through the pipe wall is given by,

asbestos

23

steel

12

21

)/ln()/ln(

)(2

k

rr

k

rr

TTLQ

where L is the length of pipe, T1 and T2 are the respective temperatures at

the inner and outer surfaces of the insulated pipe, r1 and r2 are the

respective inner and outer radius of the steel pipe, and r3 is the outer

radius of the insulated pipe. ii) Latent heat of vapourization of steam could be taken as 2200 kJ/kg.

Questions from PM3125 / Feb 2010 past paper

Page 38: R. Shanthini 19 May 2010 1 Content of Lectures 10 to 12: Heat transfer: Mathematical problems on heat transfer PM3125: Lectures 10 to 12

R. Shanthini 19 May 2010

38

1) A fermentor containing 40 m3 medium at 25oC is sterilized by an isothermal heat source, which is saturated steam at 4.5 bar absolute pressure. Determine the time taken to heat the medium to 122oC. [06 marks]

Additional Data: - The density of the medium is 1000 kg/m3.- You may use the heat exchange equation

where U (= 695 W/m2.K) is the overall heat transfer coefficient, A (= 40 m2) is the heat transfer area, t is the time taken, M is the mass of the medium, c (= 4.187 kJ/kg.K) is the heat capacity of the medium, T0 is the initial temperature of the medium, T is the final temperature of the medium and TH is the temperature of the heat source.

- You may use the steam table (provided) if necessary.

H

H

TT

TTMcUAt 0ln

Questions from PM3125 / Sept 2010 past paper

Page 39: R. Shanthini 19 May 2010 1 Content of Lectures 10 to 12: Heat transfer: Mathematical problems on heat transfer PM3125: Lectures 10 to 12

R. Shanthini 19 May 2010

39

2) Fluid A enters the shell-side of a shell-and-tube heat exchanger at a flow rate of 3.60 kg/s at 200oC and it leaves the heat exchanger at 150oC. Fluid B enters the tube-side of the exchanger at a flow rate of 2.88 kg/s at 100oC. Heat capacities of fluids A and B may be taken as 3.0 and 2.5 kJ/kg.K, respectively.

a) Show that fluid B leaves the heat exchanger at 175oC under the ideal conditions of no heat loss to the environment. [04 marks]

b) Should the heat exchanger be operated at parallel flow mode or at counter-current flow mode? Explain your answer. [04 marks]

c) Show that the log mean temperature driving force is about 36. [04

marks]d) Determine the overall heat transfer area required assuming an overall

heat transfer coefficient of 414 W/m2.K. [06 marks]

e) If the length of the tube-bundle of the heat exchanger is limited to 2 m, determine the number of 25 mm internal diameter tubes required to construct the tube-bundle. [04 marks]

Questions from PM3125 / Sept 2010 past paper

Page 40: R. Shanthini 19 May 2010 1 Content of Lectures 10 to 12: Heat transfer: Mathematical problems on heat transfer PM3125: Lectures 10 to 12

R. Shanthini 19 May 2010

40

Critical Radius of Insulation

Ti – ToQ.

[ln(ro/ri)] /2πkPL

=

ri

ro

Ti

Tor Pipe

+ [ln(r/ro)] /2πkIL + 1/hairA

Insulation

Pipe resistance could be neglected

A = 2 π r L

Additional material not used in the lectures.

Page 41: R. Shanthini 19 May 2010 1 Content of Lectures 10 to 12: Heat transfer: Mathematical problems on heat transfer PM3125: Lectures 10 to 12

R. Shanthini 19 May 2010

41

Ti – ToQ.

= [ln(r/ro)] /2πkIL + 1/(hair 2πrL)

2π L ( Ti – To) =

[ln(r/ro)] /kI + 1/(hair r)Convective resistance

Insulationresistance

Increasing r decreases convective resistance and increases heat transfer.

Increasing r increases insulation resistance and decreases heat transfer.

Critical Radius of Insulation Additional material not used in the lectures.

Page 42: R. Shanthini 19 May 2010 1 Content of Lectures 10 to 12: Heat transfer: Mathematical problems on heat transfer PM3125: Lectures 10 to 12

R. Shanthini 19 May 2010

42

= 0 at the critical radius of insulation,

which leads to rcr = kI / hair

If the outer radius of the pipe (ro) < rcr and if insulation is added to the pipe, heat losses will first

increase and go through a maximum at the insulation radius of rcr and then decrease.

dQ.

/dr

If the outer radius of the pipe (ro) > rcr and if insulation is added to the pipe, heat losses will

continue to decrease.

Critical Radius of Insulation Additional material not used in the lectures.