initial commit / add all books

10 files changed, 240 insertions, 0 deletions

diff --git a/2471/CH1/EX1.1/Ex1_1.sce b/2471/CH1/EX1.1/Ex1_1.sce
new file mode 100755
index 000000000..1958fb9df
--- /dev/null
+++ b/2471/CH1/EX1.1/Ex1_1.sce
@@ -0,0 +1,16 @@
+
+clc;
+// Example 1.1
+printf('Example 1.1\n\n');
+printf(' Page No. 08\n\n');
+// Solution 
+
+// Given
+m1= 40*10^3;// fuel oil in gallons per year
+ga= 4.545*10^-3;// m^3
+m= m1*ga;// fuel oil in m^3 per year
+Cv1= 175*10^3;// Btu per gallons
+Bt= .2321*10^6;// J per m^3
+Cv= Cv1*Bt;// in J per year per m^3 
+q=m*Cv;// in J per year
+printf(' Heat available is %3.2e J per year\n',q)
diff --git a/2471/CH1/EX1.10/Ex1_10.sce b/2471/CH1/EX1.10/Ex1_10.sce
new file mode 100755
index 000000000..69380bdaf
--- /dev/null
+++ b/2471/CH1/EX1.10/Ex1_10.sce
@@ -0,0 +1,57 @@
+clear ;

+clc;

+// Example 1.10

+printf('Example 1.10\n\n');

+printf('Page No. 19\n\n');

+// given

+

+//Monthly Energy Usage

+qunty = [15*10^3 4*10^3 90*10^3]

+cost = [4950 720  2250]// in Pound

+common_basis1  = [738 117 90]// in 10^6 Wh

+common_basis= [2655 421 324]// converted into 10^9 Joules

+unit_cost = cost ./ common_basis1// in Pound per 10^6 Wh

+p= 80;// production in tonnes

+EI = ((sum(common_basis))/p)*10^9;

+CI = sum(cost)/80;

+printf('Monthly energy index is %3.2e J per tonne \n',EI)

+printf('Monthly cost index is %.0f Pound per tonne \n\n',CI)//Deviation in answer is due to calculation error for sum of cost in the book

+

+// Boiler House Energy Audit

+qunty_b = [15000 10000]

+Com_basis_b_1 = [2655 36]// in 10^9 J

+Com_basis_b = [738 10]// in 10^6 Wh

+Cost_b = [4950 250]// in Pound

+b_output = 571*10^6;// in Wh

+EI_b = (b_output/(sum(Com_basis_b)*10^6));

+CI_b = (sum(Cost_b)/b_output)*10^3;// Pound converted into p

+printf('Energy index  for boiler is %.3f \n',EI_b)

+printf('Cost index for boiler is %3.2e p per Wh\n \n',CI_b)

+

+//Power House Energy Audit

+P_gen = 200*10^6;// Power generated in Wh

+Com_basis_p_1 = [14.4 2055 -1000]// in 10^9 J

+Com_basis_p = [4.0 571 -278]// in 10^6 Wh

+Cost_p = [100 5196 -2530]// in Pound

+CI_p = (sum(Cost_p)/P_gen)*10^3;// Pound converted into p

+printf('Cost index for power house is %3.2e p per Wh\n\n',CI_p)//Deviation in answer is due to wrong calculation in the book

+

+//Space Heating Energy Audit

+deg_days = 260;// Number of degree-days

+Com_basis_s_1 = [36 100 105]// in 10^9 J

+Com_basis_s = [10.0 27.8 29.2]// in 10^6 Wh

+Cost_s = [250 253 179]// in Pound

+EI_s = ((sum(Com_basis_s)*10^6)/deg_days)

+CI_s = (sum(Cost_s)/deg_days)

+printf('Energy index  for space heating is %3.2e Wh per degree-day\n',EI_s)

+printf('Cost index for space heating is %3.2f Pound per degree-day\n\n',CI_s)

+

+//Process Energy Audit

+T_pdt_output = 100;// in tonne

+Com_basis_pr_1 = [216 720 810 316]// in 10^9 J

+Com_basis_pr = [60 200 225 88]// in 10^6 Wh

+Cost_pr = [1500 2766 2047 540]// in Pound

+EI_pr = ((sum(Com_basis_pr)*10^6)/T_pdt_output);

+CI_pr = (sum(Cost_pr)/T_pdt_output);

+printf('Energy index  for Process Energy Audit is %3.2e Wh per tonne \n',EI_pr)

+printf('Cost index for Process Energy Audit is %.2f Pound per tonne \n',CI_pr)

diff --git a/2471/CH1/EX1.2/Ex1_2.sce b/2471/CH1/EX1.2/Ex1_2.sce
new file mode 100755
index 000000000..2824a4c11
--- /dev/null
+++ b/2471/CH1/EX1.2/Ex1_2.sce
@@ -0,0 +1,25 @@
+clear ;

