Аналитическое решение задачи о сходящейся ударной волне в газе в одномерном случае Текст научной статьи по специальности «Математика»

УДК 533.6.011.1 DOI: 10.14529/mmph170407

АНАЛИТИЧЕСКОЕ РЕШЕНИЕ ЗАДАЧИ О СХОДЯЩЕЙСЯ УДАРНОЙ ВОЛНЕ В ГАЗЕ В ОДНОМЕРНОМ СЛУЧАЕ

В.Ф. Куропатенко1,2 , Ф.Г. Магазов , Е.С. Шестаковская2

Российский федеральный ядерный центр - Всероссийский научно-исследовательский институт технической физики им. академика Е.И. Забабахина, г. Снежинск, Российская Федерация 2 Южно-Уральский государственный университет, г. Челябинск, Российская Федерация E-mail: [email protected]

Построено аналитическое решение задачи о сходящейся ударной волне в сосуде с непроницаемой стенкой, описывающее случаи плоской, цилиндрической и сферической симметрии. На границе сосуда задана отрицательная скорость, а скорость холодного идеального газа равна нулю. В начальный момент времени из этой точки начнет распространяться ударная волна к центру симметрии. Граница сосуда будет двигаться по определенному закону, согласованному с движением ударной волны. В эйлеровых переменных она движется, но в лагранжевых переменных её траектория является вертикальной линией. Получены уравнения, определяющие структуру течения газа между фронтом ударной волны и границей как функции времени и лагранжевой координаты, а также зависимость энтропии от скорости ударной волны. Для всех случаев симметрии найдены показатели автомо-дельности и соответствующие им значения безразмерных координат для широкого диапазона показателей адиабаты. Задача решена в лагранжевых координатах и принципиально отличается от ранее известных постановок задачи о схождении автомодельной ударной волны к центру симметрии и её отражении от центра, которые построены для бесконечной области в эйлеровых координатах.

Ключевые слова: ударная волна; плоская симметрия; цилиндрическая симметрия; сферическая симметрия; идеальный газ; аналитическое решение.

Введение

Развитие теории размерности и подобия величин механики сплошной среды позволило в разные годы ХХ века получить автомодельные решения задачи о фокусировке ударной волны в бесконечном идеальном газе [1-9]. Задачи, имеющие аналитическое решение, еще называют «эталонными» задачами, так как они используются для верификации математических моделей динамических процессов механики сплошных сред. Однако реальные тела имеют конечные размеры. В данной работе рассмотрена задача о сходящейся ударной волне в сосуде с непроницаемой стенкой, имеющая точное аналитическое решение, описывающее три типа одномерных движений - плоские, цилиндрически и сферически симметричные движения.

Постановка задачи

Рассматривается сосуд с непроницаемой стенкой, в котором находится газ массой M0 и начальными при t = t0 параметрами газа р0 = const, U0 = 0, P0 = 0, E0 = 0, где p - плотность, U -скорость, P - давление, E - удельная внутренняя энергия. Лагранжевой координатой является масса M . Второй независимой переменной является время t. В точке t0, M0 задана скорость U1 < 0 . Таким образом, в этой точке задан сильный разрыв, который при t > t0 распространяется к центру симметрии и в момент f фокусируется в точку M = 0 . Граница сосуда при t > t0 движется в переменных r, t, но в переменных M, t её траектория является вертикальной линией. Вообще говоря, все траектории частиц являются вертикальными линиями, вдоль которых сохраняется то значение энтропии, которое возникло на ударной волне. Параметры газа между ударной волной и границей определяются системой законов сохранения Эйлера-Гельмгольца. Уравнение состояния используются в двух формах:

Куропатенко В.Ф., Магазов Ф.Г., Шестаковская Е.С.

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

Р = (у- 1)рЕ, Р = Р(я)рг, (1)

где Р (5 ) - функция от энтропии. Ударная волна

Законы сохранения на ударной волне (УВ) при и0 = 0, Р0 = 0, Е0 = 0, Р0 = 0 имеют вид [10]:

р^ (О - )-р0 О = 0, (2)

Р0Оик -Рк = 0, (3)

Р0О\Е„ + 2и2 I -Р^и^ = 0.

(4)

