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

1Московский государственный университет имени М.В. Ломоносова

Механико-математический факультет

ша1|«и —

Тагирова Рената Рифовна Кумуляция массы плотных газовых

слоев при их ускорении в нормальном направлении

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

2 6 НОЯ 2009

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

Москва — 2009 ^

Работа выполнена на кафедре аэромеханики и газовой динамики Механико-математического факультета Московского государственного университета имени М.В. Ломоносова и в лаборатории теоретического моделирования Института космических исследований Российской академии наук

Научный руководитель:

Официальные оппоненты:

Ведущая организация:

доктор физико-математических наук, профессор К.В. Краснобаев

доктор физико-математических наук, профессор А.Н. Голубятников

кандидат физико-математических наук М. К. Ермаков

Московский физико-технический институт, факультет аэрофизики и космических исследований

Защита диссертации состоится 18 декабря 2009 г. в 15 часов на заседании диссертационного совета Д.501.001.89 при Московском государственном университете имени М.В. Ломоносова, по адресу: Российская Федерация, 119991, ГСП-1, Москва, Ленинские горы, Главное здание МГУ, Механико-математический факультет, аудитория 46-24

С диссертацией можно ознакомиться в Научной библиотеке Механико-математического факультета МГУ им. М.В. Ломоносова.

Автореферат разослан ноября 2009 г.

Ученый секретарь

диссертационного совета Д.501.001.89 доктор физико-математических наук

А.Н. Осипцов

Общая характеристика работы

Актуальность темы

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

В работах А.Н.Голубятникова, С.И.Зоненко и Г.Г.Черного1 было установлено, что если свойства среды таковы, что имеется максимум инкремента нарастания возмущений, то могут возникать достаточно регулярные "пальцеобразные" структуры, обладающие свойством накопления массы и импульса вещества слоя. Данные эффекты, являющиеся одним из проявлений неустойчивости Рэлея-Тейлора, находят широкое применение не только в технических устройствах, основанных на взрывном метании металлических оболочек, но и при решении фундаментальных научных проблем.

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

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

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

В космических условиях анализ влияния неустойчивости Рихтмайера-

1ГЬлубятников А.Н., ЗоненкоС.И., ЧерныйГ.Г. Новые модели и задачи теории кумуляции. Успехи механики. 2005. Т.З. №1. С.31-93.

2RichtmyerR.D. Taylor instability in shock acceleration of compressible fluids. Communs. Pure and Ap-pl.Math. 1960. V.13. №2. P.297-319.

3МешковЕ.Е. Неустойчивость границы раздела двух газов, ускоряемой ударной волной. Изв. АН СССР. МЖГ. 1969. >5. С.151-158.

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

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

Цели работы

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

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

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

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

Научная новизна

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

2. Исследована структура течений при слиянии частей деформированного слоя. Найдено присутствие комплекса ударных волн. Их конфигурация зависит от времени и на разных стадиях

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

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

4. Развита методика двумерных расчетов проявления неустойчивости Рихтмайера-Мешкова с учетом процессов радиационного охлаждения. Выявлена значительная роль функции охлаждения в распределении плотности газа за ударной волной.

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

Достоверность результатов

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

Научная и практическая значимость

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

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

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

Апробация работы

Основные положения и результаты, вошедшие в диссертацию, докладывались и обсуждались на семинарах кафедры аэромеханики и газовой динамики (МГУ, Москва); на семинарах лаборатории теоретического моделирования (ИКИ РАН, Москва); на семинаре кафедры волновой и газовой динамики (МГУ, Москва).

Результаты, представленные в диссертации, докладывались автором на следующих И научных конференциях: конференциях "Ломоносовские чтения" (МГУ, Москва, 2006, 2007); конференциях молодых ученых "Фундаментальные и прикладные космические исследования" ИКИ РАН (Москва, 2007-2009); XV школы-семинаре «Современные проблемы аэрогидродинамики» (Сочи, 2007); конференции-конкурсе молодых ученых НИИ механики МГУ (Москва, 2007); XXXII Академических чтениях по космонавтике (МГТУ им. Н.Э.Баумана, Москва, 2008); XV Международной конференции студентов, аспирантов и молодых ученых "Ломоносов" (МГУ, Москва, 2008); Всероссийской конференции "Астрофизика высоких энергий сегодня и завтра" (ИКИ РАН, Москва, 2008); Всероссийском конкурсе инновационных проектов "У.М.Н.И.К."(МГУ, Москва, 2008).

В соавторстве с К.В. Краснобаевым и Г.Ю. Котовой за циклы работ "Моделирование кумулятивных эффектов при ускоренном движении плотных газовых оболочек" и "Развитие моделей формирования неоднородностей структуры плотных газовых оболочек в космических условиях" автор удостоен премии в номинации "Лучшая научная работа

Института" по результатам конкурсов ИКИ РАН в 2007 и 2008 гг. Публикации

Основные результаты диссертационной работы изложены в 16-ти печатных работах, 3 из них опубликованы в журналах, которые входят в перечень ВАК. Список работ приведен в конце автореферата.

Структура и объем диссертации

Диссертационная работа состоит из введения, четырех глав, заключения и списка литературы (114 наименований). Общий объем диссертации 134 страниц.

Краткое содержание работы

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

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

В общем случае используется одножидкостная модель идеальной многокомпонентной среды. Среда является смесью совершенных газов с показателем адиабаты 7 = 5/3. Скорости и температуры Т всех компонент считаются одинаковыми. Концентрация и масса атомов г-го сорта далее обозначены соответственно через щтллщ.

Уравнения неразрывности, движения, притока тепла и уравнение состояния имеют вид

ар ¿V ,

р 16гЬ[р)+ р<11уу = р = ?щкт' р = ^Ш{П{

\ ' "' « г

Плотность р, давление р и скорость V определены как функции пространственных декартовых координат (х, у) и времени ¿. Величина

к - постоянная Больцмана. Функция Ф(п;, Т) - это количество энергии, теряемой путем излучения единицей объема среды в единицу времени. Величина Ф(щ,Т) ~ р2А(Т), где вид эффективности охлаждения А(Т) известен и рассмотрен в работах М. Кафатоса и С. Фалле.

Система уравнений решается численным методом сквозного счета, основанным на монотонной разностной схеме типа TVD второго порядка аппроксимации. В разделе 1.2 описан этот метод и рассмотрены примеры тестовых расчетов.

Вторая глава посвящена исследованию развития двумерных адиабатических деформаций плотного слоя газа, ускоряемого под действием разности давлений с обеих его сторон. В разделе 2.1 дан краткий обзор основных приближенных теорий. Рассматривается построенная Г.Г. Черным с соавторами1 инерционная модель устойчивости ускоренно движущихся тонких плотных оболочек, предсказывающая образование "пальцеобразных" уплотнений и эффектов кумуляции массы и импульса. Представлены приближенные теории столкновения двух симметричных плоских потоков вещества, описывающие условия возникновения ударных волн в потоках и критерии существования кумулятивных струй. Изложение теорий основано на результатах работ М.А. Лаврентьева, Э. Orra, а также на материалах монографии по физике взрыва под редакцией Л.П. Орленко. Отмеченные теории используются в дальнейшем для интерпретации расчетных данных.

