Осесимметричные волны в пузырьковой жидкости

Баязитова, Алина Разифовна

Осесимметричные волны в пузырьковой жидкости тема автореферата и диссертации по механике, 01.02.05 ВАК РФ

Баязитова, Алина Разифовна АВТОР

кандидата физико-математических наук УЧЕНАЯ СТЕПЕНЬ

Тюмень МЕСТО ЗАЩИТЫ

2006 ГОД ЗАЩИТЫ

01.02.05 КОД ВАК РФ

Диссертация по механике на тему «Осесимметричные волны в пузырьковой жидкости»

Автореферат диссертации на тему "Осесимметричные волны в пузырьковой жидкости"

На правах рукописи

БАЯЗИТОВА АЛИНА РАЗИФОВНА

ОСЕСИММЕТРИЧНЫЕ ВОЛНЫ В ПУЗЫРЬКОВОЙ жидкости

01.02.05 - Механика жидкости, газа и плазмы

АВТОРЕФЕРАТ диссертации на соискание ученой степени кандидата физико-математических наук

Тюмень 2006

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

Научный руководитель: член-корреспондент АН РБ, доктор

физико-математических наук, профессор Шагапов Владислав Шайхулагзамович

Научный консультант: кандидат физико-математических наук, доцент

Официальные оппоненты: доктор физико-математических наук, профессор

Защита диссертации состоится « 17 » мая 2006 г. в 13 час, на заседании диссертационного совета ДМ 212.274.09 в Тюменском государственном университете по адресу: 625003, г. Тюмень, ул. Перекопская, д. 15А

С диссертацией можно ознакомиться в библиотеке Тюменского государственного университета по адресу: 625003, г. Тюмень, ул. Семакова, д. 10.

Автореферат разослан « Ж». апреля 2006 г.

Гималтдинов Ильяс Кадирович

Лежнин Сергей Иванович,

кандидат физико-математических наук Санников Иван Николаевич

Ведущая организация: Башкирский государственный университет

Ученый секретарь диссертационного совета, кандидат физико-математических наук

ШМ-

ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ

Актуальность темы. Значительный интерес исследователей к проблемам и задачам механики пузырьковых сред обусловлен широким распространением таких систем в природе и их интенсивным использованием в современной технике. Газожидкостные среды широко используются в различных отраслях промышленности (химической, нефтеперерабатывающей, металлургической и др.) для ведения технологических процессов (например, жидкофазное окисление углеводородов) или как элемент системы обеспечения взрывобезопасности (жидкостные предохранительные затворы). Кроме того, пузырьковые жидкости могут найти применение в литотрипсии (разрушении почечных камней), в качестве активного элемента в сазере (акустическом аналоге лазера). К настоящему времени одномерные волны экспериментально и теоретически подробно изучены, и активно ведутся исследования по динамике двумерных волн (A.A. Губайдуллин, В.Е. Накоряков, В.Е.Донцов, В.К. Кедринский, Г.Г.Лазарева, В.Г. Гасенко, С.А. Ждан, В.Ю. Ляпидевский, Н.К. Вахитова, Y. Matsumoto и др.).

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

Цели работы. Теоретическое исследование динамики волн в пузырьковой жидкости в двумерной осесимметричной постановке задачи. Изучение эволюции волн давления в пузырьковой жидкости при неоднородном по сечению распределении пузырьков. Исследование распространения волн давления в трубе, содержащей сферические пузырьковые кластеры. Изучение динамики детонационных волн в трубчатом пузырьковом кластере. Анализ процесса выхода детонационных волн из трубчатого газожидкостного кластера в однородную пузырьковую область.

Научная новизна. В диссертации поставлен и решен ряд новых важных задач. Изучено влияние неоднородности распределения пузырьков по сечению трубы на динамику нелинейных волн. Объяснены механизмы усиления волн газожидкостным кластером. Рассмотрены особенности структуры ударной волны в газожидкостной среде кластерной структуры, в зависимости от числа, размеров и взаимного расположения кластеров. Исследована динами-

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

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

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

Апробация работы. Основные результаты работы были представлены и обсуждались на следующих конференциях и научных школах:

- на Международной научной конференции «Спектральная теория дифференциальных операторов и родственные проблемы» (Стерлитамак, 2003);

- на девятой Всероссийской научной конференции студентов-физиков и молодых ученых (Красноярск, 2003);

- на третьей всероссийской молодежной научной школе-конференции (Казань, 2003);

- на Всероссийской научно-теоретической конференции «ЭВТ в обучении и моделировании» (Бирск, 2004);

- на XIII сессии Российского акустического общества (Москва, 2004);

- на Всероссийской научной конференции «Современные проблемы физики и математики», посвященной 50-летию физико-математического факультета (Стерлитамак, 2004);

- на десятой Всероссийской научной конференции студентов-физиков и молодых ученых (Москва, 2004);

- на XV сессии Российского акустического общества (Нижний Новгород, 2004);

-на V Региональной школе-конференции для студентов, аспирантов и молодых ученых по математике и физике (Уфа, 2005),

- на XVI сессии Российского акустического общества (Москва, 2005) -на IV Региональной научно-методической конференции «ЭВТ в обучении и моделировании» (Бирск, 2005).

Кроме того, результаты диссертационной работы докладывались и обсуждались на научных семинарах кафедры прикладной математики и механики СГПА под руководством профессора В.Ш. Шагапова, кафедры математического анализа СГПА под руководством профессоров К.Б. Сабитова и И.А. Калиева, кафедры теоретической физики СГПА под руководством профессора А.И. Филиппова, на научном семинаре Института математики с вычислительным центром УНЦ РАН под руководством профессоров М.Д. Рамазанова, Н.Д. Морозкина и получили положительную оценку.

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

Структура и объем работы. Диссертация состоит из введения, четырех глав, заключения и списка литературы. Работа изложена на 134 страницах, иллюстрирована 64 рисунками. Список литературы состоит из 111 наименований.

КРАТКОЕ СОДЕРЖАНИЕ РАБОТЫ

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

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

В § 1.1 представлен обзор теоретических и экспериментальных работ по распространению нелинейных волн давления в пузырьковой жидкости.

В § 1.2 приведен обзор теоретических и экспериментальных работ по динамике детонационных волн в пузырьковой жидкости.

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

вует массообмен между пузырьками и жидкостью, температура жидкости постоянна, отсутствуют дробление и слипание пузырьков) имеет вид [1]

М+У(/уЗ) = 0, 0 = /,*), ~- + У(пи) = 0, с1а

3 2 4 м'я Рх~Р1 Ре ~ Р1

dt 2 а р? " piCta

wA =-

, о 4 з

p = p/+pg, at +ag = 1, pt - p, a,, ag=~mia ,

где p,°, at, и (uz, v,), ph n, a, w,y,q~ соответственно плотность, объемное содержание г'-й фазы, скорость (составляющие скорости по координатам гиг) и давление несущей жидкости, число и радиус пузырьков, скорость радиального движения, показатель адиабаты для газа, интенсивность теплообмена или тепловой поток от газа к жидкости, отнесенный к единице площади межфазной поверхности.

Интенсивность межфазного теплообмена принята в виде [1]

где Г0 = const - температура жидкости, Nu - число Нуссельта. Числа Нуссельта и Пекле заданы в виде:

Ре >100 ,ч То a\w\

Nu = ] Ре = 12(Г-1)т

Ре <100 \т8-Цкя

