CC BY
56
ГАЗОВАЯ ДИНАМИКА ПРИ ДВИЖЕНИИ ШАРА В НЕПОДВИЖНОЙ ГАЗОВОЙ СРЕДЕ
Леонид И. Гречихин
Минский государственный высший авиационный колледж, Минск, Республика Беларусь
СЮ!: 10.5937/vojtehg62-5368
ОБЛАСТЬ: механика, аэродинамика, газодинамика ВИД СТАТЬИ: оригинальная научная статья
Краткое содержание:
Разработана молекулярно-кинетическая модель газодинамического обтекания шара. Установлены три области разных механизмов взаимодействия движущегося шара с неподвижной газовой средой. Показано, что срывное течение является решающим при малых скоростях движения шара, а при скоростях движения близких к скорости звука решающую роль выполняет ударное взаимодействие шара на молекулы окружающей среды в передней полусфере. Срывное течение за шаром создает конус с неподвижной газовой средой. Установлен закон, определяющий размер срывного конуса.
Ключевые слова: срывное течение, газодинамика, срывной конус, закон срывного течения.
Введение
Беспилотные летательные аппараты (БЛА) микро и мини обладают малыми размерами и совершают полет при скоростях заведомо меньших скорости звука. Поэтому для таких летательных аппаратов срывное течение играет заметную роль. В авиации срывное течение изучалось преимущественно экспериментально. Срывное течение порождает возникновение вихрей. Подробные исследования возникновения вихрей в срывном течении выполнены Теодором фон Карманом (фон Карман, 2001, стр.208). Анализируя экспериментальные факты, им были обоснованы условия воз-никновения вихрей. Поэтому шлейф вихрей за движущимся телом был назван «Вихревая дорожка Кармана». Подробное описание эксперимен-тальных фактов возникновения срывного течения дано в работе (Беляев, Червяков, 2010, стр.38). Полный гидродинамический анализ обтекания шара газовым потоком в континуальном режиме выполнен в работе (Симаков, 2013, стр.16). Общее теоретическое обоснование срывного течения без конкретного приложения описано в работе (Гречихини др. 2012, стр.285). Получат-
e-mail: [email protected]
ся, что в ряде случаях возникает срывное течение, но причины его образования остаются не достаточно ясны. В этой связи возникает цель: разработать физическую модель срывного течения путем применения молекулярно-кинетической теории и получить лобовое сопротивления шара при его движении в неподвижной среде и на этой основе обосновать причину возникновения срывного течения. Поставленная цель может быть достигнута путем решения следующих задач:
• рассмотреть динамику обтекания шара с применением молеку-лярно-кинетической теории;
• создать модель газодинамического обтекания шара с учетом формирования срывного течения;
• разработать программное обеспечение для проведения компьютерного моделирования;
• провести конкретный расчет газодинамического обтекания шара, совершающего движение в неподвижной атмосфере путем компьютерного моделирования.
Рассмотрим последовательно поставленные задачи.
Молекулярно-кинетическое обоснование обтекания газовым потоком шара
Пусть шар (рис.1) движется в воздушной среде со скоростью V вдоль оси Х-в. Выделенный элемент поверхности на шаре
AS = 2 пШ,
(.
где Д/ =Ax/cos a, R - радиус шара и а = arctg
R V1 - х2/ R2 )
(1)
х
Рис.1 - Общая схема _ газодинамического обтекания шара S/ika 1 - Opsta sema gasodinamickog opstrujavanja kugle Figure 1 - General scheme of gas dynamic flow around the sphere
1
AFyd n =~r~ = 2npv2xsinaAx (4)
Скорость движения воздуха по нормали и вдоль поверхности шара соответственно равны
vn = v sin а и vT = v cosa (2)
Элемент поверхности в виде шарового слоя AS, нормального по отношению к направлению движения возмущает массу воздуха
Am « 2п pxAxvAt
Здесь At = Ax/vcos2 а
Масса возмущенного воз-духа Am приобретает импульс нормально поверхности
Ap = Amvn = Amv sin а (3)
Изменение импульса происходит за время At. Тогда нормально действующая сила на элемент поверхности AS составит
Ap
д7
а вдоль оси Х-в
AFyd х = 2npv2xsin2 aAx (5)
Так как вдоль поверхности выделенного шарового слоя воздух дви-жется со скоростью отличной от нуля, то возникают дополнительно еще две силы: центробежная сила и сила, обусловленная действием закона Бернулли. Обе эти силы нормально поверхности шарового слоя соответственно равны
Ar2
AF„ n = 2npv2 x-;
Ц,п R (6)
AFB n = npxv2 cos aAx
Результирующая сила, которая действует на движущийся шар вдоль оси X-в
R
Fx =J (dFgx - dF^ - dFB 4). (7)
Газодинамический анализ срывного течения
За шаром возникает срывное течение. В соответствии с законом сохранения массы, какое количество массы возмущает передняя полусфера, а это
М = р0пЯ2уМ , (8)
то такое же количество массы газа должно возвратиться в заднюю полусферу.
Возмущенный поток массы в передней полусфере составляет
м
Ат =-2— = р0 V ,
жЯ2А( 0
(9)
а в задней полусфере заполнение зоны разрежения происходит со скоростью звука, т.е.,
Ат' = ,
(10)
гДе =
ук БТ
т„
- скорость звука и у = СР / Су - отношение
теплоемкостей при постоянном давлении и постоянном объеме; кБ -постоянная Больцмана, та - средняя масса молекул воздуха.
За шаром каждая молекула воздуха относительно движущегося тела вдоль оси У-ов совершает движение со скоростью звука (узв)., а вдоль оси Х-ов скорость движения молекул воздуха равна скорости движения шара с обратным знаком (-V). Результирующая скорость заполнения зоны разрежения
VI =^1 + V2
(11)
Нормальная и тангенциальная составляющие на основании рис. 1 соответственно равны
vn = V; С0Б(в + 8) и VT = V; 8т(в - 8),
(12)
где угол в = аг^
х
Я-х/1 - х2/ Я2 у
Л 8 Г VI - V) ^
, а угол 8 = агс tg
V;)
На основании закона сохранения массы количество расталкиваемой массы передней полусферой за время прохождения фронта заполнения расстояния Я - Ь , равного (рис. 2) Аt = (Я - Ь)/Vзв составляет
1
Ат = pvпЯ
2 Я - Ь
(13)
V
зв
X
Рис. 2 - Характер срывного течения за движущимся шаром Slika 2 - Karakter otcepljene struje na aerodinamickoj izlaznoj ivici kugle koja se krece Figure 2 - Character of the shear flow on the aerodynamic outward contour of the moving sphere
Эта же масса содержится в объеме
АУ = пЯ2к - 6 жИ((Я2 + 3Ь2 + к2)
На основании (13) и (14) следует, что к = Rtgа. Из рис. 2 получаем следующее равенство
л к R
Аt = — = —. у У„„
(14)
(15)
Отсюда
h = ^
v„„
(16)
Приравнивая (15) и (16) получаем для угла, при котором масса воздуха заполняет тыльную область за шаром, приобретает значение
5 = arctg
Г v ^
(17)
v у
До критического угла ар = arc tg
fh 1 ^
течение возду-
Vl - h2/ R7
ха частично компенсирует лобовое взаимодействие шара со средой. После критического угла формируется срывное течение. По мере движения потока воздуха под углом ар возникающая центробежная сила постепенно компенсируется газодинамической силой, т.е.,
h v12 2 1 2 2
р ASv1 cos(a + 5) —- sin2 (а - 5) = — p v1 sin2 (а - 5) AS (18) v R 2
Отсюда
Rv v
cos(a + 5) =-= 0,5-^-. (19)
2hv1 v1
С учетом (16) получаем v1 = v3B a + 5 = 60o, т.е., независимо от скорости движения, радиуса шара и скорости звука конус срывного течения обладает углом раскрыва 2а = 2(600 - 5). Это является законом для газодинамического обтекания шара. Отсюда следует интересный вывод: чем больше скорость движения шара, тем дальше от тела удаляется вершина конуса.
В тыльной стороне до критического угла масса воздуха, которая взаимодействует нормально с поверхностью AS = 2ny Al =
= 2nyJR2 -x2 Ax/cos—, определяется по (10). Эта масса ударно воздействует на элемент поверхности AS нормальной силой
AF;ó.,„ = 2pv2 cos2 — + 5)AS (20)
Центробежная сила нормально поверхности AS за время взаимодействия потока воздуха на ширине шарового слоя Al составляет
Ax
AFnn = pv2cos— + 5)sin—-5)——— AS (21)
R cos(—)
Сила газодинамического воздействия за счет закона Бернулли
AF,n = 1 pv2sin2(—-5)AS (22)
C3D
Результирующая сила, возникающая при заполнении зоны разрежения до расстояния h по оси Х-ов, равна
h
F^ = J sin в(- (n + dF^n + dF'En ) (23)
0
Скорость в точке начала отрыва потока
vcp. = v!sin«cp. (24)
В конусе срывного течения воздух покоится, т.е., пребывает без движения. Поэтому на тыльную поверхность шара действуют две силы - ударное воздействие хаотического движения молекул воздуха и сила сдерживания движущегося шара вследствие пониженного давления в конусе срывного течения. Обе силы направлены нормально поверхности шара.
Газодинамическое давление за шаром в конусе срывного течения определяется следующим образом
Pomp. = Poo eXP
mav3 sin a 2k БТ
.exp
mavx sin (a + S)
2k Б T
(25)
Сила газодинамического давления неподвижного газа в срывном конусе, которая сдерживает движение шара соответственно равна
К.д = РотрП^ - h')2 (26)
Общее лобовое сопротивление шара является суммой
P = F + F' + F' + F (27)
1 x 1 X ^ 1 рез. 1 уд.,ч ^ 1 г. д. ^ ' >
Получается, что обтекание газовым потоком тел шарообразной формы достаточно сложный процесс. Поэтому возникает необходимость в создании программного обеспечения для расчета динамики обтекания шара с учетом всех процессов взаимодействия движущегося шара с окружающим неподвижным газом.
Программное обеспечение длякомпьютерного моделирования.
Программное обеспечение разработано в соответствии с алгоритмом, приведенным на рис. 3.
CID
По разработанной программе выполнен конкретный расчет для шара диаметром 1 м в зависимости от скорости его движения. Результаты приведены в таблице 1.
Таблица 1 - Зависимость лобового сопротивления шара (H) диаметром 1 м от
скорости его движения Tabela 1 - Zavisnost ceonog otpora kugle precnika 1 m od brzine njenog kretanja Table 1 - Dependence of the frontal resistance of the sphere with radius 1m on its velocity
Параметры Скорость движения , м/с
0,5 1 5 10 50 100 200 250
Fлоб. 0,20 0,79 19,8 79,2 1980 7920 31700 49500
Fтыл. -0,69 -2,76 -68,8 -274 -6100 -17000 -17100 -11000
Fсрывн. 0,64 2,55 63,2 249 5250 15000 22300 19600
Fрез. 0,15 0,59 14,2 54,0 1130 5860 36900 58200
P. v2 S, H 0,25 1,01 25,3 101 2530 10100 40500 63300
Комьютерное моделирование
Разработанное программное обеспечение позволяет осуществить анализ разных режимов движущегося шара в зависимости от скорости движения в неподвижной атмосфере, радиуса шара, состояния атмосферы и высоты полета. В качестве примера в таблице приведены результаты расчета разных действующих сил в сравнении с газодинамическим напором на мидель движущегося шара.
Выводы и предложения
Таким образом, молекулярно-кинетический анализ газодинамического обтекания шара свидетельствует о следующем. При движении шара в неподвижной атмосфере возникают три области разного механизма взаимодействия с молекулами окружающей среды:
1 - воздействие шара своей передней полусферой путем ударного действия на молекулы окружающей атмосферы со скоростью движения шара, а возникающие тангенциальные скорости движения воздуха вокруг шара формируют центробежные силы и газодинамические силы, обусловленные действием закона Бернулли.
2 - заполнение зоны разрежения до образования срывного течения происходит со скоростью звука, где возникают силы ударного действия потока молекул воздуха нормальной составляющей, а тангенциальная составляющая потока создает центробежную силу, а также газодинамическую силу в соответствии с законом Бернулли.
3 - образующийся срывной конус неподвижного воздуха создает силу газодинамического давления, обусловленного разностью давлений окружающей среды и давлением в срывном конусе.
<Е>
В таблице приведены результирующие силы, возникающие в со трех областях разного механизма взаимодействия движущегося шара с окружающей атмосферой. При сравнительно малых скоростях движения шара передняя полусфера ударным своим действием не компенсирует обратного действия потока воздуха в задней полусфере до образования срывного течения. Эту роль выполняет область срывного течения, т.е., парадокс корпускулярной теории Ньютона имеет место, но только частично. §
По мере приближения скорости движения шара к скорости звука, наоборот, ударное действие передней полусферы в основном определяет лобовое сопротивление шара. Полеты беспилотных летательных | аппаратов реализуются со скоростями не более 100 м/с. Следовательно, газовая динамика полета этих аппаратов определяется в основном срывным течением, т.е., задней полусферой. Тогда как пилотируемые большие летательные аппараты совершают полет вблизи скорости зву- <о ка, где газовая динамика полета определяется в основном механизмом взаимодействия передней полусферы с окружающим воздухом. Вот почему разработанную аэродинамику для пилотируемых летательных аппаратов нельзя применять для беспилотных летательных аппаратов. |
ш
Особо следует отметить, что для движущегося шара в неподви- J жной атмосфере имеет место закон: Сумма угла срыва потока и угла между результирующей скоростью и скоростью звука при выполнении неравенства v3e >> v равна примерно 600. Этот закон определяет размер конуса срывного течения.
Обычно полагают, что в передней полусфере давление возрастает примерно в два раза. Поэтому максимальное значение газодинамического давления составляет рхv2SMUd Эти значения в качестве сравнения приведены в таблице. Лобовое сопротивление шара меньше максимального газодинамического давления, а по мере приближения скорости движения шара к скорости звука эти давления становятся соизмеримыми.
Литература
Беляев, Е.Н., & Червяков, В.В. 2010. Срывные течения жидкостной и газовой фазы. Двигатель, 69(3), стр. 38.
Гречихин, Л.И., Лапцевич, А.А., & Куць, Н.Г. 2012. Аэродинамика летательных аппаратов.Мн: ИООО «Право и экономика», стр. 285.
Симаков, Н.Н. 2013. Расчет обтекания и сопротивления шара в ламинарном и сильно турбулентном потоке .Журнал технической физики, 83(4), стр. 16-20.
фон Карман, Т. 2001. Аэродинамика: Избранные темы в их историческом развитии.Ижевск., стр. 208.
GASNA DINAMIKA PRI KRETANJU KUGLE U NEPOKRETNOJ GASNOJ SREDINI
Leonid Ivanovic Grecihin
Drzavni visi koledz za vojno vazduhoplovstvo u Minsku, Republika Belorusija
OBLAST: masinstvo, aerodinamika, gasna dinamika VRSTA CLANKA: originalni naucni clanak
Sazetak:
U radu je razraden matematicki model gasodinamickog opstruja-vanja kugle. Postoje tri podrucja razlicitih mehanizama uzajamnog dej-stva pokretne kugle i nepokretne gasne sredine. Pokazalo se da je ot-cepljena struja presudna pri malim brzinama kretanja kugle, a pri brzi-nama kretanja kugle koje su bliske brzini zvuka presudnu ulogu ispolja-va udarno uzajamno dejstvo kugle sa sredinom u prednjoj polusferi. Otcepljena struja iza kugle cini konus sa nepokretnom gasnom sredinom. Utvrdena je zakonitost koja definise razmere tog konusa.
Kljucne reci: otcepljena struja, gasna dinamika, otcepljeni konus, za-kon otcepljene struje.
GAS DYNAMICS DURING THE SPHERE MOVING IN THE STATIONARY GASEOUS ENVIRONMENT
Leonid Ivanovich Gretchihin
Minsk State Higher Aviation College, Department of Natural Science Disciplines, Minsk, Republic of Belarus
FIELD: Mechanics, Aerodynamics, Gas Dynamics ARTICLE TYPE: Original Scientific Paper
Abstract:
This paper developed a mathematical model of the gas dynamic fluid flow around the sphere. There are three areas of different mechanisms of the interaction of moving spheres and stationary gaseous environment. It has been proved that shear flow plays a decisiverole at sphere low velocities, while at sphere velocities close to the speed of sound the critical role is exerted by the impact interaction of the sphere with the center in the front hemisphere. The shear flow behind the sphere makes a cone with stationary gaseous environment. The principle which defines the size of the cone has been established.
Key words: shear flow, gas dynamics, stalling cone, law of shear flow. Дата получения работы/Paper received on: 21. 01. 2014.
Дата получения исправленной версии работы/Manuscript corrections submitted on: 02. 04. 2014.
Дата окончательного согласования работы /Paper accepted for publishing on: 04. 04. 2014.