В разделе 2.2 моделируется невозмущенное движение плоского слоя для того, чтобы найти поле течения, которое затем в двумерной постановке принимается за основное. Рассматривается ускорение первоначально покоившегося (то есть v = О при t = 0) плоского слоя газа. Начальные условия имеют вид

х < хс : р = ро = const, р = pa = const Xc<Xi<Xc + ôty:p — pi = const, p — pi = const X > Xc + Ôq : P = P2 = const, P = P2 = const Здесь pi = ро и, следовательно, xc - координата тангенциального разрыва. Значению х — хс + ¿о, где ¿о = const - начальная толщина слоя, отвечает положение произвольного разрыва с р2 < Pi- Будем считать, что ро и р2 значительно меньше р\, то есть слой (область "1") достаточно плотный. В таком случае область "0" занята горячим газом сравнительно малой плотности, а температура в области "2" невелика. Анализ основных особенностей решения удобно проводить, вводя некоторый характерный масштаб течения £ и выражая пространственные переменные, а также

время, плотность, скорость и давление соответственно в единицах С, £/со

(со = у/ро/ро), Ро, со, Ро-

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

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

Для исследования влияния свойств среды проведены вычисления в широком диапазоне изменения плотности слоя. Типичными являются результаты расчетов с начальными величинами 6о = 0.25, ро = 1, ро = 1, р2/ро = 0.02 для достаточно плотного слоя с условиями (¡) - рг/ро = 50, Р2/Р0 = 1 и для слоя меньшей плотности (и) - рг/ро = 5, р2/ро = 0.1: Начальная амплитуда возмущений и длина волны равны А0 = —0.1 и Л = 5.

Об основных особенностях развития двумерных деформаций дает представление рис.1. Уплотнения со значительно различающимися в продольном и в поперечном направлениях масштабами формируются при начальной плотности слоя существенно (более чем на порядок) превышающей плотность окружающей среды, то есть с условиями (1).

Количественные эффекты кумуляции массы характеризуем интегральной величиной те = (М5 — М|)/М*, где

Хп 1„ Хс+Лз

Мз(у, 0 = / р(х, у, 1)йх, Мд(0) - / р(х, 0,0)вх, М*(0) = / р{х, 0,0)<1х.

Хо х0 Хс

Здесь хо, хп - координаты, ограничивающие расчетную область.

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

Рис. 1. Линии уровня плотности для 4 = 8 при условиях (¡) и (п), на каждом графике изохор слева одномерные расчеты (Ао = 0), справа - двумерные. Интегральные массы тс(у) при 4 = 4,8 для (¡) - сплошные линии и (и) - штриховые кривые.

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

Для объяснения выявленных отличий проводится анализ течений газа в слое. Согласно инерционным теориям на нелинейной стадии эволюции возмущений происходит столкновение соседних частей слоя. Качественное описание дальнейшего поведения потоков дают приближенные модели (Э. Отт, Л.П. Орленко). Предполагается, что потоки имеют одинаковые постоянные скорости и возможно выделить точку столкновения потоков. Движения относительно системы отсчета, связанной с этой точкой, могут быть разделены на два типа.

Первый тип реализуется в случае, когда скорости потоков дозвуковые. Для таких движений характерно образование в области слияния как опускающихся "вниз", так и поднимающихся "вверх" частей уплотнений (Г.Биркгоф и др., М.А.Лаврентьев). Качественно такая картина течения согласуется с расчетами при рг/ро = 5 (рис. 1, (п)), что объясняется существенно большей, чем для р\/ро = 50, скоростью звука в слое.

Второй тип движений соответствует сверхзвуковым потокам и характеризуется возникновением ударных волн в слое (Дж. Уолш и др.). В зависимости от параметров среды возможны режимы как пересечения ударных волн на оси слияния потоков (у = 0), так и его отсутствия (Л.П. Орленко). Появление ударных волн иллюстрирует рис.2(а,б) для плотного слоя, удовлетворяющего условиям (¡). Слева на графиках (а) и (б) показаны сплошными линиями изобары, а штриховыми - положение

Рис. 2. Для каждого графика о ^=6) и б (1=8): линии уровня давления (слева) и числа | Маха с векторами относительных скоростей (справа). Условия (¡): р\/ро = 50, рг/ра = 1.

слоя, ограниченного тангенциальными разрывами и и положения возникающих поверхностей разрывов (5*, 6"*, Л*). Справа от изобар изображены числа Маха, вычисленные для абсолютных скоростей V, и векторы относительных скоростей газа, рассмотренные в системе координат, связанной с движением невозмущенного тангенциального разрыва Возникающие при столкновении потоков области с резкими изменениями параметров газа (5*,¿>2,¿>2) являются скачками уплотнения, что было установлено путем анализа соотношений Рэнкина-Гюгонио. Взаимодействие ударной волны с тангенциальным разрывом £>1 должно сопровождаться появлением волны разрежения Я*, что и показывают расчеты (изобары при £ = 8).

На начальной стадии слияния частей слоя проекции скоростей V на ось у вблизи линии у = 0 малы и набегающие потоки почти параллельны оси х. Это влечет образование пересекающихся ударных волн 5* (рис. 2а). Однако по мере нарастания деформаций слоя угол между сталкивающимися частями увеличивается и возможен режим отхода ударных фронтов 5* в сторону набегающих потоков (рис. 26). В целом результаты численного моделирования показывают, что конфигурация волн сложнее, чем это предсказывается приближенными теориями.

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

______

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

Рассматриваемые в этой главе движения основаны на представлениях о развитии неустойчивостей Рэлейя-Тейлора и Рихтмайера-Мешкова. Поэтому в разделе 3.1 излагаются основополагающие результаты теорий Тейлора4 и Рихтмайера2, а также экспериментов Е.Е. Мешкова3. Дается краткий обзор современных направлений исследования неустойчивостей данного типа и физических аспектов теории, обсуждаемых в монографиях А.Ю. Демьянова, H.A. Иногамова и Э.И. Сона5,6.

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

Раздел 3.3 посвящен анализу влияния импульсного и постепенного режима ускорения плотного газа на характеристики движения деформированного тангенциального разрыва и на эффекты кумуляции массы. Моделирование двух типов ускоренного движения тангенциального разрыва D проведено в условиях, близких к реализованным в

4Taylor J. The instability of liquid surfaces when accelerated in a direction perpendicular to their planes I. Proc. Roy Soc. London. Ser. A. 1950. V. 201. № 1065. P. 192-196.

5ИногаыовН,А., Демьянов А.Ю., Сон Э.И. Гидродинамика перемешивания: периодические структуры, усиление субгармоник, инверсный каскад. М.: Изд-во МФТИ, 1999. 464 с.

"InogamovN.A. The Role of RT and RM Instabilities in Astrophysics: an Introduction. Astroph. and Space Phys. Rev. 1999. V. 10. pt. 2.

экспериментах Е.Е. Мешкова. Результаты расчетов представлены на рис. 3, где под расстоянием х понимается тот путь, который проходит невозмущенный тангенциальный разрыв. Видно, что для границ Д прошедших один и тот же путь, амплитуды смещения А в случае постепенного ускорения (кривые 4> 5, 6) могут быть больше их значений при импульсном ускорении (1).

Рис.3. Зависимость амплитуды смещения границы А/Ац от пройденного пути х/А при импульсном ускорении границы (1 - численный расчет, 2 - эксперимент3, 3 -линейная интерполяция экспериментальных данных3) и при движении слоя толщины ¿о/А = 0.075,0.2,0.5 (линии 4, 5, б).

Эффекты накопления массы в случае неустойчивости Рихтмайера-Мешкова характеризуются интегральной массой газа в слое между деформированным тангенциальным разрывом и преломленной ударной волной. При этом тпе(у) вычисляется по формулам, аналогичным случаю постепенно ускоряемого слоя (гл. 2). Только для М* используется

хс+А

выражение М*(0) = / р(х, 0,где хс - положение невозмущенного

Хс

тангенциального разрыва.

При сопоставлении те для двух режимов ускоренного движения границы найдено, что наиболее отчетливо эффекты кумуляции массы проявляются при развитии возмущений достаточно тонких слоев газа (например, для 50/А — 0.075), подвергшихся постепенному ускорению.

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

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

В разделе 4.1 приведены имеющиеся в настоящее время сведения о параметрах звездного ветра и ударных волн, генерируемых при вспышках сверхновых звезд, а также о физических процессах высвечивания в околозвездном газе. В рассматриваемых астрофизических объектах среда является полностью ионизованной оптически прозрачной водородной плазмой с малыми примесями ионов тяжелых элементов. Тепловое состояние горячей межзвездной среды определяется неупругими столкновениями частиц (ионов и электронов). В общем виде функция охлаждения представляется как Ф = прпеА(Т), где Л(Т) - эффективность охлаждения, обусловленного радиационными процессами, а пр и пе -концентрации протонов и электронов. В расчетах используется Л(Т), вычисленная в работах Д.П. Кокса, М. Кафатоса и С. Фалле. Характерной для выбранного вида кривой охлаждения является немонотонная зависимость от температуры. В частности, имеет место как возрастание эффективности охлаждения с увеличением Т [104 — 2 х 105].ЙГ), так и убывание А(Т) с ростом температуры 2х [105 — 107]К). Таким видом Л(Г), как показали расчеты, обусловлены важные особенности течения с высвечиванием.

