4.

AERODYNAMIKA PROFILOVÝCH MŘÍŽÍ

–    strana 4.3    –
–    strana 4.5    –
–    strana 4.5    –
–    Úloha 1:    –
–    strana 4.8    –
–    Úloha 2:    –
–    Úloha 3:    –
–    strana 4.10    –
–    strana 4.11    –
–    Úloha 4:    –
–    strana 4.13    –
– strana 4.14-4.21 –
strana 4.2
– autor: –
ŠKORPÍK, Jiří (LinkedIn.com/in/jiri-skorpik)
–    datum vydání:    –
Září 2022
– název: –
Aerodynamika profilových mříží
– sborník: –
– provenience: –
Brno (Česká republika)
– email: –
skorpik.jiri@email.cz

Copyright©Jiří Škorpík, 2006-2022
Všechna práva vyhrazena.

 AERODYNAMIKA PROFILOVÝCH MŘÍŽÍ
strana 4.3

Základní pojmy aerodynamiky profilových mříží

Aerodynamika profilové mříže vychází z definic aerodynamických veličin používáných při popisu aerodynamiky leteckých profilů [Škorpík, 2022]. To znamená, že u profilu zařazeného do profilové mříže lze také rozlišovat odpor a vztlak. Aerodynamické veličiny profilových mříží se také měří v aerodynamických tunelech jako v případě leteckých profilů.

Vzorce pro odpor a vztlak lopatky v lopatkové mříži
Základní aerodynamické veličiny profilové mříže lze vypočítat podle stejných vzorců jako v případě osamocených profilů, s tím rozdílem, že místo nátokové rychlosti se vychází ze střední aerodynamické rychlosti v mříži, takže odpor je ve směru této rychlosti a vztlak kolmý na tuto rychlost, viz Vzorce 1.
– 1: –
Vzorce pro odpor a vztlak lopatky v lopatkové mříži
(b) řez rovinnou lopatkovou mříží; (b) vzorce pro výpočet odporu a vztlaku přímé lopatky v rovinné lopatkové mříži; (c) vztahy mezi rychlostmi a silami působící na přímou lopatku. D [N] odpor; L [N] vztlak; F [N] výsledná síla působící na lopatku; ρ [kg·m-3] hustota pracovní tekutiny; Wm [m·s-1] střední aerodynamická rychlost v mříži; c [m] tětiva profilu; γ [°] úhel nastavení profilu v mříži; b [m] šířka lopatkové mříže; l [m] délka přímé lopatky; CD, CL [1] součinitel odporu a vztlaku profilu v profilové mříži; s [m] rozteč. θ, a-označení os v soustavě souřadnic.
Aerodynamický tunel pro měření lopatkových mříží
Aerodynamické součinitelé C profilů v mříži jsou měřeny v aerodynamických tunelech rovinných lopatkových mříží, viz Obrázek 2 pro různou kombinaci rychlostí a nátokových úhlů. Lopatková mříž v aerodynamické tunelu je tvořena několika (5 až 7) přímými lopatkami vloženými do průtokového kanálu šikmo tak, aby proud pracovní tekutiny odpovídal směru relativní rychlosti ve skutečné lopatkové mříži. Vstupní a výstupní průřez tunelu odpovídá průtočným průřezům testované lopatkové mříže. Konstrukce aerodynamického tunelu kompresorových lopatkových mříží jsou uvedeny např. v [Japikse, 1997, s. 11-7] a pro turbínové lopatkové mříže [Japikse, 1997, s. 6-22].
– 2: –
Aerodynamický tunel pro měření lopatkových mříží
Kanál aerodynamického tunelu na výtoku je tvořen pohyblivými stěnami, kterými se ovlivňuje rychlostní pole na okrajích mříže a umožňuje naklápění vstupních a výstupních kanálů pro změnu nátokového úhlu. W1, 2 [m·s-1] rychlost před a za mříži; p1, 2 [Pa] tlak před a za mříži.
 AERODYNAMIKA PROFILOVÝCH MŘÍŽÍ
