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 --- 2471/CH1/EX1.1/Ex1_1.sce | 16 +++++++++++++ 2471/CH1/EX1.10/Ex1_10.sce | 57 ++++++++++++++++++++++++++++++++++++++++++++++ 2471/CH1/EX1.2/Ex1_2.sce | 25 ++++++++++++++++++++ 2471/CH1/EX1.3/Ex1_3.sce | 22 ++++++++++++++++++ 2471/CH1/EX1.4/Ex1_4.sce | 22 ++++++++++++++++++ 2471/CH1/EX1.5/Ex1_5.sce | 21 +++++++++++++++++ 2471/CH1/EX1.6/Ex1_6.sce | 19 ++++++++++++++++ 2471/CH1/EX1.7/Ex1_7.sce | 17 ++++++++++++++ 2471/CH1/EX1.8/Ex1_8.sce | 17 ++++++++++++++ 2471/CH1/EX1.9/Ex1_9.sce | 24 +++++++++++++++++++ 10 files changed, 240 insertions(+) create mode 100755 2471/CH1/EX1.1/Ex1_1.sce create mode 100755 2471/CH1/EX1.10/Ex1_10.sce create mode 100755 2471/CH1/EX1.2/Ex1_2.sce create mode 100755 2471/CH1/EX1.3/Ex1_3.sce create mode 100755 2471/CH1/EX1.4/Ex1_4.sce create mode 100755 2471/CH1/EX1.5/Ex1_5.sce create mode 100755 2471/CH1/EX1.6/Ex1_6.sce create mode 100755 2471/CH1/EX1.7/Ex1_7.sce create mode 100755 2471/CH1/EX1.8/Ex1_8.sce create mode 100755 2471/CH1/EX1.9/Ex1_9.sce (limited to '2471/CH1') 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) + + + -- cgit