Нелинейная рефракция относительно слабых ударных волн в газах и газожидкостных средах тема автореферата и диссертации по механике, 01.02.05 ВАК РФ
Матутин, Александр Александрович
АВТОР
|
||||
кандидата физико-математических наук
УЧЕНАЯ СТЕПЕНЬ
|
||||
Саратов
МЕСТО ЗАЩИТЫ
|
||||
2007
ГОД ЗАЩИТЫ
|
|
01.02.05
КОД ВАК РФ
|
||
|
На правах рукописи
Матутин Александр Александрович
НЕЛИНЕЙНАЯ РЕФРАКЦИЯ ОТНОСИТЕЛЬНО СЛАБЫХ УДАРНЫХ ВОЛН В ГАЗАХ И ГАЗОЖИДКОСТНЫХ СРЕДАХ
01 02 05 — механика жидкости, газа и плазмы
Автореферат
диссертации на соискание ученой степени кандидата физико-математических наук
Саратов-2007
003163027
Работа выполнена на кафедре вычислительного эксперимента в механике механико-математического факультета Саратовского государственного университета им. Н Г. Чернышевского
Защита состоится 12 ноября 2007 года в 17 час. 30 мин на заседании диссертационного совета К 212.243 02 при Саратовском государственном университете им. НГ. Чернышевского по адресу 410012, г Саратов, ул. Астраханская, д 83, Саратовский государственный университет, механико-математический факультет.
С диссертацией можно ознакомиться в Научной библиотеке Саратовского государственного университета
Научный руководитель
Доктор физико-математических наук, профессор
Шиндяпин Георгий Петрович
Официальные оппоненты Доктор технических наук, профессор
Есин Александр Иванович
Кандидат физико-математических наук, доцент
Шевырев Сергей Павлович
Ведущая организания
Балтийский государственный технический
университет им Д.Ф Устинова
(БГТУ - «Военмех») (г Санкт-Петербург)
Автореферат разослан
2007 года.
Ученый секретарь диссертационного совета к ф -м. н, доцент
В.В Корнев
ОБЩАЯ ХАРАКТЕРИСТИКА РАБОТЫ Актуальность темы. Теоретическое исследование процессов нелинейной рефракции ударных волн (УВ) в газах и газожидкостных средах (ГЖС) представляет одну из фундаментальных проблем современной механики жидкости и газа Решение этой проблемы важно для понимания многих физических процессов, развития сверхзвуковой авиации, проектирования трубопроводов для транспортировки топливных смесей, движения тел в газожидкостных средах, в которых скорость звука мала и др Теоретическая важность проблемы обусловлена нелинейным характером основных уравнений и сложным разрывным характером решений, в силу чего методы и их применение при исследовании задач имеют универсальный характер и связаны с решением общих проблем выявления структуры обобщенных решений задач математической физики
В работе исследуются процессы нелинейной рефракции ударной волны (УВ), возникающие при взаимодействии ударной волны, относительно малой интенсивности (абсолютная интенсивность волн может быть велика), со свободной поверхностью, разделяющей различные газожидкостные среды, когда за фронтами УВ возникают области резких изменений параметров -области коротких волн, и ударные нагрузки многократно возрастают В этих областях процесс существенно нелинеен, что вызывает значительные трудности при анализе поставленной задачи Особое внимание уделяется нахождению и исследованию различных режимов существования нелинейной рефракции УВ, а также численному анализу поставленной задачи
Пристальное внимание к проблеме отечественных (С А Христианович, А А Гриб, Б И Заславский, Г П Шиндяпин, В К Кедринский и др) и зарубежных исследователей (L F Henderson, A. Sakurai, К Takayama, G Ben-Dor и др) на протяжении более 50 лет привело к появлению разнообразных альтернативных точек зрения на природу явления Однако в настоящее время все отчетливее вырисовывается необходимость исследования процессов рефракции в целом, для построения полной картины исследуемого явления
Настоящая работа развивает и дополняет положения асимптотической теории коротких волн, опирающейся на решение внутренних краевых задач для областей нелинейных взаимодействий (коротких волн), которые позволяют учесть влияние потока в целом на образующиеся ударно-волновые структуры, изучить характерные особенности и закономерности процессов нелинейной рефракции УВ.
Цель и задачи исследования. Целью работы является аналитическое и численное исследование различных режимов нелинейной рефракции относительно слабых УВ в газах и газожидкостных (пузырьковых) средах Исходя из этого, в работе решались следующие задачи 1
• разработка аналитических моделей нелинейной рефракции относительно слабых УВ,
• классификация режимов нелинейной рефракции УВ, анализ областей существования различных режимов нелинейной рефракции УВ;
• получение явных аналитических зависимостей для основных параметров, характеризующих возникающие ударно-волновые структуры,
• аналитическое описание границ областей существования различных режимов нелинейной рефракции УВ;
• численное решение краевых задач рефракции УВ, на основе построенных аналитических моделей,
• анализ полей давлений и скоростей при различных режимах нелинейной рефракции УВ
Методика исследований основана на использовании асимптотической теории коротких волн при постановке краевых задач для областей нелинейных взаимодействий (коротких волн), разработке и формировании аналитических моделей рефракции, соответствующих краевым задачам
Основное внимание уделяется аналитическим методам анализа, основанным на выделении структурных особенностей течений Научная новизна работы заключается в следующем
• развит новый аналитический метод исследования ударно-волновых структур и потоков при нелинейной рефракции относительно слабых УВ в газах и газожидкостных средах,
• исследованы различные режимы нелинейной рефракции УВ, найдены аналитические выражения для границ областей существования различных режимов нелинейной рефракции УВ,
• найдены области существования неклассических режимов рефракции УВ, которые расширяют физическое представление о процессах нелинейной рефракции УВ, позволяя исследовать случаи рефракции с невырожденным фронтом преломленной УВ,
• получил развитие численный метод решения краевых задач нелинейной рефракции УВ, позволяющий рассчитать поля давлений и скоростей, а также, ударно-волновые структуры, характерные для различных режимов нелинейной рефракции УВ,
• исследованы ударно-волновые структуры и течения в областях нелинейных взаимодействий. Выявлены качественные закономерности и особенности процессов нелинейной рефракции УВ, показавшие физическую адекватность разработанных аналитических и численных подходов к решению поставленной задачи
Достоверность результатов подтверждается непротиворечивостью полученных аналитических результатов с имеющимися физическими представлениями и экспериментальными данными, согласованностью результатов для рассчитанных ударно-волновых структур, полей давления и скоростей с известными результатами численного решения соответствующих краевых задач, обоснованностью используемых методов исследований
Практическая ценность работы. Результаты исследований расширяют представления о физических процессах нелинейной рефракции УВ в газах и газожидкостных средах, развивают аналитические методы и методы численного расчета течений с относительно слабыми УВ.
Результаты исследований могут быть использованы в учебном процессе при подготовке учебных пособий и чтении специальных курсов по динамике ударных волн в газах и газожидкостных средах
Апробация результатов. Результаты исследований докладывались и обсуждались на международных, всероссийских и вузовских конференциях и семинарах Среди них.
• XX, XXI Международный семинар по струйным, отрывным и нестационарным течениям Санкт-Петербург, 2004 г ; Новосибирск, 2007 г,
• XXXIV Уральский семинар по механике и процессам управления Миасс, 2004 г,
• XXIV, XXVII Российская школа по проблемам науки и технологий Миасс, 2004,2007 гг,
• XLII Международная научная конференция «Студент и научно-технический прогресс» Новосибирск, 2005 г,
• IV, V Молодежные научные школы-конференции «Лобачевские чтения» Казань, 2005; 2006 гг ;
• Международная конференция «Континуальные алгебраические логики, исчисления и нейроинформатика в науке, технике и экономике». Ульяновск, 2005; 2006 гг,
• Ежегодные научные конференции Саратовского государственного университета «Актуальные проблемы математики и механики» Саратов, СГУ 2003,2004, 2005, 2006, 2007 гг ;
• Специальные семинары кафедры вычислительного эксперимента в механике Саратовского государственного университета 2003-2007 гг
СОДЕРЖАНИЕ РАБОТЫ
Диссертация состоит из введения, трех глав, выводов по главам, заключения и списка использованной литературы, изложена с учетом рисунков на 111 страницах Работа содержит 18 рисунков Список использованной литературы включает 104 наименования
Во введении обоснована актуальность работы, сформулированы цель и задачи исследования, научная новизна, практическая значимость работы, кратко излагается содержание работы и основные результаты, выносимые на защиту
В главе 1 изложены общие основы математического моделирования процессов рефракции относительно слабых ударных волн (УВ) в газах и газожидкостных пузырьковых средах, обоснована методика исследований
В п. 1.1 рассмотрен класс задач нелинейной рефракции УВ Приведены схемы регулярной (RR), нерегулярной (NR) рефракций с волной разрежение; рефракции с УВ, замыкающей волну разрежение (RRV), рефракции с отраженной УВ, возникающей при падении плоской инициирующей УВ относительно малой интенсивности Р10 =(Pi— Ра)!В0, fi0=Poco под углом а к вертикали на свободную поверхность, разделяющую покоящиеся газожидкостные среды с газосодержаниями у+ ,у~ Отмечено, что указанные
режимы наблюдаются физически при падении УВ на свободную поверхность со стороны более плотной среды Сформулированный класс задач допускает все основные режимы рефракции, наблюдаемые при экспериментальных исследованиях.
В п. 1.2 приведен анализ применимости модели равновесной двухфазной газожидкостной среды для описания процессов рефракции относительно слабых УВ. Исследована зависимость параметров покоящейся газожидкостной среды (плотности р0, скорости звука с0), скорости УВ Ю0, числа Маха М0, термодинамических параметров Ло(у), ¿0(у) от газосодержания среды Отмечено, в частности, что значения параметра ^о(У) (А) =Р0^0^Р0с0' =1 = ¿0&р/р0), характеризующего влияние газосодержания на интенсивность УВ, для практически важного случая рефракции на поверхности, разделяющей газовую и газожидкостную среду, в основном диапазоне газосодержания Ю-7 ¿у¿Ю"4близко к значениям этого параметра для чистого газа. Поэтому в целом режимы рефракции и конфигурации УВ, возникающие при взаимодействиях газожидкостных пузырьковых средах, оказываются сходными с режимами рефракции волн на поверхности, разделяющей чистые газы — газоподобны
В п. 1.3 приведен анализ современного состояния исследований проблем нелинейной рефракции и общей теории взаимодействий относительно слабых УВ в газах и газожидкостных средах. Особое внимание уделено развитию аналитических и численных методов исследований, рассмотрено развитие асимптотической теории коротких волн в работах Гриба А А, Рыжова О С, Христиановича С А. (1955-1960), позволяющей исследовать течение в областях резких изменений параметров за фронтами взаимодействующих УВ Развитая в трудах Христиановича С А., Гриба А А, Заславского Б И, Рыжова О С, Рябинина А Г теория коротких волн охватывает неустановившиеся пространственные течения и позволяет изучить картину нелинейного взаимодействия УВ с различными препятствиями Для нелинейной системы уравнений коротких волн в работах Заславского Б И., Шиндяпина Г.П (19621971) построены достаточно широкие классы точных решений, расширяющие возможности анализа Развитие теории коротких волн связано с постановкой и анализом краевых задач нелинейных взаимодействий УВ, разработанной в работах Шиндяпина Г.П, Ковалева АД и др. (1984-1987) аналитических и численных методов решений Результаты исследований легли в основу асимптотической теории взаимодействий УВ (1984-1990), результаты которой в случае относительно слабых УВ в отличие от результатов, например, теории Уизема (1957) значительно лучше согласуются с экспериментальными данными Успехи теории коротких волн в разработке математических моделей взаимодействия и отражения УВ (Шиндяпин Г.П., 1990-1997), анализе проблемы взаимодействий УВ в условиях возникновения парадокса Неймана (Шиндяпин ГП, Гамаюнова ЕН, 2001-2004) позволяют использовать математический аппарат теории для анализа процессов нелинейной рефракции УВ
Значительную роль в развитии общей теории УВ сыграли идеи методы анализа ударно-волновых структур, предложенные в работах исследователей отечественной школы Баженовой Т В , Гвоздевой Л Г , Дулова В Г, Ускова В Н , Фонарева Л.С, Жигапко Е Ф и их учеников, а также экспериментальные исследования Заславского Б И, Сафарова Р А, Семенова А Н, Сыщиковой М.П. и др
Современное состояние теории взаимодействий УВ и теории рефракции УВ характеризуется бурным развитием. В этом направлении успешно работает ряд научных школ Австралия, Япония, США и др Ведутся интенсивные экспериментальные и численные исследования Следует, однако, отметить, что для исследуемого класса задач рефракции (падение УВ на свободную поверхность со стороны более плотной среды) экспериментальных результатов явно недостает В п 13 обсуждаются наиболее интересные результаты исследований
В работах автора (2003-2007) на базе теории коротких волн разработан аналитический подход, позволяющий сформулировать математические модели для различных режимов нелинейной рефракции (ЯЯ, N11, Ш1У, ЯУ) Приведены исследования характерных ударно-волновых структур и параметров в пространстве параметров подобия задачи, найдены аналитические выражения для границ областей существования различных режимов рефракции УВ Приведены численные и аналитические исследования полей давления и скоростей, включающие анализ течений в областях больших градиентов в окрестности точки взаимодействия падающей волны со свободной поверхностью Сравнение полученных результатов с известными экспериментальными, аналитическими и численными результатами показало физическую достоверность развитых моделей
Проведенный обзор обосновывает общую методику исследований процессов нелинейной рефракции относительно слабых УВ в газах и газожидкостных средах с помощью развитого аналитического метода исследований, основанного на выделении структуры течений (схем) при различных режимах взаимодействий
В п. 1.4 рассмотрены проблемы математического моделирования процессов нелинейной рефракции УВ Среди этих проблем нелинейность и нестационарность процессов рефракции, влияние газосодержания сред на структуру течения - режимы взаимодействия Особая сложность характерна для постановки и решения краевых задач для областей возмущений возникающих за фронтами преломленных, отраженных УВ или волнами разрежения, когда фронты УВ заранее неизвестны и находятся вместе с решениями краевых задач
В п. 1.5 дан анализ применимости модели равновесной двухфазной газожидкостной среды для описания ударно-волновых процессов Исследования взаимодействий УВ в газожидкостных пузырьковых средах (1974-1985) показывают, что эти процессы можно описать аналогично процессам взаимодействий в газах Для относительно слабых УВ используя модель равновесной газожидкостной среды, удается свести задачу к решению
нелинейного уравнения для потенциала скорости с параметром, выражающим зависимость решения от газосодержания среды В предельном случае это уравнение переходит в известное уравнение для адиабатического движения газа Лайтшпла
В п. 1.6 проведен анализ общей постановки краевых задач для описания процессов рефракции УВ в идеальных газовых или газожидкостных средах, а также краевых задач для потенциальных адиабатических течений газовых и газожидкостных сред, характерных для УВ относительно малой интенсивности Рю=(Р1-Ро)/В0, В0 = р0с^ Общие уравнения динамики при автомодельных движениях в переменных
x/c0t = !; = /• eos 8, у / c0t = r) = г sin 0 приводятся к уравнению для потенциала скоростей, которое вместе с уравнение состояния и интегралом Лагранжа составляют систему уравнений потенциальных адиабатических движений газовых и газожидкостных сред
Ф(х, у, i) = c0t /(г, 9) , — = fr , - = -i/e,
С о СQ г
с*2
с*2
fn +\fr = +Дг/е {¿~\fr +~/08]/9-2{r~fr)U
^с0>/ В0 2
G = /-rfr +~/г2 +т~2"/в ; (1)
Ро г 2г
= 2д + 1 р0c¿ | Ь PqCq d2
) Ро СУ, +cv„Y
2(а +1) ро ' 2(д + 1) р0 dp2{F(p))po
(wl
a=(l + y)& , Ь = Ц- ,
уД ' 1 + у
Здесь «,у- компоненты скорости, с- скорость звука, р- давление, р-плотность среды Зависимость решения от относительного массового газосодержания проявляется через параметры Л0, р0, с0 определяемые термодинамическими формулами (сп , ст , Я- постоянные). В предельном случае чистого газа (у -» оо) уравнения (1) переходят в известные уравнения Лайтхилла для потенциальных адиабатических течений газа.
При постановке краевых задач для уравнений (1) в качестве граничных условий выступают условия динамической совместности на фронтах УВ (индекс 1 - соответствует параметрам перед фронтом) при г = А(9)
AP=PL Ро
Р, р М/»)
РО *2-А'2 ~рос2 ' Ргм^)-/''
„ 1 + 2д ■ „ _2(1 + а) />i Pi
И[ =-; л2 =--—г- , И3 =-г- ,
2а 1 + 2а (1 + у)РоСо
условия на линиях слабого разрыва, где параметры потока непрерывны (/, =/1г, /0 =/ю, Р = 0), и условия на твердых стенках (/0 =Р0 =0) Основной сложностью решения краевых задач является нахождение фронтов УВ, положение которых заранее неизвестно.
Кроме условий на фронтах УВ для задач рефракции УВ характерным является формулировка условий на границе области возмущений с волной разрежения. Решение уравнения (1) для центрированных волн разрежения, возникающих, например, при рефракции на свободной поверхности в т А А , Пл) имеет вид
ДЬтО^-^УМ + ^о . z = 0l-тЫ/Й-Ы >
=rcos8, Ti-ii^ =rsm8 , z = tgQ , F(z) = + n^'g9 + Ф(в)/cos9 , Ф(9) = Рcosа(0 + С) , (3)
Р = ±Л/£>/(Л0-1) , а = tJ(R0 -1)/R0 Постоянные С, D в решении (3) находятся обычно при задании условий на переднем фронте центрированной волны разрежения На границе области возмущения решение (3) записывается на звуковой линии
г
Решение (3) через параметр R0 - Д0(у) зависит от газосодержания смеси
На свободной поверхности F(r,8) = 0, разделяющей среды с различным газосодержаниями у+, у~ имеют место условия (значки ± соответствуют различным областям)
~fo F
г
Условия (2)-(4) позволяют сформулировать краевые задачи для уравнения (1) для областей неоднородных потоков (областей возмущений) в задачах рефракции УВ
В п. 1.7 для решения задачи рефракции УВ в идеальных средах рассмотрено применение метода сращиваемых асимптотических разложений по малому параметру Р10 =(рх-.РоУРосо интенсивность УВ (За параметры
р0,с0 берутся их значения при у = у~) Этот метод позволяет редуцировать исходную задачу - свести решение задачи для области возмущений в целом к решению ряда задач для области (I) малых градиентов, где течение описывается с помощью линейной теории и находится аналитически, решениям задач для областей (III) значительных градиентов, примыкающих к фронтам УВ и линии слабого разрыва, где течение описывается системой нелинейных квазиодномерных уравнений, которая интегрируется аналитически, решениям в области (II) двумерных нелинейных течений со значительными градиентами параметров в окрестности точек взаимодействия УВ, где соответствующая краевая задача решается численно или аналитически приближенно
Метод позволяет существенно сузить область исследований до малой окрестности точек взаимодействия УВ - областей коротких волн Анализ общей постановки задач рефракции приводит (п 1 8) для областей (II) нелинейных взаимодействий во внутренних переменных 5, Y (X ,Y)
r= — = l + Pl0R0S, в = P^2Rq2Y, 5 = X + 1/2У2 ;
c0t
— = fr = floM. - = i/e = «/2v, p" = P^ c0 c0 r
(5)
к постановке краевой задачи, состоящей в нахождении решения ft, v системы уравнений коротких волн
2(ц-5)ц6+уу+ц = 0, jiY=v8 , ц = Р*0)=Я« , (6)
удовлетворяющих на границах области нелинейных взаимодействий условиям на фронтах УВ 8 = 5*(Y) (q„ - const)
+ Gi-H^v.-v, ц-J"«-*» (7)
Hi-i„. ^=-<^±0, av = иТу2м - №
"[0 "о Ш
асимптотическим условиям сращивания на границе с областью (I) линейного решения (5 —> —оо , -аойУ^О)
Р{1) = Ц = --{arctg^f^-arctg^A , (9)
я V-26 V-25
асимптотическим условиям сращивания на границах, с областями (III) нелинейных квазиодномерных решений (У —> ±оо, -оо ^5 ^8*) за фронтами УВ 5 = 5* (У) и линией слабого разрыва 8* = q„
д+ 1-д
5-5/
К\У) = 4
п
У2к\у)'
^ . . ^ ( I „ л (10)
Ч<ху-У ау+У
К2 (У)
Здесь в общем случае рефракции УВ = д~ при У -> -а>, - при У -><», знак «+» соответствует общему случаю, когда отраженные волны ударные, знак «-» соответствует вырождению волн в линию слабого разрыва
Кроме условий на фронтах УВ для задач рефракции характерны условия на границах области (II) с центрированной волной разрежения, описываемой решением (б)
(11)
2 У-УА л У-УА
\>(8, +^К2У-8аУ + с* , ¿-согш .
Условие на границе области возмущения с центрированной волной разрежения получаются при записи решения (11) вдоль звуковой линии (ц = 8 )
8-6^-1/2(У-УА)2 .
На свободной поверхности, разделяющей среды с различными газосодержаниями у+ , у~ имеют место условия
р+ = р~ , Л0+1/V = , 8* = = Л0-Р102у" . (12)
При решении краевых задач (5)-(12) основными трудностями являются нелинейность системы уравнений (6) и неизвестное положение ударных фронтов, удовлетворяющих дифференциальным уравнениям (7) Сложность проблемы усугубляется тем, что в условиях дозвукового характера течения (ц >5) параметры течения за образующимися ударно-волновыми структурами и сами структуры зависят от характеристик течения во всей области возмущения. Это обстоятельство побуждает при разработке схем расчета использовать решения, характеризующие поток за фронтами волн в целом Так, например, численный метод решения (см гл 3) позволяет построить итерационный процесс, состоящий в последовательном уточнении решения задачи и положения ударных фронтов. Результаты численного решения, .как показывает анализ, адекватно описывают процессы нелинейных взаимодействий УВ.
Приведенные исследования показывают, что для разработки теории нелинейных взаимодействий УВ, решения задач рефракции, необходимо развитие аналитических методов анализа краевых задач, развитие аналитических моделей рефракции УВ Эти исследования должны выявить основные ударно-волновые структуры (режимы рефракции) и области их
существования, позволить получить простые аналитические выражения для основных параметров, приближенно рассчитать поля течений в областях, прилегающих к ударно-волновым фронтам, при различных режимах нелинейной рефракции УВ
В главе 2 проведен асимптотический анализ режимов нелинейной рефракции УВ, опирающийся на рассмотренные аналитические модели ЯЯ, N11, Ю/ Эти модели отражают структуру течения в окрестности точки взаимодействия А - точки пересечения падающей волны со свободной поверхностью Методика разработки аналитических моделей рефракции УВ включает, разработку общего подхода для относительно слабых УВ на базе общих уравнений, описывающих потенциальные адиабатические течения газа, асимптотический анализ основных уравнений динамики, характеризующих соответствие, как продольных, так и поперечных вдоль свободной поверхности движений, позволяющих описать конкретные режимы рефракции, получение аналитических выражений для основных параметров, характеризующих течение и возникающие ударно-волновые структура; анализ областей существования различных режимов нелинейной рефракции УВ, общий анализ физической адекватности модели
В п. 2.1 развит общий подход для относительно слабых ударных волн (Р10 «1) к анализу режимов с волной разрежения (ИЛ, N11), опирающейся на использование общих аналитических выражений для фронтов УВ, волн разрежения и свободной поверхности (п 16 3) для потенциальных адиабатических течений для расчета параметров Ы1, те определим при известных е10 = (Р\-Ро)/Ро, а, У+, у~ > а также р0,Г0 покоящейся среды, интенсивности разрежения Е}0 ={р2 -Ро)/Ро » положения фронта преломленной УВ, переднего и заднего фронта волны разрежения, положения свободной поверхности, характеризуемых углами <э, (},, , 8 решается совместная система алгебраических уравнений Определены границы области существования типа А — перехода к режиму N11, типа В - обращения ш в О, типа С - соответствующей минимально возможному разрежению е20 (при
максимальных у") Исследования показали, что в случае раздела сред газ/ГЖС в широком диапазоне интенсивностей Ею и газосодержаний у" интенсивность преломленной УВ (е30 = е20) и изменение давления за ней (в верхней области Б) оказываются незначительными по сравнению с интенсивностью е10 и изменением давления в нижней области О Это обстоятельство позволяет сводить решение задач рефракции в этом случае к решению задач отражения УВ от свободной поверхности и не рассматривать течение в верхней области
В п. 2.2 проведен асимптотический анализ постановки задач рефракции Для относительно слабых ударных волн е«1, ё = £0(У~) ^ю уравнение выражающее равенство скоростей распространения падающей и преломленной УВ вдоль свободной поверхности «ют = сова /
(N+, N~- скорости распространения УВ) в первом приближении может быть представлено в виде
2(coso)-co -cosa-Сд )£ = (q+ cosaco -cos<DCgZ,)E ,
T-L0 r± - PoRQ n± _ „±„±2 „+ _P+ ~Po
. ¿o rj » Bo Poco . co =c0(Y ), q =-. (13)
¿5 fio A-Po
Ро Ро
Уравнение, выражающее соответствие поперечных к свободной поверхности движений в т. взаимодействия А за фронтом преломленной волны и за волной разрежения, в первом приближении может быть представлено в виде при
Р+=Р'
RqW = ±[2{XA-q+)?n +d , R =
ч 1/2
До+! V _ igra
mv=^7I, d-const (14)
е-
Исключая шиз (13) и (14) можно записать уравнение для определения величины q+, характеризующей интенсивность преломленной УВ, а также степень разрежения в случае рефракции с волной разрежения (<Г =(Р~-Ра)К,Р\-Ро))
(1-<у02г1
[U2(XA-q+)fn+df =(9--y + av-+2Cr)-C/Ë (15) R2q+ J L
В общем случае параметры рефракции УВ согласно (15), (14) определяются в пространстве параметров подобия
_« ! <\
В случае относительно малых углов a ~0(б ) в переменных теории коротких волн (5) выражение (13), (15) имеют вид
2ст =(шу2 ~av2) + ~=- — q~ ,
11 2 (1?) Ч -ZTT-Я+))3'2 + <Ча -a + .
R q+ 3 L
При анализе различных режимов рефракции с помощью (13)-(15) и (17) определяются конкретные значения ХА, d,q~.
В о. 2.3 — п. 2.6 рассмотрены различные режимы нелинейной рефракции УВ. регулярной (RR), нерегулярной (NR); регулярный с ударной волной, замыкающей область разрежения (RRV) и регулярный с отраженной ударной волной (RV) Определены соответствующие значения величин ХА, d, q~,
+
входящих в выражения (13)-(15), (17) позволяющие рассчитать основные параметры рефракции
Определены границы областей существования различных режимов в случае а ~0(б1/2)
NR Oiav 0, RR 10Sav^209, RRV 209^av В п. 2.7 исследован переход от режима регулярной рефракции с волной разрежения к режиму рефракции с отраженной ударной волной В общем случае, граница перехода в пространстве параметров подобия (16) имеет вид
(18)
1 R L
Граница (18)-переход, как для случаев су < 0 (для раздела сред типа slow-fast), так и для случаев су > 0 (для раздела сред типа fast-slow).
В п. 2.8 проведен анализ поведения характерной интенсивности д+. характеризующей интенсивность преломленной ударной волны, в пространстве параметров подобия av, , L, R (при R= 10) для различных режимов в случае рефракции на поверхности разделяющей газ/ГЖС
Рис 1 Зависимость интенсивности преломленной волны д* в случае рефракции на поверхности газ/ГЖС от параметров подобия а" ,сг,Ь при й = 1 О
На рис 1 поверхность A (av =1 0) соответствует границе перехода от режима NR к режиму RR Поверхность A' (av =2 09) соответствует границе перехода от режима RR к режиму RRV Поверхность К соответствует режиму рефракции с отраженной ударной волной для взаимодействия типа slow-fast (су <0) Кривая G', построенная согласно (18), есть граница перехода от
режима с волной разрежения к режиму с отраженной ударной волной Границы
С соответствуют предельно допустимым значениям Ь для рассматриваемых сред газ/ГЖС.
В п. 2.9 рассмотрены приложения результатов теории подобия к анализу физических параметров рефракции Для случая газ/ГЖС (Л =1.0) исследовано поведение относительной интенсивности и угла преломленной УВ со от угла падения инициирующей УВ а и газосодержания среды у при фиксированном е10 =0.1 Отмечено, что границ^ областей существования режимов регулярной и нерегулярной рефракции (а¥=10) построенной по формулам (17) при а~ согласуется с соответствующей границей
построенной с использованием общей теории п. 2.1.
В п. 2.10 для анализа физической адекватности регулярной рефракции произведены расчеты для случая двух газовых сред (С02 / С#4)
Рис 2 Сравнение теоретических результатов настоящих исследований с известными экспериментальными данными
Результаты проведенных теоретических исследований для модели ЯУ (с отраженной УВ), (рис 2), достаточно хорошо согласуются во всем диапазоне изменений углов а с известными экспериментальными результатами АЬ<1-Е1-РайаЬ (1978). Приведенные на рис. 2 результаты подтверждают наличие границы В (ш = 0), которая возникает в случае газ/газ, аналогично случаю газ/ГЖС отмеченному в п 2 1.
В целом, приведенные вп 28 — 210 результаты анализа аналитических исследований позволяют считать, что развитые аналитические модели
физически адекватно описывают процессы нелинейной рефракции для относительно слабых УВ (1 0 > ^ ^ 0 7), объясняют закономерности и особенности явления
В главе 3 развита методика численного решения краевых задач нелинейных взаимодействий относительно слабых УВ в применении к классу задач рефракции УВ для областей И.
В п. 3.1 рассмотрен метод последовательных приближений для построения решений краевых задач нелинейной рефракции УВ Для решения краевых задач в области II, сведенных в п 1 9 к первой краевой задаче с неизвестными положениями ударных фронтов, метод последовательных приближений состоит в задании начальных фронтов УВ и начального поля давления и последовательного уточнения положения ударных фронтов при решении соответствующей серии краевых задач.
В п. 3.2.1 рассмотрено применение метода прямых для построения начального приближения поля давления в различных задачах рефракции.
В п. 3.2.2 построены начальные приближения положения фронтов УВ для различных задач рефракции
В п. 3.3 для решения краевой задачи в области (II) использован неявный численный метод второго порядка точности, позволяющий при задании приближенного положения фронтов УВ п 3 2 2 и начального поля давления п 3 2 1 обеспечить достаточно быструю сходимость итераций для получения высокой точности
В п. 3.4 приведены результаты расчета полей давления для различных режимов рефракции УВ
В п. 3.4.1 приведены результаты расчета полей давления для различных режимов рефракции регулярного (av=1.0), нерегулярного (av=0.5) на поверхности разделяющей газожидкостные среды с близкими значениями газосодержания (у+ « у-) и близкими плотностями р£, Ро (R = 1 0, су = 0)
Результаты расчета показывают, что в области II происходят резкие изменения давления, как вдоль свободной поверхности, так и в направлении ортогональном ей Сравнение приведенных результатов численного расчета при относительно малых углах а~0(ё|/2), для относительно слабых УВ, с результатами численных исследований Henderson, Collela (1991) при конечных углах а ~0(1), для относительно слабых УВ (4 = 0 78), для случая газ/газ (С02 / С#4), для аналогичных режимов, показывает их качественное соответствие
В п. 3.4.2 приведены результаты расчетов полей давления для случая нерегулярного отражения от свободной поверхности соответствующих случаю газ/ГЖС (av=05) Для сравнения приведены результаты расчета поля давления для соответствующего случая при тех же значениях с помощью точных частных решений уравнений коротких волн
Сравнение показывает, что результаты расчета достаточно хорошо согласуются с аналитическими результатами вблизи фронта Маха и границы с
волной разрежения. Это обстоятельство позволяет считать, что результаты численного исследования получены с достаточной точностью.
В п. 3.4.3 приведены результаты расчета поля давления для случая регулярного отражения с возникновением УВ. замыкающей область разрежения, для случая газ/ГЖС (с^ =2 124) Результаты расчета показывают, что в области II происходят, аналогично случаю нерегулярного отражения, резкие изменения параметров
В целом, следует отметить, что развитый численный метод построения полей давления физически адекватно описывает процессы нелинейной рефракции для относительно слабых ударных волн , иллюстрирует возникновение различных режимов рефракции УВ, объясняет закономерности и особенности явления
В заключении диссертации перечислены основные результаты и общие выводы полученные в работе
В работе получены новые научные данные об ударно-волновых процессах нелинейной рефракции относительно слабых УВ в газах и газожидкостных средах
Развит аналитический подход, позволяющий получить количественные и качественные характеристики исследуемых нелинейных процессов, определить области существования параметров для различных режимов рефракции, когда интенсивность преломленной УВ сравнима с интенсивностью падающей УВ.
Часть работ по теме диссертации написана совместно с руководителем Шиндяпиным Г.П, участие которого в проведенных исследованиях состояло в постановке задач Основные результаты в совместных работах принадлежат лично автору.
Автор выражает искреннюю благодарность научному руководителю -профессору Г П Шиндяпину за поставленные задачи и большое внимание к работе.
ОСНОВНЫЕ РЕЗУЛЬТАТЫ И ВЫВОДЫ
1 Дан анализ современного состояния исследования по проблеме нелинейной рефракции ударных волн Развита аналитическая теория нелинейной рефракции ударных волн, использующая асимптотическую теорию коротких волн и модель равновесной газожидкостной среды.
2 Разработаны аналитические модели- нерегулярной, регулярной рефракции с волной разрежения, регулярной рефракции с ударной волной, замыкающей волну разрежения, с отраженной ударной волной, позволяющие исследовать ударно-волновые структуры течений различных режимов рефракции
3 Исследованы основные параметры и границы областей существования решений в пространстве безразмерных параметров подобия для различных режимов нелинейной рефракции ударных волн.
4 Развита методика численного решения краевых задач нелинейной рефракции относительно слабых ударных волн в газах и газожидкостных средах Построены картины течений для различных режимов рефракции.
нерегулярного, регулярного с волной разрежения, регулярного с ударной волной, замыкающей волну разрежения
5 Показана математическая и физическая адекватность аналитических и численных результатов путем их анализа и сравнения с известными экспериментальными данными, позволяющая рекомендовать полученные результаты для практического использования при решении проблем современного естествознания и техники
В целом приведенные результаты представляют совокупность завершенных исследований, развивающих перспективное направление — нелинейную теорию
взаимодействий ударных волн в газах и газожидкостных средах
!
ПУБЛИКАЦИИ ПО ТЕМЕ ДИССЕРТАЦИИ
1 Матутин А А Анализ режимов рефракции ударной волны с образованием волны разрежения или отраженной ударной волны / ГП Шиндяпин, А А Матутин // Известия СГУ Сб научных трудов Саратов Изд-во Сарат ун-та, 2007 Том 7 Вып 1 С 74-77 (данный журнал входит в перечень изданий ВАК рекомендованный для защит кандидатских диссертаций) (Матутин А А исследовал различные режимы рефракции в пространстве параметров подобия)
2 Матутин А А Особенности нелинейной рефракции ударной волны / Г П Шиндяпин, А А Матутин // Вестник Сарат roc технич ун-та Науч -технич журнал Саратов Изд-во Сарат гос технич ун-та, 2006 Вып 1 № 3(14) С 44-48 (Матутин А А исследовал зависимость интенсивности преломленной ударной волны от интенсивности падающей ударной волны для случая газ/ГЖС)
3 Матутин АЛ Численное исследование нелинейной рефракции ударной волны на поверхности, разделяющей газовую и газожидкостную среды / Г П Шиндяпин, А А Матутин // Известия СГУ Сб науч тр Саратов Изд-во Сарат ун-та, 2006 Том 6 Вып 1/2 С 67-74 (Матутин А А построил поля давления для режима рефракции с ударной волной, замыкающей область разрежения)
4 Матутин А А Исследование нелинейной рефракции ударной волны на поверхности, разделяющей газовую и газожидкостную среды численным методом / Г П Шиндяпин, А А Матутин // Математические методы и модели в прикладных задачах науки Сб науч тр Ульяновск Изд-во Ульяновского гос технич ун-та, 2005 Том 7 С 112-117 (Матутин А А построил поля давления для случая рефракции газ/ГЖС)
5 Матутин А А Аналитическое исследование класса задач рефракции с невырожденным фронтом преломленной ударной волны на поверхности газ/газожидкостная среда / Г П Шиндяпин, А А Матутин // Сб науч тр «Прикладная математика и механика» Ульяновск Изд-во Ульяновского гос технич ун-та, 2006 Вып 8 С 56-60 (Матутин А А рассчитал поведение основных параметров рефракции для случая с невырожденным фронтом преломленной ударной волны)
6 Матутин А А Анализ ударно-волновых структур и потоков при нелинейной рефракции ударных волн / Г П Шиндяпин, А А Матутин // Математика Механика Сб науч тр Саратов Изд-во Сарат ун-та, 2003 Вып 5 С 162-165 (Матутин А А построил с использованием точных решений поля давления для нерегулярной рефракции)
7 Матутин А А Аналитическое исследование нелинейной рефракции ударной волны на поверхности, разделяющей газовую и газо жидкостную среды / Г П Шиндяпин, А А Матутин // Механика и процессы управления Труды XXXIV Уральского семинара Екатеринбург УРО РАН, 2004 Том. 1 С 190-198 (Матутин А А рассчитал параметры рефракции на поверхности газ/ГЖС в соответствии с общей теорией постановки задач регулярной рефракции)
8 Матутин Л А Аналитическое исследование ударно-волновых структур при нелинейной рефракции ударной волны на поверхности, разделяющей газовую и газожидкостную среды / Г П Шиндяпин, А А Матутин // Тр XX Юбилейного Междунар семинара по струйным, отрывным и нестационарным течениям Санкт-Петербург, 1-4 июля 2004 г Санкт-Петербург ИПЦ СПб ГУТД, 2004 С 251-252 (Матутин А А описал основные результаты расчетов параметров рефракции в соответствии с общей теорией постановки задач рефракции)
9 Матутин А А О математических моделях и режимах рефракции ударных волн в газожидкостных пузырьковых средах / ГП Шиндяпин, А.А Матутин // Тр XXI Всероссийского семинара по струйным, отрывным и нестационарным течениям. Новосибирск, 15-18 августа 2007 г Новосибирск Параллель, 2007 С 199-201 (Матутин А А исследовал области существования и границы различных режимов рефракции в пространстве параметров подобия)
10 Матутин А А К теории нелинейной рефракции ударной волны / Г П Шиндяпин, А.А Матутин // Математика Механика Сб науч тр Саратов Изд-во Сарат ун-та, 2006 Вып 8 С 199-203 (Матутин А А рассчитал поведение параметров регулярной рефракции в соответствии с развитой асимптотической теорией)
11 Матутин А А О математических моделях и асимптотических законах подобия рефракции ударных волн в газожидкостных пузырьковых средах / Г П Шиндяпин, А А Матутин // Математические методы и модели в прикладных задачах науки Сб науч тр Ульяновск Изд-во Ульяновского гос. технич ун-та, 2006 Том 8 С 106-110 (Матутин АЛ провел расчеты параметров рефракции иллюстрирующих законы подобия процессов рефракции ударных волн)
12 Матутин А А Аналитический и численный анализ нелинейной рефракции ударной вслны / А А Матутин // Труды математического центра НИ Лобачевского Т 31 Казанское математическое общество «Лобачевские чтения — 2005» Материалы IV молодежной научной школы - конференции Казань Изд-во Казанского мат общ
2005 С 99-100
13 Матутин А А Аналитический и численный анализ нелинейной рефракции ударной волны / А А Матутин // Труды математического центра Н И Лобачевского Т. 34 Казанское математическое общество «Лобачевские чтения - 2006» Материалы V молодежной научной школы - конференции Казань Изд-во Казанского мат общ
2006 С 162-163.
14 Матутин А А Аналитическое и численное исследование процессов рефракции ударных волн / А А Матутин // Математика Механика Сб науч тр Саратов Изд-во Capar ун-та, 2005 Вып 7 С 183-186
15 Матутин А А Нелинейная рефракция ударной волны с образование ударной волны, замыкающей зону разрежения / А А Матутин // Математика Механика Сб науч тр Саратов Изд-во Сарат ун-та, 2004. Вып 6 С 198-202
Подписано к печати 04 10 07г Формат 60x84 1/16 Бумага офсетная Гарнитура Tunes New Roman Печать RISO Объем-1,2 печ л Тираж 100 экз Заказ № 139-Т
Отпечатано с готового оригинал-макета Типография Саратовского государственного университета имени Н Г Чернышевского 410012 г Саратов, ул Большая Казачья, д 112 а
Введение.
Глава 1. Математическое моделирование процессов рефракции ударных волн в газах и газожидкостных средах.
1.1. Класс задач нелинейной рефракции ударных волн.
1.2. Математическая модель газожидкостной среды.
1.3. Современное состояние исследований проблем нелинейной рефракции ударных волн.
1.4. Проблемы математического моделирования нелинейной рефракции ударных волн.
1.5. Анализ применимости модели двухфазной ГЖС для описания ударно-волновых процессов.
1.6. Постановка краевых задач для описания взаимодействий ударных волн в газовых и газожидкостных средах.
1.6.1. Дифференциальные уравнения идеальной, сжимаемой газовой и газожидкостной среды.
1.6.2. Граничные условия для областей неоднородных течений при нелинейной рефракции УВ. Условия на фронтах УВ.
1.6.3. Постановка задач для потенциальных адиабатических течений газовых и газожидкостных сред.
1.7. Метод сращиваемых асимптотических разложений решения краевых задач взаимодействия ударных волн в идеальных средах.
1.7.1. Линейные решения для областей (I) малых градиентов.
1.7.2. Нелинейные решения в областях (III), прилегающих к фронтам ударных волн и линий слабого разрыва.
1.8. Краевые задачи для областей (II) в окрестности точек взаимодействия ударных волн, областей нелинейных взаимодействий.
1.8.1. Система уравнений коротких волн.
1.8.2. Граничные условия для областей коротких волн.
1.8.3. Решение в области разрежения.
1.8.4. Координаты тройной точки.
1.8.5. Преобразование нелинейной краевой задачи к первой краевой задаче
Выводы к главе 1.
Глава 2. Асимптотический анализ режимов нелинейной рефракции ударных волн.
2.1. Аналитическое исследование режимов рефракции ударных волн. Общий подход для относительно слабых ударных волн.
2.2. Асимптотический анализ постановки задач рефракции.
2.3. Регулярная рефракция с возникновением волны разрежения (RR).
2.4. Случай нерегулярной рефракции с волной разрежения (NR).
2.5. Случай регулярной рефракции с ударной волной (RRV), замыкающей область разрежения.
2.6. Случай рефракции с отраженной ударной волной.
2.7. Переход от рефракции с волной разрежения (RR, RRV) к рефракции с отраженной ударной волной (RV).
2.8. Анализ характерной интенсивности q+ преломленной ударной волны в пространстве параметров подобия.
Актуальность темы. Теоретическое исследование процессов нелинейной рефракции ударных волн (УВ) в газах и газожидкостных средах представляет одну из фундаментальных проблем современной механики жидкости и газа. Решение этой проблемы важно для понимания многих физических процессов; развития сверхзвуковой авиации; проектирования трубопроводов, для транспортировки топливных смесей; движения тел в газожидкостных средах, в которых скорость звука мала и др. Теоретическая важность проблемы обусловлена нелинейным характером основных уравнений и сложным разрывным характером решений, в силу чего методы и их применения при исследовании задач имеют универсальный характер и связаны с решением общих проблем выявления структуры обобщенных решений задач математической физики.
В работе исследуются процессы нелинейной рефракции ударной волны (УВ), возникающие при взаимодействии ударной волны, относительно малой интенсивности (абсолютная интенсивность волн может быть велика), со свободной поверхностью, разделяющей различные газожидкостные среды, когда за фронтами УВ возникают области резких изменений параметров -области коротких волн, и ударные нагрузки многократно возрастают. В этих областях процесс существенно нелинеен, что вызывает значительные трудности при анализе поставленной задачи. Особое внимание уделяется нахождению и исследованию различных режимов существования нелинейной рефракции УВ, а также численному анализу поставленной задачи.
Пристальное внимание к проблеме отечественных (С.А. Христианович, А.А. Гриб, Б.И. Заславский, Г.П. Шиндяпин, В.К. Кедринский и др.) и зарубежных исследователей (L.F. Henderson, A. Sakurai, К. Takayama, G. Ben-Dor и др.) на протяжении более 50 лет привело к появлению разнообразных альтернативных точек зрения на природу явления. Однако, в настоящее время все отчетливее вырисовывается необходимость исследования процессов рефракции в целом, для построения полной картины исследуемого явления.
Настоящая работа развивает и дополняет положения асимптотической теории коротких волн [47, 81], опирающейся на решение внутренних краевых задач для областей нелинейных взаимодействий (коротких волн), которые позволяют учесть влияние потока в целом на образующиеся ударно-волновые структуры, изучить характерные особенности, закономерности и особые точки процессов нелинейной рефракции УВ.
Цель и задачи исследования. Целью работы является аналитическое исследование различных режимов существования нелинейной рефракции относительно слабых УВ в газах и газожидкостных (пузырьковых) средах, разработка и применение численного метода решения краевых задач при исследовании нелинейной рефракции УВ.
Исходя из этого, в работе решались следующие задачи:
• разработка аналитических моделей нелинейной рефракции относительно слабых УВ;
• классификация режимов нелинейной рефракции УВ;
• анализ областей существования различных режимов нелинейной рефракции УВ;
• получение явных аналитических зависимостей для основных параметров, характеризующих возникающие ударно-волновые структуры;
• аналитическое описание границ областей существования различных режимов нелинейной рефракции УВ;
• численное решение краевых задач рефракции УВ, на основе построенных аналитических моделей;
• анализ полей давлений и скоростей при различных режимах нелинейной рефракции УВ.
Методика исследований основана на использовании асимптотической теории коротких волн при анализе режимов нелинейной рефракции УВ, постановке краевых задач для областей нелинейных взаимодействий (коротких волн), разработке и применении численного метода решения краевых задач подобного типа.
Основное внимание уделяется аналитическим методам анализа, основанным на выделении структурных особенностей течений.
Научная новизна работы заключается в следующем:
• развит новый аналитический метод исследования ударно-волновых структур и потоков при нелинейной рефракции относительно слабых УВ в газах и газожидкостных средах;
• исследованы различные режимы нелинейной рефракции УВ, найдены аналитические выражения для границ областей существования различных режимов нелинейной рефракции УВ;
• найдены области существования неклассических режимов рефракции УВ, которые расширяют физическое представление о процессах нелинейной рефракции УВ, позволяя исследовать случаи рефракции с невырожденным фронтом преломленной УВ;
• получил развитие численный метод решения краевых задач нелинейной рефракции УВ, позволяющий рассчитать поля давлений и скоростей, а также, структуры ударно-волновых процессов нелинейной рефракции УВ;
• исследованы ударно-волновые структуры и течения в областях нелинейных взаимодействий. Выявлены качественные закономерности и особенности процессов нелинейной рефракции УВ, показавшие физическую адекватность разработанных аналитических и численных подходов к решению поставленной задачи.
Достоверность результатов подтверждается непротиворечивостью полученных аналитических результатов с имеющимися физическими представлениями и экспериментальными данными; согласованностью результатов для рассчитанных ударно-волновых структур и полей давления и скоростей с известными результатами численного решения соответствующих краевых задач; обоснованностью используемых методов исследований.
Практическая ценность работы. Результаты исследований расширяют представления о физических процессах нелинейной рефракции УВ в газах и газожидкостных средах, развивают аналитические методы и методы численного расчета течений для относительно слабых УВ. Результаты исследований могут быть использованы в учебном процессе при подготовке учебных пособий и чтении специальных курсов по динамике ударных волн в газах и газожидкостных средах. k Апробация результатов и публикации. Результаты исследований докладывались и обсуждались на международных, всероссийских и вузовских конференциях и семинарах. Среди них:
• XX, XXI Международный семинар по струйным, отрывным и нестационарным течениям. Санкт-Петербург, 2004 г.; Новосибирск, 2007 г.;
• XXXIV Уральский семинар по механике и процессам управления. Миасс, 2004 г.;
• XXIV, XXVII Российская школа по проблемам науки и технологий. Миасс, 2004; 2007 гг.;
• XLII Международная научная конференция «Студент и научно-технический прогресс». Новосибирск, 2005 г.;
• IV, V Молодежные научные школы-конференции «Лобачевские чтения». Казань, 2005; 2006 гг.;
• Международная конференция «Континуальные алгебраические логики, исчисления и нейроинформатика в науке, технике и экономике». Ульяновск, 2005; 2006 гг.;
• Ежегодные научные конференции Саратовского государственного университета «Актуальные проблемы математики и механики». Саратов, СГУ: 2003; 2004; 2005; 2006; 2007 гг.;
I • Специальные семинары кафедры вычислительного эксперимента в механике
Саратовского государственного университета. 2003-2007 гг.
Результаты исследований опубликованы в пятнадцати печатных работах.
Структура и объем работы.
Диссертация состоит из введения, трех глав, выводов по главам, заключения, списка использованной литературы, изложена с учетом рисунков k на 123 страницах. Работа содержит 19 рисунков. Список использованной
Выводы к главе III.
Получил развитие численный метод решения краевых задач нелинейной рефракции ударных волн. Разработана общая методика численного решения краевых задач нелинейной рефракции ударных волн, позволяющая построить поля давлений и скоростей в области нелинейных взаимодействий ударных волн. Построены картины течений для различных режимов рефракции ударных волн: нерегулярного, регулярного и регулярного с ударной волной, замыкающей волну разрешения.
Показано, что развитый численный метод решения краевых задач взаимодействия оказался весьма эффективным и экономичным при решении задач рефракции, позволяя получить результаты (поля давлений и положения фронтов УВ) с достаточной точностью при ограниченном числе итераций и серий краевых задач.
Результаты объясняют характерные особенности и свойства течений при различных режимах нелинейной рефракции ударных волн. Результаты численных исследований качественно и количественно подтверждаются известными результатами экспериментальных и аналитических исследований.
Заключение
В работе получены новые научные данные об ударно-волновых процессах нелинейной рефракции относительно слабых УВ в газах и газожидкостных средах. Развит аналитический подход, позволяющий получить количественные и качественные характеристики исследуемых нелинейных процессов; определить области существования параметров для различных режимов рефракции, когда интенсивность преломленной УВ сравнима с интенсивностью падающей УВ. Развита методика численного решения краевых задач нелинейной рефракции УВ, построены картины течений для различных режимов рефракции.
Основные научные результаты:
1. Развита аналитическая теория нелинейной рефракции ударных волн, использующая асимптотическую теорию коротких волн и модель равновесной газожидкостной среды.
2. Разработаны аналитические модели: нерегулярной; регулярной рефракции с волной разрежения; регулярной рефракции с ударной волной, замыкающей волну разрежения; с отраженной ударной волной, позволяющие исследовать ударно-волновые структуры течений различных режимов рефракции.
3. Проведен анализ основных параметров и границ областей существования решений в пространстве безразмерных параметров подобия для различных режимов нелинейной рефракции ударных волн.
4. Развита методика численного решения краевых задач нелинейной рефракции относительно слабых ударных волн в газах и газожидкостных средах. Построены картины течений для различных режимов рефракции: нерегулярного; регулярного с волной разрежения; регулярного с ударной волной, замыкающей волну разрежения.
5. Показана математическая и физическая адекватность аналитических и численных результатов путем их анализа и сравнения с известными экспериментальными и численными результатами исследований, позволяющая рекомендовать полученные результаты для практического использования при решении проблем современного естествознания и техники.
В целом приведенные результаты представляют совокупность завершенных исследований, развивающих перспективное направление - нелинейную теорию взаимодействий ударных волн в газах и газожидкостных средах.
1. Abd-El-Fattah A.M. Shock waves at a slow-fast gas interface. / Abd-El-Fattah A.M., Henderson L.F. // J. Fluid Mech. Vol. 89. 1978. part 1, pp. 79-85.
2. Adachi T. Mach reflection of a weak shock waves / Adachi Т., Suzuki Т., Kobayashi S. // Trans. Jap. Soc. Mech. Eng. B. 1994. V. 60. № 575. P. 2281-2286.
3. Ben-Dor G. Shock wave reflection phenomena. New-York: Springer Verlag. 1991.
4. Bleakney W. Interaction of shock waves. / Bleakney W., Taub A.H. // Rev. Mod. Phys. 21, 584. 1949.
5. Bleakney W., Taub A.H., Fletcher C.H. Rev. Mod. Phys. 23, 271. 1951.
6. Campbell I.J. Shock waves in a liquids containing gas bubbles / Campbell I.J. Pitcher A.S. //Proc. Roy. Soc. A. 243. 1958. № 1235.
7. Cole R.H. Underwater Explosions Dover, New York, 1948 and references therein.
8. Chester W. Q. J. Mech. Appl. Math. 6,440. 1953.
9. Falada A., Holt M. Phys. Fluids. 21. 1702. 1978.
10. Flores J. Shock wave interaction with the ocean surface. / Flores J., Holt M. // Phys. Fluids. 25(2). 1982. p. 238.
11. Griffith W.C. Shock waves / Griffith W.C. // J. Fluid Mech. 1981. Vol. 106. P. 81101.
12. Griffith W.C., Bleakney W. Am. J. Phys. 22, 597. 1954.
13. Gvozdeva L.G., Kazhenova T.V., Predvoditeleva O.H., Fokeev V.P. Astron. Acta. 1970.15. 503.
14. Henderson L.F. On the refraction of shock waves at a slow-fast gas interface. / Henderson L.F., Colella P., Puckett E.G. // J. Fluid Mech. Vol. 224. 1991. pp. 127.
15. Henderson L.F. Reflection of a shock wave at on air-water interface / Henderson L.F., Ma J.H., Sakurai A., Takayama K. // Fluid Dynam. Res. 1990. № 5. P. 337350.
16. Henderson L.F. On the irregular refraction of a plane shock wave at a Mach number interface. / Henderson L.F., Macpherson A.K. // J. Fluid Mech. Vol. 32. 1968. part. l.pp. 185-202.
17. Henderson L.F. Experiments on the diffraction of weak blast waves: the von Neumann paradox / Henderson L.F., Siegenthaler A. // Proc. Roy. Soc. London. Ser. A. 1980. Vol. 369. № 1739. P. 537-555.
18. Henderson L.F. The persistence of regular reflection during strong shock diffraction over rigid ramps. / Henderson L.F., Takayama K., Crutchfield W.Y. Itabashi S. //J. Fluid Mech. 2001. Vol. 431. pp. 273-296.
19. Henderson L.F. The wall-jetting effect in Mach reflection: theoretical consideration and numerical investigation. / Henderson L.F., Vasilev E.I., Ben-Dor G., Elperin T. // J. Fluid Mech. 2003. Vol. 479. pp. 259-286.
20. Holt M. Annual Reviews of Fluid Mechanics. / Holt M. //Annuals Reviews, Palo Alto. 1977. Vol. 9, p. 187.
21. Jahn R.G. The refraction of shock waves at a gaseous interface. / Jahn R.G. // J. Fluid Mech. 1.457. 1956.
22. Kameda M. Propagation of Shock Waves in Dilute Bubbly Liquids. / Kameda M., Matsumoto Y. // Trans. Jap. Soc. Mech. Eng. B. 1995. 61. № 584. P. 1236-1245.
23. Lighthill M.J. The shock strength in supersonic "conical fields". / Lighthill M.J. // Philos. Mag. 1949. Vol. 40, № 511. P. 1202-1223.
24. Neumann J. Collected Works. Oxford; London; New York; Paris: Pergamon Press. 1963.538 р.
25. Polachek H., Seeger R.J. In Fundamentals of Gas Dynamics, edited by H.W. Emmons. Princeton University Press. Princeton. NJ, 1958. 482 p.
26. Quirk J.J. On the dynamics of a shock-bubble interaction. / Quirk J.J., Kami S. // NASA Langley Research Center, Hampton. 1995.
27. Sakurai A. U.S. Army Engineering Waterways / Sakurai A. //Experiment Station
28. Technical Report. № 1-771. 1967. 28.Sasoh A. Characterization of disturbance propagation in weak shock-wave reflections. / Sasoh A, Takayama K. // J. Fluid Mech. 1994. Vol. 277. P. 331-345.
29. Smith W.R. Mutual reflection of two shock waves of arbitrary strengths. / Smith W.R. // Phys. Fluid. 1959. Vol. 2. № 5. P. 533-541.
30. Taub A.H. Phys. Rev. 72, 51. 1947.
31. White D. R. An experimental survei of the Mach reflection of shock waves /White D. R. // Proc. Second Midwestern Fluid Mech. 1952. P. 253-262.
32. Whitham G.B. A new approach to problem of Shock dynamics. Pt 1: Two dimantional problems / Whitham G.B. // J. Fluid Mech. 1957. Vol. 2. P. 145-171.
33. Багдоев А.Г. Пространственные нестационарные движения сплошной среды с ударными волнами. Ереван. Изд-во АН АрмССР. 1961. 276 с.
34. Баженова Т.В., Гвоздева Л.Г., Лагутов Ю.П. Нестационарные взаимодействия ударных и детонационных волн. М.: Наука, 1986. 207 с.
35. Булах Б.М. Нелинейные конические течения газа. М.: Наука. 1970.344 с.
36. Ван-Дайк М. Методы возмущений в механике жидкости / Пер. с англ.М.: Мир. 1967.310 с.
37. Вельмисов П.А. Асимптотические исследования нелинейных взаимодействий слабых ударных волн / Вельмисов П.А., Шиндяпин Г.П. // Аэродинамика. Межвуз. сб. Саратов: Изд-во Сарат. ун-та, 1972. Вып. 1(4). С. 78-93.
38. Волошинов А.В. Влияние газосодержания при взаимодействии ударных волн различной интенсивности в двухфазной газожидкостной среде / Волошинов А.В., Ковалев А.Д., Шиндяпин Г.П. // Изв. АНССР. МЖГ. 1986. №6. С. 172-174.
39. Волошинов А.В. К теории регулярного и маховского отражения ударных волн в двухфазной газожидкостной среде / Волошинов А.В., Ковалев А.Д., Шиндяпин Г.П. // Изв. АНСССР. МЖГ. 1985. № 2. С. 188-190.t
40. Гамаюнова E.H. Аналитическое исследование ударно-волновых структур и потоков при отражениях и взаимодействиях относительно слабых ударных волн в газах и газожидкостных средах. Дисс. на соиск. уч. степ. канд. физ.-мат. наук. 2004.
41. Гельфанд Б.Е. Преломление плоских ударных волн при взаимодействии со слоем пузырьки газа-жидкость. / Гельфанд Б.Е., Губанов А.В., Тимофеев Е.И. // Изв. АНСССР. МЖГ. 1981. №2. С. 173-176.
42. Гельфанд Б.Е. Расчет параметров нестационарных ударных волн в двухфазной среде / Гельфанд Б.Е., Губанов А.В., Тимофеев Е.И. // ФГВ. 1981. №5. С. 139-143.
43. Гриб А.А. Нерегулярное отражение плоской ударной волны в воде от свободной поверхности. / Гриб А.А., Березин А.Г. // ПМТФ. 1960. № 2. С. 34-39.
44. Гриб А.А. Теория коротких волн / Гриб А.А., Рыжов О.С., Христианович С.А. // ПМТФ. 1960. № 1. С. 63-74.
45. Гриб А.А. Об отражении плоской ударной волны от свободной поверхности. / Гриб А.А., Рябинин А.Г., Христианович С.А. // ПММ, 1956. Том. 20. Вып. 4.
46. Губайдулин А.А. Волны в жидкостях с пузырьками / Губайдулин А.А., Ивандаев А.И., Нигматулин Р.И., Хабеев И.С. // Итоги науки и техники. Механика жидкости и газа. ВИНИТИ М. 1982. Т. 17. С. 16-249.
47. Кедринский В.К. Усиление ударных волн при столкновении и фокусировке в пузырьковых средах / Кедринский В.К., Вшиков В.А., Дудникова Г.И., Шокин Ю.И. // ДАН, 1998. Т. 361, № 1. С. 41-44 .
48. Кедринский В.К. Взаимодействие и усиление волн в жидких средах с пузырьками свободного и химически активного газа / Кедринский В.К. // Информационный бюллетень РФФИ. Новосибирск. Изд-во РФФИ, 1998. Том 6. № 2. 773 с.
49. Кедринский В.К. Гидродинамика взрыва: эксперимент и модели-Новосибирск: Изд-во СО РАН. 2000,-435 с.
50. Клейнер Б.Г. Об одном классе точных частных решений уравнений коротких волн. / Клейнер Б.Г., Шиндяпин Г.П. // ПММ. Т. 34. Вып. 6. 1970. С. 11501158.
51. Ковалев А.Д. Исследование структур зон взаимодействия ударных волн. Дисс. на соиск. уч. степ, к.ф.м.н. Саратов: Сарат. ун-т. 1987.
52. Ковалев А.Д. О начальной стадии взаимодействия слабых ударных волн / Ковалев А.Д., Шиндяпин Г.П. // Аэродинамика: Межвуз. сб. Саратов: Изд-во Сарат. ун-та. 1979. Вып. 7(10). С. 116-122.
53. Ковалев А.Д. К теории гидроудара в двухфазной газожидкостной смеси. / Ковалев А.Д., Шиндяпин Г.П. // Изв. АНССР. МЖГ. 1985. № 5. С. 188-190.
54. Ковалев А.Д. Идентификация режимов течения газожидкостных смесей в горизонтальных трубопроводах / Ковалев А.Д., Терентьев С.А., Шиндяпин Г.П. // Сб. науч. тр. Исследование и измерение параметров криогенных потоков. ВНИИФТРИ. М: 1981. С. 3-12.
55. Ляхов В.И., Подлубный В.В., Титаренко В.В. Воздействие ударных волн и струй на элементы конструкций. М.: Машиностроение. 1989. 392 с.
56. Маркушин А.Г. О максимальном давлении и асимптотическом законе подобия при переходе от регулярного к маховскому взаимодействию ударных волн в газожидкостных средах /Маркушин А.Г., Шиндяпин Г.П. //
57. Аэродинамика: Межвуз. сб. Саратов: Изд-во Сарат. ун-та, 1993. Вып. 13(16). С.14-27.
58. Маркушин А.Г. Рефракция ударной волны на свободной поверхности в газожидкостной среде с образованием волны разрежения / Маркушин А.Г., Шиндяпин Г.П. // Аэродинамика: Межвуз. сб. Саратов: Изд-во Сарат. ун-та, 1991. Вып. 12(15). С. 24 -39.
59. Накоряков В.Е. Экспериментальное исследование ударных волн в жидкости с пузырьками газа / Накоряков В.Е., Покусаев Б.Г., Шрейбер И.Р., Кузнецов В.В., Малых Н.В. // Волновые процессы в двухфазных системах. Новосибирск, 1975. С. 54-97.
60. Паркин Б.Р. Ударные волны в воде с пузырьками воздуха / Паркин Б.Р., Гилмор Ф.Р., Броуд Г.Л. // Подводные и подземные взрывы. (Пер. с англ.) М.: Мир. 1974 С. 152-258.
61. Рязанцев Ю.С. Об отражении ударной волны от граничной поверхности. / Рязанцев Ю.С. //ПМТФ. 1962. № 2. С. 122-123.
62. Сагомонян А.Я. Пространственные задачи неустановившегося движения сжимаемой жидкости. М.: Изд-во МГУ. 1962. 80 с.
63. Сагомонян А.Я., Поручинов В.Б. Пространственные задачи неустановившегося движения сжимаемой среды. М.: Изд-во МГУ. 1970. 120 с.
64. Семенов А.Н. Разновидности отражения ударных волн от клина / Семенов А.Н. // Нестационарные течения газов с ударными волнами. Л.: ФТИ АН СССР. 1990. С. 162-172.
65. Титаренко В.В. Исследования взаимодействий ударных волн в идеальных газожидкостных средах асимптотическими методами / Титаренко В.В., Шиндяпин Г.П. // Нестационарные течения с ударными волнами. Л.: ФТИ АН СССР. 1990. С. 199-215.
66. Усков В.Н. Интерференция стационарных газодинамических разрывов / Усков В.Н. // Сверхзвуковые газовые струи. Новосибирск: Наука. 1983. С. 22-46.
67. Фонарев A.C. Распространение взрывных волн при частичном разрушении газонасыщенной цилиндрической оболочки / Фонарев А.С., Подлубный В.В. // Труды ЦАРИ. 1977. Вып. 1834. 23 с.
68. Христианович С.А. Ударные волны на значительном расстоянии от места взрыва. / Христианович С.А. // ПММ. 1956. Том. 20. Вып. 5. С. 599-605.
69. Христианович С.А. Механика сплошной среды. М.: Наука. 1981. 484 с.
70. Христианович С.А., Заславский Б.И. Асимптотические методы в теории слабых ударных волн. XXI Международный конгресс по прикладной механике. Н. Gortler, ed. Springer-Verlag, Berlin, 1966.
71. Христианович С.А. О нелинейном отражении слабых ударных волн. / Христианович С.А., Рыжов О.С. // ПММ. 1958. Т. 22. Вып. 5.
72. Шиндяпин Г.П. Аналитическое исследование ударно-волновых структур и потоков при отражении и взаимодействии относительно слабых ударных волн / Шиндяпин Г.П. // Аэродинамика: Межвуз. сб. науч. тр. Саратов: Изд-во Сарат. ун-та. 2001. Вып. 3. С. 31-44.
73. Шиндяпин Г.П. К теории дифракции слабых ударных волн / Шиндяпин Г.П. //Изв. АН СССР. МЖГ. 1984. № 6. С. 176-178.
74. Шиндяпин Г.П. Нелинейные взаимодействия ударных волн в газах и газожидкостных средах Учеб. пособие. Саратов: Изд-во Сарат. ун-та. 1997104 с.
75. Шиндяпин Г.П. Нерегулярное взаимодействие слабых ударных волн разной интенсивности / Шиндяпин Г.П. // ПММ 1974. Т. 38. Вып. 1. С. 105-114.
76. Шиндяпин Г.П. О режиме регулярного взаимодействия слабых ударных волн // ПММ. 1965. Т. 29. Вып. 6. С. 1128-1131.
77. Шиндяпин Г.П. О нерегулярном отражении слабых ударных волн от жесткой стенки / Шиндяпин Г.П. // ПМТФ. 1964. № 2. С. 22-28.
78. Шиндяпин Г.П. Исследование рефракции ударных волн на поверхности, разделяющей газовую и газожидкостную среды, асимптотическим методом / Шиндяпин Г.П. // Аэродинамика. Межвуз. сб. научн. тр Саратов: Изд-во Сарат. ун-та, 1995. вып. 14 (17).
79. Шиндяпин Г.П. Численное решение задачи нерегулярного отражения слабой ударной волны от жесткой стенки в идеальном газе / Шиндяпин Г.П. // ЖВМ иМФ. 1980. № 1.С. 249-254.
80. Шиндяпин Г.П., Ковалев А.Д. Математическое моделирование в задачах динамики многофазных сред. Учеб. пособие. Саратов. Изд-во Сарат. ун-та.1988.4.1.
81. Шиндяпин Г.П., Ковалев А.Д. Математическое моделирование в задачах динамики многофазных сред. Учеб. пособие. Саратов. Изд-во Сарат. ун-та.1990.4.2.
82. Матутин А.А. Аналитическое и численное исследование процессов рефракции ударных волн / Матутин А.А.// Математика. Механика: Сб. науч. тр. Саратов: Изд-во Сарат. ун-та, 2005. Вып. 7. С. 183-186.
83. Матутин А.А. Нелинейная рефракция ударной волны с образование ударной волны, замыкающей зону разрежения / Матутин А.А. // Математика. Механика: Сб. науч. тр. Саратов: Изд-во Сарат. ун-та, 2004. Вып. 6. С. 198202.