strana 4.4
Speciální aerodynamické tunely
Používají se i aerodynamické tunely pro měření zkroucených lopatek, kde se provádí měření i po délce lopatek. Dokonce se používají i aerodynamické tunely, ve kterých lze měřit aerodynamiku celého stupně při rotaci. Konstrukce takového zkušebního zařízení je uvedena např. v [Japikse, 1997, s. 6-23]. Poslední speciální skupinou aerodynamických tunelů jsou nadzvukové aerodynamické tunely, ve kterých se zkoumá i vliv efektů spojených s nadzvukovými rychlostmi.
Měřené aerodynamické veličiny a aerodynamická charakteristika profilové mříže v aerodynamickém tunelu lopatkové mříže
Naměřená aerodynamická data jsou publikována v grafické podobě, viz Obrázek 3. Na obrázku je jasně patrné, že součinitel odporu CD má nějaké minimum a zakřivení proudu Δβ nějaké maximum. To je dáno tím, že stabilní proudění lopatkovou mříží je jen v určitém intervalu nátokových úhlů i. Mimo tento interval docházi k Odtržení proudu od profilu – v případě velkého nátokového úhlu na sací straně profilů – nebo od určité záporné hodnoty nátokového úhlu na přetlakové straně profilů.
– 3: –
Aerodynamická charakteristika profilové mříže
Δβ [°] úhel zakřivení proudu; i [°] nátokový úhel nátokové rychlosti W1 – index opt označuje optimální nátokový úhel, kdy lopatky dosahují maximální hodnoty vztlaku vzhledem k hodnotě odporu, index n označuje navržený jmenovitý nátokový úhel – jedná se o takový úhel, kterému odpovídá zakřivení proudu přibližně 0,8·Δβmax – zde ještě má lopatková mříž dostatečné rezervy pro změny provozních parametrů, anichž by docházelo k odtrhávání proudu od profilu. a-stabilní oblast proudění bez odtrhávání proudu od profilu; b-oblast odtržení na sací straně lopatky; c-oblast odtrhávání na přetlakové straně lopatky.
 AERODYNAMIKA PROFILOVÝCH MŘÍŽÍ
strana 4.5
Predikce změn v aerodynamické charakteristice profilové mříže pomocí Howeelovy, Carterovy nebo Liebleinovy metody
Lopatkové kánaly v lopatkových strojích svůj tvar a velikost po délce mění, tak jak se mění profil a rozteč lopatek, které tyto kanály vymezují. Zejména u zkroucených lopatek dochází k velkým změnám nejen geometrických ale i aerodynamickým parametrům proudění, takže pro aerodynamické výpočty je nutné mít k dispozici velké množství měření. Měření aerodynamických veličin pro kombinaci tolika proměnných je časově i finančně nákladné a proto byly vyvinuty výpočtové metody predikce změny aerodynamických veličin profilů při změně některých parametrů, např. metoda Howelova, Carterova, Liebleinova...[Japikse, 1997, s. 5-14], [Dixon and Hall, 2010]. Tyto metody mají užitek především při snaze snížit potřebné množství měření při drobných změnách na profilu a profilové mříži.

Tlaková ztráta profilové mříže

Při proudění jakýmkoli lopatkovým kanálem lze indikovat tlakovou ztrátu, která vzniká v důsledku tzv. profilových ztrát zejména vnitřním tření pracovní tekutiny.

Definice tlakové ztráty v profilové mříži
V případě lopatkových kanálů, respektive lopatkových mříží je tlaková ztráta Lp definována jako rozdíl, o který se musí tlak za mříží p2 snížit, aby rychlost za mříží W2 byla stejná jako rychlost W2is při proudění mříží beze ztrát při výstupním tlaku p2is, viz Rovnice 4. Přičemž stav před mříží je předpokládán stejný pro oba případy.
– 4: –
Tlaková ztráta profilové mříže
Lp [Pa] tlaková ztráta; p2is [Pa] tlak za mříží při proudění beze ztrát; p2 [Pa] tlak za mříží při dosažení stejných rychlostí před a za mříží W1, 2.

