Экспериментальное исследование особенностей аэродинамического нагревания треугольного крыла при больших числах маха

Треугольные крылья часто являются несущим элементом высокоскоростного летательного аппарата. В связи с этим исследование теплопередачи и перехода ламинарного пограничного слоя в турбулентный на поверхности треугольных крыльев представляет большой практический интерес. В области турбулентного пограничного слоя величины тепловых потоков возрастают в два и более раз в сравнении с их ламинарным уровнем. В случае ламинарно-турбулентного перехода происходит существенная перестройка отрывного течения перед отклоненными органами управления, например, элевонами. Это в свою очередь приводит к изменению как тепловых, так и аэродинамических характеристик обтекаемого тела. Известен ряд работ [1–6] по изучению перехода на поверхности треугольных крыльев с острыми передними кромками. В некоторых работах [7, 8] исследовалось влияние радиуса затупления передних кромок на переход на наветренной стороне треугольного крыла. Было показано, что переход ламинарного пограничного слоя в турбулентный на треугольном крыле зависит от большого числа параметров. Основными из них являются число Маха набегающего потока, единичное число Рейнольдса, угол атаки, угол стреловидности и радиус затупления передних кромок крыла. Затупление вершины крыла также может вызывать образование локальных областей повышенных тепловых потоков [9]. Однако этих результатов и данных некоторых других работ недостаточно для понимания закономерностей, определяющих существование и положение областей повышенных величин тепловых потоков на поверхности крыла. В связи с этими проблемами и была выполнена настоящая экспериментальная работа.

Эксперименты проведены в аэродинамической трубе периодического действия Т-117 ЦАГИ [10]. В таблице 1 приведены основные параметры, которые были реализованы при проведении настоящих исследований.

На рис. 1 приведена схема модели, использованной в экспериментах при числе М го = 7, 5 с указанием координат сечений, выбранных для представления результатов исследований.

Модель выполнена в виде треугольного крыла с углами стреловидности передних кромок х = 75°, сферическим затуплением вершины (r = 6 мм) и цилиндрически- ми затуплениями кромок (r = 3 и 6 мм). В хвостовой части модели расположены два элевона, которые могут иметь углы отклонения 5 = 0, 10° и 20° относительно нижней поверхности крыла. Для иссле- дований при М^ = 10,5 была изготовлена модель треугольного крыла без элевонов. При этом длина модели стала равной 677 мм, а радиусы затупления обеих кромок и вершины сделаны равными 3 мм.

Таблица 1
№	а [ ° ]	5 1 [ ° ]	5 2 [ ° ]	M ~	P 0 [кг / см 2 ]	T о [К]	Re i • 10 - 6 [1 / м]	Re l • 10 - 6 [1 / м]
1	0	0	0	7,5	26,60	745	6,72	4,02
2	5	0	0	7,5	10,70	825	2,29	1,37
3	5	0	0	7,5	24,60	750	6,15	3,68
4	5	20	10	7,5	10,22	839	2,13	1,28
5	5	20	10	7,5	24,65	766	5,95	3,56
6	10	0	0	7,5	23,90	754	5,92	3,54
7	10	0	0	7,5	18,33	732	4,55	3,08
8	15	0	0	7,5	25,35	745	6,41	3,84
9	10	0	0	10,5	97,85	1100	4,8	3,25

№ — номер пуска трубы; α — угол атаки модели; δ 1 и δ 2 — углы отклонения правого и левого элевонов; M ^ — число Маха набегающего потока; P 0 и T 0 — давление и температура торможения; Re 1 — единичное число Рейнольдса, вычисленное по параметрам набегающего потока и длине, равной 1 м; Re L — число Рейнольдса, вычисленное по параметрам набегающего потока и длине крыла L

Для измерения величин тепловых потоков к поверхности крыла в данных испытаниях был применен метод термоиндикаторных покрытий [11]. В связи с этим крыло было изготовлено из малотеплопроводного материала АГ-4В. Модель, покрытая термоиндикатором с некоторой температурой плавления, после выхода трубы на режим быстро вводится в поток.

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

