Příběh uživatele
Tento příklad prezentuje experimentální a numerická data pro aerodynamické síly působící na průřez mostovky [1]. Tato data jsou široce využívána jako referenční hodnoty pro validaci CFD simulací a pro posouzení aeroelastického chování v aplikacích větrného inženýrství. Proudění kolem průřezů mostů představuje komplexní aerodynamický problém zahrnující odtržení proudění, interakci smykových vrstev a vývoj úplavu. V závislosti na geometrii a úhlu náběhu se mohou vyskytovat jak stacionární, tak nestacionární aerodynamické jevy, včetně odtrhávání vírů, gallopingu a flutteru.
Ve srovnání s jednoduchými tupými tělesy vykazují mostovky chování proudění silně závislé na geometrii. Malé geometrické detaily (např. ostrost hran, zábradlí, aerodynamické kryty) mohou významně ovlivnit rozložení tlaku a součinitele sil. Proto je přesná predikce těchto jevů pro numerické simulace obzvláště náročná. Tento příklad zásadně souvisí se Skupinou 4, protože se zabývá aeroelastickými jevy, zejména flutterem a chováním interakce tekutiny a konstrukce (FSI) průřezů mostů. Nicméně, protože plně sdružená obousměrná FSI simulace není v současné době v RWIND k dispozici, přímá reprezentace aeroelastických jevů není možná.
Proto je přijat zjednodušený přístup. V důsledku toho je tento příklad v RWIND považován za příslušející do Skupiny 1, podle Obrázek 2.2 ve WTG-Merkblatt M3, na základě vyhodnocení střední rychlosti větru a odpovídajících středních aerodynamických veličin:
- G1: Kvalitativní hodnoty s nízkými požadavky na přesnost pro použití v základním výzkumu nebo předběžném návrhu. Úsilí a požadavky na úroveň podrobnosti jsou sníženy, protože často nejsou všechny okrajové podmínky plně objasněny.
- R1: Samostatně stojící (bez okolních budov), analýza jednotlivých důležitých směrů větru.
- Z1: Statistické střední hodnoty, za předpokladu, že se týkají stacionárních procesů proudění, kde mohou být fluktuace (např. v důsledku turbulence nabíhajícího proudu) dostatečně zachyceny jinými opatřeními.
- S1: Statické účinky. Postačuje reprezentovat konstrukční model s nezbytnými mechanickými detaily, ale bez hmotových a tlumicích vlastností.
Popis
Zkoumaný případ se zaměřuje na aerodynamické chování průřezu mostovky vystavené stacionárnímu nabíhajícímu proudu při různých úhlech náběhu. Hlavním cílem této studie je stanovení stacionárních aerodynamických součinitelů sil, se zvláštním důrazem na součinitel odporu CD. Na základě stacionárních CFD simulací je vyhodnocena střední odporová síla působící na průřez a následně je převedena do bezrozměrného tvaru pomocí referenčního dynamického tlaku a charakteristických rozměrů. To umožňuje výpočet CD a konzistentní srovnání s experimentálními daty, stejně jako jeho aplikaci v konstrukčním návrhu.
Kromě součinitele odporu mohou být také hodnoceny další globální součinitele, jako je vztlak CL a moment CM, aby byla poskytnuta úplnější aerodynamická charakteristika mostovky. Primární důraz však zůstává na přesné predikci CD jako klíčového parametru pro zatížení ve směru větru.
Výchozím bodem je dynamický tlak větru, definovaný jako qo=1/2 ρU2, kde ρ je hustota vzduchu a U je střední rychlost větru. Tato veličina představuje kinetickou energii proudění vzduchu na jednotku objemu a slouží jako reference pro všechna aerodynamická zatížení.
Tabulka 1: Vstupní data průřezů mostu
| Parametr | Symbol | Hodnota | Jednotka |
|---|---|---|---|
| Rychlost volného proudu | u | 8,2 | m/s |
| Výška střechy | Href | 180 | mm |
| Hustota vzduchu – RWIND | ρ | 1,25 | kg/m³ |
| Model turbulence – RWIND | RANS K-Omega | - | - |
| Kinematická viskozita – RWIND | ν | 1,5×10⁻⁵ | m²/s |
| Řád schématu – RWIND | Druhý | - | - |
| Cílová hodnota reziduí – RWIND | 10⁻⁴ | - | - |
| Typ reziduí – RWIND | Tlak | - | - |
| Minimální počet iterací – RWIND | 800 | - | - |
| Mezní vrstva – RWIND | NL | 10 | - |
| Typ stěnové funkce – RWIND | Standardní | - | - |
Studie výpočetní sítě
Obrázek 2 představuje analýzu citlivosti sítě válcového modelu v RWIND. Vypočtený součinitel síly (Cf) mírně klesá z 0,76 při hustotě sítě 15 % na 0,71 při 25 % a dále na 0,70 při 35 %. Toto postupné snižování naznačuje, že se řešení stabilizuje s tím, jak je síť zjemňována. Malá odchylka Cf při vyšších hustotách sítě demonstruje celkovou konvergenci, což naznačuje, že další zjemňování má jen malý vliv na výsledky.
Studie výpočetní sítě musí být rovněž provedena podle následujícího odkazu:
Požadavek na přesnost WTG-Merkblatt M3
WTG-Merkblatt M3 poskytuje dvě klíčové metody pro validaci 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, za použití binárního klasifikačního přístupu (shoda nebo neshoda). Tento přístup hodnotí spolehlivost simulace výpočtem míry shody q, podobně jako spolehlivostní funkce používané v teorii spolehlivosti. Naproti tomu 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 zohlednila rozdíly v měřítku. Společně tyto metody poskytují jak kvalitativní, tak kvantitativní měřítka pro validaci simulace.
Výsledky a diskuse
Obrázek 3 ilustruje rozložení povrchového tlaku a integrované aerodynamické síly na dvourozměrné čtvercové desce. Simulace je provedena s rychlostí volného proudu 8,2 m/s a tlakové kontury se pohybují od +32,1 Pa do −45,4 Pa, což zvýrazňuje oblasti přetlaku a sání. Informační okno porovnává výsledky původní CAD geometrie s výsledky výpočetního modelu a uvádí výslednou sílu 4,8 N spolu s jejími složkami a odpovídajícím působištěm výsledných aerodynamických sil. Tento příklad demonstruje schopnost RWIND přesně vypočítat rozložení tlaku, integrovat aerodynamická zatížení a určit působiště výsledných aerodynamických sil pro účely validace a přenosu konstrukčního zatížení. Integrovaná aerodynamická zatížení získaná z výpočetního modelu vedou k celkové síle 4,8 N, včetně odporové síly 2,0 N, která je použita k výpočtu součinitele odporu CD = 0,073 na základě referenční plochy 0,648 m². Vypočtené působiště výsledných aerodynamických sil navíc poskytuje polohu výsledné aerodynamické síly působící na model. Společně tyto výsledky demonstrují konzistenci rozložení tlaku, integrovaného výpočtu sil a hodnocení aerodynamického součinitele použitého pro účely validace.
Tabulka 2 porovnává součinitele odporu získané z experimentálních měření a RWIND pro pět úhlů náběhu. Pro každý případ je uveden absolutní rozdíl (Pi − Oi), procentuální odchylka a soulad s akceptačními kritérii ±10 % a ±20 %. Výsledky ukazují vynikající shodu při −10° a 10°, s odchylkami 5,30 % a 2,56 %. Větší nesrovnalosti jsou pozorovány při −5° a 5°, zatímco případ 0° zůstává v rámci akceptačního kritéria ±20 %. Celkově 40 % případů splňuje kritérium ±10 %, zatímco 60 % splňuje kritérium ±20 %.
Tabulka 2: Porovnání součinitele odporu (Cp) mezi RWIND a experimentálními daty
| "Stupeň | Cp – Experiment (Oi) | Cp – RWIND (Pi) | Pi-Oi | Odchylka (%) | Míra shody ≤10 % | Míra shody ≤20 %" |
|---|---|---|---|---|---|---|
| -10 | 0,132 | 0,125 | -0,007 | 5,30 | 🟢 | 🟢 |
| -5 | 0,102 | 0,062 | -0,040 | 39,22 | 🔴 | 🔴 |
| 0 | 0,088 | 0,073 | -0,015 | 17,05 | 🔴 | 🟢 |
| 5 | 0,099 | 0,062 | -0,037 | 37,37 | 🔴 | 🔴 |
| 10 | 0,117 | 0,114 | -0,003 | 2,56 | 🟢 | 🟢 |
Tabulka 3 shrnuje statistické ukazatele použité k vyhodnocení shody mezi predikcemi RWIND a experimentálními referenčními daty. Validace je založena na 5 datových bodech, což vede k míře shody 40 % v rámci odchylky ±10 % a 60 % v rámci ±20 %. Celková chyba predikce je kvantifikována normalizovanou střední kvadratickou chybou (NMSE) 0,00093, střední chybou (ME) −0,0204, což ukazuje na mírné celkové podhodnocení RWINDem, střední absolutní chybou (MAE) 0,0204 a střední kvadratickou chybou (RMSE) 0,0254. Tyto metriky poskytují celkové hodnocení přesnosti predikce a vychýlení CFD modelu vzhledem k experimentálním měřením.
Tabulka 3. Statistické výkonnostní metriky validace RWIND
| Metrika | Hodnota |
|---|---|
| Počet datových bodů (N) | 5 |
| Míra shody (10 %) | 40 % |
| Míra shody (20 %) | 60 % |
| Normalizovaná střední kvadratická chyba, e² | 0,0551 |
| Střední chyba (ME) | -0,0204 |
| Střední absolutní chyba (MAE) | 0,0204 |
| Střední kvadratická chyba (RMSE) | 0,0254 |
Obrázek 4 porovnává součinitel odporu (CD) získaný z experimentálních měření a simulací RWIND pro úhly náběhu mezi −10° a +10°. Obě sady dat vykazují podobný trend ve tvaru U, s minimálním odporem při 0°. RWIND úspěšně zachycuje celkové aerodynamické chování, ale soustavně podhodnocuje součinitel odporu, zejména při středních úhlech (±5°), zatímco při extrémních úhlech (±10°) vykazuje dobrou shodu. Celkově výsledky ukazují, že RWIND spolehlivě predikuje kvalitativní trend změny odporu s úhlem náběhu, ačkoli další zlepšení v modelování turbulence a modelování u stěny by mohla zvýšit jeho kvantitativní přesnost.
