4.

AERODYNAMIKA PROFILOVÝCH MŘÍŽÍ

4.3 . . . . . . . . . . . . . .
4.5 . . . . . . . . . . . . . .
4.5 . . . . . . . . . . . . . .
4.8 . . . . . . . . . . . . . .
4.9 . . . . . . . . . . . . . .
4.10 . . . . . . . . . . . . . .
4.11 . . . . . . . . . . . . . .
4.13 . . . . . . . . . . . . . .
4.14 - 4.19 . . . . . . . . .
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.

Klíčová slovaKapitola: Základní pojmy aerodynamiky profilových mříží
4.3

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

Aerodynamika leteckých profilů

Odpor

Vztlak

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, Aerodynamika profilů]. To znamená, že u profilu zařazeného do profilové mříže lze také rozlišovat odpor a vztlak, které 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 střední rychlosti a vztlak kolmý na tuto rychlost, viz Vzorce 1.

Vzorce pro odpor a vztlak lopatky v lopatkové mříži
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.

Součinitel vztlaku

Součinitel odporu

Aerodynamický tunel

Japikse,1997

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].

Aerodynamický tunel pro měření lopatkových mříží
2: Aerodynamický tunel pro měření lopatkových mříží
Klíčová slovaKapitola: Základní pojmy aerodynamiky profilových mříží
4.4
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.

Aerodynamický tunel

Zkroucené lopatky

Nadzvuková rychlost

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.

Součinitel odporu

Zakřivení proudu

Odtržení proudu

Nátokový úhel

Aerodynamická charakteristika

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ů.

Aerodynamická charakteristika profilu v lopatkové mříži
3: Aerodynamická charakteristika profilu v lopatkové mříži
Δβ [°] ú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.
Klíčová slovaKapitola: Základní pojmy aerodynamiky profilových mříží
4.5

Lopatkový kanál

Howell, A. R.

Carter, A. D. S.

Lieblein, S.

Japikse,1997

Dixon and Hall, 2010

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 kanálem lze indikovat tlakovou ztrátu, která vzniká zejména vnitřním tření pracovní tekutiny. 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.

Tlaková ztráta profilové mříže
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.

Klíčová slovaKapitola: Profilová ztráta profilové mříže
4.6

Vnitřní tření

Mezní vrstva

Laminární proudění

Turbulentní proudění

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é turbuletní 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].

Reynoldsovo číslo

Letecký profil

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.

Drsnost povrchu

Tření v mezní vrstvě je ovliněno i kvalitou povrchu (drsnost) lopatky. Drsnost povrchu lopatky nelze příliš ovlivnit, protože pracovní tekutiny lopatkových strojů obvykle obsahují nečistoty, byť v malé koncentraci, a ty 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.

Změna drsnosti obtékaného povrchu v čase.
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ší.

Odtržení proudu

Mechanismus odtržení proudu od profilu je už 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.

Difuzorová mříž

U difuzorových profilových mříží dochází k odtržení vždy, viz také článek Proudění plynů a par difuzory [Škorpík, 2023b].

Klíčová slovaKapitola: Profilová ztráta profilové mříže
4.7

Víření

Odtoková hrana

Odtokové hrany lopatek nejsou ostré především z pevnostních důvodů, takže mezi proudy od sací a přetlakové stranách 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.

Vířením za odtokovou hranou
6: Vířením za odtokovou hranou

Rázové vlny

Kritické Machovo číslo

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].

Energetická bilance

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.