Profilová ztráta profilové mříže

Profilovými ztrátami profilové mříže se rozumí nežádoucí transformace energie (disipace energie) v pracovní tekutině, které vnikají při průtoku profilovou mříží. Mezi profilové ztráty patří zejména ztráty vnitřním třením, ztráty při odtržení proudu od profilu, ztráty při víření za odtokovou hranou a ztráty v rázové vlně při proudění vysokými rychlostmi. Suma těchto ztrát vytváří profilovou ztrátu mříže.

 AERODYNAMIKA PROFILOVÝCH MŘÍŽÍ
strana 4.6
Vnitřní tření v krátkých lopatkových kanálech
Ztráty kvůli vnitřnímu tření pracovní tekutiny jsou nevyhnutelné. Charakteristickým rysem lopatkových kanálů je jejich malý rozměr a vysoká rychlost, takže mezní vrstva se obvykle plně nevyvine a zůstává laminární, přičemž v jádru proudu mohou být turbulence, pokud před lopatkových kanálem bylo vyvinuté turbulentní proudění. Více o přechodech proudění mezi turbulentním laminárním v článku Vnitřní tření tekutiny a vývoj mezní vrstvy [Škorpík, 2023].
Vliv délky tětivy na Reynoldsovo číslo
K predikci druhu mezní vrstvy v okolí profilu je důležité tzv. Reynoldsovo číslo, které je mimo jiné přímo úměrné délce tětivy c, takže pro stejné rychlosti mají profily s malou délkou tětivy malou hodnotu Reynoldsova čísla a naopak. To je také důvod proč se uvádějí aerodynamické součinitelé leteckých profilů pro velká Reynoldosva čísla – očekává se jejich aplikace jako profily křídel letadel, jenže délka tětivy lopatek lopatkových strojů je řádově měnší než tětiva profilu křídla.
Vývoj drsnosti povrchů lopatek
Na vnitřní tření, respektive tlakovou ztrátu profilové mříže má významný vliv znečištění a drsnost povrchu lopatek. Drsnost povrchu lopatky nelze příliš ovlivnit, pokud pracovní tekutina obsahuje nečistoty, které reagují s povrchem lopatky (ulpívají na povrchu nebo ho poškozují) a určují jeho drsnost především při delším provozu, viz Obrázek 5.
– 5: –
Změna drsnosti obtékaného povrchu v čase.
(a) na počátku vysoká drsnost; (b) na počátku nízká drsnost. Ra [μm] drsnost povrchu; t [h] doba provozu. Při výrobě je nutné zvážit, jestli se vyplatí výroba velmi hladké lopatky, když se po čase její povrch stane drsnější.
Vliv odtržení proudu na funkce lopatkové mříže
Mechanismus odtržení proudu od profilu je popsán v článku Aerodynamika profilů [Škorpík, 2022]. Kromě energetické ztráty, následkem odtržení může být kolaps průtoku lopatkovou mříží a kmitání lopatek. Rozsah odtržení lze ovlivnit zejména tvarem profilu a nátokovým úhlem, respektive vhodným rozvržením změn tlaku podél profilu, viz Úloha 1. Nejcitlivější na odtrhávání proudu jsou difuzorové profilové mříže, protože současně dochází k nárůstu tlaku za profilovou mříží, viz také článek Proudění plynů a par difuzory [Škorpík, 2023b].
 AERODYNAMIKA PROFILOVÝCH MŘÍŽÍ
