10.10. ГИДРОДИНАМИЧЕСКОЕ СОПРОТИВЛЕНИЕ КРИВОЛИНЕЙНОЙ ПОВЕРХНОСТИ
Уравнение импульсов (9.5.10) можно записать в следующей форме:
(10.10.1)
Здесь
— текущее число Рейнольдса, построенное по толщине потери импульса;
— текущее числе Рейнольдса, построенное по длине контура
— формпараметр, характеризующий аэродинамическую кривизну потока;
— формпараметр, представляющий собой отношение толщин вытеснения и потери импульса;
— число Маха;
— плотность, кинематическая вязкость и скорость звука на внешней границе слоя в данном сечении
— координата, направленная вниз по потоку вдоль обвода контура; L — полная длина контура или другой его характерный размер (хорда, диаметр);
— относительное расстояние по обводу контура.
Полная длина контура рассчитывается от его передней кромки до задней или от точки разветвления потока до задней кромки тела.
Рассмотрим течение среды с постоянными физическими свойствами. Тогда коэффициент трения в общем случае будет функцией числа
и аэродинамической кривизны контура. Локальной характеристикой последнего фактора может служить формпараметр f, а интегральной — распределение скорости по контуру L, т. е. функция
(10.10.2)
Здесь
— характерная скорость, например скорость на бесконечности. В однопараметрическом приближении
Совместное решение уравнений (10.10.1) и (10.10.3) дает возможность рассчитать распределение основных параметров пограничного слоя вдоль контура.
Для ламинарного пограничного слоя имеются точные решения некоторых классов течения, характеризуемых видом функции (10.10.2), полученные Фолкнером и Скэн, Хоуартом, Гертлером и Виттингом, А. А. Дородницыным и др. Приближенные методы были предложены в работах Кармана и Польгаузена, Л. Г. Лойцянекого и др. Подробное изложение основных из этих методов дано в монографиях Л. Г. Лойцянского.
Здесь мы ограничимся приведением результирующей таблицы однопараметрического решения, полученного Н. Е. Кочиным и Л. Г. Лойцянским (табл. 10.4).
Таблица 10.4. Характеристики изотермического ламинарного пограничного слоя по приближенному однопараметрическому решению
Как видно, по этому решению отрыв изотермического ламинарного пограничного слоя происходит при значении формпараметра
(10.10.4)
а закон трения может быть выражен интерполяционной формулой
(10.10.5)
с погрешностью до 3%.
Для турбулентного пограничного слоя теоретическое определение всего сомплекса параметров отрыва впервые было сделано в однопараметрическом фиближении автором и А. И. Леонтьевым. Для изотермических условий по этому решению
(10.10.6)
зависимость (10.10.3) определяется рафиками рис. 10.6.
Экспериментальные данные Никуадзе и Фуруа хорошо подтверждают
решение в отношении величин
. Что касается значения
, достаточно точных эксперименальных данных пока не имеется и но, видимо, лежит в пределах
.
Из сопоставления формул (10.10.4) и (10.10.6) видно, что отрыв ламинарного пограничного слоя наступает при меньшем значении формпараметра
, чем отрыв слоя турбулентного. Кроме того, универсальной характеристикой грыва ламинарного пограничного слоя является произведение
, а турбулентного пограничного слоя — собственно формпараметр
.
Введем в рассмотрение функцию
десь
— формпараметр Бури—Лойцянскрго. Тогда уравнение импульсов можно записать в виде
.
(10.10.8)
Рис. 10.6. Закон трения в диффузорной области
При
; при
. Здесь
— относительное изменение коэффициента трения под влиянием неизотермичности при обтекании пластины. При
, и поэтому в области больших чисел Рейнольдса вполне приемлемо приближение
(10.10.9)
Вводя это выражение
в уравнение (10.10.8) и определяя закон трения степенной формулой, получаем
Интегрируя это уравнение для условий изотермического течения
, находим
(10.10.11)
Здесь
- безразмерная координата начала развития рассматриваемого пограничного слоя. Далее, по распределению
определяется и распределение коэффициента трения.
Более подробное изложение этой проблемы можно найти в монографиях, приведенных в списке литературы к данной главе.