Теория пограничного слоя

Теория пограничного слоя — это раздел гидроаэродинамики, изучающий тонкий слой жидкости или газа, который образуется в непосредственной близости от поверхности обтекаемого твёрдого тела (например, крыла самолёта, лопасти турбины, корпуса корабля) и в котором скорость течения изменяется от нуля (на самой поверхности) до скорости внешнего потока. Вне этого слоя влиянием вязкости можно пренебречь, и течение считается потенциальным (невязким). Теория пограничного слоя позволяет рассчитывать сопротивление трения, теплообмен и диффузию между телом и потоком, что имеет ключевое значение в авиации, судостроении, метеорологии и химической технологии.

История развития

Предпосылки возникновения

До конца XIX века гидродинамика идеальной жидкости (уравнения Эйлера) не могла объяснить возникновение лобового сопротивления тел, движущихся в жидкости или газе: согласно теории, сопротивление должно было быть нулевым (парадокс Даламбера). Опыты показывали, что сопротивление существует и зависит от вязкости среды. В 1883 году Осборн Рейнольдс экспериментально установил критерий перехода от ламинарного течения к турбулентному (число Рейнольдса), но математического описания пристеночного слоя ещё не было.

Создание теории

В 1904 году немецкий учёный Людвиг Прандтль на III Международном математическом конгрессе в Гейдельберге впервые предложил концепцию пограничного слоя. Он показал, что при больших числах Рейнольдса влияние вязкости сосредоточено в тонком слое у поверхности тела, а вне этого слоя жидкость можно считать невязкой. Прандтль вывел упрощённые уравнения Навье — Стокса для этого слоя (уравнения Прандтля), которые значительно упростили расчёты. Его работа «О движении жидкости при очень малом трении» стала фундаментом современной прикладной гидроаэродинамики.

Развитие в XX веке

В 1920–1930-е годы ученики Прандтля — Теодор фон Карман, Генрих Блазиус, Пауль Кюхеман — разработали методы расчёта ламинарного и турбулентного пограничного слоя. Блазиус в 1908 году получил точное решение для ламинарного слоя на плоской пластине. В 1930-е годы появились полуэмпирические теории турбулентности (модель Прандтля — фон Кармана). Во второй половине XX века, с развитием вычислительной техники, теория пограничного слоя стала основой для численных методов расчёта обтекания (CFD — Computational Fluid Dynamics).

Физическая сущность

Основные свойства

Пограничный слой характеризуется резким градиентом скорости поперёк течения. На поверхности тела (при y = 0) выполняется условие прилипания: скорость жидкости равна нулю. По мере удаления от стенки скорость возрастает и на внешней границе слоя достигает скорости внешнего потока (U∞). Толщина пограничного слоя δ — это расстояние от стенки, на котором скорость составляет 99 % от U∞. Для ламинарного течения δ ~ √(νx/U∞), для турбулентного — δ ~ x^(4/5) (где ν — кинематическая вязкость, x — расстояние от передней кромки).

Ламинарный и турбулентный режимы

Пограничный слой может быть ламинарным или турбулентным. Ламинарный слой — упорядоченное течение, в котором частицы движутся параллельными слоями, перенос импульса происходит за счёт молекулярной вязкости. Турбулентный слой — хаотическое, вихревое течение с интенсивным перемешиванием, перенос импульса усиливается за счёт турбулентных пульсаций. Переход ламинарного течения в турбулентное происходит при критическом числе Рейнольдса Reкр ≈ 5×10^5 (для плоской пластины) и зависит от шероховатости поверхности, уровня турбулентности набегающего потока и градиента давления.

Толщина вытеснения и толщина потери импульса

Для расчёта сопротивления и теплообмена используются интегральные характеристики:

• Толщина вытеснения δ* — расстояние, на которое смещаются линии тока внешнего потока из-за замедления жидкости в пограничном слое.
• Толщина потери импульса θ — мера уменьшения количества движения в слое по сравнению с потенциальным течением.
• Формпараметр H = δ*/θ — отношение толщины вытеснения к толщине потери импульса, характеризующее форму профиля скорости (для ламинарного слоя H ≈ 2.6, для турбулентного H ≈ 1.3–1.4).

Уравнения пограничного слоя

Уравнения Прандтля

Для двумерного стационарного течения несжимаемой жидкости с постоянной вязкостью уравнения Прандтля имеют вид:

• ∂u/∂x + ∂v/∂y = 0 (уравнение неразрывности)
• u·∂u/∂x + v·∂u/∂y = – (1/ρ)·dp/dx + ν·∂²u/∂y²

Здесь u — продольная скорость, v — поперечная скорость, p — давление, ρ — плотность, ν — кинематическая вязкость. Предполагается, что давление поперёк слоя не меняется (∂p/∂y ≈ 0), а градиент давления dp/dx определяется внешним потенциальным течением. Уравнения Прандтля параболического типа, что позволяет решать их маршевыми методами (от передней кромки вниз по потоку).

Решение Блазиуса

Для ламинарного пограничного слоя на плоской пластине (dp/dx = 0) Генрих Блазиус в 1908 году получил автомодельное решение. Вводя переменную подобия η = y·√(U∞/(νx)), он свёл уравнения к обыкновенному дифференциальному уравнению третьего порядка. Решение даёт профиль скорости, касательное напряжение на стенке τw = 0.332·ρ·U∞²/√Re_x и коэффициент трения cf = 0.664/√Re_x.