strana 4.7
–    Úloha 1:    –
Popište, která z uvedených profilových mříží (viz obrázek) bude pravděpodobně citlivější na održení proudu od profilu.
Řešení úlohy je uvedeno v Příloze 1.
Obrázek k úloze 632
Víření za odtokovou hranou lopatek
Odtokové hrany lopatek nejsou ostré především z pevnostních důvodů, takže mezi proudy od sací a přetlakové strany lopatek je mezera, ve které vznikají drobné víry kvůli odlišným rychlostí těchto proudů, viz Obrázek 6. Přičemž víry jsou typickým jevem zvyšující entropii tekutiny.
– 6: –
Vířením za odtokovou hranou
 
Rázové vlny vznikají při kritickém Machově čísle před profilovou mříží
Rychlost proudění se podél profilu mění, takže i když před mříží nedosahuje proudění rychlosti zvuku, tak v mříži této rychlosti dosáhnout může nebo ji dokonce překonat při určité nátokové rychlosti, kterou označujeme jako kritického Machovo číslo profilové mříže. Při následném poklesu rychlosti zpět na podzvukovou může, za jistých podmínek, vzniknout rázová vlna. Podrobnosti o proudění v profilových mříží vysokou rychlostí jsou uvedeny v kapitole Aerodynamika profilových mříží ve stlačitelném prostředí v [Škorpík, 2023c].
Profilové ztráty profilové mříže se projeví zvýšením entropie a tlakovou ztrátou
Profilová ztráta profilové mříže je ekvivaletní entalpii, o kterou se musí snížit izoentropický rozdíl entalpií za mříži, aby rychlost na výstupu byla stejná jako v případě proudění beze ztrát, viz Obrázek 7.
 AERODYNAMIKA PROFILOVÝCH MŘÍŽÍ
strana 4.8
– 7: –
Energetická rovnováha profilových ztrát lopatkové mříže
Energetická rovnováha profilových ztrát lopatkové mříže v h-s diagramu: (a) situace v konfuzorové mříži; (b) situace v difuzorové mříži. h [J·kg-1] entalpie; s [J·kg-1·K-1] entropie; ps [Pa] celkový tlak; Lh [J·kg-1] profilové ztráty. Index is označuje změny stavu při proudění beze ztrát.
Poměrná profilová ztráta lopatkové mříže
K vyjádření aerodynamické kvality profilové mříže se používá veličina zvaná poměrná profilová ztráta profilové mříže ξh, která je definována jako poměr mezi profilovou ztrátou mříže a celkovou disponibilní entalpii pracovní tekutiny při průtoku mříží, viz Vzorec 8. Hodnoty poměrné profilové ztráty moderních profilů nepřekračují jednotky procent. V případě jednoduchých profilů tepelných turbín běžně používaných kolem roku 1910 se pohybovaly poměrné profilové ztráty kolem 20 %.
– 8: –
Poměrná profilová ztráta lopatkové mříže
(a) konfuzorový lopatkový, nebo rovnotlakový lopatkový kanál; (b) difuzorový lopatkový kanál. ξh [1] poměrná profilová ztráta profilové mříže.
Profilové ztráty jsou součásti vnitřních ztrát
Ve stupních lopatkových strojů se mimo profilových ztrát vyskytují ještě další ztráty, které přímo nesouvisí s profilovými mřížemi. Tyto další ztráty a profilové ztráty se souhrně označují jako vnitřní ztráty, viz článek Vnitřní ztráty lopatkových strojů a jejich vliv na návrh lopatkového stroje.

Výpočet sil působící na lopatku v rovinné lopatkové mříži

Při identifikaci jednotlivých složek sil působící na lopatku v rovinné lopatkové mříže se vychází z toho, že při reálném proudění se změní tlak (klesne) na výtoku z mříže, tak aby rychlostní trojúhelníky byly stejné jako při proudění beze ztrát. Za těchto předpokladů lze sestrojit silový trojúhelník rovinné lopatkové mříže i odvodit vztahy mezi aeordynamickými veličinami profilové mříže.

 AERODYNAMIKA PROFILOVÝCH MŘÍŽÍ