где = 1&1с%р&о, с6 - коэффициенты температуропроводности, теплопроводности и теплоемкость газа соответственно.

Жидкость полагалась акустически сжимаемой, а газ калорически совершенным:

Р, = Р* + С} (А° - Ао). Ря = Р°вВТ8 ■ Во второй главе представлен переход от эйлеровых переменных к ла-гранжевым, поскольку в лагранжевых переменных первоначальные границы неоднородностей остаются неподвижными. Здесь же рассмотрено построение разностной схемы для численных расчетов. Проведено тестирование модели по данным [1].

В § 2.1 приводится переход от эйлеровой системы к лагранжевой системе координат. Математическая модель в лагранжевых переменных имеет вид

д£1_ а

ЧР1

1 ОС 0

ди2 д1

дог

~дГ дрё

а/ д1

Рщ . а8 2+ У

Рр

др[ дг дг0 дг0

да

ха/ з &'

ф/ дг дг0 дг0

г 1 (др1 дг др1 дг г0 ,1р0 [ дг0 дг0 5г0

& Ы

дг

Зи>* . 30"-1)

а?

-И'

«0

5а

аг

= = +Т4>

а/

д =

"Л 3 2 О--

4у

Н'в

'А >

"Л

рГС^

2 а

[ш = л/рё, Ре >100, N11 = 10, Ре <100,

Ем.

Ро

Ре = 12(х -1)

'О

к0 =-

_ г ( дг дг дг дг ~ г0{дг0 дг0 дг0дг0

дJ _иг + г (до2 дг ^ дг диг ди2 дг дг диг дг г г0 ^ дг0 дг0

дг$ дг0 дгп дг,

дг0 дг0

(2)

В § 2.2 описана разностная аппроксимация системы уравнений (2) для двух случаев: на границе г0= 0 задается давление или скорость среды. Система дифференциальных уравнений аппроксимирована на пространственно-временной сетке со = со/,*о)т, где

Щ= {(г„ Гу)}, г1+1 = г,+кг0, г,н/2 = г,+0.5-Л^; г = 0,1,...,М1 - 1,

0+1 = г)+Ьго, 0+1/2 = = 1,...,М2- 1,

а>т = {(4, 4+1 = 4+т, к = 0; 4 = Ат, А = 0,1,2,...)} К узлам сетки (г„ гу) отнесены сеточные функции скорости г>г!Д иГ|/ и эйлеровых переменных г,/ и г,/, к "полуцелым" точкам (Х+]/2, />+1/2) - сеточные функции всех остальных параметров.

В § 2.3 описана реализация начальных и граничных условий.

В третьей главе рассмотрена динамика нелинейных волн в жидкости при неравномерном по сечению распределении пузырьков. Для каждого из случаев распределения пузырьков (однородная пузырьковая жидкость, пузырьковое кольцо, пузырьковое ядро) произведено сравнение с экспериментами Накорякова В.Е., Донцова В.Е. [2]. Показано, что эволюция волн давления типа "ступенька" в трубе, содержащей трубчатый пузырьковый кластер, сопровождается появлением волн давления, амплитуда которых многократно превышает амплитуду первоначальной волны. Рассмотрена эволюция волн давления в трубе, содержащей различное количество газожидкостных кла-

Др = Р/-Ро>МПа

11' Г ПГ I2 Г "К ч 1

12010123 231 3 4 6 6 7 4 6 /,мс

Рис. 1 Динамика импульсного давления в трубе с кольцевым расположением пузырьков по сечению

а - схема задачи;

б - затухание волн давления по длине рабочего участка: символы 1 и 2 экспериментальные данные [2] в случае гомогенной структуры (Лр,/р0 = 2.5) и кольцевой структуры среды (Лро/рв =2.1); сплошная и штриховая линии соответствуют инициированию первоначального сигнала на границе г0 = 0 граничным давлением или жестким ударником в случае кольцевой структуры (Лр(/р0 = 2.1), пунктирная линия - расчет для гомогенной структуры среды (Ар(/р0 ~ 2.5); в- д осциллограммы давления в датчиках 01-03, расположенных соответственно на расстояниях г0=0 м 0.25 м и 0.76 м, показывающие динамику импульса при кольцевом расположении пузырьков; линии 1 - эксперимент, линии 2 и 3 соответствуют заданию на границе г0 = 0 скорости ударника или давления.

стеров. Для случая пузырьковой жидкости кластерной структуры было проведено сравнение с экспериментами Донцова В.Е. [3] и с профилем, рассчитанным по уравнению Буссинеска для газожидкостной среды.

В § 3.1 исследована эволюция волн давления в трубе при неоднородном по сечению распределении пузырьков газа.

На рис.1 представлены экспериментальные и расчетные данные по динамике волн давления в случае кольцевого расположения фреоновых пузырьков в трубе (схематическое изображение приведено на фрагменте (а) ). Параметры системы: жидкость - 50 %-ный по массе раствор глицерина с водой, газ - фреон-12, = 0.027 м, Я0 = 0.019 м, а0 = 0.53 мм, а^ = 0.01. Импульс инициировался воздействием давлением или жестким ударником, т.е. на границе г0 = 0 задавалась скорость в виде треугольника с амплитудой возмущения скорости До2 = 0.7 м/с, периодами подъема и спада скорости ударника - = 1.810"4 с и г2 = 2.2-10"4 с соответственно. Из сравнения данных, приведенных на фрагменте (б) следует, что, когда пузырьки расположены в кольцевом слое вблизи стенки, и в случае гомогенного распределения имеет место неплохое соответствие между экспериментальными и расчетными данными. Из фрагмента (б) видно, что на начальных расстояниях (г0 <0,25 м) амплитуда сформировавшегося всплеска давления во всех случаях превышает амплитуду первоначального сигнала. При этом для кольцевого режима всплеск давления в случае, когда первоначальный импульс образуется под действием граничного давления, превышает по амплитуде аналогичный всплеск давления, реализуемый, когда первоначальный импульс задается скоростью жесткого ударника. Это происходит потому, что за все время воздействия граничным давлением амплитуда и профиль давления по сечению 0<г0<Лс однородны, а в случае задания на границе г0 = 0 скорости движения среды, хотя датчик, расположенный на стенке в области пузырькового кольца фиксирует такой же по амплитуде и форме сигнал, что и в случае задания граничного давления, в целом по сечению из-за разгрузки на границе г0 = реализующееся давление неоднородно по амплитуде и по форме. Далее полученный импульс распространяется в расслоенной среде, при этом возмущения, распространяющиеся в виде предвестника вдоль осевой зоны по чистой жидкости, значительно опережают основную волну. Вследствие того, что акустический импеданс р°С1 для "чистой" жидкости, значительно выше, чем аналогичный параметр для пузырьковой смеси, возмущения, распространяющиеся вдоль центральной зоны "чистой" жидкости, быстро "съедаются" из-за разгрузки на боковой границе с пузырьковой жидкостью. Эта граница (г0 = /?0) для возмущений, распространяющихся в центральной зоне, фактически играет роль свободной поверхности. Для основных возмущений, распространяющихся по газожидкостному кольцу, граница (г0 = /?0) действует аналогично твердой стенке.

В § 3.2 рассмотрено взаимодействие волны типа "ступенька" со сферическими пузырьковыми кластерами.

На рис.2 представлены расчетные осциллограммы, иллюстрирующие динамику волн при различных радиусах кластера (сплошные линии соответствуют случаю наличия в трубе кластера, штриховые - его отсутствию). Первоначальный импульс типа "ступенька" инициируется заданием давления на границе г = 0 с амплитудой Др/ =0.3 МПа. Фрагмент (а) иллюстрирует расчетные осциллограммы, регистрируемые датчиком 01, расположенным в центре кластера, фрагмент (б) - показания датчика VI, расположенного на стенке трубы. Осциллограммы слева соответствуют радиусу кластера Лс1= 0.02 м, справа - случаю, когда радиус кластера равен радиусу трубы (ЯС, = ЯС = 0.1 м).

Видно, что когда радиус кластера мал (/?с; < Л</2), датчик, находящийся

Рис.2. Эволюция волн давления в трубе, заполненной жидкостью, содержащей пузырьковый кластер а) осциллограммы регистрируемые датчиком D1 в центре кластера; б) осциллограммы, регистрируемые датчиком D2 на стенке трубы Параметры системы: жидкость — вода, газ —воздух ао = 1 мм, <Xgo = 0.01, Rc = 0.1 м, zcl = 0.4 м.

в центре кластера (D1), иллюстрирует осцилляционный закон изменения давления. При этом амплитуда первого всплеска превышает амплитуду первоначальной волны более чем в 2 раза (р/> 15 атм). Датчик, расположенный на стенке, регистрирует сигнал с небольшими всплесками давления. В случае, когда радиус кластера равен радиусу трубы, оба датчика регистрируют значительное повышение амплитуды возмущений давления по сравнению с амплитудой инициирующей ударной волны. Анализ численных расчетов показал, что когда кластер достаточно мал (Rci < Rc/2), усиление давления внутри кластера происходит за счет двумерных и нелинейных эффектов, описанных в [4, 5]. Когда радиус кластера достаточно большой (Rd>Rc/2) и если zcl« fic,-(l+2-C//C)> где С - скорость волны в пузырьковой жидкости, увеличение амплитуды давления на участках, находящихся вблизи кластера, связано с нелинейным преломлением падающей волны на границах раздела сред, описанным в [6].

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

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

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

В § 4.1 показано, что принятая модель с учетом энерговыделения в газовой фазе содержит решение в виде стационарной волны. Для детонационных волн проведен сравнительный анализ амплитуды, формы и скорости со-литонов, полученных в двумерной осесимметричной постановке, со стационарным и экспериментальным (A.B. Пинаев, А.И. Сычев, 1987 г.) профилями и скоростями солитонов. Проведено сопоставление с результатами расчетов,

проведенных С.А. Жданом (2004 г.).

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

На рис. 3 (а) представлены эпюры давления в жидкости для динамики детонационной волны вдоль цилиндрического пузырькового кластера радиуса Я0 = 0.05 м. При расчетах на границе г = Яс задавались неотражающие условия в виде импедансного соотношения Ар1 = Лг>гр/С/, где Лр/ и Ли, - текущие значения возмущений давления и радиальной скорости для лагранжевой системы для границ расчетной области г = Де. Первоначальный импульс инициировался воздействием жесткого ударника с максимальной скоростью Ди0= 10.5 м/с и характерным временем воздействия 1,= 0.1 мс. Под воздействием граничного импульса температура внутри пузырьков в зоне вблизи поршня достигает значения Г., и возникает самоподдерживающаяся детонационная волна. Фрагмент (а) показывает, что детонационная волна распространяется вдоль трубчатого пузырькового кластера и сопровождается вол-

Рис.З. Распространение детонационной волны в трубчатом пузырьковом кластере Числовые указатели на кривых соответствуют времени в мс. Параметры системы для фрагментов (а)-(в)' жидкость -50% смесь глицерина с водой, газ - С2Н2+2 502 , а0- 1 мм, ag0 = 0.01; для фрагмента (г): жидкость - вода, газ - С2Н2+2.502, а0 = 2 мм, линии 1, 2 и 3 соответствуют ag0 = 0.005, а^ = 0.01 и ag0 = 0.02.

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

На фрагменте (б) приведены результаты численного эксперимента при тех же параметрах пузырьковой системы и инициирующего сигнала, что и на фрагменте (а), но с меньшим радиусом (/?0 = 0.025 м) пузырькового кластера. Под воздействием поршня в кластере вблизи границы поршня температура газа в пузырьках достигает значения Т*, и тем самым инициируется детонация. Но этот возникший очаг детонации не способен в дальнейшем досжать пузырьки, расположенные впереди, до температуры воспламенения Г», и детонация срывается. Видно, что на фоне детонационной волны в момент времени 0.3 мс волны давления после срыва в момент времени 1.2 мс практически незаметны. Основной причиной срыва является увеличение удельной боковой границы пузырькового кластера, через которую происходит акустическое излучение в окружающую жидкость, при уменьшении радиуса кластера.

Для иллюстрации роли акустического излучения в окружающую жидкость на фрагменте (в) представлены результаты расчетов для случая, когда пузырьковый кластер с радиусом = 0.025 м помещен в цилиндрический канал радиуса Я^ = 0.1 м с жесткими стенками. Для этой ситуации на поверхности цилиндрического канала ставится условие пг - 0. Видно, что в случае, когда акустическое излучение в окружающую жидкость экранировано жесткой стенкой, наблюдается устойчивое распространение самоподдерживающегося детонационного солитона. При этом амплитуды как детонационного солигона, так и сопровождающей волны, распространяющейся по жидкости, расположенной между кластером и стенками канала, примерно в два раза превышают значения, полученные при неотражающих условиях на границе Яс, хотя радиус пузырькового кластера в этом случае в 2 раза меньше (7?(г 0.025 м), чем для случая, приведенного на фрагменте (б) (7?о=0.05 м).

На фрагменте (г) приведена зависимость амплитуды детонационной волны от радиуса пузырькового кластера, расположенного в неограниченном объеме жидкости. Видно, что с увеличением радиуса газожидкостного кластера амплитуда детонационной волны носит немонотонный характер. Это связано с тем, что при увеличении во-первых, уменьшается удельная боковая поверхность, через которую происходит акустическое излучение в окружающую жидкость, во-вторых, уменьшается влияние волн, распространяющихся по "чистой" жидкости и поджимающих пузырьковый кластер. Начиная с некоторого значения радиуса кластера, влияние поджатая на процесс распространения детонационной волны становится несущественным.

В этом же разделе рассмотрена эволюция детонационных волн при пере-

ходе из пузырькового шнура в однородную пузырьковую жидкость. Условия при (7 = 0), 2 > 0, соответствующие исходному состоянию покоя неоднородной пузырьковой смеси в трубе, имеют вид

|(Ло < г0 < Яс)л (г < гр): а^ = 0, р = р%,

|((0 < /-0 < Ло)л (г < (г > гр}.аё = а&0,р = р^ - а§0) р& = р0,а = а0,уг = 0.

На рис. 4 представлены эпюры давления (а) и система изобар (б) в жидкости для динамики детонационной волны вдоль цилиндрического пузырькового кластера радиуса Л0 = 0 03 м и выхода волны в однородную пузырьковую жидкость. Видно, что в момент 0.4 мс детонационная волна, распространяющаяся по трубчатому кластеру, достигла зоны однородной пузырьковой жидкости и начала преломляться в эту область, и в момент 0.6 мс фронт волны имеет сферическую форму.

При уменьшении радиуса пузырькового кластера до Ло = 0.028 м под воздействием поршня в кластере также инициируется детонация, причем энергоемкости системы достаточно для того, чтобы вдоль кластера распространялась самоподдерживающаяся детонационная волна. Но при выходе со-литона в однородную пузырьковую жидкость, детонационная волна срывается, и постдетонационная волна распространяется как в неактивной пузырько-

Др/>МПа

Рис 4 Переход детонационной волны из газожидкостного шнура в однородную пузырьковую область. Параметры системы• а0 = 1.5мм, гр = 0.02м, остальные параметры такие же как для рис. 3 (фрагмент а)

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

В заключении представлены основные результаты и выводы, полученные в работе.

ОСНОВНЫЕ РЕЗУЛЬТАТЫ И ВЫВОДЫ РАБОТЫ

В работе создана и реализована численная модель динамики волн давления в пузырьковой жидкости в двумерной осесимметричной постановке. Показано, что использование предложенной разностной схемы для реализации численной модели позволяет вполне адекватно описывать волновую динамику для различных случаев расположения пузырьковой зоны. Приведено сравнение результатов расчетов с экспериментами, а также с результатами расчетов других исследователей. Применены условия неотражения от границ расчетной области, что позволило исключить влияние отраженных волн, и в то же время уменьшить время вычислений. Показана реализация принятой разностной схемы при различных способах задания первоначального сигнала: при задании на границе скорости среды или давления. Рассмотрена эволюция волн давления при неравномерном распределении пузырьков по сечению, а также волн, распространяющихся в газожидкостной среде кластерной структуры. Показано, что принятая модель с учетом энерговыделения в газовой фазе содержит решение в виде стационарной волны. Изучена динамика детонационных волн в трубчатом пузырьковом кластере, а также выход детонационной волны из кластера в однородную пузырьковую область. Исследовано влияние акустического излучения волн в зону "чистой" жидкости, окружающей газожидкостный трубчатый кластер, на параметры детонационной волны и возможность ее срыва.

По результатам исследований установлено:

1) при эволюции волны типа "ступенька" в трубе, заполненной жидкостью, содержащей трубчатый пузырьковый кластер, из-за фокусировки волны в область пузырькового кластера образуется волна, распространяющаяся вдоль кластера, амплитуды колебаний в которой могут многократно превышать амплитуду первоначальной волны;

2) при расположении пузырьков в кольцевом объеме, на участках стенки трубы, находящихся вблизи торца, через который происходит воздействие граничным давлением, из-за фокусировки возмущений в более сжимаемый пузырьковый слой, амплитуда давления может значительно превышать соответствующие значения, которые реализуются при гомогенном распределении пузырьков;

3) при воздействии жестким ударником на торец трубы, заполненной

газожидкостной смесью, в случае кольцевого режима на начальном этапе может наблюдаться сильно неоднородное "П"-образное распределение давления Но постепенно, для времен, значительно превышающих времена прохождения волн давления расстояний порядка линейных масштабов поперечной неоднородности пузырьковой системы, происходит установление на поршне более однородного давления;

4) при распространении волны типа "ступенька" в трубе, содержащей на оси пузырьковый кластер, радиус которого меньше половины радиуса трубы, происходит фокусировка волны в кластере и, как следствие этого, многократное увеличение амплитуды первоначальной волны. В случае, когда радиус кластера больше половины радиуса трубы, и кластер расположен так, что время прохождения волны через кластер было намного больше времени прохождения волны по "чистой" жидкости расстояния между кластером и границей, где поддерживается граничное давление, увеличение амплитуды давления на участках, находящихся вблизи кластера, связано с накоплением энергии, которое происходит из-за задержки волн в кластере и многократного отражения волн от кластера и торцевой границы, в которой поддерживается давление типа "ступенька";

5) амплитуда детонационного солитона, распространяющегося по трубчатому пузырьковому кластеру, расположенному в неограниченном объеме жидкости, от радиуса кластера зависит немонотонно, что обусловлено влиянием фоновой волны, излучаемой солитоном в "чистую" жидкость;

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

Цитированная литература

1. Нигматулин Р.И. Динамика многофазных сред. - М.: Наука. - Т 1 -2.-1987.-360 с.

2. Донцов В. Е., Накоряков В. Е. Волны давления в газожидкостной среде с расслоенной структурой жидкость - пузырьковая смесь // ПМТФ. 2003. Т. 44, № 4. - С. 102-107.

3. Донцов В. Е. Распространение волн давления в газожидкостной среде кластерной структуры // ПМТФ. 2005. Т. 46, № 3. - С. 50-60.

4. Нигматулин Р.К, Шагапов В.Ш., Гималтдинов И.К., Галимзянов М.Н. Двумерные волны давления в жидкости, содержащей пузырьковые зоны // ДАН. 2001. Т. 378. № 6. С. 763-768.

5. Кедринский В.К., Шокин Ю.И., Вшивков В.А., Дудникова Г.И., Лазарева Г.Г. Генерация ударных волн в жидкости сферическими пузырьковыми кластерами И ДАН. 2001. Т.381, №6, С.773-776.

6. Нигматулин Р.И., Губайдуллин A.A., АхметовА.Т., Бекишев С.А., Михайлов E.H. Экспериментальное и теоретическое моделирование эффекта аномального усиления ударных волн в высоковязких жидкостях. // ДАН. 1996. Т. 346, № 1,С. 46-50.

Публикации по теме диссертации

1. Нигматулин Р. И., Шагапов В. Ш, Гималтдинов И. К, Баязито-ва А. Р. Распространение детонационных волн вдоль трубчатого пузырькового кластера, находящегося.в жидкости // Докл. РАН. - 2005- Т. 403, №4.-С. 478-481.

2. Гималтдинов И.К., Баязитова А.Р. Уединенная детонационная волна в пузырьковой жидкости // Сб. материалов 42-ой науч. конф. 10-17 апреля 2002 г. / Отв. ред. И.Е. Карпухин. - Стерлитамак: Стерлитамак. гос. пед. ин-т, 2002 - 308 е., С.246.

3. Переведенцева (Баязитова) А.Р. Динамика волн давления в жидкости с пузырьковыми зонами в случае цилиндрической симметрии // Труды XIII сессии Российского Акустического Общества. Москва. Т.2 - С. 32-35.

4. Переведенцева (Баязитова) А Р. Динамика детонационной волны в пузырьковой жидкости в случае осевой симметрии // Девятая Всероссийская научная конференция студентов - физиков и молодых ученых. 28 марта - 3 апреля 2003 года. Сб. тезисов. Часть 1, С.396-397.

5. Гималтдинов И.К, Переведенцева (Баязитова) А.Р. Динамика детонационных волн вдоль цилиндрического объема пузырьковой смеси, находящейся в «чистой» жидкости // Спектральная теория дифференциальных операторов и родственные проблемы: Труды международной конференции

(24-28 июня 2003 г., г. Стерлитамак) / Отв. ред. К.Б. Сабитов. - Уфа: Гилем, 2003. - Т. 3. - С.78 - 84.

6. Гималтдинов И.К., БаязитоваА.Р. Об эволюции волн давления в пузырьковой трубке, находящейся в "чистой" жидкости // Труды Математического центра имени Н.И.Лобачевского. Т. 21/ Казанское математическое общество. Лобачевские чтения - 2003 // Материалы третьей всероссийской молодежной научной школы-конференции. - Казань: Издательство Казанского математического общества, 2003. - С. 79-80.

7. ГималтдиновИ.К., БаязитоваА.Р Эволюция волн давления в пузырьковой жидкости при ступенчатом распределенииЗ пузырьков // ЭВТ в обучении и моделировании: Сб. науч. трудов: в 2-х ч./ Отв. ред. С.М. Усманов // Всероссийская научно-теоретическая конференция. 21-22 мая 2004 г. - Бирск: Бирск. гос.пед. ин-т, 2004. - Ч. 1. - С. 62-65.

8. Гималтдинов И.К., Баязитова А.Р. Динамика волн давления в трубе, содержащей пузырьковую жидкость со ступенчатым распределением пузырьков по сечению // Современные проблемы физики и математики: Труды Всероссийской научной конференции (16-18 сентября 2004 г., г. Стерлитамак)/ Отв. ред. К.Б. Сабитов. - Уфа.: Гилем, 2004. - Т.2. - С. 1926.

9. БаязитоваА.Р. Динамика волн давления в трубе, заполненной жидкостью при наличии в ней газожидкостного кластера // Современные проблемы физики и математики: Труды Всероссийской научной конференции (16-18 сентября 2004 г., г. Стерлитамак)/ Отв. ред. К.Б. Сабитов. - Уфа.: Гилем, 2004. - Т.2. - С.39-43.

10. Баязитова А.Р. Динамика детонационных волн вдоль цилиндрического объема пузырьковой смеси, находящейся в "чистой" жидкости конечного и бесконечного объема // Десятая Всероссийская научная конференция студентов - физиков и молодых ученых. 1-7 апреля 2004 года. Сб. тезисов. Часть 2, С.997-998.

11. Баязитова А.Р., Гималтдинов И.К. Волны давления в трубе, содержащей газожидкостные кластеры // Труды XV сессии Российского акустического общества. Нижний Новгород, Т. 1 Физическая акустика. Распространение и дифракция волн. Геоакустика., 2004, С. 155-159.

12. Баязитова А.Р. Выход детонационной волны из газожидкостного шнура в однородную пузырьковую жидкость // V Региональная Школа-Конференция для студентов, аспирантов и молодых ученых по математике и физике: Тезисы докладов. - Уфа: РИО БашГУ, 2005. - С.28.

13. Баязитова А.Р., Гималтдинов И.К. Волны давления в трубе, содержащей газожидкостный кластер // Труды XVI сессии Российского акустического общества. Москва, Т. 1. Физическая акустика. Распространение и дифракция волн. - Москва: ГЕОС, 2005. - С. 125-129.

14. Гималтдинов И.К, БаязитоваАР. Эволюция детонационных волн в пузырьковой трубе, находящейся в объеме "чистой" жидкости // Динамика сплошной среды. Сборник научных трудов. Выпуск 123. Акустика неоднородных сред. Новосибирск, 2005, С. 12—15.

15. Баязитова А.Р Распространение детонационных волн в пузырьковом шнуре // ЭВТ в обучении и моделировании: Сб. науч. трудов: в 2-х ч./ Отв. ред. С.М. Усманов // IV Регион, научно-методич. конференция. 16-17 декабря 2005 г. - Бирск: Бирск. гос. соц.-пед. акад, 2005. - 4.1. - С. 266270.

16. Баязитова А.Р., Гималтдинов И.К., Шагапов В.Ш. Волны давления в трубе, заполненной пузырьковой смесью с неоднородным распределением по сечению // Известия РАН. Механика жидкости и газа. - 2006 - № 3 -С. 65-76.

ЛША

»-75 83

/ 1 I

Подписано в печать 11.04.2006 г. Гарнитура «Тайме». Бумага ксероксная. Формат 60х801/16. Печать оперативная. Усл.-печ. л. 1,6. Заказ № 59/06. Тираж 100 экз.

Отпечатано в типографии Стерлитамакской государственной педагогической академии: 453103, Стерлитамак, пр. Ленина, 49.

Содержание диссертации автор исследовательской работы: кандидата физико-математических наук, Баязитова, Алина Разифовна

Оглавление.

Введение.

Глава 1 . Обзор экспериментальных и теоретических исследований по динамике нелинейных и детонационных волн в пузырьковых жидкостях. Основные уравнения.

1.1. Нелинейные волны в пузырьковой жидкости.

1.1.1. Обзор экспериментальных исследований по динамике нелинейных волн.

1.1.2. Обзор теоретических исследований по динамике нелинейных волн.

1.2. Детонационные волны в пузырьковых жидкостях.

1.2.1. Обзор экспериментальных исследований по динамике детонационных волн.

1.2.2. Обзор теоретических исследований по динамике детонационных волн

1.3. Основные уравнения для смеси жидкости с газовыми пузырьками для цилиндрической симметрии.

Введение диссертация по механике, на тему "Осесимметричные волны в пузырьковой жидкости"

Задача о распространении двумерных и детонационных волн в системах жидкость - пузырьки газа интересна как в плане теоретического исследования, так и с точки зрения практических приложений. При этом интересными и важI ными являются процессы, носящие нестационарный характер и составляющие предмет изучения волновой динамики многофазных систем. Химически активные (в пузырьках находится взрывчатый газ) и пассивные (газ в пузырьках не является взрывчатым) пузырьковые среды способны воспринимать "закачиваемую" энергию, поглощая относительно слабую внешнюю импульсную нагрузку, концентрировать ее в некоторой локальной области и переизлучать с существен-# ным увеличением амплитуды. Газожидкостные среды широко используются в различных отраслях промышленности (химической, нефтеперерабатывающей, | * металлургической и др.) для ведения технологических процессов (например, | жидкофазное окисление углеводородов) или как элемент системы обеспечения взрывобезопасности (жидкостные предохранительные затворы). Кроме того, пузырьковые жидкости могут найти свое применение в литотрипсии (разрушении почечных камней), в качестве активного элемента в сазере (акустическом аналоге лазера).

Пузырьковая жидкость с горючей смесью газов (вода с пузырьками грему-; • чего газа или углеводородов с кислородом) является взрывчатым веществом

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

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

В настоящее время активно ведутся исследования по динамике двумерных волн в пузырьковой жидкости (Кедринский В.К., Лазарева Г.Г., Накоряков В.Е., Донцов В.Е., ЖданС.А., Ляпидевский В.Ю., Губайдуллин A.A., Вахитова Н.К., Matsumoto Y., и др.).

Таким образом, актуальным является разработка математической модели, адекватно описывающей реальные процессы в пузырьковой жидкости в осесим-ф метричной постановке; разработка и программная реализация алгоритмов решения задач эволюции волн давления; верификация построенной модели по экспе-t риментальным данным и путем сравнения с расчетными данными других авторов; численное исследование эволюции волн давления при неравномерном по сечению распределении пузырьков; численное исследование эволюции волн давления в трубе, содержащей газожидкостные кластеры; численное исследование динамики двумерных детонационных волн в пузырьковой жидкости.

Цели работы.

• исследование эволюции волн давления в пузырьковой жидкости при не® однородном по сечению распределении пузырьков;

• изучение распространения волн давления в трубе, содержащей сферичещ ские пузырьковые кластеры;

• исследование динамики детонационных волн в трубчатом пузырьковом кластере.

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

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

Достоверность результатов. Достоверность полученных в рамках диссертационной работы результатов диссертационной работы обоснована использованием общих законов и уравнений механики сплошной среды и согласованием полученных решений в частных случаях с результатами, известными из литературы, а также проведением сравнительных тестовых расчетов. Критерием, определяющим достоверность результатов, полученных расчетным путем, и адекватность сформулированной физико-математической модели реальному процессу, является соответствие расчетных и экспериментальных данных. Было проведено сопоставление с экспериментами для волн давления в пузырьковой жидкости в случае неоднородного по сечению распределения газа, а также для группы пузырьковых кластеров. Для детонационных волн был проведено сравнение амплитуды и формы солитонов, полученных в разработанном комплексе программ, полученных при стационарной постановке задачи и экспериментальных. Сходимость численных методов проверена на последовательности измельчающихся сеток.

Публикации. Результаты диссертации опубликованы в 16 работах.

Апробация работы. Основные результаты работы были представлены на следующих конференциях и научных школах:

- на девятой Всероссийской научной конференции студентов — физиков и молодых ученых (Красноярск, 2003);

- на третьей Всероссийской молодежной научной школе-конференции (Казань, 2003);

- на Всероссийской научно-теоретической конференции «ЭВТ в обучении и моделировании» (Бирск, 2004);

- на XIII сессии Российского акустического общества (Москва, 2004);

- на десятой Всероссийской научной конференции студентов - физиков и молодых ученых (Москва, 2004);

- на XV сессии Российского акустического общества (Нижний Новгород,

2004);

- на V Региональной школе - конференции для студентов, аспирантов и молодых ученых по математике и физике (Уфа, 2005);

- на XVI сессии Российского акустического общества (Москва, 2005)

- на IV Региональной научно-методической конференции «ЭВТ в обучении и моделировании» (Бирск, 2005).

Кроме того, результаты диссертационной работы докладывались и обсуждались на научных семинарах кафедры прикладной математики и механики СГПА под руководством профессора В.Ш. Шагапова, кафедры математического анализа СГПА под руководством профессора К.Б. Сабитова и профессора И.А. Калиева, кафедры теоретической физики СГПА под руководством профессора А.И. Филиппова, на научных семинарах Института математики с ВЦ УНЦ РАН под руководством профессоров М.Д. Рамазанова, Н.Д. Морозкина и получили положительную оценку.

Структура и объем работы.

Диссертация состоит из введения, четырех глав, заключения и списка литературы. Диссертация изложена на 134 страницах, включает библиографический список из 111 наименований работ, 64 рисунков. Рисунки, формулы и таблицы нумеруются по главам.

Заключение диссертации по теме "Механика жидкости, газа и плазмы"

Заключение

По результатам исследований установлено:

1)при эволюции волны типа "ступенька" в трубе, заполненной жидкостью, содержащей трубчатый пузырьковый кластер, из-за фокусировки волны в область пузырькового кластера образуется волна, распространяющаяся вдоль кластера, амплитуды колебаний в которой могут многократно превышать амплитуду первоначальной волны;

2) при расположении пузырьков в кольцевом объеме, на участках стенки трубы, находящихся вблизи торца, через который происходит удар, из-за фокусировки возмущений в более сжимаемый пузырьковый слой, амплитуда давления может значительно превышать соответствующие значения, которые реализуются

