MATLAB Code:
The MATLAB programming code to demonstrate the simulation results of defined objectives 1-7 is given as follows.
%==========================================================================
% Experoment on Steam Power Station
% By: Prof Walid Al-Hussaibi
% Class: 4th Year ; Date: November 2023
%--------------------------------------------------------------------------
clc;clear;close all
% Power Station Data
Cv = 8000; % Fuel calorific value
LF = 0.4 ; % Load factor of the power station
Eboiler = 0.42; % Boiler efficiency
Eturbine = 0.80; % Turbine efficiency
Eelectric = 0.90; % Boiler efficiency
Cost = 35 ; % Cost of coal in USA $ per 1 ton
Mdemand = [36000 85500 72000 68400 90750 50900 40000]; % Max Demand in kW
Week = 1:7; % The week days, i.e. [Sun Mon Tue Wed Thu Fri Sat]
%--------------------------------------------------------------------------
%----------- Objective 1:
% The plot of daily generated power units(Po)in kWh over a period of one
% week and the average.
Po = Mdemand*LF*24 % Generated power units per day in kWh Pavg
(1:7) = sum(Po)/7 % Average power per day in kWh
figure; hold on;
plot(Week,Po,'-b*','MarkerSize',8,'linewidth',1.2);hold on;
plot(Week,Pavg,'--ro','MarkerSize',6,'linewidth',1.2); hold off;
xlabel('Week Days'); ylabel('Power in kWh');grid;box legend('Total
Power','Average Power');
title('Daily Generated Power in kWh')
@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@@
%----------- Objective 2:
% The plot of daily produced heat(H)in kcal over a period of one
% week and the average. Note that 1 kWh = 860 kcal
Eoverall = Eboiler*Eturbine*Eelectric ; % Overall efficiency
H = (Po*860)/Eoverall % Heat of combustion per day in kcal Havg
(1:7) = sum(H)/7 % Average heat value per day in kcal
figure; hold on;
plot(Week,H,'-bp','MarkerSize',8,'linewidth',1.2);hold on;
plot(Week,Havg,'--ro','MarkerSize',6,'linewidth',1.2); hold off;
xlabel('Week Days'); ylabel('Heat Value in kcal');grid;box legend('Total
Heat','Average Heat');
title('Daily Produced Heat in kcal')
%----------- Objective 3:
% The plot of daily consumed coal in tons over a period of one week and the
% average.
Ccoal = H/Cv % Consumed coal per day in tons AvgCcoal
(1:7)=sum(Ccoal)/7 % Average coal consumption per day
figure; hold on;
plot(Week,Ccoal,'-bs','MarkerSize',8,'linewidth',1.2);hold on;
plot(Week,AvgCcoal,'--ro','MarkerSize',6,'linewidth',1.2); hold off;
xlabel('Week Days'); ylabel('Coal Consumption in tons');grid;box legend('Total
Consumption','Average Consumption');
title('Daily Coal Consumption in Tons')
%----------- Objective 4:
% The specific coal consumption in kg/kWh.
SpecCoalConsumption = Ccoal(1)/Po(1)
%----------- Objective 5:
% The amount of yearly fuel cost in USA $. CoalCost1year
= Cost*52*sum(Ccoal)/1000
%----------- Objective 6:
% The average load on the station in kW. AvgLoad =
sum(Po)/(24*7)
%----------- Objective 7:
% The capacity of installed generator in MW.
GeneratorCapacity = max(Mdemand) % More than this value must be considered
%-------- END of the Code -------------------------------------------------
%==========================================================================
