19.04.2014 Views

Φ Φ Φ Φ - FSB

Φ Φ Φ Φ - FSB

Φ Φ Φ Φ - FSB

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.

Mirko Tadić<br />

Termodinamika<br />

6 Otvoreni sustavi<br />

Premda je teorija termodinamike postavljena na modelu zatvorenog sustava, dobivene<br />

termodinamičke relacije mogu se uz dodatne pretpostavke primijeniti i na tehnički značajnije<br />

otvorene sustave. Njih karakterizira protok mase, ṁ (kg/s), kroz granične plohe sustava, koji<br />

je definiran kao fiksni prostor, tzv. kontrolni volumen. Na slici 6.1 prikazan je tipičan model<br />

otvorenog sustava. Granice sustava određene su stijenkom cijevi kroz koju protječe radni<br />

medij – idealni plin.<br />

Ḣ ̇ +<br />

̇<br />

1<br />

= U1<br />

p1V<br />

1<br />

1<br />

Φ snaga P Φ 2<br />

dV<br />

Ḣ ̇ +<br />

̇<br />

2<br />

= U2<br />

p2V<br />

2<br />

w 1 w 2<br />

z 1 z 2<br />

Φ toplinski tok<br />

Φ<br />

Slika 6.1 Karakterističan model otvorenog strujanja<br />

Protok mase ṁ (kg/s), odnosno količine Ṅ (kmol/s), osnovna je karakteristika otvorenih<br />

sustava. Zbog toga svi članovi u jednadžbi održanja energije (I. Zakon) imaju dimenziju J/s =<br />

W. Umjesto o toplini govorimo o toplinskom toku, Φ 12 , a umjesto o radu govorimo o snazi,<br />

P 12 . Početno stanje radnog medija je ulazno stanje (1), konačno stanje je izlazno stanje (2).<br />

Ulazno i izlazno stanje opisani su s temperaturom T ili ϑ, tlakom p i protočnim volumenom<br />

V̇ m 3 /s. Odnos protočnog volumena i protočne mase ṁ kg/s je specifični volumen v:<br />

V̇<br />

v = , m 3 /kg. (6.1)<br />

ṁ<br />

Protočni volumen je povezan s brzinom strujanja w (m/s) i površinom presjeka strujanja A<br />

(m 2 ):<br />

V ̇ = wA , m 3 /s. (6.2)<br />

Energija koju radni medij unosi, odnosno iznosi iz kontrolnog volumena, sastoji se iz dva<br />

U̇ = mu ̇ = U̇<br />

T , J/s, te<br />

dijela: unutarnje energije, srazmjerne protočnoj masi i temeperaturi, ( )<br />

snazi strujanja, srazmjernoj produktu tlaka i protočnog volumena, pV<br />

̇ , (N/m 2 )(m 3 /s) = J/s.<br />

Ova dva oblika su energije vezane sa stanjem radnog medija, a objedinjena su pojmom<br />

entalpije Ḣ , koja je također veličina stanja.<br />

Po definiciji je specifična entalpija:<br />

h = u + pv , J/kg. (6.3)<br />

Množenjem jednadžbe (6.3) s protočnom masom ṁ slijedi da je:<br />

H ̇ = mh ̇ = mu ̇ + pmv ̇ = U̇<br />

+ pV̇<br />

, J/s. (6.4)<br />

64


Mirko Tadić<br />

Termodinamika<br />

Kako je entalpija veličina stanja to se ona matematički klasificira kao totalni diferencijal, što<br />

znači da je promjena entalpije:<br />

( − )<br />

∆Ḣ = ̇ , (6.5)<br />

12<br />

m h2<br />

h1<br />

neovisna o procesu, tj. načinu promjene stanja radnog medija od ulaza do izlaza iz otvorenog<br />

sustava. Zbog toga se kaže da je entalpija konzervativno svojstvo, poput mase. Uz to, entalpija<br />

je ekstenzivna veličina, jer njena numerička vrijednost ovisi o masi.<br />

Apsolutna vrijednost specifične entalpije h nema značaja pri opisu procesa, već samo<br />

promjena njene vrijednosti tijekom procesa, ∆h 12 . Stoga je sasvim svejedno koje se toplinsko<br />

stanje smatra nultim stanjem entalpije (označeno indeksom 0), tj. stanjem pri kojem je<br />

entalpija h 0 = 0. Do izražaja dolazi samo razlika entalpija stanja (1) i (2):<br />

( h − h ) − ( h − h ) = h − h = c ( T − T ) = C ( ϑ − )<br />

∆h , J/kg. (6.6)<br />

12<br />

=<br />

2 0 1 0 2 1 p 2 1 p 2<br />

ϑ1<br />

Izborom nultog stanja entalpije s h 0 = 0, tom stanju je istovremeno pripisana određena<br />

vrijednost specifične unutarnje energije u 0 , sukladno relaciji:<br />

u = = , J/kg. (6.7)<br />

0<br />

h0<br />

+ p0v0<br />

p0v0<br />

To nema utjecaja na oblik računa, jer do izražaja dolazi samo razlika specifične unutarnje<br />

energije, ∆u 12 , za koju vrijedi:<br />

( u − u ) − ( u − u ) = u − u = c ( T − T ) = C ( ϑ − )<br />

∆u , J/kg. (6.8)<br />

12<br />

=<br />

2 0 1 0 2 1 v 2 1 v 2<br />

ϑ1<br />

Vidi se da je stvarna vrijednost u 0 potpuno nevažna.<br />

Na temelju jednadžbi (6.6) i (6.8) nastaju relacije koje su korisne za proračun otvorenih<br />

sustava:<br />

( − h ) = mc ̇ ( T − T ) = Ṅ<br />

C ( ϑ − )<br />

∆Ḣ ̇<br />

, J/s = W. (6.9)<br />

12<br />

= m h2<br />

1 p 2 1<br />

p 2<br />

ϑ1<br />