122 при гомогенном распределении пузырьков;

3) при воздействии жестким ударником на торец трубы, заполненной газожидкостной смесью, в случае кольцевого режима на начальном этапе может наблюдаться сильно неоднородное "П"-образное распределение давления. Но постепенно, для времен, значительно превышающих времена прохождения волн давления расстояний порядка линейных масштабов поперечной неоднородности пузырьковой системы, происходит установление на поршне более однородного давления;

Список источников диссертации и автореферата по механике, кандидата физико-математических наук, Баязитова, Алина Разифовна, Тюмень

1. Айдагулов Р. Р., Хабеев Н. С., Шагапов В. Ш. Структура ударной волны в жидкости с пузырьками газа с учетом нестационарного межфазного теплообмена // Прикладная математика и техническая физика. - 1977. - № 3, С. 67-74.

2. Ахмадуллин Ф.Ф. Динамика детонационных волн в неоднородной пузырьковой жидкости: Дисс. .канд. физ.-мат. наук. Тюмень. - 2005. - 102 с.

3. Баязитова А.Р. Выход детонационной волны из газожидкостного шнура в однородную пузырьковую жидкость // V Региональная Школа-Конференция для студентов, аспирантов и молодых ученых по математике и физике: Тезисы докладов. Уфа: РИО БашГУ, 2005 - С.28.