Energetická bilance profilových ztrát lopatkové mříže
7: Energetická bilance profilových ztrát lopatkové mříže
(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.
Klíčová slovaKapitola: Profilová ztráta profilové mříže
4.8

Poměrná ztráta

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 profilových mříží nepřekračují jednotky procent.

Poměrná profilová ztráta lopatkové mříže
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.

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 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

Ideální proudění vs. Proudění se ztrátami

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ýstupu z mříže, tak aby rychlosti a rychlostní trojúhelníky byly stejné jako při proudění beze ztrát. To znamená, že podle věty o změně hybnosti tekutiny budou stejné i složky síly působící ve směru-θ (tzv. obvodový směr), protože ty jsou funkcí pouze změn složek rychlostí v tomto směru, a naopak se zvýší složky sil působící ve směru-a (tzv. axiální směr). Z těchto předpokladů lze zakreslit jednotlivé složky síly působící na lopatku v rovinné lopatkové mříži, tak jako na Obrázku 9.

Rozklad sil působící na přímou lopatku určité délky v rovinné lopatkové mříži
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.
Klíčová slovaKapitola: Výpočet sil působící na lopatku v rovinné lopatkové mříži
4.9

Součinitel vztlaku

Součinitel odporu

Profilová ztráta

Tlaková ztráta

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í.

Rovnice pro výpočet aerodynamických veličin profilové mříže
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.

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

Transformační rovnice

Nožička, 1967

Japikse,1997

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í lopatové mříže. V takovém případě se tvar lopatek nerovinných mříží musí transformovat bod po bodu na rovinné mříže. Na Obrázku 11 je příklad transformace tvaru kruhové lopatkové mříže na rovinnou pomocí transformačních rovnic. Podobným způsobem lze transformovat i diagonální lopatkovou mříž [Nožička, 1967, s. 84], [Japikse, 1997, s. 7-24]. Stejně tak se musí transformovat tlakové a rychlostní pole naměřené na rovinné lopatkové mříži zpět do kruhových souřadnic postupem odvozeným v [Nožička, 1967, s. 84]. Tyto tranformační rovnice jsou odvozeny za zjednodušujícího předpokladu potenciálního proudění.

Klíčová slovaKapitola: Přenositelnost aerodynamiky rovinné...diagonálních a radiální lopatkových mříží
4.10
Transformační rovnice tvaru kruhové profilové mříže na rovinnou
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říží

Aerodynamika leteckých profilů

Aerodynamická interference

Abbott, 1959

Při návrhu profilových mříží s málo prohnutými profily levných kusových strojů, kde si nelze dovolit aerodynamický výzkum profilové mříže, lze využít dostupné osamocené profily k sestavení profilové mříže. 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.

Součinitel vztlaku

Weiningův součinitel

Lakshminar., 1996

Wang and Kruyt, 2022

Změnu součinitele vztlaku 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 profilových mříží [Wang and Kruyt, 2022].

Klíčová slovaKapitola: Transformace aerodynamických veličin...na aerodynamické veličiny profilových mříží
4.11
Weiningův součinitel
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.

Větrné turbíny

Vrtule

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

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.

Základní profil

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.

Úlohy

Ú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
Klíčová slovaKapitola: Úlohy
4.12
Ú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.
Rychlostní trojúhelník rotoru axiálního ventilátoru
U [m·s-1] obvodová rychlost.
§1   zadání:   Δp; W1; W2; ρ §3   výpočet:   ξh
§2   výpočet:   Lh; Lp    
Popisek symbolů je v Příloze 2.
Ú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
§1   zadání:   Z; cξh §4   výpočet:   s; σ
§2   odečet:   W2; W1; β1; β2; rm §5   výpočet:   Lh; CD; CL,is; CL
§3   výpočet:   W; W; Wθm; Wm; βm    
Popisek symbolů je v Příloze 3.
Klíčová slovaKapitola: Úlohy
4.13
Ú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.
§1   zadání   Δp; Wa; ρ; rm; N §9   porovnání:   Lh §4 vs. §8
§2   odvození:   rov. pro σ §10   výpočet:   γ; CL
§3   výpočet:   U; W1; β1 §11   odečet:   ν
§4   odhad:   Lh   výpočet:   c, s, Z
§5   výpočet:   W2, β2, σ   zaokrouhlení:   Z
§6   návrh:   Re   výpočet:   s, σ, b
§7   odečet:   CLisolated, iisolated=im, CDisolatedCD §12   kontrola:   CL,is
§8   výpočet:   Wm, βm, Lh    
Popisek symbolů je v Příloze 4.

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, [online], ISSN 1804-8293. 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 efekty 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