strana 4.9
Konstrukce silového trojúhelníku proudění se ztrátami
Jestliže platí rovnost relativních rychlostí pro případ proudění beze ztráta a se ztrátami, pak, podle věty o změně hybnosti tekutiny aplikovanou na rovinnou lopatkovou mříž, zůstavá stejná i obvodové složka síly. Naopak naopak, axiální složka síly proudění se ztrátami bude jiná než při proudění beze ztrát, kvůli změně tlaku za lopatkovou mříží (podle h-s diagramu na Obrázku 7) a tlakových sil působicí na kontrolní objem lopatky, viz Obrázek 9.
– 9: –
Rozklad sil působící na přímou lopatku určité délky v rovinné lopatkové mříži
(a) konfuzorová mříž; (b) difuzorová mříž. βm [°] úhel střední aerodynamické rychlosti v mříži; F [N] síla; ε‾ [°] klouzací úhel.
Rovnice pro výpočet součinitele odporu, vztlaku, profilové a tlakové ztráty profilové mříže
Rovnice pro výpočet součinitele vztlaku a odporu a profilovou ztrátu profilové mříže (Rovnice 10) lze odvodit na základě předchozího obrázku a z tlakové ztráty. Tyto rovnice lze použít i obráceně k výpočtu tlakové ztráty na základě aerodynamických veličin. Přičemž Rovnice 10d pro součinitel vztlaku pro případ proudění beze ztrát CL,is je odvozena z požadavku na zakřivení proudu a Rovnice 10b je odvozena ze skutečných hodnot aerodynamických součinitelů–pravé strany těchto rovnic se u dobře navržené profilové mříže rovnají.
– 10: –
Rovnice pro výpočet aerodynamických veličin profilové mříže
CL,is [1] součinitel vztlaku pro případ proudění bez profilových ztrát; σ [1] hustota lopatkové mříže. Orientace kótování úhlů nátokové a výtokové rychlosti β1, β2 je stejná jako orientace úhlů βm na Obrázku 9. Odvození rovnic je v Příloze 5.
–    Úloha 2:    –
Vypočítejte tlakovou a poměrnou profilovou ztrátu profilové mříže rotoru axiálního ventilátoru. Zvýšení tlaku v rotorové řadě lopatek je 500 Pa, relativní rychlost vzduchu na vstupu do rotorové řady lopatek je 46,7 m·s-1, na výstupu 35,4 m·s-1. Hustota pracovního plynu je 1,2 kg·m-3.
Řešení úlohy je uvedeno v Příloze 2.
 AERODYNAMIKA PROFILOVÝCH MŘÍŽÍ
strana 4.10
Rychlostní trojúhelník rotoru axiálního ventilátoru
U [m·s-1] obvodová rychlost.
– Úloha 3: –
Vypočítejte součinitel vztlaku a odporu profilu lopatek rotoru Kaplanovy turbíny z Úlohy 8 v článku Základní rovnice lopatkových strojů. Výpočet proveďte na středním poloměru lopatky. Počet lopatek je 6, délka tětivy lopatek na středním poloměru je 1,8 m a poměrná profilová ztráty mříže na vyšetřovaném poloměru je 2,5 %.
Řešení úlohy je uvedeno v Příloze 3.
Lopatková mříž kaplanovy turbíny

Přenositelnost aerodynamiky rovinné lopatkové mříže na aerodynamiku diagonálních a radiální lopatkových mříží

V Úloze 2 je předvedená aplikace aerodynamických dat získaných při měření na rovinné lopatkové mříži na axiální lopatkovou mříž, ale tato data mohou být aplikována i na radiální a diagonální lopatkové mříže. Jestliže mají být zachovány aerodynamické parametry rovinné lopatkové mříže i v mříži radiální, pak je nutné transformovat její rozměry pomocí transformačních rovnic.