( − u ) = mc ̇ ( T − T ) = Ṅ<br />

C ( ϑ − )<br />

∆U̇ ̇<br />

, J/s = W. (6.10)<br />

12<br />

= m u2<br />

1 v 2 1 v 2<br />

ϑ1<br />

Primjenom jednadžbe (6.4) na stanja (1) i (2) dobijaju se dvije jednadžbe, pa se razlika stanja<br />

(2) prema stanju (1) može izraziti kao:<br />

∆ Ḣ<br />

̇ ̇ ̇ . (11)<br />

12<br />

= ∆U12<br />

+ p2V2<br />

− p1V<br />

1<br />

Promjena entalpije obuhvaća i promjenu unutarnje energije i razliku energije strujanja s<br />

istovremenim promjenama volumena. Radi jasnoće možemo promotriti diferencijalni oblik<br />

jednadžbe (6.11):<br />

( pV̇<br />

) = dU̇<br />

+ pdV̇<br />

V̇<br />

dp<br />

dH<br />

̇ = dU̇<br />

+ d<br />

+ . (6.12)<br />

65


Mirko Tadić<br />

Termodinamika<br />

Značenje dU<br />

̇ (W) jasno je od ranije, s jedinom razlikom što se preko protočne mase uvela<br />

jedinica vremena. Član pdV<br />

̇ je očito energija (po vremenu) koju radni medij uzima (daje) za<br />

porast (smanjenje) svog volumen. Dakle taj član je snaga, povezana s promjenom volumena.<br />

Posljednji član, V̇ dp , je snaga čistog strujanja bez volumenske promjene, tj. samo kao<br />

transport u smjeru pada tlaka.<br />

Uvođenjem pojma entalpije nestaje potreba da sve te efekte razmatramo pojedinačno!<br />

Stanje radnog medija na ulaznom (1) i na izlaznom (2) presjeku nije jednoznačno, jer su<br />

stanja materijalnih čestica po istom presjeku različita i po brzinama i po temperaturi!<br />

Uzimanje u obzir tih činjenica vodi jako kompliciranom proračunu. Sukladno dosadašnjem<br />

principu idealiziranja fizikalnog modela možemo pretpostaviti da je stanje materije po<br />

presjeku, tj unutar proizvoljno malog volumena, dV, jedinstveno. Time će se proračun bitno<br />

pojednostaviti, ali uz izvjesnu netočnost rezultata koja se može ustanoviti eksperimentalno.<br />

Po potrebi, uvijek se može primijeniti kompleksniji pristup.<br />

U narednim razmatranjima polazi se od pretpostavke jedinstvenog stanja po presjeku.<br />

U odnosu na fizikalni model prikazan na slici 6.1 vidljivo je da zbog toplinske i mehaničke<br />

interakcije s okolišem može doći ne samo do promjene entalpije, ∆ Ḣ12<br />

, već i do promjene<br />

kinetičke ∆ Ė k , 12<br />

i potencijalne ∆ Ė p, 12<br />

energije radnog medija. Primjena I. Zakona na takav<br />

fizikalni model rezultira jednadžbom:<br />

Φ − P = ∆Ḣ<br />

+ ∆Ė<br />

k<br />

+ ∆Ė<br />

, W. (I. Zakon za otvoreni sustav). (6.13)<br />

12 12 12 , 12 p,<br />

12<br />

Promjena kinetičke energije može se odrediti iz jednadžbe:<br />

2 2<br />

w2<br />

− w1<br />

∆E ̇ k , 12<br />

= ṁ<br />

2<br />

, W, (6.14)<br />

pri čemu se za određivanje brzina može koristiti veza s protočnim volumenom prema<br />

jednadžbi (6.2):<br />

V̇<br />

w = , m/s. (6.15)<br />

A<br />

Promjena potencijalne energije postoji, ako postoji razlika u poziciji presjeka, z,<br />

( z − )<br />

∆ ̇ ̇ , W. (6.16)<br />

E p,<br />

12<br />

= mg<br />

2<br />

z1<br />

Toplinska interakcija s okolišem, u obliku toplinskog toka, Φ 12 (W), odvija se pri konačnoj<br />

razlici temperatura radnog medija (RM) i toplinskih sudionika (TS) u okolišu, tj.<br />

ireverzibilno.<br />

Mehanička interakcija s okolišem, u obliku snage P 12 (W), prepoznaje se kao neposredni<br />

dodir radnog medija s pomičnim mehaničkim sudionikom (MS) i smatra se reverzibilnom, tj.<br />

bez gubitaka.<br />

Razmotrit će se odvojeno dvije vrste procesa radnog medija:<br />

• ravnotežne promjene stanja (politropske),<br />

• neravnotežne promjene stanja, kod kojih nisu ispunjeni uvjeti unutarnje toplinske<br />

mehaničke ravnoteže.<br />

66


Mirko Tadić<br />

Termodinamika<br />

Politrope otvorenih sustava<br />

Proces u nekim važnim tehničkim uređajima može se klasificirati kao proces idealnog<br />

radnog medija u otvorenom sustavu. Kao primjer spomenimo zračni kompresor koji je<br />

pokretan električkim motorom. Radnom mediju, zraku, pripisuje se idealno ponašanje budući<br />

da su stvarni tlak i temperatura tijekom procesa daleko od uvjeta pretvorbe u kapljevito stanje<br />

(sudionika u zraku). Kao sastav zraka obično se uzima molni sastav s 21% O 2 i 79% N 2 .<br />

Standardni sastav zraka ponešto je drugačiji. Za svojstva takvog zraka postoje numerički<br />

podaci u odgovarajućim tablicama, pa se pri proračunu zrak smatra jednim plinom.<br />

Radi lakšeg razumijevanja opisa rada stapnog kompresora prikazan je na slici 6.2 presjek<br />

cilindra. Otvaranje i zatvaranje usisnog (UV) i tlačnog (TV) ventila određeno je razlikom<br />