На телах со сферическим носовым затуплением измеренная величина теплового потока q относится обычно к величине теплового потока q о , вычисленного в критической точке носового затупления по формуле Фэя–Ридделла [12] для параметров набегающего потока в данном пуске трубы. Экспериментальные результаты представляются затем в виде набора изолиний q/q о на исследуемой поверхности модели и графиков распределения величин q/q о в заданных сечениях, что и было сделано в данной работе.

На рис. 2 показана картина плавления термоиндикатора на нижней поверхности крыла (кадр из видеофайла), полученная при угле атаки модели а = 5 ° и числах М го = 7 , 5 и Re L = 3 , 68 • 10 6 (пуск № 3). Углы отклонения элевонов 5 1 и 5 2 равны нулю.

Рис. 2

Границе плавления термоиндикатора на этом рисунке соответствует величина q/q о = 0 , 036. Внутри зоны плавления значения q/q о выше, а вне зоны — ниже, чем на границе. Видно, что в задней части крыла наблюдаются две больших клиновидных области плавления термоиндикатора (далее — турбулентные клинья), величины тепловых потоков в которых близки к соответствующим величинам в окрестности передних кромок. Образование этих областей связано с переходом ламинарного пограничного слоя в турбулентный на этих участках поверхности крыла. Заметно также слабое влияние разных величин радиуса затупления передних кромок на положение и величину областей перехода.

Рис. 3. М го = 7 , 5; а = 15 ° ; Re l = 3 , 84 • 10 6 ; δ 1 = δ 2 =0

Увеличение угла атаки крыла до а = 15 ° (пуск № 8, рис. 3) приводит к появлению в центральной части поверхности крыла ещё одной области перехода пограничного слоя, индуцированного возмущениями, приходящими от затупления вершины крыла.

То есть на этом режиме обтекания модели на большей части нижней поверхности крыла течение имеет переходный и турбулентный характер.

Об этом свидетельствует также график, приведённый на рис. 4.

Точка минимума на кривой, осредня-ющей распределение относительных величин тепловых потоков q/q о вдоль данного сечения, является точкой начала перехода пограничного слоя, а точка максимума — концом. В конце области перехода как в плоскости угла атаки ( у = 0), так и в боковых сечениях ( у = ± 52 мм), проходящих через турбулентные клинья, величины тепловых потоков выходят на один и тот же турбулентный уровень.

На рис. 5 показан график распределения величин q/q _о в сечении у = — 52 мм в зависимости от угла атаки крыла а .

С увеличением а точки начала и конца перехода ламинарного пограничного слоя в турбулентный смещаются вперед. Тепловые потоки в конце области перехода резко возрастают, градиент нарастания тепловых потоков в переходной области увеличивается, а длина области перехода уменьшается. При а = 5° и числе ReL ^ (1,28-1,37) • 106 переход на поверхности крыла не происходит и в данном сечении на рис. 5 наблюдается чисто ламинарное распределение тепловых потоков (кривая 2). При числе ReL = 3,68 • 106 в конце области перехода величина (q/qо)макс ~ 0,06 (кривая 3) и более чем в три раза превосходит со- ответствующее ламинарное значение теплового потока. При этом максимальная расчётная относительная величина теплового потока на кромке крыла с радиусом затупления r = 3 мм составляет q0 кромки/q0 = 0,197. В случае угла атаки а = 15° (кривая 5) в конце области перехода величина (q/q0)макс ~ 0,16, а расчётный тепловой поток на этой кромке крыла q0 кромки/qо = 0,28.

Рис. 4. М ^ = 7 , 5; а = 15 ° ; Re L = 3 , 84 • 10 6 ; 5 1 = 5 2 = 0 ; 1)-- у = 0 ; 2) V — у = — 52 мм ( r = 3 мм), △ — у = 52 мм ( r = 6 мм)

