From b1f5c3f8d6671b4331cef1dcebdf63b7a43a3a2b Mon Sep 17 00:00:00 2001 From: priyanka Date: Wed, 24 Jun 2015 15:03:17 +0530 Subject: initial commit / add all books --- 617/CH9/EX9.1/Example9_1.sci | 34 +++++++++++++++++++++++++++ 617/CH9/EX9.10/Example9_10.sci | 27 +++++++++++++++++++++ 617/CH9/EX9.11/Example9_11.sci | 22 ++++++++++++++++++ 617/CH9/EX9.2/Example9_2.sci | 21 +++++++++++++++++ 617/CH9/EX9.3/Example9_3.sci | 11 +++++++++ 617/CH9/EX9.4/Example9_4.sci | 13 +++++++++++ 617/CH9/EX9.5/Example9_5.sci | 13 +++++++++++ 617/CH9/EX9.6/Example9_6.sci | 17 ++++++++++++++ 617/CH9/EX9.7/Example9_7.sci | 22 ++++++++++++++++++ 617/CH9/EX9.8/Example9_8.sci | 53 ++++++++++++++++++++++++++++++++++++++++++ 617/CH9/EX9.9/Example9_9.sci | 19 +++++++++++++++ 617/CH9/EX9.9b/Example9_9b.sci | 26 +++++++++++++++++++++ 12 files changed, 278 insertions(+) create mode 100755 617/CH9/EX9.1/Example9_1.sci create mode 100755 617/CH9/EX9.10/Example9_10.sci create mode 100755 617/CH9/EX9.11/Example9_11.sci create mode 100755 617/CH9/EX9.2/Example9_2.sci create mode 100755 617/CH9/EX9.3/Example9_3.sci create mode 100755 617/CH9/EX9.4/Example9_4.sci create mode 100755 617/CH9/EX9.5/Example9_5.sci create mode 100755 617/CH9/EX9.6/Example9_6.sci create mode 100755 617/CH9/EX9.7/Example9_7.sci create mode 100755 617/CH9/EX9.8/Example9_8.sci create mode 100755 617/CH9/EX9.9/Example9_9.sci create mode 100755 617/CH9/EX9.9b/Example9_9b.sci (limited to '617/CH9') diff --git a/617/CH9/EX9.1/Example9_1.sci b/617/CH9/EX9.1/Example9_1.sci new file mode 100755 index 000000000..ab297d0df --- /dev/null +++ b/617/CH9/EX9.1/Example9_1.sci @@ -0,0 +1,34 @@ +clc(); +clear; + +// To find the temperature at the free end is made of copper iron and glass + +D = 3/48; // diameter in ft +L = 9/12; // Length of steam vessel in ft +T1 = 210; // Vessel temperature in degF +T2 = 80; // Air temperature in degF +th0 = T1-T2; // Temperature difference in degF +h = 1.44; // Assumed heat coefficient in Btu/hr-ft^2-degF +C = %pi*D; // Circumference of vessel in ft +A = %pi*D*D/4; // Area of vessel in ft^2 + +// For copper +k1 = 219; // Heat conductivity of copper in Btu/hr-ft-degF +m1 = sqrt(h*C/(k1*A)); // in /ft +th1 = th0*2/(exp(m1*L)+exp(-m1*L)); +Tl1 = round(th1+T2); // The temperaure at the free end in degF +printf("Temperature at free end of the copper rod is %d degF \n",Tl1); + +// For iron +k2 = 36; // heat conductivity of copper in Btu/hr-ft-degF +m2 = sqrt(h*C/(k2*A)); // in /ft +th2 = th0*2/(exp(m2*L)+exp(-m2*L)); +Tl2 = th2+T2; // The temperaure at the free end in degF +printf(" Temperature at free end of the iron rod is %.2f degF \n",Tl2); + +// For glass +k3 = 0.64; // Heat conductivity of copper in Btu/hr-ft-degF +m3 = sqrt(h*C/(k3*A)); // in /ft +th3 = th0*2/(exp(m3*L)+exp(-m3*L)); +Tl3 = th3+T2; // The temperaure at the free end in degF +printf(" Temperature at free end of the glass rod is %.2f degF \n",Tl3); diff --git a/617/CH9/EX9.10/Example9_10.sci b/617/CH9/EX9.10/Example9_10.sci new file mode 100755 index 000000000..7713dbf7b --- /dev/null +++ b/617/CH9/EX9.10/Example9_10.sci @@ -0,0 +1,27 @@ +clc(); +clear; + +// To calculate the terminal temperature of oil and water + +To1=160; // inlet temperature of oil in degF +Cpo=0.5; // Specific heat capacity in Btu/lb-degf +Tw1=60; // Inlet temperature of water in degF +mo=1000; // Mass flow rate of oil in lb/hr +mw=2500; // Mass flow rate of water in lb/hr +Cpw=1; // Heat capacity of water in Btu/hr +X=mo*Cpo/(mw*Cpw); // Ratio of flow rates +UA=1.15*mo*Cpo; +B=UA/mo*Cpo; + +// from the graph, we can locate the point of A and B And corresponding effectiveness ratio +E=0.86; // Effectiveness ratio +To2=To1-E*(To1-Tw1); // Outlet temperature of oil in degF +printf("The outlet temperature of oil is %d degF \n",To2); + +q=mo*Cpo*(To1-To2); // Heat transferred in Btu/hr +Tw2=Tw1+(q/(mw*Cpw)); // Outlet temperature of oil in degF +printf(" The outlet tempearture of water is %.1f degF",Tw2); + + + + \ No newline at end of file diff --git a/617/CH9/EX9.11/Example9_11.sci b/617/CH9/EX9.11/Example9_11.sci new file mode 100755 index 000000000..c249ecca3 --- /dev/null +++ b/617/CH9/EX9.11/Example9_11.sci @@ -0,0 +1,22 @@ +clc(); +clear; + +// To compute the temprature distribution +h=1; // Heat transfer coefficient in Btu/hr-ft^2-degF +x=1; // Assumed thickness in ft +k=1; // Thermal conductivity in Btu/hr-ft-degF +N=h*x/k; +t0=600; +t4=200; +t1=[500 550 550 525 525 512.5 512.5 512.5 506.2 506.2 506.2 506.2 503.1 503.1]; +t2=[450 450 450 450 425 425 425 412.5 412.5 412.5 406.3 406.3 406.3 403.1]; +t3=[350 350 325 325 325 325 312.5 312.5 312.5 306.3 306.3 303.1 303.1 303.1]; + +// Assumed temperatures in degF for points 1 2 & 3 respectively +for i=1:14 +th1(i)=t0+t2(i)-2*t1(i); +th2(i)=t1(i)+t3(i)-2*t2(i); +th3(i)=t2(i)+t4-2*t3(i); +printf("Assuming t1=%.1f degF t2=%.1fdegF t3=%.1fdegF \n th1=%.1fdegF th2=%.1fdegF th3=%.1fdegF \n \n",t1(i),t2(i),t3(i),th1(i),th2(i),th3(i)); +end +printf("This way assumption must be continued till all sink strengths are zero"); \ No newline at end of file diff --git a/617/CH9/EX9.2/Example9_2.sci b/617/CH9/EX9.2/Example9_2.sci new file mode 100755 index 000000000..3cb9cf925 --- /dev/null +++ b/617/CH9/EX9.2/Example9_2.sci @@ -0,0 +1,21 @@ +clc(); +clear; + +// To find the temperature at the free end is made of copper iron and glass + +D = 3/48; // diameter in ft +L = 9/12; // Length of steam vessel in ft +T1 = 210; // Vessel temperature in degF +T2 = 80; // Air temperature in degF +th0 = T1-T2; // Temperature difference in degF +h = 1.44; // Assumed heat coefficient in Btu/hr-ft^2-degF +C = %pi*D; // Circumference of vessel in ft +A = %pi*D*D/4; // Area of vessel in ft^2 + +k = 36; // heat conductivity of copper in Btu/hr-ft-degF +m = sqrt(h*C/(k*A)); // in /ft +q=k*A*m*th0*(exp(m*L)-exp(-m*L))/(exp(m*L)+exp(-m*L)); +// Heat loss by iron rod in Btu/hr +printf("The rate of heat loss by iron rod is %.d Btu/hr",q); + + diff --git a/617/CH9/EX9.3/Example9_3.sci b/617/CH9/EX9.3/Example9_3.sci new file mode 100755 index 000000000..0d5bf2c16 --- /dev/null +++ b/617/CH9/EX9.3/Example9_3.sci @@ -0,0 +1,11 @@ +clc(); +clear; + +// To calculate the heat transfer coefficient + +x = 3/96; // Thickness of plate in ft +k = 220; // thermal conductivity in Btu/hr-ft-degF +h1 = 480; // Inner film coefficient in Btu/hr-ft^2-degF +h2 = 1250; // Outer film coefficient in Btu/hr-ft^2-degF +U = 1/((1/h1)+(x/k)+(1/h2)); // Overall heat transer coeeficient in Btu-hr-ft^2-degF +printf("Overall heat transfer coefficient is %d Btu/hr-ft^2-degF",U); \ No newline at end of file diff --git a/617/CH9/EX9.4/Example9_4.sci b/617/CH9/EX9.4/Example9_4.sci new file mode 100755 index 000000000..5d93ead76 --- /dev/null +++ b/617/CH9/EX9.4/Example9_4.sci @@ -0,0 +1,13 @@ +clc(); +clear; + +// To calculate the overall heat transfer coefficient + +r2 = 3/96; // Outer radius in ft +x = 0.1/12; // Thickness of plate in ft +r1 = r2-x; // Outer radius in ft +k = 200; // thermal conductivity in Btu/hr-ft-degF +h1 = 280; // Inner film coefficient in Btu/hr-ft^2-degF +h2 = 2000; // Outer film coefficient in Btu/hr-ft^2-degF +U = 1/((r2/(h1*r1))+(r2*log(r2/r1)/k)+(1/h2)); // Overall heat transer coeeficient in Btu-hr-ft^2-degF +printf("Overall heat transfer coefficient is %d Btu/hr-ft^2-degF",U); \ No newline at end of file diff --git a/617/CH9/EX9.5/Example9_5.sci b/617/CH9/EX9.5/Example9_5.sci new file mode 100755 index 000000000..8e9341516 --- /dev/null +++ b/617/CH9/EX9.5/Example9_5.sci @@ -0,0 +1,13 @@ +clc(); +clear; + +// To calculate LMTD for heat exchanger + +Tc1 = 120; // Inlet cold fluid temperature in degF +Tc2 = 310; // Outlet cold fluid temperature in degF +Th1 = 500; // Inlet hot fluid temperature in degF +Th2 = 400; // Outlet hot fluid temperature in degF +delt1 = Th2-Tc1; // Maximum temperature difference in degF +delt2 = Th1-Tc2; // Minimum temperature difference in degF +LMTD = (delt1-delt2)/log(delt1/delt2); // Log mean temperature difference +printf("The log mean temperature difference is %d degF",LMTD) diff --git a/617/CH9/EX9.6/Example9_6.sci b/617/CH9/EX9.6/Example9_6.sci new file mode 100755 index 000000000..c5e9b521a --- /dev/null +++ b/617/CH9/EX9.6/Example9_6.sci @@ -0,0 +1,17 @@ +clc(); +clear; + +// To calculate temperature difference for heat exchanger + +Tc1 = 120; // Inlet cold fluid temperature in degF +Tc2 = 310; // Outlet cold fluid temperature in degF +Th1 = 500; // Inlet hot fluid temperature in degF +Th2 = 400; // Outlet hot fluid temperature in degF +K = (Tc2-Tc1)/(Th2-Tc2); // Temperature ratio +R = (Th1-Th2)/(Tc2-Tc1); // Temperature ratio +delt1 = Th2-Tc1; // Maximum temperature difference in degF +delt2 = Th1-Tc2; // Minimum temperature difference in degF +LMTD = (delt1-delt2)/log(delt1/delt2); // Log mean temperature difference +f = 0.99; // Correction factor as seen from figure +LMTDc = round(LMTD*f); // Corrected log mean temperature difference +printf("Log mean temperature difference is %d degF",LMTDc); \ No newline at end of file diff --git a/617/CH9/EX9.7/Example9_7.sci b/617/CH9/EX9.7/Example9_7.sci new file mode 100755 index 000000000..800473333 --- /dev/null +++ b/617/CH9/EX9.7/Example9_7.sci @@ -0,0 +1,22 @@ +clc(); +clear; +// To calculate the outside tube area for a single-pass steam condenser + +Do=1/12; // Outside diameter of the condenser in ft +Di=0.902/12; // Outside diameter of the condenser in ft +Ts=81.7; // Steam temperature in degF +Tw1=61.4; // Water inlet temperature in degF +Tw2=69.9; // Water outlet temperature in degF +k=63; // Thermal conductivity in Btu/hr-ft-degF +v=7; // average velocity in ft/sec +h1=1270; // water side film coefficient i Btu/hr-ft^2-degF +h2=1000; // Steam side film coefficient in Btu/hr-ft^2-degF + +U=1/((Do/(Di*h1))+(Do*log(Do/Di)/(2*k))+(1/h2)); // Heat transfer coefficient +LMTD=((Ts-Tw1)-(Ts-Tw2))/log((Ts-Tw1)/(Ts-Tw2)); // Log mean temperature diff. +m=731300; // Saturated steam to be handled in lb/hr +L=1097.4-49.7; // Change in enthalpy in Btu/lb +q=m*L; // Heat required in Btu/hr +A=q/(U*LMTD); // Area of condenser in ft^2 +printf("The area of steam condenser is %d ft^2",A); + diff --git a/617/CH9/EX9.8/Example9_8.sci b/617/CH9/EX9.8/Example9_8.sci new file mode 100755 index 000000000..3d9d73c4c --- /dev/null +++ b/617/CH9/EX9.8/Example9_8.sci @@ -0,0 +1,53 @@ +clc(); +clear; + +// To calculate overall heat transfer coefficient for heat exchanger + +Tc1 = 139.7; // Inlet cold fluid temperature in degF +Tc2 = 59.5; // Outlet cold fluid temperature in degF +Th1 = 108.7; // Inlet hot fluid temperature in degF +Th2 = 97.2; // Outlet hot fluid temperature in degF +delt1 = Tc1-Th2; // Maximum temperature difference in degF +delt2 = Th1-Tc2; // Minimum temperature difference in degF +LMTD = round((delt1-delt2)/log(delt1/delt2)); +printf(" \n The log mean temperature difference is %d degF",LMTD); + +m = 18210; // Flow rate through tubes +q = m*(Th2-Tc2); // Heat loss in Btu/hr +A = 48.1; // Area in ft^2 +U = q/(A*LMTD); // Overall heat transfer coefficient +printf(" \n The overall heat transfer coefficient is %d Btu/hr-ft^2-degF \n",U); + + +// To calcalute using equations estabilished by correlation +Ts = 113; // Average tube temperature in degF +Tf = (123.9+Ts)/2; // Film temperature in degF +// At this temperature thermal properties are considered +p1 = 61.7/32.2; // Density in slug/ft^3 +u1 = 1.38/32.2; // Viscosity in slug/ft-hr +Cp1 = 1*32.2; // Btu/slug/ft +k1 = 0.366; // Thermal conductivity in Btu/hr-ft-degF +D1 = 0.375/12; // Diameter in ft +v1 = 7610; // Velocity in ft/sec +Nre1 = v1*D1*p1/u1; // Reynolds number +Npr1 = u1*Cp1/k1; // Prandtls number +Nnu1 = 0.33*Nre1^0.6*Npr1^(1/3); // Nusselt number +h1 = Nnu1*k1/D1; // Heat transfer coefficient +printf(" \n The outer heat transfer coefficient is %d Btu/hr-ft^2-degF ",h1); + +// Taking the thermal properties at 78.3 degF +p2 = 62.2/32.2; // Density in slug/ft^3 +u2 = 2.13/32.2; // Viscosity in slug/ft-hr +Cp2 = 1*32.2; // Heat capacity in Btu/slug/ft +k2 = 0.348; // Thermal conductivity in Btu/hr-ft-degF +D2 = 0.277/12; // Diameter in ft +v2 = 7140; // Velocity in ft/sec +Nre2 = v2*D2*p2/u2; // Reynolds number +Npr2 = u2*Cp2/k2; // Prandtls number +Nnu2 = 0.023*Nre2^0.8*Npr2^(0.4); // Nusselt number +h2 = Nnu2*k2/D2; // Heat transfer coefficient +printf(" \n The inner heat transfer coefficient is %d Btu/hr-ft^2-degF",h2); + +k3 = 58; +U1 = 1/((D1/(D2*h2))+(D1*log(D1/D2)/(2*k3))+(1/h1)); // Heat transfer coefficient +printf(" \n The overall heat transfer coefficient accordind to estabilished correlation is %d Btu/hr-ft^2-degF \n",U1); \ No newline at end of file diff --git a/617/CH9/EX9.9/Example9_9.sci b/617/CH9/EX9.9/Example9_9.sci new file mode 100755 index 000000000..8e3f36ba8 --- /dev/null +++ b/617/CH9/EX9.9/Example9_9.sci @@ -0,0 +1,19 @@ +clc(); +clear; + +// To determine the value of product of overall heat transfer and the total area + +To1=140; // inlet temperature of oil in degF +To2=90; // Outlet temperature of oil in degf +Cpo=0.5; // Specific heat capacity in Btu/lb-degf +Tw1=60; // Inlet tempearture of water in degF +Tw2=80; // Outlet temperature of water in degF +mo=2000; // Mass flow rate of oil in lb/hr +q=mo*Cpo*(To1-To2); // Heat transferred in Btu/hr +Cpw=1; // Heat capacity of water in Btu/hr +mw=q/(Cpw*(Tw2-Tw1)); // Mass flow rate in lb/hr +E1=(Tw1-Tw2)/(Tw1-To2); // Effective ratio + +// Seeing the effective ratio and mass flow rate ratio, from the graph we get UA +UA=1.15*mo*Cpo; +printf("The product of overall heat transfer and total area is %d Btu/hr-degF",UA); diff --git a/617/CH9/EX9.9b/Example9_9b.sci b/617/CH9/EX9.9b/Example9_9b.sci new file mode 100755 index 000000000..2918fb05d --- /dev/null +++ b/617/CH9/EX9.9b/Example9_9b.sci @@ -0,0 +1,26 @@ +clc(); +clear; + +// To calculate the temperature of surface and centre plane + +t=2; // Thickness of wall in ft +To=100; // Initial temperature of wall in degF +Tg=1000; // Temperature of hot gases exposed in degF +k=8; // Thermal conductivity in Btu/hr-ft-degF +p=162; // density in lb/ft^-3 +Cp=0.3; // Heat capacity in Btu/lb-degF +h=1.6; // Heat transfer coefficient in Btu/hr-ft^-2-degF +a=k/(p*Cp); // Thermal diffusivity + +// Considering the values of a and 4at/L^2 and hl/2k, the value of Phis, Phic and Si can be obtained +Phis=0.37; +Phic=0.41; +Si=0.62; + +Ta=Tg+(To-Tg)*Phis; // Temperature of surface in degF +printf("The temperature of surface is %d degF \n ",Ta); +Tc=Tg+(To-Tg)*Phic; // Temperature of center plane in degF +printf("The temperature of surface is %d degF \n ",Tc); +A=10; // area of wall through which heat is absorbed +q=p*Cp*t*A*Si*(To-Tg); // Heat absorbed in Btu/hr +printf("The heat absorbed by wall is %d Btu",q); \ No newline at end of file -- cgit