tlaka zraka u cilindru i tlakova u usisnom (p u ), odnosno tlačnom (p t ) vodu.<br />

tlačni vod<br />

štetni prostor<br />

cilindar<br />

p t<br />

TV<br />

STAPAJICA<br />

p u<br />

usisni vod<br />

UV<br />

STAP<br />

s<br />

stapaj<br />

Slika 6.2 Stapni kompresor<br />

Realni i teorijski proces zraka u kompresoru prikazani su u p-V dijagramu na slici 6.3.<br />

Na realni proces utječu dva faktora:<br />

• postojanje štetnog prostora u kojem uvijek ostaje dio komprimiranog zraka, koji se pri<br />

usisu miješa s usisanim zrakom,<br />

• inertnost usisnog i tlačnog ventila zbog koje pri usisu nastaje potlak, a pri ispuhu<br />

pretlak – u odnosu na tlakove u usisnom i tlačnom vodu.<br />

Teorijski proces aproksimira usis kao punjenje cilindra sa zrakom čiji se tlak i temperatura pri<br />

tome ne mijenjaju, tj. zrak je istog stanja kao u usisnom vodu. Stoga je sekvenca usisa čisti<br />

transport vanjskog zraka, koji je uzrokovan pokretom stapa. Stanje plina nije pod utjecajem<br />

bilo kakve mehaničke i toplinske interakcije iz okoliša, pa zbog toga usis nije politropska<br />

promjena. Sve točke na liniji usisa su istog stanja kao i točka označena s (1).<br />

Nakon završetka usisa, stap kreće u suprotnom smjeru, a oba ventila su zatvorena.<br />

Smanjivanjem volumena zrak se komprimira pri čemu raste tlak. Sekvenca kompresije smatra<br />

se ravnotežnom promjenom, tj. politropom, pV n = konst.<br />

Njen eksponent n ovisi o izmjeni<br />

topline, tj. hlađenju tijekom kompresije, Φ 12 . Kompresija završava u trenutku kada je tlak u<br />

cilindru dostigao vrijednost tlaka u tlačnom vodu.<br />

Rad ventila kompresora nije dirigiran nikakvim drugim mehanizmom, već samo razlikama<br />

tlaka, koje se u teorijskom procesu smatraju zanemarivo malenim, a reakcije ventila<br />

trenutnim.<br />

67


Mirko Tadić<br />

Termodinamika<br />

p<br />

2<br />

Realni proces<br />

Teorijski proces<br />

p 1<br />

p 2<br />

V 1<br />

dp<br />

{ P } okret<br />

P = 12<br />

n0<br />

12<br />

{ P } okret<br />

= −∫V( p)<br />

12<br />

dp<br />

2<br />

1<br />

V 2<br />

n = konst.<br />

V<br />

1<br />

V, m 3 /okret<br />

motor<br />

ispuh<br />

usis<br />

kompresija<br />

ZRAK<br />

Slika 6.3 Realni i teorijski proces u kompresoru<br />

Tlačni ventil se podiže s ispušnog otvora i zrak se istiskuje u tlačni vod, prema<br />

odgovarajućem spremniku komprimiranog zraka. Sekvenca ispuha nije politropska promjena,<br />

već transport zraka stanja (2), kakvo je dostignuto na kraju kompresije.<br />

Procesi koji su prikazani na slici 6.3 odnose se na jedan puni okret vratila kompresora. Zbog<br />

toga su na apscisi volumeni po okretu, V (m 3 /okret), pa se za zadani broj okretaja n 0 (okret/s)<br />

mora izračunati protočni volumen V̇ (m 3 /s) pomoću relacije:<br />

V̇ = n 0<br />

V = mv ̇ , (m 3 /s). (6.17)<br />

Snaga motora za pogon kompresora troši se tijekom kompresije i ispuha, a dio snage se<br />

dobiva tijekom usisa efektom inercije. Mehanička interakcija prema okolišu, tj. zatvorenom<br />

prostoru kućišta u kome se nalazi vratilo kompresora (na slici 6.3 opisan kao zrak), se<br />

poništava budući da se stap giba jednako u dva suprotna smjera, a gubici zbog trenja<br />

zanemaruju.<br />

Razultantni utrošak snage za pogon kompresora pri teorijskom procesu prikazan je u p-V<br />

dijagramu kao površina omeđena linijama usisa, kompresije i ispuha. Ta se površina može<br />

prikazati kao zbroj niza malih površina koje imaju smisao diferencijalno male utrošene snage,<br />

Vdp = dP. Stoga za utrošak snage po jednom okretu vratila vrijedi jednadžba:<br />

2<br />

= −∫<br />

1<br />

{ } V ( p)<br />

P okret<br />

, J/okret. (6.18)<br />

12<br />

dp<br />

Negativan predznak u jednadžbi (6.18) uveden je iz slijedećih razloga. Matematički se dp<br />

uvijek izražava prema + p smjeru, pa bi produkt Vdp imao pozitivan smisao. Kako postojeći<br />

dogovor o smislu mehaničkog rada (snage), propisuje negativan predznak svakom utrošenom<br />

radu (snazi) to se taj smisao mora osigurati i u jednadžbi (6.18).<br />

68


Mirko Tadić<br />

Termodinamika<br />

Za ukupni utrošak snage pri n 0 (okret/s) vrijedi relacija:<br />

2<br />

{ P } = −n<br />

V ( p) dp = − V ( p) dp = − mv ̇ ( p) dp = −m<br />

v( p)<br />

P12 = n<br />

̇ ̇<br />

0 12 okret 0<br />

dp , W. (6.19)<br />

∫<br />

1<br />

2<br />

∫<br />

1<br />

Lako je uočljivo, prama slici 6.3, da je površina koja odgovara utrošku snage jednaka površini<br />

koja nastaje projekcijom politropske kompresije (1)-(2) na ordinatnu os p.<br />

Za politropu vrijedi jednadžba pv n = konst.<br />

