Аналитические формулы повышенной точности для расчета распределения давления на поверхности выпуклых затупленных тел вращения произвольного очертания

УДК 533.6.011.31.5:532.582.33

Аналитические формулы повышенной точности для расчета распределения давления на поверхности выпуклых затупленных тел вращения произвольного очертания

© В.П. Котенев1, В.А. Сысенко2

1 МГТУ им. Н.Э. Баумана, Москва, 105005, Россия

2 ВПК «НПО машиностроения», г. Реутов, Московская область, 143966, Россия

Разработаны аналитические формулы для быстрого и точного расчета давления на участке поверхности тел вращения произвольного очертания, обтекаемых сверхзвуковым потоком газа. Рассмотрены примеры применения метода для пространственных течений газа.

Ключевые слова: сверхзвуковой поток, трехмерные течения газа, звуковая точка.

Введение. При применении численных алгоритмов, основанных на методе установления по времени или других итерационных подходах [1-4], успех решения задачи во многом зависит от близости начальных данных к стационарному решению. В настоящее время существует определенный разрыв между возможностями численных методов и точностью приближенных решений, что значительно увеличивает время счета или делает невозможным численное интегрирование уравнений газовой динамики. Так, использование формулы Ньютона [5] для начального распределения давления на участках, где угол встречи потока с поверхностью тела мал, часто приводит к срыву решения.

Между тем простые аналитические подходы имеют и большое самостоятельное значение [6].

В данной статье, являющейся развитием работ [7-9], получены аналитические зависимости для практически точного определения давления на выпуклых затупленных телах вращения произвольного очертания. Предложен также подход к применению этих зависимостей в трехмерных течениях.

Получение аналитической зависимости для определения давления. Рассмотрим вывод простой аналитической зависимости, дающей практически точное распределение давления на затупленных телах вращения произвольного очертания при их сверхзвуковом обтекании.

1. Обезразмеривание. Будем использовать безразмерные параметры. Давление Р отнесем к давлению в точке торможения Р0, которое определяется по известной формуле Рэлея:

у+1

Р = |Х±! М2

у-1

У 2М2

У-1 Р

1 п

СО ?

где у — показатель адиабаты (для совершенного газа у = 1,4); МС — число Маха набегающего потока; РС — давление газа в набегающем потоке.

2. Определение положения звуковой точки на сфере. Для чисел Маха в набегающем потоке более 2,5 положение звуковой точки, расчет которой приводится в [10], является функцией от отношения плотностей к на прямом скачке уплотнения:

а« = 90 -

1 2

к = у---

у + 1 (у + 1) МСС'

где а» — угол в градусах между осью тела и вектором скорости в звуковой точке.

Для умеренных сверхзвуковых чисел Маха хороший результат для определения положения звуковой точки на сфере дает использование классической формулы Ньютона:

а» = агсБт

1

1 + уМСС

1 -

1+умСС

где р* =

у +1

давление в звуковой точке, отнесенное к

давлению в точке торможения.

3. Вычисление давления на поверхности сферы. Рассмотрим функцию т, характеризующую термодинамическое состояние газа и представляющую собой отношение температуры движущегося газа к температуре изоэнтропически заторможенного потока. Величина р = ту/(у-1) есть отношение статического давления к давлению торможения, т. е. искомая величина. Как известно, функция т является вогнутой вниз плавно убывающей функцией от скорости, которая при движении газа вдоль сферы является, в свою очередь, монотонно возрастающей функцией от угла а между вектором скорости набегающего потока и вектором скорости в любой точке ее поверхности сферы. Максимальное значение т достигается в точке торможения и равно единице. Поэтому

/ , Чч 1 - tg2 (ц/2) будем считать, что т = cos (ц(а)) = --2 ^ ^. Найдем функцию

т

при а = а*

1 - tg2

= 1,

1 - tg2 (ц(а** )/2)_ 2

1 + tg2 (ц(а** )/2) у +1

Из выражений (3) и (4) следует, что tg(ц(л/2) /2) = 0 и tg (ц(а** )/2) =

у-1 Г и, -. Предположим, что tg — 1 = a а + b, тогда

У + 3 12,

tg

V 2 J

V

у-1 а-л/2

у + 3 а» - л /2

и

1 -

т =

у-1

у + 3

