Vliv modelování mezní vrstvy na přesnost aerodynamických CFD simulací

1. Úvod

V mnoha inženýrských CFD simulacích, zejména v aplikacích vnější aerodynamiky a větrného inženýrství, se turbulence v blízkosti stěny běžně modeluje pomocí standardních stěnových funkcí, spíše než přímým řešením viskózní podvrstvy. Typickým příkladem je funkce „nutkWallFunction“ používaná v OpenFOAM (a aplikovaná v programu RWIND), která odhaduje turbulentní viskozitu u stěny na základě logaritmického zákona stěny a empirických předpokladů modelování turbulence. Tato modelovací strategie umožňuje provádět simulace s relativně hrubými sítěmi, což vede k významnému snížení výpočetních nákladů a doby simulace, což je praktické pro rozsáhlé inženýrské problémy, jako jsou simulace budov, mostů a městského větru.

Zjednodušení zavedená standardními stěnovými funkcemi však také kladou významná omezení na přesnost předpovědi proudění blízko stěny. Tato omezení mohou ovlivnit veličiny, jako je smykové napětí na stěně, vývoj mezní vrstvy a lokální rozdělení tlaku, a proto je třeba je pečlivě zvážit při interpretaci aerodynamických výsledků, hodnocení přesnosti modelu nebo provádění validačních studií na základě experimentálních dat.

Tento článek pojednává o klíčových důsledcích a technických implikacích použití standardních přístupů stěnových funkcí v aerodynamických CFD simulacích za předpokladu velmi vysokých hodnot y⁺, které leží mimo doporučený nebo přijatelný rozsah. Zdůrazňuje také situace, ve kterých je tato modelovací strategie vhodná, stejně jako případy, kdy její omezení mohou ovlivnit spolehlivost výsledků.

2. Základní informace: stěnové funkce v aerodynamické CFD

V blízkosti pevných stěn obsahují turbulentní mezní vrstvy několik oblastí:

• Viskózní podvrstva (nízká 𝑦+)

• Vyrovnávací vrstva

• Logaritmická oblast (vysoká 𝑦+)

1 Struktura turbulentní mezní vrstvy ve stěnových jednotkách vykazující viskózní, vyrovnávací, logaritmické a vnější oblasti

Standardní stěnové funkce předpokládají:

• První výpočetní buňka leží v logaritmické oblasti (obvykle 𝑦+>30).

• Rychlostní profil blízko stěny se řídí logaritmickým zákonem stěny.

• Turbulentní viskozita a napětí ve smyku se odhadují pomocí empirických vztahů namísto řešení malých konstrukcí proudění.

Tato modelovací strategie se široce používá v

• externí aerodynamice,

• simulacích větrného inženýrství,

• simulacích velkých environmentálních toků,

protože umožňuje stabilní řešení s relativně hrubými sítěmi.

3. Výhody současného přístupu s vysokou hodnotou stěnové funkce y⁺

Použití standardních stěnových funkcí s vysokými hodnotami 𝑦+ nabízí několik praktických výhod pro inženýrské CFD simulace. Tím, že se vyhnete nutnosti řešit viskózní podvrstvu, umožňuje tento přístup použití hrubších sítí, což vede ke snížení výpočetních nákladů a zrychlení simulací, přičemž stále zachycuje celkové chování proudění relevantní pro mnoho velkých aerodynamických studií.

• Vysoká výpočetní účinnost

Umožňuje použití relativně hrubých sítí blízko stěn, což výrazně snižuje počet výpočetních buněk a dobu simulace.

• Snížené výpočetní náklady

Vyhýbá se extrémně jemné síti potřebné k vyřešení viskózní podvrstvy (𝑦+≈1), která by dramaticky zvýšila požadavky na CPU a paměť.

• Vhodné pro velké výpočetní domény

Účinné pro simulace zahrnující rozsáhlá prostředí, jako je městská zástavba, modely terénu a komplexní skupiny budov.

• Numerická robustnost a stabilita

Standardní přístupy založené na stěnových funkcích jsou obecně méně citlivé na nepravidelnosti sítě a mají tendenci produkovat stabilní řešení pro složité geometrie.

• Efektivní pro parametrické studie

Umožňuje rychlé porovnání více alternativních návrhů, orientací budov nebo podmínek prostředí.

• Vhodné pro rané fáze posouzení

Poskytuje rychlý aerodynamický přehled během fází koncepčního posouzení, kdy není vyžadováno podrobné rozlišení mezní vrstvy.

• Účinné pro zachycení globálního chování proudění

Přesně reprezentuje velkoplošné proudové vzorce, jako jsou úplavy, zóny zrychlení proudění a obecné rozložení větru.

• Praktický technický kompromis

Vyvažuje výpočetní účinnost a přijatelnou technickou přesnost, když je hlavním zájmem celkové aerodynamické chování spíše než detailní fyzika v oblasti blízko stěny.

4. Hlavní nevýhody v aerodynamických simulacích

4.1 Snížené rozlišení v blízkosti stěny

Standardní stěnové funkce neřeší viskózní podvrstvu. V důsledku toho:

• Nejsou zachyceny podrobné gradienty rychlosti v blízkosti stěny.

• Turbulence v blízkosti stěny jsou aproximovány, nikoli vypočítány.

• Lokální fyzikální vlastnosti proudění mohou být příliš zjednodušeny.

Toto omezení má přímý vliv na aerodynamickou přesnost, pokud převládají účinky vázané na stěny.

4.2 Nepřesné napětí ve smyku na stěně a plášťové tření

Stěnové funkce odhadují napětí ve smyku na stěně z empirických logaritmických vztahů. To může způsobit:

• Chyby v plášťovém tření,

• Sníženou přesnost v předpovědi celkového odporu pro aerodynamická tělesa,