Раздел 4.2 посвящен численному моделированию развития возмущений тангенциального разрыва, возникающего при торможении высокоскоростного звездного ветра межзвездной средой. При взаимодействии звездного ветра с окружающей плазмой возникает нестационарное течение и вокруг звезды образуется двухслойная оболочка. Внешний слой (область "2") ограничен ударной волной 5е в ионизованном межзвездном газе, а внутренний ("1") - ударной волной 5* в звездном ветре. Обычно плотность р\ <С ръ а давление р\ = р2, поэтому для температуры имеем Т\ Тг. Тангенциальный разрыв И между слоями отделяет горячую плазму звездного ветра от менее нагретой межзвездной среды. При падении на поверхность Б ударной волны со стороны горячей плазмы создаются условия для развития неустойчивости Рихтмайера-Мешкова. Происхождение таких ударных волн может быть обусловлено нестационарными процессами в звездном ветре, например, изменениями

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

Рассмотрено взаимодействие плоской ударной волны 5 с неподвижным тангенциальным разрывом, положение которого в произвольный момент времени t задается равенством хц — %о{у> О- Для определенности при Ь = 0 положено хв(0) = 0, а концентрация и температура плазмы, следуя работе С.Фалле7, равны Пю = рю/тн = 0.0735 см-3, Тю = 2.16 х 107К и п2о = Р2о/тпн = 2.57см-3, Гго = 6.14 х 105К соответственно в невозмущенных прохождением ударной волны областях "1" и "2" (га#-масса атома водорода). Скорость падающей ударной волны V является задаваемым параметром задачи и поэтому она лишь косвенно отражает нестационарные процессы в звездном ветре. В расчетах V выбиралась так, чтобы увеличение давления за фронтом ударной волны было порядка исходного давления в среде. В приводимом примере V = 1.15 х 108см/с. В качестве характерных величин задачи определены скорость преломленной ударной волны С/5, время охлаждения ¿с, концентрация пщ, температура Тю и пространственный масштаб движения С = С/д • ¿с-

Рис.4. Два типа распределения плотности после взаимодействия ударной волны с плоским тангенциальным разрывом. Тип а - температура газа перед преломленной ударной волной Str понижается со временем и тип б - температура остается постоянной. Сплошные линии - движение с высвечиванием, штриховые - адиабатическое течение. Время t = 2,5 - соответственно для кривых 1, 2.

7FalleS.A.E.G. A numerical calculation of the effect of stellar winds on the interstellar medium. Astron. Astrophys.J. 1375. V. 43. P. 323-336.

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

процесс охлаждения газа приводит к увеличению значения плотности, а также к уменьшению скорости преломленной волны 5гг- Характерна также немонотонная зависимость р от координаты х. Причем наибольшее значение плотности достигается не вблизи тангенциального разрыва £), а в области перед ним. Это обусловлено особенностями кривой охлаждения. Температура в межзвездной среде в приводимом примере ~ 106 К. Для ударных волн умеренной интенсивности темп охлаждения при переходе через скачок не меняется по порядку величины. Это означает, что из-за охлаждения со временем заметно уменьшаются температура и давление перед фронтом волны и за ним. Но в рассматриваемом диапазоне температур А(Т) возрастает при уменьшении Т. Следовательно, частицы газа, прошедшие через ударную волну в более поздние моменты времени, охлаждаются быстрее, чем те, что находятся вблизи тангенциального разрыва. Наибольшее сжатие происходит ближе к ударной волне, чем к границе £>. Факт влияния переменности темпа охлаждения на распределение р подтверждается расчетом на рис. 46, где перед скачком "¿V температура принималась постоянной. Видно, что теперь плотность монотонно возрастает по направлению к разрыву Б. Выявленная для случая, представленного на рис. 4а, особенность распределения плотности за фронтом ударной волны с высвечиванием ранее не отмечалась в литературе.

В расчетах принималось, что невозмущенная прохождением ударной волны среда и находящийся в ней тангенциальный разрыв покоятся. Применительно к течению звездного ветра это предположение справедливо, если длины волн возмущений Л малы по сравнению с расстояниями между поверхностями £¿,1), 5е и характерное время нарастания возмущений í¿nst невелико по сравнению со временем изменения параметров звездного ветра. В работе выбор таких значений А и был проведен с использованием расчетов динамики звездного ветра, выполненных С. Фалле7. Поскольку эффекты радиационного охлаждения наиболее существенны, когда порядка Ъ, то соответствующим

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

поверхности £>, так что хц{у, 0) = А0 соз(2пу/Х), А0 = 0.4, Л = 4.

Слева от графиков (рис. 5) приведены также изохоры в случае одномерного течения = ££)(£), хц(0) = 0. Образующиеся области

