Коли в'язка рідина тече по твердій поверхні, вона не може ковзати по стінці — шар рідини у безпосередньому контакті повинен відповідати швидкості поверхні (умова прилипання). Ця умова породжує тонку область — прикордонний шар, де швидкість переходить від нуля на стінці до зовнішнього значення U∞ і де в'язкі сили порівнянні з силами інерції. У 1908 році Пол Блаузіус аналітично розв'язав рівняння Нав'є–Стокса для ламінарного обтікання плоскої пластини через змінну подібності η = y√(U∞/νx), отримавши ОДР f‴ + ½ff″ = 0. Цей результат лежить в основі проектування крил літаків, розрахунків теплообміну та оцінки опору суден і підводних човнів.
Налаштовуйте швидкість зовнішнього потоку U∞ і довжину пластини L, перемикайтеся між водою та повітрям, щоб побачити вплив кінематичної в'язкості ν на товщину прикордонного шару δ₉₉ ≈ 5,0√(νx/U∞). Полотно показує криву δ(x), профілі швидкості у чотирьох перерізах і коефіцієнт тертя C̄f = 1,328/√Re_L. Усі показники оновлюються в реальному часі.
Що таке прикордонний шар Блаузіуса і хто його вивів?
Прикордонний шар Блаузіуса — точний аналітичний розв'язок рівнянь Нав'є–Стокса для ламінарного обтікання напівнескінченної плоскої пластини без градієнта тиску, виведений Полом Блаузіусом у 1908 році у його дисертації під керівництвом Прандтля. Він увів змінну подібності η, зводячи рівняння в частинних похідних до ОДР f‴ + ½ff″ = 0. Параметр стрільби f″(0) ≈ 0,4696 визначається чисельно і повністю характеризує профіль швидкості u/U∞ = f′(η).
Як товщина прикордонного шару δ зростає вздовж пластини?
99%-ва товщина шару δ₉₉ зростає як δ₉₉ ≈ 5,0√(νx/U∞) = 5,0x/√Re_x. Ця залежність від √x виникає з балансу між в'язкою дифузією імпульсу та конвекцією вниз по потоку. Подвоєння довжини пластини збільшує δ на виході у √2 ≈ 1,41 рази. Для повітря при 15 м/с над пластиною 0,5 м δ₉₉ на задній кромці становить близько 5 мм.
Що таке коефіцієнт тертя і як він розраховується?
Місцевий коефіцієнт тертя C_f(x) = τ_w/(½ρU∞²) = 0,664/√Re_x, де τ_w — напруга зсуву на стінці. Усереднений по пластині C̄_f = 1,328/√Re_L. Для плоскої пластини у повітрі при Re_L = 10⁶ C̄_f ≈ 1,33×10⁻³. Ця ламінарна формула перестає бути вірною вище Re_x ≈ 5×10⁵, де відбувається перехід до турбулентності.
Перехід відбувається, коли малі збурення підсилюються замість того, щоб загасати. Критичне Re_x,crit для плоскої пластини зазвичай становить 3×10⁵–10⁶ залежно від інтенсивності турбулентності зовнішнього потоку, шорсткості поверхні та градієнта тиску. Хвилі Толлміна–Шліхтінга зростають експоненційно і вироджуються у турбулентні плями. Перехід різко збільшує тертя шкіри, але також покращує теплообмін.
Товщина шару масштабується як δ ∝ √(ν/U∞). Кінематична в'язкість повітря при 20 °C ν ≈ 1,5×10⁻⁵ м²/с, тоді як у води ν ≈ 1×10⁻⁶ м²/с — приблизно у 15 разів менше. Тому прикордонний шар води у √15 ≈ 3,9 рази тонший за повітряний при тій самій швидкості. Це означає, що вода чинить більшу напругу зсуву на одиницю площі.
Товщина витиснення δ* ≈ 1,72√(νx/U∞) відображає, на скільки прикордонний шар «відсуває» зовнішні лінії течії від стінки. Для крила вона додається до аеродинамічної форми профілю: товщина витиснення 1 мм на хорді 200 мм — це ~0,5% ефективного збільшення товщини, що враховується при розрахунку профільного опору.
Несприятливий градієнт тиску (dP/dx > 0, потік гальмується) потовщує шар, зменшує імпульс поблизу стінки і може спричинити відрив, коли напруга зсуву на стінці обнуляється. Відрив відбувається на задній частині крил при великих кутах атаки і спричиняє зрив потоку — критичне явище в авіації. Сприятливий градієнт (прискорення) тонша і стабілізує шар.
Теорія прикордонного шару лежить в основі розрахунку опору і теплообміну практично будь-якого транспортного засобу. Крила літаків оптимізуються для відкладення переходу; субмарини контролюють шар для зниження акустичного шуму; лопатки газових турбін розраховуються за кореляціями теплообміну прикордонного шару; теплозахист космічних апаратів при вході в атмосферу розміровується за прогнозами аеродинамічного нагрівання.
Товщина імпульсу θ ≈ 0,664√(νx/U∞) пов'язана з опором через інтегральне рівняння Кармана: τ_w = ρU∞² dθ/dx. Повний опір однієї сторони пластини шириною b і довжиною L дорівнює ρU∞²bθ(L) — це елегантний результат, що дозволяє уникнути інтегрування місцевого тертя. Для профілю Блаузіуса θ = 0,133δ₉₉.
Тепловий прикордонний шар розвивається аналогічно шару швидкостей, з товщиною δ_T ≈ δ × Pr^{-1/3}, де Pr — число Прандтля. Для повітря Pr ≈ 0,71 (теплові та швидкісні шари майже однакові), для води Pr ≈ 7 (тепловий шар удвічі тонший). Місцеве число Нуссельта Nu_x = 0,332 Re_x^{1/2} Pr^{1/3} дає коефіцієнт теплопередачі безпосередньо з параметрів прикордонного шару.
На передній кромці розв'язок Блаузіуса передбачає сингулярність місцевого коефіцієнта тертя C_f ∝ x^{-1/2} → ∞. Фізично прикордонний шар на самому кінчику має нульову товщину — рідина вперше контактує зі стінкою саме тут. На практиці реальні пластини мають кінцеву товщину передньої кромки, де розв'язок Блаузіуса суворо не застосовний, але ця поправка незначна для більшості інженерних завдань.
Біля будь-якої твердої поверхні рухома рідина сповільнюється через тертя, утворюючи тонку область — примежовий шар, де швидкість зростає від нуля на стінці до повної швидкості набігаючого потоку. Ця сторінка відтворює розв'язок Блазіуса — точний аналітичний опис ламінарного примежового шару над плоскою пластиною без градієнта тиску. Спостерігайте, як шар потовщується далі за течією, як профіль швидкості вигинається назовні, і як зміна рідини змінює баланс між в'язкими та інерційними силами, який описує число Рейнольдса. Та сама фізика визначає опір крил літаків, корпусів суден і лопаток турбін.
Полотно малює межу примежового шару δ₉₉, яка зростає як √x уздовж пластини, чотири профілі швидкості u/U∞ у послідовних перерізах та частинки, що рухаються з місцевою швидкістю Блазіуса. Живі показники виводять ReL, δ₉₉, середній коефіцієнт тертя C̄f і напругу зсуву на стінці τw на задній кромці.
Перетягніть повзунок «Набігаючий потік U∞», щоб змінити швидкість течії, і «Довжина пластини L», щоб змінити відстань, яку долає рідина. Перемикайте кнопки «Вода» і «Повітря», щоб змінити кінематичну в'язкість ν, яка визначає, наскільки швидко зростає примежовий шар і скільки тертя він створює.
Людвіг Прандтль увів поняття примежового шару 1904 року; його учень Пол Блазіус точно розв'язав його для ламінарної течії 1908 року. Шкіра акул і ямки на м'ячах для гольфу навмисно використовують поведінку примежового шару, прискорюючи перехід до тоншого турбулентного шару, що знижує загальний опір.
В'язка дифузія розносить імпульс від стінки, тоді як рідина водночас рухається вниз за течією, що дає δ₉₉ ∝ √x. Тому подвоєння довжини пластини збільшує товщину на задній кромці лише приблизно у 1,4 раза, а не вдвічі.
Змодельований тут ламінарний профіль Блазіуса гладкий і самоподібний, описаний як u/U∞ = f′(η). Турбулентний профіль повніший біля стінки, ближчий до степеневого закону однієї сьомої, оскільки турбулентне перемішування переносить імпульс до поверхні набагато ефективніше, ніж сама лише в'язкість.
Кінематична в'язкість повітря приблизно у п'ятнадцять разів більша, ніж у води. Оскільки товщина масштабується як δ ∝ √ν, примежовий шар повітря за тієї самої швидкості набігаючого потоку виявляється приблизно у чотири рази товщим за водяний.
ReL = U∞L/ν порівнює інерційні та в'язкі сили. Низькі значення підтримують ламінарну течію, як і змодельовано тут, тоді як реальні течії зазвичай переходять до турбулентності, коли місцеве число Рейнольдса перевищує приблизно 500 000.
Умова прилипання змушує рідину біля пластини мати нульову швидкість відносно неї, створюючи крутий градієнт швидкості та напругу зсуву на стінці τw. Усереднення цієї напруги вздовж пластини дає коефіцієнт тертя C̄f = 1,328/√ReL.