Textbook of Heat Transfer

Scilab Textbook Companion for
Textbook Of Heat Transfer

by Sukhatme, S. P.
1
Created by
Tushar Goel
B.Tech
Chemical Engineering
Indian Institute of Technology (BHU) , Varanasi
College Teacher
Dr. Prakash Kotecha
Cross-Checked by
August 23, 2012
1
Funded by a grant from the National Mission on Education through ICT,
http://spoken-tutorial.org/NMEICT-Intro. This Textbook Companion and Scilab
codes written in it can be downloaded from the Textbook Companion Project
section at the website http://scilab.in
Book Description
Title: Textbook Of Heat Transfer
Author: Sukhatme, S. P.
Publisher: Universities Press
Edition: 4
Year: 2005
ISBN: 9788173715440
1
Scilab numbering policy used in this document and the relation to the
above book.
Exa Example (Solved example)
Eqn Equation (Particular equation of the above book)
AP Appendix to Example(Scilab Code that is an Appednix to a particular
Example of the above book)
For example, Exa 3.51 means solved example 3.51 of this book. Sec 2.3 means
a scilab code whose theory is explained in Section 2.3 of the book.
2
Contents
List of Scilab Codes 4
1 Introduction 7
2 Heat Conduction in Solids 11
3 Thermal Radiation 35
4 Principles of Fluid Flow 52
5 Heat Transfer by Forced Convection 60
6 Heat Transfer by Natural convection 78
7 Heat Exchangers 88
8 Condensation and boiling 98
9 Mass Transfer 109
3
List of Scilab Codes
Exa 1.1 Viscosity in SI system . . . . . . . . . . . . . . . . . . 7
Exa 1.2 Useful heat gain and thermal eciency . . . . . . . . . 8
Exa 1.3 Exit velocity and Temperature . . . . . . . . . . . . . 9
Exa 2.1 Heat ow rate . . . . . . . . . . . . . . . . . . . . . . 11
Exa 2.2 Heat ow rate . . . . . . . . . . . . . . . . . . . . . . 12
Exa 2.3 Engineers decision . . . . . . . . . . . . . . . . . . . . 13
Exa 2.4 Thickness of insulation . . . . . . . . . . . . . . . . . . 14
Exa 2.5 Heat loss rate . . . . . . . . . . . . . . . . . . . . . . . 16
Exa 2.6 Critical radius . . . . . . . . . . . . . . . . . . . . . . 17
Exa 2.7 Maximum temperature . . . . . . . . . . . . . . . . . 18
Exa 2.8 Steady state temperature . . . . . . . . . . . . . . . . 19
Exa 2.9 Time taken by the rod to heat up . . . . . . . . . . . 20
Exa 2.10.i Heat transfer coecient at the centre . . . . . . . . . 21
Exa 2.10.ii heat transfer coecient at the surface . . . . . . . . . 23
Exa 2.11.a Time taken by the centre of ball . . . . . . . . . . . . 24
Exa 2.11.b time taken by the centre of ball to reach temperature . 25
Exa 2.12 Temperature at the centre of the brick . . . . . . . . . 26
Exa 2.13.a Temperature at the copper n tip . . . . . . . . . . . 28
Exa 2.13.b Temperature at the steel n tip . . . . . . . . . . . . . 29
Exa 2.13.c Temperature at the teon n tip . . . . . . . . . . . . 30
Exa 2.14 Heat loss rate . . . . . . . . . . . . . . . . . . . . . . . 31
Exa 2.15 Decrease in thermal resistance . . . . . . . . . . . . . 32
Exa 2.16 Overall heat transfer coecient . . . . . . . . . . . . . 33
Exa 3.1 Monochromatic emissive power . . . . . . . . . . . . . 35
Exa 3.2 Heat ux . . . . . . . . . . . . . . . . . . . . . . . . . 36
Exa 3.3 Absorbed radiant ux and absorptivity and reectivity 37
Exa 3.4.a Total intensity in normal direction . . . . . . . . . . . 38
Exa 3.4.b Ratio of radiant ux to the emissive power . . . . . . 39
4
Exa 3.5 Rate of incident radiation . . . . . . . . . . . . . . . . 40
Exa 3.6 Shape factor F12 . . . . . . . . . . . . . . . . . . . . . 40
Exa 3.7 Shape factor . . . . . . . . . . . . . . . . . . . . . . . 42
Exa 3.8 Shape factor F12 . . . . . . . . . . . . . . . . . . . . . 43
Exa 3.9 Shape factor . . . . . . . . . . . . . . . . . . . . . . . 44
Exa 3.10 Net radiative heat transfer . . . . . . . . . . . . . . . 44
Exa 3.11 steady state heat ux . . . . . . . . . . . . . . . . . . 45
Exa 3.12 Rate of heat loss . . . . . . . . . . . . . . . . . . . . . 47
Exa 3.13 Rate of nitrogen evaporation . . . . . . . . . . . . . . 47
Exa 3.14 Rate of energy loss from satellite . . . . . . . . . . . . 48
Exa 3.15 Net radiative heat transfer . . . . . . . . . . . . . . . 49
Exa 4.1 Pressure drop in smooth pipe . . . . . . . . . . . . . . 52
Exa 4.2.a Pressure drop and maximum velocity calculation . . . 53
Exa 4.2.b Pressure drop and maximum velocity calculation . . . 54
Exa 4.3 Pressure drop and power needed . . . . . . . . . . . . 56
Exa 4.4 Thickness of velocity boundary layer . . . . . . . . . . 57
Exa 4.5 Drag coecient and drag force . . . . . . . . . . . . . 58
Exa 5.1.a Local heat transfer coecient . . . . . . . . . . . . . . 60
Exa 5.1.b Wall temperature . . . . . . . . . . . . . . . . . . . . 61
Exa 5.2 ratio of thermal entrance length to entrance length . . 62
Exa 5.3.i Length of tube . . . . . . . . . . . . . . . . . . . . . . 63
Exa 5.3.ii Exit water temperature . . . . . . . . . . . . . . . . . 64
Exa 5.4 Length of tube over which temperature rise occurs . . 66
Exa 5.5 Rate of heat transfer to the plate . . . . . . . . . . . . 68
Exa 5.6.i Heat transfer rate . . . . . . . . . . . . . . . . . . . . 69
Exa 5.6.ii Average wall tempeature . . . . . . . . . . . . . . . . 70
Exa 5.7.i Pressure drop . . . . . . . . . . . . . . . . . . . . . . . 72
Exa 5.7.ii Exit temperature of air . . . . . . . . . . . . . . . . . 73
Exa 5.7.iii Heat transfer rate . . . . . . . . . . . . . . . . . . . . 75
Exa 6.1 Average nusselt number . . . . . . . . . . . . . . . . . 78
Exa 6.2 Reduce the equation . . . . . . . . . . . . . . . . . . . 80
Exa 6.3 Time for cooling of plate . . . . . . . . . . . . . . . . 81
Exa 6.4 True air temperature . . . . . . . . . . . . . . . . . . . 83
Exa 6.5 Rate of heat ow by natural convection . . . . . . . . 85
Exa 6.6 Average Heat transfer coecient . . . . . . . . . . . . 86
Exa 7.1 Heat transfer coecient . . . . . . . . . . . . . . . . . 88
Exa 7.2 Area of heat exchanger . . . . . . . . . . . . . . . . . 89
Exa 7.3 Mean temperature dierence . . . . . . . . . . . . . . 90
5
Exa 7.4.a Area of heat exchanger . . . . . . . . . . . . . . . . . 91
Exa 7.4.b Exit temperature of hot and cold streams . . . . . . . 92
Exa 7.5 Exit Temperature . . . . . . . . . . . . . . . . . . . . 94
Exa 8.1 Average Heat Transfer Coecient . . . . . . . . . . . 98
Exa 8.2 Average heat transfer coecient and lm Reynolds num-
ber . . . . . . . . . . . . . . . . . . . . . . . . . . . . 99
Exa 8.3 Length of the tube . . . . . . . . . . . . . . . . . . . . 101
Exa 8.4 boiling regions . . . . . . . . . . . . . . . . . . . . . . 103
Exa 8.5 Initial heat transfer rate . . . . . . . . . . . . . . . . . 107
Exa 9.1 Composition on molar basis . . . . . . . . . . . . . . . 109
Exa 9.2 Diusion coecient of napthalene . . . . . . . . . . . 110
Exa 9.3.a Rate of hydrogen diusion . . . . . . . . . . . . . . . . 110
Exa 9.3.b Rate of hydrogen diusion . . . . . . . . . . . . . . . . 111
Exa 9.4.a Rate of loss of ammonia . . . . . . . . . . . . . . . . . 112
Exa 9.4.b Rate at which air enters the tank . . . . . . . . . . . . 113
Exa 9.5 Rate of evaporation . . . . . . . . . . . . . . . . . . . 114
Exa 9.6 Rate of evaporation . . . . . . . . . . . . . . . . . . . 115
Exa 9.7.a Mass transfer coecient Colburn anology . . . . . . . 116
Exa 9.7.b Mass transfer coecient Gnielinski equation . . . . . . 117
Exa 9.7.c To show mass ux of water vapour is small . . . . . . 119
Exa 9.8 Mass fraction . . . . . . . . . . . . . . . . . . . . . . . 120
6
Chapter 1
Introduction
Scilab code Exa 1.1 Viscosity in SI system
1 clear;
2 clc;
3 // A Textbook on HEAT TRANSFER by S P SUKHATME
4 // Chapter 1
5 // I nt r o duc t i o n
6
7
8 // Example 1 . 1
9 // Page 5
10 // Gi ven t hat t he v i s c o s i t y o f wat er at 100 de gr e e
Ce l s i u s i s 28. 8 106 kgf s /m2 i n MKS system ,
e xpr e s s t h i s val ue i n SI system .
11 printf(Example 1 . 1 , Page 5 \n \n)
12
13 // So l ut i o n :
14
15 // at 100 de gr e e Ce l s i u s
16 v1=28.8 * 10^ -6; // [ kgf s /m 2 ]
17 v2=28.8 * 10^-6 * 9.8; // [ N s /m 2 ]
18 printf( Vi s c o s i t y o f wat er at 100 de gr e e c e l s i u s i n
t he SI system i s %e N. s /m2 ( or kg/m s ) ,v2)
7
Scilab code Exa 1.2 Useful heat gain and thermal eciency
1 clear all;
2 clc;
3 // Textbook o f Heat Tr ans f e r ( 4 th Edi t i on ) ) , S P
Sukhatme
4 // Chapter 1 I nt r o duc t i o n
5
6 // Example 1 . 2
7 // Page 14
8 printf(Example 1 . 2 , Page 14 \n \n)
9 // So l ut i o n :
10 i=950; // r a d i a t i o n f l u x [W/m 2 ]
11 A=1.5; // ar e a [m 2 ]
12 T_i =61; // i n l e t t e mpe r at ur e
13 T_o =69; // o u t l e t t e mpe r at ur e
14 mdot =1.5; // [ kg/min ] , mass f l ow r a t e
15 Mdot =1.5/60; // [ kg/ s e c ]
16 Q_conductn =50; // [W]
17 t=0.95; // t r a n s mi s s i v i t y
18 a=0.97; // a b s o p t i v i t y
19 // f rom appendi x t a bl e A. 1 at 65 de gr e e C
20 C_p= 4183 ; // [ J/kg K]
21 // Usi ng Equat i on 1 . 4 . 1 5 , assumi ng t hat t he f l ow
t hr ough t he t ube s i s s t e ady and one di me ns i onal .
22 // i n t h i s c as e (dW/ dt ) s h a f t = 0
23 // assumi ng (dW/ dt ) s he a r i s n e g l i g i b l e
24 // eqn ( 1 . 4 . 1 5 ) r e duc e s t o
25 q=Mdot*C_p*(T_o -T_i);
26
27 // l e t n be t her mal e f f i c i e n c y
28 n=q/(i*A);
29 n_percent=n*100;
30
8
31
32 // e quat i on 1 . 4 . 1 3 y i e l d s dQ/ dt = 0
33 Q_re_radiated =(i*A*t*a)-Q_conductn -q; // [W]
34
35
36 printf( Us e f ul heat gai n r a t e i s %f W \n,q);
37 printf( Thermal e f f i c i e n c y i s %e i . e . %f per c e nt \n
,n,n_percent);
38 printf(The r a t e at whi ch ener gy i s l o s t by re
r a d i a t i o n and c o nve c t i o n i s %f W,Q_re_radiated)
Scilab code Exa 1.3 Exit velocity and Temperature
1 clear all;
2 clc;
4 // Chapter 1
5 // I nt r o duc t i o n
6
7
8 // Example 1 . 3
9 // Page 16
10 printf(Example 1 . 3 , Page 16\n\n);
11
12 // So l ut i o n :
13 // Gi ven
14 v_i =10; // [m/ s ]
15 q=1000; // [W]
16 d_i =0.04; // [m]
17 d_o =0.06; // [m]
18
19 // From appendi x t a bl e A. 2
20 rho1 =0.946; // [ kg/m 3 ] at 100 de gr e e C
21 C_p =1009; // [ J/kg K]
22
9
23 mdot=rho1*(%pi/4)*(d_i ^2)*v_i; // [ kg/ s ]
24
25
26 // I n t h i s c as e (dW/ dt ) s h a f t =0 and ( z o z i )=0
27 // From eqn 1 . 4 . 1 5 , q=mdot ( h oh i )
28 // Let dh = ( h oh i )
29 dh=q/mdot; // [ J/kg ]
30 // Let T o be t he o u t l e t t e mpe r at ur e
31 T_o=dh/C_p +100;
32
33 rho2 =0.773; // [ kg/m 3 ] at T o = 183. 4 de gr e e C
34 // From eqn 1 . 4 . 6
35 v_o=mdot/(rho2*(%pi/4)*(d_o)^2); // [m/ s ]
36
37 dKE_kg =(v_o^2-v_i^2)/2; // [ J/kg ]
38
39
40 printf( Exi t Temperature i s %f de gr e e C \n,T_o);
41 printf( Exi t v e l o c i t y i s %f m/ s \n,v_o);
42 printf(Change i n Ki ne t i c Energy per kg = %f J/kg,
dKE_kg);
10
Chapter 2
Heat Conduction in Solids
Scilab code Exa 2.1 Heat ow rate
1 clear all;
2 clc;
3
5 // Chapter 2
6 // Heat Conducti on i n S o l i d s
7
8 // Example 2 . 1
9 // Page 27
10 printf(Example 2 . 1 , Page 27 \n\n)
11
12 d_i =0.02; // [m] i nne r r a di us
13 d_o =0.04; // [m] out e r r a di us
14 r_i=d_i/2; // [m] i nne r r a di us
15 r_o=d_o/2; // [m] out e r r a di us
16 k=0.58; // [ w/m K] t her mal c o nduc t i v i t y o f tube
ma t e r i a l
17 t_i =70; // [ de gr e e C]
18 t_o =100; // [ de gr e e C]
19 l=1; // [m] per uni t l e ng t h
20 // us i ng e quat i on 2 . 1 . 5
11
21 q=l*2*( %pi)*k*(t_i -t_o)/log(r_o/r_i);
22 printf( Heat f l ow per uni t l e ng t h i s %f W/m,q);
Scilab code Exa 2.2 Heat ow rate
1 clear;
2 clc;
3
5 // Chapter 2
7
8 // Example 2 . 2
9 // Page 31
11
12 d_i =0.02; // [m] i nne r r a di us
13 d_o =0.04; // [m] out e r r a di us
14 r_i=d_i/2; // [m] i nne r r a di us
15 r_o=d_o/2; // [m] out e r r a di us
16 k=0.58; // [ w/m K] t her mal c o nduc t i v i t y o f tube
ma t e r i a l
17 t_i =70; // [ de gr e e C]
18 t_o =100; // [ de gr e e C]
19 l=1; // [m] per uni t l e ng t h
20
21 // t her mal r e s i s t a n c e o f tube per uni t l e ng t h
22 R_th_tube =(log(r_o/r_i))/(2* %pi*k*l); // [ K/W]
23
24 // f rom t a bl e 1 . 3 , heat t r a n s f e r coe f f i c i e n t f o r
c onde ns i ng steam may be t aken as
25 h=5000; // [W/m2 K]
26 // t her mal r e s i s t a n c e o f c onde ns i ng steam per uni t
l e ng t h
27 R_th_cond =1/( %pi*d_o*l*h);
12
28
29 // s i n c e R t h t ube i s much l e s s than R t h cond , we
can assume out e r s u r f a c e t o be at 100 de gr e e C
30 // hence heat f l ow r a t e per uni t meter i s
31 q=l*2*( %pi)*k*(t_i -100)/log(r_o/r_i);
32
33 printf( Thermal r e s i s t a n c e o f tube per uni t l e ng t h
i s %f K/W\n,R_th_tube);
34 printf( Thermal r e s i s t a n c e o f c onde ns i ng steam per
uni t l e ng t h i s %f K/W\n,R_th_cond);
35 printf( Heat f l ow per uni t l e ng t h i s %f W/m,q);
Scilab code Exa 2.3 Engineers decision
1 clear all;
2 clc;
3
5 // Chapter 2
7
8 // Example 2 . 3
9 // Page 31
11
12 h_w =140; // heat t r a n s f e r c o e f f i c i e n t on wat er s i de ,
[W/m2 K]
13 h_o =150; // heat t r a n s f e r c o e f f i c i e n t on o i l s i de , [
W/m2 K]
14 k=30; // t her mal c o nduc t i v i t y [W/m K]
15 r_o =0.01; // i nne r di ame t e r o f GI pi pe on i n s i d e
16 r_i =0.008; // out e r di ame t e r GI pi pe on i n s i d e
17 l=1; // [m] , per uni t l e ng t h
18
19 // Thermal r e s i s t a n c e o f i nne r GI pi pe
13
20 R_inner_GI=log((r_o/r_i))/(2* %pi*k*l);
21
22
23 // Thermal r e s i s t a n c e on t he o i l s i d e per uni t
l e ng t h
24 R_oilside =1/( h_o*%pi *2*r_i*l);
25
26
27 // Thermal r e s i s t a n c e on c ol d wat er s i d e per uni t
l e ng t h
28 R_waterside =1/( h_w*%pi *2*r_o*l);
29
30
31 // we s e e t her mal r e s i s t a n c e o f i nne r GI pi pe
c o n t r i b u t e s l e s s than 0 . 5 pe r c e nt t o t he t o t a l
r e s i s t a n c e
32
33
34 printf( Thermal r e s i s t a n c e o f i nne r GI pi pe = %f K/W
\n,R_inner_GI);
35 printf( Thermal r e s i s t a n c e on t he o i l s i d e per uni t
l e ng t h = %f K/W \n,R_oilside);
36 printf( Thermal r e s i s t a n c e on c ol d wat er s i d e per
uni t l e ng t h = %f K/W \n,R_waterside);
37 printf(So , Engi ne er i nc har ge has made a bad
d e c i s i o n );
Scilab code Exa 2.4 Thickness of insulation
1 clear all;
2 clc;
3
5 // Chapter 2
14
7
8 // Example 2 . 4
9 // Page 32
11
12 Ti = 300; // I n t e r n a l temp o f hot gas i n
de gr e e Ce l s i u s
13 OD = 0.1; // Outer di ame t e r o f l ong met al
pi pe i n met er s
14 ID = 0.04; // I n t e r n a l di amt e r e o f l ong met al
pi pe i n met er s
15 ki = 0.052; // t her mal c o nduc t i v i t y o f mi ne r al
wood i n W/mK
16 To = 50; // Outer s u r f a c e t e mpe r at ur e i n
de gr e e c e l s i u s
17 hi = 29; // heat t r a n s f e r c o e f f i c i e n t i n
t he i nne r s i d e i n W/m2 K
18 ho = 12; // heat t r a n s f e r c o e f f i c i e n t i n
t he out e r pi pe W/m2 K
19
20 // De t e r mi nat i on o f t h i c k n e s s o f i n s u l a t i o n
21 function[f] = thickness(r)
22 f = r*(10.344 + 271.15* log(r*(0.05) ^-1)) -11.75
23 funcprot (0);
24 endfunction
25 r = 0.082;
26 while 1
27 rnew = r - thickness(r)/diff(thickness(r));
28 if rnew == r then
29 r3 = rnew;
30 break;
31 end
32 r = rnew;
33 end
34 t = r3 - OD/2;
35 printf(\n Thi c kne s s o f i n s u l a t i o n = %f cm,t*100);
36 // Heat l o s s per uni t l e ng t h
37 q = 600*(22/7)*r3;
15
38 printf(\n Heat l o s s per uni t l e ng t h = %. 1 f W/m,q);
Scilab code Exa 2.5 Heat loss rate
1 clear all;
2 clc;
3
5 // Chapter 2
7
8 // Example 2 . 5
9 // Page 34
11
12 Ti = 90; //Temp on i nne r
s i d e i n de gr e e c e l s i u s
13 To = 30; //Temp on out e r
s i d e i n de gr e e c e l s i u s
14 hi = 500; // heat t r a n s f e r
c o e f f c i e n t i n W/m2 K
15 ho = 10; // heat t r a n s f e r
c o e f f c i e n t i n W/m2 K
16 ID = 0.016; // I n t e r n a l
di ame t e r i n met er s
17 t = [0 0.5 1 2 3 4 5]; // I n s u l a t i o n
t h i c k n e s s i n cm
18 OD = 0.02; // Outer di ame t e r
i n met er s
19 r3 = OD/2 + t/100; // r a di us a f t e r
i n s u l a t i o n i n met er s
20
21 i=1;
22 printf(\n I n s u l a t i o n t h i c k n e s s (cm) r 3 (m)
heat l o s s r a t e per meter (W/m) );
16
23 while i<=7
24 ql(i) = [2*( %pi)*(ID/2)*(Ti-To)]/[(1/ hi)
+(0.008/0.2)*log(r3(i)/0.01) + (0.008/ r3(i))
*(1/ho)];
25 printf(\n %. 1 f %. 3 f
%. 1 f ,t(i),r3(i),ql(i));
26 i = i+1;
27 end
28 plot(t,ql);
29 xtitle(, I n s u l a t i o n t h i c k n e s s (cm) , Heat l o s s r a t e
per uni t l engt h ,W/m);
30 printf(\n The maxima i n t he cur ve i s at r 3 = 0. 02
m);
Scilab code Exa 2.6 Critical radius
1 clear all;
2 clc;
3
5 // Chapter 2
7
8 // Example 2 . 6
9 // Page 36
11
12 h_natural = 10; // heat t r a n s f e r c o e f f i c i e n t f o r
na t ur a l c o nve c t i o n i n W/m2 K
13 h_forced = 50; // heat t r a n s f e r c o e f f i c i e n t f o r
f o r c e d c o nve c t i o n i n W/m2 K
14 // f o r a s be s t o s
15 k1 = 0.2; // t her mal c o nduc t i v i t y i n W/m K
16 // f o r mi ne r al wool
17 k2 = 0.05; // t her mal c o nduc t i v i t y i n W/m K
17
18 printf(\n c r i t i c a l r a di us o f i n s u l a t i o n i n cm);
19 printf(\n h = 10
h = 50 );
20 printf(\n As be s t os %. 1 f
%. 1 f ,k1 *100/ h_natural ,k1
*100/ h_forced);
21 printf(\n Mi ner al wool %. 1 f
%. 1 f ,k2 *100/ h_natural ,k2
*100/ h_forced);
Scilab code Exa 2.7 Maximum temperature
1 clear all;
2 clc;
3
5 // Chapter 2
7
8 // Example 2 . 7
9 // Page 43
11
12 H = 5 ; // Hei ght , [m]
13 L = 10 ; // Length , [m]
14 t = 1 ; // t hi c kne s s , [m]
15 b = t/2;
16 k = 1.05 ; // [W/m K]
17 q = 58 ; // [W/m 3 ]
18 T = 35 ; // [ C]
19 h = 11.6 ; // Heat t r a n s f e r c o e f f i c i e n t , [W/m2 K]
20 // Su b s t i t u t i ng t he va l ue s i n e quat i on 2 . 5 . 6
21 T_max = T + q*b*(b/(2*k)+1/h);
22 printf(Maximum Temperature = %f de gr e e C,T_max);
18
Scilab code Exa 2.8 Steady state temperature
1 clear all;
2 clc;
3
5 // Chapter 2
7
8 // Example 2 . 8
9 // Page 47
11
12 // The bar wi l l have two di me ns i onal v a r i a t i o n i n
t e mpe r at ur e
13 // t he d i f f e r e n t i a l e quat i on i s s ubj e c t t o boundary
c o n d i t i o n s
14 x1 = 0; // [ cm]
15 Tx1 = 30; // [ C]
16 x2 = 5; // [ cm]
17 Tx2 = 30; // [ C]
18 y1 = 0; // [ cm]
19 Ty1 = 30; // [ C]
20 y2 = 10; // [ cm]
21 Ty2 = 130; // [ C]
22 // s u b s t i t u t i n g t he t a = T30 and us i ng eqn 2 . 6 . 1 1
23 // put t i ng x = 2 . 5cm and y = 5cm i n i n f i n i t e
summati on s e r i e s
24
25
26 n = 1;
27 x1 = (1-cos(%pi*n))/(sinh (2*%pi*n))*sin(n^%pi/2)*
sinh(n*%pi);
28
19
29 n = 3;
sinh(n*%pi);
31
32 n = 5;
sinh(n*%pi);
34
35 x = x1+x3+x5;
36
37 T = x*100+30;
38 printf( Steady s t a t e t e mpe r at ur e = %f C,T);
Scilab code Exa 2.9 Time taken by the rod to heat up
1 clear all;
2 clc;
3
5 // Chapter 2
7
8 // Example 2 . 9
9 // Page 51
11
12 k = 330; // t her mal c o nduc t i v i t y i n W/m K
13 a = 95*10^( -6); // t her mal e xpans i on c o e f f i c i e n t
14 R = 0.01; // r a di us i n met er s
15 To = 77; // t e mpe r at ur e i n k e l v i n s
16 Tf = 273+50; // t e mpe r at ur e i n k e l v i n s
17 theta1 = To - Tf;
18 T = 273+10; // t e mpe r at ur e i n k e l v i n s
19 theta = T - Tf;
20 h = 20; // heat t r a n s f e r c o e f f i c i e n t i n W
20
/m2 K
21 printf(\n Theta1 = %d K,theta1);
22 printf(\n Theta = %d K ,theta);
23 printf(\n v/A = %. 3 f m,R/2);
24 printf(\n k/a = %. 4 f 10( 6) J/m3 K ,(k/a)*10^( -6))
;
25
26 time = (k/a)*(R/2)/h*log(theta1/theta);
27
28 printf(\n Time t aken by t he rod t o heat up = %. 1 f
s e c s ,time);
29 Bi = h*R/k;
30 printf(\n Bi ot number Bi = %. 2 f 10( 4) ,Bi *10^4);
31 printf(\n Si nc e Bi ot number i s much l e s s than 0 . 1 ,
t h e r e f o r e as s umpt i on t hat i n t e r n a l t e mpe r at ur e
g r a di e n t s ar e n e g l i g i b l e i s a good one );
Scilab code Exa 2.10.i Heat transfer coecient at the centre
1 clear all;
2 clc;
3
5 // Chapter 2
7
8 // Example 2 . 1 0 ( i )
9 // Page 58
10 printf(Example 2 . 1 0 ( i ) , Page 58 \n\n)
11
12 // Cent r e o f t he s l a b
13 // Gi ven data
14 b = 0.005 ; // [m]
15 t = 5*60; // ti me , [ s e c ]
16 Th = 200 ; // [ C]
21
17 Tw = 20 ; // [ C]
18 h = 150 ; // [W/m2 K]
19 rho = 2200 ; // [ kg/m 3 ]
20 Cp = 1050 ; // [ J/kg K]
21 k = 0.4 ; // [W/m K]
22 // Usi ng c ha r t s i n f i g 2. 18 and 2. 19 and eqn 2 . 7 . 1 9
and 2 . 7 . 2 0
23
24 theta = Th - Tw;
25 Biot_no = h*b/k;
26 a = k/(rho*Cp); // al pha
27 Fourier_no = a*t/b^2;
28
29 // From f i g 2 . 1 8 , r a t i o = t he t a x b0 / t he t a o
30 ratio_b0 = 0.12;
32 ratio_b1 = 0.48;
33
34 // The r e f or e
35 theta_x_b0 = theta*ratio_b0; // [ C]
36 T_x_b0 = theta_x_b0 + Tw ; // [ C]
38 T_x_b1 = theta_x_b1 + Tw ; // [ C]
39
40 // From Tabl e 2 . 2 f o r Bi = 1. 875
41 lambda_1_b = 1.0498;
42 x = 2*sin(lambda_1_b)/[ lambda_1_b +(sin(lambda_1_b))
*(cos(lambda_1_b))];
43
44 // From eqn 2 . 7 . 2 0
45 theta_x_b0 = theta*x*(exp((- lambda_1_b ^2)*Fourier_no
));
46 T_x_b0 = theta_x_b0 + Tw;
47 printf( Temperature at b=0 i s %f de gr e e C\n,T_x_b0)
;
22
Scilab code Exa 2.10.ii heat transfer coecient at the surface
1 clear all;
2 clc;
3
5 // Chapter 2
7
8 // Example 2 . 1 0 ( i i )
9 // Page 58
10 printf(Example 2 . 1 0 ( i i ) , Page 58 \n\n)
11
12 // ( i i ) Sur f a c e o f t he s l a b
13
14 b = 0.005 ; // [m]
15 t = 5*60; // ti me , [ s e c ]
16 Th = 200 ; // [ C]
17 Tw = 20 ; // [ C]
18 h = 150 ; // [W/m2 K]
19 rho = 2200 ; // [ kg/m 3 ]
20 Cp = 1050 ; // [ J/kg K]
21 k = 0.4 ; // [W/m K]
22 // Usi ng c ha r t s i n f i g 2. 18 and 2. 19 and eqn 2 . 7 . 1 9
and 2 . 7 . 2 0
23 theta = Th - Tw;
24 Biot_no = h*b/k;
25 a = k/(rho*Cp); // al pha
26 Fourier_no = a*t/b^2;
27
29 ratio_b0 = 0.12;
31 ratio_b1 = 0.48;
23
32
35 T_x_b0 = theta_x_b0 + Tw ; // [ C]
37 T_x_b1 = theta_x_b1 + Tw ; // [ C]
38
39 // From Tabl e 2 . 2 f o r Bi = 1. 875
40 lambda_1_b = 1.0498;
41 x = 2*sin(lambda_1_b)/[ lambda_1_b +(sin(lambda_1_b))
*(cos(lambda_1_b))];
42
43 // From 2 . 7 . 1 9
44 theta_x_b1 = theta_x_b0 *(cos(lambda_1_b *1));
45 T_x_b1 = theta_x_b1 + Tw;
46 printf( Temperature at b=1 i s %f de gr e e C\n,T_x_b1)
;
Scilab code Exa 2.11.a Time taken by the centre of ball
1 clear all;
2 clc;
3
5 // Chapter 2
7
8 // Example 2 . 1 1 ( a )
9 // Page 65
10 printf(Example 2 . 1 1 ( a ) , Page 65 \n\n)
11
12 D = 0.05 ; // [m]
13 To = 450 ; // [ de gr e e C]
14 Tf = 90 ; // [ de gr e e C]
15 T = 150 ; // [ de gr e e c ]
24
16 h = 115 ; // [W/m2 K]
17 rho = 8000 ; // [ kg/m 3 ]
18 Cp = 0.42*1000 ; // [ J/kg K]
19 k = 46 ; // [W/m K]
20 R = D/2;
21
22 // ( a )
23 // From eqn 2 . 7 . 3 f o r a s phe r e
24 t1 = rho*Cp*R/(3*h)*log((To-Tf)/(T-Tf)); // [ s e c ]
25 t1_min = t1/60 ; // [ min ]
26 printf(Time t aken by t he c e nt r e o f t he b a l l t o
r e ac h 150 de gr e e C i f i n t e r n a l g r a di e n t s ar e
ne g l e c t e d i s %f s e c onds i . e . %f mi nut es \n,t1 ,
t1_min);
Scilab code Exa 2.11.b time taken by the centre of ball to reach temper-
ature
1 clear all;
2 clc;
3
5 // Chapter 2
7
8 // Example 2 . 1 1 ( b)
9 // Page 65
10 printf(Example 2 . 1 1 ( b) , Page 65 \n\n)
11
12 D = 0.05 ; // [m]
13 To = 450 ; // [ de gr e e C]
14 Tf = 90 ; // [ de gr e e C]
15 T = 150 ; // [ de gr e e c ]
16 h = 115 ; // [W/m2 K]
17 rho = 8000 ; // [ kg/m 3 ]
25
18 Cp = 0.42*1000 ; // [ J/kg K]
19 k = 46 ; // [W/m K]
20 R = D/2;
21
22 // ( b)
23 // l e t r a t i o = t he t a R 0 / t he t a o
24 ratio = (T-Tf)/(To - Tf);
25 Bi = h*R/k;
26 // From Tabl e 2 . 5
27 lambda_1_R = 0.430;
28 x = 2*[sin(lambda_1_R) - lambda_1_R*cos(lambda_1_R)
]/[ lambda_1_R - sin(lambda_1_R)*cos(lambda_1_R)];
29
30 // Su b s t i t u t i ng i n e qua t t i o n 2 . 7 . 2 9 , we have an
e quat i on i n v a r i a b l e y(=at /R2)
31 // So l vi ng
32 function[eqn] = parameter(y)
33 eqn = ratio - x*exp(-( lambda_1_R ^2)*(y));
34 funcprot (0);
35 endfunction
36
37 y = 5; // ( i n i t i a l gues s , assumed val ue f o r f s o l v e
f unc t i o n )
38 Y = fsolve(y,parameter);
39
40 a = k/(Cp*rho); // al pha
41 t2 = Y*(R^2)/(a); // [ s e c ]
42 t2_min = t2/60; // [ min ]
43 printf(Time t aken by t he c e nt r e o f t he b a l l t o
r e ac h 150 de gr e e C i f i n t e r n a l t e mpe r at ur e
g r a di e n t s ar e not ne g l e c t e d i s %f s e c onds i . e . %f
mi nut es ,t2 ,t2_min);
Scilab code Exa 2.12 Temperature at the centre of the brick
26
1 clear ;
2 clc;
3
5 // Chapter 2
7
8 // Example 2. 12
9 // Page 67
10 printf(Example 2 . 1 2 , Page 67 \n\n)
11
12 a = 0.12 ; // [m]
13
14 T = 400 ; // [ C]
15 To = 25 ; // [ C]
16 t = 100/60 ; // [ hour ]
17 h = 10 ; // [W/m2 K]
18 k = 1.0 ; // [W/m K]
19 alpha = 3.33*10^ -3 ; // [m2/h ]
20 // us i ng f i g 2. 18 and eqn 2 . 7 . 2 0
21
22 x1 = h*a/k ;
23 x2 = k/(h*a);
24 x3 = alpha*t/a^2;
25
26 // Let r a t i o x = t he t a / t he t a o f o r x d i r e c t i o n , f rom
f i g 2. 18
27 ratio_x = 0.82 ;
28
29 // Si mi l a r l y , f o r y d i r e c t i o n
30 ratio_y = 0.41;
31
32 // Si mi l a r l y , f o r z d i r e c t i o n
33 ratio_z = 0.30;
34
36 total_ratio = ratio_x*ratio_y*ratio_z ;
37
27
38 T_centre = To + total_ratio *(T-To) ; // [ de gr e e C]
39 printf( Temperature at t he c e nt r e o f t he br i c k = %f
de gr e e C \n\n,T_centre);
40
41 // Al t e r n a t i v e l y
42 printf( Al t e r na t i v e l y , o bt a i ni ng Bi ot number and
va l ue s o f l ambda 1 b and us i ng eqn 2 . 7 . 2 0 , we ge t
\n)
43
44 ratio_x = 1.1310* exp ( -(0.9036^2) *0.385);
45 ratio_y = 1.0701* exp ( -(0.6533^2) *2.220);
46 ratio_z = 1.0580* exp ( -(0.5932^2) *3.469);
47 ratio = ratio_x*ratio_y*ratio_z;
48
49 T_centre = To + total_ratio *(T-To) ; // [ de gr e e C]
50 printf( Temperature at t he c e nt r e o f t he br i c k = %f
de gr e e C \n,T_centre);
Scilab code Exa 2.13.a Temperature at the copper n tip
1 clear all;
2 clc;
3
5 // Chapter 2
7
8 // Example 2 . 1 3 ( a )
9 // Page 73
10 printf(Example 2 . 1 3 ( a ) , Page 73 \n\n)
11
12 D = 0.003 ; // [m]
13 L = 0.03 ; // [m]
14 h = 10 ; // [W/m 2 ]
15 Tf = 20 ; // [ C]
28
16 T1 = 120 ; // [ C]
17
18 // ( a ) Copper f i n
19 k = 350 ; // [W/m K]
20
21 // For a c i r c u l a r c r o s s s e c t i o n
22 m = [4*h/(k*D)]^(1/2);
23 mL = m*0.03 ;
24 // T at x = L
25 T = Tf + (T1-Tf)/cosh(m*L);
26 printf(mL = %f \n,mL);
27 printf( Temperature at t he t i p o f f i n made o f copper
i s %f de gr e e C \n,T);
Scilab code Exa 2.13.b Temperature at the steel n tip
1 clear all;
2 clc;
3
5 // Chapter 2
7
8 // Example 2 . 1 3 ( b)
9 // Page 73
10 printf(Example 2 . 1 3 ( b) , Page 73 \n\n)
11
12 D = 0.003 ; // [m]
13 L = 0.03 ; // [m]
14 h = 10 ; // [W/m 2 ]
15 Tf = 20 ; // [ C]
16 T1 = 120 ; // [ C]
17
18
19 // ( b) S t a i n l e s s s t e e l f i n
29
20 k = 15 ; // [W/m K]
21
23 m = [4*h/(k*D)]^(1/2);
24 mL = m*0.03 ;
25 // T at x = L
28 printf( Temperature at t he t i p o f f i n made o f s t e e l
Scilab code Exa 2.13.c Temperature at the teon n tip
1 clear all;
2 clc;
3
5 // Chapter 2
7
8 // Example 2 . 1 3 ( c )
9 // Page 73
10 printf(Example 2 . 1 3 ( c ) , Page 73 \n\n)
11
12 D = 0.003 ; // [m]
13 L = 0.03 ; // [m]
14 h = 10 ; // [W/m 2 ]
15 Tf = 20 ; // [ C]
16 T1 = 120 ; // [ C]
17
18 // ( c ) Te f l on f i n
19 k = 0.35 ; // [W/m K]
20
22 m = [4*h/(k*D)]^(1/2);
30
23 mL = m*0.03 ;
24 // T at x = L
27 printf( Temperature at t he t i p o f f i n made o f t e f l o n
Scilab code Exa 2.14 Heat loss rate
1 clear all;
2 clc;
3
5 // Chapter 2
7
8 // Example 2. 14
9 // Page 74
11
12 L = 0.02 ; // [m]
13 t = 0.002 ; // [m]
14 b = 0.2 ; // [m]
15 theta1 = 200 ; // [ C]
16 h = 15 ; // [W/m2 K]
17 k = 45 ; // [W/m K]
18
19 Bi = h*(t/2)/k ;
20
21 // We have
22 P = 2*(b+t); // [m]
23 A = b*t ; // [m 2 ]
25 mL = ([(h*P)/(A*k)]^(1/2))*L;
26
31
27 // From e quat i on 2 . 8 . 6 , f i n e f f e c t i v e n e s s n
28 n = tanh(mL)/mL;
29 printf( Fi n Ef f e c t i v e n e s s = %f \n,n);
30
31 q_loss = n*h*40.4*2*10^ -4*200; // [W]
32 printf( Heat l o s s r a t e f rom f i n s u r f a c e = %f W,
q_loss);
Scilab code Exa 2.15 Decrease in thermal resistance
1 clear all;
2 clc;
3
5 // Chapter 2
7
8 // Example 2 . 1
9 // Page 74
11
12
13 // Fi nd De c r e as e i n t her mal Re s i s t a nc e
14 // Fi nd I n c r e a s e i n heat t r a n s f e r r a t e
15
16 h = 15 ; // [W/m 2 .K]
17 k = 300; // [W/m. K]
18 T = 200; // [ C]
19 Tsurr = 30; // [ C]
20 d = .01; // [m]
21 L = .1; // [m]
22 A = .5*.5 // [m 2 ]
23 n = 100 //Number o f Pi ns
24
25 Bi = h*d/2/k; // Bi ot Number
32
26 // Val ue o f Bi ot Number i s much l e s s than . 1
27 //Thus us i ng e quat i on 2 . 8 . 6
28 mL = (h*4/k/d)^.5*L;
29 zi = tanh(mL)/mL;
30 Res1 = 1/h/A; // Thermal r e s i s t a n c e wi t hout
f i ns , [ K/W]
31 Res2 = 1/(h*(A - n*%pi/4*d^2 + zi*(n*%pi*d*L)));//
Thermal r e s i s t a n c e wi th f i ns , [ K/W]
32
33 delRes = Res1 -Res2; // [ K/W]
34 // I n c r e a s e i n heat t r a n s f e r r a t e
35 q = (T-Tsurr)/Res2 - (T-Tsurr)/Res1; // [W]
36
37 printf(\n\n De c r e as e i n t her mal r e s i s t a n e at
s u r f a c e %. 4 f K/W. \ n I n c r e a s e i n heat t r a n s f e r
r a t e %. 1 f W,delRes ,q)
38 //END
Scilab code Exa 2.16 Overall heat transfer coecient
1 clear all;
2 clc;
3
5 // Chapter 2
7
8 // Example 2. 16
9 // Page 75
11
12 // The o r e t i c a l Probl em
13
14 printf( \n\n Thi s i s a The o r e t i c a l Problem , does not
i nv o l v e any mat hemat i cal comput at i on . );
33
15 //END
34
Chapter 3
Thermal Radiation
Scilab code Exa 3.1 Monochromatic emissive power
1 clear all;
2 clc;
3
5 // Chapter 3
6 // Thermal Radi at i on
7
8 // Example 3 . 1
9 // Page 114
10 printf(Example 3 . 1 , Page 114 \n\n);
11
12 T = 5779 ; // [ Temperature , i n Kel vi n ]
13 // From Wein s law , eqn 3 . 2 . 8
14 lambda_m = 0.00290/T ; // [m]
15 // Su b s t i t u t i ng t h i s val ue i n pl ank s law , we ge t
16 e = 2*(%pi)*0.596*(10^ -16) /(((0.5018*10^ -6) ^5)*(exp
(0.014387/0.00290) -1)) ; // [W/m2 m]
17
18 e_bl_max= e / 10^6 ;
19
20 printf( Val ue o f e mi s s i v i t y on sun s u r f a c e i s %f W/m
35
2 um \n,e_bl_max); // [W/m2 um]
21
22 e_earth = e_bl_max *((0.695*10^6) /(1.496*10^8))^2 ;
23
24 printf(The val ue o f e mmi s s i vi t y on e a r t hs s u r f a c e
i s %f W/m2 um, e_earth)
Scilab code Exa 3.2 Heat ux
1 clear all;
2 clc;
3
5 // Chapter 3
7
8 // Example 3 . 2
9 // Page 115
11
12 // Heat e mi s s i o n
13 Stefan_constt = 5.67*10^( -8); // (W/m 2 .K4)
14 T = 1500; // t e mpe r at ur e i s
i n k e l v i n s
15 eb = (Stefan_constt)*(T^(4)); // ener gy
r a di a t e d by bl ackbody
16 // e mi s s i o n i n 0 . 3um t o 1um
17 e = 0.9; // e mi s s i v i t y
18 lamda1 = 1; // wavel engt h i s i n um
19 lamda2 = 0.3; // wavel engt h i s i n um
20 D0_1 =0.5*(0.01972+0.00779); //From t a bl e 3 . 1
page 114
21 D0_2 =0; //From t a bl e 3 . 1 page
114
22 q = e*(D0_1 -D0_2)*Stefan_constt*T^(4);// i n W/m2
36
23 printf(\n wavel engt h temp = %d um K ,1*1500);
24 printf(\n wavel engt h temp at 0 . 3um = %d um K
,0.3*1500);
25 printf(\n\n Requi r ed heat f l ux , q = %d W/m2 ,q);
Scilab code Exa 3.3 Absorbed radiant ux and absorptivity and reectiv-
ity
1 clear all;
2 clc;
3
5 // Chapter 3
7
8 // Example 3 . 3
9 // Page 119
11
12
13 a0_2 =1; // a b s o r p t i v i t y
14 a2_4 =1; // a b s o r p t i v i t y
15 a4_6 =0.5; // a b s o r p t i v i t y
16 a6_8 =0.5; // a b s o r p t i v i t y
17 a8_ =0; // a b s o r p t i v i t y
18 H0_2 =0; // I r r a d i a t i o n i n W/m2 um
22 H8_ =750; // I r r a d i a t i o n i n W/m2 um
23 Absorbed_radiant_flux =1*0*(2 -0) +1*750*(4 -2)
+0.5*750*(8 -4) +0;
24 H = 750*(8 -2); // I nc i de nt f l u x
25 a = Absorbed_radiant_flux/H;
26 p = 1-a; // Si nc e t he s u r f a c e i s opaque
37
27 printf(\n Absorbed r a di a nt f l u x = %d W/m2 ,
Absorbed_radiant_flux);
28 printf(\n I nc i de nt f l u x = %d W/m2 ,H);
29 printf(\n Abs o r pt i vi t y = %. 3 f ,a);
30 printf(\n Si nc e t he s u r f a c e i s opaque r e f l e c t i v i t y
= %. 3 f ,p);
Scilab code Exa 3.4.a Total intensity in normal direction
1 clear all;
2 clc;
3
5 // Chapter 3
7
8 // Example 3 . 4 ( a )
9 // Page 123
10 printf(Example 3 . 4 ( a ) , Page 123 \n\n)
11
12
13 e = 0.08; // e mi s s i v i t y
14 T = 800; // t emper at ur e , [ K]
15
16 Stefan_constt = 5.67*10^( -8); // [W/m 2 .K 4 ]
17 // From St e f an Bol tzmann law , e quat i on 3 . 2 . 1 0
18 q = e*Stefan_constt*T^4; // [W/m 2 ]
19 printf(\n Energy e mi t t e d = %. 1 f W/m2 ,q);
20
21 // ( a )
23 in = (q/(%pi));
24 printf(\n Energy e mi t t e d normal t o t he s u r f a c e = %
. 1 f W/m2 s r ,in);
38
Scilab code Exa 3.4.b Ratio of radiant ux to the emissive power
1 clear all;
2 clc;
3
5 // Chapter 3
7
8 // Example 3 . 4 ( b)
9 // Page 123
10 printf(Example 3 . 4 ( b) , Page 123 \n\n)
11
12
13 e = 0.08; // e mi s s i v i t y
14 T = 800; // t emper at ur e , [ K]
15
16 Stefan_constt = 5.67*10^( -8); // [W/m 2 .K 4 ]
17 // From St e f an Bol tzmann law , e quat i on 3 . 2 . 1 0
18 q = e*Stefan_constt*T^4; // [W/m 2 ]
19 in = (q/(%pi));
20
21 // ( b)
22 // Radi ant f l u x e mi t t e d i n t he cone 0 <= pz i <= 50
degr ee , 0 <= t he t a <= 2 pi
23 q_cone =2*( %pi)*in*(-cos (100*( %pi /180))+cos(0))/4;
24
25 printf(\n Radi ant f l u x e mi t t e d i n t he cone =%. 1 f W/
m2 ,q_cone);
26
27 Ratio = q_cone/q;
28 printf(\n Rat i o = %. 3 f ,Ratio);
39
Scilab code Exa 3.5 Rate of incident radiation
1 clear all;
2 clc;
3
5 // Chapter 3
7
8 // Example 3 . 5
9 // Page 124
11
12 l1 = 0.5 ; // wavel ength , [ um]
13 l2 = 1.5 ; // wavel ength , [ um]
14 l3 = 2.5 ; // wavel ength , [ um]
15 l4 = 3.5 ; // wavel ength , [ um]
16 H1 = 2500 ; // [W/m2 um]
17 H2 = 4000 ; // [W/m2 um]
18 H3 = 2500 ; // [W/m2 um]
19
20 // Si nc e t he i r r i d i a t i o n i s d i f f u s e , t he s p e c t r a l
i n t e n s i t y i s gi ve n by eqn 3 . 4 . 1 4 and 3 . 4 . 8
21 // I n t e g r a t i n g i l ambda over t he d i r e c t i o n s o f t he
s p e c i f i e d s o l i d angl e and us i ng f i g 3. 12
22
23
24 flux = 3/4*[H1*(l2 -l1)+H2*(l3 -l2)+H3*(l4-l3)];
25 printf( Rate at whi ch r a d i a t i o n i s i n c i d e n t on t he
s u r f a c e = %f W/m2 ,flux);
Scilab code Exa 3.6 Shape factor F12
40
1 clear all;
2 clc;
3
5 // Chapter 3
7
8 // Example 3 . 6
9 // Page 132
11
12 // Thi s i s a t h e o r e t i c a l probl em wi th no nume r i c al
data
13 printf( Thi s i s a t h e o r e t i c a l probl em wi th no
nume r i c al data \n);
14
15 // Cons i de r i ng an e l e me nt ar y r i ng dA2 o f wi dth dr at
an a r b i t a r y r a di us r , we have
16 // r = htanB1
17 // dA2 = 2%pi r dr
18 // dA2 = 2%pi ( h2) tan ( B1) s e c 2( B1) dB1
19 // B2 = B1 , s i n c e s u r f a c e s at e p a r a l l e l , and
20 // L = h/ c os ( B1)
21 // Su b s t i t u t i ng i n eqn 3 . 6 . 7
22 // F12 = s i n 2( a )
23
24
25 printf( Cons i de r i ng an e l e me nt ar y r i ng dA2 o f wi dth
dr at an a r b i t a r y r a di us r , we have \n);
26 printf( r = htanB1 \n);
27 printf(dA2 = 2 pi r dr \n);
28 printf(dA2 = 2 pi ( h2) tan ( B1) s e c 2( B1) dB1 \n);
29 printf(B2 = B1 , s i n c e s u r f a c e s at e p a r a l l e l , and \n
);
30 printf(L = h/ c os ( B1) \n);
31 printf( Su b s t i t u t i n g i n eqn 3 . 6 . 7 \n);
32 printf(F12 = s i n 2( a ) \n);
41
Scilab code Exa 3.7 Shape factor
1 clear all;
2 clc;
3
5 // Chapter 3
7
8 // Example 3 . 7
9 // Page 134
11
12 // Thi s i s a t h e o r e t i c a l probl em wi th no nume r i c al
data
13 printf( Thi s i s a t h e o r e t i c a l probl em wi th no
nume r i c al data \n);
14
15
16 // Cons i de r i ng an e l e me nt ar y c i r c u l a r r i ng on t he
s u r f a c e o f t he s pher e s s u r f a c e at any a r b i t a r y
a ng l r B,
17 // we have B1 = B, B2 = 0 , L = R and dA 2 = 2%pi (R
2) ( s i n (B) ) dB
18 // The r e f or e , f rom e quat i on 3 . 6 . 7
19 // F12 = s i n 2( a )
20
21 printf( Cons i de r i ng an e l e me nt ar y c i r c u l a r r i ng on
t he s u r f a c e o f t he s phe r e s u r f a c e at any a r b i t a r y
a ng l r B \n);
22 printf(we have B1 = B, B2 = 0 , L = R and dA 2 = 2
pi (R2) ( s i n (B) ) dB \n);
23 printf( The r e f or e , f rom e quat i on 3 . 6 . 7 \n);
24 printf(F12 = s i n 2( a ) );
42
Scilab code Exa 3.8 Shape factor F12
1 clear all;
2 clc;
3
5 // Chapter 3
7
8 // Example 3 . 8
9 // Page 135
11
12 // From eqn 3 . 7 . 5 or f i g 3. 19
13 F65 = 0.22;
14 F64 = 0.16;
15 F35 = 0.32;
16 F34 = 0.27;
17 A1 = 3; // [m 2 ]
18 A3 = 3; // [m 2 ]
19 A6 = 6; // [m 2 ]
20
21 // Usi ng a d d i t i v e and r e c i p r o c a l r e l a t i o n s
22 // We have F12 = F16 F13
23
24 F61 = F65 - F64 ;
25 F31 = F35 - F34 ;
26
27 F16 = A6/A1*F61 ;
28 F13 = A3/A1*F31 ;
29
30 F12 = F16 - F13;
31
32 printf(F 12 = %f ,F12);
43
Scilab code Exa 3.9 Shape factor
1 clear all;
2 clc;
3
5 // Chapter 3
7
8 // Example 3 . 9
9 // Page 136
11
12 // Thi s i s a t h e o r e t i c a l probl em , does not i nv o l v e
any nume r i c al comput at i on
13 printf( Thi s i s a t h e o r e t i c a l probl em , does not
i nv o l v e any nume r i c al comput at i on \n);
14 // Denot i ng ar e a o f c o n i c a l s u r f a c e by A1
15 // Cons i de r i ng an i magi nar y f l a t s u r f a c e A2 c l o s i n g
t he c o n i c a l c a v i t y
16
17 F22 = 0 ; // Fl at s u r f a c e
18
19 // f rom eqn 3 . 7 . 2 , we have F11 + F12 = 1 and F22 +
F21 = 1
20 F21 = 1 - F22 ;
21
22 // F12 = A2/A1F21 ;
23 // F11 = 1 F12 ;
24 // F11 = 1 s i n ( a )
Scilab code Exa 3.10 Net radiative heat transfer
44
1 clear all;
2 clc;
3
5 // Chapter 3
7
8 // Example 3. 10
9 // Page 138
11
12 sigma = 5.670*10^ -8 ;
13 T1 = 473 ; // [ K]
14 T2 = 373 ; // [ K]
15 A1 = 1*2 ; // area , [m 2 ]
16 X = 0.25;
17 Y = 0.5 ;
18 // From eqn 3 . 7 . 4
19 F12 = (2/( %pi*X*Y))*[log ((((1+X^2) *(1+Y^2))/(1+X^2+Y
^2))^(1/2)) + Y*((1+X^2) ^(1/2))*atan(Y/((1+X^2)
^(1/2))) + X*((1+Y^2) ^(1/2))*atan(X/((1+Y^2)
^(1/2))) - Y*atan(Y) - X*atan(X) ];
20
21
22 q1 = sigma*A1*(T1^4-T2^4)*[(1-F12 ^2)/(2-2* F12)];
23
24 printf(Net r a d i a t i v e heat t r a n s f e r f rom t he s u r f a c e
= %f W \n,q1);
Scilab code Exa 3.11 steady state heat ux
1 clear all;
2 clc;
3
45
5 // Chapter 3
7
8 // Example 3. 11
9 // Page 141
11
12 // Al l modes o f heat t r a n s f e r ar e i nvo l ve d
13 // l e t s t e ady s t a t e heat f l u x f l o wi ng t hr ough t he
c ompos i t e s l a b be ( q/a )
14 h1 = 20; // [W/m2 K]
15 w1 = 0.2; // [m]
16 k1 = 1; // [W/m K]
17 e1 = 0.5; // e mmi s i vi t y at s u r f c e 1
18 e2 = 0.4; // e mmi s i vi t y at s u r f c e 2
19 w2 = 0.3; // [m]
20 k2 = 0.5; // [W/m K]
21 h2 = 10; // [W/m2 K]
22 T1 = 473; // [ Kel vi n ]
23 T2 = 273+40; // [ Kel vi n ]
24 stefan_cnst = 5.67e-08; // [W/m2 K 4 ]
25
26 // For r e s i s t a n c e s 1 and 2
27 function[f]= temperature(T)
28 f(1) = (T1-T(1))/(1/h1 + w1/k1) - (T(2) - T2)/(
w2/k2 + 1/h2);
29 f(2) = stefan_cnst *(T(1)^4 - T(2)^4) /(1/e1 + 1/
e2 -1) - (T(2) - T2)/(w2/k2 + 1/h2);
30 funcprot (0);
31 endfunction
32
33 T = [10 10]; // assumed i n i t i a l va l ue s f o r f s o l v e
f unc t i o n
34 y = fsolve(T,temperature);
35
36 printf(\n Steady s t a t e heat f l u x q/A = %. 1 f W/m2
,(T1-y(1))/(1/h1 + w1/k1));
46
Scilab code Exa 3.12 Rate of heat loss
1 clear all;
2 clc;
3
5 // Chapter 3
7
8 // Example 3. 12
9 // Page 145
11
12 D = 0.02 ; // [m]
13 T1 = 1000+273 ; // [ K]
14 T2 = 27+273 ; // [ K]
15 s = 5.670*10^ -8 ; // s t e f a n s c ons t ant
16 // Assumi ng t he openi ng i s c l o s e d by an i magi nar y
s u r f a c e at t empe r at ur e T1
17 // Usi ng e quat i on 3 . 1 0 . 3 , we ge t
18 q = s*1*%pi *((D/2)^2)*(T1^4-T2^4); // [W]
19
20 printf( Rate at whi ch heat i s l o s t by r a d i a t i o n = %f
W,q);
Scilab code Exa 3.13 Rate of nitrogen evaporation
1 clear all;
2 clc;
3
5 // Chapter 3
47
7
8 // Example 3. 13
9 // Page 146
11
12 D = 0.32 ; // [m]
13 D_s = 0.36 ; // [m]
14 e = 0.02 ; // e mi s s i v i t y
15 l = 201 ; // [ kJ/kg ]
16 rho = 800 ; // [ kg/m 3 ]
17 s = 5.670*10^ -8 ;
18
19 T2 = 303 ; // [ K]
20 T1 = 77 ; // [ K]
21
22 // From e quat i on 3 . 1 0 . 1
23 q1 = s*4*%pi *((D/2)^2)*(T1^4-T2^4) /[1/e+((D/D_s)^2)
*(1/e-1)]; // [W]
24
25 evap = abs(q1)*3600*24/(l*1000); // [ kg/day ]
26 mass = 4/3* %pi *((D/2)^3)*rho;
27 boiloff = evap/mass *100 ; // pe r c e nt
28
29 T_drop = (abs(q1))/(4* %pi *((D/2)^2))*(1/100); // [ C]
30
31 printf( Rate at whi ch ni t r o g e n e va po r a t e s = %f kg/
day \n,evap)
32 printf( Boi l o f f r a t e = %f pe r c e nt \n,boiloff);
33 printf( Temperature drop between l i q u i d Ni t r oge n and
i nne r s u r f a c e = %f C,T_drop);
Scilab code Exa 3.14 Rate of energy loss from satellite
1 clear all;
48
2 clc;
3
5 // Chapter 3
7
8 // Example 3. 14
9 // Page 147
11
12 D = 1 ; // [m]
13 r = 6250 ; // [ km]
14 D_surf = 300 ; // [ km]
15 s = 5.670*10^ -8;
16 e = 0.3 ;
17 Tc = -18+273 ; // [ K]
18 T_surf = 27+273 ; // [ K]
19
20 // Rate o f e mi s s i no o f r a di a nt ener gy f rom t he two
f a c e s o f s a t e l l i t e d i s c
21 r_emission = 2*e*%pi*((D/2)^2)*s*Tc^4; // [W]
22
23 // A2F21 = A1F12
24 sina = (r/(r+D_surf));
25 F12 = sina ^2;
26
27 // Rate at whi ch t he s a t e l l i t e r e c e i v e s and abs or bs
ener gy comi ng f rom e ar t h
28 r_receive = e*s*(%pi *((D/2)^2))*F12*T_surf ^4; // [W]
29
30 r_loss = r_emission - r_receive; // [W]
31
32 printf(Net Rate at whi ch ener gy i s l e a v i ng t he
s a t e l l i t e = %f W,r_loss);
49
Scilab code Exa 3.15 Net radiative heat transfer
1 clear all;
2 clc;
3
5 // Chapter 3
7
8 // Example 3. 15
9 // Page 151
11
12 // From exampl e 3. 10
13 F12 = 0.0363;
14 F11 = 0;
15 F13 = 1-F11 -F12;
16 // S i mi l a r l y
17 F21 = 0.0363;
18 F22 = 0;
19 F23 = 0.9637;
20
21 // Now, F31 = A1/A3F13
22 F31 = 2/24* F13;
24 F32 = F31;
25 F33 = 1-F31 -F32;
26
27 // Su b s t i t u t i ng i nt o e quat i on 3 . 1 1 . 6 , 3 . 1 1 . 7 ,
3 . 1 1 . 8 , we have f ( 1 ) , f ( 2 ) , f ( 3 )
28
29 function[f]=flux(B)
30 f(1)= B(1) - 0.4*0.0363*B(2) - 0.4*0.9637*B(3) -
0.6*(473^4) *(5.670*10^ -8);
31 f(2)= -0.4*0.0363*B(1) + B(2) - 0.4*0.9637*B(3)
- 0.6*(5.670*10^ -8) *(373^4);
32 f(3)= 0.0803*B(1) + 0.0803*B(2) - 0.1606*B(3);
33 funcprot (0);
50
34 endfunction
35
36 B = [0 0 0];
37 y = fsolve(B,flux);
38 printf(\n B1 = %. 1 f W/m2 ,y(1));
39 printf(\n B2 = %. 1 f W/m2 ,y(2));
40 printf(\n B3 = %. 1 f W/m2 \n,y(3));
41
43 H1 = 0.0363*y(2) + 0.9637*y(3) ; // [W/m 2 ]
44 // and
45 q1 = 2*(y(1) - H1) ; // [W]
46
47 printf(Net r a d i a t i v e heat t r a n s f e r = %f W,q1);
51
Chapter 4
Principles of Fluid Flow
Scilab code Exa 4.1 Pressure drop in smooth pipe
1 clear;
2 clc;
3
5 // Chapter 4
6 // Pr i n c i p l e s o f Fl ui d Flow
7
8 // Example 4 . 1
9 // Page 172
11
12 L = 3 ; // Length , [m]
13 D = 0.01 ; // ID , [m]
14 V = 0.2 ; // Average Ve l oc i t y , [m/ s ]
15
16 // From Tabl e A. 1 at 10 de gr e e C
17 rho =999.7 ; // [ kg/m 3 ]
18 v=1.306 * 10^-6 ; // [m2/ s ]
19
20 Re_D =0.2*0.01/(1.306*10^ -6) ;
21
52
22 // t h i s val ue i s l e s s than t he t r a n s i t i o n Reynol ds
number 2300.
23 // Hence f l ow i s l ami nar . From eqn 4 . 4 . 1 9
24 f = 16/ Re_D;
25
26 // f rom eqn 4 . 4 . 1 7
27 delta_p = 4*f*(L/D)*(rho*V^2) /2;
28
29 // s i n c e f l ow i s l ami nar
30 V_max = 2*V;
31
32 printf( Pr e s s ur e drop i s %f Pa \n,delta_p);
33 printf( Maximum v e l o c i t y i s %f m/ s ,V_max);
Scilab code Exa 4.2.a Pressure drop and maximum velocity calculation
1 clear;
2 clc;
3
5 // Chapter 4
7
8 // Example 4 . 2 ( a )
9 // Page 180
10 printf(Example 4 . 2 ( a ) , Page 180 \n\n)
11
12 L = 3 ; // [m]
13 D = 0.01 ; // [m]
14 V = 0.2 ; // [m/ s ]
15
16 // ( a )
17 printf( ( a ) I f t he t empe r at ur e o f wat er i s i n c r e a s e d
t o 80 de gr e e C \n);
18
53
19
20 // Pr o p e r t i e s o f wat er at 80 de gr e e C
21 rho = 971.8 ; // [ kg/m 3 ]
22 v = 0.365 * 10^-6 ; // [m2/ s ]
23
24 Re_D = D*V/v;
25
26 // f l ow i s t u r b i l e n t , s o f rom eqn 4 . 6 . 4 a
27
28 f=0.079*( Re_D)^( -0.25);
29 delta_p = (4*f*L*rho*V^2)/(D*2); // [ Pa ]
30 printf( Pr e s s ur e drop i s %f Pa \n,delta_p);
31
32 // f rom eqn 4 . 4 . 1 6
33
34 // x = ( T w/p) 0. 5 = ( ( f /2) 0 . 5 ) V ;
35 x = ((f/2) ^0.5)*V ;
36 y_plus = 0.005*x/(0.365*10^ -6);
37
38 // f rom eqn 4 . 6 . 1 c & 4 . 6 . 2
39
40 V_max = x*(2.5* log(y_plus) + 5.5) ; // [m/ s ]
41 ratio = V_max/V;
42 printf(V max = %f m/ s \n,V_max);
43 printf(V max/ V bar = %f \n\n,ratio);
Scilab code Exa 4.2.b Pressure drop and maximum velocity calculation
1 clear;
2 clc;
3
5 // Chapter 4
7
54
8 // Example 4 . 2 ( b)
9 // Page 180
10 printf(Example 4 . 2 ( b) , Page 180 \n\n)
11
12 L = 3 ; // [m]
13 D = 0.01 ; // [m]
14 V = 0.2 ; // [m/ s ]
15
16 // ( b)
17
18 V1=0.7;
19 v1 = 1.306 * 10^-6 ; // [m2/ s ]
20
21 printf( ( b) I f t he v e l o c i t y i s i n c r e a s e d t o 0 . 7 \n)
;
22 // i f v e l o c i t y o f wat er i s 0 . 7 m/ s
23 V1=0.7; // [m/ s ]
24 Re_D1=V1*D/(1.306*10^ -6);
25 printf( Reynol ds no i s %f \n,Re_D1);
26
27 // f l ow i s agai n t ur bul e nt
28 f1 = 0.079*( Re_D1)^( -0.25);
29
30 delta_p1 = (4*f1*L*999.7*0.7^2) /(0.01*2); // [ Pa ]
31 printf( Pr e s s ur e drop i s %f Pa \n,delta_p1);
32
33 // x1 = ( T w/p) 0. 5 = ( ( f 1 /2) 0 . 5 ) V ;
34 x1 = ((f1/2) ^0.5)*V1 ;
35
36 y1_plus = 0.005* x1/(v1);
37 printf(y+ at c e nt r e l i n e = %f \n,y1_plus);
38
39 V_max1 = x1 *(2.5* log(y1_plus) + 5.5) ; // [m/ s ]
40 printf(V max i s %f m/ s \n,V_max1);
41
42 ratio1 = V_max1/V1;
43 printf(Vmax/Vbar = %f ,ratio1);
55
Scilab code Exa 4.3 Pressure drop and power needed
1 clear;
2 clc;
3
5 // Chapter 4
7
8 // Example 4 . 3
9 // Page 181
11 P = 80 * 10^3 ; // [ Pa ]
12 L = 10 ; // [m]
13 V_bar = 1.9 ; // [m/ s ]
14 l = 0.25 ; // [m]
15 b = 0.15 ; // [m]
16
17 // Ful l y de ve l ope d f l ow
18
19 // From Tabl e A. 2 , f o r a i r at ! atm pr e s s ur e and 25
de gr e e C
20 rho = 1.185 ; // [ kg/m 3 ]
21 mew = 18.35 * 10^-6 ; // [ kg/m s ]
22
23 // At 80 kPa and 25 de gr e e C
24 rho1 = rho *(80/101.3) ; // [ kg/m 3 ]
25
26 // For gi ve n duct r =(b/a )
27 r = b/l;
28
29 D_e = (4*l/2*b/2)/(l/2 + b/2); // [m]
30
31 // From eqn 4 . 6 . 7
56
32
33 D_l = [2/3 + 11/24*0.6*(2 -0.6) ]*D_e ; // [m]
34
35 // Reynol ds no based on D l
36
37 Re = rho1*D_l*V_bar/mew;
38 printf( Reynol ds no = %f \n,Re);
39
40 f = 0.079*( Re^ -0.25) ;
41 printf( f = %f \n,f);
42
43 // From eqn 4 . 4 . 1 7
44
45 delta_P = 4*f*(L/D_l)*(rho1*(V_bar ^2) /2);
46 printf( Pr e s s ur e drop = %f Pa \n,delta_P);
47
48 power = delta_P *( V_bar*l*b)
49 printf(Power r e q ui r e d = %f W,power);
Scilab code Exa 4.4 Thickness of velocity boundary layer
1 clear;
2 clc;
3
5 // Chapter 4
7
8 // Example 4 . 4
9 // Page 189
11
12 l = 2 ; // [m]
13 b = 1 ; // [m]
14 V = 1 ; // [m/ s ]
57
15
16 // From t a bl e A. 2
17
18 rho = 1.060 ; // [ kg/m 3 ]
19 v = 18.97 * 10^-6 ; // [m2/ s ]
20
21 // At x = 1 . 5m
22 x = 1.5 ; // [m]
23 Re = V*x/v; // Reynol ds number
24
25 // From eqn 4 . 8 . 1 2
26
27 d = 5*x/(Re ^(1/2))*1000 ; // [mm]
28 printf( Thi c kne s s o f Boundary l a y e r at x = 1 . 5 i s %f
mm \n,d)
29
30 Re_l = V*l/v;
31
32 // From eqn 4 . 8 . 1 9 and 4 . 8 . 1 6
33
34 c_f = 1.328* Re_l ^ -(1/2); // drag c o e f f i c i e n t
35 printf(Drag Co e f f i c i e n t c f = %f \n,c_f);
36
37 F_d = 0.00409*(1/2)*rho *(2*l*b)*1^2;
38 printf(Drag For ce F D = %f N,F_d);
Scilab code Exa 4.5 Drag coecient and drag force
1 clear;
2 clc;
3
5 // Chapter 4
7
58
8 // Example 4 . 5
9 // Page 195
11
12 l = 2 ; // [m]
13 v = 4 ; // [m/ s ]
14
15 // From Tabl e A. 2
16
17 mew = 18.1*10^ -6; // [ N s /m 2 ]
18 rho = 1.205*1.5; // [ kg/m 3 ]
19
20 Re_l = rho*v*l/mew;
21 // Boundary l a y e r i s pa r t l y l ami nar and pa r t l y
t ur bul e nt , we s h a l l us e eqn 4 . 1 0 . 4
22 Cf = 0.074*(7.989*10^5) ^( -0.2) - 1050/ Re_l ;
23 printf(Drag c o e f f i c i e e n t i s %f \n,Cf)
24
25 D_f= Cf *1/2* rho*l*v^2;
26 printf(Drag f o r c e per meter wi dth = %f N \n,D_f);
27
28 // f rom eqn 4 . 1 0 . 1
29
30 x = 3*10^5 * (18.1*10^ -6) /(1.808*4);
31 printf( Val ue o f x c i s %f m,x);
59
Chapter 5
Heat Transfer by Forced
Convection
Scilab code Exa 5.1.a Local heat transfer coecient
1 clear;
2 clc;
3
5 // Chapter 5
6 // Heat Tr ans f e r by Forced Convect i on
7
8
9 // Example 5 . 1 ( a )
10 // Page 209
11 printf(Example 5 . 1 ( a ) \n\n)
12
13 D = 0.015 ; // [m]
14 Q = 0.05 ; // [m3/h ]
15 H = 1000 ; // [W/m 2 ]
16 T_b = 40 ; // [ de gr e e C]
17
18 // From t a bl e A. 1 , p r o p e r t i e s at 40 de gr e e C
19 k = 0.634 ; // [W/m K]
60
20 v = 0.659*10^ -6 ; // [m2/ s ]
21
22 V_bar = 4*Q/((%pi)*D^2);
23
24 Re_D = V_bar*D/v;
25
26 // The r e f or e , Lami nar Flow , f rom eqn 5 . 2 . 8
27
28 h = 4.364*k/D; // [W/m2 K]
29
30 printf( ( a ) Loc al heat t r a n s f e r c o e f f i c i e n t i s %f W/
m2 K \n,h);
Scilab code Exa 5.1.b Wall temperature
1 clear;
2 clc;
3
5 // Chapter 5
7
8
9 // Example 5 . 1 ( b)
10 // Page 209
11 printf(Example 5 . 1 ( b) \n\n)
12
13 D = 0.015 ; // [m]
14 Q = 0.05 ; // [m3/h ]
15 H = 1000 ; // [W/m 2 ]
16 T_b = 40 ; // [ de gr e e C]
17
18 // From t a bl e A. 1 , p r o p e r t i e s at 40 de gr e e C
19 k = 0.634 ; // [W/m K]
20 v = 0.659*10^ -6 ; // [m2/ s ]
61
21
22 V_bar = 4*Q/((%pi)*D^2);
23
24 Re_D = V_bar*D/v;
25
26 // The r e f or e , Lami nar Flow , f rom eqn 5 . 2 . 8
27
28 h = 4.364*k/D;
29
30 // From t he d e f i n i t i o n o f h i n eqn 5 . 2 . 3 , t he l o c a l
wal t o bul k mean t empe r at ur e d i f f e r e n c e i s gi ve n
by
31
32 T_w = H/h + T_b;
33
34 printf( ( b) Wal l Temperature Tw = %f de gr e e C,T_w);
Scilab code Exa 5.2 ratio of thermal entrance length to entrance length
1 clear;
2 clc;
3
5 // Chapter 5
7
8
9 // Example 5 . 2
10 // Page 213
12
13 // From eqn 5 . 2 . 1 2 and 4 . 4 . 2 0
14 // Let r = Lth/Le
15 // r = 0. 04305 Pr / 0 . 0 5 7 5 ;
16
62
17 function[T]=r(Pr)
18 T = 0.04305* Pr /0.0575
19 endfunction
20
21 // For Pr = 0. 01
22 r1 = r(0.01);
23 // For Pr = 0 . 1
24 r2 = r(1);
25 // For Pr = 100
26 r3 = r(100);
27
28 printf(Lth/Le at Pr = 0. 01 i s %f \n,r1);
29 printf(Lth/Le at Pr = 1 i s %f \n,r2);
30 printf(Lth/Le at Pr = 100 i s %f ,r3);
Scilab code Exa 5.3.i Length of tube
1 clear;
2 clc;
3
5 // Chapter 5
7
8
9 // Example 5 . 3 ( i )
10 // Page 215
11 printf(Example 5 . 3 ( i ) , Page 215 \n\n)
12
13 D = 0.015 ; // [m]
14 V = 1 ; // [m/ s ]
15 Tw = 90 ; // [ de gr e e C]
16 Tmi = 50 ; // [ de gr e e C]
17 Tmo = 65 ; // [ de gr e e C]
18
63
19 // ( i )
20
22 k = 0.656 ; // [W/m K]
23 rho = 984.4 ; // [ kg/m 3 ]
24 v = 0.497 * 10^-6 ; // [m2/ s ]
25 Cp = 4178 ; // [ J/kg K]
26 Pr = 3.12 ;
27 rho_in = 988.1 ; // [ kg/m 3 ]
28
29 m_dot = %pi*(D^2)*rho_in*V/4 ; // [ kg/ s ]
30
31 Re = 4* m_dot/(%pi*D*rho*v) ;
32
33 // Usi ng eqn 5 . 3 . 2 and 4 . 6 . 4 a
34 f = 0.079*( Re)^-0.25 ;
35
36 Nu = (f/2)*(Re -1000)*Pr /[1+12.7*(f/2) ^(1/2) *((Pr
^(2/3)) -1)];
37 h = Nu*k/D; // [W/m2 K]
38
39 // From t he ener gy equat i on , e x t r a c t i n g t he val ue o f
L
40 L = m_dot*Cp*(Tmo -Tmi)*[log((Tw-Tmi)/(Tw-Tmo))]/[((
Tw -Tmi) -(Tw-Tmo))*h*D*%pi]; // [m]
41
42 printf(The l e ng t h o f tube i f t he e x i t wat er
t e mpe r at ur e i s 65 de gr e e C = %f m\n,L);
Scilab code Exa 5.3.ii Exit water temperature
1 clear;
2 clc;
3
64
5 // Chapter 5
7
8
9 // Example 5 . 3 ( i )
10 // Page 215
11 printf(Example 5 . 3 ( i i ) , Page 215 \n\n)
12
13 D = 0.015 ; // [m]
14 V = 1 ; // [m/ s ]
15 Tw = 90 ; // [ de gr e e C]
16 Tmi = 50 ; // [ de gr e e C]
17 Tmo = 65 ; // [ de gr e e C]
18
20 k = 0.656 ; // [W/m K]
21 rho = 984.4 ; // [ kg/m 3 ]
22 v = 0.497 * 10^-6 ; // [m2/ s ]
23 Cp = 4178 ; // [ J/kg K]
24 Pr = 3.12 ;
25 rho_in = 988.1 ; // [ kg/m 3 ]
26
27 m_dot = %pi*(D^2)*rho_in*V/4 ; // [ kg/ s ]
28
29 Re = 4* m_dot/(%pi*D*rho*v) ;
30
31 // Usi ng eqn 5 . 3 . 2 and 4 . 6 . 4 a
32 f = 0.079*( Re)^-0.25 ;
33
34 Nu = (f/2)*(Re -1000)*Pr /[1+12.7*(f/2) ^(1/2) *((Pr
^(2/3)) -1)];
35 h = Nu*k/D; // [W/m2 K]
36
37 // From t he ener gy equat i on , e x t r a c t i n g t he val ue o f
L
38 L = m_dot*Cp*(Tmo -Tmi)*[log((Tw-Tmi)/(Tw-Tmo))]/[((
Tw -Tmi) -(Tw-Tmo))*h*D*%pi]; // [m]
39
65
40 // ( i i )
41 printf(\ nTr i a l and e r r o r method \n);
42
43 // Tr i a l 1
44 printf( Tr i a l 1\n);
45 printf(Assumed val ue o f Tmo = 70 de gr e e C\n);
46 T_mo = 70 ; // [ de gr e e C]
47 T_b = 60 ; // [ de gr e e C]
48
49 k1 = 0.659 ; // [W/m K]
50 rho1 = 983.2 ; // [ kg/m 3 ]
51 v1 = 0.478 * 10^-6 ; // [m2/ s ]
52 Cp1 = 4179 ; // [ J/kg K]
53 Pr1 = 2.98 ;
54
55 Re1 = 4*m_dot/(%pi*D*rho1*v1);
56
57 // From Bl a s i us eqn ( 4 . 6 . 4 a ) , we ge t
58 f1 = 0.005928;
59
60 // Su b s t i t u t i ng t h i s val ue i nt o t he Gn i e l i n s k i Eqn
61 Nu_d = 154.97;
62 h = Nu_d*k1/D ; // [W/m2 K]
63
64 // f rom eqn 5 . 3 . 3 , we ge t
65 Tmo1 = 73.4 ; // [ de gr e e C]
66 printf( Val ue o f Tmo obt ai ne d = 73. 4 de gr e e C\n);
67
68 // Tr i a l 2
69 printf( Tr i a l 2\n);
70 printf(Assume Tmo = 73. 4 de gr e e C\n);
71 printf( Val ue o f Tmo obt ai ne d = 73. 6 de gr e e C whi ch
i s i n r e a s o na bl y c l o s e agr eement wi th assumed
val ue . \ n)
66
Scilab code Exa 5.4 Length of tube over which temperature rise occurs
1 clear;
2 clc;
3
5 // Chapter 5
7
8
9 // Example 5 . 4
10 // Page 219
12
13 D_i = 0.05 ; // [m]
14 m = 300 ; // [ kg/min ]
15 m1 = m/60 ; // [ kg/ s e c ]
16 rho = 846.7 ; // [ kg/m 3 ]
17 k = 68.34 ; // [W/m K]
18 c = 1274; // [ J/kg K]
19 v = 0.2937*10^ -6 ; // [m2/ s ]
20 Pr = 0.00468 ;
21
22 Re_D = 4*m1/(%pi*D_i*rho*v);
23
24 // Assumi ng both t e mpe r at ur e and v e l o c i t y p r o f i l e
ar e f u l l y de ve l ope d over t he l e ng t h o f tube
25 // us i ng eqn 5 . 3 . 6
26 Nu_D = 6.3 + 0.0167*( Re_D ^0.85) *(Pr ^0.93);
27
28 h = Nu_D*k/D_i;
29
30 // Equat i ng t he heat t r a n s f e r r e d t hr ough t he wal l o f
t he tube t o t he change o f e nt hal py pf sodi um
31 L = 300/60*1274*(500 -400) /(h*%pi*D_i *30)
32
33 printf( Length o f tube over whi ch t he t e mpe r at ur e
r i s e o c c ur s = %f m,L)
67
Scilab code Exa 5.5 Rate of heat transfer to the plate
1 clear;
2 clc;
3
5 // Chapter 5
7
8
9 // Example 5 . 5
10 // Page 231
11 printf(Example 5 . 5 , Page 231 \n)
12
13 V = 15 ; // [m/ s ]
14 s=0.2 ; // [m]
15 T_m = (20+60) /2; // [ de gr e e C]
16 // Pr o p e r t i e s at mean temp = 40 de gr e e C
17 v = 16.96*10^ -6; // [m2/ s ]
18 rho = 1.128 ; // [ kg/m 3 ]
19 k = 0.0276; // [W/m K]
20 Pr = 0.699;
21 A=s^2;
22 Re_L = V*0.2/v;
23 // Thi s i s l e s s than 3105 , hence t he boundary
l a y e r may be assumed t o be l ami nar over t he
e n t i r e l e ng t h .
24 // f rom eqn 4 . 8 . 1 9
25
26 Cf = 1.328/( Re_L)^0.5
27 Fd = 2*Cf*1/2* rho*A*V^2;
28
29 // From eqn 5 . 5 . 1 0
30 Nu_l = 0.664*( Pr ^(1/3))*(Re_L ^(1/2));
68
31
32 h = Nu_l*k/s;
33 // The r e f or e r a t e o f heat t r a n s f e r q i s
34 q = 2*A*h*(60 -20);// [W]
35
36 // With a t ur bul e nt boundary l a y e r f rom l e a di ng edge
, t he drag c o e f f i c i e n t i s gi ve n by eqn 4 . 1 0 . 4
37 Cf1 = 0.074*( Re_L)^( -0.2);
38 Fd1 = 2*Cf1 *1/2* rho*A*V^2; // [ N]
39
40 // f rom eqn 5 . 8 . 3 wi th C1 = 0
41 Nu_l1 = 0.0366*(0.699^(1/3))*(Re_L ^(0.8));
42
43 h1 = Nu_l1*k/s; // [W/m2 K]
44 q1 = 2*A*h1*(60 -20);
45
46 printf( For Lami nar Boundary Layer \n);
47 printf( Rate o f Heat t r a n s f e r = %f W\n,q);
48 printf(Drag f o r c e = %f N \n \n,Fd)
49
50 printf( For Tur bul ent Boundary Layer f rom t he
l e a di ng edge \n);
51 printf( Rate o f Heat t r a n s f e r = %f W\n,q1);
52 printf(Drag f o r c e = %f N\n,Fd1)
Scilab code Exa 5.6.i Heat transfer rate
1 clear;
2 clc;
3
5 // Chapter 5
7
8
69
9 // Example 5 . 6 ( i )
10 // Page 235
11 printf(Example 5 . 6 ( i ) , Page 235 \n\n)
12
13 D = 0.075 ; // [m]
14 V = 1.2 ; // [m/ s ]
15 T_air = 20 ; // [ de gr e e C]
16 T_surface = 100 ; // [ de gr e e C]
17 T_m = (T_air+T_surface)/2;
18
19 v = 18.97*10^ -6 ; // [m2/ s ]
20 k = 0.0290 ; // [W/m K]
21 Pr = 0.696 ;
22
23 Re_D = V*D/v;
24
25 Nu = 0.3 + [(0.62*( Re_D ^(1/2))*(Pr ^(1/3)))
/[(1+((0.4/ Pr)^(2/3)))^(1/4) ]]*([1+(( Re_D /282000)
^(5/8))]^(4/5)) ;
26
27 h = Nu*k/D ; // [W/m2 K]
28
29 flux = h*( T_surface - T_air); // [W/m 2 ]
30 q = flux*%pi*D*1; // [W/m]
31
32 printf( Heat t r a n s f e r r a t e per uni t l e ng t h = %f W/m\
n,q);
Scilab code Exa 5.6.ii Average wall tempeature
1 clear;
2 clc;
3
5 // Chapter 5
70
7
8
9 // Example 5 . 6 ( i i )
10 // Page 235
11 printf(Example 5 . 6 ( i i ) , Page 235 \n\n)
12
13 D = 0.075 ; // [m]
14 V = 1.2 ; // [m/ s ]
15 T_air = 20 ; // [ de gr e e C]
16 T_surface = 100 ; // [ de gr e e C]
17 T_m = (T_air+T_surface)/2;
18
19 v = 18.97*10^ -6 ; // [m2/ s ]
20 k = 0.0290 ; // [W/m K]
21 Pr = 0.696 ;
22
23 Re_D = V*D/v;
24 Nu = 0.3 + [(0.62*( Re_D ^0.5) *(Pr ^(1/3)))/[(1+((0.4/
Pr)^(2/3)))^(1/4) ]]*[1+( Re_D /282000) ^(5/8) ]^(5/8)
;
25 h = Nu*k/D ; // [W/m2 K]
26 flux = h*( T_surface - T_air); // [W/m 2 ]
27
28 // ( i i ) Usi ng Tr i a l and e r r o r method
29 T_avg = 1500/ flux*( T_surface - T_air);
30
31 T_assumd = 130 ; // [ de gr e e C]
32 Tm= 75 ; // [ de gr e e C]
33
34 v1 = 20.56*10^ -6 ; // [m2/ s ]
35 k1 = 0.0301 ; // [W/m K]
36 Pr1 = 0.693 ;
37
38 Re_D1 = V*D/v1;
39
40
41 // Usi ng eqn 5 . 9 . 8
71
42 Nu1 = 33.99;
43 h = Nu1*k1/D;
45 T_diff = 1500/h; // [ de gr e e C]
46 T_avg_calc = 129.9 ; // [ de gr e e C]
47 printf(Assumed ave r age wal l t e mper at ur e = %f de gr e e
C\n,T_assumd);
48 printf( Cal c ul at e d ave r age wal l Temperature = %f
de gr e e C\n,T_avg_calc);
49 printf(Hence , Average wal l Temperature = %f de gr e e C
,T_avg_calc);
Scilab code Exa 5.7.i Pressure drop
1 clear;
2 clc;
3
5 // Chapter 5
7
8
9 // Example 5 . 7 ( i )
10 // Page 241
11 printf(Example 5 . 7 ( i ) , Page 241 \n \n);
12
13 // Gi ven data
14 D = 0.0125 ; // [m]
15 ST = 1.5*D ;
16 SL = 1.5*D ;
17 V_inf = 2 ; // [m/ s ]
18
19 N = 5;
20 Tw = 70; // [ de gr e e C]
21 Tmi = 30; // [ de gr e e C]
72
22 L = 1; // [m]
23 // Pr o p e r t i e s o f a i r at 30 de gr e e C
24 rho = 1.165 ; // [ kg/m 3 ]
25 v = 16.00 *10^-6 ; // [m2/ s ]
26 Cp = 1.005 ; // [ kJ/kg K]
27 k = 0.0267 ; // [W/m K]
28 Pr = 0.701;
29
30 // From eqn 5 . 1 0 . 2
31 Vmax = ST/(SL-D)*V_inf ; // [m/ s ]
32 Re = Vmax*D/v ;
33
34 // From f i g 5. 15
35 f = 0.37/4;
36 // Al so , tube ar r angement i s s quar e
37 X = 1;
38 // From eqn 5 . 1 0 . 6
39 delta_P = 4*f*N*X*(rho*Vmax ^2)/2 ; // [ N/m 2 ]
40
41 printf( ( i ) Pr e s s ur e drop o f a i r a c r o s s t he bank i s
%f N/m2 \n,delta_P);
Scilab code Exa 5.7.ii Exit temperature of air
1 clear;
2 clc;
3
5 // Chapter 5
7
8
9 // Example 5 . 7 ( i i )
10 // Page 241
11 printf(Example 5 . 7 ( i i ) , Page 241 \n \n);
73
12
13 D = 0.0125 ; // [m]
14 ST = 1.5*D ;
15 SL = 1.5*D ;
16 V_inf = 2 ; // [m/ s ]
17 N = 5;
18 Tw = 70; // [ de gr e e C]
19 Tmi = 30; // [ de gr e e C]
20 L = 1; // [m]
21
22 rho = 1.165 ; // [ kg/m 3 ]
23 v = 16.00 *10^-6 ; // [m2/ s ]
24 Cp = 1.005*1000 ; // [ J/kg K]
25 k = 0.0267 ; // [W/m K]
26 Pr = 0.701;
27
28 // From eqn 5 . 1 0 . 2
30 Re = Vmax*D/v ;
31
32 // From f i g 5. 15
33 f = 0.37/4;
35 X = 1;
36 // From eqn 5 . 1 0 . 6
38
39 // At 70 de gr e e C
40 Pr1 = 0.694 ;
41 // From t a bl e 5 . 4 and 5 . 5
42
43 C1 = 0.27;
44 m = 0.63;
45 C2 = 0.93;
46
47 // Su b s t i t ut i ng i n Eqn 5 . 1 0 . 5
48 Nu = C1*C2*(Re^m)*(Pr ^0.36) *(Pr/Pr1)^(1/4);
49 h = Nu*k/D; // [W/m2 K]
74
50
51 // For 1 m l ong tube
52 m_dot = rho *(10*1.5*D*1)*2; // [ kg/ s ]
53
54 // Su b s t i t u t i ng m dot i n 5 . 3 . 4 and s o l vi ng , we ge t
55 function[f]=temp(Tmo)
56 f(1) = h*(%pi*D*L)*50*[(Tw -Tmi) -(Tw-Tmo(1))]/[
log((Tw -Tmi)/(Tw -Tmo(1)))]-m_dot*Cp*(Tmo(1)-
Tmi) ;
57 // h( %pi DL) N( (TwTmi ) (TwTmo) ) / l o g [ ( TwTmi
) /(TwTmo) ] m dot Cp(Tmo Tmi ) ;
58 funcprot (0);
59 endfunction
60
61 Tmo = 40; // I n i t i a l assumed val ue f o r f s o l v e
f unc t i o n
62 y = fsolve(Tmo ,temp);
63 printf(Tmo = %f \n,y);
64
65 printf( ( i i ) Exi t t e mpe r at ur e o f a i r = %f de gr e e C \
n,y);
Scilab code Exa 5.7.iii Heat transfer rate
1 clear;
2 clc;
3
5 // Chapter 5
7
8
9 // Example 5 . 7 ( i i i )
10 // Page 241
11 printf(Example 5 . 7 ( i i i ) , Page 241 \n \n);
75
12
13 D = 0.0125 ; // [m]
14 ST = 1.5*D ;
15 SL = 1.5*D ;
16 V_inf = 2 ; // [m/ s ]
17 N = 5;
18 Tw = 70; // [ de gr e e C]
19 Tmi = 30; // [ de gr e e C]
20 L = 1; // [m]
21
22 rho = 1.165 ; // [ kg/m 3 ]
23 v = 16.00 *10^-6 ; // [m2/ s ]
24 Cp = 1.005*1000 ; // [ J/kg K]
25 k = 0.0267 ; // [W/m K]
26 Pr = 0.701;
27
28 // From eqn 5 . 1 0 . 2
30 Re = Vmax*D/v ;
31
32 // From f i g 5. 15
33 f = 0.37/4;
35 X = 1;
36 // From eqn 5 . 1 0 . 6
38
39 // At 70 de gr e e C
40 Pr1 = 0.694 ;
41 // From t a bl e 5 . 4 and 5 . 5
42
43 C1 = 0.27;
44 m = 0.63;
45 C2 = 0.93;
46
47 // Su b s t i t ut i ng i n Eqn 5 . 1 0 . 5
48 Nu = C1*C2*(Re^m)*(Pr ^0.36) *(Pr/Pr1)^(1/4);
49 h = Nu*k/D; // [W/m2 K]
76
50
51 // For 1 m l ong tube
52 m_dot = rho *(10*1.5*D*1)*2; // [ kg/ s ]
53
54 // Su b s t i t u t i ng m dot i n 5 . 3 . 4 and s o l vi ng , we ge t
55 function[f]=temp(Tmo)
56 f(1) = h*(%pi*D*L)*50*[(Tw -Tmi) -(Tw-Tmo(1))]/[
log((Tw -Tmi)/(Tw -Tmo(1)))]-m_dot*Cp*(Tmo(1)-
Tmi) ;
57 // h( %pi DL) N( (TwTmi ) (TwTmo) ) / l o g [ ( TwTmi
) /(TwTmo) ] m dot Cp(Tmo Tmi ) ;
58 funcprot (0);
59 endfunction
60
61 Tmo = 40; // I n i t i a l assumed val ue f o r f s o l v e
f unc t i o n
62 y = fsolve(Tmo ,temp);
63
64 // Heat t r a n s f e r r a t e q
65 q = h*(%pi*D*L)*50*((Tw-Tmi)-(Tw-y))/(log((Tw-Tmi)/(
Tw -y)));
66
67 printf( ( i i i ) Heat t r a n s f e r r a t e per uni t l e ng t h t o
a i r = %f W,q);
77
Chapter 6
Heat Transfer by Natural
convection
Scilab code Exa 6.1 Average nusselt number
1 clear;
2 clc;
3
5 // Chapter 6
6 // Heat Tr ans f e r by Nat ur al Convect i on
7
8
9 // Example 6 . 1
10 // Page 258
11 printf(Example 6 . 1 , Page 258 \n \n);
12
13 H = 0.5 ; // [m]
14 T_h = 100; // [ de gr e e C]
15 T_l = 40; // [ de gr e e C]
16
17 v = 20.02*10^ -6 ; // [m/ s ]
18 Pr = 0.694;
19 k = 0.0297; // [W/m K]
78
20
21 T = (T_h+T_l)/2 + 273 ; // [ K]
22 printf(Mean f i l m t e mpe r at ur e = %f K \n,T);
23 B = 1/T;
24
25 Gr = 9.81*B*((T_h -T_l)*H^3)/(v^2);
26 Ra = Gr*Pr;
27
28 // ( a )
29 // Exact a n a l y s i s Equat i on 6 . 2 . 1 7
30 disp( ( a ) );
31 printf( Exact a n a l y s i s \n);
32 Nu_a = 0.64*( Gr ^(1/4))*(Pr^0.5) *((0.861+ Pr)^( -1/4));
33 printf(Nu L = %f \n,Nu_a);
34
35 // ( b)
36 // I n t e g r a l method Equat i on 6 . 2 . 2 9
37 disp( ( b) );
38 printf( I n t e g r a l method \n);
39 Nu_b = 0.68*( Gr ^(1/4))*(Pr^0.5) *((0.952+ Pr)^( -1/4));
40 printf(Nu L = %f \n,Nu_b);
41
42 // ( c )
43 // McAdams c o r r e l a t i o n Equat i on 6 . 2 . 3 0
44 disp( ( c ) );
45 printf(McAdams c o r r e l a t i o n \n);
46 Nu_c = 0.59*( Ra)^(1/4);
47 printf(Nu L = %f \n,Nu_c);
48
49 // ( d)
50 // Chur c hi l l and Chu c o r r e l a t i o n Equat i on 6 . 2 . 3 1
51 disp( ( d) )
52 printf( Chur c hi l l and Chu c o r r e l a t i o n \n);
53 Nu_d = 0.68 + 0.670*( Ra ^(1/4))/[1+(0.492/ Pr)^(9/16)
]^(4/9);
54 printf(Nu L = %f \n,Nu_d);
79
Scilab code Exa 6.2 Reduce the equation
1 clear;
2 clc;
3
5 // Chapter 6
7
8
9 // Example 6 . 2
10 // Page 259
12
13 Tm = 150 ; // [ de gr e e C]
15 v = 28.95*10^ -6 ; // [m2/ s ]
16 Pr = 0.683;
17 k = 0.0357 ; // [W/m K]
18
19 B = 1/(273+ Tm); // [ K1]
20
21 // f rom eqn 6 . 2 . 3 0
22 printf( Equat i on 6 . 2 . 3 0 \n h = k/L 0 . 5 9 [ 9 . 8 1 B(Tw
Ti nf ) ( L3) 0 . 6 8 3 / ( v 2) ] ( 1 / 4 ) \n)
23 // h = k/L 0 . 5 9 [ 9 . 8 1 B(TwTi nf ) ( L3) 0 . 6 8 3 / ( v 2)
] ( 1 / 4 ) ;
24 // s i mp l i f y i n g we ge t
25 // h = 1 . 3 8 [ ( TwTi nf ) /L] ( 1 / 4 )
26 printf( Reduces t o h = 1 . 3 8 [ ( TwTi nf ) /L] ( 1 / 4 ) \n)
27
28
29 // From eqn 6 . 2 . 3 3
30 // hL/k = 0 . 1 0 [ 9 . 8 1 B(TwTi nf ) ( L3) 0 . 6 8 3 / ( v 2)
80
] ( 1 / 3 ) ;
31 printf( Equat i on 6 . 2 . 3 3 \n hL/k = 0 . 1 0 [ 9 . 8 1 B(Tw
Ti nf ) ( L3) 0 . 6 8 3 / ( v 2) ] ( 1 / 3 ) \n);
32 // s i mp l i f y i n g
33 // h = 0 . 9 5 (TwTi nf ) 1/3
34 printf( Reduces t o h = 0 . 9 5 (TwTi nf ) 1/3 \n);
35
36 printf( where h i s e xpr e s s e d i n W/m2 K, (TwTi nf )
i n C and L i n met r es \n);
Scilab code Exa 6.3 Time for cooling of plate
1 clear;
2 clc;
3
5 // Chapter 6
7
8
9 // Example 6 . 3
10 // Page 260
12
13 s = 0.2 ; // [m]
14 d = 0.005 ; // [m]
15 rho = 7900 ; // [ kg/m 3 ]
16 Cp = 460 ; // [ J/kg K]
17
18 T_air = 20 ; // [ C]
19 // For 430 C t o 330 C
20 T_avg = 380 ; // [ C]
21 Tm = (T_avg + T_air)/2 ; // [ C]
22
23
81
24 v = 34.85*10^ -6 ; // [m2/ s ]
25 Pr = 0.680 ;
26 k = 0.0393 ; // [W/m K]
27
28 Re = 9.81*1/(273+ Tm)*(T_avg -T_air)*(s^3)/(v^2)*Pr;
29
30 // From eqn 6 . 2 . 3 1
31 Nu = 0.68 + 0.670*( Re ^(1/4))/[1+(0.492/ Pr)^(4/9)
]^(4/9);
32
33 h = Nu*k/s; // [W/m2 K]
34 t1 = rho*s*s*d*Cp/((s^2)*2*h)*log ((430 - T_air)/(330 -
T_air)); // [ s ]
35 printf(Time r e q ui r e d f o r t he pl a t e t o c o o l f rom 430
C t o 330 C i s %f s \n,t1);
36
37 // f o r 330 t o 230
38 h2 = 7.348 ; // [W/m2 K]
39 t2 = rho*s*s*d*Cp/((s^2)*2*h2)*log ((330 - T_air)/(230-
T_air)); // [ s ]
C t o 230 C i s %f s \n,t2);
41
42 // f o r 230 t o 130
43 h3 = 6.780; // [W/m2 K]
44 t3 = rho*s*s*d*Cp/((s^2)*2*h3)*log ((230 - T_air)/(130-
T_air)); // [ s ]
C t o 130 C i s %f s \n,t3);
46
47 // Tot al ti me
48
49 time = t1+t2+t3;
50 minute = time /60;
51 printf(Hence , ti me r e q ui r e d f o r t he pl a t e t o c o o l
f rom 430 C t o 130 C \n = %f s \n = %f min,time ,
minute);
82
Scilab code Exa 6.4 True air temperature
1 clear;
2 clc;
3
5 // Chapter 6
7
8
9 // Example 6 . 4
10 // Page 264
12
13 D = 0.006 ; // [m]
14 e = 0.1 ;
15 Ti = 800 ; // [ C]
16 Ta = 1000 ; // [ C]
17 // Rate at whi ch heat gai ne d = net r a di a nt heat ,
g i v e s h( Ta800) = 1306. 0 ; // [W/m 2 ]
18
19 // Usi ng t r i a l and e r r o r method
20 // Tr i a l 1
21 printf( Tr i a l 1 \n);
22 // Let Ta = 1000 de gr e e C
23 printf( Let Ta = 10000 C \n);
24
25 Tm = (Ta+Ti)/2;
27 v = 155.1*10^ -6 ; // [m2/ s ]
28 k = 0.0763 ; // [W/m K]
29 Pr = 0.717 ;
30
31 Gr = 9.81*1/1173*(200*D^3)/(v^2);
83
32 Ra = Gr*Pr ;
33
34 // From eqn 6 . 3 . 2
35 Nu = 0.36 + 0.518*( Ra ^(1/4))/[1+(0.559/ Pr)^(9/16)
]^(4/9);
36 h = Nu*k/D;
37 x = h*(Ta-Ti); // [W/m 2 ]
38 printf( Val ue o f h( Ta800) = %f W/m2 , whi ch i s much
l a r g e r than t he r e q ui r e d val ue o f 1306 W/m2 \n
,x);
39
40 // Tr i a l 2
41 printf(\ nTr i a l 2 \n);
42 // Let Ta = 900
44 Ra2 = 6.42 ;
45 Nu2 = 0.9841 ;
46 h2 = 12.15 ;
47 x2 = h2 *(900 -800);
48 printf( Val ue o f h( Ta800) = %f W/m2 , whi ch i s a
l i t t l e l e s s than t he r e q ui r e d val ue o f 1306 W/m2
\n,x2);
49
50 // Tr i a l 3
51 printf(\ nTr i a l 3 \n);
52 // Let Ta = 910
54 Ra3 = 6.93 ;
55 Nu3 = 0.9963 ;
56 h3 = 12.33 ;
57 x3 = h3 *(910 -800);
58 printf( Val ue o f h( Ta800) = %f W/m2 \ nThi s val ue
i s l i t t l e more than t he r e q ui r e d val ue o f 1306 W/
m2 \n,x3);
59 // I n t e r p o l a t i o n
60 T = 900 + (910 -900) *(1306 -x2)/(x3-x2);
61 printf(\nThe c o r r e c t val ue o f Ta obt ai ne d by
i n t e r p o l a t i o n i s %f C,T);
84
Scilab code Exa 6.5 Rate of heat ow by natural convection
1 clear;
2 clc;
3
5 // Chapter 6
7
8
9 // Example 6 . 5
10 // Page 269
12
13 T_p = 75 ; // Temperature o f abs or be r pl a t e , de gr e e
C
14 T_c = 55 ; // Temperature o f g l a s s c ove r , de gr e e C
15 L = 0.025 ; // [m]
16
17 H = 2 ; // [m]
18 Y = 70 ; // de gr e e
19
20 a = 19/180* %pi ; // [ Radi ans ]
21
22 r = H/L ;
23
24 T_avg = (T_p+T_c)/2+273 ; // [ K]
25 // Pr o p e r t i e s at 65 de gr e e C
26 k = 0.0294 ; // [W/m K]
27 v = 19.50*10^ -6 ; // [m2/ s ]
28 Pr = 0.695 ;
29
30 Ra = 9.81*(1/ T_avg)*(T_p -T_c)*(L^3)/(v^2)*Pr*cos(a);
31
85
32 // From eqn 6 . 4 . 3
33 Nu = 0.229*( Ra)^0.252;
34
35 h = Nu*k/L ; // [W/m2 K]
36
37 Rate = h*2*1*(T_p -T_c); // [W]
38
39 printf( Heat t r a n s f e r r a t e = %f W,Rate);
Scilab code Exa 6.6 Average Heat transfer coecient
1 clear;
2 clc;
3
5 // Chapter 6
7
8
9 // Example 6 . 6
10 // Page 270
12
13 T_air = 30 ; // [ C]
14 D = 0.04 ; // [m]
15 T_s = 70 ; // s u r f a c e t emper at ur e , [ C]
16 V = 0.3 ; // [m/ s ]
17
18 Tm = (T_air + T_s)/2 ; // [ C]
19 // Pr o p e r t i e s at Tm
20 v = 17.95*10^ -6 ; // [m2/ s ]
21 Pr = 0.698 ;
22 k = 0.0283 ; // [W/m K]
23
24 Gr = 9.81*1/323*( T_s -T_air)*(D^3)/v^2;
86
25 Re = V*D/v ;
26 X = Gr/Re^2 ;
27 printf( Si nc e Gr/Re2 = %f i s > 0 . 2 , we have a
combi ned c o nve c t i o n s i t u a t i o n . \n\n,X);
28
29 // From Eqn 5 . 9 . 8
30 Nu_forced = 0.3 + 0.62*( Re^0.5)*(Pr ^(1/3))/[[1+(0.4/
Pr)^(2/3) ]^(1/4) ]*[1+( Re /282000) ^(5/8) ]^(4/5);
31
32 // Su b s t i t u t i ng i n Eqn 6 . 5 . 1
33 Nu = Nu_forced *[1+6.275*(X)^(7/4) ]^(1/7);
34 h = Nu*(k/D);
35 printf(The Average heat t r a n s f e r c o e f f i c i e n t = %f W
/m2 K,h);
87
Chapter 7
Heat Exchangers
Scilab code Exa 7.1 Heat transfer coecient
1 clear;
2 clc;
3
5 // Chapter 7
6 // Heat Exchanger s
7
8
9 // Example 7 . 1
10 // Page 285
12
13 h = 2000 ; // [W/m2 K]
14 // From Tabl e 7 . 1
15 U_f = 0.0001 ; // f o u l i n g f a c t o r , m2K/W
16 h_f = 1/[1/h+U_f];
17 printf( Heat t r a n s f e r c o e f f i c i e n t i n c l u d i n g t he
e f f e c t o f f oul ung = %f W/m2 K \n,h_f);
18
19 p = (h-h_f)/h*100;
20 printf( Pe r c e nt age r e duc t i o n = %f \n,p);
88
Scilab code Exa 7.2 Area of heat exchanger
1 clear;
2 clc;
3
5 // Chapter 7
7
8
9 // Example 7 . 2
10 // Page 294
12
13 m = 1000 ; // [ kg/h ]
14 Thi = 50 ; // [ C]
15 The = 40 ; // [ C]
16 Tci = 35 ; // [ C]
17 Tce = 40 ; // [ C]
18 U = 1000 ; // OHTC, W/m2 K
19
20 // Usi ng Eqn 7 . 5 . 2 5
21 q = m/3600*4174*( Thi -The) ; // [W]
22
23 delta_T = ((Thi -Tce)-(The -Tci))/log((Thi -Tce)/(The -
Tci)); // [ C]
24 printf( de l t a T = %f \n\n,delta_T);
25
26 // T1 = Th and T2 = Tc
27 R = (Thi -The)/(Tce -Tci) ;
28 S = (Tce -Tci)/(Thi -Tci) ;
29 // From f i g 7 . 1 5 ,
30 F =0.91 ;
31
89
32 printf( Taki ng T1 = Th and T2 = Tc \n)
33 printf(R = %f , S = %f \n,R,S);
34 printf(Hence , F = %f \n \n,F);
35
36 // Al t e r na t i v e l y , t aki ng T1 = Tc and T2 = Th
37 R = (Tci -Tce)/(The -Thi);
38 S = (The -Thi)/(Tci -Thi);
39
40 // Agai n f rom f i g 7 . 1 5 ,
41 F =0.91 ;
42
43 printf( Taki ng T1 = Tc and T2 = Th \n)
44 printf(R = %f , S = %f \n,R,S);
45 printf(Hence , F = %f \n,F);
46
47 A = q/(U*F*delta_T);
48 printf(\nArea = %f m2 ,A);
Scilab code Exa 7.3 Mean temperature dierence
1 clear;
2 clc;
3
5 // Chapter 7
7
8
9 // Example 7 . 3
10 // Page 295
12
13 // Because o f change o f phase , Thi = The
14 Thi = 100 ; // [ C] , Sat ur at e d steam
15 The = 100 ; // [ C] , Condensed steam
90
16 Tci = 30 ; // [ C] , Cool i ng wat er i n l e t
17 Tce = 70 ; // [ C] , c o o l i n g wat er o u t l e t
18
19 R = (Thi -The)/(Tce -Tci) ;
20 S = (Tce -Tci)/(Thi -Tci) ;
21
22 // From f i g 7. 16
23 F = 1;
24
25 // For c ount e r f l ow ar r angement
26 Tm_counter = ((Thi -Tce)-(The -Tci))/log((Thi -Tce)/(
The -Tci)); // [ C]
28 Tm = F*Tm_counter ;
29 printf(Mean Temperaature Di f f e r e n c e = %f C,Tm)
Scilab code Exa 7.4.a Area of heat exchanger
1 clear;
2 clc;
3
5 // Chapter 7
7
8
9 // Example 7 . 4 ( a )
10 // Page 302
11 printf(Example 7 . 4 ( a ) , Page 302 \n \n);
12
13 // ( a )
14 printf( ( a ) \n);
15 // Usi ng Mean Temperature Di f f e r e n c e approach
16 m_hot = 10 ; // [ kg/min ]
17 m_cold = 25 ; // [ kg/min ]
91
18 hh = 1600 ; // [W/m2 K] , Heat t r a n s f e r c o e f f i c i e n t
on hot s i d e
19 hc = 1600 ; // [W/m2 K] , Heat t r a n s f e r c o e f f i c i e n t
on c ol d s i d e
20
21 Thi = 70 ; // [ C]
22 Tci = 25 ; // [ C]
23 The = 50 ; // [ C]
24
25 // Heat Tr ans f e r Rate , q
26 q = m_hot /60*4179*( Thi -The); // [W]
27
28 // Heat gai ne d by c ol d wat er = heat l o s t by t he hot
wat er
29 Tce = 25 + q*1/( m_cold /60*4174); // [ C]
30
31 // Usi ng e quat i on 7 . 5 . 1 3
32 Tm = ((Thi -Tci)-(The -Tce))/log((Thi -Tci)/(The -Tce));
// [ C]
33 printf(Mean Temperature Di f f e r e n c e = %f C \n,Tm);
34
35 U = 1/(1/hh + 1/hc); // [W/m2 K]
36 A = q/(U*Tm); // Area , [m 2 ]
37 printf( Area o f Heat Exchanger = %f m2 \n,A);
Scilab code Exa 7.4.b Exit temperature of hot and cold streams
1 clear;
2 clc;
3
5 // Chapter 7
7
8
92
9 // Example 7 . 4 ( b)
10 // Page 302
11 printf(Example 7 . 4 ( b) , Page 302 \n \n);
12
13 // Usi ng Mean Temperature Di f f e r e n c e approach
14 m_hot = 10 ; // [ kg/min ]
15 m_cold = 25 ; // [ kg/min ]
16 hh = 1600 ; // [W/m2 K] , Heat t r a n s f e r c o e f f i c i e n t
on hot s i d e
17 hc = 1600 ; // [W/m2 K] , Heat t r a n s f e r c o e f f i c i e n t
on c ol d s i d e
18
19 Thi = 70 ; // [ C]
20 Tci = 25 ; // [ C]
21 The = 50 ; // [ C]
22
23 // Heat Tr ans f e r Rate , q
24 q = m_hot /60*4179*( Thi -The); // [W]
25
26 // Heat gai ne d by c ol d wat er = heat l o s t by t he hot
wat er
27 Tce = 25 + q*1/( m_cold /60*4174); // [ C]
28
29 // Usi ng e quat i on 7 . 5 . 1 3
30 Tm = ((Thi -Tci)-(The -Tce))/log((Thi -Tci)/(The -Tce));
// [ C]
31 U = 1/(1/hh + 1/hc); // [W/m2 K]
32 A = q/(U*Tm); // Area , [m 2 ]
33
34 m_hot = 20 ; // [ kg/min ]
35 // Flow r a t e on hot s i d e i . e . hh i s doubl ed
36 hh = 1600*2^0.8 ; // [W/m2 K]
37 U = 1/(1/hh + 1/hc); // [W/m2 K]
38 m_hC_ph = m_hot /60*4179 ; // [W/K]
39 m_cC_pc = m_cold /60*4174 ; // [W/K]
41 C = m_hC_ph/m_cC_pc ;
42 NTU = U*A/m_hC_ph ;
93
43 printf(NTU = %f \n,NTU);
44
45 // From e quat i on 7 . 6 . 8
46 e = [1 - exp(-(1+C)*NTU)]/(1+C) ;
47
48 // The r e f or e ( Thi The ) /( Thi Tci ) = e , we ge t
49 The = Thi - e*(Thi - Tci); // [ C]
50
51 // Equat i ng t he heat l o s t by wat er t o heat gai ne d by
c ol d wat er , we ge t
52 Tce = Tci + [m_hC_ph *(Thi -The)]/ m_cC_pc;
53 printf( Exi t t e mpe r at ur e o f c ol d and hot s t r eam ar e
%f C and %f C r e s p e c t i v e l y . ,Tce ,The);
Scilab code Exa 7.5 Exit Temperature
1 clear;
2 clc;
3
5 // Chapter 7
7
8
9 // Example 7 . 5
10 // Page 304
12
13 mc = 2000 ; // [ kg/h ]
14 Tce = 40 ; // [ C]
15 Tci = 15 ; // [ C]
16 Thi = 80 ; // [ C]
17 U = 50 ; // OHTC, [W/m2 K]
18 A = 10 ; // Area , [m 2 ]
19
94
20 // Usi ng e f f e c t i v e NTU method
21 // Assumi ng m cC pc = (mC p ) s
22 NTU = U*A/(mc *1005/3600);
23 e = (Tce -Tci)/(Thi -Tci);
24 // From f i g 7 . 2 3 , no val ue o f C i s f ound
c o r r e s po ndi ng t o t he above val ue s , hence
as s umpt i on was wrong .
25 // So , m hC ph must be e qual t o (mC p ) s ,
pr oc e e di ng by t r a i l and e r r o r method
26
27
28 printf(m h( kg/h NTU C e
T he (C) T he (C) ( Heat Bal ance ) );
29
30 mh = rand (1:5);
31 NTU = rand (1:5);
32 The = rand (1:5);
33 The2 = rand (1:5);
34
35 mh(1) = 200
36 NTU (1) = U*A/(mh(1) *1.161);
37 // Cor r es pondi ng Val ues o f C and e f rom f i g 7. 23
38 C = .416;
39 e = .78;
40 //From Equat i on 7 . 6 . 2 Page 297
41 The (1) = Thi - e*(Thi -Tci)
42 //From Heat Bal ance
43 The2 (1) = Thi - mc *1005/3600*( Tce -Tci)/(mh(1)
*1.161);
44 printf(\n\n %i %. 3 f %. 3 f %. 3 f
%. 2 f %. 2 f ,mh(1),NTU (1),C,e,The (1),The2 (1))
;
45
46 mh(2) = 250
47 NTU (2) = U*A/(mh(2) *1.161);
49 C = .520;
50 e = .69;
95
*1.161);
55 printf(\n\n %i %. 3 f %. 3 f %. 3 f
;
56
57 mh(3) = 300
58 NTU (3) = U*A/(mh(3) *1.161);
60 C = .624;
61 e = .625;
*1.161);
66 printf(\n\n %i %. 3 f %. 3 f %. 3 f
;
67
68 mh(4) = 350
69 NTU (4) = U*A/(mh(4) *1.161);
71 C = .728;
72 e = .57;
*1.161);
77 printf(\n\n %i %. 3 f %. 3 f %. 3 f
;
78
79 mh(5) = 400
96
80 NTU (5) = U*A/(mh(5) *1.161);
82 C = .832;
83 e = .51;
*1.161);
88 printf(\n\n %i %. 3 f %. 3 f %. 3 f
;
89
90 clf();
91 plot(mh,The ,mh ,The2 ,[295 295 200] ,[0 39.2 39.2])
92 xtitle( The vs mh , mh ( kg/ hr ) , The (C) );
93 printf(\n\n From t he pl ot , val ue o f mh = 295 kg/ hr
and c o r r e s po ndi ng l y The = 39. 2 C)
97
Chapter 8
Condensation and boiling
Scilab code Exa 8.1 Average Heat Transfer Coecient
1 clear;
2 clc;
3
5 // Chapter 8
6 // Condens at i on and Bo i l i ng
7
8
9 // Example 8 . 1
10 // Page 318
12 Ts = 80 ; // [ C]
13 Tw = 70 ; // [ C]
14 L = 1 ; // [m]
15 g = 9.8 ; // [m/ s 2 ]
16
17 // Assumi ng c onde ns at e f i l m i s l ami nar and Re < 30
18 Tm = (Ts + Tw)/2 ;
20 rho = 978.8 ; // [ kg/m 3 ]
21 k = 0.672 ; // [W/m K]
98
22 u = 381 *10^ -6 ; // [ kg/m s ]
23 v = u/rho ;
24 // At 80 C,
25 lambda = 2309 ; // [ kJ/kg ]
26 // Su b s t i t u t i ng i n eqn 8 . 3 . 9 , we ge t
27 h = 0.943*[( lambda *1000*( rho^2)*g*(k^3))/((Ts-Tw)*u*
L)]^(1/4); // [W/m2 K]
28
29 rate = h*L*(Ts-Tw)/( lambda *1000); // [ kg/m s ]
30 Re = 4*rate/u;
31 printf( Assumi ng c onde ns at e f i l m i s l ami nar and Re <
30 \n);
32 printf(h = %f W/m2 K\n,h);
33 printf(Re L = %f \n,Re);
34 printf( I n i t i a l as s umpt i on was wrong , Now
c o n s i d e r i n g t he e f f e c t o f r i p p l e s , we ge t \n);
35
36 // Su b s t i t u t i ng h = Re ( lambda 1000) u/( 4L( TsTw) )
, i n eqn 8 . 3 . 1 2
37 Re = [[[4*L*(Ts -Tw)*k/( lambda *1000*u)*(g/(v^2))
^(1/3) ]+5.2]/1.08]^(1/1.22);
38 // From eqn 8 . 3 . 1 2
39 h = [Re /(1.08*( Re ^1.22) -5.2)]*k*((g/v^2) ^(1/3)); //
[W/m2 K]
40 m = h*L*10/( lambda *1000); // r a t e o f c onde ns at i on ,
[ kg/m s ]
41
42 printf(Re = %f \n,Re);
43 printf( Heat Tr ans f e r Co f f i c i e n t = %f W/m2 K \n,h)
;
44 printf( Rate o f c onde ns at i on = %f kg/m s ,m);
Scilab code Exa 8.2 Average heat transfer coecient and lm Reynolds
number
99
1 clear;
2 clc;
3
5 // Chapter 8
7
8
9 // Example 8 . 2
10 // Page 321
12
13 Ts = 262 ; // [ K]
14 D = 0.022 ; // [m]
15 Tw = 258 ; // [ K]
16
17 Tm = (Ts+Tw)/2;
19 rho = 1324 ; // [ kg/m 3 ]
20 k = 0.1008 ; // [W/m K]
21 v = 1.90*10^ -7 // [m2/ s ] ;
22 lambda = 215.1*10^3 ; // [ J/kg ]
23 g = 9.81 ; // [m/ s 2 ]
24 u = v*rho ; // Vi s c o s i t y
25
26 // From eqn 8 . 4 . 1
27 h = 0.725*[ lambda *(rho ^2)*g*(k^3) /((Ts -Tw)*u*D)
]^(1/4);
28
29 rate = h*%pi*D*(Ts -Tw) /lambda ; // [ kg/ s m]
30 Re = 4*rate/u ;
31
32 printf( Heat t r a n s f e r c o e f f i c i e n t = %f W/m2 K\n,h)
;
33 printf( Condens at i on r a t e per uni t l e ng t h = %f kg/ s
m \n,rate);
34 printf( Fi l m Reynol ds number = %f \n,Re);
100
Scilab code Exa 8.3 Length of the tube
1 clear;
2 clc;
3
4 // A TeTwtbook on HEAT TRANSFER by S P SUKHATME
5 // Chapter 8
7
8
9 // ETwample 8 . 3
10 // Page 322
12
13 m = 25/60 ; // [ kg/ s e c ]
14 ID = 0.025 ; // [m]
15 OD = 0.029 ; // [m]
16 Tci = 30 ; // [ C]
17 Tce = 70 ; // [ C]
18 g = 9.8 ; // [m/ s 2 ]
19
20 Ts = 100 ; // [ C]
21 // Assumi ng 5 . 3 . 2 i s val i d , p r o p e r t i e s at 50 C
23 rho = 988.1 ; // [ kg/m 3 ]
24 k = 0.648 ; // [W/m K]
25 v = 0.556*10^ -6 // [m2/ s ] ;
26 Pr = 3.54 ;
27 Re = 4*m/(%pi*ID*rho*v);
28 // From eqn 4 . 6 . 4 a
29 f = 0.005635;
30 // From eqn 5 . 3 . 2
31 Nu = 198.39 ;
32 h = Nu*k/ID ;
101
33
34 // Assumi ng ave r age wal l t e mpe r at ur e = 90 C
35 Tw = 90 ; // [ C]
36 Tm = (Tw+Ts)/2;
39 rho = 961.9 ; // [ kg/m 3 ]
40 k = 0.682 ; // [W/m K]
41 u = 298.6*10^ -6 ; // [ kg/m s ]
42 lambda = 2257*10^3 ; // [ J/kg ]
43
44 h = 0.725*[ lambda *(rho ^2)*g*(k^3) /((Ts -Tw)*u*OD)
]^(1/4);
45 // Equat i ng t he heat f l ow f rom t he c onde ns i ng steam
t o t he tube wal l , t o t he heat f l ow f rom t he tube
wal l t o t he f l o wi ng wat er .
46 // So l vi ng t he s i mp l i f i e d e quat i on
47 function[f] =temp(Tw)
48 f=(100 -Tw)^(3/4) -8.3096/[ log((Tw-Tci)/(Tw-Tce))
];
49 funcprot (0);
50 endfunction
51
52 T=fsolve(Tw,temp);
53 printf( Temperature obt ai ne d f rom t r i a l and e r r o r =
%f C \n,T);
54
56 hc = 21338.77/(100 -T)^(1/4); // [W/m2 K]
57 printf( h c = %f W/m2 K \n,hc);
58
59 // Now, e quat i ng t he heat f l o wi ng f rom t he
c onde ns i ng steam t o t he tube wal l t o t he heat
gai ne d by t he water , we have
60 function[g] =lngth(l)
61 g=hc*(%pi*OD*l)*(100-T)-m*4174*(Tce -Tci);
62 funcprot (0);
63 endfunction
102
64
65 l = 0; // ( i n i t i a l gues s , assumed val ue f o r f s o l v e
f unc t i o n )
66 L = fsolve(l,lngth);
67 printf(\ nLength o f t he tube = %f m \n,L);
Scilab code Exa 8.4 boiling regions
1 clear;
2 clc;
3
4 // Pr o p e r t i e s at (Tw+Ts ) /2 = 100. 5 de gr e e c e l s i u s
5 deltaT1 = 1; // i n de gr e e c e l s i u s
6 p1 = 7.55e-4; // [ K( 1) p1 i s c o e f f i c i e n t
o f c u b i c a l e xpans i on
7 v1 = 0.294e-6; // [m2/ s e c ] v i s c o s i t y
at 100. 5 de gr e e c e l s i u s
8 k1 = 0.683; // [W/mk ] t her mal
c o nduc t i v i t y
9 Pr1 = 1.74; // Pr andt l number
10 g = 9.81; // a c c e l e r a t i o n due t o
g r a v i t y
11 L = 0.14e-2; // di ame t e r i n met er s
12 // Pr o p e r t i e s at (Tw+Ts ) /2 =102. 5
15 v2 = 0.289e-6; // [m2/ s e c ] v i s c o s i t y at
102. 5 de gr e e c e l s i u s
16 k2 = 0.684; // [W/mk ] t her mal
18 // Pr o p e r t i e s at (Tw+Ts ) /2 =105
103
21 v3 = 0.284e-6; // [m2/ s e c ] v i s c o s i t y at
105 de gr e e c e l s i u s
22 k3 = 0.684; // [W/mk ] t her mal
24
25 function[Ra]= Rayleigh_no(p,deltaT ,v,Pr)
26 Ra = [(p*g*deltaT*L^3)/(v^2)]*Pr;
27 funcprot (0);
28 endfunction
29
30 function[q] = flux(k,deltaT ,Rai ,v)
31 q=(k/L)*( deltaT)*{0.36+(0.518* Rai ^(1/4))
/[1+(0.559/v)^(9/16) ]^(4/9) };
32 funcprot (0);
33 endfunction
34
35 Ra = Rayleigh_no(p1,deltaT1 ,v1,Pr1);
36 q1 = flux(k1,deltaT1 ,Ra ,Pr1);
37 printf(\n q/A = %. 1 f W/m2 at (TwTs )=1,q1);
40 printf(\n q/A = %. 1 f W/m2 at (TwTs )=5,q2);
43 printf(\n q/A = %. 1 f W/m2 at (TwTs ) =10,q3);
44
45 //At 100 de gr e e c e l s i u s
46 Cpl = 4.220; // [ kJ/kg ]
47 lamda = 2257; // [ kJ/kg ]
48 ul = 282.4e-6; // v i s c o s i t y i s i n kg/ms e c
49 sigma = 589e-4; // Sur f a c e t e ns i o n i s i n N/m
50 pl = 958.4; // de ns i t y i n kg/m3
51 pv =0.598; // de ns i t y o f vapour i n kg/m3
52 deltap = pl-pv;
53 Prl = 1.75; // Pr andt l no . o f l i q u i d
54 Ksf = 0.013;
104
55 function[q1]= heat_flux(deltaT)
56 q1 =141.32* deltaT ^3;
57 funcprot (0);
58 endfunction
59
60 printf(\n q/A at del t aT = 5 de gr e e c e l s i u s = %. 1 f W
/m2 ,heat_flux (5));
61 printf(\nq/A at del t aT = 10 de gr e e c e l s i u s = %. 1 f W
62 printf(\n q/A at del t aT =20 de gr e e c e l s i u s = %. 1 f W
63 // qi = [ h e a t f l u x ( 5 ) , h e a t f l u x ( 10) , h e a t f l u x ( 20) ] ;
64 q = [q1 q2 q3];
65 i=1;
66 while i<=10
67 T(i)=i;
68 ql(i) = heat_flux(i);
69 i=i+1;
70 end
71 plot2d ([1 5 10],q);
72 plot2d(T,ql);
73 xtitle( Bo i l i ng cur ve , (Tw Ts ) de gr e e c e l s i u s ,
Heat f l ux , ( q/A)W/m2 );
74 L1 = (L/2)*[g*(pl-pv)/sigma ]^(1/2);
75 printf(\n Peak heat f l u x L = %. 3 f ,L1);
76 f_L = 0.89+2.27* exp ( -3.44*L1 ^(0.5));
77 printf(\n f ( l ) = %. 4 f ,f_L);
78 q2 = f_L*{( %pi /24)*lamda *10^(3)*pv ^(0.5) *[sigma*g*(
pl -pv)]^(0.25) };
79 printf(\n q/A = %. 3 e W/m2 ,q2);
80
81 Tn = poly ([0], Tn );
82 Tn1 = roots (141.32* Tn^3 - q2);
83 printf(\n TwTs = %. 1 f de gr e e c e l s i u s ,Tn1(3));
84
85
86
87 printf(\n\n Minimum heat f l u x );
105
88 q3 = 0.09* lamda *10^3* pv*[sigma*g*(pl -pv)/(pl+pv)^(2)
]^(0.25);
89 printf(\n q/A = %d W/m2 ,q3);
90 printf(\n\n St a bl e f i l m b o i l i n g );
91 Ts1 = 140; // s u r f a c e t e mper at ur e i n de gr e e
c e l s i u s
c e l s i u s
c e l s i u s
94 Twm1 = (140+100) /2; //Mean f i l m t e mper at ur e
95 // p r o p e r t i e s o f steam at 120 de gr e e c e l s i u s and
1. 013 bar
96 kv = 0.02558; // t her mal c o nduc t i v i t y i n W/mK
97 pv1 = 0.5654; // vapor de ns i t y i n kg/m3
98 uv =13.185*10^( -6); // v i s c o s i t y o f vapour i n kg/m
s e c
99 lamda1 = (2716.1 -419.1) *10^(3);// Lat ent heat o f
f u s i o n i n J/kg
100 hc = 0.62*[( kv^3)*pv*(pl -pv)*g*lamda1 /(L*uv
*(140 -100))]^(0.25);
101 printf(\n hc = %. 2 f W/m2 ,hc);
102 qrad = 5.67*10^( -8) *(413^4 - 373^4) /[(1/0.9) +1 -1];
103 printf(\n q/A due t o r a d i a t i o n = %. 2 f W/m2 ,qrad);
104 hr = qrad /(413 -373);
105 printf(\n hr = %. 2 f W/m2 K ,hr);
106
107 printf(\n Si nc e hr<hc );
108 printf(\n The t o t a l heat t r a n s f e r c o e f f i c i e n t );
109 h = hc + 0.75*hr;
110 printf( h = %. 2 f W/m2 K,h);
111 printf(\n Tot al heat f l u x = %. 3 f W/m2 K,h
*(140 -100));
112
113 hc_200 = 0.62*[( kv^3)*pv*(pl -pv)*g*lamda1 /(L*uv
*(200 -100))]^(0.25);
114 qrad1 = 5.67*10^( -8) *(473^4 - 373^4) /[(1/0.9) +1-1];
115 hr_200 = qrad1 /(200 -100);
106
116 printf(\n\n hc = %. 2 f W/m2 ,hc_200);
117 printf(\n hr = %. 2 f W/m2 K,hr_200);
118 printf(\n q/A due t o r a d i a t i o n = %. 2 f W/m2 ,qrad1)
;
119 h_200 = hc_200 +0.75* hr_200;
120 printf(\n Tot al heat f l u x = %d W/m2 ,h_200 *100);
121 hc_600 = 0.62*[( kv^3)*pv*(pl -pv)*g*lamda1 /(L*uv
*(600 -100))]^(0.25);
122 qrad2 = 5.67*10^( -8) *(873^4 - 373^4) /[(1/0.9) +1-1];
123 hr_600 = qrad1 /(600 -100)
124 printf(\n\n hc = %. 2 f W/m2 ,hc_600);
125 printf(\n hr = %. 2 f W/m2 K,hr_600);
126 printf(\n q/A due t o r a d i a t i o n = %. 2 f W/m2 ,qrad2)
;
Scilab code Exa 8.5 Initial heat transfer rate
1 clear;
2 clc;
3
5 // Chapter 8
7
8
9 // Example 8 . 5
10 // Page 337
12
13 D = 0.02 ; // [m]
14 l = 0.15 ; // [m]
15 T = 500+273 ; // [ K]
16 Tc = -196+273 ; // [ K]
17 e = 0.4;
18 s = 5.670*10^ -8;
107
19 // Fi l m b o i l i n g wi l l occur , hence eqn 8 . 7 . 9 i s
a p p l i c a b l e
20 Tm = (T+Tc)/2;
21
22 // Pr o p e r t i e s
23 k = 0.0349 ; // [W/m K]
24 rho = 0.80 ; // [ kg/m 3 ]
25 u = 23*10^ -6 ; // [ kg/m s ]
26
27 Cp_avg = 1.048 ; // [ kJ/kg J ]
28 rho_liq = 800 ; // [ kg/m 3 ]
29 latent = 201*10^3 ; // [ J/kg ]
30
31 lambda = [latent + Cp_avg *(Tm-Tc)*1000]; // [ J/kg ]
32 h_c = 0.62*[((k^3)*rho *799.2*9.81* lambda)/(D*u*(T-Tc
))]^(1/4); // [W/m2 K]
33
34 // Taki ng t he e mi s s i v i t y o f l i q u i d s u r f a c e t o be
uni t y and us i ng e quat i on 3 . 9 . 1 , t he exchange o f
r a di a nt heat f l u x
35 flux = s*(T^4-Tc^4) /(1/e+1/1 -1); // [W/m 2 ]
36 h_r = flux/(T-Tc);
37
38 // Si nc e h r < h c , t o t a l heat t r a n s f e r c o e f f i c i e n t
i s det er mi ned f rom eqn 8 . 7 . 1 1
39 h = h_c +3/4* h_r; // [W/m2 K]
40
41 flux_i = h*(T-Tc);
42 Rate = flux_i*%pi*D*l; // [W]
43
44 printf( I n i t i a l heat f l u x = %f W/m2 \n,flux_i);
45 printf( I n i t i a l heat t r a n s f e r r a t e = %f W,Rate);
108
Chapter 9
Mass Transfer
Scilab code Exa 9.1 Composition on molar basis
1 clear;
2 clc;
3
5 // Chapter 9
6 // Mass Tr ans f e r
7
8
9 // Example 9 . 1
10 // Page 349
12
13 w_a = 0.76 ;
14 w_b = 0.24 ;
15 m_a = 28 ; // [ kg/kg mol e ]
16 m_b = 32 ; // [ kg/kg mol e ]
17
18 x_a = (w_a/m_a)/(w_a/m_a+w_b/m_b);
19 x_b = (w_b/m_b)/(w_a/m_a+w_b/m_b);
20 printf(The mol ar f r a c t i o n s ar e gi ve n by \n);
21 printf( x a = %f \n,x_a);
109
22 printf( x b = %f ,x_b);
Scilab code Exa 9.2 Diusion coecient of napthalene
1 clear;
2 clc;
3
5 // Chapter 9
7
8
9 // Example 9 . 2
10 // Page 350
12
13 // From Tabl e 9 . 1 at 1 atm and 25 C
14 Dab = 0.62*10^ -5 ; // [m2/ s ]
15 // The r e f or e at 2 atm and 50 C
16 Dab2 = Dab *(1/2) *(323/298) ^1.5 ;
17 printf(Dab at 2 atm & 50 C = %e m2/ s ,Dab2);
Scilab code Exa 9.3.a Rate of hydrogen diusion
1 clear;
2 clc;
3
5 // Chapter 9
7
8
9 // Example 9 . 3 ( a )
110
10 // Page 352
12
13 t = 0.04 ; // [m]
14 A = 2 ; // [m 2 ]
15 rho1 = 0.10 ;
16 rho2 = 0.01 ;
17 D_400 = 1.6*10^ -11 ; // at 400K [m2/ s ]
18
19 // Mass Di f f u s i o n i n s o l i d s o l ut i o n , assumi ng Fi c ks
l aw i s v a l i d & s t e ady s t a t e and one di me ns i onal
d i f f u s i o n
20
21 // Subt i t ut i ng t he va l ue s i n eqn 9 . 3 . 3 , At 400 K
22
23 m_400 = A*D_400 *(rho1 -rho2)/t; // [ kg/ s ]
24 printf( Rate o f d i f f u s i o n o f Hydrogen at 400 K = %e
kg/ s \n,m_400);
Scilab code Exa 9.3.b Rate of hydrogen diusion
1 clear;
2 clc;
3
5 // Chapter 9
7
8
9 // Example 9 . 3 ( b)
10 // Page 352
12
13 t = 0.04 ; // [m]
14 A = 2 ; // [m 2 ]
111
15 rho1 = 0.10 ;
16 rho2 = 0.01 ;
17 D_1200 = 3.5*10^ -8 ; // at 1200 k [m2/ s ]
18
19 // Mass Di f f u s i o n i n s o l i d s o l ut i o n , assumi ng Fi c ks
l aw i s v a l i d & s t e ady s t a t e and one di me ns i onal
d i f f u s i o n
20
21 // At 1200 K
22 // From eqn 9 . 3 . 3
23
24 m_1200 = A*D_1200 *(rho1 -rho2)/t ;
25 printf( ( b) Rate o f d i f f u s i o n o f Hydrogen at 1200 K
= %e kg/ s \n,m_1200);
Scilab code Exa 9.4.a Rate of loss of ammonia
1 clear;
2 clc;
3
5 // Chapter 9
7
8
9 // Example 9 . 4 ( a )
10 // Page 356
12
13 L = 1 ; // [m]
14 D = 0.005 ; // [m]
15 Pa1 = 1 ; // [ atm ]
16 Pa2 = 0 ;
17 R = 8314 ;
18 T = 298 ; // [ K]
112
19
20 // Assumi ng Equi mol al c ount e r d i f f u s i o n
21 // From Tabl e 9 . 1
22 Dab = 2.80*10^ -5 ; // [m2/ s ]
23 // Subs t i t ui ng i n eqn 9 . 4 . 1 2
24 Na = -[Dab/(R*T)*(Pa2 -Pa1)*(1.014*10^5)/L]*( %pi*(D
/2) ^2);
25 R_NH3 = Na*17 ; // [ kg/ s ]
26
27 printf(Na = Nb = %e ( kg mol e ) /m2 s \n,Na);
28 printf( Rate at whi ch ammonia i s l o s t t hr ough t he
tube = %e kg/ s \n,R_NH3);
Scilab code Exa 9.4.b Rate at which air enters the tank
1 clear;
2 clc;
3
5 // Chapter 9
7
8
9 // Example 9 . 4 ( b)
10 // Page 356
12
13 L = 1 ; // [m]
14 D = 0.005 ; // [m]
15 Pa1 = 1 ; // [ atm ]
16 Pa2 = 0 ;
17 R = 8314 ;
18 T = 298 ; // [ K]
19
20 // Si nc e t he tank i s l a r g e and t he pr e s s ur e and
113
t e mpe r at ur e at t he two ends o f t he same tube ar e
same , we ar e assumi ng Equi mol al c ount e r d i f f u s i o n
21 // From Tabl e 9 . 1
22 Dab = 2.80*10^ -5 ; // [m2/ s ]
23 // Subs t i t ui ng i n eqn 9 . 4 . 1 2
24 Na = -[Dab/(R*T)*(Pa2 -Pa1)*(1.014*10^5)/L]*( %pi*(D
/2) ^2);
25
26 // Si nc e e qui mol al c ount e r d i f f u s i o n i s t aki ng pl a c e
27 Nb = - Na ;
28 // t h e r e f o r e r a t e at whi ch a i r e nt e r s t he tank
29 R_air = abs(Nb)*29 ; // [ kg/ s ]
30
31 printf( Rate at whi ch a i r e nt e r s t he tank = %e kg/ s
,R_air);
Scilab code Exa 9.5 Rate of evaporation
1 clear;
2 clc;
3
5 // Chapter 9
7
8
9 // Example 9 . 5
10 // Page 359
12
13 // Evapor at i on o f water , one di me ns i onal
14 T_w = 20+273 ; // [ K]
15 D = 0.04 ; // [m]
16 h = 0.20 ; // [m]
17 h_w = 0.03 ; // [m]
114
18
19 P = 1.014*10^5; // [ Pa ]
20 R = 8314 ; // [ J/kg mol e K]
21 P_sat = 0.02339 ; // [ bar ]
22 x_a1 = P_sat /1.014 ; // mol e f r a c t i o n at l i q vap
i n t e r f a c e
23 x_a2 = 0 ; // mol e f r a c t i o n at open top
24 c = P/(R*T_w);
25 // From Tabl e 9 . 2
26 Dab = 2.422*10^ -5 ; // [m2/ s ]
27
28 // Su b s t i t u t i ng above va l ue s i n eqn 9 . 4 . 1 8
29 flux = 0.041626* Dab /0.17* log((1-0)/(1-x_a1)); // [ kg
mol e /m2 s ]
30 rate = flux *18*( %pi/4)*(D^2);
31
32 printf( Rate o f e vapor at i on o f wat er = %e kg/ s ,rate
);
Scilab code Exa 9.6 Rate of evaporation
1 clear;
2 clc;
3
5 // Chapter 9
7
8
9 // Example 9 . 6
10 // Page 364
12
13 l = 1; // l engt h , [m]
14 w = 0.25; // wi dth , [m]
115
15 T = 293 ; // Temperature , [ K]
16 rho_infinity = 0; // [ kg/m 3 ]
17 R = 8314; // [ J/ kg K]
18
20 v = 15.06*10^ -6; // [m2/ s ]
21 // From Tabl e 9 . 2
22 Dab = 2.4224*10^ -5; // [m2/ s ]
23 Re = 2.5/v;
24 Sc = v/Dab;
25 // Si nc e Re > 3105 , we may assume l ami nar boundary
l a y e r
26 Sh = 0.664* Sc ^(1/3)*Re ^(1/2); // Sherwood number
27 h = Sh*Dab;
28
29 p_aw = 2339; // Sa t ur a t i o n pr e s s ur e o f wat er at 20
de gr e e C. [ N/m 2 ]
30 rho_aw = p_aw/(R/18*T); // [ kg/m 3 ]
31 rho_a_inf = 0 ; // s i n c e a i r i n t he f r e e s t r eam i s
dry
32 m_h = h*(2*l*w)*(rho_aw -rho_infinity);
33 printf( Rate o f e vapor at i on f rom pl a t e = %e kg/ s ,
m_h);
Scilab code Exa 9.7.a Mass transfer coecient Colburn anology
1 clear;
2 clc;
3
5 // Chapter 9
7
8
9 // Example 9 . 7 ( a )
116
10 // Page 366
12
13 D = 0.04 ; // [m]
14 V = 1.9 ; // [m/ s ]
15
16 // ( a ) Col burn anol ogy and Gn i e l i n s k i e quat i on
18 v = 15.718*10^ -6 ; // [m2/ s ]
19 rho = 1.177 ; // [ kg/m 3 ]
20 Pr = 0.7015 ;
21 Cp = 1005 ; // [ J/kg K]
22 k = 0.02646 ; // [W/m K]
23 // From Tabl e 9 . 2
24 Dab = 2.54 * 10^-5 ; // [m2/ s ]
25 Sc = v/Dab ;
26 Re = V*D/v;
27 // The f l ow i s t ur bul e nt and eqn 9 . 6 . 5 may be
a ppl i e d
28 // l e t r = h/h m
29 r = rho*Cp*((Sc/Pr)^(2/3));
30 // From Bl a s i us e quat i on 4 . 6 . 4 a
31 f = 0.079* Re^( -0.25);
32 // Su b s t i t u t i ng t h i s val ue i nt o Gn i e l i n s k i e quat i on
5 . 3 . 2
33 Nu = [(f/2)*(Re -1000)*Pr ]/[1+12.7*((f/2) ^(1/2))*((Pr
^(2/3)) -1)];
34 h = Nu*k/D;
35 h_m = h/r; // [m/ s ]
36
37 printf(h m us i ng Col burn anol ogy and Gn i e l i n s k i
e quat i on = %f \n,h_m);
Scilab code Exa 9.7.b Mass transfer coecient Gnielinski equation
117
1 clear;
2 clc;
3
5 // Chapter 9
7
8
9 // Example 9 . 7 ( b)
10 // Page 366
12
13 D = 0.04 ; // [m]
14 V = 1.9 ; // [m/ s ]
15
16 // ( b) mess t r a n s f e r c o r r e l a t i o n e q ui v a l e nt t o t he
Gl e i l i n s k i e quat i on
17
19 v = 15.718*10^ -6 ; // [m2/ s ]
20 rho = 1.177 ; // [ kg/m 3 ]
21 Pr = 0.7015 ;
22 Cp = 1005 ; // [ J/kg K]
23 k = 0.02646 ; // [W/m K]
24 // From Tabl e 9 . 2
25 Dab = 2.54 * 10^-5 ; // [m2/ s ]
26 Sc = v/Dab ;
27 Re = V*D/v;
28
30 f = 0.079* Re^( -0.25);
31
32 // Su b s t i t u t i ng i n eqn 9 . 6 . 7
33 Sh_D = [(f/2)*(Re -1000)*Sc ]/[1+12.7*((f/2))*((Sc
^(2/3)) -1)];
34 h_m1 = Sh_D*Dab/D;
35
36 printf( ( b) h m = %f \n,h_m1);
118
Scilab code Exa 9.7.c To show mass ux of water vapour is small
1 clear;
2 clc;
3
5 // Chapter 9
7
8
9 // Example 9 . 7 ( c )
10 // Page 366
11 printf(Example 9 . 7 ( c ) , Page 366 \n \n);
12
13 D = 0.04 ; // [m]
14 V = 1.9 ; // [m/ s ]
15
16 // ( c ) To show t hat mass f l u x o f wat er i s ver y s mal l
compared t o t he mass f l u x o f a i r f l o wi ng i n t he
pi pe
18 v = 15.718*10^ -6 ; // [m2/ s ]
19 rho = 1.177 ; // [ kg/m 3 ]
20 Pr = 0.7015 ;
21 Cp = 1005 ; // [ J/kg K]
22 k = 0.02646 ; // [W/m K]
23 // From Tabl e 9 . 2
24 Dab = 2.54 * 10^-5 ; // [m2/ s ]
25 Sc = v/Dab ;
26 Re = V*D/v;
27 // The f l ow i s t ur bul e nt and eqn 9 . 6 . 5 may be
a ppl i e d
28 // l e t r = h/h m
29 r = rho*Cp*((Sc/Pr)^(2/3));
119
31 f = 0.079* Re^( -0.25);
32
33 // From steam t a bl e
34 rho_aw = 1/38.77 ; // [ kg/m 3 ]
35 // l e t X = ( m a/A) max
36 X = f*rho_aw; // [ kg/m2 s ]
37
38 // l e t Y = mass f l u x o f a i r i n pi pe = (m/A)
39 Y = rho*V ; // [ kg/m2 s ]
40 ratio = X/Y ;
41 percent = ratio *100;
42
43 printf( ( c ) ( m a/A) max /( m a/A) = %f pe r c e nt Thus ,
mass f l u x o f wat er i s ver y s mal l compared t o t he
mass f l u x o f a i r f l o wi ng i n t he pi pe . ,percent )
;
Scilab code Exa 9.8 Mass fraction
1 clear;
2 clc;
3
5 // Chapter 9
7
8
9 // Example 9 . 8
10 // Page 369
12
13 V = 0.5 ; // [m/ s ]
14 T_h = 30 ; // [ C]
15 T_c = 26 ; // [ C]
120
16 Tm = (T_h+T_c)/2;
18 rho = 1.173 ; // [ kg/m 3 ]
19 Cp = 1005 ; // [ J/kg K]
20 k = 0.02654 ; // [W/m K]
21
22 alpha = k/(rho*Cp); // [m2/ s ]
23
24 // From Tabl e 9 . 2 at 301 K
25 Dab = 2.5584*10^ -5 ; // [m2/ s ]
26 lambda = 2439.2*10^3 ; // [ J/kg ]
27
28 // Su b s t i t u t i ng i n e quat i on 9 . 7 . 5
29 // l e t d i f f e r e n c e = rho awr ho a i n f i n i t y
30 difference = rho*Cp*(( alpha/Dab)^(2/3))*(T_h -T_c)/
lambda;
31
32 // From steam t a bl e
33 Psat = 3363;
34 rho_aw = Psat /(8314/18*299);
35 rho_inf = rho_aw - difference;
36 x = rho_inf/rho; // mol e f r a c t i o n o f wat er vapour i n
a i r s t r eam
37
38 PP = rho_inf *8314/18*303; // Pa r t i a l pr e s s ur e o f
wat er vapour i n a i r s t r eam
39 // From steam t a bl e p a r t i a l pr e s s ur e o f wat er vapour
at 30 C
40 PP_30 = 4246 ; // [ N/m 2 ]
41
42 rel_H = PP/PP_30;
43 percent = rel_H *100;
44
45 printf( Re l a t i v e humi di t y = %f i . e . %f pe r c e nt ,
rel_H ,percent);
121

Textbook of Heat Transfer - Sukhatme, S. P.

Uploaded by

Copyright:

Available Formats

Textbook of Heat Transfer - Sukhatme, S. P.

Uploaded by

Document Information

Original Description:

Original Title

Copyright

Available Formats

Share this document

Share or Embed Document

Sharing Options

Did you find this document useful?

Is this content inappropriate?

Copyright:

Available Formats

Textbook of Heat Transfer - Sukhatme, S. P.

Uploaded by

Copyright:

Available Formats

Scilab Textbook Companion for

You might also like