CC BY
83
Вычислительные технологии
Том 2, № 2, 1997
РОЛЬ КАВИТАЦИОННЫХ ЭФФЕКТОВ В МЕХАНИЗМАХ РАЗРУШЕНИЯ И В КРУПНОМАСШТАБНЫХ ВЗРЫВНЫХ
ПРОЦЕССАХ**
В.К.Кедринский Институт гидpодинамики им. М. А. Лавpентьева СО РАН
Иовосибщск, Россия
В. А. Вшивков, Г. И. ДудниковА, Ю.И. Шокин Институт вычислительных технологий СО РАН Иовосибщск, Россия
Представлен краткий обзор экспериментальных результатов, математических моделей и численных исследований волновых процессов и динамики структуры жидких сред с микронеоднородностями при импульсном нагружении. Показано, что присутствие микронеоднородностей может существенно изменить состояние среды, в которой начинает развиваться кавитация, структуру прикладываемого волнового поля, а также механизмы пртекающих в жидких средах физических процессов.
1. Введение
Известно, что реальные жидкости содержат микронеоднородности типа микропузырьков свободного газа, твердых частиц или их комбинаций (рис.1), которые играют роль ядер кавитации [1]. При приложении растягивающих напряжений на ядрах начинают развиваться парогазовые пузырьки: в жидкости формируется пузырьковый кластер, форма и динамика плотности которого определяются спектром начальных размеров ядер. Кавитация существенно изменяет не только состояние среды в поле растягивающих напряжений, но и параметры и структуру самого прикладываемого поля. Согласно современным представлениям, основанным на многочисленных экспериментальных исследованиях, плотность микропузырьков свободного газа в чистой жидкости имеет порядок 103-104 см-3 при Ко ~ 1.5 мкм. Теоретический спектр плотностей ядер кавитации, представленный на рис. 1, где сравниваются экспериментальные данные ряда авторов, имеет максимум распределения Мо ~ 105-106 см-3, который определен экспериментально по трекам дифракционных пятен и соответствует общей плотности ядер любой природы.
*© В.К.Кедринский, В. А. Вшивков, Г. И. Дудникова, Ю.И. Шокин, 1997.
^Работа выполнена при поддержке Российского фонда фундаментальных исследований, грант №96-0219369.
Рис. 1. Стpуктуpа микpонеодноpодностей; 1 — пузьфек свободного газа, 2 стpуктуpы. Спектp плотности ядеp кавитации любой пpиpоды.
комбинационные
Состояние реальных жидкостей с микронеоднородностями носит явно двухфазный характер, что позволило использовать соответствующую математическую модель для описания волновых процессов в жидкостях, кавитирующих за фронтом мощных волн разрежения [2, 3]:
1р + р V = 0,
оЬ
¿V V»
17 + =0, оЬ р
р = (1 - к)р1, к = кс
Я
Яс
_ 1 ЯЯ 3 , ¿^Я ч о 1 г 1
ЯИё + 2(Л)2 = р- 1»'-»]•
(1)
Как показывает опыт, неоднородность состояния среды и физика кавитационных явлений имеют прямое отношение к таким проблемам, как инициирование детонации в жидких ВВ, прочность жидкости при динамическом нагружении, и к хорошо известным, но не полностью изученным эффектам крупномасштабных взрывов контейнеров с горючим в результате аварий при транспортировке. Было показано [4, 5], что приведенная выше двухфазная система, дополненная кинетической моделью адиабатического взрыва химически активной газовой смеси в пузырьках [6]:
= Л^Те~Еа/вт(а - п)2, о'-Г =
(2)
способна описать эти эффекты. Здесь Еа — энергия активации, В — газовая постоянная, д — теплота, выделяемая в каждом элементарном акте реакции, — теплоемкость в единице объема, а — начальная концентрация исходной компоненты, п — число молекул, формирующихся в единице объема в процессе реакции.
В статье рассмотрены некоторые подходы к решению указанных проблем, для которых кавитационные эффекты могут иметь принципиальное значение. При этом нас будут интересовать два разных аспекта поведения жидких сред с микронеоднородностями при импульсном нагружении: 1 — гидродинамический (механизм кавитационного разрушения жидкости в интенсивных волнах разрежения) и 2 — волновой (возможность возбуждения
3
режима типа пузырьковой детонации при взаимодействиии ударных волн и волн разрежения). Сюда можно отнести, в частности, кавитационные процессы в легкоиспаряющихся и воспламеняющихся жидкостях, содержащихся под давлением в контейнерах, когда последние внезапно разгерметизируются, например, в результате удара. Этот процесс имеет прямое отношение к так называемым паровым взрывам, с которыми связывают катастрофические разрушения железнодорожных контейнеров и накопительных емкостей, содержащих сжиженные под давлением газы и воспламеняющиеся жидкости.
2.Структура потока при импульсном разрушении жидкости
Гидродинамическая стадия развития взрывного процесса рассматривалась в ряде работ, среди которых можно отметить результаты [7-9]. Чаще всего приводятся два механизма формирования условий для последующего развития взрыва [9]:
а) возможное извержение сжатых жидкостей в виде двухфазных кавитирующих струй и их последующее распыление, в результате которого формируется газо-капельное облако (возникает так называемый объемно-детонирующий заряд);
б) при частичной разгерметизации волны разрежения распространяются внутрь контейнера с жидким горючим и вызывают быстрые процессы типа кипения (давление внутри сосуда резко увеличивается).
Процесс образования капельного облака в результате кавитационного разрушения жидкого объема при внезапном снятии давления моделируется в экспериментах с гидродинамическими ударными трубами путем неограниченного инерционного развития зоны кавитации, возникающей при отражении ударной волны от свободной поверхности жидкости. Оказалось, что здесь имеет место целый комплекс явлений, который включает стадии зарождения и формирования кавитационных кластеров, их развитие в пеноподобную структуру и, наконец, переход в газо-капельную фазу [1].
В результате исследований первых двух стадий был сформулирован механизм развития пузырьковых кластеров, согласно которому в интенсивной фазе разрежения жидкость становится монодисперсной средой, и обнаружено, что крутизна фронта волны разрежения т играет основную роль в определении предельных растягивающих напряжений, допускаемых кавитирующей жидкостью. Например, для объемной концентрации пузырьков ко = 10-10, их радиусов Ко = 10-4см и давления р = -30 МПа в падающей волне разрежения с вертикальным фронтом время релаксации составит всего ¿сау ~ 8.3 • 10-8 с, а максимальные растягивающие напряжения в кавитирующей жидкости не превысят р* ~ -3.8 МПа, если крутизна фронта волны разрежения равна т = 1 мкс.
Показано, что волновые процессы и динамика структуры кавитирующей жидкости характеризуются рядом критериев, среди которых отметим время релаксации кавитиру-ющей среды к состоянию "насыпной плотности"пузырьков (к ~ 0.5 — 0.7 )
~ (рг/р — 1)ё-1,
где ё — скорость деформации; рг, р — плотность жидкости и средняя плотность кавитирующей среды. Очевидно, для состояния "насыпной плотности пузырьков"^ ~ ё-1, что дает соответственно 2 мс для скорости деформации ё = 500 с-1 и 0.7 мс для ё = 1330 с-1. Эти оценки вполне коррелируют с экспериментальными данными, полученными при измерении динамики средней плотности кавитационной зоны в результате компьютерной
Рис. 2. Схема экспеpимента с каплей.
%
Рис. 3. Ячеистая стpуктуpа pазpушающейся капли пpи ми^осекундном удаpе.
обработки рентгеновских негативов зоны после оцифровки на лазерном сканере РЕЛС [10] и соответствуют принятой модели инверсии структуры течения с промежуточным состоянием типа пены [1].
Однако исследование физического механизма перехода из пенной фазы в капельную потребовало разработки специальной методики, согласно которой в качестве исследуемого объекта использовалась капля жидкости, а в качестве источника возбуждения кавитации — ультракороткие ударные волны. Экспериментальная установка (рис. 2) представляет собой силовой элемент электромагнитной ударной трубки, на диафрагму которой помещается капля дистиллированной воды радиусом в несколько миллиметров. При разрядке батареи конденсаторов С на плоскую спиральную катушку Я, помещенную между проводящими мембраной и диском Б, в зазоре между мембраной и катушкой возникает импульсное магнитное поле, толкающее мембрану. При ударе мембраны по капле в последней формируется ударная волна с амплитудой 5-6 МПа и длительностью порядка 3-4 мкс. Параметры волны оценивались по динамике свободной поверхности тонкого слоя жидкости на мембране, записанной емкостным датчиком.
Эксперименты показали, что процесс разделяется на три основные стадии [11]. Вначале формируется плотная кавитационная зона ячеистой структуры, основу которой составляют крупные кластеры (порядка 1 мм), объединяющие кавитационные зародыши, растущие в поле растягивающих напряжений. При этом размер капли практически не меняется. Согласно теоретическим исследованиям структуры и параметров волны разрежения, выполненным в рамках модели (1), растягивающие напряжения в кавитирующей жидкости релаксируют очень быстро и вторая стадия процесса разрушения характеризуется инерционным развитием зоны кавитации: капля переходит в своеобразную пенную фазу с четко выделяющимися границами в виде сетки жидких жгутов (рис. 3). Элементы жидкой сетки при дальнейшем растяжении постепенно разделяются на отдельные фрагменты,
фрагменты — на струйки (третья стадия), которые в силу неустойчивости разрушаются на отдельные капли (см. правый кадр рис. 3 — увеличенный фрагмент полной структуры, момент времени около 1.5 мс после удара). Характерно, что после разрушения даже в капельном состоянии поток сохраняет ячеистый рисунок структуры в течение длительного времени.
Естественно, знания механизма инверсии двухфазного состояния как перехода кавитирующей жидкости в газо-капельную систему недостаточно для того, чтобы ответить на вопрос о связи физических параметров возникшей системы с начальными условиями: параметрами ударной волны и состоянием микронеоднородностей жидкой капли. Это задача другого исследования.
ЗКавитационные эффекты в химически активных пузырьковых средах
Кроме названных выше двух механизмов генерации высоких давлений возможен и третий. При заполнении контейнера в жидкости может образоваться большое число пузырьков, содержащих смесь воздуха с парами горючего, способную воспламеняться при адиабатическом сжатии. Если контейнер, движущийся с большой скоростью, сталкивается с препятствием, то в результате удара в жидкости возникает система волн сжатия, которые могут возбудить волну пузырьковой детонации. Наоборот, при разгерметизации контейнера в жидкости возникает система волн разрежения, результат взаимодействия которых, как будет показано ниже, может привести к неожиданным эффектам.
Явления типа пузырьковой детонации исследовались давно, в частности, при попытках понять механизм инициирования жидких взрывчатых веществ (ВВ). Так, F. Bowden and A. Yoffe [12] указали на тепловую природу инициирования взрыва жидких ВВ ударом, подчеркнув, что небольшие газовые пузырьки, которые обычно присутствуют в жидкостях, являются возможными тепловыми источниками взрыва (модель Ьо1-8ро1э). K. Andreev [13] и C. Johansson [14] предположили, что присутствие небольших частичек или газовой фазы ВВ в пузырьках — наиболее реальная причина инициирования детонации.
А. Дубовик и В. Боболев [15] исследовали ту же идею и возможность влияния капелек или паров ВВ в газовой атмосфере пузырька на развитие механизма hot-spots, оценив величину адиабатического периода индукции и критический размер, при котором капелька успеет быстро прогреться за время схлопывания пузырька. Рассматривался также механизм дробного удара, когда детонационный процесс инициируется последовательностью ударных волн: первая волна отражается от свободной поверхности, возникает кавитация, затем вторая волна сжимает пузырьки, образуя систему горячих точек, и возбуждает детонацию.
A. Campbell et al. [16] экспериментально доказали, что механизм ударного инициирования определяется взаимодействием волны с микронеоднородностями, приводящим к формированию горячих точек и к последующему усилению волны по мере распространения до параметров, когда выделяемой за фронтом энергии окажется достаточно, чтобы возникла детонация. Таким образом, можно констатировать, что история развития исследований чувствительности, например, жидких ВВ к удару показывает, что эта важнейшая характеристика безопасности зависит от количественных параметров и природы микро-неоднородностей в жидкостях, способных к взрывным процессам.
Рассмотрим некоторые сценарии развития волновых процессов и возможности усиле-
ния волн в пузырьковых средах, которые могли бы объяснить механизмы упомянутых выше аварийных ситуаций. Волновые процессы рассчитывались в рамках моделей (1)-(2) по программам, разработанным в ИВТ СО РАН. Основная система уравнений для средних давления р, плотности р и массовой скорости и в лагранжевой системе координат для цилиндрической симметрии, позволяющей исследовать фокусировку волн, имеет вид :
ди д
1 / х(г, г) Ро \ г
др дх
д? Ж _и' р _ РО
1 1 (х(г, ¿)
дх дг'
р
2
1 + р°£° )п - 1],
прО 1 — к'
1 — корв , дг
Б.
Кинетические уравнения :
„дБ 3 2 Т С2 „ Б
% +2 Б 2 = С — У — С в — Р
дТ дЖ 1ЧАТ в3(1 — Т) 0. 1ЧТБ
П— + 8(7 — 1)^и 1 Т ; — 3(7 — 1) —,
дг
дг
(3)
Ж _ а (1 — Ж )2,
Ре = С4(7 — 1)
С _ рдоТоВ с2
в |Б | |1 — Т |
Ми
>/Р£, Ре > 100, 10, Ре < 100.
2а
Сз _
С4 _
12Ко^/ро/ро
V
РоМ КоРо Вол/РоРо
Здесь Ж — доля прореагировавших компонент газовой смеси, Т — температура в газовой смеси, в _ К/Ко — относительный радиус пузырьков, Б — их радиальная скорость, а — энергия активации, £ — константа, зависящая от начальных параметров среды и состава смеси, п — теплота реакции, В — универсальная газовая постоянная, М — молярная масса. При 8 _ 0 учитывается теплообмен.
Расчет проводился для водородо-кислородной смеси (а _ 51.5, (/Ко _ 2.7• 109, п _ 70). Из-за отсутствия экспериментальных данных по микронеоднородностям в жидких ВВ и горючих жидкостях в расчетах приходилось ориентироваться на данные для простых жидкостей типа свежей или дистиллированной воды. Безразмерные переменные были определены как г _ Га/Ко, х _ хЛ/Ко, и _ ил/ л/рО/Роо, где рО, рОО — начальные давление и плотность жидкости, г _ глсо/Ко, р _ 'Ра/'Ро, Т _ ТЛ/То, р _ р^/роо, а индексом й отмечается размерная переменная.
Для расчета использована неявная разностная схема [17]; уравнения газодинамики сведены к одному нелинейному уравнению второго порядка по давлению Р, которое решалось методом итераций Ньютона. Общая разностная схема имеет вид
Qг-l/2^Р{t}l —
2Н2
Qi-l/2 + ^+1/2 +--2" %%
1 (др
т2 "грт+1(к) V дР
т+1(к)'
Ар(к) +
+Qi+1/2Ар(+)1 + —2"^^-т+ЙТкг — [хт+1/2 — х 1-1/2 + т (и™1/2 — и™1/2)] +
Т рг Т
2
2
г
к
+дт/2(С+1(к) - рГ+1(к)) - ^-1/2(рГ1(к) - Р—+1(к)) = 0,
здесь к — число итераций, а
(Х-1/2)2
п _ ^ г—1/2) „ _ РоГ.
Чг-1/2 — -2-'
Рог2—1/^ * (жГ)2' Давление и плотность рассчитывались по схеме
Рт+1(к+1) — рГ+1(к) + Рт+1(к+1) — р(рГ+1('+1))-
Условие выполнения итераций
тах | Ар(к) | < е.
г
Кинетические уравнения решаются в каждой точке интервала по времени методом Рунге— Кутта—Мерсона с переменным шагом по времени
гт <г< ¿т+1 — гт + т.
Формирование волны пузырьковой детонации (плоский одномерный случай) показано на рис. 4. Здесь N — доля прореагировавших компонент газовой смеси, Р — среднее давление в волне пузырьковой детонации, Т — температура в продуктах детонации газовой смеси, К/Ко — относительный радиус пузырьков, Б — их радиальная скорость.
Известно, что по мере распространения волны по пузырьковой среде происходит ее затухание в силу потерь на сообщение жидкому компоненту кинетической энергии и увеличение внутренней энергии газовой фазы в пузырьках при их сжатии. Если пузырьки наполнены взрывчатой газовой смесью, то при ее адиабатическом нагреве возможно возникновение реакции с большим выделением энергии. В окружающую жидкость будет излучаться волна сжатия, которая компенсирует потерю энергии падающей волны, в результате чего в среде может установиться самоподдерживающийся режим, получивший название пузырьковой детонации [18]. Экспериментально этот эффект подробно исследован А. Сычевым и А. Пинаевым [19]. Расчет показывает, что предложенная выше модель (3) вполне адекватно описывает процесс формирования и распространения волны пузырьковой детонации. Ее амплитуда при рассматриваемой объемной концентрации газовой фазы (ко — 5 • 10—3) неустойчива и устанавливается в интервале 80-100 МПа. Столкновение таких волн приводит к дополнительному их усилению (амплитуда "подска-кивает"примерно до 300 МПа). На рис. 5 показан процесс формирования и столкновения двух волн пузырьковой детонации.
При анализе возможных сценариев развития волновых процессов возникает вопрос о пороговых значениях интенсивности волн, способных к возбуждению детонационного процесса в пузырьковых активных средах. Частично ответ на него дан на рис. 6, где демонстрируется эффект столкновения слабых ударных волн (амплитуды 0.4 МПа). Сами по себе такие волны не способны возбудить детонацию, однако в результате их столкновения необходимые условия могут возникнуть. Как видно из рисунка, процесс с химической реакцией в газовой фазе, возбужденный в плоскости столкновения, по мере распространения к периферии сопровождается формированием детонационной волны, амплитуда которой уже на расстоянии 5 см от центра достигает 60 МПа.
Рис. 4. Фоpмиpование волны пузырьковой детонации; к0 = 5 • 10—3, Яо = 0.2 см, р = 5 МПа
Рис. 5. Усиление ударных волн при столкновении в химически активной пузырьковой системе; к0 = 5 • 10—2, Я0 = 0.2 см, Дг = 1 МПа.
Рис. 6. Инициирование пузырьковой детонации при столкновении слабых ударных волн; = 0.4 МПа, к0 = 0.02, Я0 = 0.1 см.
2 см
О 1 2 см
Рис. 7. Инициирование пузырьковой детонации при столкновении волн разрежения; = —10 • ехр(—1/т) МПа, т = 2 мкс, ко = 0.001, Яо = 0.01 см.
100 з
Рис. 8. Динамика амплитуды ударной волны при фокусировке; к0 = 210 2, К0 = 0.01 см, Рь = 1 МПа, г = 5 см.
Неожиданный эффект был открыт при численном анализе столкновения волн разрежения (рис. 7) [5]. Эта постановка моделирует естественный процесс импульсного нагруже-ния образца жидкого ВВ со свободной границей (в плоской постановке таких границ две): возбужденная при ударе волна сжатия взаимодействует со свободной поверхностью, от которой по образцу должна распространяться волна разрежения. Расчет демонстрирует возможность инициирования волны пузырьковой детонации в центре образца в результате столкновения волн разрежения. Эффект основан на особенностях трансформации волны разрежения по мере ее распространении в среде с микронеоднородностями и возбуждения ею кавитации: за фронтом волны возникает пульсация с интенсивной положительной фазой, которая при столкновении волн и отвечает, в конечном счете, за разогрев газовой смеси до температур воспламенения. Заметим, что интенсивности положительной фазы одиночной волны разрежения для этого, как правило, оказывается недостаточно. Столкновение волн в центре образца шириной 2 см произошло около 8.5 мкс от начала распространения, и уже к 17.5 мкс к свободной поверхности подошла детонационная волна с амплитудой около 40 МПа (см. рис. 5).
В реальной обстановке не исключена и возможность фокусировки ударных волн в пузырьковых пассивных и активных средах, которая также может привести к значительному повышению давления в среде и стать причиной разрушения контейнеров. Процесс фокусировки волн в пузырьковых средах осложняется тонкой структурой формирующегося в них волнового поля. Естественно, что в первую очередь фокусируется предвестник. Отражаясь в фокусе, он будет взаимодействовать с основной падающей волной и искажать ее
Рис. 9. Фокусировка волны пузырьковой детонации; ко = 10-2, Яа = 0.2 см, Рь = 3 МПа.
профиль.
На рис. 8 показана динамика усиления амплитуды падающей волны в результате ее фокусировки на ось в пассивной среде. Расчет показал, что амплитуда на оси достигает 350 МПа.
Еще больший эффект наблюдается, если в процессе фокусировки в среде за падающей волной возбуждается пузырьковая детонация. На рис. 9 показано распределение давления вдоль г для трех характерных моментов времени. При ¿=32 мкс падающая волна достигает оси и ее амплитуда возрастает в 5 раз. В этот момент в окрестности границы (см. на рис. 9 справа) возникает волна пузырьковой детонации, амплитуда которой к ¿=60 мкс (на расстоянии 2 см от оси) достигает уже 40 МПа. Расчет показывает, что в окрестности оси амплитуда достигает весьма значительных величин — порядка 103 МПа. Нет сомнений, что возникновение таких напряжений может привести к катастрофическим последствиям.
4. Заключение
Необходимо отметить, что возникающие при фокусировке высокие давления являются следствием инерционных эффектов, сопровождающих процесс пульсации пузырьков, и пересжатия продуктов детонации до высоких степеней. Температура последних достигает (1.5—2) 104 К, что требует введения в систему (3) дополнительных условий на изменение состояния продуктов детонации.
Другая проблема связана с недостаточностью учитываемого теплообмена для потери температуры продуктами детонации в процессе первой пульсации, что соответствовало бы экспериментальным данным о формировании действительно уединенной волны: в расчетах в структуре волны остается высокотемпературный "хвост", который как мостик связывает детонационную волну с основным возмущением, бегущим по прореагировавшей системе.
В то же время приведенные в работе экспериментальные данные, постановки и результаты расчета показали, что многие физические явления, возникающие при импульсном нагружении жидкости, а также механизмы их развития могут рассматриваться с точки зрения волновых процессов и динамики структуры кавитирующих сред.
Исследуя процессы импульсного диспергирования жидкости необходимо учитывать, что на последней стадии кавитационного разрушения в интенсивных волнах разрежения жидкий объем трансформируется в ячеистую структуру типа жидкой сетки, разрушение которой на отдельные струи и затем на капли составляет основу процесса инверсии двухфазного состояния жидкости и формирования газо-капельного облака.
В проблеме инициирования детонации жидких ВВ и крупномасштабных взрывов в аварийных ситуациях на транспорте особую роль играют механизмы, связанные с влиянием микронеоднородностей типа пузырьков инертного или химически активного газов. Применение модели пузырьковой детонации в качестве механизма hot-spot для инициирования детонации жидких ВВ и крупномасштабных взрывных процессов как результата взаимодействия ударных волн и волн разрежения, а также их фокусировки в активных пузырьковых средах расширяет спектр возможных механизмов, управляющих развитием рассмотренных явлений.
Авторы выражают признательность Ю. И. Молородову за помощь при проведении численных расчетов.
Список литературы
[1] Кедринский В. К. Нелинейные проблемы кавитационного разрушения жидкости при взрывном нагружении (обзор). ПМТФ, №3, 1993, 74-91.
[2] КЕДРИНСКИЙ В. К. Динамика кавитационной зоны при подводном взрыве вблизи свободной поверхности. ПМТФ, №5, 1975, 68-78.
[3] KEDRINSKII V. K. Negative pressure profile in cavitation zone at underwater explosion near free surface. Acta Astronáutica, 3, No. 7-8, 1976, 623-632.
[4] KEDRINSKII V. K. The role of cavitation effects in the mechanisms of destruction and exp1osive processes. Int. J. Shock Waves, 7, No. 2, 1997, 63-76.
[5] Кедринский В. К., Вшивков В. А., Дудниковл Г. И., Шокин Ю. И. Взаимодействие волн в химически активных пузырьковых средах. Докл. РАН, 349, №2, 1996, 185-188.
[6] TODES O. M. Adiabatic term explosion. J. Phys. Chem., 4, No. 1, 1933, 71.
[7] HILL L., STERTEVANT B. An experimental study of evaporation waves in a superheated liquid. In "Proc. of IUTAM Symp. on Adiabatic Waves in Liquid-Vapour Systems", Springer, Gottingen, 1989, 25-37.
[8] Chaves H., Lang H., Meier G., Speckmann H. Lecture Notes in Physics: Flow in Real Fluids. Springer, Berlin, 1985.
[9] Barbone R., Frost D., Makris A., Nerenberg J. Explosive boiling of a depressurized volatile liquid. In "Proc. IUTAM Symp on Waves in liquid/gas and liquid/vapour two-phase systems, Kluwer Academic Publishers, Kyoto, 315-324.
[10] БЕРНГАРДТ А. Р., КЕДРИНСКИй В. К., Плльчиков Е. И. Эволюция внутренней структуры зоны разрушения жидкости при импульсном нагружении. ПМТФ, №2, 1995, 99-105.
[11] КЕДРИНСКИй В. К., Бесов А. С., Гутник И. Э. Инверсия двухфазного состояния жидкости при импульсном нагружении. Докл. РАН, 352, №4, 1997, 477-479.
[12] Bowden F. P., Yoffe A. D. Fast Reaction in Solid. London Butterworths Scientific Publications, 1958.
[13] Andreev K. K. Some consideration on the mechanism of initiation of detonation in explosive. In "Proc. Royal Soc.", London, 1958, 195 A 246(1245), 257-267.
[14] Johansson C. H. The initiation of liquid explosives by shock and the importance of liquid break-up. In "Proc. Royal Soc.", London, 1958, A 246(1245), 160-166.
[15] ДУБОВИК А. В., БОБОЛЕВ В. К. Чувствительность жидких систем к удару. Наука, М., 1978.
[16] Campbell A.W., Davis W.C., Travis J.R. Phys. of Fluids, 4, No. 4, 1961, 498.
[17] САМАРСКИЙ A., Попов Ю. Разностные схемы в газовой динамике. Наука, M., 1975.
[18] Hasegawa T, Fujiwara T. Detonation in oxyhydrogen bubbled liquids. In "Proc. of 19th Int. Symp. on Combustion, Hafia, 1982, 675-683.
[19] ПИНАЕВ А., СыЧЕВ А. Влияние физико-химических свойств газа и жидкости на параметры и условия возникновения детонационных волн в системах "жидкость — газовые пузырьки". ФГВ, 23, №6, 1987, 76-84.
Поступила в редакцию 15 сентября 1995 г.