+clc;

+// Example 1.2

+printf('Example 1.2\n\n');

+printf('Page No. 09\n\n');

+// Solution 

+

+// Given

+Eo1= 1.775*10^9;// Annular energy consumption of oil in Btu

+Btu= 1055;// 1 Btu = 1055 Joules

+Eo= Eo1*Btu;// Annular energy consumption of oil in Joules

+Eg1= 5*10^3;// Annular energy consumption of gas in Therms

+Th= 1055*10^5;//  1 Th = 1055*10^3 Joules

+Eg= Eg1*Th;// Annular energy consumption of gas in Joules

+Ee1= 995*10^3;// Annular energy consumption of electricity in KWh

+KWh= 3.6*10^6;// 1 KWh = 3.6*10^6 Joules

+Ee= Ee1*KWh;// Annular energy consumption of electricity in Joules

+Et= ( Eo + Eg + Ee);// Total energy consumption

+P1= (Eo/Et)*100; // percentage of oil consumption

+P2= (Eg/Et)*100; // percentage of gas consumption

+P3= (Ee/Et)*100; // percentage of electricity consumption

+printf('percentage of oil consumption is %3.1f \n',P1)

+printf('percentage of gas consumption is %3.1f \n',P2)

+printf('percentage of electricity consumption is %3.1f \n',P3)

+

diff --git a/2471/CH1/EX1.3/Ex1_3.sce b/2471/CH1/EX1.3/Ex1_3.sce
new file mode 100755
index 000000000..e35a408b4
--- /dev/null
+++ b/2471/CH1/EX1.3/Ex1_3.sce
@@ -0,0 +1,22 @@
+clear ;

+clc;

+// Example 1.3

+printf('Example 1.3\n\n');

+printf('Page No. 10\n\n');

+// Solution 

+

+// Given

+Et = 100*10^3;// total energy production in tonnes per annum

+Eo= 0.520*10^9;// oil consumption in Wh

+Eg= 0.146*10^9;// gas consumption in Wh

+Ee= 0.995*10^9;// electricity consumption in Wh

+Io= Eo/Et;

+Ig= Eg/Et;

+Ie= Ee/Et;

+Et1= Eo + Eg + Ee;// total energy consumption

+It= Et1/Et;

+printf('oil energy index is %3.0f Wh per tonne \n',Io)

+printf('gas energy index is %3.0f Wh per tonne \n',Ig)

+printf('electricity energy index is %3.0f Wh per tonne \n',Ie)

+printf('total energy index is %3.0f Wh per tonne ',It)

+

diff --git a/2471/CH1/EX1.4/Ex1_4.sce b/2471/CH1/EX1.4/Ex1_4.sce
new file mode 100755
index 000000000..5b515e949
--- /dev/null
+++ b/2471/CH1/EX1.4/Ex1_4.sce
@@ -0,0 +1,22 @@
+clear ;

+clc;

+// Example 1.4

+printf('Example 1.4\n\n');

+printf('Page No. 10\n\n');

+// Solution 

+

+// Given

+mc= 1.5*10^3;// coke consumption in tonnes

+mg= 18*10^3;// gas consumption in therms

+me= 1*10^9;// electricity consumption in Wh

+Cc1= 72;// cost of coke in Pound per tonne

+Cg1= 0.20;// cost of gas in Pound per therm

+Ce1= 2.25*10^-5 ;// cost of electricity in Pound per Wh

+Cc= mc*Cc1;//in Pound

+Cg= mg*Cg1;//in Pound

+Ce= me*Ce1;//in Pound

+Ct= Cc+Cg+Ce;//in Pound

+printf('cost of coke consumption is %.0f Pound \n',Cc)

+printf('cost of gas consumption is %.0f Pound \n',Cg)

+printf('cost of electricity consumption is %.0f Pound \n',Ce)

+printf('total cost is %3.0f Pound \n',Ct)

diff --git a/2471/CH1/EX1.5/Ex1_5.sce b/2471/CH1/EX1.5/Ex1_5.sce
new file mode 100755
index 000000000..3b862ccdf
--- /dev/null
+++ b/2471/CH1/EX1.5/Ex1_5.sce
@@ -0,0 +1,21 @@
+clear ;

+clc;

+// Example 1.5

+printf('Example 1.5\n\n');

+printf('Page No. 11\n\n');

+// Solution 

+

+// Given

+Cc= 108.0*10^3;// cost of coke in Pound

+Cg= 3.6*10^3;// cost of gas in Pound

+Ce= 22.5*10^3;// cost of electricity in Pound

+Ct= Cc+Cg+Ce;// total cost of fuel in Pound

+E= 15*10^3;// total production in tonnes per year

+Ic= Cc/E;//Pound per tonne

+Ig= Cg/E;//Pound per tonne

+Ie= Ce/E;//Pound per tonne

+It= Ct/E;//Pound per tonne