между преломленной ударной волной и границей Г> значительно более плотные, чем при А(Т) = 0. Причем толщина сжатого ударной волной слоя газа существенно уменьшается из-за радиационных потерь энергии. В то же время глубина проникновения более горячей плазмы в холодную ("пузырь") и холодной в горячую ("струя") не слишком отличаются от имеющих место при адиабатическом течении (А.Н. Алешин и др.). Достаточно слабо выражены и различия в эффектах накопления массы в областях газа, отстающих от основного движения границы £>. Видно, что интегральные характеристики тпе(у), полученные для адиабатического течения и при учете высвечивания, отличаются мало.

Рис.5. Изохоры и функция накопления массы тпс(у) при t=5 для движения с высвечиванием (а) и для адиабатического течения (б).

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

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

изначально в межзвездной среде. Динамика такого взаимодействия и доминирующий механизм деформации контактной границы зависят от параметров ударной волны, геометрии и характерного размера конденсации Ь. В работе установлен такой диапазон изменения Л и Ь, в пределах которого развитие неустойчивости Рихтмайера-Мешкова произойдет раньше, чем ударная волна пересечет облако. Цель расчетов состоит в определении таких масштабов возмущений, для которых форма облака существенно изменится, прежде чем проявятся ранее исследовавшиеся более крупномасштабные эффекты (ускорение конденсаций как целого объекта; фокусировка огибающих уплотнение ударных волн; неустойчивость Кельвина-Гельмгольца). В разделе исследуется эволюция деформаций тангенциального разрыва, определяются их форма и скорость роста в зависимости от длины волны и от времени высвечивания. Для фиксированного момента времени найдены наибольшие размеры уплотнений.

Рис.6. Зависимость роста амплитуды от длины волны начального возмущения для движения с высвечиванием (сплошные линии) и для адиабатического течения (штриховые).

Конкретные расчеты проведены при условиях, характерных для остатка вспышки сверхновой G74.0-8.5 "Cygnus Loop" (H.A.Левенсон и др.). При этом параметры неоднородного межзвездного газа, по которому распространяется ударная волна, выбирались согласно данным С. Орландо и др.8 Задача о развитии возмущений рассматривалась в той же

8OrlandoS., PeresG., RealeF. et al. Crushing of interstellar gas clouds in supernova remnants - I. The role of thermal conduction and radiative losses. Astron. Astrophys. 2005. V. 444. P. 505-519.

постановке, что и в разд. 4.2. Следуя оценкам С. Орландо и др., в расчетах принято пю = 0.1см~3, Тш = 10\ п20 = 1см-3, Т20 = 103К, V = 3.88 х 107см/с. В соответствии с оценками ~ если Ь ~ 1пк

(парсек). Отсюда Ьс < Ь/Пв, когда Ь > 1пк.

Особенности нарастания двумерных деформаций тангенциального разрыва исследовались для различных длин волн. Использовались следующие начальные условия для возмущений: ггд(у, 0) = Аосоз(27гу/А), А0 - 0.125, Л = (0.375 - 1.5). Амплитуда возмущения А определялась как разность между глубинами проникновения горячего газа в холодный и холодного в горячий. Найденный в расчетах вид Л(А) в фиксированные моменты времени представлен на рис.6. Как и в экспериментах А.Н.Алешина и др., графики Л(А) имеют максимум, смещающийся со временем в область более длинных волн. Однако при Л ф 0 этот максимум значительно более резкий, чем в случае Л = 0. Физически это объясняется тем, что с ростом длины волны уменьшается инкремент нарастания возмущений. Следовательно, для достаточно длинных волн ипвь > поэтому прежде чем А(А) увеличится, возникнет слабо искривленный слой очень плотной плазмы, испытывающий малое смещение после действия ударной волны. Тогда в соответствии с теорией Рихтмайера должна быть малой также амплитуда А(А).

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

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

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

2. Количественно исследован режим постепенного ускорения плотного

газового слоя под действием разности давлений с обеих его

сторон. Установлены критерии формирования достаточно регулярных

"пальцеобразных" уплотнений со значительно различающимися в

продольном и в поперечном направлениях масштабами. Выявлены

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

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

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

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

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

7. Найдены такие масштабы возмущений, для которых морфология неоднородной среды вследствие неустойчивости Рихтмайера-Мешкова существенно изменится, прежде чем проявятся более крупномасштабные эффекты, ранее анализировавшиеся в литературе.

Работы автора по теме диссертации

[1] КотоваГ.Ю., Краснобаев К. В., ТагироваР.Р. Сравнительный анализ механизмов резонансного усиления возмущений ионизационно-ударного фронта // Ломоносовские чтения. Секция механики. 18-28 апреля 2006 г. Тезисы конференции. М.: Изд-во Моск. Ун-та, 2006. С. 94.

[2] Краснобаев К.В., ТагироваР.Р. Кумуляция массы и импульса в ускоренно движущейся сжимаемой среде //IV конференция молодых ученых "Фундаментальные и прикладные космические исследования". 12 апреля 2007г. Тезисы докладов. М.: ИКИ РАН, 2007. С. 44.

[3] Краснобаев К.В., ТагироваР.Р. Накопление массы и импульса в ускоренно движущейся сжимаемой среде // Ломоносовские чтения. Секция механики. 16-25 апреля 2007г. Тезисы конференции. М.: Изд-во Моск. Ун-та, 2007. С. 95.

[4] Краснобаев К. В., ТагироваР.Р. Двумерное неустановившееся ускоренное движение плотного слоя газа // Тезисы докладов XV школы-семинара "Современные проблемы аэрогидродинамики". 5-15 сентября 2007 г. М: Изд-во Моск. Ун-та, 2007. С. 63.

[5] Тагировa. P.P. Влияние переменного ускорения на развитие возмущений поверхности тангенциального разрыва // Тезисы докладов XV школы-семинара "Современные проблемы аэрогидродинамики". 5-15 сентября 2007 г. М: Изд-во Моск. Ун-та, 2007. С. 91-92.

[6] Тагирова. P.P. Деформации поверхности тангенциального разрыва в газах, движущихся с ускорением // Труды конф.-конкурса молодых ученых НИИ механики МГУ. 10-12 октября 2007 г. Под ред. акад. Г.Г.Черного и проф. В.А.Самсонова. М.: Изд-во Моск. Ун-та, 2009. С. 206-212.

[7] КотоваГ.Ю., Краснобаев К. В., ТагироваР.Р. Двумерные неустановившиеся движения фотоиспаряемых газовых оболочек // Сб. Проблемы современной механики: к 85-летию со дня рождения академика Г.Г.Черного. / Под ред. A.A.Бармина - М.: Изд-во Моск. Ун-та, 2008. С. 190-206.

[8] Краснобаев К. В., ТагироваР.Р. Исследования развития неустойчивости тангенциального разрыва под действием импульсного и постепенного ускорения // Актуальные проблемы российской космонавтики: Труды XXXII Академических чтений по космонавтике. 29 января-1 февраля 2008г. Под ред. А.К.Медведевой. М.: Комиссия РАН, 2008, С. 185-186.

[9] Тагирова. P.P. Ускоренное движение излучающей деформированной газовой оболочки // V конференция молодых ученых "Фундаментальные и прикладные космические исследования". 9 апреля 2008г. Тезисы докладов. М.: ИКИ РАН, 2008. С. 41.

