|
CRESTEREA EFICIENTEI SI REDUCEREA POLUARII IN TRANSPORTURI PRIN UTILIZAREA GRUPURILOR COGENERATIVE BINARE - EXPERIMENTE
1. Obiectivele cercetarii experimentale
Analizele preliminare realizate indica faptul ca radiatorul instalatiei de condensare este componenta grupului binar avand gabaritul si greutatea cele mai mari, suprafata de schimb de caldura a acestuia fiind de cel putin doua ori mai mare decat suprafata totala a schimbatoarelor de caldura ale CRGA. Sporul de greutate si gabarit datorat instalatiei de condensare alaturi de pierderile de apa pe care le presupune functionarea acesteia constituie principalele impedimente in procesul de adaptare a grupurilor binare ca sisteme de propulsie terestra. Din aceste considerente, instalatia de condensare constituie componenta grupului binar analizata in cadrul experimentelor.
In cadrul cercetarii experimentale, care a vizat instalatiile de condensare cu schimbator de caldura de suprafata (schema Heller), au fost urmarite urmatoarele obiective:
intensificarea schimbului de caldura (cresterea coeficientului global de schimb de caldura) in radiator in scopul reducerii suprafetei de schimb de caldura a acestuia; metodele experimentate au fost urmatoarele:
cresterea vitezei de circulatie a aerului prin radiator;
umidificarea aerului prin pulverizare de apa in curentul de aer, inaintea radiatorului;
stabilirea unei variante constructive care sa permita recuperarea, cu pierderi minime, a debitului de apa antrenat de curentul de aer (in cazul sistemului de racire cu picaturi);
evaluarea debitului de apa de pierderi atat in cazul sistemului de racire cu picaturi;
evaluarea consumului de energie necesar antrenarii aerului de catre ventilator.
2. Instalatia experimentala
In cadrul cercetarii experimentale s-a urmarit studierea procesului de cedare in mediul ambiant a caldurii latente pe care apa din instalatia de condensare a grupului binar o preia, in procesul de condensare, de la aburul evacuat din turbina. Realizarea obiectivelor prezentate in subcap. 1. a impus desfasurarea experimentelor in doua etape, fiecareia dintre acestea corespunzandu-i o varianta specifica de instalatie experimentala. In ambele variante, curentul de aer necesar racirii apei a fost furnizat prin intermediul unui tunel aerodinamic prevazut cu ventilator de 17 kW, la care au fost conectate instalatiile experimentale caracteristice etapelor. Variatia debitului/vitezei aerului se poate obtine prin actionarea unui siber.
Cele doua etape de experimentare si variantele corespunzatoare de instalatie sunt prezentate in continuare.
2.1. Instalatie pentru analiza schimbului de caldura intensificat
prin cresterea vitezei de circulatie a aerului
In aceasta etapa experimentala s-a urmarit analizarea influentei vitezei de circulatie a aerului asupra coeficientului global de schimb de caldura al unui radiator si asupra consumului de energie datorat caderii de presiune pe acesta. Varianta constructiva de instalatie care a permis desfasurarea experimentelor in aceasta etapa a cercetarii este prezentata in fig. 1, din care lipsesc dispozitivul de pulverizare a apei (DP), dispozitivul de masurare a pierderilor de apa (DMCA) si tubul Anubar (TA), conectat cu tubul manometric diferential TU2.
Temperaturile aerului inainte si dupa radiator sunt masurate cu ajutorul termocuplurilor TC1, respectiv TC2. Tubul manometric diferential TU1 permite determinarea caderii de presiune a aerului pe radiator. Debitmetrul volumetric D permite determinarea debitului de apa. Incalzirea apei care urmeaza a fi racita in radiatorul Rad se realizeaza in vasul VI, cu ajutorul a trei termoplonjoare (TP1, TP2, TP3) fiecare avand puterea termica de 1 kW. Circulatia apei este asigurata de pompa P, alimentata de la o sursa de curent continuu. Temperaturile apei la intrarea/iesirea in/din radiator se determina cu ajutorul termocuplurilor TC4, tespectiv TC1.
Suprafata totala de schimb de caldura a radiatorului este iar sectiunea de trecere a acestuia este .
2.2. Instalatie pentru analiza schimbului de caldura in cazul
umidificarii aerului prin pulverizare de apa in curentul de aer
Intensificarea schimbului de caldura prin umidificarea aerului nu constituie o noutate, activitati de cercetare in aceasta privinta fiind desfasurate, in cadrul catedrei de Masini Termice si Termotehnica a Universitatii Tehnice "Gh. Asachi" din Iasi, inca dinainte de anul 1980 [29]. Studiul actual reprezinta o continuare a cercetarilor desfasurate anterior in catedra respectiva.
Schema instalatiei corespunzatoare acestei etape este prezentata in fig. 1. In fig. 2 este prezentata o fotografie a instalatiei. Componentele caracteristice acestei variante constructive sunt dispozitivul de pulverizare a apei (DP), dispozitivul de masurare a pierderilor de apa (DMCA) si tubul Anubar (TA), conectat cu tubul manometric diferential TU2.
Ansamblul TA - TU2 permite determinarea vitezei si, ulterior, a debitului de aer. Pentru citirea mai precisa a indicatiilor tubului manometric TU2, acesta este inclinat sub un unghi α = 8˚ fata de planul orizontal.
Fig.1
Fig. 2.
Dispozitivul de masurare a consumului de apa (a debitului de apa pulverizat de DP) este compus din vasele V1 (gradat), V2 si V3. La inceputul fiecarei determinari experimentale se noteaza nivelul apei in vasul V1 (hi). Cu ajutorul robinetului r1 se regleaza debitul de apa introdus in V2 astfel incat acesta sa fie cel putin egal cu debitul de apa pulverizata de DP. Ceea ce se urmareste practic este mentinerea nivelului apei in V2 la valoarea maxima (pentru mentinerea constanta a presiunii apei la aspiratie), reducand la minim cantitatea de apa scursa din V2 in V3. Aceasta minimizare este necesara pentru
cresterea la maxim a intervalului de timp necesar umplerii lui V2, deci a intervalului de timp in care se poate determina cantitatea de apa pulverizata, ceea ce echivaleaza cu obtinerea unei precizii maxime de determinare a debitului de apa pulverizata. La finalul fiecarei masuratori, apa scursa din V2 in V3 se recupereaza prin deschiderea robinetului r2 si se reintroduce in V1. Rezulta, astfel, nivelul final al apei in V1 (hf). Debitul volumic de apa aspirata se determina in functie de diferenta si de intervalul de timp τp (cronometrat) in care nivelul apei din vasul V1 variaza de la hi la hf.
3. Etalonari
Instrumentele de masura care au necesitat etalonare sunt termocuplurile si tubul Anubar. Pentru stabilirera procedurilor de etalonare a acestor doua dispozitive, care sunt prezentate in continuare, au fost utile informatiile din [76] si [39].
3.1. Etalonarea termocuplurilor
Termocuplurile utilizate sunt de tipul cromel-alumel. In cadrul procedurii de etalonare s-a folosit drept etalon un termometru de precizie ridicata, avand valoarea unei diviziuni de 0,1˚. Cele patru termocupluri si termometrul etalon s-au montat in capacul suport al unui vas plin cu apa (VE) astfel incat rezervorul termometrului si jonctiunile de lucru ale termocuplurilor sa fie foarte apropiate, evitand insa atingerea lor (a se vedea fig. 3). Aceasta distanta minima se impune pentru a se asigura preluarea informatiilor privind temperatura dintr-o zona cu gradienti termici cat mai mici - conditie absolut necesara pentru o etalonare corecta.
Fig. 3.
Cablurile de prelungire (compensatie) ale termocuplurior sunt conectate la aparatul de masura AMT, prevazut cu termostat. Acesta mentine constanta temperatura jonctiunilor de referinta ale termocuplurilor, eliminand erorile datorate variatiei temperaturii mediului ambiant. Pentru a verifica eficienta aparatului de masura termostatat, etalonarea s-a facut de doua ori, la doua temperaturi diferite ale mediului ambiant. Citirile au fost efectuate atat la cresterea cat si la descresterea temperaturii. Astfel, au fost trasate in total cate patru curbe de etalonare pentru fiecare termocuplu.
Incalzirea apei s-a realizat utilizand un termoplonjor (TP) cu puterea termica de 500 W. Pentru trasarea curbelor de etalonare la descresterea temperaturii, racirea s-a produs liber pana la aproximativ 30˚C. Temperaturile sub aceasta valoare au fost obtinute prin introducere de gheata in vasul VE.
Tabelele 1 si 2 prezinta citirile corespunzatoare celor doua temperaturi ale mediului ambiant. In ambele tabele, in prima coloana sunt date valorile indicate de termometrul etalon, iar coloanele 2 - 5, respectiv 6 - 9 prezinta indicatiile date de termocuplurile TC1.TC4 la cresterea, respectiv descresterea temperaturii.
Tabelul 1. Citiri efectuate la temperatura mediului ambiant
Nr.
crt.
tTE
Cresterea temperaturii
Descresterea temperaturii
TC1
TC2
TC3
TC4
TC1
TC2
TC3
TC4
1.
5
5,9
5,7
5,7
5,8
7
6,6
6,6
6,8
2.
10
10,7
10,5
10,6
10,7
11,9
11,6
11,6
11,8
3.
15
15,6
15,4
15,5
15,6
17,2
16,8
16,7
16,9
20
20,6
20,5
20,5
20,6
22,2
21,8
21,7
21,9
5.
25
25,8
25,6
25,6
25,7
26,6
26,4
26,4
26,7
6.
30
30,8
30,6
30,7
30,8
31,8
31,5
31,5
31,8
7.
35
35,8
35,6
35,7
35,8
36,9
36,6
36,6
36,9
8.
40
40,8
40,6
40,7
40,8
41,9
41,6
41,6
41,9
9.
45
46,1
46,1
46,1
46,2
46,9
46,7
46,8
47,3
10.
50
51
51
51
51,2
51,5
51,3
51,4
51,7
11.
55
56,1
56,1
56,1
56,2
56,5
56,3
56,3
56,6
12.
60
60,9
60,9
61
61,2
61,6
61,4
61,4
61,7
13.
65
65,9
65,9
66,1
66,3
66,9
66,8
66,9
67,2
1
70
70,9
70,8
70,9
71,2
71,3
71,2
71,3
71,6
15.
75
76,2
76
76
76,3
76,2
76,1
76,3
76,6
16.
80
81
81
81,1
81,6
81,2
81,2
81,2
81,7
17.
85
86,1
86,2
86,4
86,7
85,6
85,7
85,7
86,3
18.
90
91,4
91,2
91,4
91,7
91,4
91,2
91,4
91,7
Tabelul 2. Citiri efectuate la temperatura mediului ambiant
Nr.
crt.
tTE
Cresterea temperaturii
Descresterea temperaturii
TC1
TC2
TC3
TC4
TC1
TC2
TC3
TC4
1.
5
7,5
7,2
7,2
7,4
5,8
5,6
5,8
6,1
2.
10
10,6
10,4
10,4
10,6
10,8
10,6
10,7
11
3.
15
15,6
15,4
15,4
15,6
15,8
15,6
15,6
16
20
20,7
20,4
20,4
20,7
20,9
20,7
20,7
21
5.
25
25,6
25,4
25,4
25,7
26
25,8
25,8
26,1
6.
30
30,7
30,5
30,5
30,7
31
30,8
30,8
31,1
7.
35
35,8
35,6
35,6
35,8
36
35,8
35,8
36,1
8.
40
40,8
40,6
40,6
40,8
41,1
40,9
40,9
41,2
9.
45
46,1
46
46
46,3
46,1
46,1
45,9
46,2
10.
50
51,2
51
51
51,2
51
51
51
51,2
11.
55
56,4
56,1
56,1
56,3
56,1
56
56
56,4
12.
60
61,5
61,3
61,2
61,5
61,3
61,2
61,2
61,5
13.
65
66,6
66,4
66,4
66,7
66,3
66,2
66,4
66,7
1
70
71,8
71,5
71,5
71,8
71,5
71,4
71,5
71,8
15.
75
77
76,8
76,8
77
76,7
76,6
76,7
77,1
16.
80
82,1
81,8
81,8
82,1
81,4
81,8
81,8
81,6
17.
85
86,8
86,7
86,7
87
86,7
86,4
87,5
87,9
18.
90
92,5
92,2
92,1
92,5
91,6
91,7
91,7
91,9
Curbele de etalonarea a termocuplurilor, corespunzatoare valorilor din tabelele 1 si 2, sunt date in fig. Se observa ca toate cele 16 curbe practic se suprapun, ceea ce confirma eficienta termostatului aparatului de masura si permite stabilirea unei ecuatii comune de calcul a temperaturilor pentru cele patru termocupluri. Pentru stabilirea ecuatiei s-au luat in considerare valorile medii aritmetice (TCm) ale indicatiilor celor patru termocupluri. Aceste valori si curba generala de etalonare sunt prezentate in fig. 5.
Formula dedusa pentru calculul a temperaturii este urmatoarea:
. (1)
Fig.
Fig. 5.
Curba generata de valorile calculate ale temperaturii este cea reprezentata cu linie punctata in fig. 5. Relatia (1) aproximeaza curba tTE = f(TCm) cu o eroare maxima de 04%, atinsa atunci cand termometrul etalon indica 5˚C (valoarea minima masurata in procesul de etalonare). De altfel, precizia mai scazuta in determinarea temperaturii este o caracteristica a temperaturilor joase. Aceasta se datoreaza faptului ca, pe masura ce temperatura scade spre 0˚C, o zecime de grad reprezinta un procent tot mai important din valoarea temperaturii, deci o eroare de o zecime de grad determina o eroare procentuala tot mai mare. Pentru experimentele efectuate, domeniul de interes este cel al temperaturilor de peste 20˚C. Pentru aceste valori, eroarea maxima de aproximare a curbei tTE = f(TCm) este de 1,3%; intre valorile temperaturii determinate cu ajutorul termocuplelor si valorile indicate de termometrul etalon este o diferenta maxima de 4,8%.
3.2. Etalonarea tubului Anubar
Tubul Anubar (TA) este alcatuit din doua tuburi, dintre care unul prevazut cu patru prize de presiune totala si montat perpendicular pe directia de curgere a aerului iar celalalt este dispus paralel cu directia curentului, cu orificiul spre aval (reprezentand priza de presiune statica).Viteza de curgere a aerului se resimte printr-o presiune diferentiala, citita la un tub manometric diferential (TU2). Pentru cresterea preciziei de masurare, TU2 a fost inclinat sub un unghi α = 8˚ fata de orizontala.
Aparatul de masura al vitezei aerului utilizat drept etalon este un anemometru cu palete axiale. Aparatul permite determinarea vitezei aerului cu precizie de pana la 1%.
Etalonarea tubului Anubar presupune stabilirea ecuatiilor matematice caracteristice curbelor de variatie a vitezelor de interes ale aerului in functie de diferenta de coloana de apa Δh2 (masurata cu TU2). Sunt doua viteze de interes si anume viteza in sectiunea tronsonului in care se afla plasat tubul Anubar (wTA) si viteza prin radiator (wm rad). Pentru etalonare se porneste ventilatorului instalatiei si se deschide treptat siberul; pentru fiecare pozitie a siberului se noteaza diferentele de coloana de apa Δh1 (masurata cu TU1) si Δh2, valoarea vitezei masurata de anemometru (wA) si valoarile indicate de TC2 si TC3. Pe tot parcursul etalonarii se verifica temperatura mediului ambiant (tamb) cu termometrul etalon TE, utilizat si la etalonarea termocuplurilor. Valorile masurate ale acestor parametrii sunt date in coloanele 1 - 6 ale tabelului 3.
Tabelul 3.
Nr.
crt.
Valori masurate
Valori calculate
Δh1
[mm H2O]
Δh2
[mm H2O]
wA
[m/s]
TC2
TC3
tamb
[˚C]
t2
[˚C]
tTA
[˚C]
wm TA
[m/s]
wc TA
[m/s]
wm rad
[m/s]
wc rad
[m/s]
1.
20,5
4
11,32
25,6
25,7
24
24,45
24,46
4,23
4,3
15,43
15,69
2.
25
4,5
12,4
25,6
25,8
24
24,45
24,56
4,63
4,59
16,9
16,75
3.
36
7
15,5
25,5
25,7
24
24,35
24,46
5,78
5,65
21,1
20,63
41
9,5
17,2
25,6
25,7
24
24,45
24,46
6,41
6,44
23,4
23,52
5.
48,5
12,5
19,27
25,6
25,8
24
24,45
24,56
7,18
7,26
26,2
26,52
6.
53,5
14,5
20,8
25,6
25,8
24
24,45
24,56
7,75
7,76
28,31
28,38
7.
59,5
17,5
22,7
25,6
25,8
24
24,45
24,56
8,45
8,49
30,89
31,06
8.
68
19,85
24,43
25,5
25,7
24
24,35
24,46
9,08
9,08
33,21
33,23
9.
79
23
26,69
25,6
25,8
24
24,45
24,56
9,91
9,78
36,28
35,78
10.
83
27
28,55
25,6
25,8
24
24,45
24,56
10,6
10,69
38,8
39,13
11.
89,5
30
30,73
25,6
25,7
24
24,45
24,46
11,4
11,37
41,73
41,62
Determinarea vitezelor wTA si presupune cunoasterea densitatilor aerului in sectiunea de trecere a anemometrului (ρA), in sectiunea tronsonului in care se afla plasat tubul Anubar (ρTA) si a densitatii medii a aerului in radiator (). Pentru a putea calcula aceste densitati este necesara cunoasterea presiunilor si a temperaturilor aerului inainte si dupa radiator (p3, t3, respectiv p2, t2) precum si in sectiunea in care se afla plasat tubul Anubar (pTA, tTA).
Pe traseul aerului, suprapresiunea maxima pe care o poate realiza ventilatorul (precizata de producator) este de 300 mm H2O (0,0294 bar). Aceasta suprapresiune, care reprezinta valoarea teoretica maxima posibila a caderii de presiune pe radiator, poate determina o variatie a densitatii de pana la 3%. In aceste conditii, pot fi facute urmatoarele estimari:
; ,
unde Δp1 reprezinta caderea de presiune pe radiator. Deoarece intre radiator si tubul Anubar nu exista rezistente locale iar pierderile liniare de presiune sunt practic nule, se admite.
Valorile temperaturilor tTA si t2 se calculeaza cu relatia (1), in functie de indicatiile termocuplurilor TC3 (plasat in acest caz in amonte de tubul Anubar, intr-o sectiune apropiata) si TC2. Urmarind valorile temperaturilor tTA si t2, date in tabelul 3, se constata o variatie foarte mica (maxim 0,21˚C) a temperaturii aerului intre sectiunile in care se masoara cele doua temperaturi. In acest caz, se poate face aproximarea t3 = t2.
Se pot determina acum, din ecuatiile de stare, densitatile ρA, ρ3 si ρm aer:
, (2)
, (3)
respectiv
, (4)
in care:
pA - presiunea in sectiunea de trecere a anemometrului; pA = p0 = 1,013 bar;
Raer - constanta aerului; ;
Tamb - temperatura absoluta a mediului ambiant;
TTA - temperatura absoluta in sectiunea in care se afla plasat tubul Anubar;
pm rad - presiunea medie a aerului in radiator; ;
Tm rad - temperatura absoluta medie a aerului; ,
unde T3 si T2 sunt temperaturile absolute inainte/dupa radiator radiator.
Cunoscand sectiunea de trecere a anemometrului si sectiunea tronsonului in care este plasat tubul Anubar (311 mm 150 mm), din ecuatia de continuitate aplicata intre sectiunile SA si STA
, (5)
rezulta
. (6)
In mod similar, se determina:
. (7)
Valorile rezultate ale vitezei sunt date in coloanele cu numarul 9 ale tabelelor 3 si In coloana 11 a tabelului 3 sunt date valorile vitezei .
Fig. 6.
Curbele de etalonare si wm rad = f(Δh2) sunt prezentate in fig. 6. Ecuatiile de etalonare a tubului Anubar, stabilite pe baza acestor curbe, sunt urmatoarele:
, (8)
. (9)
Valorile vitezelor calculate cu aceste relatii sunt date in coloanele cu numarul 8 si 10 ale tabelului 3. Curbele generate de aceste valori - si - sunt prezentate de asemenea in fig. 6, fiind trasate cu linii punctate.
Se observa profilul liniar al curbelor atunci cand Δh2 > 10 mm H2O. De aceea, prin extrapolare, se poate accepta validitatea relatiilor (8) si (9) si pentru viteze de peste 30 m/s, dar care nu depasesc exagerat aceasta valoare.
In cadrul experimentelor, diferenta de temperatura (intre temperaturile in sectiunea in care intereseaza valoarea vitezei, respectiv in sectiunea in care se preia informatia privind valoarea vitezei) nu mai poate fi neglijata. In acest caz, viteza aerului in radiator se corecteaza cu coeficientul
, (10)
viteza finala fiind
(11)
Metoda de cercetare experimentala
1. Metoda caracteristica primei etape de cercetare experimentala
In aceasta etapa a cercetarii s-a analizat posibilitatea intensificarii transferului de caldura in cazul unui radiator prin cresterea vitezei aerului de racire. In acest scop, s-a urmarit trasarea unei curbe krad = f(wrad), unde krad reprezinta coeficientul global de schimb de caldura al radiatorului, wrad fiind viteza aerului prin radiator. Pe langa efectul dorit al imbunatatirii schimbului de caldura, cresterea vitezei aerului prezinta si dezavantajul cresterii consumului energetic al ventilatorul instalatiei datorat cresterii caderii de presiune pe radiator (erad). Pentru a analiza aceasta variatieeste necesara trasarea curbei erad = f(wrad).
Modificarea vitezei aerului prin radiator s-a facut prin actionarea siberului de pe aspiratia ventilatorului. Pentru fiecare regim de viteza, determinarea parametrilor krad si erad a presupus masurarea temperaturilor apei la intrarea/iesirea in/din radiator (t4, respectiv t1), a temperaturilor aerului inainte/dupa radiator (t3, respectiv t2) si a debitului de apa. Temperaturile t1.t4 au fost determinate cu ajutorul termocuplurilor TC1.TC4 (fig. 1), aplicand relatia (1). Debitul volumic de apa s-a determinat simplu, cu relatia
, (12)
in care Vi si Vf sunt valorile citite pe cadranul debitmetrului volumetric D la inceputul, respectiv sfarsitul fiecarei experiente, iar τ este timpul scurs intre cele doua citiri, exprimat in secunde.
In continuare, au fost determinate entalpiile apei si aerului in punctele de masurare, densitatile aerului in punctele de masurare si densitatile medii (corespunzatoare temperaturilor si presiunilor medii pe circuitele de apa si aer).
Procedura de calcul a densitatii medii (ρm apa) si a entalpiilor apei (i4 - intrarea in radiator, i1 - iesirea din radiator) este cea utilizata si in cazul subrutinei APA, descrisa in subcap. 3.3.2. Deoarece P este o pompa de circulatie (variatiile densitatii si entalpiei apei la variatii mici ale presiunii sunt nesemnificative - sub 0,05% atunci cand presiunea variaza de la 1 bar la 2 bar), in toate calculele efectuate s-au adoptat .
Asa cum s-a precizat in subcap. 3.2, variatia densitatii, pe traseul aerului, este de maxim 3%. Variatia entalpiei este de maximum 0,6 %. In aceste conditii, pentru determinarea entalpiilor aerului s-a utilizat formula
, (13)
din [62] in care valorile coeficientilor an sunt urmatoarele:
Aplicand ecuatia de stare in punctul 2, se determina densitatea aerului dupa radiator:
. (14)
Densitatea aerului inainte de radiator si densitatea medie a aerului se calculeaza cu relatiile (3) si (4). Se mentioneaza ca, de aceasta data, nu este valabila egalitatea .
In continuare, se determina cantitatea de caldura cedata si diferenta de temperatura medie logaritmica, relatiile de calcul fiind urmatoarele:
, (15)
respectiv
. (16)
Semnificatiile marimilor care intervin in aceste relatii de calcul sunt urmatoarele:
- debitul masic de apa prin radiator, in kg/s;
ρm apa - densitatea medie a apei, in kg/m3;
;
.
Se poate determina acum coeficientul global de schimb de caldura:
. (17)
Calcului vitezei aerului prin radiator presupune cunoasterea debitului masic de aer. Acesta se calculeaza cu formula
. (18)
Viteza aerului prin radiator rezulta din formula
. (19)
Puterea consumata de ventilator datorita radiatorului se calculeaza cu relatia
, (20)
in care ηv este randamentul ventilatorului; in cazul de fata, .
Consumul energetic datorat caderii de presiune pe radiator este dat de raportul
(21)
Cunoscand valorile parametrilor krad si erad corespunzatoare fiecarui regim de viteza, se pot trasa curbele propuse - krad = f(wrad), respectiv erad = f(wrad).
2. Metoda de cercetare experimentala caracteristica etapei a doua
In aceasta etapa s-au analizat influentele vitezei aerului prin radiator (wrad) si a concentratiei de apa pulverizata in aer, definita ca
[kg apa/kg aer uscat], (22)
asupra coeficientului global de schimb de caldura al radiatorului. Parametrii si din relatia (22) reprezina debitul de apa pulverizata, respectiv debitul de aer uscat, in kg/s.
Pentru aprecierea schimbului de caldura, s-a urmarit trasarea curbelor krad = f(wrad) la diferite valori ale lui ψ, respectiv krad = f(ψ) pentru valori diferite ale lui wrad. S-a avut in vedere, de asemenea, trasarea curbei erad = f(wrad).
Pentru un anumit debit de aer uscat, modificarea lui ψ se obtine prin modificarea nivelului h0 al apei aspirate de dispozitivul de pulverizare DP (a se vedea fig. 1).
Parametrii caracteristici circuitului de apa (, , i4, i1, ρm apa, Qc), Δtlg si krad au fost determinati in acelasi mod ca si in prima etapa de experimentare. Trebuie facuta mentiunea ca, in aceasta a doua etapa, t2 si t3 sunt temperaturi ale aerului umed.
Pe langa Δh1 si t1.t4, masurate si in prima etapa, in etapa a doua se masoara si diferentele de coloana de apa Δh2 si Δh4; de asemenea, se cronometreaza intervalul de timp τp in care se produce diferenta de nivel Δh4 a apei din vasul V1. Acesti trei parametrii permit stabilirea consumului de apa al instalatiei (a debitului de apa pulverizata) si a vitezei wrad.
Modul de determinare a debitului volumic de apa pulverizata este prezentat in subcap. 2.2. Valorile acestui parametru se obtin cu ajutorul formulei
, (23)
unde este volumul corespunzator unui milimetru de coloana de apa din vasul V1.
Debitul masic de apa pulverizata este dat de expresia
, (24)
in care reprezinta densitatea apei pulverizate, corespunzatoare unei temperaturi temperaturii medii de 26˚C, mentinuta aproximativ constanta pe parcursul experimentelor.
Datorita diferentelor de temperatura mai mici de un grad (maxim 0,761 in tabelul 3 - valoare corespunzatoare unei variatii de densitate sub 0,3%), densitatea aerului in sectiunea STA se exprima astfel:
. (25)
Cunoscand si viteza wTA, calculata cu formula (8), se pot determina debitul volumic de aer in sectiunea STA
(26)
si debitul masic de aer uscat:
. (27)
Este acum posibila determinarea parametrului ψ, dat de relatia (22). Valoarea maxima a lui ψ obtinuta in cadrul experimentelor este 3,24∙10-3 kgapa/kgaer, ceea ce indica faptul ca, din punct de vedere al debitului volumic (de care depinde viteza), prezenta apei pulverizate practic nu conteaza. Aceasta afirmatie este sustinuta de faptul ca oprirea pulverizarii nu a determinat modificarea valorilor diferentelor de coloana de apa Δh1 (care determina wrad) si Δh2 (in functie de care se calculeaza wTA) la nici nici unul din regimurile de functionare analizate. Prin urmare, wrad se determina in functie de Δh1, ca viteza a aerului uscat. Relatiile de calcul utilizate sunt (9), (10) si (11).
Puterea consumata de ventilator datorita radiator se calculeaza cu formula
, (28)
in care este densitatea aerului uscat corespunzatoare valorilor medii pm rad si Tm rad; se calculeaza cu relatia (4). Consumul energetic erad este dat de raportul (21).
5. Rezultatele experientelor si interpretarea lor
5.1. Rezultate obtinute in prima etapa de experimentare
Atat valorile masurate direct cat si cele calculate in prima etapa a cercetarii experimentale sunt date in tabelul Desi prezentarea parametrilor masurati si calculati in aceasta etapa precum si a metodicii de calcul a fost facuta in subcapitolele 2.1 si 1, pentru urmarirea cu usurinta a datelor experimentale se impune o prezentare succinta a parametrilor prezentati in tabelul
Parametrii masurati, in ordinea din tabel, sunt:
Δh1 - caderea de presiune pe radiator [mm H2O];
TC1 - indicatia termocuplului plasat la iesirea din radiator, pe circuitul de apa;
TC2 - indicatia termocuplului plasat dupa radiator, pe circuitul de aer;
TC3 - indicatia termocuplului plasat inainte de radiator, pe circuitul de aer;
TC4 - indicatia termocuplului plasat la intrarea in radiator, pe circuitul de apa,;
Vi - valoarea indicata de debitmetrul volumetric de apa la inceputul experientei [m3];
Tabelul
Nr.
crt.
Parametrii masurati
Parametrii calculati
Δh1
[mm H2O]
TC1
TC2
TC3
TC4
Vi
[m3]
Vf
[m3]
τ
[s]
t1
[˚C]
t2
[˚C]
t3
[˚C]
t4
[˚C]
[kg/s]
Qc
[kW]
krad
[W/m2∙K]
kc rad
[W/m2∙K]
p3
[bar]
i2
[kJ/kg]
i3
[kJ/kg]
[kg/s]
wrad
[m/s]
Pc
[kW]
erad
[kJ/kg]
1.
12
45,1
43
28
51,7
95,144
95,221
743
43,84
41,75
26,76
50,44
0,103
2,823
91,07
92,31
1,014
314,5
299,6
0,189
12,89
0,024
0,128
2.
21,5
43
41,1
27,7
49,6
94,992
95,052
572
41,75
39,85
26,46
48,34
0,104
2,859
99,1
99,37
1,015
312,6
299,3
0,215
14,56
0,049
0,229
3.
33,5
41,4
39,4
27,7
47,9
94,922
94,985
599
40,15
38,15
26,46
46,64
0,104
2,825
103,97
106,87
1,016
310,9
299,3
0,243
16,42
0,086
0,355
45
39,7
38,4
27,5
46,4
94,844
94,917
703
38,45
37,15
26,26
45,14
0,103
2,876
115,86
112,60
1,017
309,9
299,1
0,266
17,9
0,126
0,476
5.
58,5
38,4
37,4
27,4
45,1
94,669
94,751
794
37,15
36,15
26,16
43,84
0,102
2,861
124
118,08
1,019
309
299
0,288
19,37
0,178
0,617
6.
79
37,2
36,1
26,9
44
94,574
94,659
827
35,95
34,85
25,66
42,75
0,102
2,891
127,38
124,73
1,021
307,7
298,5
0,317
21,19
0,263
0,83
7.
86,5
36,6
35,4
27,1
43,2
94,456
94,56
984
35,35
34,15
25,86
41,95
0,105
2,886
134,2
132,43
1,021
307
298,7
0,35
23,42
0,318
0,907
8.
100
35,7
34,6
26,8
42,5
94,354
94,439
825
34,45
33,36
25,56
41,25
0,102
2,899
138,38
137,69
1,023
306,2
298,4
0,375
24,98
0,392
1,046
9.
114,5
35,4
34,3
27,2
42
94,269
94,348
751
34,15
33,06
25,96
40,75
0,105
2,873
144,95
144,83
1,024
305,9
298,8
0,408
27,18
0,488
1,197
10.
128
34,8
33,9
27,4
41,5
94,195
94,257
591
33,56
32,66
26,16
40,25
0,104
2,909
155,53
153,77
1,026
305,5
299
0,451
30,03
0,603
1,337
11.
144
34,4
33,6
27,3
41,2
94,109
94,18
673
33,16
32,36
26,06
39,95
0,105
2,969
162,05
158,50
1,027
305,2
298,9
0,475
31,58
0,714
1,502
12.
160,5
34,1
33,2
27,5
40,8
94,043
94,102
565
32,86
31,96
26,26
39,55
0,104
2,896
163,83
165,70
1,029
304,8
299,1
0,512
34,01
0,857
1,673
13.
179,5
33,5
32,9
27,5
40,3
93,991
94,039
467
32,26
31,66
26,26
39,05
0,102
2,893
173,8
170,94
1,031
304,5
299,1
0,54
35,82
1,009
1,868
1
196
33,3
32,7
27,6
40
93,909
93,97
584
32,06
31,46
26,36
38,75
0,104
2,897
179,64
176,94
1,032
304,3
299,2
0,573
37,94
1,167
2,038
15.
211,5
33
32,4
27,6
39,7
93,842
93,904
591
31,76
31,16
26,36
38,45
0,104
2,91
185,16
183,87
1,034
304
299,2
0,611
40,44
1,342
2,196
16.
225
32,8
32,2
27,6
39,7
93,575
93,639
623
31,56
30,96
26,36
38,45
0,102
2,935
187,39
189,49
1,035
303,8
299,2
0,643
42,52
1,501
2,334
Fig. 7.
Fig. 8.
Vf - valoarea indicata de debitmetrul volumetric de apa la sfarsitul experientei [m3];
τ - timpul scurs intre indicatiile Vi si Vf ale debitmetrului volumetric de apa [s].
Parametrii calculati sunt urmatorii:
t1.t4 - temperaturile corespunzatoare indicatiilor termocuplurilor TC1.TC4 [˚C];
- debitul masic de apa [kg/s];
Qc - cantitatea de caldura cedata de apa [kW];
krad - coeficientul global de schimb de caldura [W/m2∙K];
kc rad - coeficientul global de schimb de caldura dat de subrutina RADIATOR [W/m2∙K];
p3 - presiunea aerului inainte de radiator [bar];
i2 - entalpia aerului inainte de radiator [kJ/kg];
i3 - entalpia aerului dupa radiator [kg/m3];
- debitul masic de aer [kg/s];
wrad - viteza aerului prin radiator [m/s];
Pc - puterea consumata de ventilator datorita caderii de presiune pe radiator [kW];
erad - puterea specifica consumata de ventilator [kJ/kg aer].
Pe baza datelor din tabelul 4 au fost trasate graficele prezentate in fig. 7 si 8. Curbele de variatie reprezentate in aceste figuri au evolutiile preconizate: cresterea vitezei wrad conduce la imbunatatirea schimbului de caldura (krad creste - fig. 7) dar si la cresterea caderii de presiune pe radiator (Δh1 creste), ceea ce determina cresterea consumului energetic al ventilatorului (Pc si erad cresc - fig. 8). Valorile lui kc rad, calculate de subrutina RADIATOR in functie de viteza si temperaturile masurate, corespunzatoare fiecarui regim, sunt foarte apropiate de valorile lui krad obtinute experimental. Diferentele, calculate cu relatia
, (29)
sunt de maxim 5%. Diferenta medie, pentru toate cele 16 regimuri experimentate, este de 1,62%.
In cazul configuratiei 6, in punctul caracterizat de randamentul maxim obtinut in cadrul cercetarii analitice, valoarea debitului specific de aer necesar racirii apei din radiator este de de energie produsa. Regimul optim din punct de vedere al schimbului de caldura (descris de valoarea maxima ), dintre cele analizate experimental, este caracterizat de viteza maxima si (de asemenea, valoare maxima). Rezulta ca, la acest regim, radiatorul determina un consum de 1 kW pentru fiecare 30 kW produsi de grupul binar, ceea ce inseamna o scadere a randamentului de 3,3%. In cazul regimului caracterizat de viteza minima , consumul specific de putere datorat radiatorului este . In acest caz, radiatorul determina un consum de 1 kW la 547,7 kW produsi de grup, ceea ce inseamna reducerea randamentului cu 0,18%.
Cele doua regimuri discutate reprezinta situatiile extreme experimentate. Ceea ce rezulta in urma analizei este ca o crestere de 3,3 ori a vitezei (de la 12,89 m/s la 42,52 m/s) determina o crestere de 2,1 ori a lui krad (de la 91,07 W/m2∙K la 187,39 W/m2∙K) cu pretul cresterii consumului energetic datorat radiatorului de 18,3 ori, de la 0,128 kJ/kg aer (corespunzand unei reduceri a randamentului global cu 0,18%) la 2,33 kJ/kg aer (corespunzand unei reduceri a randamentului global cu 3,3%).
5.2. Rezultate obtinute in etapa a doua de experimentare
Valorile masurate si cele calculate in aceasta etapa sunt prezentate in tabelul 5. Parametrii masurati, in ordinea din tabel, au urmatoarele semnificatii:
Δh1 - caderea de presiune pe radiator [mm H2O];
Δh2 - caderea de presiune pe tubul Anubar [mm H2O];
Δh4 - diferenta de coloana de apa corespunzatoare consumului de apa pulverizata [mm H2O];
τp - timpul de pulverizare corespunzator diferentei de coloana de apa Δh4 [s];
TC1 - indicatia termocuplului plasat la iesirea din radiator, pe circuitul de apa,;
TC2 - indicatia termocuplului plasat dupa radiator, pe circuitul de aer;
TC3 - indicatia termocuplului plasat inainte de radiator, pe circuitul de aer;
TC4 - indicatia termocuplului plasat la intrarea in radiator, pe circuitul de apa,;
tamb - temperatura mediului ambiant [˚C];
Vi - valoarea indicata de debitmetrul volumetric de apa la inceputul experientei [m3];
Vf - valoarea indicata de debitmetrul volumetric de apa la sfarsitul experientei [m3];
τ - timpul scurs intre indicatiile Vi si Vf ale debitmetrului volumetric de apa [s].
Tabelul 5.
Nr.
crt.
Parametrii marurati
Parametrii calculati
Δh1
Δh2
Δh4
τp
[s]
TC1
TC2
TC3
TC4
tamb
[˚C]
Vi
[m3]
Vf
[m3]
τ
[s]
t1
[˚C]
t2
[˚C]
t3
[˚C]
t4
[˚C]
[kg/s]
∙104
[kg/s]
Qc
[kW]
krad
[W/m2∙K]
p3
[bar]
[kg/s]
ψ∙104
wrad
[m/s]
Pc
[kW]
erad
[kJ/kg]
1.
98
10,5
45
1189
42,2
38,2
27
50,1
28,1
103,061
103,15
967
40,95
36,95
25,76
48,84
0,091
1,893
3,003
89,31
1,023
0,371
5,1
24,55
0,381
1,031
2.
40
647
38,9
34,7
26,9
46,9
28
103,171
103,215
475
37,65
33,46
25,66
45,64
0,092
3,092
3,063
101,52
0,371
8,33
3.
50
595
39,1
35,2
27,5
47,2
28,5
103,226
103,27
477
37,85
33,95
26,26
45,94
0,091
4,202
3,088
104,96
0,37
11,35
70
689
37,9
34,6
27,9
45,9
29,1
103,282
103,322
437
36,65
33,36
26,66
44,64
0,091
5,081
3,028
114,14
0,37
13,75
5.
95
683
36,4
34,9
27,3
44,3
28,7
103,359
103,39
337
35,15
33,66
26,06
43,05
0,091
6,956
3,006
130,32
0,37
18,79
6.
120
764
35,8
38,6
29,1
43,7
30,5
103,401
103,455
596
34,55
37,35
27,86
42,45
0,09
7,855
2,961
202,49
0,368
21,35
7.
130
719
35,3
39
28,4
43,3
29,8
103,47
103,519
534
34,05
37,75
27,16
42,05
0,091
9,042
3,037
221,51
0,369
24,52
8.
150
721
34,6
38,6
28
42,8
29,5
103,529
103,575
520
33,36
37,35
26,76
41,55
0,088
10,404
3,002
226,67
0,369
28,18
9.
160
671
34,9
38,5
28,5
43
30
103,59
103,667
845
33,66
37,25
27,26
41,75
0,091
11,924
3,054
227,02
0,369
32,36
10.
127
14,5
20
739
41,4
37,6
27,3
49,5
28,3
102,486
102,538
569
40,15
36,35
26,06
48,24
0,091
1,353
3,058
94,55
1,025
0,429
3,15
28,38
0,572
1,334
11.
40
712
40,4
38,9
26,9
48,1
28
102,544
102,595
541
39,15
37,65
25,66
46,84
0,093
2,809
2,999
107,26
0,43
6,54
12.
60
756
38,7
36,5
27,3
46,9
28,4
102,598
102,65
576
37,45
35,25
26,06
45,64
0,09
3,969
3,06
112,65
0,429
9,25
13.
70
672
37
34,6
26,8
45,4
27,9
102,658
102,721
716
35,75
33,36
25,56
44,14
0,087
5,209
3,056
116,76
0,43
12,12
1
80
621
35,6
34,8
27
44
28,2
102,726
102,776
581
34,35
33,56
25,76
42,75
0,085
6,442
2,99
134,79
0,429
15
15.
115
791
34,5
37,3
27,6
42,9
28,8
102,786
102,838
610
33,26
36,05
26,36
41,65
0,085
7,27
2,963
190,79
0,429
16,97
16.
120
744
34,2
37,5
28,7
42,3
30
102,849
102,894
501
32,96
36,25
27,46
41,05
0,089
8,066
3,011
234,79
0,427
18,9
17.
130
689
34,3
37,4
29
42,5
30,3
102,905
102,951
511
33,06
36,15
27,76
41,25
0,089
9,435
3,055
235,59
0,426
22,13
18.
110
531
33,8
37,3
28,2
41,9
29,7
102,966
103,02
600
32,56
36,05
26,96
40,65
0,089
10,359
3,017
237,99
0,427
24,25
19.
159,5
21,5
30
911
41,3
37,3
26,9
49,3
27,6
101,562
101,627
692
40,05
36,05
25,66
48,04
0,093
1,647
3,104
94,56
1,029
0,524
3,14
34,51
0,873
1,676
20.
50
739
40,2
38,4
27,3
48
28
101,641
101,69
533
38,95
37,15
26,06
46,74
0,091
3,383
2,963
106,39
0,524
6,46
21.
80
788
38,7
35
29,1
46,5
29,9
101,704
101,757
575
37,45
33,76
27,86
45,24
0,091
5,077
2,972
113,27
0,52
9,76
22.
130
1045
37,2
34,7
27,8
45,1
28,7
100,958
101,055
1078
35,95
33,46
26,56
43,84
0,089
6,221
2,94
119,17
0,522
11,91
23.
110
766
36,3
35,5
29,1
44
30
101,067
101,119
598
35,05
34,25
27,86
42,75
0,086
7,181
2,77
141,8
0,52
13,8
2
120
770
35,1
37,8
28,5
42,9
29,5
101,129
101,178
528
33,85
36,55
27,26
41,65
0,092
7,793
2,995
206,43
0,521
14,96
25.
130
749
33,9
37,1
28,9
41,5
29,8
101,188
101,239
534
32,66
35,85
27,66
40,25
0,095
8,679
3,004
256,66
0,521
16,67
26.
140
732
34,9
38,5
29,8
42,6
30,8
101,26
101,308
498
33,66
37,25
28,56
41,35
0,096
9,564
3,071
268,77
0,519
18,43
27.
160
748
34,8
38,4
29,5
42,5
30,7
101,32
101,371
538
33,56
37,15
28,26
41,25
0,094
10,697
3,02
259,11
0,519
20,61
28.
160
678
34,6
38,2
29,5
42,3
30,8
101,384
101,444
645
33,36
36,95
28,26
41,05
0,092
11,801
2,964
259,41
0,519
22,75
Tabelul 5 - continuare
Nr.
crt.
Parametrii marurati
Parametrii calculati
Δh1
Δh2
Δh4
τp
[s]
TC1
TC2
TC3
TC4
tamb
[˚C]
Vi
[m3]
Vf
[m3]
τ
[s]
t1
[˚C]
t2
[˚C]
t3
[˚C]
t4
[˚C]
[kg/s]
∙104
[kg/s]
Qc
[kW]
krad
[W/m2∙K]
p3
[bar]
[kg/s]
ψ∙104
wrad
[m/s]
Pc
[kW]
erad
[kJ/kg]
29.
186
24
35
761
41
39,5
27,5
48,7
28,2
101,562
101,627
692
39,75
38,25
26,26
47,44
0,093
2,3
2,988
106,86
1,031
0,557
4,13
36,62
1,082
1,952
30.
55
830
39,5
37,8
27,4
47,4
28,2
101,641
101,69
533
38,25
36,55
26,16
46,14
0,091
3,314
3,002
111,43
0,557
5,95
31.
85
877
38
36
27,6
46,1
28,4
101,704
101,757
575
36,75
34,75
26,36
44,84
0,091
4,847
3,087
120,77
0,556
8,71
32.
85
730
36,5
34,5
27,9
44,7
28,7
100,958
101,055
1078
35,25
33,26
26,66
43,45
0,089
5,823
3,052
130,5
0,556
10,48
33.
95
713
35,2
36,2
28,5
43,5
29,3
101,067
101,119
598
33,95
34,95
27,26
42,25
0,086
6,663
2,986
171,16
0,555
12,01
3
115
742
34,5
38
29,6
42,4
30,6
101,129
101,178
528
33,26
36,75
28,36
41,15
0,092
7,75
3,034
261,88
0,552
14,03
35.
135
748
34
37,9
29,1
42
30,1
101,188
101,239
534
32,76
36,65
27,86
40,75
0,095
9,025
3,162
282,52
0,553
16,32
36.
145
745
33,8
37,3
29,2
41,5
30,3
101,26
101,308
498
32,56
36,05
27,96
40,25
0,096
9,733
3,072
280,14
0,553
17,61
37.
145
683
34
37,6
29,5
41,7
30,5
101,32
101,371
538
32,76
36,35
28,26
40,45
0,094
10,616
3,021
281,92
0,552
19,22
38.
155
676
33,9
37,7
29,3
41,9
30,4
101,384
101,444
645
32,66
36,45
28,06
40,65
0,092
11,466
3,08
280,88
0,553
20,75
39.
210,5
30
45
809
38,3
35,2
27,2
46,1
27,8
100,147
100,227
843
37,05
33,95
25,96
44,84
0,094
2,782
3,06
111,56
1,034
0,635
4,38
41,62
1,391
2,193
40.
65
795
36,5
33,9
26,4
44,3
27,2
100,254
100,305
541
35,25
32,66
25,16
43,05
0,094
4,089
3,041
118,99
0,636
6,43
41.
85
780
35,1
33,2
25,8
42,9
26,5
100,319
100,378
637
33,85
31,96
24,56
41,65
0,092
5,449
2,989
126,19
0,637
8,55
42.
105
827
33,7
33,1
26,3
41,6
27,1
100,396
100,448
563
32,46
31,86
25,06
40,35
0,092
6,349
3,02
152,57
0,636
9,98
43.
115
794
32,9
35,4
27,1
40,5
27,9
100,461
100,515
573
31,66
34,15
25,86
39,25
0,094
7,243
2,965
218,46
0,634
11,42
4
135
842
32,3
35,6
28,1
40,1
28,9
100,537
100,609
779
31,06
34,35
26,86
38,85
0,092
8,018
2,985
275,26
0,632
12,68
45.
125
711
32,3
35,5
28,1
39,8
29
100,611
100,663
535
31,06
34,25
26,86
38,55
0,097
8,792
3,018
284,8
0,632
13,91
46.
135
719
31,7
35
27,6
39,4
28,6
100,682
100,731
517
30,46
33,76
26,36
38,15
0,094
9,389
3,022
285,25
0,633
14,84
47.
145
731
32,5
35,9
27,9
40
28,9
100,744
100,806
625
31,26
34,65
26,66
38,75
0,099
9,919
3,08
284,27
0,632
15,69
48.
221
33
55
844
41,3
42,1
26,3
48,9
26,9
103,683
103,767
879
40,05
40,85
25,06
47,64
0,095
3,259
3
116,08
1,035
0,675
4,83
44,09
1,555
2,312
49.
85
889
39,9
40,5
26,9
47,6
27,5
103,775
103,839
667
38,65
39,25
25,66
46,34
0,095
4,781
3,053
125,52
0,673
7,1
50.
100
885
38,8
41
25,8
46,5
26,6
103,857
103,927
749
37,55
39,75
24,56
45,24
0,093
5,651
2,975
136,82
0,675
8,37
51.
110
874
35,1
36,8
27,1
42,7
27,9
103,941
104,023
859
33,85
35,55
25,86
41,45
0,095
6,294
3,002
174,53
0,673
9,36
52.
120
854
33,3
36,9
28,1
41,1
28,8
104,045
104,115
741
32,06
35,65
26,86
39,85
0,094
7,027
3,05
261,21
0,671
10,48
53.
130
842
32,4
35,9
28
40
28,8
104,132
104,195
657
31,16
34,65
26,76
38,75
0,095
7,721
3,017
284,81
0,671
11,51
5
140
838
31,6
35,1
27,2
39,2
28,2
104,205
104,254
501
30,36
33,85
25,96
37,95
0,097
8,354
3,078
290,55
0,672
12,43
55.
150
822
31,8
35,3
27,7
39,5
28,6
104,271
104,329
617
30,56
34,05
26,46
38,25
0,094
9,125
2,997
289,63
0,671
13,6
56.
160
807
31,5
35
27,3
39,3
28,3
104,338
104,395
609
30,26
33,76
26,06
38,05
0,093
9,915
3,023
285,26
0,672
14,76
Tabelul 5 - continuare
Nr.
crt.
Parametrii marurati
Parametrii calculati
Δh1
Δh2
Δh4
τp
[s]
TC1
TC2
TC3
TC4
tamb
[˚C]
Vi
[m3]
Vf
[m3]
τ
[s]
t1
[˚C]
t2
[˚C]
t3
[˚C]
t4
[˚C]
[kg/s]
∙104
[kg/s]
Qc
[kW]
krad
[W/m2∙K]
p3
[bar]
[kg/s]
ψ∙104
wrad
[m/s]
Pc
[kW]
erad
[kJ/kg]
57.
231
36,5
45
865
38,4
36,4
27,1
46,1
27,6
99,266
99,347
845
37,15
35,15
25,86
44,84
0,095
2,602
3,051
116,76
1,036
0,718
3,63
46,96
1,725
2,405
58.
65
805
36,9
35,6
26,5
44,5
27,3
99,375
99,441
693
35,65
34,35
25,26
43,25
0,095
4,038
2,994
124,65
0,718
5,62
59.
85
862
34,8
33,5
25,8
42,7
26,7
99,497
99,564
722
33,56
32,26
24,56
41,45
0,092
4,931
3,033
133,68
0,72
6,85
60.
105
832
33,6
33,9
26
41,2
26,8
99,572
99,667
1003
32,36
32,66
24,76
39,95
0,094
6,311
2,98
160,39
0,719
8,77
61.
115
797
32,1
35,8
27,3
39,9
28
99,675
99,741
699
30,86
34,55
26,06
38,65
0,094
7,216
3,049
275,36
0,717
10,07
62.
125
840
31,7
35,4
27,5
39,5
28,1
99,755
99,816
645
30,46
34,15
26,26
38,25
0,094
7,442
3,055
295,17
0,716
10,39
63.
135
838
31,9
35,3
27,7
39,4
28,5
99,9
99,952
525
30,66
34,05
26,46
38,15
0,099
8,056
3,076
297,24
0,715
11,26
6
145
823
31,7
35,6
27,4
39,7
28,4
99,962
100,016
584
30,46
34,35
26,16
38,45
0,092
8,81
3,063
292,5
0,716
12,31
65.
155
806
31,9
35,5
27,5
39,6
28,5
100,031
100,081
531
30,66
34,25
26,26
38,35
0,094
9,617
3,002
283,39
0,715
13,44
66.
251,5
40,5
50
1197
40
40
27,3
48,1
27,9
98,507
98,555
536
38,75
38,75
26,06
46,84
0,089
2,089
2,997
117,52
1,038
0,768
2,72
50,22
2,011
2,616
67.
50
1073
38,7
37,4
26,8
46,3
27,5
98,282
98,351
708
37,45
36,15
25,56
45,04
0,097
2,33
3,062
118,89
0,769
3,03
68.
60
1059
37,2
35,8
26,3
45
27,2
98,584
98,636
551
35,95
34,55
25,06
43,75
0,094
2,833
3,044
121,75
0,77
3,68
69.
70
862
35,8
33,8
27,1
43,5
27,9
98,667
98,704
397
34,55
32,56
25,86
42,25
0,093
4,061
2,969
129,53
0,768
5,29
70.
80
822
34,3
32,7
26,3
42,1
27,2
98,737
98,796
632
33,06
31,46
25,06
40,85
0,093
4,867
3,013
139,19
0,77
6,32
71.
90
823
32,9
32,6
26,4
40,8
27,2
98,805
98,877
776
31,66
31,36
25,16
39,55
0,092
5,469
3,034
166,29
0,77
7,1
72.
110
872
31,7
35,3
27,1
39,5
27,9
98,901
98,945
471
30,46
34,05
25,86
38,25
0,093
6,308
3,017
275,18
0,768
8,21
73.
130
937
31,3
35
27,2
39,1
28,1
98,972
99,033
639
30,06
33,75
25,96
37,85
0,095
6,938
3,084
301,59
0,767
9,04
7
150
992
31,7
35,4
27,4
39,5
28,3
99,085
99,155
732
30,46
34,15
26,16
38,25
0,095
7,562
3,089
294,95
0,767
9,86
75.
261
43
45
1303
36
34,7
25,1
44,1
25,9
104,395
104,447
577
34,75
33,46
23,86
42,85
0,090
1,727
3,02
119,51
1,039
0,805
2,15
52,26
2,167
2,699
76.
50
1137
35,6
33,2
25,4
43,3
26,1
104,457
104,531
763
34,35
31,96
24,16
42,05
0,096
2,199
3,09
122,06
0,805
2,73
77.
60
1127
34,5
31,1
26,1
42
26,8
104,543
104,599
578
33,26
29,86
24,86
40,75
0,096
2,662
3,007
125,66
0,803
3,32
78.
70
981
33,1
30,8
24,5
41
25,3
104,609
104,651
445
31,86
29,56
23,27
39,75
0,094
3,568
3,086
131,99
0,807
4,42
79.
80
957
32,7
30,5
25,1
40,5
25,9
104,662
104,727
687
31,46
29,26
23,86
39,25
0,094
4,18
3,055
140,09
0,805
5,19
80.
90
931
32,4
32,6
25,8
40
26,6
104,734
104,812
827
31,16
31,36
24,56
38,75
0,094
4,834
2,968
170,2
0,803
6,02
81.
110
944
31,6
35,3
26,8
39,4
27,7
104,821
104,861
427
30,36
34,05
25,56
38,15
0,093
5,827
3,026
273,22
0,8
7,28
82.
130
1032
30,7
34,4
26,6
38,5
27,6
104,874
104,938
685
29,46
33,16
25,36
37,25
0,093
6,299
3,018
302,3
0,8
7,87
83.
150
1066
30,9
34,4
26,7
38,5
27,7
104,949
105,031
850
29,66
33,16
25,46
37,25
0,096
7,037
3,037
293,44
0,8
8,79
8
150
956
30,6
34,6
25,8
38,7
27,7
105,047
105,105
631
29,36
33,36
24,56
37,45
0,091
7,846
3,084
278,46
0,8
9,81
Tabelul 5 - continuare
Nr.
crt.
Parametrii marurati
Parametrii calculati
Δh1
Δh2
Δh4
τp
[s]
TC1
TC2
TC3
TC4
tamb
[˚C]
Vi
[m3]
Vf
[m3]
τ
[s]
t1
[˚C]
t2
[˚C]
t3
[˚C]
t4
[˚C]
[kg/s]
∙104
[kg/s]
Qc
[kW]
krad
[W/m2∙K]
p3
[bar]
[kg/s]
ψ∙104
wrad
[m/s]
Pc
[kW]
erad
[kJ/kg]
85.
272,5
45,5
30
1300
40,1
41,9
25,9
47,9
26,7
97,447
97,529
885
38,85
40,65
24,66
46,64
0,092
1,154
2,987
125,82
1,04
0,835
1,38
54,29
2,366
2,836
86.
50
1308
39,3
40,6
26,8
47,1
27,6
97,7
97,776
825
38,05
39,35
25,56
45,84
0,091
1,912
2,97
129,8
0,832
2,3
87.
50
1018
37,2
37,1
26,5
45,1
27,2
97,557
97,639
867
35,95
35,85
25,26
43,84
0,094
2,456
3,09
133,45
0,833
2,95
88.
70
1060
35,7
36,5
25,6
44
26,5
97,905
97,955
564
34,45
35,25
24,36
42,75
0,088
3,302
3,044
139,72
0,835
3,95
89.
80
1078
34,9
36,1
24,9
42,7
25,7
97,813
97,876
664
33,66
34,85
23,67
41,45
0,094
3,711
3,062
150,07
0,838
4,43
90.
100
1109
32,5
35
25,8
41,2
26,5
98,011
98,054
512
31,26
33,76
24,56
39,95
0,084
4,509
3,025
188,14
0,835
5,4
91.
100
900
30,9
35
26,8
39,1
27,7
98,084
98,127
471
29,66
33,76
25,56
37,85
0,091
5,556
3,1
303,24
0,832
6,68
92.
100
789
30,5
34,6
26,3
38,9
27,3
98,138
98,188
571
29,26
33,36
25,06
37,65
0,087
6,338
3,047
287,46
0,833
7,61
93.
290
48
50
1269
36,3
36,1
26,2
44
27
96,67
96,738
717
35,05
34,85
24,96
42,75
0,094
1,97
3,021
135,26
1,041
0,867
2,27
56,31
2,603
3,005
9
100
1504
35,1
34,9
26,4
42,7
27,2
96,754
96,835
833
33,85
33,66
25,16
41,45
0,097
3,325
3,058
148,77
0,866
3,84
95.
100
1239
33,9
34,6
26,5
41,6
27,3
96,864
96,934
737
32,66
33,36
25,26
40,35
0,094
4,036
3,027
168,62
0,866
4,66
96.
100
1072
33,1
35
26,3
40,5
27,2
96,976
97,034
577
31,86
33,76
25,06
39,25
0,1
4,665
3,079
201,53
0,866
5,39
97.
100
992
31,8
34,4
26,1
39,3
27
97,176
97,229
540
30,56
33,16
24,86
38,05
0,098
5,041
3,048
231,07
0,867
5,82
98.
100
885
30,5
33,9
26,4
38
27,2
97,095
97,157
625
29,26
32,66
25,16
36,75
0,099
5,651
3,082
308,64
0,866
6,53
99.
100
792
30,4
33,9
26,2
38,1
27
97,266
97,31
455
29,16
32,66
24,96
36,85
0,096
6,314
3,084
308,9
0,867
7,29
Fig. 9.
Fig. 10.
Parametrii calculati sunt urmatorii:
t1.t4 - temperaturile corespunzatoare indicatiilor termocuplurilor TC1.TC4 [˚C];
- debitul masic de apa prin radiator [kg/s];
- debitul masic de apa pulverizata [kg/s];
Qc - cantitatea de caldura cedata de apa [kW];
krad - coeficientul global de schimb de caldura [W/m2∙K];
p3 - presiunea aerului inainte de radiator [bar];
- debitul masic de aer uscat [kg/s];
ψ - concentratia de apa pulverizata in aer [kg apa/kg aer uscat];
wrad - viteza aerului prin radiator [m/s];
Pc - puterea consumata de ventilator datorita caderii de presiune pe radiator [kW];
erad - puterea specifica consumata de ventilator [kJ/kg aer].
Pe baza datelor din tabelul 5 au fost reprezentate curbele de variatie , trasate pentru valori multiple ale vitezei wrad (fig. 9). Se observa ca, indiferent de viteza, exista o valoare optima a lui ψ (ψopt), careia ii corespunde o valoare optima a lui krad (krad opt). La viteze mici, ψopt este valoarea limita de la care cresterea importanta a lui ψ determina o variatie nesemnificativa a lui krad. La viteze mari, ψopt corespunde valorii maxime krad max.
Fig. 9 indica, de asemenea, ca prin cresterea lui wrad se obtin cresterea lui krad max si scaderea lui ψopt. Aceasta evolutie este avantajoasa deoarece inseamna intensificarea schimbului de caldura odata cu reducerea consumului de apa. Astfel, valoarea maxima , obtinuta in etapa a doua, corespunde vitezei maxime a aerului atinse in radiator si optimului minim . Rezulta ca, pentru grupul binar de randamentul maxim, caracterizat de un debit specific de aer necesar racirii apei din radiator de de energie produsa, consumul specific de apa este 1,04∙10-5 kg apa/kJ de energie produsa. Pentru un grup binar de 150 kW, rezulta un consum efectiv de apa pentru pulverizare de . Tinand cont de avantajele tehnico-economice datorate intensificarii schimbului de caldura prin radiator (economie de material, masa si gabarit reduse ale radiatorului, deci ale instalatiei) acest consum de apa este nesemnificativ. Din pacate, datorita vitezei mari, caderea de presiune pe radiator este mare (290 mm H2O), ceea ce determina un consum energetic ridicat al ventilatorului: (fig. 10). In cazul grupului cu performante optime rezulta o reducere a randamentului cu 4,3%.
Consumul energetic minim () corespunde, evident, vitezei minime ; in acest caz, punctul optim de pe curba de viteza constanta (fig. 9) este descris de si . Reducerea de randament a grupului binar optim ca performante, impusa de radiator, scade la 1,47%, iar consumul specific de apa pentru pulverizare devine 3,64∙10-5 kg apa/kJ de energie produsa. Pentru un grup binar de 150 kW, se ajunge la un consum efectiv de apa pentru pulverizare de . Deoarece krad obtinut inacest caz este de 1,6 ori mai mare decat cel determinat in prima etapa, in absenta pulverizarii apei si pentru aceeasi viteza wrad (cand ), se poate afirma ca si acest ultim consum de apa este acceptabil.
6. Concluzii asupra experientelor efectuate
Diferentele mici intre valorile coeficientului global de schimb de caldura determinate experimental (krad) si cele date de subrutina RADIATOR (kc rad) valideaza atat procedura de calcul a suprafetelor de schimb de caldura dezvoltata in cadrul programelor din pachetul BINAR cat si metodica de experimentare caracteristica primei etape.
Rezultatele obtinute in cea de a doua etapa de experimentare confirma faptul ca pulverizarea apei in curentul de aer permite cresterea semnificativa a coeficientului global de schimb de caldura (de 1,5 1,65 ori pentru aceeasi valoare a lui waer). Rezultatul este reducerea considerabila a suprafetei de schimb de caldura a radiatorului, care compenseaza, atat din punct de vedere tehnic (reducerea masei si a gabaritului instalatiei) cat si economic (costuri mai mici, datorita economiei de material), neajunsul pe care il reprezinta consumul de apa. Pierderile energetice datorate caderii de presiune pe radiator sunt similare celor determinate in absenta pulverizarii, in prima etapa de cercetare experimentala.
Capitolul 5 - CONCLUZII
Se reaminteste faptul ca, pana in prezent, grupurile binare sunt folosite in domeniul transporturilor doar pentru propulsia navala. In cadrul grantului, cercetarile s-au concentrat asupra posibilitatilor de utilizare a grupurilor binare in domeniul transportului terestru - feroviar si auto.
Rezultatele activitatii de cercetare teoretica si experimentala desfasurata in cadrul grantului conduc la formularea urmatoarelor concluzii:
7. Diferentele mici intre valorile coeficientului global de schimb de caldura determinate experimental si cele date de subrutina RADIATOR valideaza atat procedura de calcul a suprafetelor de schimb de caldura utilizata in programele de calcul BINAR, cat si metodica de experimentare caracteristica primei etape.
8. Utilizarea sistemului de condensare Heller, in conditiile analizei teoretice efectuate, permite utilizarea grupurilor binare pentru propulsia feroviara, precum si pentru propulsia autovehiculelor grele.
9. Cresterea vitezei wrad de la 8 m/s (valoare considerata in calculele efectuate cu programele din pachetul BINAR) la 24,55 m/s (valoarea minima obtinuta in etapa a doua de experimentare) determina cresterea lui krad de la 67,2 W/m2∙K la 136,5 W/m2∙K, ceea ce conduce la reducerea suprafetei specifice SR sp de la 0,446 m2/kW la 0,228 m2/kW. Masa si volumul radiatorului instalatiei de condensare a unui grup binar de 100 kW s-ar reduce de la 174 kg / 0,4 m3 (v. subcap. 3.5) la 89 kg / 0,2 m3, valori care creeaza inca dificultati, dar pot fi acceptate in cazul unui automobil. Dupa cum s-a aratat in subcap. 5.2, viteza wrad = 24,55 m/s determina reducerea cu 1,47% a randamentului global al grupului binar, care devine 45,7%. Dar, chiar si acest randament este superior motoarelor diesel de ultima generatie, de putere egala.
10. In cazul pulverizarii apei in curentul de aer de racire, pierderile energetice datorate caderii de presiune pe radiator sunt similare celor determinate in absenta pulverizarii, in prima etapa de cercetare experimentala.
11. Rezultatele obtinute in cea de a doua etapa de experimentare confirma faptul ca pulverizarea apei in curentul de aer permite cresterea semnificativa a coeficientului global de schimb de caldura (de 1,5 1,65 ori pentru aceeasi valoare a lui waer). Astfel, daca pe langa cresterea vitezei aerului la 24,55 m/s se apeleaza si la pulverizarea apei in curentul de aer, se ajunge la situatia in care krad = 221,51 W/m2∙K si consumul de apa este 3,64∙10-3 kg/s = 13,1 kg/h. In acest caz, SR sp = 0,135 m2/kW; pentru un grup binar de 100 kW rezulta reduceri ale masei si volumului radiatorului pana la 52,7 kg / 0,12 m3. Asadar, pulverizarea permite reducerea masei si a gabaritului radiatorului cu aproximativ 40%. Metoda presupune, insa, un consum de apa care nu poate fi neglijat. Daca se are in vedere o autonomie de 5 ore de functionare a grupului binar, se impune utilizarea unui rezervor de apa de aproximativ 65 litri, care anuleaza avantajul datorat reducerii masei totale a radiatorului (masa metalica + apa continuta in radiator) si chiar introduce un spor de greutate de 36,3 kg. Pe ansamblu, metoda este avantajoasa atat din punct de vedere economic - prin economia de material realizata - cat si, mai ales, din punct de vedere al incadrarii in limitele de gabarit admise pentru sistemul de propulsie al unui automobil: avantajul reducerii considerabile a gabaritului radiatorului este mult mai important decat neajunsul pe care il reprezinta prezenta rezervorului de apa de adaos, amplasarea acestuia pe automobil fiind mult mai putin pretentioasa si mai usor de realizat. De altfel, consumul de apa constituie un neajuns doar in prima jumatate a intervalului de timp care reprezinta autonomia automobilului: pe masura ce automobilul se deplaseaza, masa sa se reduce treptat pe masura ce apa este pulverizata; dupa pulverizarea a 36,3 kg de apa (diferenta masica intre variantele fara/cu pulverizarea apei), automobilul devine la fel de greu ca si in cazul in care pulverizarea lipseste. Din acest punct si pana la golirea rezervorului de apa avantajul masic al variantei cu pulverizarea apei este din ce in ce mai important.
12. Varianta de instalatie de condensare cu pulverizarea apei in curentul de aer de racire nu poate fi adoptata in cazul in care automobilul este exploatat in zone cu resurse limitate de apa; in acest caz, solutia optima pentru reducerea gabaritului radiatorului este cresterea vitezei aerului prin radiator. Aceasta crestere este insa limitata de valorile minime ale randamentului global care justifica utilizarea grupurilor binare in locul motoarelor diesel.
Concluziile prezentate mai sus certifica aprecierea ca, prin activitatea desfasurata si prin rezultatele obtinute, au fost respectate si indeplinite obiectivele propuse la inceperea activitatii de cercetare in cadrul grantului.
PERSPECTIVE DE DEZVOLTARE A GRUPURILOR COGENERATIVE BINARE DE TRANSPORT
In toate estimarile realizate in cadrul grantului, pentru parametrii grupului binar care influenteaza randamentul si consumul specific de combustibil (temperatura dupa camera de ardere a TMG - TT; parametrii aburului viu - pabM, TabM; randamentul comprimarii - ηc; randamentele interioare ale turbinelor cu gaze - hi1 hi2; randamentul interior al turbinei cu abur - hab), nu au fost adoptate valori care sa conduca la dificultati tehnologice. Se prognozeaza ca valorile limita actuale ale acestor parametrii vor creste in viitorul apropiat, datorita materialelor de ultima generatie aparute si a noilor tehnologii de fabricatie. Aceasta va conduce la cresterea randamentului global η si la reducerea corespunzatoare a consumului specific de combustibil, cu efecte semnificative asupra avantajelor economice. Asfel, in cazul in care si (valorile tuturor celorlalti parametrii fiind neschimbate), valorile indicatorilor de performanta corespunzatori schemei termodinamice 6 se imbunatatesc considerabil. Asa cum s-a aratat in subcap. 3.6, atat performantele maxime ale grupului binar cat si gabaritul minim al radiatorului instalatiei de condensare se obtin atunci cand tai, ε si Δtrec au valorile minime luate in calcul, anume , si . In acest caz, in punctul de perfomante maxime se obtin valorile:
; ;
; .
Se observa ca, pe langa imbunatatirea performantelor (cresterea randamentului cu 3,3% si reducerea consumului de combustibil cu 15 g/kWh), rezulta si reducerea gabaritului radiatorului instalatiei de condensare (SR sp se reduce cu 14%), cu efecte directe asupra gabaritului intregului grup binar.
1. Abbot J.W., Baham G.J., COGAS - a new look at Naval Propulsion, Naval Engineering, October 1974
2. Aldea M. s.a., Cazane de abur si recipiente sub presiune - indrumator, Ed. Tehnica, Bucuresti, 1972
3. Alexe F., Cenusa V., Marculescu C., Ene M., Optimizarea centralelor electrice de medie putere cu ciclu combinat abur-gaze, Producerea, transportul si distributia energiei electrice si termice, nr.10, 2000
4. Artenov G.A, Boikov V.P., Lisov V.I., Sudovie gazoturbinie ustanovki, L. Sudostroenie, 1978
5. Bajenov G.V., Romanov V.I., Lisov V.I., Sudovaia gazoturbinnaia ustanovka, M - 25, Movskoi flot, 1976 nr. 9
6. Bannister R.L., Newby R.A., Yang W.C., Final report on the Development of a Hydrogen-Fueled Combustion Turbine Cycle for Power generation, Transactions of the ASME, vol. 121, January 1999
7. Balanescu D.T., Dragomir-Stanciu D., Possibilities to Use the Binary Cogenerative Units as Terrestrial Propulsion Systems, Proceedings of the 35th International Scientific Symposium of METRA, Vol. II, Bucuresti, 27-28 mai 2004
8. Balanescu D.T., Homutescu M.V., Small Scale Combined Cycle Units - Clean Power Systems for Terrestrial Propulsion, Proceedings of the 3rd International Conference on Electrical and Power Engineering, Tom L (LIV) Fasc. 5C, Iasi, oct. 2004
9. Balanescu D.T., Dragomir-Stanciu D., The Performances of a Small Scale Combined Cycle Unit with Low Water Consumption, XXXVI. Kraftwerkstechnisches Kolloquium, Tagungsband II, Dresden, 19- 20 oct. 2004
10. Balanescu D.T., Dragomir-Stanciu D., Binary Cogenerative Units for Terrestrial Propulsion: A Theoretical estimation of the Performances and Size, Proceedings of the 5th International Colloquium "FUELS", Esslingen, ian. 2005
11. Balanescu D.T., Analysis of the Gas Turbine Inlet Temperature Influence over the Performances of a Small Scale Combined Cycle Mobile Unit Based on a Gas Turboengine with Recuperator, Proceedings of the 1st International Conference on Thermal Engines and Environmental Engineering, Fasc. IV, Galati, iunie 2005
12. Balanescu D.T., The Influence of the Temperature Difference on Recuperator over the Performances of a Small Scale Combined Cycle Unit for Terrestrial Propulsion, Proceedings of the 1st International Conference on Thermal Engines and Environmental Engineering, Fasc. IV, Galati, iunie 2005
13. Balanescu D.T., Manolache G., Homutescu V.M., Small Scale Combined Cycle Mobile Unit with Postcombustion Chamber and Based on a Gas Turboengine with Heat Exchanger: Performance Estimation, Proceedings of the 5th International Conference on Electromechanical and Power Systems - SIELMEN 2005, Vol. I, Chisinau, oct. 2005.
14. Balanescu D.T., Homutescu V.M., Small Scale Combined Cycle Mobile Unit with postcombustion and Based on a Regenerative Gas Cycle, Proceedings of the International Conference Trans & MOTAUTO'05+, Veliko-Tarnovo, nov. 2005
15. Bettocchi R., Cantore G., Montenegro G.N., A. Peretto, Performance Optimization In Hybrid Plants with Two Pressure Levels Combined Cycles, Integrated with Geotermal Sources, Florence World Energy Research Symposium, 1994
16. Boissenin Y., Combined Cycle Power Plants:A Practical Guide to the Right Choice, 3rd International Symposium on Turbomachinery, Combined-Cycles Technologies and Cogeneration, Nisa, 1989
17. Boland O., A Comparative Evaluation Of Advanced Combined Cycle Alternatives, ASME Journal of Engineering for Gas Turbines and Power, vol. 113, pp 190-197
18. Bong D., Combined Cycle Advantages, http://www.visionengineer.com/mech/ccadvantages.shtml
19. Burducea C. s a., Centrale nuclearoelectrice de putere mare, Ed. Tehnica, Bucuresti, 1974
20. Cantemir C.G., Hubert C., Rizzoni G., Ursescu G., High-power High-speed Road Train System, International Truck and Bus Meeting and Exhibition, SAE 2003-01-3380, Ft. Worth, Texas, November 2003
21. Cantuniar C., Turbomasini termice, Ed. Matrix Rom, Bucuresti, 1998
22. Carabogdan G. s.a., Instalatii termice industriale, Ed. Tehnica, Bucuresti, 1978
23. Catina J.L., Fortune H.J., Soroka G.E., Repowering Chesterfield 1 and 2 with combined cycle, Transactions of ASME, 110, 214-299, 1988
24. Carlanescu C., Ursescu D., Manea I., Turbomotoare de aviatie - Aplicatii industriale, Ed. Didactica si Pedagogica, Bucuresti, 1997
25. Cerri G., Parametric Analysis of Combined Gas, Steam Cycles, Journal of Engineering for Gas Turbines and Power, January 1987
26. Chiesa P., Consonni S., Shift Reactors and Physical Absorption for Low - CO2 Emissions IGCC, ASME Journal of Engineering for Gas Turbines and Power, vol. 121, April 1999
27. Cohen H., Rogers G.F.C., Saravanamuttoo H.I.H., Gas Turbine Theory, 4th Edition, Longman Group Limited 1996
28. Consonni S., Larson E.D., Biomass-Gasifier/Aeroderivative Gas Turbine Combined Cycles. Part B: Performance Calculations and Economic Assessment, ASME Journal of Engineering for Gas Turbines and Power, vol. 118, 1996
29. Costin D., Contributii la studiul schimbatoarelor de caldura cu pelicula de lichid, Teza de doctorat, Iasi, 1980
30. Creta G., Curs de turbine cu abur, vol. I, Ed. I.P.Timisoara, 1975
31. Creta G., Turbine cu abur si gaze, Ed. Tehnica, Bucuresti, 1996
32. Davidescu A., Muncica H., Schimbul de caldura in instalatiile industriale, Ed.Tehnica, Bucuresti, 1971
33. Dechamps P.J., Advanced Combined Cycle Alternatives with the Latest Gas Turbines, ASME Journal of Engineering for Gas Turbines and Power, vol. 120, January 1998
34. DeMoss T.B., They're Here (almost): The 60% Efficient Combined Cycle, Power Engineering, July 1996
35. Dinculescu C. s.a., Analiza termodinamica a schemelor centralelor electrice, Ed. Acad. RSR, Bucuresti, 1967
36. Finckh H.H., Pfost H., Development Potential of Combined-Cycle (GUD) Power Plants with and without Supplemetary Firing, Journal of Engineering for Gas Turbines and Power, October 1992
37. Frutschi H. U., Control Methods for Cogeneration with Gas Turbines and Combined Cycles, International Symposium of Turbomachinery, Combined-Cycles Technologies and Cogeneration, Montreaux, 1987
38. Gheorghiu St., Cazane de abur, Ed. Didactica si Pedagogica, Bucuresti, 1966
39. Giurca V., Bazele cercetarii experimentale a masinilor termice, vol. I, II, Ed. Rotaprint, Iasi, 1992
40. Grecu T., Turbine si turbocompresoare, vol. I, Ed. Didactica si Pedagogica, Bucuresti, 1965
41. Grecu T., Cirdu M., Nicolau I., Turbine cu abur, Ed. Tehnica, Bucuresti, 1976
42. Gurevici A.M, Gazoturbinnii naduv sudovih parogeneratorov, Ed. Sudostroenie, 1978
43. Gusso R., Pucci M., Sabella D., Combined Cycles with or without Cogeneration for Small Plants, 3rd International Symposium of Turbomachinery, Combined-Cycles Technologies and Cogeneration, Nisa, 1989
44. Horlock J.H., The Evaporative Gas Turbine [EGT] Cycle, Transactions of ASME, vol. 120, April 1998.
45. Iordache N., Stefan I., Vasilescu M., Generatoare de abur si apa fierbinte recuperatoare, Oficiul de informare documentara pentru aprovizionarea tehnico-materiala si controlul gospodaririi fondurilor fixe, Bucuresti, 1988
46. Ishigai S., Steam Power Engineering - Thermal and Hydraulic Design Principles, Cambridge University Press 1999
47. Korpela T., Characterisation of the Combined Cycle Steam Process, VGB Power Tech, decembrie 2002
48. Korpela T., Combined Cycle Performance Prediction, VGB Power Tech, ianuarie 2003
49. Kuehn E., Combined Cycles Leads Efficiency Race, Power Engineering, May 1995
50. Kurzon A.G., Iodovin B.S., Sudovie kombinirovannie energetiskie ustanovki, Ed. "Sudostroenie", Leningrad, 1981
51. Kurzon A.G., Gazoturbinie ustanovski marskih sudov, Ed. Transport 1967
52. Leca A. s.a., Procese si instalatii termice in centrale nuclearoelectrice, Ed. Didactica si Pedagogica, Bucuresti, 1979
53. Lisov V.I., Glavnii gazoturbinnii agregat sudna "Capitan Smirnov", Sudostroenie, 1980, nr. 3
54. Lozza G., Bottoming Steam Cycles for Combined Gas-Steam Power Plants: A Theoretical Estimation of Steam Turbine Performance and Cycle Analysis, Proceedings, 1990 ASME Cogen-Turbo, ASME, New York, pp. 83-92
55. Mihaescu L., Prisecaru T., Oprea I., Cazane si Turbine. Notiuni de baza, Ed. Matrix Rom, Bucuresti,1999
56. Miler M.C., Instalatii cogenerative cu componente recuperate, teza de doctorat, Iasi, 2004
57. Molinari G., Calapso F., La fornitura di potenza nelle navi da crociera. Aparati Diesel-elettrici o Coges?, La Termotecnica, Luglio/Agosto 2004
58. Montenegro G.N., Bianchi M., Peretto A., Sistemi Energetici e Loro Componenti - Considerazioni Teoriche e Valutazioni Numeriche, Ed. Pitagora, Bologna, 2004
59. Neaga C., Tratat de generatoare de abur, vol. I, Ed. Agir, Bucuresti, 2001
60. 85)Neaga C., Epure Al., Indrumar - Calculul termic al generatoarelor de abur, Ed. Tehnica, Bucuresti, 1988
61. Pavlov K.F., Procese si aparate in ingineria chimica (exercitii si probleme), Ed. Tehnica, Bucuresti, 1981
62. Pop M.G. s.a., Indrumar - Tabele, nomograme si formule termotehnice, vol. I, II, III, Ed. Tehnica, Bucuresti, 1987
63. Poulton M.L., Alternative Enginees for Road Vehicles, Computational Mechanics Publications, 1994
64. Pimsner V. s.a., Energetica turbomotoarelor cu ardere interna, Ed. Academiei RSR, Bucuresti, 1964
65. Rice I.G., The Combined Reheat Gas Turbine/Steam Turbine Cycle, Part I: A Critical Analysis of the Combined reheat Gas Turbine/Steam Turbine Cycle, Journal of Engineering for Power, January 1980
66. Rice I.G., The Combined Reheat Gas Turbine/Steam Turbine Cycle, Part II: The LM 5000 Gas Generator Applied to the Combined Reheat Gas Turbine/Steam Turbine Cycle, Journal of Engineering for Power, January 1980
67. Romanov F.V.I., Lisov V.I., Pazvitie sudovih gazoturbinih ustanovok, Sudostroenie 1979 nr. 4
68. Rossiter A.P., Criteria for Integration of Combined Cycle Cogeneration Systems in the Process Industries, Heat Recovery Systems & CHP, no. 1, 1990
69. Ryder G., The Advantages of Combined Cycle Plants: A "New Generation" Technology, November 1997, http://www.probeinternational.org/pi/documents/Mekong/psr2.htm
70. Schröder K. s.a., Centrale termoelectrice de putere mare, Vol. II - Bazele constructiei centralelor termoelectrice, Ed. Tehnica, Bucuresti - 1965
71. Shortlidge C.C., Control System for a 373 kW, Intercooled, Too Spol Gas Turbine Engine Powering a Hybrid Electric World Sports Car Class Vehicle, ASME Journal of Engineering for Gas Turbines and Power, vol. 120, January 1998
72. Smith A.J., McMichael T.C., The Evolution of Combined Cycle Cogeneration Technology, 3-rd International Symposium of Turbomachinery, Combined-Cycles Technologies and Cogeneration, Nisa, 1989
73. Snoeck J., Repowering Vilvoorde, Seminaire d'Energetique TERM, Tractebel Engineering, 19 martie 2004
74. Stodola A., Turbines a vapeur et a gaz, Ed. Dunod, Paris, 1925
75. Stefanescu D., Marinescu M., Termotehnica, Ed. Didactica si Pedagogica, Bucuresti, 1983
76. Töle H., Masurari in instalatiile termice, Ed. Tehnica, Bucuresti, 1972
77. Tsunoda A., Shimoda H., Takaishi T., Mitsubishi Lean-Burn Gas Engine with World's Highest Thermal Efficiency, Mitsubishi Heavy Industries, Ltd. Technical Review, Vol.40 No.4, Aug. 2003
78. Ungureanu C., Generatoare de abur pentru instalatii energetice clasice si nucleare, Ed. Didactica si Pedagogica, Bucuresti, 1977
79. Ursescu D., Tardea I., Turbomotoare cu abur si gaze, vol. I - Procese, Ed. I.P. Iasi, 1981
80. Ursescu D., Jugureanu E., Tardea I., Instalatie binara gaze - abur, Creatia tehnica si fiabilitatea in constructia de masini, Iasi, 22 - 23 decembrie 1978
81. Ursescu G., Contributii la utilizarea recuperatoarelor de caldura pe turbomotoarele cu gaze de aviatie, transport terestru si naval, teza de doctorat, Iasi, 2003
82. Ursescu D., Ursescu G., Dragomir D., Zubcu V., Perspective de utilizare a grupurilor cogenerative gaze - abur ca sisteme de propulsie pentru automobile, Conferinta Internationala TURBO'98, Bucuresti, iulie 1998
83. Ursescu D., Homutescu C., Poeata N., Stabilirea procedeelor de calcul a ciclurilor celor mai uzuale turbine cu gaze in domeniul temperaturii 800 - 2500 K, Contract de cercetare, beneficiar - I.N.M.T. - Bucuresti, 1996
84. Valenti M., Reaching for 60 Percent, Mechanical Engineering-April 2002
85. Valentine H., A Combined Cycle Locomotive?, April 2001, http://www.internationalsteam.co.uk/trains/newsteam/modern17.htm
86. Vladea I., Tratat de termotehnica tehnica, Ed.Tehnica, Bucuresti, 1979
87. ***, GE Completes H Testing, Power Engineering International, July 2003
88. ***, GE Marine Engines - A Gas Turbine-Based System Powers Newly-Launched Brilliance Of The Seas Cruise Ship, http://www.geae.com/aboutgeae/presscenter/marine/marine_20020731.html
89. ***, http://www.transportal.no/Press/3107ge.html
90. ***, http://www.geae.com/engines/marine/lm2500.html
91. ***, http://www.pratt-whitney.com/prod_ind_st18.asp
92. ***, Proceedings of the Institute of Marine Engineering. Science and Technology, Part A: Journal of Marine Engineering, London, July 2003
93. ***, Sistemul de conversie termodinamica cu turbine cu gaze si cicluri combinate abur-gaze, Energetica (43), iulie - august 1995
94. ***, Top Plants - Baglan Bay Power Station, Cardiff, Walles, Power, iulie/august 2003
95. ***, Turbine Renewal: Shaping an Emerging Gas-Fired Power Source, Oak Ridge National Laboratory Review, Vol.33, No.1, 2000