5 files changed, 165 insertions, 0 deletions
diff --git a/1328/CH19/EX19.1/19_1.sce b/1328/CH19/EX19.1/19_1.sce
new file mode 100644
index 000000000..4d0d9851d
--- /dev/null
+++ b/1328/CH19/EX19.1/19_1.sce
@@ -0,0 +1,74 @@
+printf("\t example 19.1 \n");
+// For orientation purposes, one can make an estimate of the number of tubes required in the radiant section by assuming avg flux is 12000 Btu/(hr)*(ft^2)
+// from Fig.19.14 it can be seen that with a tube temperature of 800"F, an exit-gas temperature of l730°F will be required to effect such a flux.
+printf("\t approxiate values are mentioned in the book \n");
+Q=50000000; // Btu/hr
+QF=(Q/0.75); // efficiency of tank is 75%
+printf("\t heat liberated by the fuel : %.3e Btu/hr \n",QF);
+w1=(QF/17130); // heating value of fuel is 17130Btu/lb
+printf("\t fuel quantity : %.2e lb/hr \n",w1);
+w2=(w1*17.44); // lb of fuel fired with 17.44lb of air
+printf("\t air required : %.2e lb/hr \n",w2);
+w3=(w1*0.3); // 0.3 lb of air is used for atomizing lb of fuel
+printf("\t steam for atomizing : %.2e lb/hr \n",w3);
+QA=(w2*82); // heating value at 400F is 82Btu/lb
+printf("\t QA is : %.2e Btu/hr \n",QA);
+printf("\t QS is negligible \n");
+QW=(0.02*QF);
+printf("\t QW is : %.2e Btu/hr \n",QW);
+Qnet=(QF+QA-QW);
+printf("\t Qnet is : %.2e Btu/hr \n",Qnet);
+//Heat out m gases at 1730°F, 25 per cent excess air, 476 Btu/lb of flue gas
+QG=(476*(w1+w2+w3));
+printf("\t QG is : %.2e Btu/hr \n",QG);
+Q1=(Qnet-QG);
+printf("\t Q1 is : %.2e Btu/hr \n",Q1); // calculation mistake in book
+A=(3.14*38.5*(5/12)); // area of tube
+printf("\t area of tube is : %.1f ft^2 \n",A);
+Nt=(Q1/(12000*A)); // 12000 is avg flux
+printf("\t estimated number of tubes : %.0f \n",Nt);
+// The layout of the cross section of the furnace may be as shown m Fig. 19.16.
+// center to center distance is 8(1/2)in
+Acp=(8.5*38.5/12);
+printf("\t cold plane surface per tube : %.1f ft^2 \n",Acp); // calculation mistake in book
+a=0.937; // a=alpha, from fig 19.11 as Ratio of center-to-center/OD is 1.7
+Acp1=(Acp*a);
+printf("\t Acp1 is : %.0f ft^2 \n",Acp1);
+Acpt=(Acp1*Nt);
+printf("\t total cold plane surface is : %.1e ft^2 \n",Acpt);
+A1=(2*20.46*14.92); // from fig 19.16
+printf("\t surface of end walls : %.0f ft^2 \n",A1);
+A2=(38.5*14.92); // from fig 19.16
+printf("\t surface of side wall : %.0f ft^2 \n",A2);
+A3=(38.5*9.79); // from fig 19.16
+printf("\t surface of bridge walls : %.0f ft^2 \n",A3);
+A4=(2*20.46*38.5); // from fig 19.16
+printf("\t surface of floor and arch : %.0f ft^2 \n",A4);
+AT=(A1+A2+A3+A4);
+printf("\t AT is : %.0f ft^2 \n",AT);
+AR=(AT-Acpt);
+printf("\t AR is : %.0f ft^2 \n",AR);
+Ar=(AR/Acpt);
+printf("\t ratio of areas is : %.2f \n",Ar);
+printf("\t dimension ratio is 3:2:1 \n");
+L=((2/3)*(38.5*20.46*14.92)^(1/3));
+printf("\t length is : %.0f ft \n",L);
+printf("\t gas emissivity \n");
+// From the analysis of the fuel, the steam quantity, and the assumption that the humidity of the air is 50 per cent of saturation at 60F, the partial pressures of CO2 and H2O in the combustion gases with 25 per cent excess air are
+pCO2=0.1084;
+pH2O=0.1248
+pCO2L=1.63; // pCO2L=(pCO2*L)
+pH2OL=1.87;
+P=((pCO2)/(pCO2+pH2O));
+printf("\t percentage correction at P : %.3f \n",P);
+Pt=pCO2L+pH2OL;
+printf("\t Pt is : %.2f \n",Pt);
+// %correction estimated to be 8%
+eG=(((6500+14500)-(650+1950))/(39000-4400))*((100-8)/100); // values from fig 19.12 and 19.13, eq 19.5
+printf("\t eG is : %.3f \n",eG);
+f=0.635; // from fig 19.15 as (AR/Acpt)=1.09 and eG=0.496
+printf("\t overall exchange factor : %.3f \n",f);
+Z=(Q1/(Acpt*f));
+printf("\t Z is : %.2e \n",Z);
+printf("\t TG required (at Ts = 800F) = 1670F compared with 1730°F assumed in heat balance) \n");
+// end
diff --git a/1328/CH19/EX19.2/19_2.sce b/1328/CH19/EX19.2/19_2.sce
new file mode 100644
index 000000000..79bf17152
--- /dev/null
+++ b/1328/CH19/EX19.2/19_2.sce
@@ -0,0 +1,9 @@
+printf("\t example 19.2 \n");
+QF=50000000;
+G=22.36;
+Acpt=1500;
+printf("\t approxiate values are mentioned in the book \n");
+Q=(QF/(1+(G/4200)*(QF/Acpt)^(1/2))); // eq 19.15
+printf("\t Q is : %.2e Btu/hr \n",Q);
+printf("\t The radiant-section average rate will be 8350 Btu/(hr) (ft2), and the exit-flue-gas temperature 1540°F by heat balance. \n");
+// end
diff --git a/1328/CH19/EX19.3/19_3.sce b/1328/CH19/EX19.3/19_3.sce
new file mode 100644
index 000000000..70aeb18eb
--- /dev/null
+++ b/1328/CH19/EX19.3/19_3.sce
@@ -0,0 +1,17 @@
+printf("\t example 19.3 \n");
+Qr=1.5; // Qr=(QF2/QF1)
+Cr=1.5; // Cr=(CR2/CR1)
+Gr=140/125; // Gr=(G2/G1)
+Qr1=0.38; // Qr1=(Q1/QF1)
+printf("\t approxiate values are mentioned in the book \n");
+a1=1.63; // a1=(G1*(CR1/27)^(1/2)), from eq 19.17
+printf("\t a1 is : %.2f \n",a1);
+a2=1.37*(a1); // a2=(G2*(CR2/27)^(1/2))
+printf("\t a2 is : %.2f \n",a2);
+Qr2=(1/(1+a2)); // Qr2=(Q2/QF2),from eq 19.15
+printf("\t Qr2 is : %.2f \n",Qr2);
+Q21=(Qr2/Qr1)*(Qr); // Q21=(Q2/Q1)
+printf("\t ratio of heats is : %.2f \n",Q21);
+printf("\t Hence the radiant absorption will be increased only 22 per cent for an increase of 50 per cent in the heat liberated. \n");
+// end
+
diff --git a/1328/CH19/EX19.4/19_4.sce b/1328/CH19/EX19.4/19_4.sce
new file mode 100644
index 000000000..2c7e195d1
--- /dev/null
+++ b/1328/CH19/EX19.4/19_4.sce
@@ -0,0 +1,35 @@
+printf("\t example 19.4 \n");
+eS=0.9; // assumed
+TG=1500;
+TS=650;
+pCO2=0.1084;
+pH2O=0.1248;
+printf("\t approxiate values are mentioned in the book \n");
+L=(0.4*8.5)-(0.567*5);  // table 19.1
+printf("\t L is : %.3f ft \n",L);
+pH2OL=0.1248*L;
+pCO2L=0.1084*L;
+printf("\t pH2OL is : %.4f atm-ft \n",pH2OL);
+printf("\t pCO2L is : %.4f atm-ft \n",pCO2L);
+qH2O=1050; // at TG, from fig 19.12 ana 19.13
+qCO2=1700; // at TG, from fig 19.12 ana 19.13
+qTG=(qH2O+qCO2);
+printf("\t qTG is : %.0f \n",qTG);
+qsH2O=165; // at TS, from fig 19.12 ana 19.13
+qsCO2=160; // at TS, from fig 19.12 ana 19.13
+qTS=(qsH2O+qsCO2);
+printf("\t qTG is : %.0f \n",qTS);
+q=(0.9*(qTG-qTS)); // q=(QRC/A)
+printf("\t q is : %.1f \n",q);
+P=((pCO2)/(pCO2+pH2O));
+printf("\t percentage correction at P : %.3f \n",P);
+Pt=pCO2L+pH2OL;
+printf("\t Pt is : %.4f \n",Pt);
+// %correction estimated to be 2%
+q1=(q*0.98); // // q1=(QRC/A)
+printf("\t q1 is : %.2e \n",q1);
+hr=(q1/(TG-TS));
+printf("\t radiation coefficient is : %.2f Btu/(hr)*(ft^2)*(F) \n",hr);
+//end
+
+
diff --git a/1328/CH19/EX19.5/19_5.sce b/1328/CH19/EX19.5/19_5.sce
new file mode 100644
index 000000000..1cf416bcf
--- /dev/null
+++ b/1328/CH19/EX19.5/19_5.sce
@@ -0,0 +1,30 @@
+printf("\t example 19.5 \n");
+Q=500000;
+printf("\t approxiate values are mentioned in the book \n");
+a=(3.5+(3.14*4*(120/360)))/(2); // a=(alpha*Acp) from fig 19.17
+AR=(3+3.6+3);
+printf("\t a is : %.2f ft^2/ft \n",a);
+printf("\t AR is : %.1f ft^2/ft \n",AR);
+// Arbitrarily neglecting end wa.lls and also .the side wall refractory over 3'0" above the floor
+R=(AR/a);
+printf("\t ratio of two areas is : %.2f \n",R);
+eG=0.265;
+TG=1174; // F
+TS=500; // F
+f=0.56; // from fig 19.15 as (AR/Acpt)=2.49 and eG=0.265
+q=15300; // at TG and TS,q=(Q/(a*f))
+// However, the convection coefficient is small, 1.0 ± Btu/(hr)(ft2)("F), and AR/a is not 2.0 as in the assumptions for the Lobo and Evans equation.
+q1=(q)-(7*(TG-TS)); // q1=(Q/(a*f))
+printf("\t q1 is : %.2e Btu/(hr)*(ft^2) \n",q1);
+q2=(q1*f); // q2=(Q/(a))
+printf("\t q2 is : %.2e Btu/(hr)*(ft^2) \n",q2);
+printf("\t convection rate basis \n");
+q3=(1*(TG-TS)*(4.2/a)); // q2=(Q/(a))
+printf("\t q3 is : %.1e Btu/(hr)*(ft^2) \n",q3); // calculation mistake in book
+qt=(q2+q3); // qt=(Q/(a))
+printf("\t qt is : %.2e Btu/(hr)*(ft^2) \n",qt);
+ar=(Q/qt);
+printf("\t required a is : %.0f ft^2 \n",ar);
+L=(ar/a);
+printf("\t length required is : %.1f ft \n",L);
+// end