Индексом « w » обозначены величины на ударной волне, О - скорость ударной волны. Преобразуем эти уравнения к виду, содержащему зависимости им,, рм,, Рч>, Рм, от скорости ударной волны в лагранжевых координатах. Лагранжева координата Мм, ударной волны в одномерном случае связана с её эйлеровой координатой ^ уравнением

где / =

1, плоская симметрия

2, цилиндрическая симметрия, 5/ =

3, сферическая симметрия

Mw =^Р0<,

л

1, при т = 1 2р, при /л = 2 . 4р, при / = 3

(5)

Скорость ударной волны в лагранжевых координатах есть изменение М^, со временем

Ж = М = ГГО

(6)

Заменим эйлерову координату ударной волны её лагранжевой координатой. Для этого выразим ^ из (5) и подставим в (6)

* = (М Г^ (зр )1Л р. (7)

Выразив в (7) О через Ж и Мм! и подставив в (2)-(4), получим с помощью (1) величины на ударной волне, содержащие Ж и М^,:

и, =-

V +1

р =тггр0;

\-1/ л

у+1

(5/Р0 )-1/ '{/М^ )(1-т)/Ж;

Pw =-^р0М-2)/\-2/т (М )2(1-/)/л Ж2. у+1

Из (1), (8) и (10) следует выражение для РМ!

Рw =

2