na osnovu koje se diferenciranjem dobiva:<br />

n<br />

v dp + npv<br />

iz čega slijedi<br />

n<br />

2<br />

∫<br />

1<br />

− 1<br />

dv = 0 , (6.20)<br />

vdp = −npdv . (6.21)<br />

Za cijelu politropu između (1) i (2) vrijedi:<br />

2<br />

1<br />

2<br />

( p) dp −n<br />

p( v) dv = −nw12<br />

∫ v = ∫<br />

, (6.22)<br />

1<br />

gdje je w 12 (J/kg) specifični mehanički rad politrope čija se površina dobija projekcijom na<br />

apscisnu os. Sukladno tome, može se jednadžba (6.19) pisati u obliku:<br />

2<br />

12<br />

−ṁ v<br />

= W<br />

P<br />

1<br />

( p) dp = nmw ̇<br />

12<br />

n ̇<br />

12<br />

= ∫ , W. (6.23)<br />

2<br />

∫<br />

1<br />

P<br />

P<br />

3<br />

= −∫<br />

1<br />

2<br />

= −∫<br />

1<br />

p<br />

N/m 2<br />

3<br />

p<br />

2<br />

3<br />

−∫ V̇<br />

p dp = P12<br />

= nẆ<br />

12<br />

= n<br />

V̇<br />

1<br />

( p)<br />

n<br />

2 13<br />

p 2<br />

V̇<br />

( p)<br />

n 12<br />

1<br />

13<br />

dp<br />

12<br />

dp<br />

p 1<br />

V 2 =V 3<br />

2<br />

( ) ∫ p( V̇<br />

)<br />

V 1<br />

Ẇ<br />

2<br />

= ∫<br />

1<br />

( ̇)<br />

̇<br />

12<br />

p V dV<br />

V, m 3 /okret<br />

1<br />

dV̇<br />

Slika 6.4 Prikazi snage i mehaničkog rada politropa u p – V dijagramu<br />

Uočimo da se u jednadžbi (6.23) nalazi mehanička snaga politrope Ẇ<br />

12<br />

, zbog utjecaja protoka<br />

mase ṁ (kg/s).<br />

Ako se politropi pridruže usis i ispuh, tada se množenjem snage, Ẇ 12<br />

, s eksponentom<br />

politrope n dobiva snaga, P<br />

12<br />

, takvog otvorenog politropskog procesa.<br />

69


Mirko Tadić<br />

Termodinamika<br />

To znači da se za otvorene politropske procese mogu koristiti sve ranije izvedene relacije za<br />

mehanički rad zatvorenih sustava, ali se pri tome umjesto mase m (kg) treba uvrstiti protočna<br />

masa ṁ (kg/s), kako bi se dobila snaga Ẇ 12<br />

. Jedini izuzetak je otvoreni sustav s izohorom<br />

( V ̇ = konst.<br />

) kod koje je eksponent n = ± ∞, a mehanička snaga Ẇ 0 , tako da je umnožak<br />

nW ̇ 12<br />

neodređen. U tom slučaju treba koristiti jednadžbu (6.23) koja omogućava dobijanje<br />

rezultata:<br />

2<br />

2<br />

12<br />

= −m∫<br />

vdp = −mv<br />

̇ dp = −V̇<br />

∫<br />

2<br />

p1<br />

1<br />

1<br />

P<br />

( p − )<br />

̇ , W, (izohorni otvoreni proces). (6.24)<br />

Budući da su usis i ispuh samo transportne sekvence procesa, to pri tome nema izmjene<br />

topline. Jedina izmjena topline postoji, eventualno, samo tijekom politropske promjene koja<br />

se odvija pri zatvorenim ventilima. Stoga se toplinski tok, Φ 12 (W), računa potpuno isto kao<br />

kod zatvorenih sustava, tj. vrijede iste jednadžbe uz zamjenu mase m s protočnom masom ṁ .<br />

Isto vrijedi i za proračun promjene entropije, ∆ Ṡ12<br />

(J/K), radne tvari od usisnog stanja (1) do<br />

ispušnog stanja (2). Transportne sekvence procesa, usis i ispuh, pretpostavljene su bez<br />

gubitaka (∆S = 0).<br />

Teorijski gubitak snage uslijed nepovratnosti izmjene topline (samo pri politropi) računa se<br />

prema referentnom prirodnom toplinskom spremniku temperature T 0 = konst. pomoću<br />

relacije:<br />

∆P<br />

= T ∆ ̇<br />

0<br />

, W, (teorijski gubitak snage). (6.25)<br />

S s<br />

12<br />

=<br />

Otvoreni procesi specijalnih politropa<br />

Ostajući pri aproksimaciji procesa u otvorenim sustavima, kao kombinaciji ravnotežne<br />

promjene (politrope) idealnog plina s dvije transportne sekvence, usisom (ulazom) i ispuhom<br />

(izlazom), promotrit ćemo procese sa specijalnim politropama.<br />

Postoje i sličnosti i razlike između politropskih procesa u otvorenim i onih u zatvorenim<br />

sustavima. Kod otvorenih sustava uvodi se putem protočne mase ṁ (kg/s) utjecaj vremena i<br />

ponavljanje procesa, dok se kod zatvorenih sustava proces s konstantnom masom m (kg)<br />

odvija jednokratno, bez ponavljanja.<br />

U oba slučaja politropska promjena se odvija potpuno identično, jer se i kod otvorenih<br />

sustava ta sekvenca odvija pri zatvorenim ventilima, kao što je to slučaj kod zatvorenih<br />

sustava. Stoga je početno i konačno stanje radnog medija identično u oba slučaja. Sukladno<br />

usvojenim pretpostavkama kod otvorenih sustava ta stanja nisu pod utjecajem transportnih<br />

sekvenci, usisa i ispuha.<br />

Promjena energije radnog medija: unutarnje, kinetičke i potencijalne, tijekom politropskog<br />

