autor:
ŠKORPÍK, Jiří, ORCID: 0000-0002-3034-1696
datum vydání:
Červen 2022
název:
Aerodynamika profilů
sborník:
provenience: Brno (Česká republika)
email: skorpik.jiri@email.cz
Copyright©Jiří Škorpík, 2023 |
Co je a k čemu je aerodynamikaAerodynamika vyšetřuje silové účinky proudění na obtékaná tělesa nebo kanál, ve kterém se nachází. Aerodynamika profilu je speciální oblast aerodynamiky, která zkoumá 2D obtékání těles, které mají proměnné rozměry pouze ve dvou směrech – jedná se o profilovanou desku, viz Obrázek 1. Profilem je nazýván řez tímto tělesem. Popis geometrie profilů a pravidel pro jejich kreslení je uveden například v článku Tvary lopatek a průtočných částí lopatkových strojů [Škorpík, 2023]. 1: Aerodynamický profil-příklad P-rovina řezu; Ψ-proudnice ve vyšetřované rovině; A-znázornění vyšetřovaného profilu. l [m] délka profilové desky; x, y, z [m] souřadnice.Nestlačitelná tekutinaRychlost zvukuMachovo čísloPři nízkých rychlostech se plyn chová podobně jako kapalina, která je nestlačitelná, ale s rostoucí rychlostí se začíná projevovat vliv nízké rychlosti zvuku v plynech, které jsou v kapalinách mnohem vyšší. Rychlost zvuku je totiž rychlostí šíření tlakových poruch v tekutině. Navíc při nadzvukových rychlostech proudění může docházet k jevům souvisejících se skokovými změnami stavových veličin plynu, ke kterým při nižších rychlostech nedochází. Nicméně stlačitelnost plynu se na proudění významně projevuje až od rychlosti cca 0,3 Mach [Kadrnožka, 1991, s. 27]. Z toho důvodu je v tomto článku popsána teorie proudění nestlačitelného proudění a teorie stlačitelného proudění v článku Machovo číslo a efekty při proudění vysokými rychlostmi. Odpor profiluPři obtékaní jakéhokoliv tělesa proudem tekutiny působí tato tekutina na těleso silou. Tato síla vzniká kvůli změně hybnosti tekutiny v důsledku vloženého tělesa, změny tlaku v okolí tělesa a tření tekutiny o plochu tělesa případě kvůli dalším efektům. Složka síla ve směru nátokové rychlosti (směr volného proudu před tělesem) se nazývá odpor. |
Třecí odporTvarový odporPříspěvek k odporu od třecích sil se nazývá zkráceně třecí odpor, a od tlakových tlakový odpor nebo tvarový odpor, který vzniká vždy, protože tlakové diference Δp mezi čelní a zádní plochou profilu vzniknou i při proudění bez tření. Výsledný odpor profilu je jen přibližně součtem obou dvou odporů, viz Rovnice 2(a), protože se vzájemně ovlivňují. Základní profilNa Obrázku 2(b) je příklad vzniku odporu při obtékaní základního profilu (symterický profil). V tomto případě veškeré síly ve směru osy-z mají nulovou výslednici, ve směru osy-x nikoliv. Tvarový odpor v tomto směru vznikne tak, že na nátokové straně profilu je vyšší tlak než na odtokové. 2: Princip vzniku sil od tlaku a tření pří proudění tekutiny kolem profilu (a) rovnice síly působící na profil od proudu tekutiny; (b) příklad vzniku sil působící na profil od proudu tekutiny při obtékání symetrického profilu. V∞ [m·s-1] nátoková rychlost (rychlost před obtékaným profilem); p [Pa] statický tlak u profilu; S [m2] povrch profilové desky; F [N] síla (index p od tlaku, index F od tření); D=Fx [N] odpor; x [m] souřadnice ve směru proudění.Vztlak profiluNátokový úhelV případě, že nátoková rychlost nemá směr osy profilu nebo je profil nesymetrický, pak má výsledná síla F působící na profil dvě na sebe kolmé složky. Složka ve směru nátokové rychlosti se nazývá opět odporem a složka kolmá na směr proudění vztlakem, viz Obrázek 3. 3: Definice vztlaku CL-střední čára profilu (osa profilu) (ang. camaber line); ε [°] klouzavý úhel; c [m] tětiva profilu; i [°] nátokový úhel; L=Fy [N] vztlak. |
Klouzavý úhelÚhel mezi vztlakem a odporem profilu nazýváme klouzavým úhlem. Název pochází z klouzavého letu letounu (ustálený bezmotorový let), kdy síla F působí proti směru gravitačního zrychlení a ε označuje sklon vztlaku vůči gravitačnímu zrychlení, viz Obrázek 4. 4: Princip klouzavého letu U [m·s-1] rychlost letounu; m [kg] hmotnost letounu; g [m·s-2] gravitační zrychlení.Průběh tlaku kolem profiluPřetlaková stranaSací stranaVztlakJestliže při obtékání profilu vznikne vztlak L, pak to znamená, že tlaková diference mezi dvěma body profilu, které leží proti sobě kolmo na směr nátokové rychlosti, je různá od nuly. Na Obrázku 5 je průběh tlaku kolem profilu, u kterého vzniká vztlak. Proto rozlišujeme přetlakovou stranu a sací stranu profilu, tedy stranu profilu v průměru s vyšším tlakem a stranu v průměru s nižším tlakem. Na sací straně profilu je průměrně nižší tlak a vyšší rychlost, protože proudění na této straně má delší trajektorie v okolí profilu než na přetlakové straně profilu. 5: Změna tlaku podél osamoceného profilu V okolí profilu je nejvyšší tlak u nátokové hrany profilu, jejíž stěna je téměř kolmá na nátokovou rychlost a proto v těchto místech dojde k prudkému snížení hybnosti ve směru nátoku a tedy poklesu dynamického tlaku a nárůstu tlaku statického. (a) průběh tlaku na sací straně; (b) průběh tlaku na přetlakové straně.Základní profilI symetrický neboli základní profil může vytvořit přetlakovou a sací stranu, pokud je vložen do proudu šikmo, respektive nátokový úhel je různý od nuly. |
XfoilZměnu tlaku po profilu se měří pomocí otvorů v profilu. Konkrétní hodnoty tlakového součinitele z měření jsou uvedeny např. v [Kousal, 1980, s. 142]. Existují ale i velmi přesné numerické modely profilu, které dokaží měření nahradit, například software Xfoil. Tlakový součinitelPři proudění se tedy část kinetické energie tekutiny transformuje na tlakovou energii, takže velikost tlaku v okolí profilu se může maximálně zvýšit o hodnotu dynamického tlaku před profilem podle vzorce pmax=p∞ + 0,5·ρ·V2∞. Pro popis změn statického a dynamického profilu byla zavedena veličina tlakový součinitel (Vzorec 6). Tlakový součinitel profilu CP popisuje jak se mění statický tlak v okolí profilu na úkor dynamického před profilem. Tlakový součinitel profilu může dosahovat maximální hodnoty 1 a jestliže je záporný, pak to znamená, že ve vyšetřovaném bodě poklesl tlak pod tlak před profilem p∞, respektive rychlost proudění v tomo dodě je větší než nátoková rychlost (obrázek u Vzorce 6). 6: Tlakový součinitel (a) průběh tlakového součinitele na sací straně; (b) průběh tlakového součinitele na přetlakové straně. CP [1] tlakový součinitel profilu [Nožička, 1967, s. 27]; ρ [kg·m-3] hustota; V [m·s-1] rychlost ve vyšetřovaném místě v okolí profilu; p∞ [Pa] tlak před profilem.Průběh rychlosti proudění kolem profiluRychlostní sondaPrůběh rychlosti kolem profilu lze měřit přímo rychlostní sondou a nebo vypočítat z průběhu podílu rychlosti ve vyšetřovaném bodě V a nátokové rychlosti V∞, viz Obrázek 7, který lze stanovit z průběhu tlaků podle Rovnice 6. 7: Průběh rychlosti podél osamoceného profilu (a) situace na sací straně; (b) situace na přetlakové straně. |
Projevy vnitřního tření v mezní vrstvě profiluTření tekutiny o plochy profilové desky nevytváří pouze odpor, ale také tzv. mezní vrstvu7., ve které dochází ke ztrátě tlakové a kinetické energie tekutiny. Tato disipace energie způsobuje jednak odklon proudu za profilem a jednak může způsobyt odtržení proudu od profilu: Deviační úhelProtože rychlosti na sací a přetlakové straně profilu jsou různé, je různá i jejich vnitřní tření v mezní vrstvě a tedy i hybnost na odtokové hraně, kde by se oba proudy slévají. To znamená, že směr proudění těsně za profilem nemusí mít směr daný tvarem profilu, respektive směrem střední čáry profilu, protože převáží směr proudu s vyšší hybností. Odklon tohoto směru od střední čáry profilu CL se nazývá deviační úhel, Obrázek 8. 8: Deviační úhel δ [°] deviační úhel.Odtržení prouduV blízkosti odtokových hran, kde má proudění v mezní vrstvě nejnižší energie, může dojít k odtržení proudu (kolaps mezní vrstvy v důsledku průniku okolní tekutiny do mezní vrstvy) od profilu a tedy i ke změně sil působící na profil – dokonce i ke ztrátě potřebného vztlaku, viz Obrázek 9. K nějakému odtržení dojde téměř vždy. Velikost odtržení se zvětšuje s nakloněním profilu vůči směru nátokové rychlosti, respektive se zakřivením profilu. Samotný proces odtržení je často přetržitý (nestacionární), po odtržení se může mezní vrstva opět stabilizovat (sníží se tlak v oblasti odtokové hrany) a záhy odtrhnout. 9: Hlavní toky pracovní tekutiny při odtržení proudu od profilu |
Výpočet sil působící na profil od proudu tekutinyOdporVztlakSoučinitel vztlakuSoučinitel odporuPolára profiluVztahy pro výpočet odporu a vztlaku profilu lze odvodit z definice tlakového součinitele, viz Vzorce 10. Odpor i vztlak profilové desky jsou přímo úměrné její velikosti, dynamického tlaku před deskou aerodynamickými součinitely zvanými součinitel odporu CD a součinitel vztlaku CL, které zohledňují tvar profilu desky, nátokový úhel a charakter proudění viz Vzorce 10. Pro použití těchto vzorců je nutné mít naměřeny hodnoty součinitele odporu a vztlaku vyšetřovaného profilu a to pro různé nátokové úhly i a Reynoldsova čísla Re. Tato aerodynamická data se publikují ve formě tabulek nebo grafu. 10: Rovnice vztlaku a odporu profilu (a) vztahy pro výpočet odporu a vztlaku profilové desky; (b) graf závislosti L(D, i), tzv. polára profilu. CD [1] součinitel odporu; CL [1] součinitel vztlaku. Odvození vzorců pro odpor a vztlak profilu je v Příloze 1.Aerodynamická dataV současné době již existují rozsáhlé katalogy naměřených aerodynamických dat pro různé tvary profilu nebo projektilů, začít lze třeba zde: [Abbott and Doenhoff, 1959], [Hošek, 1949, s. 390], [Kneubuehl, 2004, s. 76]. Existují i neveřejné, soukromé katalogy firem. U profilů se zároveň měří i další aerodynamické charakteristiky (např. profilový moment [Hošek, 1949, s. 278], [Abbott and Doenhoff, 1959, s. 4], který je důležitý pro návrh polohy těžiště letounu a pro pevnostní výpočet délky křídla/lopatky lopatkových strojů na krut).
OdkazyŠKORPÍK, Jiří, 2023, Tvary lopatek a průtočných částí lopatkových strojů, Transformační technologie, Brno, [on-line], ISSN 1804-8293. https://turbomachinery.education/tvary-lopatek-a-prutocnych-casti-lopatkovych-stroju.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.
HOŠEK, Josef, 1949, Aerodynamika vysokých rychlostí, Naše vojsko, Praha.
KADRNOŽKA, Jaroslav, 1991, Teorie lopatkových strojů, Vysoké učení technické, Brno, ISBN 80-214-0275-X.
KNEUBUEHL, Beat, 2004, Balistika střely, přesnost střelby, účinek, Naše vojsko, Praha, ISBN 80-206-0749-8.
KOUSAL, Milan, 1980, Spalovací turbíny, Nakladatelství technické literatury n. p., Praha.
NOŽIČKA, Jiří, 1967, Analogové metody v proudění, Academia, Praha.
©Jiří Škorpík, LICENCE
|