VYSOKÉ UČENÍ TECHNICKÉ V BRNĚ BRNO UNIVERSITY OF TECHNOLOGY
FAKULTA STROJNÍHO INŽENÝRSTVÍ ENERGETICKÝ ÚSTAV FACULTY OF MECHANICAL ENGINEERING ENERGY INSTITUTE
NÁVRH KOMPRESORU TURBODMYCHADLA COMPRESSOR TURBOCHARGER DESIGN
DIPLOMOVÁ PRÁCE MASTER’S THESIS
AUTOR PRÁCE
Bc. DAVID DOMANSKÝ
AUTHOR
VEDOUCÍ PRÁCE SUPERVISOR
BRNO 2011
doc. Ing. JAN FIEDLER, Dr.
-2-
Vysoké učení technické v Brně, Fakulta strojního inženýrství Energetický ústav Akademický rok: 2010/2011
ZADÁNÍ DIPLOMOVÉ PRÁCE student(ka): Bc. David Domanský který/která studuje v magisterském navazujícím studijním programu obor: Energetické inženýrství (2301T035) Ředitel ústavu Vám v souladu se zákonem č.111/1998 o vysokých školách a se Studijním a zkušebním řádem VUT v Brně určuje následující téma diplomové práce: Návrh kompresoru turbodmychadla v anglickém jazyce: Compressor Turbocharger Design Stručná charakteristika problematiky úkolu: Mez pumpovani u kompresoru. Studie proveditelnosti numerickeho modelovaní meze pumpovaní u předem zvoleného kompresoru pomoci dostupných CFD nástrojů.
Cíle diplomové práce: -Studium problematiky přeplňování spalovacích motorů a meze pumpování kompresoru -Návrhový výpočet předem zvoleného kompresoru s využitím CFD nástrojů a sestavení kompresorové mapy -Diskuse možnosti konstrukčních úprav kompresoru s cílem posunutí meze pumpování k nižším průtokům -Apikace CFD výpočtu na vlastní konstrukční úpravu kompresoru a srovnání se stávajícím řešením
Seznam odborné literatury: Kousal, M.: Spalovací turbiny, SNTL Praha 1980 Firemní literatura
Vedoucí diplomové práce: doc. Ing. Jan Fiedler, Dr. Termín odevzdání diplomové práce je stanoven časovým plánem akademického roku 2010/2011. V Brně, dne 20.10.2010 L.S.
_______________________________ doc. Ing. Zdeněk Skála, CSc. Ředitel ústavu
_______________________________ prof. RNDr. Miroslav Doupovec, CSc. Děkan fakulty
ABSTRAKT Diplomová práce se zabývá kompresory turbodmychadel používaných pro přeplňování spalovacích motorů. Cílem práce je studie proveditelnosti numerického modelovaní meze pumpování v kompresoru. V úvodní části je stručně popsán princip a účel přeplňování spalovacích motorů a přehled kompresorů, které se pro přeplňování používají. Hlavní část práce se zabývá nestabilitami proudění v kompresoru a jejich dopadem na pracovní oblast kompresoru. Uvedeny jsou zde některé možnosti rozšíření této pracovní oblasti. Praktická část obsahuje zjednodušený výpočet stupně radiálního kompresoru a analýzu výpočtu proudění pomocí CFD simulace.
ABSTRACT The master's thesis deals with compressor turbochargers used for supercharging of combustion engines. The aim of the thesis is feasibility study of compressor surge limit simulation. First of all the principle and purpose of supercharged combustion engines is described and the compressors' types using for supercharging are mentioned as well. The main part of the thesis is focused on flow instabilities and their influence on compressor operating zone. There are mentioned some possibilities of operating zone extension. The practical part comprises the simplified computation of a one stage of centrifugal compressor and flow analysis by means of CFD simulation.
KLÍČOVÁ SLOVA Turbodmychadlo, radiální kompresor, pumpování, oběžné kolo, CFD, charakteristika kompresoru
KEY WORDS Turbocharger, Centrifugal Compressor, Surge, Impeller, CFD, Compressor Map
-5-
-6-
BIBLIOGRAFIKÁ CITACE DOMANSKÝ, D. Návrh kompresoru turbodmychadla. Brno: Vysoké učení technické v Brně, Fakulta strojního inženýrství, 2011. 82 s. Vedoucí diplomové práce doc. Ing. Jan Fiedler, Dr.
-7-
-8-
PROHLÁŠENÍ O PŮVODNOSTI Prohlašuji, že jsem diplomovou práci vypracoval samostatně za použití literatury uvedené v seznamu použitých zdrojů. V Brně dne: ………………….
Podpis: …………………
-9-
- 10 -
PODĚKOVÁNÍ Děkuji vedoucímu diplomové práce doc. Ing. Janu Fiedlerovi, Dr. a konzultantovi Ing. Milanu Nejedlému, PhD. za cenné rady a podněty, které mi byly poskytnuty při řešení této práce.
- 11 -
- 12 -
OBSAH 1.
ÚVOD .......................................................................................................................... - 15 -
2.
PRINCIP PŘEPLŇOVÁNÍ SPALOVACÍCH MOTORŮ........................................... - 16 -
3.
KOMPRESORY POUŽÍVANÉ PRO PŘEPLŇOVÁNÍ ............................................. - 18 -
4.
3.1
Mechanicky poháněné kompresory ...................................................................... - 18 -
3.2
Turbodmychadla ................................................................................................... - 20 -
RADIÁLNÍ KOMPRESOR ......................................................................................... - 21 4.1
Radiální kompresorový stupeň ............................................................................. - 21 -
4.2
Pracovní proces v radiálním stupni kompresoru................................................... - 22 -
5.
CHARAKTERISTIKA RADIÁLNÍHO KOMPRESORU .......................................... - 23 -
6.
ROTUJÍCÍ ODTRŽENÍ PROUDU .............................................................................. - 26 -
7.
AERODYNAMICKÉ UCPÁNÍ V ČÁSTECH KOMPRESORU ............................... - 27 -
8.
9.
7.1
Ucpání v sacím hrdle ............................................................................................ - 27 -
7.2
Ucpání v oběžném kole ......................................................................................... - 27 -
7.3
Ucpání v difuzoru ................................................................................................. - 28 -
PUMPOVÁNÍ KOMPRESORU .................................................................................. - 29 8.1
Iniciace pumpování kompresoru ........................................................................... - 29 -
8.2
Rozdělení typu pumpování ................................................................................... - 30 -
8.3
Rozdělení meze pumpování .................................................................................. - 31 -
DETEKCE MEZE PUMPOVÁNÍ ............................................................................... - 33 9.1
Detekce meze pumpování dle teploty ................................................................... - 33 -
9.2
Detekce meze pumpování dle tlaku ...................................................................... - 35 -
9.3
Detekce meze pumpování dle vstupního rychlostního profilu ............................. - 36 -
10. DISKUSE NĚKTERÝCH NÁVRHŮ NA ROZŠÍŘENÍ STABILNÍ OBLASTI KOMPRESORU ........................................................................................................... - 38 10.1 Diskuse vlivu tvaru lopatek kompresoru .............................................................. - 38 10.2 Diskuse vlivu obtokového kanálu na vstupu do oběžného kola ........................... - 38 10.3 Diskuse vlivu lopatek s výřezy ............................................................................. - 39 10.4 Diskuse vlivu bezlopatkového difuzoru ............................................................... - 41 10.5 Diskuse vlivu natáčivých vstupních a difuzorových lopatek ................................ - 42 -
- 13 -
11. ÚVODNÍ NÁVRH GEOMETRIE KOMPRESORU ................................................... - 44 11.1 Základní rozměry oběžného kola .......................................................................... - 44 11.2 Vytvoření modelu oběžného kola v BladeGen...................................................... - 45 12. 1-D VÝPOČET STUPNĚ KOMPRESORU PO STŘEDNÍM PROUDOVÉM VLÁKNU ..................................................................................................................... - 49 12.1 Stav proudu vzduchu na vstupu do oběžného kola ............................................... - 50 12.2 Stav proudu vzduchu na výstupu z oběžného kola................................................ - 52 12.3 Návrh ztrátového modelu pro stanovení izoentropické účinnosti ......................... - 54 12.4 Výsledná charakteristika kompresoru pro 1-D výpočet ........................................ - 57 13. 2-D VÝPOČET STUPNĚ KOMPRESORU POMOCÍ VistaTF .................................. - 60 13.1 Příprava modelu části kola v programu DesignModeler ....................................... - 60 13.2 Sumář nastavení výpočtu ve VistaTF .................................................................... - 61 13.3 Editace konfiguračních souborů ............................................................................ - 62 13.5 Výpočet otáčkové větve charakteristiky................................................................ - 63 13.6 Výpočet nestabilní oblasti ..................................................................................... - 64 13.7 Výpočet ucpání kompresoru .................................................................................. - 68 13.4 Výsledná charakteristika kompresoru pro 2-D výpočet ........................................ - 70 14. NÁVRH ÚPRAVY GEOMETRIE KOMPRESORU PRO ROZŠÍŘENÍ STABILNÍ OBLASTI...................................................................................................................... - 72 14.1 Úprava geometrie kompresoru v programu BladeGen.......................................... - 72 14.2 Výsledná otáčková větev upraveného oběžného kola ........................................... - 74 15. ZÁVĚR .......................................................................................................................... - 75 SEZNAM POUŽITÝCH ZDROJŮ ..................................................................................... - 76 SEZNAM HLAVNÍCH POUŽITÝCH OZNAČENÍ .......................................................... - 80 -
- 14 -
1.
ÚVOD
Spalovací motory mají za sebou již poměrně dlouho historii. Již od doby vynálezu spalovacího motoru bylo snahou inženýrů zvyšovat jeho výkon. Jedním ze způsobů, kterým lze zvýšení výkonu dosáhnout, je zvětšení objemu spalovacího prostoru válce. Tento způsob ovšem vedle požadovaného zvýšení výkonu s sebou přináší nevýhodu velkého zastavěného prostoru motorem, větší hmotnosti a s tím spojené ceny. Odlišnou cestou, jak dosáhnout zvýšení výkonu motoru při zachování stávající velikosti, je zvýšení jeho účinnosti. Způsobem, kterým toho lze docílit, je použití turbodmychadel využívajících pro svoji činnost energii obsaženou ve výfukových spalinách. Primární funkcí turbodmychadla je přivádění většího množství vzduchu do spalovacího prostoru válce a v důsledku toho přivádění většího množství paliva. Tím dochází ke zvýšení efektivního tlaku na píst a ke zvýšení motoru bez potřeby zvyšování otáček nebo zvětšování zdvihového objemu motoru. Vzhledem k tomu, že je pro pohon turbíny využito odpadního tepla spalin, dochází také ke zvýšení termické účinnosti a ke snížení měrné spotřeby paliva. Konstrukce turbodmychadel pro motory malých výkonů, které jsou používané v osobních automobilech, představuje v současné době nejvíce kombinaci jednostupňové radiálně axiální turbíny a jednostupňového radiálního kompresoru. Vývoj každého lopatkového stroje je dlouhým a náročný procesem. Jistou úsporu jak časovou, tak ekonomickou představuje využití CFD simulací proudění ve stupni lopatkového stroje. Na základě výpočtů provedených v počítačových programech lze provést optimalizaci geometrie a funkčnosti stroje již v raných projekčních fázích a odpadá tak zdlouhavá a finančně náročná výroba prototypů. Výsledky výpočtů zároveň umožňují analýzu veličin, které by jinak byly experimentálně neměřitelné.
- 15 -
2.
PRINCIP PŘEPLŇOVÁNÍ SPALOVACÍCH MOTORŮ
Přeplňování je velice efektivní způsob, kterým lze ovlivňovat výkonové, ekonomické, provozní a konstrukční parametry spalovacích motorů. Princip přeplňování a vliv na užitečný výkon motoru lze jednoduše odvodit ze vztahu 2.1. ܲ ൌ ܳሶ ή ߟ ൌ ܪ௨ ή ݉ሶ ή ߟ
2.1
Hodnota spodní výhřevnosti paliva ܪ௨ pro benzin Natural 95 je 43,50 MJ×kg-1 a pro motorovou naftu 42,5 MJ×kg-1 [12]. Ke spálení určitého množství paliva ve válci, je potřeba do válce přivést tomu odpovídající množství spalovacího vzduchu. Vzájemný vztah obou hmotnostních toků je dán výrazem 2.2. ݉ሶ ൌ
݉ሶ௦ ߣ ή ܮ௧
2.2
Součinitel přebytku vzduchu ve válci ߣ vyjadřuje tzv. bohatost směsi. Při ߣ ൏ ͳ se jedná o směs bohatou (nedostatek vzduchu). V tomto režim dochází k podstechiometrickému spalování a k tvorbě produktů nedokonalého spalování, zejména CO. Z tohoto důvodu a vzhledem k nevyužití veškerého potenciálu paliva, je práce v tomto režimu nežádoucí. Při ߣ ൌ ͳ jde o stechiometrické složení směsi, při kterém pracuje převážné většina dnešních vozidlových zážehových motorů. Při ߣ ͳ je směs chudá, s přebytkem vzduchu, s kterou pracují všechny zážehové motory [13]. Hmotnostní tok vzduchu, zůstávajícího ve válci pro spalování, je pomocí parametrů motoru vyjádřen vztahem 2.3. ݉ሶ௦ ൌ ܸ ή
݊ெ ߝ ሺߙ െ ͳሻ ή ߩ௦ ή ߟ ή Ͳ ή ݅ ߝെͳ
2.3
V tomto vztahu představuje kompresní poměr ߝ poměr objemu válce při dolní úvrati pístu k objemu válce při horní úvrati. Moderní zážehové motory mají rozsah kompresního poměru od 8 do 11, vznětové pak od 12 do 24. Motory s turbodmychadly mají obvykle menší kompresní poměr než motory atmosférické z důvodů omezení materiálovými vlastnostmi, technologií a kvalitou paliva. Plnící účinnost ߟ je definována vztahem 2.4 ߟ ൌ
ܸ ή ܶ ή ଵ ݉ ൌ ܸு ή ߩ ܸு ή ܶଵ ή
2.4
Jedná se o poměr hmotnosti náplně válce ݉ při stavu na konci sání ܶଵ , ଵ k hmotnosti, která by byla ve válci při ideálním sacím zdvihu při stavu plnícího vzduchu ܶ , . Na obr. 2.1 je p-V diagram průběhu výměny náplně válce čtyřdobého přeplňovaného motoru, z kterého je zřejmý význam uvedených veličin.
- 16 -
Obr. 2.1 Průběh výměny náplně ve válci [8] Kombinací předcházejících vztahů je potom možné vztah pro efektivní výkon motoru popsat rovnicí 2.5. ܲ ൌ ܸ ή
݊ெ ୳ ߝ ሺߙ െ ͳሻ ή ή ߩ௦ ή ߟ ή ߟ ή Ͳ ή ݅ ߣ ή ܮ௧ ߝെͳ
2.5
Z rovnice vyplývá, že na výkon motoru má přímý vliv hustota vzduchu ߩ௦ . Principem přeplňování je tedy dodání větší hmotnosti vzduchu do spalovacího prostoru válce zvýšením hustoty ߩ௦ . Ze stavové rovnice plynu je hustota vzduchu vyjádřena vztahem 2.6. ߩ௦ ൌ
௦ ܴ ή ܶ௦
2.6
Hustotu je tedy možné zvýšit zvýšením tlaku a také snížením teploty vzduchu v sání motoru. Zvýšení tlaku je realizováno stlačením vzduchu v kompresoru. Pro snížení teploty se používá chladič plnicího vzduchu, který je umístěn za kompresorem. [4], [8]
- 17 -
3.
KOMPRESORY POUŽÍVANÉ PRO PŘEPLŇOVÁNÍ
Kompresory používané pro přeplňování spalovacích motorů je možná rozdělit do dvou skupin na kompresory mechanicky poháněné od motoru a turbokompresory.
3.1
Mechanicky poháněné kompresory
Při tomto způsobu přeplňování je kompresor spojen mechanickou vazbou přímo s klikovou hřídelí motoru. Kompresor tak odebírá část výkonu motoru pro pokrytí vlastní spotřeby energie, což je hlavní nevýhodou tohoto pohonu. Mechanicky poháněné kompresory tak výrazně snižují mechanickou účinnost přeplňovaného motoru, k čemuž je nutné přičíst ještě ztráty v pružných členech a spojkách, které jsou zařazeny do pohonu pro případ náhlé změny otáček motoru. Mezi další nevýhody patří vysoká cena, větší hmotnost a hluk. Hlavní výhodou je velice rychlá odezva na změny výkonu. Protože je kompresor mechanicky spojen s klikovou hřídelí, každá změna otáček motoru se okamžitě projeví na otáčkách kompresoru. [8], [14] Pístové kompresory Tyto kompresory jsou nejstarším typem, který se používá pro přeplňování. Jejich výhodou je, že se s otáčkami téměř nemění izoentropická účinnost. Nevýhodou je složitá konstrukce, velké rozměry, při práci se dostává olej do stlačovaného vzduchu. Kvůli špatnému vyvážení setrvačných sil se nehodí pro vysoké otáčky. Uplatnění těchto kompresorů se nachází u velkých dvoudobých motorů, které mají velmi nízké otáčky. Rootsovo dmychadlo V současné době velice rozšířené dmychadlo používané například automobilkou Mercedes-Benz v řadách C, E a SLK. Dmychadlo je tvořeno párem dvoj nebo třízubých rotorů. Zvláštností oproti jiným kompresorům je, že ke kompresy vzduchu, nejprve dopravovaného v prostoru mezi zuby, dochází až při spojení mezizubového prostoru s výtlačným potrubím. Tento způsob komprese s sebou ale přináší nevýhodu značného hluku a tlakových pulzací ve výtlaku. Pro omezení tlakových pulzací se proto používá rotorů se zakroucenými zuby (obr. 3.1). Další nevýhodou je malé dosažitelné poměrné stlačení při nízkých otáčkách motoru, způsobené ztrátou mezerou mezi zuby rotoru a skříní. Výhodou je poměrně malá závislost izoentropické účinnosti na otáčkách, možnost práce při vysokých otáčkách, jednoduchá konstrukce a s tím spojená nízká cena. [8], [15]
- 18 -
Obr. 3.1 Rootsovo dmychadlo [27] Lysholmovo šroubové dmychadlo Obdobně jako Rootsovo se i dmychadlo Lysholmovo skládá ze dvou rotorů. Rotory mají různý příčný průřez. Jednotlivé tvary průřezů se vinou po šroubovici po délce rotoru. Vzduch je na jedné straně nasáván do mezirotorových prostorů, kde se při pootočení uzavírá a následně při axiálním pohybu dochází k jeho stlačení. Lysholmovo dmychadlo má malou změnu účinnosti s otáčkami a lze s ním dosahovat vyšších stlačení než s Rootsovým dmychadlem. Nevýhodou je náročnost výroby a z důvodů velkých poměrných stlačení i citlivost na zvětšení vůlí mezi rotory. Tento typ kompresoru používá například automobilka Mazda.
Obr 3.2 Lysholmovo dmychadlo[27] - 19 -
Proudové kompresory Proudové kompresory jsou další možností při mechanickém přeplňování. Dle směru vystupujícího proudu je lze rozdělit na kompresory axiální a radiální. Axiální kompresory se pro přeplňování v současné době již téměř nepoužívají z důvodu relativně malého stlačení vzduchu ve stupni. Malé poměrné stlačení lze řešit použitím vícestupňového kompresoru, nicméně takové řešení je prostorově náročné a vylučuje proto použití v osobních automobilech. Radiální kompresory nabízejí oproti axiálním výhodnější parametry, a proto jsou při mechanickém přeplňování využívány daleko častěji. Stupeň radiálního kompresoru je podrobněji popsán v následujících kapitolách.[8]
3.2
Turbodmychadla
Historie přeplňování turbodmychadly je spjata se jménem a patenty švýcarského inženýra Alfreda Büchiho. První patent popisující princip přeplňovaného spalovacího motoru turbodmychadlem obdržel Büchi již v roce 1905. Na uvedení do praxe ovšem čekal až do roku 1925, kdy tato technologie našla využití ve velkých lodních motorech. V automobilovém průmyslu začala být tato technologie zajímavá až v polovině 50. let, když již bylo dosaženo dostatečné kompaktnosti, spolehlivosti a cenové přijatelnosti motorů s turbodmychadlem. Kompresor turbodmychadla je poháněn turbínou, která je uložena na společné hřídeli (obr. 3.3). Zdrojem energie pro pohon turbíny jsou výfukové plyny, jež vychází se spalovacího prostoru válce. Tento způsob spolupráce má oproti mechanicky poháněným kompresorům výhodu v tom, že kompresor neodebírá výkon motoru a ke svému chodu využívá odpadního tepla spalin. V případě dobrého návrhu a volby turbodmychadla připojeného k motoru lze dosáhnout ke zvýšení celkové termické účinnosti okolo 5 %.
Obr. 3.3 Schematický nákres turbodmychadla [28] - 20 -
4.
RADIÁLNÍ KOMPRESOR
K přeplňování spalovacích motorů se nejvíce používají kompresory radiální. Svou pozici si získali zejména díky relativně nízkým výrobním nákladům, jednoduchou konstrukcí, malou hmotností a schopností v jednom stupni realizovat velké stlačení vzduchu.
4.1
Radiální kompresorový stupeň
Radiální kompresorový stupeň je tvořen (obr. 4.1) stacionární vstupní částí (0 - 1), oběžným radiálním kolem (1 – 2) a na ně navazujícím statorem. Stator se obvykle dělí na bezlopatkový (2 – 3) a lopatkový difuzor (3 - 5). V případě jednostupňového kompresoru následuje za difuzorem výstupní spirální skříň. Ve vstupní části stupně dochází k urychlení případně k usměrnění proudu. Vzduch následně vstupuje do oběžného kola, ve kterém se energie vzduchu přivádí a zvyšuje se jeho tlak a rychlost. V difuzoru se poté snižuje kinetická energie proudu na výstupu z oběžného kola na energii tlakovou. Výstupní spirální skříň nakonec proud sbírá a svádí ho k jednomu nebo několika výstupním hrdlům. Zároveň dochází ve spirální skříni k dalšímu snižování kinetické energie proudu na energii tlakovou. [2], [29]
Obr. 4.1 Schéma stupně radiálního kompresoru s vyznačením vztažných rovin [9]
- 21 -
4.2
Pracovní proces v radiálním stupni kompresoru
Pracovní proces v jednostupňovém radiálním kompresoru lze znázornit v i – s diagramu, tak jak je na obr. 4.2. Bod 01 představuje vzduchu vstupujícího do kompresoru z okolí. Urychlování vzduchu ve vstupní části způsobí pokles tlaku z ଵ na ଵ v případě izoentropického proudění a z na ଵpokud jsou uvažovány ztráty ve vstupní části. Změna v entalpii odpovídá vzrůstu kinetické energie (ܿଵଶ Τʹ). Energie přivedená vzduchu oběžným kolem je označena spojnicí bodu 1 a 2. Tomu odpovídající izoentropický proces je označen 1 – 2iz. Pokud by se uvažovalo úplná izoentropická změna kinetické energie vystupující z oběžného kola na energii tlakovou, byla by velikost výstupního tlaku ଶ . Protože je ale proces v difuzoru se ztrátami (2 - 5) a na konci difuzoru zůstává proud s určitou kinetickou energií (ܿହଶ Τʹ), je skutečný statický tlak za difuzorem ହ . [9]
Obr. 4.2 i - s diagram radiálního kompresorového stupně [9]
- 22 -
5.
CHARAKTERISTIKA RADIÁLNÍHO KOMPRESORU
U vozidlových spalovacích motorů pracuje kompresor turbodmychadla velice často v režimech, které se výrazně odlišují od návrhového bodu kompresoru. Jedná se především o různé otáčky kompresoru, změně hmotnostního toku a odlišné termodynamické stavy vzduchu v sání kompresoru. V závislosti na těchto podmínkách se mění parametry, které charakterizují provozní režim kompresoru – poměrné stlačení a účinnost. [3] Látkové vlastnosti stlačovaného média (tepelná vodivost, měrná tepelná kapacita, vazkost, plynová konstanta) mají obvykle na provozní režim kompresoru nepatrný vliv, a proto se při sestavování kompresorové charakteristiky zanedbávají. I přes toto zjednodušení je stav pracovní látky za kompresorem při změně provozního režimu funkcí čtyř proměnných 5.1 a 5.2. ଶ ൌ ݂ଵ ሺଵ ǡ ܶଵ ǡ ݊ǡ ݉ሶሻ
5.1
ߟ ൌ ݂ଶ ሺଵ ǡ ܶଵ ǡ ݊ǡ ݉ሶሻ
5.2
Tyto funkční závislosti je pro radiální kompresory velice obtížné stanovit výpočtem. Z tohoto důvodu se hodnoty funkcí získávají především experimentální cestou, nicméně v posledních letech je zřejmý sílící význam CFD výpočtů. Pro vynesení funkcí do rovinného grafu je nutné zvolit dvě veličiny jako konstantní parametr grafu. Funkce poté přejdou do tvaru 5.3 a 5.4. ଶ ሺଵ ǡ ܶଵ ሻ ൌ ݂ଵǡ ሺ݊ǡ ݉ሶሻ
5.3
ߟ ሺଵ ǡ ܶଵ ሻ ൌ ݂ଶǡ ሺ݊ǡ ݉ሶሻ
5.4
Výsledné funkce je poté možné zobrazit v grafu, přičemž tento graf se nazývá absolutní charakteristika kompresoru (obr. 4.1). [2]
Obr. 4.1 Absolutní charakteristika kompresoru [3] Tato charakteristika vyjadřuje závislost poměrného stlačení P୩ výstupního a vstupního tlaku na hmotnostním průtoku ݉ሶ. Je nutné zdůraznit, že se jedná o poměrné stlačení celkových - 23 -
tlaků, tedy P୩ ൌ ଶ Ȁଵ . Celkový tlak je určen statickou ௦௧ a dynamickou složkou ௗ , viz vztah 5.5. ൌ ௦௧ ௗ ൌ ௦௧
ܿଶ ߩ ʹ
5.5
V charakteristice jsou také vyneseny křivky účinnosti ߟ௭ a čáry konstantních otáček ݊ ൌ ݇ݐݏ݊, které se nazývají větve charakteristiky. Absolutní charakteristika je velice názorná, ale závislost na konkrétních hodnotách termodynamického stavu v sání ଵ ǡ ܶଵ , podstatně snižuje její použitelnost v praxi. Bylo by nutné proměřovat velké množství takovýchto charakteristik pro různé hodnoty vstupních tlaků a teplot, aby bylo dosaženo reprezentativního zobrazení provozních vlastností kompresoru. [3]
Vzhledem k této nevýhodě je užitečné sestavit charakteristiku, kterou lze použít i při změně stavu v sání kompresoru. Takováto charakteristika se nazývá univerzální charakteristika kompresoru (obr. 4.2).
Obr. 4.2 Univerzální charakteristika kompresoru [3] Aby bylo možné takovou charakteristiku vytvořit, je nutné splnit podmínky fyzikální podobnosti proudění v kompresoru. Tato podobnost je dána rovností kriteriálních parametrů podobnosti: geometrická podobnost, Reynoldsovo číslo Re, Machovo číslo Ma, Freudovo číslo Fr a izoentropický exponent k. Geometrická podobnost je pro daný kompresor přirozeně splněna. Vzhledem již výše zmíněnému zjednodušení, při kterém lze zanedbat vliv látkových vlastností stlačovaného média není vyžadována rovnost Re, Fr a k. Jediným rozhodujícím kritériem je rovnost Machova čísla v odpovídajících průřezech kompresoru [29]. Při respektování rovnosti Machova čísla je poté charakteristika vyjádřena funkcemi 5.6 a 5.7. ߎ ൌ ݂ଵ ቆ ߟ௭ ൌ ݂ଶ ቆ
݉ሶ ή ඥܶଵ ݊ ǡ ቇ ଵ ඥܶଵ ݉ሶ ή ඥܶଵ ݊ ǡ ቇ ଵ ඥܶଵ
- 24 -
5.6
5.7
Podrobné odvození komplexních parametrů „bezrozměrového“ průtoku a „bezrozměrových" otáček lze najít v [2]. Univerzálních charakteristika má všeobecnější platnost než absolutní charakteristika, nicméně použití bezrozměrných parametrů přináší i jednu nevýhodu. Bezrozměrové parametry jsou čísla, která nedávají přímou představu o technických hodnotách n a ݉ሶ. Proto se někdy charakteristika vyjadřuje pomocí tzv. redukovaných hodnot nr a ݉ሶ , které jsou definovány vztahy 5.8 a 5.9. ܶଵ ݉ሶ ൌ ݉ሶ ή ඨ ή ܶ ଵ ܶ ݊ ൌ ݊ ή ඨ ܶଵ
5.8
5.9
kde a ܶ jsou standardní vstupní podmínky, jejichž hodnoty jsou obvykle brány podle norem ISO. Charakteristika se potom znázorňuje v souřadnicích ݉ሶ a ߎ jako soubor křivek ݊ [4]. Její příklad je na obr. 4.3.
Obr. 4.3 Charakteristika kompresoru s redukovaným průtokem a otáčkami [5] Charakteristika kompresoru tedy zobrazuje stabilní oblast práce kompresoru při různých provozních režimech. Tato charakteristika je z levé strany oddělena od nestabilní oblasti tzv. mezí pumpování. Vlevo od této meze, tedy ve směru se zmenšujícím se hmotnostním tokem až po nulový hmotnostní průtok, je provoz kompresoru nestabilní. Příčinou jsou dva úkazy nestability a to rotující odtržení proudu a pumpování. Pravá strana pracovní oblasti kompresoru je také omezena a to tzv. aerodynamickým ucpáním kompresoru. - 25 -
6.
ROTUJÍCÍ ODTRŽENÍ PROUDU
Rotující odtržení proudu je nestabilita, která zasahuje hmotnostní proud v jedné nebo více zónách, v kterých dochází ke kmitání. Kmitající zóny rotují po obvodu rotoru, přičemž počet zón se při dalším postupu do oblasti nestability charakteristiky kompresoru mění. Jednotlivé zóny postihují nejprve proud na vnější straně lopatkové mříže, postupně se rozšiřují po jejím obvodu a zachvacují proud po délce lopatek, schematicky na obr. 5.1. V závislosti na velikosti jednotlivých zón se jejich obvodová rychlost rotace pohybuje mezi 20 – 80 % rychlosti rotoru, přičemž větší zóny rotují s menší rychlostí než zóny malé. Frekvence kmitání není závislá na velikosti a tvaru potrubí, na které je kompresor připojen. Tento režim proudění se tedy od pumpování kompresoru liší především tím, že nezasahuje celý proudu ale pouze jeho část a zároveň hmotnostní tok rotorem zůstává stále kladný. [8], [7], [16]
Obr. 5.1 Schematické znázornění zón kmitání Vznik rotujícího odtržení proudu lze popsat na kompresorové mříži obr. 5.2. Se snižujícím se průtokem dochází ke zvětšování kladného úhlu náběhu i, dochází k odtrhávání proudu na podtlakové straně lopatky. Dojde-li v kanálu a k odtržení proudění dříve než v sousedních kanálech, což je způsobeno například nerovnoměrností proudu před mříží, či drobnou odchylkou v geometrii lopatek, ucpe se vírem částečně kanál a, zatímco v sousedních kanálech se průtok zvýší. Tím dojde ke zvětšení úhlu náběhu pro kanál c a ke zmenšení pro kanál b. To vede k odtržení proudu v c a postupné stabilizaci proudění v kanálu b. Odtržení proudu se tedy postupně posouvá proti směru hodinových ručiček a dochází k rotaci po obvodu rotoru kompresoru. [4]
Obr. 5.2 Iniciace vzniku rotujícího odtržení proudu [9]
- 26 -
7.
AERODYNAMICKÉ UCPÁNÍ V ČÁSTECH KOMPRESORU
Aerodynamické ucpání nastane v případě, kdy průtoková rychlost v části kompresoru dosáhne rychlosti zvuku. Kompresor v tomto případě dosáhl svého průtočného maxima a již není schopen pro dané podmínky zpracovat větší hmotností tok. V radiálním kompresoru nastává ܽܯൌ ͳ nejčastěji v sacím hrdle kompresoru, v hrdle oběžného kola, při vstupu do bezlopatkového difuzoru nebo v nejužším místě kanálu lopatkového difuzoru. Pro případ jednorozměrného izoentropického proudění ideálního plynu jsou stanoveny vztahy pro výpočet hmotnostního toku, při kterém dojde k ucpání.
7.1
Ucpání v sacím hrdle
V sání kompresoru dojde k zahlcení v případě, že je překročen hmotnostní tok vzduchu stanovený vztahem 7.1. ାଵ
ʹ ଶήሺିଵሻ ݉ሶ ൌ ܣ௧ ή ߩ ή ܽ ή ൬ ൰ ߢͳ
(7.1)
Vzhledem k tomu, že ߩ a ܽ jsou podmínky v sání kompresoru, je ucpání v sacím hrdle závislé pouze na podmínkách okolního prostředí a není ovlivněno procesy probíhajícími v oběžném kole.
7.2
Ucpání v oběžném kole
V oběžném kole může dosáhnout proudění rychlosti zvuku nejčastěji v hrdle, což je místo s nejmenší plochou průřezu průtočného kanálu. Toto místo se nachází velice blízko vstupu do oběžného kola (obr. 7.1).
Obr. 7.1 3D model lopatkování oběžného kola s vyznačeným hrdlem - 27 -
Na rozdíl od sání kompresoru nastává ucpání v oběžném kole v případě, že rychlosti zvuku dosáhne relativní rychlost proudu. Maximální hmotnostní tok, který dokáže oběžné kolo zpracovat je potom definován vztahem 7.2. ۇ ݉ሶ ൌ ܣଵ௧ ή ߩଵ ή ܽଵ ή ۈ ۉ
ାଵ ଶ ଶήሺିଵሻ ܷଵ ଶ ۊ ܽଵ
ʹ ሺߢ െ ͳሻ ή
ۋ
ߢͳ
ی
(7.2)
Z tohoto vztahu je zřejmé, že zvýšením otáček oběžného kola se zvýší i hmotnostní tok vzduchu, s kterým může kompresor pracovat, aniž by došlo k zahlcení.
7.3
Ucpání v difuzoru
Vztah pro hmotnostní tok, při kterém dojde k zahlcení ve statické části kompresoru, byl již uveden pro sání kompresoru (7.1). V difuzoru není ovšem uvažováno se stavem vzduchu na vstupu, ale se stavem na výstupu oběžného kola. Stav na výstupu z kola je ale ovlivněn vlastními procesy, které v kole probíhají. Proto vztah 7.1 přechází do tvaru 7.3. ݉ሶ ൌ ܣସ ή ߩଵ ή ܽଵ ή
ͳ ሺߢ െ ͳሻ ή ߟ ή ߪ ή
ܷଶଶ ଶ ൨ ܽଵ
ܷଶ ͳ ሺߢ െ ͳሻ ή ߪ ή ଶଶ ൨ ܽଵ
ǡହ
ାଵ
ʹ ଶήሺିଵሻ ή൬ ൰ ߢͳ
(7.3)
Ve vztahu 7.3 je proces v difuzoru uvažován jako izoentropický a v oběžném kole jako neizoentropický. Zároveň platí pouze pro oběžné kolo s radiálními lopatkami. Z výše uvedených poznatků vyplývá, že hmotnostní tok kompresorem, při kterém dojde k ucpání, závisí na geometrii kompresoru, vlastnostech vzduchu na vstupu do kompresoru a otáčkách oběžného kola. Při provozu lze tedy ucpání kompresoru ovlivnit pouze změnou otáček oběžného kola, přičemž je ale nutné respektovat to, že může dojít k ucpání ve statických částech kompresoru. [9], [10]
- 28 -
8.
PUMPOVÁNÍ KOMPRESORU
Pumpování kompresoru je stav proudění, při kterém se vyskytují výkyvy toku v celém kompresním systému. Výkyvy přitom mohou být tak velké, že se během krátké části pumpovního cyklu obrátí směr proudění. Oscilační zatěžování a odlehčování kompresoru při pumpování také vyvíjí střídavé zatížení ložisek podporujících kompresorovou hřídel, což může mít fatální důsledky na celý kompresní systém. Při pumpování se vytvářejí mechanické vibrace a je zřetelně slyšitelné.
8.1
Iniciace pumpování kompresoru
Iniciaci pumpování lze vysvětlit na kompresorové mapě obr. 8.1a, na které je zobrazena kompresorová charakteristika a charakteristika odběru (škrtícího ventilu) pro různé polohy otevření. Zpočátku je hmotnostní tok kompresoru a stlačení stabilní a odpovídá hmotnostnímu toku a stlačení na škrtícím ventilu ்ܻଵ . V bodě (*), je tedy tok kompresorem absolutně stabilní. To znamená, že při malém snížení hmotnostního toku dojde ke zvýšení hodnoty stlačení kompresoru a ke snížení hodnoty stlačení na škrtícím ventilu. Tok se potom bude urychlovat a zvětšovat, dokud nebude obnovena původní rovnováha.
Obr. 8.1 a) Stabilita kompresoru b) Časová závislost tlakového a průtokového součinitele [6]
V případě, že je ventil uzavřen například na polohu ்ܻଷ , stává se tok nestabilní. Malé snížení hmotnostního průtoku zapříčiní mírné snížení stlačení na kompresoru a prudké na ventilu. Ve výtlačném potrubí mezi kompresorem a ventilem je poté určitý objem vzduchu, jehož tlak je větší než tlak vzduchu dodávaný kompresorem po snížení hmotnostního toku. Kompresor již není schopen pracovat proti zvýšenému tlaku a dochází ke zpomalování proudění nebo dokonce k obrácení směru toku v kompresoru. Vývoj cyklu po iniciaci pumpování je schematicky znázorněn na obr. 8.1b, pro polohy škrtícího ventilu ்ܻସ . Cyklus začíná v místě největšího stlačení a v kompresorové charakteristice se pohybuje ve smyslu proti směru hodinových ručiček. V důsledku snížení hmotnostního toku, poklesne tlak ve výtlačném potrubí. Při určitém tlaku je kompresor schopný pracovat proti tlaku ve výtlačném potrubí a průtok kompresorem se obnoví.
- 29 -
Následně se tlak v kompresním systému zvyšuje do bodu největšího stlačení a cyklus se opakuje.
8.2
Rozdělení typu pumpování
V závislosti na kompresním systému a provozních podmínkách je možné rozlišovat pumpování mírné, klasické, hluboké a modifikované. Mírné – v tomto případě je frekvence oscilací blízká tzv. Helmholtzově frekvenci, tj. rezonanční frekvenci potrubí se vzduchem připojeného ke kompresoru. V tomto režimu nedochází k obrácení směru proudění plynu v kompresoru. Obr. 8.2a, označení 1. Klasické – je charakteristické většími tlakovými oscilacemi s menší frekvencí, které závisí na rotační rychlosti a provozních podmínkách. Za těchto podmínek se dynamika systému stává nelineární, což se projevuje většími harmonickými kmity. Cyklus pro systémy s malým objemem je na obr. 8.2a označen 2. Pro systémy s velkým objemem jsou harmonické kmity více zřejmé a na krátkou chvíli může dojít i ke změně směru hmotnostního toku. Na obrázku označen 3. Hluboké – při tomto typu pumpování dochází k velkým oscilacím tlaku a hmotnostního toku, který se stává v části cyklu negativní, obr. 8.2b. Při obráceném průtoku se okamžitý provozní bod pohybuje v kompresorové charakteristice dolů po křivce nazývané jako záporná charakteristika toku. Tato křivka charakterizuje odpor, který klade lopatková mříž obrácenému toku. V kladném směru toku provozní bod těsně sleduje charakteristiku kompresoru v ustáleném stavu. Frekvence při hlubokém pumpování závisí na plnění a vyprazdňování výtlačného potrubí. Modifikované – pumpování může nastat také v kombinaci s rotačním odtržením proudu. V takovém případě je hmotnostní tok kompresorem nestálý a axiálně nesymetrický, což je způsobeno superpozicí rotačního odtržení a oscilacemi hmotnostního toku v celé lopatkové mříži kompresoru. [6]
Obr. 8.2 a) Mírné a klasické pumpování b) Hluboké pumpování [6]
- 30 -
8.3
Rozdělení meze pumpování
Stacionární mez pumpování Stacionární mez pumpování je v kompresorové charakteristice čárou, která odděluje stabilní oblast práce kompresoru od nestabilní. V celé řadě zdrojů (např. [9]) se mez pumpování definuje jako místo, ve kterém se křivka konstantních otáček stává horizontální, neboli první derivace poměrného stlačení podle hmotnostního toku je nulová. ݀ȫ݇ ݀݉ሶ
ൌͲ
(8.1)
Z kompresorové charakteristiky na obr. 4.3 vyplývá, že je tento přístup do značné míry idealizovaný a konzervativní. Skutečná mez pumpování je v tomto případě posunuta především při vysokých otáčkách směrem k nižším průtokům než stanovuje kritérium 8.1. K určení skutečné meze pumpování je proto zapotřebí použít především experimentálního měření (kap. 9). Poloha stacionární meze není stálá, je závislá na provozních podmínkách a velikosti objemu zapojeného před a za kompresorem. Při zvětšování objemu zapojeného za kompresorem se stacionární mez pumpování posouvá v kompresorové charakteristice směrem vpravo, ke zvětšujícímu se hmotnostnímu toku. [18] Vliv provozních podmínek je patrný obzvláště v automobilech, které pracují v různých klimatických pásmech, teplotách prostředí a nadmořských výškách. Tyto podmínky mají vliv na posun meze pumpování směrem k nižším, ale i vyšším hmotnostním průtokům.
Pokud dojde ke snížení tlaku na vstupu z hodnoty na hodnotu , dojde k posunu meze pumpování směrem k vyšším průtokům (obr. 8.3a). K posunu k vyšším průtokům dojde i v případě pokud se zvýší teplota v sání kompresoru z hodnoty ܶ na hodnotu ܶ (obr. 8.3b). [16]
Obr. 8.3 a) Vliv tlaku v sání b) Vliv teploty v sání [16] Dynamická mez pumpování Dynamická mez pumpování je stejně tak jako statická mez pumpování dělicí čárou od nestabilní části kompresorové charakteristiky vpravo a stabilní vlevo. Rozdíl je - 31 -
v nestacionárním průtoku způsobený připojením kompresoru k přeplňovanému motoru. Motor budí odběry jednotlivých válců kmity podle otáček motoru a počtu válců s různou frekvencí a podle zdvihového objemu válců nebo poměrného stlačení plnícího vzduchu. Tyto nestacionární podmínky ovlivňují stacionární charakteristiku kompresoru tak, že dochází k posuvu meze pumpování v obou směrech vodorovné osy kompresorové charakteristiky a k poklesu stlačení plnícího vzduchu. Dynamickou mez pumpování je nutné sledovat především u motorů s nízkými otáčkami a s velkým zdvihovým objemem válců, kde je vzhledem k zdvihovému objemu válců mezi motorem a kompresorem relativně malý objem. V tomto případě dochází k znatelnému posunu meze pumpování směrem vpravo. U vysokootáčkových motorů s malým zdvihovým objemem válců a relativně velkým objemem mezi kompresorem a motorem je posuv meze pumpování směrem vpravo zanedbatelný, přičemž je možné pozorovat i posuv vlevo směrem k nižším průtokům. [18]
- 32 -
9.
DETEKCE MEZE PUMPOVÁNÍ
Při stanovování stabilní oblasti práce kompresoru je pro určení meze pumpování možné využít úkazů projevujících se v její blízkosti.
9.1
Detekce meze pumpování dle teploty
Jedním z příznaků přibližující se meze pumpování je zvyšující se teplota v sacím potrubí před oběžným kolem kompresoru. V [1] byly prezentovány výsledky experimentu zabývajícím se měřením tohoto jevu. Na obrázku 9.1 je charakteristika kompresoru pro čtyři různé otáčky oběžného kola, na které je pomocí izoterem zobrazena diference mezi teplotou vzduchu na vstupu do kompresoru a teplotou před oběžným kolem.
Obr. 9.1 Charakteristika kompresoru s vyznačením teplotních diferencí [1] Tato charakteristika byla získána pomocí teplotních čidel, jejichž rozmístění je schematicky znázorněno na obr. 9.2. Nejbližší teplotní sonda měří teplotu v 2 mm vzdálenosti a nejvzdálenější je umístěna ve vzdálenosti 93 mm. Graf naměřených teplot pro každou sondu při konstantních otáčkách oběžného kola a různých hmotnostních průtocích je na obr. 9.3. Z grafu je zřejmý výrazně se zvyšující teplotní gradient při nízkém hmotnostním průtoku a blížící se mezi pumpování. Jsou zde také patrné vyšší hodnoty teplot naměřené sondou ve vzdálenosti 19 mm od oběžného kola, než sondou umístěnou nejblíže lopatkové mříže. Příčinou je větší zpětný tok po vnějším obvodu oběžného kola a tedy vyšší teplota vzduchu.
- 33 -
Obr. 9.2 Rozmístění jednotlivých teplotních sond pře oběžným kolem [1]
Obr. 9.3 Zvýšení teploty v různé vzdálenosti od oběžného kola [1] Z obrázků 9.1 a 9.3 vyplývá, že je teplota jasným indikátorem nestabilit vyskytujících se v proudění kompresorem při přibližování se k mezi pumpování. Přičemž na blížící se pumpování kompresoru může upozornit výrazný nárůst gradientu teploty. [1]
- 34 -
9.2
Detekce meze pumpování dle tlaku
Další z indikátorů oznamující blížící se pumpování kompresoru je tlak v difuzoru. Na obr. 9.5a je zobrazena iniciace pumpování. Zavíráním škrtícího ventilu ்ܻ , klesá hmotnostní průtok a při určité hodnotě je iniciováno pumpování, což je zřejmé z grafu poměrného stlačení (obr. 9.5b). Iniciace pumpování je v tomto případě definována jako bod, ve kterém začíná narůstat amplituda průběhu tlaku a tlak začíná oscilovat na jiné frekvenci (obr 9.5c). Ihned po iniciaci následuje oblast rozvoje pumpování (na obr. 9.5c označena B). Tato oblast se vyznačuje menšími tlakovými amplitudami, než se rozšíří v plně vyvinuté pumpování – oblast C. [6]
Obr. 9.4 a) Škrtící ventil b) Poměrné stlačení v závislosti na čase c) Detailní pohled na poměrné stlačení [6] Při měření tlakových oscilací je nutné respektovat odlišný průběh tlaků v různých částech kompresoru. Na obr. 9.6 jsou zobrazeny oscilogramy průběhu tlaků kompresoru s lopatkovým difuzorem na konci oběžného kola, v mezeře mezi kolem a difuzorem, v hrdle difuzoru a na konci difuzoru. Z porovnání průběhu tlaků v různých částech kompresoru vyplývá, že nejdříve dochází k tlakovým oscilacím při iniciaci pumpování v hrdle difuzoru.
- 35 -
Obr. 9.5 Průběh tlaků v různých místech kompresoru [11]
9.3
Detekce meze pumpování dle vstupního rychlostního profilu
Významným indikátorem blížícího se pumpování kompresoru je vstupní rychlostní profil proudění. Jak se kompresor se snižujícím se hmotnostním tokem přibližuje k hranici pumpování, dochází na vstupní části oběžného kola na vnějším průměru k vývoji recirkulační oblasti proudu. Při normálních provozních podmínkách je rychlostní profil na vstupu do oběžného kola v axiálním směru téměř rovnoměrný. V blízkosti meze pumpování je pozorovatelné na vnějším průměru záběrníku snižování axiální rychlosti proudu mnohem výrazněji než na středním průměru. S dalším snižováním průtoku dochází v této oblasti ke změně celkového směru proudění, než nastane vlastní pumpování kompresoru (obr. 9.1). [17]
- 36 -
Obr. 9.6 Rychlost proudu na vnějším průměru záběrníku [17] Rychlostní profil proudění v normálním provozním režimu a v režimu v blízkosti meze pumpování je na obrázku 9.2. Proudění v normálním režimu je zřetelně vyrovnané s nezáporným rychlostním profilem po celém průřezu kanálu oběžného kola. V blízkosti pumpování je na vnějším průměru oběžného kola vyvinutá recirkulační oblast, ve které dochází k obrácení směru proudění. Tento závěr je zároveň v souladu s kap. 9.1.
Obr. 9.7 Rychlostní profily proudění při normálních provozních podmínkách a v blízkosti meze pumpování [20]
- 37 -
10.
DISKUSE NĚKTERÝCH NÁVRHŮ NA ROZŠÍŘENÍ STABILNÍ OBLASTI KOMPRESORU
Pro kompresory turbodmychadel v automobilovém průmyslu je důležitá schopnost pracovat v široké stabilní pracovní oblasti charakteristiky. Velikost pracovní oblasti lze ovlivnit celou řadou aktivních a pasivních zařízení a také vlastním tvarem lopatek kompresoru. Mezi pasivní prvky patří například obtokový kanál, ve kterém se přepouští vzduch mezi záběrníkem a vstupem do kompresoru, nebo také připouštění vzduchu do oblasti bezlopatkového difuzoru. Aktivními prvky jsou především natáčivé lopatky na vstupu do kompresoru a natáčivé lopatky difuzoru. Z důvodů požadavků na malé rozměry a nízkou náročnost výroby turbodmychadel je dávána přednost především pasivním prvkům.
10.1 Diskuse vlivu tvaru lopatek kompresoru Nejvýznamnější vliv na průběh kompresorové charakteristiky má tvar lopatek na výstupu z oběžného kola. Ze tří možných typů oběžných kol – stupeň s dopředu zahnutými lopatkami, s radiálními lopatkami a s lopatkami dozadu zahnutými, nabízí největší šířku stabilní oblasti kolo s lopatkami dozadu zahnutými (ߚଶ ͻͲι). Takovýto tvar lopatek je příznivý pro proudění vzduchu, při kterém dochází k menším ztrátám a tím ke zvýšení účinnosti kompresoru. Oproti kolu s radiálními lopatkami je také daleko méně rozvinutý sekundární vír uvnitř kola, což se projevuje větším součinitelem skluzu a větší rovnoměrností na výstupu z oběžného kola. Nevýhodou je nižší dosahovaná hodnota stlačení oproti lopatkám radiálním nebo dopředu zahnutým. Proto je nutné pro dosažení požadovaného stlačení zvýšit rychlost kompresoru, což zvyšuje nároky na pevnost oběžného kola. [23]
Obr. 10.1 Kompresor s dozadu zahnutými lopatkami [9]
10.2 Diskuse vlivu obtokového kanálu na vstupu do oběžného kola Dalším způsobem rozšíření pracovní oblasti je pomocí obtokového kanálu. Jedná se o kanál vedený skříní kompresoru spojující prostor nad záběrníkem s prostorem před oběžným kolem. Takovéto spojení umožňuje v závislosti na tlakovém spádu přetékání vzduchu mezi těmito dvěma místy (obr. 10.2). - 38 -
V případě nízkých otáček a malého hmotnostního toku je tlak v záběrníku větší než ve vstupu do kompresoru a dochází k přetékání směrem ze záběrníku do přívodu. Důsledkem toho se před oběžným kolem vytvoří vír, který zmenší efektivní průtočnou plochu. Tím se zvýší osová rychlost proudu, což příznivě ovlivní rychlostní trojúhelník před oběžným kolem. Zmenší se úhel náběhu a tím se stabilizuje chod kompresoru. Při vysokých otáčkách a velkých hmotnostních tocích je tlak v přívodu do kompresoru větší než v záběrníku. V tomto případě dojde k odsávání vzduchu z místa před oběžným kolem a přepuštění vzduchu do záběrníku, což příznivě ovlivní proudění v oběžném kole. V oblasti návrhového bodu jsou tlaky téměř vyrovnané, a tudíž nedochází k žádnému přetékání obtokovým kanálem. [23]
Obr. 10.2 Schéma proudění v obtokovém kanále [24] Použití přepouštěcího kanálu má výrazný vliv na rozšíření stabilní oblasti kompresoru především směrem k nižším průtokům a tedy posunu meze pumpování směrem vlevo. Dle experimentu popsaného v [24] bylo dosaženo posunu meze pumpování oproti původní charakteristice o 18 %. Nepříznivě bylo ovšem ovlivněno poměrné stlačení drobným poklesem především při vysokých otáčkách a snížení účinnosti kompresoru přibližně o 1,5 %. Snížení účinnosti je způsobeno třecími ztrátami v obtokovém kanálu a ztrátami při míšení proudění před oběžným kolem.
10.3 Diskuse vlivu lopatek s výřezy V experimentu prezentovaném v [25] je popsán pozitivní vliv na rozšíření stabilní oblasti kompresoru pomocí lopatek s výřezy. Výřezy jsou umístěny v 30 % celkové osové vzdálenosti od výstupní hrany lopatky a představují 0,2 mm široké otvory v lopatkách, vedoucí od náboje až po vnější hranu oběžného kola (obr. 10.3). - 39 -
Umístěním tohoto prvku do místa, ve kterém dochází k odtrhávání proudění na sací straně lopatky, způsobí přetékání vzduchu z tlakové strany lopatky na stranu sací. To má za následek omezení růstu mezní vrstvy a snížení skluzu.
Obr. 10.3 Lopatky s výřezy [25] Srovnání výsledků měření poměrného stlačení pro čtyři různé rychlosti kompresoru s neupravenými lopatkami a lopatkami upravenými výřezy, je uvedeno na obrázku 10.4. Ve všech případech došlo k rozšíření stabilní oblasti k větším hmotnostním průtokům. Rozšíření směrem k nižším hmotnostním průtokům a tedy posunutí meze pumpování lze pozorovat pouze u vysokých otáček oběžného kola. Celkové rozšíření pracovní oblasti je pro otáčky 10 000 min-1 5,3 %, pro střední otáčky 15 000 a 20 000 min-1 je rozšíření zhruba 31 % a pro nejvyšší otáčky 60 %. Rozdíl ve stlačení lze považovat za zanedbatelný vzhledem k možným nepřesnostem v měření.
Obr. 10.4 Kompresorová charakteristika upraveného a neupraveného oběžného kola [25] - 40 -
Srovnání izoentropická účinnosti pro čtyři různé rychlosti upraveného a neupraveného oběžného kola je na kompresorové charakteristice obrázku 10.5. Při nízkých otáčkách dochází při hmotnostním toku 0,03 kg/s až 0,06 kg/s k poklesu účinnosti přibližně o 9 %. K poklesu o 10 % dochází v celém rozsahu při 15 000 min-1. Odlišný trend je zachycen při vyšších otáčkách, kde při určitém hmotnostním průtoku dochází dokonce ke zvýšení účinnosti zhruba o 5 %.
Obr. 10.5 Srovnání izoentropická účinnosti upraveného a neupraveného oběžného kola [25]
10.4 Diskuse vlivu bezlopatkového difuzoru Dalším důležitým faktorem, který výrazným způsobem ovlivňuje charakteristiku kompresoru je použití lopatkového nebo bezlopatkového difuzoru. Stabilní oblast kompresoru s lopatkovým difuzorem je užší než s kompresorem s difuzorem lopatkovým. Příčinou je odtrhávání proudění od lopatek difuzoru v případě, že je vstupní úhel proudu významně odlišný od úhlu návrhového. [9] Ze srovnání charakteristik na obr. 10.6 vyplývá, že při použití lopatkového difuzoru dochází v tomto případě k většímu stlačení ve stupni pro všechny rychlosti v místech nejvyšší účinnosti. Větší je i dosahovaná účinnost kompresoru a to přibližně o 6 %. Při použití bezlopatkového difuzoru se ovšem stabilní oblast rozšiřuje nejen směrem k větším hmotnostním průtokům v celé oblasti charakteristiky, ale také dochází k posunu meze pumpování k nižším hmotnostním průtokům. Za povšimnutí také stojí srovnání účinnostních ostrovů, kde při použití lopatkového difuzoru dochází k posunu oblasti s nejvyšší účinností k nižším hmotnostním průtokům a nižšímu poměrnému stlačení.
- 41 -
Obr. 10.6 Porovnání pracovních oblastí kompresoru s lopatkovým difuzorem (modře) a bezlopatkovým difuzorem (černě) [9]
10.5 Diskuse vlivu natáčivých vstupních a difuzorových lopatek Aktivním prvkem rozšiřujícím stabilní oblast radiálního kompresoru jsou natáčivé vstupní lopatky a natáčivé difuzorové lopatky. V obou případech se mohou lopatky za běhu kompresoru natáčet a ovlivňovat tak rychlostní trojúhelníky proudění při různých pracovních režimech. Natáčivé vstupní lopatky (VIGV) se umisťují před oběžné kolo kompresoru a ovlivňují tak proudění vstupující do záběrníku. Lopatky nastavují vstupní rychlost ܿଵ, čímž se změní i relativní rychlost ݓଵ, tak aby došlo k minimalizaci úhlu náběhu a tím k ustálení proudění. Nevýhodou je, že při odklonu rychlosti ܿଵ od osového směru ve smyslu otáčení kola dojde ke snížení práce přenesené kolem. Příklad charakteristiky kompresoru s různým natočením vstupních lopatek je na obr. 10.7. Na obdobném principu pracují i natáčivé difuzorové lopatky, které se svým natočením snaží minimalizovat úhel náběhu a zabránit tak odtrhávání proudění od profilu lopatky. Příklad charakteristiky kompresoru s různými úhly natočení difuzorových lopatek je na obr. 10.7. Vzhledem ke skutečnosti, že je to právě difuzor, který má značný vliv na stabilní rozsah práce kompresoru, lze natáčením lopatek dosáhnout posunu stabilní oblasti ve velkém rozsahu hmotnostního průtoku. [9]
- 42 -
Obr. 10.7 Charakteristiky kompresorů s natáčivými vstupními a difuzorovými lopatkami [9]
- 43 -
11.
ÚVODNÍ NÁVRH GEOMETRIE KOMPRESORU
Základem pro správný návrhový výpočet vlastností kompresoru je správně definovaná a zpracovaná geometrie. K tomuto účelu byl použit program Bladegen systému ANSYS.
11.1 Základní rozměry oběžného kola Pro zpracování návrhového výpočtu kompresorové charakteristiky byly zadány základní rozměry oběžného kola, které jsou uvedeny v tab. 11.1. Výstupní úhel lopatek ߚଶǡ je zvolen 90 °, tedy radiální lopatky. Pro oběžné kolo je dále uvažováno zjednodušení zanedbáním mezery mezi kolem a skříní. Průměr náboje Vnější průměr Průměr kola Axiální délka Výstupní šířka kola Počet lopatek
ܦଵ
[mm]
ܦଵ௦
[mm]
ܦଶ
[mm]
ܮଵ
[mm]
ܾଶ
[mm]
ݖ
[-]
Tab. 11.1 Sumář základních rozměrů oběžného kola
Obr 11.1 Základní rozměry oběžného kola
- 44 -
ͻǡͷ
͵ͷǡ Ͷ
ͳͶǡͻ ͵ǡͻ ͻ
11.2 Vytvoření modelu oběžného kola v BladeGen Základní rozměry z tab. 11.1 byly použity při vytváření počítačového modelu oběžného kola kompresoru. K vytvoření modelu byla zvolena komponenta BladeGen programu ANSYS. Jedná se o specializovaný program pro tvorbu 3D modelů lopatkových strojů jako jsou kompresory, turbíny, ventilátory apod., a to jak axiálních, tak i radiálních oběžných kol. BladeGen je zároveň velice dobře propojený s dalšími komponentami systému ANSYS, což umožňuje použití vytvořeného modelu v dalších výpočtech například pevnostní analýzy nebo při výpočtech proudění ve stupni lopatkového stroje. Výstupem jsou soubory typu *.bgd, nicméně program umožňuje export geometrie i do jiných formátů jako je Pro/Engineer a další CAD systémy.
Obr. 11.1 Volba základních rozměrů oběžného kola Mimo základní rozměry je třeba stanovit také rozložení úhlů po délce lopatky. Oběžné kolo je rozděleno na pět vrstev, které kopírují zakřivený tvar kanálu a vedou od vstupní po výstupní hranu kola (obr. 11.2).
Obr. 11.2 Rozdělení oběžného kola na vrstvy - 45 -
Na každé této vrstvě je definováno zakřivení lopatky úhlem ߚ. Je třeba upozornit na to, že v komponentě BladeGen je úhel ߚ určen americkým způsobem, jako úhel mezi relativní rychlostí a osou kola (obr. 11.3).
Obr. 11.3 Definice úhlového systému v programu BladeGen
Pro radiální lopatky je proto výstupní úhel ߚ stanoven na 0 °. Dále je zvolen úhel ߚ na vstupu do oběžného kola na vnějším průměru ܦଵ௦ dle [10] na 59 °. Zbylý tvar lopatky je vytvořen jako plynulý přechod mezi úhly na vstupní a výstupní hraně. Rozložení úhlů v každé vrstvě je na obr. 11.4, kde horizontální osa označená % M (LE to TE), je vzdálenost na vrstvě mezi vstupní a výstupní hranou vyjádřená v procentech.
Obr. 11.4 Rozložení ߚ úhlů lopatky na vrstvách modelu oběžného kola Další požadovanou hodnotou je tloušťka lopatky. Protože pevnostní analýza přesahuje rámec této práce, je volena tloušťka konstantní po celé lopatce 1 mm s ohledem na [26], kde je stanovena minimální tloušťka pro malé kompresory na 0,35 mm. Rozložení tloušťky ve všech vrstvách oběžného kola je na obr. 11.5.
- 46 -
Obr. 11.5 Rozložení tloušťky lopatky na vrstvách modelu oběžného kola Výsledný 3D model oběžného kola radiálního kompresoru s radiálními lopatkami vytvořený v programu BladeGen je na obr. 11.6.
Obr. 11.6 Model oběžného kola radiálního kompresoru
- 47 -
Z rozměrů oběžného kola a jeho modelu je možné vypočítat, případně odměřit další potřebné parametry. Vypočítané parametry Vstupní průtočná plocha
ܦଵ௦ ଶ ܦଵ ଶ ܣଵ ൌ ߨ ή ቆ൬ ൰ െ ൬ ൰ ቇ ʹ ʹ ܣଵ ൌ ߨ ή
ǡଷହ ଶ ൬ቀ ቁ ଶ
െ
ǡଽହ ଶ ቀ ቁ ൰ ଶ
(11.1) ൌ ͻǡʹͶͷ ή ͳͲିସ ሾ݉ଶ ሿ
Výstupní průtočná plocha ܣଶ ൌ ߨ ή ܦଶ ή ܾଶ
(11.2)
ܣଶ ൌ ߨ ή ͲǡͲͶ ή ͲǡͲͲ͵ͻ ൌ ͷǡ͵ሾ݉ଶ ሿ
Hydraulický průměr
dle [30] se hydraulický průměr určí ze vztahu 11.3. ܦ௬ௗ
ܦଵ௦ െ ܦଵ ή ܾଵ ߨ ή ܦଶ ή ܾଶ ʹ ൌ ݖή ܾଵ ߨ ή ܦଵ ݖή ܾଶ ߨ ή ܦଶ
ܦ௬ௗ ൌ
ߨή ߨή
ͲǡͲ͵ͷ െ ͲǡͲͲͻͷ ή ͲǡͲͳ͵ͳ ߨ ή ͲǡͲͶ ή ͲǡͲͲ͵ͻ ʹ ൌ ͺǡͲͶͺ ή ͳͲିଷ ሾ݉݉ଶ ሿ ͻ ή ͲǡͲͲ͵ͻ ߨ ή ͲǡͲͶ ͻ ή ͲǡͲͳ͵ͳ ߨ ή ͲǡͲʹ
(11.3)
TRIM Trim je hodnota běžně se používající při popisu kompresoru, dle [37] je definována vztahem 11.4. ܶ ݉݅ݎൌ ܶ ݉݅ݎൌ
ଶ ܦଵ௦ ܦଶଶ
(11.4)
ͲǡͲ͵ͷଶ ൌ ͷͻǡͺͻሾെሿ ͲǡͲͲͶଶ
Parametry odměřené z 3-D modelu v programu BladeGen Délka středního průtočného vlákna ܮ ൌ ʹͲǡͲͺ݉݉
Hrdlo oběžného kola
ܣଵ௧ ൌ ͲͶǡͶ݉݉ଶ
- 48 -
12.
1-D VÝPOČET STUPNĚ KOMPRESORU PO STŘEDNÍM PROUDOVÉM VLÁKNU
Proudění ve stupni radiálního kompresoru je výrazně prostorového charakteru a tedy velice komplikované. Kompletní a přesná analýza proudění v oběžném kole představuje mnoho problémů, pro jejichž řešení je vyžadováno specializovaných programů a výkonných počítačů. Následující výpočet si z výše uvedeného důvodu neklade ambice na stanovení skutečných poměrů na výstupu z oběžného kola kompresoru. Jeho význam spočívá v odhadu pracovní oblasti kompresoru v různých provozních režimech. Výpočet stupně po středním proudovém vláknu představuje jednodimenzionální způsob řešení, při kterém jsou užity vztahy termodynamiky ideálního plynu a energetické bilance. Dále jsou uvažovány střední výpočtové rychlosti stanovené z rovnice kontinuity a z příslušných kinematických vztahů v daném průřezu. Pro odhad účinnosti kompresoru je použit také jednoduchý ztrátový model. Výpočet je proveden pro oběžné kolo, jakožto hlavní část kompresoru, udávající trend v tvaru kompresorové charakteristiky. Výpočet je zpracován v programu Microsoft Excel pro pět různých rychlostí rotace oběžného kola. Maximální rychlost rotace je volena s ohledem na pevnostní omezení oběžného kola. Pro hliníkové oběžné kolo se maximální obvodová rychlost dle [26] pohybuje v rozmezí ͵ͲͲȂ ͷͲͲ݉ ή ି ݏଵ.
V následujících kapitolách je popsán postup výpočtu pro příklad provozního režimu kompresoru uvedeného v tab. 12.1. Otáčky oběžného kola Hmotnostní tok kompresorem Teplota okolí Atmosférický tlak
݊
ሾ݉݅݊ିଵ ሿ
݉ሶ
ሾ݇݃ ή ି ݏଵ ሿ
ܶଵ
ሾܭሿ
ଵ
ሾܲܽሿ
Tab. 12.1 Provozní podmínky kompresoru pro 1-D výpočet
- 49 -
130000 0,11 288 101325
12.1 Stav proudu vzduchu na vstupu do oběžného kola Fyzikální vlastnosti vzduchu Poissonova konstanta ߢ ൌ ͳǡͶሾെሿ
Měrná plynová konstanta
ܴ ൌ ʹͺǡͲͷሾ ܬή ݇݃ିଵ ή ି ܭଵ ሿ
Měrná tepelná kapacita
ܿ ൌ ͳͲͲͷሾ ܬή ݇݃ିଵ ή ି ܭଵ ሿ
Vstupní rychlostní trojúhelník na středním proudovém vláknu Při výpočtu vstupního rychlostního trojúhelníku proudu vzduchu lze uvažovat, že obvodová složka celkové rychlosti je nulová [9]. Proto platí vztah 12.1 ܥଵ ൌ ܥଵ
(12.1)
Obvodová rychlost na vstupu ଶ ଶ ܦଵ௦ ܦଵ ݊ ݑଵ ൌ ߨ ή ඨ ή ʹ Ͳ
(12.2)
ͲǡͲ͵ͷଶ ͲǡͲͲͻͷଶ ͳ͵ͲͲͲͲ ݑଵ ൌ ߨ ή ඨ ή ൌ ͳǡ͵Ͷሾ݉ ή ି ݏଵ ሿ ʹ Ͳ
Z rovnice kontinuity se stanoví celková rychlost ܿଵ ൌ
݉ሶ Ͳǡͳͳ ൌ ൌ ͻǡͲͺሾ݉ ή ି ݏଵ ሿ ߩଵ ή ܣଵ ͳǡʹʹͷ ή ͻǡʹͷ ή ͳͲିସ
kde
ߩଵ ൌ
ଵ ͳͲͳ͵ʹͷ ൌ ൌ ͳǡʹʹͷሾ݇݃ ή ݉ିଷ ሿ ܴ ή ܶଵ ʹͺǡͲͷ ή ʹͺͺ
(12.3)
(12.4)
Relativní rychlost na vstupu
ݓଵ ൌ ටܿଵଶ ݑଵଶ ൌ ඥͻǡͲʹଶ ͳǡ͵Ͷଶ ൌ ʹͲʹǡͳͺሾ݉ ή ି ݏଵ ሿ
Obr. 12.1 Vstupní rychlostní trojúhelník na středním proudovém vláknu - 50 -
(12.5)
Vstupní rychlostní trojúhelník na vnějším průměru oběžného kola Obvodová rychlost na vstupu
ݑଵ௦ ൌ ߨ ή ܦଵ௦ ή
݊ ͳ͵ͲͲͲͲ ൌ ߨ ή ͲǡͲ͵ͷ ή ൌ ʹͶʹǡʹͲሾ݉ ή ି ݏଵ ሿ Ͳ Ͳ
(12.6)
Relativní rychlost na vstupu
ଶ ݓଵ௦ ൌ ටܿଵଶ ݑଵ௦ ൌ ඥͻǡͲʹଶ ʹͶʹǡʹͲଶ ൌ ʹͲǡͻ͵ሾ݉ ή ି ݏଵ ሿ
(12.7)
Vstupní rychlostní trojúhelník na patním průměru lopatky Obvodová rychlost na vstupu ݑଵ ൌ ߨ ή ܦଵ ή
݊ ͳ͵ͲͲͲͲ ൌ ߨ ή ͲǡͲͲͻͷ ή ൌ Ͷǡሾ݉ ή ି ݏଵ ሿ Ͳ Ͳ
(12.8)
Relativní rychlost na vstupu
ଶ ݓଵ ൌ ටܿଵଶ ݑଵ ൌ ඥͻǡͲʹଶ Ͷǡଶ ൌ ͳͳǡͶͶሾ݉ ή ି ݏଵ ሿ
(12.9)
Termodynamické vlastnosti proudu na vstupu Rychlost zvuku na vstupu ܽଵ ൌ ඥߢ ή ܴ ή ܶଵ ൌ ඥͳǡͶ ή ʹͺǡͲͷ ή ʹͺͺ ൌ ͵ͶͲǡʹͳሾ݉ ή ି ݏଵ ሿ
(12.10)
Machovo číslo pro absolutní rychlost ܽܯଵ ൌ
ܿଵ ͻǡͲʹ ൌ ൌ ͲǡʹͺͷሾȂ ሿ ܽଵ ͵ͶͲǡʹͳ
(12.11)
Machovo číslo pro relativní rychlost ܽܯଵ ൌ
ݓଵ ʹͲʹǡͳͺ ൌ ൌ ͲǡͷͻͶሾȂ ሿ ܽଵ ͵ͶͲǡʹͳ
(12.11)
Celková teplota na vstupu ܶଵ ൌ ܶଵ ή ൬ͳ
ߢെͳ ͳǡͶ െ ͳ ή ܽܯଵଶ ൰ ൌ ʹͺͺ ή ൬ͳ ή Ͳǡʹͺͷଶ ൰ ൌ ʹͻʹǡͻሾܭሿ ʹ ʹ
Celkový tlak na vstupu ଵ ଵ
(12.12)
ିଵ ߢെͳ ൌ ଵ ή ൬ͳ ή ܽܯଵଶ ൰ ʹ
ଵǡସ
ଵǡସିଵ ͳǡͶ െ ͳ ൌ ͳͲͳ͵ʹͷ ή ൬ͳ ή Ͳǡʹͺͷଶ ൰ ൌ ͳͲʹͳͺǡͻͺሾܲܽሿ ʹ
- 51 -
(12.13)
Hustota vzduchu na vstupu z celkových stavů ߩଵ ൌ
ଵ ͳͲʹͳͺǡͻͺ ൌ ൌ ͳǡʹሾ݇݃ ή ݉ିଷ ሿ ܴ ή ܶଵ ʹͺǡͲͷ ή ʹͻʹǡͻ
(12.14)
Statická entalpie vzduchu na vstupu
݅ଵ ൌ ܶଵ ή ܿ ൌ ʹͺͺ ή ͳͲͲͷ ൌ ʹͺͻͶͶͲሾ ܬή ݇݃ିଵ ሿ
(12.15)
12.2 Stav proudu vzduchu na výstupu z oběžného kola Ideální výstupní rychlostní trojúhelník Při výpočtu výstupního rychlostního trojúhelníku pro kompresor s radiálními lopatkami lze dle [8] a [19] uvažovat, že se výstupní radiální složka celkové rychlosti rovná celkové vstupní rychlosti. ܿଶǡ ൌ ܿଵ
ܿଶǡ ൌ ͻǡͲͺሾ݉ ή ି ݏଵ ሿ
(12.16)
݊ ͳ͵ͲͲͲͲ ൌ ߨ ή ͲǡͲͶ ή ൌ ͵ͳ͵ǡͳͳሾ݉ ή ି ݏଵ ሿ Ͳ Ͳ
Obvodová rychlost na výstupu ݑଶǡ ൌ ߨ ή ܦଶ ή
(12.17)
Pro radiální lopatky (ߚଶ ൌ ͻͲι) se výstupní relativní rychlost rovná radiální složce celkové rychlosti. ݓଶǡ ൌ ܿଶǡ
ݓଶǡ ൌ ͻǡͲͺሾ݉ ή ି ݏଵ ሿ
(12.18)
Ideální výstupní celková rychlost
ଶ ଶ ܿଶǡ ൌ ටݓଶǡ ݑଶǡ ൌ ඥͻǡͲͺଶ ͵ͳ͵ǡͳͳଶ ൌ ͵ʹǡͺʹሾ݉ ή ି ݏଵ ሿ
(12.19)
Výstupní rychlostní trojúhelník ovlivněný skluzem V radiálním kompresorovém stupni dochází vlivem rotace oběžného kola k vytvoření protisměrného víru v mezilopatkovém kanále, v důsledku kterého dochází k odklonu výstupní relativní rychlosti proti smyslu otáčení oběžného kola. Dochází tak ke zmenšení obvodové složky celkové rychlosti ܿଶ௨ (Obr. 12.2). Tento jev se popisuje součinitelem skluzu Ɋ, který je pro radiální lopatky dle [9] vyjádřen vztahem 12.19. Ɋൌ
ܿଶ௨ Ͳǡ͵ ή ߨ ൌͳെ ሾെሿ ݑଶ ݖ
Ɋൌͳെ
Ͳǡ͵ ή ߨ ൌ Ͳǡͺሾെሿ ͻ
(12.20)
Obvodová výstupní rychlost ݑଶ ൌ ݑଶǡ
ݑଶ ൌ ͵ͳ͵ǡͳͳሾ݉ ή ି ݏଵ ሿ
(12.21)
- 52 -
Obvodová složka výstupní celkové rychlosti
ܿଶ௨ ൌ Ɋ ή ݑଶ ൌ Ͳǡͺ ή ͵ͳ͵ǡͳͳ ൌ ʹͶͶǡʹሾ݉ ή ି ݏଵ ሿ
(12.22)
Radiální složky výstupní celkové rychlosti ܿଶ ൌ ܿଶǡ
ܿଶ ൌ ͻǡͲͺሾ݉ ή ି ݏଵ ሿ
(12.23)
Celková rychlost na výstupu
ଶ ଶ ܿଶ ൌ ටܿଶ௨ ܿଶ ൌ ඥʹͶͶǡʹଶ ͻǡͲͺଶ ൌ ʹʹǡͺͶሾ݉ ή ି ݏଵ ሿ
(12.24)
Relativní rychlost na výstupu z obr. 12.2 ଶ ଶ ଶ ݓଶ ൌ ටݓଶ௨ ݓଶ ൌ ටሺݑଶ െ ܿଶ௨ ሻଶ ܿଶ
ݓଶ ൌ ඥሺ͵ͳ͵ǡͳͳ െ ʹͶͶǡʹሻଶ ͻǡͲͺଶ ൌ ͳͳͻǡͲʹሾ݉ ή ି ݏଵ ሿ
(12.25)
Úhel celkové výstupní rychlosti k radiálnímu směru ߙଶ ൌ
ିଵ
ܿଶ ͻǡͲͺ ൌ
ିଵ ൌ ͺǡ͵ʹሾιሿ ܿଶ ʹʹǡͺͶ
(12.26)
Obr. 12.2 Ideální výstupní rychlostní trojúhelník a skutečný rychlostní trojúhelník ovlivněný skluzem Termodynamické vlastnosti proudu na výstupu Statická výstupní entalpie vyjádřená z energetické bilance na rotorové lopatkové řadě
݅ଵ
ݓଵଶ ݑଵଶ ݓଶଶ ݑଶଶ െ ൌ ݅ଶ െ ʹ ʹ ʹ ʹ
݅ଶ ൌ ݅ଵ
ݓଵଶ ݑଵଶ ݓଶଶ ݑଶଶ െ െ ʹ ʹ ʹ ʹ
݅ଶ ൌ ʹͺͻͶͶͲ
(12.27)
ʹͲʹǡͳͺଶ ͳǡ͵Ͷଶ ͳͳͻǡͲʹଶ ͵ͳ͵ǡͳͳଶ െ െ ൌ ͵͵ͲͺͺǡͻͶሾ ܬή ݇݃ିଵ ሿ ʹ ʹ ʹ ʹ
Výstupní statická teplota ܶଶ ൌ
݅ଶ ͵͵ͲͺͺǡͻͶ ൌ ൌ ͵͵ͶǡͶʹሾܭሿ ܿ ͳͲͲͷ
(12.28)
- 53 -
Rychlost zvuku na výstupu
ܽଶ ൌ ඥߢ ή ܴ ή ܶଶ ൌ ඥͳǡͶ ή ʹͺǡͲͷ ή ͵͵ͶǡͶʹ ൌ ͵ǡͲሾ݉ ή ି ݏଵ ሿ ܿଶ ʹʹǡͺͶ ൌ ൌ ͲǡͳͲሾെሿ ܽଶ ͵ǡͲ
(12.29)
Machovo číslo na výstupu pro absolutní rychlost ܽܯଶ ൌ
(12.30)
ݓଶ ͳͳͻǡͲʹ ൌ ൌ Ͳǡ͵ʹͶሾെሿ ܽଶ ͵ǡͲ
Machovo číslo na výstupu pro relativní rychlost ܽܯଶ ൌ
(12.31)
݉ሶ Ͳǡͳͳ ൌ ൌ ʹǡͲͳͲͶሾ݇݃ ή ݉ିଷ ሿ ܣଶ ή ܿଶ ͲǡͲͲͲͷ͵ ή ͻǡͲͺ
Hustota vzduchu na výstupu ߩଶ ൌ
ଶ ൌ ߩଶ ή ܴ ή ܶଶ ൌ ʹǡͲͳͲͶ ή ʹͺǡͲͷ ή ͵͵ͶǡͶʹ ൌ ͳͻʹͻͻͷǡ͵ͺሾܲܽሿ
(12.32)
Statický tlak na výstupu
(12.33)
Pro určení celkového výstupního tlaku ଶ je zapotřebí stanovit izoentropickou účinnost oběžného kola. Této problematice je věnována následující kapitola.
12.3 Návrh ztrátového modelu pro stanovení izoentropické účinnosti V následující části jsou uvedeny vztahy pro výpočet změny entalpie vlivem ztrát v oběžném kole. Pokud není uvedeno jinak, je vycházeno z literatury [21] a [33]. Ztráty na vstupu do oběžného kola
Ztráty na vstupu do oběžného kola vznikají, pokud je úhel vstupní relativní rychlosti ߚଵ odlišný od vstupního úhlu lopatky ߚଵǡ (obr. 12.3). Vstupní ztráta určuje množství energie, která je spotřebována na přizpůsobení se směru proudění vstupnímu tvaru (úhlu) lopatky.
Obr. 12.3 Rozdíl vstupního úhlu proudu a vstupního úhlu lopatky [9] - 54 -
ଶ ݓ௨ ଶ ൌ Ͳǡͷ ή ൫ݑଵ െ ܿଵ௨ െ ܿ݃ݐሺߚଵǡ ሻ ή ܿଵ ൯ ʹ ൌ Ͳǡͷ ή ሺͳǡ͵Ͷ െ Ͳ െ ܿ݃ݐሺͶͻιሻ ή ͻǡͲͺሻଶ ൌ Ͷ͵ʹͲǡͳሾ ܬή ݇݃ିଵ ሿ
ο݅ ൌ ο݅
(12.34)
Ztráta při průtoku oběžným kolem
Ztráta při průtoku oběžným kolem zahrnuje ztráty třením v mezních vrstvách, sekundárním prouděním a vířením při odtržení proudění [29]. Pro stanovení množství ztracené energie platí vztah 12.35. ο݅ௗ ൌ ͲǡͲͷ ή ܦଶ ή ݑଶଶ ൌ ͲǡͲͷ ή ͲǡͷͶଶ ή ͵ͳ͵ǡͳͳଶ ൌ ʹͳͲǡ͵ሾ ܬή ݇݃ିଵ ሿ
kde
ܦ ൌ ͳ െ
(12.35)
Ͳǡͷ ή Ɋ ή ݑଶଶ ݑଶଶ
ݓଶ ݓଵ௦ ݓଵ௦ ݖ ܦଵ௦ ʹ ή ܦଵ௦ ݓଶ ή ቆቀߨቁ ή ቀͳ െ ܦଶ ቁ ܦଶ ቇ
Ͳǡͷ ή Ͳǡͺ ή ͵ͳ͵ǡͳͳ ͳͳͻǡͲʹ ͵ͳ͵ǡͳͳଶ ܦ ൌ ͳ െ ʹͲǡͻ͵ ʹͲǡͻ͵ ͻ ͲǡͲ͵ͷ ʹ ή ͲǡͲ͵ͷ ͳͳͻǡͲʹ ή ቆቀߨቁ ή ቀͳ െ ͲǡͲͶ ቁ ͲǡͲͶ ቇ ଶ
(12.36)
ܦ ൌ ͲǡͷͶሾെሿ
Třecí ztráty na povrchu mezilopatkového kanálu Tyto ztráty vyjadřují míru ztracené energie způsobenou adhezními silami mezi povrchem kanálu, tvořeného povrchem lopatek, povrchem náboje a povrchem skříně, jež je v kontaktu se vzduchem. Ztracená energie je vyjádřena vztahem 12.37. ο݅௦
ܥ ή ܮ ή ݉ሶଶ ൌ ʹ ή ܦ௬ௗ ή ߩଵଶ ή ܣଵଶ ή ଶ ߚଵ
ο݅௦ ൌ
ǡͲͻ ή ͳͲିଷ ή ͲǡͲʹͲͶ ή Ͳǡͳͳଶ ൌ ͳ͵ͺǡͶ͵ሾ ܬή ݇݃ିଵ ሿ ʹ ή ͺǡͲͶͺ ή ͳͲିଷ ή ͳǡʹʹͷଶ ή ͻǡʹͶͷଶ ή ͳͲି଼ ή ଶ Ͷͻι
(12.37)
Ztrátový koeficient povrchového tření
ܥ ൌ Ͷ ή ݂ ൌ Ͷ ή ͳǡͷ ή ͳͲିଷ ൌ ǡͲͻ ή ͳͲିଷ ሾെሿ
(12.38)
kde
݂ ൌ Ͳǡ͵ͳͶ ή ܴ݁ଶିǡଶହ ൌ Ͳǡ͵ͳͶ ή ͳ͵ͶͶͻͻͶͺʹିǡଶହ ൌ ͳǡͷ ή ͳͲିଷ ሾെሿ
(12.39)
Reynoldsovo číslo na výstupu je pro lopatkové stroje určeno vztahem 12.40 ܴ݁ଶ ൌ
ݑଶ ή ݎଶ ͵ͳ͵ǡͳͳ ή ͲǡͲʹ͵ ൌ ൌ ͳ͵ͶͶͻͻͶͺʹሾെሿ ߥଶ ͳǡͲͳʹ ή ͳͲିହ
- 55 -
(12.40)
Dynamickou viskozitu lze určit pro teplotu ܶଶ ze Sutherlandova vztahu 12.41. [32] ܶ ܥ ܶଶ ߟଶ ൌ ߤ ή ൬ ൰ήቆ ቇ ܶଶ ܥ ܶ
ଵǡହ
ʹͻͳǡͳͷ ͳʹͲ ͵͵ͶǡͶʹ ଵǡହ ߟଶ ൌ ͳǡͺʹ ή ͳͲିହ ή ൬ ൰ή൬ ൰ ൌ ʹǡͲ͵ͷ ή ͳͲିହ ሾܲܽ ή ݏሿ ͵͵ͶǡͶʹ ͳʹͲ ʹͻͳǡͳͷ
kde
(12.41)
ߟ - dynamická viskozita při ܶ (ߟ ൌ ͳǡͺʹ ή ͳͲିହ ሾܲܽ ή ݏሿሻ ܶ – referenční teplota (ܶ ൌ ʹͻͳǡͳͷሾܭሿሻ
– ܥSutherlandova konstanta pro vzduch ( ܥൌ ͳʹͲሾܭሿ)
Kinematická viskozita
ߥଶ ൌ
ߟଶ ʹǡͲ͵ͷ ή ͳͲିହ ൌ ൌ ͳǡͲͳʹ ή ͳͲିହ ሾ݉ଶ ή ି ݏଵ ሿ ʹǡͲͳͲͶ ߩଶ
(12.41)
Směšovací ztráty
Jsou způsobené nevyrovnaným rychlostním a tlakovým polem těsně za výstupem z oběžného kola [29]. Výpočet této ztráty vychází z teorie dvouzónového proudění v oběžném kole (obr. 12.4), které je blíže popsáno v [20]. Množství ztracené energie je vyjádřeno vztahem 12.42. ο݄௫ ൌ ο݄௫ ൌ
ͳ ͳ െ ߝ െ ܾ ܿ כଶଶ ή ൬ ൰ή ͳ ଶ ߙଶ ͳെߝ ʹ
ͳ ͳ െ Ͳǡͷ͵ʹ െ ͳ ʹʹǡͺͶଶ ή ൬ ൰ή ൌ Ͳͺǡ͵ሾ ܬή ݇݃ିଵ ሿ ͳ ଶ ͺǡ͵ʹι ͳ െ Ͳǡͷ͵ʹ ʹ
(12.42)
Ve vztahu 12.42 je ܾ כpoměr výstupní plochy kola a vstupní plochy difuzoru, z důvodu neuvažování difuzoru je voleno ܾ כൌ ͳሾെሿ
Obr. 12.4 Schéma dvouzónového proudění na výstupu z oběžného kola [21] - 56 -
Poměr sekundárního průtočného průřezu k celkovému průtočnému průřezu ߝ je dle [36] určen vztahem 12.43. ߝൌ ߝൌ
െͲǡͲ ඥͲǡͲଶ Ͷ ή Ͳǡͻ͵ ή ߯ ʹ ή Ͳǡͻ͵
(12.43)
െͲǡͲ ඥͲǡͲଶ Ͷ ή Ͳǡͻ͵ ή Ͳǡ͵ ൌ Ͳǡͷ͵ʹሾെሿ ʹ ή Ͳǡͻ͵
kde ߯je dle [20] pro malé kompresory ߯ ൌ Ͳǡ͵ሾെሿ
Izoentropická účinnost oběžného kola
Izoentropickou účinnost radiálního kompresoru lze vypočítat dle [19] ze vztahu 12.44.
ߟ௭ ൌ
ȟ୧ୢ െ ȟߟ െ ȟߟ െ ȟߟ௩ ȟ୧ୢ ȟ୪୭ୱୱ
(12.44)
Pro oběžné kolo s radiálními lopatkami bez uvažování mezery mezi kolem a skříní přechází vztah do tvaru 12.45. ߟ௭ ൌ ߟ௭ ൌ
Ɋ ή ݑଶଶ െ ȟߟ Ɋ ή ݑଶଶ ȟ୪୭ୱୱ
Ͳǡͺ ή ͵ͳ͵ǡͳͳଶ െ ͲǡͲ͵ ൌ Ͳǡͺʹሾെሿ Ͳǡͺ ή ͵ͳ͵ǡͳͳଶ ͳʹͻǡʹ
(12.45)
Ve vztahu 12.45 představuje člen ȟ୪୭ୱୱ sumu uvažovaných ztrát v oběžném kole vypočtených v předchozích krocích. Suma ztrát je poté určena rovnicí 12.46. ȟ୪୭ୱୱ ൌ ȟ୧୬ୡ ȟୠ୪ୢ ȟୱ ȟ୫୧୶
ȟ୪୭ୱୱ ൌ Ͷ͵ʹͲǡͳ ʹͳͲǡ͵ ͳ͵ͺǡͶ͵ Ͳͺǡ͵ ൌ ͳʹͻǡʹሾ ܬή ݇݃ିଵ ሿ
(12.46)
Člen ȟߟ vyjadřuje pokles účinnosti vlivem zpětného proudění v oběžném kole a je dle [9] jako typickou hodnotu možné uvažovat 3 % pokles. ȟߟ ؐ ͲǡͲ͵ሾെሿ
12.4 Výsledná charakteristika kompresoru pro 1-D výpočet Výsledný tlakový poměr S využitím výsledků z předešlých kapitol lze stanovit velikost celkového tlaku za oběžným kolem ଶ . Dle [19] je tento tlak definován vztahem 12.47. ଶ ଶ
ߟ௭ ή Ɋ ή ݑଶଶ ିଵ ൌ ଵ ή ቆͳ ቇ ܶଵ ή ܿ
ଵǡସ
Ͳǡͺʹ ή Ͳǡͺ ή ͵ͳ͵ǡͳͳଶ ଵǡସିଵ ൌ ͳͲʹͳͺǡͻͺ ή ቆͳ ൌ ʹͳʹͲͺͷǡͺʹሾܲܽሿ ቇ ʹͻʹǡͻ ή ͳͲͲͷ
- 57 -
(12.47)
Poměr výstupního a vstupního celkového tlaku ȫ ൌ
ଶ ʹͳʹͲͺͷǡͺʹ ൌ ଵ ͳͲʹͳͺǡͻͺ
(12.48)
Limitace charakteristiky aerodynamickým ucpáním K aerodynamickému ucpání kompresoru dojde v případě, kdy je přikročen kritický hmotnostní tok, viz kap. 7.
݉ሶ௧ ݉ሶ௧
ۇ ൌ ܣଵ ή ߩଵ ή ܽଵ ή ۈ ۉ
ʹ ሺߛ െ ͳሻ ή ߛͳ
ାଵ ଶ ଶήሺିଵሻ ܷଵ ଶ ۊ ܽଵ
ۋ ی
ଵǡସାଵ
Ͷǡଶ ଶήሺଵǡସିଵሻ ʹ ሺͳǡͶ െ ͳሻ ή ͵Ͷʹǡͻଶ ൌ ǡͲͶͶ ή ͳͲିସ ή ͳǡʹ ή ͵Ͷʹǡͻ ή ൮ ൲ ͳǡͶ ͳ
(12.49)
݉ሶ௧ ൌ Ͳǡͳͺʹ͵ሾ݇݃ ή ି ݏଵ ሿ
݉ሶ ൏ ݉ሶ௧ ֜ provozní bod leží vlevo od kritického průtoku, takže nedochází k zahlcení kompresoru. Limitace charakteristiky mezí pumpování V souladu s kap. 8.3 je nestabilní část pracovní oblasti práce kompresoru oddělena od stabilní části mezí pumpování. Dělícím místem mezi stabilní a nestabilní oblastí je uvažováno místo nejvyššího dosaženého stlačení při daných otáčkách. Mez pumpování poté tvoří spojnici vrcholů křivek větví charakteristiky.
- 58 -
2,5
n = 150000 min-1
2,25
2
Πk = p02/p01 [-]
n = 130000 min-1
1,75
n = 110000 min-1
1,5
n = 90000 min-1 1,25
n = 70000 min-1 T1 = 288 [K] p1 = 0,1013 [MPa] 1 0
0,05
0,1
0,15
m [kg/s] Obr. 12.1 Kompresorová charakteristika pro 1-D výpočet
- 59 -
0,2
13.
2-D VÝPOČET STUPNĚ KOMPRESORU POMOCÍ VistaTF
Jak již bylo zmíněno v předchozí kapitole, přesného výpočtu proudění v kompresoru lze docílit pouze s využitím specializovaných programů. Vedoucí pozici mezi těmito programy zaujímá systém ANSYS, který má v této oblasti již dlouhou historii. Obsahuje programy pro 3-D výpočet proudění jako je například CFX, ale také pro zjednodušený případ 2-D výpočtu proudění v programu VistaTF. A právě pro program VistaTF ve spolupráci s firmou Honeywell zpracováno posouzení možnosti výpočtu kompresorové charakteristiky. VistaTF je program, který vyplňuje mezeru mezi jednoduchým 1-D výpočtem a komplexním 3-D výpočtem proudění. Jeho výhodou je především rychlost a nenáročnost výpočtů, které lze realizovat na běžných počítačích, bez potřeby spouštění výpočtů v počítačových clusterech. V programu VistaTF lze analyzovat proudění v různých typech lopatkových strojů a to jak samostatný stupeň, tak i vícestupňové stroje. Pro účely této práce je rozhodující možnost výpočtu proudění v radiálním kompresoru. Návrhový výpočet kompresorové charakteristiky v programu VistaTF obnáší spolupráci s několika dalšími programy systému ANSYS. Geometrii kompresoru, tak jak byla vytvořena v programu BladeGen (kap. 11), nelze přímo převést do VistaTF. Model je nutné nejprve importovat z BladeGen do programu DesignModeler. Zde je již poté možno vytvořit geometrii kompatibilní s VistaTF. VistaTF poté po nastavení počátečních podmínek provede výpočet pro danou geometrii. Po provedení výpočtu se výsledek zobrazí v programu CFD-Post, který slouží pro vyhodnocování výsledků. Metodika řešení je schematicky zobrazena na obr. 13.1.
BladeGen
DesignModeler
VistaTF
CFD-Post
Obr. 13.1 Schéma řešení výpočtu proudění ve VistaTF
13.1 Příprava modelu části kola v programu DesignModeler Mezikrok mezi převodem geometrie z BladeGen do VistaTF je import geometrie do komponenty DesignModeler. K importu slouží funkce Load BGD, nacházející se na panelu nástrojů BladeEditor. Po vybrání souboru typu *.bgd vytvořeného v BladeGen, se načtou všechna potřebná data pro tvorbu geometrie. Model lopatkování oběžného kola se následně vytvoří volbou funkce Generate. Model vygenerovaný v prostředí DesignModeler je na obr. 13.2.
- 60 -
Obr. 13.2 3-D Model lopatkování radiálního oběžného kola Velice užitečnou funkcí programu DesignModeler je možnost výpočtu umístění a plochy hrdla oběžného kola. K tomuto slouží funkce ThroatArea, jejíž výsledek je na obr. 7.1. Nicméně nejdůležitější funkcí je VistaTFExport. Po zvolení této funkce a vybrání geometrie lopatkování, je umožněno pomocí funkce Generate vytvořit data potřebná pro správný chod programu VistaTF.
13.2 Sumář nastavení výpočtu ve VistaTF Počáteční nastavení výpočtu je dostupné z prostředí Workbench systému ANSYS. Po výběru komponenty Vista TF je možné provést úpravu nastavení v položce Setup. Souhrn nabízených možností nastavení je v tab. 13.1. Počet proudových vláken Maximum iterací
Smysl rotace stroje Rychlost otáčení Volba proudění Hmotnostní tok Vstupní celkový tlak Vstupní celková teplota
Nastavení počtu proudových vláken použitých při výpočtu. Je možné volit mezi 9 a 27. Pro výpočet je doporučená hodnota 17 vláken Definuje maximální počet iteračních smyček při výpočtu. Pro výpočet jednoduchého stupně radiálního kompresoru je doporučená hodnota 500 iterací Určení směru rotace oběžného kola Nastavení počtu otáček kola za jednotku času Volba okrajových podmínek výpočtu. Lze vybrat hmotnostní tok, nebo hmotnostní tok zároveň se stlačením (ଶ Ȁଵ), případně hmotnostní tok s rozdílem tlaků (ଶ െ ଵ ) Nastavení hmotnosti toku tekutiny na vstupu do oběžného kola Zadání celkového tlaku na vstupu do oběžného kola, který je možné spočítat z rovnice 12.13 Zadání celkové teploty na vstupu do oběžného kola, kterou je možné spočítat z rovnice 12.12 - 61 -
Vstupní úhel Referenční průměr Polytropická účinnost Volba tekutiny Měrná tepelná kapacita
Nastavení úhlu celkové rychlosti, pod kterým vstupuje tekutina do oběžného kola. Úhel se zadává od osového směru, přičemž kladná hodnota znamená po směru rotace kola Hodnota referenčního průměru je dle [31] na výstupu oběžného kola, tedy průměr ܦଶ Nastavení polytropické účinnosti kompresoru Je možné nastavit výpočet pro ideální plyn, reálný plyn nebo kapalinu Nastavení měrné tepelné kapacity tekutiny
Izoentropický koeficient Zadání izoentropického koeficientu Dynamická viskozita Počáteční ܿ Τݑ
Zadání dynamické viskozity. Pokud je zadána 0, program spočítá dynamickou viskozitu podle vlastních vztahů, přičemž je tato možnost nabízena pouze pro ideální plyn počáteční odhad poměru osové rychlosti ܿ a charakteristické rychlosti ݑ , který je dle programové dokumentace vyjádřen vztahem: ܿ ܿ ൌ ݑ ܦ ʹ ή݊
Tab. 13.1 Výčet základního nastavení výpočtu ve VistaTF
13.3 Editace konfiguračních souborů Mimo základní nastavení přístupné z prostředí Workbench je možné upravit nastavení výpočtu i editací konfiguračních souborů. Konfigurační soubory jsou celkem tři. Zvlášť upravují nastavení pro kontrolní data (*.cont), aerodynamická data (*.aert) a korelační data (*.cort). Popis parametrů, jež lze v těchto souborech měnit je v [34]. Zajímavou volbu nastavení představují zejména parametry i_expert v konfiguračním souboru *.cont a i_loss v souboru *.cort. i_expert – Změnou této hodnoty z 0 na 1 umožní programu využít všechny modely pro určení aerodynamického ucpání oběžného kola. Při nastavení této hodnoty je vyžadován přesný výpočet plochy hrdla, což je popsáno v kapitole 13.1. Program v tomto případě při zadání většího hmotnostního průtoku, než je kolo schopné zpracovat, nekonverguje. i_loss – Nabízejí se tři možnosti této hodnoty: · i_loss = 0 – · i_loss = 1 – · i_loss = 2 -
V tomto případě nejsou při výpočtu uvažovány ztráty v oběžném kole. Ztráty v kole jsou definovány pomocí polytropické účinnosti zadané v základním nastavení. K výpočtu ztrát jsou použity modely ztrát obsažené v programu, přičemž je možné tyto jednotlivé modely dále modifikovat.
Pro výpočet charakteristiky kompresoru je vhodné použití volby i_expert = 1, proto bude dále použita při všech výpočtech. Nastavení proměnné i_loss na hodnotu 2 se zdá být také nejlepší volbou pro výpočet charakteristiky. V [34] se ovšem uvádí, že je pro získání adekvátních - 62 -
výsledků žádoucí upravit ztrátové modely pomocí multiplikačních koeficientů, které lze také nastavit v konfiguračním souboru *.cort. Dokumentace k úpravě těchto koeficientů není bohužel dostupná a jedinou možností získání potřebných dat je experimentální měření. Z tohoto důvodu je pro další výpočet zvolena hodnota i_loss = 1.
13.5 Výpočet otáčkové větve charakteristiky Při nastavování výpočtu programu VistaTF, lze vhodně využít výsledky získané z výpočtu po středním proudovém vláknu z kap. 12. Při volbě hmotnostního toku kompresorem při zadaných otáčkách lze s pomocí obr. 12.1 zvolit hodnotu ze stabilní části kompresorové charakteristiky. S využitím hodnot z výpočtu v kap. 12 je také zadána celková teplota a celkový tlak vzduchu na vstupu. Nastavení VistaTF je na obr. 13.3.
Obr. 13.3 Nastavení VistaTF Po provedení výpočtu se výsledek zobrazí v programu CFD-Post. Hodnoty, které se nemění po délce proudového vlákna, jako je např. poměrné stlačení nebo izoentropická účinnost oběžného kola, se vypíší do tabulky (obr. 13.4a). Hodnoty veličiny, které se po délce proudového vlákna mění, jsou zobrazeny pomocí barevných kontur na řezu mezilopatkového kanálu oběžného kola. Příkladem může být rozložení meridiální rychlosti na obr. 13.4b. Vypočtené hodnoty je nakonec možné zanést do kompresorové mapy.
- 63 -
Obr. 13.4 a) Tabulka výsledných hodnot b) Rozložení velikosti meridiální rychlosti proudění
13.6 Výpočet nestabilní oblasti Důležitým bodem pro vykreslení kompresorové charakteristiky je určení meze pumpování. Dle kap. 9 je důležitým vodítkem při stanovení nestabilního provozu kompresoru zpětné proudění na vnějším průměru záběrníku a s tím spojené vytvoření recirkulační oblasti před oběžným kolem. Tento závěr potvrzují i výsledky 3D výpočtu proudění, které byly poskytnuty firmou Honeywell. Na obr. 13.5 je zachycen režim proudění v blízkosti meze pumpování. Je zde patrná oblast se zvýšenou entropií před oběžným kolem, značící nestabilní proudění. Na obr 13.6 je rozložení rychlostí v prostorovém modelu mezilopatkového kanálu při stejném režimu proudění.
Obr. 13.5 Recirkulační oblast před oběžným kolem [35] - 64 -
Obr. 13.6 Rozložení rychlostí v kanálu oběžného kola [35] VistaTF dle [34] neumožňuje vzhledem k použitému matematickému modelu výpočet zpětného proudění, a proto nelze pozorovat v blízkosti meze pumpování tvoření zpětného proudění na vnějším průměru záběrníku, ani oblast zvýšené entropie před oběžným kolem. Nicméně i přes tuto obtíž je možné z vypočtených dat určit alespoň přibližnou polohu meze pumpování. V [22] je definován kritický úhel náběhu na středním proudovém vláknu, při jehož překročení lze předpokládat, že dochází k utržení proudění v lopatkové mříži. Přičemž při prvotním návrhu kompresoru je možné takový stav považovat za pumpování. Kritický úhel náběhu je funkcí Machova čísla vypočteného z relativní vstupní rychlosti na středním proudovém vláknu. ଶ ݅௧ ൌ ͶͶǡͷͷ͵ െ Ͷͻǡͷ ή ܽܯଵ ͳͶǡͳ ή ܽܯଵ
(13.1)
Rozhodovací kritérium vyjádřené vztahem 13.1 lze aplikovat na výstupní data, jelikož VistaTF umožňuje výpočet úhlu náběhu (obr. 13.8) a rozložení Machova čísla pro relativní rychlost na proudovém vláknu (obr. 13.9). Příklad výpočtu pro ݊ ൌ ͳ͵ͲͲͲͲሾ݉݅݊ିଵ ሿ.
nestabilní
režim
je proveden pro
݉ሶ ൌ ͲǡͲͺሾ݇݃ ή ି ݏଵ ሿ při
Hodnoty úhlů ߚ odečtené z grafu na obr. 13.7. Systém značení úhlů viz obr. 11.3. ߚଵǡ ൌ െͶሾιሿ ߚଵ ൌ െሾιሿ
- 65 -
Machovo číslo pro relativní rychlost odečtené z grafu na obr. 13.8 ܽܯଵ ൌ Ͳǡͳሾെሿ Úhel náběhu ݅ ൌ ȟߚ ൌ ȁߚଵ െ ߚଵ ȁ ൌ ȁሺെሻ െ ሺെͶሻȁ ൌ ʹͲሾιሿ
(13.2)
݅௧ ൌ ͶͶǡͷͷ͵ െ Ͷͻǡͷ ή Ͳǡͳ ͳͶǡͳ ή Ͳǡͳଶ ൌ ͳͻǡ͵ሾιሿ Vzhledem ke skutečnosti, že ݅ ݅௧ lze uvažovat, že v tomto režimu dochází k utržení proudění a kompresor je v nestabilním režimu. Pro úplnost je pro tento provozní bod vykresleno i rozložení entropie (obr. 13.9) a meridiální rychlosti (obr. 13.10).
Obr. 13.7 Rozložení úhlů β proudu (beta_fl mean) a lopatky na středním proudovém vláknu (beta_bl mean)
- 66 -
Obr. 13.8 Rozložení Machova čísla relativní rychlosti na středním proudovém vláknu
Obr. 13.9 Rozložení entropie v oběžném kole - 67 -
Obr. 13.10 Rozložení meridiální rychlosti v oběžném kole
13.7 Výpočet ucpání kompresoru Stejně tak jak je důležité stanovit omezení kompresorové charakteristiky z levé strany mezí pumpování, je nezbytné určit omezení pravé strany maximálním hmotnostním tokem, který je kompresor schopný zpracovat. Dle programové dokumentace [34] není matematický model programu vhodný pro výpočet aerodynamického ucpání, protože dochází k průměrování proudění v obvodovém směru a není tak možné identifikovat proudění s vysokým ܽܯna sací straně lopatky. I přes takové omezení obsahuje model schopný stanovit určitý odhad ucpání kompresoru. Tento model je přesnější než 1-D přiblížení, protože bere v úvahu změnu Machových čísel po výšce lopatky. Výstupem výpočtu je parametr Choke Ratio, který představuje poměr aktuálního průtoku a maximálního možného průtoku v daném místě. Proto maximální možnou hodnotou, kterou může tento parametr nabývat je Choke Ratio = 1. Na obr. 13.11 je výsledek výpočtu, který značí ucpání lopatkové mříže.
- 68 -
Obr. 13.11 Rozložení parametru Choke Ratio po řezu kanálu oběžného kola
- 69 -
13.4 Výsledná charakteristika kompresoru pro 2-D výpočet Výše uvedené poznatky jsou aplikovány při návrhovém výpočtu kompresorové charakteristiky na obr. 13.12 a při stanovení izoentropické účinnosti (obr. 13.13). Průběh izoentropické účinnosti je pro všechny počítané rychlosti oběžného kola bez očekávaného trendu klesání účinnosti s blížící se mezí pumpování a ucpáním. K výraznému poklesu izoentropické účinnosti dochází pouze v místech vypočítaného ucpání oběžného kola, kde již parametr Choke Ratio dosahuje hodnoty 1. Tento trend poukazuje na silnou závislost výpočtu na hodnotě odhadované polytropické účinnosti zadávané v základním nastavení výpočtu (tab. 13.1).
n = 150000 min-1 2,4 mez pumpováni aerodynamické ucpání 2,2
2
Πk = p02/p01 [-]
n = 130000 min-1 1,8
n = 110000 min-1 1,6
n = 90000 min-1 1,4
n = 70000 min-1 1,2
T1 = 288 [K] p1 = 0,1013 [MPa] 1 0
0,05
0,1
0,15
0,2
m [kg/s] Obr. 13.12 Výsledná charakteristika kompresoru pro 2-D výpočet ve VistaTF - 70 -
0,25
0,9
0,85
ηiz [-]
0,8
0,75 n = 150000 min-1 n = 130000 min-1 n = 110000 min-1 n = 90000 min-1 n = 70000 min-1 0,7
0,65
T1 = 288 [K] p1 = 0,1013 [MPa] 0,6 0
0,05
0,1
0,15
0,2
m [kg/s] Obr. 13.13 Izoentropická účinnost oběžného kola dle 2-D výpočtu ve VistaTF
- 71 -
0,25
14.
NÁVRH ÚPRAVY GEOMETRIE KOMPRESORU PRO ROZŠÍŘENÍ STABILNÍ OBLASTI
Jednou z možností jak rozšířit stabilní pracovní oblast kompresoru je dle kap. 10 použití dozadu zahnutých lopatek oběžného kola. Návrh úpravy oběžného kola tady spočívá ve změně úhlu ߚଶǡ .
14.1 Úprava geometrie kompresoru v programu BladeGen Úprava výstupního úhlu lopatky ߚଶǡ je aplikován na model oběžného kola, jehož postup tvorby je popsán v kap. 11. Výstupní úhel ߚଶǡ je zvolen ߚଶǡ ൌ െ͵ͷሾιሿ vzhledem systému značení na obr. 11.3. Výsledné rozložení úhlu ߚ lopatky na vrstvách modelu je na obr. 14.1. Rozložení tloušťky lopatky na vrstvách modelu zůstává nezměněno a tloušťka lopatky zůstává konstantní ͳሾ݉݉ሿ (Obr. 14.2). Výsledný 3–D model oběžného kola je na obr. 14.3.
Obr. 14.1 Rozložení ߚ úhlů lopatky na vrstvách modelu oběžného kola
- 72 -
Obr. 14.2 Rozložení tloušťky lopatky na vrstvách modelu oběžného kola
Obr. 14.3 3-D model oběžného kola s návrhem úpravy na rozšíření stabilní oblasti
- 73 -
14.2 Výsledná otáčková větev upraveného oběžného kola V programu VistaTF je na upraveném modelu oběžného kola s dozadu zahnutými lopatkami proveden výpočet otáčkové větve charakteristiky. Otáčky oběžného kola jsou voleny ݊ ൌ ͳ͵ͲͲͲͲ݉݅݊ିଵ. Vypočtená otáčková větev je na obr. 14.4 porovnána s otáčkovou větví oběžného kola s radiálními lopatkami. Dle [29] lze u oběžných kol s dozadu zahnutými lopatkami očekávat oproti lopatkám radiálním pokles stlačení a posun ucpání kompresoru směrem k menším průtokům. Výpočet těmto tvrzením vyhovuje, přičemž dochází i k mírnému posunu meze pumpování k nižším průtokům. Skutečný hmotnostní průtok, při kterém dojde k pumpování, ovšem nelze z důvodů softwarových limitací jednoznačně stanovit (kap. 13.6).
2,4 dozadu zahnuté lopatky radiální lopatky 2,2
Πk = p02/p01 [-]
2
1,8
1,6
1,4
1,2
T1 = 288 [K] p1 = 0,1013 [MPa]
1 0
0,05
0,1
0,15
0,2
m [kg/s]
Obr. 14.4 Srovnání otáčkové větve oběžného kola s dozadu zahnutými lopatkami a oběžného kola s lopatkami radiálními - 74 -
15. ZÁVĚR Cílem práce bylo provést studii proveditelnosti numerického modelování meze pumpování kompresoru pomocí dostupného CFD nástroje. Práce byla směřována pro firmu Honeywell, ze které vzniknul požadavek na zkušební návrhový výpočet oběžného kola programem VistaTF systému ANSYS. Teoretická část práce se věnuje studii pracovní oblasti kompresoru a jevům, které stabilní pracovní oblast limitují. Jedná se především o pumpování, utržení proudu a aerodynamické ucpání kompresoru. Před vlastním výpočtem v programu VistaTF byl v programu Excel 2007 zpracován jednodimenzionální výpočet oběžného kola, který vychází ze vztahů termodynamiky ideálního plynu, energetické bilance, rovnice kontinuity a středních výpočtových rychlostí. Při výpočtu byl uvažován i jednoduchý ztrátový model. Výsledné hodnoty byly použity k odhadu pracovních oblastí kompresoru při různých provozních režimech. Výpočet v programu VistaTF vyžaduje spolupráci s několika dalšími programy systému ANSYS. Model oběžného kola byl vytvořen v programu BladeGen, ze kterého byl následně importován do programu DesignModeler. DesignModeler již umožňuje vytvoření geometrických dat, se kterými dokáže VistaTF pracovat. Následně byl zpracován výpočet výstupních vlastností proudění při různých provozních podmínkách oběžného kola. Během práce v programu byly ovšem objeveny výrazné softwarové limitace vylučující komplexní výpočet proudění v oběžném kole. V práci jsou také diskutovány možnosti rozšíření stabilní pracovní oblasti. Jednou z možností je použití dozadu zahnutých lopatek. Tato úprava byla aplikována na geometrii oběžného kolo z předcházejícího výpočtu. Nicméně vzhledem k již výše zmíněným softwarovým limitacím VistaTF nebylo dostatečně možné prokázat vliv na rozšíření stabilní oblasti oběžného kola.
- 75 -
SEZNAM POUŽITÝCH ZDROJŮ [1]
ANDERSEN, Johannes; LINDSTÖM, Fredrik; WESTIN, Fredrik. Surge Definitions for Radial Compressors in Automotive Turbochargers. SAE International Journal of Engines [online]. 2009, vol. 1, n. 1, [cit. 2010-10-29]. Dostupný z WWW:
. ISSN 0148-7191.
[2]
KOUSAL, Milan. Spalovací turbíny. 2., zcela přeprac. vydání. Praha : Státní nakladatelství technické literatury, 1980. 623 s.
[3]
KADRNOŽKA, Jaroslav. Lopatkové stroje. Brno : Akademické nakladatelství CERM, 2003. 177 s. ISBN 80-7204-297-1.
[4]
MACEK, Jan. Spalovací turbiny, turbodmychadla a ventilátory : přeplňování spalovacích motorů. Vyd. 4. Praha : Nakladatelství ČVUT, 2006. 206 s. ISBN 8001-03529-8.
[5]
Garrett : Products [online]. c2009-2010 [cit. 2010-12-18]. Product Catalog. Dostupné z WWW: .
[6]
MEULEMAN, Corina H.J. Measurement and Unsteady Flow Modelling of Centrifugal Compressor Surge [online]. Eindhoven : Technische Universiteit Eindhoven, 2002. 134 s. Dizertační práce. Technische Universiteit Eindhoven. Dostupné z WWW: . ISBN 90-3862564-2.
[7]
WILLEMS, Franciscus P.T. Modeling and Bounded Feedback Stabilization of Centrifugal Compressor Surge [online]. Eindhoven : Technische Universiteit Eindhoven, 2000. 106 s. Dizertační práce. Technische Universiteit Eindhoven. Dostupné z WWW: . ISBN 90-3862931-1.
[8]
HOFMANN, Karel. Turbodmychadla, vozidlové turbiny a ventilátory : přeplňování spalovacích motorů. 2., nezm. vyd. Praha : SNTL, 1985. 134 s.
[9]
WATSON, N.; JANOTA, M.S. Turbocharging the Internal Combustion Engine. London : The Macmillan Press, 1982. 623 s. ISBN 0333242904, ISBN 9780333242902.
[10]
DIXON, S.L.; HALL, C.A. Fluid Mechanics and Thermodynamics of Turbomachinery. 6th edition. Burlington : Butterworth-Heinemann, 2010. 459 s. ISBN 9781856177931, ISBN 978-1856177931.
[11]
PAMPREEN, Ronald C. Compressor Surge and Stall. 1st edition. Norwich (Vermont) : Concepts ETI, 1993. 457 s. ISBN 0-933283-05-9.
[12]
VLK, František. Paliva a maziva motorových vozidel. 1. vyd. Brno : František Vlk, 2006. 376 s. ISBN 80-239-6461-5. - 76 -
[13]
SCHOLZ, C. Motorová paliva [online]. Liberec : Technická univerzita v Liberci, 2009. 12 s. Učební Text. Technická univerzita v Liberci. Dostupné z WWW: .
[14]
PULKRABEK, Willard W. Engineering Fundamentals of the Internal Combustion Engine. New Jersey : Prentice Hall, 1997. 426 s. ISBN 0135708540.
[15]
HIERETH, Hermann; PRENNINGER, Peter. Charging the Internal Combustion Engine [online]. Wien : Springer-Verlag, 2007 [cit. 2011-03-09]. Dostupné z WWW: . ISBN 978-3-21133033-3.
[16]
MCMILLAN, Gregory K. Centrifugal and Axial Compressor Control [online]. [s.l.] : Instrument Society of America, 1983 [cit. 2011-05-04]. Dostupné z WWW: . ISBN 0-87664-744-1.
[17]
MCKEE, R.J.; EDLUN, C.E.; PANTERMUEHL, P.J. Development on an Active Surge Control System [online]. 2000 [cit. 2011-05-20]. Dostupné z WWW: .
[18]
HOFMANN, Karel. Alternativní pohony [online]. Brno : Vysoké učení technické v Brně, Fakulta strojního inženýrství, [2008] [cit. 2010-24-10]. Dostupné z WWW: .
[19]
JIANG, Wei; KHAN, Jamil; DOUGAL, Roger A. Dynamic centrifugal compressor model for system simulation.Journal of Power Sources [online]. 2006, vol. 158, n. 2, [cit. 2010-11-15]. Dostupný z WWW: .
[20]
JAPIKSE, David. Centrifugal Compressor Design and Performance. Norwich (Vermont) : Concepts ETI, 1996. 387 s. ISBN 0-933283-03-2.
[21]
BOTHA, Barend W.; MOOLMAN, Adriaan. Determing the Impact of the Different Losses on Centrifugal Compressor Design. In R & J Journal [online]. [s.l.] : SA Mechanical Engineer, 2005 [cit. 2011-03-20]. Dostupné z WWW: .
[22]
LI, Shuqi; ZHANG, Junyue; XING, Weidong. Two-Zone Modeling Prediction Method of Centrifugal Compressor Performance. In XU, Jianzhong, et al. Fluid Machinery and Fluid Mechanics : 4th International Symposium (4th ISFMFE) [online]. [s.l.] : Springer-Verlag, 2009 [cit. 2011-04-02]. Dostupné z WWW: . ISBN 978-3540-89748-4.
- 77 -
[23]
JÍLEK, Adolf. Proudění v částech odstředivého kompresoru [online]. Brno : Univerzita obrany, 2007. 38 s. Písemná práce ke státní doktorské zkoušce. Univerzita obrany, Fakulta vojenských technologií, Katedra letecké a raketové techniky. Dostupné z WWW: .
[24]
HUNZIKER, R.; DICKMANN, H-P; EMMRICH, R. Numerical and experimental investigation of a centrifugal compressor with an inducer casing bleed system. In Proceedings of the Institution of Mechanical Engineers, Part A : Journal of Power and Energy [online]. London : SAGE Publications, 2001 [cit. 2011-05-07]. Dostupné z WWW: .
[25]
SAADELDIN ELDALIL, Khaled M. Experimental Study of Centrifugal Compressor Characteristics with Slotted Impeller blades. [online]. 2010 [cit. 2011-04-21]. Dostupné z WWW: .
[26]
CAME, P.M.; ROBINSON, C.J. Centrifugal compressor design. In Proceedings of the Institution of Mechanical Engineers, Part C : Journal of Mechanical Engineering Science [online]. London : SAGE Publications, 1999 [cit. 2011-04-10]. Dostupné z WWW: .
[27]
HARRIS, William. HowStuffWorks [online]. 2006 [cit. 2011-02-24]. How Superchargers Work. Dostupné z WWW: .
[28]
GARRETT, T.K.; NEWTON, K.; STEEDS, W. The Motor Vehicle. 13th edition. Oxford : Butterworth-Heinemann, 2001. 1214 s. ISBN 0750644494.
[29]
KADRNOŽKA, Jaroslav. Tepelné turbiny a turbokompresory : základy teorie a výpočtů. 1. vyd. Brno : CERM, 2004. 308 s. ISBN 80-720-4346-3.
[30]
LI, Wen-Guang. Blade Exit Angle Effects on Performance of a Standard Industrial Centrifugal Oil Pump. [online]. 200? [cit. 2011-04-11]. Dostupné z WWW: .
[31]
CASEY, Michael; ROBINSON, Chris. A New Streamline Curvature Throughflow Method for Radial Turbomachinery. In ASME Conference Proceedings [online]. [s.l.] : ASME, 2008 [cit. 2011-04-20]. Dostupné z WWW: .
[32]
Viscosity. In Wikipedia : the free encyclopedia [online]. St. Petersburg (Florida) : Wikipedia Foundation, 24 February 2002, last modified on 1 April 2011 [cit. 201103-27]. Dostupné z WWW: .
- 78 -
[33]
GRAVDAHL, J.T. Modeling and Control of Surge and Rotating Stall in Compressors [online]. Trondheim (Norway) : Norwegian University of Science and Technology, 1998. 141 s. Dizertační práce. Norwegian University of Science and Technology, Department of Engineering Cybernetics. Dostupné z WWW: .
[34]
ANSYS, Canonsburg (Pennsylvania). TurboSystem [online]. 2010 [cit. 2011-01-05]. Dostupné z WWW: .
[35]
NEJEDLÝ, Milan. Entropy in Compressor Design. 2011. 3 s.
[36]
OH, H.W.; YOON, E.S.; CHUNG, M.K. Systematic two-zone modelling for performance prediction of centrifugal compressors. In Proceedings of the Institution of Mechanical Engineers, Part A : Journal of Power and Energy[online]. London : SAGE Publications, 2002 [cit. 2011-04-28]. Dostupné z WWW: .
[37]
TurboByGarrett [online]. 2009-2010 [cit. 2011-03-12]. Turbo Tech 102. Dostupné z WWW: .
- 79 -
SEZNAM HLAVNÍCH POUŽITÝCH OZNAČENÍ ܽ ܣ ܾ ܾכ ܿ ܿ ܥ ܥ ܦ ܦ ݂ ܪ௨ ݅
[m×s-1] [m2] [m] [-] [m×s-1] [J×kg-1×K-1] [m] [-] [m] [-] [-] [kJ×kg-1] [°] [-]
ܶ ܶ௦ ݑ ܸ ݓ ݓ௨
[J×kg-1] [m] [m] [kg×kg-1] [kg×s-1] [kg×s-1] [kg×s-1] [-] [min-1] [min-1] [Pa] [kW] [kJ×s-1] [J×kg-1×K-1] [-] [K] [K] [m×s-1] [m3] [m×s-1] [m×s-1]
்ܻ ݖ
[-] [-]
ߙ
[°] [-]
ߚ
[-]
ܮଵ ܮ ܮ௧ ݉ሶ ݉ሶ ݉ሶ௦ ܽܯ ݊ ݊ெ ܲ ܳሶ ܴ ܴ݁
rychlost zvuku průřez šířka kola poměr výstupní šířky kola a vstupní šířky difuzoru celková rychlost proudu měrná tepelná kapacita Sutherlandova konstanta ztrátový koeficient povrchového tření průměr difuzní faktor třecí faktor spodní výhřevnost paliva úhel náběhu počet otáček motoru na jeden pracovní oběh (i = 1 pro dvoudobý motor, i = 2 pro čtyřdobý), (2.5) entalpie axiální délka oběžného kola délka středního průtočného vlákna hmotnost vzduchu teoreticky potřebná pro spálení 1 kg paliva hmotnostní tok hmotnostní tok paliva přiváděného do pracovního prostoru válce hmotnostní tok spalovacího vzduchu zůstávajícího ve válci Machovo číslo otáčky oběžného kola otáčky motoru tlak užitečný výkon motoru tepelný tok přiváděný palivem do pracovního prostoru válce měrná plynová konstanta Reynoldosovo číslo teplota celková teplota vzduchu v sání motoru obvodová rychlost zdvihový objem všech válců motoru relativní rychlost rozdíl mezi obvodovou složkou ideální a skutečné relativní rychlosti poloha ventilu počet lopatek úhel celkové rychlosti od radiálního směru součinitel propláchnutí spalovacího prostoru (ߙ ൌ ͳ pro plně vypláchnutý spalovací prostor, ߙ ൌ Ͳ pro nevypláchnutý), (2.5) úhel relativní rychlosti od obvodového směru - 80 -
ο݅ௗ ο݅ ο݅௦௦ ο݅௫ ο݅௦ ȟߟ ȟߟ ȟߟ௩ ߝ
[-] [-] [-] [-] [-] [-] [-] [-] [-] [-]
ߟ
[Pa×s] [-] [-] [-] [-] [-] [-] [-] [m2×s-1] [-] [kg×m-3] [-] [-] [-]
ߟ ߟ ߟ ߢ
ߣ Ɋ ߥ ߎ ߩ ߮ ߯ ߰
Použité indexy Ͳ ͳ ʹ ͵ Ͷ ͷ ͲͲ Ͳͳ Ͳʹ ܾ ܿ ݀ ݄ ݄݀ݕ ݅݀ ݅ݖ ݇ݐ݅ݎ ݉ ݇
ztráty při průtoku oběžným kolem ztráty na vstupu do oběžného kola suma ztrát směšovací ztráty třecí ztráty na povrchu mezilopatkového kanálu pokles účinnosti vlivem zpětného proudění pokles účinnosti vlivem mezery mezi kolem a skříní pokles účinnosti vlivem ztrát ve spirální skříni kompresní poměr poměr sekundárního průtočného průřezu k celkovému průřezu, (12.43) dynamická viskozita (12.41) účinnost užitečná účinnost motoru účinnost stupně kompresoru plnící účinnost Poissonova konstanta součinitel přebytku vzduchu ve válci skluz kinematická viskozita poměrné stlačení celkových tlaků hustota proudícího vzduchu průtokový součinitel podíl sekundárního hmotnostního toku tlakový součinitel
-
řez před vstupními lopatkami vstup do oběžného kola výstup z oběžného kola vstup do lopatkového difuzoru hrdlo lopatkového difuzoru výstup z lopatkového difuzoru celkový stav před vstupními lopatkami celkový stav na vstupu do oběžného kola celkový stav na výstupu z oběžného kola lopatka celkový dynamický vnitřní průměr hydraulický ideální izoentropický kritický meridiální oběžné kolo - 81 -
ݎ ݈݁ݎ ݂݁ݎ ݏ ݐݏ ݄ݐ ݑ
-
redukovaná hodnota radiální vztažené k relativní rychlosti referenční vnější průměr sání motoru (kap.2) statický hrdlo obvodový
- 82 -