|
– autor: –
ŠKORPÍK, Jiří (LinkedIn.com/in/jiri-skorpik)
– datum vydání: –
Leden 2023; Leden 2025 (2. vydání)
– název: –
Aerodynamika větrných turbín
– web: –
– provenience: – Brno (Česká republika)
– email: – skorpik.jiri@email.cz
Copyright©Jiří Škorpík, 2023-2025 |
Rozdělení větrných turbín podle principu transformace kinetické energie větru na práciPodle orientace osy rotace se rozdělují větrné turbíny na horizontální (HAWT– Horizontal axis wind turbine) a vertikální osou rotace (VAWT– Vertical axis wind turbine). Turbíny s horizontální osou rotace pracující na principu lopatkového stroje, ve kterém vzduch koná práci, když protéká lopatkovými kanály rotoru, viz Obrázek 1283a. Turbíny s vertikální osou rotace pracují na principu proměnlivého odporu lopatek během jedné otáčky, takže při stejném směru větru kladou v jednotlivých fázích otáčky rotoru různý odpor, viz Obrázek 1283b. – 1283: – Základní typy rotorů větrných turbín ![]() (a) čelní pohled na větrnou turbínu s horizontální osou rotace a válcový řez A-A její lopatkové mříže; (b) čelní pohled na větrnou turbínu s vertikální osou rotace a řez A-A jejího rotoru (na obrázku je tzv. Darrieova turbína, ale existují i jiné koncepty turbín typu VAWT); (c), (d) detail změny orientace lopatek a jejich aerodynamický odpor při otáčení rotoru – v poloze (c) má lopatka větší odpor než v poloze (d) a tak vzniká moment sil na rotoru, který způsobuje otáčení turbíny. AS [m2] opsaná plocha rotoru (swept area); r [m] vyšetřovaný poloměr rotoru; rt [m] poloměr rotoru; U [m·s-1] obvodová rychlost na vyšetřovaném poloměru; V [m·s-1] absolutní rychlost proudu vzduchu; Vθ [m·s-1] obvodová složka absolutní rychlosti; W [m·s-1] relativní rychlost proudu vzduchu. Index 1 označuje stav před rotorem, index 2 za rotorem.
|
Aerodynamické součinitelé větrných turbínAerodynamické veličiny rotorů větrných turbín poskytují informace pro návrh, výpočet rotoru i o jeho zatížení. Mezi aerodynamické veličiny větrné turbíny patří axiální a obvodový součinitel, které popisují změny rychlostí před a za rotorem, a výkonový a náporový součinitel, které popisují energetické, respektive silové zatížení rotoru.
– 311: – ![]() ST-proudová trubice rotoru (stream tube); SWA-oblast zpomaleného proudu vzduchu (slowed wind area); AS-okolní proudění (ambient stream). A [m2] průtočná plocha; a [1] axiální součinitel (axial induction factor). Index 0 označuje stav na vstupu do proudové trubice, respektive rychlost větru.
– 305: – ![]() a' [1] obvodový součinitel (angular induction factor); Ω [rad · s-1] úhlová rychlost osového víru za rotorem; ω [rad · s-1] úhlová rychlost rotoru. Záporné znaménko v čitateli neguje zápornou hodnotu rychlosti V2θ, respektive úhlovou rychlost víru za turbínou Ω, které jsou záporné (směřují proti kladnému smyslu obvodového směru). |
– 260: – ![]() (a) parametry pro celý rotor; (b) parametry na vyšetřovaném poloměru (lokální hodnota daného součinitele). CP [1] výkonový součinitel; CT [1] náporový součinitel; Pi [W] výkon větru přenesený na rotor větrné turbíny; Pwind [W] kinetický výkon větru; Pd [Pa] dynamický tlak větru; T [N] nápor větru na rotor. ρ [kg·m-3] hustota vzduchu. Rovnice pro kinetický výkon větru Pwind je odvozena v Příloze 260. Porovnání ideálního a reálného rotoru větrné turbínyParametry ideálního rotoru větrné turbíny se označují jako Betzův limit. Reálné rotory těchto parametrů nemohou dosáhnout kvůli ztrátám, které nutně při transformaci kinetické energie větru na práci rotoru vznikají. Navíc u reálných rotorů dochází ke zborcení proudové trubice, což zvyšuje hodnotu náporového součinitele oproti ideálnímu stavu. |
– 326: –
Parametry ideálního rotoru větrné turbíny
a, a', CP, CT [1]. Odvození hodnot je uvedeno v Příloze 326. Jiný (novější) postup odvození je v [Tyagi and Schmitz, 2025].
|
|
– 799: – Výkonový součinitel rotoru větrné turbíny a ztráty podle Betze ![]() LCP [1] ztráty na výkonovém součiniteli; Ut [m·s-1] obvodová rychlost na špici lopatky; λ [1] součinitel rychloběžnosti; λstart [1] součinitel rychloběžnosti, při kterém je větrná turbína schopna samostatného chodu. Index opt označuje optimální parametry. Data z [Hau, 2006, s. 98].
– 615: –
Data z [Hau, 2006, s. 98]
|
|
– 599: – Turbulizátor lopatky větrné turbíny (vortex generators) ![]() (a) detail turbulizátoru na přetlakové straně lopatky; (b) turbulizátor od společnosti 3M s lepicí ploškou (foto 3m.com), který lze použít k dodatečnému vybavení lopatky turbíny (například při použití turbíny v lokalitě s nižší rychlostí větru než na jakou byla navržena) – všimněte si profilu žeber turbulizátoru.
|
– 357: – Náporový součinitel rotoru větrné turbíny ![]() (a) rovnice odvozená pro ideální rotor; (b) aproximace naměřených hodnot náporového součinitele rotoru s počátkem v a=0,3262. Za axiální součinitel a je dosazena střední hodnota po délce lopatek. Odvození rovnice (a) je uvedeno v Příloze 357. Návrh tvaru lopatky větrné turbíny pomocí aerodynamických součinitelůNávrh lopatky se provádí analytickým 2D výpočtem, tj. rozdělení lopatky na elementární stupně, proto též anglický název Blade element method (BEM), viz Obrázek 642, s. 10.10. Předpokladem je, že proudění rotorem nemá radiální složku. Návrh začíná výběrem vhodného profilu a odhadem hodnot axiálního a obvodového součinitele na vyšetřovaném poloměru. Z těchto predikcí lze vypočítat, pomocí teorie leteckého profilu, viz článek Aerodynamika profilů [Škorpík, 2022], rychlosti v rychlostních trojúhelnících, úhel nastavení profilu, délku tětivy i síly působící na vyšetřovaný element lopatky. Odtud lze stanovit i lokální hodnoty výkonového a náporového součinitele. Ověření přesnosti odhadu hodnot axiálního a obvodového součinitele na vyšetřovaném poloměru lze pomocí energetické a silové rovnováhy počítaného elementu lopatky. Celkový výkon turbíny je pak součtem výkonů všech elementů lopatky. |
|
– 642: – Elementární stupeň větrné turbíny ![]() Při 2D výpočtu elementárního stupně větrné turbíny očekáváme, že pro každý vyšetřovaný poloměr lopatky bude vycházet rozdílná délka tětivy profilu c a úhel nastavení profilu γ. U špice lopatky bývají obvykle aerodynamické prvky pro snížení okrajových ztrát případně aerodynamickou regulaci turbíny. c [m] délka tětivy; γ [°] úhel nastavení profilu v mříži; Δr [m] výška elementu lopatky.
|
– 1065: – ![]() Wm [m·s-1] střední aerodynamická rychlost; βm [°] úhel střední aerodynamické rychlosti; Δβ [°] zakřivení proudu; im [°] nátokový úhel střední aerodynamické rychlosti (úhel mezi střední aerodynamickou rychlostí a tětivou). Index m označuje střední aerodynamickou rychlost. Odvození tohoto vzorce je uvedeno v Příloze 1065.
|
|
– 172: – Větrná elektrárna Smith-Putnam ![]() Větrná elektrárna Smith-Putnam z roku 1941 byla konstruována s ohledem na co nejmenší výrobní náklady. Přímé lopatky byly klasické žebrované konstrukce. Vnitřní žebra lopatek byly z oceli a vnější ocelový plášť z nerezové oceli. Lopatky byly dlouhé 20 m a každá vážila 8 tun. Průměr rotoru 53,3 m při jmenovitém výkonu 1 250 kW. Regulace výkonu elektrárny byla pomocí natáčení lopatek. – Úloha 166: –
Vypočítejte základní úhly lopatky rotoru větrné turbíny o průměru 140 m a třemi lopatkami. Výpočet proveďte na poloměru 10 m pro profil NACA 63-209. Hodnoty axiálního a obvodového součinitele zatím pouze odhadněte. Řešení úlohy je uvedeno v Příloze 166.
– 1084: – Rozklad síly působící na lopatku na vyšetřovaném poloměru ![]() θ-obvodový směr; a-axiální směr (osový). D [N] odpor; F [N] síla působící na element lopatky od proudu vzduchu; L [N] vztlak; ε‾ [°] klouzací úhel.
– Úloha 153: –
Vypočítejte délku tětivy a jednotkové síly působící na element lopatky rotoru z Úlohy 166. Výpočet proveďte pro jmenovitou rychlost větru 6,8 m·s-1 a hustotu vzduchu 1,2 kg·m-3. Řešení úlohy je uvedeno v Příloze 153.
|
– 358: – ![]() m [kg·s-1] hmotnostní tok; wE [kg·s-1] Eulerova práce na vyšetřovaném poloměru. (a) výkonový součinitel podle definice Eulerovy práce (Euler work); (b) výkonový součinitel podle teorie leteckého profilu (Airfoil theory); (c) podmínka energetické rovnováhy. Odvození rovnic je uvedeno v Příloze 358.
– 1089: – ![]() (a) podmínka správného návrhu hodnot axiálního a obvodového součinitele s Pradtlovou korekcí; (b) Prandtlův součinitel. CF [1] Prandtltův součinitel. Součin wE·CF představuje střední hodnotu Eulerovy práce na vyšetřovaném poloměru. Prandtlův součinitel CF byl původně odvozen pro vrtule, nicměně se ho používá při výpočtu větrných turbín i v současnosti [Wilson et al., 1976], [Hansen, 2008, s. 52].
– 1088: – ![]() |
|
(a) náporový součinitel odovzený z věty o změně hybnosti (Theorem of momentum change) – vzorec je odovzen již v Příloze 1088 a je totožný se Vzorcem 357, s. 10.9 a jeho platnost je potvrzena příbližně jen do hodnoty a=1/3; (b) náporový součinitel podle teorie leteckého profilu; (c) podmínka silové rovnováhy. Odvození rovnice (b) je uvedeno v Příloze 1088.
– Úloha 900: –
Proveďte výpočet výkonového a náporového součinitele v oblasti elementu lopatky větrné turbíny řešené v Úloze 166 a 153. Vypočítejte také otáčky rotoru. Řešení úlohy je uvedeno v Příloze 900.
– 448: – ![]() Parametry lopatky rotoru větrné turbíny uvedené v Úlohách 166, 153 a 900 s Prandtlovou korekcí. Parametry jsou počítány pro případ, kdy veličiny CL, CD a im jsou konstatní. r [m]; c [m]; γ [m]; CP, CT [1]; CF [1]. Data pro konstrukci grafů jsou uvedena v Příloze 448.
Výkonová charakteristika větrné elektrárnyObvyklou výkonovou charakteristikou větrných elektren bývá závislost jejího výkonu P na rychlosti větru V0 pro provozní otáčky, viz Obrázek 1098, s. 10.16. Tuto charakteristiku lze zkonstruovat z bezrozměrné charakteristiky větrných turbín CP-λ [Manwell et al., 2002, s. 129] nebo přesněji výpočtem výkonu rotoru větrné turbíny pro jednotlivé rychlosti větru (v Úloze 461, s. 10.15 je uveden výpočet výkonu pro optimální rychlost větru) a současným započítáním ztrát v mechanismech, elektrickém generátoru a vlastní spotřeby větrné elektrárny. |
|
– Úloha 461: –
Vypočítejte očekávaný výkon rotoru větrné turbíny z řešené v Úlohách 166, 153 a 900. Konečné parametry lopatek jsou uvedeny na Obrázku 448, s. 10.14. Vypočítejte také výkonový a náporový součinitel tohoto rotoru. Řešení úlohy je uvedeno v Příloze 461.
Provoz větrných elektrárenRegulace větrné elektrárny reaguje na změny rychlosti větru s ohledem na výkonovou charakteristiku větrné elektrárny a její bezpečnost. Regulace větrné turbíny začíná jejím startem a lze ji provádět čtyřmi základními způsoby a to aerodynamickou regulací, natáčením lopatek, změnou otáček a bočením.
|
|
– 988: – ![]() (a) výkonová charakteristika asynchronního generátoru; (b) výsuvně-natáčívá aerodynamická brzda na špici lopatky v aktivním stavu. T [N·m] kroutící moment; N [s-1] otáčky. Index shutdown označuje stav, kdy je nutné větrnou turbínu odstavit brzdou; index G označuje generátor; index n označuje jmenovité parametry (nominal).
– 1098: – Výkonová charakteristika větrné elektrárny ![]() modře-větrná elektrárna s aerodynamickou regulací; černě-větrná elektrárna s regulací natáčením lopatek.
|
– 1095: – Schéma zapojení multifrekvenčního generátoru větrné elektrárny ![]() (a) charakteristika multifrek. generátoru; (b) elektrárny s multifrek. generátorem musí obsahovat také výkonovou elektroniku pro změnu frekvence dodávaného proudu. 1-multifrek. generátor; 2-usměrňovač střídavého proudu na proud stejnosměrný; 3-střídač stejnosměrného proudu; 4-připojení elektrárny k přenosové soustavě. AC-střídavý proud; DC-stejnosměrný proud.
|
|
– 526: – ![]() Konstrukce lopatky větrné turbínyKonstrukce lopatky větrné turbíny je většinou odvislá od velikosti turbíny. Nejmenší turbíny mají lopatky z jednoho kusu materiálu, tzv. plné lopatky. Větší lopatky jsou duté s žebrovanou konstrukcí uvnitř podobné konstrukcím křídel letadel. Lopatky velkých větrných turbín jsou skořepinové konstrukce s několika vrstvami a s výztuhami, viz Obrázek 596. – 596: – Jednotlivé vrstvy skořepinové lopatky ![]()
– 911: – Výroba lopatky větrné turbíny ![]() Na obrázku je fáze nanášení vrstev skelných vláken před zalitím polyesterem v továrně společnosti LM Wind Power [Thomsen, 2004]. |
– 218: – Testování odolnosti lopatky při úderu bleskem ![]() Testování odolnosti lopatky v továrně společnosti LM Wind Power (Nizozemsko) [Thomsen, 2004].
|
|
– 419: – Lopatka větrné turbíny po renovaci svého povrchu – první po dvaceti letech provozu ![]() Lopatka o délce 28 m větrné turbíny Vestas V52.
OdkazyŠKORPÍK, Jiří, 2022, Aerodynamika profilů, fluid-dynamics.education, Brno, https://fluid-dynamics.education/aerodynamika-profilu.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.
ANDERSON, M. B., 1980, A vortex-wake analysis of a horizontal axis wind turbine and a comparison with modified blade element theory, Proceedings of the Third International Symposium on Wind Energy Systems, Copenhagen, BHRA Fluid Engineering, paper no. H1 357–374.
DIXON, S., HALL, C., 2010, Fluid Mechanics and Thermodynamics of Turbomachinery, Elsevier, Oxford, ISBN 978-1-85617-793-1.
HANSEN, Martin, 2008, Aerodynamics of wind turbines, Earthscan Ltd., London, ISBN 978-1-84407-438-9.
HAU, Erich, 2006, Wind Turbines – fundamentals, technologies, Applications, Economics, Springer, Berlin, Heidelberg, New York, ISBN-10-3-540-24240-6.
LIEW, Jaime, HECK, Kirby, HOWLAND, Michael, 2024, Unified momentum model for rotor aerodynamics across operating regimes, Nature Communications 15, 6658, https://doi.org/10.1038/s41467-024-50756-5.
MANWELL, J.F., MCGOWAN, J.G. and ROGERS A.L., 2002, Wind Energy Explained – Theory, Design and Application, John Wiley & Sons Ltd, ISBNS: 0-471-49972-2.
MILLER, Rudolf, HOCHRAINER, A., LÖHNER, K., PETERMANN, H., 1972, Energietechnik und Kraftmaschinen, Rowohlt taschenbuch verlag GmbH, Hamburg, ISBN 3-499-19042-7.
SHARPE, D. J., 1990, Wind turbine aerodynamics, Wind Energy Conversion Systems, Englewood Cliffs.
STIESDAL, Henrik, 1999, The wind turbine components and operation, Bonus Energy A/S, Brande. Dostupné z http://www.windmission.dk.
THOMSEN, Troels, 2004, Reliability of large rotor blades, AusWIND 2004 in July 28 – 30, Launceston, Tasmania.
TYAGI, Divya, SCHMITZ, Sven, 2025, Glauert's optimum rotor disk revisited – a calculus of variations solution and exact integrals for thrust and bending moment coefficients, Wind Energy Science, 10, 451–460, https://doi.org/10.5194/wes-10-451-2025.
WILSON, R., LISSAMAN, P., WALKER, Stel, 1976, Applied Aerodynamics of Wind Power Machines, Oregon State University.
|