[10] Тагирова. P.P. Развитие возмущений ускоренно движущегося газового слоя // Материалы докладов XV Международной конференции студентов, аспирантов и молодых ученых "Ломоносов". Подсекция "Механика". Апрель 2008. М.: СП "Мысль", 2008. С. 15-16.

[11] Краснобаев К.В., Тагирова P.P. Моделирование неустойчивости ускоренно движущейся газовой оболочки // Изв. РАН. МЖГ. 2008. №5. С. 163-172.

[12] Тагирова. P.P. Взаимодействие ударной волны от сверхновой с неоднородностями в межзвездной среде // Всероссийская конференция "Астрофизика высоких энергий сегодня и завтра". 24-26 декабря 2008. Аннотации докладов. ИКИ РАН. 2008. С. 36.

[13] Тагирова. P.P. Моделирование двумерных деформаций ускоренно движущихся контактных поверхностей в сжимаемых средах // Сборник инновационных проектов "У.М.Н.И.К.". М: Изд-во Моск. Унта, 2009. С. 78-81.

[14] Тагирова. P.P. Роль переменности темпа высвечивания в распределении параметров плазмы за ударной волной // VI конференция молодых ученых. "Фундаментальные и прикладные космические исследования". 9-10 апреля 2009 г. Тезисы докладов. М.: ИКИ РАН, 2009. С. 40.

[15] Краснобаев К.В., Тагирова P.P. О проявлении неустойчивости Рихтмайера-Мешкова в неоднородной межзвездной среде с высвечиванием // Письма в Астрон. ж. 2009. Т. 35. № 5. С. 364371.

[16] КопговаГ.Ю., Краснобаев К.В., Тагирова P.P. Неустановившиеся движения плотных оболочек, порождаемых ионизационно-ударным фронтом в межзвездной среде / / Вести. МГУ. Сер.1. Математика, механика. 2009. № 3. С. 56-62.

Подписано в печать . //. 03 Формат 60x90 1/16. Усл. печ. л. (5 Тираж {00 экз. Заказ 52

Отпечатано с оригинал-макета на типографском оборудовании механико-математического факультета МГУ имени М. В. Ломоносова

Введение.

1 Модель среды, численный метод расчета течения

1.1 Основные уравнения, их приведение к безразмерному виду

1.2 Численный метод решения. Результаты тестовых расчетов

1.3 Выводы

2 Двумерные адиабатические деформации слоя, ускоряемого под действием разности давлений с обеих его сторон

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

2.2 Одномерное ускоренное движение плоского слоя газа

2.3 Численное моделирование эволюции двумерных возмущений параметров слоя. Эффекты накопления массы. Структура поля течения

2.4 Выводы

3 Нарастание двумерных адиабатических возмущений при импульсном и постепенном ускорении границы.

3.1 Элементы теории неустойчивости Рэлея-Тейлора и Рихтмайера-Мешкова.

3.2 Моделирование неустойчивости Рихтмайера-Мешкова с помощью TVD-схемы Лакса-Фридрихса.

3.3 Анализ движений возмущенных границ при различных режимах ускорения

3.4 Выводы

4 Развитие неустойчивости Рихтмайера-Мешкова в среде с высвечиванием

4.1 Физические процессы в межзвездной среде.

4.2 Развитие деформаций тангенциального разрыва, возникающего при истечении звездного ветра.

4.3 Взаимодействие ударной волны от сверхновой с неоднородной межзвездной средой

4.4 Выводы

Проблеме эволюции возмущений ускоренно движущейся границы, разделяющей среды различной плотности, посвящены многочисленные экспериментальные и теоретические исследования [111, 109, 46, 3, 23, 14, 48, 86]. Основополагающей является классическая работа Дж. Тейлора [111], в которой в линейной постановке было установлено, что неустойчивость тангенциального разрыва развивается при действии ускорения, перпендикулярно направленного к разрыву, от менее плотной среды к более плотной (неустойчивость Рэ лея-Тейлор а). В дальнейшем значительное внимание уделялось изучению роли неустойчивости Рэлея-Тейлора в процессах перемешивания жидкостей и газов, вопросам стабилизации возмущений, анализу влияния на рост возмущений физических свойств среды и ускорения, зависящего от времени [11, 12, 19, 49, 7, 18, 23, 50, 10, 86]. В настоящее время значительный интерес представляет исследование ускоренного движения деформированных тонких слоев плотного вещества под действием разности давлений с обеих сторон слоя. Этот интерес обусловлен тем, что, как было показано в работах А.Н.Голубятникова, С.И.Зоненко и Г.Г.Черного [21, 14, 15], если свойства среды таковы, что имеется максимум инкремента нарастания возмущений, то могут возникать достаточно регулярные структуры "пальцеобразной" формы с эффектами накопления массы и импульса вещества слоя. Данные эффекты находят широкое применение не только в технических устройствах, основанных на взрывном метании металлических пластин и оболочек [17, 52], но и при решении фундаментальных научных проблем. Например, проблемы устойчивости газовых оболочек представляют интерес в космической газовой динамике при исследовании морфологических особенностей неоднородной структуры межзвездной среды [73, 96, 89, 32, 74, 83, 101, 25, 26]. Именно приложениям к этому кругу задач в основном и посвящена настоящая работа. При этом рассматриваются плоские двумерные неустановившиеся движения идеального газа.

Во многих астрофизических объектах наблюдается ускоренное движение плотных холодных слоев или оболочек: в расширяющихся областях газа, нагретого излучением горячей звезды; в окрестности "межзвездных пузырей", выдуваемых звездным ветром; в остатках сверхновых звезд. Особенностью таких оболочек является то, что в них часто присутствуют неоднородности плотности с сильно различающимися в продольном и в поперечном направлениях размерами. Как показывают наблюдения, массы таких неоднородностей могут быть порядка масс звезд и планет. Следовательно, моделирование развития подобных возмущений является необходимой составной частью общей теории возникновения самогравитирующих конденсаций [113, 5, 79, 24, 55, 89, 100, 103].

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

Однако уплотнения могут возникать и под действием иных механизмов развития возмущений границы между газами различной плотности. В частности, в космических условиях существенна неустойчивость Рихтмайера-Мешкова [109, 46], обусловленная импульсным ускорением тангенциального разрыва под влиянием падающей на разрыв ударной волны. При этом проявление неустойчивости может косвенно свидетельствовать об условиях возникновения и формирования слоев [73, 113, 110, 71, 32, 102, 74]. Поэтому еще одна цель диссертационной работы состоит в определении влияния режимов ускоренного движения на процесс образования' и морфологию уплотнений. Здесь принимаются во внимание две характерные для космических условий особенности движений. Первая состоит в том, что с тангенциальным разрывом могут взаимодействовать ударные волны как умеренной, так и очень большой интенсивности (динамическое давление на несколько порядков превышает тепловое) [114, 98, 94, 93, 107]. Вторая особенность связана с тем, что на динамику газа существенно влияют процессы радиационного охлаждения [5, 59, 95]. Применительно к явлениям в звездных ветрах это обстоятельство отмечалось уже в ранних работах [80, 24]. Важная роль высвечивания установлена также в случае взаимодействия сильных ударных волн от вспышек сверхновых с неоднородностями в межзвездной среде [99, 81, 102, 97]. С учетом указанных особенностей в диссертации моделируется развитие неустойчивости Рихтмайера-Мешкова в случаях, когда тангенциальный разрыв отделяет горячую плазму звездного ветра от менее нагретого межзвездного газа и когда ударная волна от вспышки сверхновой звезды взаимодействует с границей области плотного вещества (конденсации), содержащегося изначально в межзвездной среде.

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

