Project 2 - Rankine cycle Simulator

Aim : To simulate rankine cycle

Process

1-2 Isentropic compresion in pump

2-3 Constant pressure heat addition in boiler

3-4 Isentropic expansion in turbine

4-1 Constant pressure heat rejection in condeser

clear all
close all
clc

disp('RANKINE CYCLE SIMULATOR');
disp('1-2 = isentropic expansion in the turbine');
disp('2-3 = constant pressure heat rejection');
disp('3-4 = isentropic compression in the pump');
disp('4-1 = contant pressure heat addition in the boiler');

%to type inputs and get the text given in the problem

t1 = input('enter the temperature at the turbine inlet(in degree celcius)');

p1 = input('enter the pressure at turbine inlet(i bar)');
p2 = input('enter the pressure at the condensor(in bar)');

fprintf('Results');
disp(' ');

%stage 1 turbine inlet
s1 = XSteam('s_pT',p1,t1);
h1 = XSteam('h_pT',p1,t1);

% stage 2 condensor inlet
sf2 = XSteam('sL_p',p2);
sg2 = XSteam('sV_p',p2);
hf2 = XSteam('hL_p',p2);
hg2 = XSteam('hV_p',p2);

s2=s1;
t2 = XSteam('Tsat_p',p2);

%split the formula
x2 = ((s2-sf2)/(sg2-sf2));
h2 = hf2 + (x2*(hg2-hf2));

%stage 3 pump inlet
p3=p2;
t3 = XSteam('Tsat_p',p2);
s3 = XSteam('sL_p',p3);
h3 = XSteam('hL_p',p3);

cp4 = XSteam('cpL_p',p3);
cv4 = XSteam('cvL_p',p3);
gamma = (cp4/cv4);
g = (gamma-1)/gamma;

%stage 4 boiler inlet
p4=p1;
s4=s3;
h4 = XSteam('h_ps',p4,s4);
t4 = t3*((p4/p3)^g);

%turbine work
wt = h1-h2;
%pump work
wp = h4-h3;
% net work
w_net = wt-wp;

%heat in and out of the system
q_in = h1-h4;
q_out = h2-h3;

n=w_net/q_in;
np = n*100;

SSC=3600/w_net;
tf_sat = XSteam('Tsat_p',p1);
sf_sat = XSteam('sL_p',p1);
hf_sat = XSteam('hL_p',p1);
sg_sat = XSteam('sV_p',p1);
hg_sat = XSteam('hV_p',p1);

%plot figure 1
figure(1);
hold on
temp= linspace(0,400,1000);
for i=1:length(temp)
    aa(i) =XSteam('sL_T',temp(i));
end
for j=1:length(temp)
    bb(j) = XSteam('sL_T',temp(j));
end
plot(aa,temp,'Color','b')
plot(bb,temp,'Color','b')

temp=XSteam('Tsat_p',p1);
t0=linspace(temp,t1,300);
for i=1:length(t0)
    c11(i)=XSteam('s_pt',p1,t0(i));
end

t5 =XSteam('Tsat_p',p4);
s5 =XSteam('sL_T',t5);

plot([s1 s2],[t1 t2],"Color",'r');
plot([s2 s3],[t2 t3],"Color",'r');
plot([s3 s4],[t3 t4],"Color",'r');
plot([s4 s5],[t4 t5],"Color",'r');
plot([sf_sat sg_sat],[tf_sat tf_sat],"Color",'r');
plot(c11,t0,'Color','r');

title('T-S diagram')
xlabel('entropy(KJ/kgK)')
ylabel('Temperature(K)')

text(s1,t1,'1');
text(s2,t2,'2');
text(s3,t3,'3');
text(s4,t4,'4');
text(s5,t5,'5');
text(sg_sat,tf_sat,'1s');

%plot figure 2
figure(2);
hold on

temp1=linspace(0,400,1000);
for i=1:length(temp1)
    cc(i)=XSteam('hV_T',temp1(i));
end
for j=1:length(temp1);
    dd(j)=XSteam('hV_T',temp1(j));
end

plot(aa,cc,"Color",'b');
plot(bb,dd,"Color",'b');

plot([s1 s2],[h1 h2],"Color",'r');
plot([s2 s3],[h2 h3],"Color",'r');
plot([s3 s4],[h3 h4],"Color",'r');
plot([s4 s1],[h4 h1],"Color",'r');

title('H-S Diagram')
xlabel('entropy(KJ/kgK)')
ylabel('enthalpy(KJ/Kg)')

text(s1,h1,'1');
text(s2,h2,'2');
text(s3,h3,'3');
text(s4,h4,'4');

%printing the output
fprintf('at stage point 1')
fprintf('P1 is : %f',p1);
fprintf('T1 is : %f deg celcius',t1);
fprintf('H1 is : %f KJ/kg',h1);
fprintf('S1 is : %f KJ/KgK',s1);

disp(' ');
fprintf('at stage point 2')
fprintf('P2 is : %f',p2);
fprintf('T2 is : %f deg celcius',t2);
fprintf('H2 is : %f KJ/kg',h2);
fprintf('S2 is : %f KJ/KgK',s2);
fprintf('X is : %f rn',x2);

disp(' ');
fprintf('at stage point 3')
fprintf('P3 is : %f',p3);
fprintf('T3 is : %f deg celcius',t3);
fprintf('H3 is : %f KJ/kg',h3);
fprintf('S3 is : %f KJ/KgK',s3);

disp(' ');
fprintf('at stage point 4')
fprintf('P4 is : %f',p4);
fprintf('T4 is : %f deg celcius',t4);
fprintf('H4 is : %f KJ/kg',h4);
fprintf('S4 is : %f KJ/KgK',s4);

disp(' ');
fprintf('wt : %f KJ/Kg',wt);
fprintf('wp : %f KJ/Kg',wp);
fprintf('Wnet is %f KJ/Kg ',w_net);
fprintf('Nterm : %f %%n ',n*100)
fprintf('SSC : %f KJ/ KWh',SSC);