You also want an ePaper? Increase the reach of your titles
YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.
2 <strong>Teoretične</strong> <strong>osnove</strong><br />
Brez obširnejše razlage so nanizani glavni fizikalni zakoni iz mehanike tekočin in<br />
termodinamike, ki so potrebni za razumevanje učne snovi v zvezi z energetskimi<br />
stroji in napravami. Na koncu poglavja je dodanih še nekaj osnov o podobnostnih<br />
kriterijih.<br />
2.1 Mehanika tekočin<br />
2.1.1 Gostota<br />
Plini in kapljevine se med seboj razlikujejo po gostoti, stisljivosti itn., imajo pa sicer<br />
mnogo skupnih lastnosti. Z eno besedo se imenujejo tekočine ali s tujko fluidi. Velja:<br />
dV<br />
V<br />
ϱ = dm<br />
dV<br />
= dp<br />
E = dϱ<br />
ϱ<br />
ϱ =<br />
ϱ 0<br />
1 + ∆p<br />
E<br />
(2.1.1)<br />
(2.1.2)<br />
(2.1.3)<br />
pri čemer je ∆p / (N/m 2 ) tlačna razlika in E / (N/m 2 ) modul elastičnosti, preglednica<br />
2.1.<br />
Za mnoge praktične primere zadostujejo nekatere poenostavitve, npr. za kapljevine<br />
vzamemo, da so praktično nestisljive, za pline pa, da se vedejo kot idealni plini, gostota<br />
plinov se spreminja po plinski enačbi. Pri majhnih spremembah tlakov pogosto<br />
zadošča, da tudi pline obravnavamo kot nestisljive, npr. tok zraka skozi vetrnice in<br />
ventilatorje.<br />
9
10 2 TEORETIČNE OSNOVE<br />
Tabela 2.1: Modul elastičnosti za nekatere snovi<br />
Snov Modul elastičnosti<br />
T / 0 ◦ C E / (N/cm 2 )<br />
Benzol 1, 2 · 10 5<br />
Voda 2, 1 · 10 5<br />
Živo srebro 2, 9 · 10 6<br />
Jeklo 2, 1 · 10 7<br />
2.1.2 Hidrostatični tlak<br />
Kapljevina lahko sprejema samo tlačne sile. Če zanemarimo njeno težo, deluje hidrostatični<br />
tlak v notranjosti kapljevine enakomerno na vse strani (Pascalov zakon).<br />
Dejansko je v večini tehničnih primerov potrebno upoštevati težo kapljevine, hidrostatični<br />
tlak se namreč povečuje sorazmerno z globino:<br />
p = p 0 + ϱ · g · ∆H (2.1.4)<br />
pri tem je p 0 tlak okolice in ϱ gostota kapljevine.<br />
2.1.3 Vzgon<br />
Na telo, ki je potopljeno v tekočino, deluje sila vzgona. Ta sila vzgona je enaka teži<br />
izpodrinjene tekočine. Če telo plava, potem je sila vzgona v ravnotežju s silo teže:<br />
teža plavajočega telesa je enaka teži izpodrinjene tekočine (Arhimedov zakon):<br />
F V zg = ϱ · g · V (2.1.5)<br />
pri tem je V prostornina izpodrinjene tekočine.<br />
2.1.4 Površinska napetost<br />
Če povečujemo površino neke tekočine, npr. milnega mehurčka, je treba dovajati<br />
delo:<br />
∆W σ = σ · ∆A (2.1.6)<br />
To delo je sorazmerno povečanju površine mehurčka. Sorazmernostni faktor σ / (N/m<br />
= J/m 2 ) se imenuje površinska napetost in je pri 20 ◦ C za vodo σ = 72, 75 · 10 −3<br />
N/m, za veliko večino organskih tekočin pa σ = (20–40) · 10 −3 N/m.
2.1 MEHANIKA TEKOČIN 11<br />
Zgled. Sile v mehurčku<br />
Zaradi površinske napetosti je v mehurčku večji tlak, kot pa v okolici. Za kroglasti<br />
mehurček, ki je potopljen v kapljevini, je to mogoče enostavno izračunati: če<br />
prerežemo mehurček s premerom d, potem mora biti sila zaradi površinske napetosti<br />
po obsegu kroga (ϱ · π · d) v ravnotežju s tlačno silo, ki deluje na površino kroga<br />
(∆p σ · π · d 2 /4) :<br />
σ · π · d = ∆p σ · π · d2<br />
4<br />
V urejeni obliki:<br />
∆p σ = 4 · σ<br />
d<br />
(2.1.7)<br />
(2.1.8)<br />
Povečanje tlaka ∆p σ v mehurčku je premo sorazmerno s površinsko napetostjo σ in<br />
obratno sorazmerno s premerom mehurčka d.<br />
2.1.5 Kontinuitetna enačba<br />
Za stacionarni tok tekočine velja:<br />
∫<br />
∫<br />
ṁ 1 = ṁ 2 = ϱ 1 · v m1 · dA 1 =<br />
A 1<br />
ϱ 2 · v m2 · dA 2<br />
A 2<br />
(2.1.9)<br />
pri tem je ṁ masni tok, ϱ gostota, v m srednja (krajevno povprečena) hitrost tekočine<br />
in A prerez toka tekočine. Pri majhnih prerezih in zanemarljivih krajevnih spremembah<br />
hitrosti velja poenostavljen zapis kontinuitetne enačbe:<br />
ṁ 1 = ṁ 2 = ϱ 1 · v m1 · A 1 = ϱ 2 · v m2 · A 2 (2.1.10)<br />
2.1.6 Termična enačba stanja<br />
Veličine stanja p, T in ϱ neke tekočine so med seboj odvisne, povezuje jih termična<br />
enačba stanja. Za realne tekočine je lahko zveza ϱ = ϱ(p, T ) zelo komplicirana, v<br />
večini primerov pa so mogoče poenostavitve. Za idealni plin velja plinska enačba:<br />
p<br />
ϱ · R · T = 1 (2.1.11)<br />
Za realne pline se pogosto uporablja korigirana plinska enačba:<br />
p<br />
ϱ · R · T = ξ (2.1.12)
12 2 TEORETIČNE OSNOVE<br />
Faktor stisljivosti ξ je določen eksperimentalno, odvisen je od snovi, njene temperature<br />
in tlaka. Navadno je napisan v obliki polinoma, ki je primeren za računalniško<br />
obdelavo. Za praktični izračun plinov in par, npr. vodne pare, se uporabljajo tudi tabele<br />
in grafični diagrami. Lastnosti realnih plinov se razlikujejo od idealnih tem bolj,<br />
čim višji je tlak in čim nižja je temperatura, drugače napisano: čim večja je gostota<br />
plina. Za vodo in druge nestisljive snovi smemo predpostaviti:<br />
ϱ = konst. (2.1.13)<br />
2.1.7 Energijski izrek<br />
Energijski izrek je za mehaniko tekočin mogoče izpeljati iz osnovnega Newtonovega<br />
zakona:<br />
F =<br />
d(m · v)<br />
dt<br />
(2.1.14)<br />
Obe strani enačbe množimo s prirastkom dolžine dL, katere smer ustreza smeri sile<br />
F , in predpostavimo m = konst.:<br />
F · dL = m · dv · dL (2.1.15)<br />
dt<br />
Na levi strani dobimo delo dW = F · dL, na desni strani izraz, ki pomeni prirastek<br />
kinetične energije: m · (dL/dt) · dv = m · v · dv. Od tod sledi energijski izrek<br />
mehanike, ki pravi, da je delo enako razliki kinetične energije.<br />
dW = m · v · dv (2.1.16)<br />
V integralni obliki:<br />
∆W = m ( )<br />
2 · v2 2 − v1<br />
2<br />
(2.1.17)<br />
Če zanemarimo sile med posameznimi atomi, velja energijski izrek tudi za skupek<br />
atomov neke snovi, pri tem je kinetična energija vsota kinetičnih energij vseh atomov.<br />
To kinetično energijo sestavljajo:<br />
• energija, ki jo imajo atomi zaradi svojega neurejenega gibanja (pri idealnih<br />
plinih ustreza to njihovi notranji energiji dU) in<br />
• energija, ki ustreza hitrosti težišča celotne plinske mase.
2.1 MEHANIKA TEKOČIN 13<br />
V splošnem ni mogoče določiti, kolikšen del dovedenega dela povečuje notranjo<br />
energijo in kolikšen del makroskopsko kinetično energijo. V mnogih primerih se<br />
vse dovedeno delo porablja samo za povečevanje makroskopske kinetične energije.<br />
Natančno velja to le za nestisljive idealne snovi brez trenja. Ustrezna oblika energijskega<br />
izreka za ta posebni primer je Bernoullijeva enačba:<br />
m · p1 + m · g · H 1 + m · v2 1<br />
ϱ 1 2 = m · p2 + m · g · H 2 + m · v2 2<br />
ϱ 2 2<br />
(2.1.18)<br />
Enačba pove, da se je energija, ki rezultira iz tlaka tekočine in njene zemeljske privlačnosti,<br />
brez izgub spremenila v kinetično energijo.<br />
Zgled. Hitrost iztoka kapljevine iz posode<br />
S kakšno hitrostjo izteka olje iz narisanega cevovoda, če zanemarimo trenje? Razlika<br />
višin ∆H = 6 m, slika 2.1. Tlak v posodi je enak tlaku na iztoku p 1 = p 2 , začetna<br />
hitrost (hitrost nižanja gladine) je zanemarljiva: v 1 = 0. Iz Bernoullijeve enačbe<br />
dobimo:<br />
ϱ · g · ∆H = ϱ 2 · v2 2<br />
Od tod je mogoče izračunati izstopno hitrost olja:<br />
v 2 = √ 2 · g · ∆H = 10,6 m/s<br />
Hitrost je enaka, kot če bi olje prosto padalo s te višine.<br />
Slika 2.1: Iztok kapljevine iz posode
14 2 TEORETIČNE OSNOVE<br />
Tabela 2.2: Primerjava energijskega in impulznega izreka<br />
Skalarni energijski izrek<br />
Newtonov zakon: F = m · dv<br />
dt<br />
dL<br />
F · dL = m · dv<br />
dt · dL<br />
dW = m · dL<br />
dt · dv = m · v · dv<br />
∆W = m 2 · (v<br />
2 2 − )<br />
v2 1<br />
Vektorski impulzni izrek<br />
dt<br />
F · dt = m · dv<br />
dI = m · dv<br />
∆I = m · (v 2 − v 1 )<br />
∆ ˙ I = ṁ · (v 2 − v 1 ) = F<br />
2.1.8 Impulzni izrek<br />
Tudi impulzni izrek za mehaniko tekočin je mogoče izpeljati iz osnovnega Newtonovega<br />
zakona, enačba 2.1.14. Obe strani množimo s prirastkom časa dt in predpostavimo<br />
m = konst.:<br />
F · dt = m · dv (2.1.19)<br />
Na levi strani dobimo impulz sile dI = F · dt, na desni strani izraz, ki pomeni<br />
prirastek hitrosti m · dv. Od tod sledi impulzni izrek, ki pravi, da je impulz sile enak<br />
spremembi gibalne količine:<br />
dI = m · dv (2.1.20)<br />
V tehniki imamo navadno opravka z masnim tokom ṁ in ne z maso m; za stacionarni<br />
tok velja v integralni obliki:<br />
∆ ˙ I = ṁ · (v 2 − v 1 ) = F (2.1.21)<br />
Tok tekočine spreminja pod vplivom zunanje sile F hitrost za ∆v. Ker so sile in<br />
hitrosti vektorji, je smer sile enaka smeri spremembe hitrosti. Čeprav je impulz vektorska<br />
veličina in imamo zato na razpolago tri enačbe, je te enačbe navadno lažje<br />
uporabljati kot pa skalarni energijskih izrek, preglednica 2.2.<br />
Impulzni izrek, uporabljen za tekočino v kanalu med lopatjem, nam da Eulerjevo<br />
enačbo, ki je podrobneje obrazložena pri turbinskih strojih.<br />
Stanje tekočine je v splošnem opisano z gostoto ϱ, temperaturo T , tlakom p in vektorjem<br />
hitrosti ⃗v(v x , v y , v z ) . Če poznamo snovne lastnosti tekočine (dinamična viskoznost<br />
η, toplotna prevodnost λ, izobarna in izohorna specifična toplota c p in c v
2.1 MEHANIKA TEKOČIN 15<br />
itd.), je v splošnem na razpolago dovolj enačb za določitev šestih neznank: ϱ, T ,<br />
p, ⃗v(v x , v y , v z ). Prav tako imamo tudi za reševanje problemov na razpolago šest<br />
enačb: kontinuitetno enačbo, termično enačbo stanja, energijski izrek in tri enačbe<br />
impulznega izreka. Načelno je tako mogoče za dane razmere (znana geometrija,<br />
sile in robni pogoji) določiti vsako stanje tekočine. Praktične težave, ki nastanejo<br />
pri računanju, pa so matematične narave: Navier-Stokesove enačbe (kontinuietna,<br />
gibalna in energijska enačba v splošni obliki) so namreč nelinearne parcialne diferencialne<br />
enačbe, ki so analitično nerešljive.<br />
2.1.9 Trenje v toku<br />
Doslej še ni bilo omenjeno, da moramo pri tekočinah upoštevati tudi notranje trenje,<br />
to je, da imamo v tekočini strižne napetosti τ, ki so odvisne od hitrosti te tekočine.<br />
Če položimo na mokra tla ploščo s ploščino A, lahko ugotovimo, da potrebujemo<br />
za premikanje plošče neko silo F , slika 2.2. Ta sila je odvisna od ploščine A, od<br />
hitrosti pomika v, od razdalje med ploščo in tlemi B in od sorazmernostnega faktorja<br />
η. Sorazmernostni faktor je specifičen za vsako tekočino in ga imenujemo dinamična<br />
viskoznost:<br />
F = η · A · v<br />
B<br />
(2.1.22)<br />
splošneje:<br />
F<br />
A = τ = η · dv<br />
dy<br />
(2.1.23)<br />
Da izločimo vpliv mase tekočine, delimo dinamično viskoznost z gostoto. To veličino<br />
imenujemo kinematična viskoznost: ν = η/ϱ.<br />
Slika 2.2: Strižne sile v tekočini
16 2 TEORETIČNE OSNOVE<br />
laminarni tok<br />
turbulentni tok<br />
Slika 2.3: Hitrostni profil v cevi okroglega prereza<br />
V osnovi imamo opravka z dvema vrstama toka. Laminarni tok imenujemo tok, kjer<br />
potekajo tokovnice po vsem tokovnem prerezu vzporedno. Menjava energije med<br />
tokovnicami je majhna in je izključno posledica trenja tekočine. V tehniki imamo<br />
večinoma opravka s turbulentnim tokom: delci tekočine zadevajo drug ob drugega,<br />
zaradi tega nastane nepravilno gibanje v vzdolžni in prečni smeri toka. Menjava<br />
energije v tekočini je večja kot pri laminarnem toku in je posledica trenja in medsebojnih<br />
elastičnih trkov delcev tekočine. Zaradi večje izmenjave energije znotraj toka<br />
je hitrost toka po celotnem prerezu zelo enakomerna, z izjemo tanke mejne plasti ob<br />
steni δ, slika 2.3. S kakšnim tokom imamo opravka, nam pove Reynoldsovo število<br />
Re, ki je eden od brezdimenzijskih kriterijev podobnosti; tok v okroglem cevovodu<br />
je npr. turbulenten, če je Re > 2300. Več o tem je napisano v poglavju o zakonih<br />
podobnosti.<br />
2.1.10 Tlačne izgube v cevovodih in armaturah<br />
Pri toku tekočin skozi cevi, odcepe, armature itd. pride zaradi notranjega trenja v<br />
tekočini do izgube tlaka. Za ravno cev velja enačba:<br />
∆p = λ · L<br />
d · ϱ · v2<br />
2<br />
(2.1.24)<br />
Izguba specifične tlačne energije ∆p/ϱ je premo sorazmerna z dolžino cevovoda L,<br />
s kinetično energijo v 2 /2 in obratno sorazmerna s premerom cevovoda d. Sorazmernostni<br />
faktor λ se imenuje koeficient tekočinskega trenja in ga je mogoče za la-
2.1 MEHANIKA TEKOČIN 17<br />
minarni tok izpeljati matematično iz Navier-Stokesove diferencialne enačbe, Hagen-<br />
Poiseuillov zakon:<br />
λ = 64<br />
(2.1.25)<br />
Re<br />
kjer je Re Reynoldsovo število.<br />
Za turbuletni tok pa je koeficient tekočinskega trenja λ mogoče določiti samo eksperimentalno,<br />
saj se ga zaradi naključne narave turbulence ne da analitično določiti.<br />
Za praktično uporabo je primerna Colebrookova poenostavitev Prandtlove enačbe, ki<br />
velja za celotno turbulentno področje in za hidravlično gladke cevi. Koeficient trenja<br />
je odvisen samo od Reynoldsovega števila:<br />
λ = (<br />
0,309 ) 2<br />
(2.1.26)<br />
lg Re<br />
7<br />
V splošnem je koeficient tekočinskega trenja odvisen tudi od hrapavosti stene. Dobre<br />
rezultate daje zopet Colebrookova enačba z upoštevanjem absolutne hrapavosti<br />
stene k:<br />
1<br />
λ = [ ( )]<br />
2 · lg 2,51<br />
2<br />
(2.1.27)<br />
Re·√λ + 0,27 · k<br />
d<br />
Vrednosti za koeficient tekočinskega trenja λ so navadno prikazane grafično v Moodyjevem<br />
diagramu, slika 2.4. Diagram prikazuje koeficient trenja za različna Reynoldsova<br />
števila za gladke in za hrapave cevi. Namesto absolutne hrapavosti k je v<br />
diagramu vrisana brezdimenzijska relativna hrapavost k/d. Pri tem je merilo hidravlično<br />
gladke cevi določeno z debelino mejne plasti δ in absolutno hrapavostjo stene<br />
k. Absolutna hrapavost je za neko cev konstanta, medtem ko se debelina mejne plasti<br />
zmanjšuje z naraščanjem Reynoldsovega števila: neka cev se pri majhnih pretokih<br />
vede kot hidravlično gladka, pri velikih pa kot hrapava. Vrednosti srednjih absolutnih<br />
hrapavosti k za najpogosteje uporabljene cevi so zbrane na sliki 2.5.<br />
Izgubo tlaka v kolenih, lokih, odcepih, iztokih, priključkih, ventilih, zasunih itd.<br />
lahko določimo samo z meritvami. Podobno kot pri enačbi za ravne cevi velja:<br />
∆p = ζ · ϱ · v2<br />
(2.1.28)<br />
2<br />
Pri tem je koeficient izgub ζ navadno tabeliran v priročnikih za posamezne vrste<br />
armatur. Za gladke okrogle cevi dobimo medsebojno zvezo:<br />
ζ = λ · L<br />
d<br />
(2.1.29)
18 2 TEORETIČNE OSNOVE<br />
Slika 2.4: Moodyjev diagram; koeficient tekočinskega trenja λ za ravne okrogle cevi<br />
v odvisnosti od Reynoldsovega števila Re in relativne hrapavosti k/d
2.1 MEHANIKA TEKOČIN 19<br />
Slika 2.5: Srednje absolutne hrapavosti k najpogosteje uporabljanih cevi<br />
Zgled. Tlačne izgube v cevovodu<br />
Kakšne so tlačne izgube v hidravlično gladkem cevovodu iz jeklene valjane cevi zunanjega<br />
premera 323,9 mm in debeline stene 10,0 mm ter dolžine 260 m, če teče<br />
skozi cev 780 m 3 /h vode pri 20 ◦ C, dinamična viskoznost vode η = 10 3 kg/(m · s)?<br />
Predhodno je treba izračunati pomožne veličine:<br />
Tok vode ˙V = 0,217 m 3 /s<br />
Notranji premer cevovoda d = 0,3039 m<br />
Prerez cevovoda A = 0,0725 m 2<br />
Hitrost vode v cevovodu v = 2,99 m/s<br />
Reynoldsovo število Re = 907 ′ 760<br />
Koeficient trenja λ = 0,0118<br />
Izguba tlaka v cevovodu:<br />
∆p = λ · L<br />
d · ϱ · v2<br />
2<br />
= 45 ′ 115 Pa = 0,45 bar
20 2 TEORETIČNE OSNOVE<br />
2.2 Termodinamika in prenos toplote<br />
2.2.1 Zakon o ohranitvi energije<br />
Termodinamika razširja mehaniko tekočin z uvedbo novih oblik energije. V nekem<br />
določenem in izoliranem sistemu je vsota vseh energij konstantna; napisati smemo,<br />
da ostane v sistemu konstantna tudi vsota vseh eksergij in anergij. Če se torej v<br />
sistemu poveča energija ene vrste, se mora zaradi tega v tem istem sistemu zmanjšati<br />
energija neke druge vrste. Pod imenom sistem je treba razumeti prostor ali pa količino<br />
snovi, za katero si lahko predstavljamo, da je omejena s stenami. Vse, kar je zunaj teh<br />
dejanskih ali pa namišljenih sten, je okolica. Tako lahko v splošnem ločimo različne<br />
termodinamične sisteme:<br />
• Odprti sistemi: preko meje sistema poteka izmenjava snovi in energije (primer:<br />
motor z notranjim zgorevanjem).<br />
• Energijsko odprti (diatermični) sistemi: preko meje sistema je možna samo<br />
izmenjava energije, ne pa snovi (primer: segrevanje snovi v zaprti posodi).<br />
• Zaprti sistemi: preko meje sistema ne prehaja niti snov niti energija (primer:<br />
adiabatna kompresija v zaprti izolirani posodi).<br />
Delo in toplota sta prehodni energiji, ki delujeta samo na meji med sistemom in<br />
okolico in se uporabljata za prenos energije; pravimo, da sistem oddaja ali sprejema<br />
delo in toploto. Razlika med delom in toploto se opazi samo v okolici, ne pa tudi v<br />
sistemu samem, slika 2.6. Kadar npr. sistem opravlja delo, se to izraža v sočasnem in<br />
urejenem gibanju vseh delcev v okolici (npr. translatorno gibanje bata). Pri prenosu<br />
toplote preko meje sistema pa se to izraža kot neurejeno, termično gibanje delcev v<br />
okolici (npr. termično gibanje molekul tekočine v okolici).<br />
V tehniki imamo največkrat opravka z odprtimi sistemi, kjer v sistem periodično ali<br />
zvezno vstopa in izstopa delovna snov. Če upoštevamo vse vrste energij (tlačno,<br />
kinetično in potencialno) ter masno bilanco, potem lahko za splošni, odprti, neadiabatni,<br />
stacionarno delujoč sistem zapišemo:<br />
U 1 + p 1 · V 1 + m · v2 1<br />
2 + m · g · H 1 + Q 12 =<br />
= U 2 + p 2 · V 2 + m · v2 2<br />
2 + m · g · H 2 + W t12 (2.2.1)<br />
To je oblika energijske enačbe, ki se uporablja pri energetskih strojih, pri čemer se<br />
indeks 1 nanaša na začetek procesa in indeks 2 na konec procesa, ki poteka na meji<br />
sistema. Vsota vseh energij in v sistem dovedene toplote Q 12 je enaka vsoti vseh
2.2 TERMODINAMIKA IN PRENOS TOPLOTE 21<br />
delo<br />
toplota<br />
Slika 2.6: Razlika med delom in toploto; A - urejeno gibanje molekul, B - neurejeno<br />
gibanje molekul<br />
energij in tehničnemu delu W t12 , ki smo ga dobili na izstopu iz sistema, pri čemer<br />
velja: m = ϱ · V . Energije ni mogoče niti proizvesti niti uničiti, ampak samo spremeniti<br />
iz ene oblike v drugo.<br />
Bernoullijeva enačba je posebna oblika energijske enačbe in velja za tok nestisljive<br />
tekočine brez upoštevanja trenja, brez dovoda toplote in odvoda tehničnega dela:<br />
p 1 · V 1 + m · v2 1<br />
2 + m · g · H 1 = p 2 · V 2 + m · v2 2<br />
2 + m · g · H 2 (2.2.2)<br />
Če je sistem namenjen le opravljanju dela, ne pa tudi drugim energijskim pretvorbam,<br />
se splošni zapis prvega glavnega zakona termodinamike, enačba (2.2.1), poenostavi:<br />
Q 12 − W t12 = (U 2 + p 2 · V 2 ) − (U 2 + p 2 · V 2 ) (2.2.3)<br />
kjer sta notranja in tlačna energija (sumanda v vsakem od oklepajev) veličini stanja,<br />
zato jih lahko nadomestimo z novo veličino stanja - entalpijo:<br />
Q 12 − W t12 = m · (h 2 − h 1 ) (2.2.4)<br />
Energijo nekega stacionarno delujočega sistema poveča dovedena toplota Q 12 , zmanjša<br />
pa jo iz sistema pridobljeno tehnično delo W t12 . Pri odprtih, stacionarno delujočih<br />
sistemih, je treba delovno snov ves čas dovajati, iztrošeno snov pa odvajati. Pri tem<br />
potrebujemo stalno polnilno (−p 1 ·V 1 ) in praznilno delo (+p 2 ·V 2 ). Razlika teh dveh<br />
del je volumensko, imenovano tudi absolutno delo:<br />
W u12 = W t12 − p 1 · V 1 + p 2 · V 2 (2.2.5)
22 2 TEORETIČNE OSNOVE<br />
Slika 2.7: Energetski stroji; A in B - odprti sistem: preko meje sistema prehajata<br />
zvezno ali periodično snov in energija, C - zaprti sistem: preko meje sistema prehaja<br />
samo energija<br />
Polnjenje in praznjenje delovne snovi se opravlja neprekinjeno pri turbinskih ali<br />
pretočnih strojih in periodično pri volumenskih ali izrivnih strojih, ki imajo sesalne<br />
in tlačne ventile, slika 2.7.<br />
V zaprtem sistemu tehničnega dela W t12 nimamo, prvi glavni zakon termodinamike<br />
pa dobi obliko:<br />
Q 12 − W u12 = U 2 − U 1 (2.2.6)<br />
Notranjo energijo zaprtega sistema poveča dovedena toplota Q 12 , zmanjša pa jo<br />
iz sistema pridobljeno volumensko delo, imenovano tudi delo enkratne ekspanzije<br />
W u12 . Z enkratno ekspanzijo je ostala delovna snov v energetskem stroju zaprta, njena<br />
uporabnost je izčrpana, s tem pa je izčrpana tudi uporabnost tega stroja, slika 2.7.<br />
2.2.2 Krožni procesi<br />
V tehniki je pomembno pridobivati delo iz toplote stalno in nepretrgoma. Enkratna<br />
pridobitev dela ni zanimiva. Stalno in zvezno pridobivanje dela pa je mogoče le, če<br />
se delovna snov po končani ekspanziji vrne v začetno stanje. Da bi se to zgodilo, je<br />
treba vložiti nekaj dela in odvesti preostanek toplote v okolico. Cilj je, da je pri konstantnem<br />
dovodu toplote pridobljeno delo čim večje, vloženo pa čim manjše. Take<br />
procese imenujemo krožne procese, stroje, v katerih ti procesi potekajo, pa energetske<br />
stroje.<br />
Pri toplotnih krožnih procesih si preobrazbe v značilnih diagramih p – v ali T – s<br />
sledijo v smeri urinega kazalca, zato jih imenujemo desni krožni procesi. Toplota<br />
teče s telesa z višjo temperaturo preko delovne snovi na telo z nižjo temperaturo, pri
2.2 TERMODINAMIKA IN PRENOS TOPLOTE 23<br />
Slika 2.8: Shematski prikaz delovanja krožnega procesa<br />
tem pa se del toplote preobrazi v delo. Take krožne procese uporabljamo v toplotni<br />
tehniki za pridobivanje mehanskega dela oz. električne energije, slika 2.8A. Levi<br />
krožni procesi so nasprotni desnim, preobrazbe si sledijo v nasprotni smeri urinega<br />
kazalca. Vlagamo delo, toploto pa črpamo s telesa z nižjo temperaturo na telo z<br />
višjo. Take krožne procese uporabljamo v hladilni tehniki in v toplotnih črpalkah,<br />
slika 2.8B.<br />
Za proučevanje dejanskih krožnih procesov, ki jih srečamo v tehniški praksi, si pomagamo<br />
z ustreznimi teoretičnimi, primerjalnimi krožnimi procesi. Najpomembnejši<br />
taki procesi, ki so termodinamično gledano dobri in ki so se zaradi enostavnosti uveljavili<br />
v praksi, so zbrani v preglednicah 2.3 in 2.4. Pri tem so vzete za opis procesa<br />
vedno osnovne termodinamične preobrazbe, kot so izentropa, izoterma, izobara in<br />
izohora.<br />
Pri periodično delujočih volumenskih batnih strojih se posamezne preobrazbe v času<br />
ene periode izvajajo v enem samem stroju, npr. v valju motorja, zato je dogajanja<br />
najlažje prikazati v diagramu p – v, pri tem vse specifične energije nastopajo kot<br />
ploščine.<br />
Pri zvezno delujočih turbinskih strojih pa se posamezne preobrazbe izvajajo v več<br />
ločenih strojih in napravah, ki so med seboj povezane v postroj ali postrojenje. Ne<br />
glede na to dejstvo bo zaradi lažje primerjave privzeto, da se vse preobrazbe nanašajo<br />
na en sam stroj, dogajanja pa so prikazana v diagramu T – s, kjer spet vse specifične<br />
energije nastopajo kot površine, ali v Mollierovem diagramu h – s, kjer so specifične<br />
energije prikazane kot daljice.
24 2 TEORETIČNE OSNOVE<br />
Tabela 2.3: Termodinamične značilnosti najvažnejših teoretičnih krožnihprocesov<br />
Carnot Joule Stirling<br />
1-2 izoterma 3-4 1-2 izentropa 3-4 1-2 izoterma 3-4<br />
2-3 izentropa 4-1 2-3 izobara 4-1 2-3 izohora 4-1<br />
Q do Q do Q do<br />
Q 34 = m · R · T 3 · ln p 3<br />
p 4<br />
Q 23 = m · c p · (T 3 − T 2 ) Q 34 = m · R · T 3 · ln p 3<br />
p 4<br />
= m · R · T 3 · ln v 4<br />
v 3<br />
= m · (h 3 − h 2 ) = m · R · T 3 · ln v 1<br />
v 2<br />
Q od Q od Q od<br />
Q 12 = m · R · T 1 · ln p 2<br />
p 1<br />
Q 41 = m · c p · (T 4 − T 1 ) Q 12 = m · R · T 1 · ln p 2<br />
p 1<br />
= m · R · T 1 · ln v 1<br />
v 2<br />
= m · (h 4 − h 1 ) = m · R · T 1 · ln v 1<br />
v 2<br />
W = Q do − Q od W = Q do − Q od W = Q do − Q od<br />
W = m · R · (T 3 − T 1 ) · ln p 2<br />
p 1<br />
W = m · c p · (T 1 − T 2 + T 3 − T 4 ) W = m · R · (T 3 − T 1 ) · ln p 2<br />
p 1<br />
= m · R · (T 3 − T 1 ) · ln v 1<br />
v 2<br />
= m · (h 1 − h 2 + h 3 − h 4 ) = m · R · (T 3 − T 1 ) · ln v 1<br />
v 2<br />
η t = 1 − T 1<br />
T 3<br />
η t = 1 − T 1<br />
T 2<br />
= 1 − T 4<br />
T 3<br />
η t = 1 − T 1<br />
T 3
2.2 TERMODINAMIKA IN PRENOS TOPLOTE 25<br />
Tabela 2.4: Termodinamične značilnosti najvažnejših teoretičnih krožnihprocesov<br />
Otto Diesel Clausius-Rankine<br />
1-2 izentropa 3-4 1-2 izentropa 3-4 1-2 izentropa 3-4<br />
2-3 izohora 4-1 2-3 izobara, 4-1 izohora 2-3 izobara 4-1<br />
Q do Q do Q do<br />
Q 23 = m · c v · (T 3 − T 2 ) Q 23 = m · c p · (T 3 − T 2 ) Q 23 = m · (h 3 − h 2 )<br />
= m · (u 3 − u 2 ) = m · (h 3 − h 2 )<br />
Q od Q od Q od<br />
Q 41 = m · c v · (T 4 − T 1 ) Q 41 = m · c v · (T 4 − T 1 ) Q 41 = m · (h 4 − h 1 )<br />
= m · (u 4 − u 1 ) = m · (u 4 − u 1 )<br />
W = Q do − Q od W = Q do − Q od W = Q do − Q od<br />
W = m · c v · (T 1 − T 2 + T 3 − T 4 ) W = m ·<br />
= m ·<br />
R<br />
κ−1 · W = m · (h 1 − h 2 + h 3 − h 4 )<br />
R<br />
κ−1 · (T 1 − T 2 + T 3 − T 4 ) ·[κ · (T 3 − T 2 ) + (T 4 − T 1 )] ≈ m · (h 3 − h 4 )<br />
η t = 1 − T 4−T 1<br />
T 3 −T 2 η t = 1 − T 4−T 1<br />
κ·(T 3 −T 2 ) η t = 1 − h 4−h 1<br />
h 3 −h 2<br />
≈ 1 − h 4<br />
h 3
26 2 TEORETIČNE OSNOVE<br />
2.2.3 Delo, moč in izkoristek<br />
V tehniki je predvsem važno tehnično delo W t , to je delo, ki ga dobimo na gredi<br />
idealno delujočega ekspanzijskega stroja, pri tem se v delovni snovi zmanjša tlak od<br />
neke začetne do neke končne vrednosti.<br />
∫ p2<br />
W t = m · v · dp (2.2.7)<br />
p 1<br />
Podobno velja za kompresijski stroj. Tehnično delo W t je delo, ki ga je treba dovajati<br />
gredi idealno delujočega kompresijskega stroja, pri tem se v delovni snovi poviša tlak<br />
od neke začetne do neke končne vrednosti.<br />
To delo ekspanzije oziroma kompresije, če poteka proces v nasprotni smeri, je v<br />
diagramu diagramu p – V enako ploskvi p 1 -1-2-p 2 , v diagramu h – s pa daljici 1 -<br />
2, slika 2.9. Pri tem poteka teoretično ekspanzija ali kompresija delovne snovi po<br />
izentropi.<br />
Dejansko tehnično delo je pri ekspanzijskih strojih manjše, in pri kompresijskih strojih<br />
večje, saj stroji niso idealno delujoči, ampak se pri vseh strojih del energije zaradi<br />
več ali manj nepopolnih konstrukcij neželeno spremeni v toploto in odteka v okolico.<br />
Slika 2.9: Tehnični delo idealno delujočega ekspanzijskega oz. kompresijskega stroja
2.2 TERMODINAMIKA IN PRENOS TOPLOTE 27<br />
Delo enkratne ekspanzije oziroma kompresije W u je v diagramu p – V enako ploskvi<br />
pod krivuljo V 1 -1-2-V 2 in je povezano s spremembo prostornine:<br />
W u =<br />
∫ V2<br />
V 1<br />
p · dV (2.2.8)<br />
To ekspanzijsko ali kompresijsko delo je tisto, ki ga opravi ali dobi delovna snov, ki<br />
je v neki zaprti posodi pri eni sami preobrazbi (ekspanziji ali kompresiji).<br />
Poleg omenjenega tehničnega dela je pri toplotnih krožnih procesih važno delo krožnega<br />
procesa W ; to je pridobljeno tehnično delo, zmanjšano za vloženo tehnično<br />
delo, ki je potrebno, da se vrne delovna snov v začetno stanje. Pri Jouleovem krožnem<br />
procesu je npr. delo krožnega procesa enako delu plinske turbine W tE , zmanjšano za<br />
delo kompresorja W tK , ki stiska delovno snov na prvotni tlak:<br />
W = W tE − W tK (2.2.9)<br />
Moč P je delo W opravljeno v časovni enoti; velja:<br />
P = dW dt<br />
= Ẇ (2.2.10)<br />
Pri vsaki pretvorbi ene vrste energije v drugo imamo opravka z izgubami: del dovedene<br />
energije se ne spremeni v želeno obliko, npr. v toploto namesto v delo. Izgube<br />
imajo dva različna izvora in so posledica:<br />
• narave pretvorbe energij v krožnem procesu in<br />
• nepopolnih konstrukcij strojev in naprav.<br />
Merilo za učinkovitost preobrazbe energije je izkoristek. Ta je vedno definiran kot<br />
razmerje med izkoriščeno in vloženo energijo (delom, toploto) ali med izkoriščenim<br />
in vloženim energijskim tokom (močjo, toplotnim tokom). Pogosto so izkoristki definirani<br />
glede na dane merilne možnosti, in sicer tako, da nazorno pokažejo slabosti<br />
posameznih sklopov stroja ali naprave. V nadaljevanju je navedenih nekaj definicij<br />
za glavne vrste izkoristkov, ki so pomembni za energetske stroje. Dejansko delo<br />
ali dejanska moč, ki jo dobimo iz energetskega pogonskega stroja, je torej manjša;<br />
dejanska moč, ki jo moramo vložiti v energetski delovni stroj, pa večja.<br />
Pri toplotnih pogonskih strojih imamo opraviti tako s krožnim procesom delovne<br />
snovi kot tudi s strojem, zato obravnavamo izgube ločeno glede na izvor.
28 2 TEORETIČNE OSNOVE<br />
Termični izkoristek toplotnega (pogonskega) stroja<br />
Termični izkoristek se nanaša zgolj na nepovračljivosti pretvorb energije v krožnem<br />
procesu. Definiran je kot razmerje med delom ali močjo povračljivo delujočega<br />
krožnega procesa – torej idealnega stroja – in toploto ali toplotnim tokom, ki preide<br />
na delovno snov v krožnem procesu. Najboljši termični izkoristek ima povračljiv<br />
Carnotov krožni proces. Ta je tem večji, čim višja je temperatura dovoda toplote in<br />
čim nižja je temperatura odvoda toplote. Praktično Carnotovega krožnega procesa ni<br />
mogoče uresničiti, zato se mu skušamo pri dejanskih krožnih procesih čim bolj približati,<br />
slika 2.10. Pojem termičnega izkoristka ima smisel samo pri toplotnih krožnih<br />
procesih; upošteva dejstvo, da je toplota le delno pretvorljiva v delo.<br />
Q<br />
η t =<br />
˙ do − Q ˙ od<br />
= Ṗ<br />
(2.2.11)<br />
Q˙<br />
do Q do<br />
Indicirani in notranji izkoristek<br />
Pri volumenskih batnih strojih govorimo o indiciranem izkoristku, izračunamo ali<br />
določimo ga iz indikatorskega diagrama p - v. Pri turbinskih strojih govorimo o notranjem<br />
izkoristku, izračunamo ali določimo ga iz diagramov T - s ali h - s in je<br />
primerljiv z indiciranim. Oba, indicirani in notranji izkoristek se nanašata na nepopolne<br />
izvedbe konstrukcij strojev.<br />
Slika 2.10: Carnotov in carnotiziran krožni proces; A - izobaren dovod toplote,<br />
B - izentropna ekspanzija, C - izobaren odvod toplote, Č - izentropna kompresija
2.2 TERMODINAMIKA IN PRENOS TOPLOTE 29<br />
Indicirani izkoristek je pri pogonskih strojih razmerje med dejansko močjo na gredi<br />
stroja in močjo idealno delujočega stroja. Pove, koliko energijskega toka prehaja z<br />
delovne snovi na bat volumenskega stroja oz. koliko energijskega toka prehaja v<br />
gonilniku na lopatice turbinskega stroja. Nasprotno je pri delovnih strojih.<br />
Pogonski stroj:<br />
η i = P i<br />
P<br />
(2.2.12)<br />
Delovni stroj:<br />
η i = P P i<br />
(2.2.13)<br />
Zmnožek termičnega in notranjega izkoristka:<br />
η t · η i =<br />
P i<br />
˙ Q do<br />
(2.2.14)<br />
je razmerje med močjo, ki je na razpolago na gredi toplotnega pogonskega stroja, in<br />
toplotnim tokom, ki je bil doveden v krožni proces.<br />
Mehanski izkoristek<br />
Mehanski izkoristek je pri pogonskih strojih razmerje moči na gredi med strojem in<br />
generatorjem, pri delovnih strojih pa razmerje moči na gredi med elektromotorjem in<br />
strojem. Upošteva izgube zaradi trenja v ležajih, drsnikih in vodilih.<br />
Pogonski stroj:<br />
η m = P e<br />
P i<br />
(2.2.15)<br />
Delovni stroj:<br />
η m = P i<br />
P e<br />
(2.2.16)<br />
Dejanski (efektivni) izkoristek<br />
Dejanski izkoristek je zmnožek posameznih izkoristkov. Pri pogonskih strojih je to<br />
razmerje med dejansko potrebnim in teoretično izračunanim energijskim tokom. Pri<br />
delovnih strojih je nasprotno: dejanski izkoristek stroja je razmerje med teoretičnim<br />
in dejansko vloženim energijskim tokom.<br />
Pogonski stroj:<br />
η e =<br />
Ṗ · Pi<br />
Q do P · Pe = P e<br />
(2.2.17)<br />
P i Q˙<br />
do
30 2 TEORETIČNE OSNOVE<br />
Delovni stroj:<br />
η e = P P i<br />
· Pi<br />
P e<br />
= P P e<br />
(2.2.18)<br />
Termični izkoristek krožnega procesa je značilnost toplotnih pogonskih strojev, postrojev<br />
ali postrojenj, medtem ko pri drugih strojih, npr. aerohidravličnih, nima pravega<br />
smisla, saj pretvorba energije v njih ni vezana na krožni proces (toploto). Termični<br />
izkoristek η t , ki ima približne vrednosti od 0,5 do 0,6, najmočneje vpliva na<br />
dejanski izkoristek vsakega toplotnega pogonskega stroja, slika 2.11.<br />
Če imamo poleg naštetih izgub še druge, npr. izgube zaradi nepopolne toplotne izolacije<br />
η I , izgube zaradi sevanja η S , izgube zaradi zobniškega ali jermenskega prenosa<br />
η R itd., velja:<br />
η e = η t · η i · η m · η I · η S · η R · · · (2.2.19)<br />
Pogonski stroj:<br />
P > P i > P e (2.2.20)<br />
Delovni stroj:<br />
P e > P i > P (2.2.21)<br />
Slika 2.11: Shematski prikaz najvažnejših izgub pri pogonskih strojih
2.2 TERMODINAMIKA IN PRENOS TOPLOTE 31<br />
2.2.4 Tok tekočin skozi šobe<br />
Za pravilno delovanje energetskih strojev, posebno turbinskih, je pomembno razumevanje<br />
toka tekočine skozi kanale. Če v kanalu med dvema turbinskima lopaticama<br />
hitrost tekočine narašča, statični tlak pa pada, imenujemo tak kanal šoba (konfuzor);<br />
v nasprotnem primeru govorimo o difuzorju. Tok tekočine skozi šobo ima odločilno<br />
vlogo pri zasnovi vseh vrst turbinskih strojev (plinskih, parnih in vodnih turbin, turbokompresorjev<br />
itd.) ter različnih vrst raketnih pogonov. Hitrostne razmere na vstopu<br />
in na izstopu iz kanala je mogoče določiti z energijsko enačbo. Za izentropni tok<br />
tekočine (Q 12 = 0), brez opravljanja dela (W t12 = 0) in neupoštevanja potencialne<br />
energije (H 1 = H 2 ) lahko energijsko enačbo (2.1.1) poenostavimo:<br />
h 1 + v2 1<br />
2 = h 2 + v2 2<br />
2<br />
(2.2.22)<br />
Vsota tlačne, notranje in kinetične energije se pri plinskih tokovih pogosto označuje<br />
s totalno entalpijo, slika 2.12:<br />
h tot = h + v2<br />
2<br />
(2.2.23)<br />
Slika 2.12: Pretvorba energije pri izentropni ekspanziji
32 2 TEORETIČNE OSNOVE<br />
Če pospešujemo plin v šobi iz skoraj mirujočega stanja (v1 2 ≪ v2 2 ) na neko določeno<br />
izstopno hitrost v 2 , velja energijska enačba v poenostavljeni obliki (primer: turbinska<br />
stopnja parne turbine):<br />
√<br />
v 2 = 2 · (h 1 − h 2 ) (2.2.24)<br />
Za idealni plin se enačba (2.2.24) za izstopno hitrost poenostavi (primer: turbinska<br />
stopnja plinske turbine):<br />
√<br />
v 2 = 2 · c p · (T 1 − T 2 ) (2.2.25)<br />
Za kapljevine privzamemo, da so nestisljive (ϱ 1 = ϱ 2 = ϱ), notranja energija ostane<br />
praktično konstantna (U 2 = U 1 ); v tem primeru se energijska enačba nadalje poenostavi,<br />
slika 2.12 (primer: šoba Peltonove turbine):<br />
√<br />
v 2 = 2 · p1 − p 2<br />
(2.2.26)<br />
ϱ<br />
Če je tlačna energija pred šobo posledica hidrostatičnega tlaka, se hitrost v 2 lahko<br />
izrazi tudi z višinsko razliko, kot je to prikazano v razdelku 2.1.7.<br />
V nadaljevanju se bomo omejili na obravnavo izentropnega toka idealnega plina brez<br />
trenja skozi šobo. Enačbo (2.2.25) preoblikujemo z upoštevanjem:<br />
• izobarne specifične toplote:<br />
c p = R ·<br />
κ<br />
κ − 1<br />
• termične enačbe stanja idealnega plina:<br />
(2.2.27)<br />
T 1 = p 1<br />
ϱ 1 · R<br />
(2.2.28)<br />
• izentropne ekspanzije:<br />
( ) 1<br />
ϱ 1 p2 κ<br />
=<br />
ϱ 2 p 1<br />
( ) κ−1<br />
T 2 p2 κ<br />
=<br />
T 1 p 1<br />
(2.2.29)<br />
(2.2.30)<br />
• kontinuitetne enačbe:<br />
ṁ = ϱ 1 · v 1 · A 1 = ϱ 2 · v 2 · A 2 (2.2.31)
2.2 TERMODINAMIKA IN PRENOS TOPLOTE 33<br />
Pri pogoju, da je vstopna hitrost majhna (v1 2 ≪ v2 2 ), izrazimo hitrost in gostoto masnega<br />
toka na izstopu iz šobe:<br />
v 2 =<br />
=<br />
√<br />
2 · c p · T 1 ·<br />
√ 2 ·<br />
κ<br />
κ − 1 · p1<br />
ϱ 1<br />
·<br />
(<br />
1 − T )<br />
2<br />
=<br />
T 1<br />
[<br />
1 −<br />
( ) 1<br />
ṁ<br />
A = ϱ p2 κ<br />
2 · v 2 = ϱ 1 · ·<br />
p 1 = √ 2 · ϱ 1 · p 1 · √ κ<br />
κ − 1 ·<br />
( ) κ−1 ]<br />
p2 κ<br />
p 1<br />
[<br />
( )<br />
√ κ 2 ·<br />
κ − 1 · p1<br />
κ−1 ]<br />
p2 κ<br />
· 1 −<br />
=<br />
ϱ 1 p 1<br />
[ (p2 ) 2 ( ) κ+1 ]<br />
κ p2 κ<br />
− =<br />
p 1 p 1<br />
(2.2.32)<br />
= √ 2 · ϱ 1 · p 1 · ψ (2.2.33)<br />
Slika 2.13: Pretočna funkcija v odvisnosti od tlačnega razmerja
34 2 TEORETIČNE OSNOVE<br />
Konvergentna šoba<br />
V enačbi (2.2.33) vsebuje prvi koren veličine stanja pred šobo (indeks 1), drugi pa<br />
vrednosti, ki so odvisne samo od plina in od razmerja tlakov za in šobo in pred njo<br />
(p 2 /p 1 ). Vrednost drugega korena lahko izračunamo za poljubno mesto v šobi, če<br />
poznamo potek lokalnega tlaka v njej. Imenujemo ga pretočna funkcija:<br />
ψ(p/p 1 ) = √ κ<br />
[ ( ) 2 ( ) κ+1 ]<br />
p κ p κ<br />
− (2.2.34)<br />
κ − 1 p 1 p 1<br />
Za dani plin je v enačbi (2.2.34) eksponent izentrope znan, pretočna funkcija je torej<br />
odvisna le od razmerja lokalnega tlaka v šobi in tlaka pred njo (p/p 1 ). Pri konvergentni<br />
šobi je tlačno razmerje na vstopu (p 1 /p 1 ) = 1, nato zaradi pospeševanja toka plina<br />
pada. Funkcija ima dve ničli: na vstopu pri (p/p 1 ) = 1 in na izstopu (p/p 1 ) = 0,<br />
slika 2.13.<br />
Kontinuitetno enačbo za tok idealnega plina v šobi brez trenja lahko zapišemo glede<br />
na enačbo (2.2.33) v obliki:<br />
ṁ = A · √2<br />
· ϱ 1 · p 1 · ψ (2.2.35)<br />
Da je zadoščeno kontinuitetni enačbi, se mora z zmanjševanjem prereza povečevati<br />
pretočna funkcija, dokler ne doseže svoje največje vrednosti, pri kateri mora biti<br />
prerez najmanjši; to pa je lahko pri konvergentni šobi samo na izstopu. Največjo<br />
vrednost pretočne funkcije dobimo pri pogoju ∂ψ/∂(p 2 /p 1 ) = 0:<br />
1<br />
2 · κ ·<br />
p 1<br />
)kr<br />
2 · κ<br />
κ − 1 ·<br />
[ ( ) 2−κ<br />
2<br />
κ ·<br />
p2 κ<br />
− κ + 1 ( ) 1<br />
p2 κ<br />
·<br />
p 1 κ p 1<br />
]<br />
= 0 (2.2.36)<br />
Prvi del zmnožka ne more biti enak nič, pač pa je lahko nič razlika v oglatem oklepaju.<br />
Po preureditvi dobimo:<br />
( ( ) κ<br />
p2 2 κ−1<br />
=<br />
(2.2.37)<br />
κ + 1<br />
Tabela 2.5: Kritična tlačna razmerja in pripadajoče vrednosti pretočne funkcije za<br />
različne idealne pline<br />
Vrsta plina κ (p 2 /p 1 ) kr ψ max<br />
Enoatomni plini 1,667 0,487 0,514<br />
Dvoatomni plini 1,400 0,528 0,484<br />
Triatomni plini 1,300 0,546 0,472<br />
1,135 0,577 0,449
2.2 TERMODINAMIKA IN PRENOS TOPLOTE 35<br />
To je tlačno razmerje, pri katerem doseže pretočna funkcija ψ svojo največjo vrednost,<br />
imenujemo jo kritično ali Lavalovo razmerje. Največjo vrednost pretočne<br />
funkcije je mogoče izračunati tako, da vstavimo enačbo (2.2.37) v enačbo (2.2.32).<br />
Pri kritičnem tlačnem razmerju doseže tudi hitrost v 2 svojo največjo vrednost, imenujemo<br />
jo kritična ali Lavalova hitrost:<br />
√<br />
[ ( )]<br />
κ<br />
v 2 = v kr = 2 ·<br />
κ − 1 · p1 2<br />
1 −<br />
ϱ 1 κ + 1<br />
√ √<br />
√<br />
2 · κ<br />
=<br />
κ + 1 · p1 2 · κ<br />
=<br />
ϱ 1 κ + 1 · R · T 1 = 2 · κ − 1<br />
κ + 1 · c p · T 1 (2.2.38)<br />
Mogoče je dokazati, da je kritična ali Lavalova hitrost enaka zvočni hitrosti toka.<br />
Enačba (2.2.38) je v napisani obliki prikladna, saj je kritična hitrost določena z veličinami<br />
stanja plina pred šobo, ki so navadno znane, slika 2.13 in preglednica 2.5.<br />
Z upoštevanjem enačb (2.2.30) in (2.2.37) dobimo kritično temperaturo pri kritični<br />
hitrosti v kr :<br />
T kr = 2<br />
κ + 1 · T 1 (2.2.39)<br />
Če povežemo enačbi (2.2.38) in (2.2.39), dobimo kritično hitrost, izraženo z lokalnimi<br />
veličinami stanja:<br />
v kr = √ κ · R · T kr (2.2.40)<br />
S kakšnimi hitrostmi imamo opravka, nam pove Machovo število, ki je definirano kot<br />
razmerje med dejansko in kritično hitrostjo plina skozi šobo:<br />
Ma = v<br />
v kr<br />
(2.2.41)<br />
Tudi Machovo število je – enako, kot že preje omenjeno Reynoldsovo število – eden<br />
od važnejših brezdimenzijskih kriterijev. Ti bodo opisani v posebnem poglavju. Tlak<br />
v najmanjšem pretočnem prerezu šobe A min ne pade pod kritičnega, tudi tedaj ne,<br />
če imamo za šobo vakuum, izstopna hitrost pa ne more biti večja od zvočne hitrosti,<br />
Ma = 1.
36 2 TEORETIČNE OSNOVE<br />
Konvergentno-divergentna šoba<br />
Gostota masnega toka (ṁ/A) doseže na izstopu iz konvergetne šobe svojo največjo<br />
vrednost, enačba (2.2.35):<br />
(ṁ )<br />
= √ 2 · ϱ 1 · p 1 · ψ max (2.2.42)<br />
A<br />
max<br />
ali drugače napisano: pri kritičnem tlačnem razmerju je masni tok plina skozi najožji<br />
ali kritični prerez šobe A min :<br />
ṁ<br />
A min<br />
= √ 2 · ϱ 1 · p 1 · ψ max (2.2.43)<br />
Glede na enačbo (2.2.43) in lastnost pretočne funkcije ψ, ki doseže največjo vrednost<br />
pri kritični hitrosti, lahko dobimo večjo hitrost od kritične le, če se začne prerez šobe<br />
povečevati. Pri nadaljnji ekspanziji prevlada namreč zmanjševanje specifične gostote<br />
ϱ (večanje specifične prostornine) plina proti povečevanju hitrosti plina. Šobo, ki<br />
zadostuje temu pogoju, imenujemo konvergento-divergentna ali Lavalova šoba. Sestavljena<br />
je iz treh delov:<br />
• konfuzor; podkritično območje, kjer se prerez zmanjšuje, hitrost pa narašča;<br />
• grlo; najožje mesto, kjer imamo zvočno hitrost in<br />
• difuzor; nadkritično območje, kjer se prerez povečuje, hitrost pa narašča.<br />
Izkustvo je pokazalo, da kot v difuzorskem delu šobe ne sme biti večji kot 10–12 ◦ ,<br />
sicer se tok odlepi od sten dufuzorja.<br />
Krivulje Fanno<br />
Povezava kontinuitetne (2.1.9) in energijske enačbe (2.2.23) nam da enačbo:<br />
h tot = h + v2<br />
2 = h + 1 (ṁ ) 2<br />
2 · ϱ 2 ·<br />
(2.2.44)<br />
A<br />
Ob poznanju totalne entalpije lahko za poljubno gostoto ϱ z enačbo (2.2.44) izračunamo<br />
specifično entalpijo pri konstantni vrednosti (ṁ/A) in dobimo krivuljo Fanno.<br />
Drugače napisano, krivulja Fanno podaja vsa možna stanja plina (p, T , ϱ, h, v) pri<br />
toku skozi določen prerez s konstantno gostoto masnega toka (ṁ/A).<br />
Na sliki 2.14 so v diagramu h - s narisane krivulje Fanno za konvergentno-divergentno<br />
šobo za tri različne gostote masnega toka. Vsaka krivulja ima pri ds = 0<br />
navpično tangento. Glede na drugi stavek termodinamike velja za to točko:<br />
T · ds = dh − 1 · dp (2.2.45)<br />
ϱ
2.2 TERMODINAMIKA IN PRENOS TOPLOTE 37<br />
Slika 2.14: Krivulje Fanno za tri različne prereze šobe pri konstantni gostoti masnega<br />
toka<br />
Z odvajanjem enačbe (2.2.44) po gostoti ϱ dobimo:<br />
dh − 1 (ṁ ) 2<br />
ϱ 3 · dϱ = 0 (2.2.46)<br />
A<br />
S povezavo enačb (2.1.9), (2.2.45) in (2.2.46) dobimo:<br />
(ṁ<br />
A<br />
) 2<br />
= 1 ϱ 2 · dp<br />
dϱ<br />
(2.2.47)<br />
in nadalje:<br />
v kr =<br />
√<br />
dp<br />
dϱ<br />
(2.2.48)<br />
Mogoče je dokazati, da je enačba (2.2.48) identična z enačbo (2.2.40): v točki krivulje<br />
Fanno z navpično tangento imamo torej kritično ali zvočno hitrost.
38 2 TEORETIČNE OSNOVE<br />
Slika 2.15: Labirintno tesnjenje v diagramu h – s<br />
Zgled. Labirintno tesnjenje<br />
Labirintno tesnjenje omogoča tesnjenje rotirajoče gredi glede na mirujoče ohišje turbinskega<br />
stroja pri poljubno veliki tlačni razliki med notranjostjo stroja in okolico,<br />
npr.: pri kompresorjih ali pri plinskih in parnih turbinah.<br />
Dogajanje je mogoče spremljati v diagramu h – s. Skozi labirintne reže enakega<br />
prereza A teče delovna snov, pri tem ekspandira, zato se močno poveča njena hitrost<br />
in s tem kinetična energija. V prostoru za režo se kinetična energija delovne snovi<br />
zaradi vrtinčenja preobrazi v toploto in pri konstantnem tlaku dovaja delovni snovi.<br />
Postopek se ponovi pri vsakem naslednjem labirintu, pri čemer se tlak od labirinta do<br />
labirinta zmanjšuje. Proces je adiabaten, zato ostane totalna entalpija nespremenjena.<br />
Ekspanzija delovne snovi se konča na krivulji Fanno, kot to prikazuje slika 2.15. Pri<br />
danem začetnem tlaku p 1 , znanem prerezu rež A in številu labirintov je treba s poskušanjem<br />
več krivulj doseči končni tlak p 2 . Pri premajhnem številu labirintov lahko<br />
doseže delovna snov zvočno hitrost, kar ni priporočljivo. V nasprotju z navadnim<br />
tesnjenjem dosežemo pri labirintnih tesnilkah tesnjenje le, če teče skozi vse labirinte<br />
določena količina delovne snovi.
2.2 TERMODINAMIKA IN PRENOS TOPLOTE 39<br />
Tlačni skoki<br />
Konvergentno-divergentna šoba zahteva za vsako ekspanzijo posebno konstrukcijo.<br />
Za izbrani protitlak p O je potreben pravilno izračunan izstopni prerez šobe A 2 . In<br />
nasprotno: vsaki šobi je za izbrani vhodni tlak p 1 potreben točno izračunan izhodni<br />
tlak p 2 . Slika 2.16 prikazuje v brezdimenzijski obliki potek tlakov vzdolž Lavalove<br />
šobe za različne protitlake p 2 .<br />
A<br />
B<br />
C<br />
Č<br />
D<br />
E<br />
p 1 > p O > p A<br />
Pri prevelikem tlačnem razmerju (p 2 /p 1 ) je tok delovne snovi podoben toku<br />
v konvergentni šobi, kjer je tlak vedno manjši od kritičnega. V tem tlačnem<br />
območju tudi v grlu šobe ni dosežena kritična hitrost.<br />
p O = p B<br />
Pri zmanjševanju tlačnega razmerja (p 2 /p 1 ) pride v grlu šobe do kritičnega<br />
tlaka, tok doseže kritično hitrost, ki v difuzorju preide v podkritično.<br />
p B > p O > p C<br />
Pri nadaljnjem zmanjševanju tlačnega razmerja (p 2 /p 1 ) doseže tok v grlu šobe<br />
kritično hitrost, ki nato v difuzorju preide v nadkritično. Ker je protitlak p 2 še<br />
vedno premajhen, nastane v difuzorju nezvezni tlačni skok iz nadkritičnega v<br />
podkritično območje. Na izstopu iz šobe je hitrost podkritična.<br />
Nezvezni tlačni skok je pojav, pri katerem se hitrost delovne snovi v trenutku<br />
sprevrže iz nadkritične v podkritično. V diagramu h – s se tak tlačni skok izrazi<br />
kot tlačna sprememba iz nadtlačne v podtlačno krivuljo Fanno.<br />
p O = pČ<br />
Tlačno razmerje (p 2 /p 1 ) je toliko zmanjšano, da se tlačni skok premakne na<br />
izstop difuzorja. Na izstopu iz šobe je tako hitrost kritična.<br />
pČ > p O > p D<br />
Tlačno razmerje (p 2 /p 1 ) se zmanjšuje še naprej. V tem območju, ki je nadkritično,<br />
nastanejo t. i. poševni tlačni skoki. Poševni tlačni skok preide v<br />
nezveznega, kakor hitro je v območju p B > p O > pČ. Hitrost na izstopu iz<br />
šobe je kritična.<br />
p O = p E<br />
Tlačno razmerje (p 2 /p 1 ) je enako izračunanemu. Tok ekspandira točno po<br />
izračunani krivulji, zato ni kompresijskih skokov. Hitrost na izstopu iz šobe je<br />
nadkritična.
40 2 TEORETIČNE OSNOVE<br />
F<br />
p E > p O > p F<br />
Tlačno razmerje (p 2 /p 1 ) je manjše od izračunanega, dejanski protitlak je nižji<br />
od računskega (doslej je bil dejanski tlak, razen v primeru E, vedno višji od<br />
računskega). Tok dokončno ekspandira po izstopi iz šobe. Hitrost na izstopu<br />
iz šobe je nadkritična.<br />
Slika 2.16: Tlačne razmere vzdolž Lavalove šobe
2.2 TERMODINAMIKA IN PRENOS TOPLOTE 41<br />
2.2.5 Prenos toplote<br />
Prenos toplote je načelno mogoč na tri načine: s prevodom toplote (kondukcijo), s<br />
konvekcijo in s sevanjem.<br />
Prevod toplote<br />
Prenos toplote je snovnega (npr. molekulskega) izvora in temelji na izmenjavi impulza<br />
med sosednjimi delci snovi, ki nihajo okrog svojih stalnih ravnotežnih leg.<br />
Toplota se razširja po notranjosti telesa, s tem da se prenaša z molekul, ki imajo več<br />
energije (so toplejše) na molekule, ki imajo manj energije (so hladnejše). Pri tem je<br />
toplotni tok na enoto površine premo sorazmeren s temperaturno razliko in obratno<br />
sorazmeren z razdaljo med površinama telesa, Fourierjev zakon:<br />
˙Q<br />
A = λ · ∆T<br />
∆y<br />
(2.2.49)<br />
Enačba velja za trdna telesa in tekočine. Sorazmernostni faktor je toplotna prevodnost<br />
λ in je v določenem območju veljavnosti enačbe konstantna vrednost. V splošnem<br />
pa je toplotna prevodnost funkcija snovi in temperature.<br />
Podobna zakonitost velja tudi za strižno silo, deljeno s ploščino, torej za strižno napetost<br />
v toku, primerjaj razdelek 2.1.9, Newtonov zakon:<br />
F<br />
A = τ S = η · ∆v<br />
∆y<br />
(2.2.50)<br />
in za spremembo koncentracije snovi, deljeno s ploščino, npr. prenos trdnih delcev<br />
ṁ 1 v dimnih plinih, Fickov zakon:<br />
(ṁ )<br />
= D · ∆ϱ i<br />
(2.2.51)<br />
A ∆y<br />
S<br />
Ta analogija velja pri tekočinah samo v primerih, kadar ne prevladuje konvekcija zaradi<br />
vzgonskih in zunanjih sil, npr. za bližino stene (indeks S), kjer je tok laminaren.<br />
Za spremembo temperature snovi je bistvena prevodnost λ, za spremembo hitrosti<br />
viskoznost η, in za spremembo koncentracije snovi difuzija D.<br />
Konvekcija in konvektivni prestop toplote<br />
Medtem ko pri prevodu toplota ”pronica” skozi materijo, se pri konvekciji toplota<br />
prenaša po prostoru z materijo. Transport materije je v tem primeru makroskopski in<br />
zaradi tega mogoč samo pri tekočinah.
42 2 TEORETIČNE OSNOVE<br />
Glede na silo, ki povzroča gibanje tekočine, razlikujemo naravno in prisilno konvekcijo.<br />
O naravni konvekciji govorimo, če se vzpostavi tok pod vplivom vzgonskih sil<br />
zaradi razlik gostot. Vzrok za to je lahko razlika temperatur ali razlika koncentracij.<br />
Pri prisilni konvekciji povzroča tok zunanja sila, npr. sila, ki jo ustvari razlika tlakov.<br />
O konvektivnem prestopu toplote ali snovi govorimo, če gre za toplotni ali snovni<br />
tok med tekočino in steno ali če gre za toplotni ali snovni tok neposredno med dvema<br />
tekočinama. Pri konvektivnem prestopu se ogrejejo delci tekočine, ki pridejo v stik<br />
s toplejšo steno, potujejo zaradi vzgonskih ali zunanjih sil od stene proč, prenašajo<br />
toploto, ki so jo sprejeli, na hladnejšo tekočino in končno na hladnejšo steno.<br />
Konvektivni prestop toplote in snovi opisujeta enačbi:<br />
˙Q<br />
A = v · ϱ · c p · ∆T (2.2.52)<br />
ṁ i<br />
A = v · ∆ϱ i (2.2.53)<br />
ki pa praktično nista uporabni, saj se hitrost tekočine v krajevno močno spreminja in<br />
jo zato ni mogoče določiti.<br />
Čeprav se dogajanje v sredini toka tekočine razlikuje od dogajanja ob steni (ali na<br />
meji med dvema tekočinama), je prestop toplote (snovi) mogoče zadovoljivo zajeti z<br />
veličino α (β), ki upošteva dogajanje ob steni (na meji) in v sredini toka.<br />
˙Q<br />
A<br />
= α · ∆T (2.2.54)<br />
ṁ i<br />
A = β · ∆ϱ i (2.2.55)<br />
Sevanje<br />
Sevanje je bistveno drugačen način prenosa toplote: ne potrebuje nobenega materialnega<br />
prevodnika. Pri prenosu toplote s sevanjem je treba omeniti tri faze:<br />
• pretvorbo termične notranje energije sevajočega telesa v elektromagnetno valovanje<br />
z valovnimi dolžinami od 0,8 µm do 300 µm;<br />
• širjenje elektromagnetnega valovanja po prostoru;<br />
• absorpcijo elektromagnetnega valovanja obsevanega telesa, ponovna pretvorba<br />
v termično notranjo energijo.
2.2 TERMODINAMIKA IN PRENOS TOPLOTE 43<br />
Podobnega prenosa snovi ni.<br />
Vsako telo s temperaturo večjo od 0 K seva. Količina oddane (sprejete) toplote je odvisna<br />
od emisivnosti (absorptivnosti) površine telesa, površine same in temperature.<br />
Emisivnost telesa je enaka njegovi absorptivnosti. Telesu, ki pri določeni temperaturi<br />
na svoji površini emitira (absorbira) največ sevalne energije, pravimo črno telo.<br />
Gostota toplotnega toka je premo sorazmerna četrti potenci temperature, Stefanov<br />
zakon:<br />
˙Q<br />
A = ε · σ · T 4 (2.2.56)<br />
Sorazmernosti faktor je Stefanova konstanta σ = 5, 67 · 10 −8 W/(m 2· K 4 ) za sevanje<br />
črnega telesa; ε < 1 je emisijski koeficient, ki pomeni razmerje energije, ki jo oddaja<br />
površina nekega telesa, in energije, ki jo oddaja črno telo. Sevalna moč črnega telesa<br />
pri temperaturi okolice 300 K je potemtakem:<br />
˙Q<br />
A = 5, 67 · 10−8 · 300 4 = 459 W/m 2 (2.2.57)<br />
Pri prenosu toplote med dvema različno ogretima površinama je treba upoštevati ne<br />
samo sevanje tople površine, ampak tudi sevanje hladne.<br />
Prehod toplote<br />
Pod izrazom prehod toplote razumemo prenos toplote s tekočine z višjo temperaturo<br />
skozi steno, za katero pa ni nujno, da je prisotna, na tekočino z nižjo temperaturo,<br />
v splošnem torej zaporedoma: prestop toplote (α 1 ) - prevod toplote (λ) - prestop<br />
toplote (α 2 ), slika 2.17.<br />
˙Q 1 = α 1 · A 1 · (T 1 − T S1 ) (2.2.58)<br />
˙Q S = λ δ · A S · (T S1 − T S2 ) (2.2.59)<br />
˙Q 2 = α 2 · A 2 · (T S2 − T S2 ) (2.2.60)<br />
Vsa toplota ali toplotni tok prehaja s tekočine 1 skozi steno v tekočino 2, zato velja:<br />
˙Q = ˙Q 1 = ˙Q S = ˙Q 2 (2.2.61)<br />
Če je stena zakrivljena, je potrebno ploščino površine A S nadomestiti s srednjo ploščino<br />
površine po enačbi:<br />
A mS = A 1 − A 2<br />
ln A 1<br />
A 2<br />
(2.2.62)
44 2 TEORETIČNE OSNOVE<br />
Slika 2.17: Prehod toplote<br />
Za tanke, zakrivljene stene velja poenostavitev:<br />
A mS ≈ A 1 ≈ A S ≈ A 2 ≈ A (2.2.63)<br />
Če napisane enačbe preuredimo tako, da iz njih izpadejo vrednosti T S1 in T S2 , ki so<br />
le težko merljive, in predpostavimo tanko ravno steno, dobimo:<br />
( 1<br />
(T 1 − T S1 ) + (T S1 − T S2 ) + (T S2 − T 2 ) = + δ α 1 λ + 1 ) ˙Q<br />
· (2.2.64)<br />
α 2 A<br />
[<br />
]<br />
˙Q<br />
A = 1<br />
· (T 1 − T 2 ) (2.2.65)<br />
1/α 1 + δ/λ + 1/α 2<br />
Izraz v oglatem oklepaju se imenuje toplotna prehodnost k, pri čemer so vrednosti α<br />
in λ eksperimentalno določene in so za vsako tekočino drugačne; pogosto so napisane<br />
v brezdimenzijski obliki, primerjaj poglavje 2.3. Lastnost koeficienta toplotne<br />
prehodnosti je, da ne more biti večji od najmanjše toplotne prestopnosti min(α 1 , α 2 ).<br />
To pomeni, da je treba za izboljšanje toplotne prehodnosti k, npr. prenosnika toplote,<br />
izboljšati mesto, kjer je toplotna prestopnost α 1 najslabša. Če poznamo številčno<br />
vrednost za toplotno prehodnost k, lahko izračunamo toplotni tok po enačbi:<br />
˙Q = k · A · (T 1 − T 2 ) (2.2.66)<br />
To je ena od najpogosteje uporabljenih enačb pri prenosu toplote. Toplotni tok<br />
˙Q je torej premo sorazmeren s toplotno prehodnostjo k, s ploščino stene A med<br />
tekočinama in z razliko temperatur ∆T med toplo in hladno tekočino.
2.2 TERMODINAMIKA IN PRENOS TOPLOTE 45<br />
2.2.6 Goriva in zgorevanje<br />
Goriva<br />
Goriva so snovi, ki pri višjih temperaturah kemično reagirajo s kisikom (oksidirajo),<br />
pri tem spremenijo kemično sestavo in sočasno oddajo toploto. Vsebujejo sestavine<br />
(z malimi črkami so označeni masni deleži posameznih sestavin goriva):<br />
w c kg ogljika / kg goriva<br />
w h kg vodika / kg goriva<br />
kg žvepla / kg goriva<br />
w s<br />
w o<br />
kg kisika / kg goriva<br />
w v kg vode / kg goriva<br />
w p kg pepela / kg goriva<br />
w c + w h + w s + w o + w v + w p = 1 (2.2.67)<br />
Ogljik C, vodik H 2 in žveplo S so gorljive snovi, voda in pepel so nepotreben in<br />
odvečen balast. Slabše vrste goriv, npr. trboveljski rjavi premog in velenjski lignit,<br />
vsebujejo do 50 % negorljivih snovi in vlage.<br />
Pri stanju okolice je kemična sestava goriv trajno stabilna. Reakcija s kisikom poteka<br />
vedno pri bistveno višjih temperaturah, kot jo ima okolica, zato so goriva pomembni<br />
nosilci kemično vezane notranje energije. Nasprotno od drugih virov energije lahko<br />
kemično vezano energijo goriv pretvorimo v toploto tam, kjer jo potrebujemo, in<br />
takrat, kadar jo potrebujemo.<br />
Zgorevalna toplota in kurilnost<br />
Pomembna veličina za oceno goriva je zgorevalna toplota goriva H s . To je pri izobarnem<br />
zgorevanju vsa sproščena toplota, ki je enaka razliki entalpij udeleženih snovi<br />
pred zgorevanjem in po njem. Pri tem predpostavimo, da so v procesu zgorevanja<br />
vse udeležene snovi ohlajene na 0 ◦ C pri tlaku okolice 1,01325 bar (t. i. normalno<br />
stanje).<br />
Za praktično uporabo pa je najpomembnejša kurilnost goriva H i , preglednica 2.6.<br />
To je tisti del zgorevalne toplote, ki jo dobimo, če ne izkoristimo uparjalne toplote<br />
vodne pare v dimnih plinih, kar približno ustreza primeru, ko vse produkte zgorevanja<br />
ohladimo do temperature rosišča. Pri tem je izvor vode v dimnih plinih nepomemben;<br />
voda lahko nastane iz vodika v gorivu ali pa je prišla v proces zgorevanja kot vlaga v<br />
gorivu ali zraku. Za procese zgorevanja, kjer je v dimnih prisotna voda, velja:<br />
H s − H i = w V,D · r (2.2.68)<br />
pri čemer je w V,D količina vode, nastale pri zgorevanju 1 kg goriva in r kondenzacijska<br />
(uparjalna) toplota vode, r =2,5 MJ/kg pri 0 ◦ C.
46 2 TEORETIČNE OSNOVE<br />
Tabela 2.6: Kurilnosti nekaterih važnejših goriv<br />
Gorivo<br />
Kurilnost goriva<br />
H i / (MJ/kg)<br />
Antracit 31,82<br />
Koks 29,31<br />
Črni premog 27,21<br />
Rajvi premog 17,52<br />
Lignit 12,31<br />
Tekoča goriva 41,87<br />
Mazut 39,80<br />
Zemeljski plin 34,10<br />
Kapljeviti naftni plin 46,00<br />
Drugi plini ≈20,00<br />
V večini današnjih tehničnih naprav se izkorišča samo kurilnost goriva, kajti kondenzirana<br />
vodna para (H 2 O) reagira z žveplovim dvokisom (SO 2 ), ki je v manjših<br />
količinah pogosto v produktih zgorevanja (navadno jih imenujemo dimni plini), in<br />
tvori žveplasto (H 2 SO 3 ) in nato žvepleno kislino (H 2 SO 4 ). Kislini sta močno korozivni<br />
in povzročata škodo na konstrukcijskih materialih.<br />
Če poznamo kemično sestavo goriva, lahko z zgorevalnimi toplotami posameznih<br />
sestavin ocenimo kurilnost tega goriva. Za trdna in tekoča goriva velja empirična<br />
enačba:<br />
(<br />
H i = 33, 9 · w c + 121, 4 ·<br />
w h − w o<br />
8<br />
)<br />
+ 10, 5 · w s − 2, 5 · w v (2.2.69)<br />
Podobne enačbe je mogoče najti v strokovni literaturi tudi za plinasta goriva.<br />
Zgorevanje<br />
Zgorevanje je oksidacija goriva pri visokih temperaturah. Proces je eksotermen, pri<br />
njem se torej sprošča toplota. Pri zgorevanju sproščena kemična energija poveča<br />
notranjo energijo dimnih plinov in s tem njihovo temperaturo.<br />
Zgorevanje je popolno, če se vse molekule goriva spojijo s kisikom – oksidirajo;<br />
v dimnih plinih torej ni več sledi goriva ali produktov delne oksidacije, npr. CO.<br />
Najvišjo temperaturo dimnih plinov dobimo, če:
2.2 TERMODINAMIKA IN PRENOS TOPLOTE 47<br />
• med zgorevanjem toplote dimnih plinov ne odvajamo – imamo adiabatno zgorevanje;<br />
• je zgorevanje popolno;<br />
• je zgorevanje stehiometrično, brez presežka kisika.<br />
Dejanska temperatura dimnih plinov v kuriščih je občutno nižja zaradi spontanega<br />
odvoda toplote skozi stene, nepopolnega zgorevanja in presežka zraka. Sam proces<br />
zgorevanja je hiter in nepovračljiv, zato je – eksergijsko gledano – neugoden.<br />
Zgorevalne procese lahko opišemo s stehiometričnimi enačbami, najvažnejše so:<br />
C + O 2 = CO 2 + 406,1 MJ/kmol<br />
1,0 kg C + 2,7 kg O 2 = 3,7 kg CO 2 + 33,9 MJ/kg<br />
CO + 1 2 O 2 = CO 2 + 282,7 MJ/kmol<br />
1,0 kg CO + 0,56 kg O 2 = 1,56 kg CO 2 + 10,1 MJ/kg<br />
H 2 + 1 2 O 2 = H 2 O (para) + 241,9 MJ/kmol<br />
1,0 kg H 2 + 8,0 kg H 2 O = 9,0 kg H 2 O + 121,0 MJ/kg<br />
S + O 2 = SO 2 + 296,7 MJ/kmol<br />
1,0 kg S 2 + 1,0 kg O 2 = 2,0 kg SO 2 + 9,3 MJ/kg<br />
(2.2.70)<br />
(2.2.71)<br />
(2.2.72)<br />
(2.2.73)<br />
Pri tem so zaokrožene molske mase M glavnih udeležencev zgorevanja:<br />
M C = 12 kg/kmol<br />
M CO = 28 kg/kmol<br />
M H2 = 2 kg/kmol (2.2.74)<br />
M S<br />
= 32 kg/kmol<br />
M O2 = 32 kg/kmol<br />
Produkti zgorevanja so torej CO 2 , H 2 O in pogosto tudi SO 2 . Pri nepopolnem zgorevanju<br />
najdemo v zgorelih dimnih plinih še CO, nezgorele aromate in C (saje), pri<br />
prevelikem presežku zraka pa SO 3 in O 2 . Dušik N 2 iz zraka teče skozi proces zgorevanja<br />
praktično nespremenjen, pri visokih temperaturah se deloma veže s kisikom v<br />
dušikove okside NO, NO 2 , N 2 O 5 itd. (skupna oznaka NO x ). Pri temperaturah nad<br />
800 ◦ C se začne v večji meri pojavljati disociacija plinov, med seboj začneta reagirati<br />
tudi inertna plina dušik in kisik.
48 2 TEORETIČNE OSNOVE<br />
Tabela 2.7: Razmernik zraka<br />
Naprava<br />
Gorilniki plinskih turbin<br />
Motorji z notranjim zgorevanjem<br />
Vročevodni in parni kotli<br />
Razmernik zraka λ<br />
2−3<br />
1,0−1,1<br />
1,05−1,6<br />
V tehničnih napravah jemljemo kisik za zgorevanje iz zraka; na zemeljski površini<br />
sestavljata zrak prostorninska deleža kisika ϕ O2 ≈ 21 % in dušika ϕ N2 ≈ 79 %. Navedena<br />
prostorninska deleža veljata za suh zrak. Navadno imamo opravka z vlažnim<br />
zrakom, zato je potrebno pri natančnejših računih upoštevati tudi ustrezni delež vode<br />
v zraku. Najmanjša količina zraka v Zmin , ki je potrebna pri zgorevalnem procesu 1<br />
kg goriva, je potemtakem:<br />
v Zmin = v O min<br />
0, 21<br />
m 3 /kg (2.2.75)<br />
kjer je v Omin minimalna količina kisika za popolno zgorevanje izražena s prostornino<br />
pri temperaturi 0 ◦ C in tlaku 1,01325 bar (normalni kubični meter). Dejanska količina<br />
zraka, ki je dovedena v zgorevalni prostor, mora biti v praksi večja, da je zgorevanje<br />
zanesljivo popolno:<br />
λ =<br />
v Z<br />
v Zmin<br />
> 1 (2.2.76)<br />
Razmernik zraka λ je definiran kot količnik med dejansko in teoretično potrebno<br />
količino zraka. Odvisen je od namena zgorevanja, vrste goriva in konstrukcije zgorevalnega<br />
prostora, preglednica 2.7.<br />
Za pravilno zgorevanje je treba razmernik zraka stalno nadzorovati, to je še posebej<br />
pomembno pri vseh vrstah kotlov. Pri premajhnem razmerniku zraka pride do nepopolnega<br />
zgorevanja; v dimnih plinih se pojavijo CO, aromati in saje, v pepelu pa<br />
ostajajo ostanki nezgorelega goriva. Pri prevelikem razmerniku zraka se popolnost<br />
zgorevanja sicer izboljša, toda pri tem se mora segrevati tudi večja količina zraka,<br />
pretežno sestavljenega iz dušika N 2 , ki pri procesu zgorevanja ne sodeluje, njegova<br />
toplota pa neizkoriščena odteka v okolico. Posledica prevelikega razmernika zraka<br />
pa je tudi nižja temperatura zgorevanja, ki vpliva na slabši izkoristek zgorevanja.
2.2 TERMODINAMIKA IN PRENOS TOPLOTE 49<br />
Tabela 2.8: Največji možni prostorninski delež ogljikovega dvokisa CO 2 v dimnih<br />
plinih nekaterih goriv<br />
Gorivo ϕ CO2 max<br />
/ %<br />
Zemeljski plin 11,0−13<br />
Bencin ≈15,5<br />
Kurilno olje 15,5−17,5<br />
Rjavi premog 18,5−19,7<br />
Lignit<br />
18,5−20,5<br />
Les ≈20,5<br />
Ogljik (oglje) 21<br />
Razmernik zraka kontroliramo s kemično analizo plinov: merimo količino kisika v<br />
zgorelih dimnih plinih. Še lažje je določiti razmernik zraka iz izmerjene količine<br />
ogljikovega dvokisa CO 2 v suhih dimnih plinih, kjer velja za večino trdih in tekočih<br />
goriv:<br />
λ = ϕ CO 2 max<br />
ϕ CO2<br />
(2.2.77)<br />
ϕ CO2 max<br />
je največji možni prostorninski delež ϕ CO2 v suhih dimnih plinih, ki nastane<br />
pri popolnem zgorevanju brez presežka zraka in je značilna veličina za vsako gorivo,<br />
preglednica 2.8. Kolikor večji je razmernik zraka pri zgorevanju določenega goriva,<br />
toliko manjši je dejanski prostorninski delež ϕ CO2 v dimnih plinih.<br />
Dimni – izpušni plini<br />
Za produkte zgorevanja se pri kotlih uporablja izraz dimni plini, pri motorjih z notranjim<br />
zgorevanjem in gorilnikih plinskih turbin pa izraz izpušni plini. Kurilnost H i<br />
je odvisna od kemične sestave goriva. Na dlani je misel, da mora obstajati odvisnost<br />
med kurilnostjo na eni strani ter zgorevalnim zrakom in zgorelimi plini na drugi. Na<br />
osnovi meritev je nastalo več empiričnih enačb za trda, kapljevita in plinasta goriva,<br />
ki to pričakovanje potrjujejo. Preglednica 2.9 prikazuje take zveze, ki jih je praksa<br />
dobro potrdila.<br />
Nadalje se je pokazalo, da se v splošnem ni treba ozirati na kemično sestavo goriva<br />
in da je srednja specifična toplota dimnih plinov c pD / (kJ/(m 3· K)) približno enaka<br />
za vse dimne pline, ne glede na vrsto goriva. To velja le kot groba ocena. Ker se<br />
dimni plini vedejo kot idealni plini, mora biti zato za vsa goriva enaka tudi odvisnost
50 2 TEORETIČNE OSNOVE<br />
Tabela 2.9: Izkustveni podatki za minimalno količino zraka in minimalno količino<br />
dimnih plinov (λ = 1), izraženih v normalnih kubičnih metrih na kilogram goriva<br />
v odvisnosti od kurilnosti H i / (MJ/kg) trdega ali kapljevitega oz. MJ/m 3 plinastega<br />
goriva<br />
Vrsta goriva Kurilnost Minimalna količina Minimalna količina<br />
H i zraka dimnih plinov<br />
v Z min / (m 3 /kg)<br />
v D min / (m 3 /kg)<br />
oz. (m 3 /m 3 ) oz. (m 3 /m 3 )<br />
Trdo MJ/kg 0,241 · H i + 0,50 0,213 · H i + 1,65<br />
Kapljevito MJ/kg 0,203 · H i + 2,00 0,265 · H i<br />
Plinasto bogato MJ/m 3 0,260 · H i + 0,25 0,272 · H i + 0,25<br />
Plinasto revno MJ/m 3 0,209 · H i 0,173 · H i + 1,00<br />
specifične entalpije dimnih plinov h D / (kJ/m 3 ) od temperature T D pri pogoju, da<br />
imajo ti dimni plini enak presežek zraka:<br />
h D = H i<br />
v D<br />
(2.2.78)<br />
kjer je v D / (m 3 /kg) volumen dimnih plinov, izražen na kilogram mase goriva.<br />
Na sliki 2.18 je prikazana entalpija dimnih plinov v odvisnosti od njihove temperature.<br />
Vsa toplota, ki se sprosti pri zgorevanju 1 kg goriva s kurilnostjo H i , je vsebovana<br />
v dimnih plinih, pri tem je bila temperatura vseh udeležencev zgorevanja na<br />
začetku enaka temperaturi okolice T O . Velja:<br />
H i = m D · c pD · (T D − T O ) (2.2.79)<br />
kjer je m D / (kg/kg) je masa dimnih plinov, ki se sprosti pri zgorevanju 1 kg goriva,<br />
c pD / (kJ/(kg K)) srednja specifična toplota dimnih plinov in T D temperatura dimnih<br />
plinov.<br />
Pri gorilnikih plinskih turbin so temperaturne razlike manjše, zato je dopustno vstaviti<br />
v račune srednjo specifično toploto. Pri motorjih z notranjim zgorevanjem pa<br />
so temperaturne razlike in tlaki tako veliki, da je treba srednjo specifično toploto<br />
izračunati po odsekih. Dejanska količina zraka v dimnih plinih, izražena v normalnih<br />
kubičnih metrih, pa se računa po enačbi:<br />
v Z,D = v Z − v Z min<br />
v D<br />
(2.2.80)
2.2 TERMODINAMIKA IN PRENOS TOPLOTE 51<br />
Slika 2.18: h D -T diagram za dimne pline<br />
Vpliv dimnih plinov na okolico<br />
Pri vsakem zgorevanju nastanejo v dimnih plinih spojine, ki so za okolico škodljive,<br />
predvsem: ogljikov dvokis CO 2 in dušikovi oksidi NO x , pogosto tudi žveplov dvokis<br />
SO 2 . To velja za vse industrijske procese, za termoelektrarne, za motorje z notranjim<br />
zgorevanjem in tudi za vse zgorevalne procese v široki porabi.<br />
Ogljikov dvokis CO 2 je eden glavnih povzročiteljev tvorjenja učinkov tople grede<br />
v ozračju. Drugi plini, ki povzročajo enak učinek, so še metan, dušikovi oksidi in<br />
freoni. Ti plini, za katere se je udomačilo ime ”toplogredni plini”, tvorijo v zgornjih<br />
plasteh zemeljskega ozračja sloj, ki ima enak učinek kot steklo v rastlinjakih:<br />
absorbirajo dolgovalovno energijo sončnega sevanja, se pri tem segrejejo in del te<br />
akumulirane toplote vračajo na Zemljo. Ta toplota povzroča počasno povečevanje<br />
povprečne temperature Zemljine atmosfere.<br />
Dušikovi oksidi NO x , od katerih je najnevarnejši NO 2 , posredno povzročajo tvorjenje<br />
ozona O 3 in sodelujejo pri tvorjenju toplotne grede. Njihov škodljivi vpliv na okolico<br />
še ni popolnoma pojasnjen.<br />
Žveplov dvokis SO 2 , reagira z vlago v dimnih plinih in zraku ter tvori žvepleno kislino.<br />
Iz ozračja se s časom te kapljice izločijo in padajo na zemljo kot kisli dež. Ta<br />
uničuje predvsem iglaste gozdove, škodi tudi drugemu rastlinju in povečuje kislost<br />
zemeljske površine.
52 2 TEORETIČNE OSNOVE<br />
2.3 Podobnost in dimenzijska analiza<br />
2.3.1 Kriteriji podobnosti<br />
Pri veliki večini energetskih strojev in naprav imamo opravka s turbulentnimi tokovi,<br />
ki analitično niso rešljivi. Na drugi strani pa je bistven čim bolj natančen preračun<br />
novega stroja ali naprave, kajti čim večji je stroj, ki ne deluje pravilno, tem dražje<br />
in časovno daljše so kasnejše spremembe in popravki. Te težave premosti v mnogih<br />
primerih izdelava primernega, praviloma pomanjšanega modela, ki pa se mora<br />
skladati z dejansko izvedbo v vseh bistvenih kriterijih podobnosti. Ti so izraženi kot<br />
brezdimenzijska števila, ki jih v splošnem lahko določimo na dva načina:<br />
• z zapisom enačb v brezdimenzijski obliki, npr. Navier-Stokesovih enačb, in<br />
• z dimenzijsko analizo.<br />
Dva fizikalna pojava sta si podobna, če je vrednost njunih karakterističnih brezdimenzijskih<br />
števil enaka, s tem je zadoščeno kriterijem podobnosti. Praktična uporabnost<br />
metode je v tem, da lahko rezultate laboratorijskih meritev na modelu prenesemo na<br />
realne naprave brez zahtevnega eksperimentalnega preverjanja. V nadaljevanju so<br />
obravnavani samo nekateri najpomembnejši kriteriji podobnosti s področja prenosa<br />
impulza, toplote in snovi.<br />
Model (indeks M) in izvedba (indeks I) si morata biti:<br />
• geometrijsko podobna (podobnost dolžinskih dimenzij)<br />
• kinematično podobna (podobnost vektorjev hitrosti in pospeškov)<br />
• dinamično podobna (podobnost vektorjev sil)<br />
• termično podobna in<br />
• snovno podobna<br />
V praksi se izkaže, da ni mogoče zadostiti vsem kriterijem podobnosti med modelom<br />
in izvedbo, zato se navadno zadovoljimo z ujemanjem tistih kriterijev podobnosti, ki<br />
imajo največji vpliv na opazovani pojav.<br />
Geometrijska podobnost<br />
Geometrijska podobnost med modelom in izvedbo je osnovni kriterij pri uporabi teorije<br />
podobnosti. Za njeno izpolnitev morajo biti karakteristične dimenzije modela in<br />
izvedbe v določenem medsebojnem razmerju:
2.3 PODOBNOST IN DIMENZIJSKA ANALIZA 53<br />
L M<br />
L I<br />
= konst.<br />
2<br />
L M<br />
2<br />
L I<br />
= konst. (2.3.1)<br />
3<br />
L M<br />
3<br />
L I<br />
= konst.<br />
Primer za geometrijsko podobnost je brezdimenzijsko število π, razmerje med obsegom<br />
kroga O in njegovim premerom d:<br />
π = O d<br />
(2.3.2)<br />
Kinematična podobnost<br />
Kinematična podobnost zahteva, da so si vektorji hitrosti in pospeškov modela in<br />
izvedbe med seboj proporcionalni. S kriteriji kinematične podobnosti se navadno ni<br />
treba posebej ukvarjati, saj jim je avtomatično zadoščeno z izpolnjevanjem dinamične<br />
podobnosti.<br />
Dinamična podobnost<br />
Za dosego podobnosti med modelom in izvedbo mora biti zadoščeno kriterijem dinamične<br />
podobnosti, kjer gre za razmerje sil, ki delujejo na tekočino pri modelu in<br />
Slika 2.19: Dinamična podobnost tokov na modelu in izvedbi
54 2 TEORETIČNE OSNOVE<br />
izvedbi, sočasno pa morata biti izpolnjeni tudi geometrijska ter kinematična podobnost,<br />
kar pomeni, da morajo za oba sistema veljati enake enačbe.<br />
Na masni delec tekočine delujejo v splošnem masne, tlačne in viskozne sile, ki so v<br />
ravnotežju z vztrajnostno silo. Model in izvedba sta si dinamično podobna, če so si<br />
posamezne sile v poljubni točki opazovanega pojava v sorazmerju, slika 2.19.<br />
Z dimenzijsko analizo lahko zapišemo za posamezne sile:<br />
vztrajnostna sila<br />
viskozna sila<br />
F v = m · a = ϱ · V · ∆v<br />
∆t<br />
∼ ϱ · L 3 ·<br />
v<br />
L/v = ϱ · v2 · L 2 (2.3.3)<br />
F η = η · A · ∆vx<br />
∆y<br />
∼ η · v · L (2.3.4)<br />
masna sila teže<br />
F g = m · g ∼ g · ϱ · L 3 (2.3.5)<br />
masna sila vzgona<br />
tlačna sila<br />
F β = β · g · ϱ · V · ∆T ∼ g · ϱ · β · L 3 · ∆T (2.3.6)<br />
F p = ∆p · A ∼ ∆p · L 2 (2.3.7)<br />
tlačna sila z upoštevanjem stisljivosti tekočin v povezavi z enačbo (2.1.2)<br />
F ξ = ∆p · A = ∆ϱ<br />
ϱ · E · A ∼ E · L2 (2.3.8)<br />
Kriterije dinamične podobnosti zapišemo kot razmerja med posameznimi silami.<br />
Razmerje med vztrajnostno in viskozno silo:<br />
F v<br />
= ϱ · v2 · L 2<br />
F η η · v · L<br />
= ϱ · v · L<br />
η<br />
= v · L<br />
ν<br />
= Re (2.3.9)<br />
kjer je kriterij podobnosti Re Reynoldsovo število. Za podobnost tokov dveh različnih<br />
tekočin, kjer imata prevladujoč vpliv vztrajnostna in viskozna sila, je dinamična podobnost<br />
izpolnjena z enakostjo Reynoldsovih števil:<br />
[ ] [ ]<br />
Fv Fv<br />
= ⇒ Re M = Re I (2.3.10)<br />
F η F η<br />
M<br />
I
2.3 PODOBNOST IN DIMENZIJSKA ANALIZA 55<br />
Razmerje med vztrajnostno in masno silo teže je:<br />
F v<br />
= ϱ · v2 · L 2<br />
F g g · ϱ · L 3 = v2<br />
g · L = F r (2.3.11)<br />
kjer je kriterij podobnosti F r Froudovo število. Za podobnost tokov dveh različnih<br />
tekočin, kjer imata prevladujoč vpliv vztrajnostna sila in masna sila teže, je dinamična<br />
podobnost izpolnjena z enakostjo Froudovih števil:<br />
[ ] [ ]<br />
Fv Fv<br />
= ⇒ F r M = F r I (2.3.12)<br />
F g F g<br />
M<br />
I<br />
Razmerje med masno silo vzgona in viskozno silo je:<br />
F β<br />
= g · ϱ · β · L3 · ∆T<br />
F η η · v · L<br />
= g · β · L3 · ∆T ν<br />
ν 2 ·<br />
v · L = Gr<br />
Re<br />
(2.3.13)<br />
kjer je kriterij podobnosti kvocient med Gr Grashofovim številom in Re Reynoldsovim<br />
številom. Za podobnost tokov dveh različnih tekočin, kjer imata prevladujoč<br />
vpliv masna sila vzgona in viskozna sila, je dinamična podobnost izpolnjena z enakostjo<br />
Grashofovih števil:<br />
[<br />
g · β · L3 · ∆T<br />
ν 2 ]<br />
M<br />
=<br />
[ ]<br />
g · β · L3 · ∆T<br />
ν 2<br />
Razmerje med tlačno in vztrajnostno je:<br />
F p<br />
F v<br />
=<br />
∆p · L2<br />
ϱ · v 2 · L 2 =<br />
I<br />
⇒ Gr M = Gr I (2.3.14)<br />
∆p = Eu (2.3.15)<br />
ϱ · v2 kjer je kriterij podobnosti Eu Eulerjevo število (v ameriški literaturi pogosto navedeno<br />
kot ”pressure coefficient”). Za podobnost tokov dveh različnih tekočin, kjer<br />
imata prevladujoč vpliv masna sila teže in vztrajnostna sila, je dinamična podobnost<br />
izpolnjena z enakostjo Eulerjevih števil:<br />
[ Fp<br />
F v<br />
]M<br />
=<br />
[ Fp<br />
F v<br />
]I<br />
⇒ Eu M = Eu I (2.3.16)<br />
Razmerje med vztrajnostno in tlačno silo je pri stisljivih tekočinah (v povezavi z<br />
enačbama 2.1.2 in 2.2.47 za kritično hitrost):<br />
F v<br />
= ϱ · v2 · L 2<br />
F ξ E · L 2 = v2<br />
E/ϱ = v2<br />
vkr<br />
2<br />
= Ma 2 (2.3.17)
56 2 TEORETIČNE OSNOVE<br />
kjer je kriterij podobnosti Ma Machovo število. Za podobnost tokov dveh različnih<br />
tekočin, kjer imata prevladujoč vpliv vztrajnostna sila in sila zaradi kompresije ali<br />
ekspanzije tekočine, je dinamična podobnost izpolnjena z enakostjo Machovih števil:<br />
[ ] [ ]<br />
Fv Fv<br />
= ⇒ Ma M = Ma I (2.3.18)<br />
F ξ F ξ<br />
M<br />
I<br />
Pri vrednostih Ma < 0, 3 lahko vplive stisljivosti tekočine na tokovno polje zanemarimo.<br />
Medsebojno povezavo najvažnejših brezdimenzijskih števil, ki pridejo v poštev pri<br />
prenosu impulza v tekočinah, prikazuje slika 2.20.<br />
Navadno lahko zadostimo samo enemu kriteriju podobnosti med modelom in izvedbo.<br />
Kadar je pomembno poznati več vplivov, ki so značilni za tok tekočine, je<br />
treba meritve ponavljati, tako da je vsakokrat zadoščeno tistemu kriteriju, ki bistveno<br />
določa tok tekočine. Take meritve so navadno obsežne in drage.<br />
Slika 2.20: Brezdimenzijska števila pri stacionarnem toku tekočine
2.3 PODOBNOST IN DIMENZIJSKA ANALIZA 57<br />
Zgled. Ventil za vodik<br />
Ugotoviti je treba, ali ustreza ventil, ki je vgrajen v cevovodu za zrak, tudi toku<br />
vodika. Tehnični podatki:<br />
zrak: T = 20 ◦ C vodik: T = 40 ◦ C<br />
p = 1 bar<br />
p = 8 bar<br />
ν = 1,5 · 10 5 m 2 /s ν = 1,45 · 10 5 m 2 /s<br />
v = 10 m/s<br />
v = 10 m/s<br />
Na tok skozi ventil odločilno vplivata vztrajnostna F v in viskozna sila F η , medtem<br />
ko smemo vpliv masne sile teže F g in silo zaradi razlike tlakov v cevovodu F p zanemariti.<br />
Oba tokova skozi ventil sta si podobna, če je izpolnjen kriterij Re Z = Re H .<br />
[ ] [ ]<br />
v · d v · d<br />
=<br />
ν ν<br />
Z<br />
H<br />
Od tod je mogoče takoj izračunati hitrost vodika v cevovodu: v H = 9,7 m/s. Ventil<br />
ustreza spremenjenim razmeram, saj je hitrost vodika v cevovodu v normalnih mejah.<br />
Termična podobnost<br />
Termična podobnost je težje uporabljiva, saj poleg kriterijev dinamične podobnosti<br />
zahteva še izpolnjevanje dodatnih kriterijev, ki so značilni za prenos toplote. Poleg<br />
veličin, ki nastopajo pri dinamični podobnosti, imamo pri termični podobnosti tokov<br />
še dve novi veličini: toplotni tok ˙Q in temperaturno razliko ∆T . Zveza med toplotnim<br />
tokom in drugimi veličinami, ki so bistvene za prenos toplote, je določena s<br />
sorazmernostnimi faktorji, ki pomenijo določeno snovno lastnost, npr. toplotna prevodnost<br />
λ, toplotna prestopnost α, specifična toplota pri konstantnem tlaku c p itd.<br />
Podobno kot pri dinamični podobnosti uporabimo dimenzijsko analizo za zapis toplotnega<br />
toka tudi pri termični podobnosti:<br />
prevod (kondukcija) toplote<br />
konvekcija<br />
˙Q λ = λ · A · ∆T<br />
∆y<br />
∼ λ · L · ∆T (2.3.19)<br />
˙Q K = ṁ · c p · ∆T = ϱ · v · A · c p · ∆T ∼ ϱ · c p · v · L 2 · ∆T (2.3.20)<br />
konvektivni prestop toplote<br />
sevanje<br />
˙Q α = α · A · ∆T ∼ α · L 2 · ∆T (2.3.21)<br />
˙Q σ = ε · σ · A · ∆T 4 ∼ ε · σ · L 2 · ∆T 4 (2.3.22)
58 2 TEORETIČNE OSNOVE<br />
Kriterije termične podobnosti lahko zapišemo kot razmerja med posameznimi toplotnimi<br />
tokovi, ki morajo biti za primer dinamične podobnosti tokov enaka za model in<br />
izvedbo.<br />
Razmerje med konvektivnim prestopom in prevodom toplote je:<br />
˙Q α<br />
˙Q λ<br />
= α · L2 · ∆T<br />
λ · L · ∆T<br />
= α · L<br />
λ<br />
= Nu (2.3.23)<br />
kjer je kriterij podobnosti Nu Nußeltovo število. Pri obravnavi prenosa toplote z ene<br />
snovi na drugo, kjer imata prevladujoč vpliv konvektivni prestop in prevod toplote,<br />
je termična podobnost izpolnjena z enakostjo Nußeltovih števil:<br />
[ ] ˙Q α<br />
˙Q λ<br />
M<br />
=<br />
[ ] ˙Q α<br />
˙Q λ<br />
I<br />
⇒ Nu M = Nu I (2.3.24)<br />
Razmerje med konvekcijo in kondukcijo (prevodom toplote) je:<br />
˙Q K<br />
˙Q λ<br />
= ϱ · c p · v · L 2 · ∆T<br />
λ · L · ∆T<br />
= ϱ · c p · v · L<br />
λ<br />
= P e (2.3.25)<br />
kjer je kriterij podobnosti P e Pecletovo število. Podobno kot je Reynoldsovo število<br />
razmerje med turbulentim in laminarnim prenosom impulza v tekočini, je Pecletovo<br />
število razmerje med konvekcijo in kondukcijo. Pri prenosu toplote, kjer imata prevladujoč<br />
vpliv na porazdelitev temperatur konvekcija in kondukcija, je termična podobnost<br />
izpolnjena z enakostjo Pecletovih števil:<br />
[ ] ˙Q K<br />
˙Q λ<br />
M<br />
=<br />
[ ] ˙Q K<br />
˙Q λ<br />
I<br />
⇒ P e M = P e I (2.3.26)<br />
Pogosto se namesto Pecletovega števila uporablja Prandtlovo število P r, ki je kvocient<br />
Pecletovega in Reynoldsovega števila:<br />
P r = ˙Q K<br />
· Fη = η · c p<br />
˙Q λ F v λ<br />
= P e<br />
Re<br />
(2.3.27)<br />
Iz zgornjega izraza je razvidno, da je Prandtlovo število funkcija snovnih lastnosti<br />
tekočine P r = P r(p, T ). Vrednosti za snovne lastnosti so navadno navedene v termodinamičnih<br />
tabelah ali diagramih. Prandtlovo število je pomemben pokazatelj pri<br />
opisu prenosa toplote, ker vsebuje informacijo o razmerju med debelino hidravlične<br />
in termične mejne plasti. Tako pomeni neposredno povezavo med hitrostnim in temperaturnim<br />
poljem toka tekočine.
2.3 PODOBNOST IN DIMENZIJSKA ANALIZA 59<br />
Razmerje med konvektivnim prestopom toplote in konvekcijo je:<br />
˙Q α<br />
˙Q K<br />
=<br />
α · L 2 · ∆T<br />
ϱ · v · c p · L 2 · ∆T =<br />
α<br />
ϱ · v · c p<br />
= St (2.3.28)<br />
kjer je podobnostni kriterij St Stantonovo število. To število je mogoče izraziti tudi s<br />
kombinacijo že znanih brezdimenzijskih števil:<br />
St = Nu<br />
P e =<br />
Nu<br />
Re · P r<br />
Razmerje med sevanjem in prevodom toplote je:<br />
˙Q σ<br />
˙Q λ<br />
= ε · σ · L2 · ∆T 4<br />
λ · L · ∆T<br />
= ε · σ · L · ∆T 3<br />
λ<br />
(2.3.29)<br />
= Sf (2.3.30)<br />
kjer je kriterij podobnosti Sf Stefanovo število. Pri prenosu toplote, kjer imata prevladujoč<br />
vpliv na porazdelitev temperatur sevanje in prevodnost, je termična podobnost<br />
izpolnjena z enakostjo Stefanovih števil:<br />
[ ] ˙Q σ<br />
˙Q λ<br />
M<br />
=<br />
[ ] ˙Q σ<br />
˙Q λ<br />
I<br />
⇒ Sf M = Sf I (2.3.31)<br />
Medsebojno povezavo najvažnejših brezdimenzijskih števil, ki pridejo v poštev pri<br />
prenosu toplote, prikazuje slika 2.21.<br />
Slika 2.21: Brezdimenzijska števila pri stacionarnem prenosu toplote
60 2 TEORETIČNE OSNOVE<br />
Slika 2.22: Podobnost modela in izvedbe pri prisilni konvekciji<br />
Na sliki 2.22 sta narisani dve tokovni polji in v teh poljih dva valja. Z vidika prisilne<br />
konvekcije sta si toka podobna, če so izpolnjeni naslednji pogoji:<br />
• geometrijska podobnost: velikost valjev in velikost območja opazovanja sta si<br />
proporcionalni;<br />
• kinematična podobnost: porazdelitvi hitrosti na mejah območja opazovanja sta<br />
si podobni po velikosti in smeri;<br />
• dinamična podobnost: enakost Reynoldsovih števil v območju opazovanja,<br />
Re M = Re I ;<br />
• termična podobnost: enakost Prandtlovih in Nußeltovih števil v območju opazovanja,<br />
P r M = P r I in Nu M = Nu I . Enakost Prandtlovih števil kaže na<br />
proporcionalnost temperaturnega polja, enakost Nußeltovih pa na proporcionalnost<br />
prestopa toplote.<br />
Iz omenjenega primera za prisilno konvekcijo je razvidno, da se kriterij termične<br />
podobnosti lahko izrazi s funkcijsko povezavo Nu = Nu(Re, P r). Pri naravni konvekciji<br />
na tok tekočine bistveno vpliva masna sila vzgona, zato namesto Reynoldsovega<br />
števila v kriteriju termične podobnosti nastopa Grashofovo število: Nu =<br />
Nu(Gr, P r).<br />
Eksperimentalne raziskave so pokazale, da je mogoče empirične enačbe za števila<br />
Nu prikazati v potenčni obliki:<br />
Nu = K · Gr m · P r n naravna konvekcija (2.3.32)<br />
Nu = K · Re m · P r n prisilna konvekcija (2.3.33)
2.3 PODOBNOST IN DIMENZIJSKA ANALIZA 61<br />
Konstanta K in eksponenta m in n so določeni z meritvami za vsako izvedbo prenosnika<br />
toplote posebej. Če na prenos toplote vplivajo še druge okoliščine toka, npr.<br />
geometrija telesa, stisljivost tekočine itd., je treba enačbe ustrezno razširiti in z meritvami<br />
določiti nove konstante. Empirične enačbe so zbrane v tehniških priročnikih in<br />
drugi literaturi, ki obravnava prenos toplote. Pri enačbah je navedeno tudi območje<br />
veljavnosti takih enačb. Primer praktične enačbe za prestop toplote na steno okrogle<br />
cevi za stacionarni turbuletni tok:<br />
Nu = α · d<br />
λ<br />
= 0,0235 · (Re 0,8 − 230) · (1,8 · P r 0,3 − 0,8) ·<br />
·<br />
[ ( ] ( ) d 3 P r 0,14<br />
1 + ·<br />
L)2<br />
P r S<br />
(2.3.34)<br />
Enačba velja v območju Re > 2300 in 0,7 < P r < 1000. Člen v oglatem oklepaju<br />
upošteva vpliv dolžine cevi L: na začetku cevi imamo namreč nerazvito mejno<br />
plast, prestop toplote je zato boljši. Mejna plast postaja z naraščajočo dolžino cevi<br />
debelejša, prestop toplote se slabša. Zadnji člen v enačbi upošteva vpliv spremembe<br />
snovnih lastnosti zaradi razlike temperatur tekočine na sredini cevi in ob steni. Omenjeni<br />
vpliv je majhen in je ga treba upoštevati le, kadar so te temperaturne razlike<br />
velike.<br />
Zgled. Prestop toplote iz tekočine na steno cevi<br />
Določiti je treba toplotno prestopnost α za kotlovsko cev, v kateri teče vroča voda, za<br />
naslednje tehnične podatke:<br />
T = 150 ◦ C<br />
p = 150 bar<br />
d = 0,038 m<br />
L = 20 m<br />
v = 3,0 m/s<br />
Iz tabel za lastnosti vode in vodne pare povzamemo:<br />
ϱ = 925,1 kg/m 3<br />
c p = 4,263 kW/(kg K)<br />
λ = 691,8·10 6 kW/(m K)<br />
η = 186,1·10 6 kg/(m s)
62 2 TEORETIČNE OSNOVE<br />
Iz navedenih podatkov izračunamo manjkajoče vrednosti:<br />
P r = η · c p<br />
λ = 186,1 · 10−6 · 4,263<br />
691,8 · 10 −6 = 1,147<br />
( P r<br />
P r S<br />
) 0,14<br />
≈ 1<br />
Re = ϱ · v · d<br />
η<br />
=<br />
Enačba (2.3.34) za prestop toplote:<br />
925,1 · 3,0 · 0,038<br />
186,1 · 10 −6 = 566 ′ 692<br />
Nu = 0,0235 · (566 ′ 692 0,8 − 230) · (1,8 · 1,147 0,3 − 0,8) ·<br />
[ ( ) 2<br />
]<br />
0,038 3<br />
· 1 +<br />
= 1022<br />
20<br />
Povprečna toplotna prestopnost z vode na steno cevi je:<br />
α = λ d<br />
· Nu =<br />
691,8 · 10−6<br />
0,038<br />
· 1022 = 18,6 kW/(m 2 · K)<br />
Snovna podobnost<br />
Pogosto imamo opraviti z dvo- ali večfaznimi sistemi, torej sistemi, v katerih nastopajo<br />
snovi (tekočine) različnih gostot ali agregatnih stanj. Primer take naprave je<br />
hladilni stolp. Na meji med tekočinama – med zrakom in hladilno vodo – poteka poleg<br />
intenzivnega prenosa toplote tudi izmenjava snovi v obeh smereh: voda prehaja<br />
na zrak in v manjši meri zrak v hladilno vodo. Pri snovni podobnosti imamo dve<br />
novi veličini: masni tok i-te komponente ṁ i , ki prestopa iz ene tekočine v drugo,<br />
in razlika gostot i-te komponente ∆ϱ i (razlika gostot je večkrat definirana tudi kot<br />
razlika parcialnih tlakov ali kot razlika koncentracije i-te komponente v eni in drugi<br />
tekočini). Mehanizmi prenosa snovi so popolnoma primerljivi z mehanizmi prenosa<br />
toplote. Tudi pri prenosu snovi je zveza med masnim tokom i-te komponente in<br />
drugimi veličinami določena s sorazmernostnimi faktorji, ki predstavljajo določeno<br />
snovno lastnost, npr. difuzija D, snovna prestopnost β.<br />
Snovno podobnost je mogoče izraziti podobno kot termično podobnost:<br />
difuzija<br />
ṁ i,D = D i · A · ∆ϱ i<br />
∆y<br />
∼ ϱ · D · L (2.3.35)
2.3 PODOBNOST IN DIMENZIJSKA ANALIZA 63<br />
Difuzijski prenos snovi, ki je posledica razlike gostot i-te komponente, ki prehaja<br />
iz ene tekočine na drugo, je primerljiv s prevodom toplote, ki je posledica razlike<br />
temperatur.<br />
snovna konvekcija<br />
ṁ i,K = v · A · ∆ϱ i ∼ ϱ · v · L 2 (2.3.36)<br />
Prenos snovi zaradi konvekcije je posledica gibanja, enako kot je to pri prenosu toplote.<br />
konvektivni prestop snovi<br />
ṁ i,β = β · A · ∆ϱ i ∼ ϱ · β · L 2 (2.3.37)<br />
Tudi kriterije snovne podobnosti lahko zapišemo – enako kot pri termični podobnosti<br />
– kot razmerja med posameznimi snovnimi tokovi, ki morajo biti za primer dinamične<br />
podobnosti tokov enaka za model in izvedbo.<br />
ṁ i,β<br />
= ϱ · β · L2<br />
ṁ i,D ϱ · D · L = β · L = Sh (2.3.38)<br />
D<br />
kjer je kriterij podobnosti Sh Sherwoodovo število - analogno kot Nußeltovo število<br />
pri konvektivnem prestopu toplote. Pri obravnavi prenosa snovi i-te komponente iz<br />
ene snovi na drugo, npr. hlapenje vode v zrak, kjer imata prevladujoč vpliv konvektivni<br />
prehod in difuzija snovi, je snovna podobnost izpolnjena z enakostjo Sherwoodovih<br />
števil:<br />
[ ]<br />
ṁi,β<br />
ṁ i,D<br />
M<br />
=<br />
[ ]<br />
ṁi,β<br />
ṁ i,D<br />
I<br />
⇒ Sh M = Sh I (2.3.39)<br />
Razmerje med konvekcijo in difuzijo, viskozno in vztrajnostno silo je:<br />
ṁ i,K<br />
· Fη = ϱ · v · L2<br />
ṁ i,D F v ϱ · D · L · η · v · L<br />
ϱ · v 2 · L 2 = ν D<br />
= Sc (2.3.40)<br />
kjer je kriterij podobnosti Sc Schmidtovo število - analogno kot Prandtlovo število<br />
pri termični podobnosti. Tudi Schmidtovo število je snovna lastnost in pomeni pomembno<br />
povezavo med hitrostnim in koncentracijskim poljem dveh sistemov.<br />
Razmerje med konvektivnim prestopom snovi in konvekcijo je:<br />
ṁ i,β<br />
ṁ i,K<br />
= ϱ · β · L2<br />
ϱ · v · L 2 = β v = St II (2.3.41)
64 2 TEORETIČNE OSNOVE<br />
kjer je kriterij podobnosti St II Stantonovo II število - analogno kot Stantonovo število<br />
pri konvektivnem prestopu toplote. Stantonovo II število je mogoče izraziti tudi s<br />
kombinacijo že znanih brezdimenzijskih števil:<br />
St II =<br />
Sh<br />
Re · Sc<br />
(2.3.42)<br />
Preurejena enačba (2.3.34) za prestop snovi na steno okrogle cevi za stacionarni turbuletni<br />
tok:<br />
Sh = α · d<br />
λ<br />
= 0,0235 · (Re 0,8 − 230) · (1,8 · Sc 0,3 − 0,8)<br />
·<br />
[ ( ] ( ) d 3 Sc 0,14<br />
1 + ·<br />
L)2<br />
Sc S<br />
(2.3.43)<br />
Enačba velja pri enakih robnih pogojih kot enačba (2.3.34): Re > 2300 in 0, 7 <<br />
Sc < 1000. Izvrednotena in grafično prikazana je na sliki 2.23.<br />
Slika 2.23: Pregledni diagram za določitev toplotne in snovne prestopnosti za stacionarni<br />
turbuletni tok v okrogli cevi
2.3 PODOBNOST IN DIMENZIJSKA ANALIZA 65<br />
2.3.2 Dimenzijska analiza<br />
V primerjavi s posameznim fizikalnim pojavom, kot je npr. tok tekočine v cevi,<br />
je delovanje stroja mnogo bolj zapleteno, zato v večini primerov ne poznamo vseh<br />
enačb, ki bi omogočale določitev brezdimenzijskih števil. V tem primeru s pridom<br />
uporabimo dimenzijsko analizo, imenovano tudi Buckinghamov ali Π-teorem: vsako<br />
dimenzijsko pravilno enačbo je namreč mogoče zapisati kot povezavo brezdimenzijskih<br />
števil. Dimenzijsko pravilna enačba:<br />
f(x 1 , x 2 , . . . , x n ) = 0 (2.3.44)<br />
kjer so x 1 , x 2 , . . . , x n poljubne dimenzijske veličine, pomembne za opis nekega fizikalnega<br />
pojava, je torej mogoče zapisati tudi v brezdimenzijski obliki:<br />
F (Π 1 , Π 2 , . . . , Π m ) = 0 (2.3.45)<br />
Število karakterističnih brezdimenzijskih števil je m = n − i, pri tem je n število<br />
fizikalnih veličin in i število osnovnih mer. Tako se zmanjša razsežnost opazovanega<br />
fizikalnega pojava na m brezdimenzijskih produktov n fizikalnih veličin. To<br />
omogoča pri raziskavah strojev in naprav bistveno zmanjšanje eksperimentalnega in<br />
teoretičnega dela. Razen tega lahko izmerjene ali izračunane karakteristike različnih<br />
strojev ali naprav, ki so podane v brezdimenzijski obliki, neposredno medsebojno<br />
primerjamo, ker so neodvisne od uporabljenega merskega sistema. Navadno uporabljamo<br />
mednarodni SI-sistem.<br />
Za mehanske probleme je število osnovnih mer i = 3, npr. dolžina L, čas T in masa<br />
M, pri termičnih pa i = 4, potrebna je še temperatura Θ. Posamezne veličine, ki<br />
so značilne za določen fizikalni pojav ali delovanje posameznega stroja ali naprave,<br />
lahko torej zapišemo z osnovnimi merami, najvažnejše so zbrane v preglednici 2.10.<br />
Oznake za osnovne mere so v tem poglavju prilagojene ustaljeni praksi in ne odgovarjajo<br />
siceršnjim oznakam v knjigi.<br />
Za uporabo dimenzijske analize je treba za opazovani fizikalni pojav poznati vse<br />
bistvene fizikalne veličine, saj sicer ne dobimo ustrezne množice brezdimenzijskih<br />
števil. V nadaljevanju je uporaba dimenzijske analize prikazana na treh značilnih<br />
primerih iz energetskega strojništva.
66 2 TEORETIČNE OSNOVE<br />
Tabela 2.10: Mere nekaterih napogosteje uporabljenih fizikalnih veličin v mednarodnem<br />
merskem sistemu SI<br />
Fizikalna veličina Oznaka Osnovna mera<br />
Dolžina<br />
L<br />
Višina H L<br />
Premer<br />
d<br />
Absolutna hrapavost<br />
k<br />
Ploščina A L 2<br />
Prostornina, volumen V L 3<br />
Čas t T<br />
Hitrost v L · T −1<br />
Pospešek<br />
a<br />
Zemeljski pospešek<br />
g<br />
L · T −2<br />
Volumenski tok ˙V L3 · T −1<br />
Vrtilna frekvenca<br />
f<br />
Kotna hitrost ω T −1<br />
Vrtilna frekvenca<br />
n<br />
Kinematična viskoznost ν L 2 · T −1<br />
Snovna difuzivnost<br />
D<br />
Masa m M<br />
Gostota ϱ M · L −3<br />
Dinamična viskoznost η M · L −1 · T −1<br />
Sila F M · L · T −2<br />
Tlak<br />
p<br />
Napetost (mehanska) σ, τ<br />
M · L −1 · T −2<br />
Masni tok ṁ M · T −1<br />
Delo, energija<br />
W<br />
Toplota Q M · L 2 · T −2<br />
Moment sile<br />
M<br />
Moč<br />
P<br />
Toplotni tok<br />
˙Q<br />
M · L 2 · T −3<br />
Temperatura T Θ<br />
Plinska konstanta<br />
R<br />
Izobarna spec. toplota c p L 2 · T −2 · Θ −1<br />
Izohorna spec. Toplota c v<br />
Toplotna prevodnost λ M · L · T −3 · Θ −1
2.3 PODOBNOST IN DIMENZIJSKA ANALIZA 67<br />
Zgled. Tlačne izgube v okrogli cevi<br />
V preglednici 2.11 so podane najvažnejše fizikalne veličine, ki vplivajo na tok tekočine.<br />
Na osnovi teh veličin je mogoče z dimenzijsko analizo določiti najvplivnejša<br />
karakteristična brezdimenzijska števila, ne da bi pri tem uporabili znane enačbe za<br />
tok viskozne tekočine v okrogli cevi.<br />
Tabela 2.11: Fizikalne veličine za tok viskozne tekočine v cevi<br />
Fizikalna Zmanjšanje tlaka, Dinamična Absolutna Premer Hitrost Gostota<br />
veličina deljeno z dolžino viskoznost hrapavost<br />
Eksponent a b c č d e<br />
Oznaka ∆p/L η k d v ϱ<br />
Mera M · L −2 · T −2 M · L −1 · T −1 L L L · T −1 L · T −3<br />
Tok v cevi je podan kot funkcija šestih spremenljivk n = 6:<br />
f( ∆p , η, k, d, v, ϱ) = 0 (2.3.46)<br />
L<br />
Število mer obravnavanega primera je i = 3, fizikalni pojav lahko zato zapišemo kot<br />
funkcijo treh brezdimenzijskih števil m = n − i = 3:<br />
F (Π 1 , Π 2 , Π 3 ) = 0 (2.3.47)<br />
Splošno karakteristično brezdimenzijsko število obravnavanega problema je:<br />
Π . = (M · L −2 · T −2 ) a · (M · L −1 · T −1 ) b ·<br />
· (L) c · (L)č · (L · T −1 ) d · (M · L −3 ) e = (2.3.48)<br />
= M a+b+e · L −2a−b+c+č+d−3e · T −2a−b−d<br />
Ob pogoju, da je Π brezdimenzijsko število, mora biti vsota eksponentov vsake mere<br />
enaka nič, zato dobimo za vsako brezdimenzijsko število Π sistem treh linearnih<br />
enačb s šestimi neznankami. Izbira treh eksponentov je poljubna, preostale tri pa<br />
izračunamo. Število trojic karakterističnih brezdimenzijskih števil je neskončno veliko,<br />
v praksi pa vzamemo vedno kombinacijo Π 1 , Π 2 , Π 3 , ki imajo fizikalni pomen.<br />
V preglednici 2.12 so za vsako karakteristično brezdimenzijsko število podane izbrane<br />
vrednosti eksponentov a, b in c ter izračunane vrednosti preostalih treh eksponentov<br />
po enačbah:<br />
a + b + e = 0<br />
−2 · a − b + c + č + d − 3 · e = 0<br />
−2 · a − b − d = 0
68 2 TEORETIČNE OSNOVE<br />
Če za Π 1 izberemo vrednosti eksponentov: a = 1, b = 0 in c = 0, lahko izračunamo<br />
preostale tri eksponente:<br />
⎫<br />
1 + 0 + e = 0 ⎪⎬ č = 1<br />
−2 − 0 + 0 + č + d − 3 · e = 0<br />
−2 − 0 − d = 0<br />
⎪⎭ ⇔ d = −2<br />
(2.3.49)<br />
e = −1<br />
Enako izračunamo vrednosti eksponentov za brezdimenzijska števila Π 2 in Π 3 .<br />
Tabela 2.12: Matriki izbranih in izračunanih eksponentov za tok viskozne tekočine v<br />
cevi<br />
Izbrani eksponenti<br />
a b c<br />
Π 1 1 0 0<br />
Π 2 0 1 0<br />
Π 3 0 0 1<br />
Izračunani eksponenti<br />
č d e<br />
Π 1 1 −2 −1<br />
Π 2 −1 −1 −1<br />
Π 3 −1 0 0<br />
Karakteristična brezdimenzijska števila Π 1 , Π 2 in Π 3 , ki opisujejo viskozni tok v<br />
cevi, so tako določena:<br />
Π 1 =<br />
( ) ∆p a<br />
· (η) b · (k) c · (d)č · (v) d · (ϱ) e = ∆p<br />
L<br />
ϱ · v 2 · d<br />
L = Eu · d<br />
L (2.3.50)<br />
Brezdimenzijsko število Π 1 je za poznano geometrijo cevi sorazmerno Eulerjevem<br />
številu.<br />
( ) ∆p a<br />
Π 2 = · (η) b · (k) c · (d)č · (v) d · (ϱ) e η<br />
=<br />
L<br />
ϱ · v · d = 1 (2.3.51)<br />
Re<br />
Brezdimenzijsko število Π 2 je obratno sorazmerno Reynoldsovem številu.<br />
Π 3 =<br />
( ) ∆p a<br />
· (η) b · (k) c · (d)č · (v) d · (ϱ) e = k L<br />
d<br />
(2.3.52)<br />
Brezdimenzijsko število Π 3 je poznano kot relativna hrapavost.
2.3 PODOBNOST IN DIMENZIJSKA ANALIZA 69<br />
Tok viskozne tekočine v cevi je funkcija treh karakterističnih brezdimenzijskih števil:<br />
( ∆p<br />
F (Π 1 , Π 2 , Π 3 ) = F<br />
ϱ · v 2 · d<br />
L , 1<br />
Re , k )<br />
= 0 (2.3.53)<br />
d<br />
Izguba tlaka je navadno prikazana kot funkcijska odvisnost od preostalih dveh brezdimenzijskih<br />
števil Π 2 in Π 3 :<br />
∆p = φ(Π 2 , Π 3 ) · ϱ · v 2 · L<br />
d = λ 2 · ϱ · v2 · L<br />
d<br />
(2.3.54)<br />
kjer funkcijo φ(Π 2 , Π 3 ) izrazimo s koeficientom tekočinskega trenja λ, ustrezno<br />
enačbi (2.1.24):<br />
φ(Π 2 , Π 3 ) = φ(Re, k/d) =<br />
λ(Re, k/d)<br />
2<br />
(2.3.55)<br />
Zgled. Brezdimenzijske karakteristike hidravličnih turbinskih strojev<br />
V preglednici 2.13 so podane najvažnejše fizikalne veličine, ki vplivajo na delovanje<br />
hidravličnih turbinskih strojev, kot so turbinske črpalke in vodne turbine. Enako kot v<br />
prejšnjem zgledu bodo z dimenzijsko analizo določena najvplivnejša karakteristična<br />
brezdimenzijska števila. Za izračun privzamemo nestisljivost tekočine in konstantno<br />
temperaturo toka.<br />
Tabela 2.13: Fizikalne veličine za hidravlične turbinske stroje<br />
Fizikalna Moč na Spec. Izkoristek Dinamična Volumenski Gostota Vrtilna Premer<br />
veličina gredi energija viskoznost tok frekvenca rotorja<br />
Eksponent a b c č d e f g<br />
Oznaka P gH η νϱ ˙V ϱ n d<br />
Mera M · L 2 · T −3 L 2 · T −2 - M · L −1 · T −1 L 3 · T −1 M · L −3 T −1 L<br />
Karakteristika hidravličnega turbinskega stroja je podana kot funkcija osmih spremenljivk<br />
n = 8:<br />
f(P, gH, η, νϱ, ˙V , ϱ, n, d) = 0 (2.3.56)<br />
Število mer obravnavanega primera je i = 3, primer lahko zato zapišemo kot funkcijo<br />
petih brezdimenzijskih števil m = n − i = 5:<br />
F (Π 1 , Π 2 , Π 3 , Π 4 , Π 5 ) = 0 (2.3.57)
70 2 TEORETIČNE OSNOVE<br />
Tabela 2.14: Matriki izbranih in izračunanih eksponentov za hidravlične turbinske<br />
stroje<br />
Izbrani eksponenti<br />
Izračunani eksponenti<br />
a b c č d<br />
Π 1 1 0 0 0 0<br />
Π 2 0 1 0 0 0<br />
Π 3 0 0 1 0 0<br />
Π 4 0 0 0 1 0<br />
Π 5 0 0 0 0 1<br />
e f g<br />
Π 1 −1 −3 −5<br />
Π 2 0 −2 −2<br />
Π 3 0 0 0<br />
Π 4 −1 −1 −2<br />
Π 5 0 −1 −3<br />
Splošno karakteristično brezdimenzijsko število obravnavanega problema je:<br />
Π . = (M · L 2 · T −3 ) a · (L 2 · T −2 ) b · (1) c · (M · L −1 · T −1 )č ·<br />
· (L 3 · T −1 ) d · (M · L −3 ) e · (T −1 ) f · (L) g = (2.3.58)<br />
= M a+č+e · L 2a+2b−č+3d−3e+g · T −3a−2b−č−d−f<br />
V preglednici 2.14 so za vsako karakteristično brezdimenzijsko število podane izbrane<br />
vrednosti eksponentov a, b, c, č in d ter izračunane vrednosti preostalih treh<br />
eksponentov po že znanem načinu.<br />
Karakteristična brezdimenzijska števila Π 1 , Π 2 , Π 3 , Π 4 , Π 5 in Π 6 , ki opisujejo delovanje<br />
hidravličnih turbinskih strojev, so tako določena:<br />
Π 1 = (P ) a · (gH) b · (η) c · (νϱ)č · ( ˙V ) d · (ϱ) e · (n) f · (d) g =<br />
P<br />
=<br />
ϱ · n 3 · d 5 (2.3.59)<br />
Brezdimenzijsko število Π 1 je poznano kot močnostno število.<br />
Π 2 = (P ) a · (gH) b · (η) c · (νϱ)č · ( ˙V ) d · (ϱ) e · (n) f · (d) g =<br />
gH<br />
=<br />
n 2 · d 2 (2.3.60)<br />
Brezdimenzijsko število Π 2 je poznano kot energijsko število.<br />
Π 3 = (P ) a · (gH) b · (η) c · (νϱ)č · ( ˙V ) d · (ϱ) e · (n) f · (d) g =<br />
= η (2.3.61)<br />
Brezdimenzijsko število Π 3 je preprosto izkoristek turbinskega stroja.<br />
Π 4 = (P ) a · (gH) b · (η) c · (νϱ)č · ( ˙V ) d · (ϱ) e · (n) f · (d) g =<br />
νϱ<br />
=<br />
ϱ · n · d 2 = ν<br />
(n · d) · d = 1<br />
(2.3.62)<br />
Re u
2.3 PODOBNOST IN DIMENZIJSKA ANALIZA 71<br />
Brezdimenzijsko število Π 4 je obratno sorazmerno z obodnim Reynoldsovim številom<br />
(indeks u). Izkazalo se je, da je tok v hidravličnih turbinskih strojih vedno močno<br />
turbulenten, zato viskoznost tekočine ne vpliva bistveno na delovanje stroja. Obodno<br />
Reynoldsovo število lahko v večini tehničnih primerov zanemarimo.<br />
Π 5 = (P ) a · (gH) b · (η) c · (νϱ)č · ( ˙V ) d · (ϱ) e · (n) f · (d) g =<br />
˙V<br />
=<br />
n · d 3 (2.3.63)<br />
Brezdimenzijsko število Π 5 je poznano kot pretočno število.<br />
Karakteristika hidravličnih turbinskih strojev je podana kot funkcija štirih brezdimenzijskih<br />
števil:<br />
F<br />
(<br />
P<br />
ϱ · n 3 · d 5 ,<br />
gH<br />
n 2 · d 2 , η,<br />
)<br />
˙V<br />
n · d 3 = 0 (2.3.64)<br />
Karakteristike hidravličnih turbinskih strojev so navadno prikazane kot funkcijska<br />
odvisnost pretočnega števila od močnostnega in energijskega števila ter od izkoristka<br />
turbinskega stroja:<br />
P<br />
ϱ · n 3 · d 5 = F ∗ ( ˙V<br />
n · d 3 )<br />
gH<br />
n 2 · d 2 = F ∗∗ ( ˙V<br />
n · d 3 )<br />
η = F ∗∗∗ ( ˙V<br />
n · d 3 )<br />
(2.3.65)<br />
(2.3.66)<br />
(2.3.67)<br />
Zgled. Brezdimenzijske karakteristike toplotnih turbinskih strojev<br />
V preglednici 2.15 so podane najvažnejše fizikalne veličine, ki vplivajo na delovanje<br />
toplotnih turbinskih strojev, kot so turbinski kompresorji, parne in plinske turbine.<br />
Tudi pri toplotnih turbinskih strojih je – enako kot pri hidravličnih strojih –<br />
tok skozi stroj močno turbulenten, zato viskoznost tekočine ne vpliva bistveno na<br />
delovanje stroja. Obodno Reynoldsovo število lahko zanemarimo. Nasprotno od<br />
hidravličnih turbinskih strojev moramo pri toplotnih turbinskih strojih upoštevati stisljivost<br />
tekočine: volumenski tok je potrebno nadomestiti z masnim. Pri dimenzijski<br />
analizi je temperaturo primerno zapisati kot zmnožek plinske konstante s temperaturo<br />
(R T ), da se izognemo pisanju mere za temperaturo Θ.
72 2 TEORETIČNE OSNOVE<br />
Tabela 2.15: Fizikalne veličine za toplotne turbinske stroje<br />
Fizikalna Masni Vrtilna Izstopni Izstopna Eksponent Premer Vstopni Vstopna<br />
veličina tok frekv. tlak temperatura izentrope rotorja tlak temperatura<br />
Eksponent a b c č d e f g<br />
Oznaka ṁ n p 2 R T 2 κ d p 1 R T 1<br />
Mera M · T −1 T −1 M · L −1 · T −2 L 2 · T −2 - L M · L −1 · T −2 L 2 · T −2<br />
Karakteristika toplotnega turbinskega stroja je podana kot funkcija osmih spremenljivk<br />
n = 8:<br />
f(ṁ, n, p 2 , RT 2 , κ, d, p 1 , RT 1 ) = 0 (2.3.68)<br />
Število mer obravnavanega primera je i = 3, primer lahko zato zapišemo kot funkcijo<br />
petih brezdimenzijskih števil m = n − i = 5:<br />
F (Π 1 , Π 2 , Π 3 , Π 4 , Π 5 ) = 0 (2.3.69)<br />
Splošno karakteristično brezdimenzijsko število obravnavanega problema je:<br />
Π . = (M · T −1 ) a · (T −1 ) b · (M · L −1 · T −2 ) c · (L 2 · T −2 )č ·<br />
· (1) d · (L) e · (M · L −1 · T −2 ) f · (L 2 · T −2 ) g = (2.3.70)<br />
= M a+c+f · L −c+2č+e−f+2g · T −a−b−2c−2č−2f−2g<br />
Tabela 2.16: Matriki izbranih in izračunanih eksponentov za toplotne turbinske stroje<br />
Izbrani eksponenti<br />
a b c č d<br />
Π 1 1 0 0 0 0<br />
Π 2 0 1 0 0 0<br />
Π 3 0 0 1 0 0<br />
Π 4 0 0 0 1 0<br />
Π 5 0 0 0 0 1<br />
Izračunani eksponenti<br />
e f g<br />
Π 1 −2 −1 0,5<br />
Π 2 1 0 −0,5<br />
Π 3 0 −1 0<br />
Π 4 0 0 −1<br />
Π 5 0 0 0
2.3 PODOBNOST IN DIMENZIJSKA ANALIZA 73<br />
V preglednici 2.16 so za vsako karakteristično brezdimenzijsko število podane izbrane<br />
vrednosti eksponentov a, b, c, č in d ter izračunane vrednosti preostalih treh eksponentov<br />
po že znanem načinu.<br />
Karakteristična brezdimenzijska števila Π 1 , Π 2 , Π 3 , Π 4 in Π 5 , ki opisujejo delovanje<br />
toplotnih turbinskih strojev, so tako določena:<br />
Π 1 = (ṁ) a · (n) b · (p 2 ) c · (RT 2 )č · (κ) d · (d) e · (p 1 ) f · (RT 1 ) g =<br />
= ṁ · √RT<br />
1<br />
p 1 · d 2 (2.3.71)<br />
Brezdimenzijsko število Π 1 je poznano kot pretočno število.<br />
Π 2 = (ṁ) a · (n) b · (p 2 ) c · (RT 2 )č · (κ) d · (d) e · (p 1 ) f · (RT 1 ) g =<br />
= n · d √ RT1<br />
(2.3.72)<br />
Brezdimenzijsko število Π 2 je poznano kot brezdimenzijska vrtilna frekvenca.<br />
Π 3 = (ṁ) a · (n) b · (p 2 ) c · (RT 2 )č · (κ) d · (d) e · (p 1 ) f · (RT 1 ) g =<br />
= p 2<br />
p 1<br />
(2.3.73)<br />
Brezdimenzijsko število Π 3 pomeni povečanje (zmanjšanje) tlaka pri kompresiji (ekspanziji)<br />
v kompresorju (turbini) med vstopom in izstopom turbinskega stroja in je<br />
primerljivo z energijskim številom pri hidravličnih turbinskih strojih.<br />
Π 4 = (ṁ) a · (n) b · (p 2 ) c · (RT 2 )č · (κ) d · (d) e · (p 1 ) f · (RT 1 ) g =<br />
= RT 2<br />
RT 1<br />
= T 2<br />
T 1<br />
(2.3.74)<br />
Brezdimenzijsko število Π 4 pomeni povečanje (zmanjšanje) temperature pri kompresiji<br />
(ekspanziji) v kompresorju (turbini) med vstopom in izstopom turbinskega stroja<br />
in je, fizikalno gledano, povezano z notranjim izkoristkom turbinskega stroja.<br />
Π 5 = (ṁ) a · (n) b · (p 2 ) c · (RT 2 )č · (κ) d · (d) e · (p 1 ) f · (RT 1 ) g =<br />
= κ (2.3.75)<br />
Brezdimenzijsko število Π 5 je preprosto eksponent izentrope. Karakteristika toplotnega<br />
turbinskega stroja je podana kot funkcija petih brezdimenzijskih števil:<br />
(ṁ √ · RT1 n · d<br />
F<br />
p 1 · d 2 , √ , p 2<br />
, T )<br />
2<br />
, κ = 0 (2.3.76)<br />
RT1 p 1 T 1
74 2 TEORETIČNE OSNOVE<br />
Karakteristike toplotnih turbinskih strojev so navadno prikazane kot funkcijska odvisnost<br />
pretočnega števila od brezdimenzijske vrtilne frekvence, razmerja tlakov in<br />
temperatur ter od eksponenta izentrope.<br />
(ṁ √ )<br />
n · d<br />
· R<br />
√ = F ∗ T1<br />
R T1 p 1 · d 2 (2.3.77)<br />
(ṁ √ )<br />
p 2<br />
· R<br />
= F ∗∗ T1<br />
p 1 p 1 · d 2 (2.3.78)<br />
(ṁ √ )<br />
T 2<br />
· R<br />
= F ∗∗∗ T1<br />
T 1 p 1 · d 2 (2.3.79)<br />
(ṁ √ )<br />
· R<br />
η = F ∗∗∗ T1<br />
p 1 · d 2 (2.3.80)<br />
(ṁ √ )<br />
· R<br />
κ = F ∗∗∗∗ T1<br />
p 1 · d 2 (2.3.81)<br />
Pretočno število je odvisno od razmerja temperatur, mogoče pa je dokazati, da lahko<br />
razmerje temperatur nadomestimo z izkoristkom. Pri večini praktičnih primerov je<br />
odvisnost pretočnega števila od eksponenta izentrope majhna in je navadno vsebovana<br />
v izkoristku.