+printf(' coke cost index is %3.2f Pound per tonne \n',Ic)

+printf(' gas cost index is %3.2f Pound per tonne\n',Ig)

+printf(' electricity cost index is %3.2f Pound per tonne\n',Ie)

+printf(' total cost index is %3.2f Pound per tonne\n',It)

diff --git a/2471/CH1/EX1.6/Ex1_6.sce b/2471/CH1/EX1.6/Ex1_6.sce
new file mode 100755
index 000000000..2cc062b72
--- /dev/null
+++ b/2471/CH1/EX1.6/Ex1_6.sce
@@ -0,0 +1,19 @@
+clear ;
+clc;
+// Example 1.6
+printf('Example 1.6\n\n');
+printf('Page No. 11\n\n');
+// Solution 
+
+// Given
+G1= 11.72*10^3;// hourly consumption of gas in therms
+th= 34.13;// in Watts
+G= G1*th;// hourly consumption of gas in Watts
+O1= 4.32*10^9;// hourly consumption of oil in Joules
+J= .278*10^-3;// in Watts
+O= O1*J;// hourly consumption of oil in Watts
+E= 500*10^3;// hourly consumption of electricity in Watts
+// Pie Chart Representation : one input argument  x=[G O E]
+pie([G O E],["gas" "oil" "electricity"]);// Please see the graphics window
+printf('The Pie chart is plotted in the figure');
+
diff --git a/2471/CH1/EX1.7/Ex1_7.sce b/2471/CH1/EX1.7/Ex1_7.sce
new file mode 100755
index 000000000..1dc8ea9fe
--- /dev/null
+++ b/2471/CH1/EX1.7/Ex1_7.sce
@@ -0,0 +1,17 @@
+close

+clear ;

+clc;

+// Example 1.7

+printf('Example 1.7\n\n');

+printf('Page No. 12\n\n');

+// Solution 

+

+// Given

+O= 150*10^3;// energy consumption in office heating in Watts

+L= 120*10^3;// energy consumption in lighting in Watts

+B= 90*10^3;// energy consumption in boiler house in Watts

+P= 180*10^3;// energy consumption in process in Watts

+// Pie Chart Representation : one input argument  x=[O L B P]

+pie([O L B P],["office heating" "lighting" "boiler heating" "process"]);// Please see the graphics window

+printf('The Pie chart is plotted in the figure');

+

diff --git a/2471/CH1/EX1.8/Ex1_8.sce b/2471/CH1/EX1.8/Ex1_8.sce
new file mode 100755
index 000000000..a56f83edc
--- /dev/null
+++ b/2471/CH1/EX1.8/Ex1_8.sce
@@ -0,0 +1,17 @@
+clear ;

+clc;

+// Example 1.8

+printf('Example 1.8\n\n');

+printf('Page No. 16\n\n');

+// given

+

+qunty= [40 10000 400 90000]

+unit_price= [29 0.33 0.18 0.025]

+cost= (unit_price .* qunty)// in Pound

+common_basis= [310 492 11.7 90]// in 10^6 Wh

+per_unit_cost= (unit_price .* qunty) ./ common_basis// Pound per 10^6 Wh

+p= 150;// production in tonnes

+EI= sum(common_basis)*10^6/150

+CI= sum(unit_price .* qunty)/150

+printf('energy index is %3.2f Wh per tonne \n',EI)

+printf('cost index  is %3.2f Wh per tonne \n',CI)

diff --git a/2471/CH1/EX1.9/Ex1_9.sce b/2471/CH1/EX1.9/Ex1_9.sce
new file mode 100755
index 000000000..e7d6fd5d7
--- /dev/null
+++ b/2471/CH1/EX1.9/Ex1_9.sce
@@ -0,0 +1,24 @@
+clear ;

+clc;

+// Example 1.9

+printf('Example 1.9\n\n');

+printf('Page No. 17\n\n');

+//given

+

+p= [50, 55, 65, 50, 95, 90, 85, 80, 60, 90, 70, 110, 60, 105];// weakly production in tonnes

+s= [0.4, 0.35, 0.45, .31, 0.51,0.55, 0.45, 0.5, 0.4, 0.51, 0.4, 0.6, 0.45, 0.55];// weakly steam consumption in 10^6 kg

+coefs = regress(p,s);

+new_p = 0:120

+new_s = coefs(1) + coefs(2)*new_p;

+plot(p,s,'r*');

+set(gca(),"auto_clear","off")

+

+plot(new_p,new_s);// please see the corresponding graph in graphic window

+xtitle('weakly steam consumption-production','weakly output (tonnes)','steam consumption/week (10^6 kg)')

+l = legend([_('Given data'); _('Fitting function')],2);

+

+in= coefs(1)*10^6;// intercept of graph in kg/weak

+printf('At zero output the steam consumption is %3.0f in kg/weak \n',in)

+

+

+