KD PDF
KD PDF
KD PDF
(Rankine Cycle)
Previous
Cycle (thermodynamic cycle) is a sequence of processes that begins
and ends at the same state.
A power cycle (or a heat engine) is an installation with a working fluid,
that experiments a thermodynamic cycle, and it produces work as a
consequence of the heat exchanged by the working fluid.
Working fluid is the substance that makes the cycle passing through
the different devices of the installation.
(Steam=vapor water)
Contents
Introduction
Rankine cycle
Principal Work and Heat Transfers. Thermal efficiency
How to improve thermal efficiency of Rankine cycle?
Modifications of Rankine cycle:
Reheat
Regenerative steam power cycle
Introduction
Steam power plants are the most conventional source
of electric power. These power plants are variations of
a thermodynamic cycle in which water is the working
fluid. In this kind of cycle the water is in liquid phase
in a part of the cycle and it is in vapor phase in another
one.
On the other hand, steam turbine are the responsible
of more than 70% of power electricity generated in the
world. They are present in many types of power plants,
such as nuclear plants, conventional thermal plants,
combined cycles, biomass plants and solar power
plants.
Introduction
Introduction
In this figure we have an schematic draw of a very simplified steam
power plant. (Source: Moran Shapiro)
Rankine cycle
T-s diagram for Rankine cycle
Diagrama T-s pel cicle de Rankine ideal
T
1
Qcaldera
4
3
2
Qcondensador
s3 = s4
s1 = s2
= h1 h 2
cond
= h
3) Pump
wbomba = h4 h3
4) Boiler
q boiler
= h1 h 4
Thermal efficiency
Thermal efficiency:
wturb w pump
wnet
=
=
qboyler
qboyler
Thermal efficiency can be also written as an alternative form :
T c ( s1 s 4 )
Tc
=1
= 1
T m ( s1 s 4 )
Tm
Tm (average heat-addition temperature)
Tc (temperature at which condenser operates)
1
Qcaldera
Tm
4
3
2
Qcondensador
s3 = s4
s1 = s2
T
1
Tm
4'
3'
s 3 ' = s 4'
4
3
2
2'
s3 = s4
s 1 = s 2 = s2'
Problems
Decreases the quality at the exit of the turbine (X2<x2). The presence of
more than 10% liquid water (quality 90%) can cause erosion of the blades.
Air leaks in the condenser (if the condenser pressure decreases a lot)
T
1
1'
Tm'
Tm
4'
4
3
s 3 = s 4 = s 4'
2'
s 1' = s 2'
s1 = s2
Problems
Decreases the quality at the exit of the turbine (Quality 90%)
Nowadays the maximum boiler pressure achieved is 30 MPa (>22,1 MPa)
supercritical cycle
1'
1
Tm'
Tm
4
s 3 = s 4 = s 4'
2'
s 1 = s 2 s 1' = s 2'
Exercises
1.- A Rankine cycle has steam entering the turbine at 30 bar and 400C,
and a condenser pressure of 0,05 bar. Turbine and pump are
considered isentropic (adiabatic and reversible process). Determine:
a) the net work of the cycle; b) the heat added to the cycle; c)
thermal efficiency of the cycle. (source: M. J. Moran and H.N.
Shapiro, Fundamentals of Engineering Thermodynamics, 2nd
Edition (1992) John Wiley & Sons, Inc. (USA)).
2.- Solve Prob.1 adding an isentropic turbine efficiency of 90%.
3.- Make a MATLAB program to study the improvement of the thermal
efficiency of vapor power plants with the boiler pressure
4.- Make a MATLAB program to study the improvement of the thermal
efficiency of vapor power plants with the turbine inlet temperature.
5.- Make a MATLAB program to study the improvement of the thermal
efficiency of vapor power plants with the decrease of the
condenser pressure.
clear;
% Prob. 1 Ideal Rankine Cycle
P1=30; % pressure bar
T1=400; % temperature Celsius
P2=0.05; % bar
% Previously we calculate thermodynamical properties in diff.
states
h1=XSteam('h_pT',P1,T1)
s1=XSteam('s_pT',P1,T1)
s2=s1
% Detemine in which state is water
sl=XSteam('sL_p',P2)
sg=XSteam('sV_p',P2)
if s2<=sg
x2=(s2-sl)/(sg-sl);
disp ('liquid-vapor mixture:')
disp ('x2=')
disp (x2)
h2=XSteam('hL_p',P2)+x2*(XSteam('hV_p',P2)XSteam('hL_p',P2))
else
disp ('superheated vapor')
h2=XSteam('h_ps',P2,s2)
end
Prob. 1 cont.
h3=XSteam('hL_p',P2)
s3=XSteam('sL_p',P2)
s4=s3;
h4=XSteam('h_ps',P1,s4)
% a) the net work of the cycle
W_turb=h1-h2
W_pump=h4-h3
w_net=W_turb-W_pump
% b) The heat added to the cycle
Q_add=h1-h4
%c) Thermal efficiency
Eff=w_net/Q_add
Prob. 1 cont.
fprintf('Solutions: \n')
fprintf('a)net work: %6.3fkJ/kg \n',w_net)
fprintf('b)Heat added: %6.3fkJ/kg \n',Q_add)
fprintf('c)Thermal efficiency: %6.3f \n',Eff)
if s2s<=sg
x2s=(s2s-sl)/(sg-sl);
disp ('liquid-vapor mixture:')
disp ('x2s=')
disp (x2s)
h2s=XSteam('h_px',P2,x2s)
else
disp ('superheated vapor')
h2s=XSteam('h_ps',P2,s2s)
end
W_s=h1-h2s
% IST=W_12/W_s
fun=@(x)IST-(h1-x)/(h1-h2s)
h2=fsolve(fun,1000)
%Detemine in which state is water in 2
hl=XSteam('hL_p',P2)
hg=XSteam('hV_p',P2)
if h2<=hg
x2=(h2-hl)/(hg-hl);
disp ('liquid-vapor mixture:')
disp ('x2=')
disp (x2)
s2=XSteam('sL_p',P2)+x2*(XSteam('sV_p',P2)-XSteam('sL_p',P2))
else
disp ('superheated vapor')
s2=XSteam('s_ph',P2,h2)
end
h3=XSteam('hL_p',P2)
s3=XSteam('sL_p',P2)
s4=s3;
h4=XSteam('h_ps',P1,s4)
Prob.3
clear
hold off;
vp=linspace(80,130,20)
for i=1:length(vp)
vr(i)=rendiment_rankine2(vp(i))
end
plot(vp,vr,'o-');
title Thermal efficiency vs. boiler pressure (bar)'
xlabel('P1 (bar)')
ylabel('\eta')
Plot Prob. 3
Thermal efficiency vs. boiler pressure (bar)
0.425
0.42
0.415
0.41
0.405
0.4
80
85
90
95
100
105
110
P1 (bar)
115
120
125
130
Tmr i
2
Tm
6
5
s5 = s6
Source: Moran-Shapiro
W and Qabs
s1 = s2
s3 = s4 s
Excercise
Water is the working fluid in a vapor power cycle with reheat.
Steam at 80 bar, 480C enters the first-stage turbine and
expands to 24 bar. Then, it is reheated to 480C before
entering the second turbine stage, where it expands to the
condenser pressure of 0.08 bar. The mass flow rate of steam is
73 kg/s. Each turbine stage operates with an isentropic
efficiency of 88%. The pump operates with an isentropic
efficiency of 92%. Determine for the cycle:
a) The net power developed, in kW
b) The thermal efficiency
c) the rate of heat transfer to cooling water passing to the
condenser, in kW
clear;
% Prob. 8 Reheat
P1=80; % bar
T1=480; %C
P2=24; % bar
P3=P2;
T3=480; % C
P4=0.08; % bar
m=73; % kg/s
RIT=0.88; % isentropic efficiency for each turbine stage
RITP=0.92; % isentropic pump efficiency
% Determination of thermodynamic properties of each state of cycle
h1=XSteam('h_pT',P1,T1);
s1=XSteam('s_pT',P1,T1);
s2s=s1;
% Detemine in which state is water at 2s state
sl=XSteam('sL_p',P2)
sg=XSteam('sV_p',P2)
if s2s<=sg
x2s=(s2s-sl)/(sg-sl)
disp ('liquid-vapor mixture:')
disp ('x2s=')
disp (x2s)
h2s=XSteam('h_px',P2,x2s)
disp ('superheated vapor')
h2s=XSteam('h_ps',P2,s2s);
end
W_12s=h1-h2s; % in kJ/kg
W_12=RIT*W_12s % actual work first-stage turbine
h2=h1-W_12;
s2=XSteam('s_ph',P2,h2);
h3=XSteam('h_pT',P3,T3);
s3=XSteam('s_pT',P3,T3);
s4s=s3
% Detemine in which state is water at 4s state
sl=XSteam('sL_p',P4)
sg=XSteam('sV_p',P4)
if s4s<=sg
x4s=(s4s-sl)/(sg-sl);
disp ('liquid-vapor mixture:')
disp ('x4s=')
disp (x4s)
h4s=XSteam('h_px',P4,x4s);
else
disp ('superheated vapor')
h4s=XSteam('h_ps',P4,s4s)
end
W_34s=h3-h4s;
W_34=RIT*W_34s % actual work second-stage turbine
h4=h3-W_34;
s4=XSteam('s_ph',P4,h4);
h5=XSteam('hL_p',P4);
s5=XSteam('sL_p',P4);
s6s=s5;
h6s=XSteam('h_ps',P1,s6s);
W_56s=h6s-h5
W_56=W_56s/RITP
h6=h5+W_56;
s6=XSteam('s_ph',P1,h6);
% a)
Power_turb=m*W_12+m*W_34
Power_pump=m*W_56
Power_net=Power_turb-Power_pump
% b)
Q_add=m*(h1-h6)+m*(h3-h2)
ThEf=Power_net/Q_add
% c)
Q_cond=m*(h4-h5) % in absolute value
fprintf('Solutions: \n')
fprintf('a)net power: %6.3fkW \n',Power_net)
fprintf('b)Thermal efficiency: %6.8f \n',ThEf)
fprintf('c)The rate of heat transfer to cooling
water (condenser): %6.3fkW \n',Q_cond)
Source: Moran-Shapiro
2
6
a
5
4
s 4 = s 5 s6 = s 7
s1 = s2 = s3
In a steam power plant there are several feedwater heater (open and closer). The number of
fedwater heaters used is based on economic considerations, since incremental increases in
thermal efficiency achieved with each additional heater must justify the added capital costs.
The maximum number of heaters used is about 10. A steam power plant used to have one
reheat
(maximum
2).
Exercise
A power plant operates on a regenerative vapor power cycle
with one open feedwater heater. Steam enters the first turbine
stage at 120bar and 520C and expands to 10 bar, where some of
the steam is extracted and diverted to the open feedwater
operating at 10 bar. The remaining steam expands through the
second turbine stage to the condenser pressure of 0.06 bar.
Saturated liquid exists the open feedwater heater at 10 bar. For
isentropic processes in the turbines and pumps, determine:
a) The fraction of mass flow rate extracted from 1st turbine
stage (sol: y=0.2336)
b) The net power, considering a mass flow rate if 106 kg/h
(sol: 3325.1 MW)
c) The thermal efficiency (sol: 0.4554)
Source: Moran-Shapiro
clear;
P1=120
T1=520
P2=10
P3=0.06
m1=10^6; %kg/h
m=10^6/3600; %kg/s
% Calculate state 2
s1=XSteam('s_pT',P1,T1);
h1=XSteam('h_pT',P1,T1);
s2=s1;
sg2=XSteam('sV_p',P2);
sl2=XSteam('sL_p',P2);
if s2<=sg2
x2=(s2-sl2)/(sg2-sl2);
h2=XSteam('hL_p',P2)+x2*(XSteam('hV_p',P2)-XSteam('hL_p',P2));
else h2=XSteam('h_ps',P2,s2);
end
% Calculate state 3
s3=s2;
sg3=XSteam('sV_p',P3);
sl3=XSteam('sL_p',P3);
if s3<=sg3
x3=(s3-sl3)/(sg3-sl3);
h3=XSteam('hL_p',P3)+x3*(XSteam('hV_p',P3)-XSteam('hL_p',P3));
else h3=XSteam('h_ps',P3,s3);
end
h4=XSteam('hL_p',P3);
s4=XSteam('sL_p',P3);
s5=s4
h5=XSteam('h_ps',P2,s5);
h6=XSteam('hL_p',P2);
s6=XSteam('sL_p',P2);
s7=s6;
h7=XSteam('h_ps',P1,s7);
% a) fraction of the total flow extracted;
fun=@(y)h6-h5*(1-y)-h2*y
y=fsolve(fun,0.1)
% b) net power developed;
W_turb=m*(h1-h2)+m*(1-y)*(h2-h3);
W_pump=m*(1-y)*(h5-h4)+m*(h7-h6);
W_net=W_turb-W_pump
%c) Thermal efficiency;
eff=W_net/(m*(h1-h7))
fprintf('Solutions: \n')
fprintf('a)extracted fraction y: %6.3f \n',y)
fprintf('b)net power: %6.3fkW \n',W_net)
fprintf('c)Thermal efficiency: %6.3f \n',eff)
Exercise
From the last exercise, make a graph of the thermal efficiency versus extraction
pressure, P2.
0.455
0.45
0.445
0.44
0.435
0.43
0.425
10
20
30
40
50
60
Pext (bar
70
80
90
100
function [eff]=Regenerative(P2)
% P2 is the extraction pressure
This is a function.
P1=120
This program is called Regen_func
T1=520
P3=0.06
m1=10^6; %kg/h
m=10^6/3600; %kg/s
% Calculate state 2
s1=XSteam('s_pT',P1,T1);
h1=XSteam('h_pT',P1,T1);
s2=s1;
sg2=XSteam('sV_p',P2);
sl2=XSteam('sL_p',P2);
if s2<=sg2
x2=(s2-sl2)/(sg2-sl2);
h2=XSteam('hL_p',P2)+x2*(XSteam('hV_p',P2)-XSteam('hL_p',P2));
else h2=XSteam('h_ps',P2,s2);
end
% Calculate state 3
s3=s2;
sg3=XSteam('sV_p',P3);
sl3=XSteam('sL_p',P3);
if s3<=sg3
x3=(s3-sl3)/(sg3-sl3);
h3=XSteam('hL_p',P3)+x3*(XSteam('hV_p',P3)-XSteam('hL_p',P3));
else h3=XSteam('h_ps',P3,s3);
end
h4=XSteam('hL_p',P3);
s4=XSteam('sL_p',P3);
s5=s4
h5=XSteam('h_ps',P2,s5);
h6=XSteam('hL_p',P2);
s6=XSteam('sL_p',P2);
s7=s6;
h7=XSteam('h_ps',P1,s7);
% a) fraction of the total flow extracted;
fun=@(y)h6-h5*(1-y)-h2*y
y=fsolve(fun,0.1)
% b) net power developped;
W_turb=m*(h1-h2)+m*(1-y)*(h2-h3);
W_pump=m*(1-y)*(h5-h4)+m*(h7-h6);
W_net=W_turb-W_pump
%c) Thermal efficiency;
eff=W_net/(m*(h1-h7))
end
Bibliography
M. J. Moran and H.N. Shapiro, Fundamentals
of Engineering Thermodynamics, 2nd Edition
(1992) John Wiley & Sons, Inc. (USA)
J.B. Jones and R.E. Dugan, Engineering
Thermodynamics, (1996) Prentice Hall Inc.