4. Галимзянов М. Н., Гималтдинов И. К., Шагапов В. Ш. Двумерные волны даления в жидкости, содержащей пузырьки // Изв. РАН. МЖГ.-2002.- № 2. -С. 139-147.

5. Галимзянов М.Н. Динамика двумерных волн в пузырьковой жидкости: Дисс.к-та физ.-мат. наук. Уфа - 2004. - 110 с.

6. Гималтдинов И.К., Баязитова А.Р. Уединенная детонационная волна в пузырьковой жидкости // Сб. материалов 42-ой науч. конф. 10-17 апреля 2002 г. /ф Отв. ред. И.Е. Карпухин. Стерлитамак: Стерлитамак. гос. пед. ин-т, 2002, С.246.

7. Губайдуллин А. А., Ивандаев А. И., Нигматулин Р. И. Исследование нестационарных ударных волн в газожидкостных смесях пузырьковой структуры // Прикладная математика и техническая физика. 1978.- № 2. - С. 78-86.

8. Донцов В. Е., Накоряков В. Е. Волны давления в газожидкостной среде с расслоенной структурой жидкость—пузырьковая смесь // Прикладная математика и техническая физика. 2003. Т. 44, № 4. - С. 102-107.

9. Донцов В.Е., Накоряков В.Е. Структура волн давления в жидкости с газожидкостными кластерами. XV сессия Российского акустического общества. Т. 1 Физическая акустика. Распространение и дифракция волн. Геоакустика. С. 1551. Ф 159.