procesa računa se po istom principu, tj. s istim jednadžbama u oba slučaja. Pod utjecajem<br />

protočne mase ṁ (kg/s), umjesto mase m (kg), nastaje samo razlika u dimenziji tih energija.<br />

Jednako tako, umjesto topline, Q 12 (J), sada govorimo o toplinskom toku, Φ 12 (J/s). Kako<br />

ćemo zanemariti izmjenu topline pri transportnim sekvencama usisa i ispuha to je grafički<br />

70


Mirko Tadić<br />

Termodinamika<br />

prikaz politropskog procesa u ravnini T-s identičan onom kod zatvorenih sustava. Ukratko,<br />

usis i ispuh se ne vide u dijagramu.<br />

Jedina razlika očituje se u promjeni grafičkog prikaza u ravnini p-V, budući da se transportne<br />

sekvence (usis i ispuh) kod otvorenih sustava odražavaju na ukupni mehanički efekt, tj. snagu.<br />

Zbog ponavljanja procesa poništava se mehanička interakcija prema okolišnjem zraku, koja<br />

kod jednokratnih procesa u zatvorenim sustavima dolazi do izražaja kao rad W 0 .<br />

Grafički prikazi otvorenih procesa u p-V i T-s dijagramu dani su u nastavku za posebne<br />

politropske procese.<br />

Izohora (V = konst.; n = ± ∞)<br />

p<br />

N/m 2<br />

n = ±∞<br />

ispuh 2<br />

2<br />

p 2 T 2<br />

ϑ 2<br />

2<br />

{ P } okret<br />

= −∫ V dp<br />

OS<br />

12<br />

1<br />

1<br />

p 1<br />

1<br />

T<br />

usis<br />

1<br />

p 1<br />

2<br />

ϑ 1<br />

q 12<br />

= ∫Tds<br />

1<br />

V, m 3 V1 = V /okret<br />

2<br />

s 1 s 2<br />

Izobara (p = konst.; n = 0)<br />

2<br />

12<br />

= n0<br />

12<br />

0∫<br />

1<br />

p2<br />

1<br />

{ P } = −n<br />

V dp = V̇<br />

( p )<br />

P<br />

okret<br />

−<br />

T<br />

K<br />

OS<br />

v 1 = v 2<br />

n = ±∞<br />

p 2<br />

( T )<br />

Φ = ̇ ̇<br />

v<br />

−<br />

12<br />

mq12<br />

= mc<br />

2<br />

T1<br />

Slika 6.5 Otvoreni izohorni proces u p-V i T-s dijagramu<br />

s, J/(kg K)<br />

p<br />

N/m 2<br />

p<br />

T 2<br />

1 2<br />

s 1<br />

T<br />

K<br />

v 2<br />

V 1 V 2 V, m 3 /okret<br />

s 2<br />

OS<br />

OS<br />

2<br />

usis n = 0<br />

ϑ 1<br />

ispuh<br />

ϑ 2 v 1<br />

T 1<br />

1<br />

q 12 > 0<br />

p 1 = p 2<br />

P<br />

12<br />

= 0<br />

Φ = ṁ<br />

q = mc ̇<br />

p<br />

( T − )<br />

12 12<br />

2<br />

T1<br />

s, J/(kg K)<br />

Slika 6.6 Otvoreni izobarni proces u p-V i T-s dijagramu<br />

71


Mirko Tadić<br />

Termodinamika<br />

Izentropa (S = konst.; n = κ)<br />

p<br />

T<br />

N/m 2<br />

K<br />

v<br />

usis<br />

1<br />

1<br />

p 1<br />

p 1 1<br />

MS<br />

dp<br />

ϑ 1<br />

T 1<br />

2<br />

n = κ<br />

2<br />

{ P } okret ∫ V<br />

12<br />

= −κ ( p)dp<br />

2<br />

1<br />

p 2<br />

ispuh ϑ 2<br />

T 2<br />

v 2<br />

p 2<br />

V 1<br />

V<br />

V 2<br />

κ<br />

P12<br />

= n<br />

̇<br />

0 12 okret 1<br />

−<br />

κ −1<br />

{ P } = mR( T T )<br />

V, m 3 /okret<br />

2<br />

s 1 = s 2<br />

q 12 = 0<br />

s, J/(kg K)<br />

Izoterma (T = konst.; n = 1)<br />

Slika 6.7 Otvoreni izentropski proces u p-V i T-s dijagramu<br />

p<br />

ϑ T<br />

N/m 2 1<br />

K<br />

usis<br />

p 1 1<br />

OS<br />

MS<br />

dp<br />

v 2<br />

p 1 p<br />

1<br />

2 2<br />

T 1 = T 2<br />

2<br />

{ P } okret ∫ V<br />

12<br />

= − ( p)dp<br />

1<br />

2<br />

p q 2 12 > 0<br />

ispuh<br />

ϑ 2<br />

n =1<br />

v 1<br />

s 1<br />

V 1 V<br />

V2<br />

p<br />

P12<br />

= n0{ P12<br />

} = mRTln ̇<br />

okret<br />

p<br />

1<br />

2<br />

V, m 3 /okret<br />

Φ<br />

= ̇ ̇<br />

s 2<br />

( s − )<br />

12<br />

mq12<br />

= mT<br />

2<br />

s1<br />

s, J/(kg K)<br />

Slika 6.8 Otvoreni izotermni proces u p-V i T-s dijagramu<br />

Maksimalni rad i eksergija<br />

Već smo upoznali način određivanja teorijskog gubitka na radu u zatvorenim, odnosno<br />

snazi u otvorenim sustavima, pomoću kojih se ocjenjuje termodinamička valjanost procesa u<br />

odnosu na nepovratnost izmjene topline.<br />

72


Mirko Tadić<br />

Termodinamika<br />

Preostaje nam da razmotrimo kriterij vrednovanja zadanog stanja radnog medija, polazeći od<br />

