close

Вход

Забыли?

вход по аккаунту

?

Применение тангенциального выдува струй для снижения сопротивления сверхкритических профилей при больших дозвуковых скоростях..pdf

код для вставкиСкачать
УЧЕНЫЕ
Том XXXX
ЗАПИСКИ
ЦАГИ
2009
№1
УДК 629.735.33.015.3.062.4
ПРИМЕНЕНИЕ ТАНГЕНЦИАЛЬНОГО ВЫДУВА СТРУЙ
ДЛЯ СНИЖЕНИЯ СОПРОТИВЛЕНИЯ СВЕРХКРИТИЧЕСКИХ
ПРОФИЛЕЙ ПРИ БОЛЬШИХ ДОЗВУКОВЫХ СКОРОСТЯХ
В. Д. БОКСЕР, А. В. ВОЛКОВ, А. В. ПЕТРОВ
Рассматривается проблема создания высокоскоростных сверхкритических крыловых
профилей с пониженным уровнем сопротивления при больших дозвуковых скоростях. Предлагаемая концепция заключается в использовании тангенциального выдува высоконапорной
струи малой интенсивности для ликвидации волнового отрыва, обусловленного скачком уплотнения, замыкающим местную сверхзвуковую зону на профиле. Выполнены расчетные
оценки возможности уменьшения сопротивления сверхкритических крыловых профилей второго поколения с низким уровнем волнового сопротивления в условиях турбулентного обтекания путем тангенциального выдува струй на верхней и нижней поверхностях при значени6
ях су = 0.55 и Re = 40 · 10 в диапазоне чисел М = 0.7 ÷ 0.88.
Ключевые слова: профиль, тангенциальный выдув, распределение давления, сопротивление.
Создание сверхкритических крыловых профилей [1] позволило, с одной стороны, существенно продвинуться по скорости по сравнению с обычными профилями без снижения аэродинамического качества. С другой стороны, использование сверхкритических профилей дало возможность значительно увеличить их максимальную относительную толщину при фиксированном
числе Маха, что позволило разработать утолщенные крылья большого удлинения (λ = 9 — 11) и
малой стреловидности (χ1/4 = 25 — 30°) [2], успешно реализованные в компоновках современных
отечественных пассажирских и транспортных околозвуковых самолетов (Ту-204, Ил-96, Ан-124,
Ан-70). Проведенные в ЦАГИ подробные физические исследования особенностей околозвукового обтекания сверхкритических профилей [3], развития местной сверхзвуковой зоны на их поверхностях [4], возникновения и развития отрыва потока на режимах М > Мкр [5] показывают,
что дальнейшее продвижение по числам Маха сопровождается существенным ростом сопротивления, обусловленного отрывом из-под скачка уплотнения (волновым отрывом).
Одним из возможных путей снижения сопротивления профилей при больших околозвуковых скоростях (М > 0.8) является использование тангенциального выдува струй сжатого воздуха
с целью ликвидации волнового отрыва.
Идея тангенциального выдува струй для ликвидации диффузорных отрывов на крыле и закрылках при малых скоростях (М = 0.15 — 0.3) известна [6]. Созданы полуэмпирические методы
расчета коэффициента импульса струи, потребного для ликвидации отрыва при малых скоростях
[7 — 10]. Из литературы известен пример экспериментального исследования тангенциального
выдува струи на тонком симметричном профиле RAE-102 ( c = 6%) при числе М = 0.75 [11].
На перспективных пассажирских и транспортных околозвуковых самолетах с утолщенными
сверхкритическими стреловидными крыльями ( c = 16 — 13 — 11% и более) обеспечение экономичного крейсерского полета с высоким уровнем аэродинамического качества (Kmax ≈ 18 — 20)
при числах Мкрейс = 0.85 — 0.9 представляет сложную проблему вследствие возникновения интенсивных волновых отрывов на крыле.
8
Целью настоящей работы являются
предварительные расчетные оценки возможности
снижения
сопротивления
сверхкритических профилей большой относительной толщины ( c > 10%) при околозвуковых скоростях (М > 0.8) за счет
подавления волнового отрыва путем использования тангенциального выдува
струй.
Постановка задачи. Метод расчета. Рассматривается обтекание крылового
профиля вязким сжимаемым потоком газа
при наличии развитой местной сверхзвуковой зоны, замыкаемой скачком уплотнения. Вблизи точки отрыва потока X отр ,
Рис. 1. Расчетная схема обтекания профиля с тангенциальным
выдувом струи
вызванного скачком уплотнения, осуществляется выдув из щелевого сопла струи сжатого воздуха с интенсивностью, достаточной для
перемещения точки отрыва к задней кромке крыла (рис. 1).
Для расчета вязкого отрывного околозвукового обтекания профиля с интенсивным скачком
уплотнения использована программа VISTRAN, (А. В. Волков, С. В. Ляпунов, 1994 г.), базирующаяся на зональном методе [12]. В данном методе область невязкого течения рассчитывается
с помощью уравнения полного потенциала скорости, а турбулентный пограничный слой — с помощью метода Грина. Согласование решений во внешней невязкой и внутренней вязкой области
1 d (ρuδ*)
, устанавливающего связь
течения производится посредством условия Лайтхилла: v =
ρ dS
между нормальной скоростью на поверхности профиля v и продольной скоростью u c величиной
толщины вытеснения пограничного слоя δ*, получаемой из решения вязкой задачи.
Расчет турбулентного пограничного слоя на поверхности профиля осуществляется методом
Грина путем решения системы обыкновенных дифференциальных уравнений:
dθ/dS = f(θ, H, u), dH/dS = g(θ, H, u),
где θ — толщина потери импульса, Н = δ*/θ — формпараметр пограничного слоя, S — направление вдоль контура профиля.
Настоящая программа позволяет рассчитывать режимы вязкого отрывного околозвукового
обтекания профилей с интенсивными скачками уплотнения. Пример расчета такого режима
в сравнении с экспериментом [13] представлен на рис. 2.
Для приближенного учета влияния струи на обтекание
профиля использован полуэмпирический метод [7], согласно
которому коэффициент импульса струи, потребный для ликвидации отрыва пограничного слоя, может быть определен по
формуле:
cμ R =
2Δθ
⎛ Vср ⎞
0.85 ⎜1 −
⎟
Vc ⎠
⎝
2
,
(1)
где Δθ = Δθ b — относительная разность толщины потери
импульса на задней кромке профиля ( θ1 ) и в точке отрыва пограничного слоя ( θ0 ) при отсутствии выдува струи (см. Рис. 2. Сравнение расчетного и экспериментального распределения давления на
рис. 1); Vср — средняя скорость на обдуваемой части профи- профиле
NACA 64A010 при значениях
ля; Vс — скорость истечения струи из щелевого сопла
( Vс ≈ 3Vср ); b — хорда профиля.
M = 0.8, α = 6.2°, Re = 2 · 106:
◦ ◦ ◦ — эксперимент [13] сy = 0.62;
—⎯ — расчет сy = 0.618
9
Формула (1) позволяет получить удовлетворительное согласование расчетных и экспериментальных значений коэффициента cμ R для профилей с тангенциальным выдувом до- и сверхзвуковых струй на верхнюю поверхность отклоненных закрылков при больших углах их отклонения ( δз ≈ 50D ), т. е. в условиях развитого отрыва пограничного слоя в дозвуковом потоке [7, 8].
Данное соотношение не требует детального расчета обтекания профиля с тангенциальным выдувом, а лишь интегрально учитывает величину импульса струи, потребного для восстановления
безотрывного обтекания, путем уменьшения толщины потери импульса пограничного слоя в точке выдува струи θ0 на определенную величину η.
С целью учета тангенциального выдува струи программа VISTRAN была модифицирована
путем введения следующих параметров: θs = θ0 η и H s = H 0 η , где s — точка отрыва пограничного слоя; η — переменный параметр; θ0 , H0 — толщина потери импульса и формпараметр в точке
отрыва при отсутствии выдува струи.
Условиями выбора параметра η являются δ1∗ = δ∗0 , θ1 = θ0 , c f 1 ≈ 0, где δ1∗ , θ1 — соответственно толщины вытеснения и потери импульса пограничного слоя, c f 1 — коэффициент трения
на задней кромке профиля при наличии выдува.
В процессе расчета из полного сопротивления (сх) выделялось волновое сопротивление
( cxв ) профиля, рассчитываемое согласно [14]:
cхв ≈
А(М) В(М ск − 1)
,
kc(М cк − 1)
(2),
где функции А(М), В(М ск − 1) , c(М cк − 1) приведены в [14]; М cк — число Маха перед скачком
уплотнения на поверхности; k — кривизна поверхности у скачка уплотнения.
Рис. 3. Сравнение расчетного и экспериментального распределения давления
на профиле RAE 102 с тангенциальным выдувом струи при значениях
M = 0.75, α = 6°, Re ≈ 1.5 · 106
10
Пример согласования результатов расчетов с экспериментальными данными [11] приведен
на рис. 3. Удовлетворительное согласование расчетного и экспериментального распределения
давления по профилю как при отсутствии выдува, так и с выдувом свидетельствует о возможности применения предложенного приближенного метода для качественного анализа воздействия
тангенциального выдува на аэродинамические характеристики профилей при больших дозвуковых скоростях.
Анализ результатов расчета. Расчетные исследования по влиянию тангенциального выдува струй на сопротивление сверхкритических крыловых профилей с максимальными относительными толщинами c = 10 — 14% проведены в условиях турбулентного обтекания при фиксированном крейсерском значении коэффициента подъемной силы су = 0.55 в диапазоне чисел
6
М = 0.7 — 0.88 и постоянном натурном числе Re = 40 · 10 .
На рис. 4 показаны форма толстого сверхкритического профиля ( c = 14%) с уплощенной
верхней и искривленной нижней поверхностями и распределение давления при его докритическом (М = 0.6) и закритическом (М = 0.83 ) обтекании. Для такого типа профилей характерной
является «подрезка» в хвостовой части нижней поверхности, компенсирующая потерю подъемной силы, вызванную уплощением верхней поверхности [1, 2]. Как видно из рис. 4, переход
от числа М = 0.6 к числу М = 0.83 сопровождается возникновением местных сверхзвуковых зон
на обеих поверхностях. При этом волновое сопротивление ( cxв ) данного профиля ( c = 14%) и
сверхкритических профилей других максимальных относительных толщин ( c = 10 — 12%) малó
по сравнению с полным сопротивлением cx (рис. 5). Так, например, при числе М = 0.85 рассматриваемые профили с относительными толщинами c = 10 — 14% имеют cxв ≈ 0.005 — 0.015 при
cx ≈ 0.028 — 0.064.
Рис. 4. Распределение давления по сверхкритическому профилю с максимальной относительной толщиной c = 14% при докритическом
(M = 0.6) и закритическом (M = 0.83) обтекании
11
Рис. 5. Изменение коэффициентов полного (сx) и волнового (схв) сопротивления по числам Маха
для сверхкритических профилей c различными максимальными относительными толщинами
На рис. 6 показано изменение положения точки отрыва пограничного слоя ( xотр ) на верхней и нижней поверхностях исследуемых профилей по числам Маха. Зависимости xотр (М) носят сложный характер, обусловленный воз-
Рис. 6. Изменение положения точки отрыва
пограничного слоя по числам Маха для
сверхкритических профилей с различными
максимальными относительными толщинами
аэродинамическое качество K эф =
cx
никновением скачков уплотнения и их перемещением
вдоль верхней и нижней поверхностей профиля с увеличением числа М.
С целью ликвидации волнового отрыва в данной работе предложено использовать тангенциальный выдув тонкой высокоскоростной струи (до Мс = 2.7) в окрестности
точки отрыва пограничного слоя, возникающего на профиле при отсутствии выдува.
Влияние тангенциального выдува на подавление волнового отрыва на профилях большой относительной толщины ( c = 14%) при фиксированном угле атаки α = 3.22º и
числе М = 0.83 показано на рис. 7, а. Видно существенное
(примерно вдвое) увеличение подъемной силы (за счет
расширения местной сверхзвуковой зоны на верхней поверхности и увеличения давления в хвостовой части нижней поверхности), сопровождающееся ростом сопротивления профиля (примерно на 35%). При этом эффективное
cу
возрастает приблизительно на 30%.
+ cμ R
На рис. 7, б показано, что и при фиксированном значении c у = 0.55 тангенциальный выдув
струи непосредственно перед точкой отрыва исходного течения приводит к полной ликвидации
отрыва, смещению скачка уплотнения вниз по потоку и восстановлению давления ( c p1 ≈ 0.25) на
уплощенной верхней поверхности. На нижней поверхности при этом происходит существенное
восстановление давления в хвостовой части ( x > 0.65) за счет ликвидации отрыва
( xотр. н. п = 0.64 — без выдува струи, xотр = 1 — при наличии выдува струи в точке x = 0.64).
12
а)
б)
Рис. 7. Влияние тангенциального выдува струи на распределение давления по сверхкритическому
профилю c = 14% при фиксированном угле атаки (а) и при фиксированном значении коэффициента
подъемной силы (б)
Сравнение зависимостей cx(M), приведенных на рис. 8, показывает, что тангенциальный
выдув уменьшает полное сопротивление профиля с относительной толщиной c = 14% во всем
исследованном диапазоне чисел М = 0.7 — 0.87. Учет величины cμ R , рассчитанной по формуле (1),
несколько снижает общий положительный эффект от выдува струи (см. зависимости cx + cμ R
и cμ R на рис. 8). Отметим, что волновое сопротивление, рассчитанное по зависимости (2), при
выдуве струи возрастает (см. рис. 8) за счет расширения местной сверхзвуковой зоны на верхней
поверхности (см. рис. 7, б). Так, например, при числе М = 0.83 прирост волнового сопротивления
за счет ликвидации отрыва составил Δ cxв ≈ 0.0087 или 87% от уровня cxв при отрывном обтекании ( cxв отр ≈ 0.01). Факт снижения темпа или прекращения роста волнового сопротивления при
числах М > 0.84 на исследованных профилях большой относительной толщины ( c ≥ 12%) обусловлен повышенной кривизной хвостового участка верхней поверхности, где располагаются
скачки уплотнения при полностью безотрывном обтекании, поскольку величина cxв ∼ 1/k (см. (2)).
13
Рис. 8. Влияние тангенциального выдува струи на полное (сx) и волновое (схв) сопротивление
профиля при различных числах Маха
Рис. 9. Изменение относительного выигрыша в полном сопротивлении сверхкритических профилей с различными максимальными относительными толщинами за счет тангенциального
выдува струй с ростом числа Маха
cxвыд
( M ) , где cxвыд — коэффициент сопротивления
cx
профиля при наличии выдува, наглядно демонстрируют количественный эффект от тангенциального выдува для исследуемых профилей. Расчеты указывают на наличие для каждой максимальПредставленные на рис. 9 зависимости
ной относительной толщины профиля характерного числа Маха ( М∗выд ), при котором эффективность выдува максимальна. Как видно из рис. 9, снижение сопротивления на режиме М∗выд для
рассматриваемых профилей ( c = 10 — 14% ) достигает 18 — 35% (М = 0.83 — 0.78). При дальнейшем увеличении числа Маха (М > М∗выд ) эффективность выдува снижается вследствие заметного увеличения при этом волнового сопротивления (см. рис. 8).
14
Рис. 10. Изменение «вязкой» компоненты полного сопротивления ( сх − схв ) по числам Маха
для сверхкритических профилей с различными максимальными относительными толщинами
при наличии и отсутствии тангенциального выдува струи
Рис. 11. Изменение относительного вклада ( cxотр c x ) сопротивления, обусловленного отрывом
пограничного слоя ( cxотр ), по числам Маха для сверхкритических профилей с различными
максимальными относительными толщинами
Поскольку основным назначением тангенциального выдува струи является ликвидация отрыва при околозвуковом обтекании профилей, представляют интерес зависимости «вязкой» компоненты сопротивления ( сх − схв ) от числа М, приведенные на рис. 10. Видно заметное снижение
величины ( сх − схв ) за счет выдува струи вследствие ликвидации сопротивления, обусловленного
отрывом потока ( cxотр ).
На рис. 11 показан относительный вклад сопротивления cxотр cx , обусловленного отрывом
пограничного слоя ( cxотр ) и представляющего собой разность величин ( сх − схв ) при отсутствии и
наличии выдува струи, в полное сопротивление профиля для исследованных профилей в диапазоне чисел М = 0.7 — 0.88. Видно, что максимальные величины cxотр max ≈ 28 — 43% достигаются
15
при числах М = 0.83 — 0.8. Заметное снижение величины cxотр при числах М > 0.82 у исследованных профилей обусловлено, по всей вероятности, существенным уменьшением протяженности области отрыва на верхней поверхности при этих режимах обтекания.
Заключение. Проведенные расчетные оценки показали, что применение тангенциального
выдува высокоскоростных струй на сверхкритических крыловых профилях средних и больших
максимальных относительных толщин ( c = 10 — 14%) может обеспечить ликвидацию волнового
отрыва при достаточно малой интенсивности выдува и, как следствие, снижение сопротивления
примерно на 20 — 30% при больших дозвуковых скоростях (М ≥ 0.8). Указанные величины, полученные с использованием приближенной расчетной модели, следует рассматривать как предварительные. Эти результаты должны быть уточнены с применением более совершенных численных методов и подтверждены экспериментально.
ЛИТЕРАТУРА
1. С е р е б р и й с к и й Я. М., Н и к о л а е в а К. С., Б о к с е р В. Д. Новые типы сверхкритических околозвуковых профилей // ТВФ. 1970. № 5.
2. П а в л о в е ц Г. А., Б о к с е р В. Д., Л я п у н о в С. В. Аэродинамика крыловых
профилей. — В кн.: Аэродинамика и динамика полета магистральных самолетов / Под ред.
Г. С. Бюшгенса. — Москва — Пекин: Изд. ЦАГИ, Авиаиздат. КНР, 1995.
3. Б о к с е р В. Д. Некоторые особенности околозвукового обтекания профилей // Ученые записки ЦАГИ. 1980. Т. XI, № 2.
4. Б о к с е р В. Д. Экспериментальное исследование высоты местной сверхзвуковой
зоны и волнового сопротивления при околозвуковом обтекании профиля // Ученые записки
ЦАГИ. 1981. Т. XII, № 6.
5. Б о к с е р В. Д. Развитие отрыва и его влияние на аэродинамику сверхкритических
профилей при околозвуковых скоростях // Ученые записки ЦАГИ, 1988. Т. XIX, № 5.
6. L a c h m a n n G. V. Boundary layer and flow control. — Pergamon, London, 1961.
7. T h o m a s F. Untersuchungen uber die Erhohung des Auftriebes von Tragflugen mittels
Grenzschichtbeeinflussing durch Ausblasen // Z. Flugwiss, V. 10, № 2. 1962.
8. П е т р о в А. В., Ш е л о м о в с к а я В. В. Метод расчета коэффициента импульса
струи, потребного для ликвидации отрыва потока на профиле крыла // Труды ЦАГИ, 1979,
вып. 1977.
9. П е т р о в А. В. Расчет аэродинамических характеристик крыльев со струйной механизацией // Труды ЦАГИ, 1984, вып. 2235.
10. P e t r o v A. V. Separation flow over a high-lift wings and active flow control. High-lift
and separation control / CEAS European Forum, RAES Proc.: London, 1995.
11. M a l m u t h N. D., M u r p h y W. D., S h a n k a r V., C o l e J. D., C u m b e r b a t c h E. Studies of upper surface blown airfoils in incompressible and transonic flows //
AIAA-80-1270, 1980.
12. W o l k o v A. V. a n d L y a p u n o v S. V. Numerical prediction of transonic viscous
separated flow past an airfoil // Theor. and Comp. fluid dynam. 1994. V. 6, № 1.
13. J o h n s o n D. A., B a c h a l o W. D. Transonic flow past a symmetrical airfoil — inviscid and turbulent flow propeties // AIAA J. 1980. V. 3, № 18.
14. Б о к с е р В. Д., С е р е б р и й с к и й Я. М. Приближенный метод определения
волнового сопротивления профиля при наличии местной сверхзвуковой зоны // Ученые записки ЦАГИ. 1978. Т. IX, № 5.
_________________
Рукопись поступила 14/VI 2007 г.
16
1/--страниц
Пожаловаться на содержимое документа