Transformační rovnice pro převod rozměrů rovinné lopatkové mříže na radiální mříž
Transformační rovnice pro převod rozměrů rovinné lopatkové mříže na mříž radiální a obráceně při zachování aeorodynamickýchz vlastností lze odvodit pro potenciální proudění lopatkovou mříží, viz Obrázek 11. Transformace tvaru profilu se provádí bod po bodu. Podobným způsobem lze transformovat i diagonální lopatkovou mříž [Nožička, 1967, s. 84], [Japikse, 1997, s. 7-24]. Stejným postupem lze transformovat tlakové a rychlostní pole naměřené na rovinné lopatkové mříži do kruhových souřadnic postupem odvozeným v [Nožička, 1967, s. 84].
 AERODYNAMIKA PROFILOVÝCH MŘÍŽÍ
strana 4.11
– 11: –
Transformační rovnice tvaru kruhové profilové mříže na rovinnou
r [m] poloměr lopatkové mříže.

Transformace aerodynamických veličin leteckých profilů na aerodynamické veličiny profilových mříží

Jestliže by byl výzkum nového typu profilové mříže nákladný, pak lze profilovou mříž sestavit z leteckých profilů alespoň v případech požadavku na malé prohnutí profilů. Výhodou je, že jsou snadno dostupná aerodynmická data z měření velkého množství osamocených profilů, viz například [Abbott and Doenhoff, 1959]. Nevýhodou je, že lze očekávat, že aerodynamické součinitele leteckého profilu budou mít jinou hodnotu, než tyto součinitele stejného profilu, ale zařazeného v profilové mříži. To je dáno jednak aerodynamickou interferencí sousedících profilů a jednak rozdílnou definicí aerodynamických součinitelů – v případě leteckých profilů jsou vztažena k nátokové rychlosti, v případě profilových mříží ke střední aerodynamické rychlosti v mříži. Nicméně existují ověřené postupy pro přepočet těchto součinitelů.

Přepočet součinitele vztlaku leteckého profilu v profilové mříži pomocí Weiningova součinitele
Změnu součinitele vztlaku leteckého profilu po vložení do profilové mříže lze přibližně stanovit pomocí Weiningova součinitele. Weiningův součitel je definován jako podíl součinitele vztlaku profilové mříže složené z plochých desek CL a součinitele vztlaku osamocené ploché desky CLisolated, který je označován písmenem K, viz Vzorec 12 [Lakshminarayana, 1996, s. 212]. Tento poměr je odvozený pro ploché desky, ale pro přibližnou změnu aerodynamických součinitelů málo prohnutých profilů se běžně používá při návrhu i profilových mříží [Wang and Kruyt, 2022].
– 12: –
Weiningův součinitel
CLisolated [1] součinitel vztlaku osamocené desky. Vzorec má platnost přibližně do hodnoty 1/σ= 0,7, od této hodnoty lze vycházet z grafu uvedeného v [Lakshminarayana, 1996, s. 212]. Od hodnoty 1/σ=2,5 je součinitel K blízký hodnotě 1.
 AERODYNAMIKA PROFILOVÝCH MŘÍŽÍ