Структура диссертации

Диссертационная работа состоит из введения, четырех глав, заключения и списка литературы.

4.4 Выводы

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

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

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

Результаты, полученные в данной главе, изложены в [62, 63, 65, 64, 66, 38].

Заключение

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

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

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

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

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

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

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

7. Найдены такие масштабы возмущений, для которых морфология неоднородной среды вследствие неустойчивости Рихтмайера-Мешкова существенно изменится, прежде чем проявятся более крупномасштабные эффекты, ранее анализировавшиеся в литературе.

1. Алешин А.Н., ЛазареваЕ.В., Зайцев С.Г. и др. Исследование линейной, нелинейной и переходной стадий развития неустойчивости Рихтмайера-Мсшкова // Докл. АН. СССР. 1990. Т. 310. №5. С. 11051108.

2. Алешин А.Н., Демченко В.В., Зайцев С.Г. и др. Взаимодействие ударного фронта с волнообразным контактным разрывом // Изв. АН СССР. МЖГ. 1992. №5. С. 168-174.

3. Анисимов С.И., Прохоров A.M., Фортов В.Е. Применение мощных лазеров для исследования вещества при сверхвысоких давлениях / / Успехи физ. наук. 1984. Т. 142. №3. С. 395-434.

4. Баранов В.В., Ермаков М.К., Лебедев М.Г. Трехкомпонентная газодинамическая модель взаимодействия солнечного ветра с межзвездной средой. Изв. АН СССР. МЖГ. 1982. № 5. С. 122-128.

5. Баранов В.В., Краснобаев К.В. Гидродинамическая теория космической плазмы. М.: Наука, 1977. 335 с.

6. Баранов В.В., Краснобаев К.В., Куликовский А.Г. Модель взаимодействия солнечного ветра с межзвездной средой // Докл. АН. СССР. 1970. Т. 194. №1. С. 41-44.

7. Бахрах С.М., Симонов Г. П. Аналитическое и численное исследование Рэлей-Тейлоровской неустойчивости тонкого жидкого слоя // ВАНТ, серия "Математическое моделирование физических процессов". 1995. №3. С. 39-46.

8. Бочкарев Н.Г. Основы физики межзвездной среды. М.: Изд-во МГУ, 1991. 352 с.

9. Бранд Дэю. Солнечный ветер. Введение в проблему. М.: Мир, 1973.

10. ВолковН.Б., МайерА.Е., ЯловецА.П. Нелинейная динамика контактной границы сплошных срсд с различной плотностью // ЖТФ. 2003. Т. 73. №3. С. 1-9.

11. ГерценштейнС.Я., Чернявский В. М. О нелинейном развитии двумерных и трехмерных возмущений при неустойчивости Рэлея-Тейлора//Изв. АН СССР. МЖГ. 1985. №2. С.38-46.

12. ГерценштейнС.Я., Чернявский В. М., ШтемлерЮ.М. О неустойчивости Рэлея-Тейлора при больших значениях времени //Изв. АН СССР. МЖГ. 1989. № 5. С. 8-18.

13. Годунов С.К., Забродин А.В., Иванов М.Я. и др. Численное решение многомерных задач газовой динамики. М.: Наука, 1976. 400с.

14. Голубятников А.Н., Зоненко С.И., Черный Г. Г. Новые модели и задачи теории кумуляции // Успехи механики. 2005. Т. 3. №1. С. 31-93.

15. Голубятников А.Н., Зоненко С.И., Черный Г.Г. Новые модели оболочек, метаемых взрывом // Прикладная матем. и мех. 2007. Т. 71. С. 727-743.

16. Демченко В.В., Сергеев М. А. Гидродинамическая неустойчивость при высокоскоростном ударе // Матем. мод-ие. 2002. Т. 14. №10. С. 87-94.

17. Зайцев С.Г., КривецВ.В., Титов С.Н. и др. Развитие неустойчивости Рэлея-Тэйлора в сжимаемых средах // Изв. АН. МЖГ. 1999. №3. С. 1625.

18. Зайцев С.Г., ЛебоИ.Г., Розанов В.Б. и др. Гидродинамическая неустойчивость области контакта газовых сред, движущихся ускоренно // Изв. АН. МЖГ. 1991. №6. С. 15-21.

19. Зельдович Я.В., Райзер Ю.П. Физика ударных волн и высокотемпературных явлений. М.: Наука. Изд. 2-ое, 1963. 688 с.

20. Зоненко С.И., ЧерныйГ.Г. Новый вид кумуляции энергии и импульса метаемых взрывом пластин и оболочек // Докл. РАН. 2003. Т. 390. № 1. С. 46-50.

21. ИмшенникВ.С. Вспышки сверхновых и историческая хронология // Успехи физ. наук. 2000. Т. 170. №5. С. 553-557.

22. Иногамов Н.А., Демьянов А.Ю., СонЭ.И. Гидродинамика перемешивания: периодические структуры, усиление субгармоник, инверсный каскад. М.: Изд-во МФТИ, 1999. 464 с.

23. КапланС.А., Пикелънер С.Б. Физика межзвездной среды. М.: Наука, 1979.591с.

24. КотоваГ.Ю., Краснобаев К. В. Нелинейные деформации ускоренно движущейся излучающей оболочки // Химическая физика. 2008. Т. 27. №5. С. 81-86.

25. КотоваГ.Ю., Краснобаев К.В. Ускорение сферической нейтральной оболочки, формируемой ионизационно-ударным фронтом в неоднородной межзвездной среде // Письма в Астрон. журн. 2009. Т. 35. №3. С. 189-198.

26. КотоваГ.Ю., Краснобаев К. В., Тагирова P.P. Сравнительный анализ механизмов резонансного усиления возмущений ионизационно-ударного фронта // Ломоносовские чтения. Секция механики. 18-28 апреля 2006 г. Тезисы конф. М.: Изд-во Моск. Ун-та. 2006. С. 94.

27. КотоваГ.Ю., Краснобаев К. В., Тагирова P.P. Двумерные неустановившиеся движения фотоиспаряемых газовых оболочек // Сборник: Проблемы современной механики / Под ред. А.А. Бармина. М.: Изд-во МГУ, 2008. С. 190-206.

28. КотоваГ.Ю., Краснобаев К.В., Тагирова Р.Р. Неустановившиеся движения плотных оболочек, порождаемых ионизационно-ударным фронтом в межзвездной среде // Вестн. МГУ. Сер.1. Математика, механика. 2009. №3. С. 56-62.

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

30. Краснобаев К.В. Неустойчивость газового слоя, ускоряемого "реактивным" механизмом / / Нелинейные задачи теориигидродинамической устойчивости. Материалы университетской школы. Изд-во МГУ, Москва. 1976. С. 39-40.