Интегральные соотношения

Для инженерных расчётов широко используется интегральное соотношение Кармана (1921): dθ/dx + (2θ + δ*)·(1/U∞)·dU∞/dx = cf/2

Это уравнение связывает толщину потери импульса θ, толщину вытеснения δ*, коэффициент трения cf и градиент давления. Для его замыкания необходимы дополнительные эмпирические или полуэмпирические соотношения (например, закон стенки для турбулентного слоя).

Отрыв пограничного слоя

Механизм отрыва

При обтекании тел с положительным градиентом давления (диффузорная часть, задняя часть крыла) скорость в пограничном слое может упасть до нуля, а затем стать отрицательной (обратное течение). В точке отрыва касательное напряжение на стенке обращается в ноль. За точкой отрыва образуется зона возвратного течения и вихревая область, что резко увеличивает лобовое сопротивление (кризис обтекания). Отрыв — одна из главных причин срыва потока на крыле самолёта и потери подъёмной силы.

Управление отрывом

Для предотвращения или затягивания отрыва применяются:

• Вихрегенераторы — небольшие пластинки на поверхности, создающие вихри, которые подмешивают энергию из внешнего потока в пограничный слой.
• Отсос пограничного слоя — удаление замедленной жидкости через щели или перфорацию.
• Сдув пограничного слоя — подача струй воздуха или газа вдоль поверхности для ускорения течения.
• Использование профилей с ламинаризацией — специальная форма крыла, обеспечивающая длительный ламинарный участок.

Применение

Авиация и ракетная техника

Теория пограничного слоя лежит в основе расчёта аэродинамического сопротивления крыла, фюзеляжа, оперения. Для сверхзвуковых летательных аппаратов учитывается сжимаемость газа и образование ударных волн, которые могут взаимодействовать с пограничным слоем (взаимодействие ударной волны с пограничным слоем). Теплообмен в пограничном слое критичен для проектирования теплозащиты космических аппаратов при входе в атмосферу.

Судостроение

Расчёт сопротивления трения корпуса судна основан на теории турбулентного пограничного слоя. Для снижения сопротивления применяются специальные покрытия (например, имитирующие кожу акул — риблеты) и системы подачи пузырьков воздуха (воздушная смазка).

Энергетика

В турбомашинах (газовые и паровые турбины, компрессоры) пограничный слой на лопатках определяет потери энергии и эффективность. Отрыв пограничного слоя в диффузорах и на лопатках — одна из основных причин снижения КПД.

Метеорология и океанология

Пограничный слой атмосферы (планетарный пограничный слой) — нижний слой атмосферы (высотой от нескольких десятков метров до 1–2 км), в котором трение о поверхность Земли существенно влияет на ветер. Аналогично существует пограничный слой океана (экмановский слой), определяющий перенос тепла и импульса между океаном и атмосферой.

Современное состояние и численные методы

Модели турбулентности

Для расчёта турбулентного пограничного слоя используются различные модели замыкания:

• Модели с нулевыми уравнениями (модель Прандтля — длина пути смешения).
• Модели с одним уравнением (модель Спаларта — Аллмараса).
• Модели с двумя уравнениями (k-ε, k-ω, SST Ментера).
• Модели напряжений Рейнольдса (RSM) — наиболее точные, но требующие больших вычислительных ресурсов.

Прямое численное моделирование (DNS)

При малых числах Рейнольдса возможно прямое численное решение полных уравнений Навье — Стокса без дополнительных моделей турбулентности. DNS позволяет детально изучать структуру турбулентного пограничного слоя, но требует огромных вычислительных мощностей и пока неприменима для инженерных расчётов реальных объектов.

Метод LES (Large Eddy Simulation)

Моделирование крупных вихрей — компромисс между DNS и RANS (Reynolds-Averaged Navier-Stokes). Крупные вихри рассчитываются напрямую, а мелкомасштабные — моделируются. LES широко используется в аэродинамике и гидродинамике для задач с отрывом и взаимодействием.

Интересные факты

• Пограничный слой может быть «невидимым»: на крыле самолёта его толщина составляет от нескольких миллиметров до нескольких сантиметров, но именно он определяет 50–70 % общего лобового сопротивления на дозвуковых скоростях.
• В 1960-е годы советские учёные (Г. И. Петров, В. Я. Нейланд) разработали теорию взаимодействия пограничного слоя с ударными волнами, что позволило создать гиперзвуковые летательные аппараты (например, «Буран»).
• В природе пограничный слой играет роль в полёте птиц, плавании рыб и даже в движении крови по сосудам (в микрососудах толщина пограничного слоя сравнима с диаметром сосуда).

Источники

• Прандтль Л. «Гидроаэродинамика». — М.: Издательство иностранной литературы, 1951.
• Шлихтинг Г. «Теория пограничного слоя». — М.: Наука, 1974.
• Ландау Л. Д., Лифшиц Е. М. «Гидродинамика». — М.: Наука, 1986.
• Лойцянский Л. Г. «Механика жидкости и газа». — М.: Дрофа, 2003.
• Белов И. А., Кудрявцев Н. А. «Теплопередача и сопротивление пакетов труб». — Л.: Энергия, 1971.
• Cebeci T., Cousteix J. «Modeling and Computation of Boundary-Layer Flows». — Horizons Publishing, 2005.