Příběh uživatele
V tomto příkladu spočítáme průměrný součinitel síly a součinitel momentu pro Tokijská polytechnická univerzita (TPU) jako experimentální příklad (M0/S0), například ty, které se vztahují na proces statických výpočtů na základě WTG-Merkblatt-M3.
Tato část výsledků patří do „skupiny 1“ podle „obrázku 2.2“ ve WTG-Merkblatt-M3:
- G2: Absolutní hodnoty se středními požadavky na přesnost. Oblast použití může zahrnovat parametry nebo předběžné studie, pokud jsou plánovány pozdější studie s vyšší přesností (např. zkouška ve větrném tunelu třídy G3).
- R2: Samostatné, všechny relevantní směry větru s dostatečně jemným rozlišením směru.
- Z1: Statistické střední hodnoty, pokud se týkají stacionárních procesů proudění, kde lze kolísání (např. v důsledku přibližující se turbulence proudění) dostatečně zachytit jinými opatřeními.
- S1: Statické účinky. Stačí znázornit statický model s nezbytnými mechanickými detaily, ale bez hmoty a tlumicích vlastností.
Popis
Toto srovnání, které zkoumá aerodynamické chování v různých rovinách, má za cíl poskytnout komplexnější pochopení rozložení sil. Tyto poznatky jsou zásadní pro hodnocení celkové stability modelu, aerodynamického výkonu a odezvy konstrukce za podmínek proudění. Průměrné profily rychlosti a intenzity turbulence proudění reprodukované ve větrném tunelu odpovídají profilům terénní kategorie IV podle normy Japonského architektonického institutu (AIJ). Exponent 0,25 uvedený v některých číslech byl po dotazu klasifikován jako typografická chyba. Předpokládaný model je znázorněn na obrázku 1:
Předpoklady, z nichž vychází analýza a simulace, jsou shrnuty v tabulce 1, která nabízí jasný přehled parametrů a podmínek zohledněných v celé studii.
Tabulka 1: Vstupní data 3D obdélníkového modelu
| Parametr | Symbol | Hodnota | Jednotka |
|---|---|---|---|
| Referenční rychlost větru | UH | 11 | m/s |
| Výška střechy | Href | 0,4 | m |
| Exponent profilu | α | 0,25 | - |
| Kategorie terénu | - | IV | - |
| Hustota vzduchu – RWIND | ρ | 1,25 | kg/m³ |
| Turbulentní model – RWIND | RANS K-Omega | - | - |
| Kinematická viskozita – RWIND | ν | 1,5×10⁻⁵ | m²/s |
| Pořadí schématu – RWIND | Druhé | - | - |
| Reziduální cílová hodnota – RWIND | 10⁻⁴ | - | - |
| Reziduální typ – RWIND | Tlak | - | - |
| Minimální počet iterací – RWIND | 800 | - | - |
| Mezní vrstva – RWIND | NL | 10 | - |
| Typ stěnové funkce – RWIND | Vylepšená / Smíšená | - | - |
Studie výpočetní sítě
Obrázek 2 shrnuje studii citlivosti sítě a ukazuje, že s nárůstem hustoty sítě z 20 % na 55 % se součinitel síly (Cf) zvyšuje z 0,95 na 1,05 a poté se stabilizuje. To naznačuje, že simulace dosahuje nezávislosti na síti při 55 %, což znamená, že další zahuštění sítě již významně nemění výsledky. Obrázek zdůrazňuje význam kvality sítě pro zajištění spolehlivých výsledků simulace větru.
Studii vytvoření sítě je nutné provést podle následujícího odkazu:
Obrázek 3 představuje diagram, který porovnává výsledky RWIND s experimentálními daty, pokud jde o vykreslení normalizovaných rychlostních profilů podle výšky. Výsledky vykazují velmi dobrou shodu.
Požadavky na přesnost WTG-Merkblatt
WTG-Merkblatt M3 poskytuje dvě klíčové metody pro ověření výsledků simulace. Metoda „Hit Rate“ vyhodnocuje, kolik simulovaných hodnot Pi správně odpovídá referenčním hodnotám Oi v rámci definované tolerance, a to pomocí binární klasifikace (hit nebo miss). Tento přístup hodnotí spolehlivost simulace výpočtem míry úspěšnosti q, podobně jako funkce spolehlivosti používané v teorii spolehlivosti. Naopak metoda „normalizované střední kvadratické chyby“ (e2) nabízí podrobnější hodnocení přesnosti kvantifikací průměrné kvadratické odchylky mezi simulovanými a referenčními hodnotami, normalizované tak, aby zohledňovala rozdíly ve stupnici. Společně tyto metody poskytují jak kvalitativní, tak kvantitativní měřítka pro ověření simulace.
Výsledky a diskuse
Tabulka 3 ukazuje velmi dobrou shodu mezi výsledky RWIND a referenčními údaji WTG pro normalizované hodnoty rychlosti. Všechny odchylky jsou v přijatelném rozmezí (pod 10 %), takže je dosaženo „míry úspěšnosti“ „q = 100 %“ a „normalizovaná střední kvadratická chyba“ „e2 = 0,00001“ je extrémně nízká. Tyto výsledky potvrzují, že program RWIND přesně reprodukuje referenční profil větru a splňuje přísná validační kritéria.
Tabulka 3: Metrika validace pro rychlostní profil
| WTG - u / uref | RWIND - u / uref | Odchylka (%) | n |
|---|---|---|---|
| 0,394 | 0,390 | 1,045 | 1,00 |
| 0,478 | 0,473 | 1,093 | 1,00 |
| 0,566 | 0,564 | 0,291 | 1,00 |
| 0,629 | 0,629 | 0,081 | 1,00 |
| 0,675 | 0,674 | 0,183 | 1,00 |
| 0,713 | 0,712 | 0,072 | 1,00 |
| 0,750 | 0,753 | 0,422 | 1,00 |
| 0,784 | 0,783 | 0,098 | 1,00 |
| 0,822 | 0,819 | 0,362 | 1,00 |
| 0,869 | 0,865 | 0,463 | 1,00 |
| 0,897 | 0,901 | 0,496 | 1,00 |
| 0,939 | 0,940 | 0,113 | 1,00 |
| 1,010 | 1,005 | 0,493 | 1,00 |
| 1,065 | 1,070 | 0,540 | 1,00 |
| 1,123 | 1,124 | 0,044 | 1,00 |
| 1,161 | 1,162 | 0,121 | 1,00 |
| 1,195 | 1,195 | 0,051 | 1,00 |
| 1,237 | 1,240 | 0,263 | 1,00 |
| 1,266 | 1,270 | 0,291 | 1,00 |
| 1,299 | 1,302 | 0,223 | 1,00 |
Obrázek 4 ilustruje, jak jsou v programu RWIND definovány zóny plochy pro přesný výpočet součinitelů síly na každé ploše modelu budovy vystavené působení větru. 3D model s jednoduchým obdélníkovým tvarem má každé ploše (přední, zadní, levé, pravé a horní) přiřazenu samostatnou zónu označenou barevným kódem. Tyto zóny umožňují programu RWIND samostatně analyzovat a vypočítat aerodynamické síly a rozložení tlaku nezávisle pro každý povrch. Místo umístění četných bodových sond jsme se rozhodli definovat povrchové zóny, protože program RWIND může přímo určit součinitel síly pro každou definovanou zónu.
Ve středu obrázku 4 je otevřeno okno „Upravit zónu“ pro zónu č. 2 – přední jako příklad. Tento podrobný výstup obsahuje klíčové aerodynamické parametry. Promítnutá plocha této zóny ve směru větru je 0,04 m². Vypočtená odporová síla na této ploše je 1,7 N a odpovídající součinitel síly (Cx) je 0,56. V následujícím vzorci je jako příklad vypočítán součinitel síly větru (aerodynamický součinitel síly) pro konkrétní tlakovou zónu. Tuto hodnotu lze poté porovnat s hodnotou zobrazenou na záložce Info v dialogovém okně Upravit data modelu.
Obrázek 5 zobrazuje výsledky následného zpracování simulace větru provedené pro vysokopodlažní budovu pomocí modelu turbulence k-omega. Na pravé straně obrázku je srovnávací tabulka, která obsahuje koeficienty síly získané z CFD simulace spolu s referenčními hodnotami z pokynů WTG. Přední plocha má součinitel síly 0,562, zatímco referenční hodnota WTG je 0,540, což odpovídá odchylce 4,07 %. Zadní plocha, která je vystavená podtlaku v důsledku jevů v úplavu, vykazuje mírně větší odchylku 6,09 %. Boční plochy také vykazují malé rozdíly s odchylkami v rozmezí 3,68 % až 5,38 %. Je pozoruhodné, že celkový součinitel síly ve směru větru, který představuje globální odporový účinek, vykazuje pouze odchylku 1,95 %, což podtrhuje přesnost a spolehlivost CFD simulace při zachycení celkového aerodynamického chování konstrukce.
Skutečné hodnoty síly v newtonech pro každou stranu budovy jsou uvedeny pod touto tabulkou. Přední plocha je vystavena kladné odporové síle 1,7 N, zatímco zadní, pravá a levá plocha jsou vystaveny záporným silám v důsledku sání a bočního odporu. Celková čistá síla ve směru větru činí 3,2 N, což odpovídá celkové síle uvedené v souhrnu sítě. V simulaci byla použita referenční rychlost větru 11 m/s a plocha čela 0,04 m².
Tabulka 4 ukazuje, že program RWIND přesně replikuje referenční průměrné součinitele síly WTG na všech plochách, s odchylkami v rozmezí od 1,95 % do 6,09 % a ‚‘'q=100 % jako míra shody‘‘. Nízká normalizovaná střední kvadratická chyba e2=0,0015 potvrzuje výraznou shodu mezi simulací a měřeními, čímž účinně splňuje validační standardy.
Metrika validace pro součinitele síly mezi WTG a RWIND
| Součinitel síly | Cf,průměr – WTG | Cf – RWIND – k-omega | Odchylka (%) | n |
|---|---|---|---|---|
| Přední plocha | 0,540 | 0,562 | 4,07 | 1,00 |
| Zadní plocha | –0,528 | –0,496 | 6,09 | 1,00 |
| Pravá boční plocha | –0,829 | –0,860 | 3,68 | 1,00 |
| Levá boční plocha | –0,847 | –0,893 | 5,38 | 1,00 |
| Celková síla ve směru větru | 1,070 | 1,058 | 1,95 | 1,00 |
| Celková síla ve směru bočního větru | 0,018 | 0,020 | 10,00 | 1,00 |
| Moment ve směru větru | 0,582 | 0,520 | 10,65 | 1,00 |
| Moment ve směru bočního větru | 0,010 | 0,0099 | 1,00 | 1,00 |