1

10. MODELAREA TURBINEI CU ABUR I A

GRUPULUI AFERENT N REGIM STAIONAR
10.1. Calculul parametrilor ciclului pentru determinarea unui regim
de referin
Schema termic dup care se va realiza calculul este prezentat n figura
10.1, de mai jos, iar datele constructive i parametrii principali sunt dup datele
tehnice cunoscute ale grupului echipat cu turbine de 330 MW, i anume:
puterea la borne P
B
= 330 MW
presiunea aburului viu p
0
= 192 bar
temperatura aburului viu t
0
= 540
0
C
temperatura aburului supranclzit t
0
= 540
0
C

Pentru calcul s-au considerat grupuri de 330 MW de producie
romneasc (toate agregatele de producie romneasc).
Aceast schem nu ia n consideraie pierderile de abur la labirini i tijele
ventilului, apa de injecie pentru reglarea temperaturii aburului viu i
supranclzit intermediar, alte recuperri de cldur.
Fa de aceasta, calculele din lucrare s-au efectuat n urmtoarele ipoteze:
pierderi de presiune pe conducte i n prenclzitoare constante;
diferenele finale de temperatur pe prenclzitoare constante;
re-pomparea condensatului la PJP 2;
randamentele pompelor constante;
se recupereaz n circuit 85% din pierderile la generator i mecanice;
cota procentual de recuperri de abur a fost raportat la debitul de
abur viu i considerat la presiunea prizei respective.
2

Fig. 10.1. Schema termic de calcul
Stabilirea limitelor de variaie a presiunii de condensare
S-au considerat urmtoarele valori:
- minimul de iarn pentru temperatura turnului: t
turn
= 1518
0
C ;
- maximul de var pentru temperatura turnului: t
turn
= 3033
0
C ;
- minimul de iarn pentru temperatura rului: t
ru
= 45
0
C;
- maximul de var pentru temperatura rului: t
ru
= 2831
0
C .
innd cont de cota sub care se gsete, n sistemul de rcire, apa din ru i
cea din circuitul turnului de rcire se calculeaz o medie de temperatur a apei
de rcire:
t = c t
ru
+ (1 - c) t
turn

unde: c - cota de ru , c = 50% .
Pentru minimul de iarn: t = 0,55 + (1 - 0,5)18 t = 11,5
0
C
3
Deci: t
c
= t + t + t
unde: t = 812
0
C; t = 36
0
C .
Atunci:
t
c1
= 11,5 + 8 + 3 t
c1
= 22,5
0
C
t
c2
= 11,5 + 12 + 6 t
c2
= 29,5
0
C .
Prin urmare vom avea presiunile:
p
c1
= 0,0272 bar
p
c2
= 0,0412 bar .
Alegem o izobar exact: p
c min
= 0,03 bar , pentru care avem temperatura
t
c min
= 24,10
0
C .
Se presupune t = 4
0
C i atunci vom avea :
t = t
c min
- t - t t = 24,103 - 11,5 - 3 t = 9,603
0
C .

10.2. Stabilirea unui regim de referin
Avnd la baz caracteristicile schemei termice a grupului de 330 MW,
prezentat n ,,Instruciuni de exploatare pentru turbina FIC - 330 MW, s-au
stabilit debitele i presiunile prin turbin, pentru schema i notaiile din figura
10.1.
Din bibliografia ICEMENERG , pentru o putere de 330 MW , cu notaiile
din tabelul din figura 10.1, rezult urmtoarele valori :
admisie n CIP : p
3
= 182,4 bar ;
priza 7 (ieirea din CIP): p
5
=48,7 bar ;
admisie n CMP : p
8
=44,705 bar ;
priza 6 : p
9
=30,701bar ;
priza 5 : p
10
=18,497 bar ;
priza 4 : p
11
= 8,6835 bar ;
priza 3 (ieire CMP) : p
12
= 4,6204 bar ;
4
priza 2 : p
14
= 1,2684 bar ;
priza 1 : p
15
= 0,31 bar ;
ieire CJP : p
16
= p
17
= 0,0373 bar .
Pentru luarea n considerare a pierderilor de presiune pe conducte i a
pierderilor la etanri , care n datele iniiale erau neglijate, s-a fcut astfel :
- cunoscnd presiunea la intrarea n turbin, respectiv punctul 3 din
schem p
3
= 182,4 bar , am estimat presiunile regimului de referin , nainte de
intrarea n ventilele de nchidere rapid VIR i de reglare ale corpului CMP;
astfel c aici am estimat o pierdere de presiune p 2 % .
- p
2
=p
3
/0,98 = 182,4/0.98 = 186,11 p
2
= 186,11 bar
De asemenea la ieirea din cazan s-a estimat o pierdere de cca.: p 34 %.
p
1
= p
2
/0,969 = 182,4/0,969 =192 p
1
= 192 bar .
Alte predeterminri ale presiunilor necunoscute nc:
Cunoscnd presiunea la ieirea din CIP, p
5
=48,7 bar, (presiune
corespunztoare alimentrii prenclzitorului 7 din schem), am estimat
presiunea la intrare n supranclzitorul intermediar de :
p
6
= p
5
0,971
p
6
= 48,7 0,971 = 47,302
p
6
= 47,302 bar .
n supranclzitorul intermediar considerat o pierdere de presiune p=
0,968 bar , astfel nct : p
7
= p
6
0,968
p
7

= 47,302 0,968 = 45,804
p
7
= 45,804 bar .
Presiunea treptei de reglare , respectiv p
4
, am determinat-o din condiia
unei cderi de entalpie pe treapt de aproximativ 60 kJ, concomitent cu o
temperatur la ieire din treapt de cca. 500
0
C, astfel nct p
4
= 151,88 bar, ceea
ce corespunde i cu presiunile n treapta de reglare prezentate de ICSITEE.
5
Estimez pierderile pe conductele de abur supranclzit intermediar CMP
inclusiv n ventilele de nchidere rapid MP i ventilele de reglare (VR):
p = p
8
/p
7
=44,705 / 45,804
p = 0,976
Pentru corpul de medie presiune se cunosc :
Are patru prize i anume: priza care alimenteaz prenclzitorul 3, respectiv
ieire din CMP la p
12
= 4,6204 bar . Celelalte prize sunt: cea care alimenteaz
prenclzitorul 4 : p
11
= 8,6835 bar , priza care alimenteaz prenclzitorul 5 : p
10

= 18,497 bar , priza care alimenteaz prenclzitorul 6 : p
9
=30,701 bar (din
datele ICEMENERG).
Deoarece ntre corpurile CMP i CJP am considerat o pierdere de presiune
de cca. p = 2% p
13
= p
12
0,98 = 4,6204 0,98
p
13
= 4,5304 bar .
Presiunea de condensare , respectiv ieirea din CJP a avut o valoare de p
c

= 0,0373 bar , astfel c :
p
16
= p
17
= p
c
= 0,0373 bar (din datele ICEMENERG).
Prizele n CJP , corp n dublu flux i dou corpuri n schem, sunt deci
prizele care alimenteaz prenclzitoarele 1 i 2, avnd valoarea :
- priza care alimenteaz prenclzitorul 1 :
p
15
= 0,31 bar
- priza care alimenteaz prenclzitorul 2 :
p
14
= 1,2684 bar (din datele ICEMENERG).
Deoarece regimul de referin l-am considerat cu TPA , alimentarea acesteia
se face de la ieirea din corpul CMP. Am considerat p 3% pierderi de
presiune pe conducta turbinei principale (turbina TPA) , astfel c la intrarea n
turbopomp:
p
18
= p
12
0,97 = 4.6204 0,97
p
18
= 4,4818 bar
6
Presiunea la ieirea din turbopomp am considerat-o aceeai cu presiunea
de condensare, i pentru regimul de referin: p
19
= 0,0373 bar .
Ca urmare a acestor estimri i de asemeni ca urmare a datelor luate
conform ICEMENERG , presiunile de referin sunt :
p
1
= 192 bar ; p
7
= 45,804 bar ; p
13
= 4,5304 bar ;
p
2
= 186,11 bar ; p
8
= 44,705 bar ; p
14
= 1,2684 bar ;
p
3
= 182,4 bar ; p
9
= 30,701 bar ; p
15
= 0,31 bar ;
p
4
= 151,88 bar ; p
10
= 18,497 bar ; p
16
= 0,0373 bar ;
p
5
= 48,7 bar ; p
11
= 8,6835 bar ; p
17
= 0,0373 bar ;
p
6
= 47,302 bar ; p
12
= 4,6204 bar ; p
18
= 4,4818 bar ;
p
19
= 0,0373 bar .
Pentru estimarea debitelor prin turbin am cunoscut tot din datele
respective, valorile debitelor de abur la prize, inclusiv TPA, date mai jos :
debitul la priza 1 : d
51
= 11,7 kg/s ;
debitul la priza 2 : d
52
= 13,06 kg/s ;
debitul la priza 3 : d
53
= 16 kg/s ;
debitul la priza 4 : d
54
= 10,06 kg/s ;
debitul la priza 5 : d
55
= 12,083 kg/s ;
debitul la priza 6 : d
56
= 14 kg/s ;
debitul la priza 7 : d
57
= 17,42 kg/s ;
debitul la TPA : d
18
= 15,17 kg/s .
De asemenea , am considerat debitul de abur intrat n turbin (din datele
ICEMENERG), respectiv la admisie n CIP :
d
1
= 273,33 kg/s (=
[ ]
[ ]
[ ] s kg
s
h t
/ 33 , 273
3600
10 / 984
3
=

) .

Necesitatea recalculrii unui regim de referin
7
Calculele termodinamice ale proprietilor ap - abur , nu sunt indicate ,
astfel c este strict necesar ca toate evalurile s fie efectuate pentru comparaie
cu acelai set de proprieti ap - abur .
Un alt motiv care a stat la baza necesitii recalculrii schemei l constituie
valoarea entalpiei la prize care nu este cunoscut , i care presupune
determinarea randamentelor interne prin turbin cu relaii care s fie valabile n
toate regimurile .
Datele decalcul ICEMENERG nu iau n considerare pierderile de abur la
labirini i n tijele ventilului , i nici apa de injecie pentru reglarea temperaturii
aburului viu i a aburului supranclzit intermediar , ceea ce constituie unul din
motivele recalculrii unui regim de referin pentru a putea estima celelalte
regimuri de funcionare ce vor fi analizate.
Estimarea pierderilor la labirini i la ventilele de nchidere rapid (VIR) a
fost considerat ca o cot procentual din debitul de abur viu , ceea ce pe de alt
parte respect i formula lui MARTIN (din R. Spencer i [Link]) :
d =
1
1
v
p
C ,
unde: C - o constant; p
1
,v
1
: sunt presiunea, respectiv volumul
specific la intrarea n pachetul de labirini; respect aceasta deoarece
conform relaiei simplificate a lui Stodola:

