autor:
ŠKORPÍK, Jiří – LinkedIn.com/in/jiri-skorpik
datum vydání:
Září 2024
název:
Termodynamika turbokompresorů
sborník:
provenience: Brno (Česká republika)
email: skorpik.jiri@email.cz
Copyright©Jiří Škorpík, 2006-2024 |
Komprese v turbokompresoruKompresní poměrCharakteristickým rysem komprese v turbokompresoru je plynulost děje transformace práce na tlakovou a vnitřní energii pracovního plynu. Základním požadavkem je zvýšení tlaku, respektive dosažení požadovaného kompresního poměru (viz Rovnice 1), při minimálním zvýšení teploty, která významně narůstá zejména pokud není komprese chlazená. 1: Kompresní poměr kompresoru p [Pa] tlak; ε [1] kompresní poměr. Index i označuje stav na vstupu do turbokompresoru, index e označuje stav na výstupu z turbokompresoru.Při vyšetřování průběhu komprese je třeba odlišovat jednostupňovou a vícestupňovou kompresi. Obě komprese vyžadují trochu jiný přístup k pochopení příčin ztrát a přístup ke stanovení opatření k jejich snížení. Adiabatická kompreseVýpočtový model adiabatické komprese se používá v případech, při kterých se neočekává významný vliv sdílení tepla s okolím turbokompresoru. Izoentropická kompreseVýstupní rychlostIdeální vnitřní práce adiabatické komprese odpovídá práci izoentropické komprese. Při ideální kompresi se obvykle očekává stejná výstupní rychlost jako při reálné kompresi. To znamená, že reálný stroj musí mít o něco větší průtočné průřezy než stroj ideální, protože ztrátové teplo měrný objem pracovního plynu zvětšuje. Přídavné ztrátyZtrátové teploVnitřní práceh-s diagramT-s diagramZtrátyCharakteristickým rysem reálné komprese jsou tzv. přídavné ztráty Δ. Tyto ztráty jsou ekvivalentní práci, kterou je nutné přivést navíc kvůli zvýšení měrného objemu pracovního plynu způsobeného ztrátovým teplem. Na Obrázku 2 je příklad komprese v kompresoru nebo v jeho stupni v h-s a T-s diagramu. Zatímco v h-s diagramu lze odlišit pouze ztráty jako celek, tak v T-s diagramu lze odlišit jednotlivé typy ztrát. |
2: Vnitřní práce turbokompresoru při adibatické kompresi v h-s a T-s diagramu Lq [J·kg-1] ztrátové teplo ([Škorpík, 2024]); ΔeK [J·kg-1] rozdíl kinetické energie mezi vstupem a výstupem (obvykle nevýznamný rozdíl); h [J·kg-1] entalpie; Δh [J·kg-1] rozdíl entalpii; Δhis [J·kg-1] rozdíl entalpii při izoentropické kompresi; Lw [J·kg-1] vnitřní ztráty v kompresoru (práce přívedená do stupně navíc oproti izoentropické kompresi); s [J·kg-1·K-1] entropie; T [K] absolutní teplota; V [m·s-1] rychlost; v [m3·kg-1] měrný objem; wi [J·kg-1] vnitřní práce; wis [J·kg-1] vnitřní práce při izoentropické kompresi; Δ [J·kg-1] přídavné ztráty. Index is označuje stavy při izoentropické kompresi, index s celkový stav. T-s diagram je sestrojen při nevýznamném rozdílu kinetických energií. Rovnice jsou odvozeny v Příloze 5.Vícestupňová kompreseVnitřní účinnostVnitřní účinnost stupněSoučinitel přídavných ztrátTypickou vlastností vnitřní účinnosti vícestupňové komprese ηi je, že je menší než průměrná vnitřní účinnost jednotlivých stupňů ηj, viz Obrázek 3, na kterém je příklad h-s diagramu komprese turbokompresoru se Z stupni. Příčinou jsou přídavné ztráty. Je tedy zřejmé, že vnitřní ztráty ve stupni kompresoru zhoršují účinnost v následujícím stupni. Poměr průměrné hodnoty vnitřní účinnosti jednotlivých stupňů η,j ku vnitřní účinnosti měřené mezi prvním a posledním stupněm ηi se nazývá součinitel přídavných ztrát 1+f, viz Úloha 1. 3: Vícestupňová adiabatická komprese v kompresoru |
1+f [1] součinitel přídavných ztrát (preheat factor); Z [-] počet stupňů; Δj [J·kg-1] přídavné ztráty jednoho stupně; ηi [1] vnitřní účinnost komprese mezi body 1-Z. Index j označuje j-tý stupeň. Rovnice jsou odvozeny pro předpoklad, že všechny stupně zpracovávají stejný entalpický spád a komprese je adiabatická. Pro přehlednost není v obrázku zakreslena kinetická energie absolutní rychlosti. Rovnice jsou odvozeny v Příloze 6.
Polytropická kompreseVratná polytropická kompreseOhřevV některých případech kompresi ovlivňuje sdílení tepla s okolím kompresoru. Například při záměrném chlazení kompresoru, nebo, když je komprimován kryogenní plyn, který je ohříván okolním prostředí. V takových případech je komprese podobná polytropické kompresi – porovnávací ideální kompresí je v tomto případě vratná polytropická komprese. Polytropickou kompresi popisují Rovnice 4. Tyto rovnice lze odvodit z rovnice prvního zákona termodynamiky. 4: Vnitřní práce kompresoru pro případ q>0 q [J·kg-1] teplo sdílené s okolím; Lu [J·kg-1] disipované množství energie, respektive úhrn různých druhů energií transformovaných na vnitřní energii plynu v průběhu komprese. Index pol označuje vratnou polytropickou kompresi. Na obrázku je případ q>0 (přívod tepla – jestliže he, pol-he,is>0, pak se jedná o součet přivedeného tepla a přídavných ztrát, které vznikl v důsledku přívodu tepla). T-s diagram je sestrojen při nevýznamném rozdílu kinetických energií.ChlazeníPro úplnost jsou na Obrázku 5 uvedeny případy s odvodem tepla při kompresi. |
5: Vnitřní práce kompresoru pro případ chlazené komprese (a) případ pro Lw<-q; (b) případ kdy Te=Ti (zdánlivě izotermická komprese – teplota chladícího média musí být v tomto případě nižší než teplota pracovního plynu na vstupu do kompresoru Ti). Jestliže he,pol-he,is<0, pak se jedná o součet odvedeného tepla a úspory práce v důsledku chlazení komprese. T-s diagram je sestrojen při nevýznamném rozdílu kinetických energií.Energetická bilanceVnitřní účinnostPři vytváření energetických bilancí polytropické komprese je nutné definovat práci při vratné polytropické kompresi wpol. Často se za wpol používá práce vratná izotermická komprese, zejména je-li komprese chlazená, viz Úloha 2. Zejména při kompresi s přívodem tepla se práce komprese porovnává s prací při izoentropické kompresi wis. Při uvádění účinností je nutné uvést jaký děj byl vybrán jako porovnávací, aby hodnota vnitřní účinnosti měla vypovídající hodnotu. Provedení chlazení turbokompresorůChlazení komprese je nejefektivnější způsob snížení vnitřní práce kompresoru, přičemž existuje několik možností jak toho dosáhnout. Komprimovaný plyn v průběhu komprese lze kontinuálně chladit dvěma způsoby a to povrchově nebo vstřikování chladící kapaliny do komprimovaného plynu. Chlazení lze ale provést i přetržitě po stupních pomocí tzv. mezichlazení. Každé chlazení ale generuje nový druh ztrát, takže efektivní chlazení lze provést jen za určitých podmínek. Povrchové chlazeníRadiální stupeňPovrchové neboli plášťové chlazení lze provést u dvouplášťových kompresorů, přičemž mezi plášti proudí chladící kapalina, která ochlazuje pracovní plyn uvnitř, viz Obrázek 6. Kompresory s povrchovým chlazením jsou složité a drahé – jsou nutné rozvodové kanálky a v dělící rovině hrozí únik chladící kapaliny do komprimovaného plynu a obráceně. |
6: Turbokompresor s jedenácti radiálními stupni a s povrchovým chlazením Vyrobil Demag.VzduchPovrchové chlazení je málo účinný způsob a tak se používá u kompresorů s malým stlačením v jednom stupni, u jednostupňových turbokompresorů a dmychadel je to jediný způsob chlazení komprese. Na druhou stranu, jeho nízká účinnost umožňuje jeho použití pro chlazení i vlhkého vzduchu obsahující prach, s tím, že teplota povrchu kompresoru neklesne pod rosný bod vzduchu. Chlazení vstřikováním chladící kapalinyVstřikováním chladící kapaliny do proudu pracovního plynu se plyn ochlazuje v důsledku odpařování chladící kapaliny. Rychlost odpaření a tedy i ochlazení závisí, mimo jiné, na vzájemné teplosměnné ploše chladící kapaliny a pracovního plynu, proto jsou vstřikovací trysky (Obrázek 7) konstruovány tak, aby měly co největší rozptyl. I tak k odpaření je nutný určitý úsek a proto jsou pro chlazení vstřikováním vhodnější radiální stupně (chladící kapalina se vstřikuje v místě za statorovými lopatkami směrem do vratného kanálu k dalšímu stupni) v případě axiálních stupňů by se musela v místě vstřiku zvětšit mezera mezi stupni. |
7: Princip vstřikovací trysky chladící kapaliny kompresoru VzduchČpavekMnožství chladící kapaliny závisí na tlaku, požadované teplotě a složení výsledné směsi po vychlazení. Například, pokud je komprimovaným plynem vzduch, tak lze vstříknout pouze takové množství chladící vody, aby po odpaření byla relativní vlhkost vzduchu menší než 100 %, v opačném případě zůstanou ve vzduchu kapičky vody. Při kompresi čpavku se používá kapalný čpavek, při kompresi nitrózní plynů se používá slabý roztok kyseliny dusičné apod. VlhkostNevýhodou tohoto způsobu chlazení je, že na výtlaku kompresoru plyn obsahuje jistou vlhkost. To znamená, že využití takového plynu je omezeno pouze na aplikace, kde vlhkost v plynu není překážkou k jeho využití. Jedná se zejména o takové procesy, při kterých může dojít ke kondenzaci par v plynu obsažených. Například při použití v pneumatických pohonech apod. Vnější chlazení (Mezichlazení)Vícestupňový turbokompresorSdílení tepla s okolím nemusím být rozloženo rovnoměrně po celé kompresi, například při použití vnějšího chlazení, respektive mezichlazení. To spočívá ve vyvedení komprimovaného plynu za vybranými stupni kompresoru mimo kompresor do rekuperačního výměníku tepla (nejčastěji tvořený žebrovanými trubkami), kde se pomocí chladící kapaliny (obvykle voda) plyn ochladí. Například v případě Obrázku 8, kde je provedeno mezichlazení u třístupňového turbokompresoru, lze rozdělit celou kompresi na 3 samostatné komprese a vnitřní práci kompresoru vypočítat z rozdílů entalpií a odvedeného tepla, viz Úloha 3. Výhodou vhodně navrženého mezichlazení je i to, že práce jednotlivých stupňů a jejich pracovní podmínky jsou si podobné (rychlostní trojúhelníky, teploty apod.), je ale nutné počítat s tím, že mezi stupni se zmenšuje měrný objem plynu, proto první stupeň po mezichlazení bude mít menší vstupní průřezy než výstupní průřezy předchozího stupně. |
8: Princip komprese s vnějším chlazením C-mezichladiče plynu; e*-konečný stav pracovního plynu na výstupu z kompresoru v případě komprese bez chlazení. pC1, pC2 [Pa] tlak před vstupem do mezichladičů; w*i [J·kg-1] vnitřní práce kompresoru pro případ komprese bez chlazení.Radiální stupeňNa Obrázku 9 je příklad konstrukce turbokompresoru se sedmi radiálními stupni a se dvěma mezichladiči, přičemž první mezichladič je umístěný za druhým stupněm, druhý za čtvrtým stupněm. Všimněte si, že před demontáží vrchní skříně kompresoru je nutné demontovat i chladiče (viz příčný řez). Pro jednodušší demontáž kompresoru existují i jiná konstrukční uspořádání – napojení vstupu i výstupu mezichladiče přes spodní skříň. V případě vícetělesových turbokompresorů se může mezichlazení instalovat v propoji mezi jednotlivými tělesy. 9: Turbokompresor se sedmi radiálními stupni a se dvěma mezichladiči Vyrobil Escher Wyss. |
Tento způsob chlazení doprovází větší konstrukční i investiční náročnost (mimo kompresor je nutné pořídit zařízení pro chlazení), a proto se obvykle provádí jen od určité velikosti turbokompresoru nebo k tomu musí být další než ekonomický důvod, například bezpečnost (u hořlavých plynů), stabilita plynu (molekuly při vyšší teplotě mohou disociovat, atd.). Efektivita chlazeníVnitřní práce kompresoru vs. Tlaková ztrátaKomprese vzduchuKomprese metanuKomprese páryKomprese heliaChlazení nemusí vždy znamenat výrazné snížení příkonu. Jakékoliv chlazení ovlivňuje termodynamiku komprese (chlazení zvyšuje tlakovu ztrátu, třetím tekutiny o teplosměnné plochy a tvořením vírů při vstřikování chladící kapaliny apod.), takže vždy existuje hranice efektivity chlazení. Jak plyne z porovnání práce izoentropické komprese s prací izotermické komprese, která odpovídá kompresi při dokonalém chlazením v grafu na Obrázku 10. Z grafu, respektive rovnice je například zřejmé, že pokud by navýšení práce na kompresi vzduchu kvůli tlakovým ztrátám bylo 10 %, tak by chlazení mělo pozitivní význam až při kompresních poměrech 2, při kompresi metanu dokonce až při 5 % apod. Při reálných kompresí je úspora práce mnohem menší, takže se vyplácí chlazení od vyšších kompresních poměrů než ukazuje graf. 10: Rozdíl mezi prací izoentropické a izotermické komprese wit [J·kg-1] práce izotermické komprese; Δwi [%] maximální teoretická úspora práce díky chlazení, Δwi=(wis·w-1it-1)100; κ [1] poměr tepelných kapacit pracovního plynu (κ=1,13 například CH4, κ=1,22 například C2H4, κ=1,33 například pára H2O, κ=1,4 například vzduch, κ=1,67 například He). Odvození rovnice je uvedeno v Příloze 7. |
Termodynamický návrh stupně turbokompresoruJapikse,1997Základními typy kompresorových stupňů jsou přetlakový axiální stupeň a radiální. Ojediněle se lze setkat s diagonálním provedením jednostupňového kompresoru. Většina stacionárních aplikací spadá do oblasti radiálních stupňů, až při požadavku na vyšší průtok cca od 15 m3·s-1 se kompresor konstruuje jako axiální, protože od takových průtoků už je účinnost axiálních turbokompresorů vyšší než radiálních. Na druhou stranu u radiálních stupňů se lépe realizuje mezichlazení, což více než kompenzuje nižší účinnost jednotlivých stupňů. Doporučené hodnoty podobnostních součinitelů pro návrh jednotlivých stupňů jsou uvedeny v [Japikse, 1997, s. 1-3]. h-s diagramyVentilační ztráta rotoruNa Obrázku 11 je h-s diagram axiálního stupně kompresoru na vyšetřovaném poloměru (Obrázku 11a) a radiálního stupně kompresoru (Obrázku 11b). Na každém vyšetřovaném poloměru může být vliv sdíleného tepla qE na Eulerovu práci různý – záleží na typu chlazení a teplu z ventilační ztráty rotoru. Modře je vyznačena energetická bilance celého stupně. 11: h-s diagram stupně kompresoru na poloměru r (a) axiální stupeň; (b) radiální stupeň. Lh [J·kg-1] profilové ztráty; qE [J·kg-1] sdílené teplo s okolím vyšetřované proudnice; ∑ L [J·kg-1] vnitřní ztráty stupně, součet všech ztrát ve stupni. Index 1 označuje stav před rotorovou řadou lopatek.Stupeň reakceStupně reakce R radiálních stupňů (centrifugální) jsou větší než nula, protože při nulovém stupni reakce by došlo, vlivem odstředivých sil, ke zvýšení relativní rychlosti W2, což je dobře patrné z h-s diagramu na Obrázku 11b. |
Eulerova práceOdtržení prouduZtráta zpětným prouděnímNa Obrázku 12 je očekávaný průběh Eulerovy práce wE,ref axiálního stupně kompresoru a Eulerova práce wE,is při proudění beze ztrát. Rozdíl mezi těmito pracemi jsou profilové ztráty na daném poloměru lopatky. Profilové ztráty jsou u okrajů lopatek největší a proto v těchto místech nelze dosáhnout potřebného zvýšení tlaku a naopak lze zde očekávat odtrhávání proudu a dokonce i ztrátu zpětným prouděním. Rozdíly v Eulerových prací mezi jádrem proudu a okraji lopatek jsou ještě větší v případě přímých lopatek axiálních stupňů, a proto se příliš u kompresorů nepoužívají. 12: Porovnání Eulerových prácí axiálního stupně kompresoru při reálném a izoentropickém proudění r [m] poloměr stupně; wE,is [J·kg-1] průběh Eulerovy práce při proudění beze ztrát; wE,ref [J·kg-1] navrhovaný lineární (konstantní) průběh Eulerovy práce částečně respektující průměrné ztráty stupně; Lw,m [J·kg-1] průměrné profilové ztráty stupně. Index h označuje patní poloměr lopatek, index t poloměr u špic lopatek.Lopatky turbokompresorůAxiální stupeňKuželový stupeňAxiální nebo kuželové stupně kompresorů obsahují zkroucené lopatky, ale existuje i nemálo případu s přímými lopatkami. h-s diagram pro radiální stupeň lze využít i při konstrukci kuželového stupně kompresoru, ve kterém kvůli snižujícímu se měrnému objemu klesá poloměr stupně, respektive délka lopatek, viz Obrázek 13 (výpočet kuželového stupně je uveden v článku Vnitřní ztráty lopatkových strojů a jejich vliv na návrh lopatkového stroje). Kvůli tenkým lopatkám nelze u kompresorových stupňů použít bandáže. |
13: Rychlostní trojúhelník axiálního kompresorového stupně R-rotorová řada lopatek; S-statorová řada lopatek. l [m] délka lopatek; rm [m] střední poloměr; U [m·s-1] obvodová rychlost; V [m·s-1] absolutní rychlost; W [m·s-1] relativní rychlost; α [°] úhel absolutní rychlosti; β [°] úhel relativní rychlosti. Rychlostní trojúhelník je nakreslen pro střední poloměr a stupeň reakce 0,5.Změna relativní vlhkosti vzduchu při kompresiRelativní vlhkostKondenzacePři kompresi vlhkého vzduchu se zvyšuje tlak plynů i tlak páry obsažený ve vzduchu. Při adiabatické kompresi bude obsah páry na konci komprese vždy v přehřátém stavu, a to i v případě komprese sytého vzduchu. Znamená to, že relativní vlhkost na konci komprese bude vždy menší než na počátku, a proto ke kondenzaci páry ve vzduchu nemůže docházet. Nicméně při vyšším tlaku vzroste i teplota kondenzace vlhkosti ve vzduchu z původní absolutní teploty Ti,C na teplotu Te,C, viz Obrázek 14. To znamená, že pokud klesne teplota komprimovaného vzduchu pod teplotu Te,C, začne vlhkost ve vzduchu kondenzovat. 14: T-s diagram komprese páry ve vzduchu p [Pa] parciální tlak vodní páry ve vzduchu; TC [K] absolutní teplota kondenzace páry ve vzduchu při tlaku na počátku komprese (index i) a konci komprese (index e); x [1] suchost páry. Na obrázku je případ izoentropické komprese. |
KondenzaceKe kondenzaci páry v komprimovaném vzduchu může dojít například v potrubí během jeho distribuce ke spotřebičům nebo při jeho ochlazování v mezichladičích kompresoru, či zásobnících komprimovaného vzduchu apod. V těchto případech se obvykle počítá s tím, že komprimovaný vlhký vzduch bude ochlazen na teplotu okolí, tj. teplotu na sání kompresoru Ti. Úkolem konstruktéra či projektanta tedy je stanovit, jestli při této teplotě dojde k vyloučení kondenzátu a v jakém množství podle Rovnice 15. Tato rovnice byla odvozena za předpokladu, že se vlhký vzduch vychladí na teplotu na sání, jestliže bude výsledná teplota při chlazení menší bude množství vyloučeného kondenzátu větší. 15: Množství vyloučeného kondenzátu z komprimovaného a ochlazeného vlhkého vzduchu mC [kg] množství vyloučeného kondenzátu z komprimovaného a ochlazeného vlhkého vzduchu zpět na teplotu ti (záporná hodnota znamená, že relativní vlhkost vzduchu na konci komprese a po vychlazení φe bude menší než 1, a proto nebude docházet ke kondenzaci); Vi [m3] objem zkomprimovaného vzduchu měřený na sání; v''i [m3·kg-1] měrný objem syté páry při teplotě na sání ti; φ [1] relativní vlhkost vzduchu. Odvození této rovnice je uvedeno v Příloze 8.Měrný objem sytých par v Rovnici 15 je funkcí teploty v''=f(t), proto je možné sestrojit nomogram pro určení množství vyloučeného kondenzátu z komprimovaného a ochlazeného vlhkého vzduchu jako funkci teploty na sání, viz Nomogram 16. |
16: Nomogram pro přibližné stanovení množství vyloučeného kondenzátu z komprimovaného a ochlazeného vlhkého vzduchu mC [g·m-3]; ε [1]; ti [°C] <>teplota; φi [%] |
ÚlohyÚloha 1:
Turbokompresor nasává vzduch o teplotě 15 °C a tlaku 0,1013 MPa, na výstupu z turbokompresoru má vzduch teplotu 293 °C a tlak 0,802 MPa. Určete přídavné ztráty, součinitel přídavných ztrát a vnitřní účinnost ηi. Turbokompresor má 12 stupňů. Komprese je nechlazená, respektive uvažujte adiabatickou kompresi. Řešení úlohy je uvedeno v Příloze 1.
Užití lineární aproximace termodynamických změn v T-s diagramu k přibližnému určení velikosti přídavných ztrát při kompresi T [K]; t [°C]; s [J·K-1·kg-1]
Popisek symbolů je uveden v Příloze 1.
Úloha 2:
Stanovte vnitřní izoentropickou, polytropickou a izotermickou účinnost turbokompresoru, který komprimuje suchý vzduch. Vstupní teplota vzduchu je 14,34 °C, výstupní je 480,6 °C. Tlak na vstupu je atmosférický a kompresní poměr je 23. Vnitřní příkon turbokompresoru je 12,6 MW. Turbokompresor je vybaven plášťovým chlazením o výkonu 0,8 MW. Řešení úlohy je uvedeno v Příloze 2.
Užití lineární aproximace vratné polytropické komprese v T-s diagramu k přibližnému určení stavu epol T [K]; t [°C]; s [J·K-1·kg-1]; q [kJ·kg-1]
Popisek symbolů je uveden v Příloze 2.
|
Úloha 3:
Stanovte vnitřní účinnost turbokompresoru, který komprimuje suchý vzduch. Vstupní teplota vzduchu je 14,34 °C, výstupní je 156,6 °C. Tlak na vstupu je atmosférický a kompresní poměr je 23. Vnitřní příkon turbokompresoru je 10,6 MW. Turbokompresor je vybaven dvěma mezichladiči na tlakových úrovní 0,7 MPa a 1,4 MPa. Chladící výkon chladičů je 6,5 MW. Řešení úlohy je uvedeno v Příloze 3.
Popisek symbolů je uveden v Příloze 3.
Úloha 4:
Proveďte základní návrh rozměrů rotoru jednostupňového radiálního turbokompresoru s axiálním vstupem (bez předřazených lopatek). Lopatky rotoru mají na výstupu radiální směr, viz obrázek. Parametry suchého vzduchu na vstupu do rotoru jsou: 101,33 kPa, 15 °C. Požadovaný tlak ze statorové řady lopatek je 0,44 MPa. Požadovaný hmotnostní tok je 0,7225 kg·s-1. Řešení úlohy je uvedeno v Příloze 4.
ΔWθ [m·s-1] odchylka obvodové složky relativní rychlosti na výstupu z rotoru způsobená protiběžným vírem (skluz). |
Popisek symbolů je uveden v Příloze 4.
OdkazyŠKORPÍK, Jiří, 2023, Proudění plynů a par tryskami, fluid-dynamics.education, Brno, https://fluid-dynamics.education/proudeni-plynu-a-par-tryskami.html.
ŠKORPÍK, Jiří, 2024, Technická termomechanika, engineering-sciences.education, Brno, https://engineering-sciences.education/technicka-termomechanika.html.
DIXON, S., HALL, S., 2010, Fluid Mechanics and Thermodynamics of Turbomachinery, Elsevier, Burlington, ISBN 978-1-85617-793-1.
JAPIKSE, David, 1997, Introduction to turbomachinery, Oxford University Press, Oxford, ISBN 0-933283-10-5.
KADRNOŽKA, Jaroslav, 2004, Tepelné turbíny a turbokompresory, Akademické nakladatelství CERM, s.r.o., Brno, ISBN 80-7204-346-3.
©Jiří Škorpík, LICENCE
|