• Nepřesný průběh napětí ve smyku podél ploch.

Pro analýzu aerodynamického výkonu mohou být tyto odchylky významné.

4.3 Omezení při predikci separace proudění

Stěnové funkce předpokládají rovnovážné chování mezní vrstvy, které se stává neplatným v následujících případech:

• Silné nepříznivé tlakové gradienty

• Rychlé zpomalení proudění

• Odtržení způsobené složitou geometrií

Mezi typické důsledky patří:

• Nesprávná predikce začátku separace proudění

• Chyby v umístění opětovného uchycení

• Nepřesná struktura proudění

Vzhledem k tomu, že separace silně ovlivňuje aerodynamické síly, představuje to významnou nevýhodu.

4.4 Snížená přesnost v oblastech proudění a recirkulace

Vzhledem k tomu, že turbulence v blízkosti stěny není plně vyřešena:

• Vývoj smykové vrstvy může být modelován nepřesně

• Recirkulační zóny mohou mít nesprávnou velikost

• Vírové struktury mohou být příliš rozptýlené

To má vliv na odhad aerodynamického zatížení a vizualizaci proudění.

4.5 Citlivost na umístění y⁺

Platnost stěnové funkce závisí na správném umístění první buňky sítě blízko stěny:

• Vysoká hodnota 𝑦+ → logaritmický zákon platný

• Nízká hodnota 𝑦+ → předpoklady neplatí

Střední hodnota 𝑦+ (pufrovací oblast) → zvyšuje se nejistota modelování

Špatné posouzení sítě může proto způsobit skryté chyby, i když se simulace jeví jako numericky stabilní.

4.6 Omezená schopnost pro přechodové proudění

Standardní stěnové funkce předpokládají plně turbulentní podmínky proudění. V důsledku toho:

• Přechod z laminárního do turbulentního proudění nelze zachytit přesně.

• Vývoj mezní vrstvy může být zkreslen.

To je důležité pro aerodynamické aplikace, které zahrnují:

• Nízkou intenzitu turbulence

• Hladké povrchy

• Proudění kolem aerodynamického profilu

4.7 Snížená přesnost pro lokální tlaková maxima

Empirické modelování stěn může vyhlazovat tlakové gradienty v blízkosti ploch, což vede k:

• Podhodnocení nebo nadhodnocení součinitelů tlaku pro maximální hodnoty

• Snížené přesnosti při posuzování lokálního zatížení

• Potenciálním chybám ve výpočtech zatížení pro statické výpočty

4.8 Závislost na empirických předpokladech logaritmického zákona

Stěnové funkce se opírají o zjednodušené empirické modely, které předpokládají:

• Plně vyvinuté turbulentní mezní vrstvy

• Plynulý vývoj proudění

• Mírné tlakové gradienty

Pokud se skutečné aerodynamické podmínky od těchto předpokladů odchylují, přesnost klesá.

5. Technický kompromis: přesnost vs. efektivita

Navzdory nevýhodám se standardní stěnové funkce nadále běžně používají v technickém CFD, protože poskytují několik důležitých praktických výhod:

• Vylepšená numerická stabilita

• Snížené požadavky na síť

• Nižší výpočetní náklady

• Větší přístupnost pro technické aplikace, což umožňuje provádět simulace v reálném čase a s omezenými zdroji

V praktických inženýrských pracovních postupech jsou simulace s extrémně vysokým rozlišením, které plně řeší viskózní podvrstvu, často výpočetně velmi nákladné a pro běžné konstrukční úkoly ne vždy proveditelné. Proto přístupy založené na stěnových funkcích představují nejen kompromis mezi přesností a náklady, ale také praktickou strategii modelování, která zpřístupňuje aerodynamické simulace pro každodenní inženýrské aplikace. Jak je také zdůrazněno v článku, vysoce detailní a výpočetně náročné simulace často ještě nejsou praktické v typických inženýrských projektech, kde jsou časové a výpočetní zdroje omezené.

Z těchto důvodů jsou standardní přístupy založené na stěnových funkcích často přijatelné pro:

• Vizualizaci proudění

• Studie v rané fázi posouzení

• Koncepční aerodynamické hodnocení

• Topografické nebo terénní studie větru

• Odhad globálního zatížení

Mohou však být nevhodné pro aplikace, kde je rozhodující podrobná fyzika proudění blízko stěn, jako například:

• Výpočty maximálního tlaku

• Posouzení lokálního tlaku na fasádu

• Nestacionární simulace zahrnující podrobné kolísání proudění

• Lokální účinky v blízkosti ostrých hran nebo rohů

• Přesná analýza separace a opětovného uchycení

Tabulka 1 klasifikuje aerodynamické aplikace CFD podle jejich vhodnosti při použití přístupu s vysokým standardem 𝑦+ stěnové funkce, kde není rozlišena viskózní podvrstva. Doporučují se aplikace zaměřené na globální chování proudění, některé analýzy požadující střední přesnost se doporučují podmíněně, zatímco studie závislé na podrobné fyzice blízko stěn se nedoporučují.

Tabulka 1: Použitelnost aerodynamických CFD simulací RWIND za předpokladu vysokého y⁺ standardní stěnové funkce

6. Závěr

Standardní přístupy založené na stěnových funkcích, jako je nutkWallFunction, představují praktickou a efektivní metodu pro modelování turbulentních mezních vrstev v aerodynamických simulacích. Nicméně závislost na empirických logaritmických aproximacích s sebou nese významná omezení, včetně snížené přesnosti fyziky blízko stěn, méně spolehlivé predikce separace a citlivosti na strukturu sítě. Inženýři by měli tyto nevýhody brát v úvahu při interpretaci výsledků simulací, provádění validace CFD nebo porovnávání simulací s experimentálními daty.