činjenice da sva moguća početna stanja nisu jednako vrijedna u odnosu na mogućnost<br />

dobijanja korisnog rada, odnosno snage.<br />

Očito je, da kriterij ne smije imati nikakve realne karakteristike koje bi utjecale na rezultat<br />

procjene. Stoga je logično da se kao osnova kriterija uzme idealni teorijski proces radnog<br />

medija u kojem nema nepovratne izmjene topline. Mehanička interakcija tijekom procesa<br />

smatra se reverzibilnom, tj. bez gubitaka, kao i u svim prethodnim razmatranjima. Tijekom<br />

procesa radni medij prolazi kroz stanja u unutarnjoj ravnoteži.<br />

Radni sustav uključuje samo nužne sudionike: sam radni medij, neopisanog mehaničkog<br />

sudionika čije su osobine idealizirane (bez trenja), te postojeći prirodni toplinski spremnik<br />

(okolišnji zrak) temperature T 0 = konst., koji je ujedno i mehanički sudionik tlaka p 0 = konst..<br />

Sustav je potpuno izoliran od drugih utjecaja. Pretpostavit ćemo da radni medij zadanog<br />

stanja (T 1 , p 1 , V 1 ) nije u ravnoteži s okolišem, jer u protivnom slučaju ne bi bio moguć<br />

nastanak procesa niti teorijski. Neravoteža s okolišem pruža mogućnost provođenja procesa i<br />

dobivanja korisnog mehaničkog rada, odnosno snage. Proces traje do uspostave ravnoteže s<br />

prirodnim spremnikom (okolišnjim zrakom), kada radni medij postiže temperaturu T 0 i tlak<br />

p 0 .<br />

Stanje plina u potpunosti je zadano podacima temperature T 1 , tlaka p 1 i volumena V 1 , odnosno<br />

V̇ 1, na osnovu čega se iz jednadžbe stanja može izračunati masa m, odnosno protočna masa<br />

ṁ , plina:<br />

p1V<br />

1<br />

m = , kg,<br />

RT1<br />

(zatvoreni sustav), (6.26)<br />

p1V̇<br />

1<br />

ṁ = , kg/s,<br />

RT<br />

(otvoreni sustav). (6.27)<br />

1<br />

Da bi se izbjegli gubici na radu, odnosno snazi, koji su spomenuti na početku, tijekom procesa<br />

ne smije doći do nepovratne izmjene topline. Teorijski zamišljena povratna izmjena topline<br />

mogla bi biti samo pri istoj temperaturi radnog medija i toplinskog spremnika. Ta je<br />

temperatura zadana zadavanjem stanja prirodnog toplinskog spremnika, okolišnjeg zraka.<br />

Tamperatura T 0 je upravo ta povratna izoterma po kojoj smije teći proces plina. Naravno,<br />

proces uravnotežavanja može koristiti i bilo koju politropu bez izmjene topline, tj. izentropu,<br />

S = konst., dS = 0. Promjena entropije takvog izoliranog sustava jednaka je nuli, ∆S s = 0 ili<br />

Ṡ = 0 , kao i teorijski gubitak na radu, ∆W = 0, odnosno snazi, ∆P = 0.<br />

∆ s<br />

S postavljenim uvjetom reverzibilnosti izmjene topline uklonjena je poslijednja realna<br />

karakteristika, pa takav proces ima samo teorijski smisao. Pri takvom procesu teorijski se<br />

dobiva najveći mogući mehanički rad, odnosno najveća moguća snaga (eksergija), s obzirom<br />

na zadana stanja radnog medija i prirodnog spremnika.<br />

Maksimalni rad – zatvoreni sustav<br />

Teorijski idealan proces u zatvorenom sustavu smije biti samo kombinacija bilo koje<br />

izentrope, pri kojoj nema izmjene topline, dok se izmjena topline smije odvijati samo tijekom<br />

povratne izoterme T 0 , koja je zadana temperaturom prirodnog spremnika, okolišnjim zrakom.<br />

U posebnim slučajevima bit će dovoljna samo jedna od tih promjena.<br />

73


Mirko Tadić<br />

Termodinamika<br />

Na slici 6.10 prikazan je jedan slučaj takvog idealnog procesa za zadano stanje (1) radnog<br />

medija i zadano stanje (0) prirodnog spremnika, okolišnjeg zraka temperature T 0 i tlaka p 0 .<br />

p<br />

p 1<br />

1<br />

p 1<br />

1<br />

T 1<br />

ϑ 1<br />

v 0 p 0<br />

2<br />

0<br />

T 0<br />

2<br />

W max<br />

0<br />

p q 0<br />

ϑ o<br />

20 > 0<br />

n = κ<br />

p 2<br />

W 0<br />

s 1 = s 2<br />

v<br />

T<br />

1<br />

V V, m 3 1<br />

V0 s 0<br />

s, J/(kg K)<br />

Slika 6.9 Idealni proces zatvorenog sustava u ravnini p-V i T-s<br />

Dogod je temperatura radnog medija različita od temperature T 0 smije se koristiti samo<br />

izentropski proces bez izmjene topline. Tek kada plin postigne tu temperaturu (stanje 2 u<br />

dijagramima) može se zamisliti teorijska izmjena topline, pri istim temperaturama plina i<br />

toplinskog spremnika, okolišnjeg zraka. Proces se po toj povratnoj izotermi odvija sve dok se<br />

ne uspostavi potpuna ravnoteža plina i okoliša, tj. plin poprimi stanje s temperaturom T 0 i<br />

tlakom p 0 (stanje 0 u dijagramima).<br />

Položaj stanja (2) u ravnini p-V može se procijeniti ako se izračuna pripadni tlak p 2 :<br />

κ<br />

⎛ T<br />

κ−1<br />

0<br />

⎞<br />

p<br />

2<br />

= p1<br />

⎜<br />

⎟ ,<br />

T1<br />

(za sjecište izentrope 1-2 i izoterme 2-0). (6.28)<br />