Рис. 5. М ^ = 7 , 5 ; 5 1 = 5 2 = 0 ; у = — 52 мм; 1) — — а = 0 , Re l = 4 , 02 • 10 6 ; 2) B — а = 5 ° , Re l = (1 , 28 - 1 , 37) • 10 6 ; 3) V — а = 5 ° , Re l = (3 , 56 - 3 , 68) • 10 6 ; 4) △ — а = 10 ° , Re l = 3 , 54 • 10 6 ; 5) M — а = 15 ° , Re l = 3 , 84 • 10 6

Рис. 6. M ^ = 7 , 5 ; а = 5 ° ; а) Re l = (1 , 28 - 1 , 37) - 10 6 ; 1) 5 1 = 5 2 = 0 , V — y = — 109 мм,-- y = 109 мм; 2) ! — 5 2 = 10 ° , y = — 109 мм;

3) 7 — 5 1 = 20 ° , y = 109 мм; б) Re l  =  (3 , 56 - 3 , 68) - 10 6 ;

1) 5 1 = 52 = 0, V — y = —109 мм,--y = 109 мм; 2) ! — 52 = 10°, y = —109 мм; 3) 7 — 5 1 = 20°, y = 109 мм

На рис. 6а, б приведены графики распределения относительных величин тепловых потоков, полученных при а = 5° в сечениях у = ± 109 мм, проходящих через середину элевонов (рис. 1) для случаев неот-клонённого элевона 5 1 = 52 = 0 и при отклонении его на 10° и 20° вниз. При малом числе Рейнольдса (рис. 6а) перед отклонённым на угол 20° элевоном существу- ет отрыв ламинарного пограничного слоя с резко выраженным пиком теплового потока в точке присоединения (x/xмакс ~ 1, 17) на элевоне. Величина (q/q0)макс ~ 0,43 в точке присоединения и примерно в 20 раз превышает величину q/q0 в этой точке на неотклонённом элевоне.

При числе ReL = 3,56 - 106 в области существования турбулентного клина на поверхности крыла области отрыва пограничного слоя перед отклонённым на 5 1 = 20° элевоном практически нет. На обоих графиках вертикальная линия на оси x/xмакс справа обозначает заднюю кромку элевона.

Распределения величин q/q ₀ на нижней поверхности крыла в поперечном сечении, расположенном на расстоянии x = 500 мм от его вершины, приведены на рис. 7.

При а = 5 ° и малом числе Re L в этом сечении, как уже отмечалось выше, существует ламинарный пограничный слой. При больших числах Рейнольдса мы видим образование областей переходного и турбулентного течения, которые при а = 15 ° распространяются от y/y макс — - 0 , 88 до y/y макс — 0 , 85. Координата y _макс соответствует точке пересечения поперечного сечения с образующими передних кромок крыла (рис. 1).

На рис. 8 показан график зависимости чисел Рейнольдса начала и конца перехода Re _t иRe _T , определённых в сечениях y = ± 52 мм, от угла атаки а . Числа Re t и Re T вычислялись по параметрам набегающего потока и расстояниям до точек начала x t и конца перехода x T соответственно.

Наиболее сильное уменьшение величин чисел Рейнольдса Re _t происходит в диапазоне а = 0—10 ° . При дальнейшем увеличении угла атаки до а = 15 ° величина числа Re t практически не изменяется. Минимальное значение числа Рейнольдса конца перехода реализуется при а = 15 ° и составляет Re T — 1 , 4 • 10 6 . Здесь же нанесены соответствующие данные, полученные в работе [5] при числе М ^ = 8 на крыле с острыми передними кромками. При нулевом угле атаки число Рейнольдса начала перехода на этом крыле примерно в 2 раза больше, чем на крыле с затупленными кромками. С увеличением α величины Re t , измеренные в обеих работах, сближаются. Длина области перехода в работе [5] меньше, чем на крыле с затупленными кромками. По-видимому, эти расхождения обусловлены различием в параметрах набегающего потока, геометрии кромок и положении измерительных сечений.