10. Донцов В. Е. Распространение волн давления в газожидкостной среде кла-^ стерной структуры // Прикладная математика и техническая физика. 2005.1. Т. 46, № 3. С. 50-60.

11. Ждан С. А., Ляпидевский В. Ю. Детонация в двухслойной пузырьковой среде // Физика горения и взрыва. -2002. -Т. 38, № 1. С. 123-128.

12. Ждан С. А. О стационарной детонации в пузырьковой среде // Физика горения и взрыва. -2002.-Т. 38, № 3. С. 85-95.

13. Ждан С. А. Детонация столба химически активной пузырьковой среды в• жидкости // Физика горения и взрыва. -2003. -Т. 39, № 4. С. 107-112.

14. Иорданский С. В. Об уравнениях движения жидкости, содержащей пузырьки газа // Прикладная математика и техническая физика. 1960.-№3. -С. 102-110.

15. Ильгамов М.А., Гильманов А.Н. Неотражающие условия на границах расчетной области. М.: ФИЗМАТЛИТ, 2003. - 240 с.

16. Кедринский В. К. Ударные волны в жидкости с пузырьками газа // Физика горения и взрыва. 1980.-Т. 16, № 5. - С. 14-25.

17. Кедринский В. К. Гидродинамика взрыва: эксперимент и модели. -Новосибирск: Издательство СО РАН, 2000. 435 с.

