CC BY
13
УДК 533.16
Моделирование расхода газа через ламинарный пограничный слой на поверхности полусферы в сверхзвуковом воздушном потоке
© В.В. Горский1'2, В. А. Сысенко2
1 МГТУ им. Н.Э. Баумана, Москва, 105005, Россия
2
ОАО «ВПК «НПО машиностроения», г. Реутов Московской обл., 143966, Россия
Приведены результаты оценки точности для инженерной методики расчета массового расхода газа через ламинарный пограничный слой на полусфере из работы [1]. Предложена аналогичная инженерная методика повышенной точности.
Ключевые слова: массовый расход, пограничный слой, воздушный поток, газовый поток, термохимическое равновесие, тепломассоперенос.
Введение. Расчет массового расхода газа через ламинарный пограничный слой широко используют в теории пограничного слоя. В процессе численного решения уравнений пограничного слоя этот расход находят изопределяющего его из интегрального выражения [1]. Однако при использовании различных интегральных решений уравнений пограничного слоя информация о расходе газа может быть получена только путем введения определенных допущений о профилях газодинамических функций в пограничном слое.
Для полусферы, в частности, в работе [1] приведена методика расчета расхода газа через ламинарный пограничный слой, полученная способом эффективной длины и скорректированная по результатам строгих численных решений уравнений пограничного слоя. Однако информация о точности этой методики в литературе отсутствует.
Решению этой задачи и рассмотрению возможности построения более корректной методики этого типа на базе систематических решений уравнений ламинарного пограничного слоя и посвящена данная работа.
Расчетно-теоретические исследования, результаты которых приведены далее, получены в широком диапазоне изменения определяющих факторов: числа Маха Мот в набегающем воздушном потоке,
давления торможения Ро этого потока и энтальпийного фактора , под которым понимают отношение энтальпии воздуха при температуре «стенки» к энтальпии торможения газового потока. Под термином «стенка» здесь и далее понимают поверхность тела, обтекаемую набегающим на него газовым потоком.
Исследования выполнены для многокомпонентной газовой смеси в пограничном слое на абсолютно каталитичной «стенке», находящейся в состоянии термохимического равновесия с расчетом диффузионного тепломассопереноса в рамках уравнений Стефана - Максвелла.
В процессе проведения исследований принято:
• химический состав газовой смеси ограничен набором химических веществ О, 02, К, N2, N0, Аг, образованных из химических элементов О, К, Аг;
• в нормальных условиях воздух характеризуется следующим мольным химическим составом [2]:
к о2 =0,2095; к к2 =0,7808; к Аг =0,0097; • область изменения определяющих параметров задана в виде
Ыме[4,25], Р0 е[0,001,ртах], Я е
п п
И,тт' И, тах
• для расчета переносных свойств многокомпонентной газовой смеси использован метод Гиршфельдера [3];
• использованы потенциальная функция межмолекулярного взаимодействия Леннарда - Джонса [3] и методика расчета параметров этой функции из [4], базирующаяся на современных расчетно-теоретических данных по вязкости воздуха, находящегося в состоянии термохимического равновесия [5, 6];
• использованы уточненные данные по параметрам идеального газа, набегающего на «стенку», полученные методом сплайновой аппроксимации [7];
• давление торможения измеряется в МПа;
• Ртах = 10 при Мот >10 и ртах = 1 в противном случае;
• Я шп — значение энтальпийного фактора, соответствующее
температуре «стенки», равной примерно 300 К;
• Значение Я тах минимально в диапазоне между 0,7 и значением энтальпийного фактора, соответствующим температуре кипения атомарного углерода.
Результаты расчетно-теоретических исследований. Массовый расход газа через пограничный слой ехр относится к числу важных характеристик пограничного слоя и определяется интегральным выражением вида
Уе
ехр = 2 кг | рпёу, (1)
0
где у — координата, отсчитываемая от «стенки» в направлении внешней нормали к ней; г — удаление образующей сферы от оси,
проходящей через ее «критическую точку» в направлении вектора скорости набегающего газового потока; р, u — плотность газовой смеси и тангенциальная составляющая вектора скорости в пограничном слое. Индекс е относится к внешней границе пограничного слоя. При использовании переменных ц Лиза - Дородницына
£( s) = Jpe^euer2ds'; ц( s,y) = UU~ Jpdy' и безразмерной функции тока
о о
f, определяемой выражением f = u/ue, формула (1) принимает вид
exp = 2nfe^. (2)
Здесь j — криволинейная координата, отсчитываемая вдоль образующей сферы от ее «критической точки», т.е. от точки торможения набегающего на сферу газового потока; fe — значение безразмерной функции тока на условной внешней границе пограничного слоя, на которой координата це удовлетворяет условию 1 - f < 0,005 при
ц>це. Индекс ц означает частную производную по этой координате.
Численное интегрирование уравнений пограничного слоя в данной работе проведено на неравномерной сетке по координате ц, в которой шаг между узлами этой координаты увеличивается по определенному закону по мере удаления от «стенки».
В работе [1] приведена следующая приближенная формула для расчета массового расхода газа через тонкий ламинарный пограничный слой на непроницаемой «стенке» для полусферы:
exp (J) = PeMeUeXeff (J). (3)
Здесь — коэффициент динамической вязкости газа; xeff (s) —
длина цилиндра, обтекаемого газовым потоком с параметрами, соответствующими точке полусферы с координатой s, для которого выполняются условия [8]:
• равенства удельного теплового потока в его концевом сечении аналогичному потоку на полусфере в точке с координатой s;
• равенства интегрального теплового потока к поверхности цилиндра аналогичному потоку на поверхности полусферы, ограниченной сверху сечением с координатой s.
На рис. 1 приведены результаты сопоставления расхода газа через пограничный слой по поверхности полусферы, рассчитанные по формуле (3), с аналогичными результатами численных решений уравнений пограничного слоя. Здесь Sexp = explit/expnum; aeff = л/2 -s —
эффективный угол атаки, рад; aeff Ит = max (aeff;0,01) — эффектив-
ный угол атаки, ограниченный снизу значением, равным 0,01; ехрпит
— толщина потери импульса, найденная в результате численного решения уравнений пограничного слоя.
"ехр 1,20
°'7°0 2 4 6 8 10 12 14 16 18 20 22 24 26 (^0^
Рис. 1. Зависимость отношения массового расхода газа через ламинарный пограничный слой, полученного по формуле (3) из работы [1] и в результате численного решения уравнений пограничного слоя, от эффективного угла атаки: Мм = 25 (ооо); 15 (□□□); 8 (000); 4 (АДА)
Как следует из рассмотрения представленных данных, погрешность инженерного расчета по формуле (3) не превышает примерно 25 %, что является приемлемым для многих технических приложений. В то же время актуально проведение исследований, направленных на создание аналогичной инженерной методики расчета расхода газа через пограничный слой, характеризуемой существенно меньшей погрешностью, чем формула (3).
Как показали результаты проведенных исследований, наиболее просто эта задача решается путем:
• введения поправочной функции к формуле (3), т.е. расчета расхода газа через пограничный слой по формуле
exp
(aeff ) = eXPlit (aeff ) Fexp ( Rh , aeff,
); (4) введения функции ^(aeff Иш), определенной на области |0
вида
С(а
ей, Кт
\ 1 • = — min
I 3
6,ctg (а
й 11т
1 1 •
1--min
12
6,Ctg (aeff, lim )
расчета поправочной функции по формуле
F =
exp
0,8 + 10,0056 + 0,2789Rh - 0,3961Rh +1,0883R^3
)c(a
eff, lim
-1
На рис. 2 представлено сопоставление расходов газа через ламинарный пограничный слой на поверхности полусферы, рассчитанных по формулам (3) и (4), с аналогичными данными, полученными в рамках численного решения уравнений пограничного слоя. Здесь
е«ф, к = ^к ^ш; к=^ ^ num.
1,20 а
'ехр,пшп
Рис. 2. Сопоставление массового расхода газа через ламинарный пограничный слой на полусфере, рассчитанного по формулам (3) и (4), с аналогичными данными, полученными в рамках численного решения уравнений пограничного слоя: 1 — k = lit; 2 — k = eng
Из представленных на этом рисунке данных следует, что переход в инженерных расчетах от использования формулы (3) к применению формулы (4) позволяет качественно снизить погрешность. Так, в частности, максимальная погрешность вычислений снижается при этом с 25 до 9 %.
Выводы. Установлено, что погрешность широко используемой на практике инженерной формулы, предназначенной для расчета толщины потери импульса в ламинарном пограничном слое на полусфере, сопряжена с внесением в расчет погрешностей, достигающих 25 %.
Предложена модификация этой формулы, применение которой позволяет снизить максимальную погрешность вычислений до 9 %.
ЛИТЕРАТУРА
[1] Землянский Б.А., Лунев В.В., Власов В.И. и др. Конвективный теплообмен летательных аппаратов. Москва, Физматлит, 2014, 377 с.
[2] Предводителев А.С., Ступоченко Е.В., Плешанов А.С. и др. Таблицы термодинамических функций воздуха (для температур от 200 до 6000 К и давлений от 0,00001 до 100 атмосфер). Москва, ВЦ АН СССР, 1962, 268 с.
[3] Гиршфельдер Дж., Кертис Ч., Берд Р. Молекулярная теория газов и жидкостей. Москва, Изд-во иностранной литературы, 1961, 929 с.
[4] Горский В.В., Федоров С.Н. Об одном подходе к расчету вязкости диссоциированных газовых смесей, образованных из кислорода, азота и углерода. Инженерно-физический журнал, 2007, т. 80, № 5, с. 97-101.
[5] Соколова И. А. Коэффициенты переноса и интегралы столкновений воздуха и его компонент. В кн.: Физическая кинетика, Аэрофизические исследования. Новосибирск, Институт теоретической и прикладной механики СО АН СССР. Сб. трудов № 4, 1974, с. 39-104.
[6] Capitelli M., Colonna G., Gorse C., D'Angola A. Transport properties of high temperature air in local thermodynamic equilibrium. The European Physical Journal, 2000, no. 11, рp. 279-289.
[7] Горский В.В. Метод сплайновой аппроксимации. Журнал вычислительной математики и математической физики РАН, 2007, т. 47, № 6, с. 939-943.
[8] Авдуевский В.С., Галицейский Б.М., Глебов Г.А. Основы теплопередачи в авиационной и ракетно-космической технике. В.К. Кошкин, ред. Москва, Машиностроение, 1975, 624 с.
Статья поступила в редакцию 05.11.2014
Ссылку на эту статью просим оформлять следующим образом: Горский В.В., Сысенко В.А. Моделирование расхода газа через ламинарный пограничный слой на поверхности полусферы в сверхзвуковом воздушном потоке. Математическое моделирование и численные методы, 2014, № 4, с. 88-94.
Горский Валерий Владимирович родился в 1939 г., окончил Московский авиационный институт им. Серго Ордоникидзе в 1963 г. Главный научный сотрудник ОАО «ВПК «НПО машиностроения». Д-р техн. наук, профессор кафедры «Вычислительная математика и математическая физика» МГТУ им. Н.Э. Баумана. Автор более 100 научных работ.
Сысенко Валентина Алексеевна родилась в 1978 г., окончила МГТУ им. Н.Э. Баумана в 2001 г. Канд. техн. наук, старший научный сотрудник отдела аэродинамики ОАО «ВПК «НПО машиностроения». Автор 12 научных работ в области прикладной математики. e-mail: [email protected]
B.B. ropcKuu, B.A. CbicenKO
Simulation of gas flow through the laminar boundary layer on the hemisphere surface in a supersonic air flow
© V.V. Gorskiy1,2, V.A. Sysenko2
1 Bauman Moscow State Technical University, Moscow, 105005, Russia
2 JSC "MIC "NPO Mashinostroyenia", Reutov-town, Moscow Region, 143966, Russia
The article presents estimated accuracy of the engineering design procedure of the mass flow rate of gas through the laminar boundary layer on a hemisphere of [1]. A similar engineering method of extra accuracy is proposed.
Keywords: mass flow, boundary layer, air flow, gas flow, thermochemical equilibrium, heat and mass transfer.
REFERENCES
[1] Zemlyanskiy B.A., Lunev V.V., Vlasov V.I. Konvektivnyi teploobmen le-tatelnykh apparatov [Thermal convection of Aircraft]. Moscow, Fizmatlit Publ., 2014, 377 p.
[2] Predvoditelev A.S., Stupochenko E.V., Pleshanov A.S. Tablitsy termodinamich-eskikh funktsiy vozdukha (dlya temperature ot 200 do 6000 K i davleniya ot 0,00001 do 100 atm.) [Tables of air thermodynamic functions (for temperatures between 200 and 6000 K and pressures from 0.00001 to 100 atm)]. Moscow, Vychislitelnyi Tsentr AN SSSR Publ., 1962, 268 p.
[3] Girshfelder G., Kertis Ch., Berd R. Molekulyarnaya teoriya gazov i zhidkostey [Molecular theory of gases and liquids]. Moscow, Inostrannaya Literatura Publ., 1961, 929 p.
[4] Gorskiy V.V., Fedorov S.N. Inzhenerno-Fizicheskii Zhurnal — Journal of Engineering Thermophysics, 2007, vol. 80, no. 5, pp. 97-101.
[5] Sokolova I.A. Koeffitsienty perenosa i integraly stolknoveniy vozdukha i ego component [Transfer coefficients and collision integrals of air and its components]. In: Fizicheskaya kinetika. Aerofizicheckie issledovaniya [Physical kinetics. Aerophysical research]. Novosibirsk, Proceedings of the Institute of Theoretical and Applied Mechanics SO AS USSR, no. 4, 1974, pp. 39-104.
[6] Capitelli M., Colonna G., Gorse C., D'Angola A. Transport properties of high temperature air in local thermodynamic equilibrium. The European Physical Journal, 2000, no. 11, pp. 279-289.
[7] Gorskiy V.V. Zhurnal vychislitelnoy matematiki i vychislitelnoy fiziki RAN — Journal of Computational Mathematics and Computational Physics RAS, 2007, vol. 47, no. 6, pp. 939-943.
[8] Avduevskiy V.S., Galitseyskiy B.M., Glebov G.A. Osnovy teploperedachi v avi-atsionnoy i raketno-kosmicheskoy tekhnike [Fundamentals of heat transfer in the aviation and aerospace technology]. Koshkin V.K., ed. Moscow, Mashi-nostroenie Publ., 1975, 624 p.
Gorskiy V.V. (b.1939) graduated from Moscow Aviation Institute named after Sergo Ordzhonikidze in 1963. Senior staff scientist at the JSC "MIC "NPO Mashinostroyenia". Dr. Sci. (Eng.), professor of the Computational Mathematics and Mathematical Physics Department at Bauman Moscow State Technical University. Author of more than 133 publications. Scientific interests: ablation heat protection, high temperature heat-mass transfer, numerical methods of solving equations of mathematical physics.
Sysenko V.A. (b.1978) graduated from Bauman Moscow State Technical University in 2001. Ph.D., senior scientist of the Aerodynamics Department at the JSC "MIC "NPO Mashinostroyenia". Author of 12 publications in the field of applied mathematics. Scientific interests: high temperature heat-mass transfer, numerical methods of solving equations of mathematical physics. e-mail: [email protected]