(

V-1

Р0

(г-(т-2)/т) 5 -2/

(М )2(1~Л)/Ж2.

у+1\у+1

В точке t = ^ Mw = M0,Uw = и^,Рм, = Рм,0,Рw = Рм выражения (9)-(11) принимают вид:

^ = и 0

( Ж ^

\ Ж0

М,

\(т-1)/т

,Р = Р0

( Ж >

2 (М \2(т-1)/

\ Ж0)

Гш\2 г М ^2(m-1Vm

>К = Р0

Ж

\ Ж0 J

По аналогии с [11, 12] зададим траекторию ударной волны в виде

Mw = М0 р^)п,

(8)

(9)

(10)

(11)

. (12) (13)

где р = (^ -1)- t0). Продифференцировав Мм, по t, получим выражение для скорости

0

ударной волны в лагранжевых координатах

2

W = W0j

n—1

где

Wo =-

M 0n

tf - to

Исключив в (13) и (14) функцию от времени, получим зависимость W от Mw

W_ W

/ \(n—l]/n

V M 0 J

С помощью соотношения (16) исключим W в (12):

Uw = Uwo

(M \(n—m)lmn

V M o J

P = P

(M "\2(n—m)mn

wo

V Mo J

Fw = Fw0

(M ^2(n—m)mn

V Mo J

(14)

(15)

(16) (17)

Поскольку Мк и ^ связаны уравнением траектории ударной волны, то из (13) и (17) следуют

зависимости Uw, Pw, Fw от t:

w w w

и = и Л"-т)/т Р = Р {А"-т)1т р = р (А(18) Величина ^ вдоль траектории частицы с координатой М№ постоянна. Следовательно, зави-

симость энтропии от массы между ударной волной и границей газа имеет вид

F = F

( M ^2(n—m)mn

w0

Из (5) следует зависимость rw от Mw

V M0 J

( m

w

V M 0 J

(19)

(20)

Выбор траектории УВ в форме (13) определяет характер зависимостей W и D от Mw . При n < 1 W ® —¥ , если Mw ® 0, при n = 1 W = W0 = const и при n > 1 уменьшается при уменьшении Mw . Эти типы зависимостей приведены на рис. 1. Зависимость D от Mw следует из (7) и (20)

D

(M \(n—m)mn

D

V M 0 J

(21)

Из (21) следует, что при n <m D ®—¥, если Mw ® 0, при n = m D = D0 = const и при n > m D ® 0 при уменьшении Mw (рис. 2). Фактическая реализация одного из этих режимов возможна лишь при организации соответствующего течения газа между фронтом УВ и наружной границей.

Течение газа за ударной волной

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

— и = 0,

dt

M

dt

dU dt

d( rm—1u)

dr| + = 0,

M

+ s,,r

M

dM

.1 4ft)

dM

= 0.

(22)

(23)

(24)

Эти уравнения содержат три искомых функции г, р и и . Величина ^ определяется на ударной волне и зависит только от М (17).

rw = r0

Куропатенко В.Ф., Магазов Ф.Г., Шестаковская Е.С.

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

Перейдём в (22)-(24) к новым искомым функциям

Я = гт, С = гт-1и.

0

Ж0

р

0

О0

(25)

Рис. 1. Схематическое изображение возможных зависимостей W(Mw)

Рис. 2. Схематическое изображение возможных зависимостей й(М^)

После перехода к функциям Я и С уравнения (22)-(24) примут вид:

дЯ 1 -/С = 0,

д )м др1 2 дС п

— | + 5лр2-= 0,

дt )м т дМ '

'С + 5

М

Из (17), (19) и (25) следуют зависимости Я^ и СМ! от Мм,:

М,

2С 2 Я= 0.

К = Я0

М

Cw = С0

( м 1(п-1)/

V М0 У

(26)

(27)

(28) (29)

Уравнения (26)-(28) являются основными для отыскания Я, С и рв области интегрирования

Mw < М < М0, to < t < t/ .

Перейдём от переменных t, М к переменным t, X(t, М). С помощью уравнений для производных

д 1 Г дХ 1 ( д 1 Г д 1Г дХ

11 = (11 +(

дt

М

дt

преобразуем уравнения (26)-(28)

X

дЯ дt

V дХ У, V дt

М

М

\ дХ у, V дМ

+

(дЯ 1

X

V дХу?

§| = 0^

М

Эр дt

+

др 1 ( дХ

V да V Э,

+ /

(дС 1

М

V эХу,

К

дМ

=0,

'дС' . дt

(дС

Л

У(

дХ1 (л- 1С

дt

М

Я

+ 5/Я

2(м-1)/т

Р

( ЭР 1

^ | +ГРрг-1 дМ 1 7 И

др1( дХ

\ эху, V дм

(30)

(31) = 0. (32)

0

0

0

0

1

1

п

I

X

Зависимость М) зададим так, чтобы на ударной волне было Х = 1- Из (13) следует, что проще всего взять такую зависимость в виде

Х =М <Р~" - (33)

М0

Для разделения переменных представим Я, р и С в виде произведений функций от времени на функции от Х

Я = ак ^)Т(Х) р=ар^)5(Х) С = а (^(Х). (34)

Поскольку на ударной волне Х = 1, то значения Т1= Т(1), 81 = 8(1), Z1 = Z (1) должны быть постоянными. Из (29) и (33) следует зависимость Я^ (^)

К = ЯоРп - (35)

Сравнив эту зависимость с (34) на ударной волне, получим выражение для аЯ

ая (г) = ЯоРпТ1"1. (36)

Аналогично для ар и аС получаем соотношения

ар=Ро (у+1^ ас (() = Сорп-1 ZГ1. (37)

Подставив (34)-(37) в (30)-(32) и воспользовавшись (15), получим три уравнения для Т, 8 и

Z:

ХТ ' = А1, (38)

81B1Z '-Х1х8' = 0, (39)

-Х ^ + <СуХХ8' = С2, (40)

^ 81

где штрих означает дифференцирование по X. Коэффициенты уравнений (38)-(40) А1, В1, С1, С2 с помощью (5), (6), (12) и (23) преобразуются к виду

2ZT 282 8^-1 Т 2('и-1)/тх~((т-2)п+2т)/тп

А =Т-(7+1)Z1, В = (у-1)8]_2 ' С1 = 8у-1 Т12(т-1^т ,

2(т-1)Z2 Т1 (п -1)Z 2(п-т)8

m(r+1)Z2т nzl ~L mn d

Уравнения (38)-(40) образуют относительно T', S', Z' систему линейных неоднородных уравне-

С2 М(г+1)^т С1

-(40) образуют отн ний. Определитель системы равен

А = ВС1Х-Х2.

Если А Ф 0, то решение системы (37)-(39) существует и имеет вид

т' = А, 8'= В С 81, ^ = ХС^ . (42)

X А А

Интегрирование системы уравнений (42) начинается в точке Х= 1 (на ударной волне). Расчёты показывают, что существует промежуток значений п таких, что определитель в ноль не обращается. При некотором значении п* определитель обращается в ноль при Х = Х* ■ В этой точке решение существует, если С2 тоже обращается в ноль. Каждому значению у соответствует одно значение п* (таблицу). В этой же таблице приведены значения Х*, при которых одновременно

А(Х*) = 0, С2(Х*) = 0.

Куропатенко В.Ф., Магазов Ф.Г., Аналитическое решение задачи о сходящейся

Шестаковская Е.С. ударной волне в газе в одномерном случае

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

Значения n и X*, соответствующие различным значениям показателя адиабаты у

m у 1,1 1,2 4/3 1,4 5/3

1 n* 1,0 1,0 1,0 1,0 1,0

X* 4,691312 3,464715 2,828906 2,646222 2,236450

2 n* 1,770501 1,722331 1,684516 1,670651 1,631252

X* 6,768540 4,873946 3,896265 3,616019 2,990161

3 n* 2,387916 2,271434 2,183068 2,151532 2,065135

X* 7,959997 5,717071 4,559431 4,227062 3,481885

Статья выполнена при поддержке Правительства РФ (Постановление № 211 от 16.03.2013 г.), соглашение № 02.А03.21.0011.

Литература

1. Guderley, G. Starke kugelige und zylindrische Verdichtungsstobe in der Nahe des Kugelmittelpunktes bzw. der Zylinderachse / G. Guderley // Luftfartforschung. - 1942. - Т. 19, № 9. - С. 302312.

2. Седов, Л.И. О неустановившихся движениях сжимаемой жидкости / Л.И. Седов // Доклады Академии наук СССР. - 1945. - Т. 47, № 2. - С. 94-96.

3. Станюкович, К.П. Автомодельные решения уравнений гидромеханики, обладающих центральной симметрии / К.П. Станюкович // Доклады Академии наук СССР. - 1945. - Т. 48, № 5. -С.331-333.

4. Брушлинский, К.В. Об автомодельных решениях некоторых задач газовой динамики / К.В. Брушлинский, Я.М. Каждан // Успехи математических наук. - 1963. - Т. 18, № 2. - С. 3-23.

5. Седов, Л.И. Методы подобия и размерности в механике / Л.И. Седов. - Москва: Гостехиз-дат, 1954. - 328 с.

6. Сидоров, А.Ф. Процессы безударного конического сжатия и разлета газа / А.Ф. Сидоров, О.Б. Хайруллина // Прикладная математика и механика. - 1994. - Т. 58, № 4. - С. 81-92.

7. Сидоров, А.Ф. Об оптимальном безударном сжатии газовых слоев / А.Ф. Сидоров // Доклады Академии наук СССР. - 1990. - Т. 313, № 2. - С. 283-287.

8. Крайко, А.Н. Сферически и цилиндрически симметричное нестационарное сжатие идеального газа / А.Н. Крайко // VII Забабахинские научные чтения, Снежинск, 8-12 сентября 2003. http://www.vniitf.ru/rig/konfer/7zst/reports/s6/6-1.pdf

9. Крайко, А.Н. Быстрое цилиндрически и сферически симметричное сильное сжатие идеального газа / А.Н. Крайко // Прикладная математика и механика. - 2007. - Т. 71, № 5. - С. 744760.

10. Куропатенко, В.Ф. Модели механики сплошных сред / В.Ф. Куропатенко. - Челябинск: Челябинский государственный университет, 2007. - 302 с.

11. Куропатенко, В.Ф. Динамическое сжатие холодного газового шара / В.Ф. Куропатенко, Е.С. Шестаковская, М.Н. Якимова // Доклады академии наук. - 2015. - Т. 461, № 5. - С. 530-532.

12. Kuropatenko, V.F. Analytical solution of the problem of a shock wave in the collapsing gas in Lagrangian coordinates. / V.F. Kuropatenko, E.S. Shestakovskaya // AIP Conference Proceedings. -2016. - Vol. 1770, Issue 1. - P. 030069.

Поступила в редакцию 18 сентября 2017 г.

Bulletin of the South Ural State University Series "Mathematics. Mechanics. Physics" _2017, vol. 9, no. 4, pp. 52-58

DOI: 10.14529/mmph170407

ANALYTICAL SOLUTION OF THE PROBLEM OF A CONVERGENT SHOCK IN GAS FOR ONE-DIMENSIONAL CASE

V.F. Kuropatenko1'2 , F.G. Magazov2, E.S. Shestakovskaya

2

1 Russian Federal Nuclear Center - Zababakhin All-Russian Scientific Research Institute of Technical Physics, Snezhinsk, Russian Federation

2 South Ural State University, Chelyabinsk, Russian Federation E-mail: [email protected]

The analytical solution of the problem of a convergent shock in the vessel with an impermeable wall is constructed for the cases of planar, cylindrical and spherical symmetry. The negative velocity is set at the vessel boundary. The velocity of cold ideal gas is zero. At the initial time the shock spreads from this point into the center of symmetry. The boundary moves under the particular law which conforms to the movement of the shock. In Euler variables it moves but in Lagrange variables its trajectory is a vertical line. Equations that determine the structure of the gas flow between the shock front and the boundary as a function of time and the Lagrange coordinate as well as the dependence of the entropy on the shock wave velocity are obtained. Self-similar coefficients and corresponding critical values of self-similar coordinates were found for a wide range of adiabatic index. Thus, the problem is solved for Lagrange coordinates. It is fundamentally different from previously known formulations of the problem of the self-convergence of the self-similar shock to the center of symmetry and its reflection from the center which has been constructed for the infinite area in Euler coordinates.

Keywords: shock wave; planar symmetry; cylindrical symmetry; spherical symmetry; ideal gas; analytical solution.

References

1. Guderley G. Starke kugelige und zylindrische Verdichtungsstobe in der Nahe des Kugelmittelpunktes bzw. der Zylinderachse. Luftfartforschung, 1942, Vol. 19, no. 9, pp. 302-312.

2. Sedov L.I. Doklady Akademii nauk SSSR, 1945, Vol. 47, no. 2, pp. 94-96. (in Russ.).

3. Stanyukovich K.P. Doklady Akademii nauk SSSR, 1945, Vol. 48, no. 5, pp. 331-333. (in Russ.).

4. Brushlinskiy K.V., Kazhdan Ya.M. Russian Mathematical Surveys, 1963, Vol. 18, no. 2, pp. 1-22. DOI: 10.1070/RM1963v018n02ABEH001133

5. Sedov L.I. Metody podobiya i razmernosti v mekhanike. Moscow, Gostekhizdat Publ., 1954, 328 p. (in Russ.).

6. Sidorov A.F., Khayrullina O.B. Prikladnaya matematika i mekhanika, 1994, Vol. 58, no. 4, pp. 81-92. (in Russ.).

7. Sidorov A.F. Doklady Akademii nauk SSSR, 1990, Vol. 313, no. 2, pp. 283-287. (in Russ.).

8. Krayko A.N. Sfericheski i tsilindricheski simmetrichnoe nestatsionarnoe szhatie ideal'nogo gaza (Spherically and cylindrically symmetric nonstationary compression of an ideal gas). VII Zaba-bakhinskie nauchnye chteniya (Proc. VII Zababakhin Scientific Talks Conference), Snezhinsk, 2003. http://www.vniitf.ru/rig/konfer/7zst/reports/s6/6-1.pdf (in Russ.).

9. Krayko A.N. Prikladnaya matematika i mekhanika, 2007, Vol. 71, no. 5, pp. 744-760. (in Russ.).

10. Kuropatenko V.F. Modeli mekhaniki sploshnykh sred [Models of continuum mechanics]. Chelyabinsk: Chelyabinskiy gosudarstvennyy universitet Publ., 2007, 302 p. (in Russ.).

11. Kuropatenko V.F., Shestakovskaya E.S., Yakimova M.N. Dynamic compression of a cold gas sphere. Doklady Physics, 2015, Vol. 60, Issue 4, pp. 180-182. DOI: 10.1134/S1028335815040059

12. Kuropatenko V.F., Shestakovskaya E.S. Analytical solution of the problem of a shock wave in the collapsing gas in Lagrangian coordinates. AIP Conference Proceedings, 2016, Vol. 1770, Issue 1, p. 030069. DOI: 10.1063/1.4964011

Received September 18, 2017