18. Кедринский В. К., Шокин Ю. И., Вшивков В. А., Дудникова Г. И., Лазаре-Ф ва Г. Г. Генерация ударных волн в жидкости сферическими пузырьковыми кластерами // Докл. РАН.- Т.381. № 6, 2001, С.773-776.

19. Когарко Б. С. Об одной модели кавитирующей жидкости // ДАН СССР. -1961. Т. 137, №6. -С. 1331-1333.

20. Когарко Б. С. Одномерное неустановившееся движение жидкости с возникновением и развитием кавитации // ДАН СССР. 1964. Т. 155, № 4. - С. 779782.

21. Кочин Н.Е., Кибель И.А., Розе Н.В. Теоретическая гидромеханика, М.: Физматгиз. 1963. Часть 1., - 584 с.

22. ЗБ.Кузнецов В. В., Донцов В. Е. Ударные волны умеренной амплитуды в двухфазной среде. Гидродинамические течения и волновые процессы // Под ред. В. Е. Накорякова. Новосибирск.: ИТФ СОАН СССР. 1983., С. 29-34.

23. Кузнецов Н.М., Копотев В.А. Структура волны и условие Чепмена-Жугепри гетерогенной детонации в жидкостях с пузырьками газа // Докл. АН СССР. -1989. Т. 304, №4. -С.850-853.