На рис. 9 показано распределение относительных величин тепловых потоков в сечениях y = 0 и y = — 52 мм для двух чисел Маха М ^ = 7 , 5 и М ^ = 10 , 5 и угла атаки крыла а = 10 ° . Число Re L в обоих случаях было приблизительно одинаковым.

y/y макс

Рис. 7. x = 500 мм; 1) а = 0 , Re L = 4 , 02 • 10 6 ; 2) ! — а = 5 ° , Re l = (1 , 28 - 1 , 37) • 10 6 ; 3) V — а = 5 ° , Re l = (3 , 56 - 3 , 68) • 10 6 ; 4) △ — а = 10 ° , Re l = 3 , 54 • 10 6 ; 5) M — а = 15 ° , Re l = 3 , 84 • 10 6

Рис. 8. 1) — — Re _t ( z = 52 мм, r = 3 мм), настоящие данные; 2) △ — Re t ( z = — 52 мм, r = 6 мм), настоящие данные; 3) — Re t ( z = 52 мм, r = 3 мм), 7 — Re t ( z = — 52 мм, r = 6 мм), настоящие данные; 4) V — Re t ( z = 21 , 5 мм, r = 0 , М ^ = 8 ) [5]; 5) ! — Re t ( z = 21 , 5 мм, r = 0 , М ^ = 8 ) [5]

Рис. 9. a = 10 ° ; 1)-- y = 0 ; М ^ = 7 , 5 ; Re l = 3 , 54 • 10 6 ; V — y = 0 ; М ^ = 10 , 5 ; Re l = 3 , 25 • 10 6 ; 2) — y = — 52 мм; М ^ = 7 , 5 ; Re l = 3 , 54 . 10 6 ; 3) — y = — 52 мм; М ^ = 10 , 5 ; Re l = 3 , 25 • 10 6

Рис. 10. а = 10 ° , x = 400 мм; а) 1)--М ^ = 7 , 5 ; Re l = 3 , 54 • 10 6 , r = 6 мм; 2) △ —

М ^ = 7 , 5 ; Re l = 3 , 08 • 10 6 , r = 3 мм; 3) V

М ^ = 10 , 5 ; Re l = 3 , 25 • 10 6 , r = 3 мм

В плоскости угла атаки крыла данные для двух вариантов обтекания практически совпадают. Это же наблюдается и в ламинарной области бокового сечения y = — 52 ( x/x макс <  0 , 2). Точка начала перехода в этом сечении с ростом числа М ^ смещается назад. При этом максимальная величина теплового потока в конце области перехода составляет здесь ( q/q о ) макс — 0 , 084.

В поперечном сечении x = 400 мм (рис. 10 а) при М ^ = 10 , 5 и радиусе затупления вершины крыла r = 3 мм появляются два ярко выраженных пика тепловых потоков вблизи плоскости симметрии крыла, индуцированных затуплением его вершины. Они хорошо видны на рис. 10 б в виде двух полос плавления термоиндикатора, идущих от окрестности вершины крыла. На этом рисунке приведён кадр видеозаписи, полученной в пуске трубы при М ^ = 10 , 5, Re L = 3 , 25 • 10 6 и радиусе затупления вершины крыла r = 3 мм. Границе изменения цвета термоиндикатора на этом рисунке соответствует величина q/q о = 0 , 034. Горизонтальной линией показано сечение x = 400 мм.