(

-л/2 ^

а** - л /2

1 +

у-1

у + 3

(

•-л/2 ^

а** - л /2

или

т = cos

2 arctg

у-1 л/2-а

'у + 3 л/2-а*

Формулы (5) или (6) наряду с (1), (2) позволяют найти давление

а] = т (а) . Входящие в них углы содержатся в числителе и

знаменателе дроби, поэтому их можно брать как в радианах, так и в градусах.

4. Определение положения звуковой точки для произвольного выпуклого тела вращения. Для затупленного тела вращения определим положение звуковой точки а*, воспользовавшись методом работы [11], применение которого можно также найти в [12, 13].

5. Расчет давления на поверхности затупленного тела вращения. Распределение давления на поверхности затупленного тела вращения будем искать в виде

2

2

Р =

1 -

У -1 у+ 3

с

о-л/2 о* — л/2

2л X (о)^-V у—1

1 +

у -1

(

у+ 3

о-л/2 о* - л / 2

2

V у

Так как в звуковой точке должно выполняться равенство

Г

р* =

Л

у +1

у-1

, то X (о*) = 1. Будем считать функцию Х(о) линейной.

Пусть вторая точка Х0, через которую проходит прямая Х(о), лежит в небольшой окрестности о* и Х( х0) = Х0. Тогда Х(о) имеет вид

Х(о) = 1 +

1 - Х 0 _ \ (о - о*).

О* — Л"0

Положение точки Х0 определим эмпирически в зависимости от положения звуковой точки на поверхности тела из условия: если о» >о*, то X) = 1,15 о*; если о** <о*, то х0 = 0,85 о*.

В небольшой окрестности звуковой точки для расчета давления с высокой точностью можно воспользоваться модификацией формулы Ньютона:

если о > о*, то Р = б1п о + -

2

СОБ2 о*

-соб о;

2

Б1П2 о

если о<о*, то Р = (р* - Р„ )—г--ъ Рю.

Б1П2 о*

Величину Х0 будем искать из равенства давления, рассчитываемого по формуле (7), давлению, рассчитываемому по формулам (9) или (10) в точке Х0. Таким образом, в любой точке тела можно вычислить значение Х(о), а значит, и давление по формуле (7).

Анализ результатов. Для примера приведем результаты расчетов по предложенному методу для двух типов эллипсоидов с соотношением полуосей Ыа = 1/2 (рис. 1) и Ыа = 3/2 (рис. 2) и для параболоидов вида г2 = 2рх (где 2 — безразмерная величина, отнесенная к радиусу кривизны тела в критической точке) с параметрами р = 0,25 (рис. 3, а) и р = 1 (рис. 3, б-г и рис. 4, а, б) при различных значениях числа Маха. Как видно из сравнения расчетных данных с табличными [14, 15], применение формулы (7) позволяет получить результат с максимальной относительной погрешностью, не превышающей 5 % для всех вариантов. В то же время на некоторых режимах формула Ньютона имеет максимальную относительную погрешность до 100 %.

у

Р 0,9

0,7

0,5

0,3

0,1 0

Р 0,9

0,7

0,5

0,3

0,1 О

Р 0,9

0,7

0,5

0,3

0,1 О

в

Р 0,9

0,7 0,5 0,3 0,1

О 10 20 30 40 50 60 70 о, град

г

Рис. 1. Распределение давления Р на эллипсоиде с соотношением полуосей Ыа = 1/2 при числе Маха М = 2 (а), М = 4 (б), М = 6 (в) и М = 20 (г),

полученное по формуле Ньютона (---), по формуле (7) (...), по табличным

данным [14]

х-

У

*

9

**

А

а

у

/

✓ г

/

/

У /

V

Л*' У*

б

у /

/ г

* г

/

>

г/

У

л* /

у

у'

л--'

Р 0,9

0,7

0,5

0,3

0,1

Р 0,9

0,7

0,5

0,3

0,1

Р 0,9

0,7 0,5 0,3 ОД

0 10 20 30 40 50 60 70 о, град в

Рис. 2. Распределение давления Р на эллипсоиде с соотношением полуосей Ыа = 3/2 при числе Маха М = 4 (а), М = 6 (б), М = 20 (в), полученное

по формуле Ньютона (---), по формуле (7) (...), по табличным данным

[14] (...)'

у*

У 'У СГ*

А * У

У ✓ г"

/

(У

у

* 3>

«о»

а

..г» -Г"*

У* у .. г

4 у

/ ✓ Г*

/

/

б

/ г

/ ✓

✓

✓

/

У 'у

✓

•••• *

Р 0,9

0,7

0,5

0,3

19 29 39 49 59 69 о, град а

Р 0,9

0,7 0,5 0,3

35 45 55 65 75 а, град б

Р 0,9

0,7 0,5 0,3

35 45 55 65 75 а, град в

Р 0,9

0,7

0,5

0,3

35 45 55 65 75 а, град г

Рис. 3. Распределение давления Р на параболоидах с параметром р = 0,25 при числе Маха М = 2 (а) и с параметром р = 1 при числе Маха M = 4 (б), М = 6 (в),

М = 8 (г), полученное по формуле Ньютона (---), по формуле (7) (...),

по табличным данным [14]

К- - *****

•Л

«•V

< г*

г--"" .—1

с

> •г"

/

Р 0,9

0,7

0,5

0,3

г-г

Г*

А

35

45

55

65

75

а, град

Р 0,9

0,7

0,5

0,3

У

У

35

45

55

65

75

с, град

б

Рис. 4. Распределение давления Р на параболоидах с параметром р = 1 при

числе Маха М = 10 (а) и М = 20 (б), полученное по формуле Ньютона (---),

по формуле

Эти формулы можно применять и в трехмерных течениях, используя правило «местного тела», когда давление считается одинаковым в точках с одинаковыми углами встречи потока с поверхностью. Опробуем этот метод для расчета давления в различных сечениях параболоида вращения г = 1,125 г2 (рис. 5), обтекаемого сверхзвуковым потоком идеального газа под ненулевым углом атаки. В соответствии с правилом «местных тел» давление на параболоиде в трехмерном течении определяется по формуле

С

Р =

1 -

У-1 у + 3

с

агссоБ

-к

V

У(с)т=т

о» - ж/2

1 +

У-1 у + 3

агссоз(Г, Г,, я)-71 о - ж/2

У У

где Уао — вектор скорости набегающего потока; п — внешняя единичная нормаль к поверхности параболоида. Угол а* описывает положение звуковой точки на параболоиде в осесимметричном течении.

2

Результаты расчета давления на указанном параболоиде, обтекаемого под углами атаки а = 5° и а = 10°, в различных сечениях z = const, а также табличные данные [14] в зависимости от центрального угла 3 представлены на рис. 6 (а = 5°) и на рис. 7 (а = 10°). На графиках видно, что в сечениях z < 5 метод «местных тел» позволяет более точно рассчитать давление по сравнению с формулой Ньютона. Метод Ньютона дает максимальную относительную погрешность в указанных сечениях 20 %, а формула «местных тел» — 12 % для а = 10°, z = 4 (см. рис. 7, в), в остальных сечениях — не более 7 %.

В качестве «эквивалентного тела» можно рассматривать сферу. Тогда формула для расчета давления упрощается. На некоторых из этих графиков приведено распределение давления на параболоиде, вычисленного по методу «местных сфер», когда давление на теле в

точке, где угол arccos^-F^/F^, л) совпадает с углом ^-о на сфере,

полагалось равным давлению на сфере с данным углом встречи потока с поверхностью:

(

P =

1 -

Y-1 Y + 3

arccos^/F^, и)-я с.. - л/2

У V

1 +

Y-1 Y + 3

arccos

- л/2

- Л

У У

Анализ результатов показывает, что расчет давления по правилу «местных сфер» также дает более точные результаты по сравнению с формулой Ньютона, однако несколько уступает по точности (« 5 %) расчетам по правилу «местных тел».

Р 0,9

0,7 0,5

а

Р 0,9

0,7 0,5

б

Р ( 0,9

0,7 0,5

в

Р 0,9

0,7 0,5

0 20 40 60 80 100 120 140 3,град

г

Рис. 6. Распределение давления Р на параболоидах, обтекаемых под углом а = 5°, в сечениях г = 1 (а), г = 2 (б), г = 5 (в), г = 10 (г) при числе Маха

М = 4, полученное по формуле Ньютона (---), по формуле (11) (»),

по табличным данным [14]

'41 > _

- - — ' г-.,

* «

ч г,.,

N _ <

ч

ч ••1 •ч г-.. мч

•ч

41 Гв#.

ч' >

ч 1

••• 'Ч

— - " -

ч > Г».

V

-Л • -—»

Р 0,9

0,7 0,5 0,3

Р 0,9

0,7 0,5 0,3

Р 0,9

0,7 0,5 0,3

в

Р 0,9

0,7 0,5 0,3

0 20 40 60 80 100 120 140 9, град

г

Рис. 7. Распределение давления Р на параболоидах, обтекаемых под углом а = 10°, в сечениях г = 2 (а), г = 3 (б), г = 4 (в), г = 5 (г) при числе Маха

М = 4, полученное по формуле Ньютона (---), по формуле (11) (—),

по формуле

**ч V

ч к

\ > *>>

ч >

ч V

Ч . —.. ....

»»» . • - - - -

ч

«

N ч,

V \ 1* *ч 1.

ч. > Ч

ч Ч "•V

\

' - - "и

Вывод. Предложенные в данной работе формулы (7) и (11) в большинстве случаев дают лучшие результаты, чем метод Ньютона. Как и ожидалось, при наличии угла атаки с наветренной стороны формула Ньютона применима только тогда, когда угол встречи потока с поверхностью тела достаточно большой. Таким образом, данный метод может быть использован для получения начального распределения, а также для быстрой оценки давления на поверхности выпуклых затупленных тел вращения. В дальнейшем предполагается опробовать этот метод для более широкого класса поверхностей.

ЛИТЕРАТУРА

[1] Ковалев В.А. Гетерогенные каталитические процессы в аэротермодинамике. Москва, Физматлит, 2002, 224 с.

[2] Белоцерковский О.М., Андрущенко В.А., Шевелев Ю.Д. Динамика пространственных вихревых течений в неоднородной атмосфере. Москва, ЯНУС-К, 2000, 456 с.

[3] Котенев В.П., Сахаров В.И., Тирский Г.А. О расчете сверхзвукового обтекания затупленных тел химически неравновесным потоком газа. Журнал вычислительной математики и математической физики, 1987, т. 27, № 6, с. 411-415.

[4] Gross A., Fasel H.F. High-Order-Accurate Numerical Method for Complex Flow. AIAA Journal, 2008, vol. 46, no. 1, pp. 204-214.

[5] Черный Г.Г. Газовая динамика. Москва, Наука, 1986, 424 с.

[6] Калугин В.Т. Аэрогазодинамика органов управления полетом летательных аппаратов. Москва, Изд-во МГТУ им. Н.Э. Баумана, 2004, 688 с.

[7] Котенев В.П. Уравнения двумерных течений газа в динамических переменных. Информационные технологии, 2007, № 1, с. 37-41.

[8] Котенев В.П. Приближенный метод расчета давления на поверхности затупленных тел. Информационные технологии, 2007, № 2, с. 76-80.

[9] Димитриенко Ю.И., Котенев В.П., Захаров А.А. Метод ленточных адаптивных сеток для численного моделирования в газовой динамике. Москва, Физматлит, 2011, 280 с.

[10] Лунев В.В. Течение реальных газов с большими скоростями. Москва, Физматлит, 2007, 760 с.

[11] Котенев В. П. Определение положения звуковой точки на поверхности выпуклого затупленного тела. Вестник МГТУ им. Н.Э. Баумана. Сер. Естествознание, 2011, спец. вып. «Математическое моделирование», с. 150-153.

[12] Котенев В.П., Сысенко В.А. Уточненный метод быстрой оценки давления на поверхности гладких затупленных тел. Вестник МГТУ им. Н.Э. Баумана,, 2012, № 3, с. 64-75.

[13] Котенев В.П., Сысенко В.А. Метод быстрой оценки параметров на поверхности затупленных тел, обтекаемых сверхзвуковым потоком, с учетом равновесных физико-химических превращений. Инженерный журнал: наука и инновации, 2013, вып. 7. URL: http://engjournal.ru/ catalog/mathmodel/aero/

[14] Любимов А.Н., Русанов В.В. Течения газа около тупых тел. Москва, Наука, 1970, 379 с.

[15] Белоцерковский О.М. Расчет обтекания осесимметричных тел с отошедшей ударной волной: расчетные формулы и таблицы полей течений. Москва, ВЦ АН СССР, 1961, 56 с.

Статья поступила в редакцию 06.03.2014

Ссылку на эту статью просим оформлять следующим образом: Котенев В.П., Сысенко В.А. Аналитические формулы повышенной точности для расчета распределения давления на поверхности выпуклых затупленных тел вращения произвольного очертания. Математическое моделирование и численные методы, 2014, № 1, с. 68-81.

Котенев Владимир Пантелеевич родился в 1956 г., окончил МГУ им. М.В. Ломоносова в 1978 г. Д-р техн. наук, профессор кафедры вычислительной математики и математической физики МГТУ им. Н.Э. Баумана, начальник отдела аэродинамики в ВПК «НПО машиностроения». Автор более 40 научных работ в области прикладной математики, численных и аналитических методов исследования течения газа при обтекании поверхности летательных аппаратов. e-mail: rotvp@mail.ru

Сысенко Валентина Алексеевна родилась в 1978 г., окончила МГТУ им. Н.Э. Баумана в 2001 г. Канд. техн. наук, старший научный сотрудник в отделе аэродинамики в ВПК «НПО машиностроения». Автор 12 научных работ в области прикладной математики. e-mail: dv-sys@yandex.ru

Analytical formula with improved accuracy for calculating pressure distribution on the surface of convex, blunt rotation bodies of arbitrary shape

© V.P. Kotenev1, V.A. Sysenko2

1 Bauman Moscow State Technical University, Moscow, 105005, Russia

2 JSC MIC NPO Mashinostroyenia, Reutov, 143966, Russia

The authors developed the analytical formula for fast and accurate calculation of pressure on the surface of rotation bodies with arbitrary shape, which were flown by supersonic gas. The paper provides examples of applying the method for three-dimensional flows of gas.

Keywords: supersonic gas flow, three-dimensional gas flows, sonic point.

REFERENCES

[1] Kovalev V.A. Geterogennye kataliticheskie protsessy v aerotermodinamike [Heterogeneous catalytic processes in aerothermodynamics]. Moscow, Fizmatlit Publ., 2002, 224 p.

[2] Belotserkovsky O.M., Andrushchenko V.A., Shevelev Yu.D. Dinamika prostranstvennykh vikhrevykh techenii v neodnorodnoi atmosphere [Dynamics of spatial vortex flows in an inhomogeneous atmosphere]. Moscow, YaNUS-K Publ., 2000, 456 p.

[3] Kotenev V.P., Sakharov V.I., Tirsky G.A. Zhurnal vychislitel'noi matematiki i matematicheskoi fiziki — J. Computational Mathematics and Mathematical Physics, 1987, vol. 27, no. 6, pp. 411-415.

[4] Gross A., Fasel H. F. AIAA J.., vol. 46, no. 1, January 2008, pp. 204-214.

[5] Chernyi G.G. Gazovaya dinamika [Gas dynamics]. Moscow, Nauka Publ., 1986, 424 p.

[6] Kalugin V.T. Aerogazodinamika organov upravleniya poletom letatel'nykh apparatov [Aerogasdynamics of aircraft flight controls]. Moscow, Bauman MSTU Publ., 2004, 688 p.

[7] Kotenev V.P. Informatsionnye tekhnologii — Information technologies, 2007, no. 1, pp. 37-41.

[8] Kotenev V.P. Informatsionnye tekhnologii — Information technologies, 2007, no. 2, pp. 76-80.

[9] Dimitrienko Yu.I., Kotenev V.P., Zakharov A.A. Metod lentochnykh adaptivnykh setok dlya chislennogo modelirovaniya v gazovoi dinamike [Adaptive mesh tape method for numerical simulation of gas dynamics]. Moscow, Fizmatlit Publ., 2011, 280 p.

[10] Lunev V.V. Techenie real'nykh gazov s bol'shimi skorostyami [Real gases flow at high speeds]. Moscow, Fizmatlit Publ., 2007, 760 p.

[11] Kotenev V.P. Vestnik MGTU im. N.E. Baumana. Seriya Estestvennye nauki. Spetsial'nyi vypusk «Matematicheskoe modelirovanie» — Herald of Bauman Moscow State Technical University, Natural Science Series, 2011, special issue "Mathematical Modeling", pp. 150-153.

[12] Kotenev V.P., Sysenko V.A. Vestnik MGTU im. N.E. Baumana. Seriya — Herald of Bauman Moscow State Technical University, Series , 2012, no. 3, pp. 64-75.

[13] Kotenev V.P., Sysenko V.A. Inzhenernyi zhurnal: nauka i innovatsii — Engineering Journal: Science and Innovation, 2013, issue 7. Available at: http://engjournal.ru/catalog/mathmodel/aero/

[14] Lyubimov A.., Rusanov V.V. Techeniya gaza okolo tupykh tel [Gas flow around blunt bodies]. Moscow, Nauka Publ., 1970, 379 p.

[15] Belotserkovsky O.M. Raschet obtekaniya osesimmetrichnykh tel s otoshedshei udarnoi volnoi. (Raschetnye formuly i tablitsy polei techeniya). [Calculation of flow axial-symmetric bodies with detached shock wave. (Calculation formulas and tables of the flow fields)]. Moscow, CC AS USSR, 1961, 56 p.

Kotenev V.P. (b. 1956) graduated from Lomonosov Moscow State University in 1978. Doctor of Science (Technical), head of research department of Scientific-Industrial Association mashinostroenie. Author of about 40 publication in the field of applied mathematics, numerical and analytical methods to investigate gas streams near by the surface of flying vehicles. e-mail: rotvp@mail.ru

Sysenko V.A (b. 1978) graduated from Bauman Moscow State University in 2001. Candidate of Science (Technical), worker of Scientific-Industrial Association mashinostroenie. Author of 12 publication in the field of applied mathematics. e-mail: dv-sys@yandex.ru