Aerodynamické interakce: numerická a experimentální validace interakce vrtulí a ramen multikoptér

Přehled problému: proč interakce vrtulí a ramen určuje výkon multikoptér

U multikoptér proudění z vrtule interaguje s blízkou strukturou – zejména s rameny a kryty motorů. Vzniká tím komplexní kombinace efektů: lokální odtržení, zesílení vírů, periodické tlaky a změny indukované rychlosti, které přímo ovlivňují tah, účinnost a hluk. Vzhledem k tomu, že typická vzdálenost hrany ramene od konce vrtule je často < 0,2–0,35R (R = poloměr vrtule), jedná se o zónu silné aerodynamické interakce, kde lineární aproximace selhávají a je nezbytná kombinace CFD simulací a validace v aerodynamickém tunelu.

Základní bezrozměrné ukazatele a definice

• Koeficient tahu: CT = T / (u03C1 n^2 D^4), kde T je tah, u03C1 hustota, n otáčky [1/s], D průměr vrtule.
• Koeficient výkonu: CP = P / (u03C1 n^3 D^5).
• Indukované číslo: u03BB = vi / (n D) pro visení; v_i je indukovaná rychlost.
• Solídnost (solidity): u03C3 = B c / (u03C0 R) (B = počet listů, c = střední délka listu).
• Advance ratio (poměr postupu): J = V / (n D); pro visení je J ≈ 0, v reálném letu J ≈ 0,1–0,6.

Mechanismy interakce: jaké jevy sledovat

• Clipped inflow (omezený přítok): rameno narušuje přítok k základu listu; vzniká asymetrie zatížení a pokles lokálního Cl.
• Tip-vortex & arm impingement (náraz vírů na rameno): špičkové víry periodicky narážejí na hrany ramene → tlakové pulzace, hlukové tóny.
• Wake re-ingestion (znovunasávání proudu): šíření toku od horních vrtulí na spodní (koaxiální/stacked uspořádání) a zpětné nasávání kolem ramen.
• Bluff-body shedding (uvolňování vírů od tělesa): při hranaté geometrii ramene (Re ~ 10⁵–10⁶) vzniká Kármánova ulička, která interaguje se smykovými vrstvami listů.
• Swirl & kontra-rotace: rotující proud modifikuje úhel náběhu downstream vrtule/ramene; kontra-rotace redukuje netto víření, ale mění spektrum tlakových fluktuací.

CFD přístup: výběr modelu, domény a sítě

• Geometrie a varianty: parametrizace ramene (profil kruh/ovál/airfoil-like), hrany (radius/srážení), vzdálenost od roviny rotace (z/R), vysunutí ramene od osy (r/R), kryt motoru (nacelle), ochranný kruh (duct/guard).
• Rotace listů:
  • MRF (Multiple Reference Frame): rychlá stacionární aproximace pro screening.
  • Sliding Mesh / Overset: časově přesné řešení interakce s ramenem, zachycuje periodické děje (doporučeno pro validaci).
• Modely turbulence: k-ω SST pro RANS základ; DES/DDES nebo (S)LES pro zachycení špičkových vírů a interakčního spektra. Při (S)LES cílit na y⁺ ≈ 1–2 na listech a ramenech; při RANS postačí y⁺ ≈ 30–100 s wall functions, ale hrozí podmodeování odtržení.
• Doména a okrajové podmínky: hranice minimálně 5–10D od vrtule; inlet s uniformním nebo profilovaným proudem (pro dopředný let), outlet s p=0; stěny ramen a listů no-slip; rotace zadaná přes RPM.
• Diskretizace a časování: pro sliding mesh zvolit krok Δθ ≈ 1–3° na krok; konvergence hodnocena na periodické stacionaritě (CT, CP) přes 3–5 otáček.
• Síť: hexa/prism vrstvy na listech (20–30 proudových vrstev), lokální zjemnění kolem špiček a hran ramene; celkem 30–80 mil. buněk pro detailní LES jednoho rotoru s ramenem.

Metodika nastavení: krok za krokem

1. Izolovaný rotor (baseline): vytvořte CT–RPM–P mapu v MRF; validujte proti dostupným datasheetům nebo jednoduchému tahoměrnému stojanu.
2. Přidání ramene: nejprve stacionární MRF pro 5–10 variant trubky/profilu a polohy; vyberte 2–3 nejlepší řešení.
3. Časově přesný model: sliding mesh pro vybrané varianty; analyzujte spektrum tlaku na ramenech a listech (FFT).
4. Parametrická studie: z/R ∈ [−0,3; +0,3], r/R ∈ [0,6; 0,95], tloušťka ramene ∈ [0,03; 0,12], hrana = {ostrý, srážený, fillet R/t ∈ [0,2; 1,5]}.
5. Dopředný let: nastavte J = 0,2–0,4 a úhel náklonu 5–15°, sledujte posun vírů a crossflow na ramenech.

Validace v aerodynamickém tunelu: návrh experimentu

