24.10.2014 Views

Teoretične osnove

Teoretične osnove

Teoretične osnove

SHOW MORE
SHOW LESS

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.

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!