24. Лазарева Г.Г. Численное моделирование усиления ударных волн в пузырьковых средах: Дисс.канд. физ.-мат. наук. Новосибирск. 2003. 102 с.

25. Ландау Л.Д., ЛифшицЕ.М. Теоретическая физика: Учебное пособие. В 10 т. Т. V. Ч. 1. Статистическая физика. 3-е изд., перераб. и допол. — М.: Наука. 1976.-584 с.

26. Ляпидевский В.Ю. О скорости пузырьковой детонации // Физика горения и взрыва. 1990. - № 4. С. 139-140.

27. Ляпидевский В.Ю. Структура детонационных волн в многокомпонентных пузырьковых средах // Физика горения и взрыва 1997.-Т. 33, № 3. - С. 104-113.

28. Ляхов Г. М. Ударные волны в многокомпонентных средах // Изв. АН СССР. Сер. механика и маш. 1959. - № 1. - С. 20-25.

29. Малых Н. В., Огородников И. А. О применении уравнения Клейна- Гордона для описания структуры импульсов сжатия в жидкости с пузырьками газа // В сб.: Динамика сплош. среды. Вып. 29. Новосибирск.-1977. - С. 143-148.

30. Малых Н. В., Огородников И. А. Скорость и затухание импульсов большой амплитуды в слое жидкости с пузырьками газа // В сб.: Переход ламинарн. по-граничн. слоя в турбулент. Двухфазные потоки. Новосибирск. -1978. С. 38-51.

31. Мейз Дж. Теория и задачи механики сплошных сред. М.: Мир. 1974. -318 с.

32. Милн-Томсон JI.M. Теоретическая гидродинамика М.: Мир, 1964. — 660 с.

33. Накоряков В. Е., Донцов В. Е. Взаимодействие ударной волны со сферическим пузырьковым кластером в жидкости //.Докл. РАН. 2003. -Т. 391, №2. -С. 199-202.

34. Накоряков В.Е., Донцов В.Е., Волны давления в расслоенной среде жидкость газожидкостная среда // ДАН. - 2002. - Т. 386, № 1, С. 48-50

35. Накоряков В. Е., Покусаев Б. Г., Шрейбер И. Р. Распространение волн в газо- и парожидкостных средах// Новосибирск: ИФТ СОАН СССР. 1983, 237 с.

36. Накоряков В. Е., Покусаев Б. Г., Шрейбер И. Р. Волновая динамика газо- и парожидкостных сред. М.: Энергоатомиздат. 1990.

37. Нигматулин Р. И. Мелкомасштабные течения и поверхностные эффекты в Щ гидромеханике многофазных сред // ПММ. 1971.- № 3. - С. 541-563.

38. Нигматулин Р. И. Основы механики гетерогенных сред. М.: Наука. 1978. -4 336 с.

39. Нигматулин Р. И. Эффекты и их математическое описание при распространении волн в пузырьковых средах // В сб.: Избр. вопр. соврем, мех., Т. 4.1. М.- 1981.-С. 65-89.

40. Нигматулин Р. И. Динамика многофазных сред. М.: Наука. Т1, 2, 1987. -360 с.

41. Нигматулин Р.И., Губайдуллин A.A., Ахметов А.Т., Бекишев С.А., Михай-Ф лов E.H. Экспериментальные и теоретические моделирования эффекта аномального усиления ударных волн в высоковязких жидкостях// ДАН, 1996, Т. 346, №1, С 46-50.

42. Нигматулин Р. И., Хабеев Н. С., Шагапов В. Ш. Об ударных волнах в жидкости с пузырьками газа // Докл. АН СССР. 1974, Т. 214, № 4. - С. 779-782.

43. Нигматулин Р. И., Шагапов В. Ш., ВахитоваН.К. Проявление сжимаемости несушей фазы при распространении волны в пузырьковой среде // Докл. РАН. 1989.- Т. 304, № 5. - С. 1077-1088.

44. Нигматулин Р. И., Шагапов В. Ш., Гималтдинов И. К., Баязитова А. Р. Распространение детонационных волн вдоль трубчатого пузырькового кластера, находящегося в жидкости // Докл. РАН. 2005.- Т. 403, № 4. - С. 1-4.

45. Нигматулин Р. И., Шагапов В. Ш., Гималтдинов И. К., Галимзянов М. Н. Двумерные волны давления в жидкости, содержащей пузырьковые зоны // Докл. РАН. 2001.- Т. 378, № 6. - С. 763-767.

46. Николаев Ю.А., Топчиан М.Е. Расчет равновесных течений в детонационных волнах в газах // Физика горения и взрыва. -1977.-Т.13, № 3. С. 393-404.

47. Переведенцева А.Р. Динамика волн давления в жидкости с пузырьковыми зонами в случае цилиндрической симметрии // Труды XIII сессии Российского Акустического Общества. Москва. Т.2 С. 32-35.

48. Переведенцева А.Р. Динамика детонационной волны в пузырьковой жидкости в случае осевой симметрии // Девятая Всероссийская научная конференция студентов-физиков и молодых ученых. 28 марта-3 апреля 2003 года. Сб. тезисов. Часть 1, с.396-397.

49. Пинаев А. В., Сычев А. И. Обнаружение и исследование самоподдерживающихся режимов детонации в системах жидкое горючее-пузырьки окислителя // Докл. АН СССР. 1986. -Т. 290, № 3. - С. 611-615.

50. Пинаев А. В., Сычев А. И. Влияние физико-химических свойств газа и жидкости на параметры и условия существования волн детонации в системах жидкость-пузырьки газа // Прикладная математика и техническая физика. -1987.- № 6. С. 76-84.

51. Пинаев А. В., Сычев А. И. Обнаружение и исследование самоподдерживающихся режимов детонации в системах жидкое горючее пузырьки окислителя // Докл. АН СССР.-1986. -Т. 290, № 3.- С. 611-615.

52. Пинаев А. В., Сычев А. И. Влияние физико химических свойств газа и жидкости на параметры и условия возникновения детонационных волн в системах жидкость - пузырьки газа // Физика горения и взрыва.- 1987. -Т. 23, № 6. -С. 76-84.

53. Пинаев А. В. Передача пузырьковой детонации через слой инертной жидкости// Физика горения и взрыва.- 2004. -Т. 40, № 2. -С. 105-110.