v
p
v
p
d
d
nom
nom
nom
= ,
unde: d - debitul scpat prin labirini la un alt regim; d
nom
- debitul
scpat prin labirini la regimul de referin .
Deci modificarea raportului p/v aproximeaz debitul (conform lui Rizhkin).
Deci rezult n final (din Rizhkin), pierderile la labirini i ventile de nchidere
rapid (adoptate la pachetul de 330 MW).
b
1
: scpri la tija ventilului b
1
= 0,0028 .
b
2
: scpri pe partea de nalt presiune a CIP - ului b
2
= 0,0072 .
8
b
3
: scpri la labirinii prii de ieire din CIP b
3
= 0,0053 .
b
4
: scpri la labirinii de etanare dinspre partea de nalt
presiune a CMP - ului b
4
= 0,0016 .
b
5
: scpri la labirinii de ieire ai CMP - ului b
5
= 0,0003 .
Aburul care intr pentru etanarea de atmosfer (prin camerele de ventilare
de la sfrit), se consider o recuperare n condensator cu o cot: b
6
= =0,00053 .
Pentru reglarea temperaturii aburului am considerat urmtoarele :
pentru reglarea temperaturii aburului viu se injecteaz ap de
alimentare dup prenclzire n valoare de: a
inj v
= 0,1% = 0,001;
pentru reglarea temperaturii aburului supranclzit intermediar am
considerat o injecie: a
inj SI
= 0,08% = 0,0008 .
innd seama de cele de mai sus , au rezultat urmtoarele debite prin
turbin (debite care constituie punctul de pornire pentru calculul regimului de
referin i anume n mod special al randamentelor: = f(debite prin turbin):
d
1
= 274,1 kg/s ; d
7
= 252,72 kg/s ; d
13
= 184,89 kg/s ;
d
2
= 274,1 kg/s ; d
8

= 252,72 kg/s ; d
14
= 184,89 kg/s ;
d
3
= 273,33 kg/s ; d
9
= 252,28 kg/s ; d
15
= 171,83 kg/s ;
d
4
= 271,37 kg/s ; d
10
= 238,28 kg/s ; d
16
= 160,13 kg/s ;
d
5
= 269,92 kg/s ; d
11
= 226,2 kg/s ; d
17
= 160,28 kg/s ;
d
6
= 252,5 kg/s ; d
12
= 216,14 kg/s .
Alte date necesare calculului au fost:
temperatura aburului viu t
0
= 540
0
C i
temperatura aburului supranclzit t
SI
= 540
0
C .
10.3. Calculul destinderii aburului n turbin
Punctul 1 : p
1
= 192 bar ; t
1
= 540
0
C
( ) ( ) kJ/kg 1 , 3376 1 , 3376 540 ; 192 ;
1
0
1 1 1
= = = = h C bar h t p h h
( ) ( ) K kJ/kg 3393 , 6 3393 , 6 / 1 , 3376 ; 192 ;
1 1 1 1
= = = = s kg kJ bar s h p s s
9
( ) ( ) /kg m 0169472 , 0 0169472 , 0 540 ; 192 ;
3
1
0
1 1 1
= = = = v C bar v t p v v
titlul aburului: x = 1 .
Pe conducta cazan turbin am estimat o pierdere de entalpie:
h
caz-tb
3 kJ/kg.

Punctul 2 : h
2
= h
1
- h
caz-tb
= 3376,1 - 3 h
2
= 3373,1 kJ/kg ;
p
2
= 186,11 bar ;
( ) ( ) kJ/kg 67 , 536 67 , 536 / 1 , 3373 ; 11 , 186 ;
2 2 2 2
= = = = t kg kJ bar t h p t t
( ) K kJ/kg 3476 , 6 3476 , 6 ;
2 2 2 2
= = = s h p s s
( ) = = 017454 , 0 ;
2 2 2
h p v v v
2
=0,017454 m
3
/kg
titlul aburului : x = 1.

Punctul 3 : p
3
=182,4 bar ; h
3
=h
2
= 3373,1 kJ/kg (are loc o laminare n
ventile);
( ) ( ) C t kg kJ bar t h p t t
0
3 3 3 3
22 , 535 22 , 535 / 1 , 3373 ; 4 , 182 ; = = = =
( ) K kJ/kg 3553 , 6 3553 , 6 ;
3 3 3 3
= = = s h p s s

( ) ( )
/kg m 017813 , 0
017813 , 0 / 1 , 3373 ; 4 , 182 ;
3
3
3 3 3
=
= = =
v
kg kJ bar v h p v v

Cderea teoretic pe treapta de reglare :
( )
3 4 3 4 3
; s p h h h h h
t ttr
= = ,
unde: p
4
=151,88 bar i s
4t
= s
3
= 6,3553 kJ/kgK (proces izentropic).
( ) kg kJ h K kg kJ bar h h
ttr ttr
/ 17 , 56 17 , 56 / 3553 , 6 ; 88 , 151 1 , 3373 = = =
(cum s-a aproximat i la calculul presiunii p
4
).
Randamentul treptei de reglare (TR treapt de reglare cu o singur
coroan i cu un diametru d
TR
=865 mm constructiv); deci
TR
l-am calculat cu
relaia lui Schegliiev:

|
|

\
|

1
5
1
4
/ TR
10
10 2
83 , 0
v
p
G
K
c u
,
10
unde: p
1
, v
1
sunt presiunea respectiv volumul specific la intrarea n
TR;
G,[kg/s]; p,[bar]; v,[m
3
/kg] ).
Unde K
u/ c
=
( )
( )
(

optim
fictiv
c u
c u
f
1
1
/
/
, (raportul vitezelor dat n diagrama 6.9); n
care :
( )
( )

= =
1 2
cos
/
1
1
x
opt optim
K
x c u , unde s-a considerat :

1
= 13
0
14
0
; ( cos
1max
=1)
: grad de reaciune al paletelor: = 0,05 ;
K
u
= 0,95
Atunci avem : ( )
( )
. 5 , 0
05 , 0 1 2
95 , 0 1
/
1
=

= =
opt optim
x c u
Iar
( ) ( ) 405386 , 0
17 , 56 2000
865 , 0 50
/
2000
50
/
1 1
=

fictiv
tTR
TR
fictiv
c u
h
d
c u
Diagrama 6.9, [7], am transformat-o ntr-o funcie parabolic de gradul trei ,
astfel c: K
u/c
= C
1
x
3
+ C
2
x
2
+ C
3
x + C
4
,
unde : C = [0,3604 ; -2,119 ; 3,1523 ; -0,03936] .
Notez cu
( )
( )
optim
fictiv
c u
c u
x
1
1
/
/
= , necunoscuta ecuaiei de gradul trei.
81077 , 0 81077 , 0
5 . 0
405386 , 0
= = = x x .
Atunci vom avea c :
K
u/c
=0,36040,81077
3
-2,1190,81077
2
+3,15230,81077 - 0,03936
K
u/c
= 0,94175.
Facem nlocuirile n relaia randamentului :
11

7596 , 0
7596 , 0
017813 , 0
10 4 , 182
333 , 273
10 2
83 , 0 94175 , 0
5
4
=
=
|
|

\
|

TR
TR

Entalpia real la ieirea din treapt va fi:
h
4
=h
3
-
TR
h
tTR
h
4
= 3373,1 - 0,759656,17 = 3330,5
h
4
= 3330,5 kJ/kg.
Astfel nct:

Punctul 4:t
4
=t(p
4
;h
4
)=t(151,88bar;3330,5kJ/kg) t
4
=507,55
0
C;
s
4
=s(p
4
;h
4
)=s(151,88bar;3330,5kJ/kg)s
4
= 6,3726kJ/kgK
v
4
=v(p
4
;h
4
)=v(151,88bar;3330,5kJ/kg) v
4
=0,0208 m
3
/kg.

Estimarea randamentului pe treptele de presiune CIP :
estimarea randamentului a fost fcut cu relaia :
|

\
|
|

\
|

+
|
|
|
|

\
|

\
|
=
tTP
vs tTP
med med
CIP
H
H
v D
n

1
20000
600
1
4 , 0
50
1
925 , 0
2
;
unde: n - frecvena (= 50Hz)

vs
= 5 , 3 5 , 3 sin
1
1
2
=
vs
z

(z - numrul de trepte de presiune).
Atunci : 011 , 0 011 , 0
13 , 317
5 , 3
= = =
tTP
vs
tYP
vs
H H


Diametrul mediu D
med
, n aceast situaie este debitul care intr n corp ,
deoarece nu sunt prize n CIP :
Volumul specific mediu v
med
: [ ] kg m v v v v v
e i med
/
3
5 4
= = .
La nceput ntr-o prim iteraie v
e
l iau volumul teoretic de la ieire:
v
e
= v
5t
= v(p
5
; s
4
) = v(48,7 bar ; 6,3726 kJkgK) = 0,0507
12
v
e
= 0,0507 m
3
/kg.
Atunci vom avea pentru volumul specific mediu valoarea :

kg m v
v v v v v v
m
m t e i m
/ 0325 , 0
0325 , 0 0507 , 0 0208 , 0
3
5 4
=
= = = =

Cderea teoretic pe treapta de presiune :

( ) ( )
kg kJ H
K kg kJ bar h h s p h h H
tTP
tTP
/ 13 , 317
/ 3726 , 6 ; 7 , 48 ;
4 4 5 4
=
= =

Prin urmare randamentul treptei de presiune dup o prim estimare va fi :
( )
( ) I a aproximati
CIP
TP
CIP
TP
8466 , 0
8466 , 0 011 , 0 1
20000
600 13 , 317
1
0325 , 0 33 , 273
4 , 0
925 , 0
=
=
|

\
|

+
|

\
|

Cu acest randament determinm entalpia real la ieirea din CIP :

kg kJ h
h h h h
CIP
TP tTP
/ 3062
8466 , 0 13 , 317 5 , 3330
5
5 4 5
=
= =

Cu aceast valoare se recalculeaz volumul specific la ieire :
v
5
=v(p
5
; h
5
) = v(48,7bar ; 3062kJ/kg)
v
5
= 0,053 m
3
/kg
Se recalculeaz volumul specific mediu :
053 , 0 0208 , 0
,
5 4
= =
m m
v v v v v
m
=0,0332m
3
/kg (iteraia II)
i deci noul randament va avea valoarea :
( ) 8478 , 0 011 , 0 1
20000
600 13 , 317
1
0332 , 0 33 , 273
4 , 0
925 , 0 =
|

\
|

+
|

\
|

=
CIP
TP

8478 , 0 =
CIP
TP
(iteraia II)
Pentru acest randament entalpia va avea valoarea :

kg kJ h
h h h h
CIP
TP tTP
/ 6 , 3061
8478 , 0 13 , 317 5 , 3330
5
5 4 5
=
= =

Volumul specific real la ieire va avea valoare :
v
5
=v(p
5
; h
5
) = v(48,7bar ; 3061,6kJ/kg)
13
v
5
= 0,053 m
3
/kg ,
observm c-am obinut acelai volum specific i deci vom avea i acelai
volum specific mediu : v
m
= 0,0332 m
3
/kg (iteraia III).

Deci noul randament va avea valoarea :
( )
8478 , 0
011 , 0 1
20000
600 13 , 317
1
0332 , 0 33 , 273
4 , 0
925 , 0
=
=
|

\
|

+
|

\
|

=
CIP
TP

8478 , 0 =
CIP
TP
(iteraia III)
Cum n aceast ultim iteraie valoarea randamentului a fost identic cu cea
rezultat din iteraia anterioar , o vom considera valoare final .
Entalpia la ieirea din CIP va avea valoarea :

kg kJ h
h h h h
CIP
TP tTP
/ 6 , 3061
8478 , 0 13 , 317 5 , 3330
5
5 4 5
=
= =
;
aceeai cu valoarea corespunztoare iteraiei II .
Punctul 5 : h
5
= 3061,6 kJ/kg ;
p
5
= 48,7 bar ;
t
5
= t(p
5
; h
5
) = t(48,7bar ; 3061,6kJ/kg) t
5
= 345,99
0
C ;
v
5
= v(p
5
; h
5
) = v(48,7bar ; 3061,6kJ/kg) v
5
= 0,053 m
3
/kg;
s
5
= t(p
5
; t
5
) = t(48,7bar ; 345,99
0
C) s
5
= 6,4502 kJ/kgK ;
Se calculeaz randamentul global al CMP - ului:
8345 , 0
13 , 317 17 , 56
8478 , 0 13 , 317 7596 , 0 17 , 56
=
+
+
=

+
+
=
CMP
CMP
tTP tTR
CIP
TP tTP TR tTR CMP
H H
H H

Parametrii la intrarea n supranclzitorul intermediar :
Punctul 6 : p
6
= 47,302 bar ;
de asemenea se consider o pierdere de entalpie de 2 kJ/kg :
14
h
6
= h
5
- 2 = 3061,6 - 2 h
6
= 3059,6 kJ/kg ;
t
6
= t(p
6
; h
6
) = t(47,302bar ; 3059,6kJ/kg) t
6
= 343,92
0
C ;
s
6
=s(p
6
;h
6
)=s(47,302bar ; 3059,6kJ/kg) s
6
= 6,4591kJ/kgK;
v
6
=t(p
6
;s
6
)=v(47,302bar; 6,4591kJ/kgK) v
6
= 0,0545m
3
/kg;

Punctul 7 (ieirea din prenclzitorul intermediar) :
p
7
= 45,804 bar ;
t
7
= 540
0
C (luat din datele iniiale) ;
h
7
= h(p
7
; t
7
) = h(45,804bar ; 540
0
C) h
7
= 3528,1 kJ/kg ;
s
7
= s(p
7
; t
7
) = s(45,804bar ; 540
0
C) s
7
= 7,1356 kJ/kgK ;
v
7
= v(p
7
; t
7
) = v(45,804bar ; 540
0
C) v
7
= 0,078475 m
3
/kg ;

Punctul 8 (intrarea n CMP) :
p
8
= 44,705 bar ;
se va considera o pierdere de entalpie egal cu 2 kJ/kg :
h
8
= h
7
- 2 = 3528,1- 2 h
8
= 3526,1 kJ/kg ;
t
8
= t(p
8
; h
8
) = t(44,705bar ; 3526,1kJ/kg) t
8
= 538,68
0
C ;
s
8
= s(p
8
;h
8
) = s(44,705bar;3526,1kJ/kg) s
8
= 7,1441kJ/kgK;
v
8
= v(s
8
;t
8
)= v(7,1441kJ/kgK;538,68
0
C) v
8
= 0,0813 m
3
/kg;

Randamentul mediu al CMP - ului se calculeaz conform Troyanovski:
|

\
|
|

\
|

+ |

\
|

=
tTP
vs tTP
m m
CMP
H
H
v d

1
20000
700
1
2 , 0
92 , 0 ;
unde : d
m
=
e i m
e i
v v v si
d d
=
+
2
.
Debitul mediu la ieire va avea valoarea :
s kg d
d d
d
m m
/ 43 , 234 43 , 234
2
14 , 216 72 , 252
2
12 8
= =
+
=
+
= ;
15
Pentru calculul volumului specific mediu, trebuie s calculm volumul
specific teoretic la ieire:
( ) ( ) 8425 , 0 / 1441 , 7 ; 6204 , 4 ;
8 12 12
= = = K kg kJ bar v s p v v
t

v
12t
= 0,8425 m
3
/kg .
Ca urmare la nceput ntr-o prim aproximaie vom avea pentru volumul
specific mediu valoarea : 2618 , 0 8425 , 0 0813 , 0
12 8
= = =
t m
v v v
v
m
= 0,2618 m
3
/kg .
Cderea de entalpie pe CMP va avea valoarea :
( )
( ) 43 , 799 / 1441 , 7 ; 6204 , 4 1 , 3526
;
8 12 8
= =
= =
K kg kJ bar h
s p h h h
CMP
t

kg kJ h
CMP
t
/ 43 , 799 = .

Deci randamentul va avea valoarea:

9162 , 0
43 , 799
5 , 3
1
20000
700 43 , 799
1
2618 , 0 43 , 234
2 , 0
92 , 0
=
|

\
|

|

\
|

+
|

\
|

=
CMP
CMP

Cu acest randament calculez entalpia real la ieire :
7 , 2793 9162 , 0 43 , 799 1 , 3526
8 12
= = =
CMP CMP
t
h h h
h
12
= 2793,7 kJ/kg ;
pentru care volumul specific are valoarea:

( ) ( )
kg m v
kg kJ bar v h p v v
/ 4916 , 0
4916 , 0 / 7 , 2793 ; 6204 , 4 ;
3
12
12 12 12
=
= = =

Trebuie recalculat volumul specific mediu:
1999 , 0 4916 , 0 0813 , 0
12 8
= = = v v v
m

v
m
= 0,1999 m
3
/kg .

Cu aceast valoare calculm un nou randament :
16
9163 , 0
43 , 799
5 , 3
1
20000
700 43 , 799
1
1999 , 0 43 , 234
2 , 0
92 , 0
=
|

\
|

|

\
|

+
|

\
|

=
CMP
CMP

Diferena ntre ultimele dou randamente fiind foarte mic , aceast ultim
valoare va fi considerat final i se va trece la calculul definitiv al entalpiei,
volumului specific i volumului specific corespunztoare punctului 12 :
582 , 2793 9163 , 0 43 , 799 1 , 3526
8 12
= = =
CMP CMP
t
h h h
h
12
= 2793,582 kJ/kg ;

( ) ( )
; / 4916 , 0
4916 , 0 / 582 , 2793 ; 6204 , 4 ;
3
12
12 12 12
kg m v
kg kJ bar v h p v v
=
= = =


( ) ( )
; / 3046 , 7
3046 , 7 / 582 , 2793 ; 6204 , 4 ;
12
12 12 12
K kg kJ s
kg kJ bar s h p s s
=
= = =

Avnd randamentul pe CMP i presiunile la prize , am calculat entalpiile
dup o relaie ce ine seama de destinderea real n turbin, destindere ce se
desfoar dup un arc de cerc al crui centru a fost calculat cu metodologia
urmtoare:

( )
( ) ( )
2
0
10 177 , 0
371
650

+ +
)
`

+
=
B B
B
S h h R
h h
S ;
(relaie valabil pentru zona supranclzit ).
Fa de entalpia calculat n mod obinuit ca:
( ) [ ]
( ) [ ]
kg kJ h
K kg kJ bar h h
s p h h h h
CMP
/ 1 , 3069
9163 , 0 / 1441 , 7 ; 6835 , 8 1 , 3526 1 , 3526
;
11
11
8 11 8 8 11
=
=
=

Dar conform literaturii: h
real
< 7 15 kJ/kg fa de valorile determinate ca
mai sus ; h
11real
< h
11
.

La p
10
= 18,497 bar ( ) [ ] =
CMP
s p h h h h
8 10 8 8 10
;
( ) [ ]
. / 3260
9163 , 0 / 1441 , 7 ; 497 , 18 1 , 3526 1 , 3526
10
10
kg kJ h
K kg kJ bar h h
=
=

17

La p
9
= 30,701 bar ( ) [ ] =
CMP
s p h h h h
8 9 8 8 9
;
( ) [ ]
. / 9 , 3406
9163 , 0 / 1441 , 7 ; 701 , 30 1 , 3526 1 , 3526
9
9
kg kJ h
K kg kJ bar h h
=
=

Fa de aceste valori , innd seama de destinderea real sub form de arc de
cerc, se obin:
h
11
= 3061,9 kJ/kg (< 7 15 kJ/kg fa de cel anterior)
s
11
= s(p
11
;h
11
) =s(8,6835bar;3061,9kJ/kg) s
11
= 7,2045kJ/kgK
h
10
= 3253,8 kJ/kg
s
10
= s(p
10
;h
10
) =s(18,497bar;3253,8kJ/kg) s
10
= 7,1706kJ/kgK
h
9
= 3403,2 kJ/kg
s
9
=s(p
9
;h
9
) = s(30,701bar ; 3403,2kJ/kg) s
9
= 7,1537 kJ/kgK

Punctul 13 (punctul de intrare n CJP) :
Am considerat o pierdere de presiune ntre ieire din CMP i
intrare n CJP:
p 2% , astfel nct :
p
13
= 4,5304 bar
Pe conducta CMP CJP am considerat o pierdere de entalpie de cca. 1
kJ/kg, atunci:
h
13
= h
12
- 1= 2793,582-1 h
13
= 2792,582 kJ/kg.
t
13
= t(p
13
; h
13
) = t(4,5304bar ; 2792,582kJ/kg) t
13
= 239,34
0
C;
s
13
= s(p
13
;h
13
) = s(4,5304bar;2792,582kJ/kg)s
13
= 7,1616kJ/kgK;
v
13
= v(p
13
;h
13
) =v(4,5304bar;2792,582kJ/kg) v
13
= 0,4368 m
3
/kg.
Titlul aburului are valoarea: x
abur
= 1 .

Estimarea randamentului pe CJP :
18

CJP
t
umed
umed
tTP CJP
H
H
K
H
|

\
|

+ =
10000
400
1 87 , 0 ;
unde:
- ( ) ( ) K kg kJ bar h h s p h h H
CJP
t
= = / 1616 , 7 ; 0373 , 0 ;
13 13 16 13

kg kJ H H
CJP
t
CJP
t
/ 33 , 700 3 , 2092 582 , 7292 = =
- coeficientul K
umed
utilizeaz relaia din Troyanovski, Filipov:
( )
CJP
t
um
t
hum ev umed
H
H x x
K |

\
|
+
=
2
1 1 8 , 0 1
16 13
,
;
n care
ev,um
= 0,1.
Pentru calculul cderii umede
um
t
H am calculat presiunea de saturaie a
interseciei destinderii teoretice cu curba limit :
p
sat
(s
v
) = 3,1928 bar
h
umed
= h(vap. saturai) - h
16t
= h - h
16t
, unde h este entalpia aburului pe
curba limit : h(p
sat
= 3,1925bar) = 2729,1 h = 2729,1 kJ/kg , astfel c
destinderea teoretic n zona umed va fi :

um
t
H = h - h
16t
= 2729,1-2092,3
um
t
H = 636,86 kJ/kg .
Titlul destinderii teoretice este x
16
= 0,8113 .
Atunci fcnd nlocuirile obinem:
( ) 9382 , 0
33 , 700
86 , 636
2
8113 , 0 1
1 1 , 0 1 8 , 0 1 =
(

+
=
umed umed
K K .
- |

\
|

+
|
|

\
|

1
1 , 0
1
2
10
2
0
3
z flux
c
umed
n
v G
H

;
2
2
2
2
2 2
) sec (
l
l
d
esapare de totala tiunea l d
z
z

=
=

;
Pentru aceast turbin cunosc n
flux
= 4 i lungimea paletei l
p
= 860 mm.
19

|
|
|
|
|

\
|

+
|
|

\
|

=

1
1 , 0
1 10 5
2
2
16 16
4
p
flux flux flux
umed
l
s s n
v D
H

;
unde: s
fl ux
= 6,1 (seciune pe flux , valoare luat din culegerea de
probleme de Centrale Nuclearo-Electrice .
Debitul are valoarea stabilit la subcapitolul anterior : D
16
= 160,13 kg/s .
h
16
=h
13
-h
t
CJP

CJP
, unde
CJP
=
0
K
umed
=
umed
tTP
K
H
|

\
|

+
10000
400
1 87 , 0

CJP
=0,8408 ;
atunci : h
16
= 2792,582 - 700,330,8408 h
16
= 2203,8 kJ/kg ;
v
16
= v(h
16
;p
16
) = v(2203,8kJ/kg ; 0,0373bar) v
16
= 31,8613 m
3
/kg.
Prin urmare:
kg kJ H
H
umed
umed
/ 205 , 23
1
86 , 0
1 , 6
1 , 0
1
1 , 6 4
8613 , 31 13 , 160
10 5
2
2
4
=

|
|
|
|

\
|

+
|

\
|

Recalculez: =

=
33 , 700
205 , 23
8408 , 0
, CJP
CJP
t
umed
CJP
CJP
H
H

8076 , 0 =
CJP
.
Cu noul
CJP
se recalculeaz entalpia real de ieire din turbin:
h
16
= h
13
- h
t
CJP

CJP
= 2792,582 - 700,330,8076 h
16
= 2227 kJ/kg
i noul volum specific va avea valoarea : v
16
= v(h
16
; p
16
) = v(2227kJ/kg ;
0,0373bar)
v
16
= 32,2155 m
3
/kg ;
noul titlu al aburului va avea valoarea :
x
16
= 0,8667.
Cu aceste valori se reia calculul :
20
( ) 9563 , 0
33 , 700
86 , 636
2
8667 , 0 1
1 1 , 0 1 8 , 0 1 =
(

+
=
umed umed
K K .
Se recalculeaz :
CJP
=
0
K
umed
= 0,8961 0,9563 =0,857
CJP
=0,857.
Se recalculeaz H
umed
:
kg kJ H
H
umed
umed
/ 72 , 23
1
86 , 0
1 , 6
1 , 0
1
1 , 6 4
2155 , 32 13 , 160
10 5
2
2
4
=

|
|
|
|

\
|

+
|

\
|

Atunci randamentul va avea valoarea:
=

=
33 , 700
72 , 23
857 , 0
, CJP
CJP
t
umed
CJP
CJP
H
H

8231 , 0 =
CJP
, (iteraia II) .

Refacem calculul cu noul
CJP
:
h
16
= h
13
- h
t
CJP

CJP
=2792,582-700,330,8231 h
16
= 2216,1 kJ/kg;
v
16
=v(h
16
;p
16
)=v(2216,1kJ/kg;0,0373bar) v
16
= 32,0499 m
3
/kg;
x
16
= x(h
16
; p
16
) = v(2216,1kJ/kg ; 0,0373bar) x
16
= 0,8622 ;
Se recalculeaz : K
umed
= 0.9549 ;
CJP
= 0,8557 ; H
umed
= 23,48 kJ/kg ;
noul randament va avea valoarea :
CJP
= 0.8222 (iteraia III).

Se reia din nou calculul i va rezulta : K
umed
= 0.955 ;
CJP
= 0,8558;
H
umed
= 23,5 kJ/kg ;
CJP
= 0.8222 (iteraia IV).
Aceste ultime valori rmn definitive deoarece se poate observa c
randamentul are aceeai valoare.

Prin urmare pentru punctul 16 avem urmtoarele valori :
h
16
=h
13
-h
t
CJP

CJP
=2792,582-700,330,8222h
16
= 2216,89 kJ/kg;
v
16
= v(h
16
;p
16
) =v(2216,89kJ/kg;0,0373bar)v
16
=32,0595 m
3
/kg;
21
x
16
= x(h
16
; p
16
) = x(2216,89kJ/kg ; 0,0373bar) x
16
= 0,8625 ;
s
16
=s(h
16
;p
16
) = s(2216,89kJ/kg;0,0373bar)s
16
=7,3853kJ/kgK;
t
16
= t
saturaie
=27,78
0
C.
Corpul CJP are dou prize , ale cror entalpii sunt:
- la p
15
= 0,31 bar
h
15
= h
13
- ( ) [ ] [ = 582 , 2792 582 , 2792 ;
15 13 15 13
h s p h h
CJP

( )] kg kJ h K kg kJ bar h / 5 , 2434 8222 , 0 / 9716 , 6 ; 31 , 0
15
= ;
la p
14
=1,2684 bar
h
14
= h
13
- ( ) [ ] [ = 582 , 2792 582 , 2792 ;
14 13 14 13
h s p h h
CJP

( )] kg kJ h K kg kJ bar h / 5 , 2608 8222 , 0 / 9716 , 6 ; 2684 , 1
14
= .
Fa de aceste valori , valorile entalpiilor reale conform destinderii dup un
arc de cerc sunt : h
15
= 2429,3 kJ/kg (valoare final) ;
h
14
= 2602,6 kJ/kg (valoare final) .
Cu aceste valori entropiile sunt :
s
14
= s(h
14
;p
14
)= s(2602,6kJ/kg;1,2684bar) s
14
= 7,061kJ/kgK;
s
15
= s(h
15
;p
15
) = s(2429,3kJ/kg;0,31bar) s
15
= 7,1823 kJ/kgK;
Viteza aburului la ieirea din turbin :
s m w w
s n
v d
w
ab ab
flux flux
ab
/ 4 , 210
2 , 6 4
0595 , 32 13 , 160
16 16
=

= ,
Pierderile la eapare sunt de: 39,42 kJ/kg; atunci entalpia real innd cont
de aceste pierderi va avea valoarea :
h
17
= h
16
+ h , (intr efectiv n condensator)
h
17
= 2216,89 +39,42 h
17
= 2256,2 kJ/kg.
v
17
= v(h
17
; p
17
) = v(2256,2kJ/kg ; 0,0373bar) v
17
= 32,66m
3
/kg;
x
17
= x(h
17
; p
17
) = x(2256,2kJ/kg ; 0,0373bar) x
17
= 0,8788 .

Calculul randamentului turbopompei de alimentare :
h
18
= h
12
= 2793,58 kJ/kg ,
22
p
18
= p
12
0,97 p
12
= 4,6204 0,97 p
12
= 4,4818 bar .
s
18
= s(h
18
;p
18
) = s(2793,58kJ/kg; 4,4818bar) s
18
= 7,2658 kJ/kgK;
v
18
= v(h
18
; p
18
) = v(2793,58kJ/kg ; 4,4818bar) v
18
= 0,5149 m
3
/kg;
Abur supranclzit cu x
ab
= 1 .
h
tTPA
= h
18
- h(s
18
; p
19
) h
tTPA
= 2793,6 - h(7,2658m
3
/kg; 0,0373bar)
h
tTPA
= 751,5 kJ/kg .
Calculul randamentului TPA l-am efectuat cu aceeai relaie de la calculul
randamentului CJP - ului :
TPA
t
umed
umed
TPA
t
CJP
H
H
K
H

|
|

\
|

+ =
10000
400
1 87 , 0 ;
unde vom nota:
|

\
|

+ = |

\
|

+ =
10000
400 5 , 751
1 87 , 0
10000
400
1 87 , 0
0 0

TPA
t
H

9006 , 0
0
= .
Presiunea la care destinderea teoretic intersecteaz curba limit : p
sat
=
1,3211 bar , creia i corespunde o entalpie a vaporilor:
h([Link]) = h = 2687,3 kJ/kg .
H
t
um
= h - h
19,t

h
19,t
=h(s
18
;p
19
)h
19,t
=h(7,2658m
3
/kg;0,0373bar) h
19,t
= 2180,8kJ/kg.
Atunci : H
t
um
= 2687,3 - 2180,8 H
t
um
= 506,54 kJ/kg .
Titlul teoretic la ieire: x
19,t
= x(s
18
;p
19
) x
19,t
= x(7,2658m
3
/kg ;
0,0373bar) x
19,t
= x
19
= 0,8477 .
( )
TPA
t
um
t
hum ev umed
H
H x x
K |

\
|
+
=
2
1 1 8 , 0 1
19 18
,

( )
. 9581 , 0
5 , 751
54 , 506
2
8477 , 0 1
1 1 , 0 1 8 , 0 1
=

+
=
umed
umed
K
K

Atunci avem:
TPA
=
0
K
umed

TPA
= 0,9006 0,9581
23

TPA
= 0,8629.
Pentru TPA am considerat n = 1 un singur flux , cu o palet lp = 0,508 m i
o seciune pe flux : s
flux
= 2,43 ; astfel c pierderile prin umiditate :
|
|
|
|
|

\
|

+
|
|

\
|

=

1
1 , 0
1 10 5
2
2
19 19 4
p
flux flux flux
umed
l
s s n
v D
H

;
unde : D
19
= D
18
= 15,17 kg/s ,
v
19
= v(s
18
; p
19
) v
19
= v(7,2658m
3
/kg ; 0,0373bar) v
19
= 31,51m
3
/kg.
kg kJ H
H
umed
umed
/ 3165 , 20
1
508 , 0
43 , 2
1 , 0
1
43 , 2 1
51 , 31 17 , 15
10 5
2
2
4
=

|
|
|
|

\
|

+
|

\
|

Noul:

TPA
=
TPA
- 8358 , 0
5 , 751
3165 , 20
8623 , 0
=
=

TPA
TPA
TPA
t
umed
H
H

.
n urma recalculrii rezult :
TPA
0,8358 (iteraia II) ; calculul va fi
considerat final deoarece nu exist diferen ntre cele dou rezultate, i deci ca
urmare vom avea , procednd ca i n cazul CJP:
h
19
= 2304,8 kJ/kg ;
v
19
= v(h
19
; p
19
) = v(2304,8kJ/kg ; 0,0373bar) v
19
= 33,39 m
3
/kg ;
s
19
= s(h
19
; p
19
) = s(2304,8kJ/kg ; 0,0373bar) s
19
= 7,6757kJ/kgK;
x
19
= x(s
19
;p
19
) x
19
= x(7,6757kJ/kgK/kg;0,0373bar) x
19
= 0,8983;
t
19
= t
17
= 27,78
0
C .
n urma calculului destinderii prin turbin, au rezultat pentru punctele din
turbin valorile prezentate n tabelul 10.1, ce urmeaz.
Ceea ce intereseaz n primul rnd sunt entalpiile pentru a putea face
bilanul pe circuitul de prenclzire regenerativ.
24

Tabel 10.1. Valorile parametrilor termodinamici
Punct

p
[ bar ]
t
[ C ]
h
[ kJ/kg ]
s
[ kJ/kg ]
1 192 540 3376.14 274.1
2 186.11 536.67 3373.1 274.1
3 182.4 535.22 3373.1 273.33
4 151.88 507.55 3330.5 271.37
5 48.7 345.99 3061.6 269.92
6 47.302 343.92 3059.6 252.5
7 45.804 540 3528.1 252.72
8 44.705 538.68 3526.1 252.72
9 30.701 479.96 3403.2 252.28
10 18.497 406.47 3253.8 238.28
11 8.6835 309.61 3061.9 226.2
12 4.6204 239.98 2793.58 216.14
13 4.5304 239.34 2792.58 184.89
14 1.2684 121.12 2602.6 184.89
15 0.31 70.10 2429.3 171.83
16 0.0373 27.78 2216.89 160.13
17 0.0373 27.78 2256.2 160.28
18 4.4818 239.73 2793.58 15.04
19 0.0373 27.78 2304.8 15.04
10.4. Calculul circuitului de prenclzire regenerativ
Schema circuitului de prenclzire regenerativ pentru PIP - uri i PJP - uri
este ce prezentat n figura 1.
25
Din calculul destinderii n turbin se cunosc parametrii termodinamici ai
aburului la prize, respectiv : p , t , h , s , v , x. Aceste valori sunt luate pentru
regimul de referin i necesitatea calculului de prenclzire deriv din motivele
pe care le-am artat anterior i anume :
ne luarea n considerare a pierderilor la labirini ;
neprecizarea explicit de ctre GCSI a tabelelor de ap - abur, dup
care au fost efectuate calculele ;
nu se ine seama de o serie de recuperri de cldur n curent .
S-au mai fcut urmtoarele ipoteze i anume :
Pierderile de presiune pe conductele turbin - prenclzitoare sunt
constante i egale cu 3% ;
Pierderile de presiune pe parte de abur n prenclzitor,
desupranclzitor ca i rcitoare am considerat-o de 1% ;
Pierderile de presiune pe parte de ap le-am considerat de 0,8bar
pentru PIP - uri ; 0,2bar pentru DS ; 0,5bar pentru PJP - uri i 0,1bar
pentru rcitoarele de condensat ;
Pentru temperatura de saturaie a aburului la intrarea n prenclzitor i
apei la ieire le-am considerat de 5,7
0
C pentru PIP-uri i de 2,8
0
C pentru
PJP-uri , innd seama de caracteristica constructiv a acestor
prenclzitoare i anume : oel la PIP-uri (rezisten termic mai mare i
evi mai groase) i pentru PJP-uri am considerat evi din OLC obinuit
(cu grosimi de evi mai mici i rezisten termic mai mic);
Diferena de temperatur ntre ieirea condensatorului din
prenclzitoare i intrarea apei la cele cu rcitor de condensat ca aparat
distinct sau nglobat ;
n cazul PJP - ului 2 , am considerat c ieirea condensatului se face la
saturaie ;
Randamentele principalelor pompe le-am considerat constante cu
ncrcarea (debitul) i anume : = 0,82 pentru pompe de alimentare ; =
26
0,76 pentru pompe de condensat i = 0,68 pentru pompe de condensat
secundar;
Randamentul mecanic al cuplajului pomp-TPA l-am considerat
mec
=
0,98 ;
Pierderile de presiune n staia de tratare chimic le-am considerat de 4
bar i o presiune maxim admis de staia de tratare de 6 bar, ceea ce mi-a
permis dimensionarea PC treapta I, iar pentru PC treapta II am fcut
ipoteza: pierderile de presiune n capul de alimentare ale degazorului sunt
de 20% ;
Am considerat c 85% din pierderile la generator i pierderi mecanice
transmise hidrogenului respectiv uleiului, se recupereaz n circuit i
anume ntre pompele PC treapta I i PC trepta II .

Schema i procentul de recuperare sunt justificate i de indicaiile din
,,CEI ca i din ,,Spencer i Cotton . Aceste recuperri de cldur, pentru
uurina calcului, au fost raportate la debitul de abur viu (proporional unor
energii specifice) justificabil prin aceea c att pierderile n generator ct i cele
mecanice depind de puterea brut a grupului, putere brut care n esen depinde
de debitul intrat n turbin ; cu att mai mult cu ct ipotezele de calcul au fost cu
parametrii constani ai aburului viu la intrare i temperatura de supranclzire
constant (am variat doar presiunea de condensare).
Valorile acestor energii specifice recuperate (valori raportate la debit) sunt:
pentru pierderi generator 16,9381 kJ/kg i pentru pierderi mecanice sunt 4,067
kJ/kg .
Am luat n considerare recuperarea aburului , att sub form de cantitate de
cldur considerat raportat la debit de 7,281 kJ/kg i sub form de debit
recuperat n condensator ca o cot de 2,43310
-3
(=d
20
) din debitul de abur viu .
n schema de prenclzire am inut seama de asemenea , de cantitatea de cldur
27
echivalent de la condensatorul de abur ejector cu o cot q
k
= 7,281 (extrage abur
de la prenclzitoare).
Recuperarea pierderilor la labirini de nalt presiune respectiv de la intrarea
n CIP i n CMP se face conform schemei grupului de 330 MW n degazor i
anume: d
49
=
1 4 2
2
1
7
2
d b b |

\
|
+ . Entalpia corespunztoare a fost considerat
ca medie ponderat a entalpiilor debitelor respective , deci :
h
49
=
4 2
8 4 3 2
2
1
7
2
2
1
7
2
b b
h b h b
+
+
.
O parte din aburul de la labirinii CIP - ului se recupereaz n PJP-ul 3 ,
formnd debitul :
1 3 2 48
2
1
7
4
d b b d |

\
|
+ = . Similar ca mai sus, entalpia va
fi (ca medie ponderat): h
48
=
3 2
5 3 3 2
2
1
7
4
2
1
7
4
b b
h b h b
+
+
.
Surplusul de abur de la regulatorul abur-labirini se recupereaz n PJP-ul 1
i anume :
1 6 5 4 3 2 50
2
1
2
1
7
1
d b b b b b d |

\
|
+ + + = respectiv:
h
50
=
5 4 3 2
12 5 8 4 5 3 3 2
2
1
2
1
7
1
2
1
2
1
7
1
b b b b
h b h b h b h b
+ + +
+ + +
.
Toate aceste debite recuperate am considerat c au presiunea prizei
respective (vezi justificare i de la cotele de debit).

Observaie : Cotele procentuale recuperate n degazor , PJP 3 iPJP1au fost
considerate dup debitele date de ICEMENERG i dup cotele prezentate de
Rizhkin (pentru grupuri de 800 MW).
Am considerat c din de-supranclzitoare iese abur supranclzit cu 20
0
C
(deci t
sat
+ dt
x
; dt
x
= 20
0
C ).
28

Scopul final
Scopul final al acestui calcul este acela de a determina ct mai exact
debitele de abur necesare la prize pentru a calcula corect debitele ce circul prin
fiecare seciune de turbin , precum i volumele specifice la intrare n grupul de
trepte i presiunile, care se constituie date de referin pentru regimuri
nenominale dar staionare ale centralei.
Pentru scopul actualei cercetri, aceste regimuri nenominale constituie de
fapt variaia presiunii de condensare.
De asemeni, pentru comparaie, am efectuat calculul funcionrii cu EPA .
Acest calcul se justific prin aceea c condiiile reale de funcionare a grupului
de 330 MW impun restricii datorit puterii calorifice inferioare a
combustibilului , astfel c sarcina uzual este de aproximativ 200 250 MW .

Calculul propriu-zis al circuitului :
Bilanul pe prenclzitorul 7 :
47 58 37 37 36 36 58 58
h d d h h d h d + = +
unde : h
58
= h(p
58
; t
58
) .
t
58
= t
sat
(p
58
) + dt
x
; dt
x
= 20
0
C (ipotez)
p
58
= 0,970,99p
5

p
58
= 0,970,9948,7
p
58
= 46,9893 bar (conform
ipotezei)
t
58
= t
sat
(46,9893 bar)+ 20 t
58
= 260,17 + 20
t
58
= 280,17
0
C
h
58
= h(46,9893 bar ; 280,17
0
C)
h
58
= 2870,28 kJ/kg
h
37
= h(p
37
; t
37
)
29
unde : p
37
= p
al
+ p
DS6
+ p
DS7
; dar p
al
= 1,33 p
0
p
al
= 1,33 192
p
al
= 255,36 bar.
Atunci presiunea p
37
va avea valoarea : p
37
= 255,36 + 0,2 + 0,2
p
37
= 255,76bar.
t
37
= t
sat
(p
58
) - dt
x
; dt
x
= 5,7
0
C
t
37
= 260,17 - 5,7 t
37
= 254,4
0
C.
nlocuind n relaia entalpiei obinem:
h
37
= h(255,76bar ; 254,4
0
C) h
37
= 1107,4 kJ/kg .
h
36
= h(p
36
; t
36
)
unde: p
36
= p
37
+ p
PIP
p
36
= 255,76 + 0,8 p
36
= 258,5bar .
t
36
= t
sat
(p
59
) - dt
x
.
p
59
= 0,97 0,99 p
9
p
59
= 0,97 0,99 30,701
p
59
= 29,48bar .
t
36
= t
sat
(29,48bar) - 5,7
t
36
= 233,7 - 5,7 t
36
= 228
0
C .
i atunci:
h
36
= h(258,8bar ; 228
0
C) h
36
= 985,7 kJ/kg .
h
47
= h(p
47
; t
47
)
unde : p
47
= 0,99 p
58
p
47
= 0,99 46,9893 p
47
= 46,5bar .
t
47
= t
36
+ 2 t
47
= 228 + 2 t
47
= 230
0
C .
h
47
= h(46,5bar ; 230
0
C) h
47
= 993,6 kJ/kg .
Dar : d
36
= d
37
= d
1
= (= d
33
= = d
39
) = 274,1 kg/s .
Din bilanul pe prenclzitorul 7 avem relaia :

=
6 , 993 28 , 2870
1 , 274 7 , 985 1 , 274 4 , 1107
58
47 58
36 36 37 37
58
d
h h
d h d h
d
d
58
= 17,768 kg/s
Atunci: a
58
= 0,064847
Dar : d
47
= d
58
= 17 768 kg/s
30
a
47
= 0,064867.

Bilanul pe de-supranclzitorul DS 7 :
d
37
h
37
+ d
57
h
57
= d
38
h
38
+ d
58
h
58

h
38
=
38
58 58 57 57 37 37
d
h d h d h d +
;
unde : d
38
= d
37
= 274,1 kg/s ,
atunci relaia de mai sus devine:
( )
58 57
38
57
37 38
h h
d
d
h h + = ;
unde : d
57
= d
58
=17,768 kg/s ,
h
57
= 3062,7 kJ/kg (luat din destindere)
( ) 28 , 2870 7 , 3062
1 , 274
768 , 17
4 , 1107
38
+ = h
h
38
= 1119,87 kJ/kg .
p
38
= p
al
+ p
DS6
p
38
= 255,36 + 0,2 p
38
= 257,5bar
t
38
= t(p
38
; h
38
) = t(257,5bar ; 1119,87kJ/kg) t
38
= 257
0
C .

Bilanul pe prenclzitorul 6 :
d
35
h
35
+ d
47
h
47
+ d
59
h
59
=
= d
36
h
36
+ (d
59
+ d
47
) h
46

unde : h
59
= 2866,06 kJ/kg (din destindere)
h
47
= 993,62 kJ/kg , d
47
= 17,768 kg/s ,
d
35
= d
34
= 274,1 kg/s ,
p
35
= p
36
+ p
PIP
p
35
= 258,5 + 0,8
p
35
= 259,3 bar ,
t
34
= t
sat
(p
60
) - dt
x
, dt
x
= 5,7
0
C ,
p
60
= 0,97 (0,99)
4
p
10

p
60
= 0,97 (0,99)
4
18,497
31
p
60
= 17,24 bar .
Atunci temperatura va avea valoarea:
t
34
= t
sat
(17,24bar) - 5,7 t
34
= 207,5 -5,7
t
34
= 201,35
0
C ,
t
35
= t
34
+ t
DS
t
35
= 201,35 + 4 t
35
= 205,35
0
C ,
h
35
= h(p
35
; t
35
) h
35
= h(259,3bar ; 205,35
0
C) h
35
= 886,2 kJ/kg ,
h
46
= h(p
46
; t
46
) , unde : p
46
= 0,99 p
59
p
46
= 0,99 0,97 0,99 p
9
p
46

= 0,99
2
0,97 30,701 p
46
= 29,5 bar ,
t
46
= t
35
+ 2 t
46
= 205,35 + 2 t
46
= 207,35
0
C ,
h
46
= h(29,5bar ; 207,35
0
C) h
46
= 890,2 kJ/kg .
nlocuim n bilan:

46 59
35 35 47 47 46 47 36 36
59
h h
h d h d h d h d
d

+
=

2 , 890 06 , 2866
2 , 886 1 , 274 62 , 993 768 , 17 2 , 890 768 , 17 7 , 985 1 , 274
59

+
= d
d
59
= 12,864 kg/s , atunci: a
59
= 0,046947 ;
Dar d
56
= d
59
= 12,864 kg/s , atunci : a
56
= 0,046947 ;
d
46
= d
59
+ d
47
d
46
= 12,864 + 17,768 d
46
= 30,63 kg/s ,
atunci : a
46
= 0,111794 .

Bilan pe de-supranclzitorul 6 :
d
38
h
38
+ d
56
h
56
= d
39
h
39
+ d
56

h
59

h
39
= ( )
59 56
38
56
38
h h
d
d
h + ,
h
56
= 3408,18 kJ/kg (din destindere)


32
( ) 06 , 2866 18 , 3408
1 , 274
864 , 12
87 , 1119
39
+ = h
h
39
= 1145,3 kJ/kg .
t
39
= t(p
39
; h
39
)
t
39
= t(255,36bar; 1145,3kJ/kg) t
39
= 262,44
0
C.
Ceea ce nseamn o cretere de temperatur pe DS 6 de :
t
39
- t
38
= 262,44 - 257 t
39
- t
38
= 514
0
C = t
DS6
;
ceea ce concord cu valorile cunoscute (din ipoteze).
Presiunea efectiv de intrare n cazan dup capul de alimentare este :
p
40
= 0,97 p
39

p
40
= 0,97 255,36 p
40
= 249 bar.
Entalpia h
40
= h
39
(un proces de laminare n ventile)
h
40
= 1145,3 kJ/kg.
d
40
= d
39
- a
injo
d
1

d
40
= 274,1 - 0,001 274,1 d
40
= 273,73 kg/s.

Bilanul pe prenclzitorul 5 :
d
60
h
60
+ d
46
h
46
+ d
33
h
33
= d
34
h
34
+ (d
60
+ d
46
) h
45
,
unde d
60
+ d
46
= d
45
.
45 60
46 46 33 33 45 46 34 34
60
h h
h d h d h d h d
d

+
=
h
60
= 2851,1 kJ/kg (din destindere)
t
34
= 201,35
0
C
p
34
= p
35
+ p
DS
p
34
= 259,5 bar
h
34
= h(p
34
; t
34
) h
34
= 868,6
kJ/kg
d
46
= d
59
+ d
47
d
46
= 12,864 +
17,768
33
d
46
= 30,63 kg/s
h
46
= 890,2 kJ/kg .


Punctul ,,33 :calculul parametrilor apei la intrarea n prenclzitoarele de
nalt presiune :
Parametrii apei la ieirea din degazor :
p
32
= 0,99 p
54
p
32
= 0,99 0,97p
11

p
32
= 0,99 0,97 8,6835 p
32
= 8,46bar
t
32
= t
sat
(p
32
) t
32
= t
sat
(8,46bar) t
32
= 172,7
0
C
h
32
= h
sat
(p
32
) h
32
= h
sat
(8,46bar) h
32
= 731,2 kJ/kg .
Creterea de entalpie n pomp :
( )
PA
pompa
v v p p
h

33 32 32 33
100
= ,
unde:
PA
= 0,82 (conform ipotezelor)
p
33
= p
34
+dp
PIP
p
33
= 259,5 + 0,8 p
33
= 260,3 bar
v
32
= v
sat
(p
32
) v
32
= v
sat
(8,46bar) v
32
= 0,001118m
3
/kg
v
33
= v(p
33
;t
33
) v
33
= v(260,3bar;177,7
0
C) v
33
= 0,0011048 m
3
/kg
unde: t
33
= t
32
+ dt
pomp
,
se estimeaz dt
pomp
= 5
0
C t
33
= 172,7+ 5
t
33
=177,7
0
C.
Atunci creterea de entalpie n pomp va fi :
( )
. / 2 , 34
82 , 0
0011048 , 0 001118 , 0 100 46 , 8 3 , 260
kg kJ h
h
pompa
pompa
=

=

h
33
= h
32
+ h
pomp
h
33
= 731,2 + 34,2 h
33
= 765,4 kJ/kg .
Se recalculeaz t
33
= t(p
33
; h
33
)
t
33
= t(260,3bar ; 765,4kJ/kg) t
33
= 177,5
0
C .
34
Se poate observa c diferena ntre valoarea anterioar i cea obinut acum
este nesemnificativ , deci aceast ultim valoare va fi considerat final . Nu
este nevoie s refacem calculul pentru volumul specific (am estimat corect
creterea de temperatur pe pomp).
t
45
= t
33
+ dt
PIP
t
45
= 177,5 + 2 t
45
= 179,5
0
C
p
45
= 0,99 p
60
p
45
= 0,99 17,24 p
45
= 17,8 bar
h
45
= h(p
45
; t
45
) h
45
= h(17,8bar ; 179,5
0
C) h
45
= 765,3 kJ/kg
Din bilan se va calcula debitul masic d
60
:
s kg d
d
/ 734 , 11
3 , 765 1 , 2851
2 , 890 63 , 30 4 , 765 1 , 274 3 , 765 63 , 30 6 , 868 1 , 274
60
60
=

+
=

Atunci : a
60
= 0,042825. Dar : d
55
= d
60
= 11,734 kg/s .

Bilanul pe desupranclzitorul DS 5:
d
55
h
55
+ d
34
h
34
= d
35
h
35
+ d
60
h
60

( )
60 55
35
60
34 35
h h
d
d
h h + =
h
55
= 3262,7 kJ/kg .
nlocuind n relaie se obine:
( )
. / 2 , 886
1 , 2851 7 , 3262
1 , 274
734 , 11
6 , 868
35
35
kg kJ h
h
=
+ =

t
35
= t(p
35
; h
35
) t
35
= t(259,3bar ; 886,2kJ/kg) t
35
= 205,37
0
C.
Prin urmare creterea de temperatur pe DS 5 este :
t
DS5
= t
35
-t
34
t
DS5
=205,37-201,35
t
DS5
= 4,03
0
C, valoare apropiat de
cea pe care am estimat-o anterior , prin
urmare nu este necesar s relum calculele.
Bilanul pe degazor
35
Din bilanul pe degazor:
d
31
= d
32
- d
54
- d
49
- d
45

Din bilanul energetic :
d
54
h
54
+ d
49
h
49
+ d
45
h
45
+
+(d
32
- d
54
-d
49
- d
45
) h
31
= d
32
h
32


( )
31 54
31 45 49 32 45 45 49 49 32 32
54
h h
h d d d h d h d h d
d

= ;
unde : h
45
= 765,3 kJ/kg ,
p
54
= 0,97 p
11
p
54
= 0,97 8,6835 p
54
= 8,4229 bar,
d
45
= d
60
+ d
46
= d
60
+ d
59
+ d
58
= d
60
+ d
56
+ d
58

d
45
=11,734+ +12,864+17,768
d
45
= 42,37 kg/s ,
a
45
= 0,15457862 ,

0016 , 0
2
1
0072 , 0
7
2
1 , 3526 0016 , 0
2
1
1 , 3373 0072 , 0
7
2
0016 , 0
2
1
7
2
2
1
7
2
h
2
8 4 3 2
49
+
+
=
+
+
=
b
h b h b

s kg / 99 , 3415 =
49
h ,
d
49
= |

\
|
+ = |

\
|
+ 1 , 274 0016 , 0
2
1
0072 , 0
7
2
2
1
7
2
49 1 4 2
d d b b
s kg d / 783 , 0
49
= ,
d
32
= d
33
+ a
inj SI
d
1
,
unde a
inj SI
reprezint debitul de injecie n supranclzitorul intermediar.
d
32
= 274,1+274,1 0,0008 d
32
= 274,22 kg/s ,
h
11
= 3073,96 kJ/kg (din destindere) ,
t
31
= t
sat
(p
53
)-dt
PJP
; p
53
=0,97p
12
p
53
=0,974,6204 p
53
=4,4817bar
atunci :
t
31
= t
sat
(4,4817bar)-2,8 t
31
=148,51-2,8 t
31
=145,71
0
C ,
p
31
= 1,2 p
54
p
31
= 1,2 8,4229 p
31
= 10,25 bar ,
36
h
31
= h(p
31
;t
31
) h
31
= h(10,25bar;145,71
0
C) h
31
= 613,9 kJ/kg.
Prin nlocuire se va obine :
, / 577 , 9
9 , 613 96 , 3073
9 , 613 37 , 42 9 , 613 783 , 0
9 , 613 96 , 3073
9 , 613 22 , 274 3 , 765 37 , 42 99 , 3415 783 , 0 22 , 274 2 , 731
54
54
s kg d
d
=

+
+
+

=

a
54
= 0,034952 ( a
pa
) .

[Link] pe prenclzitorul 3:
Bilanul masic :
d
44
= d
48
+ d
53

Bilanul energetic :
d
53
h
53
+ d
48
h
48
+ d
30
h
30
=
= d
30
h
31
+ (d
48
+ d
53
) h
44

44 53
30 30 48 48 31 30 44 48
53
h h
h d h d h d h d
d

+
= ;
Din datele iniiale:
. / 52 , 3251
0053 , 0
2
1
0072 , 0
7
4
6 , 3061 0053 , 0
2
1
1 , 3373 0072 , 0
7
4
2
1
7
4
2
1
7
4
. / 854 , 1
1 , 274 0053 , 0
2
1
0072 , 0
7
4
2
1
7
4
48
3 2
5 3 3 2
48
48
1 3 2 48
kg kJ h
b b
h b h b
h
s kg d
d b b d
=
+
+
=
+
+
=
=
|

\
|
+ = |

\
|
+ =

p
53
= 0,97 p
12

p
53
= 0,97 4,6204 p
53
= 4,57 bar.
h
53
= 2940,4 kJ/kg (din destindere) .
t
53
= t(p
53
; h
53
) t
53
= t(4,57bar ; 2940,4kJ/kg) t
53
= 239,73
0
C.
t
30
= t
29
+ dt
PCs
,
37
unde : dt
PCs
= contribuie asupra temperaturii a pompelor de condensat
secundar , estimez dt
PCs
= 0,4
0
C ,
t
29
= t
sat
(p
52
)- dt
PJP
, p
52
= 0,97 p
14
p
52
= 0,97 1,2684
p
52
= 1,2303 bar ,
t
29
= t
sat
(1,2303bar) - 2,8 t
29
= 103,05
0
C ,
t
30
= 103,05 + 0,4 t
30
= 103,45
0
C ,
p
30
= p
31
+ dp
PJP
+ dp
rc

p
30
= 10,25+ 0,5+ 0,1 p
30
= 10,85 bar ,
h
30
= h(p
30
; t
30
) h
30
= h(10,85bar ; 103,45
0
C) h
30
= 434,03 kJ/kg ,
p
44
= 0,99 p
53
p
44
= 0,99 0,97 p
12

p
44
= 0,99 0,97 4,6204 p
44
=4,5265bar
t
44
= t
30
+ dt
PJP
t
44
= 103,45+ 5 t
44
= 108,45
0
C ,
h
44
= h(p
44
; t
44
) h
44
= h(4,5265bar ; 108,45
0
C) h
44
= 454,7 kJ/kg ,
Prin nlocuire vom avea :
, / 94 , 13
7 , 454 4 , 2904
03 , 434 49 , 221 52 , 3251 854 , 1 9 , 613 49 , 221 7 , 454 854 , 1
53
53
s kg d
d
=

+
=

a
p3
= 0,05087,
d
44
= d
48
+ d
53
d
44
= 1,854+ 13,94 d
44
= 15,794 kg/s .

Bilanul pe prenclzitorul 2 :
Din bilanul masic n punctul de
recirculare :
d
29
= d
28
= d
30
- d
48
- d
53
- d
52
.
Bilanul energetic:
d
52
h
52
+ (d
30
- d
48
- d
53
- d
52
) h
28

+
+ d
44
h
44
= (d
30
- d
48
- d
53
- d
52
) h
29

38
+ (d
52
+ d
44
)h
42

( ) ( )
42 29 28 52
44 44 28 53 48 30 42 44 29 53 48 30
52
h h h h
h d h d d d h d h d d d
d
+
+
=
unde : h
29
= h(p
29
; t
29
) h
29
= h(10,85bar ; 103,05
0
C) h
29
=
=432,54 kJ/kg ,
d
44
= 15,794 kg/s (parametru calculat) ,
p
42
= 0,99 p
52
p
42
= 0,99 0,97 p
14
p
42
= 0,99 0,97 1,26842
p
42
=1,218bar.
Deoarece acest prenclzitor nu are rcitor de condensat pentru c nu se
justific subrcirea deoarece condensatul se re-pompeaz n circuitul de
prenclzire dup PJP 2 , condensul iese saturat . Prin urmare :
t
42
= t
sat
(p
42
) t
42
= t
sat
(1,218 bar) t
42
= 105,56
0
C ,
p
28
= p
29
+ dp
PJP
p
28
= 10,85+ 0,5 p
28
= 11,35 bar ,
t
28
= t
sat
(p
51
) - dt
PJP
, unde : p
51
= 0,97 p
15
p
51
= 0,97 0,31
p
51
= =0,3036bar,
t
28
= t
sat
(0,3036bar) - 2,8 t
28
= 66,6
0
C ,
h
28
= h(p
28
; t
28
) h
28
= h(11,35bar ; 66,6
0
C) h
28
= 279,5 kJ/kg .
Prin nlocuire se va obine :
( )
( )
, / 899 , 12
56 , 442 54 , 432 5 , 279 63 , 2715
5 , 279 94 , 13 854 , 1 49 , 221
56 , 442 54 , 432 5 , 279 63 , 2715
7 , 454 794 , 15 56 , 442 794 , 15 54 , 432 94 , 13 854 , 1 49 , 221
52
52
s kg d
d
=
+

+
+
=

a
52
= a
p2
= 0,047076 .
Deci pot calcula debitul :
d
42
= d
44
+ d
52
d
42
= 15,794+ 12,899 d
42
= 28,7 kg/s ,
dar : d
42
= d
43
d
43
= 28,7 kg/s ,
unde d
43
reprezint debitul dup pompa de condensat secundar .
d
28
= d
29
= d
30
- d
42
d
28
= d
29
= 221,49 - 28,7 d
28
= d
29
= 192,8 kg/s .
Creterea de entalpie n pompa de condensat secundar .
39
Bilanul pe pomp :
( )
PCS
PCs
v v p p
h

42 43 42 43
100
= ;
unde
PCs
= =0,68 (din ipoteze) .
p
42
= 0,99 p
52
p
42
= 0.99 0,97 p
14
p
42
= 1,218 bar .
Dar p
43
= p
30
= p
29
= 10,85 bar .
v
42
= v
sat
(p
42
) v
42
= v
sat
(1,218bar) v
42
= 0,001048 m
3
/kg .
Pentru calculul volumului specific de ap dup pomp , v
43
, estimez pentru
nceput o cretere de temperatur n pomp de : t
PCs
= =0,3
0
C; atunci vom avea
pentru volumul specific valoarea :
v
43
= v(p
43
; t
43
) v
43
= v(10,85bar ; 105,86
0
C) v
43
= 0,0010795 m
3
/kg,
unde t
43
= t
42
+ t
PCs
t
43
= 105,56+ 0,3 t
43
= 105,86
0
C .
( )
. / 507 , 1
68 , 0
001048 , 0 0010795 , 0 100 218 , 1 85 , 10
kg kJ h h
PCs PCs
=

=

Prin urmare creterea de temperatur n PCse va avea valoarea :
t
43
- t
42
= 105,86 - 105,56 = 0,3 ;
se observ identitatea cu valoarea considerat anterior , i atunci nu mai
sunt necesare iteraii pentru calculul creterii de entalpie pe pompa de condensat
secundar .

Concluzii :
Pentru punctul de recirculare: bilanul energetic n acest punct conduce la:
. / 03 , 434
49 , 221
89 , 443 7 , 28 8 , 192 54 , 432
30
30
31
43 43 29 29
30
kg kJ h
h
d
h d h d
h
=

+
=
+
=

t
30
= t(p
30
; h
30
) t
30
= t(10,85bar ; 434,03kJ/kg) t
30
= 103,4
0
C , aceeai
cu valoarea propus anterior , deci nu mai este necesar s reiau calculul , altfel
era nevoie s reiterez de la bilanul pe PJP- ul 3 i PJP- ul 2 .
40

Bilanul pe prenclzitorul 1 :
d
50
h
50
+ d
51
h
51
+ d
27
h
27
=
= d
27
h
28
+ (d
50
+ d
51
)h
41

41 51
50 50 27 27 28 27 41 50
51
h h
h d h d h d h d
d

+
= .
Pentru d
50
am estimat anterior (n ipoteze) ; unde d
50
reprezint surplusul de
abur de la regulatorul abur - labirini .
)
. / 4 , 3199
0003 , 0
2
0016 , 0
2
0053 , 0
7
0072 , 0
7 , 2793 0003 , 0 1 , 3526
2
0016 , 0
6 , 3061
2
0053 , 0
1 , 3373
7
0072 , 0
2 2 7
2 2 7
. / 614 , 1 1 , 274 00053 , 0 0003 , 0
2
0016 , 0
2
0053 , 0
7
0072 , 0
2 2 7
50
50
5
4 3 2
12 5 8
4
5
3
3
2
50
50
50 1 6 5
4 3 2
50
kg kJ h
h
b
b b b
h d h
b
h
b
h
b
h
s kg d
d d b b
b b b
d
=

+ + +
+ + +
=

+ + +
+ + +
=
= +

\
|
+ + + =
|

\
|
+ + + =

h
51
= 2530,6 kJ/kg (din destindere) , unde t
28
= 66,6
0
C , (calculat anterior) .
p
27
= p
28
+ p
PJP1
+ p
rc.
p
27
= 11,35+ 0,5+ 0,1 p
27
= 11,95 bar ; unde
p
27
reprezint presiunea apei la ieirea din treapta 2 a pompei de condens
principal .
Pentru calculul temperaturii este necesar s se calculeze bilanul pe
recuperatoare i s se in seama de creterea de entalpie n PC
tr1
i PC
tr2
.
Din datele anterioare avem parametrii condensului la ieirea din
condensator , i anume :
41
Punctul 21 : p
21
= 0,0373 bar ;
t
21
= t
sat
(p
21
) t
21
= t
sat
(0,0373bar) t
21
= 27,78
0
C;
h
21
=h
sat
(p
21
)h
21
=h
sat
(0,0373bar)h
21
=116,42kJ/kg.
Calculul creterii de entalpie n PC trepta 1 :
Cum prin ipotez am considerat presiunea maxim la care rezist
rezervoarele staiei de tratare chimic este 6 bar , rezult c am impus p
25
= 6 bar
(la intrarea n staie).
innd seama de pierderile de presiune n rcitoarele de hidrogen i ulei , de
0,2 bar , rezult c presiunea la ieirea treapta 1 de pompare de condensat
principal (3 x 55%) :
p
22
= p
25
+ 0,2 p
22
= 6 + 0,2 p
22
= 6,2 bar .
( )
1 ,
21 22 21 22
1 ,
100
tr PC
tr PC
v v p p
h

= ; unde
PC,tr1
= 0,76 (din ipoteze) .
v
21
= v
sat
(p
21
) v
21
= v
sat
(0,0373bar) v
21
= 0,0010039 m
3
/kg .
Estimez o cretere de temperatur dup PC treapta 1 de: t = 0,1
0
C.
t
22
= t
21
+ t t
22
= 27,78+ 0,1 t
22
= 27,88
0
C .
v
22
= v(t
22
; p
22
) v
22
= v(27,88
0
C ; 6,2bar) v
22
= 0,0010033 m
3
/kg .
Atunci prin nlocuire avem :
( )
. / 65 , 0
76 , 0
0010033 , 0 0010039 , 0 100 0373 , 0 2 , 6
1 ,
1 ,
kg kJ h
h
tr PC
tr PC
=

=

Se trece la iteraia 2 : t
22
= t(p
22
; h
22
); unde h
22
= h
21
+ h
PC,tr1
h
22
=
=116,42+ 0,65 h
22
= 117,09 kJ/kg .
t
22
= t(6,2bar ; 117,09kJ/kg) t
22
= 27,82
0
C ,valoarea este foarte
apropiat de cea corespunztoare iteraiei anterioare , deci va rmne ca valoare
final.
Creterea de temperatur pe PC treapta 1 va fi :
42
t = t
22
- t
21
t = 27,82 - 27,78 t = 0,04
0
C , prin urmare nu se mai
reitereaz .
Pentru a calcula temperatura la intrare n staia de tratare chimic a apei ,
fac bilan global pe recuperatoarele de cldur , corespunztor pierderile
mecanice , n generator i condensului de la abur labirini (CAE):
( )
22
20 20 1 1 9
22 25
85 , 0
d
h d d q d p p
h h
CAE dm d
+ +
+ = ; unde ,,0,85
apare n ipoteza c 85% din pierderi se recupereaz .
d
20
- am luat n considerare recuperarea aburului din etanri sub form de
debit recuperat n condensator ca o cot de 2,433 x 10
-3
din d
1
.
d
20
= (2,433 x 10
-3
) d
1
d
20
= 2,433 10
-3
274,1 d
20
= 0,667 kg/s .
t
20
= t
22
+ dt
2
t
20
= 27,82+ 5 t
20
= 32,82
0
C .
h
20
= h
sat
(t
20
) h
20
= h
sat
(32,82
0
C) h
20
= 137,47 kJ/kg .
Ca urmare se poate calcula h
25
:
( )
803 , 192
47 , 137 667 , 0 1 , 274 281 , 7 1 , 274 067 , 4 9381 , 16 85 , 0
09 , 117
25
+ +
+ = h
h
25
= 152,3 kJ/kg .
t
25
= t(p
28
; h
25
) t
25
= t(6bar ; 152,3kJ/kg) t
25
= 36,27
0
C .
Recuperatoarele : rcire ulei i rcire generator , sunt puse n paralel cu
condensatorul abur-ejector n schem i mprirea debitului de condens prin
acestea se face din ipoteza c la ieirea din cele dou entalpia este constant i
egal cu h
25
.
h
25
= h
23
= h
24
= 152,3 kJ/kg .
( )
( )
. / 8 , 138
09 , 117 3 , 152
067 , 4 9381 , 16 274 85 , 0
85 , 0
24
24
22 25
1
24
s kg d
d
h h
p p d
d
dm dg
=

+
=

+
=

d
25
= d
22
= d
21
= 192,803 kg/s .
d
23
= d
25
- d
24
d
23
= 192,803 - 138,8 d
23
= 54 kg/s .
43
n staia de tratare chimic a apei n-am considerat pierderi de temperatur
ci numai de presiune :
p
26
= p
25
- p
ST
p
26
= 6 - 4 p
26
= 2 bar .
t
25
= t
26
= 36,27
0
C .
h
26
= h(p
26
; t
26
) h
26
= h(2bar ; 36,27
0
C) h
26
= 151,97 kJ/kg .
v
26
= v(p
26
; t
26
) v
26
= v(2bar ; 36,27
0
C) v
26
= 0,0010063 m
3
/kg .
Creterea de entalpie n pompa de condensat principal, treapta 2 :
( )
2 ,
26 27 26 27
2 ,
100
tr PC
tr PC
v v p p
h

= ;
unde :
PC,tr2
= 0,76 (din ipoteze).
Estimez o cretere de temperatur n pomp : t
PC,tr2
= 0,1
0
C .
t
27
= t
26
+ t
PC,tr2
t
27
= 36,27 + 0,1 t
27
= 36,37
0
C .
v
27
= v(p
27
;t
27
) v
27
= v(11,95bar;36,37
0
C) v
27
= 0,0010059 m
3
/kg .
nlocuind obinem:

( )

=
76 , 0
0010063 , 0 0010059 , 0 100 2 95 , 11
2 ,tr PC
h
h
PC,tr2
= 1,32 kJ/kg .
h
27
= h
26
+ h
PC2
h
27
= 151,97 + 1,32 h
27
= 153,29 kJ/kg .
t
27
= t(p
27
; h
27
) t
27
= t(11,95bar ; 153,29kJ/kg) t
27
= 36,37
0
C .
Creterea de temperatur pe PC treapta 2 : t
PC,tr2
= t
27
- t
26
t
PC,tr2
=
36,37-36,27 t
PC,tr2
= 0,1
0
C , deci am estimat bine , astfel nct nu mai este
nevoie s reiau calculul .
Punctul 41: p
41
=0,99p
51
p
41
=0,990,97p
15
p
41
=0,990,970,31
p
41
= 0,03006 bar .
t
41
= t
27
+ dt
2
t
41
= 36,37+ 5 t
41
= 41,37
0
C .
h
41
= h(p
41
; t
41
) h
41
= h(0,03006bar ; 41,37
0
C)
h
41
= 173,13 kJ/kg .
44
Am stabilit toi parametrii apei , aburului i condensului intrat i ieit n
PJP 1 , deci pot calcula debitul necesar la priza 1 din schem:
. / 826 , 8
13 , 173 7 , 2530
4 , 3199 614 , 1 29 , 153 803 , 192 48 , 279 803 , 192 13 , 173 614 , 1
51
51
s kg d
d
=

+
=

a
p1
= 0,032212 .
d
41
= d
50
+ d
51
d
41
= 1,164 + 8,826 d
41
= 9,99 kg/s .
n concluzie , n urma acestui calcul au rezultat urmtoarele valori pentru
debitele de la priz :
pentru prenclzitorul 1- priza 1 :

=
=
032212 , 0
/ 826 , 8
51
51
a
s kg d

pentru prenclzitorul 2 - priza 2 :

=
=
047076 , 0
/ 899 , 12
52
52
a
s kg d

pentru prenclzitorul 3 - priza 3 :

=
=
05087 , 0
/ 938 , 13
53
53
a
s kg d

pentru prenclzitorul 4 - priza 4 :

=
=
034952 , 0
/ 577 , 9
54
54
a
s kg d

pentru prenclzitorul 5 - priza 5 :

=
=
042825 , 0
/ 734 , 11
55
55
a
s kg d

pentru prenclzitorul 6 - priza 6 :

=
=
046947 , 0
/ 864 , 12
56
56
a
s kg d

pentru prenclzitorul 7 - priza 7 :

=
=
064847 , 0
/ 768 , 17
57
57
a
s kg d

Calculul debitului necesar la turbopomp :
( )
( )
TPA mec
h h
h h d
d
, 19 18
32 33 32
18


= ;
unde : h
18
= h
12
= 2940,39 kJ/kg ,
p
18
= p
19
0.97 p
18
= 4,6204 0,97 p
18
= 4,4818 bar .
Din destinderea pe TPA am calculat
mec,TPA
= 0,8358 .
45

( ) [ ] [
( )] . / 4 , 2305 8358 , 0 / 2658 , 7 ; 0373 , 0
6 , 2793 6 , 2793 ;
19
19 18 19 18 18 19
kg kJ h K kg kJ bar h
h s p h h h h
TPA
=
= =


mec,TPA
= 0,98 (din ipoteze) .
nlocuind n relaia de mai sus vom avea :

( )
( )
s kg d d / 044 , 15
98 , 0 4 , 2305 39 , 2940
17 , 731 31 , 765 219 , 274
18 18
=


= .
Cum n calculul anterior am considerat randamente i debite prin turbin
calculate cu valoarea din d
18
, rezult c nu apar diferene care s justifice
recalcularea circuitului .

Recalcularea debitelor prin turbin, pe care le consider debite de referin:
d
2
= d
1
= 274 kg/s ,
d
3
= d
2
- d
1
b
1
d
3
= 274 - 274 0,0028 d
3
= 273,233 kg/s ,
d
4
= d
3
- d
1
b
2
d
4
= 273,233 - 274 0,0072 d
3
= 271,260 kg/s ,
d
5
= d
4
- d
1
b
3
d
3
= 271,260 - 274 0,0053 d
3
= 269,808 kg/s ,
d
6
= d
5
- d
57
d
3
= 269,808 - 17,768 d
6
= 252,040 kg/s ,
d
7
= d
6
+ d
1
a
injSI
d
7
= 252,04 - 274 0,0016 d
7
= 252,259 kg/s,
d
8
= d
7
= 252,259 kg/s ,
d
9
= d
8
- d
1
b
4
d
9
= 252,259 - 274 0,0016 d
9
= 251,82 kg/s ,
d
10
= d
9
- d
56
d
10
= 251,82 - 12,864 d
3
= 238,957 kg/s ,
d
11
= d
10
- d
55
d
11
= 238,957 - 11,734 d
3
= 227,223 kg/s ,
d
12
= d
11
- d
54
d
12
= 227,223 - 9,577 d
3
= 217,646 kg/s ,
d
13
= d
12
- d
1
b
5
- d
53
- d
18

d
13
= 217,646 - 274 0,0003-13,938- 15,044 d
13
= 188,581 kg/s ,
d
14
= d
13
= 188,581 kg/s ,
d
15
= d
14
- d
52
d
15
= 188,581 - 12,899 d
15
= 175,683 kg/s ,
d
16
= d
15
- d
51
d
16
= 175,683 - 8,826 d
16
= 166,857 kg/s ,
d
17
= d
16
+ d
1
b
6
d
17
= 166,857 + 274 0,00053 d
17
= 167,002 kg/s .
46

Prin urmare, cum debitele nu difer cu prea mult de cele anterioare,
consider acestea ca debite ale regimului de referin pentru calculele anterioare
i de asemenea volumele de referin sunt:
debitele de referin
d
r
: [274 ; 274 ; 273,233 ; 271,260 ; 269,808 ; 252,04 ; 252,259 ;
251,82 ; 238,95; 227,223 ; 217,646 ; 188,581 ; 188,581 ; 175,683 ;
166,857 ; 167,002 ; 15,044;14,52] .
volumele specifice de referint
v
r
: [1,694710
-2
; 1,745210
-2
; 1,781110
-2
; 2,079410
-2
; 5,302310
-2

; 5,450910
-2
;7,947510
-2
; 8,133810
-2
; 1,100410
-1
; 1,653010
-1
;
3,030210
-1
; 5,020710
-1
; 5,112910
-1
; 1,410710
0
; 4,856510
0
;
3,341310
1
; 3,409810
1
; 5,176210
1
; 3,341310
1
] .

10.5. Programul de calcul al blocului de 330 MW
% CAZAN DE 1035 t/h PE LIGNIT
etafcaz=0.91;Hatm=760; depF=12; taer0=40; etacaz=0.875;
paerF=(Hatm-depF)/760; Dmc0=4*100/3.6; g=9.80665; cpt=0.442;
VF=(43.9-11.8)*15.36*15.36+(11.8-1.2)*15.36^2-2*(11.8-1.2).^2/sqrt(3)*15.36;

% COMBustibil Solid, lignit de Oltenia-esantion mediu
xc=19.5; xh=1.76; xn=1.9; xo=8; xs=0.64; xa=25.2; xw=43;
Alfa0=1.25; xaer=10; tcomb=180; hcomb=1.13*(tcomb+273.15);
% Date la regimul de referinta

% Dm0, dinjip [kg/s] - debit masic abur viu
Dm0=274; ainjip1=0.0005; ainjip2=0.0005; Dms20=Dm0*(1-ainjip1);
Dmal0=(1-ainjip1-ainjip2)*Dm0;
% Dmcip0, ainjsi, Dmsi0 [kg/s] - debit masic abur CIP, injectie SI, dupa SI
47
Dmcip0=251.8; ainjsi=1e-3; Dmsi0=Dmcip0*(1+ainjsi);
% pab0, tab0 [bar] -presiune, temperatura abur viu = iesire S3
pab0=196; tab0=540;
% pabmp0, tabmp0, [bar] -presiune, temperatura iesire SI2 -spre CMP
pabmp0=45.8; tabmp0=540;

% DATE FUNCTIONALE
% Date functionale ECOnomizor
% p1eco0, t1eco0, [bar] -presiune, temperatura intrare ECO
p1eco0=246; t1eco0=250;
h1eco0=hfpt(p1eco0,t1eco0); [v1eco0,alf1eco,bet1eco]=vfpt(p1eco0,t1eco0);
% p2eco0, t2eco0 [bar] -presiune, temperatura iesire eco = intrare VAP
p2eco0=234; t2eco0=370; peco0=(p1eco0+p2eco0)/2;
h2eco0=hfpt(p2eco0,t2eco0); [v2eco0,alf2eco,bet2eco]=vfpt(p2eco0,t2eco0);
% Qeco0 [kW] -caldura transmisa supraincalzitorului eco in regim stationar
Qeco0=Dmal0*(h2eco0-h1eco0);
% tteco0 [C] - temperatura tevilor vaporizatorului in regim stationar
teco0=(t1eco0+t2eco0)/2; kglc=43.67; tteco0=Qeco0/kglc/Dmal0.^0.8+teco0;

% Date functionale VAPorizator
% p1v0, t1v0, p2v0, t2v0 [bar,C] -presiune, temp intrare, iesire VAP
p1v0=234; t1v0=370; [v1v0,alf1v,bet1v]=vfpt(p1v0,t1v0);
p2v0=224; t2v0=425; [v2v0,alf2v,bet2v]=vfpt(p2v0,t2v0);
% Qv0 [kW] -caldura transmisa vaporizatorului in regim stationar
h1v0=hfpt(p1v0,t1v0); h2v0=hfpt(p2v0,t2v0);
Qv0=Dmal0*(h2v0-h1v0);
% ttv0 [C] - temperatura tevilor vaporizatorului in regim stationar
tv0=(t1v0+t2v0)/2; kglv=0.4443; ttv0=(Qv0/kglv).^(1/3)+tv0;

% Date functionale S1
% p1s10, t1s10, [bar] -presiune, temperatura intrare S1 =iesire VAP
p1s10=224; t1s10=425;
h1s10=hfpt(p1s10,t1s10); [v1s10,alf1s1,bet1s1]=vfpt(p1s10,t1s10);
% p2s10, t2s10 [bar] -presiune, temperatura iesire S1 = intrare S2
48
p2s10=212; t2s10=450;
h2s10=hfpt(p2s10,t2s10); [v2s10,alf2s1,bet2s1]=vfpt(p2s10,t2s10);
% Qs10 [kW] -caldura transmisa supraincalzitorului S1 in regim stationar
Qs10=Dmal0*(h2s10-h1s10);
% tts10 [C] - temperatura tevilor vaporizatorului in regim stationar
ts10=(t1s10+t2s10)/2; kglc=43.67; tts10=Qs10/kglc/Dmal0.^0.8+ts10;

% Date functionale S2
% p1s20, t1s20, [bar] -presiune, temperatura intrare S2 =iesire S1
p1s20=212; t1s20=450;
%h1s20=hfpt(p1s20,t1s20);
h1s20=(h2s10*Dmal0+h1eco0*(Dms20-Dmal0))/Dms20;
t1s20=tfph(p1s20,h1s20); [v1s20,alf1s2,bet1s2]=vfpt(p1s20,t1s20);
% p2s20, t2s20 [bar] -presiune, temperatura iesire S2 = intrare S3
p2s20=204; t2s20=490;
h2s20=hfpt(p2s20,t2s20); [v2s20,alf2s2,bet2s2]=vfpt(p2s20,t2s20);
% Qs20 [kW] -caldura transmisa supraincalzitorului S2 in regim stationar
Dms20=(1-ainjip2)*Dm0; Qs20=Dms20*(h2s20-h1s20);
% tts20 [C] - temperatura tevilor vaporizatorului in regim stationar
ts20=(t1s20+t2s20)/2; kglc=43.67; tts20=Qs20/kglc/Dms20.^0.8+ts20;

% Date functionale S3
% p1s30, t1s30, [bar] -presiune, temperatura intrare S3 =iesire S2
p1s30=204; t1s30=490;
%h1s30=hfpt(p1s30,t1s30);
h1s30=(h2s20*Dms20+h1eco0*(Dm0-Dms20))/Dm0;
t1s30=tfph(p1s30,h1s30); [v1s30,alf1s3,bet1s3]=vfpt(p1s30,t1s30);
% p2s30, t2s30 [bar] -presiune, temperatura abur viu = iesire S3
p2s30=pab0; t2s30=tab0;
h2s30=hfpt(p2s30,t2s30); [v2s30,alf2s3,bet2s3]=vfpt(p2s30,t2s30);
% Qs30 [kW] -caldura transmisa supraincalzitorului S3 in regim stationar
Qs30=Dm0*(h2s30-h1s30);
% tts30 [C] - temperatura tevilor vaporizatorului in regim stationar
ts30=(t1s30+t2s30)/2; kglc=43.67; tts30=Qs30/kglc/Dm0.^0.8+ts30;
49

% Date functionale SI_1
% p1si10, t1si10, [bar, C] -presiune, temperatura intrare SI1=iesire CIP
p1si10=47.4; t1si10=348;
h1si10=hfpt(p1si10,t1si10); [v1si10,alf1si1,bet1si1]=vfpt(p1si10,t1si10);
% p2si10, t2si10, [bar] -presiune, temperatura iesire SI1
p2si10=46.3; t2si10=470;
h2si10=hfpt(p2si10,t2si10); [v2si10,alf2si1,bet2si1]=vfpt(p2si10,t2si10);
% Qsi10 [kW] -caldura transmisa SI_2 in regim stationar
Qsi10=Dmcip0*(h2si10-h1si10);
% ttsi10 [C] - temperatura tevilor vaporizatorului in regim stationar
tsi10=(t1si10+t2si10)/2; kglc=43.67; ttsi10=Qsi10/kglc/Dmcip0.^0.8+tsi10;
% Date functionale SI_2
% p1si20, t1si20, [bar, C] -presiune, temperatura intrare SI2=iesire SI1
p1si20=46.3; t1si20=470;
%h1si20=hfpt(p1si20,t1si20);
hinjmp=hlsatp(8); h1si20=(h2si10*Dmcip0+hinjmp*(Dmsi0-Dmcip0))/Dmsi0;
t1si20=tfph(p1si20,h1si20);[v1si20,alf1si2,bet1si2]=vfpt(p1si20,t1si20);
% p2si20, t2si20, [bar] -presiune, temperatura iesire SI2 -spre CMP
p2si20=pabmp0; t2si20=tabmp0;
h2si20=hfpt(p2si20,t2si20); [v2si20,alf2si2,bet2si2]=vfpt(p2si20,t2si20);
% Qsi20 [kW] -caldura transmisa SI_2 in regim stationar
Qsi20=Dmsi0*(h2si20-h1si20);
% ttsi20 [C] - temperatura tevilor vaporizatorului in regim stationar
tsi20=(t1si20+t2si20)/2; kglc=43.67; ttsi20=Qsi20/kglc/Dmsi0.^0.8+tsi20;

% Qcaz0 [kW] - caldura totala
Qcaz0=Qeco0+Qv0+Qs10+Qs20+Qs30+Qsi10+Qsi20;

% Date functionale PAR
% t2apar0, t1apar0 [C] - temperatura aerului la iesirea, intrarea din PAR
t2apar0=297.5; t1apar0=taer0; tapar0=(t1apar0+t2apar0)/2;
% h2apar0 [kJ/kg] -entalpia aerului umed preincalzit necesar arderii
haert0=(haert(0)+0.00161*xaer*hh2o(0))/(1+0.00161*xaer);
50
h2apar0=(haert(t2apar0)+0.00161*xaer*hh2o(t2apar0))/(1+0.00161*xaer)-haert0;
h1apar0=(haert(t1apar0)+0.00161*xaer*hh2o(t1apar0))/(1+0.00161*xaer)-haert0;
hapar0=(h1apar0+h2apar0)/2;
% t2gpar0, t1gpar0 [C] - temperatura gazelor la iesirea, intrarea din PAR
t2gpar0=140; t1gpar0=311;
% fsa, arec, fraer - fractie scapari, cota recirculare, fractie volum aer
fsa=0.1; arec=0.06; fraer=0.4;
% npar0 [rot/s] - turatie PAR
npar0=3/60; kglpar=20.867;
% p1apar0, p2apar0, p1gpar0, p2gpar0 [ Pa ] - presiuni
p1apar0=101325+456*g; p2apar0=101325+309*g;
p1gpar0=101325-351*g; p2gpar0=101325-207*g;


% DATE CONSTRUCTIVE
% Date constructive ECOnomizor
% Veco [m^3] -Volumul ECOnomizorului
Seco=20948; dieco=31.8e-3; deeco=40.8e-3; nteco=950;
Leco=Seco/nteco/pi/deeco; Veco=Seco*dieco^2/4/deeco; Veco=79;
%Aeco [m^2] - aria de curgere
Aeco=nteco*pi*dieco.^2/4;
% fecoe[1/m^4] -coeficient echivalent de pierderi
fecoe=1e+5*(p1eco0-p2eco0)/Dmal0^2/v2eco0;
% Mteco [kg] -masa tevilor ECOnomizorului
Mteco=Seco*(deeco^2-dieco^2)/4/deeco*7858;

% Date constructive VAPorizator
% Vv [m^3] -Volumul vaporizatorului
Vv=101;
% Sv [m^2] -suprafata de schimb de caldura vaporizator (MIT-I)
Sv=3825; dve=40.8e-3;
% dv [m] - diametrul interior tevi vaporizator
dv=31.8e-3;
% Lv [m] - lungimea tevilor vaporizator
51
Lv=43.9-11.235+(11.235-1.2)/sin(pi/3);
% fve[1/m^4] -coeficient echivalent de pierderi in vaporizator
fve=1e+5*(p1v0-p2v0)/Dmal0^2/v2v0;
%Av [m^2] - aria de curgere a vaporizatorului
Av=4*256*pi*dv.^2/4;
% ffv [-] - coeficient de frecare tevi vaporizator
anm=44.38*exp(-0.01688*p1v0); bnm=0.7878-0.4762e-3*p1v0; rugr=0.1e-3/dv;
rnm=1+anm*(1./2).^bnm; ffvmonof=1/4/(log10(3.72*dv/rugr))^2;
fcj=1.36+0.00725*p1v0+0.0738*(Dmal0/Av/1000)-0.0077*p1v0*(Dmal0/Av/1000);
ffv=ffvmonof*rnm*fcj;
% zetav [-] - coeficient pierderi locale tevi vaporizator
zetav=1;
% Mtv [kg] -masa tevilor vaporizatorului
Mtv=4*256*pi/4*(dve^2-dv^2)*Lv*7858;

% Date constructive S1
% Vs1 [m^3] -Volumul supraincalzitorului S1
Ss1=4734; Ls1=40;
Vs1=pi*26.9e-3^2/4*1364*11.9+pi*31.8e-3^2/4*684*28.1+ ...
pi*12/4*(175e-3^2+163e-3^2)*40+pi*44.8e-3^2/4*4*40+ ...
pi*12/4*88.9e-3^2*5*12; Vs1=16.3;
%As1 [m^2] - aria de curgere
As1=1364*pi*26.9e-3.^2/4;
% fs1e[1/m^4] -coeficient echivalent de pierderi
fs1e=1e+5*(p1s10-p2s10)/Dmal0^2/v2s10;
% Mts1 [kg] -masa tevilor S1
Mts1=(pi*(35.9e-3^2-26.9e-3^2)/4*1364*11.9+ ...
pi*(40.8e-3^2-31.8e-3^2)/4*684*28.1+ ...
pi*12/4*(235e-3^2-175e-3^2+211e-3^2-163e-3^2)*40+ ...
pi*(57.4e-3^2-44.8e-3^2)/4*4*40+ ...
pi*12/4*(106.5e-3^2-88.9e-3^2)*5*12)*7858;

% Date constructive S2
% Vs2 [m^3] -Volumul supraincalzitorului S2
52
Ss2=5669; dis2=31.8e-3; des2=40.8e-3; nts2=760; Ls2=Ss2/nts2/pi/des2;
Vs2=pi*dis2^2/4*nts2*Ls2; Vs2=19.5;
%As2 [m^2] - aria de curgere
As2=nts2*pi*dis2.^2/4;
% fs2e[1/m^4] -coeficient echivalent de pierderi
fs2e=1e+5*(p1s20-p2s20)/Dms20^2/v2s20;
% Mts2 [kg] -masa tevilor S2
Mts2=nts2*pi/4*(des2^2-dis2^2)*Ls2*7858;

% Date constructive S3
% Vs3 [m^3] -Volumul supraincalzitorului S3
Ss3=1827; dis3=31.8e-3; des3=44.4e-3; nts3=736; Ls3=Ss3/nts3/pi/des3;
Vs3=pi*dis3^2/4*nts3*Ls3; Vs3=6.2;
%As3 [m^2] - aria de curgere
As3=nts3*pi*dis3.^2/4;
% fs3e[1/m^4] -coeficient echivalent de pierderi
fs3e=1e+5*(p1s30-p2s30)/Dm0^2/v2s30;
% Mts3 [kg] -masa tevilor
Mts3=nts3*pi/4*(des3^2-dis3^2)*Ls3*7858;

% Date constructive SI_1
% Vsi1 [m^3] -Volumul supraincalzitorului intermediar SI_1
Ssi1=7889; disi1=51e-3; desi1=57.4e-3; ntsi1=950; Lsi1=Ssi1/ntsi1/pi/disi1;
Vsi1=pi*disi1^2/4*ntsi1*Lsi1; Vsi1=104;
%Asi1 [m^2] - aria de curgere
Asi1=ntsi1*pi*disi1.^2/4;
% fsi1e[1/m^4] -coeficient echivalent de pierderi
fsi1e=1e+5*(p1si10-p2si10)/Dmcip0^2/v2si10;
% Mtsi1 [kg] -masa tevilor SI_1
Mtsi1=ntsi1*pi/4*(desi1^2-disi1^2)*Lsi1*7858;

% Date constructive SI_2
% Vsi2 [m^3] -Volumul supraincalzitorului intermediar SI_2
Ssi2=2857; disi2=38e-3; desi2=44.4e-3; ntsi2=1128; Lsi2=Ssi2/ntsi2/pi/disi2;
53
Vsi2=pi*disi2^2/4*ntsi2*Lsi2; Vsi2=38;
%Asi2 [m^2] - aria de curgere
Asi2=ntsi2*pi*disi2.^2/4;
% fsi2e[1/m^4] -coeficient echivalent de pierderi
fsi2e=1e+5*(p1si20-p2si20)/Dmsi0^2/v2si20;
% Mtsi2 [kg] -masa tevilor SI_2
Mtsi2=ntsi2*pi/4*(desi2^2-disi2^2)*Lsi2*7858;

% Preincalzitor de Aer Regenerativ PAR
% Mrot, romt [kg] -masa, densitate metal in rotatie
Mrot=230000; romt=7800;
% dpar, hpar [m] - diametrul, inaltime PAR
dpar=10.608; hpar=2.3;
% Vlt [m^3] - volumul liber total
Vlt=pi*dpar.^2/4*hpar-Mrot/romt;

% Calcul ARDERE
% % hi [ kJ/kg ] - Puterea calorifica inferioara
hi=4.18.*(81.*xc+246.*(xh-xo./8)+25.*xs-6.*xw);
% Mtaer, Mtaerum [kg/kg] - masa teoretica aer uscat, respectiv umed
MO2=31.999; MC=12.01; MH2=2.016; MS=32.066;
Mtaer=(MO2/MC*xc+MO2/2/MH2*xh+MO2/MS*xs-xo)/100/0.2320;
Mtaerum=Mtaer*(1+0.00161*xaer);
% Mmaerum [kg/kmol] - masa molara aer umed
Maer=28.967; Raer=0.28703; MH2O=18.015; RH2O=0.46153; R=8.314473915;
Mmaerum=(Maer+0.00161*xaer*MH2O)/(1+0.00161*xaer); Raerum=R/Mmaerum;
% Mase molare compusi gazosi
% MgCO2 [kg/kg] - masa de CO2 produsa
MCO2=44.01; RCO2=0.188921; MgCO2=MCO2*xc/MC/100;
% MgSO2 [kg/kg] - masa de SO2 produsa
MSO2=64.063; RSO2 =0.129785967; MgSO2=MSO2/MS*xs/100;
% MgH2O [kg/kg] - masa de vapori de H2O produsa si intrata aer
MgH2O=(MH2O/MH2*xh+xw)/100+Alfa0*(Mtaerum-Mtaer);
% Mexaer [kg/kg] -masa de aer uscat in exces
54
Mexaer=(Alfa0-1)*Mtaer;
% MgN2 [kg/kg] - masa de N2 din gaze
MN2=28.013; RN2=0.2968; MgN2=0.768*Alfa0*Mtaer+xn/100;
% MgO2 [kg/kg] - masa de O2 din gaze
MO2=31.999; RO2=0.25984; MgO2=0.232*(Alfa0-1)*Mtaer;
% Mgaze [kg/kg] - masa gazelor fara continut de cenusa
Mgaze=MgCO2+MgSO2+MgH2O+MgN2+MgO2;
% Mmgaze [kg/kmol] - Masa molara aparenta a gazelor
Mmgaze=Mgaze/(MgCO2/MCO2+MgSO2/MSO2+MgH2O/MH2O+MgN2/MN2+
MgO2/MO2);
% Rgaze [kJ/(kg K)] - constanta de gaz perfect a gazelor produse
R=101325*22.414/273.15/1000; Rgaze=R/Mmgaze;
% Manorg [kg/kg] - masa de anorganic-cenusa
Manorg=xa/100;
% Mzgfocar [kg/kg] - masa de zgura ramasa in focar
xzgfocar=15/100; Mzgfocar=xzgfocar*Manorg; Mcenusa=Manorg-Mzgfocar;
% Compozitie cenusa Rovinari
MSiO2=60.089; MAl2O3= 75.2727; MFe2O3= 133.0327; MCaO=56.8;
MMgO=40.32;MSO3=80.058; xSiO2=0.487; xAl2O3=0.2801; xFe2O3=0.0981;
xCaO=0.0545; xMgO=0.0245; xSO3=0.0407;
Mmanorg=xSiO2*MSiO2+xAl2O3*MAl2O3+xFe2O3*MFe2O3+xCaO*MCaO+ ...
xMgO*MMgO+xSO3*MSO3;
% Mgaze [kg/kg] - Masa gazelor cu cenusa impure
Mgaze=MgCO2+MgSO2+MgH2O+MgN2+MgO2+Mcenusa;
% Mmgaze [kg/kmol] - Masa molara aparenta a gazelor continind si cenusa
Mmgaze=Mgaze/(MgCO2/MCO2+MgSO2/MSO2+MgH2O/MH2O+MgN2/MN2+ ...
MgO2/MO2+(Manorg-Mzgfocar)/Mmanorg);
% Rgaze [kJ/(kg K)] - constanta de gaz perfect a gazelor cu cenusa
Rgaze=R/Mmgaze;

% hgaze0 [kJ/kg] - entalpia initiala-adiabatica- a gazelor
hgaze0=etafcaz*(hcomb+hi+h2apar0*Mtaerum*Alfa0)/Mgaze;
tf1=1500; dtf=10; ind=1;
hgF0=(MgCO2*hco2(tf1)+MgSO2*hso2(tf1)+MgH2O*hh2o(tf1)+ ...
55
MgN2*hn2(tf1)+MgO2*ho2(tf1)+Mcenusa*hcenusa(tf1))/Mgaze;
while abs(hgaze0-hgF0) > 0.5
hgF0=(MgCO2*hco2(tf1)+MgSO2*hso2(tf1)+MgH2O*hh2o(tf1)+ ...
MgN2*hn2(tf1)+MgO2*ho2(tf1)+Mcenusa*hcenusa(tf1))/Mgaze;
if hgF0 < hgaze0;
if ind==1; dtf=dtf/2; end; ind=0; tf1=tf1+dtf;
else
if ind~=1; dtf=dtf/2; end; ind=1; tf1=tf1-dtf;end
end
% tgF0 [grd.C] -temperatura adiabatica a gazelor din focar
tgF0=tf1;
% hgF0 [kJ/kg] -entalpia gazelor la temperatura adiabatica
hgF0;
% rogF0 [kg/m^3] -densitatea gazelor la temperatura adiabatica
rogF0=paerF/1.01325/22.414*273.15/(273.15+tgF0)*Mmgaze;

% Entalpii gaze si temperaturi pe traseul gazelor
% Iesire ECOnomiazor estimat
tgeeco0=t1gpar0;
hgeeco0=(MgCO2*hco2(tgeeco0)+MgSO2*hso2(tgeeco0)+MgH2O*hh2o(tgeeco0)+ ...
MgN2*hn2(tgeeco0)+MgO2*ho2(tgeeco0)+Mcenusa*hcenusa(tgeeco0))/Mgaze;
% Dmc0 [kg/s] - Debit Masic de Combustibil necesar
Qpierd=Qcaz0*(1-etacaz/etafcaz);
Dmc0=(Qcaz0+Qpierd)/Mgaze/(hgF0-hgeeco0);
% Dmgaze [kg/s] - debit masic de gaze, in regim stationar
Dmgaze0=Dmc0*Mgaze;

% Iesire Vaporizator
hg2v0=hgF0-Qv0/Dmc0/Mgaze-Qpierd/Dmc0/Mgaze; tf=1100; dtf=10; ind=1;
hgaz=(MgCO2*hco2(tf)+MgSO2*hso2(tf)+MgH2O*hh2o(tf)+ ...
MgN2*hn2(tf)+MgO2*ho2(tf)+Mcenusa*hcenusa(tf))/Mgaze;
while abs(hgaz-hg2v0) > 0.5
hgaz=(MgCO2*hco2(tf)+MgSO2*hso2(tf)+MgH2O*hh2o(tf)+ ...
MgN2*hn2(tf)+MgO2*ho2(tf)+Mcenusa*hcenusa(tf))/Mgaze;
56
if hgaz < hg2v0;
if ind==1; dtf=dtf/2; end; ind=0; tf=tf+dtf;
else
if ind~=1; dtf=dtf/2; end; ind=1; tf=tf-dtf;end; end; tg2v0=tf;
% tgv0 [C] - temperatura medie de calcul a gazelor
tgv0=(tgF0+tg2v0)/2;
hgv0=(MgCO2*hco2(tgv0)+MgSO2*hso2(tgv0)+MgH2O*hh2o(tgv0)+ ...
MgN2*hn2(tgv0)+MgO2*ho2(tgv0)+Mcenusa*hcenusa(tgv0))/Mgaze;

% Iesire Supraincalzitor primar nr.3, S3
hg2s30=hg2v0-Qs30/Dmc0/Mgaze-(33.4^2-29^2)/33.4^2*Qs10/Dmc0/Mgaze;
tf=1000; dtf=10; ind=1;
hgaz=(MgCO2*hco2(tf)+MgSO2*hso2(tf)+MgH2O*hh2o(tf)+ ...
MgN2*hn2(tf)+MgO2*ho2(tf)+Mcenusa*hcenusa(tf))/Mgaze;
while abs(hgaz-hg2s30) > 0.5
hgaz=(MgCO2*hco2(tf)+MgSO2*hso2(tf)+MgH2O*hh2o(tf)+ ...
MgN2*hn2(tf)+MgO2*ho2(tf)+Mcenusa*hcenusa(tf))/Mgaze;
if hgaz < hg2s30;
if ind==1; dtf=dtf/2; end; ind=0; tf=tf+dtf;
else
if ind~=1; dtf=dtf/2; end; ind=1; tf=tf-dtf;end; end; tg2s30=tf;
% tgs30 [C] - temperatura medie de calcul a gazelor
tgs30=(tg2v0+tg2s30)/2;
hgs30=(MgCO2*hco2(tgs30)+MgSO2*hso2(tgs30)+MgH2O*hh2o(tgs30)+ ...
MgN2*hn2(tgs30)+MgO2*ho2(tgs30)+Mcenusa*hcenusa(tgs30))/Mgaze;

% Iesire Supraincalzitor Intermediar nr.2, SI_2
hg2si20=hg2s30-Qsi20/Dmc0/Mgaze-(29^2-25^2)/33.4^2*Qs10/Dmc0/Mgaze;
tf=1000; dtf=10; ind=1;
hgaz=(MgCO2*hco2(tf)+MgSO2*hso2(tf)+MgH2O*hh2o(tf)+ ...
MgN2*hn2(tf)+MgO2*ho2(tf)+Mcenusa*hcenusa(tf))/Mgaze;
while abs(hgaz-hg2si20) > 0.5
hgaz=(MgCO2*hco2(tf)+MgSO2*hso2(tf)+MgH2O*hh2o(tf)+ ...
MgN2*hn2(tf)+MgO2*ho2(tf)+Mcenusa*hcenusa(tf))/Mgaze;
57
if hgaz < hg2si20;
if ind==1; dtf=dtf/2; end; ind=0; tf=tf+dtf;
else
if ind~=1; dtf=dtf/2; end; ind=1; tf=tf-dtf;end; end; tg2si20=tf;
% tgsi20 [C] - temperatura medie de calcul a gazelor
tgsi20=(tg2s30+tg2si20)/2;
hgsi20=(MgCO2*hco2(tgsi20)+MgSO2*hso2(tgsi20)+MgH2O*hh2o(tgsi20)+ ...
MgN2*hn2(tgsi20)+MgO2*ho2(tgsi20)+Mcenusa*hcenusa(tgsi20))/Mgaze;

% Iesire Supraincalzitor primar nr.2, S2
hg2s20=hg2si20-Qs20/Dmc0/Mgaze-(25^2-20^2)/33.4^2*Qs10/Dmc0/Mgaze;
tf=1000; dtf=10; ind=1;
hgaz=(MgCO2*hco2(tf)+MgSO2*hso2(tf)+MgH2O*hh2o(tf)+ ...
MgN2*hn2(tf)+MgO2*ho2(tf)+Mcenusa*hcenusa(tf))/Mgaze;
while abs(hgaz-hg2s20) > 0.5
hgaz=(MgCO2*hco2(tf)+MgSO2*hso2(tf)+MgH2O*hh2o(tf)+ ...
MgN2*hn2(tf)+MgO2*ho2(tf)+Mcenusa*hcenusa(tf))/Mgaze;
if hgaz < hg2s20;
if ind==1; dtf=dtf/2; end; ind=0; tf=tf+dtf;
else
if ind~=1; dtf=dtf/2; end; ind=1; tf=tf-dtf;end; end; tg2s20=tf;
% tgs20 [C] - temperatura medie de calcul a gazelor
tgs20=(tg2si20+tg2s20)/2;
hgs20=(MgCO2*hco2(tgs20)+MgSO2*hso2(tgs20)+MgH2O*hh2o(tgs20)+ ...
MgN2*hn2(tgs20)+MgO2*ho2(tgs20)+Mcenusa*hcenusa(tgs20))/Mgaze;

% Iesire Supraincalzitor Intermediar nr.1, SI_1
hg2si10=hg2s20-Qsi10/Dmc0/Mgaze-(20^2-10^2)/33.4^2*Qs10/Dmc0/Mgaze;
tf=1000; dtf=10; ind=1;
hgaz=(MgCO2*hco2(tf)+MgSO2*hso2(tf)+MgH2O*hh2o(tf)+ ...
MgN2*hn2(tf)+MgO2*ho2(tf)+Mcenusa*hcenusa(tf))/Mgaze;
while abs(hgaz-hg2si10) > 0.5
hgaz=(MgCO2*hco2(tf)+MgSO2*hso2(tf)+MgH2O*hh2o(tf)+ ...
MgN2*hn2(tf)+MgO2*ho2(tf)+Mcenusa*hcenusa(tf))/Mgaze;
58
if hgaz < hg2si10;
if ind==1; dtf=dtf/2; end; ind=0; tf=tf+dtf;
else
if ind~=1; dtf=dtf/2; end; ind=1; tf=tf-dtf;end; end; tg2si10=tf;
% tgsi10 [C] - temperatura medie de calcul a gazelor
tgsi10=(tg2s20+tg2si10)/2;
hgsi10=(MgCO2*hco2(tgsi10)+MgSO2*hso2(tgsi10)+MgH2O*hh2o(tgsi10)+ ...
MgN2*hn2(tgsi10)+MgO2*ho2(tgsi10)+Mcenusa*hcenusa(tgsi10))/Mgaze;

% Iesire ECOnomizor
ces1=1-(20^2-10^2)/33.4^2-(25^2-20^2)/33.4^2-(29^2-25^2)/33.4^2-(33.4^2-29^2)/33.4^2;
hg2eco0=hg2si10-Qeco0/Dmc0/Mgaze-ces1*Qs10/Dmc0/Mgaze;
tf=1000; dtf=10; ind=1;
hgaz=(MgCO2*hco2(tf)+MgSO2*hso2(tf)+MgH2O*hh2o(tf)+ ...
MgN2*hn2(tf)+MgO2*ho2(tf)+Mcenusa*hcenusa(tf))/Mgaze;
while abs(hgaz-hg2eco0) > 0.5
hgaz=(MgCO2*hco2(tf)+MgSO2*hso2(tf)+MgH2O*hh2o(tf)+ ...
MgN2*hn2(tf)+MgO2*ho2(tf)+Mcenusa*hcenusa(tf))/Mgaze;
if hgaz < hg2eco0;
if ind==1; dtf=dtf/2; end; ind=0; tf=tf+dtf;
else
if ind~=1; dtf=dtf/2; end; ind=1; tf=tf-dtf;end; end; tg2eco0=tf;
% tgeco0 [C] - temperatura medie de calcul a gazelor
tgeco0=(tg2si10+tg2eco0)/2;
hgeco0=(MgCO2*hco2(tgeco0)+MgSO2*hso2(tgeco0)+MgH2O*hh2o(tgeco0)+ ...
MgN2*hn2(tgeco0)+MgO2*ho2(tgeco0)+Mcenusa*hcenusa(tgeco0))/Mgaze;
% Supraincalzitor primar nr.1, S1
tgs10=(tg2v0+tg2eco0)/2; hgs10=(hg2v0+hg2eco0)/2;
% Preincalzitor de Aer Regenerativ PAR
h2gpar0=(MgCO2*hco2(t2gpar0)+MgSO2*hso2(t2gpar0)+MgH2O*hh2o(t2gpar0)+ ...
MgN2*hn2(t2gpar0)+MgO2*ho2(t2gpar0)+Mcenusa*hcenusa(t2gpar0))/Mgaze;
t1gpar0=tg2eco0;
h1gpar0=(MgCO2*hco2(t1gpar0)+MgSO2*hso2(t1gpar0)+MgH2O*hh2o(t1gpar0)+ ...
MgN2*hn2(t1gpar0)+MgO2*ho2(t1gpar0)+Mcenusa*hcenusa(t1gpar0))/Mgaze;
59
tgpar0=(t1gpar0+t2gpar0)/2; hgpar0=(h1gpar0+h2gpar0)/2;
10.6. Rezultatele simulrii funcionrii grupului de 330 MW
Schema termic de calcul este cea prezentat n figura 10.1, iar datele
constructive i parametrii principali sunt dup antreprenor. Aceast schem nu ia
n consideraie pierderile de abur la labirini i tijele ventilului, apa de injecie
pentru reglarea temperaturii aburului viu i supranclzit intermediar, alte
recuperri de cldur.
Fa de aceasta, calculele din lucrare s-au efectuat n urmtoarele ipoteze:
pierderi de presiune pe conducte i n prenclzitoare constante;
diferenele finale de temperatur pe prenclzitoare constante;
repomparea condensatului la PJP 2;
randamentele pompelor constante;
se recupereaz n circuit circa 85 % din pierderile la generator
(16,9381 kJ/kg) i mecanice (4,067 kJ/kg);
cota procentual de recuperri de abur a fost raportat la debitul de
abur viu i considerat la presiunea prizei respective, fiind corelat i cu
datele grupului de 800 MW.
Comparnd valorile date de furnizor cu cele ale regimului de referin
calculat, pentru a obine puterea electric de 330 MW, din tabelul 1 se observ
c eroarea relativ este acceptabil (<2 %), regimul constituind referin pentru
calculul altor regimuri staionare nenominale.
Tab. 1. Comparaia valorilor calculate cu cele date de furnizor
Nr.
punct
Debit
iniial
[ kg/s ]
Debit
calculat
[ kg/s ]
Eroare
debit
[ % ]
Presiune
iniial
[ bar ]
Presiune
calculat
[ bar ]
Eroare
presiune
[ %]
1. 274.0000 274.0000 0 192.0000 192.0000 0
2. 274.0000 274.0000 0 186.1056 186.1056 0
3. 273.3330 273.2328 0.0368 182.4021 182.4021 0
60
4. 273.3330 271.2600 0.7584 151.9553 151.9553 0
5. 273.3330 269.8078 1.2897 48.8371 48.9787 -0.2899
6. 255.833 252.0875 0.6890 47.5899 47.5899 0
7. 255.833 252.3067 0.6890 46.1023 46.1023 0
8. 255.9722 252.3067 1.4320 44.9145 45.0089 -0.2102
9. 255.8333 251.8683 1.5498 30.8909 31.1153 -0.7264
10. 241.9722 239.0556 1.2201 18.6326 18.7857 -0.8217
11. 229.8889 227.2969 1.1275 8.745 8.8306 -0.9788
12. 219.8333 217.7151 0.9635 4.709 4.7207 -0.2463
13. 188.6667 188.6062 0.0321 4.6091 4.6287 -0.4252
14. 188.6667 188.6062 0.0321 1.2503 1.2566 -0.5039
15. 175.6111 175.6936 -0.0470 0.3130 0.3130 0
16. 164.4444 165.8597 -0.8607 0.0373 0.0373 0
10.7. ncadrarea programului turbinei de 330 MW n modulele
cazanului
n lucrarea de fat se elaboreaz modelele matematice ale proceselor
staionare i dinamice care urmeaz s fie implementate pe un simulator n timp
real, destinat personalului de exploatare.
Cum cazanul de abur cu strbatere forat nu poate funciona fr turbin i
serviciile interne care fac parte din blocul de 330 MW, suplimentar temei
iniiale, s-a preconizat i elaborarea unor modele staionare i dinamice ale
turbinei, care s determine valorile parametrilor de schimb ntre turnin i cazan.
n figura 10.2 se prezint schema de ncadrare a turbinei n blocul energetic
de 330 MW, cu mrimile care sunt necesare funcionrii cazanului modelat.
Aceste mrimi sunt :
debitul de abur primar preluat de turbin la ncrcarea / descrcarea
grupului; parametrii temperatur i presiune ai acestui debit de abur sunt
determinai n modulul cazan;
debitul de abur la ieirea din corpul de nalt presiune al turbinei, DSI,
cu presiunea PSI i temperatura TSI; valori calculate n modelul turbinei;
61
debitul de abur trimis de la supranclzitorul intermediar 2 al cazanului
la corpul de medie presiune al turbinei, DSE, cu presiunea TSE (valori
calculate n modulul cazanului);
temperatura apei de alimentare cazan, Tal i debitul apai de alimentare,
Wal, calculate n modulul turbinei.
Pe baza analizei fcute anterior a rezultat c fiecare din marimile de schimb
menionate mai sus i care aparin modulului turbinei se pot calcula pentru
regimul staionar, cu mare precizie, cu ajutorul unor relaii polinomiale.
Notnd cu y mrimea ce dorim s fie calculat, cu c
1
... c
5
coeficienii i cu
D debitul de abur absorbit de corpul de nalt presiune al turbinei, se obine:
5 4 3 2 1
2 3 4
c D c D c D c D c y + + + + =
Menionm c, aceast relaie permite s se calculeze debitele de abur
schimbate ntre supranclzirea intermediar cazan-turbin, presiuni i
temperaturi abur i ap, putere electric, putere termic, etc.
Constantele c
1
... c
5
difer pentru fiecare mrime calculat i sunt redate n
tabelul 10.1 i tabelul 10.2.