⎝<br />

⎠<br />

Teorijski rad plina je, po definiciji, površina ispod linija procesa (1-2-0). Kako se proces (1-2)<br />

odvija po zakonu izentrope, a proces (2-0) po zakonu izoterme, to se mehanički rad može<br />

izračunati samo odvojeno za svaku sekvencu. Stoga je<br />

W , J. (mehanički rad plina pri idealnom procesu). (6.29)<br />

1− 2−0<br />

= W1<br />

−2<br />

+ W2−0<br />

No, W 1-2-0 se može izraziti i na drugi način, kako će se pokazati u nastavku.<br />

U promatranom primjeru ovaj je rad pozitivan tj. predaje se prisutnim mehaničkim<br />

sudionicima. Dio tog rada troši se na potiskivanje okolišnjeg zraka u iznosu W 0 , a ostatak se<br />

može iskoristiti. Zbog idelnog procesa taj koristana rad je ujedno i najveći mogući rad, W max ,<br />

koji se teorijski može dobiti od plina zadanog stanja – u odnosu na zadano stanje prirodnog<br />

spremnika, okolišnjeg zraka (T 0 , p 0 ).<br />

Odnos teorijskog rada plina W 1-2-0 , rada okoline W 0 i maksimalnog korisnog rada W max , uvijek<br />

je određen jednadžbom oblika:<br />

W max<br />

= W 1 − 2−0<br />

+ W0<br />

, (6.30)<br />

Rad prema okolini može se izraziti, pomoću definicije mehaničkog rada, na slijedeći način:<br />

74


Mirko Tadić<br />

Termodinamika<br />

0,ok<br />

W0 p0<br />

dV0<br />

= p0<br />

0,ok<br />

−<br />

1,ok<br />

1,ok<br />

( V V )<br />

= ∫ , J, (rad prema okolini). (6.31)<br />

Razliku volumena okolišnjeg zraka moramo zamijeniti s istom takvom razlikom volumena<br />

radnog medija, ali suprotnog smisla:<br />

V<br />

( V )<br />

0, ok<br />

V1,<br />

ok<br />

= −<br />

0<br />

−V1<br />

− , (6.32)<br />

pa umjesto jednadžbe (6.31) možemo pisati:<br />

W<br />

( V − )<br />

0<br />

p0<br />

1<br />

V0<br />

= , J, (rad prema okolini). (6.33)<br />

Za određivanje teorijskog rada plina W 1-2-0 , tijekom procesa (1-2-0), iskoristit ćemo jednadžbu<br />

I. Zakona:<br />

Q − W = ∆U<br />

, J, (6.34)<br />

1−2−0<br />

1−2−0<br />

1−2−0<br />

koju možemo uskladiti s procesom u dvije sekvence, (1-2) i (2-0):<br />

Q + Q −W<br />

= ∆U<br />

+ ∆U<br />

, J. (6.35)<br />

1−2<br />

2−0<br />

1−2−0<br />

1−2<br />

2−0<br />

Na sekvenci izentrope (1-2) nema izmjene topline, Q 1-2 = 0, dok se toplina na sekvenci (2-0)<br />

odvija pri konstantnoj temperaturi, pa se može izračunati primjenom II. Zakona:<br />

( S − S ) = T ( S − )<br />

Q<br />

−<br />

, J, (jer vrijedi: S 2 = S 1 ). (6.36)<br />

2 0<br />

= T0∆S<br />

2−0<br />

= T0<br />

0 2 0 0<br />

S1<br />

Promjena unutarnje energije radnog medija tijekom procesa (1-2-0) iznosi:<br />

( U<br />

2<br />

−U1) + ( U<br />

0<br />

−U<br />

2<br />

) = U<br />

0<br />

−<br />

1<br />

∆U<br />

−<br />

, J, (jer vrijedi: U 2 = U 0 ). (6.37)<br />

1 2<br />

+ ∆U<br />

2−0<br />

=<br />

U<br />

Rezultate jednadžbi (6.35) i (6.37) možemo uvrstiti u jednadžbu (6.35) i preoblikovati u izraz<br />

kojim se izražava W 1-2-0 :<br />

( S − )<br />

W<br />

−<br />

, J. (6.38)<br />

1 2−0<br />

= U1<br />

−U<br />

0<br />

− T0<br />

1<br />

S0<br />

Konačno, uvrštavanjem jednadžbi (6.38) i (6.33) u jednadžbu (6.29) dobivamo izraz za<br />

izračunavanje maksimalnog rada za proces između zadanih stanja 1-0:<br />

( S − S ) + p ( V − )<br />

W max, 1 0<br />

= U1<br />

−U<br />

0<br />

− T0<br />

1 0 0 1<br />

V0<br />

−<br />

, J, (maksimalni rad). (6.39)<br />

Sve veličine na desnoj strani jednadžbe se odnose na radni medij, a pojedini članovi se<br />

računaju prema slijedećim relacijama:<br />

( T )<br />

U −<br />

v<br />

− , J, (6.40)<br />

1<br />

U<br />

0<br />

= mc<br />

1<br />

T0<br />

75


Mirko Tadić<br />

Termodinamika<br />

⎛ T p ⎞<br />

T<br />

p<br />

, J, (6.41)<br />

0<br />

1<br />

1<br />

( S − ) =<br />

⎜ − ⎟ 1<br />

S0<br />

mT0<br />

c ln R ln<br />

⎝ T0<br />

p0<br />

⎠<br />

( V )<br />

p − , J, pomoću:<br />

0 1<br />

V0<br />

mRT<br />

1<br />

V<br />

1<br />

= i<br />

p1<br />

mRT<br />

0<br />

V<br />

0<br />

= . (6.42)<br />

p0<br />

U dijagramu p-V uvijek se W max,1-0 može prikazati površinom, koja je omeđena s linijama<br />

idealnog procesa (1-2-0) i s dvije pomoćne linije (nisu linije procesa): izohorom V 1 i izobarom<br />

