Хвиля 53: перехід від ламінарного до турбулентного режиму, ісламська геометрія та граничний шар Блазіуса

Хвиля 53 одним поглядом

〜 Перехід від ламінарного до турбулентного режиму

У 1883 році Осборн Рейнольдс уводив тонку нитку фарбника в трубу і змінював швидкість течії. Нижче критичної швидкості фарбник утворював ідеально пряму нитку — ламінарну течію. Вище певного порогу нитка раптово розпадалася і перемішувалась по всьому перерізу труби — турбулентна течія. Безрозмірний параметр, що керує цим переходом, — число Рейнольдса:

Ламінарний профіль швидкості

Для стаціонарної, повністю розвиненої ламінарної течії в круглій трубі (течія Хагена — Пуазейля) профіль швидкості є точно параболічним:

u(r) = u_max (1 − r²/R²)

де r — радіальна відстань від осі труби, а R — радіус труби. Максимальна швидкість на осі вдвічі перевищує середню по перерізу: u_max = 2 Ū. Усі елементи рідини рухаються прямими, паралельними шляхами — змішування не відбувається.

Турбулентний профіль швидкості

У повністю турбулентній течії в трубі (Re ≫ 10 000) усереднений за часом профіль швидкості значно пласкіший. Емпіричний степеневий закон 1/n дає гарне наближення:

u(r)/u_c ≈ (1 − r/R)^1/n

де n ≈ 7 для Re ~ 10⁵. Поверх цього середнього профілю накладаються турбулентні флуктуації швидкості u′(t) і v′(t). Їхня середньоквадратична величина відносно середньої — це турбулентна інтенсивність I = u_rms/Ū, зазвичай 5–15% у розвиненій турбулентності.

Симуляція використовує частинки-трасери: кольорові потоки входять зліва. У ламінарному режимі їхні траєкторії залишаються прямими, а кольорові смуги ніколи не змішуються. Коли Re перетинає поріг переходу, флуктуації v-компоненти розсіюють частинки по каналу, відтворюючи змішування фарбника, яке спостерігав Рейнольдс у своєму оригінальному досліді.

✦ Ісламські геометричні візерунки

Ісламське геометричне мистецтво — яке зустрічається в мечетях, палацах і рукописах від Андалусії до Центральної Азії — будується з невеликого набору взаємопов'язаних зірчастих багатокутників, розеток і сполучних плиток. Візерунки демонструють ідеальну трансляційну і (локально) обертальну симетрію, не залишаючи ані пропусків, ані накладень.

Зірчасті багатокутники {n/k}

Правильний зірчастий багатокутник {n/k} будується шляхом розміщення n рівновіддалених точок на колі та з'єднання кожної k-ї точки. Цю нотацію запровадив Шлефлі. Ключові приклади:

• {6/2} — зірка Давида (гексаграма); gcd(6,2) = 2 утворює два окремі рівносторонні трикутники.
• {8/3} — восьмипроменева зірка (октаграма); gcd(8,3) = 1 утворює єдину замкнену 8-променеву зірку — поширену в османському кахельному мистецтві та візерунках Гіріх.
• {10/4} — десятипроменева зірка Гіріх; канонічний елемент системи плиток Гіріх, знайдений у святині Дарб-і Імам, Ісфахан (1453 р. н.е.). Решта форм-проміжків (бантик, п'ятикутник, видовжений шестикутник, ромб) заповнюють площину без накладень.
• {12/5} — дванадцятипроменева османська зірка; поширена в турецькій кахлі Ізнік і марокканській мозаїці зелідж.

Логіка мозаїки

Генератор розміщує центри зірок або на квадратній решітці (з кроком d), або на шестикутній решітці (dx = d, dy = d√3⁄2 зі зсувом кожного другого ряду). Для кожного центру малюються шляхи багатокутника {n/k} і, за бажанням, заповнюються. Менший правильний n-кутник у центрі (радіус ≈ 0.45 R) створює внутрішню зону заповнення зі зниженою непрозорістю, імітуючи контраст між тілом зірки та навколишніми сполучними плитками.

Доступно шість історичних колірних палітр: Ісламське золото (бурштин і охра), Альгамбра (лісовий зелений), Османський синій (кобальт і небесний), Марокканська теракота (обпалена сієна), Могольський білий (майже білий на чорному) та Неонова ніч (фіолетовий на майже чорному). Готовий візерунок можна експортувати як PNG без втрат.

〰 Граничний шар Блазіуса

Коли в'язка рідина обтікає тверду поверхню, умова прилипання змушує швидкість рідини на стінці дорівнювати нулю. Розвивається тонкий шар — граничний шар — у якому швидкість зростає від нуля до значення зовнішнього потоку U∞. Поза цим шаром течія практично невязка.

Рівняння Блазіуса

У 1907–1908 роках Генріх Блазіус (учень Прандтля) показав, що для стаціонарної ламінарної течії над напівнескінченною плоскою пластиною рівняння Нав'є — Стокса зводяться, за наближеннями граничного шару, до єдиного нелінійного ОДР третього порядку через змінну подібності η = y √(U∞ / νx):

f‴ + ½ f f″ = 0   за умов   f(0) = f′(0) = 0, f′(∞) = 1

Тут f′(η) = u/U∞ — безрозмірний профіль швидкості. Розв'язок універсальний — усі ламінарні граничні шари над плоскою пластиною збігаються на цій єдиній кривій, якщо їх побудувати як функцію η.

Метод стрільби та РК4

Гранична умова f′(∞) = 1 не задана при η = 0, тому пряме інтегрування неможливе. Натомість використовується метод стрільби: задається початкове значення α = f″(0) ≈ 0.4696 (відома константа пристінного зсуву Блазіуса), інтегрування проводиться вперед методом РК4 з кроком Δη = 0.02, і перевіряється, що f′ → 1 при η → 12. Отримана таблиця розв'язку використовується під час виконання через лінійну інтерполяцію для обчислення u/U∞ у будь-якій точці (x, y) на полотні.

Ключові результати

• Товщина граничного шару: δ₉₉(x) ≈ 5.0 √(νx/U∞)
• Місцеве поверхневе тертя: C_f(x) = 0.664 / √Re_x
• Середнє поверхневе тертя по довжині пластини L: C̄_f = 1.328 / √Re_L
• Дотичне напруження на стінці: τ_w(x) = 0.332 ρ U∞² / √Re_x

На полотні знизу показана плоска пластина з оболонкою граничного шару δ₉₉(x) (помаранчева заливка), профілями швидкості u/U∞ у чотирьох точках x (білі криві) та стрілками зовнішнього потоку над шаром. Перемикання між водою (ν = 1 × 10⁻⁶ м²/с) і повітрям (ν = 1.5 × 10⁻⁵ м²/с) миттєво оновлює товщину граничного шару та значення поверхневого тертя.

Питання чи пропозиції? Звертайтеся через сторінку контактів.