diff options
Diffstat (limited to '914/CH11')
| -rwxr-xr-x | 914/CH11/EX11.1/ex11_1.sce | 29 | ||||
| -rwxr-xr-x | 914/CH11/EX11.10/ex11_10.sce | 62 | ||||
| -rwxr-xr-x | 914/CH11/EX11.2/ex11_2.sce | 22 | ||||
| -rwxr-xr-x | 914/CH11/EX11.3/ex11_3.sce | 26 | ||||
| -rwxr-xr-x | 914/CH11/EX11.5/ex11_5.sce | 58 | ||||
| -rwxr-xr-x | 914/CH11/EX11.6/ex11_6.sce | 23 | ||||
| -rwxr-xr-x | 914/CH11/EX11.7/ex11_7.sce | 28 | ||||
| -rwxr-xr-x | 914/CH11/EX11.8/ex11_8.sce | 23 | ||||
| -rwxr-xr-x | 914/CH11/EX11.9/ex11_9.sce | 15 |
9 files changed, 286 insertions, 0 deletions
diff --git a/914/CH11/EX11.1/ex11_1.sce b/914/CH11/EX11.1/ex11_1.sce new file mode 100755 index 000000000..17592bd85 --- /dev/null +++ b/914/CH11/EX11.1/ex11_1.sce @@ -0,0 +1,29 @@ +clc;
+warning("off");
+printf("\n\n example11.1 - pg497");
+// given
+K_drywall=0.28; //[Btu/ft*degF] - thermal conductivity of dry wall
+K_fibreglass=0.024; //[Btu/ft*degF] - thermal conductivity of fibre glass
+K_concrete=0.5; //[Btu/ft*degF] - thermal conductivity of concrete
+T4=0; //[degF]
+T1=65; //[degF]
+deltaT=T4-T1; //[degF]
+a=1; //[ft^2] - assuming area of 1 ft^2
+deltax1=0.5/12; //[ft]
+deltax2=3.625/12; //[ft]
+deltax3=6/12; //[ft]
+R1=deltax1/(K_drywall*a); //[h*degF/Btu]
+R2=deltax2/(K_fibreglass*a); //[h*degF/Btu]
+R3=deltax3/(K_concrete*a); //[h*degF/Btu]
+qx=deltaT/(R1+R2+R3);
+q12=-qx;
+q23=-qx;
+q34=-qx;
+deltaT1=(-q12)*deltax1*(1/(K_drywall*a));
+T2=T1+deltaT1;
+deltaT2=(-q23)*deltax2*(1/(K_fibreglass*a));
+T3=T2+deltaT2;
+deltaT3=(-q34)*deltax3*(1/(K_concrete*a));
+T4=T3+deltaT3;
+printf("\n\n T1 = %f degF\n T2 = %f degF\n T3 = %f degF\n T4 = %f degF",T1,T2,T3,T4);
+
diff --git a/914/CH11/EX11.10/ex11_10.sce b/914/CH11/EX11.10/ex11_10.sce new file mode 100755 index 000000000..3a9884c47 --- /dev/null +++ b/914/CH11/EX11.10/ex11_10.sce @@ -0,0 +1,62 @@ +clc;
+warning("off");
+printf("\n\n example11.10 - pg544");
+// given
+Ui=325; //[W/m^2*K] - overall heat transfer coefficient
+Thi=120; //[degC] - inlet temperature of hydrocarbon
+Tho=65; //[degC] - outlet temperature of hydrocarbon
+Tci=15; //[degC] - inlet temperature of water
+Tco=50; //[degC] - outlet temperture of water
+cp=4184; //[J/kg*K] - heat capacity of water
+ch=4184*0.45; //[J/kg*K] - heat capacity of hydrocarbon
+wc=1.2; //[kg/sec] - mass flow rate of water
+wh=((wc*cp)*(Tco-Tci))/((ch)*(Thi-Tho));
+qtotal=wc*cp*(Tco-Tci);
+// (a) - parallel double pipe
+F=1;
+Thi=120; //[degC] - inlet temperature of hydrocarbon
+Tho=65; //[degC] - outlet temperature of hydrocarbon
+Tci=15; //[degC] - inlet temperature of water
+Tco=50; //[degC] - outlet temperture of water
+deltaT1=Thi-Tci;
+deltaT2=Tho-Tco;
+LMTD=(deltaT2-deltaT1)/(log(deltaT2/deltaT1));
+Ai=qtotal/((Ui*LMTD));
+printf("\n\n (a) parallel double pipe \n Ai = %f m^2",Ai);
+// (b) - counter flow
+F=1;
+Thi=120; //[degC] - inlet temperature of hydrocarbon
+Tho=65; //[degC] - outlet temperature of hydrocarbon
+Tco=15; //[degC] - inlet temperature of water
+Tci=50; //[degC] - outlet temperture of water
+deltaT1=Thi-Tci;
+deltaT2=Tho-Tco;
+LMTD=(deltaT2-deltaT1)/(log(deltaT2/deltaT1));
+Ai=qtotal/((Ui*LMTD));
+printf("\n\n (b) counter flow \n Ai = %f m^2",Ai);
+// (c) - 1-2 shell and tube
+Thi=120; //[degC] - inlet temperature of hydrocarbon
+Tho=65; //[degC] - outlet temperature of hydrocarbon
+Tci=15; //[degC] - inlet temperature of water
+Tco=50; //[degC] - outlet temperture of water
+Z=(Thi-Tho)/(Tco-Tci);
+nh=(Tco-Tci)/(Thi-Tci);
+deltaT1=Thi-Tco;
+deltaT2=Tho-Tci;
+F=0.92;
+LMTD=(F*(deltaT2-deltaT1))/(log(deltaT2/deltaT1));
+Ai=qtotal/((Ui*LMTD));
+printf("\n\n (c) 1-2 shell and tube \n Ai = %f m^2",Ai);
+// (d) - 2-4 shell and tube
+Thi=120; //[degC] - inlet temperature of hydrocarbon
+Tho=65; //[degC] - outlet temperature of hydrocarbon
+Tci=15; //[degC] - inlet temperature of water
+Tco=50; //[degC] - outlet temperture of water
+Z=(Thi-Tho)/(Tco-Tci);
+nh=(Tco-Tci)/(Thi-Tci);
+F=0.975;
+LMTD=(F*(deltaT2-deltaT1))/(log(deltaT2/deltaT1));
+Ai=qtotal/((Ui*LMTD));
+printf("\n\n (d) 2-4 shell and tube \n Ai = %f m^2",Ai);
+
+
diff --git a/914/CH11/EX11.2/ex11_2.sce b/914/CH11/EX11.2/ex11_2.sce new file mode 100755 index 000000000..7883d0d4b --- /dev/null +++ b/914/CH11/EX11.2/ex11_2.sce @@ -0,0 +1,22 @@ +clc;
+warning("off");
+printf("\n\n example11.2 - pg501");
+// given
+r1=(2.067/2)/(12); //[ft]
+r2=r1+0.154/12; //[ft]
+r3=r2+3/12; //[ft]
+L=1; //[ft]
+Ka=26; //[Btu/h*ft*degF]
+Kb=0.04; //[Btu/h*ft*degF]
+T1=50; //[degF]
+Ra=(log(r2/r1))/(2*%pi*L*Ka);
+Rb=(log(r3/r2))/(2*%pi*L*Kb);
+R=Ra+Rb;
+deltaT=-18; //[degF] - driving force
+Qr=-(deltaT/(R));
+disp(Qr);
+deltaT1=(-Qr)*(Ra);
+T2=T1+deltaT1;
+printf("\n\n The interface temperature is \n T2 = %f degF",T2);
+
+
diff --git a/914/CH11/EX11.3/ex11_3.sce b/914/CH11/EX11.3/ex11_3.sce new file mode 100755 index 000000000..dcd62a9f3 --- /dev/null +++ b/914/CH11/EX11.3/ex11_3.sce @@ -0,0 +1,26 @@ +clc;
+warning("off");
+printf("\n\n example11.3 - pg502");
+// given
+Ra=8.502*10^-4; //[h*degF*Btu^-1]
+Rb=5.014; //[h*degF*Btu^-1]
+r1=(2.067/2)/(12); //[ft]
+r2=r1+0.154/12; //[ft]
+r3=r2+3/12; //[ft]
+d1=2*r1;
+d0=2*r3;
+h0=25; //[Btu/h*ft^2*degF]
+h1=840; //[Btu/h*ft^2*degF]
+L=1; //[ft] - considering 1 feet length
+R0=1/(h0*%pi*d0*L);
+R1=1/(h1*%pi*d1*L);
+R=R0+R1+Ra+Rb;
+disp(R);
+deltaT=-400; //[degF]
+Qr=-(deltaT)/R;
+disp(Qr);
+// the heat loss calculated above is the heat loss per foot.therefore for 500 ft
+L=500;
+Qr=Qr*L;
+printf("\n\n the heat loss for a 500 feet pipe is \n qr = %e Btu/h",Qr);
+
diff --git a/914/CH11/EX11.5/ex11_5.sce b/914/CH11/EX11.5/ex11_5.sce new file mode 100755 index 000000000..fa9d0b66a --- /dev/null +++ b/914/CH11/EX11.5/ex11_5.sce @@ -0,0 +1,58 @@ +clc;
+warning("off");
+printf("\n\n example11.5 - pg521");
+// given
+Nre=50000;
+d=0.04; //[m] - diameter of pipe
+// physical properties of water
+T1=293.15; //[K]
+T2=303.15; //[K]
+T3=313.15; //[K]
+p1=999; //[kg/m^3] - density of water at temperature T1
+p2=996.0; //[kg/m^3] - density of water at temperature T2
+p3=992.1; //[kg/m^3] - density of water at temperature T3
+mu1=1.001; //[cP] - viscosity of water at temperature T1
+mu2=0.800; //[cP] - viscosity of water at temperature T2
+mu3=0.654; //[cP] - viscosity of water at temperature T3
+k1=0.63; //[W/m*K] - thermal conductivity of water at temperature T1
+k2=0.618; //[W/m*K] - thermal conductivity of water at temperature T2
+k3=0.632; //[W/m*K] - thermal conductivity of water at temperature T3
+cp1=4182; //[J/kg*K] - heat capacity of water at temperature T1
+cp2=4178; //[J/kg*K] - heat capacity of water at temperature T2
+cp3=4179; //[J/kg*K] - heat capacity of water at temperature T3
+Npr1=6.94; // prandtl no. at temperature T1
+Npr2=5.41; // prandtl no. at temperature T2
+Npr3=4.32; // prandtl no. at temperature T3
+// (a) Dittus -Boelter-this correction evalutes all properties at the mean bulk temperature,which is T1
+kmb=0.603
+h=(kmb/d)*0.023*((Nre)^(0.8))*((Npr1)^0.4);
+printf("\n\n (a) Dittus -Boelter\n the heat transfer coefficient is \n h = %f W/m^2*K = %f Btu/ft^2*h^-1*degF",h,h*0.17611);
+// (b) Seider Tate-this correlation evaluates all the properties save muw at the mean bulk temperature
+h=(kmb/d)*(0.027)*((Nre)^0.8)*((Npr1)^(1/3))*((mu1/mu3)^0.14);
+printf("\n\n (b) Seider Tate\n the heat transfer coefficient is \n h = %f W/m^2*K = %f Btu/ft^2*h^-1*degF",h,h*0.17611);
+// (c) Sleicher-Rouse equation
+a=0.88-(0.24/(4+Npr3));
+b=(1/3)+0.5*exp((-0.6)*Npr3);
+Nref=Nre*(mu1/mu2)*(p2/p1);
+Nnu=5+0.015*((Nref)^a)*((Npr3)^b);
+h=Nnu*(kmb/d);
+printf("\n\n (c) Sleicher-Rouse equation\n the heat transfer coefficient is \n h = %f W/m^2*K = %f Btu/ft^2*h^-1*degF",h,h*0.17611);
+// (d) Colbum Analogy- the j factor for heat transfer is calculated
+jh=0.023*((Nref)^(-0.2));
+Nst=jh*((Npr2)^(-2/3));
+U=(Nre*mu1*10^-3)/(d*p1);
+h=Nst*(p1*cp1*U);
+printf("\n\n (d) Colbum Analogy\n the heat transfer coefficient is \n h = %f W/m^2*K = %f Btu/ft^2*h^-1*degF",h,h*0.17611);
+// (e) Friend-Metzner
+f=0.005227;
+Nnu=((Nre)*(Npr1)*(f/2)*((mu1/mu3)^0.14))/(1.20+((11.8)*((f/2)^(1/2))*(Npr1-1)*((Npr1)^(-1/3))));
+h=Nnu*(kmb/d);
+printf("\n\n (e) Friend-Metzner\n the heat transfer coefficient is \n h = %f W/m^2*K = %f Btu/ft^2*h^-1*degF",h,h*0.17611);
+// (f) Numerical analysis
+Nnu=320;
+h=Nnu*(kmb/d);
+printf("\n\n (f) Numerical analysis\n the heat transfer coefficient is \n h = %f W/m^2*K = %f Btu/ft^2*h^-1*degF",h,h*0.17611);
+
+
+
+
diff --git a/914/CH11/EX11.6/ex11_6.sce b/914/CH11/EX11.6/ex11_6.sce new file mode 100755 index 000000000..fae0d7b69 --- /dev/null +++ b/914/CH11/EX11.6/ex11_6.sce @@ -0,0 +1,23 @@ +clc;
+warning("off");
+printf("\n\n example11.6 - pg525");
+// given
+Tw=680; //[K] - temperature at the wall
+Tb=640; //[K] - temperature at the bulk
+Tf=(Tw+Tb)/2; //[K]
+Nre=50000;
+vmb=2.88*10^-7;
+vf=2.84*10^-7;
+Nref=Nre*(vmb/vf);
+k=27.48;
+d=0.04;
+// from table 11.3 the prandtl no. is
+Npr=8.74*10^-3
+// constant heat flow
+Nnu=6.3+(0.0167)*((Nref)^0.85)*((Npr)^0.93);
+h=Nnu*(k/d);
+printf("\n\n constant heat flow\n h = %f W/m^2*K = %f Btu/ft^2*h*degF",h,h*0.17611);
+// constant wall temperature
+Nnu=4.8+0.0156*((Nref)^0.85)*((Npr)^0.93);
+h=Nnu*(k/d);
+printf("\n\n constant wall temperature\n h = %f W/m^2*K = %f Btu/ft^2*h*degF",h,h*0.17611);
diff --git a/914/CH11/EX11.7/ex11_7.sce b/914/CH11/EX11.7/ex11_7.sce new file mode 100755 index 000000000..bda1cbc8e --- /dev/null +++ b/914/CH11/EX11.7/ex11_7.sce @@ -0,0 +1,28 @@ +clc;
+warning("off");
+printf("\n\n example11.7 - pg536");
+// given
+di=0.620; //[inch] - internal diameter
+d0=0.750; //[inch] - outer diameter
+Ai=0.1623; //[ft^2/ft]
+Ao=0.1963; //[ft^2/ft]
+wc=12*(471.3/0.9425);
+cp=1; //[Btu/lbm*degF] - heat capacity of water
+Tco=110;
+Tci=50;
+qtotal=wc*cp*(Tco-Tci);
+deltaH_coldwater=3.6*10^5;
+deltaH_vapourization=1179.7-269.59;
+wh=deltaH_coldwater/deltaH_vapourization;
+hi=80; //[Btu/h*ft^2*degF]
+ho=500; //[Btu/h*ft^2*degF]
+km=26; //[Btu/h*ft*degF]
+Ui=1/((1/hi)+((Ai*log(d0/di))/(2*%pi*km))+(Ai/(Ao*ho)));
+disp(Ui)
+deltaT1=300-50;
+deltaT2=300-110;
+LMTD=(deltaT1-deltaT2)/(log(deltaT1/deltaT2));
+A=qtotal/(Ui*LMTD);
+L=A/Ai;
+printf("\n\n the length of the heat exchanger is \n L = %f ft",L);
+
diff --git a/914/CH11/EX11.8/ex11_8.sce b/914/CH11/EX11.8/ex11_8.sce new file mode 100755 index 000000000..9d6218535 --- /dev/null +++ b/914/CH11/EX11.8/ex11_8.sce @@ -0,0 +1,23 @@ +clc;
+warning("off");
+printf("\n\n example11.8 - pg537");
+// given
+L=30; //[ft] - length
+Ai=0.1623*L;
+di=0.620; //[inch] - internal diameter
+d0=0.750; //[inch] - outer diameter
+Ao=0.1963*L; //[ft^2/ft]
+wc=12*(471.3/0.9425);
+cp=1; //[Btu/lbm*degF] - heat capacity of water
+deltaH_coldwater=3.6*10^5;
+deltaH_vapourization=1179.7-269.59;
+wh=deltaH_coldwater/deltaH_vapourization;
+hi=80; //[Btu/h*ft^2*degF]
+ho=500; //[Btu/h*ft^2*degF]
+km=26; //[Btu/h*ft*degF]
+Ui=1/((1/hi)+(((Ai/L)*log(d0/di))/(2*%pi*km))+(Ai/(Ao*ho)));
+deltaT1=300-50;
+deltaT=deltaT1/(exp((Ui*Ai)/(wc*cp)));
+Tsat=300;
+Tc2=Tsat-deltaT;
+printf("\n\n Therefore, the outlet temperature of the cold fluid is \n Tc2 = %f degF",Tc2);
diff --git a/914/CH11/EX11.9/ex11_9.sce b/914/CH11/EX11.9/ex11_9.sce new file mode 100755 index 000000000..44c318004 --- /dev/null +++ b/914/CH11/EX11.9/ex11_9.sce @@ -0,0 +1,15 @@ +clc;
+warning("off");
+printf("\n\n example11.9 - pg538");
+// given
+Ai=4.869;
+wc=6000;
+cp=1;
+Rf=0.002;
+Uclean=69.685;
+Udirty=1/(Rf+(1/Uclean));
+deltaT1=300-50;
+deltaT2=deltaT1/(exp((Udirty*Ai)/(wc*cp)));
+Th2=300;
+Tc2=Th2-deltaT2;
+printf("\n\n the outlet temperature is \n Tc2 = %f degF",Tc2);
|