CC BY
169
2009
НАУЧНЫЙ ВЕСТНИК МГТУ ГА серия Аэромеханика, прочность, поддержание летной годности ВС
№141
УДК: 533.601 (075.8)
АЭРОДИНАМИЧЕСКИЕ ХАРАКТЕРИСТИКИ ГЛА ДКИХ ТЕЛ ВРАЩЕНИЯ В ДОЗВУКОВОМ ПОТОКЕ
ЛЕ КУОК ДИНЬ, Н.В. СЕМЕНЧИКОВ, ЧАН ВАН ВЬЕТ, О.В. ЯКОВЛЕВСКИЙ Статья представлена доктором технических наук, профессором Ципенко В.Г.
Приводятся результаты численных исследований обтекания гладких тел вращения турбулентным потоком вязкого воздуха малых скоростей. Получены распределения давления в меридиональных сечениях тела, коэффициенты сил и моментов, действующих на тело. Показано влияние удлинения тела на его аэродинамические характеристики.
Ключевые слова: аэродинамические характеристики, турбулентность, крыло, обтекание, число Рейнольдса.
Исследование обтекания тел вращения различной формы проводится в течение многих лет. К настоящему времени накоплен обширный материал по их аэродинамическим характеристикам и параметрам обтекания [1-3]. Тем не менее, особенности изменения аэродинамических характеристик гладких тел различного удлинения в зависимости от изменения различных определяющих их параметров рассмотрены еще недостаточно.
В данной работе представлены результаты численных исследований обтекания турбулентным потоком вязкого воздуха гладких тел вращения при углах атаки, изменявшихся в пределах от 0 до 90°. Исследования проведены при скоростях невозмущенного потока У¥ = 20...50 м/с. При этом число Рейнольдса, подсчитанное по длине тела, составляло Яе = (1,6...4)х106.
Рассматривались гладкие тела вращения различного удлинения, координаты контура которых в плоскости "ХОУ", связанной с телом системы координат с началом в вершине тела, определялись соотношением [4]:
у = 0,972£[ х (1 - х )(1,5 - х )]1/2,
где х = х/Ь, Ь - длина тела; Б - диаметр миделевого сечения тела.
Для решения задачи использовались осредненные по Рейнольдсу уравнения Навье-Стокса, записанные для пространственного случая течения [5] и замкнутые моделью турбулентности "кв". На теле ставилось условие прилипания, на входе в счетную область задавалась скорость невозмущенного потока. Задача решалась численно. Расчетная сетка имела блочную структуру. Структура счетной области выбиралась из условия минимизации количества ячеек, форма которых ухудшала сходимость и аппроксимацию. Для достоверного определения газодинамических параметров сетка была адаптирована под особенности течения и имела сгущение у поверхности тела, где имеется развитый пограничный слой.
В результате расчетов были получены параметры потока в расчетных точках, найдены в связанной системе координат коэффициенты сил и моментов, действующих на рассматриваемые тела. Коэффициенты сил и момента тангажа были отнесены к '2/3. Коэффициенты момента тангажа были подсчитаны относительно центра объема тела и отнесены к ', где W - объем тела.
На рис. 1 - 3 представлены некоторые результаты расчетов, позволяющие составить представление о влиянии на величины коэффициентов сил и момента тангажа тел углов атаки, удлинения тела и чисел Рейнольдса.
На рис. 1 приведены типичные зависимости Сх = Да), Су = Да) и ш2 = Да), полученные для тел различного удлинения при скорости невозмущенного потока Уда = 30 м/с, что соответствовало числу Яе = 2,4х106. Видно, что величины коэффициента продольной силы нелинейно изменяются по углу атаки.
Рис. 1. Зависимости Сх = :Т(а), Су = :Т(а) и ш2 = Г(а) для тел различных удлинений X при скорости невозмущенного потока Уда = 30 м/с
в
Характер указанного изменения сильно зависит от удлинения тела. При этом в наибольшей степени изменение положения тела по углу атаки влияет на коэффициенты Сх при удлинении тела X = 4,5, которое примерно соответствует удлинениям малоразмерных дирижаблей. Как видим, рост угла атаки вызывает падение величин указанного коэффициента, причем они при а ~ 32° становятся отрицательными. При угле атаки а = 80° на кривой Сх = f(a) этого тела наблюдается характерный минимум. При а >80° величины коэффициента Сх снова растут. Отмеченный минимум величин Сх имеет место и в случае больших, чем 4,5, удлинений тела. Но он наблюдается при угле а = 70° и при этом также сильно уменьшается абсолютная величина коэффициента Cxmin (хотя он по-прежнему остается отрицательным).
Зависимости коэффициента нормальной силы от угла атаки также являются нелинейными. При углах атаки а < 20° величины коэффициента нормальной силы рассматриваемых тел слабо зависят от удлинения тела. Но при больших, чем а = 20°, углах атаки наблюдается большое различие в величинах этого коэффициента в зависимости от удлинения тела. При а = Const в диапазоне углов атаки а = 20°...50° наибольшими значениями коэффициентов нормальной силы
обладает тело с удлинением X = 4,5, а при а > 50° наибольшие значения коэффициента Су имеет тело с удлинением X = 12.
Зависимость коэффициента момента тангажа от угла атаки, показанная на рис. 1в, является типичной для неоперенных тел. Как видим, изменение удлинения тела в пределах от 4,5 до 8 при углах атаки а < 60° слабо влияет на величины указанного коэффициента. Существенное влияние удлинения начинает проявляться при X > 8, причем максимум величин коэффициентов момента тангажа всех рассмотренных тел имеет место при углах атаки а = 45°...50°.
На рис. 2 приведены зависимости Сх = Да), Су = Да) и т2 = У(а) для тела с удлинением X = 4,5 при различных скоростях набегающего на тело потока, а на рис. 3 - аналогичные зависимости для тела с удлинением X = 8.
0,15 0,1 0,05
о
-0,05 -0,1 -0,15 -0,2 -0,25 -0,3
а На
а
Су 1,2 1
0,8 0,6 0,4
0.2 -(— -----------------------------------------------------------
0 ---г- 1 I -1----------1------------
о 20 40 60 80 100
а На
б
Шг
1,2 1
0,8 0,6
0,4 0,2 0
0 20 40 60 80 100
а На
в
Рис. 2. Зависимости Сх = Д(а), Су = Да) и ш2 = Да) для тела с удлинением X = 4,5 при различных скоростях невозмущенного потока
Шх
аКа
в)
у=20
У=30
-~т4г~- У=50
Рис. 3. Зависимости Сх = Д(а), Су = Д(а) и ш2 = Д(а) для тела с удлинением X = 8 при различных скоростях невозмущенного потока
б
Видно, что изменение скорости набегающего на тело потока влияет по-разному на величины коэффициентов сил и момента тангажа в зависимости от удлинения тела. Так, из приведенных на графиках примеров видно, что изменение скорости набегающего потока может практически не отражаться на величинах коэффициентов продольной силы и момента тангажа тела с удлинением X = 4,5. И в то же время в случае тела с удлинением X = 8 эти коэффициенты в большей степени подвержены влиянию указанного фактора. Как показывают расчеты, у тела с удлинением X = 12 коэффициенты момента тангажа снова слабо зависят от величины скорости
невозмущенного потока. Отметим, что для всех рассмотренных тел ее рост при всех углах а = Const вызывает падение величин коэффициента нормальной силы, причем в большей степени для тела с удлинением X = 4,5.
ЛИТЕРАТУРА
1. Краснов Н.Ф. Аэродинамика тел вращения. - М.: Машиностроение, 1964.
2. Петров К.П. Аэродинамика элементов летательных аппаратов. - М.: Машиностроение, 1985.
3. Аэродинамика летательных аппаратов; Под ред. Г. А. Колесникова. - М.: Машиностроение, 1993.
4. Фомина Н.Н. Атлас форм корпусов дирижаблей. // Труды ЦАГИ, вып. 238, 1935.
5. Бондарев Е.Н., Дубасов В.Т., Рыжов Ю.А., Свирщевский С.Б., Семенчиков Н.В. Аэрогидромеханика. - М.: Машиностроение, 1993.
AERODYNAMIC CHARACTERISTICS OF A SMOOTH BODIES OF REVOLUTION
Le Quok Dinh, Semenchikov N.V., Chang Wang Viet, Yakovlevsky O.V.
The results of numerical researches of a viscous air turbulent flow around a smooth bodies of revolution are resulted at small and large angles of attack. The coefficients of forces and pitch moment are received. The appreciable influence of body aspect ratio on the aerodynamic characteristics of a bodies is found.
Сведения об авторах
Ле Куок Динь, 1986 г.р., студент - дипломник кафедры аэродинамики летательных аппаратов МАИ (ГТУ) им. С. Орджоникидзе, область научных интересов - аэродинамика летательных аппаратов.
Семенчиков Николай Витальевич, 1941 г.р., окончил МАИ им. С. Орджоникидзе (1964), кандидат технических наук, профессор кафедры аэродинамики летательных аппаратов МАИ (ГТУ) им. С. Орджоникидзе, автор 120 научных работ, область научных интересов - вихревые, струйные и отрывные течения, аэродинамика летательных аппаратов.
Чан Ван Вьет, 1986 г.р., студент - дипломник кафедры аэродинамики летательных аппаратов МАИ (ГТУ) им. С. Орджоникидзе, область научных интересов - аэродинамика летательных аппаратов.
Яковлевский Олег Васильевич, 1932 г.р., окончил МФТИ (1955), кандидат технических наук, профессор кафедры аэродинамики летательных аппаратов МАИ (ГТУ) им. С. Орджоникидзе, автор более 120 научных работ, область научных интересов - турбулентные струйные течения, аэроакустика, аэродинамика летательных аппаратов и промышленная аэродинамика.