31. Краснобаев К. В. Неустойчивость тонкой фотоиспаряемой околозвездной оболочки// Письма в Астрон. ж. 2004. Т. 30. №7. С. 500-505.

32. Краснобаев К.В., ТагироваР.Р. Кумуляция массы и импульса в ускоренно движущейся сжимаемой среде //IV конференция молодых ученых "Фундаментальные и прикладные космические исследования". 12 апреля 2007г. Тезисы докладов. М.: ИКИ РАН, 2007. С. 44.

33. Краснобаев К.В., ТагироваР.Р. Накопление массы и импульса в ускоренно движущейся сжимаемой среде // Ломоносовские чтения. Секция механики. 16-25 апреля 2007г. Тезисы конференции. М.: Изд-во Моск. Ун-та, 2007. С. 95.

34. Краснобаев К.В., ТагироваР.Р. Двумерное неустановившееся ускоренное движение плотного слоя газа // Тезисы докл. XV школы-семинара "Современные проблемы аэрогидродинамики". 5-15 сентября 2007г. М: Изд-во Моск. Ун-та, 2007. С. 63.

35. Краснобаев К.В., ТагироваР.Р. Моделирование неустойчивости ускоренно движущейся газовой оболочки // Изв. РАН. МЖГ. 2008. №5. С. 163-172.

36. Краснобаев К.В., ТагироваР.Р. О проявлении неустойчивости Рихтмайера-Мешкова в неоднородной межзвездной среде с высвечиванием // Письма в Астрон. ж. 2009. Т. 35. №5. С. 364-371.

37. Куликовский А.Г., Погорелое Н.В., Семенов А.Ю. Математические вопросы численного решения гиперболических систем уравнений. М.: Физматлит, 2001. 608 с.

38. ЛаврентьевМ.А. Кумулятивный заряд и принципы его работы // Успехи мат. наук. 1957. Т. 12. №4. С. 41-56.

39. ЛаврентьевМ.А., ШабатБ.В. Проблемы гидродинамики и их математические модели. М.: Наука, 1973. 416 с.

40. Ландау Л.Д., ЛифшицЕ.М. Механика сплошных сред. М.: Гостех-издат, 1954. 788 с.

41. ЛенгК. Астрофизические формулы. T.l. М.: Мир, 1978, 448 с.

42. Лозинская Т.А. Сверхновые звезды и звездный ветер. Взаимодействие с ветром галактики. М.: Наука, 1986.

43. ЛойцянскийЛ.Г. Механика жидкости и газа. М.: Наука, 1970. 904с.

44. Мешков Е.Е. Неустойчивость границы раздела двух газов, ускоряемой ударной волной // Изв. АН СССР. МЖГ. 1969. №5. С. 151-158.

45. Мешков Е.Е. Исследования гидродинамических неустойчивостей в лабораторных экспериментах. Саров: ФГУП "РФЯЦ-ВНИИЭФ", 2006. 139 с.

46. Неуважаев В.Е., Паршуков И.Э. Изучение устойчивости границ раздела жидкостей при совместном действии импульсного и постоянного ускорений // Мат. моделирование. 1993. Т. 5. №2. С. 16-24.

47. Опарин A.M., Иногамов Н.А., Демьянов А.Ю. О спектральных и статистических свойствах Рэлей-Тейлоровского перемешивания // Письма в ЖЭТФ. 2000. Т. 72, №10, С. 704-710.

48. Овсянников Л.В. Лекции по основам газовой динамики. Москва-Ижевск: Институт компьютерных исследований. Изд. 2-е, 2003. 336с.

49. Физика взрыва /Под ред. Л.П. Орленко Т.2. М.: Физматлит, 2002. 648 с.

50. Пикельнер С.Б., Щеглов П. В. Движение газа в диффузных туманностях и звездный ветер // Астрон. ж. 1968. Т. 45. С. 953—961.

51. ПилюгинН.Н., ТирскийГ.А. Динамика ионизированного излучающего газа. М.: Изд-во МГУ, 1989. 309с.

52. Поттпаш С. Планетарные туманности. М.: Мир, 1987. 351с.

53. Седов Л.И. Методы подобия и размерности в механике. М.: Наука, Изд. 10-е. 1987.432 с.

54. Смирнов Н.Н. Нестационарное формирование кумулятивной струи в плотной среде // Вестник МГУ. Сер.1. Математика, механика, 1985, №6. С. 53-62.

55. Соболев В.В. Курс теоретической астрофизики. М.: Наука. Изд. 3-е. 1985. 504 с.

56. СпитцерЛ. мл. Физические процессы в межзвездной среде. М.: Мир, 1981.

57. ТагироваР.Р. Влияние переменного ускорения на развитие возмущений поверхности тангенциального разрыва / / Тезисы докладов XV школы-семинара "Современные проблемы аэрогидродинамики". 5-15 сентября 2007 г. М: Изд-во Моск. Унта, 2007. С. 91-92.

58. ТагироваР.Р. Ускоренное движение излучающей деформированной газовой оболочки //V конф. молодых ученых "Фундаментальные и прикладные космические исследования". 9 апреля 2008г. Тезисы докладов. М.: ИКИ РАН, 2008. С. 41.

59. ТагироваP.P. Развитие возмущений ускоренно движущегося газового слоя // Материалы докладов XV Международной конференции студентов, аспирантов и молодых ученых "Ломоносов". Подсекция Механика. Апрель 2008. М.: СП "Мысль", 2008. С. 15-16.

60. ТагироваР.Р. Взаимодействие ударной волны от сверхновой с неоднородностями в межзвездной среде // Всероссийская конференция "Астрофизика высоких энергий сегодня и завтра". 24-26 декабря 2008. Аннот. докладов. ИКИ РАН. 2008. С. 36.

61. ТагироваР.Р. Моделирование двумерных деформаций ускоренно движущихся контактных поверхностей в сжимаемых средах // Сборник инновационных проектов «У.М.Н.И.К.». М: Изд-во Моск. Унта, 2009. С.78-81.

