CC BY
19
44
ВМУ. Серия 3. ФИЗИКА. АСТРОНОМИЯ. 2017. № 6
РАДИОФИЗИКА, ЭЛЕКТРОНИКА, АКУСТИКА
Показатель автомодельности сильных сходящихся цилиндрических ударных волн в газе с однородной плотностью
, Н. Н. Сысоев1, С. А. Шутеев1,а, С. Т. Белякин2,6
Московский государственный университет имени М. В. Ломоносова, физический факультет, 1 Центр гидрофизических исследований; 2 кафедра общей физики. Россия, 119991, Москва, Ленинские горы, д. 1, стр. 2. E-mail: a shuteev@phys.msu.ru, 6 nesu@phys.msu.ru Статья поступила 05.10.2016, подписана в печать 14.10.2016.
Предложена модель схождения сильных цилиндрических ударных волн (УВ) в газе с однородной плотностью. Дифференциальные уравнения с частными производными этой модели сведены к обыкновенным дифференциальным уравнениям, из которых определен закон схождения таких УВ и зависимость их показателя автомодельности а от эффективных показателей адиабаты газа перед (y) и за (7eff) фронтами УВ: а = f(7,Yeff). Показано, что эта зависимость для цилиндрических УВ в пределах ошибки измерений согласуется с опытными данными.
Ключевые слова: сходящиеся сильные цилиндрические ударные волны в газе, показатель автомодельности, дифференциальное уравнение. УДК: 537.523. PACS: 52.35.Tc.
Юсупалиев У.
Введение
Исследования УВ интересны как с точки зрения фундаментальной науки, так и с точки зрения практических задач, связанных, например, с термоядерным синтезом или медицинскими приложениями (раздробление почечных камней и др.). Впервые УВ исследовал Гудерлей [1] (в 1942 г.), а затем Ландау и Станюкович [2, 3] (в 1945 г.). Закон схождения УВ Гудерлеем был задан в виде
R(t)= R0(1 - t/t0)a, а Ландау и Станюковичем — в виде R(t)= Asw(-t)a,
(1)
(2)
где а — показатель автомодельности, Л0 — начальный радиус, Аэ-да — размерный параметр, в (1) и в (2) моменты кумуляции волн равны ¡0 и 0 соответственно. Система уравнений задачи состояла из уравнений непрерывности, Эйлера и адиабатично-сти, решения которых Ландау и Станюкович искали в виде
о = О-у(а р=роО(£), с2 = (а-)22(а (3)
где £ = г/Я(£) — автомодельная переменная, г — координата в цилиндрической системе координат, V — радиальная скорость газа относительно системы координат, связанной с неподвижным газом, находящимся в объеме цилиндра с радиусом Л, р — плотность газа, с — скорость звука в газе, р0 — плотность газа перед фронтом УВ. Искомые величины V, р и с относятся к области г ~ Л за фронтом УВ вблизи оси цилиндра.
В указанных работах численным методом были определены значения величины а лишь для некоторых значений показателя адиабаты для невозмущенного газа 7, а именно 1.4 [1-3], 1, 1.67, 3 и то [2] (таблица). Начиная с этих работ определение значений величины а является одной из основных задач теоретического и экспериментального исследований закономерностей распространения УВ в газе с однородной плотностью.
Впоследствии задача Гудерлея-Ландау-Станю-ковича была исследована различными приближенными аналитическими и численными методами авторами работ [4-14]. Основные результаты указанных работ состояли в определении значения величины а для широкого диапазона величины 7 (см. таблицу).
Задача о сильных сходящихся УВ была также рассмотрена приближенными методами [15-17]. На основе анализа большого числа частных и предельных случаев для величины а сходящихся УВ предложено следующее приближенное выражение [18]:
n + Г
n = 1 + 2 +
7
П7~ V7-T
(4)
Экспериментальному определению величины а для сходящихся цилиндрических УВ посвящены работы [19-24]. В них было показано, что для такой УВ значение а в пределах ошибки измерений совпадает с предсказываемой авторами работ [1-3] величиной.
Несмотря на давнюю историю исследований сходящихся УВ в газе, до сих пор не определена аналитическая зависимость а = [ (7,7ей). Понятно,
n
а=
Значения коэффициента автомодельности а для сильных сходящихся цилиндрических ударных волн в газе с однородной плотностью
Y ^ 1 y = 1.4 y = 1.67 y = 3.0 Y ^
Guderley (1942) [1] 0.835
Landau, Stanyukovich (1955) [2, 3] 1.0 0.834 0.83 0.81 0.5
Butler (1954) [4] 0.835
Lazarus, Richtmyer (1977) [10] y = 1.1-6.0 0.885-0.75 Y = 10 y = 50 0.74 0.73 y = 100 0.73 0.727
Fujimoto, Mishkin (1978) [5] 0.828
Chester (1954), Whitham (1957) [15, 16] 0.835 0.815 0.772 0.707
Chisnell (1998) [18] 1.0 0.83 0.815 0.773 0.707
Nakamura (1983) [6] 0.834
Valiyev, Kraiko (2011) [14] y = 1.001-1.91 0.980-0.804 y = 2.37-3.0 0.789-0.776
По формуле (25) 1.0 0.857-0.87 0.80-0.82 0.67-0.70 0.5
что иметь дело с такой зависимостью, согласующейся с опытными данными, всегда предпочтительнее, так как она может помочь установить новые, ранее неизвестные закономерности исследуемого явления. В связи с этим данное сообщение посвящено установлению именно такой зависимости и сравнению ее с опытными данными для сходящихся цилиндрических УВ.
1. Экспериментальное определение значений величины а для сходящейся цилиндрической УВ
УВ генерировалась с помощью электрического взрыва проводников — импульсного электрического разряда. В качестве проводников использовалась металлизированная лавсановая пленка толщиной 5 мкм и шириной 0.15 м. Эта пленка закреплялась на внутренней поверхности текстолитовой трубы (толщиной 0.01 м) с внутренним диаметром 2Л0 = 0.1 м. Использовался как диаметрально расположенный подвод электрической энергии, так и односторонний, когда оба подводящих электрода расположены с одной стороны кольца на расстоянии 2 мм. Второй вариант расположения электродов оказался предпочтительным, так как позволял генерировать ударную волну необходимой формы. Такое разрядное устройство помещалось в цилиндрическую камеру из прозрачного огрстекла диаметром 0.2 м, имеющую систему откачки и напуска газа (Аг, СО2 (7 = 1.3), воздух при давлении р0 = 105 Па). С помощью такого устройства удалось генерировать УВ с числом Маха М = 6-15.
Генератор УВ в качестве полезной нагрузки входил в состав разрядного контура, который имел батарею конденсаторов емкости С = 30 мкФ (с зарядным напряжением и0 = 5-25 кВ) и управляющий воздушный разрядник. Для регистрации характе-
ристик УВ была использована теневая установка. При этом для достижения необходимого пространственно-временного разрешения УВ (0.5 мм, 10-7 с) применялся сверхскоростной фоторегистратор СФР-2М, работающий как в покадровом режиме, так и в режиме фоторазвертки во времени. Работа разряда и фоторегистратора была синхронизирована между собой с помощью блока управления. Изучалась динамика схождения цилиндрической УВ к ее оси симметрии (0z). Это достигалось регистрацией теневой картины схождения УВ в перпендикулярном к оси направлении.
Предварительные опыты показали, что эффект кумуляции сходящейся УВ проявляется в относительно малой области (10-20 мм) в окрестности оси 0z. Поэтому для обнаружения усиления сходящейся УВ особое внимание уделялось пространственно-временному разрешению волны именно в этой области.
На рис. 1 приведена последовательность кадров фоторегистрации процесса схождения фронта цилиндрической УВ. Следует подчеркнуть, что процесс схождения фронта УВ сильно зависит от начальных возмущений: наблюдаемое нами симметричное ее схождение имеет место не во всех случаях.
Из рисунка определялась зависимость фронта УВ R(t) для различных режимов ее генератора: начальное число Маха M изменялось путем варьирования напряжения U0. Для вычисления значений показателя а из зависимости R(t) для различных чисел М за начало отсчета времени t= 0 (следуя [1-3]) был принят момент кумуляции УВ.
Экспериментальная зависимость а = f (7, Yeff) приведена на рис. 2, где приведены также опытные данные работ [21, 22, 24].
ВМУ. Серия 3. ФИЗИКА. АСТРОНОМИЯ. 2017. № 6 Фронт ударной волны
ооооо
о
г
Рис. 1. Процесс схождения цилиндрической ударной волны к оси симметрии 0г и отражения от нее. Момент кумуляции УВ на оси 0г принят за t = 0. Длительность экспозиции кадра сверхскоростного фоторегистратора равна и время между метками равны 2 мкс. Радиус кольца Л0 = 5 см
а 1.00 0.95 0.90 0.85 0.80 0.75 0.70
А-[21] ■-[22] 0 — [24]
• — Данные нашей работы у — Данные нашей работы (слабо ионизированный газ)
Решение системы уравнений (5)-(7) будем искать, в отличие от (3), в виде
V = ЩЫО, р = р0ё (£),
(9)
где в качестве масштаба плотности выбрана плотность газа р0 перед фронтом УВ, а масштаб скорости Я^) = В отличие от работ [1-3], в данной постановке задачи радиус Я^) фронта УВ является искомой функцией.
Подставляя (8) и (9) в уравнения (5)-(7) и приведя их к безразмерному виду, получим следующую систему:
1.0
1.5
2.0 У
Рис. 2. Экспериментальные и теоретические зависимости коэффициента автомодельности а сильных сходящихся цилиндрических волн от показателя адиабаты 7: 1 — построена по формулам (4), 2 — построена по формуле (25)
2. Модель схождения УВ
Для установления зависимости а = / (7, ) используем те же уравнения, что и в работах [1-3]:
ё (и - £) + и' + и = 0, (10)
ЯЛ - ¿и + 7е„и' +Оиё = 0, (11)
ЛЛ - ¿и + и' (£ - и) = 0, (12)
в которой штрих означает дифференцирование по £, а точка — по t.
Используя (10) и разделяя переменные в (11) и (12), получим систему уравнений с константой разделения С:
др др до ро -£ + ^ + р— + ?-= 0, дt дг дг г
до д V = др
д +0 дГ = - РдГ'
— 1п рр"7ей + V — 1п рр-7ей = 0,
(5)
(6) (7)
ля
Л
= с,
си' ' п и' + и
£ и - 7еИи - ^^
= с,
где р — давление газа. Воспользуемся одним из основных свойств сильных УВ — равенством внутренней энергии и кинетической энергии единицы массы газа за их фронтами [25], откуда для давления получим
(7еИ - 1^2
£ и - и' - п{и' + и) = С.
(13)
(14)
(15)
Граничное условие для функции и(£) при £ = 1 формулируется следующим образом:
и(1) =
2
Л т^
V уму'
(16)
р =
= п pv .
(8)
7 + 1 ^ 7М2 которое определено из формулы, полученной в [3].
Решением уравнения (13) при C = 1с начальным условием
Я(0) = Й0, = 0, (17)
будет функция
R(t)= Ro 1 --
t
1/(1—С)
. to) ■
А при С = 1с другим начальным условием
R(0)=0, (-«(0) = 0,
(18)
(19)
уравнение (11) имеет следующие решения: при t < 0
= Asw • Ы)1/(1-С), (20)
а при t ^ 0
Asw • ^/(1-С). (21)
Решения (18) и (20) с показателем автомодельности
1 (22)
1С
совпадают с законами (1) (самоподобное решение Гудерлея) и (2) (самоподобное решение Ландау-Станюковича) в газе с однородной плотностью соответственно.
Константу С в (13) можно получить из (15), зная ы(1) и и'(1). Последнее найдем из (14) и (15) при £= 1:
2 (1 - М )(т - 3 + М) ",(1)=^ (; +1)2 . (23)
Далее, получив выражение для константы С, с учетом (20) определим искомое выражение для величины а:
= 1-
Y + 1
Y + 1 - 2(l - M) (2 + +
+ 4 1-
YeffA 1 + Yeff yM2 J Y + 1
г
(24)
Для сильных сходящихся цилиндрических УВ (М ^ 1) из (24) следует, что выражение для а = [(т, 7ей) имеет следующий вид:
(7 + 1)2
1
2 y + 1 +(Yeff + 1)(y - 1)'
(25)
3. Сравнение результатов модели с опытными данными
На рис. 2 представлены кривая 1, построенная по формуле (4), и кривая 2, построенная по формуле (25) при 7ей = 0.97 (учет ионизации газа за фронтом сильной УВ). На рисунке, кроме наших
опытных данных, также приведены экспериментальные данные работ [21, 22, 24]. Видно, что кривая 2 согласуется с опытными данными для сходящихся сильных цилиндрических УВ, созданных различными способами, чего нельзя сказать о приближенной формуле (4).
Заключение
Таким образом, на основе предложенной нами модели установлена ранее неизвестная зависимость а = f (y, Yeff) для сходящихся сильных цилиндрических УВ, согласующаяся с опытными данными.
Список литературы
1. Guderley G. Luftfahrtforschung. 1942. 19(9). P. 302.
2. Ландау Л.Д., Станюкович К.П. // ДАН СССР. 1945. 46. № 9. С. 399.
3. Stanyukovich K. Unsteady Motion of Continous Media. Oxford: Pergamon, 1960.
4. Butler D. Converging Spherical and Cylindrical Shocks. Report N 54. N.Y.: Burgess Hill, 1954.
5. Fujimoto Y., Mishkin E. // Phys. Fluids. 1978. 21. P. 1933.
6. Nakamura Y. // Phys. Fluids. 1983. 26. P. 1234.
7. Hafner P. // J. Appl. Math. 1988. 48. P. 1244.
8. Van DykeM., Guttman A. // J. Fluid Mech. 1982. 120. P. 451.
9. Ponchaut N., Hornung H.G., Mouton D.I. // J. Fluid. Mech. 2006. 560. P. 103.
10. Lazarus R., Richtmyer R. // Los Alamos Scientific Laboratory of the University of California. Los Alamos, NM, 1977.
11. Lazarus R. // SIAM J. Numer. Anal. 1981. 18. P. 316.
12. Ponchaut N., Hornung H.G., Mouton D.I. // J. Fluid. Mech. 2006. 560. P. 103.
13. Hornung H.G., Pullin D.I., Ponchaut N. // Acta Mech. 2008. 201. P. 31.
14. Валиев Х.Ф., Крайко А.Н. // ПММ. 2011. 75(2). С. 314. (Valiyev Kh.F, Kraiko A.N. // J. Appl. Math. Mech. 2011. 75(2). P. 218.)
15. Chester W. // Philos. Mag. 1954. 45. P. 1293.
16. Chisnell R. // J. Fluid Mech. 1957. 2. P. 286.
17. Whitham G. // J. Fluid Mech. 1957. 2. P. 145.
18. Chisnell R. // J. Fluid Mech. 1998. 354. P. 357.
19. Matsuo H., Nakamura Y. // J. Appl. Phys. 1980. 51. P. 3126.
20. Matsuo H, Nakamura Y. // J. Appl. Phys. 1981. 52. P. 4503.
21. Baronets P. // Fluid Dynam. 1984. 19. P. 503.
22. Baronets P. // Fluid Dynam. 1994, 29. P. 129.
23. Takayama K., Kleine H., Gronig H. // Exp. Fluids. 1987. 5. P. 315.
24. Kjellander M., Tillmark N., Apazidisa N. // Phys. Fluids. 2011. 23. P. 116103.
25. Чу К., Гросс Р. Ударные волны в физике плазмы // Физика высокотемпературной плазмы. М.: Мир, 1972. С. 262.
v
48
ВМУ. Серия 3. ФИЗИКА. АСТРОНОМИЯ. 2017. № 6
The self-similarity index of the convergence of strong cylindrical shock waves in a gas with a uniform density
, N.N. Sysoev1, S.A. Shuteev1a, S.T. Belyakin2b
1 Gydrophysical Research Center; 2 Department of General Physics,
Faculty of Physics, Lomonosov Moscow State University. Moscow 119991, Russia.
E-mail: a shuteev@phys.msu.ru, b nesu@phys.msu.ru.
A model of the convergence of cylindrical shock waves (SWs) in a gas with a uniform density has been considered. The partial differential equations of this model have been reduced to ordinary differential equations, from which the law of convergence of such shock waves and the dependence а = f(y, yeff of their self-similarity index а on the heat-capacity ratio in front of the shock wave (yeff) and behind the shock wave front (yeff) of the gas have been found. This dependence for cylindrical shock waves has been shown to agree with the experimental data within the measurement error.
Keywords: converging strong cylindrical shock waves in gas, self-similarity index, differential equation. PACS: 52.35.Tc. Received 5 October 2016.
English version: Moscow University Physics Bulletin. 2017. 72, No. 6. Pp. 539-543.
Сведения об авторах
1. Юсупалиев Усен — канд. физ.-мат. наук, ст. науч. сотрудник.
2. Сысоев Николай Николаевич — доктор физ.-мат. наук, профессор, директор; тел.: (495) 939-17-97, e-mail: nesu@phys.msu.ru.
3. Шутеев Сергей Александрович — кандидат физ.-мат. наук, ст. науч. сотрудник; тел.: (495) 939-17-97, e-mail: shuteev@phys.msu.ru.
4. Белякин Сергей Тихонович — инженер, тел.: (495) 939-51-56, e-mail: nesu@phys.msu.ru.
U. Yusupaliev