54. Самарский А. А., Попов Ю. П. Разностные схемы газовой динамики. М.: Наука. 1975.-352 с.

55. СычевА. И. Воспламенение систем жидкость-пузырьки газа ударной волной // Физика горения и взрыва. -1985.- № 1. -С. 130-134.

56. Сычев А. И., Пинаев А. В. Самоподдерживающаяся детонация в жидкостях с пузырьками активного газа // Прикладная математика и техническая физика. -1986. № 1.

57. Сычев А. И. Пузырьковая детонация в полидисперсных средах // Физика горения и взрыва. -1997. -Т. 33, № З.-С. 114-119.

58. Сычев А. И. Переход волны пузырьковой детонации в жидкость.// Физика горения и взрыва. -2002. -Т. 38, № 2. -С. 99-103.

59. Сычев А. И. Структура волны пузырьковой детонации// Физика горения и взрыва. -1994. -Т. 30, № 4. -С. 119-124.

60. Сычев А. И. Влияние размера пузырьков на характеристики волн детонации.// Физика горения и взрыва. -1995. -Т. 31, № 5. -С. 83-91.

61. Троцюк A.B., Фомин П.А. Модель пузырьковой детонации // Физика горения и взрыва. 1992. Т.28, №4. - С. 129-136

62. Численные методы в задачах физики взрыва и удара: Учебник для втузов./ A.B. Бабкин, В.И. Колпаков, В.Н. Охитин, В.В. Селиванов. М.: Изд-во МГТУ им. Н.Э.Баумана, 2000. - 516 с.

63. Шагапов В.Ш., Вахитова Н.К. Волны в пузырьковой системе при наличии химических реакций в газовой фазе. // Физика горения и взрыва. 1989. № 6. -С. 14-22.

64. Шагапов В.Ш., Абдрашитов Д.В. Структура волн детонации в пузырьковой жидкости. // Физика горения и взрыва. 1992.- №.4 - С. 89-95.

65. Biesheuvel A., van Wijngaarden L. Two fhase flow equations for a dilute dispersion of gas bubbles in liquid//J. Fluid Mech. 1983, № 148.-P. 301-318.

66. Beylich A.E., Gulhan A. Waves in reactive bubbly liquids // Proc. IUTAM Symp on Adiabatic Waves in liquid Vapor Systems. Gettingen. -FRG. -1989. -P. 3948.

67. Campbell J., Pitcher A. S. Shock waves in a liquid containing gas bubbles // Proc. Rog. Soc. London. 1959. -A 234.- №1235. - P. 534-545.

68. Derzho O. G., Malykh N. Y. Formation of strong pressure pulses reflected from water—bubble layers // Arch. Mech. 1942. -V. 4-5. - P. 463-473.

69. Drumheller P. S., Bedford A. A theory of bubbly fluids // J. Acoust. Soc. Amer. 1979.- V. 66, № 1. - P. 197-208.

70. Enever K. J. Further investigation of the theory of shock waves in bubbly liquids // Bull. etud. et rech. 1977. -A, № 2. - P. 29-46.

71. Fujiwara K. et al. New methods for generation cylindrical imploding shock // Proc. 19th Intern. Symp. on Shock Waves. Marseille, France, 26-30 July, 1993. V. 4. P. 81-86.

72. Isuzukawa K., Horiuchi M. Experimental and numerical studies of blast wave focusing in water // Proc. 19th Intern. Symp. on Shock Waves. Marseille, France, 26-30 July, 1993. V.3.P. 347-350.

73. KamedaM., ShimauraN., Higashino F., Matsumoto Y.// Phys.Fluids. 1998. V.10. № 10, P. 2661-2668.

74. Nagoya H. et al. Underwater shock wave propagation and focusing in inho-mogeneous media tube // Proc. 19th Intern. Symp. on Shock Waves. Marseille, France, 26-30 July, 1993. V. 3. P. 439-444.

75. Hasegava T., Fujiwara T. Detonation in oxyhydrogen bubbled liquids // Proc. 19th Intern. Symp. on Combustion. Haifa, 1982.

76. Kashinsky O.N., Radnin V.V. // Intern. J. Multiphase Flow. 1999. V. 25. № 1. P. 109-138.

77. Kashinsky O.N., Timkin L.S. // Experim. in Fluid. 1999. V. 26. № 2. P. 305314.

78. Kedrinskii V.K. The role of cavitation effects in the mechanisms of destruction and explosive processes // J. Shock Waves. 1997. V. 7, N 2. P. 63-76.

79. Kedrinskii V.K., Mader Ch. Accidential detonation in bubbly liquids // Proc. 16th Intern. Symp. on Shock Tube and Waves / H. Groenig (Ed.). 1987. P. 371-376.

80. Kedrinskii V.K., Mader Ch. On the velocity of bubble detonation // Proc. 13th Intern. Symp. on Nonlinear Acoustics. Bergen, Norway, 1993. P. 442-447.

81. Kuwahara M. et al. The problems of focused shock waves effect on biological tissues // Proc. 18th Intern. Symp. on Shock Waves. Sendai, Japan, 21-26 July, 1991. V. 1. P. 41-48.

82. Martin C. S., Padmanablan M., WiggertB. G. Pressure waves propagation in two-phase bubbly air-water mixtures // Proc. 2nd Int. Conf. press Surges. London, 1976, Cranfield, 1977. - P. 1-16.

83. Martin C. S., PadmanablanM., WiggertB. G. Pressure waves propagation in two-phase bubbly air-water mixtures // Bull. Dir. etud. et. resh. 1977, A, №2. -P. 47-66.

84. Miksis M. J., TingL. Effects of bubbly layers on wave propogation // J. Acoust. Soc. Amer. 1989, V. 86, № 6. - P. 2349-2385.

85. Noordzij L. Shock waves in bubble-liquid mixture // Phys. Comm. Twente Univ. Techn.-1971, V. 3, №1.-P. 51.

86. Noordzij L., Wijngaarden L.van. Relaxation effects, caused by relative motion, on shock waves in gas-bubble/liquid mixtures // J. Fluid Mech. 1974, V. 66, № l.-P. 115-143.

87. Parkin B. R., GilmoreF. R., BrodeH. L. Shock waves in bubbly water. Memorandum RM-2795-PR. Abridged. - 1961.

88. Rath H. J. Unsteady pressure waves and shock waves in elastic tubes containing bubbly air-water mixtures // Acta Mech. 1981. V. 38, № 1-2. - P. 1-17.

89. Stuka C. et al. Nonlinear transmission of focused shock waves in nondegassed water//Ibid. P. 445-448.

90. Scarinci T., BassinX., Lee J., Frost D. Propagation of a reactive waves in a bubbly liquid // Proc. 18th ISSW / K. Takayama (Ed.). V. 1. P. 481-484.

91. Takayama K. High pressure generation by shock wave focusing in ellipsoidal cavity I I Proc. Intern. Workshop on Shock Wave Focusing / K. Takayama (Ed.). Sen-dai, Japan, 1989. P. 217-226.

92. Van Wijngaarden L. On the equations of motion for mixtures of liquid and gas bubbles. J. Fluid Mech. 1968, v. 33, № 3. . p. 465-473.

93. Watanabe M. et al. Shock wave focusing in a vertical annular shock tube // Proc. 19th Intern. Symp. on Shock Waves. Marseille, France, 26-30 July, 1993. V. 4. P. 99-104.