62. УиземДж. Линейные и нелинейные волны. М.: Мир, 1977. 622 с. (

63. ЧерныйГ.Г. Газовая динамика. М.: Наука, 1988. 424 с.

64. Шкловский И. С. Сверхновые звезды и связанные с ними проблемы. М: Наука, 1976.

65. AllerL.H. The abundances of the elements // NY: Interscience publication. 1961.

66. BedogniR., Woodward P. R. Shock wave interactions with interstellar clouds // Astron. Astrophys. 1990. V.231. P. 481-498.

67. BirkhoffG., MacDougallD.P., PughE.M. TaylorG.F. Explosives with lined cavities // J. Appl. Phys. 1948. V. 19. №6. P. 563-582.

68. CapriottiE.R. The structure and evolution of planetary nebulae // Astro-phys. J. 1973. V. 179. P. 495-516.

69. CapriottiE.R., Kendall A. D. The origin and physical properties of the cometary knots in NGC 7293 // Astrophys. J. 2006. V. 642. P. 923-932.

70. ChevalierR.A., Self-similar solutions for the interaction of stellar ejecta with an external medium // Astrophys. J. 1982. V. 258. P. 790-797.

71. CoxD.P., DaltabuitE. Radiative cooling of a low-density plasma // Astrophys. J. 1971. V. 167. P. 113-117.

72. CoxD.P., Tucker W.H. Ionization equilibrium and radiative cooling of a low-density plasma // Astrophys. J. 1969. V. 157. P. 1157-1167.

73. Dalgarno A., McCrayR.A. Heating and ionization of HI regions // Ann. Revs. Astron. Astrophys. 1972. V. 10. P. 375-426.

74. ElmegreenB.GLada C.J. Sequential formation of subgroups in OB associations // Astrophys. J. 1977. V. 214. P. 725-741.

75. Falle S.A.E.G. A numerical calculation of the effect of stellar winds on the interstellar medium // Astron. Astrophys.J. 1975. V.43. P. 323-336.

76. FragileP.C., MurrayD.S., AnninosP. Radiative shock-induced collapse of intergalactic clouds // Astrophys. J. 2004. V. 604. P. 74-87.

77. Garcia-Segura G., Franco J. From ultracompact to extended HII regions // Astrophys. J. 1996. V. 469. P. 171-188.

78. GolubiatnikovA.N., Doroshenko T.A. On the dynamics of conducting envelope of magnetic star // Gravitation & Cosmology. 2006. V. 12. №2-3 (46-47). P. 140-143.

79. Graham J.R., LevensonN.A., Hester J. C. et al. An X-ray and optical study of the interaction of the Cygnus Loop supernova remnant with an interstellar cloud // Astrophys. J. 1995. V.444. P. 787-795.

80. HartenA. On a class of high resolution total-variation-stable finite-difference schemes// SIAM J. Numer. Anal. 1984. V.21. №1. P. 1-23.

81. HearnN.C., PlewaT., DrakeR.T. et al. Flash code simulations of Rayleigh-Taylor and Richtmyer-Meshkov instabilities in laser-driven experiments // Astrophys. Space Sci. 2007. V. 307. P. 227-231.

82. HoffmannK.A., ChiangS.T. Computational fluid dynamics. V. 1. KS, USA: Engineering Education Systems. 4th ed. 2000. 486 p.

83. Physics of the heliospheric. boundaries / Eds. V.V. Izmodenov, R. Kallcnbach // ESA Publication, Scientific Reports of International Space Science Institute (Bern, Switzerland). 2006. V. 5. P. 69-141.

84. InogamovN.A. The Role of RT and RM instabilities in astrophysics: an introduction // Astroph. and Space Phys. Rev. 1999. V. 10. pt. 2.

85. Jones T.W., RudnickL., JunB. et al. 1051 Ergs: The evolution of shell supernova remnants // Publ. of the Astron. Soc. of the Pacific. 1998. V. 110. P.125-151.

86. KafatosM. Time-dependent radiative cooling of a hot low-density cosmic gas // Astrophys. J. 1973. V. 182. P. 433-447.

87. Kane J., Arnett D., Remington B. A. et al. Two-dimensional versus three-dimensional supernova hydrodynamic instability growth // Astrophys. J., 2000. V. 528. P. 989-994.

88. Klein R.I., BudilK.S., Perry T.S. etal. The interaction of supernova remnants with interstellar clouds: experiments on the nova laser // Astrophys. J. 2003. V. 583. P. 245-259.

89. Klein R.I., McKee C.F., ColellaP. On the hydrodynnamic interaction of shock waves with interstellar clouds. I. Nonradiative shocks in small clouds// Astrophys. J. 1994. V.420. P. 213-236.

90. Lebedev S.V., Chittenden J.P., Beg F.N. et al. Laboratory astrophysics and collimated stellar outflows: the production of radiatively cooled hypersonic plasma jets // Astrophys. J. 2002. V. 564. P. 113-119.

91. LeflochB., LazareJJBCastetsA. Cometary globules. III. Triggered star formation in 1С 1848 // Astron. Astrophys. 1997. V. 324. P. 249-262.

92. Levenson N.A., Graham J.R. Environmental impact on the southeast limb of the Cygnus Loop // Astrophys. J. 2005. V. 622. P. 366-376.

93. McKee C. F., HollenbachD. J.; SeabG. C. et al. The structure of the time-dependent interstellar shocks and grain destruction in the interstellar medium // Astrophys. J. 1987. V. 318. P. 674-701.

94. Mellema G., Kurk J.D., Rottgering H.J. Evolution of clouds in radio galaxy cocoons // Astron. Astrophys. 2002. V.395. P.L13-L16.

95. MihalasD., Weibel-Mihalas B. Foundations of radiation hydrodynamics. NY: Dover publications, Inc., 1999. 718p.

96. MizutaA., Kane J. 0., PoundM.W. et al. Formation of pillars at the boundaries between HII regions and molecular clouds. // Astrophys. J. 2006. V. 647. P. 1151-1158.

97. Orlando S., Peres G., RealeF. et al. Crushing of interstellar gas clouds in supernova remnants I. The role of thermal conduction and radiative losses // Astron. Astrophys. 2005. V.444. P. 505-519.

98. OsterbrockD.E., Ferland G.J. Astrophysics of gaseous nebulae and active galactic nuclei. 2, Sausalito, C.A.: University Science Books, 2006. 480 p.

99. OttE. Nonlinear evolution of the Rayleigh-Taylor instability of a thin layer // Phys. Rev. Lett. 1972. V. 29. №21. P. 1429-1432.

100. Parker E. Dynamics of the interplanetary gas and magnetic // Astrophys. J. 1958. V. 128. №3. P. 664-676.

101. Patnaude D. J., FesenR. A., Raymond J. C., et al. An isolated, recently shocked ISM cloud in the Cygnus Loop supernova remnant // Astronomical J. 2002. V. 124. P. 2118-2134.

102. Poludnenko A.Y., DannenbergK.K., DrakeR.P. et al. A Laboratory investigation of supersonic clumpy flows: experimental design and theoretical analysis // Astrophys. J. 2004. V. 604. №1. P. 213-221.

103. Raymond J. C., CoxD.P., Smith B.W. Radiative cooling of a low-density plasma // Astroph. J. 1976. V. 204. P. 290-292.

104. Richtmyer R.D. Taylor instability in shock acceleration of compressible fluids // Communs. Pure and Appl.Math. 1960. V. 13. №2. P. 297-319.

105. Spitzer L.Jr. Acoustic waves in supernova remnants // Astrophys. J. 1982. V. 262. №1. P. 315-321.

106. Taylor J. The instability of liquid surfaces when accelerated in a direction perpendicular to their planes I // Proc. Roy Soc. London. Ser. A. 1950. V. 201. №1065. P. 192-196.

107. Walsh J. M., ShrefflerR.G., Willig F. G. Limiting conditon for jet formation in high velocity collisions. // J. Appl. Phys. 1953. V.24. №3. P. 349359; см. также сб. "Механика", вып. 4, Иностранная литература, 1954.

108. Woodward P.R. Shock-driven implosion of interstellar gas clouds and star formation // Astrophys. J. 1976. V. 207. pt. 1. P. 484-501.

109. Woodward P.R., ColellaP. The numerical simulation of two-dimensional fluid flow with strong shocks // J. Сотр. Phys. 1984. V. 54. P. 115-173.