• Konfigurace: montáž rotoru a ramene v měřítku 1:1 nebo 1:2 na náboji s minimálním narušením držáku (sting/balance). Zajistěte stejná RPM a geometrii jako v CFD.
• Měření:
  • 6-komponentní váha (tah, krouticí moment, boční síly) → CT a CP.
  • Pitot/rake nebo PIV v rovině 0,1–0,3D pod rotory → profil indukovaných rychlostí.
  • Tlakové senzory na hranách ramene (rozteč 0,05–0,1D) → tlakové pulzace a dominantní frekvence.
  • Mikrofonové pole (akustika) → tonalita na BPF (blade-passing frequency) a jejích harmonikách.
• Protokol: pro každý variant 3 opakování, 5 úrovní RPM (např. 2000–5000 rpm), zaznamenávání teploty a tlaku (korekce hustoty).
• Korekce tunelu: korekce blokáže (< 5 % je ideál), interference s držákem; při otevřených zkušebních sekcích kontrolujte recirkulaci.

Porovnání CFD a tunelových dat: metriky a tolerance

• Primární: ΔC_T a ΔC_P oproti baseline rotor-only; cíl ±3–5 % pro RANS, ±2–3 % pro LES u průměrných hodnot.
• Spektra: porovnejte hlavní pík na BPF = B·RPM/60 a 1×–4× harmoniky; tolerance ±5 % ve frekvenci, ±2 dB v SPL pro relativní porovnání.
• Průvodiče/kinematika: PIV vs. CFD – RMSE indukované rychlosti v rovině pod rotorem < 10 % maxima.
• Nejistota měření: vyhodnoťte dle ASME PTC 19.1 – kombinovaná standardní nejistota cíleně < 1,5 % pro tah a < 2,5 % pro moment.

Výsledky typických studií: co se obvykle potvrdí

• Vliv vzdálenosti ramene od roviny rotace (z/R): posun ramene pod rovinu rotoru o 0,1–0,2R snižuje interakční ztráty CT o 1–3 procentní body ve srovnání s nulovou vzdáleností.
• Štíhlejší, zaoblená ramena: t/D ≤ 0,06 a fillet R/t ≥ 0,5 snižují tlakové pulzace na 0,6–0,8× oproti hranatým profilům.
• Fairing před hranou ramene (leading-edge cap): redukce CP o 2–4 % při zachování CT (tedy vyšší „figure of merit”).
• Kontra-rotace sousedních vrtulí: vyrovnává víření, zlepšuje rovnoměrnost indukce, ale může přidat sekundární píky v akustice – kompromis podle priorit (účinnost vs. hluk).

Akustika: interakční hluk a jak jej modelovat

• CAA workflow: CFD (časově přesné tlaky) → akustická analýza (Ffowcs Williams–Hawkings / integrační formulace) → mapy vyzařování.
• Dominantní zdroj: rozptyl špičkových vírů na hranách ramene; sekundární odtržení na ramenech při dopředném letu.
• Snížení hluku: přiblížení ramene k „tiché zóně” mezi listy (fázové sladění), zaoblení hran, mikro-zoubky na trailing edge vrtule (serrations) – testovat dle spektra.

Doporučení pro návrh ramene a rozmístění vrtulí

• Poloha ramene: z/R = −0,1 až −0,2 (rameno mírně pod rovinou rotoru); r/R ≈ 0,75–0,85, aby se minimalizovalo impingement špičkového víru.
• Profil ramene: eliptický nebo plytký airfoil-like (t/c ≈ 0,12–0,16), s plynulým náběhem; vyhněte se ostrým hranám a kvadratickým profilům.
• Fillet a přechody: štědré radiusy na napojení na držák motoru a trup; R/t ≥ 0,5.
• Guard/duct: pokud k dispozici, zajistěte dostatečný vůli (tip clearance) 1,5–3 % D a plynulý vstupní náběh; jinak hrozí pokles CT.
• Konfigurace rotace: sousední vrtule v kontra-rotaci; pro X-frame zvažte fázování listů tak, aby se BPF píky nesčítaly.

Příklad parametrického CFD porovnání (výstupní tabulky)

Kontrolní seznam pro CFD–tunelový program

• Definovat baseline rotor-only (CT, CP) a testovací RPM/J.
• Parametrizovat geometrii ramene a polohu (z/R, r/R, fillet).
• Spustit MRF screening → top 3 varianty.
• Pro top varianty sliding mesh s konvergencí za 3–5 otáček.
• Navrhnout tunelový test (váha, PIV/mikrofony, tlakové senzory) s korekcemi blokáže.
• Vyhodnotit ΔCT, ΔCP, spektrum @BPF, RMSE proudění; doplnit nejistoty.
• Iterovat geometrii (fillet, leading-edge cap, z/R) a potvrdit zlepšení ≥ 1–2 p.b. v CT/CP nebo −2 dB SPL.

Typické chyby a jak se jim vyhnout

• Malá CFD doména: příliš blízké hranice falešně tlumí víry → nadhodnocení CT.
• Nedostatečný časový krok: aliasing BPF → nepřesná spektra a pulzace momentu.
• Ignorování tunelových korekcí: bez korekcí blokáže/interference jsou ΔCT v rádu procent zavádějící.
• Nezkoordinované RPM a napájení: u elektrických pohonů se mění krouticí moment s teplotou; stabilizujte napětí a zaznamenávejte teplotu vinutí.

Rozšíření: dopředný let, zemní efekt a vícerotorová interakce

• Dopředný let (J>0): boční příkopy zvyšují asymetrii zatížení