На рис. 10а при М^ = 10,5 (кривая 3) также видны небольшие максимумы q/qо, связанные с началом образования турбулентных клиньев в этой области поверхности модели, которые хорошо видны на рис. 10б. При числе М^ = 7,5 и r = 3 мм координаты точек (q/qо)макс в областях перехода не изменяются, то есть турбулентные клинья образуются в тех же местах поверхности крыла, что и при М^ = 10,5 (кривые 2 и 3). В центральной части крыла уже существует турбулентное течение. Увеличение радиуса затупления вершины крыла до r = 6 мм (кривая 1) приводит к тому, что пики q/qо в центральной части крыла отсутствуют, а области перехода более сильно развиты и расположены дальше от плоскости симметрии, чем при r = 3 мм.

При числах М ^ = 7 , 5 и 10 , 5 с помощью метода термоиндикаторных покрытий было выполнено исследование аэродинамического нагревания модели треугольного крыла с углом стреловидности передних кромок х = 75 ° . Целью исследований было как получение данных по величинам тепловых потоков на поверхности крыла, так и определение положения областей перехода ламинарного пограничного слоя в турбулентный.

Представленные результаты показали, что при а = 5°, числе М ^ = 7,5 и пониженных числах Rel — 1,3 • 106 на нижней поверхности крыла практически полностью реализовывалось ламинарное течение. При этом на элевоне с углом отклонения 20° наблюдается типичный для областей присоединения ламинарного пограничного слоя локальный пик теплового потока, максимальная величина которого примерно в 20 раз превышает величину теплового потока в этой точке на неоткло-нённом элевоне.

При числах Re L = (3 , 25-4 , 02) • 10 6 и числе М ^ = 7 , 5 распределение тепловых потоков на нижней поверхности модели имеет существенно неравномерный характер. Это связано с образованием здесь на всех исследованных углах атаки модели двух локальных областей ламинарно-турбулентного перехода (турбулентных клиньев) у каждой передней кромки. Турбулентные клинья, вызванные переходом вблизи передних кромок крыла, существуют при всех а . Увеличение угла атаки приводит к резкому возрастанию величин тепловых потоков в области турбулентных клиньев и смещению вперед точек начала и конца перехода. Величины ( q/q о ) _макс в конце области перехода приблизительно в 3 раза превышают соответствующее ламинарное значение в точке начала перехода.

При а = 5 ° перед большей частью поверхности отклонённых элевонов на крыле расположены турбулентные клинья. Судя по тому, что в них перед задней кромкой крыла не наблюдаются зоны малых величин тепловых потоков, можно сказать, что в областях переходного и турбулентного течений отрыва пограничного слоя перед элевонами практически нет. Максимальные величины тепловых потоков на элевоне с углом отклонения 5 2 = 10 ° приблизительно в 4 раза больше, чем на неот-клонённом, а при 5 1 = 20 ° увеличиваются примерно на порядок.

Наиболее сильное влияние увеличения угла атаки на уменьшение числа Re _t в сечениях у = ± 52 мм происходит в диапазоне а = 0-10 ° (приблизительно от Re t и 1 , 1 • 10 6 до Re t и 0 , 5 • 10 6 ). При дальнейшем увеличении угла атаки до а = 15 ° величина числа Re t практически не изменяется. Минимальное значение числа Рейнольдса конца перехода реализуется при а = 15 ° и составляет Re T и 1 , 4 • 10 6 .

В случае а = 10 ° на нижней поверхности крыла в его центральной части вблизи задней кромки появляется ещё одна область ламинарно-турбулентного перехода, индуцированная возмущениями от носового затупления. При а = 15° эта об- ласть развивается в ярко выраженный турбулентный клин, конец перехода пограничного слоя в котором расположен в точке с координатой xT/хмакс и 0,6.

Кроме того, при М ^ = 10 , 5 и а = 10 ° в окрестности плоскости симметрии крыла под воздействием возмущений, индуцированных затуплением вершины крыла, возникают резко выраженные локальные пики тепловых потоков.

Величина радиуса затупления вершины крыла имеет определяющее значение на положение областей перехода на поверхности крыла в поперечном направлении.

Работа выполнена при частичной финансовой поддержке ВЦП РН-ПВШ 2.1.1.5904.