p 0 .<br />

Radi ilustracije prikazan je na slici 6.10 proces dobijanja maksimalnog rada u slučaju kada je<br />

zadano stanje plina (T 1 , p 1 , V 1 ) ispod stanja prirodnog spremnika, okolišnjeg zraka (T 0 , p 0 ).<br />

p<br />

p 0<br />

p 0<br />

2<br />

q v 20 < 0 1 p 1<br />

p 2<br />

v 0<br />

0<br />

T<br />

0<br />

T 0<br />

W max<br />

ϑ o<br />

2<br />

ϑ 1<br />

V V, m 3 0<br />

s 0<br />

n = κ<br />

p 1<br />

1<br />

T 1<br />

1<br />

V1 s 1 = s 2<br />

s, J/(kg K)<br />

Slika 6.10 Idealni proces s plinom stanja ispod okolišnjeg<br />

Za teorijski proces (1-2-0) potrebno je utrošiti rad W 1-2-0 (površina ispod procesa – bez oznake<br />

u dijagramu). Pri tome okolina pomaže dajući rad W 0 = p 0 (V 1 – V 0 ). Razlika je koristan<br />

(pozitivan) rad W max .<br />

Minimalni rad<br />

U slučaju da je početno stanje radnog medija u ravnoteži s prirodnim spremnikom<br />

(stanje s oznakom 0 u dijagramima) može se po istom principu odrediti teorijski minimalan<br />

rad, W min , da bi se radni medij doveo u stanje (1). Idealni proces je identičan, ali se odvija u<br />

suprotnom smjeru, pa je konačni rezultat za W min,0-1 jednak rezultatu za W max,1-0 , samo<br />

suprotnog, negativnog predznaka.<br />

( S − S ) + p ( V − )<br />

W min, 0 1<br />

= U<br />

0<br />

−U1<br />

− T0<br />

0 1 0 0<br />

V1<br />

−<br />

, J, (minimalni rad). (6.43)<br />

Eksergija – otvoreni sustav<br />

Ako raspolažemo sa stalnim dotokom ṁ (kg/s) radnog medija zadanog stanja, (T 1 , p 1 ),<br />

može se primjenom istog principa idealnog procesa odrediti eksergija (= maksimalna snaga),<br />

koja pripada tom stanju u odnosu na referentno stanje prirodnog spremnika, okolišnjeg zraka,<br />

76


Mirko Tadić<br />

Termodinamika<br />

(T 0 , p 0 ). Zbog ponavljanja, idealnom procesu pripadaju transportne sekvence usisa i ispuha.<br />

Primjer eksergijskog procesa prikazan je na slici 6.10.<br />

p<br />

p<br />

T<br />

1<br />

usis 1<br />

p 1<br />

1<br />

T 1<br />

ϑ 1<br />

v 0 p 0<br />

{ E} okret<br />

2<br />

0<br />

T 0<br />

2<br />

p 2<br />

ϑ 0<br />

0<br />

p q 0 20 > 0<br />

ispuh<br />

n = κ<br />

V, m 3 /okret<br />

s 1 = s 2<br />

V 1<br />

v 1<br />

V0 s 0<br />

s, J/(kg K)<br />

Slika 6.11 Eksergijski proces (1-0) u p-V i T-s ravnini<br />

Budući da se radi o otvorenom sustavu grafički prikaz eksergije (po okretu) u ravnini p-V<br />

predstavlja površinu, dobivenu projekcijom idealnog procesa (1-2-0) na ordinatnu os p.<br />

Prema I. Zakonu za otvorene sustave je:<br />

Φ1 −2−0<br />

− P1<br />

−2−0<br />

= ∆Ḣ 1−2−0<br />

, W, (6.44)<br />

odnosno<br />

Φ + Φ − P = ∆Ḣ + ∆Ḣ<br />

, W. (6.45)<br />

1−2<br />

2−0<br />

1−2−0<br />

1−2<br />

2−0<br />

Za pojedine sekvence idealnog procesa vrijedi:<br />

0 Φ = mT ̇ s − s = mT ̇ s − s = T S − , W, (6.46)<br />

Φ , ( ) ( ) ( )<br />

1 −2<br />

=<br />

2−0<br />

0 0 2 0 0 1 0 0<br />

S1<br />

∆Ḣ<br />

1 2<br />

= Ḣ<br />

2<br />

− Ḣ<br />

− 1<br />

, ∆H<br />

̇ 2 − 0<br />

= Ḣ<br />

0<br />

− Ḣ<br />

2<br />

= 0 , W, ( Ḣ 2<br />

= Ḣ<br />

). (6.47)<br />

0<br />

Snaga idealnog procesa je eksergija stanja (1) u odnosu na referentni prirodni spremnik stanja<br />

(0):<br />

P ≡ E , W, (6.48)<br />

1−2−0<br />

1−0<br />

Uvrštavanjem jednadžbi (6.46), (6.47) i (6.48) u jednadžbu (6.44) dobiva se:<br />

E<br />

( S − )<br />

̇ ̇<br />

, W, (eksergija). (6.49)<br />

1−0<br />

= H1<br />

− H<br />

0<br />

− T0<br />

1<br />

S0<br />

Računanje pojedinih članova vrši se pomoću jednadžbi:<br />

( T )<br />

Ḣ ̇ ̇<br />

p<br />

− , (6.50)<br />

T<br />

1<br />

− H<br />

0<br />

= mc<br />

1<br />

T0<br />

0<br />

⎛ T p ⎞<br />

̇ , (6.51)<br />

1<br />

1<br />

( S − ) =<br />

⎜ − ⎟ 1<br />

S0<br />

mT0<br />

c<br />

p<br />

ln R ln<br />

⎝ T0<br />

p0<br />

⎠<br />

koje zahtijevaju poznavanje samo svojstava početnog (1) i konačnog (0) stanja.<br />

77

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

Saved successfully!

Ooh no, something went wrong!