strana 4.12
Přepočet součinitele vztlaku leteckého profilu u větrných turbín a vrtulí
Mimo Weinigova součinitele existují i jiné způsoby jak určit součinitel vztlaku profilu v mříži, jestliže jsou známy jeho parametry při jeho měření jako osamoceného profilu. Tyto způsoby se používají při výpočtu lopatkových mříží větrných turbín a vrtulí, více v článku Aerodynamika větrných turbín.
Součinitel odporu leteckých profilů není aerodynamickou interferencí s okolními lopatkami příliš ovlivněn
Jestliže nebude docházet k odtrhávání proudu od profilu, lze očekávat, že sounitel odporu profilu v mříži CD nebude tak citlivý na aerodynamickou interferenci jako součinitel vztlaku CL. Takže aerodynamický výpočet takto navržené profilové mříže lze začít tak, že z hodnot součinitele vztlaku a odporu profilu lze podle Rovnice 10 pro součinitel vztlaku CL a navržený úhel střední aerodynamické rychlosti βm vypočítat součinitel vztlaku CL,is atd., viz Úloha 4.
Součinitel odporu základního profilu zůstává přibližně stejný i po prohnutí
V případě profilů s vytvořených prohnutím základního profilu (obvykle profily s prohnutím větším než u leteckých profilů) lze při výpočtu sil na lopatku vycházet z vlastností výchozího základního profilu. Výpočet vychází z předpokladu, že po prohnutí se příliš nezmění součinitel odporu CD profilu, takže na základě výpočtu odporu profilu v mříží D lze stanovit i další hodnoty sil podle Obrázku 9.
– Úloha 4: –
Navrhněte parametry nátokové a výtokové rychlosti a zakřivení proudu a rozteč difuzorové profilové mříže sestavené z leteckého profilu typu NACA 65-410, jestliže je tato profilová mříž určena pro rotor axiálního ventilátoru s požadovaným zvýšením tlaku 500 Pa a nachází se na středním kvadratickém poloměru lopatek. Hustota pracovního plynu je 1,2 kg·m-3, střední kvadratický poloměr lopatek 1050 mm a otáčky ventilátoru jsou 325 min-1, hodnota axiální rychlosti před rotorem je 30 m·s-1.
Řešení úlohy je uvedeno v Příloze 4.
Lopatková mříž axiálního ventilátoru
rm [m] střední poloměr lopatek.
 AERODYNAMIKA PROFILOVÝCH MŘÍŽÍ
strana 4.13

Odkazy

ŠKORPÍK, Jiří, 2022, Aerodynamika profilů, Transformační technologie, Brno, [online], ISSN 1804-8293. https://fluid-dynamics.education/aerodynamika-profilu.html.
ŠKORPÍK, Jiří, 2023, Vnitřní tření tekutiny a vývoj mezní vrstvy, fluid-dynamics.education, Brno, https://fluid-dynamics.education/vnitrni-treni-tekutiny-a-vyvoj-mezni-vrstvy.html.
ŠKORPÍK, Jiří, 2023b, Proudění plynů a par difuzory, Transformační technologie, Brno, [online], ISSN 1804-8293. https://fluid-dynamics.education/proudeni-plynu-a-par-difuzory.html.
ŠKORPÍK, Jiří, 2023c, MACHOVO ČÍSLO A JEVY PŘI PROUDĚNÍ VYSOKÝMI RYCHLOSTMI, Transformační technologie, Brno, [on-line], ISSN 1804-8293. https://fluid-dynamics.education/machovo-cislo-a-efekty-pri-proudeni-vysokymi-rychlostmi.html.
ABBOTT, Ira, DOENHOFF, Albert, 1959, Theory of wing sections, including a summary of airfoil data, Dover publications, inc., New York, ISBN-10:0-486-60586-8.
DIXON, S., HALL, C., 2010, Fluid Mechanics and Thermodynamics of Turbomachinery, Elsevier, Oxford, ISBN 978-1-85617-793-1.
JAPIKSE, David, 1997, Introduction to turbomachinery, Oxford University Press, Oxford, ISBN 0-933283-10-5.
LAKSHMINARAYANA, Budugur, 1996, Fluid Dynamics and Heat Transfer of Turbomachinery, John Wiley & Sons, Toronto, ISBN 0-471-85546-4.
NOŽIČKA, Jiří, 1967, Analogové metody v proudění, Academia, Praha.
WANG, Jie, KRUYT, Niels, P., 2022, "Design for High Efficiency of Low-Pressure Axial Fans With Small Hub-to-Tip Diameter Ratio by the Vortex Distribution Method" J. Fluids Eng., 144(8), https://doi.org/10.1115/1.4053555.
©Jiří Škorpík, LICENCE